Next Article in Journal
Actinomyces spp. Prosthetic Vascular Graft Infection (PVGI): A Multicenter Case-Series and Narrative Review of the Literature
Next Article in Special Issue
The Win–Win Effects of an Invasive Plant Biochar on a Soil–Crop System: Controlling a Bacterial Soilborne Disease and Stabilizing the Soil Microbial Community Network
Previous Article in Journal
Biological Control of the Plant Pathogens
Previous Article in Special Issue
The Combined Applications of Microbial Inoculants and Organic Fertilizer Improve Plant Growth under Unfavorable Soil Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Generation Ecosystem Selection of Rhizosphere Microbial Communities Associated with Plant Genotype and Biomass in Arabidopsis thaliana

1
Department of Biology, University of Massachusetts, Boston, MA 02125, USA
2
Institute of Plant Biology, Biological Research Centre, 6726 Szeged, Hungary
*
Author to whom correspondence should be addressed.
Microorganisms 2023, 11(12), 2932; https://doi.org/10.3390/microorganisms11122932
Submission received: 2 November 2023 / Revised: 1 December 2023 / Accepted: 1 December 2023 / Published: 6 December 2023

Abstract

:
The role of the microbiome in shaping the host’s phenotype has emerged as a critical area of investigation, with implications in ecology, evolution, and host health. The complex and dynamic interactions involving plants and their diverse rhizospheres’ microbial communities are influenced by a multitude of factors, including but not limited to soil type, environment, and plant genotype. Understanding the impact of these factors on microbial community assembly is key to yielding host-specific and robust benefits for plants, yet it remains challenging. Here, we conducted an artificial ecosystem selection experiment for eight generations of Arabidopsis thaliana Ler and Cvi to select soil microbiomes associated with a higher or lower biomass of the host. This resulted in divergent microbial communities shaped by a complex interplay between random environmental variations, plant genotypes, and biomass selection pressures. In the initial phases of the experiment, the genotype and the biomass selection treatment had modest but significant impacts. Over time, the plant genotype and biomass treatments gained more influence, explaining ~40% of the variation in the microbial community’s composition. Furthermore, a genotype-specific association of plant-growth-promoting rhizobacterial taxa, Labraceae with Ler and Rhizobiaceae with Cvi, was observed under selection for high biomass.

1. Introduction

The conventional understanding of the host phenotype involves genetics and the environment shaping observable traits. Yet, the last two decades have underscored the microbiome’s significance in shaping the host phenotype, driven by extensive research on its role in ecology, evolution, and host health. The plant microbiome represents a rich source of functional diversity that is not encoded within the host genome. The interactions between the plant host and its microbiome are dynamic and reciprocal, with the plant shaping its immediate environment by exuding specific metabolites, thereby promoting the growth of specific microbial taxa, while the microbiome in turn influences plant health and growth [1]. A multitude of biotic and abiotic factors influence the dynamic nature of the microbiome, including but not limited to plant root exudates [2,3,4]; soil type [5,6,7,8,9]; environment [10]; and various aspects of the plant, including species [11], genotype [4,12,13,14,15] and developmental stage [16,17].
Previous studies in Arabidopsis thaliana have shown that the microbial diversity of the soil decreased with proximity to the endophytic compartment, as well as that microbes found in the bulk soil were not enriched in the endophytic compartment, and vice versa [18,19,20]. Furthermore, the microbiomes of the greenhouse-grown plants were similar to those of the field-grown plants, even across different plant compartments—rhizosphere, woody stem, and endophytic compartment—indicating an active role of the plant host in creating and maintaining an environment where certain microbes have improved fitness [21]. The soil environment was a stronger predictor of the rhizosphere’s microbiome structure, though plant genotype has been shown to have a weak influence in some studies [18,20,22]. A follow-up study confirmed that the variation in microbiome communities depends more on the environment than on different Arabidopsis varieties and sister species. Though there were some differences between different host species, the core microbiome’s interaction with the plant host at the root endophytic zone was largely consistent and reproducible [23].
Unfortunately, most prior studies were conducted in a single plant growth cycle, precluding an analysis of the temporal stability of the plant–microbiome association. Observations from plant–soil feedback studies highlight that plant–microbiome interactions can be altered by successive growth cycles of a given plant in the same soil [24,25]. Prior studies have employed plant-mediated selection on the soil ecosystem over multiple generations, and this has resulted in consistent effects on plant characteristics such as biomass, flowering time, and germination [26,27]. By selecting for soils where hosts show the desired phenotype, such as high biomass or altered flowering time, it becomes possible to enrich microbes that modulate host traits. Additionally, studies have shown that plants’ responses to abiotic stress, such as drought and salt stress, could be influenced in two ways. The first of these is through generational selection, where plants with desirable stress responses are chosen over multiple generations [28,29,30]. The second occurs by introducing beneficial microbes that are associated with plants and have previously lived in similar drought or salt conditions [31]. Mueller et al. [28] used generational selection to create microbial communities that could improve plant seed production by up to 205%. A soil inoculum obtained from drought-exposed soils improved wheat biomass under drought conditions [31]. Older studies, such as that of Swenson et al. [27], have also shown similar results under optimal growing conditions without necessarily exploring the microbiome. A recent study conducted on the rhizosphere microbiome of wild and domesticated tomato plants over multiple generations demonstrated an escalating influence of the host’s genotype on the microbiome community [32]. A separate study on the phyllosphere microbiome of various tomato genotypes indicated a declining impact of host genotype across successive generations [33].
Despite several earlier efforts, the development of a robust and host-specific microbial community with long-lasting beneficial effects on the plant host remains a challenge. Using an artificial ecosystem selection experiment, we sought to advance our understanding of how the plant genotype, environment, and biomass selection treatment (henceforth referred to as biomass treatment) impact the assembly of host–microbial communities. In this study, we hypothesized that: (i) Genotype, biomass treatment, and environment would play a critical role in microbial community assembly; (ii) generational selection of the soil ecosystem would lead to phenotypes with high and low biomasses; (iii) the influence of plant genotype over multiple generations would result in genotype and biomass-treatment-specific microbiomes. This study enhances our understanding of the temporal dynamics involved in microbial community assembly in the rhizosphere and sets the stage for targeted methodologies establishing resilient plant–microbiome interactions, with potential applications in sustainable agriculture and ecosystem health.

2. Materials and Methods

2.1. Multi-Generation Selection of Soil Ecosystem

Arabidopsis thaliana Cvi and Ler accessions were cultivated using custom-made “rhizotubes” (Stuart Morey, Univ. of Massachusetts, Boston, unpublished). The tubes were equipped with a black polyethylene sleeve, which effectively blocked light penetration. This design also facilitated the convenient removal of the plants from the pot, granting full access to the root system (Supplementary Figure S1).
The potting soil used was PRO-MIX PGX, a commercial mixture comprising 80–90% sphagnum peat moss and small quantities of perlite. It was autoclaved twice for 40 min, with a 48 h interval between each sterilization. The sterile potting soil was then sifted through a 3 mm sieve and combined with field soil at a 6:1 ratio for the first generation. The field soil, obtained from the Center for Agricultural Research in Waltham, Massachusetts, was composed of 44% sand, 49% silt, and 7% clay, which is representative of agricultural and grassland ecosystems.
The mixture was homogenized using a custom-made cement mixer and attached to sterile bins to avoid cross-contamination. Distilled water was gradually added to achieve a final ratio of 1:2 (water to soil). The soil was incubated at room temperature for two days following inoculation before the seeds were planted. All seeds used in the experiment were obtained from a single parent plant. Before planting, the seeds were surface-sterilized by treatment with a solution of 50% bleach v/v and a drop of Tween 20 for 10 min, then rinsed ten times with sterile distilled water. We placed 3–4 seeds in the center of each pot and kept only one seedling per pot after emergence. All plants were grown at 22 °C day/18 °C night with a 12/12 h day–night cycle in a controlled growth chamber. The relative humidity ranged from 35–60%, and the light intensity was 96 µE. Fertilizer was not applied to the plants.
In the first generation of the experiment, 100 plants of each accession (Cvi and Ler) were grown separately in individual pots. The plants were arranged in a randomized block design to avoid batch effects. All plants were harvested 35 days after germination. The above-ground portion of the plant, the rhizosphere, and the bulk soil from each rhizotube were separated, ensuring that there was no cross-contamination (Supplementary Figure S1). The above-ground parts of the plant, the stem, and the leaves were dried at 70 °C for 4 days. All of the plants were weighed individually on a closed weighing scale accurate to 1 mg. The root–soil complex (comprising the rhizosphere and the endosphere, henceforth referred to as the rhizosphere) was obtained by shaking the excess soil off the root and placing the root and the remaining attached soil into a sterile 5 mL tube. Bulk soil samples were transferred to a sterile Ziplock bag. Tubes and bags were immediately transferred to dry ice and then stored at −80 °C for DNA analysis.
All subsequent generations were inoculated with a mixture of previously collected bulk soil from the top five and bottom five plants based on above-ground biomass, then combined with sterile potting soil at a ratio of 1:14 (50 gm inoculum:700 gm sterile potting soil). This resulted in two different treatments for each genotype: high-biomass Ler, low-biomass Ler, high-biomass Cvi, and low-biomass Cvi. This process was repeated for the remaining seven generations. Starting with the third generation, plants were grown in uninoculated sterile potting soil as a control for random environmental variation (REV). In generation 6, all plants died for unknown reasons. The experiment was then continued using soil from generation 5. Soil sourced from the top six to ten plants and the bottom six to ten plants based on above-ground biomass was used as the inoculum to replant generation 6.

2.2. DNA Extraction and 16S rRNA Amplicon Library Prep

Microbial DNA was isolated from the frozen rhizosphere samples using the Machery-Nagel Nucleospin Soil DNA extraction kit (MACHEREY-NAGEL Inc., Allentown, PA, USA). Approximately 0.1 g of the rhizosphere soil sample was used for DNA isolation. All samples were diluted to 5 ng ul−1 with PCR-grade water. The 16S rRNA gene was amplified from the isolated DNA samples in triplicate in 96-well PCR plates. The PCR primers used for the 16S rRNA V4 region were 515F (5′-GTGYCAGCMGCCGCGGTAA-3′) [34] and 806R (5′-GGACTACNVGGGTWTCTAAT-3′) [35] for downstream paired-end Illumina (Illumina, Inc., San Diego, CA, USA) barcoded sequencing [36]. The PCR cycling conditions were as follows: 94 °C for 3 min; 25 cycles of 94 °C for 45 s; 50 °C for 60 s; 72 °C for 90 s; and final elongation at 72 °C for 10 min. The triplicate amplified samples were pooled, then purified and normalized using the SequalPrep™ Normalization Plate Kit (Invitrogen Corporation, Carlsbad, ON, Canada). Finally, multiplexed paired-end sequencing was carried out on the Illumina MiSeq (SY-410-1003) platform using earth microbiome project (EMP) primers [37].

2.3. Sequence Data Analysis and Statistics

The paired-end sequences obtained from the Illumina MiSeq were demultiplexed with QIIME 2 and converted into individual sequence fastq files for each sample [36,38]. The rest of the sequence processing was carried out in R (Version 4.2.2) using the DADA2 package (Version 1.24.0) [39]. The reads were processed in R using the following command in DADA2: ‘filterAndTrim(fnFs, filtFs, fnRs, filtRs, truncLen = c (145,145), minLen = 50, maxN = 0, maxEE = c (2,2), truncQ = 2, rm.phix = TRUE, compress = TRUE, multithread = TRUE)’. De novo OTU (operational taxonomic units) picking was performed using DADA2 (37), which resolves amplicon sequence variants (ASVs) to single-nucleotide differences. Chimeras were removed. Taxonomy was assigned using a 100% cut-off rate for species-level identification with the (2022) GTDB 16S rRNA reference database [40]. The phylogenetic tree was constructed using IQTree [41] with 1000 ultrafast bootstraps, and Modelfinder was used to identify the best model. In this study, the best-fit model was SYM + R10, chosen according to BIC. The data comprising the OTU table, phylogenetic tree, taxonomy table, and sample meta-data were then parsed through R studio using the phyloseq package [42]. Here, the processing removed OTUs identified as mitochondrial or chloroplast, which had less than five reads across all samples, and samples that had fewer than 8000 reads. The packages phytools, phyloseq, microbiome, and ggplot2 were used for further data analysis. All sequence data and metadata can be downloaded from PRJNA1010286.

2.4. Diversity Analyses

Alpha diversity was computed using two metrics, Observed and the Shannon diversity index. Beta diversity was calculated using weighted and unweighted UniFrac, and subsequently, principle coordinate analysis (PCoA) plots were constructed using the first two axes, which explained the most variance. Linear mixed effect models, via the lmerTest package in R, were used to determine changes in alpha and beta diversity over time (lmer (DiveristyMetric, Genotype*Generation*Biomass Treatment, (1|Generation)), with generation as a fixed and random effect. Pairwise comparisons were determined with emmeans and adjusted p-values were obtained with Tukey’s test.

2.5. Adonis

PERMANOVA was carried out for each generation using the adonis2 test. Both weighted and unweighted UniFrac distance matrices were used to account for the abundance, presence/absence, and phylogeny of the microbiome. The following model was used ‘adonis2 (distance.matrix ~ Biomass Treatment + Genotype, data = meta, permutations = 999)’.

2.6. Neutral Model

To determine the impact of neutral processes like drift and dispersal, or by deterministic selective forces, such as plant genotype and biomass treatments, on microbial community assembly, we carried out an analysis that was described by Burns et al. [43] to fit microbiome abundance data to the neutral model for prokaryotes from Sloan et al. [44]. The following design was used: sncm.fit (spp = generation (n)_Biomass Treatment, pool = generation (n), stats = T), where n is the nth generation and ‘Biomass Treatment’ represents the high- or low-biomass samples.

2.7. Differential Abundance

There are several challenges in estimating differentially abundant (DA) taxa in microbiome data. These include high variability in their abundance, zero-inflated data, and the compositional nature of the data. Nearing et al. [45] compared several methods across 38 datasets and found that different methods often identified different sets of DA taxa. The DA taxa were estimated using ANCOMBC2 [46], a conservative and robust approach. ANCOMBC2 incorporated bias correction, effectively addressing sampling-specific and sequencing biases present in the data. This feature ensures that the analysis is not skewed by any systematic errors introduced during the sampling or sequencing processes. It also conducts a sensitivity analysis for the pseudo-count addition to assess the impacts of different pseudo-count values on zero counts for each taxon. The analysis was run on Ler and Cvi samples, from OTU to family taxonomic levels, as follows: (fixed effect = (Biomass Treatment + Genotype + Generation), group = Biomass Treatment).

2.8. Spearman Correlation for Differentially Abundant OTUs

The ‘cor’ function from the stats package (version 4.2.3) was employed to calculate Spearman’s correlation coefficient between differentially abundant OTUs in the high-biomass-treatment samples of Ler. This was carried out for samples from both the early phase (generation 1 and generation 2) and the late phase (generation 6, generation 7, and generation 8) of the experiment. The results were then visualized using the corrplot package (version 0.92). Correlations were computed separately for the early and late phases of the experiment, enabling us to observe changes in the correlation coefficients at different stages of the experiment.

3. Results

3.1. Changes in Above-Ground Plant Biomass

The ecosystem selection experiment was carried out with respect to the above-ground plant biomass in two different accessions of Arabidopsis thaliana—Ler and Cvi—for eight generations. The two selected treatments were characterized by high and low above-ground biomass. The absolute values for the mean biomass of plants (n = 50) from each genotype and treatment changed substantially through the course of the eight generations of the experiment (Supplementary Tables S10 and S11). Despite the stochastic fluctuations in biomass caused by random environmental variations (REVs), the high-biomass-treatment plants were always the same as or greater than the low-biomass-treatment plants for both genotypes (Figure 1a,b). Due to the significant drop in biomass between generations 1 and 2, an uninoculated control in sterile potting soil was planted in all subsequent generations to serve as a reference for REV. The uninoculated control shows that the drastic drop in biomass after generations 1 and 2 was not a result of the inoculated soil, but was due to other subtle, but crucial, environmental factors that could not be held constant. The above-ground biomass for each generation in terms of deviations from the mean was plotted to give a clearer representation of the difference in phenotype seen in every generation (Figure 2a,b). Significant differences between the low- and high-biomass treatments become apparent from generation 4 onwards. Plants of all treatments in generation 6 died, resulting in a sharp dip, which necessitated repeating that generation (described in Section 2).

3.2. Microbial Community Composition

After preprocessing the 16S rRNA metagenomic amplicon reads in DADA 2, a total of 1,897,732 reads with an average of 12,822 reads per sample were obtained. Singletons and chimeras were removed during pre-processing. The most well-represented phyla were Proteobacteria (36.2%), Bacteroidetes (10.3%), Planctomycetes (7.8%), and Actinobacteria (7.6%) (Supplementary Figure S2). The most well-represented classes identified in the dataset were Alphaproteobacteria (24.7%), Gammaproteobacteria (11.4%), Bacteroidia (10.1%), and Verrucomicrobiae (6.2%) (Figure 3). This relative abundance is typical of microbiome data, heavily weighted by a few abundant groups and a long tail of rare taxa.

3.3. Trends in Alpha Diversity

We used two metrics to assess alpha diversity (Figure 4). Observed (F7,160 = 333.054, p < 0.001) estimated species richness and showed a significant decrease in the number of OTUs (1477 to ~580) between generation 1 and generation 5, with less pronounced changes in diversity after generation 5. The Shannon diversity index (F7,160 = 82.23, p < 0.001), which was more sensitive to the difference in abundance, exhibited a comparable trend. Both metrics demonstrated a rapid decline in diversity between generations 1 and 5, with a less pronounced decrease after generation 5. A negative exponential function was fitted to the Observed metric (R2 = 0.81, p < 0.01) and the Shannon metric (R2 = 0.58, p < 0.01). The Observed metric displayed a better fit, indicating the increased stability of the microbiome in later generations. In addition, a linear mixed-effect model showed significant interactions between genotype and biomass treatment (F = 8.71, p < 0.01); genotype and generation (F =3, p < 0.01); and genotype, biomass treatment, and generation (F = 3, p < 0.01) for Shannon (Supplementary Table S5). No significant interaction terms were found for Observed, which did not consider the abundance of OTUs in its calculations.

3.4. Beta Diversity and Principal Coordinate Analysis

Principal coordinate analysis (PCoA) of the weighted UniFrac distances was conducted. They were plotted with the first two axes that captured the most variance in the data, from 73% in generation 1 to more than 90% in generation 8 (Figure 5). The UniFrac distance considers the evolutionary relationships between taxa, which makes it more biologically meaningful than other distance metrics that do not consider the evolutionary history. Weighted UniFrac is more sensitive to differences in the abundance of taxa. The microbiomes associated with the biomass treatments and genotypes were very similar in generation 1, but began to diverge by generation 2. By generation 3, the microbiomes of the Ler and Cvi genotypes had diverged, and by generation 5, the high- and low-biomass treatments appeared to have further distinguished the microbial communities. The divergence of the microbial community in generations 4 and 5 caused by the high/low biomass treatments also aligned with the onset of notable differences in the above-ground plant biomass within these generations.
The complex interplay between genotype and biomass treatment was modeled using a PERMANOVA test with both weighted and unweighted UniFrac distances to account for differences in abundance and presence or absence of OTUs, respectively. The resulting R2 values for genotype and biomass treatment and residuals as a proxy for REVs were plotted (Figure 6a,b). For both the weighted and the unweighted UniFrac, the variance in the microbial community explained by both genotype and biomass treatment increased significantly over the eight generations. Despite observing a decrease in the variance accounted for by the residuals, which could potentially signify a decline in the influence of REVs, they accounted for ~50 percent of the variability within the microbial community in generation 8. For all metrics, weighted UniFrac exhibited more fluctuations compared to the unweighted UniFrac. The weighted UniFrac analyses found that genotype (F1,39 = 1.55, p < 0.05) and biomass treatment (F1,39 = 1.33, p < 0.001) explained 5.7% and 13.7%, respectively, of the dissimilarity between microbiomes in generation 1. By generation 8, the genotype (F1,14 = 5.14, p < 0.05) and biomass treatment (F1,14 = 5.48, p < 0.001) explained 24% and 22%, respectively, of the dissimilarity between microbiomes (Supplementary Tables S6 and S7). For unweighted UniFrac, the influence of genotype on the microbial community increased from 3.3% in generation 1 (F1,39 = 1.33, p < 0.01) to 26% in generation 8 (F1,39 = 5.13, p < 0.001), and for biomass, it increased from 3.8% in generation 1 (F1,14 = 1.55, p < 0.001) to 13% in generation 8 (F1,14 = 2.58, p < 0.05) (Supplementary Tables S8 and S9).

3.5. Sloan Neutral Model

Many forces can alter a microbial community’s structure and dynamics. For our study, these can be divided into selective forces, such as the plant genotype being colonized and the biomass selection process, or neutral, stochastic forces like dispersal that are innate in any experiment or environment. The Sloan neutral model was fit to the data to assess the importance of selective versus neutral drivers of change in the experiment. If neutral processes were the driving force in the microbial community assembly, the null hypothesis would be that all generations of the experiment would fit the model equally well. The model fit to neutrality represented by the R2 value decreased from 0.795 in generation 1 to 0.59 in generation 8, demonstrating the increasing importance of selective drivers over the course of the experiment (Figure 7).

3.6. Differential Abundance and Correlation Analyses

Differentially abundant (DA) taxa that distinguished the high-biomass treatments from the low-biomass treatments were determined using ANCOMBC2 at the family level. Taxa lacking classification at the family level were denoted by the preceding identifiable taxonomic tier. The OTUs that were present in higher (red) and lower (blue) abundances in the high biomass treatment were identified (Figure 8a,b). Among the DA taxa in this study, several are known to benefit plants. Bacteria within the class Bacilli promote plant growth [47]. Members of the class Gemmatimonadetes have been shown to exhibit a positive association with vegetation restoration, plant richness, and soil nutrients [48]. Cytophagacea are chemoorganotrophs, important for remineralizing organic materials into micronutrients. They could support both mycelial growth and plant nutrition [49].
Enterobacteriaceae, Paenibacillaceae, and JACDCH01 were all present in higher abundances in the high-biomass treatments of both Ler and Cvi. Most Paenibacillaceae members predominantly inhabited soil, frequently in close association with plant roots [50]. These rhizobacteria are known to play a significant role in enhancing plant growth, and they possess potential applications in agriculture. Enterobacteriaceae has been reported to enhance plant growth [51,52,53].
Additionally, Order_NS11-12g, BJHT01, JACDCH01, and UBA6156 were more abundant in Cvi under high-biomass conditions, but less abundant in Ler under high-biomass conditions. The majority of these were uncultivated or candidate taxa. Two alternative families in the order Rhizobiales appeared to be differentially associated with the two plant genotypes: Labraceae with Ler and Rhizobiaceae with Cvi [54]. These could indicate genotype-specific interactions with different members of the microbial community.
The relative abundances of all OTUs that exhibited higher abundance (enriched) in samples from the high-biomass treatment in Ler and Cvi are depicted in Figure 8c and Figure 8d, respectively. These OTUs consistently exhibited an increasing trend of relative abundance across successive generations. A linear model was employed to test this trend, resulting in R-squared values of 0.57 (p < 0.01) for Ler and 0.73 (p < 0.01) for Cvi. Furthermore, OTUs that were lower in abundance (depleted) in samples from the high-biomass treatment of Ler and Cvi are depicted in Figure 8e and Figure 8f, respectively. These OTUs persisted in low abundance across all generations, suggesting a lack of selective pressure on this group of microbes in the samples which received the high-biomass treatment.
We then conducted a Spearman’s correlation analysis using 58 OTUs identified as enriched in high-biomass-treatment samples of Ler (Figure 9a,b). This analysis was carried out for samples from both the early phases (generations 1 and 2) and the late phases (generations 6, 7, and 8) of the experiment. Our observations indicated that, during the early phases, the associations between OTUs were comparatively weaker than during the later phases, suggesting a trend of OTUs being selected in groups rather than individually. Additionally, we presented specific examples of OTUs that showed significant differences in the number of pairwise positive and negative associations between the early and late phases of the experiment (Figure 9c). This visualization highlights trends spanning both the early and late phases of the experiment, with certain OTUs demonstrating an increase in positive pairwise associations and a decrease in negative pairwise associations as generations progressed. Notable examples include OTU1606_Chitinophagaceae, OTU162_Bdellovibrionaceae, OTU94_Class_Bacilli, and OTU359_Order_AKYH767. Conversely, others displayed the opposite trend, such as OTU763_Fimbriimonadaceae, OTU142_Sphingobacteriaceae, and OTU1005_Labraceae.

4. Discussion

We applied artificial ecosystem selection to the eight generations in Arabidopsis thaliana Ler and Cvi to select soil microbiomes associated with the higher or lower biomass of the host. In contrast to some previous studies [26,32], we did not fertilize plants, thus maintaining nutrient limitation and thereby promoting the interaction of the plant with the microbiome [55]. Over the course of eight generations, a response to the selection, most apparent after generations 4 and 5, was evident in both the gradually shifting microbiome composition and the plant biomass. The microbiome selection process noticeably influenced the plant biomass despite large phenotypic variation from one generation to the next. Stochasticity in growth through generations may be attributed to random environmental variations (REVs), as indicated by Swenson et al. [27]. This was observed in the uninoculated sterile potting soil reference, which showed similar patterns of variability to the inoculated treatments across generations (Figure 1). Ecological variability among generations, even under the conditions of a controlled growth chamber, is a common characteristic of similar studies [56], and may be due to minute fluctuations in the growth chamber conditions or the batch and age of the potting soil.
We observed a gradual decline in microbial species richness during the initial generations of the experiment. This was evident from the observed OTU counts, as well as from the Chao1 and Shannon diversity indices. This pattern is consistent with findings from other studies [31,32,33], highlighting the impact of selection pressures as the microbial communities adapted to the host plant’s environment. However, in our study, which continued for about twice the number of generations as these earlier studies, we observed a stabilization of richness and alpha diversity in the later generations. This was further substantiated by the stronger fit to the negative exponential function when compared to a linear model, indicating constant and proportional changes in alpha diversity and a trend of the rhizosphere microbiome towards reduced diversity loss, accompanied by heightened stability, over eight generations. Previous studies have suggested that the stability in alpha diversity metrics in later generations may possibly be due to the retention of microbiome community members through robust selection and improved fitness with the host [57,58,59]. Over time, less fit individuals are outcompeted by microbes that have better fitness in the host environment [59,60].
To further understand the complex interplay between the forces of directed selection (plant genotype and biomass treatment), we conducted a PCoA. Over successive generations, the results demonstrated a strengthening effect of the plant genotypes and biomass treatment on the microbial community (Figure 5). This result was modeled using PERMANOVA for both weighted and unweighted UniFrac distances. Both the weighted and unweighted UniFrac showed marked increases in the proportion of variance explained by genotype and biomass treatment, with the unweighted UniFrac exhibiting larger increases from generation to generation (Figure 6). Both metrics considered phylogenetic relatedness. Weighted UniFrac considers the absolute abundance of OTUs, and is generally used to study changes in microbial community structure. Unweighted UniFrac considers only the presence/absence of OTUs, and is generally used to study changes in microbial community composition. The results suggest that the variance in the microbiome’s composition (presence/absence of OTUs), as explained by genotype, consistently increases, whereas the variance in the microbiome structure (considers the abundance of OTUs), as explained by genotype, greatly fluctuates with changes in abundant taxa.
Here, we would like to highlight that pronounced shifts in biomass between the high/low biomass treatments, alpha diversity, and beta diversity all occur around generations 4 and 5 (Figure 1 and Figure 2). This marks the point at which the changes in alpha diversity stabilize (Figure 4) and the divergence between microbial communities by high/low biomass treatments in the PCoA plots becomes more pronounced (Figure 5). It is possible that, over the first four to five generations, the initial microbial community underwent a period of restructuring before it stabilized and formed four distinct communities under selection by high/low biomass treatments and by plant genotype. We propose that a complex interaction between plants and their associated microbiome ensued, where differences in root exudation patterns between the two genotypes presumably established associations with microbes in the early generations. Simultaneously, during this period, the biomass treatment likely promoted host-specific, microbe-mediated interactions that modulated plant biomass. This restructuring of the microbial community was driven by both structural (abundance of OTUs) and compositional (presence/absence of OTUs) effects in the microbiome (Figure 6a,b). Differences in the abundance of taxa driven by the biomass treatment and genotype selection pressures were a larger contributor to the generation-to-generation variation during the experiment.
Changes in the assembly of microbial communities can arise from either selective pressure, like the plant’s genotype or the biomass treatment, as observed in this experiment, or stochastic processes, like minute changes in the growth chamber humidity or microenvironment. To gain deeper insights into the influence of selection on microbial community assembly, a neutral model was fit to the data [43,44]. Interestingly, this revealed a progressive decline in fit to neutrality, indicated by decreasing R2 values and increasing AIC values (Figure 7). This suggests an increasing influence of selective forces such as biomass treatment and genotype over multiple generations of ecosystem selection. This has also previously been observed by Morella et al. [33].
Previous studies in Arabidopsis thaliana have often detected only a weak genotypic effect on the rhizosphere microbiome [18,23,61]. However, a majority of these studies are focused on a single growth cycle. Insights from plant–soil feedback research emphasize that the interactions between plants and their microbiomes can be modified by consecutive growth cycles of the same plant species in its respective soil [24,25]. Recent research has presented conflicting findings on the influence of the host genotype [32,33,62]. The different experimental designs, communities being assessed, and the length of these experiments, alongside the high stochastic environmental variation that is clearly associated with these microbial studies, likely account for the conflicting findings. In our investigation, we found that during the initial stages, genotype and biomass treatments had modest, but significant, impacts. Over time, the plant genotype and biomass treatments exerted an increasing influence. Together, these selective forces explained ~40% of the variation in the microbial community’s composition in these later generations.
A key aim of microbiome engineering has been to develop host-specific microbial communities that impart lasting beneficial effects to the plant host [26,28,31,63]. Herein, we show the enrichment of some common taxa in high-biomass treatments, but also some genotype-specific changes. Despite starting from the same soil, only three common families were enriched in the high-biomass treatment samples of both the Ler and Cvi genotypes. In addition, well-known plant-growth-promoting rhizobacteria were enriched in a genotype-specific manner for both the Ler and Cvi samples that underwent the high-biomass treatment. In the order Rhizobiales, Labraceae was enriched in Ler, while Rhizobiaceae was enriched in Cvi [54]. Moreover, in high-biomass conditions, Cvi had increased levels of the taxa Order_NS11-12g, BJHT01, JACDCH01, and UBA6156, while Ler had decreased levels of these. These data suggest genotype-specific plant genetic control of the rhizosphere microbiome.
To further explore how the abundance of these OTUs changed across successive generations, we grouped all the DA OTUs that were enriched or depleted in the samples with the high-biomass treatment and calculated the sum of their relative abundances. The OTUs enriched in the samples of high biomass treatment displayed an increasing trend of relative abundance over subsequent generations (Figure 8c,d). A linear model was fitted to these data with high R2 values, reinforcing the notion that the genotype and biomass treatments exerted selective pressure on entire microbial groups rather than individual species. Conversely, the OTUs that were depleted in the samples which underwent high-biomass treatment consistently remained in a state of low relative abundance throughout the entirety of the experiment (Figure 8e,f). This observation may indicate that the selective pressures exerted by the genotype and biomass treatments do not influence this particular group of microbes.
The results of the correlation analysis revealed a strong directional shift in the relationships between the OTUs in the later phase of the experiment compared to the early phase of the experiment. This analysis provides further empirical evidence that the OTUs were selected in groups. This is particularly evident for members of Chitinophagaeae, which showed a strong increase in pairwise positive associations and a concomitant decrease in negative associations. Strengthening of microbial associations has been observed in previous succession studies on different plant species, including Arabidopsis thaliana [8,64,65]. Luo et al. also found stronger positive associations in the rhizosphere, highlighting the importance of this zone in promoting stronger microbial associations [64].
Taken collectively, these findings suggest that the presence of diverse selection pressures influences the rhizosphere microbiome. Initially, it seems that neutral processes seem to play a major role in determining the structure and composition of the microbiome (Figure 7). As generations progress, we observed an increase in the strength of directional selection, possibly driven by genotype and biomass treatment, resulting in the enrichment of specific groups of microbes, including some PGPRs (Figure 8c,d). The persistence of depleted microbes in the high-biomass-treatment samples suggests the possibility of negative frequency-dependent selection preserving them at a low abundance. On the other hand, it is possible that this group did not experience any observable selection pressure.
Multiple studies have shown variations in the structure of the rhizosphere microbiome, even within closely related plant genotypes, highlighting the importance of genotype-specific root exudates in forming associations with the corresponding microbiome [66,67,68,69,70,71,72]. While previous studies have observed qualitative distinctions in the root exudate profiles of Ler and Cvi, these profiles lack quantitative characterization [67]. Currently, our understanding is limited to some characterization of root exudates for the Ler [67,68,73], and to the best of our knowledge, there has been no in-depth characterization of Cvi’s root exudate profile. Closer associations in the rhizosphere are formed between plants and bacterial partners capable of metabolizing specific compounds. In a previous Arabidopsis thaliana study, it was observed that rhizosphere-enriched microbes, such as Terracoccus sp. 273, excelled in utilizing SA as a carbon source. Conversely, microbes like Mitsuaria sp. were unable to do so and were found in lower abundance [61]. Other plants, such as maize, release benzoxazinoids (BXs) in their root exudates, and there is a negative correlation in terms of abundance between Glomeromycota members and BXs [24,74,75]. It is important to emphasize that microbes are also able to reciprocally influence the root exudation patterns of plants. Pea plants treated with Pseudomonas aeruginosa PJHU15 and Bacillus subtilis BHHU100 showed a significant increase in gallic acid and other phenolics compared to control plants [76]. Furthermore, in response to Pseudomonas syringae pv tomato infection, Arabidopsis thaliana roots shift their exudation patterns, releasing more malic acid. Malic acid is a chemoattractant, selectively inducing Bacillus subtilis to bind to the infected plant’s roots and, thus, improving plant disease resistance [77]. Future explorations should prioritize profiling root exudates, along with comparative metagenomics, metatranscriptomics, and metabolomics studies. This will play a role in decoding chemical crosstalk, which could eventually foster reproducible interactions between plants and microbes, leading to an overall enhancement in plant fitness.
The relationship between plants and microorganisms in the rhizosphere involves complex and diverse interactions which influence crucial ecological and physiological processes. These interactions can be mutualistic, competitive, or antagonistic. Many biotic and abiotic factors act in concert to influence the dynamic rhizosphere microbiome. Gaining insight into how these factors influence the assembly of microbial communities is crucial for obtaining targeted and long-lasting advantages for plants. Nevertheless, establishing a strong and host-specific microbial community that consistently provides beneficial effects continues to be a persistent challenge. We have shown that, despite stochastic fluctuation due to REVs, it is possible to select for microbial communities that impact biomass in a genotype-specific manner within four generations, addressing this key challenge in microbial community engineering [31,63]. The rhizosphere microbiome that evolves under plant-mediated selection has the potential to offer improved survivability and efficacy when applied as an inoculum to the plant [78,79]. Multi-species ecosystem selection provides a perspective on microbiome assembly that indicates the holobiont, and the evolving relationship between the host and the microbiome contributes to emerging properties beyond those predicted by the host’s genotype or the initial composition of the microbial community. This study enhances our understanding of the temporal dynamics involved in microbial assembly in the rhizosphere, and has implications for sustainable agriculture, evolution, and ecology.

5. Conclusions

The detrimental impact of agrochemicals has emphasized the importance of adopting sustainable agricultural practices [80,81]. Microbial formulations of PGPRs as inoculants may boost soil fertility, promote favorable plant phenotypes, and reduce the harmful effects of agrochemicals. However, numerous challenges remain in terms of transitioning from the lab to the field—in particular, the survivability and efficacy of these individual strains [82,83,84]. Ecosystem selection overcomes these challenges by selecting for entire groups of microbes rather than individuals, mitigating the potential loss of function due to the extinction of essential microbes. We utilized multi-generational ecosystem selection with respect to above-ground plant biomass for eight generations in two common inbred genotypes of Arabidopsis thaliana: Ler and Cvi. The composition of the resulting rhizosphere microbiome was shaped by a complex interplay between environmental factors, plant genotype, and biomass treatment. In the first few generations of the experiment, the genotype and biomass treatments played a modest role in shaping the microbial community. Over time, the plant genotype and biomass treatments exerted increasing influence, explaining ~40% of the difference in the microbial community composition. We successfully cultivated microbiomes that were associated with either increased or reduced plant biomass, and an enrichment of known plant-growth-promoting rhizobacteria (PGPR) was observed with the high-biomass treatment. The microbiomes identified in this study alter plant biomass and growth, presumably by altering the availability of nutrients in the soil or by mitigating biotic/abiotic stressors.
Translating the findings of such studies from a lab to a field setting requires further optimization of the inoculum transfer [28,85]. For instance, Wright et al. made use of artificial selection to select for active degraders of chitin by timing the harvest window to collect inoculum at the same time that the desired trait (chitinase activity) was at its peak and not at a pre-defined incubation time [86]. However, the dynamic nature of rhizosphere microbiomes, where member abundances constantly shift based on the host plant’s genotype and developmental stage, as well as environmental factors, adds complexity to the challenge of identifying the optimal time point for transfer in plant microbiome studies. Future research should focus on the meticulous refinement of transfer timing, storage methods, and application techniques for the selected microbial inoculum to achieve sustainable agricultural practices and improved soil health through ecosystem selection.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11122932/s1, Figure S1: Rhizotube schematic; Figure S2: Taxonomic distribution of the microbial community represented in terms of the relative abundance at the phylum level. Table S1: Statistical comparisons of biomass for Ler; Table S2: Statistical comparisons of biomass for Cvi; Table S3: ANOVA for Observed alpha diversity index; Table S4: ANOVA for Shannon alpha diversity index; Table S5: Linear mixed-effect model for Shannon diversity index; Table S6: PERMANOVA-weighted UniFrac distances for generation 1; Supplementary Table S7: PERMANOVA-weighted UniFrac distances for generation 8; Supplementary Table S8: PERMANOVA-unweighted UniFrac distances for generation 1; Table S9: PERMANOVA-unweighted UniFrac distances for generation 8. Table S10. Mean and Standard deviation for biomass of Ler for all generations. The raw value of mean biomass of Ler are referred to as ‘mean_BM’ and the standard deviation is referred to as ‘sd_bm’. Table S11. Mean and Standard deviation for biomass of Cvi for all generations.

Author Contributions

N.S.: conceptualization, methodology, data analysis and visualization, data collection, writing—original draft, writing—reviewing and editing. P.S.: methodology, data analysis, writing—original draft, writing—reviewing and editing. T.C.M.: data collection. R.K.: resources, funding, writing—reviewing and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was made possible in part by funding from NSF (IOS 0847691).

Data Availability Statement

All sequence data and metadata can be downloaded from PRJNA1010286.

Acknowledgments

We would like to express our sincere gratitude to the following individuals for their advice at various stages of this project: David Sloan Wilson, Jenny Kao-Kniffin, Deepa Agashe, Gergely Maróti, Douglas Woodhams, Emeritus Michael Shiaris, Karen Ricciardi, Patrick Kearns, Kourosh Zarringhalam, David Degras-Valabregue, Liam Revell, Jay T. Lennon, Sneha Berry, Molly Bletz, Brandon LaBumbard, Cassandra Lyons, Haohui Duan, and Michael Schmelz. This work was made possible in part by funding from NSF (IOS 0847691).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gouda, S.; Kerry, R.G.; Das, G.; Paramithiotis, S.; Shin, H.S.; Patra, J.K. Revitalization of plant growth promoting rhizobacteria for sustainable development in agriculture. Microbiol. Res. 2018, 206, 131–140. [Google Scholar] [CrossRef]
  2. Jones, D.L.; Nguyen, C.; Finlay, R.D. Carbon flow in the rhizosphere: Carbon trading at the soil-root interface. Plant Soil 2009, 321, 5–33. [Google Scholar] [CrossRef]
  3. Dennis, P.G.; Miller, A.J.; Hirsch, P.R. Are root exudates more important than other sources of rhizodeposits in structuring rhizosphere bacterial communities? FEMS Microbiol. Ecol. 2010, 72, 313–327. [Google Scholar] [CrossRef]
  4. Hartmann, A.; Schmid, M.; van Tuinen, D.; Berg, G. Plant-driven selection of microbes. Plant Soil 2009, 321, 235–257. [Google Scholar] [CrossRef]
  5. Chaparro, J.M.; Sheflin, A.M.; Manter, D.K.; Vivanco, J.M. Manipulating the soil microbiome to increase soil health and plant fertility. Biol. Fertil. Soils 2012, 48, 489–499. [Google Scholar] [CrossRef]
  6. Marschner, P.; Crowley, D.; Yang, C.H. Development of specific rhizosphere bacterial communities in relation to plant species, nutrition and soil type. Plant Soil 2004, 261, 199–208. [Google Scholar] [CrossRef]
  7. Garbeva, P.; van Veen, J.A.; van Elsas, J.D. Microbial diversity in soil: Selection microbial populations by plant and soil type and implications for disease suppressiveness. Annu. Rev. Phytopathol. 2004, 42, 243–270. [Google Scholar] [CrossRef] [PubMed]
  8. Tkacz, A.; Cheema, J.; Chandra, G.; Grant, A.; Poole, P.S. Stability and succession of the rhizosphere microbiota depends upon plant type and soil composition. ISME J. 2015, 9, 2349–2359. [Google Scholar] [CrossRef]
  9. Berg, G.; Smalla, K. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 2009, 68, 1–13. [Google Scholar] [CrossRef]
  10. Jin, H.; Yang, X.; Liu, R.; Yan, Z.; Li, X.; Li, X.; Su, A.; Zhao, Y.; Qin, B. Bacterial community structure associated with the rhizosphere soils and roots of Stellera chamaejasme L. along a Tibetan elevation gradient. Ann. Microbiol. 2018, 68, 273–286. [Google Scholar] [CrossRef]
  11. Barbosa Lima, A.; Cannavan, F.S.; Navarrete, A.A.; Teixeira, W.G.; Kuramae, E.E.; Tsai, S.M. Amazonian Dark Earth and Plant Species from the Amazon Region Contribute to Shape Rhizosphere Bacterial Communities. Microb. Ecol. 2015, 69, 855–866. [Google Scholar] [CrossRef]
  12. Marques, J.M.; da Silva, T.F.; Vollu, R.E.; Blank, A.F.; Ding, G.C.; Seldin, L.; Smalla, K. Plant age and genotype affect the bacterial community composition in the tuber rhizosphere of field-grown sweet potato plants. FEMS Microbiol. Ecol. 2014, 88, 424–435. [Google Scholar] [CrossRef]
  13. Bodenhausen, N.; Bortfeld-Miller, M.; Ackermann, M.; Vorholt, J.A. A Synthetic Community Approach Reveals Plant Genotypes Affecting the Phyllosphere Microbiota. PLoS Genet. 2014, 10, e1004283. [Google Scholar] [CrossRef]
  14. Zancarini, A.; Mougel, C.; Voisin, A.S.; Prudent, M.; Salon, C.; Munier-Jolain, N. Soil Nitrogen Availability and Plant Genotype Modify the Nutrition Strategies of M. truncatula and the Associated Rhizosphere Microbial Communities. PLoS ONE 2012, 7, e47096. [Google Scholar] [CrossRef]
  15. Wei, S.; Liu, B.; Ni, K.; Ma, L.; Shi, Y.; Leng, Y.; Zheng, S.; Gao, S.; Yang, X.; Ruan, J. Rhizosphere Microbial Community Shows a Greater Response Than Soil Properties to Tea (Camellia sinensis L.) Cultivars. Agronomy 2023, 13, 221. [Google Scholar] [CrossRef]
  16. Yuan, J.; Chaparro, J.M.; Manter, D.K.; Zhang, R.; Vivanco, J.M.; Shen, Q. Roots from distinct plant developmental stages are capable of rapidly selecting their own microbiome without the influence of environmental and soil edaphic factors. Soil Biol. Biochem. 2015, 89, 206–209. [Google Scholar] [CrossRef]
  17. Houlden, A.; Timms-Wilson, T.M.; Day, M.J.; Bailey, M.J. Influence of plant developmental stage on microbial community structure and activity in the rhizosphere of three field crops. FEMS Microbiol. Ecol. 2008, 65, 193–201. [Google Scholar] [CrossRef]
  18. Lundberg, D.S.; Lebeis, S.L.; Paredes, S.H.; Yourstone, S.; Gehring, J.; Malfatti, S.; Tremblay, J.; Engelbrektson, A.; Kunin, V.; Rio, T.G.D.; et al. Defining the core Arabidopsis thaliana root microbiome. Nature 2013, 488, 86–90. [Google Scholar] [CrossRef] [PubMed]
  19. Bulgarelli, D.; Rott, M.; Schlaeppi, K.; Ver, E.; Themaat, L.V.; Ahmadinejad, N.; Assenza, F.; Rauf, P.; Huettel, B.; Reinhardt, R.; et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 2012, 488, 91–95. [Google Scholar] [CrossRef] [PubMed]
  20. Coleman-Derr, D.; Desgarennes, D.; Fonseca-Garcia, C.; Gross, S.; Clingenpeel, S.; Woyke, T.; North, G.; Visel, A.; Partida-Martinez, L.P.; Tringe, S.G. Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 2016, 209, 798–811. [Google Scholar] [CrossRef] [PubMed]
  21. 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]
  22. Wagner, M.R.; Lundberg, D.S.; Del Rio, T.G.; Tringe, S.G.; Dangl, J.L.; Mitchell-Olds, T. Host genotype and age shape the leaf and root microbiomes of a wild perennial plant. Nat Commun. 2016, 7, 12151. [Google Scholar] [CrossRef] [PubMed]
  23. Schlaeppi, K.; Dombrowski, N.; Oter, R.G.; Ver Loren Van Themaat, E.; Schulze-Lefert, P. Quantitative divergence of the bacterial root microbiota in Arabidopsis thaliana relatives. Proc. Natl. Acad. Sci. USA 2014, 111, 585–592. [Google Scholar] [CrossRef] [PubMed]
  24. Hu, L.; Robert, C.A.M.; Cadot, S.; Zhang, X.; Ye, M.; Li, B.; Manzo, D.; Chervet, N.; Steinger, T.; Van Der Heijden, M.G.; et al. Root exudate metabolites drive plant-soil feedbacks on growth and defense by shaping the rhizosphere microbiota. Nat. Commun. 2018, 9, 2738. [Google Scholar] [CrossRef] [PubMed]
  25. Van der Putten, W.H.; Bardgett, R.D.; Bever, J.D.; Bezemer, T.M.; Casper, B.B.; Fukami, T.; Kardol, P.; Klironomos, J.N.; Kulmatiski, A.; Schweitzer, J.A.; et al. Plant-soil feedbacks: The past, the present and future challenges. J. Ecol. 2013, 101, 265–276. [Google Scholar] [CrossRef]
  26. Panke-buisse, K.; Poole, A.C.; Goodrich, J.K.; Ley, R.E.; Kao-kniffin, J. Selection on soil microbiomes reveals reproducible impacts on plant function. ISME J. 2014, 9, 980–989. [Google Scholar] [CrossRef] [PubMed]
  27. Swenson, W.; Wilson, D.S.; Elias, R. Artificial ecosystem selection. Proc. Natl. Acad. Sci. USA 2000, 97, 9110–9114. [Google Scholar] [CrossRef] [PubMed]
  28. Mueller, U.G.; Juenger, T.E.; Kardish, M.R.; Carlson, A.L.; Burns, K.M.; Edwards, J.A.; Smith, C.C.; Fang, C.C.; Des Marais, D.L. Artificial Selection on Microbiomes To Breed Microbiomes That Confer Salt Tolerance to Plants. mSystems 2021, 6, e01125-21. [Google Scholar] [CrossRef]
  29. Lau, J.A.; Lennon, J.T. Rapid responses of soil microorganisms improve plant fitness in novel environments. Proc. Natl. Acad. Sci. USA 2012, 109, 14058–14062. [Google Scholar] [CrossRef]
  30. terHorst, C.P.; Lennon, J.T.; Lau, J.A. The relative importance of rapid evolution for plant-microbe interactions depends on ecological context. Proc. R. Soc. B 2014, 281, 20140028. [Google Scholar] [CrossRef]
  31. Jochum, M.D.; McWilliams, K.L.; Pierson, E.A.; Jo, Y.K. Host-mediated microbiome engineering (HMME) of drought tolerance in the wheat rhizosphere. PLoS ONE 2019, 14, e0225933. [Google Scholar] [CrossRef]
  32. Cordovez, V.; Rotoni, C.; Dini-Andreote, F.; Oyserman, B.; Carrión, V.J.; Raaijmakers, J.M. Successive plant growth amplifies genotype-specific assembly of the tomato rhizosphere microbiome. Sci. Total Environ. 2021, 772, 144825. [Google Scholar] [CrossRef] [PubMed]
  33. Morella, N.M.; Weng, F.C.H.; Joubert, P.M.; Metcalf, C.J.E.; Lindow, S.; Koskella, B. Successive passaging of a plant-associated microbiome reveals robust habitat and host genotype-dependent selection. Proc. Natl. Acad. Sci. USA 2020, 117, 1148–1159. [Google Scholar] [CrossRef]
  34. Parada, A.E.; Needham, D.M.; Fuhrman, J.A. Every base matters: Assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 2016, 18, 1403–1414. [Google Scholar] [CrossRef] [PubMed]
  35. Apprill, A.; Mcnally, S.; Parsons, R.; Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 2015, 75, 129–137. [Google Scholar] [CrossRef]
  36. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef] [PubMed]
  37. Thompson, L.R.; Sanders, J.G.; McDonald, D.; Amir, A.; Ladau, J.; Locey, K.J.; Prill, R.J.; Tripathi, A.; Gibbons, S.M.; Ackermann, G.; et al. A communal catalogue reveals Earth’s multiscale microbial diversity. Nature 2017, 551, 457–463. [Google Scholar] [CrossRef]
  38. 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]
  39. 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]
  40. Chaumeil, P.A.; Mussig, A.J.; Hugenholtz, P.; Parks, D.H. GTDB-Tk v2: Memory friendly classification with the genome taxonomy database. Bioinformatics 2022, 38, 5315–5316. [Google Scholar] [CrossRef]
  41. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; Von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [PubMed]
  42. Bunn, A.G. A dendrochronology program library in R (dplR). Dendrochronologia 2008, 26, 115–124. [Google Scholar] [CrossRef]
  43. Burns, A.R.; Stephens, W.Z.; Stagaman, K.; Wong, S.; Rawls, J.F.; Guillemin, K.; Bohannan, B.J. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 2016, 10, 655–664. [Google Scholar] [CrossRef] [PubMed]
  44. Sloan, W.T.; Woodcock, S.; Lunn, M.; Head, I.M.; Curtis, T.P. Modeling taxa-abundance distributions in microbial communities using environmental sequence data. Microb. Ecol. 2007, 53, 443–455. [Google Scholar] [CrossRef] [PubMed]
  45. Nearing, J.T.; Douglas, G.M.; Hayes, M.G.; MacDonald, J.; Desai, D.K.; Allward, N.; Jones, C.M.; Wright, R.J.; Dhanani, A.S.; Comeau, A.M.; et al. Microbiome differential abundance methods produce different results across 38 datasets. Nat. Commun. 2022, 13, 342. [Google Scholar] [CrossRef] [PubMed]
  46. Lin, H.; Peddada, S.D. Analysis of compositions of microbiomes with bias correction. Nat. Commun. 2020, 11, 3514. [Google Scholar] [CrossRef]
  47. Radhakrishnan, R.; Hashem, A.; Abd Allah, E.F. Bacillus: A biological tool for crop improvement through bio-molecular changes in adverse environments. Front. Physiol. 2017, 8, 667. [Google Scholar] [CrossRef]
  48. Mujakić, I.; Piwosz, K.; Koblížek, M. Phylum Gemmatimonadota and Its Role in the Environment. Microorganisms 2022, 10, 151. [Google Scholar] [CrossRef]
  49. Dworkin, M.; Falkow, S.; Rosenberg, E.; Schleifer, K.H.; Stackebrandt, E. The Prokaryotes: A Handbook on the Biology of Bacteria: Archaea. Bacteria: Firmicutes, Actinomycetes; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  50. Grady, E.N.; MacDonald, J.; Liu, L.; Richman, A.; Yuan, Z.C. Current knowledge and perspectives of Paenibacillus: A review. Microb. Cell Factories 2016, 15, 203. [Google Scholar] [CrossRef]
  51. Ramesh, A.; Sharma, S.K.; Sharma, M.P.; Yadav, N.; Joshi, O.P. Plant Growth-Promoting Traits in Enterobacter cloacae subsp. dissolvens MDSR9 Isolated from Soybean Rhizosphere and its Impact on Growth and Nutrition of Soybean and Wheat Upon Inoculation. Agric. Res. 2014, 3, 53–66. [Google Scholar] [CrossRef]
  52. Kämpfer, P.; Ruppel, S.; Remus, R. Enterobacter radicincitans sp. nov., a plant growth promoting species of the family Enterobacteriaceae. Syst. Appl. Microbiol. 2005, 28, 213–221. [Google Scholar] [CrossRef]
  53. Khalifa, A.Y.Z.; Alsyeeh, A.M.; Almalki, M.A.; Saleh, F.A. Characterization of the plant growth promoting bacterium, Enterobacter cloacae MSR1, isolated from roots of non-nodulating Medicago sativa. Saudi J. Biol. Sci. 2016, 23, 79–86. [Google Scholar] [CrossRef]
  54. Harkes, P.; van Steenbrugge, J.J.M.; van den Elsen, S.J.J.; Suleiman, A.K.A.; de Haan, J.J.; Holterman, M.H.M.; Helder, J. Shifts in the Active Rhizobiome Paralleling Low Meloidogyne chitwoodi Densities in Fields Under Prolonged Organic Soil Management. Front. Plant Sci. 2020, 10, 1697. [Google Scholar] [CrossRef] [PubMed]
  55. Meier, M.A.; Xu, G.; Lopez-Guerrero, M.G.; Li, G.; Smith, C.; Sigmon, B.; Herr, J.R.; Alfano, J.R.; Ge, Y.; Schnable, J.C.; et al. Association analyses of host genetics, root-colonizing microbes, and plant phenotypes under different nitrogen conditions in maize. eLife 2022, 11, e75790. [Google Scholar] [CrossRef] [PubMed]
  56. Porter, A.S.; Evans-Fitz Gerald, C.; McElwain, J.C.; Yiotis, C.; Elliott-Kingston, C. How well do you know your growth chambers? Testing for chamber effect using plant traits. Plant Methods 2015, 11, 44. [Google Scholar] [CrossRef] [PubMed]
  57. Gruber-Dorninger, C.; Pester, M.; Kitzinger, K.; Savio, D.F.; Loy, A.; Rattei, T.; Wagner, M.; Daims, H. Functionally relevant diversity of closely related Nitrospira in activated sludge. ISME J. 2015, 9, 643–655. [Google Scholar] [CrossRef] [PubMed]
  58. Lacombe, S.; Bradley, R.L.; Hamel, C.; Beaulieu, C. Do tree-based intercropping systems increase the diversity and stability of soil microbial communities? Agric. Ecosyst. Environ. 2009, 131, 25–31. [Google Scholar] [CrossRef]
  59. Nottingham, A.T.; Scott, J.J.; Saltonstall, K.; Broders, K.; Montero-Sanchez, M.; Püspök, J.; Bååth, E.; Meir, P. Microbial diversity declines in warmed tropical soil and respiration rise exceed predictions as communities adapt. Nat. Microbiol. 2022, 7, 1650–1660. [Google Scholar] [CrossRef]
  60. Zegeye, E.K.; Brislawn, C.J.; Farris, Y.; Fansler, S.J.; Hofmockel, K.S.; Jansson, J.K.; Wright, A.T.; Graham, E.B.; Naylor, D.; McClure, R.S.; et al. Selection, Succession, and Stabilization of Soil Microbial Consortia. mSystems 2019, 4, e00055-19. [Google Scholar] [CrossRef]
  61. Lebeis, S.L.; Paredes, S.H.; Lundberg, D.S.; Breakfield, N.; Gehring, J.; McDonald, M.; Malfatti, S.; del Rio, T.G.; Jones, C.D.; Tringe, S.G.; et al. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science 2015, 349, 860–864. [Google Scholar] [CrossRef]
  62. Brown, S.P.; Grillo, M.A.; Podowski, J.C.; Heath, K.D. Soil origin and plant genotype structure distinct microbiome compartments in the model legume Medicago truncatula. Microbiome 2020, 8, 139. [Google Scholar] [CrossRef] [PubMed]
  63. Basu, A.; Prasad, P.; Das, S.N.; Kalam, S.; Sayyed, R.Z.; Reddy, M.S.; El Enshasy, H. Plant Growth Promoting Rhizobacteria (PGPR) as Green Bioinoculants: Recent Developments, Constraints, and Prospects. Sustainability 2021, 13, 1140. [Google Scholar] [CrossRef]
  64. Luo, J.; Guo, X.; Tao, Q.; Li, J.; Liu, Y.; Du, Y.; Liu, Y.; Liang, Y.; Li, T. Succession of the composition and co-occurrence networks of rhizosphere microbiota is linked to Cd/Zn hyperaccumulation. Soil Biol. Biochem. 2021, 153, 108120. [Google Scholar] [CrossRef]
  65. Krishna, M.; Gupta, S.; Delgado-Baquerizo, M.; Morriën, E.; Garkoti, S.C.; Chaturvedi, R.; Ahmad, S. Successional trajectory of bacterial communities in soil are shaped by plant-driven changes during secondary succession. Sci. Rep. 2020, 10, 9864. [Google Scholar] [CrossRef] [PubMed]
  66. Micallef, S.; Colόn-Carmona, A. Genetic and Developmental Control of Rhizosphere Bacterial Communities. In Molecular Microbial Ecology of the Rhizosphere, 1st ed.; De Bruijn, F.J., Ed.; Wiley: New York, NY, USA, 2013; pp. 257–263. Available online: https://onlinelibrary.wiley.com/doi/10.1002/9781118297674.ch24 (accessed on 19 November 2023).
  67. Micallef, S.A.; Channer, S.; Shiaris, M.P.; Colón-Carmona, A. Plant age and genotype impact the progression of bacterial community succession in the Arabidopsis rhizosphere. Plant Signal. Behav. 2009, 4, 777–780. [Google Scholar] [CrossRef] [PubMed]
  68. Monchgesang, S.; Strehmel, N.; Schmidt, S.; Westphal, L.; Taruttis, F.; Muller, E.; Herklotz, S.; Neumann, S.; Scheel, D. Natural variation of root exudates in Arabidopsis thaliana-linking metabolomic and genomic data. Sci. Rep. 2016, 6, 29033. [Google Scholar] [CrossRef] [PubMed]
  69. Matsushima, C.; Shenton, M.; Kitahara, A.; Wasaki, J.; Oikawa, A.; Cheng, W.; Ikeo, K.; Tawaraya, K. Multiple analysis of root exudates and microbiome in rice (Oryza sativa) under low P conditions. Arch. Microbiol. 2021, 203, 5599–5611. [Google Scholar]
  70. Huang, X.F.; Chaparro, J.M.; Reardon, K.F.; Zhang, R.; Shen, Q.; Vivanco, J.M. Rhizosphere interactions: Root exudates, microbes, and microbial communities. Botany 2014, 92, 267–275. [Google Scholar] [CrossRef]
  71. Sun, H.; Jiang, S.; Jiang, C.; Wu, C.; Gao, M.; Wang, Q. A review of root exudates and rhizosphere microbiome for crop production. Env. Sci. Pollut. Res. 2021, 28, 54497–54510. [Google Scholar] [CrossRef]
  72. Li, Z.; Guo, W.; Mo, C.; Tang, R.; He, L.; Du, L.; Li, M.; Wu, H.; Tang, X.; Huang, Z.; et al. Root Metabolism and Effects of Root Exudates on the Growth of Ralstonia solanacearum and Fusarium moniliforme Were Significantly Different between the Two Genotypes of Peanuts. Genes 2023, 14, 528. [Google Scholar] [CrossRef]
  73. Badri, D.V.; Vivanco, J.M. Regulation and function of root exudates. Plant Cell Environ. 2009, 32, 666–681. [Google Scholar] [CrossRef]
  74. Cotton, T.E.A.; Pétriacq, P.; Cameron, D.D.; Meselmani, M.A.; Schwarzenbacher, R.; Rolfe, S.A.; Ton, J. Metabolic regulation of the maize rhizobiome by benzoxazinoids. ISME J. 2019, 13, 1647–1658. [Google Scholar] [CrossRef]
  75. Cadot, S.; Guan, H.; Bigalke, M.; Walser, J.C.; Jander, G.; Erb, M.; van der Heijden, M.G.; Schlaeppi, K. Specific and conserved patterns of microbiota-structuring by maize benzoxazinoids in the field. Microbiome 2021, 9, 103. [Google Scholar] [CrossRef]
  76. Korenblum, E.; Dong, Y.; Szymanski, J.; Panda, S.; Jozwiak, A.; Massalha, H.; Meir, S.; Rogachev, I.; Aharoni, A. Rhizosphere microbiome mediates systemic root metabolite exudation by root-to-root signaling. Proc. Natl. Acad. Sci. USA 2020, 117, 3874–3883. [Google Scholar] [CrossRef] [PubMed]
  77. Rudrappa, T.; Czymmek, K.J.; Paré, P.W.; Bais, H.P. Root-Secreted Malic Acid Recruits Beneficial Soil Bacteria. Plant Physiol. 2008, 148, 1547–1556. [Google Scholar] [CrossRef] [PubMed]
  78. Trivedi, P.; Leach, J.E.; Tringe, S.G.; Sa, T.; Singh, B.K. Plant–microbiome interactions: From community assembly to plant health. Nat. Rev. Microbiol. 2020, 18, 607–621. [Google Scholar] [CrossRef]
  79. Gu, S.; Wei, Z.; Shao, Z.; Friman, V.P.; Cao, K.; Yang, T.; Kramer, J.; Wang, X.; Li, M.; Mei, X.; et al. Competition for iron drives phytopathogen control by natural rhizosphere microbiomes. Nat. Microbiol. 2020, 5, 1002–1010. [Google Scholar] [CrossRef] [PubMed]
  80. Bai, Y.C.; Chang, Y.Y.; Hussain, M.; Lu, B.; Zhang, J.P.; Song, X.B.; Lei, X.S.; Pei, D. Soil Chemical and Microbiological Properties Are Changed by Long-Term Chemical Fertilizers That Limit Ecosystem Functioning. Microorganisms 2020, 8, 694. [Google Scholar] [CrossRef]
  81. Xu, Q.; Ling, N.; Chen, H.; Duan, Y.; Wang, S.; Shen, Q.; Vandenkoornhuyse, P. Long-Term Chemical-Only Fertilization Induces a Diversity Decline and Deep Selection on the Soil Bacteria. mSystems 2020, 5, e00337-20. [Google Scholar] [CrossRef]
  82. Cardinale, M.; Ratering, S.; Suarez, C.; Zapata Montoya, A.M.; Geissler-Plaum, R.; Schnell, S. Paradox of plant growth promotion potential of rhizobacteria and their actual promotion effect on growth of barley (Hordeum vulgare L.) under salt stress. Microbiol. Res. 2015, 181, 22–32. [Google Scholar] [CrossRef]
  83. Meena, M.; Swapnil, P.; Divyanshu, K.; Kumar, S.; Harish; Tripathi, Y.N.; Zehra, A.; Marwal, A.; Upadhyay, R.S. PGPR-mediated induction of systemic resistance and physiochemical alterations in plants against the pathogens: Current perspectives. J. Basic Microbiol. 2020, 60, 828–861. [Google Scholar] [CrossRef] [PubMed]
  84. Kamilova, F.; Okon, Y.; de Weert, S.; Hora, K. Commercialization of Microbes: Manufacturing, Inoculation, Best Practice for Objective Field Testing, and Registration. In Principles of Plant-Microbe Interactions: Microbes for Sustainable Agriculture; Lugtenberg, B., Ed.; Springer International Publishing: Cham, Switzerland, 2015; pp. 319–327. [Google Scholar] [CrossRef]
  85. Dubey, S.; Sharma, S. Rhizospheric Engineering by Plant-Mediated Indirect Selection of Microbiome for Agricultural Sustainability. Crit. Rev. Plant Sci. 2021, 40, 379–397. [Google Scholar] [CrossRef]
  86. Wright, R.J.; Gibson, M.I.; Christie-Oleza, J.A. Understanding microbial community dynamics to improve optimal microbiome selection. Microbiome 2019, 7, 85. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Ecosystem selection leads to changes in above-ground plant biomass. Above-ground biomass for genotypes Ler (a) and Cvi (b) of Arabidopsis thaliana. Dark lines represent the high-biomass treatment, and light lines represent the low-biomass treatment. Each point represents the mean of n = 50 microcosms. The mean biomass of the selection lines differed for several generations of both genotypes (t-test p-value * < 0.05; *** < 0.001).
Figure 1. Ecosystem selection leads to changes in above-ground plant biomass. Above-ground biomass for genotypes Ler (a) and Cvi (b) of Arabidopsis thaliana. Dark lines represent the high-biomass treatment, and light lines represent the low-biomass treatment. Each point represents the mean of n = 50 microcosms. The mean biomass of the selection lines differed for several generations of both genotypes (t-test p-value * < 0.05; *** < 0.001).
Microorganisms 11 02932 g001
Figure 2. Above-ground biomass represented as deviation from the mean for each generation. By plotting deviations from the mean, the influence of REVs on above-ground biomass is mitigated for genotypes Ler (a) and Cvi (b). The dark lines denote high-biomass treatments, while the light lines represent low-biomass treatments. Each data point represents the average of n = 50 microcosms (t-test p-value * < 0.05; **** < 0.001).
Figure 2. Above-ground biomass represented as deviation from the mean for each generation. By plotting deviations from the mean, the influence of REVs on above-ground biomass is mitigated for genotypes Ler (a) and Cvi (b). The dark lines denote high-biomass treatments, while the light lines represent low-biomass treatments. Each data point represents the average of n = 50 microcosms (t-test p-value * < 0.05; **** < 0.001).
Microorganisms 11 02932 g002
Figure 3. Taxonomic distribution of the microbial community is represented in terms of the relative abundance at the Class level. The plot shows eight generations (1–8) of the high- and low-biomass treatments of both the Cvi (left) and Ler (right) genotypes of Arabidopsis thaliana. Ten samples were assayed for the generation 1 experiments, but only five for all subsequent generations. The most represented classes across all samples were Alphaproteobacteria (24.7%), Gammaproteobacteria (11.4%), Bacteroidia (10.1%), and Verrucomicrobiae (6.2%).
Figure 3. Taxonomic distribution of the microbial community is represented in terms of the relative abundance at the Class level. The plot shows eight generations (1–8) of the high- and low-biomass treatments of both the Cvi (left) and Ler (right) genotypes of Arabidopsis thaliana. Ten samples were assayed for the generation 1 experiments, but only five for all subsequent generations. The most represented classes across all samples were Alphaproteobacteria (24.7%), Gammaproteobacteria (11.4%), Bacteroidia (10.1%), and Verrucomicrobiae (6.2%).
Microorganisms 11 02932 g003
Figure 4. Changes in alpha diversity. Two alpha diversity indices were computed, including the Observed and Shannon metrics. Both metrics demonstrated a rapid decline in diversity between generations 1 and 5, with a less pronounced decrease after generation 5. Additionally, a negative exponential function was fit to the Observed (R2 = 0.81, p < 0.01) and Shannon (R2 = 0.58, p < 0.01) metrics, with Observed showing a better fit, indicating the increased stability of the microbiome in later generations (Wilcox test, *** < 0.001).
Figure 4. Changes in alpha diversity. Two alpha diversity indices were computed, including the Observed and Shannon metrics. Both metrics demonstrated a rapid decline in diversity between generations 1 and 5, with a less pronounced decrease after generation 5. Additionally, a negative exponential function was fit to the Observed (R2 = 0.81, p < 0.01) and Shannon (R2 = 0.58, p < 0.01) metrics, with Observed showing a better fit, indicating the increased stability of the microbiome in later generations (Wilcox test, *** < 0.001).
Microorganisms 11 02932 g004
Figure 5. Host genotype and biomass treatment shape the composition of the microbiome. Principle coordinate analysis (PCoA) was conducted using unweighted UniFrac for generations 1 through 8. Each point represents an individual sample, with microbial communities defined as Cvi, represented by red/orange circles, and Ler, represented by dark/light blue triangles. The first two coordinate axes to be plotted accounted for the highest variation in the data, ranging from 13% in generation 1 to approximately 52% in generation 8. While there was no discernible clustering in generation 1, both genotypes and high/low biomass treatments exhibited greater clustering over the course of the experiment.
Figure 5. Host genotype and biomass treatment shape the composition of the microbiome. Principle coordinate analysis (PCoA) was conducted using unweighted UniFrac for generations 1 through 8. Each point represents an individual sample, with microbial communities defined as Cvi, represented by red/orange circles, and Ler, represented by dark/light blue triangles. The first two coordinate axes to be plotted accounted for the highest variation in the data, ranging from 13% in generation 1 to approximately 52% in generation 8. While there was no discernible clustering in generation 1, both genotypes and high/low biomass treatments exhibited greater clustering over the course of the experiment.
Microorganisms 11 02932 g005
Figure 6. Variance in microbiome structure and composition, explained by genotype and biomass treatment increases over generations. The R2 values were computed using the PERMANOVA test with the adonis2 function for generations 1 through 8 for (a) weighted and (b) unweighted UniFrac distances. The resulting R2 and variance explained by residual values for the biomass treatment (brown), genotype (blue), and residuals (black) are plotted. Residuals are plotted as a proxy for random environmental variations (REVs). There was a marked increase in the influence of both genotype and biomass treatment on differences in the microbial community from generations 1 to 8. Despite observing a decrease in the variance accounted for by the residuals, which could potentially signify a decline in the influence of REVs, the model elucidated more than 50 percent of the variability within the microbial community in the 8th generation. These differences are more pronounced in the weighted vs. the unweighted UniFrac model (distance.matrix ~ biomass treatment + genotype).
Figure 6. Variance in microbiome structure and composition, explained by genotype and biomass treatment increases over generations. The R2 values were computed using the PERMANOVA test with the adonis2 function for generations 1 through 8 for (a) weighted and (b) unweighted UniFrac distances. The resulting R2 and variance explained by residual values for the biomass treatment (brown), genotype (blue), and residuals (black) are plotted. Residuals are plotted as a proxy for random environmental variations (REVs). There was a marked increase in the influence of both genotype and biomass treatment on differences in the microbial community from generations 1 to 8. Despite observing a decrease in the variance accounted for by the residuals, which could potentially signify a decline in the influence of REVs, the model elucidated more than 50 percent of the variability within the microbial community in the 8th generation. These differences are more pronounced in the weighted vs. the unweighted UniFrac model (distance.matrix ~ biomass treatment + genotype).
Microorganisms 11 02932 g006
Figure 7. Fit of the neutral model diminished as generations progressed. The occupancy (prevalence of OTU in samples), plotted against log10 abundance, is depicted using the Burns model. OTUs that are neutral, i.e., not selected for or against, are indicated in black. Green- and blue-colored OTUs signify positive and negative selection, respectively. This plot highlights the presence of numerous microbes that underwent directed selection throughout the experiment.
Figure 7. Fit of the neutral model diminished as generations progressed. The occupancy (prevalence of OTU in samples), plotted against log10 abundance, is depicted using the Burns model. OTUs that are neutral, i.e., not selected for or against, are indicated in black. Green- and blue-colored OTUs signify positive and negative selection, respectively. This plot highlights the presence of numerous microbes that underwent directed selection throughout the experiment.
Microorganisms 11 02932 g007
Figure 8. Differentially abundant OTUs. ANCOMBC2 results are illustrated as log2 fold changes at the family level for high-biomass vs. low-biomass samples. Each taxon listed in the figure represents a single OTU in either the Ler (a) or Cvi (b) microbial community, identified to the family level if possible. The scale indicates log2 fold changes, with red and blue representing positive and negative fold changes, respectively, in the high-biomass samples. Well-known plant-growth-promoting bacteria, such as Paenibacillaceae, Bacilli, Labraceae, Rhizobiaceae, andBdellovibrio were present in higher abundance. A vast majority of the other taxa identified were uncultivated, and there is little to no information on them. (c,d): Box plots of relative abundance for all OTUs that were enriched in high-biomass treatment samples of genotypes Ler and Cvi, respectively, with colors indicating different generations. These OTUs demonstrated a consistent trend of increasing relative abundance across successive generations. To test this trend (red line), a linear model was fitted to the data, resulting in an R2 value of 0.57 (p < 0.01) for Ler and 0.73 (p < 0.01) for Cvi. (e,f) are box plots of relative abundance for all OTUs that were depleted in samples from the high-biomass treatment in genotypes Ler and Cvi, respectively, and are colored according to generation.
Figure 8. Differentially abundant OTUs. ANCOMBC2 results are illustrated as log2 fold changes at the family level for high-biomass vs. low-biomass samples. Each taxon listed in the figure represents a single OTU in either the Ler (a) or Cvi (b) microbial community, identified to the family level if possible. The scale indicates log2 fold changes, with red and blue representing positive and negative fold changes, respectively, in the high-biomass samples. Well-known plant-growth-promoting bacteria, such as Paenibacillaceae, Bacilli, Labraceae, Rhizobiaceae, andBdellovibrio were present in higher abundance. A vast majority of the other taxa identified were uncultivated, and there is little to no information on them. (c,d): Box plots of relative abundance for all OTUs that were enriched in high-biomass treatment samples of genotypes Ler and Cvi, respectively, with colors indicating different generations. These OTUs demonstrated a consistent trend of increasing relative abundance across successive generations. To test this trend (red line), a linear model was fitted to the data, resulting in an R2 value of 0.57 (p < 0.01) for Ler and 0.73 (p < 0.01) for Cvi. (e,f) are box plots of relative abundance for all OTUs that were depleted in samples from the high-biomass treatment in genotypes Ler and Cvi, respectively, and are colored according to generation.
Microorganisms 11 02932 g008
Figure 9. Pairwise spearman correlation analysis for differentially abundant OTUs. Spearman correlation results are illustrated between OTUs enriched in the high-biomass-treatment samples of the Ler genotype in the early phase (generation 1 and generation 2) (a) and late phase (generation 6, generation 7, and generation 8) (b) of the experiment. The color gradient depicts the correlation coefficient’s strength for each OTU pair, from −1 (dark red) to 1 (dark blue). The legend represents a key, where each number corresponds to the OTU name and taxonomy at the family level, when available. (c): Barplot of a subset of OTUs. The bars represent OTUs, along with the number of pairwise positive or negative associations. This visualization highlights trends over the early and late phases of the experiment, with some OTUs exhibiting increases in positive pairwise associations (OTU1606_Chitinophagaceae, OTU162_Bdellovibrionaceae, OTU94_Class_Bacilli, and OTU359_Order_AKYH767) as generations progressed and decreases in negative pairwise associations, while others exhibited the opposite trend (OTU763_Fimbriimonadaceae, OTU142_Sphingobacteriaceae, and OTU1005_Labraceae).
Figure 9. Pairwise spearman correlation analysis for differentially abundant OTUs. Spearman correlation results are illustrated between OTUs enriched in the high-biomass-treatment samples of the Ler genotype in the early phase (generation 1 and generation 2) (a) and late phase (generation 6, generation 7, and generation 8) (b) of the experiment. The color gradient depicts the correlation coefficient’s strength for each OTU pair, from −1 (dark red) to 1 (dark blue). The legend represents a key, where each number corresponds to the OTU name and taxonomy at the family level, when available. (c): Barplot of a subset of OTUs. The bars represent OTUs, along with the number of pairwise positive or negative associations. This visualization highlights trends over the early and late phases of the experiment, with some OTUs exhibiting increases in positive pairwise associations (OTU1606_Chitinophagaceae, OTU162_Bdellovibrionaceae, OTU94_Class_Bacilli, and OTU359_Order_AKYH767) as generations progressed and decreases in negative pairwise associations, while others exhibited the opposite trend (OTU763_Fimbriimonadaceae, OTU142_Sphingobacteriaceae, and OTU1005_Labraceae).
Microorganisms 11 02932 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shankar, N.; Shetty, P.; Melo, T.C.; Kesseli, R. Multi-Generation Ecosystem Selection of Rhizosphere Microbial Communities Associated with Plant Genotype and Biomass in Arabidopsis thaliana. Microorganisms 2023, 11, 2932. https://doi.org/10.3390/microorganisms11122932

AMA Style

Shankar N, Shetty P, Melo TC, Kesseli R. Multi-Generation Ecosystem Selection of Rhizosphere Microbial Communities Associated with Plant Genotype and Biomass in Arabidopsis thaliana. Microorganisms. 2023; 11(12):2932. https://doi.org/10.3390/microorganisms11122932

Chicago/Turabian Style

Shankar, Nachiket, Prateek Shetty, Tatiana C. Melo, and Rick Kesseli. 2023. "Multi-Generation Ecosystem Selection of Rhizosphere Microbial Communities Associated with Plant Genotype and Biomass in Arabidopsis thaliana" Microorganisms 11, no. 12: 2932. https://doi.org/10.3390/microorganisms11122932

APA Style

Shankar, N., Shetty, P., Melo, T. C., & Kesseli, R. (2023). Multi-Generation Ecosystem Selection of Rhizosphere Microbial Communities Associated with Plant Genotype and Biomass in Arabidopsis thaliana. Microorganisms, 11(12), 2932. https://doi.org/10.3390/microorganisms11122932

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