This study examined the abundance of bacteria and nirS-type denitrifiers associated with the rhizospheres of three emergent macrophyte species (Juncus effusus, Typha latifolia, and Peltandra virginica) to gain a greater understanding of plant-microbe interactions in wetland soils. Sampling of plant and soil properties was performed during the growing season (June) and following plant senescence (November) at two tidal freshwater wetlands. Quantitative polymerase chain reaction was used to determine the abundance of bacteria (16S rRNA) and nirS-type denitrifier genes from the rhizosphere and rhizoplane of each plant species and from nearby unvegetated soils. For bacteria, there was a positive rhizosphere effect that did not differ significantly across plant species. In contrast, significant differences in the abundance of nirS-type denitrifiers were observed across the plant species. Rhizosphere abundance was ∼2-fold greater in Peltandra virginica and 4-fold greater in Typha latifolia compared to Juncus effusus. For both bacteria and nirS-type denitrifiers, plant effects were greater during the growing season, and abundance was highly correlated with soil pH, moisture, and organic matter content. Overall, these results demonstrate plant effects on the rhizosphere microbial community can be species‐specific and that there is a synergistic relationship between plant species and environmental conditions.

1. Introduction

Wetlands are ecosystems of intense biogeochemical transformations, especially in regard to nitrogen cycling. They are estimated to remove nearly 20% of anthropogenic nitrogen inputs worldwide [1, 2]. One biogeochemical pathway of particular interest is denitrification, the stepwise reduction of nitrate (NO3) or nitrite (NO2) to a N-gas (N2O or N2), which is then lost to the atmosphere. At the ecosystem level, rates of nitrogen removal have been linked with broad-scale changes in the size and composition of wetland plant communities [36]. Most prior efforts to disentangle this relationship have focused on how plant communities, sediment properties, and hydrology interact, and considerably less is known about the small-scale interactions that take place in the rhizosphere and rhizoplane of individual plants. However, such information is essential if we are to develop a predictive understanding of ecosystem nitrogen dynamics based on plant functional characteristics and plant community composition [7].

Wetland plants affect soil biogeochemistry in a variety of ways, including via direct effects on microbial communities and indirect effects mediated through changes in soil physicochemical properties [8, 9]. For example, organic carbon release via root exudates has been found to affect microbial community composition [1013], influence overall microbial biomass [1416], and enhance rates of microbial activity [1719]. Since the quality and quantity of exudates can vary based on plant species, age, and health [2022], this is one mechanism by which differences in plant community composition can influence microbial community structure and function. Similarly, microbial communities may respond to changes in redox conditions and oxygen availability that arise from differences in radial oxygen loss (ROL) associated with the varying root morphologies and gas transport systems of wetland plants [3, 23]. These sorts of plant-mediated changes in carbon and oxygen availability are especially likely to affect denitrification, as it is primarily an anoxic heterotrophic process.

The metabolic pathways of denitrification are generally well studied. One of the most important steps is nitrite reduction, wherein NO2 is converted to a gaseous state as nitric oxide (NO). This step distinguishes denitrifiers from other nitrate-metabolizing microorganisms [24] and is regulated by the nitrite reductase enzyme. Denitrifiers carry one of two variants of the gene for nitrite reductase: nirS and nirK, which have both been used extensively as molecular markers for denitrifying bacteria [2529]. These genes are functionally equivalent, exhibit differential habitat selection [24, 30, 31], and, until very recently, were believed to be mutually exclusive [32]. The research presented here focuses on nirS, as it dominates our experimental system [33]. In particular, we studied how the abundance of bacteria and nirS-type denitrifiers changes in the rhizosphere and rhizoplane of three emergent freshwater macrophytes (Juncus effusus, Peltandra virginica, and Typha latifolia) in two tidal freshwater wetlands. Our main objectives were to determine the extent to which plant species enhanced nirS abundance and to evaluate whether plant effects on the microbial community are dependent upon soil conditions or season (i.e., growing season versus plant senescence).

2. Materials and Methods

2.1. Site Description

Two tidal freshwater wetlands were selected along the James River in the coastal plain of Virginia (USA), near Richmond. The first, referred to hereafter as “Rice,” is a 28-hectare wetland on the property of Virginia Commonwealth University’s Rice Rivers Center (37°19′38″ N, 77°12′13″ W). Originally forested wetland, Rice was dammed in the 1920s and used as a recreational lake for ∼80 years. Wetland soil conditions and vegetation began to naturally return to the site following a partial breach of the earthen dam in 2007. The remainder of the dam was removed in 2010 as part of a targeted restoration effort, which fully restored tidal flow and natural hydrology to the wetland. The second site (“James”) is a relatively undisturbed 1,700 ha wetland within the James River National Wildlife Refuge (37°15′51.1″ N, 77° 8′51.8″ W). Rice and James have similar soil texture, tidal regimes, and emergent macrophyte communities. Sampling was conducted at both sites during autumn senescence (November 2011) and during the peak of the growing season the following year (June 2012).

2.2. Experimental Design

This study focused on three macrophytes species, which differ in their root morphologies (Figure 1) and gas transport systems: Juncus effusus (common rush, hereafter Juncus), Typha latifolia (broadleaf cattail, hereafter Typha), and Peltandra virginica (arrow arum, hereafter Peltandra). Specifically, Typha has a rhizomatous root system while Juncas and Peltandra do not. Juncas utilizes a diffusive gas transport system [34] whereas Peltandra and Typha have pressured systems [35, 36], which tend to transport more oxygen belowground [37]. While rates of ROL can vary based on specific growth conditions, the potential of this particular Typha species to release oxygen is generally much higher than this Juncus [38].

For the November sampling, four replicate plants of each type were collected from each site. In June, five replicates were used. Only mature plants were included (based on height: Juncus: ∼100 cm, Peltandra: ∼60 cm, and Typha: ∼300 cm). In the field, each plant and its associated soil were collected via coring (25 cm diameter) to a depth of 30 cm, and intact plant-soil cores were transported upright to the laboratory in buckets. In addition, soil-only reference cores (8 cm diameter, 30 cm deep; 4 replicates in November and 5 in June) were collected from nearby unvegetated areas (i.e., no aboveground biomass within 30 cm) and transported back to the laboratory on ice. At Rice, all plant and soil samples were collected within a 15 × 30 m area of the wetland; at James, it was necessary to consider a slightly larger area (30 × 100 m) in order to encounter all three plant types.

2.3. Soil Properties

Upon return to the lab, soil-only reference cores were extracted into air-tight plastic bags, sealed, and gently homogenized by hand. Soil redox potential and pH were immediately measured using a Hanna Combo probe (QA Supplies, Norfolk, VA, USA). Soil samples were then subdivided for future analysis of soil properties (∼35 g, stored at 4°C up to 1 week) and the associated microbial community (∼15 g, stored at −20°C up to one 1 month).

Soil moisture content was determined gravimetrically (100 ± 5°C for 72 h), and organic matter content (%) was measured as the mass lost on ignition following combustion (400°C for 12 h). Total carbon (C) and nitrogen (N) content were determined using a Perkin Elmer CHNS/O Analyzer (Waltham, MA, USA) following the grinding and acidification of samples using 10% hydrochloric acid.

2.4. Rhizosphere and Rhizoplane Sampling

Immediately upon return to the lab, the intact root ball for each plant core was shaken gently by hand so that adhering bulk soil fell away. Rhizosphere samples were then obtained from the soil that remained immediately surrounding the roots of the plant. Five subsamples were collected along the depth of the core for each plant and pooled to yield a total of ∼35 g of soil per replicate plant. As with the unvegetated reference cores, the material was archived for analysis of soil properties and microbial communities.

For the June sampling event, rhizoplane samples were also collected from each plant for microbial community analyses. A modified version of the protocol outlined by Phillips et al. [39] was performed. First, several intact live roots were cut from the plant and shaken gently by hand to remove loose soil. These roots, along with any remaining root-bound soil, were then combined to generate one aggregate sample per plant (∼30 g, minimum of 5 root strands), mixed with 200 ml of sterile phosphate-buffered saline (PBS) buffer (pH = 7.4) in a 250-ml Erlenmeyer flask, and placed on a shaker table (170 rpm for 1 h) to produce a soil slurry. An additional 40 ml of PBS buffer was used to rinse the roots and flask, and the rinse was combined with the original slurry. The rhizoplane soil was then collected by centrifugation (1100 × g, 10 min), and pelleted soils were stored frozen for molecular genetic analysis (−20°C up to 1 month).

2.5. Plant Biomass

Plant biomass was only measured during the June sampling event; all data are reported as grams dry weight (DW) per plant. For belowground biomass, roots were harvested from each core, rinsed by hand, dried (70°C for 48 h), and weighed. Aboveground biomass was estimated after sorting and drying (100 ± 5°C for 72 h) live material. Subsamples of the aboveground material were archived for subsequent analysis of C : N of the leaf tissue.

2.6. Microbial Abundance

Whole-community DNA was extracted from 0.5 g subsamples of soil using MoBio PowerSoil DNA Isolation Kit (Carlsbad, CA, USA) and then stored at −20°C. DNA extractions were verified using agarose gel electrophoresis; purity and concentration were determined using a Nanodrop ND-1000 (Thermo Scientific, Wilmington, DE, USA).

To estimate total bacterial abundance, quantitative polymerase chain reaction (qPCR) was performed targeting the 16S rRNA gene with the Eub338 and Eub517 primers [40]. Genomic DNA from Escherichia coli (Strain #11775; ATCC, Manassas, VA, USA) was used for the standard curve (average efficiency = 98.8%, r2 = 0.995). Reactions (20 μl) were performed with 1.2 ng template DNA and 0.1 µM concentration of each primer; thermal cycling conditions were 95°C for 4 min followed by 40 cycles of 30 s at 95°C, 30 s at 55.5°C, and 60 s at 72°C.

A preliminary analysis using a subset of samples indicated nirK abundance (assessed via qPCR following Mosier and Francis [41]) was very low (101–103 copies per g dry weight of soil) and showed no significant plant or site effects; consequently, our analyses focused on nirS. Those qPCRs were performed using nirS-specific primers cd3aF and R3cd [26] and genomic DNA from Paracoccus denitrificans (ATCC Strain #17741) for the standard curve (average efficiency = 93.4%, r2 = 0.991). Reactions (20 µl) included 10 ng of template DNA and 0.1 µm concentrations of each primer; thermal cycling conditions were 95°C for 4 min followed by 50 cycles of 30 s at 95°C, 30 s at 56°C, and 60 s at 72°C.

All qPCR assays used the SsoAdvanced SYBR Green qPCR Supermix (BioRad, Hercules, CA, USA) and a BioRad CFX 96 Real-Time System. Data were analysed using Bio-Rad CFX Manager (Version 2.1), and results were reported as the number of gene copies per g dry weight of soil after averaging three technical replicates per sample.

2.7. Statistical Analyses

Shapiro–Wilks tests confirmed that data for all soil parameters and plant biomass metrics were normally distributed. Log(10) transformations were necessary for 16S rRNA and nirS gene abundances prior to statistical analyses. To investigate main effects of site (Rice vs. James) and plant type (Juncus, Typha, and Peltandra), data were examined using a two-factor analysis of variance (ANOVA) applied to each sampling date (4 replicates for reference soil and each plant type in November, 5 of each in June). Whenever a significant interaction effect was obtained, separate single-factor ANOVAs were run to assess plant differences within sites, and t-tests were used to compare a given plant species across sites. Tukey’s HSD was used for all post hoc comparisons.

To quantify the relative influence that each plant type had on gene abundance, a “plant-effect ratio” was calculated as the abundance in the rhizosphere (or rhizoplane) divided by abundance in the unvegetated reference soil. The null hypothesis of no plant effect predicts a ratio of one, wherein the proportion of the community that is bacteria (or nirS denitrifiers) is the same in the rhizosphere as in the bulk soil. Ratios greater than one represent enhanced abundance of the microbial group proximal to the plant roots, and ratios less than one indicate suppression. One-sample t-tests were used to determine when the plant-effect ratios were significantly different from one.

Lastly, stepwise multiple regression ( entry = 0.10, p exclusion = 0.15, for final inclusion = 0.05) was used to quantify the relationship between soil properties (pH, redox, moisture, OM, and C : N) and microbial abundance in the rhizosphere. Because the extent and mechanism of any plant effects were expected to differ depending on season (i.e., November senescence vs. June growth), data for each sampling event were analysed separately. Equivalent analysis of rhizoplane abundance could not be conducted because comparable environmental data were not available. All statistical analyses were performed using SPSS (Version 21, IBM Corp., Armonk, NY, USA).

3. Results

3.1. Microbial Abundance in the Rhizosphere

Bacterial abundance in the rhizosphere ranged from 107 to 10816S rRNA gene copies per gDW of soil (Table 1). In November, no statistically significant differences due to either site or plant type were detected. In contrast, differences due to both factors were significant for June (two-factor ANOVA, all ). We found bacteria abundance at that time to be lowest in the unvegetated soil (mean across both sites: 1.3 (±0.4) × 107), 2- to 3-fold higher for Peltandra and Typha (3.7 (±0.7) × 107) and 7-fold higher for Juncus (8.4 (±2.0) × 107). Overall, bacterial abundances were greater at James than Rice by a factor of ∼1.5–2.

The abundance of denitrifiers in the rhizosphere was lower than that for bacteria, ranging from 106 to 107nirS gene copies per gDW of soil (Table 1). For both sampling events, a significant effect of the plant type was observed (all ). Unvegetated soils consistently had the lowest abundance of nirS denitrifiers, followed by soil samples from the rhizosphere of Juncus and Peltandra (∼5-fold greater), and then soil samples from Typha (∼10-fold greater overall). There was no site effect in November; however, in June, a significant site effect was observed (), wherein abundance was again higher at James than Rice.

3.2. Plant-Effect Ratios

Plant-effect ratios for rhizosphere samples from November (Figure 2(a)) indicate no consistent pattern of enhancement of bacteria abundance relative to the unvegetated soils (only the t-test for James-Peltandra was significant ()). In contrast, plants sampled in June had significantly higher bacterial abundance (Figure 2(a)) in their rhizospheres compared to unvegetated soils by a factor of ∼2 (all , except Rice-Typha ). A two-factor ANOVA showed no differences in plant-effect ratio associated with the site or plant type.

For denitrifiers (Figure 2(b)), plant-effect ratios were significantly greater than 1 for Peltandra and Typha for both sites and sampling events (all ). In November, denitrifier abundance was 2 to 3 times greater in the rhizosphere of these two plants compared to unvegetated soil, whereas the enrichment was 4- to 8-fold in June. For both sampling dates, the enrichment effect was significantly greater at Rice compared to James (site effect from two-factor ANOVA; November , June ).

Effect ratios were also calculated for rhizoplane samples in June. For bacteria, the rhizoplane effect ratio was always greater than the rhizosphere effect ratio; this trend is evidenced by all points being above the 1:1 line in Figure 3(a). In the case of denitrifier abundance, rhizosphere and rhizoplane effect ratios were similar (along the 1:1 line in Figure 3(b)).

3.3. Plant Characteristics

Plant characteristics were assessed during the growing season (June), and significant differences were observed due to plant type but not site (Table 2). Specifically, the C : N content of the aboveground vegetation of Juncus (45.5 ± 7.4, combining both sites) was higher than that of Peltandra (22.0 ± 2.0) and Typha (25.8 ± 2.0). Plant species also differed with respect to aboveground plant biomass (gDW per individual) wherein Peltandra (53 ± 11) and Typha (71 ± 5) had significantly less dry mass than Juncus (135 ± 35) across both sites. Similarly, belowground biomass (gDW per individual) differed across species but not sites. Typha had significantly less root mass (36 ± 6) than Peltandra (282 ± 47) and Juncus (237 ± 38). The average belowground biomass for unvegetated reference soil was 0.6 ± 0.2 mgDW root material per gDW of soil and did not differ across sites.

3.4. Soil Properties

For both sampling events, soil properties varied between sites and plant types (Table 1). In general, the lowest soil moisture values were associated with Juncus rhizosphere samples, and the wettest conditions were associated with Typha. The sole exception to this was at James in November, when soil moisture was uniformly high for all samples (averaging 81 ± 1%). Soils at James were generally wetter than soils at Rice, and both sites had higher moisture in November (74 ± 2) than June (57 ± 2).

Soil redox potential (mV) was typically higher for Juncus (group means ranged from 210 to 321, excluding Nov-James) except in November at the James site where redox values, like soil moisture values, did not differ significantly. Redox potential was consistently higher at Rice (68 ± 26) compared to James (23 ± 16), though the sites were not significantly different for either month.

In November, there were no differences in pH across plant species, but pH was slightly higher at Rice (overall mean and standard error: 7.2 ± 0.2) compared to James (6.0 ± 0.0). In June, the pattern was reversed. pH was consistently higher at James, and there were also differences in pH across the plant species. Means for the unvegetated soil were ∼1 unit lower than the rhizosphere soils.

Irrespective of sampling date and plant type, the overall mean for OM (%) was higher at James (20 ± 1) than Rice (14 ± 1). Values were ∼5% higher during plant senescence in November for both sites. Values for Juncus were lowest regardless of the site or date. Organic matter quality as assessed by C : N was fairly consistent throughout the experiment (14.7 ± 0.2). Values associated with Typha were always the lowest (group means ranged from 12.4 to 13.8) though not significantly different from the unvegetated soil (13.2–15.7). Similarly, values for Juncus were generally the highest (15.5–19.4), and this distinction was especially noticeable in June.

3.5. Multiple Regression Analysis

The environmental data were relatively poor predictors of bacteria and denitrifier abundance in November. For both datasets, a significant multiple regression model could be created using pH and C : N, but the amount of variance explained was low (df = 31; bacteria: adj-R2=0.17, ; denitrifiers: adj-R2=0.18, ). This is in sharp contrast to June, when the environmental data were excellent predictors of microbial abundance. For both bacteria and denitrifier abundance, the best model included a positive relationship with pH and OM and a negative correlation with soil moisture. The adjusted R2 values for this data set were considerably higher: 0.81 and 0.86 for bacteria and denitrifiers, respectively (both ).

4. Discussion

The influence of plant species on rhizosphere soil conditions has been well documented (e.g., Lamers et al. [42], Grayston et al. [11], and Bodelier et al. [43]); however, corresponding effects on the soil microbial community remain understudied. In this research, we examined two microbial groups, total bacteria (via 16S rRNA gene) and denitrifiers (via nirS), and found them to be variably distributed across the rhizosphere and rhizoplane microhabitats of three wetland macrophyte species. We found the effect to differ based on site, which is consistent with prior research that there is a synergistic relationship between plant species and soil type that shapes microbial abundance in the rhizosphere [5, 4446]. The physical-chemical properties of the native soil can be important directly via their impact on the microbial community or indirectly via the selective pressure they exert that determines which plants are successful in certain areas [47, 48].

4.1. Plant Species Effects

For both sites and dates, the abundance of bacteria was lowest in the unvegetated soil. Abundance was usually higher in the rhizosphere samples, but no consistent differences were observed across the macrophyte species. Plant-effect ratios indicate little-to-no rhizosphere enhancement during November and only a modest effect (∼1.5-fold increase) in June (Figure 2(a)). Moreover, the plant-effect ratio did not differ significantly based on the plant type. In contrast, enrichment for denitrifiers was more pronounced and varied across plant species. The rhizosphere associated with Typha and Peltandra had a significantly greater abundance of nirS denitrifiers compared to Juncus. This pattern followed expectations of ROL associated with each plant. In addition, the finding of an especially strong rhizosphere effect for Typha suggests a greater change in microenvironmental conditions occurs near the roots of this plant, relative to the other plant types studied. In fact, Typha was usually found in the wetter areas of the marsh and has the greatest expected levels of ROL, which would set up a strong gradient in environmental conditions associated with the rhizosphere.

The plant-soil continuum is a complex and heterogeneous environment, and high spatial variation shapes the microbial community that exists within the matrix. Typically, the rhizoplane is thought to be the epicenter of soil microbial activity due to high nutrient availability [49, 50]. In this study, the bacterial abundance in the rhizoplane was greater when compared to the rhizosphere or bulk soil (Figure 3(a)). In contrast, nirS denitrifiers did not show increased abundance in the rhizoplane.

4.2. Season and Site Effects

In general, microbial abundance estimates were higher in November than in June. This increase is coincident with increased litter availability [51] and root senescence, which could both provide a source of fresh OM to fuel microbial growth. Our results are consistent with previous findings examining annual dynamics of denitrification and soil carbon pools [5153]. While greater abundances were observed in November, plant-effect ratios were greatest in June. For example, we saw few rhizosphere plant effects for bacteria in November, but a significant increase in abundance for all three plant species at both sites in June (i.e., five out of six plant-effect ratios were significantly greater than one) (Figure 2). For denitrifiers, plant effects were evident for both sampling events, but they were much larger in June. The increased importance of vegetation during June may be due, at least in part, to the increase in plant activity during the growing season. During this time, we expect higher rates of photosynthesis, rhizodeposition, and ROL. These stimuli effectively increase plant effects in the rhizosphere compared to bulk soils.

The wetland plant species considered in this study varied in terms of aboveground and root biomass and foliar C : N content, but these parameters did not change across sites. However, the sites did differ with respect to soil physical-chemical properties. Specifically, soil moisture was higher at James, as was OM content (Table 1). For other properties (redox, pH, and C : N), within-site variation due to plant type was greater than differences across sites. The difference in OM observed between the two study sites most likely has to do with the site history. In particular, clear-cutting of the forested wetland at Rice prior to dam construction in 1920 removed a great deal of organic carbon, and the accumulation of sediments and OM decreased during impoundment. In contrast to Rice, the site at James is relatively undisturbed, although the surrounding areas were logged during the 1970s–1980s.

With regards to microbiology, differences across sites were apparent, especially in June. Both total bacterial abundance and nirS denitrifier abundance were higher at James than Rice. Moreover, we observed a difference in nirS plant effect ratios based on the site (Figure 2(b)) and documented a greater effect of plants at Rice. One possible explanation for this relates back to OM availability. In wetlands with lower OM reserves, such as Rice, plant exudates may be a more important source of carbon and energy to the bacterial community; thus, macrophyte effects may be greater. Microorganisms in the soils at James, with its consistently higher OM, may be less reliant on these root exudates.

4.3. Relationship with Soil Properties

Multiple regression analysis was applied to determine which attributes of the soil environment were good predictors of microbial abundance for each sampling event. In June, a model including soil moisture, and pH, OM content allowed us to explain 81% and 86% of the variation in bacterial and nirS denitrifier abundance, respectively. This finding is consistent with our assertion that soil OM is more limiting during June and that root exudates are an important component of the resource pool at this time. The inclusion of pH in the model could be interpreted as a reflection of root activity. Root exudates can alter soil pH by releasing low molecular weight organic acids and alcohols [54] or by releasing hydrogen or bicarbonate ions [55]. The latter is more likely to have occurred in the soils examined in this study, as the rhizosphere pH tended to be greater than the reference soil pH. Regardless of what is causing the pH fluctuations, our identification of pH as a key driver of the microbial community is consistent with several recent studies across diverse ecosystems [40, 56, 57]. When the data for November are considered, a different suite of predictors are important. In particular, we see C : N emerge in the model for both the bacterial and nirS denitrifying community. This makes sense if we hypothesize that the quality of the native soil OM (i.e., C : N) will be more important at this time of year when root activity is less. This effect could be especially important as plant litter accumulates in the soil because the C : N content of the vegetation differed by plant species (Table 2). Interestingly, the variance explained by the multiple regression models was much lower for November (17 and 18% for bacteria and nirS denitrifiers, respectively) compared to June and indicates other environmental factors not considered in this study may be more important drivers of abundance. For example, if abundance at this time is largely affected by senescing resource inputs from plants, as hypothesized above, variations in the timing and intensity of these inputs across plant species may be an important factor that is not well captured by our sampling design. These pulses of inputs have been shown to be the major determinant of seasonal variation of soil microbial biomass [58], with different plant types and species demonstrating different responses to temporal variation.

5. Conclusions

Wetland plant species were found to have varying abundances of nirS denitrifying bacteria in their rhizospheres. Unvegetated reference soils always had the lowest abundances followed by Juncus and Peltandra, with Typha consistently having the highest abundances. Peltandra’s and Typha’s higher abundance of nirS denitrifying bacteria could be a result of a combination of soil conditions, ROL associated with the various root morphologies, and differences in root exudation. With regards to the total bacterial community, only a modest rhizosphere effect was detected, and it did not differ across plant species. For both microbial groups, plant effects were greater in June compared to November, most likely due to higher activity associated with the growing season. These stimuli effectively increase plant effects in the rhizosphere compared to bulk soils. Analysis of site effects revealed a stronger influence of plants at the younger, recently restored wetland; we hypothesize that this may be due to site history and/or differences in the standing stock of soil organic matter. Collectively, our results show that microbial abundance in the rhizosphere may be determined by a synergistic relationship between plant species and soil conditions. These findings have implications for our understanding of plant-microbe relationships in saturated soils and contribute to a growing appreciation of how this relationship may impact ecosystem biogeochemistry. Moreover, the work may be of interest to applied scientists interested in engineering wetland systems, as our results suggest that the effects of plant community composition can be mediated through microbial responses.

Data Availability

All of the original data related to the current study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This research was funded in part by the Undergraduate Research Opportunities Program (UROP) Summer Fellowship from Virginia Commonwealth University as well as the VCU Rice Rivers Center Student Research Grant, both to Morina. Thanks to Jamie Gillespie, David Berrier, Rana Mehr, Aaron Porter, and Nicki Dadashian for their assistance in the field and in the lab. Thanks to Chansotheary Dang for her edits in this manuscript. We would like to acknowledge Cyrus Brame with the James River National Wildlife Refuge for assistance with site access. This study was supported in part by the VCU Rice Rivers Center (Publication # 86).