Next Article in Journal
Sprouty3, but Not Sprouty1, Expression Is Beneficial for the Malignant Potential of Osteosarcoma Cells
Previous Article in Journal
Identification and Function of Apicoplast Glutaredoxins in Neospora caninum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spike Formation Is a Turning Point Determining Wheat Root Microbiome Abundance, Structures and Functions

1
Soil, Water and Environmental Sciences, Volcani Research Center, Rishon LeZion 7505101, Israel
2
Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 7610001, Israel
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(21), 11948; https://doi.org/10.3390/ijms222111948
Submission received: 7 October 2021 / Revised: 24 October 2021 / Accepted: 29 October 2021 / Published: 4 November 2021
(This article belongs to the Section Molecular Plant Sciences)

Abstract

:
Root selection of their associated microbiome composition and activities is determined by the plant’s developmental stage and distance from the root. Total gene abundance, structure and functions of root-associated and rhizospheric microbiomes were studied throughout wheat growth season under field conditions. On the root surface, abundance of the well-known wheat colonizers Proteobacteria and Actinobacteria decreased and increased, respectively, during spike formation, whereas abundance of Bacteroidetes was independent of spike formation. Metagenomic analysis combined with functional co-occurrence networks revealed a significant impact of plant developmental stage on its microbiome during the transition from vegetative growth to spike formation. For example, gene functions related to biofilm and sensorial movement, antibiotic production and resistance and carbons and amino acids and their transporters. Genes associated with these functions were also in higher abundance in root vs. the rhizosphere microbiome. We propose that abundance of transporter-encoding genes related to carbon and amino acid, may mirror the availability and utilization of root exudates. Genes related to antibiotic resistance mechanisms were abundant during vegetative growth, while after spike formation, genes related to the biosynthesis of various antibiotics were enriched. This observation suggests that during root colonization and biofilm formation, bacteria cope with competitor’s antibiotics, whereas in the mature biofilm stage, they invest in inhibiting new colonizers. Additionally, there is higher abundance of genes related to denitrification in rhizosphere compared to root-associated microbiome during wheat growth, possibly due to competition with the plant over nitrogen in the root vicinity. We demonstrated functional and phylogenetic division in wheat root zone microbiome in both time and space: pre- and post-spike formation, and root-associated vs. rhizospheric niches. These findings shed light on the dynamics of plant–microbe and microbe–microbe interactions in the developing root zone.

1. Introduction

The rhizosphere is a hotspot for plant and soil microbiome interactions. Rhizosphere microbiome composition and function are mostly shaped by soil type, plant species and plant genotypes [1,2,3,4]. Plants select their associated bacteria by depositing specific exudates in the soil–root interface [5,6]. Roots secrete large amounts of carbon as exudates that can be divided into two major classes of metabolites. The first consists of high-molecular-weight compounds such as polysaccharides and proteins, and the second of low-molecular-weight metabolites, including amino acids, organic acids, ions, sugars, vitamins and various secondary metabolites [7,8,9]. The root tip is the major source of root exudation [10,11,12], along with root elongation zones and root hairs [13,14]. As a result, a gradient of exudates is created from the roots to the surrounding soil, with the highest concentration near the root surface [15]. Thus, the root surface provides a more nutrient-rich environment than the rhizosphere, enabling microorganisms to colonize and proliferate.
Plant roots and their microbiome (both root-associated and rhizosphere) are continuously interacting. The root-associated bacteria can manipulate plant exudates by accelerating their diffusion, as well as altering plant carbon allocation to the roots [16]. In many plant species, specific microbiome functions have been shown to be associated with the microorganism–plant interactions. These functions include motility and adhesion to roots, biofilm formation, metabolism of nitrogen, carbohydrates and vitamins and degradation of xenobiotics [4,17,18]. In Arabidopsis, it has been proposed that roots of young seedlings secrete high levels of sugars to attract colonization of diverse bacteria, whereas in older plants, roots secrete phenolic substances to select specific rhizosphere bacteria [19,20]. A plant’s developmental stage has also been shown to affect root-exudation patterns [16,21]. This developmental stage-dependent exudation was recently shown to influence rhizosphere community structure in broomcorn millet [22] and root-associated community structure in rice [23].
Although the root surface provides an environment that is much richer in nutrients, most studies have focused on the influence of exudates on the rhizospheric bacterial community [1,9,24,25,26] rather than the root-associated one. Moreover, the effects of plant growth stage on rhizosphere community functions, and particularly on the root-associated community, remain unclear. To close this gap, we followed field-grown wheat throughout the growth season from emergence to spike maturation and compared the rhizosphere bacterial community (inhabiting the soil surrounding the roots) to that of root-associated bacteria (sampled from washed roots). We describe and discuss the changes in root-associated and rhizosphere bacterial community structure and functions with stages of plant development.

2. Results

Wheat grown under commercial field conditions was sampled from sprout emergence up to week 17, when some plants produced green spikes whereas in others, the spikes were already more mature. Rhizosphere (soil sampled by shaking from the roots) and washed roots (root-associated) were collected at weeks 1, 2, 5, 9, 12 and 17. The spatial (rhizosphere vs. root-associated) and temporal (before vs. after spike formation) effects on the bacterial community structures and functional profiles were analyzed.

2.1. Influence of Niche and Plant Developmental Stage on Bacterial Community Structure and Functional Profile

After filtering out the plant host genes from those of the root and rhizosphere microbiomes, around 80% of the sequences were related to bacteria, 16% could not be assigned to any known taxonomic group, and the rest were assigned to Eukaryota and Archaea. The Eukaryota, Archaea, and unassigned sequences were filtered out from the root and rhizosphere microbiome and the bacterial sequences were further analyzed; 28% and 20% (root and rhizosphere, respectively) of those sequences were assigned to known KEGG functions.
Clustering patterns of the bacterial community structure (Figure 1A) and gene functions (Figure 1B) in the two niches (i.e., root-associated and rhizosphere) are presented in nMDS orientation; nMDS orientation of bacterial community structure was based on all annotated bacterial taxonomic levels (Figure 1A). Gene functions in the nMDS orientation (Figure 1B) were constructed based on annotation of genes with known KO groups from the KEGG database. Niche origin had a major influence on both bacterial community structure and gene functions (calculated using the Adonis function in the Vegan R package; p < 0.0001). A lesser, but still significant influence of plant age (especially before and after spike formation) on the combined root-associated and rhizosphere bacterial community structure (p < 0.005) and functions was found (p < 0.05). The influence of wheat age was more pronounced in each niche separately than in the two niches combined. Spike formation was the turning point, leading to major shifts in both root-associated and rhizosphere bacterial community structures and gene functions. Thus, throughout the study, we compared the results obtained before and after spike formation. Differences between the niches were found in each week from sprout emergence until spike formation, in the structure of the bacterial communities (Figure 1C) and their gene functions (Figure 1D). However, niche differences for both parameters decreased after spike formation.
In addition, spike formation was also the turning point, leading to significant increase in total bacteria abundance in root-associated community (Figure 2A, e.g., 16S rRNA/TEF and bacteria genes per plant genes from metagenome). Surprisingly, not only spike formation influenced community total root-associated gene function but also the specific functional gene nosZ (Figure 2B), known as sole enzyme capable to reduce N2O, a greenhouse gas [27,28]. With spike formation, abundance of nosZ gene per total root-associated bacterial community decreased, possibly indicating diminished root-associated microbiome capability to reduce N2O. Unlike abundance of nosZ per total root bacterial community, its abundance per rhizosphere community was developmental stage dependent. In the rhizosphere, nosZ abundance per total community have increased with plant age (Figure 2B). This pattern can be explained by the difference in nitrate availability according to the proximity to the root.
Structural co-occurrence networks were constructed based on bacterial phylogeny of all annotated taxonomic levels (nodes) from the metagenome (for simplicity, nodes were marked with bacterial phyla). During wheat growth, the number of interactions (edges) between bacterial communities in the rhizosphere was almost eight times higher than that between root-associated bacterial communities (5019 vs. 641 edges; Figure 3A,B). When the complete network (all time points together) was divided between pre-spike formation and post-spike formation, the number of edges increased: whereas the complete network for root-associated bacterial communities was composed of 641 edges, pre- and post-spike formation networks were composed of 1207 and 2411 edges, respectively (Figure 3A). Similar trends were observed for the rhizosphere (Figure 3B). This supports the hypothesis that spike formation has a very strong influence on the microbiome, affecting both root-associated and rhizosphere bacterial communities (Figure 1A,B and Figure 2).
Networks of functional genes co-occurrence were constructed based on annotated genes with known KO numbers from the KEGG database. These included six categories: cellular processes, environmental information processes, genetic information processes, diseases, metabolism and systems. The number of interactions in the co-occurrence networks of functional genes in the root-associated bacterial community was slightly higher compared to the rhizosphere bacterial community during wheat growth (Figure 3C,D, respectively). Five distinct functional gene clusters in the root-associated community interacted in the co-occurrence networks before spike formation, compared to only one cluster after spike formation (Figure 3C). The biggest cluster in the root-associated microbiome before spike formation was related to metabolism of carbohydrates, amino acids, nucleotides and vitamins. Another cluster was related to processes involved mainly in carbon fixation, which negatively interacted with metabolism of methane, nitrogen and secondary metabolites, and to plant–pathogen interactions. The third cluster was related to processes involved in quorum-sensing, biofilm formation and cell motility. The fourth cluster was related to processes involved in transcription and the immune response. The fifth cluster was related to processes involved in proliferation of bacterial infections, prokaryotic defense system and metabolism of secondary metabolites, including antibiotics. The last cluster found in the root-associated microbiome after spike formation was related to processes involved mainly in carbon fixation, and metabolism of methane, nitrogen and carbohydrates (Figure 3C). In the rhizosphere, on the other hand, three distinct clusters were formed after spike formation, compared to only one before spike formation (Figure 3D). The latter cluster was related to processes involving mainly in carbohydrate metabolism and transcription. The biggest cluster in the rhizosphere after spike formation was related to biosynthesis of secondary metabolites, enzymes, signal transduction, membrane transport, cell motility and amino acid and nucleotide metabolisms. The two other clusters were related to antimicrobial resistance genes, which interacted negatively with energy metabolism, and membrane trafficking, which interacted positively with transcriptional regulators.

2.2. Influence of Niche and Spike Formation on Bacterial Community Composition

At the phylum level, Proteobacteria dominated the wheat root-associated bacterial community (47–67% relative abundance), followed by Actinobacteria (23–42%) and Bacteroidetes (4–6.5%) (Figure 4A). In the rhizosphere, Actinobacteria were dominant (44–50%), followed by Proteobacteria (32–35%) and Acidobacteria (9–12%) (Figure 4A). Proteobacteria, Actinobacteria and Bacteroidetes showed a niche preference and were dependent on plant growth stage. Proteobacteria and Bacteroidetes and their affiliated groups were more abundant in the root-associated bacterial community compared to the rhizosphere, while most Actinobacteria-affiliated groups were more abundant in the rhizosphere (Figure 4B). The abundance of Proteobacteria in the roots was significantly higher before spike formation, whereas the abundance of Actinobacteria was higher after spike formation (Figure 4C). The abundance of Bacteroidetes in the root-associated community was not affected by spike formation. In the rhizosphere, the abundance of Acidobacteria was higher before spike formation, whereas the abundance of Actinobacteria and Proteobacteria was higher after spike formation.

2.3. The Effect of Niche and Spike Formation on Specific Root-Associated and Rhizosphere Gene Functions

To zoom in on functional genes that were significantly enriched in one niche over the other (root or rhizosphere) or before or after spike formation, they were divided into eight functional categories to reduce data complexity. The functional categories were related to: (1) carbons; (2) nitrogen; (3) biofilm and sensorial movement; (4) antibiotic synthesis and resistance; (5) ion transporters; (6) amino acid transporters; (7) vitamin transporters; and (8) metal transporters (Figure 5A). Overall, during wheat growth, the highest enrichment in functional gene abundance in the root-associated compared to rhizosphere bacterial community was related to two categories: biofilm and sensorial movement, and antibiotic synthesis and resistance (Figure 5B).
The biofilm and sensorial movement group included genes related to secretion systems type I, III, IV and VI, secretion systems related to pili and flagella and mannose-sensitive hemagglutinin (MSHA), biofilm formation, motility, quorum-sensing and chemotaxis. Interestingly, 96 genes from this group were abundant in the root-associated community before spike formation, and none were enriched after spike formation (Figure 6A,B). Antibiotic synthesis and resistance genes were separated into two subgroups in relation to the timing of spike formation: multidrug and resistance toward (1) beta-lactams, (2) antimicrobial peptides and (3) invading RNA elements (CRISPR systems) before spike formation, and biosynthesis of various antibiotics, including polyketide, after spike formation (Figure 6B).
The gene cluster related to carbons consisted mainly of transporters and had twice the number of significantly abundant functional genes in the roots compared to the rhizosphere (136 vs. 67 genes; Figure 5B), while amino acid transporters showed almost three times the number of genes (22 vs. 8 genes, respectively). Genes involved in carbon compound biosynthesis and degradation processes were more abundant in roots than in the rhizosphere. In addition, significantly enriched genes related to amino acid transporters (possibly related to availability of root exudates) were two times more abundant before compared to after spike formation in the roots (17 vs. 9 genes) (Figure 6B,A).
Out of the genes related to the nitrogen cycle, two gene groups (urea transporter and denitrification) showed significant niche preference. Genes related to the urea transporter are more abundant in root-associated bacterial community, while genes related to denitrification are more abundant in the rhizosphere community (Figure 5B).
In the rhizosphere, the number of significantly enriched genes assigned to one of the eight selected functional groups (i.e., related to carbons, nitrogen, biofilm and sensorial movement, antibiotic synthesis and resistance, ion transporters, amino acids, vitamins and metal transporters) was much smaller (23%) compared to root-associated bacteria (116 vs. 503 genes), respectively. The significantly enriched genes in the rhizosphere belonged to three groups: (1) carbon transporters, biosynthesis and metabolism; (2) antibiotic biosynthesis and resistance and (3) amino acid biosynthesis and metabolism (Figure 6A,C). Significantly abundant genes in the rhizosphere assigned to antibiotic biosynthesis and resistance and to amino acid biosynthesis and metabolism were more abundant before spike formation.
It is interesting to note that only nine functions (lipopolysaccharide transporter and biosynthesis, beta-lactam resistance, metabolism of methane, metabolism of alanine, aspartate and glutamate, metabolism of arginine and proline and biosynthesis of peptidoglycan, lysine and enediyne antibiotics, Figure 6B) comprised up to 45% of the significantly enriched genes in the rhizosphere. Their abundance was significantly influenced by spike formation in both rhizosphere and roots.

2.4. Links between Taxonomy and Significantly Changed Functions

Links between selected representative bacterial phyla and significantly abundant functional genes before and after spike formation in the roots are shown in Figure 7B. The phyla proteobacteria, Actinobacteria and Bacteroidetes were selected for these analyses for three reasons: (1) these three groups had the highest relative abundance among all phyla in the roots (Figure 3A); (2) Proteobacteria and Bacteroidetes were significantly enriched in the roots compared to the rhizosphere; (3) Proteobacteria and Actinobacteria abundance in the roots was significantly influenced by spike formation. proteobacteria was the dominant phylum in the roots (47–67%), and therefore, most of its significantly enriched gene subgroups were more abundant than those of Actinobacteria and Bacteroidetes, even when calculating abundance per cell (Supplemental Table S3).

3. Discussion

Roots select their associated microbiome members and affect their activity by providing a nutrient-rich environment based on distance from the root [29]. Some of the rhizodepositions subsequently diffuse from the root surface into the rhizosphere [30]. Many crops, including wheat, have been found to secrete increasing amounts of organic carbon for the first 50 days after planting, after which organic carbon deposition decreases dramatically [21]. From that point on, root senescence is occurring rapidly and eventually resulting in root decomposition. We propose, based on metagenomic analyses, that root-associated and rhizosphere community abundance, structure and their gene repertoires respond to a turning point in wheat development (spike formation) that occurs concomitantly with changes in root organic carbon depositions.
In this study, in field grown wheat, root-associated bacteria abundance per root increased after spike formation. This increase in bacteria abundance is mainly attributed to decomposition of wheat roots at the final stages of wheat growth manifested by a decrease in TEF gene copy number while total bacteria abundance remains comparable through weeks 12, 17 green and 17 yellow.
Applying two genomic approaches (DNA sequencing and qPCR) we demonstrated that there is a potential for higher abundance of genes related to denitrification in rhizosphere compared to root-associated microbiome during wheat growth, possibly due to competition with the plant over nitrogen in the root vicinity. Based on the abundance of urea transporter genes, root-associated bacteria possibly import urea, used as a fertilizer in the field, as a source of nitrogen for building amino acids and nucleotides.
Spatial-temporal dynamics of root development leads to a constant need for bacteria to follow the root and re-assemble the associated microbiome [31,32]. Root-associated bacteria have important traits that provide them with a selective advantage in colonization. These include, among others, attachment and competition mechanisms such as motility and chemotaxis, that enable sensing and reaching the root surface [33,34]. It has been shown that non-motile or reduced-motility mutants are unable to compete in root colonization [35]. Another trait is biofilm formation, which creates a better protected niche [36]. In this study, three microbiome functional clusters were shown to be specifically associated with root proximity: biofilm and sensorial movement; antibiotic production and resistance; carbon biosynthesis, degradation and transporters, and amino acid biosynthesis and transporters. The first cluster includes genes related to secretion systems (type I, II, III, IV and VI, pili, flagella and MSHA), biofilm formation, motility, quorum-sensing and chemotaxis. In the second cluster, half of the genes were associated with biosynthesis of polyketide antibiotics which are known to facilitate competition and communication between microorganisms [37]. The third gene cluster included transporters for carbon and amino acids which were significantly enriched in the roots compared to the rhizosphere. Transporters of fructose, glycerol and glucose and mannose were enriched in the root-associated bacteria compared to the rhizosphere. Indeed, those compounds have been recently identified as wheat root exudates [38]. Based on root-associated bacterial gene enrichment, it may be concluded that during root colonization and growth, bacteria sense the presence of emerging roots and move toward them from the bulk soil, express resistance genes to compete with other bacteria, then utilize the exudates, proliferate on the root surface and form a biofilm.
Interestingly, in the root-associated microbiome, the three functional gene clusters were negatively influenced at the time of spike formation. The biofilm and sensorial movement genes were abundant before spike formation, while the roots are growing and bacterial colonization is actively occurring, but none of these genes were abundant after spike formation. Moreover, genes related to antibiotic resistance mechanisms were abundant before spike formation. These include genes coding for resistance to beta-lactams and antimicrobial peptides, both provide tools to remove competing organisms, as well as CRISPR/Cmr systems that provide immune responses directed against the RNA of invading elements [39,40]. However, after spike formation, genes related to the biosynthesis of various antibiotics, including polyketides, were enriched. This observation suggests that during colonization and biofilm formation, bacteria need genes that will allow them to cope with antibiotics producers, whereas in the mature biofilm stage, they invest in inhibiting new colonizers.
Before and after spike formation, a large portion of significantly abundant genes were related to carbon and amino acid transporters, especially ABC transporters. Amino acids, important component of root exudates, are crucial for the early steps of root colonization [41] and indeed, more than twice the amount of amino acid-related genes were abundant in the root microbiome before as compared to after spike formation. Changes in gene abundance of those bacterial carbon and amino acid transporters may probably reflect the carbohydrate and amino acid quantities and types of exudates produced at different stages of plant development. For instance, glucose and mannose, glycerol, inositol, methyl galactoside and trehalose are known wheat root exudates [38], and indeed their transporters were enriched in the microbiome before spike formation. However, transporters for fructose, sorbitol and mannitol, also known wheat root exudates [38], were abundant in the root microbiome after spike formation. Chemical assessment of the actual composition of root exudates in the soil is difficult and inaccurate [42], thus, it is suggested that microbiome-abundant transporter-encoding genes may reflect the composition of the root exudates available to the bacteria in situ.
Nine functions whose gene abundance was significantly influenced by spike formation in both root-associated and rhizosphere niches. Genes related to these functions sum up to half of the total significantly enriched genes in the rhizosphere. Moreover, all of these genes are enriched in the rhizosphere before spike formation, when root exudation is expected to be at its peak [21], supporting the notion that root exudates diffuse into the surrounding soil, and contribute to creating the rhizosphere.
In this study, the wheat root community consisted mostly of proteobacteria, Actinobacteria and Bacteroidetes, known to be wheat root colonizers [43,44,45]. Interestingly, in the root-associated niche, the abundance of proteobacteria and Actinobacteria was influenced by spike formation, whereas that of Bacteroidetes was not. Proteobacteria were significantly more abundant and, per cell, had a higher abundance of significantly enriched genes before spike formation (i.e., genes related to carbon metabolism and transporters, amino acid metabolism, degradation and transporters, antibiotic resistance, quorum-sensing, chemotaxis, biofilm formation, and various types of secretion systems, denitrification and metal and ion transporters) compared to Actinobacteria. On the other hand, Actinobacteria, which were significantly more abundant after spike formation, had more significantly enriched genes at this time point per cell (i.e., genes related to carbon transporters, degradation and metabolism, amino acid transporters and biosynthesis) compared to proteobacteria. Some Actinobacteria have the ability to partially degrade plant cell wall components, such as lignocellulose and other complex molecules [46], which are characteristic of mature and remnants of wheat roots [47]. In addition, genes encoding cellobiose transporter were abundant after spike formation, suggesting consumption of cellobiose, a product of cellulose hydrolysis [48]. These changes in abundance with spike formation may indicate a major role for proteobacteria and Actinobacteria in wheat root colonization by adapting their functional profile to follow root growth and secretion patterns, as well as decay.
We demonstrated functional and phylogenetic division in the microbiome of wheat root zone in both time and space (pre- and post-spike formation, and root-associated vs. rhizospheric niches). As it is not feasible to directly examine mechanisms occurring underground under field conditions, we suggest and discuss several possible plant-microbe interactions based on traits found in the metagenome. Differences in wheat bacterial community structure and gene functions manifest the functional division between niches as well as developmental stages, reflecting plant physiology status during plant and its microbiome co-development in the field. These findings shed light on the dynamics of the fascinating plant–microbe and microbe–microbe interactions in the root zone.

4. Materials and Methods

4.1. Soil and Root Sampling

Wheat (Triticum aestivum cv. Tishrey) was cultivated in a field at the Volcani Center, Rishon LeZion, Israel. The soil was clay loam (36% clay, 22% silt and 42% sand) and the soil parameters prior to wheat sowing were: pH 7.99 ± 0.04, electrical conductivity 99 ± 2 µS/m, N-NO3 0.55 ± 0.05 mg/kg, NH4 0.36 ± 0.02 mg/kg, P-PO4 0.06 ± 0.01 mg/kg, total soluble organic carbon 4.64 ± 0.08 mg/kg, and total soluble nitrogen 0.87 ± 0.01 mg/kg. The field was amended once with cattle urine (48 m3/acre) prior to wheat sowing in November 2015. Weather conditions during wheat growth are provided in Supplemental Figure S1.
During wheat growth and development, microbiomes were sampled as follows: rhizosphere microbiome was collected from the soil dislodged from the roots by shaking. The samples were immediately placed on ice and stored at −80 °C for further analyses. For the root-associated microbiome, roots were washed three times with saline solution (0.85% NaCl) and immediately placed on ice and stored at −80 °C for further analyses. Samples were collected at the following developmental stages: emergence with one or two leaves (1st week, 16 December 2015), tillering (2nd week), jointing (5th week), heading (9th week), anthesis (12th week) and spike maturation (17th week). Root and rhizosphere samples were collected in triplicate in weeks 1, 2, 5, 9 and 12; in week 17, two subsamples were collected for each spike category (i.e., green and mature yellow, Supplemental Figure S2).

4.2. Rhizosphere and Root DNA Extraction for Sequencing and qPCR

To obtain the root-associated microbiome, wheat roots were vortexed three time with 85% saline solution, until no visible soil particles were attached to the roots. Total DNA was extracted from 0.4 g of complete root system and 0.3 g rhizosphere soil, using the Exgene Soil DNA mini-isolation kit (GeneAll, Seoul, Korea) according to the manufacturer’s instructions.

4.3. Generation of qPCR Plasmid Standards

Plasmids containing the 16S rRNA gene, translation elongation factor 1 (TEF, a plant housekeeping gene) and nosZ gene were generated as described previously [49,50,51]. Each PCR amplification product was ligated into pGEM-T Easy Vector (Promega, Madison, WI, USA) and plasmids were transformed into BioSuper Escherichia coli DH5α competent cells (Bio-Lab, Jerusalem, Israel). Circular plasmid DNAs were used as the standards to create calibration curves at 10-fold dilutions for gene quantification by real-time qPCR.

4.4. Assessment of Gene Copy Numbers by qPCR

Copy numbers of total bacterial community (16S rRNA gene), plant housekeeping gene (TEF) and selected denitrifying gene (nosZ) were assessed using selected primers (Supplemental Table S1) in wheat roots using the StepOnePlus Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). Triplicates from whole genomic DNA were diluted up to 6 ng/µL and 1 µL was used in a 20-µL final reaction volume together with 50 µM forward and revers primers and 10 µL 1X FAST or Power SYBR MasterMix (Thermo Fisher Scientific, Waltham, MA, USA). Three biological and three technical replicates were analyzed for each individual soil or root DNA sample. Reaction efficiency was monitored in each run by means of an internal standard curve (constructed plasmids) using duplicates of 10-fold dilutions of standards ranging from 108–102 copies per reaction. Efficiency was 89–95% for all target genes and runs, and R2 values were greater than 0.99. Copy numbers of the target genes were calculated based on the relative calibration curve of plasmid copy numbers. All data analyses were conducted using StepOne software v2.3 (Applied Biosystems, Waltham, MA, USA).

4.5. Shotgun Sequencing

Shotgun metagenomic libraries were prepared as described previously [52]. The final size-selected pool was sequenced in an Illumina NovaSeq instrument with an S4 flow cell, employing 2 × 150 base reads. Library preparation and pooling were performed at the University of Illinois at Chicago Sequencing Core (UICSQC), and sequencing was performed by Novogene Corporation (Chula Vista, CA, USA).
In total, we obtained 1019 Gb of information, with 33–44 million sequences per root sample and 34–53 million sequences per rhizosphere sample. These sequence data were submitted to the Sequence Read Archive (SRA) of the NCBI databases under accession number SUB8180046, BioProject: PRJNA664890.
All reads were subjected to quality control using FastQC v0.11.3 [53] and barcode trimming using Trimmomatics v0.32 [54]. Reads were mapped to the whole wheat metagenome using Bowtie2 v2.3.5.1 [55], and mapped reads (wheat genes) were filtered out from each sample. Then, short Illumina reads in triplicate from each week (1, 2, 5, 9 and 12) were assembled using SPADES v3.13.0 [56] into longer contigs, to create five wheat-root microbiome catalogues, one for each week. Six rhizosphere microbiome catalogues were created, one for each week (1, 2, 5 and 9), and two for week 12 due to its overall size. Reads with four replicates from week 17 were separated into two root microbiome catalogues and two rhizosphere microbiome catalogues. Altogether seven root-associated catalogues and eight rhizosphere catalogues were created (Supplemental Table S2). Prodigal v2.6.2 [57] was used for protein-coding gene prediction for each of those 15 catalogues. To create a non-redundant set of genes for each niche separately, we used CD-HIT-EST software v4.8.1 [58] with a similarity threshold of 95%. Those genes were used as the root-associated gene catalogue, which included 53 million partial genes, and for the rhizosphere gene catalogue, which included 91 million partial genes. These gene catalogues were searched against the non-redundant NCBI protein database using DIAMOND sensitive algorithm v0.9.24.125 [59] to assign taxonomic and functional annotations. Results were then uploaded to MEGAN Ultimate edition software v6.15.2 [60]. The lowest common ancestor (LCA) algorithm was applied (parameters used with minimum bit score of 70, minimum support of 5%, and 30% top threshold) to compute the assignment of genes to specific taxa. For functional annotation, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [61] was used. Following annotation, to generate taxonomic and functional count tables, each library was mapped to the gene catalogue with Trinity mapping software v2.8.4 [62], with Bowtie2-modified parameters (--no-unal --gbar 99999999 -k 250 --dpad 0 --mp 1,1 --np 1 --score-min L,0,-0.9 -L 20 -i S,1,0.50).

4.6. Data Analyses

All sequencing data analyses were performed using R statistical software. A non-metric multidimensional scaling (nMDS) ordination plot was constructed using R package VEGAN v.2.5-5 [63]. The data matrix was transformed using normalized count transformation. For community structure, ordination was generated using the Bray–Curtis dissimilarity matrix, and for community function of KEGG orthologous (KO) groups, ordination was generated using the Euclidean dissimilarity matrix. Permutational multivariate analysis was used to compute the variance between microbial community structure or function and experimental parameters (niche and spike), using the Adonis function in the Vegan R package [64]. Changes in structural and functional similarity between wheat root and rhizosphere microbial communities in each week were calculated using Student’s t-test in JMP 14 Pro software (SAS Institute Inc., Cary, NC, USA) and statistical significance was set at p < 0.05. Network co-occurrence analysis of bacterial communities in root-associated and rhizosphere samples were constructed using Cytoscape v.3.7.2 software’s [65] CoNet app [66]. For comparison of taxonomic changes and functional traits, differential abundance was estimated using DESeq2 [67] and was considered significant when the difference in abundance between genes had a false discovery rate (FDR)-adjusted p-value <0.05. For comparison of changes in taxonomic and functional traits, the bacterial read counts table was binned into KO groups, based on DIAMOND-MEGAN annotation.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms222111948/s1.

Author Contributions

A.U.-T., Y.H. and D.M. jointly designed the experiments and wrote the manuscript. A.U.-T. performed the experiments and analyzed the microbiome sequences. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Israel Science Foundation (Research Project 1281/13).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequence data from root and rhizosphere samples was submitted to the Sequence Read Archive (SRA) of the NCBI databases under accession number SUB8180046, BioProject: PRJNA664890.

Acknowledgments

We acknowledge the computational help of Matan Cohen.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berg, G.; Smalla, K. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 2009, 68, 1–13. [Google Scholar] [CrossRef] [Green Version]
  2. İnceoğlu, Ö.; Falcão Salles, J.; van Elsas, J.D. Soil and Cultivar Type Shape the Bacterial Community in the Potato Rhizosphere. Microb. Ecol. 2012, 63, 460–470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Philippot, L.; Raaijmakers, J.M.; Lemanceau, P.; Van Der Putten, W.H. Going back to the roots: The microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 2013, 11, 789–799. [Google Scholar] [CrossRef] [PubMed]
  4. Bulgarelli, D.; Garrido-Oter, R.; Münch, P.C.; Weiman, A.; Dröge, J.; Pan, Y.; McHardy, A.C.; Schulze-Lefert, P. Structure and Function of the Bacterial Root Microbiota in Wild and Domesticated Barley Resource Structure and Function of the Bacterial Root Microbiota in Wild and Domesticated Barley. Cell Host Microbe 2015, 17, 392–403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Jones, D.L.; Nguyen, C.; Finlay, R.D. Carbon flow in the rhizosphere: Carbon trading at the soil–root interface. Plant Soil 2009, 321, 5–33. [Google Scholar] [CrossRef]
  6. Mendes, R.; Garbeva, P.; Raaijmakers, J.M. The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. Oxford Academic. FEMS Microbiol. Rev. 2013, 37, 634–663. [Google Scholar] [CrossRef] [PubMed]
  7. Bertin, C.; Yang, X.; Weston, L.A. The role of root exudates and allelochemicals in the rhizosphere. Plant Soil. 2003, 256, 67–83. [Google Scholar] [CrossRef]
  8. Bais, H.P.; Weir, T.L.; Perry, L.G.; Gilroy, S.; Vivanco, J.M. The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol. 2006, 57, 233–266. [Google Scholar] [CrossRef] [Green Version]
  9. Badri, D.V.; Vivanco, J.M. Regulation and function of root exudates. Plant Cell Environ. 2009, 32, 666–681. [Google Scholar] [CrossRef]
  10. Jaeger, C.H.; Lindow, S.E.; Miller, W.; Clark, E.; Firestone, M.K. Mapping of sugar and amino acid availability in soil around roots with bacterial sensors of sucrose and tryptophan. Appl. Environ. Microbiol. 1999, 65, 2685–2690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Doan, T.H.; Doan, T.A.; Kangas, M.; Ernest, A.E.; Tran, D.; Wilson, C.L.; Holmes, A.E.; Doyle, E.L.; Brooks, T.L.D. A low-cost imaging method for the temporal and spatial colorimetric detection of free amines on maize root surfaces. Front. Plant Sci. 2017, 8, 1513. [Google Scholar] [CrossRef] [Green Version]
  12. Sasse, J.; Martinoia, E.; Northen, T. Feed your friends: Do plant exudates shape the root microbiome? Trends Plant Sci. 2018, 23, 25–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rovira, A.D. Zones of exudation along plant roots and spatial distribution of micro-organisms in the rhizosphere. Pestic. Sci. 1973, 4, 361–366. [Google Scholar] [CrossRef]
  14. Sauer, D.; Kuzyakov, Y.; Stahr, K. Spatial distribution of root exudates of five plant species as assessed by14C labeling. J. Plant Nutr Soil Sci. 2006, 169, 360–362. [Google Scholar] [CrossRef]
  15. Vidal, A.; Hirte, J.; Bender, S.F.; Mayer, J.; Gattinger, A.; Höschen, C.; Schädler, S.; Iqbal, T.M.; Mueller, C.W. Linking 3d soil structure and plant-microbe-soil carbon transfer in the rhizosphere. Front. Environ. Sci. 2018, 6, 9. [Google Scholar] [CrossRef]
  16. Canarini, A.; Kaiser, C.; Merchant, A.; Richter, A.; Wanek, W. Root Exudation of Primary Metabolites: Mechanisms and Their Roles in Plant Responses to Environmental Stimuli. Front. Plant Sci. 2019, 10, 157. [Google Scholar] [CrossRef] [Green Version]
  17. Yan, Y.; Kuramae, E.E.; De Hollander, M.; Klinkhamer, P.G.L.; Van Veen, J.A. Functional traits dominate the diversity-related selection of bacterial communities in the rhizosphere. ISME J. 2017, 11, 56–66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Kamutando, C.N.; Vikram, S.; Kamgan-Nkuekam, G.; Makhalanyane, T.; Greve, M.; Le Roux, J.; Richardson, D.M.; Cowan, D.A.; Valverde, A. The functional potential of the rhizospheric microbiome of an invasive tree species. Acacia Dealbata. Microb. Ecol. 2019, 77, 191–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Badri, D.V.; Chaparro, J.M.; Zhang, R.; Shen, Q.; Vivanco, J.M. Application of natural blends of phytochemicals derived from the root exudates of arabidopsis to the soil reveal that phenolic-related compounds predominantly modulate the soil microbiome. J. Biol. Chem. 2013, 288, 4502–4512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Chaparro, J.M.; Badri, D.V.; Bakker, M.G.; Sugiyama, A.; Manter, D.K.; Vivanco, J.M. Root Exudation of Phytochemicals in Arabidopsis Follows Specific Patterns That Are Developmentally Programmed and Correlate with Soil Microbial Functions. PLoS ONE 2013, 8, e55731. [Google Scholar] [CrossRef]
  21. Pausch, J.; Kuzyakov, Y. Carbon input by roots into the soil: Quantification of rhizodeposition from root to ecosystem scale. Glob. Chang. Biol. 2018, 24, 1–12. [Google Scholar] [CrossRef]
  22. Na, X.; Cao, X.; Ma, C.; Ma, S.; Xu, P.; Liu, S.; Wang, J.; Wang, H.; Chen, L.; Qiao, Z. Plant stage, not drought stress, determines the effect of cultivars on bacterial community diversity in the rhizosphere of broomcorn millet (Panicum miliaceum L.). Front. Microbiol. 2019, 10, 828. [Google Scholar] [CrossRef] [PubMed]
  23. Edwards, J.A.; Santos-Medellín, C.M.; Liechty, Z.S.; Nguyen, B.; Lurie, E.; Eason, S.; Phillips, G.; Sundaresan, V. Compositional shifts in root-associated bacterial and archaeal microbiota track the plant life cycle in field-grown rice. PLoS Biol. 2018, 16, e2003862. [Google Scholar] [CrossRef] [Green Version]
  24. Haichar, F.E.Z.; Marol, C.; Berge, O.; Rangel-Castro, J.I.; Prosser, J.; Balesdent, J.; Heulin, T.; Achouak, W. Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2008, 2, 1221–1230. [Google Scholar] [CrossRef] [PubMed]
  25. Lomond, J.S.; Tong, A.Z. Rapid analysis of dissolved methane, ethylene, acetylene and ethane using partition coefficients and headspace-gas chromatography. J. Chromatogr. Sci. 2011, 49, 469–475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Li, T.; Li, Y.; Shi, Z.; Wang, S.; Wang, Z.; Liu, Y.; Wen, X.; Mo, F.; Han, J.; Liao, Y. Crop development has more influence on shaping rhizobacteria of wheat than tillage practice and crop rotation pattern in an arid agroecosystem. Appl. Soil Ecol. 2021, 165, 104016. [Google Scholar] [CrossRef]
  27. Betlach, M.R.; Tiedje, J.M. Kinetic explanation for accumulation of nitrite, nitric oxide, and nitrous oxide during bacterial denitrification. Appl. Environ. Microbiol. 1981, 42, 1074–1084. [Google Scholar] [CrossRef] [Green Version]
  28. Thomson, A.J.; Giannopoulos, G.; Pretty, J.; Baggs, E.M.; Richardson, D.J. Biological sources and sinks of nitrous oxide and strategies to mitigate emissions. Philos Trans. R. Soc. Lond. B Biol. Sci. 2012, 367, 1157–1168. [Google Scholar] [CrossRef] [Green Version]
  29. Bazin, M.J.; Markham, P.; Scott, E.M.; Lynch, J.M. Population dynamics and rhizosphere interactions. In The Rhizosphere; Lynch, J.M., Ed.; John Wiley and Sons Ltd.: Hoboken, NJ, USA, 1990; pp. 99–127. [Google Scholar]
  30. Reinhold-Hurek, B.; Bünger, W.; Burbano, C.S.; Sabale, M.; Hurek, T. Roots shaping their microbiome: Global hotspots for microbial activity. Annu. Rev. Phytopathol. 2015, 53, 403–424. [Google Scholar] [CrossRef]
  31. Kuzyakov, Y.; Razavi, B.S. Rhizosphere size and shape: Temporal dynamics and spatial stationarity. Soil Biol. Biochem. 2019, 135, 343–360. [Google Scholar] [CrossRef]
  32. Chaparro, J.M.; Badri, D.V.; Vivanco, J.M. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2014, 8, 790–803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. De Weert, S.; Vermeiren, H.; Mulders, I.H.M.; Kuiper, I.; Hendrickx, N.; Bloemberg, G.V.; Vanderleyden, J.; De Mot, R.; Lugtenberg, B.J.J. Flagella-driven chemotaxis towards exudate components is an important trait for tomato root colonization by Pseudomonas fluorescens. Mol. Plant Microbe Interact. 2002, 15, 1173–1180. [Google Scholar] [CrossRef] [Green Version]
  34. Lareen, A.; Burton, F.; Schäfer, P. Plant root-microbe communication in shaping root microbiomes. Plant. Mol. Biol. 2016, 90, 575–587. [Google Scholar] [CrossRef] [Green Version]
  35. Harshey, R.M. Bacterial motility on a surface: Many ways to a common goal. Annu. Rev. Microbiol. 2003, 57, 249–273. [Google Scholar] [CrossRef] [PubMed]
  36. Lugtenberg, B.J.J.; Dekkers, L.; Bloemberg, G.V. Molecular determinants of rhizosphere colonization by Pseudomonas. Annu. Rev. Phytopathol. 2001, 39, 461–490. [Google Scholar] [CrossRef] [PubMed]
  37. Khosla, C. Structures and mechanisms of polyketide synthases. American Chemical Society. J. Org. Chem. 2009, 74, 6416–6420. [Google Scholar] [CrossRef] [PubMed]
  38. Iannucci, A.; Fragasso, M.; Beleggia, R.; Nigro, F.; Papa, R. Evolution of the crop rhizosphere: Impact of domestication on root exudates in tetraploid wheat (Triticum turgidum L.). Front. Plant Sci. 2017, 8, 2124. [Google Scholar] [CrossRef] [Green Version]
  39. Hall, C.W.; Mah, T.F. Molecular mechanisms of biofilm-based antibiotic resistance and tolerance in pathogenic bacteria. FEMS Microbiol. Rev. 2017, 41, 276–301. [Google Scholar] [CrossRef]
  40. Shah, S.A.; Vestergaard, G.; Garrett, R.A. CRISPR/Cas and CRISPR/Cmr immune systems of archaea. In Regulatory RNAs in Prokaryotes; Springer: Wien, Austria, 2012; pp. 163–182. [Google Scholar]
  41. Allard-Massicotte, R.; Tessier, L.; Lécuyer, F.; Lakshmanan, V.; Lucier, J.-F.; Garneau, D.; Caudwell, L.; Vlamakis, H.; Bais, H.P.; Beauregard, P.B. Bacillus subtilis early colonization of Arabidopsis thaliana roots involves multiple chemotaxis receptors. MBio 2016, 7, e01664-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Oburger, E.; Jones, D.L. Sampling root exudates—Mission impossible? Rhizosphere 2018, 6, 116–133. [Google Scholar] [CrossRef]
  43. Donn, S.; Kirkegaard, J.A.; Perera, G.; Richardson, A.E.; Watt, M. Evolution of bacterial communities in the wheat crop rhizosphere. Environ. Microbiol. 2015, 17, 610–621. [Google Scholar] [CrossRef]
  44. Mahoney, A.K.; Yin, C.; Hulbert, S.H. Community structure, species variation, and potential functions of rhizosphere-associated bacteria of different winter wheat (Triticum aestivum) cultivars. Front. Plant. Sci. 2017, 8, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Chen, S.; Waghmode, T.R.; Sun, R.; Kuramae, E.E.; Hu, C.; Liu, B. Root-associated microbiomes of wheat under the combined effect of plant development and nitrogen fertilization. Microbiome 2019, 7, 136. [Google Scholar] [CrossRef] [Green Version]
  46. Kirby, R. Actinomycetes and Lignin Degradation. Adv. Appl. Microbiol. 2005, 58, 125–168. [Google Scholar]
  47. Watt, M.; Hugenholtz, P.; White, R.; Vinall, K. Numbers and locations of native bacteria on field-grown wheat roots quantified by fluorescence in situ hybridization (FISH). Environ. Microbiol. 2006, 8, 871–884. [Google Scholar] [CrossRef] [PubMed]
  48. Karnaouri, A.; Matsakas, L.; Krikigianni, E.; Rova, U.; Christakopoulos, P. Valorization of waste forest biomass toward the production of cello-oligosaccharides with potential prebiotic activity by utilizing customized enzyme cocktails. Biotechnol. Biofuels 2019, 12, 1–19. [Google Scholar] [CrossRef] [Green Version]
  49. Usyskin-Tonne, A.; Hadar, Y.; Yermiyahu, U.; Minz, D. Elevated CO2 has a significant impact on denitrifying bacterial community in wheat roots. Soil Biol. Biochem. 2020, 142, 107697. [Google Scholar] [CrossRef]
  50. Rocha, J.; Cacace, D.; Kampouris, I.; Guilloteau, H.; Jäger, T.; Marano, R.; Karaolia, P.; Manaia, C.M.; Merlin, C.; Fatta-Kassinos, D.; et al. Inter-laboratory calibration of quantitative analyses of antibiotic resistance genes. J. Environ. Chem. Eng. 2020, 8, 102214. [Google Scholar] [CrossRef]
  51. Tovi, N.; Frenk, S.; Hadar, Y.; Minz, D. Host Specificity and Spatial Distribution Preference of Three Pseudomonas Isolates. Front. Microbiol. 2019, 9, 3263. [Google Scholar] [CrossRef]
  52. Usyskin-Tonne, A.; Hadar, Y.; Yermiyahu, U.; Minz, D. Elevated CO2 and nitrate levels increase wheat root-associated bacterial abundance and impact rhizosphere microbial community composition and function. ISME J. 2020, 15, 1073–1084. [Google Scholar] [CrossRef]
  53. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 15 September 2021).
  54. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hyatt, D.; Chen, G.-L.; LoCascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010, 11, 119. [Google Scholar] [CrossRef] [Green Version]
  58. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef]
  59. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  60. Huson, D.H.; Beier, S.; Flade, I.; Górska, A.; El-Hadidi, M.; Mitra, S.; Ruscheweyh, H.J.; Tappu, R. MEGAN community edition—Interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comput. Biol. 2016, 12, e1004957. [Google Scholar] [CrossRef] [Green Version]
  61. Clarke, K.; Ainsworth, M. A method of linking multivariate community structure to environmental variables. Mar. Ecol. Prog Ser. 1993, 92, 205–219. [Google Scholar] [CrossRef]
  62. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef]
  63. Dixon, P. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 2003, 14, 927–930. [Google Scholar] [CrossRef]
  64. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Austral. Ecol. 2001, 26, 32–46. [Google Scholar]
  65. Shannon, P. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  66. Faust, K.; Raes, J. CoNet app: Inference of biological association networks using Cytoscape. F1000Research 2016, 5, 1519. [Google Scholar] [CrossRef] [PubMed]
  67. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Shifts in root-associated and rhizosphere bacterial community structure and function with plant age. nMDS ordination plots show clustering patterns of wheat root-associated and rhizosphere microbial community structures (A) and functions (B) during wheat growth. Data matrix was transformed using normalized count transformation in the DESeq2 package, and then ordination was generated using Bray-Curtis dissimilarity and Euclidean distances, respectively. (C) Changes in structural and (D) functional similarity between wheat root-associated and rhizosphere microbial communities at each week. Different lowercase letters indicate significant difference (p < 0.05) by Student’s t-test. Weeks 1, 2, 5, 9 and 12: n = 6; week 17 green (green spike) and yellow (yellow spike): n = 2 for each. A black dot in Figure 1C, indicate the outlier at week 12.
Figure 1. Shifts in root-associated and rhizosphere bacterial community structure and function with plant age. nMDS ordination plots show clustering patterns of wheat root-associated and rhizosphere microbial community structures (A) and functions (B) during wheat growth. Data matrix was transformed using normalized count transformation in the DESeq2 package, and then ordination was generated using Bray-Curtis dissimilarity and Euclidean distances, respectively. (C) Changes in structural and (D) functional similarity between wheat root-associated and rhizosphere microbial communities at each week. Different lowercase letters indicate significant difference (p < 0.05) by Student’s t-test. Weeks 1, 2, 5, 9 and 12: n = 6; week 17 green (green spike) and yellow (yellow spike): n = 2 for each. A black dot in Figure 1C, indicate the outlier at week 12.
Ijms 22 11948 g001
Figure 2. Influence of wheat developmental stage on total bacterial community and nosZ gene abundance in roots and rhizosphere. (A) Changes in abundance of total bacterial community (16S rRNA) and its relative abundance per root tissue (e.g., per the plant housekeeping gene TEF, and per total plant genes from metagenome) during wheat growth. (B) Changes in relative abundance of denitrifying gene nosZ per total bacterial community in roots and rhizosphere as well as nitrate content in rhizosphere during wheat growth. Each bar represents the average of three biological each of three technical replicates with standard error. Different uppercase and lowercase letters indicate significant difference (p ≤ 0.05) by Student’s t test. N-NO3 = nitrate nitrogen.
Figure 2. Influence of wheat developmental stage on total bacterial community and nosZ gene abundance in roots and rhizosphere. (A) Changes in abundance of total bacterial community (16S rRNA) and its relative abundance per root tissue (e.g., per the plant housekeeping gene TEF, and per total plant genes from metagenome) during wheat growth. (B) Changes in relative abundance of denitrifying gene nosZ per total bacterial community in roots and rhizosphere as well as nitrate content in rhizosphere during wheat growth. Each bar represents the average of three biological each of three technical replicates with standard error. Different uppercase and lowercase letters indicate significant difference (p ≤ 0.05) by Student’s t test. N-NO3 = nitrate nitrogen.
Ijms 22 11948 g002
Figure 3. Network co-occurrence analysis of bacterial communities of root and rhizosphere samples constructed using CoNet, a Cytoscape plugin. Green and red links (i.e., edges) indicate significant positive and negative correlations between two nodes, respectively. In (A,B), each node represents taxa affiliated to all taxonomic levels from the metagenome, labeled at the phylum level. In (C,D), each node represents a group of functional genes, annotated with known KO groups from the KEGG database.
Figure 3. Network co-occurrence analysis of bacterial communities of root and rhizosphere samples constructed using CoNet, a Cytoscape plugin. Green and red links (i.e., edges) indicate significant positive and negative correlations between two nodes, respectively. In (A,B), each node represents taxa affiliated to all taxonomic levels from the metagenome, labeled at the phylum level. In (C,D), each node represents a group of functional genes, annotated with known KO groups from the KEGG database.
Ijms 22 11948 g003
Figure 4. Changes in root-associated and rhizosphere bacterial community structure as a function of niche (roots vs. rhizosphere) and plant age. (A) Shifts in root and rhizosphere bacterial community structure at the phylum and order levels as influenced by plant age. Significantly abundant groups in root or rhizosphere (niche) (B) and with spike formation (C). Changes in bacterial abundance were calculated using DESeq2 with cutoff FDR-adjusted p-value < 0.05. Before spike formation, n = 9; after spike formation, n = 10. Asterisk (*) indicate bacteria groups at order level, that their abundance in root and rhizosphere microbiome was influence by spike formation.
Figure 4. Changes in root-associated and rhizosphere bacterial community structure as a function of niche (roots vs. rhizosphere) and plant age. (A) Shifts in root and rhizosphere bacterial community structure at the phylum and order levels as influenced by plant age. Significantly abundant groups in root or rhizosphere (niche) (B) and with spike formation (C). Changes in bacterial abundance were calculated using DESeq2 with cutoff FDR-adjusted p-value < 0.05. Before spike formation, n = 9; after spike formation, n = 10. Asterisk (*) indicate bacteria groups at order level, that their abundance in root and rhizosphere microbiome was influence by spike formation.
Ijms 22 11948 g004
Figure 5. Spatial effect on microbiome functions. Significantly abundant functional genes in two niches (root-associated vs. rhizosphere) summarized in functional groups (A) and genes encoding for those functional groups (B).
Figure 5. Spatial effect on microbiome functions. Significantly abundant functional genes in two niches (root-associated vs. rhizosphere) summarized in functional groups (A) and genes encoding for those functional groups (B).
Ijms 22 11948 g005
Figure 6. Temporal effect on microbiome functions. Significantly abundant functional genes before and after spike formation, divided into groups (A) and in roots (B) and rhizosphere (C).
Figure 6. Temporal effect on microbiome functions. Significantly abundant functional genes before and after spike formation, divided into groups (A) and in roots (B) and rhizosphere (C).
Ijms 22 11948 g006
Figure 7. Schematic model presenting significantly abundant functional genes. (A) Schematic representation of significantly abundant functional genes before and after spike formation in roots and rhizosphere. Size of each rectangle is proportional to sum of significantly enriched functional gene counts in each subcategory, calculated based on data in Figure 6B,C. (B) Schematic representation of significantly abundant functional genes before and after spike formation in roots in the three selected, dominant bacterial phyla proteobacteria, Actinobacteria and Bacteroidetes. Size of each rectangle is proportional to normalized amount of the gene functions (using normalized count transformation in DESeq2 package), divided by relative abundance of each relevant taxa. List of all normalized counts is provided in Supplemental Table S3.
Figure 7. Schematic model presenting significantly abundant functional genes. (A) Schematic representation of significantly abundant functional genes before and after spike formation in roots and rhizosphere. Size of each rectangle is proportional to sum of significantly enriched functional gene counts in each subcategory, calculated based on data in Figure 6B,C. (B) Schematic representation of significantly abundant functional genes before and after spike formation in roots in the three selected, dominant bacterial phyla proteobacteria, Actinobacteria and Bacteroidetes. Size of each rectangle is proportional to normalized amount of the gene functions (using normalized count transformation in DESeq2 package), divided by relative abundance of each relevant taxa. List of all normalized counts is provided in Supplemental Table S3.
Ijms 22 11948 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Usyskin-Tonne, A.; Hadar, Y.; Minz, D. Spike Formation Is a Turning Point Determining Wheat Root Microbiome Abundance, Structures and Functions. Int. J. Mol. Sci. 2021, 22, 11948. https://doi.org/10.3390/ijms222111948

AMA Style

Usyskin-Tonne A, Hadar Y, Minz D. Spike Formation Is a Turning Point Determining Wheat Root Microbiome Abundance, Structures and Functions. International Journal of Molecular Sciences. 2021; 22(21):11948. https://doi.org/10.3390/ijms222111948

Chicago/Turabian Style

Usyskin-Tonne, Alla, Yitzhak Hadar, and Dror Minz. 2021. "Spike Formation Is a Turning Point Determining Wheat Root Microbiome Abundance, Structures and Functions" International Journal of Molecular Sciences 22, no. 21: 11948. https://doi.org/10.3390/ijms222111948

APA Style

Usyskin-Tonne, A., Hadar, Y., & Minz, D. (2021). Spike Formation Is a Turning Point Determining Wheat Root Microbiome Abundance, Structures and Functions. International Journal of Molecular Sciences, 22(21), 11948. https://doi.org/10.3390/ijms222111948

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

Article Metrics

Back to TopTop