# Plant Nutrient Resource Use Strategies Shape Active Rhizosphere Microbiota Through Root Exudation

Julien Guyonnet, Martin Guillemet, Audrey Dubost, Laurent Simon, Philippe

Ortet, mohamed Barakat, Thierry Heulin, Wafa Achouak, Feth El Zahar

Haichar

► To cite this version:

Julien Guyonnet, Martin Guillemet, Audrey Dubost, Laurent Simon, Philippe Ortet, et al. Plant Nutrient Resource Use Strategies Shape Active Rhizosphere Microbiota Through Root Exudation. Frontiers in Plant Science, 2018, 9, pp.1-13. 10.3389/fpls.2018.01662. cea-01968948

HAL Id: cea-01968948

https://cea.hal.science/cea-01968948

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Plant Nutrient Resource Use Strategies Shape Active Rhizosphere Microbiota Through Root Exudation

Julien P. Guyonnet¹, Martin Guillemet¹,², Audrey Dubost¹, Laurent Simon³, Philippe Ortet⁴,⁵, Mohamed Barakat⁴,⁵, Thierry Heulin⁴,⁵, Wafa Achouak⁴,⁵ and Feth el Zahar Haichar¹*

¹ Laboratoire d'Ecologie Microbienne, UMR CNRS 5557, Univ Lyon, Université Claude Bernard Lyon 1, UMR INRA 1418, Villeurbanne, France, ² Master de Biologie, École Normale Supérieure de Lyon, Université Claude Bernard Lyon 1, Université de Lyon, Lyon, France, ³ CNRS, UMR 5023 LEHNA, Univ Lyon, Université Claude Bernard Lyon 1, Université Lyon 1, ENTPE, Villeurbanne, France, ⁴ CNRS, Laboratory for Microbial Ecology of the Rhizosphere and Extreme Environment, UMR 7265 BIAM, CEA, Aix Marseille Univ, Saint-Paul-lès-Durance, France, ⁵ CNRS, FR3098 ECCOREV, Aix Marseille Univ, Aix-en-Provence, France

Plant strategies for soil nutrient uptake have the potential to strongly influence plant-microbiota interactions, due to the competition between plants and microorganisms for soil nutrient acquisition and/or conservation. In the present study, we investigate whether these plant strategies could influence rhizosphere microbial activities via root exudation, and contribute to the microbiota diversification of active bacterial communities colonizing the root-adhering soil (RAS) and inhabiting the root tissues. We applied a DNA-based stable isotope probing (DNA-SIP) approach to six grass species distributed along a gradient of plant nutrient resource strategies, from conservative species, characterized by low nitrogen (N) uptake, a long lifespans and low root exudation level, to exploitative species, characterized by high rates of photosynthesis, rapid rates of N uptake and high root exudation level. We analyzed their (i) associated microbiota composition involved in root exudate assimilation and soil organic matter (SOM) degradation by 16S-rRNA-based metabarcoding. (ii) We determine the impact of root exudation level on microbial activities (denitrification and respiration) by gas chromatography. Measurement of microbial activities revealed an increase in denitrification and respiration activities for microbial communities colonizing the RAS of exploitative species. This increase of microbial activities results probably from a higher exudation rate and more diverse metabolites by exploitative plant species. Furthermore, our results demonstrate that plant nutrient resource strategies have a role in shaping active microbiota. We present evidence demonstrating that plant nutrient use strategies shape active microbiota involved in root exudate assimilation and SOM degradation via root exudation.

**Keywords:** active bacterial community, denitrification, microbiota, microbial activities, plant nutrient use strategies, rhizosphere, root exudation, stable isotope probing (SIP)

## INTRODUCTION

Several studies have shown that plants can host up to 1000 distinct microbial species in their different organs (Bulgarelli et al., 2013; Hacquard et al., 2015). Plant-associated microbiota benefit from their host, since plants provide easily degradable carbon *via* root exudation. Indeed, up to 20% of the carbon fixed by photosynthesis is released in the rhizosphere, which constitutes nutrient

**OPEN ACCESS**

**Edited by:**

Essaid Ait Barka,
Université de Reims Champagne
Ardenne, France

**Reviewed by:**

Collin M. Timm,
Johns Hopkins University,
United States
Sameh Selim,
Institut Polytechnique LaSalle
Beauvais, France

*Correspondence:

Feth el Zahar Haichar
zahar.haichar@univ-lyon1.fr

**Specialty section:**

This article was submitted to
Plant Microbe Interactions,
a section of the journal
Frontiers in Plant Science

**Received:** 09 July 2018

**Accepted:** 26 October 2018

**Published:** 23 November 2018

**Citation:**

Guyonnet JP, Guillemet M, Dubost A, Simon L, Ortet P, Barakat M, Heulin T, Achouak W and Haichar FZ (2018) Plant Nutrient Resource Use Strategies Shape Active Rhizosphere Microbiota Through Root Exudation. *Front. Plant Sci.* 9:1662. doi: 10.3389/fpls.2018.01662
sources and signaling molecules for plant-microbe interactions (Bais et al., 2006; Haichar et al., 2014). In turn, plant-associated microbiota influence plant growth by supplying nutrients such as nitrogen (N), phosphorus, or potassium (Singh A. et al., 2004), as well as protecting plants from phytopathogens (Mendes et al., 2011; Samain et al., 2017).

Plants adopt different strategies to interact with their environment for the acquisition or conservation of nutrients (Lambers et al., 2008). These nutrient use strategies could be explained by plant functional traits (McGill et al., 2006; Violle et al., 2007). Indeed, plants with higher photosynthetic capacity and efficient N uptake are referred to as exploitative (or fast-growing) plant species (Gross et al., 2007). By contrast, plants with lower nutrient uptake and lower photosynthetic activities but higher levels of leaf and root dry mass are referred to as conservative (or slow-growing) plant species (Tilman, 1990; Ryser and Lambers, 1995; Aerts and Chapin, 1999). These contrasted strategies of nutrient management influence the input and output of C-resources, from the assimilation of carbon by photosynthesis to its release in the rhizosphere soil *via* root exudation (Guyonnet et al., 2018).

The influence of root exudation of C-resources on microbial communities inhabiting the RAS and colonizing the root systems has been revealed using a stable isotope probing (SIP) approach (Haichar et al., 2008, 2012). These studies displayed significant differences in the active communities in the rhizosphere of four plant species. Although previous studies have explored the link between plant functional traits and plant performance, less is known about the impact of plant nutrient use strategies on bacterial populations that assimilate root exudates, and how they contribute to soil organic matter (SOM) degradation by rhizosphere bacteria.

The aim of this work was to determine the impact of plant nutrient resource strategies on microbial activities of carbon (i.e., respiration) and nitrogen (i.e., denitrification) cycles in soils, as well as on the diversity of active microbiota through root exudation. Specifically, we hypothesized that (i) exploitative plants by exhibiting higher exudation rates than conservative plants will induce an increase in microbial activity (particularly denitrification activity). (ii) Exploitative plants might select a specific bacterial community through root exudation as compared to conservative plants. In addition, by exuding more carbon, exploitative species could promote SOM degradation *via* the priming effect.

To investigate this issue, six plant species in the family Poaceae were selected along a nutrient resource gradient, from exploitative to conservative species. The plants were cultivated in the same soil for 10 weeks in a growth chamber and were then submitted to ¹³CO₂ labeling for 1 week in order to applied SIP. This approach was employed to distinguish the metabolically active bacterial populations that assimilate root exudates from those growing on SOM. In addition, microbial activities (i.e., respiration and denitrification) were measured for each plant rhizosphere.

# MATERIALS AND METHODS

## Plant Growth

Six perennial C3 grass species distributed along a gradient of plant nutrient use strategies (Gross et al., 2009; Maire et al., 2009; Cantarel et al., 2015; Guyonnet et al., 2018) were examined: two conservative species, *Festuca paniculata* (FP) and *Sesleria caerulea* (SC); and four exploitative species, *Bromus erectus* (BE), *Anthoxanthum odoratum*, (AO), *Dactylis glomerata* (DG), and *Trisetum flavescens* (TF). All studied grasses belong to the *Pooideae* subfamily (Bouchenak-Khelladi et al., 2008). All species were sampled in the field (French Alps) and separated into individual tillers according to Guyonnet et al. (2017). Plants were cultivated on luvisol with no added nitrogen source and were collected at La Côte Saint-André (Isère, France), which is continuously cropped with maize (Guyonnet et al., 2017). The soil was sieved (2 mm mesh size) and 170 g were allocated to each plastic pot. The three individual plants from each species were grown in a greenhouse (13 h day at 22°C/11 h night at 18°C), with an approximate light intensity of 8-10 klux for 11 weeks. Soil moisture was manually controlled. Three pots containing bulk soil (BS) without any plants were incubated as controls under the same conditions.

## Plant ¹³C-Labeling

Continuous labeling was initiated 10 weeks after plant growth, according to a previously reported protocol (Haichar et al., 2008). Plants were placed in a growth chamber ("PHYTOTEC" facility, CEA Cadarache, France) equipped for automatic control of light, temperature, moisture, evapotranspiration, irrigation, and CO₂ concentration. The day-night period was set at 8 h/16 h, respectively; light intensity was 13.5 klux; maximum daily temperatures ranged from 20 to 22°C; air moisture was adjusted to 80%; and CO₂ concentration was maintained at 350 µl l⁻¹. CO₂ partial pressure was kept constant by injection of pure ¹³CO₂ (>99% atom ¹³C purchased from Cortec Net, Paris, France) during active photosynthesis. The isotope excess in the chamber was maintained at >95% atom ¹³C during the 7 days of labeling. Plants and BS were collected after 1 week of ¹³CO₂ labeling.

## Plant Harvesting

At the end of the labeling experiment, three plants from each population were used to measure microbial activities and to analyze active bacterial diversity. The roots of each plant were manually separated from the root-adhering soil (RAS). Roots were washed carefully with distilled water to remove any remaining soil particles and then immediately frozen in liquid N₂ and stored at -80°C for DNA extraction. RAS was then carefully separated from the rest of the fine roots and stored at 4°C for microbial activities measurements. Few grams of the RAS were lyophilized, and stored at -80°C for DNA extraction.

## Microbial Activities

Substrate-induced respiration (SIR) and the denitrifying enzymatic assay (DEA) were measured in triplicates according
to Bardon et al. (2014) using fresh soil samples (10 g equivalent dry weight of RAS from each plant and BS). SIR was determined as the linear rate of production of CO₂ during short-term (6 h) incubation under aerobic conditions. Glucose (1.2 mg C g⁻¹ dry soil) was added to the soil samples and the soil moisture was brought to 70% water holding capacity. DEA was determined as the linear rate of production of N₂O during short-term (6 h) incubation under anaerobic conditions. Acetylene gas (C₂H₂) was used to inhibit nitrous oxide reductase activity and avoid N₂ production. Glucose (0.5 mg C g⁻¹ dry soil), glutamic acid (0.5 mg C g⁻¹ dry soil) and KNO₃ (50 µg N-NO₃⁻ g⁻¹ dry soil) were added to the soil samples and the soil moisture was brought to 100% water holding capacity. Gasses (CO₂ and N₂O) were measured using gas chromatography (Micro GCR3000, SRA Instrument; Marcy-l'Étoile, France).

## DNA Extraction and Gradient Fractionation

DNA was extracted in triplicates from either 3 g of RAS using the “FastDNA™ SPIN Kit for Soil” (MP Biomedicals; United States) or from 0.5 g of fresh roots as described by Ranjard et al. (2003). For extraction from roots, DNA was derived from endophytic microorganisms as well as microorganisms firmly attached to the root surface.

Extracted DNA was quantified using a Qubit® 2.0 Fluorometer (Life Technologies; United States). Extracted DNA from RAS was fractionated by CsCl equilibrium density-centrifugation according to Haichar et al. (2008, 2012), using 4.9 mL “Optiseal Tubes” (Beckman Coulter; United States) in a VTi 65.2 rotor (Beckman Coulter) at 20°C and 45,000 rpm, for 42 h. Buoyant density was obtained by weighing, and DNA was quantified in the CsCl fraction using the Picogreen assay (Molecular Probes). DNA was purified from CsCl salts using the Geneclean turbo kit (Qbiogene; Montreal, QC, Canada).

## Isotope Analysis

The ¹³C content of DNA (total, heavy- and light-DNA fractions) was measured by IRMS (Isoprime 100, Isoprime, Ltd.) coupled with an elemental analyzer (Thermo FlashEA 1112, ThermoElectron) as described in Haichar et al. (2007). 5 µl of DNA was transferred to 3.5-mm × 5-mm “Ultra Clean” tin capsules (Elemental Microanalysis; United Kingdom) and dried for 2 h at 60°C. Tin capsules were then submitted to mass spectrometry.

## Sequencing and Bioinformatics Analysis

Five µl samples (~50 ng) from each plant's DNA in triplicates were sent for sequencing to FASTERIS (Switzerland) using the MiSeq illumina technology, comprising fractions representative of light-, heavy-, and root-DNA in triplicates for each plant. The V3–V4 domain of the 16S rRNA gene was amplified with the primers 16S rRNA Fwd primer 3'-CCTACGGNGGCWGCAG-5' and 16S rRNA Rev primer 3'-GACTACHVGGGTATCTAATCC-5'. The tagged primers and barcodes were trimmed off by FASTERIS.

We analyzed data using the FROGS software (Escudie et al., 2015). Corrected 16S rRNA gene sequences were clustered into operational taxonomic units (OTUs) using the swarm method (Mahé et al., 2014), in which the distance aggregate was set at *d* = 3. Potential chimeric formations were removed using the VSEARCH tool. OTUs with abundance greater than 0.005% were retained. The taxonomic identification of bacterial and archaeal 16S rRNA genes was performed using the SILVA database release 123 with RDP Classifier (Wang et al., 2007).

## Network Inference and Graph Display

To study bacteria inhabiting the RAS, OTUs from the RAS found in more than 25% of the samples were selected; the same procedure was performed for root-colonizing bacteria. Microbial networks were constructed using the R package phyloseq (version1.19.1). The pairwise co-occurrence of OTUs over the different samples was assessed using a Bray–Curtis distance (Williams et al., 2014). A proportion table was used instead of an abundance table, so that only the patterns of repartition would be compared and the absolute abundance of any one given OTU would not have an effect. A threshold of 0.12 was used to select significant associations. Every node, and thus every OTU, was positioned on a graph using Cytoscape (version 3.5.1). The horizontal position of every node was calculated according to the preference of the bacteria toward exploitative or conservative plants as follows: $x = \frac{\%exploitative - \%conservative}{\%exploitative + \%conservative}$ where %exploitative and %conservative represent the proportion of reads detected in the exploitative and conservative samples, respectively. Additionally, the vertical position of RAS-inhabiting bacteria conveys the inclination toward either the light (L) or heavy (H) compartment: $\frac{\%H - \%L}{\%H + \%L}$ where %H and %L are the proportion of reads detected in the light- and heavy-DNA fractions, respectively.

## Spider Diagrams

Spider diagrams were plotted in order to present the relative specificity of bacterial phyla toward their hosts. Each community was tested, including bacteria that feed on SOM (light-DNA), bacteria that feed on root exudates (heavy-DNA), and bacteria that colonize root tissues (root-DNA). In each compartment, the total number of OTUs in each phylum was calculated for all the plants and divided by the mean abundance of this particular phylum in all the plants studied. The logarithms of these ratios are displayed.

## Statistical Analyses

For each plant, one-way analysis of variance (ANOVA) and *post hoc* Tukey HSD tests were performed to examine differences in microbial activities between plant strategies. When necessary, data were log- or square-transformed prior to analysis to ensure conformity with the assumptions of normality and homogeneity of variances. When it was not possible (i.e., for exudation variables), non-parametric Kruskal–Wallis and Wilcoxon tests were used. The correlation between microbial activities and root exudation activity was tested by building a linear model with log-transformed microbial
activity data. Effects with *p*-value < 0.05 are referred to as significant.

A PCoA (principal coordinates analysis) based on a Bray-Curtis distance was used to compare bacterial diversity retrieved from heavy-, light-, and root-DNA fractions from each plant rhizosphere. For each plant and DNA fraction, the triplicates for OTUs abundance were regrouped. Then, the table of OTU abundance for each fraction from the studied plants and BS samples was treated with the phyloseq package in R, according to McMurdie and Holmes (2013).

The core microbiota were identified using QIIME (Caporaso et al., 2010) and were determined by plotting OTU abundance in the core at 5% intervals (from 50 to 100% of samples). To remove sample heterogeneity, we performed rarefaction on the number of sequences. The approach of rarefaction is to randomly sample the same number of sequences from each sample, and use this data to compare the communities at a given level of sampling effort. We defined the core microbiota of each plant as the OTUs present in 100% of the samples. To determine core microbiota, all samples from each plant were compared across all compartments (light-, heavy-, and root-DNA fractions). Any taxa that were ubiquitous across all samples were then defined as part of the core microbiota of the compartment. Venn diagrams were then constructed from these data, using ‘VennDiagram’ from the R package to show common and unique OTUs within the conservative (FP, SC), the exploitative (AO, BE, DG, TF) and conservative versus exploitative plant groups. All statistical analyses were performed using R 3.1.2 (R Core Team, 2015).

## Data Deposition

The sequence data was deposited at EMBL-ENA public database (http://www.ebi.ac.uk/ena/data/view/PRJEB25083).

# RESULTS

## Impact of Root Exudation on Microbial Activity

In order to determine the impact of root exudation levels of each plant species, measured in our previous work by Guyonnet et al. (2018), on microbial activities, respiratory activities are plotted against root exudation level (Figure 1). Respiratory activities under aerobic and anaerobic conditions that were measured in triplicate on the RAS increased in response to enhanced root exudation activity, as reported in Figure 1. SIR and DEA were log-correlated with root exudation activity (Pearson correlation for DEA and SIR = 0.6; *p* < 0.05). SIR measured on the RAS fraction of SC, FP, AO, BE, and DG plants did not significantly differ from the unique BS sample (Figure 1A). However, SIR was significantly higher for the RAS of TF as compared to SC (Tukey test, *p* < 0.05). Potential denitrification in RAS was significantly higher for AO, BE, DG, and TF than for BS (Tukey test, *p* < 0.05) (Figure 1B). The SIR and DEA activities in the RAS fraction of the four exploitative species (TF, AO, BE, DG) were significantly higher than the RAS fractions of the two conservative species (FP, SC), which did not differ from the BS.

## ¹³C Enrichment of Microbial Communities

To confirm the assimilation of ¹³C-root exudates by soil microbiota, we measured the ¹³C/¹²C ratio of DNA extracted from the RAS recovered from each plant rhizosphere. The δ¹³C of DNA extracted from the RAS of non-labeled plants was approximately -28.7‰ on average, corresponding to natural ¹³C abundance levels. Furthermore, δ¹³C values were, on average,

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

**FIGURE 1** | Correlation between (A) root exudation rate and SIR and (B) root exudation rate and DEA in the rhizosphere of six grasses, distributed on a gradient of plant resource use strategies. Specimens include *Sesleria caerulea* (SC), *Festuca paniculata* (FP), *Anthoxanthum odoratum* (AO), *Berectu erectus* (BE), *Trisetum flavescens* (TF) and *Dactylis glomerata* (DG), as well as BS (control, in the absence of plants). Data represent the means, and error bars represent SE. Significance level is at *p* < 0.05 between plants and BS.
approximately -24.3‰ for BS, -22.2 and -19.2‰ for SC and FP respectively (conservative species), -12.5‰ for BE, 0.01‰ for AO, -5.44‰ for DG, and 13.9‰ for TF (exploitative species).

DNA extracted from each RAS sample was fractionated by ultracentrifugation, and the $\delta^{13}$C of the fractions along the gradient were determined. An example of the partitioning of DNA content and $^{13}$C-labeling among gradient fractions of increasing density for DNA extracted from the RAS of DG plants is presented in **Figure 2**. The DNA density profile shows a unique peak ranging from 1.62 to 1.63, with $^{13}$C-labeling measured by IRMS reaching a maximum density between 1.64 and 1.67. Similar profiles were obtained for DNA extracted from the RAS of all other studied grasses. One fraction with a buoyant density of 1.62 g ml$^{-1}$ CsCl was selected as representative of the light-DNA fraction, representing bacterial populations mainly involved in SOM degradation. The fraction with a buoyant density of 1.66 g ml$^{-1}$ CsCl was selected as representative of the heavy-DNA from the bacterial community actively assimilating $^{13}$C-root exudates (Figure 2).

## Diversity of Active Microbiota Colonizing the Plant Rhizosphere

To investigate the impact of root exudates produced by each plant species on the diversity of active microbiota, three fractions

from each plant: the heavy-DNA fraction, the light-DNA fraction, and root-DNA were analyzed by 16S rRNA gene sequencing. Sequencing resulted in 164,640 16S rRNA bacterial sequences corresponding to 1,421 Operational Taxonomic Units per sample, clustered using the swarm method with *d* = 3. The rarefaction curves (rarefaction depth about 2744 sequences), which displayed the observed OTU richness as a function of the sequencing effort, indicate that the sequencing depth was nearly reached to completely comprehend the diversity present in BS, RAS and root tissues of all the plants (Supplementary Figure S1). Alpha-diversity and Inverted Simpson indices were not significantly different between BS and RAS for each plant, although they were significantly decreased in root compartment of the six species (Supplementary Table S1).

To assess structural similarities and/or differences in active microbiota involved in SOM degradation (light-DNA fractions) and root exudate assimilation (heavy-DNA fractions) from the RAS and root tissues (root-DNA fractions) among the plant species, a PCoA plot was generated based on Bray-Curtis distances (Supplementary Figure S2). Bray-Curtis distances provide a measure of community composition differences between samples based on OTU counts, regardless of the taxonomic assignment. Ordinations based on this metric demonstrated a strong separation (Axis 1 = 60.9%) of the bacterial community colonizing plants root tissues (root-DNA

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Incorporation of $^{13}$C-root exudates into microbial community DNA derived from the RAS of *Dactylis glomerata* plants after 1 week of $^{13}$CO$_2$ labeling. $^{12}$C (light) and $^{13}$C (heavy) DNA were separated by CsCl density gradient centrifugation. Total DNA was quantified fluorometrically (■), and $^{13}$C ($\bar{p}^{13}$C) values were measured by IRMS within gradient fractions (1–10) (♦).</figcaption></figure>
fractions) from those inhabiting the RAS involved in SOM
degradation (light-DNA fractions) or in root exudate assimilation
(heavy-DNA fractions) (**Supplementary Figure S2A**). The
PCoA generated in **Supplementary Figure S2B** shows a clear
distinction between the microbiota retrieved from light-DNA
fractions and heavy-DNA fractions from the RAS with those
obtained from the BS, indicating the impact of plant *via*
root exudates on active microbiota diversity. In addition, by
analyzing the diversity of bacteria colonizing the RAS of all
the plants, we observed a clear separation and significant
differences (p < 0.001) between heavy- and light-DNA fractions
(triangles vs. circles; **Supplementary Figure S2B**). However,
the influence of the plant species on SOM-degrading bacterial
diversity retrieved from light-DNA fractions appears to be
marginal (p = 0.06) (**Figure 2**). The PCoA ordination indicates
a strong separation between active root-colonizing microbiota
(**Supplementary Figure S2C**, Axis 1 = 50%; Axis 2 = 18.4%) of
the six plant species, revealing the high selectivity exerted by root-
colonizing microbiota (**Supplementary Figure S2C**), but without
clustering of conservative vs. exploitative species.

Several phyla were present in different amounts in the three
studied DNA fractions, such as *Proteobacteria*, *Actinobacteria*,
*Bacteroidetes*, *Firmicutes*, and *Acidobacteria* (Figure 3). For
example, in the light-DNA fraction, the *Firmicutes* phylum was
more abundant in the rhizosphere of FP as compared to the other
plants, whereas the *Verrucomicrobia* phylum was more abundant

in the RAS fraction of DG (Figure 3A). For the heavy-DNA
fraction, members of the *Firmicutes* and *Bacteroidetes* phyla were
more abundant in FP and BE, respectively, in comparison to the
RAS fraction of other plant species (Figure 3B). In addition, the
*Nitrospirae* phylum was more abundant in the RAS fraction of
AO and SC than in the RAS of other plant species. For root-
DNA, the *Acidobacteria* phylum was more abundant on FP roots
as compared to other plant species roots (Figure 3C). We also
observed that the *Nitrospirae* was more abundant on BE and
DG roots, as compared to the other plants. At this level of
investigation, two plant species appear to strongly modify the
bacterial composition that is influenced directly or indirectly by
root exudation. Specifically, the *Firmicutes* are stimulated in the
RAS fraction (SOM degraders and root exudate consumers) and
the *Acidobacteria* are stimulated in the root tissue fraction of
*F. paniculata* (conservative species), whereas the *Bacteroidetes*
and *Nitrospirae* are stimulated in the rhizosphere (root exudate
consumers only) and the root tissue fraction of *B. erectus*
(exploitative species).

Core Microbiota of Exploitative Plants
Species

Operational taxonomic units distribution among all
compartments (root tissues, heavy- and light-DNA fractions)
were analyzed to highlight the existence of a common core

A

<figure><img src="image_6.png" /></figure>

Soil organic matter degraders (light-DNA)

B

<figure><img src="image_7.png" /></figure>

Root exudates consumers (heavy-DNA)

C

<figure><img src="image_8.png" /></figure>

Root exudates consumers (root-DNA)

<figure><img src="image_9.png" /></figure>

**FIGURE 3** | Relative abundance of five major bacterial phyla in the rhizosphere of two conservative plants (*Festuca paniculata* and *Sesleria caerulea*) and four exploitative plants (*Bromus erectus*, *Anthoxanthum odoratum*, *Dactylis glomerata*, and *Trisetum flavescens*). The relative abundances were calculated as the ratio of the abundance of a phylum in a given plant divided by the mean abundance of this phylum in all plants. This was performed for every compartment: (A) bacteria feeding on SOM (light-DNA), (B) root exudate consumers (heavy-DNA), and (C) bacteria colonizing the root tissues.
microbiota among the exploitative plants species (Figure 4). As shown in Figure 4A, 322 OTUs were shared by the four exploitative plant species corresponding to the core microbiota involved in SOM degradation. The Venn diagram also reveals unique OTUs corresponding to the microbiota involved in SOM degradation found in the rhizosphere of AO (370), DG (257), TF (413), and BE (931).

In the heavy-DNA fractions corresponding to root exudate consumers in the RAS, 218 OTUs were common to the four exploitative plant species. These shared OTUs can be regarded as the core microbiota involved in root exudate assimilation in the RAS among exploitative plant species. Unique OTUs corresponding to microbiota involved in root exudate assimilation were also detected with high numbers noted in AO (500) and BE (571) plants (Figure 4A). For the bacterial community colonizing the root tissues and involved in root exudate assimilation, 120 common OTUs were identified and could be considered as the root-associated core microbiota. The shared and unique OTUs in the rhizosphere of AO, DG, TF, and BE (light-, heavy-, and root-DNA fractions) are presented in Supplementary Tables S2–S4, respectively.

## Core Microbiota of Conservative Plants Species

As for the exploitative species, the conservative species (SC and FP) shared OTUs corresponding to the core microbiota involved in SOM degradation (236 OTUs), root exudate assimilation in the RAS (28), and in the root tissues (325) (Figure 4B). The Venn diagram also reveals that unique OTUs were observed, with 732 OTUs occupying 76% of the total OTU abundance (968) for SOM degraders in the FP rhizosphere. Furthermore, in the SC rhizosphere, 789 OTUs (96% of all OTUs) were unique to root exudate consumers inhabiting the RAS, and 950 OTUs (75% of all OTUs) were unique to root exudate consumers colonizing the root tissues (Figure 4B). The shared and unique OTUs in the SC and FP rhizosphere (light-, heavy-, and root-DNA fractions) are presented in Supplementary Tables S5–S7, respectively.

## Impact of Plant Nutrient Use Strategies on Diversity of Active Microbiota

The plant nutrient strategy (conservative vs. exploitative) influenced the composition of bacterial communities, as exploitative plants selected more specific SOM degraders (1,822 OTUs), exudate consumers in the RAS (2,691) and root-inhabiting bacteria (1,750) than conservative species (334, 168, and 643 OTUs, respectively) (Figure 4C). However, some bacterial OTUs are common to the conservative and exploitative plant species, including 921 OTUs for the SOM degraders, 686 OTUs for the root exudate consumers in the RAS, and 1,032 OTUs for the root colonizers (Supplementary Tables S8–S10).

## Co-occurrence Networks Among Plant Nutrient Use Strategies

Co-occurrence interactions between bacterial taxa were examined in order to determine which taxa respond according to plant strategy (conservative vs. exploitative), carbon nature

(fresh carbon/root exudates vs. organic matter) and ecological niche (RAS and root system), thereby providing a means to examine the environmental traits underlying the variation. The

<figure><img src="image_10.png" /></figure>

**FIGURE 4** | Venn diagram displaying the shared and unique bacterial OTUs at 100% identity among light-DNA (SOM degraders), heavy-DNA, and root-DNA fractions (root exudate consumers) retrieved from (A) exploitative species [*Bromus erectus* (BE), *Anthoxanthum odoratum* (AO), *Dactylis glomerata* (DG), and *Trisetum flavescens* (TF)], (B) conservative species [*Festuca paniculata* (FP) and *Sesleria caerulea* (SC)], and (C) exploitative and conservative species [*Bromus erectus* (BE), *Anthoxanthum odoratum* (AO), *Dactylis glomerata* (DG), *Trisetum flavescens* (TF), *Festuca paniculata* (FP), and *Sesleria caerulea* (SC)].
<figure><img src="image_11.png" /></figure>

**FIGURE 5** | Microbial networks of bacteria inhabiting (A) the RAS and involved in root exudate assimilation (heavy-DNA) or SOM degradation (light-DNA). Microbial networks of bacteria colonizing (B) the root tissues of conservative plants (Festuca paniculata and Sesleria caerulea) and exploitative plants (Bromus erectus, Anthoxanthum odoratum, Dactylis glomerata, and Trisetum flavescens). Each node represents a gene rRNA sequence, whereas the color indicates the phylum. Bacteria co-occurrence patterns, represented by gray lines, were assessed using a Bray–Curtis distance based on the repartition of OTUs in different plants and compartments. The resulting networks are positioned on a graph representing the gradient from conservative to exploitative plant species on the horizontal axis, and the gradient from SOM degraders (light-DNA) to root exudate consumers (heavy-DNA) on the vertical axis in graph (A). Overlapping nodes were scattered.

microbiota network of the bacterial association is summarized in **Figure 5**. Overall, the ecological networks were markedly different among exploitative and conservative species.

## Co-occurrence Network Analysis for RAS Microbiota

The resulting microbial network of the OTUs in RAS consisted of 444 nodes (OTUs) linked by 2,555 edges. The average number of neighbors for one node is 11.509, indicating that a given bacterial species will closely share its dispersal pattern with 11.5 other species in the outlined ecological niches (**Figure 5**). The clustering coefficient represents how the nodes are embedded in their neighborhood, and thus the extent to which they tend to cluster together. This clustering coefficient was 0.35, while the characteristic path length representing the characteristic distance between two random nodes was 3.3. Overall, this indicates that the different niches outlined were dense, with highly connected nodes that clustered with each other.

The structural properties (referred as the number of nodes and connections) of heavy-DNA fraction from the conservative species network were greater than those from exploitative species network, indicating a greater connection and closer relationships between microbial taxa inhabiting the RAS of conservative species (**Figure 5A** and **Supplementary Table S11**). Indeed, the number of OTUs belonging to the different phyla in the RAS of conservative species was higher than the exploitative

species (i.e., five *Actinobacteria* and five *Bacteroidetes* for conservative species vs. two *Actinobacteria* and two *Bacteroidetes* for exploitative species). Remarkably, *Acidobacteria*, *Chloroflexi* and *Firmicutes* were found exclusively in the microbial network of conservative species, whereas, *Verrucomicrobia* and *Armatimonadetes* were found only in exploitative species microbial network (**Supplementary Table S11**).

For conservative species, there were strong co-occurrence patterns between members of the classes *Betaproteobacteria*, *Bacilli* and *Fusobacteria*, all involved in SOM degradation. By contrast, co-occurrence for root exudate consumers was observed between members of the classes *Actinobacteria*, *Sphingobacteriia*, *Clostridia*, *Bacilli*, *Alphaproteobacteria*, *Betaproteobacteria*, *Gammaproteobacteria*, *Deltaproteobacteria*, and *Verrucomicrobia* (**Supplementary Table S11**). For exploitative species, strong co-occurrence was observed between members of the classes *Flavobacteriia*, *Alphaproteobacteria*, *Gammaproteobacteria*, *Betaproteobacteria*, *Actinobacteria* and *Verrucomicrobiae*, all involved in root exudate assimilation. Notably, strong co-occurrence was observed between different classes such as *Alphaproteobacteria* (*Nitrobacter*), *Nitrospira*, *Actinobacteria*, and *Acidobacteria*. These taxa presented greater versatility, as they were found in the RAS of conservative and exploitative species, and in association with SOM degradation and root exudate assimilation (**Supplementary Table S11**).
# Co-occurrence Network Analysis for Root-Associated Microbiota

There were 244 nodes and 464 edges in the network representing the root colonizing bacteria. The average number of neighbors for a given node was 3.8, while the clustering coefficient was 0.33 with a characteristic path length of 5.9. The root network is therefore less clustered, and the nodes are less linked to each other.

On the root tissues, regardless of the plant nutrient use strategy, strong co-occurrence patterns were observed between members of the *Proteobacteria*, *Bacteroidetes*, and *Chloroflexi* phyla. However, some co-occurrences were specifically observed in the root tissues of conservative and exploitative species (Figure 5B and Supplementary Table S11). For instance, members of the *Actinobacteria* tended to co-occur exclusively on the roots of exploitative species, whereas members of the *Verrucomicrobia* appeared to be more restricted to the roots of conservative plants (Figure 5B).

# DISCUSSION

In the past decade, numerous studies have documented the importance of plant resource use strategies in community structure and ecosystem functioning (Orwin et al., 2010; Grigulis et al., 2013; Cantarel et al., 2015). However the influence of plant nutrient use strategies through root exudation on microbial activities and diversity remains unexplored. Here, we have investigated how plant nutrient resource strategies influence the active microbiota that inhabit the RAS and colonize the root system through root exudation.

# The Influence of Root Exudation Levels on Microbial Activities Among Plant Nutrient Use Strategies

Root exudates represent the main energy-labile C source for microbiota and are known to affect their structure, activity, abundance and diversity (Grayston et al., 1997; Berg, 2009). In our previous study, we demonstrated that exploitative species exude more carbon than conservative species (Guyonnet et al., 2018), probably due to their larger photosynthetic capacity (Personeni and Loiseau, 2004; De Deyn et al., 2008). Even though root exudation rate of AO and DG plant is somewhat similar to the conservative species, these plants are classified as exploitative species as it present high specific leaf area (SLA) and specific root length (SRL), but low leaf dry mass content (LDMC) and root dry mass content (RDMC) (Guyonnet et al., 2018). Plant nutrient management strategies, which influence root exudate activity (Guyonnet et al., 2018) and quality (Guyonnet et al., 2017), could therefore impact microbial functioning in the RAS. Indeed, SIR and DEA in the RAS were found to increase along with root exudation in the rhizosphere of exploitative plants. Several studies have described an increase in DEA coupled to increased C resources that are released via root exudation, e.g., sugars, organic acids, or amino acids (Jones et al., 2004; Mounier et al., 2004; Singh B.K. et al., 2004; Henry et al., 2008). Furthermore, it has been shown that artificial root exudate

supply in the soil can increase DEA (Henry et al., 2008). This increase in denitrifying activity can be explained by the fact that denitrification is a heterotrophic activity, and denitrifying microorganisms are able to use root exudates as C-source. In this study, we demonstrated that an increase in DEA is positively correlated with root exudation in RAS, indicating that exploitative grass species characterized by high root exudation activity could stimulate more DEA in the RAS. Rates of DEA and SIR did not differ between exploitative plant species, despite their differences in root exudation activity.

# Influence of Root Exudation on Plant Microbiota

In this study, SIP was applied to root exudates in the rhizosphere of six plant species presenting different nutrient and soil colonization strategies, allowing the differentiation of the active microbiota involved in root exudate assimilation from bacterial communities that degrade SOM.

Differences were observed between light- and heavy-DNA fractions derived from RAS, thereby indicating ¹³C-root exudate assimilation by certain bacterial communities. Indeed, bacterial diversity analysis revealed that the *Sphingobacteria* class is influenced by root exudation in the RAS. *Sphingobacteria* belongs to the *Bacteroidetes* phylum, and these bacteria exhibit copiotrophic attributes (Fierer et al., 2007). These taxa thrive in conditions of elevated C availability (rhizosphere), and display relatively rapid growth rates (Fierer et al., 2012). The genus *Pseudomonas* exhibited significantly increased relative abundance in the heavy-DNA fraction as compared to the light-DNA fraction and root-DNA fractions, confirming that this bacterial genus is very efficient in root exudate assimilation in the RAS.

Unlabeled but reactive populations correspond to organisms probably specialized in the degradation of stable SOM (Fontaine et al., 2003). Furthermore, microorganisms such as *Thermoleophilia* (belonging to the phylum *Actinobacteria*) and certain *Alphaproteobacteria* may produce SOM degrading enzymes, enabling growth at the expense of SOM. *Thermoleophilia* are also known to degrade complex carbon sources in soil (Coombs and Franco, 2003).

The relative abundance of the microbiota inhabiting the root tissues is higher than the abundance of the one inhabiting the RAS. By contrast, the bacterial diversity is lower before the diffusion of nutrients into the RAS, indicating the specific recognition and nutritional selection of microbiota on root tissues. In agreement with several reports (Haichar et al., 2008; El Khalloufi et al., 2016), our data, based on alpha-diversity and inverted-Simpson indices, confirm that the plant selects bacterial populations in the rhizosphere, and more specifically on the root tissues where more intimate plant-bacteria interactions occur. For this reason, the root system is considered as the most selective habitat.

Several bacterial groups that colonize the root tissues, such as *Betaproteobacteria* or *Sphingobacteria* (within the *Bacteroidetes* phylum), are known to be copiotrophic bacteria (Fierer et al., 2012). These bacteria use labile carbon for their metabolism and
exhibit a fast growth cycle (r-type strategy). In contrast, some other classes are less abundant in the roots than in the RAS, such as *Actinobacteria*, *Acidobacteria* and *Nistrospira*. These classes, known to be oligotrophic bacteria (Fierer et al., 2007; Koch et al., 2015), are more competitive in different environments with low nutrient concentration and display a low growth cycle (K-type strategy). Remarkably, *Acidobacteria* were overrepresented in the RAS but relatively depleted on root tissues, as previously reported in the rhizosphere of maize (Peiffer et al., 2013) and *Medicago sativa* (El Khalloufi et al., 2016).

*Fibrobacteria* were more abundant on the roots of *A. odoratum* and *S. caerulea* than in any other plants. In addition, this class, which appears to be more abundant in the AO rhizosphere, could be considered as generalist bacteria (Haichar et al., 2008) since it was involved in root exudation assimilation and SOM degradation. Remarkably, *Actinobacteria* colonized the root tissues of all plants in our study, even though they displayed an unequal distribution and were more abundant in the roots of *B. erectus* and *D. glomerata* plants (which presented the largest quantity of carbon exuded into the soil). Although *Actinobacteria* are more frequently located in the rhizosphere soil of many plant species (Barka et al., 2016), some of them are known to form more intimate associations with plants and to colonize their internal tissues, as previously observed for wheat (Coombs and Franco, 2003) and *Arabidopsis thaliana* (Bulgarelli et al., 2013) plants.

## Influence of Plant Nutrient Use Strategies on Active Microbiota

Overall, taxa from the exploitative and conservative plants could be separated into two general categories according to their location: generalists, which are broadly distributed among plants from each plant strategy and represent core microbiota of each fraction (SOM degradation, root exudate assimilation in RAS and on root tissues); and specialists, which are locally abundant in one species and specifically involved in SOM degradation or root exudate assimilation in RAS and on root tissues.

The plants that we examined in this study presented one of two different strategies. Specifically, plants were either conservative species characterized by a slow growth rate, low nutrient uptake capacities and a low root exudate amount, or exploitative species characterized by a fast growth rate and a high root exudate amount (Gross et al., 2007; Guyonnet et al., 2018). Nevertheless, we did observe a core microbiota specialist for SOM degradation and root exudate assimilation in RAS and on root tissues. Our results suggest that these taxa are highly competitive and are able to adapt to the two contrasting plant strategies.

By exuding more carbon into the rhizosphere, exploitative species attract more taxa on root tissues and RAS. They also stimulate more taxa involved in SOM degradation by "priming effect" mechanism (Fontaine et al., 2003). Moreover, Guyonnet et al. (2017) have noted that the composition of primary metabolites is different between exploitative and conservative species. Indeed, exploitative species exude more diverse primary metabolites (i.e., sugar, organic acids), whereas conservative species exude more amino acids. Therefore, in addition to root exudate amount, its composition can be critical

in regulating plant microbiota, as also demonstrated by Shade and Handelsman (2011) for wheat and lupin plants.

## Meta-Community Co-occurrence Network

In this study, networks were built not only to determine the links between taxa but also between taxa and plant nutrient use strategies, according to bacterial location (RAS vs. root tissues) and the carbon source (root exudates vs. SOM). Globally, the microbial network of the root is less clustered and connected in comparison to that of the RAS. Although this is related to the smaller number of nodes, it is also linked to the niches that are less delimited than in the RAS.

Some of the observed co-occurrence patterns reveal or confirm interesting ecological patterns for taxa that have not been well-studied. For instance, members of *Verrucomicrobia* tend to co-occur exclusively on the roots of conservative species, suggesting that even though they are abundant and ubiquitous in soils (Sangwan, 2005; Bergmann et al., 2011; Barberán et al., 2012), they may be able to compete for and inhabit more specific and selective niches. In the phylum *Nitrospirae*, *Nitrospira* can co-occur in the RAS of conservative and exploitative species. We observed that *Nitrospira* was associated with SOM degradation and assimilation of root exudates, suggesting that *Nitrospirae* acts as a generalist and is stimulated by plants *via* both strategies. In addition, *Nitrospira*, which can perform complete nitrification (Daim et al., 2015), could provide nitrate to plants, since both the conservative and exploitative species analyzed in previous work present high affinities for nitrate (Grassein et al., 2015).

## CONCLUSION

Our study demonstrates that plant nutrient use strategies impact microbial activities (respiration and denitrification) *via* root exudation, as well as the diversity of active microbiota involved in root exudate assimilation and SOM degradation. The differences in denitrification activity observed among plant nutrient use strategies could be explained partly by the bacterial diversity, since we detected differences between the plant strategies. These results therefore establish that plants can stimulate microbial activities *via* root exudation, by modifying the diversity of plant microbiota. In addition to bacteria, plants could stimulate other microorganisms in the RAS such as fungi that play a key role in denitrification process (Maeda et al., 2015) and SOM degradation (Buée et al., 2009), and which could also be affected by plant nutrient management strategy.

## AUTHOR CONTRIBUTIONS

FZH designed the experiments. JG conducted the experiments. LS performed the ¹³C measurement. JG, AD, MG, MB, PO, WA, TH, and FZH analyzed the results. FZH, JG, WA, and TH wrote the manuscript. All authors have read and approved the final manuscript.
## FUNDING

This work was supported by the French National programme "EC2CO-MicrobiEn, CNRS" (RhizoDen Project).

## ACKNOWLEDGMENTS

We thank Elise Lacroix for her help in plant culture and harvesting. Plant culture was performed at the greenhouse platform (FR BioEnviS, University Lyon 1).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01662/full#supplementary-material

**FIGURE S1** | Bacterial diversity in the rhizosphere of *Bromus erectus* (BE), *Anthoxanthum odoratum* (AO), *Dactylis glomerata* (DG), *Trisetum flavescens* (TF), *Festuca paniculata* (FP) and *Sesleria caerulea* (SC), as characterized by rarefaction curves of OTUs clustered using the swarm method with *d* = 3. H: heavy-DNA fraction; L: light-DNA fraction; R: root-DNA.

**FIGURE S2** | Principal coordinate analysis (PCoA) profile of the active bacterial diversity (using Bray-Curtis distances) colonizing the RAS and inhabiting the root tissues of *Bromus erectus* (BE), *Anthoxanthum odoratum* (AO), *Dactylis glomerata* (DG), *Trisetum flavescens* (TF), *Festuca paniculata* (FP), and *Sesleria caerulea* (SC). (A) PCoA of active microbiota obtained from light- and heavy-DNA fractions from the BS, and from the RAS and root-DNA fractions. (B) PCoA of active microbiota involved in SOM degradation (light-DNA fractions) and in root exudate consumption (heavy-DNA fractions). (C) PCoA of active microbiota colonizing the root tissues.

**TABLE S1** | Sample coverage, species richness and species diversity indices.

**TABLE S2** | Taxonomic affiliation of OTUs retrieved from light-DNA fractions (SOM degraders): (file 1) shared by exploitative plants (DG, TF, BE, and AO) and representing the core microbiota; (file 2) present in AO only; (file 3) present in BE only; (file 4) present in DG only; (file 5) present in TF only.

**TABLE S3** | Taxonomic affiliation of OTUs retrieved from heavy-DNA fractions (root exudate consumers): (file 1) shared by exploitative plants (DG, TF, BE, and AO)

and representing the core microbiota; (file 2) present in AO only; (file 3) present in BE only; (file 4) present in DG only; (file 5) present in TF only.

**TABLE S4** | Taxonomic affiliation of OTUs retrieved from root-DNA fractions (root exudate consumers): (file 1) shared by exploitative plants (DG, TF, BE, and AO) and representing the core microbiota; (file 2) present in AO only; (file 3) present in BE only; (file 4) present in DG only; (file 5) present in TF only.

**TABLE S5** | Taxonomic affiliation of OTUs retrieved from light-DNA fractions (SOM degraders): (file 1) shared by conservative plants (FP and SC) and representing the core microbiota; (file 2) present in FP only; (file 3) present in SC only.

**TABLE S6** | Taxonomic affiliation of OTUs retrieved from heavy-DNA fractions (root exudate consumers): (file 1) shared by conservative plants (FP and SC) and representing the core microbiota; (file 2) present in FP only; (file 3) present in SC only.

**TABLE S7** | Taxonomic affiliation of OTUs retrieved from root-DNA fractions (root exudate consumers): (file 1) shared by conservative plants (FP and SC) and representing the core microbiota; (file 2) present in FP only; (file 3) present in SC only.

**TABLE S8** | Taxonomic affiliation of OTUs retrieved from light-DNA fractions (SOM degraders): (file 1) shared by exploitative (DG, BE, AO, and TF) and conservative plants (FP and SC) and representing the core microbiota; (file 2) present in exploitative plants only; (file 3) present in conservative plants only.

**TABLE S9** | Taxonomic affiliation of OTUs retrieved from heavy-DNA fractions (root exudate consumers): (file 1) shared by exploitative (DG, BE, AO and TF) and conservative plants (FP and SC) and representing the core microbiota; (file 2) present in exploitative plants only; (file 3) present in conservative plants only.

**TABLE S10** | Taxonomic affiliation of OTUs retrieved from root-DNA fractions (root exudate consumers): (file 1) shared by exploitative (DG, BE, AO and TF) and conservative plants (FP and SC) and representing the core microbiota; (file 2) present in exploitative plants only; (file 3) present in conservative plants only.

**TABLE S11** | Table of the different bacteria present in the co-occurrence analysis (Figure 5). Each sheet displays the taxonomic affiliation of the bacterial communities uncovered for the RAS, including root exudate consumers of exploitative plant species (Exploit_heavy DNA fractions), root exudate consumers of conservative plant species (Cons_heavy DNA fractions), SOM degraders from the rhizosphere of conservative species (Cons_light DNA fractions), and the bacteria common to both strategies. The same information was prepared for the bacteria found in the root network (Figure 5), and the position on the gradient from conservative to exploitative was reported.

## REFERENCES

Aerts, R., and Chapin, I. I. I. (1999). "The mineral nutrition of wild plants revisited: a re-evaluation of processes and patterns," in *Advances in Ecological Research*, ed. D. G. Raffaelli (Cambridge, MA: Academic Press), 1-67.

Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., and Vivanco, J. M. (2006). The role of root exudates in rhizosphere interactions with plants and other organisms. *Annu. Rev. Plant Biol.* 57, 233-266. doi: 10.1146/annurev.arplant.57.032905.105159

Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. *ISME J.* 6, 343-351. doi: 10.1038/ismej.2011.119

Bardon, C., Piola, F., Bellvert, F., Haichar, F. Z., Comte, G., Meiffren, G., et al. (2014). Evidence for biological denitrification inhibition (BDI) by plant secondary metabolites. *New Phytol.* 204, 620-630. doi: 10.1111/nph.12944

Barka, E. A., Vatsa, P., Sanchez, L., Gaveau-Vaillant, N., Jacquard, C., Klenk, H.-P., et al. (2016). Taxonomy, physiology, and natural products of actinobacteria. *Microbiol. Mol. Biol. Rev.* 80, 1-43. doi: 10.1128/MMBR.00019-15

Berg, G. (2009). Plant-microbe interactions promoting plant growth and health: perspectives for controlled use of microorganisms in agriculture. *Appl. Microbiol. Biotechnol.* 84, 11-18. doi: 10.1007/s00253-009-2092-7

Bergmann, G. T., Bates, S. T., Eilers, K. G., Lauber, C. L., Caporaso, J. G., Walters, W. A., et al. (2011). The under-recognized dominance of *Verrucomicrobia* in soil bacterial communities. *Soil Biol. Biochem.* 43, 1450-1455. doi: 10.1016/j.soilbio.2011.03.012

Bouchenak-Khelladi, Y., Salamin, N., Savolainen, V., Forest, F., Bank, M., van der Chase, M. W., et al. (2008). Large multi-gene phylogenetic trees of the grasses (Poaceae): progress towards complete tribal and generic level sampling. *Mol. Phylogenet. Evol.* 47, 488-505. doi: 10.1016/j.ympev.2008.01.035

Buée, M., Boer, W. D., Martin, F., van Overbeek, L., and Jurkevitch, E. (2009). The rhizosphere zoo: an overview of plant-associated communities of microorganisms, including phages, bacteria, archaea, and fungi, and of some of their structuring factors. *Plant Soil* 321, 189-212. doi: 10.1007/s11104-009-9991-3

Bulgarelli, D., Schlaeppi, K., Spaepen, S., van Themaat, E. V. L., and 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

Cantarel, A. A. M., Pommier, T., Desclos-Theveniau, M., Diquélou, S., Dumont, M., Grassein, F., et al. (2015). Using plant traits to explain plant-microbe relationships involved in nitrogen acquisition. *Ecology* 96, 788-799. doi: 10.1890/13-2107.1
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. *Nat. Methods* 7, 335–336. doi: 10.1038/nmeth. f.303

Coombs, J. T., and Franco, C. M. M. (2003). Isolation and identification of actinobacteria from surface-sterilized wheat roots. *Appl. Environ. Microbiol.* 69, 5603–5608. doi: 10.1128/AEM.69.9.5603-5608.2003

Daims, H., Lebedeva, E. V., Pjevac, P., Han, P., Herbold, C., Albertsen, M., et al. (2015). Complete nitrification by *Nitrospira* bacteria. *Nature* 528, 504–509. doi: 10.1038/nature16461

De Deyn, G. B., Cornelissen, J. H. C., and Bardgett, R. D. (2008). Plant functional traits and soil carbon sequestration in contrasting biomes. *Ecol. Lett.* 11, 516–531. doi: 10.1111/j.1461-0248.2008.01164.x

El Khaloufi, F., Oufdou, K., Bertrand, M., Lahrouni, M., Oudra, B., Ortet, P., et al. (2016). Microbiote shift in the *Medicago sativa* rhizosphere in response to cyanotoxins extract exposure. *Sci. Total Environ.* 539, 135–142. doi: 10.1016/j.scitotenv.2015.08.127

Escudie, F., Auer, L., Bernard, M., Cauquil, L., Vidal, K., Maman, S., et al. (2015). *FROGS* find rapidly OTU with galaxy solution. *Bioinformatics* 34, 1287–1294. doi: 10.1093/bioinformatics/btx791

Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. *Ecology* 88, 1354–1364. doi: 10.1890/05-1839

Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J., Bradford, M. A., and Knight, R. (2012). Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. *ISME J.* 6, 1007–1017. doi: 10.1038/ismej.2011.159

Fontaine, S., Mariotti, A., and Abbadie, L. (2003). The priming effect of organic matter: a question of microbial competition? *Soil Biol. Biochem.* 35, 837–843. doi: 10.1016/S0038-0717(03)00123-8

Grassein, F., Lemauviel-Lavenant, S., Lavorel, S., Bahn, M., Bardgett, R. D., Desclos-Theveniau, M., et al. (2015). Relationships between functional traits and inorganic nitrogen acquisition among eight contrasting European grass species. *Ann. Bot.* 115, 107–115. doi: 10.1093/aob/mcu233

Grayston, S. J., Vaughan, D., and Jones, D. (1997). Rhizosphere carbon flow in trees, in comparison with annual plants: the importance of root exudation and its impact on microbial activity and nutrient availability. *Appl. Soil Ecol.* 5, 29–56. doi: 10.1016/S0929-1393(96)00126-6

Grigulis, K., Lavorel, S., Krainer, U., Legay, N., Baxendale, C., Dumont, M., et al. (2013). Relative contributions of plant traits and soil microbial properties to mountain grassland ecosystem services. *J. Ecol.* 101, 47–57. doi: 10.1111/1365-2745.12014

Gross, N., Kunstler, G., Liancourt, P., De Bello, F., Suding, K. N., and Lavorel, S. (2009). Linking individual response to biotic interactions with community structure: a trait-based framework. *Funct. Ecol.* 23, 1167–1178. doi: 10.1111/j.1365-2435.2009.01591.x

Gross, N., Suding, K. N., and Lavorel, S. (2007). Leaf dry matter content and lateral spread predict response to land use change for six subalpine grassland species. *J. Veg. Sci.* 18, 289–300. doi: 10.1111/j.1654-1103.2007.tb02540.x

Guyonnet, J. P., Cantarel, A. A. M., Simon, L., and Haichar, F. Z. (2018). Root exudation rate as a functional trait involved in plant nutrient-use strategy classification. *Ecol. Evol.* 8, 8573–8581. doi: 10.1002/ece3.4383

Guyonnet, J. P., Vautrin, F., Meiffren, G., Labois, C., Cantarel, A. A. M., Michalet, S., et al. (2017). The effects of plant nutritional strategy on soil microbial denitrification activity through rhizosphere primary metabolites. *FEMS Microbiol. Ecol.* 93:fix022. doi: 10.1093/femsec/fix022

Hacquard, S., Garrido-Oter, R., González, A., Spaepen, S., Ackermann, G., Lebeis, S., et al. (2015). Microbiota and host nutrition across plant and animal kingdoms. *Cell Host Microbe* 17, 603–616. doi: 10.1016/j.chom.2015.04.009

Haichar, F. Z., Achouak, W., Christen, R., Heulin, T., Marol, C., Marais, M.-F., et al. (2007). Identification of cellulolytic bacteria in soil by stable isotope probing. *Environ. Microbiol.* 9, 625–634. doi: 10.1111/j.1462-2920.2006.01182.x

Haichar, F. Z., Marol, C., Berge, O., Rangel-Castro, J. I., Prosser, J. I., Balesdent, J., et al. (2008). Plant host habitat and root exudates shape soil bacterial community structure. *ISME J.* 2, 1221–1230. doi: 10.1038/ismej.2008.80

Haichar, F. Z., Roncato, M.-A., and Achouak, W. (2012). Stable isotope probing of bacterial community structure and gene expression in the rhizosphere of

*Arabidopsis thaliana*. FEMS Microbiol. Ecol. 81, 291–302. doi: 10.1111/j.1574-6941.2012.01345.x

Haichar, F. Z., Santaella, C., Heulin, T., and Achouak, W. (2014). Root exudates mediated interactions belowground. *Soil Biol. Biochem.* 77, 69–80. doi: 10.1016/j.soilbio.2014.06.017

Henry, S., Texier, S., Hallet, S., Bru, D., Dambreville, C., Chèneby, D., et al. (2008). Disentangling the rhizosphere effect on nitrate reducers and denitrifiers: insight into the role of root exudates. *Environ. Microbiol.* 10, 3082–3092. doi: 10.1111/j.1462-2920.2008.01599.x

Jones, D. L., Hodge, A., and Kuzyakov, Y. (2004). Plant and mycorrhizal regulation of rhizodeposition. *New Phytol.* 163, 459–480. doi: 10.1111/j.1469-8137.2004.01130.x

Koch, H., Lücker, S., Albertsen, M., Kitzinger, K., Herbold, C., Spieck, E., et al. (2015). Expanded metabolic versatility of ubiquitous nitrite-oxidizing bacteria from the genus *Nitrospira*. *Proc. Natl. Acad. Sci. U.S.A.* 112, 11371–11376. doi: 10.1073/pnas.1506533112

Lambers, H., Raven, J. A., Shaver, G. R., and Smith, S. E. (2008). Plant nutrient-acquisition strategies change with soil age. *Trends Ecol. Evol.* 23, 95–103. doi: 10.1016/j.tree.2007.10.008

Maeda, K., Spor, A., Edel-Hermann, V., Heraud, C., Breuil, M.-C., Bizouard, F., et al. (2015). N2O production: a widespread trait in fungi. *Sci. Rep.* 5:9697. doi: 10.1038/srep09697

Mahé, F., Rognes, T., Quince, C., de Vargas, C., and Dunthorn, M. (2014). Swarm: robust and fast clustering method for amplicon-based studies. *PeerJ* 2:e593. doi: 10.7717/peerj.593

Maire, V., Gross, N., Da Silveira Pontes, L., Picon-Cochard, C., and Soussana, J.-F. (2009). Trade-off between root nitrogen acquisition and shoot nitrogen utilization across 13 co-occurring pasture grass species. *Funct. Ecol.* 23, 668–679. doi: 10.1111/j.1365-2435.2009.01557.x

McGill, B. J., Enquist, B. J., Weiher, E., and Westoby, M. (2006). Rebuilding community ecology from functional traits. *Trends Ecol. Evol.* 21, 178–185. doi: 10.1016/j.tree.2006.02.002

McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. *PLoS One* 8:e61217. doi: 10.1371/journal.pone.0061217

Mendes, R., Kruij, M., Bruijn, I., de Dekkers, E., Voort, M., van der Schneider, J. H. M., et al. (2011). Deciphering the rhizosphere microbiome for disease-suppressive bacteria. *Science* 332, 1097–1100. doi: 10.1126/science.1203980

Mounier, E., Hallet, S., Chèneby, D., Benizri, E., Gruet, Y., Nguyen, C., et al. (2004). Influence of maize mucilage on the diversity and activity of the denitrifying community. *Environ. Microbiol.* 6, 301–312. doi: 10.1111/j.1462-2920.2004.00571.x

Orwin, K. H., Buckland, S. M., Johnson, D., Turner, B. L., Smart, S., Oakley, S., et al. (2010). Linkages of plant traits to soil properties and the functioning of temperate grassland. *J. Ecol.* 98, 1074–1083. doi: 10.1111/j.1365-2745.2010.01679.x

Peiffer, J. A., Spor, A., Koren, O., Jin, Z., Tringe, S. G., Dangl, J. L., et al. (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. *Proc. Natl. Acad. Sci. U.S.A.* 110, 6548–6553. doi: 10.1073/pnas.1302837110

Personeni, E., and Loiseau, P. (2004). How does the nature of living and dead roots affect the residence time of carbon in the root litter continuum? *Plant Soil* 267, 129–141. doi: 10.1007/s11104-005-4656-3

Ranjard, L., Lejon, D. P. H., Mougel, C., Schehrer, L., Merdinoglu, D., and Chaussod, R. (2003). Sampling strategy in molecular microbial ecology: influence of soil sample size on DNA fingerprinting analysis of fungal and bacterial communities. *Environ. Microbiol.* 5, 1111–1120. doi: 10.1046/j.1462-2920.2003.00521.x

Ryser, P., and Lambers, H. (1995). Root and leaf attributes accounting for the performance of fast- and slow-growing grasses at different nutrient supply. *Plant Soil* 170, 251–265. doi: 10.1007/BF00010478

Samain, E., van Tuinen, D., and Selim, S. (2017). The plant-growth-promoting rhizobacterium *Paenibacillus* sp. strain B2 stimulates wheat defense mechanisms against septoria leaf blotch and root colonization by Curtobacterium plantarum. *Biol. Control* 114, 87–96. doi: 10.1016/j.biocontrol.2017.07.012

Sangwan, S. (2005). "Virtual community success: a uses and gratifications perspective," in *Proceedings of the 38th Annual Hawaii International
*Conference on System Sciences*, (Washington, DC), 193c. doi: 10.1109/HICSS. 2005.673

Shade, A., and Handelsman, J. (2011). *Beyond the Venn Diagram: The Hunt For a Core Microbiomeemi_2585*. Available at: http://handelsmanlab.yale.edu/sites/default/files/files/BeyondtheVenndiagram2011.pdf

Singh, A., Bicudo, J. R., Tinoco, A. L., Tinoco, I. F., Gates, R. S., Casey, K. D., et al. (2004). Characterization of nutrients in built-up broiler litter using trench and random walk sampling methods. *J. Appl. Poult. Res.* 13, 426-432. doi: 10.1093/japr/13.3.426

Singh, B. K., Millard, P., Whiteley, A. S., and Murrell, J. C. (2004). Unravelling rhizosphere-microbial interactions: opportunities and limitations. *Trends Microbiol.* 12, 386-393. doi: 10.1016/j.tim.2004.06.008

R Core Team (2015). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Tilman, D. (1990). Constraints and tradeoffs: toward a predictive theory of competition and succession. *Oikos* 58, 3-15. doi: 10.2307/3565355

Violle, C., Navas, M.-L., Vile, D., Kazakou, E., Fortunel, C., Hummel, I., et al. (2007). Let the concept of trait be functional. *Oikos* 116, 882-892. doi: 10.1111/j.0030-1299.2007.15559.x

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (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

Williams, R. J., Howe, A., and Hofmockel, K. S. (2014). Demonstrating microbial co-occurrence pattern analyses within and between ecosystems. *Front. Microbiol.* 5:358. doi: 10.3389/fmicb.2014.00358

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Guyonnet, Guillemet, Dubost, Simon, Ortet, Barakat, Heulin, Achouak and Haichar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.