Research Article | Ecological and Evolutionary Science

Long-Term Nutrient Enrichment of an Oligotroph-Dominated Wetland Increases Bacterial Diversity in Bulk Soils and Plant Rhizospheres

Regina B. Bledsoe, Carol Goodwillie, Ariane L. Peralta
Barbara J. Campbell, Editor
DOI: 10.1128/mSphere.00035-20

ABSTRACT

In nutrient-limited conditions, plants rely on rhizosphere microbial members to facilitate nutrient acquisition, and in return, plants provide carbon resources to these root-associated microorganisms. However, atmospheric nutrient deposition can affect plant-microbe relationships by changing soil bacterial composition and by reducing cooperation between microbial taxa and plants. To examine how long-term nutrient addition shapes rhizosphere community composition, we compared traits associated with bacterial (fast-growing copiotrophs, slow-growing oligotrophs) and plant (C3 forb, C4 grass) communities residing in a nutrient-poor wetland ecosystem. Results revealed that oligotrophic taxa dominated soil bacterial communities and that fertilization increased the presence of oligotrophs in bulk and rhizosphere communities. Additionally, bacterial species diversity was greatest in fertilized soils, particularly in bulk soils. Nutrient enrichment (fertilized versus unfertilized) and plant association (bulk versus rhizosphere) determined bacterial community composition; bacterial community structure associated with plant functional group (grass versus forb) was similar within treatments but differed between fertilization treatments. The core forb microbiome consisted of 602 unique taxa, and the core grass microbiome consisted of 372 unique taxa. Forb rhizospheres were enriched in potentially disease-suppressive bacterial taxa, and grass rhizospheres were enriched in bacterial taxa associated with complex carbon decomposition. Results from this study demonstrate that fertilization serves as a strong environmental filter on the soil microbiome, which leads to distinct rhizosphere communities and can shift plant effects on the rhizosphere microbiome. These taxonomic shifts within plant rhizospheres could have implications for plant health and ecosystem functions associated with carbon and nitrogen cycling.

IMPORTANCE Over the last century, humans have substantially altered nitrogen and phosphorus cycling. Use of synthetic fertilizer and burning of fossil fuels and biomass have increased nitrogen and phosphorus deposition, which results in unintended fertilization of historically low-nutrient ecosystems. With increased nutrient availability, plant biodiversity is expected to decline, and the abundance of copiotrophic taxa is anticipated to increase in bacterial communities. Here, we address how bacterial communities associated with different plant functional types (forb, grass) shift due to long-term nutrient enrichment. Unlike other studies, results revealed an increase in bacterial diversity, particularly of oligotrophic bacteria in fertilized plots. We observed that nutrient addition strongly determines forb and grass rhizosphere composition, which could indicate different metabolic preferences in the bacterial communities. This study highlights how long-term fertilization of oligotroph-dominated wetlands could alter diversity and metabolism of rhizosphere bacterial communities in unexpected ways.

INTRODUCTION

The soil microbiome is critical for plant health, fitness, and diversity, especially in nutrient-limited environments (1–4). In particular, within the rhizosphere, plants provide carbon (C) resources to soil microorganisms in exchange for nutrients such as nitrogen (N) and phosphorus (P). However, nutrient enrichment has been documented to disrupt plant-microbe mutualisms (2). Over the last century, agricultural fertilization and the burning of fossil fuels and biomass have indirectly led to nutrient deposition onto historically low-nutrient ecosystems (5–8). Nutrient enrichment generally causes reduced plant species diversity (9, 10) sometimes as a shift in plant functional types with an increase in grass biomass and loss of forb diversity (11–13). Fertilization has also been shown to decrease soil microbial diversity across cropland, grassland, forest, and tundra ecosystems (14–16). Despite patterns that have emerged from these bulk soil studies, it is less clear how changes in soil microbial diversity due to nutrient additions influence rhizosphere microbial community assembly and diversity. We address this knowledge gap by comparing changes in rhizosphere bacterial community composition of a grass and forb within a long-term fertilization experiment.

Both bulk soil matrix (i.e., not in contact with plant roots) properties and plant identity influence rhizosphere microbial communities. The bulk soil matrix is the reservoir of microbial diversity from which rhizosphere-associated microbial communities are selected; therefore, shifts in bulk soil microbial communities affect rhizosphere assemblages (17–19). In many cases, N, N and P, and N-P-K (nitrogen-phosphorus-potassium) fertilization decreases soil bacterial diversity (14–16). Additionally, nutrient enrichment selects for more copiotrophic (i.e., fast-growing, r-strategists) microbial heterotrophs that preferentially metabolize labile C sources versus oligotrophic (i.e., slow-growing, K-strategist) microbial species, which can metabolize complex C sources (20–23). A molecular marker to identify life history strategy (i.e., copiotroph or oligotroph) is rRNA (rrn) gene copy number (23–26). Bacterial taxa are estimated to contain 1 to 15 rRNA gene copies, with faster-growing taxa containing higher numbers of gene copies than slower-growing taxa (20, 23–27). Specifically, bacterial growth rate is limited by the transcription rate of rRNA, such that growth rate is estimated to double with doubling of rRNA gene copy number (23). Further, several studies indicate that fertilization increases the abundance of copiotrophic bacterial groups within Actinobacteria, Alphaproteobacteria, and Gammaproteobacteria and decreases abundance in oligotrophic bacterial groups within Acidobacteria, Nitrospirae, Planctomycetes, and Deltaproteobacteria of bulk soils (15, 21, 28, 29). Additionally, copiotrophic taxa within Alpha-, Beta-, and Gammaproteobacteria, Actinobacteria, Firmicutes, and Bacteroidetes are dominant members of some rhizosphere communities (17, 30, 31).

While the bulk soil environment is the primary source of rhizosphere diversity, plant species also influence rhizosphere bacterial community assembly due to variation in rhizodeposition (30–33). Rhizodeposits include nutrients, exudates, root cells, and mucilage released by plant roots (34). Plants allocate 5 to 20% of photosynthetically fixed C belowground (35–37). Some estimates suggest that up to 40% of fixed C is translocated belowground (38), and grasses are suggested to be near that upper limit with ∼30% of fixed C allocated belowground (39). These rhizodeposits also include root exudates which are composed of sugars, organic acids, phenolic compounds, and amino acids (1, 17, 40, 41). Differences in plant physiology influencing the quantity and composition of root exudates can affect rhizosphere bacterial community composition. For example, C4 grasses have higher photosynthetic rates (i.e., fix more C) and greater root biomass allocation compared to C3 plants, resulting in a greater quantity of root exudates (42, 43). C3 plant root exudates can contain a greater variety of organic acids and amino acids along with the sugars mannose, maltose, and ribose compared to C4 plant root exudates, which can contain several sugar alcohols (i.e., inositol, erythritol, and ribitol) (44). However, N fertilization has been shown to increase C assimilation in plants but decrease belowground allocation of assimilated C while increasing total C into soils as rhizodeposits (39, 45). Prior studies revealed that root exudation of organic C can be higher in both low-nutrient scenarios (46, 47) and high-nutrient scenarios (48, 49). Further, differences in soil nutrient status can change the composition (i.e., carbohydrates, organic acids, and amino acid concentrations) of root exudates (46, 50). Thus, fertilization and plant-specific rhizodeposition patterns of C3 forbs and C4 grasses are predicted to differentially affect rhizosphere bacterial community structure.

In this study, we address the following question: to what extent does long-term fertilization (N-P-K) of bulk soil shift rhizosphere bacterial communities of two plant species representing distinct functional types (i.e., a C3 forb and a C4 grass)? First, we hypothesize that nutrient addition will decrease bacterial species diversity and increase the abundance of copiotrophic taxa in all soils, especially rhizosphere soils due to increased availability of labile C from root exudates. We expect that fertilization will stimulate microbial activity of faster-growing copiotrophic species, which would outcompete slower-growing oligotrophic species and result in decreased bacterial diversity. This effect is predicted to be amplified within plant rhizospheres due to the availability of labile C substrates in root exudates, which should preferentially select for copiotrophic bacteria. Second, we hypothesize that fertilization will be the primary factor determining differences in rhizosphere communities and plant identity will secondarily influence the rhizosphere community. If bulk soil is the reservoir for the rhizosphere community, then fertilization will more strongly determine rhizosphere bacterial diversity and community composition. In addition, plant type can also affect rhizosphere communities due to differences in root exudate composition; however, fertilization effects will constrain rhizosphere effects. As a result, plant species are expected to associate with unique core microbiomes that differ between fertilization treatments.

To test these hypotheses, bulk and rhizosphere soils were sampled from two plant species (a grass and a forb) from fertilized and unfertilized plots at a long-term disturbance and fertilization experiment (established in 2003). Bacterial communities were identified using 16S rRNA amplicon sequencing which allowed binning of bacterial taxa as copiotrophic or oligotrophic by estimating the average rRNA (rrn) gene copy number. By evaluating differences in taxonomic information and 16S rRNA gene copy numbers of bulk and rhizosphere soils of two plant species with associated soil properties (i.e., ammonium, nitrate, soil pH, carbon, and moisture), we provide insight into biotic and abiotic processes that are contributing to rhizosphere bacterial community assembly.

RESULTS

Soil source and fertilization distinguish soil properties.The main effect of fertilization was significantly different in the soil physiochemical property of pH (P = 0.02), and the main effect of soil source (bulk versus rhizosphere) was significantly different in the soil physiochemical properties of pH (P < 0.001), nitrate (P < 0.0001), C percent (P = 0.03), and N percent (P = 0.04; Table 1; see also Table S1 in the supplemental material). Rhizosphere soils were more similar to each other in soil properties than to bulk soils (Table 1, Table S1, Tukey’s honestly significant difference [HSD] test, P < 0.05). Specifically, bulk soil had lower total C and N and nitrate concentrations than forb rhizospheres with grass rhizospheres having the highest values (Table 1, Table S1, Tukey’s HSD test, P < 0.05). Soil pH was lowest in rhizosphere soils compared to bulk soils but higher in fertilized soils compared to unfertilized soils within soil sources (Table 1, Table S1, Tukey HSD, P < 0.05).

View this table:
TABLE 1

Soil physiochemical properties after 12 years of fertilization and mowing disturbancea

TABLE S1

Summary of two-way ANOVA comparing soil properties among soil source (bulk, grass rhizosphere, forb rhizosphere) and fertilization treatments. Download Table S1, DOCX file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Fertilization increased soil bacterial diversity in bulk and rhizosphere soils.Chao1 bacterial richness (P < 0.0001) and Shannon H′ diversity (P < 0.0001) were higher in fertilized soils than in unfertilized soils (Fig. 1A, Table S2). In addition, the main effect of soil source influenced bacterial diversity; bulk soil bacterial diversity was significantly higher than rhizosphere soil diversity (Tukey’s HSD test, P < 0.05, Fig. 1B, Table S2). Finally, results revealed a positive relationship between Shannon H′ diversity and pH, where pH explained 71% (P = 0.0003) and 32% (P = 0.03) of the variation in bacterial diversity in unfertilized and fertilized treatments, respectively, across all soil sources (Fig. 1C).

FIG 1
FIG 1

Bacterial diversity patterns according to soil source, fertilization, and soil pH. Boxplots of bacterial diversity for Chao1 richness (A) and Shannon diversity Hʹ (B) associated with soil source (bulk, grass rhizosphere, forb rhizosphere) and fertilization treatment. (C) Linear regression of soil pH and bacterial community Shannon diversity H′ by fertilization treatment with 95% confidence intervals (fertilized, R2 = 0.32, P = 0.03; unfertilized: R2 = 0.71, P = 0.0003). Colors indicate fertilization treatment. Asterisks indicate significant differences between fertilization treatments, and lowercase letters represent significant differences between soil sources (Tukey’s HSD test, P < 0.05).

TABLE S2

Summary of two-way ANOVA comparing bacterial community Chao1 richness (A) and Shannon diversity H′ (B) metrics among soil source and fertilization treatments. Source represents bulk, grass rhizosphere, and forb rhizosphere, and treatment represents fertilized and unfertilized mowed treatments. Main effects that were significantly different (ANOVA P < 0.05) are bold. Download Table S2, DOCX file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Copiotroph-to-oligotroph ratios indicated oligotroph-dominated bacterial communities.Across all samples, we detected 9 to 30 copiotrophic taxa and 82 to 190 oligotrophic taxa at the class level. This resulted in copiotroph-to-oligotroph ratios of <0.2 within all treatment combinations. Nutrient additions significantly decreased the ratio of copiotrophs to oligotrophs in bulk soils compared to rhizosphere soils (Tukey’s HSD test, P < 0.05; Table S4; Fig. 2). Finally, there was no relationship between bacterial Shannon H′ diversity and copiotroph-to-oligotroph ratio (fertilized, R2 = −0.01, P = 0.38; unfertilized, R2 = 0.14, P = 0.13) (see Fig. S1 in the supplemental material).

FIG 2
FIG 2

Comparison of bacterial life history traits. Boxplots of copiotroph-to-oligotroph ratios (based on 16S rRNA sequences) according to soil source (bulk, grass rhizosphere, and forb rhizosphere) and fertilization treatment. Boxplots are colored according to fertilization treatment. Lowercase letters indicate significant differences among soil sources (Tukey’s HSD test, P < 0.05).

FIG S1

Linear regression of copiotroph-to-oligotroph ratio and Shannon diversity H′ by fertilization treatment. Gray confidence bands represent 95% confidence intervals. Fertilized, R2 = −0.01, P = 0.38; unfertilized, R2 = 0.14, P = 0.13. Download FIG S1, PDF file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Fertilization treatment and soil source influenced bacterial community composition.Specifically, fertilization treatment (along principal-coordinate analysis [PCoA] axis 1) explained 31.6% of variation in bacterial community composition, while soil source (primarily bulk versus rhizosphere) separated bacterial composition (along PCoA axis 2) and explained 22.5% of bacterial community variation (Fig. 3). The main effects of soil source (permutational multivariate analysis of variance [PERMANOVA], R2 = 0.23, P = 0.001) and fertilization treatment (PERMANOVA, R2 = 0.281, P = 0.001) influenced bacterial community composition (Table S3A). According to pairwise comparisons, rhizosphere bacterial community composition was similar between grass and forb rhizosphere samples within fertilization treatments (Table S3B). When examining relationships between community composition and soil characteristics, higher soil pH and moisture were correlated with fertilized bulk soils (Fig. 3). Further, higher concentrations of soil C and N were correlated with rhizosphere community composition (Fig. 3).

FIG 3
FIG 3

Ordination based on principal-coordinate analysis depicting bacterial community composition. Colors represent fertilization treatment, and symbols represent soil source (bulk soil, grass rhizosphere, and forb rhizosphere). Vectors represent soil factors that are correlated with patterns in bacterial community composition (P < 0.05) (soil pH [pH], soil gravimetric moisture percentage [moisture], total soil carbon [C%], total soil nitrogen [N%]).

TABLE S3

Summary of PERMANOVA main effects (soil source and fertilization treatment) and interaction (A) and pairwise PERMANOVA comparisons of soil sources (bulk, grass rhizosphere, forb rhizosphere) within fertilization treatments (B). Download Table S3, DOCX file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S4

Summary of two-way ANOVA comparing bacterial community copiotroph-to-oligotroph ratio among soil source and fertilization treatments. Download Table S4, DOCX file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Different bacterial taxa (OTUs) represented fertilization treatments and plant species.We compared bacterial community taxonomic shifts in unfertilized and fertilized bulk soils and then grass and forb rhizospheres, concluding with differences in microbiome structure between the two plant species. Within bulk soil samples, important indicator species for bacterial communities within unfertilized plots were from the class Alphaproteobacteria with 1 operational taxonomic unit (OTU) from the order Rhizobiales and 2 OTUs from Rhodospirillales and 3 OTUs from the class Spartobacteria (Table S5). In contrast, fertilized bulk soils were best represented by members of the class Actinobacteria with 1 OTU from the order Actinomycetales and 2 OTUs from the order Solirubrobacterales. While OTUs within Rhizobiales were identified as indicator species for bacterial communities in unfertilized bulk soils, this order was in greatest relative abundance compared to other orders within both fertilization treatments (Fig. S2).

FIG S2

Comparisons of bulk soil top OTU relative abundances (>1%) grouped by order. A single asterisk represents an indicator taxa for unfertilized treatment, and a double asterisk represents indicator taxa for fertilized plots (Table S5). Boxplots are colored according to fertilization treatment (gray for unfertilized, green for fertilized). Download FIG S2, PDF file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S5

Summary of bacterial taxa (OTUs) characteristic to each soil source and fertilization treatment based on indicator species analysis. Listed are the top OTUs that are significantly associated with each soil source and fertilization treatment group. Download Table S5, DOCX file, 0.01 MB.
Copyright © 2020 Bledsoe et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Comparisons of rhizosphere bacterial OTU presence/absence data revealed that forb (1,249 OTUs) and grass (1,019 OTUs) rhizospheres have distinct but overlapping microbiomes. Of the 1,621 total OTUs found in rhizosphere soils, 647 are broadly distributed and are observed in all plant rhizospheres and bulk soils regardless of treatment. Therefore, less than half of the forb (48%) and grass (37%) rhizosphere members were unique to that plant functional type, and broadly distributed OTUs dominate plant microbiomes especially in grasses.

Of OTUs that were represented only in the grass microbiome (n = 372), only 22 bacterial families are represented at > 0.075% relative abundance. Within those top OTUs, unfertilized grass rhizospheres were enriched in 9 families, while fertilized plots were enriched in 19 families (Fig. 4). Indicator species for unfertilized grass rhizospheres included 2 OTUs, one in the genus Singulisphaera and family Planctomycetaceae (indicator value [IndVal] = 0.38, P = 0.026) and an unclassified Spartobacteria OTU (IndVal = 0.44, P = 0.008; Table S5). Indicator species for fertilized grass rhizospheres included two OTUs, one in the genus Planctomyces and family Planctomycetaceae (IndVal = 0.42, P = 0.011) and one in the genus Actinoallomurus and family Thermomonosporaceae (IndVal = 0.36, P = 0.045; Table S5).

FIG 4
FIG 4

Comparisons of top OTU relative abundances (>0.075%) at the family level between fertilization treatments for grass rhizosphere bacterial communities. Asterisks represent indicator species present within the family (see Table S5 in the supplemental material). Colors indicate relative abundance increases from cool to warm (green yellow, orange, and red). White boxes indicate taxa present at <0.075% relative abundance.

Of the OTUs that were represented only in the forb microbiome (n = 602), only 21 bacterial families are represented at >0.1% relative abundance. Within those top OTUs, unfertilized forb rhizospheres were enriched in 10 families, while fertilized plots were enriched in 16 families (Fig. 5). Indicator species for unfertilized forb rhizospheres included two OTUs, Acidobacteria Gp1 (IndVal = 0.42, P = 0.02), and an unclassified Proteobacteria (IndVal = 0.46, P = 0.033; Table S5). Indicator species for fertilized forb rhizospheres included an OTU in Acidobacteria Gp1 (IndVal = 0.34, P = 0.041) class and an unclassified bacterial OTU (IndVal = 0.60, P = 0.017; Table S5).

FIG 5
FIG 5

Comparisons of top OTU relative abundances (>0.1%) at the family level between fertilization treatments for forb rhizosphere bacterial communities. Asterisks represent indicator species present within the family (Table S5). Colors indicate relative abundance increases from cool to warm (green yellow, orange, and red). White boxes indicate taxa present at <0.1% relative abundance.

DISCUSSION

In this study, nutrient addition increased bacterial species diversity (H′) and richness in bulk and rhizosphere soils. These results were similar to the results of O’Brien et al. (51) but contrary to our prediction and the results of other studies (14–16). Overall, bulk soils had the greatest bacterial diversity and highest pH values compared to rhizosphere soils. Since pH is known to be a strong driver of bacterial diversity, which can have a positive relationship with pH (52, 53), this increase in diversity may be due, in part, to the greater bulk soil pH compared to rhizosphere soil pH. The difference in pH between soil types is possibly due to organic acids in plant root exudates released into the rhizosphere (41); however, we did not analyze the composition of root exudates. Additionally, pH tended to be lower in unfertilized treatments, and diversity was more strongly related to pH in unfertilized soils than in fertilized soils. This may be due to the sensitivity of bacteria to acidic soils (53). The increase in bacterial diversity is likely the result of soil pH and niche differentiation due to fertilization increasing nutrient availability and rhizodeposition by plants, which introduces organic C resources for heterotrophs (17, 32). In dilution to extinction experiments, decreases in microbial diversity can result in loss of microbial functional diversity (54, 55). Therefore, increases in microbial diversity could result in increased microbial functional diversity, which could increase C cycling and promote N mining particularly in plant rhizospheres (56).

Bacterial taxa identified in rhizosphere samples are putatively involved in nutrient cycling and disease-suppressive functions. For example, fertilized forb rhizospheres were enriched in taxa from the family Streptomycetaceae, of which many produce antibiotics (57) and Sphingomonadaceae, which include taxa with disease suppression potential against fungal pathogens (58) (Fig. 5). This increase in disease-suppressive bacterial taxa suggests a potential increase in plant-pathogenic taxa within fertilized rhizospheres; however, this study did not specifically address disease suppression in soils. In contrast, fertilized grass rhizospheres were enriched with taxa putatively involved in N2 fixation (Acetobacteraceae) (59) and also Chitiniphagaceae and Conexibacteraceae, which have been implicated in decomposition of recalcitrant C sources (60, 61) (Fig. 4). Bacterial taxa in the Xanthomonadaceae family, which have previously been found in environments containing glyphosate (62), and Caulobacteraceae, which grows optimally on pesticides (63), are also more abundant in fertilized grass rhizospheres (Fig. 4). Since fertilization increased bacterial diversity and shifted composition, it is possible that fertilization has stimulated root exudation. The relative increase in complex C-degrading bacterial taxa in the grass rhizosphere could also be due to greater inputs of phenolics and terpenoids used as allelochemicals by the plant as revealed in past studies (64, 65). These differences in bacterial composition between the two plant species could be due to differences in composition of root exudates released into the rhizosphere (37); however, we did not analyze the composition of root exudates in the present study. Together, results suggest that nutrient addition enriches forb rhizospheres with putatively disease-suppressive bacteria and grass rhizospheres with taxa capable of decomposing complex C sources.

Within bulk soil bacterial members, putative nitrogen cycling taxa in the order Rhizobiales were enriched across all fertilization treatments (66, 67). This is not surprising, considering the limited amount of nitrogen in both unfertilized and fertilized soils at the study site. Despite the increase in taxa capable of N2 fixation in fertilized rhizospheres, these bacteria will acquire soil N if it is available (68). Therefore, these taxa may be less cooperative with plant associates than the same taxa from unfertilized soils, thereby reducing plant benefit (2, 69). This was not specifically tested in this study but could be an important future research topic.

Contrary to our prediction, bulk soils had a higher copiotroph-to-oligotroph ratio (based on rrn gene copy number) than rhizospheres. Characteristic of the copiotrophic life history strategy is the ability to rapidly decompose labile C sources; therefore, we expected that C-rich root exudates in the rhizosphere would support higher proportions of copiotrophic species (17). Additionally, fertilization did not increase the relative abundance of copiotrophic taxa. Rather, the observed copiotroph-to-oligotroph ratios were low in all samples with unfertilized bulk soils having the greatest proportion (0.22) and unfertilized grass rhizospheres having the lowest (0.13) copiotroph-to-oligotroph ratios. We suggest that the dominance of oligotrophs reflects the low-nutrient history of this wetland (29, 70), which is in contrast to agricultural systems that undergo regular fertilization at target rates intended to support high nutrient requirements for enhanced crop production (e.g., corn).

These results are in contrast to our first hypothesis and in agreement with our second hypothesis. Analyses of bacterial diversity and copiotroph-to-oligotroph ratios revealed an increase in bacterial diversity in response to fertilization and dominance of oligotrophs across all treatments within the study wetland. The low-nutrient history of the study site is likely the primary factor shaping bacterial community composition within the wetland. In agreement with our second hypothesis, comparisons of bulk and rhizosphere bacterial communities revealed that rhizospheres were more similar to each other than to bulk soil bacterial communities within fertilization treatments. Core plant microbiomes were predominantly composed of broadly distributed taxa; therefore, changes in bulk soil bacterial composition due to nutrient enrichment can directly alter plant microbiome composition and indirectly diminish benefits to plants if nutrient enrichment selects for more competitive bacterial taxa. These results highlight the importance of bulk soils as reservoirs of diversity for plant rhizospheres, which could have further implications for agricultural plant species in maintaining beneficial microbial communities.

Overall, this study revealed that long-term fertilization of oligotroph-dominated soils in low-nutrient wetlands increases bacterial species diversity. This increase in bacterial diversity has the potential to result in increased C and nutrient cycling that could lead to declines of wetland C storage potential. Nutrient enrichment also differentially alters plant rhizosphere composition in a way that suggests metabolic changes within soil bacterial communities. These metabolic changes could indirectly impact plant species diversity by providing an advantage to one species versus another through disease suppression or by increasing plant-available N through promotion of soil organic matter decomposition. If indirect fertilization supports rhizosphere bacterial communities that can enhance recalcitrant or labile C decomposition, wetland C storage potential could decline. Based on this study, bacterial taxonomic characterization sheds light on fertilization effects on plant-bacterial relationships. As such, nutrient enrichment effects on the metabolic diversity of bacterial communities could be even more pronounced in naturally low-nutrient ecosystems and warrants further investigation.

MATERIALS AND METHODS

Study site and experimental design.A long-term experimental site was established in 2003 to test the effects of fertilization, mowing, and the interaction on wetland plant communities. The site is located at East Carolina University’s West Research Campus in Greenville, NC, USA (35.6298N, 77.4836W). A description of the study site and experimental design can be found in Goodwillie and Franch (71) and is summarized here. This site is classified as a jurisdictional wetland but historically described as a mosaic of wet pine flatwood habitat, pine savanna, and hardwood communities. Soils were characterized as fine, kaolinitic, thermic Typic Paleaquults (Coxville series) with a fine sandy loam texture which are ultisols that are acidic and moderate to poorly drained soil types (https://soilseries.sc.egov.usda.gov/osdname.aspx). The annual mean temperature is 17.2°C, and annual precipitation is 176 cm (https://www.climate.gov/maps-data/dataset/). Treatments are replicated on eight 20 × 30 m blocks, and the nitrogen-phosphorus-potassium (N-P-K) 10-10-10 pellet fertilizer is applied three times per year (February, June, and October) for a total annual supplementation of 45.4 kg ha−1 for each nutrient. Plots are mowed by bush-hog and raked annually to simulate a fire disturbance (71).

We compared rhizosphere and bulk soil microbiomes in mowed unfertilized and fertilized plots, where herbaceous species dominated. Soil samples were collected at mowed/unfertilized and mowed/fertilized plots in four out of eight replicate blocks to reduce variability due to hydrology. Half the site is located adjacent to a ditch (drier soils) compared to away from the ditch, where soil conditions are wetter. Since this hydrologic gradient has resulted in distinct plant communities (C. Goodwillie, M.W. McCoy, and A. L. Peralta, submitted for publication), we collected samples from the wetter plots (away from the drainage ditch).

Bulk and rhizosphere soil sampling.We collected soil samples on 29 September 2015, approximately 3 months after the last fertilization treatment. Due to annual mowing and raking in sample plots, there was limited biomass accumulated in the organic horizon. We focused soil sampling and analysis on the mineral horizon. For a single composite bulk soil sample, we collected two soil cores (12-cm depth, 3.1-cm diameter) near each of the three permanently installed 1-m2 quadrats used for annual plant surveys. Each composite bulk soil sample was homogenized and passed through a 4-mm sieve, and any plant material was removed before further analysis. At each plot, rhizosphere soils were collected from the C3 forb Euthamia caroliniana (L.) Greene ex Porter & Britton and C4 grass Andropogon virginicus L. Rhizosphere soils were a composite of three root systems of the same species. Roots were gently dislodged from soil and neighboring roots and placed in a paper bag. After vigorous shaking, soil in the bag was processed for abiotic analysis. The roots were placed into 50-ml centrifuge tubes with 30 ml sterilized Nanopure water and shaken at 100 rpm for 1 h. Washed roots were removed, and the soil and water mixture was freeze dried to remove water. Freeze-dried rhizosphere samples were stored at –80°C until DNA extraction.

Soil chemical and physical characteristics.We measured gravimetric soil moisture by drying 20 to 30 g of field soil (moist soil) at 105°C for 24 h. We calculated percent moisture as the difference in weight of moist and dried soils divided by the oven-dried soil weight. Oven-dried samples were ground and measured for pH by mixing a 1:1 (soil-water) solution. A subsample of oven-dried soil was sieved with a 500-μm mesh and analyzed for total carbon and total nitrogen using an elemental analyzer (2400 CHNS Analyzer; Perkin-Elmer, Waltham, MA, USA) at the Environmental and Agricultural Testing Service laboratory (Department of Crop and Soil Sciences at NC State). Approximately 5 g of field soil (moist soil) was extracted with 45 ml of 2 M KCl, and available ammonium (NH4+) and nitrate (NO3−) ions were colorimetrically measured using a SmartChem 200 auto analyzer (Unity Scientific, Milford, MA, USA) at the East Carolina University Environmental Resources Laboratory.

Bacterial community analyses.We extracted DNA from soils using the Qiagen DNeasy PowerSoil kit. We used this DNA as the template in PCRs using barcoded primers (bacterial 515FB/806R) originally developed by the Earth Microbiome Project to target the V4 region of the bacterial 16S subunit of the rRNA gene (72). For each sample, three 50-μl PCR libraries were prepared by combining 30.75 μl molecular grade water, 5 μl Perfect Taq 10× buffer, 10 μl Perfect Taq 5× buffer, 1 μl deoxynucleoside triphosphates (dNTPs) (40 mM total, 10 mM each), 0.25 μl Perfect Taq polymerase, 1 μl forward barcoded primer (10 μM), 1 μl reverse primer (10 μM), and 1 μl DNA template (10 ng μl−1). Thermocycler conditions for PCRs were as follows: (i) initial denaturation (94°C for 3 min); (ii) 30 cycles with 1 cycle consisting of 94°C for 45 s, 50°C for 30 s, and 72°C for 90 s; (iii) final elongation (72°C, 10 min). Triplicate PCRs were combined and cleaned using the AMPure XP magnetic bead protocol (Axygen, Union City, CA, USA). Cleaned PCR product were quantified using QuantIT dsDNA BR assay (Thermo Scientific, Waltham, MA, USA) and diluted to a concentration of 10 ng μl−1 before pooling libraries in equimolar concentration of 5 ng μl−1. We sequenced pooled libraries using the Illumina MiSeq platform using paired-end reads (Illumina reagent kit v2, 500 reaction kit) at the Indiana University Center for Genomics and Bioinformatics Sequencing Facility. Sequences were processed using mothur (v1.40.1) (73) MiSeq pipeline (74). We assembled contigs from the paired-end reads, quality trimmed using a moving average quality score (minimum quality score, 35), aligned sequences to the SILVA rRNA database (v128) (75), and removed chimeric sequences using the VSEARCH algorithm (76). We created operational taxonomic units (OTUs) by first splitting sequences based on taxonomic class and then binning into OTUs based on 97% sequence similarity. The SILVA rRNA database was then used to assign taxonomic designations to OTUs.

Samples were rarefied to 43,811 reads (lowest read count among all samples) and resampled. We used vegan::diversity (77) to calculate bacterial species diversity by calculating Shannon diversity index (H´) because it accounts for species abundance and evenness and rare species (78, 79). We estimated bacterial richness using Chao1 species richness because it is nonparametric and also considers rare species (79, 80). Shannon diversity was calculated using the vegan::diversity function and Chao1 OTU richness using vegan::estimate (77). We assigned gene copy number to each OTU using RDP classifier (v2.12) (81) integrated with the rrn operon database developed by the Schimdt Laboratory at the Michigan Microbiome Project, University of Michigan (23, 27). Higher gene copy numbers (≥5) represent the copiotrophic lifestyle, and lower gene copy numbers (<5) represent the oligotrophic lifestyle (20, 24, 82). The numbers of copiotrophs and oligotrophs were summed for each soil sample to calculate the copiotroph-to-oligotroph ratio within a soil bacterial community.

Statistical analyses.All statistical analyses were performed in the R statistical environment (RStudio v1.1.383, Rv3.4.0) (83). We used a two-way model of analysis of variance (ANOVA) to compare the main effects of soil source, fertilization treatment, and the interaction to test for differences in OTU diversity and richness, copiotroph-to-oligotroph ratios, and soil parameters (soil pH, total carbon, total nitrogen, extractable ammonium and nitrate, and soil moisture). Significant interactions were compared with Tukey’s post hoc analysis using the agricolae::HSD.test R function (84). We examined diversity by visualizing bacterial community responses to fertilization and rhizosphere association using principal-coordinate analysis (PCoA) based on Bray-Curtis dissimilarity. We used permutational multivariate analysis of variance (PERMANOVA) to test for differences in bacterial community composition among treatments and within treatment using pairwise comparisons. Hypothesis testing using PERMANOVA was performed using the vegan::adonis function (77). We examined the relationship between soil parameters and bacterial Bray-Curtis dissimilarity patterns using the vegan::envfit function (77). Soil parameters with P < 0.05 were represented on the PCoA plot as vectors scaled by strength of correlation. We performed Dufrene-Legendre indicator species analysis using the labdsv::indval function (85) to identify specific community members that represented each soil source and fertilization treatment combination.

Data availability.An .Rmd file that includes annotated code and R scripts, data files, and the sequence processing batch file used in this study can be found in a public GitHub repository at https://github.com/PeraltaLab/WRC15_Rhizo. Raw sequence files for each sample can be accessed through the NCBI SRA BioProject accession number PRJNA599142.

ACKNOWLEDGMENTS

We thank M. Beamon, C. Eakins, J. LeCrone, J. Stiller, and S. Wilkinson for laboratory and field assistance. We thank J. Gill and the East Carolina University grounds crew for their efforts in maintaining the long-term ecological experiment. We also thank two anonymous reviewers for their suggestions that greatly improved this manuscript.

This work was supported by the National Science Foundation (GRFP to R.B.B. and DEB 1845845 to A.L.P.) and East Carolina University.

R.B.B., C.G., and A.L.P. conceived and designed the research. R.B.B. collected and analyzed the data. R.B.B. wrote the manuscript. All authors performed field work and edited the manuscript.

FOOTNOTES

    • Received January 13, 2020.
    • Accepted May 7, 2020.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. van der Heijden MGA,
    2. Bardgett RD,
    3. van Straalen NM
    . 2008. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol Lett 11:296−310. doi:10.1111/j.1461-0248.2007.01139.x.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Weese DJ,
    2. Heath KD,
    3. Dentinger BTM,
    4. Lau JA
    . 2015. Long-term nitrogen addition causes the evolution of less-cooperative mutualists. Evolution 69:631–642. doi:10.1111/evo.12594.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Regus JU,
    2. Wendlandt CE,
    3. Bantay RM,
    4. Gano-Cohen KA,
    5. Gleason NJ,
    6. Hollowell AC,
    7. O’Neill MR,
    8. Shahin KK,
    9. Sachs JL
    . 2017. Nitrogen deposition decreases the benefits of symbiosis in a native legume. Plant Soil 414:159–170. doi:10.1007/s11104-016-3114-8.
    OpenUrlCrossRef
  4. 4.↵
    1. Jach-Smith LC,
    2. Jackson RD
    . 2018. N addition undermines N supplied by arbuscular mycorrhizal fungi to native perennial grasses. Soil Biol Biochem 116:148–157. doi:10.1016/j.soilbio.2017.10.009.
    OpenUrlCrossRef
  5. 5.↵
    1. Fowler D,
    2. Coyle M,
    3. Skiba U,
    4. Sutton MA,
    5. Cape JN,
    6. Reis S,
    7. Sheppard LJ,
    8. Jenkins A,
    9. Grizzetti B,
    10. Galloway JN,
    11. Vitousek P,
    12. Leach A,
    13. Bouwman AF,
    14. Butterbach-Bahl K,
    15. Dentener F,
    16. Stevenson D,
    17. Amann M,
    18. Voss M
    . 2013. The global nitrogen cycle in the twenty-first century. Philos Trans R Soc Lond B Biol Sci 368:20130164. doi:10.1098/rstb.2013.0164.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Guignard MS,
    2. Leitch AR,
    3. Acquisti C,
    4. Eizaguirre C,
    5. Elser JJ,
    6. Hessen DO,
    7. Jeyasingh PD,
    8. Neiman M,
    9. Richardson AE,
    10. Soltis PS,
    11. Soltis DE,
    12. Stevens CJ,
    13. Trimmer M,
    14. Weider LJ,
    15. Woodward G,
    16. Leitch IJ
    . 2017. Impacts of nitrogen and phosphorus: from genomes to natural ecosystems and agriculture. Front Ecol Evol 5:1–9.
    OpenUrl
  7. 7.↵
    1. Wang R,
    2. Balkanski Y,
    3. Boucher O,
    4. Ciais P,
    5. Peñuelas J,
    6. Tao S
    . 2015. Significant contribution of combustion-related emissions to the atmospheric phosphorus budget. Nature Geosci 8:48–54. doi:10.1038/ngeo2324.
    OpenUrlCrossRef
  8. 8.↵
    1. Galloway JN,
    2. Dentener FJ,
    3. Capone DG,
    4. Boyer EW,
    5. Howarth RW,
    6. Seitzinger SP,
    7. Asner GP,
    8. Cleveland CC,
    9. Green PA,
    10. Holland EA,
    11. Karl DM,
    12. Michaels AF,
    13. Porter JH,
    14. Townsend AR,
    15. Vörösmarty CJ
    . 2004. Nitrogen cycles: past, present, and future. Biogeochemistry 70:153–226. doi:10.1007/s10533-004-0370-0.
    OpenUrlCrossRef
  9. 9.↵
    1. Harpole WS,
    2. Sullivan LL,
    3. Lind EM,
    4. Firn J,
    5. Adler PB,
    6. Borer ET,
    7. Chase J,
    8. Fay PA,
    9. Hautier Y,
    10. Hillebrand H,
    11. MacDougall AS,
    12. Seabloom EW,
    13. Williams R,
    14. Bakker JD,
    15. Cadotte MW,
    16. Chaneton EJ,
    17. Chu C,
    18. Cleland EE,
    19. D’Antonio C,
    20. Davies KF,
    21. Gruner DS,
    22. Hagenah N,
    23. Kirkman K,
    24. Knops JMH,
    25. La Pierre KJ,
    26. McCulley RL,
    27. Moore JL,
    28. Morgan JW,
    29. Prober SM,
    30. Risch AC,
    31. Schuetz M,
    32. Stevens CJ,
    33. Wragg PD
    . 2016. Addition of multiple limiting resources reduces grassland diversity. Nature 537:93–96. doi:10.1038/nature19324.
    OpenUrlCrossRef
  10. 10.↵
    1. WallisDeVries MF,
    2. Bobbink R
    . 2017. Nitrogen deposition impacts on biodiversity in terrestrial ecosystems: mechanisms and perspectives for restoration. Biol Conserv 212:387–389. doi:10.1016/j.biocon.2017.01.017.
    OpenUrlCrossRef
  11. 11.↵
    1. Dickson TL,
    2. Foster BL
    . 2011. Fertilization decreases plant biodiversity even when light is not limiting. Ecol Lett 14:380–388. doi:10.1111/j.1461-0248.2011.01599.x.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Song L,
    2. Bao X,
    3. Liu X,
    4. Zhang Y,
    5. Christie P,
    6. Fangmeier A,
    7. Zhang F
    . 2011. Nitrogen enrichment enhances the dominance of grasses over forbs in a temperate steppe ecosystem. Biogeosciences 8:2341–2350. doi:10.5194/bg-8-2341-2011.
    OpenUrlCrossRefWeb of Science
  13. 13.↵
    1. Stevens CJ,
    2. Dise NB,
    3. Gowing DJG,
    4. Mountford JO
    . 2006. Loss of forb diversity in relation to nitrogen deposition in the UK: regional trends and potential controls. Global Change Biol 12:1823–1833. doi:10.1111/j.1365-2486.2006.01217.x.
    OpenUrlCrossRef
  14. 14.↵
    1. Zeng J,
    2. Liu X,
    3. Song L,
    4. Lin X,
    5. Zhang H,
    6. Shen C,
    7. Chu H
    . 2016. Nitrogen fertilization directly affects soil bacterial diversity and indirectly affects bacterial community composition. Soil Biol Biochem 92:41–49. doi:10.1016/j.soilbio.2015.09.018.
    OpenUrlCrossRef
  15. 15.↵
    1. Wang C,
    2. Liu D,
    3. Bai E
    . 2018. Decreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition. Soil Biol Biochem 120:126–133. doi:10.1016/j.soilbio.2018.02.003.
    OpenUrlCrossRef
  16. 16.↵
    1. Zhou J,
    2. Jiang X,
    3. Wei D,
    4. Zhao B,
    5. Ma M,
    6. Chen S,
    7. Cao F,
    8. Shen D,
    9. Guan D,
    10. Li J
    . 2017. Consistent effects of nitrogen fertilization on soil bacterial communities in black soils for two crop seasons in China. Sci Rep 7:3267. doi:10.1038/s41598-017-03539-6.
    OpenUrlCrossRef
  17. 17.↵
    1. Bulgarelli D,
    2. Schlaeppi K,
    3. Spaepen S,
    4. Ver Loren van Themaat E,
    5. Schulze-Lefert P
    . 2013. Structure and functions of the bacterial microbiota of plants. Annu Rev Plant Biol 64:807–838. doi:10.1146/annurev-arplant-050312-120106.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Mendes LW,
    2. Kuramae EE,
    3. Navarrete AA,
    4. Van Veen JA,
    5. Tsai SM
    . 2014. Taxonomical and functional microbial community selection in soybean rhizosphere. ISME J 8:1577–1587. doi:10.1038/ismej.2014.17.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. de Ridder-Duine AS,
    2. Kowalchuk GA,
    3. Klein Gunnewiek PJA,
    4. Smant W,
    5. Van Veen JA,
    6. De Boer W
    . 2005. Rhizosphere bacterial community composition in natural stands of Carex arenaria (sand sedge) is determined by bulk soil community composition. Soil Biol Biochem 37:349–357. doi:10.1016/j.soilbio.2004.08.005.
    OpenUrlCrossRefWeb of Science
  20. 20.↵
    1. Fierer N,
    2. Bradford MA,
    3. Jackson RB
    . 2007. Toward an ecological classification of soil bacteria. Ecology 88:1354–1364. doi:10.1890/05-1839.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Leff JW,
    2. Jones SE,
    3. Prober SM,
    4. Barberán A,
    5. Borer ET,
    6. Firn JL,
    7. Harpole WS,
    8. Hobbie SE,
    9. Hofmockel KS,
    10. Knops JMH,
    11. McCulley RL,
    12. La Pierre K,
    13. Risch AC,
    14. Seabloom EW,
    15. Schütz M,
    16. Steenbock C,
    17. Stevens CJ,
    18. Fierer N
    . 2015. Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe. Proc Natl Acad Sci U S A 112:10967–10972. doi:10.1073/pnas.1508382112.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Goldfarb KC,
    2. Karaoz U,
    3. Hanson CA,
    4. Santee CA,
    5. Bradford MA,
    6. Treseder KK,
    7. Wallenstein MD,
    8. Brodie EL
    . 2011. Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance. Front Microbiol 2:94. doi:10.3389/fmicb.2011.00094.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Roller BRK,
    2. Stoddard SF,
    3. Schmidt TM
    . 2016. Exploiting rRNA operon copy number to investigate bacterial reproductive strategies. Nat Microbiol 1:16160. doi:10.1038/nmicrobiol.2016.160.
    OpenUrlCrossRef
  24. 24.↵
    1. Klappenbach JA,
    2. Dunbar JM,
    3. Schmidt TM
    . 2000. rRNA operon copy number reflects ecological strategies of bacteria. Appl Environ Microbiol 66:1328–1333. doi:10.1128/aem.66.4.1328-1333.2000.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Yano K,
    2. Wada T,
    3. Suzuki S,
    4. Tagami K,
    5. Matsumoto T,
    6. Shiwa Y,
    7. Ishige T,
    8. Kawaguchi Y,
    9. Masuda K,
    10. Akanuma G,
    11. Nanamiya H,
    12. Niki H,
    13. Yoshikawa H,
    14. Kawamura F
    . 2013. Multiple rRNA operons are essential for efficient cell growth and sporulation as well as outgrowth in Bacillus subtilis. Microbiology 159:2225–2236. doi:10.1099/mic.0.067025-0.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Stevenson BS,
    2. Schmidt TM
    . 2004. Life history implications of rRNA gene copy number in Escherichia coli. Appl Environ Microbiol 70:6670–6677. doi:10.1128/AEM.70.11.6670-6677.2004.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Stoddard SF,
    2. Smith BJ,
    3. Hein R,
    4. Roller BRK,
    5. Schmidt TM
    . 2015. rrnDB: improved tools for interpreting rRNA gene abundance in bacteria and archaea and a new foundation for future development. Nucleic Acids Res 43:D593–D598. doi:10.1093/nar/gku1201.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Francioli D,
    2. Schulz E,
    3. Lentendu G,
    4. Wubet T,
    5. Buscot F,
    6. Reitz T
    . 2016. Mineral vs. organic amendments: microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilization strategies. Front Microbiol 7:1446. doi:10.3389/fmicb.2016.01446.
    OpenUrlCrossRef
  29. 29.↵
    1. Ho A,
    2. Di Lonardo DP,
    3. Bodelier P
    . 2017. Revisiting life strategy concepts in environmental microbial ecology. FEMS Microbiol Ecol 93:fix006. doi:10.1093/femsec/fix006.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Zarraonaindia I,
    2. Owens S,
    3. Weisenhorn P,
    4. West K,
    5. Hampton-Marcell J,
    6. Lax S,
    7. Bokulich N,
    8. Mills D,
    9. Martin G,
    10. Taghavi S,
    11. van der Lelie D,
    12. Gilbert J
    . 2015. The soil microbiome influences grapevine-associated microbiota. mBio 6:e02527-14. [CrossRef] doi:10.1128/mBio.02527-14.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Matthews A,
    2. Pierce S,
    3. Hipperson H,
    4. Raymond B
    . 2019. Rhizobacterial community assembly patterns vary between crop species. Front Microbiol 10:581. doi:10.3389/fmicb.2019.00581.
    OpenUrlCrossRef
  32. 32.↵
    1. Lundberg DS,
    2. Lebeis SL,
    3. Paredes SH,
    4. Yourstone S,
    5. Gehring J,
    6. Malfatti S,
    7. Tremblay J,
    8. Engelbrektson A,
    9. Kunin V,
    10. del Rio TG,
    11. Edgar RC,
    12. Eickhorst T,
    13. Ley RE,
    14. Hugenholtz P,
    15. Tringe SG,
    16. Dangl JL
    . 2012. Defining the core Arabidopsis thaliana root microbiome. Nature 488:86–90. doi:10.1038/nature11237.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Uroz S,
    2. Buée M,
    3. Murat C,
    4. Frey-Klett P,
    5. Martin F
    . 2010. Pyrosequencing reveals a contrasted bacterial diversity between oak rhizosphere and surrounding soil. Environ Microbiol Rep 2:281–288. doi:10.1111/j.1758-2229.2009.00117.x.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Philippot L,
    2. Raaijmakers JM,
    3. Lemanceau P,
    4. Van Der Putten WH
    . 2013. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol 11:789–799. doi:10.1038/nrmicro3109.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Hütsch BW,
    2. Augustin J,
    3. Merbach W
    . 2002. Plant rhizodeposition — an important source for carbon turnover in soils. J Plant Nutr Soil Sci 165:397–407. doi:10.1002/1522-2624(200208)165:4<397::AID-JPLN397>3.0.CO;2-C.
    OpenUrlCrossRef
  36. 36.↵
    1. Jones DL,
    2. Hodge A,
    3. Kuzyakov Y
    . 2004. Plant and mycorrhizal regulation of rhizodeposition. New Phytol 163:459–480. doi:10.1111/j.1469-8137.2004.01130.x.
    OpenUrlCrossRefWeb of Science
  37. 37.↵
    1. el Zahar Haichar F,
    2. Santaella C,
    3. Heulin T,
    4. Achouak W
    . 2014. Root exudates mediated interactions belowground. Soil Biol Biochem 77:69–80. doi:10.1016/j.soilbio.2014.06.017.
    OpenUrlCrossRef
  38. 38.↵
    1. Jones DL,
    2. Nguyen C,
    3. Finlay RD
    . 2009. Carbon flow in the rhizosphere: carbon trading at the soil-root interface. Plant Soil 321:5–33. doi:10.1007/s11104-009-9925-0.
    OpenUrlCrossRefWeb of Science
  39. 39.↵
    1. Kuzyakov Y,
    2. Domanski G
    . 2000. Carbon input by plants into the soil. J Plant Nutr Soil Sci 163:421–431. doi:10.1002/1522-2624(200008)163:4<421::AID-JPLN421>3.0.CO;2-R.
    OpenUrlCrossRef
  40. 40.↵
    1. Bertin C,
    2. Yang X,
    3. Weston LA
    . 2003. The role of root exudates and allelochemicals in the rhizosphere. Plant Soil 256:67–83. doi:10.1023/A:1026290508166.
    OpenUrlCrossRefWeb of Science
  41. 41.↵
    1. Dakora FD,
    2. Phillips DA
    . 2002. Root exudates as mediators of mineral acquisition in low-nutrient environments. Plant Soil 245:35–47. doi:10.1023/A:1020809400075.
    OpenUrlCrossRefWeb of Science
  42. 42.↵
    1. Kellogg EA
    . 2013. C4 photosynthesis. Curr Biol 23:R594–R599. doi:10.1016/j.cub.2013.04.066.
    OpenUrlCrossRef
  43. 43.↵
    1. Schmitt MR,
    2. Edwards GE
    . 1981. Photosynthetic capacity and nitrogen use efficiency of maize, wheat, and rice: a comparison between C3 and C4 photosynthesis. J Exp Bot 32:459–466. doi:10.1093/jxb/32.3.459.
    OpenUrlCrossRefWeb of Science
  44. 44.↵
    1. Vranova V,
    2. Rejsek K,
    3. Skene KR,
    4. Janous D,
    5. Formanek P
    . 2013. Methods of collection of plant root exudates in relation to plant metabolism and purpose: a review. J Plant Nutr Soil Sci 176:175–199. doi:10.1002/jpln.201000360.
    OpenUrlCrossRef
  45. 45.↵
    1. Kuzyakov Y,
    2. Siniakina SV,
    3. Ruehlmann J,
    4. Domanski G,
    5. Stahr K
    . 2002. Effect of nitrogen fertilisation on below-ground carbon allocation in lettuce. J Sci Food Agric 82:1432–1441. doi:10.1002/jsfa.1202.
    OpenUrlCrossRefWeb of Science
  46. 46.↵
    1. Wu FY,
    2. Chung AKC,
    3. Tam NFY,
    4. Wong MH
    . 2012. Root exudates of wetland plants influenced by nutrient status and types of plant cultivation. Int J Phytoremediation 14:543–553. doi:10.1080/15226514.2011.604691.
    OpenUrlCrossRef
  47. 47.↵
    1. Yin H,
    2. Li Y,
    3. Xiao J,
    4. Xu Z,
    5. Cheng X,
    6. Liu Q
    . 2013. Enhanced root exudation stimulates soil nitrogen transformations in a subalpine coniferous forest under experimental warming. Glob Chang Biol 19:2158–2167. doi:10.1111/gcb.12161.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Phillips RP,
    2. Bernhardt ES,
    3. Schlesinger WH
    . 2009. Elevated CO2 increases root exudation from loblolly pine (Pinus taeda) seedlings as an N-mediated response. Tree Physiol 29:1513–1523. doi:10.1093/treephys/tpp083.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Uselman SM,
    2. Qualls RG,
    3. Thomas RB
    . 2000. Effects of increased atmospheric CO2, temperature, and soil N availability on root exudation of dissolved organic carbon by a N-fixing tree (Robinia pseudoacacia L.). Plant Soil 222:191–202. doi:10.1023/A:1004705416108.
    OpenUrlCrossRefWeb of Science
  50. 50.↵
    1. Carvalhais LC,
    2. Dennis PG,
    3. Fedoseyenko D,
    4. Hajirezaei MR,
    5. Borriss R,
    6. Von Wirén N
    . 2011. Root exudation of sugars, amino acids, and organic acids by maize as affected by nitrogen, phosphorus, potassium, and iron deficiency. J Plant Nutr Soil Sci 174:3–11. doi:10.1002/jpln.201000085.
    OpenUrlCrossRef
  51. 51.↵
    1. O’Brien SL,
    2. Gibbons SM,
    3. Owens SM,
    4. Hampton-Marcell J,
    5. Johnston ER,
    6. Jastrow JD,
    7. Gilbert JA,
    8. Meyer F,
    9. Antonopoulos DA
    . 2016. Spatial scale drives patterns in soil bacterial diversity. Environ Microbiol 18:2039–2051. doi:10.1111/1462-2920.13231.
    OpenUrlCrossRef
  52. 52.↵
    1. Rousk J,
    2. Bååth E,
    3. Brookes PC,
    4. Lauber CL,
    5. Lozupone C,
    6. Caporaso JG,
    7. Knight R,
    8. Fierer N
    . 2010. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J 4:1340–1351. doi:10.1038/ismej.2010.58.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    1. Fierer N,
    2. Jackson RB
    . 2006. The diversity and biogeography of soil bacterial communities. Proc Natl Acad Sci U S A 103:626–631. doi:10.1073/pnas.0507535103.
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    1. Philippot L,
    2. Spor A,
    3. Hénault C,
    4. Bru D,
    5. Bizouard F,
    6. Jones CM,
    7. Sarr A,
    8. Maron PA
    . 2013. Loss in microbial diversity affects nitrogen cycling in soil. ISME J 7:1609–1619. doi:10.1038/ismej.2013.34.
    OpenUrlCrossRefPubMedWeb of Science
  55. 55.↵
    1. Trivedi C,
    2. Delgado-Baquerizo M,
    3. Hamonts K,
    4. Lai K,
    5. Reich PB,
    6. Singh BK
    . 2019. Losses in microbial functional diversity reduce the rate of key soil processes. Soil Biol Biochem 135:267–274. doi:10.1016/j.soilbio.2019.05.008.
    OpenUrlCrossRef
  56. 56.↵
    1. Weidner S,
    2. Koller R,
    3. Latz E,
    4. Kowalchuk G,
    5. Bonkowski M,
    6. Scheu S,
    7. Jousset A
    . 2015. Bacterial diversity amplifies nutrient-based plant-soil feedbacks. Funct Ecol 29:1341–1349. doi:10.1111/1365-2435.12445.
    OpenUrlCrossRef
  57. 57.↵
    1. Kinkel LL,
    2. Schlatter DC,
    3. Bakker MG,
    4. Arenz BE
    . 2012. Streptomyces competition and co-evolution in relation to plant disease suppression. Res Microbiol 163:490–499. doi:10.1016/j.resmic.2012.07.005.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Chapelle E,
    2. Mendes R,
    3. Bakker P,
    4. Raaijmakers JM
    . 2016. Fungal invasion of the rhizosphere microbiome. ISME J 10:265–268. doi:10.1038/ismej.2015.82.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Saravanan VS,
    2. Madhaiyan M,
    3. Osborne J,
    4. Thangaraju M,
    5. Sa TM
    . 2008. Ecological occurrence of Gluconacetobacter diazotrophicus and nitrogen-fixing Acetobacteraceae members: their possible role in plant growth promotion. Microb Ecol 55:130–140. doi:10.1007/s00248-007-9258-6.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Deng J,
    2. Gu Y,
    3. Zhang J,
    4. Xue K,
    5. Qin Y,
    6. Yuan M,
    7. Yin H,
    8. He Z,
    9. Wu L,
    10. Schuur EAG,
    11. Tiedje JM,
    12. Zhou J
    . 2015. Shifts of tundra bacterial and archaeal communities along a permafrost thaw gradient in Alaska. Mol Ecol 24:222–234. doi:10.1111/mec.13015.
    OpenUrlCrossRef
  61. 61.↵
    1. Seki T,
    2. Matsumoto A,
    3. Shimada R,
    4. Inahashi Y
    . 2019. Conexibacter arvalis sp. nov., isolated from a cultivated field soil sample. Int J Syst Evol Biol 62:2400–2404. doi:10.1099/ijs.0.036095-0.
    OpenUrlCrossRef
  62. 62.↵
    1. Newman MM,
    2. Hoilett N,
    3. Lorenz N,
    4. Dick RP,
    5. Liles MR,
    6. Ramsier C,
    7. Kloepper JW
    . 2016. Glyphosate effects on soil rhizosphere-associated bacterial communities. Sci Total Environ 543:155–160. doi:10.1016/j.scitotenv.2015.11.008.
    OpenUrlCrossRef
  63. 63.↵
    1. Lingens F,
    2. Blecher R,
    3. Blecher H,
    4. Blobel F,
    5. Eberspacher J,
    6. Frohner C,
    7. Gorisch H,
    8. Gorisch H,
    9. Layh G
    . 1985. Phenylobacterium immobile gen. nov. sp. nov., a Gram-negative bacterium that degrades the herbicide chloridazon. Int J Syst Bacteriol 35:26–39. doi:10.1099/00207713-35-1-26.
    OpenUrlCrossRef
  64. 64.↵
    1. Rice EL
    . 1972. Allelopathic effects of Andropogon virginicus and its persistence in old fields. Am J Bot 59:752–755. doi:10.1002/j.1537-2197.1972.tb10148.x.
    OpenUrlCrossRef
  65. 65.↵
    1. Inderjit,
    2. Duke SO
    . 2003. Ecophysiological aspects of allelopathy. Planta 217:529–539. doi:10.1007/s00425-003-1054-z.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Anderson CR,
    2. Condron LM,
    3. Clough TJ,
    4. Fiers M,
    5. Stewart A,
    6. Hill RA,
    7. Sherlock RR
    . 2011. Biochar induced soil microbial community change: implications for biogeochemical cycling of carbon, nitrogen and phosphorus. Pedobiologia (Jena) 54:309–320. doi:10.1016/j.pedobi.2011.07.005.
    OpenUrlCrossRef
  67. 67.↵
    1. Rilling JI,
    2. Acuña JJ,
    3. Sadowsky MJ,
    4. Jorquera MA
    . 2018. Putative nitrogen-fixing bacteria associated with the rhizosphere and root endosphere of wheat plants grown in an andisol from southern Chile. Front Microbiol 9:2710. doi:10.3389/fmicb.2018.02710.
    OpenUrlCrossRef
  68. 68.↵
    1. Da Costa PB,
    2. Passaglia L
    . 2015. How fertilization affects the selection of plant growth promoting rhizobacteria by host plants, p 967–974. In de Bruijn FJ (ed), Biological nitrogen fixation. John Wiley & Sons, Inc, Hoboken, NJ.
  69. 69.↵
    1. Ai C,
    2. Liang G,
    3. Sun J,
    4. Wang X,
    5. He P,
    6. Zhou W,
    7. He X
    . 2015. Reduced dependence of rhizosphere microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils. Soil Biol Biochem 80:70–78. doi:10.1016/j.soilbio.2014.09.028.
    OpenUrlCrossRef
  70. 70.↵
    1. Song W,
    2. Kim M,
    3. Tripathi BM,
    4. Kim H,
    5. Adams JM
    . 2016. Predictable communities of soil bacteria in relation to nutrient concentration and successional stage in a laboratory culture experiment. Environ Microbiol 18:1740–1753. doi:10.1111/1462-2920.12879.
    OpenUrlCrossRef
  71. 71.↵
    1. Goodwillie C,
    2. Franch WR
    . 2006. An experimental study of the effects of nutrient addition and mowing on a ditched wetland plant community: results of the first year. J North Carolina Acad Sci 122:106–117.
    OpenUrl
  72. 72.↵
    1. Caporaso JG,
    2. Lauber CL,
    3. Walters WA,
    4. Berg-Lyons D,
    5. Huntley J,
    6. Fierer N,
    7. Owens SM,
    8. Betley J,
    9. Fraser L,
    10. Bauer M,
    11. Gormley N,
    12. Gilbert JA,
    13. Smith G,
    14. Knight R
    . 2012. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J 6:1621–1624. doi:10.1038/ismej.2012.8.
    OpenUrlCrossRefPubMedWeb of Science
  73. 73.↵
    1. Schloss PD,
    2. Westcott SL,
    3. Ryabin T,
    4. Hall JR,
    5. Hartmann M,
    6. Hollister EB,
    7. Lesniewski RA,
    8. Oakley BB,
    9. Parks DH,
    10. Robinson CJ,
    11. Sahl JW,
    12. Stres B,
    13. Thallinger GG,
    14. Van Horn DJ,
    15. Weber CF
    . 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537–7541. doi:10.1128/AEM.01541-09.
    OpenUrlAbstract/FREE Full Text
  74. 74.↵
    1. Kozich JJ,
    2. Westcott SL,
    3. Baxter NT,
    4. Highlander SK,
    5. Schloss PD
    . 2013. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol 79:5112–5120. doi:10.1128/AEM.01043-13.
    OpenUrlAbstract/FREE Full Text
  75. 75.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glöckner FO
    . 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590–D596. doi:10.1093/nar/gks1219.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Rognes T,
    2. Flouri T,
    3. Nichols B,
    4. Quince C,
    5. Mahé F
    . 2016. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi:10.7717/peerj.2584.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Oksanen J
    . 2019. R package ‘vegan,’ 2.5-6. https://cran.r-project.org/web/packages/vegan/index.html.
  78. 78.↵
    1. Shannon CE
    . 1948. A mathematical theory of communication. Bell Syst Tech J 27:379–423. doi:10.1002/j.1538-7305.1948.tb01338.x.
    OpenUrlCrossRef
  79. 79.↵
    1. Kim BR,
    2. Shin J,
    3. Guevarra RB,
    4. Lee JH,
    5. Kim DW,
    6. Seol KH,
    7. Lee JH,
    8. Kim HB,
    9. Isaacson RE
    . 2017. Deciphering diversity indices for a better understanding of microbial communities. J Microbiol Biotechnol 27:2089–2093. doi:10.4014/jmb.1709.09027.
    OpenUrlCrossRef
  80. 80.↵
    1. Chao A
    . 1984. Nonparametric estimation of the number of classes in a population. Scand J Stat 11:265–270.
    OpenUrlWeb of Science
  81. 81.↵
    1. Wang Q,
    2. Garrity GM,
    3. Tiedje JM,
    4. Cole JR
    . 2007. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 73:5261–5267. doi:10.1128/AEM.00062-07.
    OpenUrlAbstract/FREE Full Text
  82. 82.↵
    1. Roller BRK,
    2. Schmidt TM
    . 2015. The physiology and ecological implications of efficient growth. ISME J 9:1481–1487. doi:10.1038/ismej.2014.235.
    OpenUrlCrossRef
  83. 83.↵
    R Core Team. 2019. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
  84. 84.↵
    1. de Mendiburu F
    . 2019. R package ‘agricolae,’ 1.3–2. https://cran.r-project.org/web/packages/agricolae/index.html.
  85. 85.↵
    1. Roberts DW
    . 2019. R package ‘labdsv,’ 2.0-1. https://cran.r-project.org/web/packages/labdsv/index.html.
Back to top
Download PDF
Citation Tools
Print
Alerts
Email

KEYWORDS

copiotroph
fertilization
oligotroph
plant-microbe
rhizosphere

Related Articles

Cited By...