- 1ELKH-EKKE Lendulet Environmental Microbiome Research Group, Eszterházy Károly Catholic University, Eger, Hungary
- 2Doctoral School of Environmental Sciences, Hungarian University of Agricultural and Life Sciences, Gödöllő, Hungary
- 3Food and Wine Research Institute, Eszterházy Károly Catholic University, Eger, Hungary
- 4Institute for Viticulture and Enology, Eszterházy Károly Catholic University, Eger, Hungary
The composition, diversity and dynamics of microbial communities associated with grapevines may be influenced by various environmental factors, including terroir, vintage, and season. Among these factors, terroir stands out as a unique possible determinant of the pathobiome, the community of plant-associated pathogens. This study employed high-throughput molecular techniques, including metabarcoding and network analysis, to investigate the compositional dynamics of grapevine fungal pathobiome across three microhabitats (soil, woody tissue, and bark) using the Furmint cultivar. Samples were collected during late winter and late summer in 2020 and 2021, across three distinct terroirs in Hungary’s Tokaj wine region. Of the 123 plant pathogenic genera found, Diplodia, Phaeomoniella, and Fusarium displayed the highest richness in bark, wood, and soil, respectively. Both richness and abundance exhibited significant disparities across microhabitats, with plant pathogenic fungi known to cause grapevine trunk diseases (GTDs) demonstrating highest richness and abundance in wood and bark samples, and non-GTD pathogens prevailed soil. Abundance and richness, however, followed distinct patterns Terroir accounted for a substantial portion of the variance in fungal community composition, ranging from 14.46 to 24.67%. Season and vintage also contributed to the variation, explaining 1.84 to 2.98% and 3.67 to 6.39% of the variance, respectively. Notably, significant compositional differences in fungi between healthy and diseased grapevines were only identified in wood and bark samples. Cooccurrence networks analysis, using both unweighted and weighted metrics, revealed intricate relationships among pathogenic fungal genera. This involved mostly positive associations, potentially suggesting synergism, and a few negative relationships, potentially suggesting antagonistic interactions. In essence, the observed differences among terroirs may stem from environmental filtering due to varied edaphic and mesoclimatic conditions. Temporal weather and vine management practices could explain seasonal and vintage fungal dynamics. This study provides insights into the compositional dynamics of grapevine fungal pathobiome across different microhabitats, terroirs, seasons, and health statuses. The findings emphasize the importance of considering network-based approaches in studying microbial communities and have implications for developing improved viticultural plant health strategies.
Introduction
The complex interplay between the grapevines and their environment extends beyond the mere physicochemical properties of the soil, it highly possibly resonates within shaping the composition, diversity, and dynamics of the microbial communities associated with grapevines, which play a crucial role in grapevine health and quality (Gladstones, 1992; Van Leeuwen and Seguin, 2006; OIV, 2010; Vaudour et al., 2015; Bokulich et al., 2016; Belda et al., 2017; Hunter et al., 2020; Pauvert et al., 2020; Gobbi et al., 2022). Among these environmental factors, terroir, vintage, and season stand out as assumed key determinants of the pathobiome composition and dynamics (Bokulich et al., 2014; Pauvert et al., 2020; Bekris et al., 2021).
Terroir, encompassing the unique combination of soil, climate, and topography of a vineyard, has a profound influence on the physiology of grapevines, potentially influencing the composition and interactions of the pathobiome within grapevine microhabitats (Hunter et al., 2020; Gobbi et al., 2022). However, the specific mechanisms by which terroir shapes the pathobiome remain largely unexplored. Vintage, the year in which grapes are harvested, is another environmental factor known to influence grapevine characteristics, potentially through its impact on the pathobiome (Bekris et al., 2021). By examining the microbial communities at different stages of the growing season, we can gain insights into the dynamics of the pathobiome and its potential influence on quality. Furthermore, season, the specific time of year during which grapes are grown and harvested, can further modulate the composition of the grapevine microbiome. The unique environmental conditions of each season, such as temperature and rainfall, can influence the abundance and diversity of microbial taxa within grapevine microhabitats (Pauvert et al., 2020).
The application of network-based approaches, particularly unweighted and weighted metric analysis, provides an avenue to discern the structure and complexity of these microbial communities. By quantifying the interactions among microbial taxa and assessing their relative importance within the network, this approach helps to elucidate the intricate interplay between microbial taxa and their environment (Herren and McMahon, 2017; Banerjee et al., 2018; Yousefi et al., 2023). Unweighted network analysis, based on presence/absence of taxa, provides insight into community structures and interactions that may be driven by environmental filters or competitive exclusion. In contrast, weighted networks, accounting for the abundance of different taxa, enable detection of more subtle and potentially important interactions that are often overlooked in unweighted analysis, including those influenced by the relative abundance of taxa (Barberán et al., 2012; Yousefi et al., 2023). This methodology allows for a deeper understanding of community structures and interactions that may be driven by environmental filters or competitive exclusion (Barberán et al., 2012; Banerjee et al., 2018; Yousefi et al., 2023).
Amidst the vast vineyards and microbial communities, the pathobiome plays a critical role in plant health by maintaining a balanced ecosystem that limits the growth and spread of pathogens (Vayssier-Taussat et al., 2014; Brader et al., 2017; Bass et al., 2019; Pauvert et al., 2020). Understanding the composition and dynamics of the pathobiome is essential for developing effective strategies to manage grapevine diseases and improve vineyard sustainability. Inside the grapevine’s pathobiome, the specter of Grapevine Trunk Diseases (GTDs) poses a significant threat. GTDs, often attributed to wood endophytes with the potential to switch to a pathogenic lifestyle, present a significant threat to the current vineyard sustainability (Sakalidis et al., 2011; Laurent et al., 2020; Bettenfeld et al., 2020). The prevalence of these diseases, and currently considered the most destructive biotic factor of grapevines worldwide, exacerbated by factors such as vineyard intensification and environmental changes, makes them a focal point in the examination of the grapevine pathobiome (Songy et al., 2019; Bettenfeld et al., 2020).
Despite previous studies on microbiome terroir and seasonality in grapevines (Bokulich et al., 2014; Gilbert et al., 2014; Knight et al., 2015; Zarraonaindia et al., 2015; Gobbi et al., 2020, 2022; Geiger et al., 2022), our understanding of the extent to which microbial communities differ among terroirs and the significance of these differences at varying spatial scales is still limited (Gobbi et al., 2022). Furthermore, there is increasing evidence indicating compositional differences among microbiomes of different grapevine parts or microhabitats. This suggests the role of strong niche-based processes in community assembly (Singh et al., 2018; Del Frari et al., 2019; Martínez-Diz et al., 2019; Swift et al., 2021; Geiger et al., 2022; Molnár et al., 2022). Yet, it remains unknown whether the effect of environmental differences on microbiome communities varies among these microhabitats. For example, it is still unclear whether soil or bark microbial communities, which are more exposed to the elements, exhibit greater differences among terroirs than microbial communities in living inner woody tissues, such as the xylem and phloem.
To address these knowledge gaps, our study embarks on an exploration of the composition, diversity, and dynamics of grapevine plant pathogenic fungi. We employed high-throughput molecular techniques, including metabarcoding and network analysis, to characterize the microbial communities associated with different grapevine microhabitats, and its possible influences by terroirs, vintage, seasons, and health status. Specifically, we are looking at grapevines that are symptomatic vs. asymptomatic for Esca-type grapevine trunk disease, all of the same cultivar (Furmint) and rootstock. Previous studies have investigated the microbiome of grapevines, but our study is the first to systematically compare the pathobiome composition and dynamics across different terroirs, seasons, and health statuses. Accordingly, we have formulated the following hypotheses: (1) The grapevine fungal pathobiome differs among different microhabitats, with distinct microbial communities inhabiting soil, bark, and wood tissue; (2) The diversity, abundance, and distribution of the pathobiome are influenced by different parameters, such as terroir, seasonality, vintage aspects, and plant health status. The findings from this study are expected to provide a more understanding of environmental factors that influence the diversity and distribution of pathogenic, offering valuable insights for better understanding of the pathobiome in viticulture. By employing network analysis, this study provides a more in-depth understanding of the structure and complexity of these microbial communities, which helped in uncovering microbial interactions and dynamics within these ecosystems.
These findings have significant implications for viticulture, as they provide a deeper understanding of the factors that influence the composition and dynamics of the grapevine pathobiome. This knowledge can be used to develop targeted interventions to control Grapevine Trunk Diseases and improve overall vineyard sustainability. Moreover, our results have broader implications for microbial ecology, offering a novel perspective on the influence of environmental factors on microbial communities.
Materials and methods
Study area and sampling
Samples were collected in three vineyards in the southern Tokaj wine region, Hungary: in the Szt. Tamás dűlő near Mád (48°11ʹ29.9”N 21°17ʹ39.0ʹE), and the southern (48°06ʹ27.1”N 21°21ʹ59.2ʹE) and northern (48°08ʹ46.3”N 21°22ʹ37.2ʹE) slopes of Nagy-Kopasz Hill near the town of Tokaj. These terroirs represent different mesoclimatic and edaphic conditions. In each vineyard, we selected eight grapevines; four of which had been marked as symptomatic, exhibiting typical leaf symptoms of the Esca type of Grapevine Trunk Disease (GTD), and four symptomless neighboring plants selected for each symptomatic grapevine. Samples were collected from the same marked plants at four sampling times: in February and August of 2020 and 2021.
All sampled grapevine plants were of the cultivar Furmint, an old and popular white variety endemic to the Carpathian Basin. This is the principal cultivar of the Tokaj wine region used to produce the botrytized Aszú wines, as well as dry white wines. The vineyards, aged between 26 to 44 years, featured the same scion (V. vinifera cv. Furmint clone T85) grafted on the same rootstock (V. berlandieri x riparia Teleki Kober 5BB). All were managed identically, with cordon training and the same conventional plant protection regime.
Three different microhabitats were sampled for each plant: bark tissue, perennial woody tissue, and bulk soil. The soil samples were taken from four sampling points within 50 cm around the trunk of the selected grapevine, to a depth of 20 cm. The samples were combined into a composite sample for each of the four grapevines per terroir and microhabitat. The bark was peeled with sterile scalpel at four points of the trunk, while woody tissue was sampled from underneath the four peeled areas of the trunk with a drill that was sterilized with a 4% (v/v) bleach solution between plants.
The samples were kept in a cooler in the field and, within five hours of collecting, were transported to the laboratory and stored at −80°C. A total of 288 composite samples were collected, representing the three microhabitats in four symptomatic and four asymptomatic plants per vineyard, sampled over two seasons across two years. All samples were lyophilized for a minimum of 72 h and homogenized using a tissue lyser and steel beads prior to DNA extraction.
Molecular work
Metagenomic DNA was extracted from approximately 0.5 mL of lyophilized and homogenized soil, bark, and woody tissues from each sample using the NucleoSpin® Plant and Soil DNA isolation kit (Macherey-Nagel Gmbh & Co., Düren, Germany). This extraction was performed according to the manufacturer’s protocol for both the tissue and soil samples. The extracted DNA samples were PCR amplified using primers fITS7 (Ihrmark et al., 2012) and ITS4 (White et al., 1990). The PCR conditions were as follows: one cycle of 95°C for 5 min; then 35 cycles of 95°C for 20 s, 56°C for 30 s, and 72°C for 1.5 min; and concluded with one cycle at 72°C for 7 min. Prior to PCR, the DNA samples were normalized for concentration.
To link the sequences to the sample source, a second PCR was performed using the same primers, which this time were equipped with Illumina adaptors and sample-specific, 8-nucleotide barcodes. The normalization of DNA extracts, PCR reactions, and 250-base paired-end Illumina NovaSeq sequencing of amplicon libraries were conducted at BaseClear B.V. (Leiden, the Netherlands) and at BIOMI Kft. (Gödöllő, Hungary).
Besides DNA metabarcoding, the inner wood samples have also been used to isolate fungi on agar media to obtain more information on the pathobiome related to GTDs. The entire ITS rDNA region of ca. 800 cultures has been sequenced and have also been detected by metabarcoding (Geiger et al., unpublished data).
Bioinformatic work
Raw DNA sequences were processed with the dada2 package (Callahan et al., 2016), which is implemented in R (R Core Team, 2020). This package is designed to resolve fine-scale DNA sequence variation and offers improved elimination of artifactual sequences. Since dada2 does not involve clustering sequences into Operational Taxonomic Units (OTUs), and is robust in removing spurious data, the output of unique Amplicon Sequence Variants (ASVs) captures both intra- and interspecific genetic variation of fungi found in the samples. This allows for the exploration of strain-level differences in inter- and intraspecific interactions.
Based on the quality score profiles, all forward and reverse reads were truncated to 250 and 200 bp, respectively, and were quality filtered, with the maximum number of expected errors (maxEE) allowed in a read set to 2. The filtered reads were denoised, the two-directional reads were merged, and clustered into sequence variants, which were later subjected to chimera filtering.
Taxonomic assignments of fungi were made based on the UNITE database of reference sequences, which represent all fungal Species Hypotheses (SHs) based on a dynamic delimitation (Kõljalg et al., 2013), using USEARCH v. 11 (Edgar, 2010). We assigned ASVs identified to at least the genus level to putative functional groups using the curated reference database of FungiTraits ver. 1.2 (Põlme et al., 2020).
In this paper, we investigate plant pathogenic fungi within the main genera, highlighting those known to be associated with GTDs, according to our former literature search (Geiger et al., 2022). All sequences of fungal ASVs analyzed in this paper have been deposited in GenBank (OQ370580-OQ371285).
Statistical analyses
All statistical analyses were performed in the R environment for statistical computing (R Core Team, 2020). The ASV table was normalized for subsequent statistical analyses by rarefying the number of high-quality fungal sequences through random subsampling of the smallest library (15,246 sequences) using the rrarefy function in the vegan R package (Oksanen et al., 2013). The resulting matrix contained 1,616,808 read counts of 5,858 fungal ASVs that served as input for the subsequent analyses.
Rarefied fungal read abundance and ASV richness were compared among samples according to categorical variables, i.e., microhabitat (bark, living wood, and bulk soil), terroir (North, South and Szt. Tamás), vintage (2020 and 2021), season (late winter and late summer) and health state (symptomatic and asymptomatic), using analysis of variance (ANOVA) and Tukey’s HSD test in R. ASV richness and rarefied read abundance of plant pathogens were graphically presented as boxplots using the ggplot2 R package (Wickham, 2016).
To visualize differences in fungal community composition among the samples, we used the vegan R package (Oksanen et al., 2013) to run non-metric multidimensional scaling (NMDS) with Bray-Curtis’s distance measure on the Hellinger-transformed ASV table, using the metaMDS function in vegan with 999 permutations in R, and a secondary matrix containing variables such as sampling source, vintage, season, health, and terroir. In addition, we performed PerMANOVA (adonis) using vegan to estimate the amount of variation explained by the variables, which were tested for significance independently and by including significant variables in a combined model that accounted for any correlations among them.
Climatic variables were measured using on-site weather stations at multiple terroirs in the Tokaj wine region. Temperature and precipitation changes during the two vintages (2020 and 2021) and seasons (Late Winter and Late Summer) were assessed via the stored climatic data at met.boreas2.hu/eke/ (Local Hungarian Meteorological Service, 2023). The data were subjected to permutational multivariate analysis of variance using the adonis function in the vegan R package (Oksanen et al., 2013). The resulting table (Table 1) illustrates the proportion of variation (%) in fungal community composition explained by different climatic variables in our seasons and vintages.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-t001.jpg)
Table 1. Proportion of variation (%) of the temperature and precipitation recorded of the area of the fungal community being explained by different seasonality (Late winter and Late summer), and year/vintage (2020 and 2021) calculated with permutational multivariate analysis of variance, based on the fungal community matrix.
We inferred and visualized possible interactions among ASVs of plant pathogenic fungi in all the three different sources as a cooccurrence table and network graph, using the cooccur and visNetwork packages, respectively (Veech, 2013; Griffith et al., 2016; Almende et al., 2019) based on presence-absence matrices. The cooccur package was used to select pairs of ASVs that co-occurred significantly more often or less often than expected by chance. We then used the co-occurrence table for network visualization, indicating both positive and negative pairwise correlations.
The weighted and unweighted network analysis was performed using the R package igraph (version 0.10.1, Csardi and Nepusz, 2006). It was performed descriptive and comparative statistics on network measures, where they were statistically analyzed to assess the effects of microhabitat, terroir, season, vintage, and health status. The statistical metrics analyzed: Average degree, network density and modularity (unweighted network analysis), and degree and betweenness centrality (weighted network analysis). See material and methods section of network metrics analysis for further details.
Furthermore, we performed indicator species analysis using the multipatt function in the indicspecies package (De Cáceres et al., 2012) with 999 permutations, to identify characteristic and differential fungal ASVs for the different categorical variables: microhabitat, vintage, season, terroir, and health type.
Network metrics analysis
The network analysis was conducted using the R package igraph (version 0.10.1, Csardi and Nepusz, 2006), which provided the toolkit for our network data analysis. Raw data were transformed into adjacency matrices for subsequent statistical analysis. These matrices represented the presence or absence of ASVs for the unweighted analysis, and the abundance of different ASVs for the weighted analysis.
The function ‘graph_from_adjacency_matrix()’ from the igraph package was employed to convert these matrices into network objects. For the weighted network analysis, the ‘weighted’ parameter was set to true to account for the abundance values. The Louvain method was applied to identify clusters or modules within the networks. This method is suitable for both weighted and unweighted networks. The ‘cluster_louvain()’ function from igraph was utilized to apply the Louvain algorithm, detecting communities in each subgraph and calculating a modularity score based on these communities.
A variety of network metrics were computed for each network (weighted and unweighted). These include the average degree (using the ‘degree()’ function), network density (using the ‘graph.density()’ function), and betweenness centrality (using the ‘betweenness()’ function). Modularity, a measure of the division strength of a network into modules, was obtained from the results of the Louvain community detection.
Custom functions and loops were designed and implemented in the code to facilitate the computation of network measures across multiple networks. These functions significantly reduced the amount of repetitive code, enhancing the efficiency and readability of the analysis. Each function was designed to accept a network object as input, perform the necessary calculation, and return the result. These functions were applied to each network in a loop, with the results compiled into a data frame for subsequent analysis.
Results
Among the 3,010 ASVs assigned to functional groups, 733 ASVs represented 123 plant pathogenic genera, including fifteen with known associations with GTDs. The five most dominant plant pathogenic genera were: Phaeomoniella (80 ASVs), Devriesia (36 ASVs), Fusarium (35 ASVs), Diplodia (29 ASVs), and Alternaria (25 ASVs). Within the three different microhabitats, Diplodia showed the highest richness in bark and woody tissues, followed by Phaeomoniella, and Devriesia. In soil, Fusarium and Alternaria were the most diverse. The distribution and sequence read abundance of all the plant pathogenic group genera among the samples are shown in Table 2. The entire ITS rDNA region of ca. 800 cultures have been sequenced and more than 90% of these have also been detected by metabarcoding (Geiger et al., unpublished data).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-t002.jpg)
Table 2. Distribution and rarefied sequence read abundance of all the plant pathogenic group genera among the samples, which were within different source types (wood, bark tissues, and soil bulk) of plant pathogenic fungal main community composition and its different terroir, season, year, and health.
The variance analyses for richness and abundance values of different plant pathogens across different microhabitats (sources), terroirs, vintage, season, and health types revealed notable distinctions. For plant pathogenic fungi, richness exhibited significant differences among microhabitats/sources (p = 7.64 × 10−8), vintage (p = 1.74 × 10−14), and season (p = 2.79 × 10−11). However, the terroir (p = 0.168) and health types (p = 0.242) did not significantly impact richness. Regarding rarefied fungal read abundance values of plant pathogenic fungi, significant differences were observed among different microhabitats (p = 2×10−16), terroirs (p = 0.000168), and health types (p = 0.033922). Meanwhile, the influence of vintage (p = 0.785) and season (p = 0.363) on abundance was not significant. Additionally, total richness values were significantly different for vintage (p = 2.60 × 10−14) and season (p = 2.32 × 10−11), whereas abundance values showed significant differences with terroirs (p = 7.86 × 10−4) and health types (p = 0.033922).
The fungal richness was greatest in soil, followed by bark and wood (Figure 1A). In terms of abundance values, wood showed by far the highest amount, followed by bark and soil (Figure 1B). Plant pathogens abundance considering health types was highest in symptomatic plants (Sympt.) compared to non-symptomatic (Asympt.) (Figure 1B). Another observation was the highest richness in August (late summer) compared to February (late winter) in the seasonality parameter, while observed richness in 2020 was higher than in 2021 (Figure 1A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-g001.jpg)
Figure 1. Comparison of plant pathogenic richness (A) and rarefied fungal read abundance (B) across the four-parameters and source types, which represented as Szt (Szt. Tamás terroir), y20 (year 2020), y21 (year 2021), Feb (February), and Aug (August). Means were compared using ANOVA and Tukey’s HSD tests, with letters denoting significant differences.
Fungal community composition was visualized in a two-dimensional NMDS ordination for all the plant pathogenic community (stress value = 0.07869204), as well as each source: wood (stress value = 0.08765602), bark (stress value = 0.09239610) and soil (stress value = 0.07753184) (Figure 2). PERMANOVA, using adonis, indicated that the composition was strongly correlated with terroirs, vintage, and season, in this order of strength, for all microhabitats, while significant compositional difference of fungi between asymptomatic and symptomatic grapevines was only observed in bark and inner woody tissue (Figure 3; Table 3). Although fungi from all three microhabitats showed compositional differences among terroirs, soil communities showed the greatest effect of terroir on compositional turnover, explaining 24.67% of variance (p = 0.0001), followed by bark (16.56%, p = 0.0001) and woody tissue (14.46%, p = 0.0001). The effect of vintage showed the same above pattern as terroir, although it was much weaker, while the influence of season was greatest on bark communities (2.98%, p = 0.0002), followed by soil (1.93%, p = 0.0036), and wood (1.84%, p = 0.0337). As for vintage, it has the same tendency as terroir. Health type only explained ca. 3% of the compositional variance in woody tissue and bark (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-g002.jpg)
Figure 2. Two-dimensional non-metric multidimensional scaling (NMDS) ordination illustrating the combined community structure of plant pathogenic fungi with ellipses representing the ordination scalability of the microhabitats: wood, bark, and soil. The overall stress value for this combined analysis is 0.07869204.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-g003.jpg)
Figure 3. Separate two-dimensional non-metric multidimensional scaling (NMDS) ordinations representing the community structure of plant pathogenic fungi in each microhabitat: wood, bark, and soil. Each graph demonstrates variations according to different parameters, namely terroir (North, South, and Szt. Tamás), health type (Asymptomatic and Symptomatic), vintage (20 for year 2020, and 21 for year 2021), and season (A = August/Late summer, and F = February/Late winter). Individual stress values for each microhabitat analysis are as follows: wood = 0.08765602, bark = 0.09239610, and soil = 0.07753184.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-t003.jpg)
Table 3. Proportion of variation (%) in different source types (wood, bark tissues and soil bulk) of fungal community composition explained by different terroir, season, year, and health calculated with permutational multivariate analysis of variance, based on the fungal community matrix.
The assessment of climatic variables, temperature and precipitation, displayed significant influence on fungal community composition (Table 1). Both seasonal and vintage parameters contributed to the variation in fungal communities. The late winter season exhibited a statistically significant impact on fungal community composition, explaining 1.188% of the variance (p = 0.0012) in temperature. Similar to late winter, the late summer season significantly contributed to compositional changes, accounting for 1.179% of the variance (p = 0.0011). This underscores the sensitivity of fungal communities to temperature fluctuations during the colder and warmer months. In contrast to temperature, precipitation during late winter did not show a significant impact on fungal communities (0.239%, p = 0.7607). Late summer precipitation, on the other hand, played a more substantial role, explaining 1.074% of the variance (p = 0.0031).
The year 2020 contributed significantly to the variance in fungal communities, with temperature explaining 1.09% (p = 0.0021) and precipitation explaining 2.215% (p = 0.0001). In 2021, the impact of temperature reduced to 0.646% (p = 0.0451), and precipitation explained 0.781% of the variance (p = 0.035). Although less pronounced, these results indicate continued but potentially attenuated seasonal effects on fungal communities in the subsequent year. Cooccurrence analyses and network graphs revealed differences in network complexity of plant pathogenic fungi among microhabitats in each terroir. The size of networks of significant relationships increased from woody tissue (265 ASVs) to bark (566 ASVs) and to soil (1,669 ASVs) (Figures 4–6; Supplementary Table S1). Positive relationships were more prevalent than negative ones in all microhabitats, as evidenced by the co-occurrence analysis (Figures 4–6). In soil (Figure 6), there were 1,407 positive and 262 negative relationships between pairs of ASVs, in bark (Figure 5) there were 414 positive and 152 negative relationships, and in wood (Figure 4), there were 195 positive and 70 negative relationships, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-g004.jpg)
Figure 4. Cooccurrence comparison visual network of the microhabitat wood visualized by the method ‘visNetwork’. Nodes represented by each specific ASVs with positive cooccurrence being represented by the blue lines and red dotted lines the negatives. A follow up table with specificities (descriptive statistics) can be seen in Supplementary Table S1. Legend specifies the main important genus represented by ASVs shown in the network.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-g005.jpg)
Figure 5. Visual network comparing the cooccurrence of microhabitat in Bark ASVs using the ‘visNetwork’ method. Nodes represent individual ASVs, with blue lines indicating positive cooccurrence and red dotted lines indicating negative cooccurrence. A table with descriptive statistics can be found in Supplementary Table S1. The legend highlights the main genus represented by ASVs in the network.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1322559/fmicb-14-1322559-HTML/image_m/fmicb-14-1322559-g006.jpg)
Figure 6. Network visualization comparing the co-occurrence of microhabitat soil ASVs using the ‘visNetwork’ method. Each node represents a specific ASV, and the blue lines indicate positive co-occurrence while the red dotted lines indicate negative co-occurrence. Supplementary Table S1 provides additional descriptive statistics. The legend indicates the key genus represented by ASVs in the network with rendered imaging for size fitting of the visualization.
The network metrics analysis elucidated the structure of microbial communities across different conditions, focusing on the specific metrics of betweenness and degree for the weighted analysis, and average degree, network density, and modularity for the unweighted analysis (Supplementary Tables S2, S3). Considering the abundance of different ASVs, the weighted network analysis showed from the weighted network metrics that for the bark microhabitat during late winter in vintage year ‘21’, the southern terroir showed the highest degree (3.662650602) for asymptomatic health status. Interestingly, the bark microhabitat in the northern terroir demonstrated the highest betweenness (147.316092) for symptomatic health status (Supplementary Table S2).
However, the soil microhabitat during late summer of vintage year ‘20’, the southern terroir showed the highest degree (4.418918919) for symptomatic health status, while the soil northern terroir during late summer of the same year displayed the highest betweenness (150.5374717) for symptomatic health status (Supplementary Table S2). As for the wood microhabitat during late winter in vintage year ‘21’, the southern terroir exhibited the highest degree (3.672727273) for symptomatic health status, while the northern terroir, late winter of the same year, demonstrated the highest betweenness (120.4482759) (Supplementary Table S2) for symptomatic health status.
The unweighted network analysis, considering the presence or absence of different ASVs, revealed for the bark microhabitat during late winter of vintage year ‘21’, the southern terroir demonstrated the highest average degree (3.69047619) with a network density of 0.044463569 and a modularity of 0.266118626 for asymptomatic health status (Supplementary Table S3). In the soil microhabitat during late winter in vintage year ‘21’, the southern terroir had the highest network density (0.071969697) with an average degree of 2.303030303 and a modularity of 0.457063712 for asymptomatic health status. For the wood microhabitat during late winter in vintage year ‘21’, the southern terroir showed the highest average degree (3.593220339) with a network density of 0.061952075 and a modularity of 0.268823425 for symptomatic health status (Supplementary Table S3).
Collectively, the network structures of the fungal pathobiome differ significantly among the grapevine microhabitats, with intricate and robust network structures evident in certain conditions.
The number of indicator species was the highest in the soil, representing the genera Fusarium, Lactera, Alternaria, Coniothyrium, Dactylonectria, Plectosphaerella, Botryosphaeria etc. (Supplementary Table S4). The most significant indicators of terroir in soil were Lectera (North), Alternaria (South), and Coniothyrium (Szt. Tamás). For season, the most significant indicator was Fusarium (late summer), while none of the ASVs was an indicator for late winter. Lastly, the strongest indicators for vintage in soil were Dactylonectria (2020) and Plectosphaerella (2021). In bark samples, Diplodia, Microstroma, Devriesia, Botrytis was the most significant identified as indicators among the parameters, with an especially high presence of the genus Diplodia. In wood, significant indicators mostly belonged to grapevine trunk disease, such as Phaeomoninella, Diplodia, and Eutypella. From those, the genus Diplodia and Eutypella were strongly significant to the symptomatic of the health type parameter, and Devriesia and Phaeomoniella to the asymptomatic (Supplementary Table S2). Interestingly, in wood terroir the genus Pleurophoma and Diplodia were most significant for North and South, but in Szt. Tamás there were only indicators of the genus Phaeomoniella.
Discussion
This study delves into the intricate structure and dynamics of the grapevine pathobiome, primarily focusing on how microhabitat, terroir, vintage, season, and health status shape these complex fungal communities. In alignment with our first hypothesis, our findings illustrate a significant variation in richness and composition of plant pathogenic fungi among different microhabitats. Specifically, a decrease in the richness of plant pathogenic fungi was observed from soil to bark and wood samples, while the richness of pathogens associated with GTDs followed an increasing trend. This reinforces the concept of niche specialization within these fungal communities, demonstrating that certain microhabitats favor the proliferation of specific fungal groups. In addition, our findings suggest that selection pressure of abiotic environmental factors is buffered to a certain extent inside the grapevine trunk, as indicated by the somewhat weaker, although still significant, effects of terroir, season, and vintage on fungal community composition in wood, compared to bark and soil communities.
The soil, known to be greatly influenced by both mesoclimatic and edaphic factors worldwide (Tedersoo et al., 2014; Větrovský et al., 2019) as well as in the wider region of study in northeastern Hungary (Geml, 2019; Geml et al., 2022), exhibited the highest fungal richness, particularly during the late summer.
Adding depth to our understanding, the indicator species analysis shed light on specific fungal taxa serving as microbial fingerprints across various microhabitats. For instance, in the soil, where NMDS highlighted vintage influences, we found Fusarium indicative of late summer and terroir-specific markers such as Lactera (North) and Alternaria (South). Bark samples pointed toward Diplodia’s prominence, and in the wood, disease-associated genera like Phaeomoninella and Diplodia stood out, especially in relation to plant health status. These indicator species, in harmony with NMDS findings, emphasize the intricate connections between fungal communities, their microhabitats, and external factors.
The isolation of fungi from inner wood samples has provided additional insights into the pathobiome related to GTDs. More than 90% of the fungal genera detected through metabarcoding were also isolated from inner wood samples, confirming their presence within the grapevine (Geiger et al., unpublished data). This suggests that these genera are well-established members of the grapevine pathobiome and play a significant role in grapevine health and disease development.
Building upon this, it’s important to delve deeper into the nuances of vintage variations. In terms of vintage, distinct patterns were observed across different microhabitats. The NMDS ordination (Figure 3) illustrates a clearer separation of fungal communities between the years 2020 and 2021 in the soil samples, indicating a notable impact of vintage in this microhabitat. The year 2020 exhibited higher fungal richness compared to 2021. This difference could be attributed to the distinct climatic conditions each year presented, with 2020 having an early and extended rainy period, whereas 2021 experienced a severe heatwave and drought, as reported by the Copernicus Climate Change Service (C3S) and the European State of Climate (ESOTC).
These adverse weather conditions created significant mesoclimatic differences and could have led to a temporary decrease in overall plant pathogenic fungal diversity. This likely is the result of a selection pressure favoring taxa that can withstand low moisture content and adapt to shifts in the diffusion rates of nutrients and microbial signals. Many of these microbial signals, such as mycotoxins, volatile organic compounds (VOCs), and extracellular enzymes, are influenced by environmental moisture content (Henry et al., 2007; Smith and De Smet, 2012; Song et al., 2012; Barnard et al., 2013; Bokulich et al., 2014; Naylor et al., 2017; Santos-Medellín et al., 2017; Naylor and Coleman-Derr, 2018). Under drier conditions, the diffusion of these signals in the soil might be hindered, affecting both fungal-fungal and plant-fungal interactions. This can modify competition dynamics, nutrient acquisition strategies, and even the establishment of symbiotic relationships.
Nevertheless, the highest rarefied fungal read abundance of plant pathogens was observed in wood and in symptomatic plants, underscoring that both microhabitat and health status exert a significant influence on the diversity and abundance of fungal communities. Our second hypothesis suggested that the diversity, abundance, and distribution of the pathobiome are influenced by different parameters, such as terroir, seasonality, and vintage aspects. In support of this hypothesis, we found that compositional shifts in these fungal communities were chiefly influenced by site-specific environmental factors, often referred to as “microbial terroir.” These observations align with previous research indicating that both mesoclimatic and edaphic factors substantially shape fungal communities worldwide (Tedersoo et al., 2014; Větrovský et al., 2019). It becomes increasingly apparent that unraveling the intricate relationships between different microhabitats and their unique ‘microbial terroir’ characteristics is essential for comprehending the nuanced composition of the pathobiome.
Moreover, our investigation revealed a notable influence of vintage and seasonality on the pathobiome. In terms of richness and abundance, the soil and late summer season presented the highest fungal richness, while wood and symptomatic plants exhibited the highest abundance. These patterns suggest that both environmental conditions and plant health status can significantly influence the diversity and abundance of fungal communities. The NMDS ordination further substantiated these patterns, revealing distinct clustering of fungal communities based on terroir, vintage, and season, thereby affirming a strong influence of these factors on community composition.
The integration of climatic variables, temperature, and precipitation into our analysis provides a nuanced understanding of their impact on fungal community composition across seasons and vintages. Late winter and late summer exerted distinct influences, with temperature changes during these seasons significantly contributing to the variance. Interestingly, late summer precipitation played a more substantial role than its winter counterpart, underlining the sensitivity of fungal communities to warmer and drier conditions. Vintage variations added another layer to this narrative, with the year 2020 significantly shaping fungal communities through both temperature and precipitation, while 2021 exhibited a more tempered effect. These findings enhance our insights of the pathobiome dynamics, affirming that while it is evident that the weather of the actual year and season can influence microbial community dynamics, the mesoclimate of the particular terroirs plays a more significant role in shaping the overall diversity and composition of the microbial communities. These climatological variations were further depicted in Supplementary Figure S1 and Table 1 from Local Hungarian Meteorological Service (2023).
The uniformity of vineyard management practices, including consistent fungicide application, across all three terroirs during various seasons and vintages is a salient aspect of our study. While this uniformity minimizes direct impacts on the observed pathobiome variations, we acknowledge the potential indirect influence of fungicides on seasonal dynamics. Fungicide application introduces a dynamic element, influencing fungal community abundance and diversity, with documented short-term and long-term effects on non-target fungi (Leroux et al., 2002; Rodriguez-Morelos et al., 2021). Selective pressures on fungal populations, a known phenomenon in agricultural systems (Tranel and Wright, 2002), can result from fungicide applications, especially when timed with disease prevalence and seasonal factors. This aligns with findings in agricultural and environmental microbiology, highlighting fungicides’ capacity to shape the structure and diversity of microbial communities (Fournier et al., 2020; Zhang et al., 2020; Ma et al., 2021; Sun et al., 2023). While our study does not directly measure fungicide residues or their specific impact on fungal communities, recognizing fungicide application as a potential contributor to seasonal changes across terroirs suggests avenues for future research. Investigating the specific impacts and temporal dynamics of different fungicides on the pathobiome would enhance our understanding of factors shaping microbial communities in viticultural ecosystems.
Moreover, our network analysis suggested potential interactions, based on co-occurrence patterns, within these fungal communities. The unweighted analysis revealed differences in average degree, network density, and modularity across different parameters, pointing to the complexity and diversity of these communities. Bark samples exhibited higher average degree and modularity compared to soil and wood samples, possibly indicating a more complex network structure in this microhabitat due to its role as a critical overwintering site for various fungi and animals (Geiger et al., 2022).
Similarly, the weighted network analysis emphasized the role of the abundance data, revealing the importance of certain fungal genera such as Diplodia and Eutypella, as central ‘hubs’ in the network. These key players, often associated with GTDs, underscore the potential influence of disease dynamics on community assembly and vice versa.
To better understand these findings, it is crucial to interpret them in light of network metrics. The ‘degree’ of a node, representing the number of its connections, reflects its potential influence or importance within the network (Barabási and Oltvai, 2004). In our study, higher degrees in certain genera indicate their pivotal role in community interactions. Network density, which measures the proportion of actual connections relative to all possible connections, provides insights into the overall interconnectedness of the community (Newman, 2003). A higher density in our findings could imply a more integrated fungal community. Furthermore, ‘modularity’ assesses the degree to which the network can be partitioned into smaller, tightly connected groups or modules (Newman and Girvan, 2004). The higher modularity observed in bark samples suggests the existence of distinct microbial sub-communities or groups within that microhabitat. Biologically, more modular networks can contribute to system resilience, as disturbances might be confined within a single module, preventing them from affecting the entire system (Stouffer and Bascompte, 2011).
These observations support the notion that microbial interactions, which are effectively mapped through such network metrics, play a crucial role in the development and severity of grapevine diseases for this studied area. These interactions might involve synergistic or antagonistic relationships between co-infecting microbes (Hogan, 2006; Hosni et al., 2011; Lamichhane and Venturi, 2015).
Our findings offer valuable insights into the complex interplay between environmental factors and fungal communities in grapevines, supporting both our initial hypotheses. The metabarcoding, bioinformatics and network metrics suite of analyses conducted in this study enhances our understanding of the complex plant pathogenic fungal communities within grapevines and the environmental factors influencing them. Recognizing that the grapevine physiology and environment can strongly affect fungal behavior is crucial (Bruno and Sparapano, 2006; Freitas et al., 2009). By identifying relevant key players, understanding the possible effects of different parameters on fungal communities, and recognizing potential early indicators of disease, we are better equipped to devise effective strategies for managing impactful plant diseases in the fields. This knowledge can ultimately contribute to improved grapevine health and productivity. Building upon our research, there’s potential to not only advance the scientific understanding but also revolutionize grapevine cultivation and disease management strategies in the broader context of viticulture, creating a foundation for future studies, innovations and sustainable practices.
Looking forward, our study lays the groundwork for future investigations that could enhance our uncovering of the grapevine pathobiome. Firstly, a more in-depth exploration into the nuanced impact of vintage variations on fungal communities across different microhabitats would contribute significantly. Understanding how climatic conditions, especially during exceptional years, shape microbial interactions and disease dynamics could unravel critical insights. Additionally, delving into the direct impacts of fungicides, particularly their residues, on fungal communities could provide a more comprehensive understanding of their role in shaping seasonal variations. The network analysis opens avenues for studying microbial interactions further, emphasizing the need for experimental validations to decipher the nature of relationships indicated by network metrics. Investigating the specific responses of key genera, such as Diplodia and Eutypella, to environmental changes and their role in disease development could provide targeted strategies for disease management. Overall, this study serves as a steppingstone, and future research endeavors should capitalize on these insights to propel advancements in viticulture, offering innovative solutions for sustainable grapevine cultivation and disease control.
Conclusion
This investigation brings to light the multifaceted dynamics of grapevine pathobiomes across varied microhabitats. It accentuates the substantial impact of both environmental and temporal determinants on these dynamics. Viewing this through the concept of “microbial terroir” lets us start to comprehend the critical role that spatial differences may play in molding the microbial environment surrounding grapevines. Our research, bolstered by advanced network-based approaches, which played a pivotal role in our analytical methodology, offers deeper insights into grapevine-microbe interactions for the specific area of study. This understanding is valuable, given its broader implications in wider agricultural and ecological spheres. Furthermore, our findings can serve as a guide to refine sustainable viticulture practices, with a potential focus on tailored plant health management. While this study stands as a milestone in demystifying the intricate relationship between grapevines and their resident microbes, specifically in Hungary, it also beckons toward a future rich with possibilities for even more detailed exploration in the realm of pathobiome interplay.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) are GenBank; accession numbers OQ370580-OQ371285.
Author contributions
CL: Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Visualization, Writing – original draft, Writing – review & editing, Conceptualization, Validation. AG: Writing – review & editing, Data curation, Investigation, Validation. AM: Writing – review & editing, Investigation. KV: Resources, Writing – review & editing. GK: Writing – review & editing, Investigation. ZZ: Validation, Writing – review & editing, Investigation, Resources. JG: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing, Resources, Validation.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was financially supported by the Lendület Program (award no. 96049) of the Hungarian Academy of Sciences and the Eötvös Lóránd Research Network, awarded to JG, Stipendium Hungaricum PhD. scholarships (Tempus Foundation and Hungarian University of Agriculture and Life Sciences) to CL and GK.
Acknowledgments
We thank the owners of the Szepsy Winery for providing access and logistical support to sample grapevines in their vineyard for our study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1322559/full#supplementary-material
References
Almende, B. V., Thieurmel, B., and Robert, T. (2019). visNetwork: network visualization using ‘Vis.Js’ library. R package version 2.0.8. Available at: https://CRAN.R-project.org/package=visNetwork
Banerjee, S., Schlaeppi, K., and van der Heijden, M. G. (2018). Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol. 16, 567–576. doi: 10.1038/s41579-018-0024-1
Barabási, A. L., and Oltvai, Z. N. (2004). Network biology: understanding the cell's functional organization. Nat. Rev. Genet. 5, 101–113. doi: 10.1038/nrg1272
Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2011.119
Barnard, R. L., Osborne, C. A., and Firestone, M. K. (2013). Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 7, 2229–2241. doi: 10.1038/ismej.2013.104
Bass, D., Stentiford, G. D., Wang, H.-C., Koskella, B., and Tyler, C. R. (2019). The Pathobiome in animal and plant diseases. Trends Ecol. Evol. 34, 996–1008. doi: 10.1016/j.tree.2019.07.012
Bekris, F., Vasileiadis, S., Papadopoulou, E., Samaras, A., Testempasis, S., Gkizi, D., et al. (2021). Grapevine wood microbiome analysis identifies key fungal pathogens and potential interactions with the bacterial community implicated in grapevine trunk disease appearance. Environ. Microbiome 16:23. doi: 10.1186/s40793-021-00390-1
Bettenfeld, P., Fontaine, F., Trouvelot, S., Fernandez, O., and Courty, P.-E. (2020). Woody Plant Declines. What’s Wrong with the Microbiome? Trends in Plant Science. doi: 10.1016/j.tplants.2019.12.024
Belda, I., Zarraonaindia, I., Perisin, M., Palacios, A., and Acedo, A. (2017). From vineyard soil to wine fermentation: microbiome approximations to explain the ‘terroir’ concept. Front. Microbiol. 8, 1–12. doi: 10.3389/fmicb.2017.00821
Bokulich, N. A., Collins, T. S., Masarweh, C., Allen, G., Heymann, H., Ebeler, S. E., et al. (2016). Associations among wine grape microbiome, metabolome, and fermentation behavior suggest microbial contribution to regional wine characteristics. MBio 7, e00631–e00616. doi: 10.1128/mBio.00631-16
Bokulich, N. A., Thorngate, J. H., Richardson, P. M., and Mills, D. A. (2014). Microbial biogeography of wine grapes is conditioned by cultivar, vintage, and climate. Proc. Natl. Acad. Sci. 111, E139–E148. doi: 10.1073/pnas.1317377110
Brader, G., Compant, S., Vescio, K., Mitter, B., Trognitz, F., Ma, L.-J., et al. (2017). Ecology and genomic insights on plant-pathogenic and -nonpathogenic endophytes. Annu. Rev. Phytopathol. 55, 61–83. doi: 10.1146/annurev-phyto-080516-035641
Bruno, G., and Sparapano, L. (2006). Effects of three esca-associated fungi on Vitis vinifera L.: III. Enzymes produced by the pathogens and their role in fungus-to-plant or in fungus-to-fungus interactions. Physiol. Mol. Plant Pathol. 69, 182–194. doi: 10.1016/j.pmpp.2007.04.006
Callahan, B. J., McMurdie, P. J., and Rosen, M. J. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Csardi, G., and Nepusz, T. (2006). The igraph software package for complex network research. InterJournal, Complex Systems, 1695. Available at: http://igraph.org
De Cáceres, M., Legendre, P., Wiser, S. K., and Brotons, L. (2012). Using species combinations in Indicator value analyses. Methods Ecol. Evol. 3, 973–982. doi: 10.1111/j.2041-210X.2012.00246.x
Del Frari, G., Gobbi, A., Aggerbeck, M. R., Oliveira, H., Hansen, L. H., and Ferreira, R. B. (2019). Characterization of the wood mycobiome of Vitis vinifera in a vineyard affected by esca. Spatial distribution of fungal communities and their putative relation with leaf symptoms. Front. Plant Sci. 10:1646. doi: 10.3389/fpls.2019.00910
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Fournier, B., Santos, S., Gustavsen, J. A., Imfeld, G., Lamy, F., Mitchell, E. A. D., et al. (2020). Impact of a synthetic fungicide (fosetyl-Al and propamocarb-hydrochloride) and a biopesticide (Clonostachys rosea) on soil bacterial, fungal, and protist communities. Sci. Total Environ. 738:139635. doi: 10.1016/j.scitotenv.2020.139635
Freitas, R., Rego, C., Oliveira, H., and Ferreira, R. B. (2009). Interactions among grapevine disease-causing fungi. The role of reactive oxygen species. Phytopathol. Mediterr. 48, 117–127.
Geiger, A., Karácsony, Z., Golen, R., Váczy, K. Z., and Geml, J. (2022). The compositional turnover of grapevine-associated plant pathogenic fungal communities is greater among Intraindividual microhabitats and terroirs than among healthy and esca-diseased plants. Phytopathology 112, 1029–1035. doi: 10.1094/PHYTO-05-21-0190-R
Geml, J. (2019). Soil fungal communities reflect aspect-driven environmental structuring and vegetation types in a Pannonian forest landscape. Fungal Ecol. 39, 63–79. doi: 10.1016/j.funeco.2018.12.005
Geml, J., Leal, C. M., Nagy, R., and Sulyok, J. (2022). Abiotic environmental factors drive the diversity, compositional dynamics and habitat preference of ectomycorrhizal fungi in Pannonian forest types. Front. Microbiol. 13:1007935. doi: 10.3389/fmicb.2022.1007935
Gilbert, J. A., van der Lelie, D., and Zarraonaindia, I. (2014). Microbial terroir for wine grapes. Proc. Natl. Acad. Sci. 111, 5–6. doi: 10.1073/pnas.1320471110
Gobbi, A., Acedo, A., Imam, N., Santini, R. G., Ortiz-Álvarez, R., Ellegaard-Jensen, L., et al. (2022). A global microbiome survey of vineyard soils highlights the microbial dimension of viticultural terroirs. Commun. Biol. 5:241. doi: 10.1038/s42003-022-03202-5
Gobbi, A., Kyrkou, I., Filippi, E., Ellegaard-Jensen, L., and Hansen, L. H. (2020). Seasonal epiphytic microbial dynamics on grapevine leaves under biocontrol and copper fungicide treatments. Sci. Rep. 10, 1–13. doi: 10.1038/s41598-019-56741-z
Griffith, D. M., Veech, J. A., and Marsh, C. J. (2016). Cooccur: probabilistic species co-occurrence analysis in R. J. Stat. Softw. 69, 1–17. doi: 10.18637/jss.v069.c02
Henry, A., Doucette, W., Norton, J., and Bugbee, B. (2007). Changes in crested wheatgrass root exudation caused by flood, drought, and nutrient stress. J. Environ. Qual. 36, 904–912. doi: 10.2134/jeq2006.0425sc
Herren, C., and McMahon, K. (2017). Cohesion: a method for quantifying the connectivity of microbial communities. ISME J. 11, 2426–2438. doi: 10.1038/ismej.2017.91
Hogan, D. A. (2006). Talking to themselves: autoregulation and quorum sensing in fungi. Eukaryot. Cell 5, 613–619. doi: 10.1128/EC.5.4.613-619.2006
Hosni, T., Moretti, C., Devescovi, G., Suarez-Moreno, Z. R., Fatmi, M. B., Guarnaccia, C., et al. (2011). Sharing of quorum-sensing signals and role of interspecies communities in a bacterial plant disease. ISME J. 5, 1857–1870. doi: 10.1038/ismej.2011.65
Hunter, J. K., Tarricone, L., Volschenk, C., Giacalone, C., Melo, M. S., and Zorer, R. (2020). Grapevine physiological response to row orientation-induced spatial radiation and microclimate changes. OENO One 54, 411–433. doi: 10.20870/oeno-one.2020.54.2.3100
Ihrmark, K., Bödeker, I. T. M., Cruz-Martinez, K., Friberg, H., Kubartova, A., Schenk, J., et al. (2012). New primers to amplify the fungal ITS2 region – evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 82, 666–677. doi: 10.1111/j.1574-6941.2012.01437.x
Knight, S., Klaere, S., Fedrizzi, B., and Goddard, M. R. (2015). Regional microbial signatures positively correlate with differential wine phenotypes: evidence for a microbial aspect to terroir. Sci. Rep. 5:14233. doi: 10.1038/srep14233
Kõljalg, U., Nilsson, R. H., Abarenkov, K., Tedersoo, L., Taylor, A. F. S., Bahram, M., et al. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277. doi: 10.1111/mec.12481
Laurent, B., Marchand, M., Chancerel, E., Saint-Jean, G., Capdevielle, X., Poeydebat, C., et al. (2020). A Richer Community of Botryosphaeriaceae Within a Less Diverse Community of Fungal Endophytes in Grapevines than in Adjacent Forest Trees Revealed by a Mixed Metabarcoding Strategy. Phytobiomes Journal 4:252–267.
Lamichhane, J. R., and Venturi, V. (2015). Synergisms between microbial pathogens in plant disease complexes: a growing trend. Front. Plant Sci. 6:385. doi: 10.3389/fpls.2015.00385
Leroux, P., Fritz, R., Debieu, D., Albertini, C., Lanen, C., Bach, J., et al. (2002). Mechanisms of resistance to fungicides in field strains of Botrytis cinerea. Pest Manag. Sci. 58, 876–888. doi: 10.1002/ps.566
Local Hungarian Meteorological Service . (2023). Climatographical data retrieved from the Meteorological Stations. Available at: http://met.boreas2.hu/eke/
Ma, G., Gao, X., Nan, J., Zhang, T., Xie, X., and Cai, Q. (2021). Fungicides alter the distribution and diversity of bacterial and fungal communities in ginseng fields. Bioengineered 12, 8043–8056. doi: 10.1080/21655979.2021.1982277
Martínez-Diz, M. D. P., Andrés-Sodupe, M., Bujanda, R., Díaz-Losada, E., Eichmeier, A., and Gramaje, D. (2019). Soil-plant compartments affect fungal microbiome diversity and composition in grapevine. Fungal Ecol. 41, 234–244. doi: 10.1016/j.funeco.2019.07.003
Molnár, A., Geml, J., Geiger, A., Leal, C. M., Kgobe, G., Tóth, A. M., et al. (2022). Exploring relationships among grapevine chemical and physiological parameters and leaf and berry Mycobiome composition. Plan. Theory 11:1924. doi: 10.3390/plants11151924
Naylor, D., and Coleman-Derr, D. (2018). Drought stress and root-associated bacterial communities. Front. Plant Sci. 8:2223. doi: 10.3389/fpls.2017.02223
Naylor, D., DeGraaf, S., Purdom, E., and Coleman-Derr, D. (2017). Drought and host selection influence bacterial community dynamics in the grass root microbiome. ISME J. 11, 2691–2704. doi: 10.1038/ismej.2017.118
Newman, M. E. (2003). The structure and function of complex networks. SIAM Rev. 45, 167–256. doi: 10.1137/S003614450342480
Newman, M. E., and Girvan, M. (2004). Finding and evaluating community structure in networks. Phys. Rev. E 69:026113. doi: 10.1103/PhysRevE.69.026113
OIV . (2010). Definition of vitivinicultural “terroir”. The International Organization of Vine and Wine, Resolution OIV/VITI 333/2010.
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., et al. (2013). Vegan: community ecology package. Version 2.0-5. Available at: http://vegan.rforge.r-project.org/
Põlme, S., Abarenkov, K., Henrik Nilsson, R., Lindahl, B. D., Clemmensen, K. E., Kauserud, H., et al. (2020). FungalTraits: a user-friendly traits database of fungi and fungus-like stramenopiles. Fungal Divers. 105, 1–16. doi: 10.1007/s13225-020-00466-2
Pauvert, C., Fort, T., Calonnec, A., Faivre d’Arcier, J., Chancerel, E., Massot, M., et al. (2020). Microbial association networks give relevant insights into plant pathobiomes. bioRxiv, 2020.02.21.958033. doi: 10.1101/2020.02.21.958033
R Core Team . (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing: Vienna, Austria.
Rodriguez-Morelos, V. H., Calonne-Salmon, M., Bremhorst, V., Garcés-Ruiz, M., and Declerck, S. (2021). Fungicides with contrasting mode of action differentially affect hyphal healing mechanism in Gigaspora sp. and Rhizophagus irregularis. Front. Plant Sci. 12:642094. doi: 10.3389/fpls.2021.642094
Sakalidis, M. L., Hardy, G. E. S., and Burgess, T. I. (2011). Class III endophytes, clandestine movement amongst hosts and habitats and their potential for disease; a focus on Neofusicoccum australe. Australasian Plant Pathology, 40, 510–521. doi: 10.1007/s13313-011-0077-3
Santos-Medellín, C., Edwards, J., Liechty, Z., Nguyen, B., and Sundaresan, V. (2017). Drought stress results in a compartment-specific restructuring of the rice root-associated microbiomes. MBio 8, e00764–e00717. doi: 10.1128/mBio.00764-17
Singh, P., Santoni, S., This, P., and Péros, J.-P. (2018). Genotype-environment interaction shapes the microbial assemblage in grapevine's phyllosphere and carposphere: an NGS approach. Microorganisms 6:96. doi: 10.3390/microorganisms6040096
Smith, S., and De Smet, I. (2012). Root system architecture: insights from Arabidopsis and cereal crops. Philos. Trans. Royal Society B 367, 1441–1452. doi: 10.1098/rstb.2011.0234
Songy, A., Fernandez, O., and Clément, C. (2019). Grapevine trunk diseases under thermal and water stresses. Planta 249, 1655–1679. doi: 10.1007/s00425-019-03111-8
Song, F., Han, X., Zhu, X., and Herbert, S. J. (2012). Response to water stress of soil enzymes and root exudates from drought and non-drought tolerant corn hybrids at different growth stages. Can. J. Soil Sci. 92, 501–507. doi: 10.4141/cjss2010-057
Stouffer, D. B., and Bascompte, J. (2011). Compartmentalization increases food-web persistence. Proc. Natl. Acad. Sci. 108, 3648–3652. doi: 10.1073/pnas.1014353108
Sun, M., Wang, H., Shi, C., Li, J., Cai, L., Xiang, L., et al. (2023). Effect of azoxystrobin on tobacco leaf microbial composition and diversity. Front. Plant Sci. 13:1101039. doi: 10.3389/fpls.2022.1101039
Swift, J. F., Hall, M. E., Harris, Z. N., Kwasniewski, M. T., and Miller, A. J. (2021). Grapevine microbiota reflect diversity among compartments and complex interactions within and among root and shoot systems. Microorganisms 9:92. doi: 10.3390/microorganisms9010092
Tedersoo, L., Bahram, M., Polme, S., Koljag, U., Yorou, N. S., Wijensundera, R., et al. (2014). Global diversity and geography of soil fungi. Science 346:1256688. doi: 10.1126/science.1256688
Tranel, P. J., and Wright, T. R. (2002). Resistance of weeds to ALS-inhibiting herbicides: what have we learned? Weed Sci. 50, 700–712. doi: 10.1614/0043-1745(2002)050[0700:RROWTA]2.0.CO;2
Van Leeuwen, C., and Seguin, G. (2006). The concept of terroir in viticulture. J. Wine Res. 17, 1–10. doi: 10.1080/09571260600633135
Vaudour, E., Costantini, E., Jones, G., and Mocali, S. (2015). An overview of the recent approaches to terroir functional modelling, footprinting and zoning. Soil 1, 287–312. doi: 10.5194/soil-1-287-2015
Vayssier-Taussat, M., Albina, E., Citti, C., Cosson, J.-F., Jacques, M.-A., Lebrun, M.-H., et al. (2014). Shifting the paradigm from pathogens to pathobiome: new concepts in the light of meta-omics. Front. Cell. Infect. Microbiol. 4:29. doi: 10.3389/fcimb.2014.00029
Veech, J. A. (2013). A probabilistic model for analyzing species co-occurrence. Glob. Ecol. Biogeogr. 22, 252–260. doi: 10.1111/j.1466-8238.2012.00789.x
Větrovský, T., Kohout, P., Kopecký, M., Machac, A., Man, M., Bahnmann, B. D., et al. (2019). A meta-analysis of global fungal distribution reveals climate-driven patterns. Nat. Commun. 10:5142. doi: 10.1038/s41467-019-13164-8
White, T. M., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA for phylogenetics” in PCR protocols: A guide to methods and applications. eds. M. A. Innis, D. H. Gelfand, J. J. Sninsky, and T. J. White (Cambridge, MA: Academic Press), 315–321.
Yousefi, B., Melograna, F., Galazzo, G., van Best, N., Mommers, M., Penders, J., et al. (2023). Capturing the dynamics of microbial interactions through individual-specific networks. Front. Microbiol. 14:1170391. doi: 10.3389/fmicb.2023.1170391
Zarraonaindia, I., Owens, S. M., Weisenhorn, P., West, K., Hampton-Marcell, J., Lax, S., et al. (2015). The soil microbiome influences grapevine-associated microbiota. mBio 6, e02527–e02514. doi: 10.1128/mBio.02527-14
Keywords: DNA metabarcoding, fungi, grapevine trunk diseases, microbiome, plant pathogen
Citation: Leal CM, Geiger A, Molnár A, Váczy KZ, Kgobe G, Zsófi Z and Geml J (2024) Disentangling the effects of terroir, season, and vintage on the grapevine fungal pathobiome. Front. Microbiol. 14:1322559. doi: 10.3389/fmicb.2023.1322559
Edited by:
Fred O. Asiegbu, University of Helsinki, FinlandReviewed by:
Nagaraju Yalavarthi, Central Silk Board, IndiaPierluigi Reveglia, Spanish National Research Council (CSIC), Spain
Copyright © 2024 Leal, Geiger, Molnár, Váczy, Kgobe, Zsófi and Geml. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: József Geml, geml.jozsef@uni-eszterhazy.hu