Next Article in Journal
Pythium Damping-Off and Root Rot of Capsicum annuum L.: Impacts, Diagnosis, and Management
Previous Article in Journal
Bacillus velezensis CE 100 Inhibits Root Rot Diseases (Phytophthora spp.) and Promotes Growth of Japanese Cypress (Chamaecyparis obtusa Endlicher) Seedlings
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 4
10.3390/microorganisms9040822
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rootstocks Shape Their Microbiome—Bacterial Communities in the Rhizosphere of Different Grapevine Rootstocks

by
Leonie Dries
1,*,†,
Simone Bussotti
1,†,
Carlo Pozzi
2,
Robert Kunz
1,
Sylvia Schnell
3,
Otmar Löhnertz
1 and
Anne Vortkamp
4
1
Department of Soil Science and Plant Nutrition, Hochschule Geisenheim University, 65366 Geisenheim, Germany
2
Department of Agricultural and Environmental Sciences, University of Milano, 20133 Milano, Italy
3
Research Center for BioSystems, Land Use, and Nutrition (IFZ), Institute of Applied Microbiology, Justus-Liebig University Giessen, 35392 Giessen, Germany
4
REACH EUREGIO Start-up Center, University of Muenster, 48151 Muenster, Germany
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2021, 9(4), 822; https://doi.org/10.3390/microorganisms9040822
Submission received: 24 March 2021 / Revised: 12 April 2021 / Accepted: 12 April 2021 / Published: 13 April 2021
(This article belongs to the Section Plant Microbe Interactions)
Browse Figures
Versions Notes

Abstract

:
The microbiota associated with the rhizosphere is responsible for crucial processes. Understanding how the plant and its bacterial community interact is of great importance to face the upcoming agricultural and viticultural challenges. The composition of the bacterial communities associated with the rhizosphere of grapevines is the result of the interaction between many drivers: biogeography, edaphic factors, soil management and plant genotype. The experimental design of this study aimed to reduce the variability resulting from all factors except the genotype of the rootstock. This was made possible by investigating four ungrafted grapevine rootstock varieties of the same age, grown on the same soil under the same climatic conditions and managed identically. The bacterial communities associated with the rhizosphere of the rootstocks 1103 Paulsen, 140 Ruggeri, 161-49 Couderc and Kober 5BB were characterized with the amplicon based sequencing technique, targeting regions V4–V5 of 16S rRNA gene. Linear discriminant analysis effect Size (LEfSe) analysis was performed to determine differential abundant taxa. The four rootstocks showed similarities concerning the structure of the bacteria assemblage (richness and evenness). Nonetheless, differences were detected in the composition of the bacterial communities. Indeed, all investigated rootstocks recruited communities with distinguishable traits, thus confirming the role of rootstock genotype as driver of the bacteria composition.

1. Introduction

The continuous improvement of cultivation-independent techniques is driving a paradigm shift in the field of biology; the expression of a specific phenotype is no longer determined solely by the interaction between genotype and environment but has to take into account the host-associated microorganisms (microbiota) [1]. The microbiome (the genome of the microbiota) provides adaptability and metabolic diversity, broadening the plant capacity to overcome environmental changes and to cope with challenging conditions [2,3]. The microbes can be divided in epiphytes when they colonize the exterior surfaces or in endophytes when they colonize the interior surfaces [4]. Moreover, the distribution of the microbiota is not even throughout the plant, which can be interpreted as an assembly of several niches. The main distinction is between the above ground and the belowground part. However, agricultural productivity is based on microbial activity, much of which occurs in the soil [5]. The productivity of agricultural systems is highly dependent on the functional processes of microbial communities in the soil, more specifically in the rhizosphere [6]. In particular, the rhizosphere microbiome can positively influence plant development, overall vigor and plant growth by providing nutrients, and has generally been targeted to identify positive relationships between plants and microorganisms [7,8]. However, even if there are no physical barriers, the composition of the soil microbiota is not the same as the rhizosphere microbiota. The soil serves as the initial inoculum, but then the plant actively recruits the bacteria from the soil [9]. Thus, the rhizospheric communities are a subset of the more complex soil microbiota [10]. This selection occurs at the level of the root compartment, which consists of the rhizosphere, rhizoplane and root endosphere [11,12]. Due to this selection, the choice of the right rootstock for recruitment of the soil microbiota is crucial, which is also relevant for grapevine. Even though rootstocks were adopted for their resistance to Phylloxera, it is necessary to take into account several other factors, such as grafting compatibility, rooting and propagation ability, resistance to nematodes and pathogens, tolerance to lime, salinity, drought and nutrient uptake [13]. However, the impact of genetic variation on the composition of the microbiome of grapevine species, particularly in the root compartment, is poorly studied [14]. And it can be expected that different rootstocks differentially select their microbial communities from the surrounding soil [4]. Consequently, rootstock genotype affects the selection and recruitment of bacteria that colonize aboveground plant organs such as flowers and fruits and can alter fruit quality [4,15]. Thus, rootstocks could play a role in determining the plant-associated bacterial community by affecting the microbial terroir [16]. This study considered rootstocks derived from the breeding of the most commonly used Vitis riparia, Vitis berlandieri and Vitis rupestris [17]. The experimental design of this study aimed to reduce the variability resulting from environmental factors to unravel the effect of the genotype of the grapevine rootstock on the bacterial communities in the rhizosphere. This was made possible by using ungrafted grapevines of the same age, grown in the same soil under the same climatic conditions and managed identically.

2. Materials and Methods

2.1. Vineyard Location and Vines

The experiment site was located at Geisenheim University, Germany (49°59′N, 7°57′E; 96 m above sea level) in the wine-growing region of Rheingau. The climate of the region is temperate oceanic (Köppen-Geiger classification: Cfb) [18]. The mean annual temperature of the 1986–2010 period is 10.5 °C and the total annual precipitation averages 543.0 mm [19]. The grapevines were planted in a soil mixture (two parts loess soil with one part sand) in cylindrical tubes with a diameter of 15 cm and a height of 120 cm in the year 2015 (Figure 1). The tubes are perforated at the bottom side allowing the water to drain. The following rootstocks were selected for the investigations: 1103 Paulsen (V. berlandieri x V. rupestris, clone 1 Gm, “1103 P” when abbreviated), 140 Ruggeri (V. berlandieri x V. rupestris clone 4 Gm, “140 Ru” when abbreviated), 161-49 Couderc (V. berlandieri x V. riparia clone 3 Gm, “161-49 C” when abbreviated) and Kober 5 BB (V. berlandieri x V. riparia clone 13-3 Gm).

2.2. Rhizosphere and Soil Sampling

Rhizosphere and soil sampling was performed in July 2020. The tubes with the grapevines were extracted from the ground. Subsequently, each tube was opened using an angle grinder, which completely exposed the root system. This procedure ensured accurate selection and sampling of the rhizosphere soil, which was taken in the very close proximity to the roots. Three replicates per rootstock genotype were investigated. For each replicate, two rhizosphere soil samples (>2 g each) for microbiome analyses were collected and stored at −80 °C. Additionally, for each replicate one soil sample (>100 g) was collected for soil analyses and stored at 4 °C.

2.3. Soil Analyses

All soil analyses were performed at the Department of Soil Science and Plant Nutrition, Hochschule Geisenheim University. The soil samples were dried for three days at 30 °C and sieved to a size of 2 mm. For all soil samples (12 samples, n = 3 per rootstock) the pH level as well as the concentration of phosphorus (P), potassium (K) and magnesium (Mg) were measured as described previously by Hendgen et al. [20], according to the VDLUFA database of methods [21]. Statistical analyses were conducted with Microsoft® Excel®.

2.4. DNA Extraction, Sample Preparation and Ion Torrent Sequencing

DNA from 500 mg rhizosphere soil per sample was extracted according to Lueders et al. [22] and Mettel et al. [23] with the following modification: instead of the described sodium dodecyl sulfate (SDS) solution, TPM (50 mM TrisHCl, 1.7% Polyvinylpyrrolidone, 20 mM MgCl2) was used. Extracted DNA was solved in water and stored at −80 °C. The Ion Torrent Polymerase-Chain-Reactions (PCRs) were conducted as described in Kaplan et al. [24]. The first Ion Torrent PCR was performed with KAPA HiFi Polymerase kit (Roche Sequencing Store, Basel, Switzerland) according to the manufacturer’s protocol. During amplification of the partial sequence of the hypervariable regions (V4-V5), the primer 926R (5′-CCGYCAATTYMTTTRAGTTT-3′ [25]) and 520F (5′-AYTGGGYDTAAAGNG-3′ [26]) were used. Amplification parameters were 3 min at 95 °C followed by 35 cycles at 98 °C for 20 s, 55 °C for 30 s, 72 °C for 30 s and finally 72 °C for 5 min. The amplified PCR product was confirmed with agarose gel electrophoresis and used as the template of the second Ion Torrent PCR reaction and barcoding according to Kaplan et al. [24]. The PCR was performed with the following amplification parameters: 3 min at 95 °C followed by 10 cycles at 98 °C for 20 s, 55 °C for 30 s, 72 °C for 30 s and finally 72 °C for 7 min. The PCR products were loaded onto an agarose gel and subsequently purified using NucleoSpin® Gel and PCR Clean-up (MACHEREY-NAGEL GmbH & Co. KG, Düren, Germany). Further, this product was purified with DNA purification beads NucleoMagVR NGS clean-up kit (MACHEREY-NAGEL GmbH & Co. KG, Düren, Germany). The Ion Torrent Sequencing was done according to Kaplan et al. [24].

2.5. Bioinformatic Analysis of the Sequencing Data

Bioinformatic analysis was performed with QIIME 2 2020.8 [27]. The obtained raw sequences were demultiplexed using cutadapt [28] with no errors allowed in the barcode sequences. Quality control, sequence denoising, clustering to amplicon sequence variants (ASVs) (i.e., 100% operational taxonomic units) [29], dereplication and removal of chimera sequences were done with DADA2 [30] (via q2-dada2). Thereby the first 15 bp were deleted and all sequences were cut to a length of 330 bp. Taxonomy was assigned to the ASVs using the q2-feature classifier [31], which uses naïve Bayes machine-learning classifiers [32] to assign taxonomies based on sequence k-mer frequencies. As a reference database, a trained version of the Silva ribosomal database [33] version 137 was used. All ASVs belonging to chloroplasts and mitochondria were removed. All ASVs were aligned with mafft [34] (via q2-alignment) and used to construct a phylogeny with fasttree2 [35] (via q2-phylogeny). Alpha-diversity metrics (Faith’s Phylogenetic Diversity [36]—measures of microbiome richness, and Pielou’s Evenness) were calculated in R version 4.0.2 [37] using the packages phyloseq version 1.34.0 [38] and picante version 1.8.2 [39]. Beta-diversity metrics (weighted UniFrac [40], unweighted UniFrac [41] and Bray-Curtis dissimilarity [42]—measures of microbiome composition dissimilarity), and Principal Coordinate Analysis (PCoA) were estimated using q2-diversity after samples were rarefied to 2220 sequences per sample. The distance matrices were also analyzed by the PERMANOVA test [43], followed by the PERMDISP test [44] to ensure the homogeneity of dispersion among the four rootstocks. Both the PERMANOVA and PERMDISP tests were performed using 999 permutations. Linear discriminant analysis effect Size (LEfSe) was run to assess differences in bacteria distribution among the rootstocks [45]. The online Galaxy platform [46] was used. The threshold for the logarithmic linear discriminant analysis (LDA) score was set at 2.0 and the non-parametric factorial Kruskal Wallis sum-rank test (α = 0.05) was performed with an all-against-all strategy. Significant taxa were visualized with a bar graph and in a taxonomic cladogram.

3. Results

3.1. Soil Analyses

The main chemical characteristics of the sampled soils were measured and a Kruskal-Wallis H test was conducted. The pH levels (Table 1) ranged between 7.53 ± 0.09 and 7.60 ± 0.00 without a statistical difference (p = 0.89). The P, K and Mg contents (Table 1) were also statistically similar (respectively, p = 0.34, p = 0.63, p = 0.15).

3.2. High-Throughput Amplicon Sequencing

The Ion Torrent sequencing generated 910,225 16S rRNA gene V4-V5 amplicon sequences, 72% of which (642,114) passed the Ion Torrent Sequencing quality selection. After barcoding errors cleaning, the remaining 578,539 reads were demultiplexed to assign each read to the respective sample. Four samples were removed from the subsequent analysis due to the low number of reads (<2827 features; one sample per rootstock). The denoised and filtered reads were clustered into 2034 amplicon sequence variants (ASVs).

3.3. Diversity and Richness of Bacterial Communities

Pielou’s evenness and Faith’s phylogenetic diversity were used to calculate alpha-diversity (i.e., the diversity within the group). All measures indicated that no significant differences between the structure of the rhizosphere communities were observable (Figure 2). Pielou’s evenness reported values approaching 1, meaning that bacteria were evenly distributed among species. To calculate beta-diversity (i.e., the diversity between groups) the estimators Bray-Curtis, unweighted UniFrac and weighted UniFrac were used. The distance matrices were ordered on a PCoA (Figure 3). The permutation test with pseudo-F ratios revealed a significant effect of the factor “rootstock genotype” on beta-diversity for Bray-Curtis (p = 0.001), unweighted UniFrac (p = 0.001) and weighted UniFrac (p = 0.002) indices. PERMDISP was run to assess the homogeneity of the dispersion. The test led to a p-value of 0.131 for Bray-Curtis, 0.572 for unweighted UniFrac and 0.258 for weighted UniFrac.

3.4. Taxonomy and Distribution among Rootstocks

The taxonomy was assigned using a trained classifier against the Silva database. The 2034 different ASVs, all belonging to the bacteria domain, clustered into 36 different phyla, 97 classes, 199 orders, 257 families and 365 genera (Supplementary Materials Table S1). The rhizosphere was dominated by Acidobacteriota (35%), Proteobacteria (22%), Latescibacteriota (15%), Methylomirabilota (6%) and Gemmatimonadota (4%). The most represented classes were Vicinamibacteria (17%), Gammaproteobacteria (15%), Latescibacterota (15%), Acidobacteriae (9%), Alphaproteobacteria (7%), Methylomirabilia (6%) and Gemmatimonadetes (4%). The linear discriminant analysis effect size (LEfSe) detected 68 bacterial clades differently distributed among the rootstocks (Figure 4 and Figure 5). Significant differences (LDA score > 2.0) were found for each taxonomic level: phylum (2), class (6), order (15), family (20) and genus (25). Proteobacteria was the most heterogeneously distributed phylum (21 clades, 12 belonging to Alphaproteobacteria and 9 belonging to Gammaproteobacteria), followed by the phyla Chloroflexi (11 clades), Methylomirabilota (7 clades) and Myxococcota (7 clades). The rootstock 1103 Paulsen reported the highest number of differentially abundant taxa (32).

4. Discussion

The experimental design of this study aimed at minimizing the variability from all factors but the rootstock genotype. This was made possible by using grapevines of the same age, located in the same soil sharing also the same climatic conditions, and managed identically. The amplicon sequencing analysis revealed that rootstock rhizospheres recruited complex bacterial communities mainly composed of Acidobacteriota, Proteobacteria, Latescibacterota, Methylomirabilota and Gemmatimonadota. These results are in partial agreement with the data of other recent studies [4,14,16]. According to Berlanas et al. [14], the dominant phyla found in two vineyards located in Northeastern Spain were Proteobacteria, Actinobacteria, Acidobacteria and Bacteroidetes. Similarly, Marasco et al. [16] reported phylum Proteobacteria as prevalent, followed by Actinobacteria, Bacteroidetes, Chloroflexi and Acidobacteria in an Oltrepò Pavese (Italy) vineyard. These results highlight that there is a core of bacteria always recruited by grapevine despite the different conditions (e.g., biogeography, and climatic or edaphic factors). Among them, Proteobacteria is commonly the dominant phylum followed by Actinobacteria. Both phyla play an important role in carbon cycling and production of secondary metabolites [47]. Proteobacteria comprise organisms with a broad variety of metabolic capabilities; alpha-, beta-, gamma- and delta- Proteobacteria are the classes mostly described in soil studies. In this study, only Gammaproteobacteria and Alphaproteobacteria were found. Both classes are considered copiotrophs (r-strategist) and can easily colonize environments characterized by abundance of nutrients, such as the rhizosphere [48]. Instead, Actinobacteria are described as oligotrophic k-strategists growing slowly and with low nutritional requirements [49]. In this project, the percentage of ASVs assigned to this phylum was rather low (0.56%) in comparison to other studies. Conversely, the percentage shown by Acidobacteriota phylum was unusually high (~ 35%), and greater than Proteobacteria phylum (~22 %). Acidobacteriota phylum is one of the most widespread and abundant on the planet [50], but its abundance is reported to decrease moving from bulk soil to rhizosphere [2]. Information about Acidobacteriota functionality is scarce since there are issues in cultivating the majority of the members of this phylum [51]. A negative correlation between the abundance of Acidobacteriota and concentration of the organic carbon in soil was reported [50]. This, together with the low growth rate, may indicate this phylum as oligotrophic [50]. Additionally, Zarraonaindia et al. [4] reported high abundance of Acidobacteriota together with unusually low abundance of Actinobacteria. It may be suggested that Acidobacteriota can replace Actinobacteria when some biotic or abiotic factors hamper the normal thriving of the latter. Indeed, both phyla are described as oligotrophic k-strategist [50] and occupy the same niche. Remarkably, this study is the first to report the occurrence of the phylum Latescibacterota (previously known as WS3) among the most abundant phyla in vineyard soils. Gemmatimonadota phylum was first associated with vineyard soils by Novello et al. [52], and this study confirmed their presence in the grapevine rhizosphere. Gemmatimonadota is reported as one of the top nine phyla found in soils [53]. Nonetheless, little is known about its metabolic capabilities since most of its members are uncultivated. It was observed that the microbial communities associated with the four rootstocks did not significantly differ in evenness and richness. Indeed, all the four rootstocks showed a bacterial community with no clearly dominating species and with similar phylogenetic richness. This is consistent with the theory of a well conserved core of bacteria associated with the grapevine rhizosphere [14]. This study, however, showed that rootstocks 1103P, 140 Ru, 161-49C and Kober 5BB were able to assemble distinct bacterial microbiota at the rhizosphere level. The percentage explained by weighted (46%) and unweighted UniFrac (22%) suggests that the differences are mainly driven by a different recruitment of the same bacteria. Moreover, the LEfSe analysis detected 68 clades significantly differentially recruited by the four investigated rootstocks. These differences confirm the rootstock genotype as a substantial driver of the bacterial communities associated with the rhizosphere. The rootstock genotype may exert its influence on the rhizosphere associated communities through specific root exudates [54,55,56] and by means of its immune system. Indeed, it has been proven that the immune system can affect not only the microbial composition of the root interior surfaces but also the communities thriving near the root [57,58]. Marasco et al. [16] have investigated whether different rootstocks affect the recruitment of the bacteria from the surrounding soil comparing ungrafted and grafted grapevines of the Vitis vinifera cultivar Barbera, all cultivated in the same soil. They concluded that the rootstock can influence the diversity and richness of the bacterial communities associated with the root and that the rootstock is a factor determining the specificity of the microbiota independently of the scion cultivar. D’Amico et al. [59] have highlighted the importance of the microbiota recruited by the rootstock with a cultivation-independent technique: the same scion (V. vinifera cv Lambrusco) was grafted onto two different rootstocks (1103 Paulsen and Kober 5BB) in the same vineyard, where only the grapevines grafted on 1103 Paulsen had potassium absorption problems. The analysis of the microbiota with an amplicon-based approach has revealed that the rootstock 1103 Paulsen is not able to successfully recruit several potassium solubilizing microorganisms [59]. Thus, the knowledge of the microbiota associated with diverse rootstocks is valuable to make the right choice at the moment of the implant, to intervene when the choice is already taken, and also for breeders to select for rootstocks with a specific microbiota.

5. Conclusions

These results show that the rootstock genotype affects the composition but not the structure of the bacterial community in the rhizosphere. However, recognizing the associations between bacteria and rootstock is crucial but still not sufficient to exploit this information from a practical perspective. In conclusion, it is essential to understand the recruitment mechanisms of grapevine to microbial communities and to identify the beneficial taxa in a more detailed context. Following studies should use a whole genome sequencing approach and combine it with other meta -omics methods. For example, the use of metatranscriptomics and metaproteomics together allows the characterization of the microbial communities thriving in the rhizosphere and simultaneously the comprehension of their role. However, applying several -omics methods will inevitably lead to a better comprehension of the rhizosphere microbiome in relation to soil type, grapevine variety, rootstock and environmental conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9040822/s1, Table S1: Amplicon Sequenced Variants Table.

Author Contributions

Conceptualization, A.V., L.D. and S.B.; methodology, L.D., S.B. and R.K.; formal analysis, S.B.; investigation, L.D. and S.B.; resources S.S. and O.L.; data curation, L.D. and S.B.; writing—original draft preparation, L.D. and S.B.; writing—review and editing, L.D., S.B., C.P., S.S., R.K., O.L. and A.V.; visualization, S.B.; supervision, C.P., S.S., O.L. and A.V.; project administration, L.D., A.V. and O.L.; funding acquisition, A.V. and O.L. All authors have read and agreed to the published version of the manuscript.

Funding

We thank Forschungsring des Deutschen Weinbaus (FDW) for funding this project.

Data Availability Statement

The sequencing reads from this study are openly available in the National Center for Biotechnology Information (ncbi) under accession no. PRJNA716534. This data can be found at https://www.ncbi.nlm.nih.gov/bioproject/716534.

Acknowledgments

We appreciate and acknowledge the members of the Department of Soil Science and Plant Nutrition, Hochschule Geisenheim University, for their assistance in sample taking and soil analyses. We thank Joachim Schmid, Department of Grapevine Breeding, Hochschule Geisenheim University, for the opportunity to sample the experimental vineyard. We thank Annette Reineke and Olivia Herczynski, Department of Crop Protection, Hochschule Geisenheim University, for the opportunity to use the molecular biology laboratory and for their helpful advices. We are very grateful to Bellinda Schneider, Justus-Liebig University Gießen, for technical support with Ion Torrent Sequencing. We further thank Maximilian Hendgen for carefully proofreading the manuscript. We acknowledge support by the German Research Foundation (Deutsche Forschungsgemeinschaft DFG)—project number 432888308—and the Open Access Publishing Fund of Geisenheim University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bordenstein, S.R.; Theis, K.R. Host Biology in Light of the Microbiome: Ten Principles of Holobionts and Hologenomes. PLoS Biol. 2015, 13, e1002226. [Google Scholar] [CrossRef] [Green Version]
  2. Rosenberg, E.; Zilber-Rosenberg, I. Microbes Drive Evolution of Animals and Plants: The Hologenome Concept. mBio 2016, 7, e01395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bulgarelli, D.; Schlaeppi, K.; Spaepen, S.; Ver Loren van Themaat, E.; Schulze-Lefert, P. Structure and functions of the bacterial microbiota of plants. Annu. Rev. Plant Biol. 2013, 64, 807–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zarraonaindia, I.; Owens, S.M.; Weisenhorn, P.; West, K.; Hampton-Marcell, J.; Lax, S.; Bokulich, N.A.; Mills, D.A.; Martin, G.; Taghavi, S.; et al. The soil microbiome influences grapevine-associated microbiota. mBio 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bakker, M.G.; Manter, D.K.; Sheflin, A.M.; Weir, T.L.; Vivanco, J.M. Harnessing the rhizosphere microbiome through plant breeding and agricultural management. Plant Soil 2012, 360, 1–13. [Google Scholar] [CrossRef]
  6. Vega-Avila, A.D.; Gumiere, T.; Andrade, P.A.M.; Lima-Perim, J.E.; Durrer, A.; Baigori, M.; Vazquez, F.; Andreote, F.D. Bacterial communities in the rhizosphere of Vitis vinifera L. cultivated under distinct agricultural practices in Argentina. Antonie Van Leeuwenhoek 2015, 107, 575–588. [Google Scholar] [CrossRef] [PubMed]
  7. Ryan, P.R.; Dessaux, Y.; Thomashow, L.S.; Weller, D.M. Rhizosphere engineering and management for sustainable agriculture. Plant Soil 2009, 321, 363–383. [Google Scholar] [CrossRef]
  8. Taye, Z.M.; Helgason, B.L.; Bell, J.K.; Norris, C.E.; Vail, S.; Robinson, S.J.; Parkin, I.A.P.; Arcand, M.; Mamet, S.; Links, M.G.; et al. Core and Differentially Abundant Bacterial Taxa in the Rhizosphere of Field Grown Brassica napus Genotypes: Implications for Canola Breeding. Front. Microbiol. 2019, 10, 3007. [Google Scholar] [CrossRef] [Green Version]
  9. Gilbert, J.A.; van der Lelie, D.; Zarraonaindia, I. Microbial terroir for wine grapes. Proc. Natl. Acad. Sci. USA 2014, 111, 5–6. [Google Scholar] [CrossRef] [Green Version]
  10. Berendsen, R.L.; Pieterse, C.M.J.; Bakker, P.A.H.M. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012, 17, 478–486. [Google Scholar] [CrossRef]
  11. Edwards, J.; Johnson, C.; Santos-Medellín, C.; Lurie, E.; Podishetty, N.K.; Bhatnagar, S.; Eisen, J.A.; Sundaresan, V. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. USA 2015, 112, E911–E920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Bonito, G.; Reynolds, H.; Robeson, M.S.; Nelson, J.; Hodkinson, B.P.; Tuskan, G.; Schadt, C.W.; Vilgalys, R. Plant host and soil origin influence fungal and bacterial assemblages in the roots of woody plants. Mol. Ecol. 2014, 23, 3356–3370. [Google Scholar] [CrossRef]
  13. Granett, J.; Walker, A.M.; Kocsis, L.; Omer, A.D. Biology and management of grape phylloxera. Annu. Rev. Entomol. 2001, 387–412. [Google Scholar] [CrossRef] [PubMed]
  14. Berlanas, C.; Berbegal, M.; Elena, G.; Laidani, M.; Cibriain, J.F.; Sagües, A.; Gramaje, D. The Fungal and Bacterial Rhizosphere Microbiome Associated With Grapevine Rootstock Genotypes in Mature and Young Vineyards. Front. Microbiol. 2019, 10, 1142. [Google Scholar] [CrossRef] [PubMed]
  15. Bokulich, N.A.; Thorngate, J.H.; Richardson, P.M.; Mills, D.A. Microbial biogeography of wine grapes is conditioned by cultivar, vintage, and climate. Proc. Natl. Acad. Sci. USA 2014, 111, E139–E148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Marasco, R.; Rolli, E.; Fusi, M.; Michoud, G.; Daffonchio, D. Grapevine rootstocks shape underground bacterial microbiome and networking but not potential functionality. Microbiome 2018, 6, 3. [Google Scholar] [CrossRef] [PubMed]
  17. Dry, P.R.; Coombe, B.G.; Anderson, C.J. (Eds.) Viticulture, 2nd ed.; Winetitles Pty.: Adelaide, Australia, 2005; ISBN 0975685007. [Google Scholar]
  18. Reineke, A.; Selim, M. Elevated atmospheric CO2 concentrations alter grapevine (Vitis vinifera) systemic transcriptional response to European grapevine moth (Lobesia botrana) herbivory. Sci. Rep. 2019, 9, 2995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. DWD Climate Data Center. Historical Daily Station Observations (Temperature, Pressure, Precipitation, Sunshine Duration, etc.) for Germany; Version v006; DWD: Köln, Germany, 2018. [Google Scholar]
  20. Hendgen, M.; Hoppe, B.; Döring, J.; Friedel, M.; Kauer, R.; Frisch, M.; Dahl, A.; Kellner, H. Effects of different management regimes on microbial biodiversity in vineyard soils. Sci. Rep. 2018, 8, 9393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Bassler, R. (Ed.) Band I: Die Untersuchung von Böden. In Handbuch der Landwirtschaftlichen Versuchs- und Untersuchungsmethodik; VDLUFA-Verlag: Darmstadt, Germany, 2011. [Google Scholar]
  22. Lueders, T.; Manefield, M.; Friedrich, M.W. Enhanced sensitivity of DNA- and rRNA-based stable isotope probing by fractionation and quantitative analysis of isopycnic centrifugation gradients. Environ. Microbiol. 2004, 6, 73–78. [Google Scholar] [CrossRef]
  23. Mettel, C.; Kim, Y.; Shrestha, P.M.; Liesack, W. Extraction of mRNA from soil. Appl. Environ. Microbiol. 2010, 76, 5995–6000. [Google Scholar] [CrossRef] [Green Version]
  24. Kaplan, H.; Ratering, S.; Felix-Henningsen, P.; Schnell, S. Stability of in situ immobilization of trace metals with different amendments revealed by microbial 13C-labelled wheat root decomposition and efflux-mediated metal resistance of soil bacteria. Sci. Total Environ. 2019, 659, 1082–1089. [Google Scholar] [CrossRef]
  25. Quince, C.; Lanzen, A.; Davenport, R.J.; Turnbaugh, P.J. Removing noise from pyrosequenced amplicons. BMC Bioinf. 2011, 12, 38. [Google Scholar] [CrossRef]
  26. Claesson, M.J.; O’Sullivan, O.; Wang, Q.; Nikkilä, J.; Marchesi, J.R.; Smidt, H.; de Vos, W.M.; Ross, R.P.; O’Toole, P.W. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS ONE 2009, 4, e6669. [Google Scholar] [CrossRef] [Green Version]
  27. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  28. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011, 17, 10. [Google Scholar] [CrossRef]
  29. Callahan, B.J.; McMurdie, P.J.; Holmes, S.P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017, 2639–2643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  31. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2′s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef] [PubMed]
  32. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  33. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  34. Katoh, K.; Misawa, K.; Kuma, K.-i.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 3059–3066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree 2—Approximately Maximum-Likelihood Trees for Large Alignments. PLoS ONE 2010, 5, e9490. [Google Scholar] [CrossRef] [PubMed]
  36. Faith, D.P. Conservation evaluation and phylogenetic diversity. Biol. Conserv. 1992, 61, 1–10. [Google Scholar] [CrossRef]
  37. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  38. McMurdie, P.J.; Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  39. Kembel, S.W.; Cowan, P.D.; Helmus, M.R.; Cornwell, W.K.; Morlon, H.; Ackerly, D.D.; Blomberg, S.P.; Webb, C.O. Picante: R tools for integrating phylogenies and ecology. Bioinformatics 2010, 26, 1463–1464. [Google Scholar] [CrossRef] [Green Version]
  40. Lozupone, C.A.; Hamady, M.; Kelley, S.T.; Knight, R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 2007, 73, 1576–1585. [Google Scholar] [CrossRef] [Green Version]
  41. Lozupone, C.; Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef] [Green Version]
  42. Bray, J.R.; Curtis, J.T. An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecol. Monogr. 1957, 27, 325–349. [Google Scholar] [CrossRef]
  43. Anderson, M.J. Permutational Multivariate Analysis of Variance (PERMANOVA); Wiley & Sons: Hoboken, NJ, USA, 2017; pp. 1–15. [Google Scholar] [CrossRef]
  44. Anderson, M.J. Distance-based tests for homogeneity of multivariate dispersions. Biometrics 2006, 62, 245–253. [Google Scholar] [CrossRef]
  45. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [Green Version]
  46. Afgan, E.; Baker, D.; Batut, B.; van den Beek, M.; Bouvier, D.; Cech, M.; Chilton, J.; Clements, D.; Coraor, N.; Grüning, B.A.; et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018, 46, W537–W544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Jenkins, S.N.; Waite, I.S.; Blackburn, A.; Husband, R.; Rushton, S.P.; Manning, D.C.; O’Donnell, A.G. Actinobacterial community dynamics in long term managed grasslands. Antonie Van Leeuwenhoek 2009, 95, 319–334. [Google Scholar] [CrossRef]
  48. Fierer, N.; Bradford, M.A.; Jackson, R.B. Toward an ecological classification of soil bacteria. Ecology 2007, 1354–1364. [Google Scholar] [CrossRef] [PubMed]
  49. Fierer, N.; Schimel, J.P.; Holden, P.A. Variations in microbial community composition through two soil depth profiles. Soil Biol. Biochem. 2003, 35, 167–176. [Google Scholar] [CrossRef]
  50. Kielak, A.M.; Barreto, C.C.; Kowalchuk, G.A.; van Veen, J.A.; Kuramae, E.E. The Ecology of Acidobacteria: Moving beyond Genes and Genomes. Front. Microbiol. 2016, 7, 744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Navarrete, A.A.; Kuramae, E.E.; de Hollander, M.; Pijl, A.S.; van Veen, J.A.; Tsai, S.M. Acidobacterial community responses to agricultural management of soybean in Amazon forest soils. FEMS Microbiol. Ecol. 2013, 83, 607–621. [Google Scholar] [CrossRef] [Green Version]
  52. Novello, G.; Gamalero, E.; Bona, E.; Boatti, L.; Mignone, F.; Massa, N.; Cesaro, P.; Lingua, G.; Berta, G. The Rhizosphere Bacterial Microbiota of Vitis vinifera cv. Pinot Noir in an Integrated Pest Management Vineyard. Front. Microbiol. 2017, 8, 1528. [Google Scholar] [CrossRef]
  53. Janssen, P.H. Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 2006, 72, 1719–1728. [Google Scholar] [CrossRef] [Green Version]
  54. Matthews, A.; Pierce, S.; Hipperson, H.; Raymond, B. Rhizobacterial Community Assembly Patterns Vary Between Crop Species. Front. Microbiol. 2019, 10, 581. [Google Scholar] [CrossRef] [Green Version]
  55. Badri, D.V.; Chaparro, J.M.; Zhang, R.; Shen, Q.; Vivanco, J.M. Application of natural blends of phytochemicals derived from the root exudates of Arabidopsis to the soil reveal that phenolic-related compounds predominantly modulate the soil microbiome. J. Biol. Chem. 2013, 288, 4502–4512. [Google Scholar] [CrossRef] [Green Version]
  56. Bulgarelli, D.; Rott, M.; Schlaeppi, K.; Ver Loren van Themaat, E.; Ahmadinejad, N.; Assenza, F.; Rauf, P.; Huettel, B.; Reinhardt, R.; Schmelzer, E.; et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 2012, 488, 91–95. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, Z.; Beskrovnaya, P.; Melnyk, R.A.; Hossain, S.S.; Khorasani, S.; O’Sullivan, L.R.; Wiesmann, C.L.; Bush, J.; Richard, J.D.; Haney, C.H. A Genome-Wide Screen Identifies Genes in Rhizosphere-Associated Pseudomonas Required to Evade Plant Defenses. mBio 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  58. Lebeis, S.L.; Paredes, S.H.; Lundberg, D.S.; Breakfield, N.; Gehring, J.; McDonald, M.; Malfatti, S.; Glavina del Rio, T.; Jones, C.D.; Tringe, S.G.; et al. PLANT MICROBIOME. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science 2015, 349, 860–864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. D’Amico, F.; Candela, M.; Turroni, S.; Biagi, E.; Brigidi, P.; Bega, A.; Vancini, D.; Rampelli, S. The Rootstock Regulates Microbiome Diversity in Root and Rhizosphere Compartments of Vitis vinifera Cultivar Lambrusco. Front. Microbiol. 2018, 9, 2240. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 09 00822 g001 550
Figure 1. Experimental vineyard. The grapevine rootstocks are planted in cylindrical tubes with a diameter of 15 cm and a height of 120 cm.
Figure 1. Experimental vineyard. The grapevine rootstocks are planted in cylindrical tubes with a diameter of 15 cm and a height of 120 cm.
Microorganisms 09 00822 g001
Microorganisms 09 00822 g002 550
Figure 2. Boxplots reporting the results of alpha diversity analysis. (a) Faith’s phylogenetic plot showing similar values for the four rootstocks. (b) Pielou’s evenness plot showing values higher than 0.895 for all the rootstocks.
Figure 2. Boxplots reporting the results of alpha diversity analysis. (a) Faith’s phylogenetic plot showing similar values for the four rootstocks. (b) Pielou’s evenness plot showing values higher than 0.895 for all the rootstocks.
Microorganisms 09 00822 g002
Microorganisms 09 00822 g003 550
Figure 3. Principal Coordinate Analyses (PCoAs) reporting the distance matrices of the beta diversity indexes: (a) Bray-Curtis, (b) unweighted UniFrac and (c) weighted UniFrac.
Figure 3. Principal Coordinate Analyses (PCoAs) reporting the distance matrices of the beta diversity indexes: (a) Bray-Curtis, (b) unweighted UniFrac and (c) weighted UniFrac.
Microorganisms 09 00822 g003
Microorganisms 09 00822 g004 550
Figure 4. Histogram of the linear discriminant analysis (LDA) scores reveals the most differentially abundant taxa among the four rootstocks 1103 Paulsen, 140 Ruggeri, 161-49 Couderc and Kober 5BB.
Figure 4. Histogram of the linear discriminant analysis (LDA) scores reveals the most differentially abundant taxa among the four rootstocks 1103 Paulsen, 140 Ruggeri, 161-49 Couderc and Kober 5BB.
Microorganisms 09 00822 g004
Microorganisms 09 00822 g005 550
Figure 5. Linear discriminant analysis (LDA) effect size taxonomic cladogram comparing the bacterial diversity of the four rootstocks. Each ring represents a taxonomic level, with kingdom, phylum, class, order, family and genus appearing from the center to the periphery. Each node represents a taxonomic unit found in the dataset. For each significantly discriminating taxon detected, the corresponding node and branch region in the taxonomic cladogram is colored according to the highest ranked group for that taxon. If the taxon is not significantly different between rootstocks, the corresponding node is colored yellow.
Figure 5. Linear discriminant analysis (LDA) effect size taxonomic cladogram comparing the bacterial diversity of the four rootstocks. Each ring represents a taxonomic level, with kingdom, phylum, class, order, family and genus appearing from the center to the periphery. Each node represents a taxonomic unit found in the dataset. For each significantly discriminating taxon detected, the corresponding node and branch region in the taxonomic cladogram is colored according to the highest ranked group for that taxon. If the taxon is not significantly different between rootstocks, the corresponding node is colored yellow.
Microorganisms 09 00822 g005
Table
Table 1. Main chemical characteristics of the soil of the four different rootstocks. Phosphorus (P), potassium (K) and magnesium (Mg) values are expressed in mg kg-1.
Table 1. Main chemical characteristics of the soil of the four different rootstocks. Phosphorus (P), potassium (K) and magnesium (Mg) values are expressed in mg kg-1.
Soil Characteristics1103 P140 Ru161–49CKober 5BB
pH7.53 ± 0.127.53 ± 0.097.60 ± 0.007.57 ± 0.12
P117.72 ± 9.42194.75 ± 76.05152.60 ± 19.82140.97 ± 52.48
K165.98 ± 20.33237.90 ± 137.09193.64 ± 19.56157.68 ± 58.68
Mg116.61 ± 12.3986.45 ± 52.66152.79 ± 15.04138.72 ± 22.57
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dries, L.; Bussotti, S.; Pozzi, C.; Kunz, R.; Schnell, S.; Löhnertz, O.; Vortkamp, A. Rootstocks Shape Their Microbiome—Bacterial Communities in the Rhizosphere of Different Grapevine Rootstocks. Microorganisms 2021, 9, 822. https://doi.org/10.3390/microorganisms9040822

AMA Style

Dries L, Bussotti S, Pozzi C, Kunz R, Schnell S, Löhnertz O, Vortkamp A. Rootstocks Shape Their Microbiome—Bacterial Communities in the Rhizosphere of Different Grapevine Rootstocks. Microorganisms. 2021; 9(4):822. https://doi.org/10.3390/microorganisms9040822

Chicago/Turabian Style

Dries, Leonie, Simone Bussotti, Carlo Pozzi, Robert Kunz, Sylvia Schnell, Otmar Löhnertz, and Anne Vortkamp. 2021. "Rootstocks Shape Their Microbiome—Bacterial Communities in the Rhizosphere of Different Grapevine Rootstocks" Microorganisms 9, no. 4: 822. https://doi.org/10.3390/microorganisms9040822

APA Style

Dries, L., Bussotti, S., Pozzi, C., Kunz, R., Schnell, S., Löhnertz, O., & Vortkamp, A. (2021). Rootstocks Shape Their Microbiome—Bacterial Communities in the Rhizosphere of Different Grapevine Rootstocks. Microorganisms, 9(4), 822. https://doi.org/10.3390/microorganisms9040822

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop