Li Jing
Ai Jia-min
Liu Xiao-dong
Deng Zhen-shan- College of Life Sciences, Yan’an University, Yan’an, China
In addition to the rhizobia, other non-rhizobial endophytes (NREs) have been simultaneously isolated from the root nodules. The existence of NREs in leguminous root nodules is a universal phenomenon, and they have the potential to enhance legume survival, especially under conditions of environmental stress. However, the diversity and biogeographic patterns of microbial communities inhabiting root nodules are not well studied or understood. Here, we explored and characterized the diversity of NRE bacteria by using 16S rRNA gene high-throughput amplicon sequencing. Additionally, we compared the biogeography and co-occurrence patterns in review of the bacterial microbiota inhabiting the rhizosphere, the bulk soil and the root nodule bacterial communities associated with Sophora davidii, a native N-fixing wild leguminous shrub in hilly and gully regions of the Loess Plateau of China. The results showed the presence of a large diversity of bacteria belonging to 81 phyla, 154 classes, 333 orders, 463 families, and 732 genera inside the nodules. Proteobacteria were dominant in the nodule and rhizosphere soil samples, and Actinomycetes were dominant in the bulk soil samples. Mesorhizobium was the dominant genus in the nodules, accounting for between 60.15 and 83.74% of the bacteria. The microbial community composition of the NRE in the root nodules differed from that in the rhizosphere soil and the bulk soil of S. davidii. Moreover, we found that the biogeographic patterns and assembly process of the rhizobia and non-rhizobia communities differed in the root nodule, the rhizosphere soil and the bulk soil. Furthermore, the correlation analysis between the soil’s physical and chemical properties and the bacteria showed that available phosphorus was the predominant factor affecting the bacterial diversity within the rhizosphere soil. Finally, our results revealed that the microbial network diagram of co-occurrence patterns showed more complexes in the soil than in the root nodules. This indicates that only specific microorganisms could colonize and thrive in the rhizosphere through the selection and filtering effects of roots. In conclusion, there are significant differences in bacterial community composition in the nodules, rhizosphere and bulk soil in the hilly and gully region of the Loess Plateau, which is the result of environmental filtration. Our study improves the understanding of the biogeographic patterns and diversity of bacterial microbiota inhabiting root nodules and can help quantify and define the root nodule assemblage process of S. davidii.
Introduction
One of the unique characteristics of legumes is that they can form special root or stem nodules on plant roots. The legume root nodule is a unique niche for symbiotic nitrogen-fixing occupied by bacteria so called rhizobia. In recent years, many studies on legumes have found that root nodules are a complex niche, in addition to rhizobium, there have also been non-rhizobial root nodule endophytes (NRE) isolated from legume nodules without clear indications of their role within the host, and an absence of positive nodulation tests (De Meyer et al., 2015). The reported NREs include Agrobacterium (Beiierinck and van Delden, 1902; Cummings et al., 2009), Bacillus (Rajendran et al., 2008; Palaniappan et al., 2010), Pseudomonas (Aserse et al., 2013), Enterobacter (Benhizia et al., 2004; Ibáñez et al., 2008), Paenibacillus (Zakhia et al., 2006; Deng et al., 2011; Lai et al., 2015), Burkholderia (Raupach and Kloepper, 1998; Rasolomampianina et al., 2005; Diouf et al., 2007), Staphylococcus (Deng et al., 2011), Pantoea (Li et al., 2007; Aserse et al., 2013), Methylobacterium (Sy et al., 2001; Palaniappan et al., 2010), Herbaspirillum (Valverde et al., 2003; Weiss et al., 2012), Microvirga (Ardley et al., 2012), Paraburkholderia (Sawana et al., 2014), Rhodococcus (Ampomah and Huss-Danell, 2011), Variovorax (Hoque et al., 2011; Aserse et al., 2013), Klebsiella (Pandya et al., 2013), Arthrobacter (Barnawal et al., 2014), Brevibacterium (Xu et al., 2014; Gopalakrishnan et al., 2016), and Streptomyces (Schrey and Tarkka, 2008; Sreevidya et al., 2016). Although different bacteria have been isolated from root nodules of various leguminous plants, there are still many bacteria that cannot be obtained by pure culture, therefore, the biological significance of these associations is largely unknown.
Although many of the NREs are not capable of nitrogen fixation, they have the potential to enhance legume survival, especially under conditions of environmental stress (Martínez-Hidalgo and Hirsch, 2017). And they are often detected within nodules obtained from soil.
These bacteria are opportunistic since they can infect nodules when rhizobia induce nodule formation (Ibáñez et al., 2008; Leite et al., 2017; Sarah and Djamel, 2021), and more than one species of non-rhizobia could be found in the nodules of Trifolium repens (Olenska et al., 2020). As many of the taxa identified in nodules are common in soil, it is possible that many of them are inadvertently entrapped during the rhizobia nodulation process. This provides evidence that both the soil bacterial community and the host plant might influence the nodule microbiome. Zgadzaj et al. (2015, 2016) concluded that endophyte nodule occupancy is host-controlled and that the ability of legumes to select a broad taxonomic range of root-associated bacteria along with rhizobia is likely to contribute to plant growth and ecological performance. However, Leite et al. (2017) reported that the non-rhizobial bacterial community in cowpea nodules was influenced by the soil type rather than the plant genotype.
However, there is currently limited information regarding the entry mode of endophytic bacteria, and the biological significance of these non-rhizobial bacteria is largely unknown. In addition to the plants themselves, rhizosphere soils are another niche of root-affected soil, and they represent a special area affected by plant roots (Li et al., 2021). The microbial community structure within the rhizosphere soil also has an important relationship with plant growth and development. Therefore, the simultaneous study of the diversity of endophytic bacteria and rhizosphere bacteria is of great significance for elucidating the function of these bacteria and tapping these biological resources (Tian and Zhang, 2017). Biogeographic patterns of microorganisms play important roles in generating and maintaining biodiversity, and community assembly processes are usually studied simultaneously with microbial biogeography (Zhang K. et al., 2017; Liu et al., 2021). Co-occurrence model analysis was applied to explore the potential interactions within niche spaces shared by community members (Kara et al., 2013).
Sophora davidii (Franch.) Skeels is a perennial shrub of the leguminous family, it grows in river valley dunes and bushes on hillside roads below 2,500 m. It is a native plant in China and is widely distributed in Guizhou, Ningxia, Shaanxi, Yunnan, and Sichuan provinces (Tai et al., 2011). As a deep-rooted nitrogen-fixing wild leguminous shrub, S. davidii has an important ecological value in plant community succession, soil improvement and soil erosion control in arid areas due to their strong adaptability, strong stress resistance, and developed root system (Wu et al., 2007). Moreover, compared with herbaceous plants, nitrogen-fixing shrubs can significantly increase the contents of available nitrogen and phosphorus in soil (Song et al., 2010). At present, research on S. davidii is mainly focused on the analysis of their chemical constituents (Ma et al., 2021; Zhang et al., 2021). A recent study focused on the diversity and distribution of rhizobia associated with S. davidii (Cao et al., 2021). However, studies on non-rhizobial composition and assembly in root nodules have not been reported, and little is known about the biodiversity of endophytic bacteria residing within the roots of S. davidii.
In this study, we compared the biogeography and co-occurrence patterns of the bacterial microbiota inhabiting the rhizosphere, the bulk soil and the root nodule bacterial communities of S. davidii collected from the hilly and gully regions of the Loess Plateau of China. The aim was to address the following questions: (i) clarify the diversity and community structure of rhizobia and the NREs of S. davidii and to estimate the ecological drivers for the bacterial communities in the rhizosphere, the bulk soil and the root nodule; (ii) clarify the distinct bacterial microbial co-occurrence patterns of the three compartments; and (iii) identify the biogeographic patterns and assembly process of the rhizobia and non-rhizobia communities that vary among the compartments.
Materials and methods
Sampling sites and nodule collection
Root nodules of the S. davidii at flowering-fruiting stages were collected from May to July 2021 at six sample sites in the hilly and gully regions of the Loess Plateau of China. The geographical locations of the sampling sites are shown in Figure 1 (the corresponding geographical features are shown in Supplementary Table 1). The distance between sampling sites was at least 10 km. Each sampling site included three subsites separated by more than 1 km. Fifteen plants with a spacing of at least 30 m were randomly selected and uprooted from each subsite. Three to four undamaged and healthy root nodules of similar size were excised from the lateral roots of each plant. These nodules were brushed away from the soil debris and then immediately desiccated over silica gel in 1.5 mL microtubes (Vincent, 1970). At the same time the root nodules were collected, the rhizosphere soil and bulk soil were collected. Rhizosphere soil is collected by digging the entire plant root out of the soil with minimum injury to the roots, gently shaking the roots until the soil that is not tightly adhering to the roots is removed, and then vigorously shaking the roots to collect the soil tightly adhering to the roots into sterile bags. The soil outside the rooting area is collected as bulk soil (Wang and Zabowski, 1998). One part of the soil samples were air-dried in the shade. The soils samples physicochemical characteristics; organic matter (SOM), total nitrogen (TN), ammonium nitrogen, nitrate-nitrogen (NN), alkali hydrolysable nitrogen, available phosphorus (AP), available potassium (AK), pH, and electrical conductivity (EC) were measured using the protocols described by the United States Department of Agriculture [USDA] (1996). The remainder was stored at –80°C for high-throughput sequencing at Novogene Biotechnology Co., Ltd. (Tianjin, China).
Figure 1. Geographical distribution of the six sampling sites of S. davidii in China. The simplified map shows six sampling areas of Huangling (HL), Fuxian (FX), Baota (BT), Yanchang (YC), Zhidan (ZD), and Zichang (ZC) in Yan’an city, Shaanxi Province. Each sampling area includes three sampling points.
To eliminate other microbial interference on the nodules’ surface, a surface sterilization procedure was performed on the nodules. The healthy and unruptured root nodules were rinsed under running water, wiped dry with filter paper, surface sterilized by immersion in 95% ethanol for 30 s, rinsed with sterile distilled water three times, then immersed in 5% sodium hypochlorite for 3 min, and rinsed with sterile distilled water eight times. To verify the effect of surface sterilization, the sterile distilled water washed for the last time was added to potato dextrose agar plates. If there was no colony growth on the plate, it indicated that the surface disinfection of the root nodules was effective. All samples were stored at –80°C until DNA extraction.
DNA extraction and sequencing
Approximately 2 g of disinfected root nodules was put into a precooled mortar, and then liquid nitrogen was added and thoroughly ground into a fine powder. Total DNA was extracted using a high-efficiency Plant Genomic DNA Extraction Kit (Tiangen, DP350, China). Soil samples (0.2 g) were weighed, and a soil genomic DNA extraction kit (Tiangen, DP336, China) was used to extract the soil DNA. The DNA concentration and purity were detected and evaluated by a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, United States) and 1% agarose gel electrophoresis. The samples were diluted with sterile distilled water to 1 ng/μL. Primers 16S-341F (5′-CCTAYGGGRBGCASCAG-3′) and 16S-806R (5′-GGACTACNNGGGTATCTAAT-3′) were used for the PCR amplification of the 16S rRNA V3+V4 region. PCRs were carried out with 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM forward and reverse primers and approximately 10 ng of template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 30 s. Finally, the samples were incubated at 72°C for 5 min. The same volume of 1 × loading buffer (containing SYBR green) was mixed with the PCR products, and electrophoresis was performed on a 2% agarose gel for detection. The PCR products were mixed in equidensity ratios. Then, the PCR products were purified with a Qiagen Gel Extraction Kit (Qiagen, Germany).
Sequencing libraries were generated using a TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, United States). The library quality was assessed on the Qubit@2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on an Illumina NovaSeq 6000 platform, and 250 bp paired-end reads were generated (Zhang C. S. et al., 2017).
Microbial community analysis
The paired-end reads were assigned to samples based on their unique barcode and truncated by cutting off the barcode and primer sequence. Then, the paired-end reads were merged using FLASH (VI.2.7)1 (Magoc and Salzberg, 2011), a very fast and accurate analysis tool that merges paired-end reads when some of the reads overlap the read generated from the opposite end of the same DNA fragment. The splicing sequences were assigned raw tags. The quality filtering of the raw tags was performed under specific filtering conditions to obtain a high-quality clean tag (Bokulich et al., 2013) according to the QIIME (Version 1.9.1)2 (Caporaso et al., 2010) quality control process. The tags were compared with the reference database (Silva database)3 using the UCHIME algorithm (UCHIME Algorithm)4 (Edgar et al., 2011) to detect chimera sequences, and then the chimera sequences were removed (Haas et al., 2011). Then, the Effective Tags were finally obtained.
Sequence analysis was performed by UPARSE software (UPARSE v7.0.1001)5 (Edgar, 2013). The operational taxonomic units (OTUs) were clustered according to 97% similarity (Edgar, 2010) and assigned to each sample with the QIIME pipeline for sequence analysis (Caporaso et al., 2010).
Representative sequences for each OTU were screened for further annotation. For each representative sequence, the Silva Database (see text footnote 3; Quast et al., 2013) was used based on the Mothur algorithm to annotate taxonomic information. To study the phylogenetic relationship of different OTUs and the difference in the dominant species in different samples, multiple sequence alignments were conducted using MUSCLE software (Version 3.8.31)6 (Edgar, 2004).
Network construction
To determine the co-occurrence network pattern of the bacterial communities in root chambers and the structural differences of bacterial communities in different niches, Spearman correlation analysis was carried out among the most abundant bacterial genera in each root chamber. The data with a correlation coefficient greater than 0.6 and a p-value < 0.5 (Barberán et al., 2012) were selected to form a correlation network through pairwise comparison of genus abundance, in which each node represents a genus and each edge represents a significant correlation between the two nodes. The number of nodes and edges, average path length (APL), network diameter, average degree, graph density, clustering coefficient (CC), and modularity were calculated by igraph packages in R language, and the results were visualized by Gephi software (Bastian et al., 2009).
Statistical analysis
One-way ANOVAs and the Tukey–Kramer test were used to test the alpha diversity of the samples (including observed species, Chao l, Shannon, Simpson, ACE and Good coverage). Beta diversity [including principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA)] was used to analyze and evaluate the differences in species complexity in the samples. A least discriminant analysis (LDA) effect size (LEfSe) taxonomic cladogram was used to identify bacterial biomarkers (LDA > 4) in the different treatments using LEfSe software (Segata et al., 2011). Alpha diversity and beta diversity were calculated by QIIME and displayed with R software (Caporaso et al., 2010). The correlation between the soil’s physicochemical properties and the bacterial community was calculated and plotted by R.
Results
Physicochemical characteristics of the soil
The physicochemical characteristics of the rhizosphere soil of S. davidii are shown in Table 1. In general, the soil of each sampling site was alkaline, and the pH was between 7.85 and 8.26. The contents of AP and TN in the sample ZCS were 78.87 and 53.93 mg/kg, respectively, which were significantly higher than those of other samples. The content of AK in the sample FXS was the highest at 249.33 mg/kg. The range of soil organic matter was between 5.25 and 25.24 g/kg. The sample HLS had an organic matter content of 25.24 g/kg, which was significantly higher than that of the other samples. The contents of AP, AK and organic matter in the rhizosphere soil were generally higher than those in the bulk soil. The contents of NN and alkaline hydrolyzable nitrogen in the rhizosphere soil at different sample sites were lower than those in the corresponding bulk soil, but the content of NN was the inverse. In the Huangling and Fuxian samples, the content of TN in the rhizosphere soil was lower than that in the bulk soil, but the results were inverse in the other four sample groups. The conductivity of the sample ZDS, which was 681.67 was much higher than that of the other samples.
Table 1. Soil physicochemical property of each sample.
Microbial composition and diversity in the root nodule, the rhizosphere, and the bulk soil
Illumina MiSeq sequencing generated a total of 4,158,024 raw tags representing 54 samples, with individual reads ranging from 55,931 to 89,355 for each sample. After quality control, a total of 4,130,398 clean tags remained. After removing the chimeras, 3,217,738 effective tags (with an average length of 405–419 bp) were obtained for the generation of OTUs. The values of Q20 ranged from 97.92 to 98.5, indicating that the quality of the databases was high (Table 2). The trend of rarefaction curves suggested that the sequencing data were sufficient for each sample to represent the bacterial community (Figure 2). A total of 114,808 OTUs were detected.
Table 2. Quantitative statistics of bacterial OTUs and tags of each sample.
Figure 2. Rarefaction curves based on the sequences of the V3+V4 region of the 16S rRNA gene from samples associated with S. davidii. HLN, FXN, BTN, YCN, ZDN, and ZCN indicate root nodule samples from the six sampling sites. HLS, FXS, BTS, YCS, ZDS, and ZCS indicate the rhizosphere soil samples of the six sampling sites. HLCK, FXCK, BTCK, YCCK, ZDCK, and ZCCK indicate the bulk soil samples of the six sampling sites. N = 3.
High-throughput sequencing revealed bacterial community diversity at the phylum level in different samples (Figure 3A). A total of 81 phyla were identified, and the relative abundances of the first 10 phyla are shown in Figure 3A. Proteobacteria had an absolute advantage in root nodules, accounting for 85.77%. Proteobacteria, Actinobacteria and Acidobacteria were the dominant phyla in the rhizosphere soil, accounting for 28.93, 15.41, and 12.71%, respectively. Actinomycetes and Proteobacteria were the dominant phyla in the bulk soil samples, accounting for 32.61 and 21.93%, respectively. The proportion of Acidobacteria in the rhizosphere soil (12.71%) was significantly higher than that in the nodule soil (0.74%) and the bulk soil (3.76%). Compared with the rhizosphere soil and the root nodules, Bacteroidetes were more dominant in the bulk soil (8.30%). However, there was a certain proportion of unrecognized sequences in all of the samples (same below).
Figure 3. The relative abundance of endophytic bacteria in the three groups of samples at different taxonomic levels. (A) is the relative abundance of bacteria at the phylum level. (B) is the relative abundance of bacteria at the class level. (C) is the relative abundance of bacteria at the genus level.
At the class level (Figure 3B), Alphaproteobacteria was the most dominant in the root nodules, accounting for 81.13%. In addition, Actinobacteria, Gammaproteobacteria, Bacteroidia, and Bacilli accounted for 3.57, 4.63, 1.5, and 2.06%, respectively, and were the dominant class of bacteria (≥1%). In the rhizosphere soil samples, Alphaproteobacteria, Gammaproteobacteria, and Actinobacteria accounted for a large proportion, at 15.69, 13.19, and 8.88%, respectively. In the bulk soil, Actinobacteria (25.27%), Alphaproteobacteria (12.29%), Gammaproteobacteria (9.64%), Bacteroidia (8.25%), Acidimicrobia (6.19%), Thermoleophilia (3.41%), Chloroflexia (1.56%), and Vicinamibacteria (1.4%) were the dominant bacteria. Overall, the content of Alphaproteobacteria in the root nodules was significantly higher than that in the rhizosphere soil and the bulk soil, while the content of Actinobacteria and Gammaproteobacteria was significantly lower than that in the rhizosphere soil and the bulk soil. It can be seen that different plant-related tissues are selective for the microorganisms present in them.
A total of 732 different genera were detected at the genus level. There were significant differences in the overall composition of bacteria in the different samples. The most prominent phenomenon was that Mesorhizobium was overwhelmingly dominant in the root nodules (75.07%) compared with the rhizosphere soil and the bulk soil. In addition to the 10 genera listed in Figure 3C, other strains accounted for 92.45% in the rhizosphere soil and 78.15% in the bulk soil. The results showed that the bacterial communities in the root nodules and rhizosphere soil were diverse, but the bacterial species richness in the root nodules was lower than that in the soil samples.
In addition, we compared the corresponding number of OTUs in different grouping samples (Figure 4). The number of specific OTUs in nodule and rhizosphere soil samples was 581 and 4,427, respectively, and there were 2,759 common OTUs in total (Figure 4A). The special OTUs of the rhizosphere soil and the bulk soil samples were 2,371 and 1,353, respectively, and there were 4,815 common OTUs in total (Figure 4B).
Figure 4. Venn diagrams show the number of shared and unique OTUs in different nodules and soil samples. (A) is shared and unique OTUs between nodules and rhizosphere soils. (B) is shared and unique OTUs between rhizosphere soil and bulk soil.
The alpha diversity index of the three groups of samples is shown in Figure 5. The observed species index of the rhizosphere soil sample (2,383) was significantly higher than that of the bulk soil (1,544) and nodule sample (576; Figure 5A). The ACE, Chao1 and PD whole tree indices of the rhizosphere soil were 3,284.316, 3,152.988, and 204.567, respectively (Figures 5B–D). The ACE, Chao1, and PD whole tree indices of the root nodules were 867.9, 794.718, and 74.77 (Figures 5B–D), respectively. The ACE, Chao1 and PD whole tree indices in the bulk soil were 2,245.706, 2,152.737, and 151.514 (Figures 5B–D), respectively. The Shannon index and Simpson index of the rhizosphere soil were 9.533 and 0.995, respectively (Figures 5E,F). The Shannon and Simpson indices in the bulk soil were 7.498 and 0.954, respectively (Figures 5E,F). In the root nodules, Shannon and Simpson indices were only 2.434 and 0.405, respectively (Figures 5E,F). The above data showed that the community diversity and abundance within the rhizosphere soil were significantly higher than those of the bulk soil and the root nodules. The alpha diversity index indicates that the rhizosphere can recruit a variety of different microorganisms at the same time and has an enrichment effect on them. Nodules, on the other hand, act as filters, screening out a group of microbes that are beneficial to both themselves and their host plants.
Figure 5. Alpha diversity index of the microbial communities within the three types of S. davidii samples. (A–F) represent observed species, and the results of the ACE, Chao1, PD Whole Tree, the Shannon and Simpson index of the different samples of nodules, rhizosphere soil and bulk soil. Statistical analyses were performed by the paired Wilcoxon rank-sum test and Tukey’s test. The significance is denoted by asterisks where ⋅P < 0.05.
The rhizosphere microbial community is closely related to the composition of endophytic bacteria within the root nodules
Principal coordinate analysis, the UPGMA, and LDA effect size (LEfSe) analysis were used for beta diversity analysis to compare the microbial compositions of the different samples. The UPGMA tree showed that all the samples were divided into two different clusters: Group 1 was composed of root nodule samples, and Group 2 was composed of the rhizosphere soil and the bulk soil samples (Figure 6A). This indicates that the microflora in the root nodule is different from that in the rhizosphere soil. The PCoA results based on unweighted UniFrac distance showed that there were differences in microflora among different samples, which were 18.31% (PC1) and 7.78% (PC2), respectively (Figure 6B). The community structure of the rhizosphere soil was similar to that of the bulk soil, while there were significant differences in the bacterial community structure between the root nodules and the soil, indicating that the root nodules and the soil had their own unique bacterial community structure.
Figure 6. Beta diversity analysis. (A) is cluster pedigrees of different flora within the different samples based on unweighted UniFrac distances. (B) is principal coordinate analysis (PCoA) based on the relative abundance of bacterial OTUs.
Through LEfSe analysis, we obtained species with significant differences in abundance in the different types of samples (Figure 7). At the family level (Figure 7A), Rhizobiaceae was significantly enriched within the root nodule samples. In the rhizosphere soil samples, Vicinamibacteraceae, Pyrinomonadaceae, Gemmatimonadaceae, and Nitrosomonadaceae were significantly enriched. Hymenobacteraceae, Chitinophagaceae, Nocardioidaceae, Micrococcaceae, Oxalobacteraceae, and Sphingomonadaceae were more abundant in the bulk soils. To assess the significance of discriminatively specific bacterial taxa, LEfSe analysis was accomplished from the phylum to genus across the three sites using the default logarithmic (LDA) value of 4 (Figure 7B). The clade graph showed that 24 bacterial groups were enriched in the bulk soil sample, and 21 were enriched in the rhizosphere soil sample, followed by the nodule sample (5 bacterial groups).
Figure 7. Results of the LEfSe analysis. (A) is phylogenetic cladistic diagram of the LEfSe. The circle of radiation from the center to the periphery represents the level of classification from the phylum to the species. Each small circle at different classification levels represents a classification at that level, and the diameter of the small circle is proportional to the relative abundance. (B) is histogram of the LDA value distribution (LDA > 4).
To further demonstrate the visual clustering of bacterial communities in the above analysis, different plant compartments were examined using ANOSIM analysis (Table 3). At the OTU level, the bacterial microbiota of all plant compartments differed significantly from one another.
Table 3. Analysis of similarity.
Bacterial communities associated with the soil’s physicochemical properties
Table 4 shows that there were no significant correlations between the soil’s physicochemical properties and the bacterial diversity within the rhizosphere soil and the bulk soil. To further clarify the interactions between the soil’s physicochemical properties and the microorganisms, we analyzed the correlation between the top 35 species and the soil’s physicochemical properties at the genus level (Figure 8). In the rhizosphere soil samples, Arthrobacter, Ellin6055, Acinetobacter, Solirubrobacter, Rubrobacter, and Aeromonas possessed an extremely significant positive correlation with the AP. Microvirga was positively correlated with the TN. Reyranella was positively correlated with EC. Mesorhizobium, Dongia, Subgroup 10 and Chryseolinea possessed an extremely significant negative correlation with AP. Chryseolinea and Reyranella possessed an extremely significant negative correlation with pH. In the bulk soil samples, Acinetobacter and Rufibacter were significantly positively correlated and negatively correlated with AP, respectively. There was a very significant positive correlation between Saccharopolyspora and NN. In general, AP interacted with the bacterial diversity in both the rhizosphere soil and the bulk soil. TN had significant effects on some bacterial diversity within the rhizosphere soil but had no significant effects on the bacterial diversity within the bulk soil samples. AK, SOM, NN, pH, and EC were all related to bacterial diversity within the different types of soil.
Table 4. Correlations Analysis between bacterial alpha-diversity and soil physicochemical properties.
Figure 8. Heatmap of the correlation between the soil’s physicochemical properties and the genus-level bacteria. (A) is the correlation between physicochemical properties of rhizosphere soils and the genus-level bacteria. (B) is the correlation between physicochemical properties of bulk soils and the genus-level bacteria. * and ** indicate P < 0.05 and P < 0.01, respectively.
Rhizocompartmental bacterial community co-occurrence network analysis
We carried out Spearman correlation analysis and statistical testing on the standardized 16S rRNA gene sequencing data, calculated the correlation between various species, selected the top 100 bacteria with the highest horizontal correlation, and drew the co-expression analysis network diagram based on gephi (Figure 9). Table 5 calculates the important topological properties of different samples to describe the complex patterns of the relationships between nodes. In the root nodules, the APL was 4,465, and the diameter was 14. The CC and modularity index (MD) were 0.410 and 0.576, respectively. In the rhizosphere soil, the APL was 3,767, and the diameter was 12. The CC was 0.599, and the MD was 0.380. In the bulk soil, the APL was 3,102, and the diameter was 11. The CC was 0.560, and the MD was 0.401. A MD value > 0.4 indicates that the network has a modular structure. It can be seen that the root nodule network diagram obviously has a modular structure. In addition, there were 280 correlations founds in the root nodule samples, 370 correlations in the rhizosphere soil samples and 555 correlations in the bulk soil samples. The interaction network of microorganisms in the root nodules is simpler than that in the soil, but the correlation between bacteria is enhanced. The network diagram showed that the bacterial community in root nodules was mainly composed of Proteobacteria, Actinobacteriota, Bacteroidota, Firmicutes, Acidobacteriota, and Myxococcota, in which Proteobacteria were dominant. At the genus level, Ellin6067, Subgroup 10, and Blastococcus showed a robust correlation with other genera (Figure 9A). Compared with root nodules, there were more dominant species of bacteria in the rhizosphere soil and the bulk soil. In the rhizosphere soil, Adhaeribacter, Microlunatus and Ellin6055 showed a wide range of robust correlations (Figure 9B). In the bulk soil, the composition of the bacterial community was the most complex at the phylum level, and there was a robust correlation among the Subgroup 10, Edwardsiella, Bradyrhizobium, and Pedomicrobium (Figure 9C).
Figure 9. Microbial interaction networks in the different compartments. The larger the circle layer of the species is, the greater the connection between the circle layer is, and the higher the correlation degree with other species is, which represents the core species in the correlation. Different colors represent the taxonomic distribution of the corresponding species at the gate level. The edge color indicates a positive (red) and negative (blue) correlation. The width of the edge represents the relevant value. Only a significant correlations (r > 0.6 percent, P < 0.05) were shown. (A) is the co-occurrence network of the bacterial communities within the root nodule. (B) is the co-occurrence network of the bacterial communities within the rhizosphere soil. (C) is the co-occurrence network of bacterial communities within the bulk soil.
Table 5. Topological properties of co-occurring networks.
Discussion
In recent decades, many studies have been focused on the bacteria residing in legume root nodules. In addition to rhizobia, other NREs have also been simultaneously isolated from the root nodules, including Rahnella, Bacillus, Enterobacter, Pantoea, Agrobacterium, Pseudomonas, Cohnella, Herbaspirillum, Paenibacillus, Klebsiella, Endobacter, Chryseobacterium, Stenotrophomonas, Azospirillum, Gluconacetobacter, Variovorax, Arthrobacter, Micromonospora, Hyphomicrobium, Microbacterium, Rhodococcus, Frondihabitans, Kocuria, Providencia, and Sphingomonas, from root nodules of several legume species (De Meyer et al., 2015; Martínez-Hidalgo and Hirsch, 2017; Rahal and Chekireb, 2021). These findings demonstrated that the existence of NREs in leguminous root nodules is a universal phenomenon, and the occurrence of certain rhizobia was correlated with the presence of a particular NRE, suggesting that their presence may not be accidental. Nevertheless, most studies are observational, and few focus on the co-habitation of diverse NREs inside nodules. Moreover, most of the previous studies have mainly involved analysis of the bacterial diversity of rhizobia and legume plant–microbe interactions. Little systemic observation has been paid to the non-symbiotic endophytic bacteria within leguminous root nodules, and their colonization meaning and ecological role are thus far unknown.
A recent study showed that other bacterial genera are present in Trifolium nodules commonly called NREs, and the data presented in this study suggest the presence of a large diversity of bacteria belonging to 50 genera within the nodules of indigenous legumes in Flanders (De Meyer et al., 2015). However, although a few novel non-rhizobia endophytes belonging to various genera have been isolated, the majority of the non-rhizobia bacterial community in the root nodules of several legume species remains uncharacterized. In a previous study, traditional techniques, such as microbial isolation and cultivation under laboratory conditions, were used to investigate the cultural endophytes within root nodules. However, only a small fraction of microorganisms (less than 1%) from the natural environment and root nodules are cultivable. To better understand the diversity of the bacterial microbiota in root nodules, more efficient technology is needed. In recent years, this problem has been resolved by next-generation sequencing of the 16S ribosomal RNA (rRNA) gene as a culture-independent molecular technique. Here, we have explored and characterized the diversity of NRE bacteria associated with S. davidii, a native N-fixing wild leguminous shrub that is a potential tool species for vegetation and ecosystem rehabilitation in hilly and gully regions of the Loess Plateau of China, by using 16S rRNA gene high-throughput amplicon sequencing, which allowed the full-depth analysis of bacterial community diversity within root nodules. In this study, we showed the presence of a large diversity of bacteria belonging to 81 phyla, 154 classes, 333 orders, 463 families, and 732 genera inside nodules associated with S. davidii of indigenous legumes. To our knowledge, this is the first report of such an extensive bacterial diversity analysis in root nodules of indigenous legumes in hilly and gully regions of the Loess Plateau of China, although similar study on cultured legumes in the same region have been reported (Xiao et al., 2017; Zai et al., 2021).
Legumes are known as pioneer plants and as enhancers of nutritional status in cultivated soils, which has been explained by their capacity to recruit rhizobia, since they give an indication of host specificity and nitrogen-fixing ability (Diouf et al., 2010). Legume root nodules constitute an environmental niche for the accommodation of specific traditional rhizobia and NRE bacteria. Previous studies found that Mesorhizobium sp. X, M. waimense, and M. amorphae were the dominant and universal microsymbionts of S. davidii in the forest area of the Loess Plateau in northern Shaanxi Province (Cao et al., 2021). However, to date, there has been no evaluation of the role of NREs in pioneer legumes of S. davidii grown in poor soils on the Loess Plateau of China. In the present study, we also found the dominant genus Mesorhizobium (75.07%) in the root nodules of S. davidii, which is consistent with the results reported by Cao et al. (2021), and confirmed that Mesorhizobium species were dominant symbionts for genus Sophora (Zhao et al., 2010; Jiao et al., 2015; Tan et al., 2015; Tagelsir and Rabab, 2017). Similarly, the legumes may select the specific and appropriate NRE from environment (De Meyer et al., 2015). Taken together with our findings, the leguminous plants not only select a certain and specific NRE but also the rhizobia inside the root nodules prefer certain NREs. Nevertheless, based on the above research background, the co-habitation of diverse bacteria inside nodules raises questions regarding the entry mode of endophytic bacteria and the role of the environment, since the biogeographic patterns and assembly process of rhizobial and NRE communities vary among root nodules, rhizosphere and bulk soil compartments. In an effort to address the issue mentioned above, we compared the biogeography and co-occurrence patterns of the bacterial communities obtained from root nodules samples, rhizosphere samples and bulk soil samples of S. davidii in hilly and gully regions of the Loess Plateau of China.
Microbial composition difference in the root nodule, the rhizosphere, and the bulk soil
The physicochemical properties of the rhizosphere, which is surrounded and affected by plant roots, are different from those of the surrounding bulk soil (Dessaux et al., 2016). Previous studies have reported that the microbiomes of the rhizosphere are a part of the bulk soil community derived from root selection (Philippot et al., 2013; Mendes et al., 2014). However, the microbial community composition of the bulk soil and the rhizosphere environment compartments are distinctly different (Reinhold-Hurek et al., 2015), which revealed that its diversity decreases with proximity to the root (Bulgarelli et al., 2013; Donn et al., 2015). Based on alpha and beta diversity analyses, the bacterial richness and diversity from different habitats were compared, and we found that the microbial community composition of the rhizosphere soil differed from that of the bulk soil of S. davidii. Moreover, the alpha diversities of the microbial communities decreased from the bulk soil to the rhizosphere to the root nodule, which is consistent with the study of Xiao et al. (2017). Our results show that there were 63 phyla of endophytic bacteria within the root nodules, 69 phyla of bacteria within the rhizosphere soil, and 62 phyla of bacteria within the bulk soil. Similar to a previous study (Mendes et al., 2013), we observed in our survey that Proteobacteria was the most dominant phylum, and its proportion in the root nodules (85.77%) was much higher than that in the rhizosphere soil (28.93%). In the bulk soil, only Actinomycetes (32.61%) were enriched, suggesting that bacterial community distribution is habitat-specific (Shi et al., 2014).
Furthermore, an interesting gradient was observed in which the bacterial diversity decreased as the distance from the soil increased. In addition, the proportion of Acidobacteria in rhizosphere soil (12.71%) was significantly higher than that in the nodules (0.74%) and the bulk soil (3.76%), and this finding was consistent with previous reports in many plants and soils (Bulgarelli et al., 2012; Chen et al., 2019; Lin et al., 2019). In the present study, Proteobacteria were absolutely dominant within the root nodules, which may be the result of the differentiation of specific bacteria within different niches and the active filtration of bacteria by the host plants (Bulgarelli et al., 2013; Yuan et al., 2021). This might be because Proteobacteria thrive better under high nutrient conditions (Zhao et al., 2017) and tend to exploit more labile C sources (Fierer et al., 2007). Overall, these findings suggest that the alpha diversity with root proximity showed a decreasing trend owing to the root filtration and selection effect (Dennis et al., 2010). Additionally, previous studies have confirmed that Acidobacteria dynamics are affected by soil pH (Fan et al., 2017). The results of the soil’s physicochemical properties showed that all of the rhizosphere soil samples were alkaline, and the content of Acidobacteria was abundant, which indicated that the living environment of Acidobacteria was not limited to an acidic environment (Xiong et al., 2012).
Important abiotic factors affecting microbial community structure
In this study, a total of 732 different genera were detected, and there were significant differences in the relative abundance of these bacteria. Comparison of the bacterial communities associated with S. davidii revealed both ubiquitous and specific members in different sampling sites, including Sphingomonas, Mesorhizobium, RB41 and Arthrobacter. Moreover, the bacterial richness and diversity varied in the three different compartments, with the highest bacterial diversity occurring in the rhizosphere soil. Similar results have been reported that plant endophytic bacteria mainly come from rhizosphere soil, and the microbial community within the rhizosphere soil can affect the bacterial diversity in or related to the root nodules (Leite et al., 2017). The bacterial community in the rhizosphere soil showed higher species richness than endophytic bacteria (Tan et al., 2017; Yang et al., 2017).
The pivotal influential factors and their contributions to variations in community composition were determined by PCoA (Figure 6). Our findings showed that microbial communities in the three compartments were associated with various edaphic physical and chemical properties, indicating that they had significantly different biogeographic patterns among the three compartments and that they had their own unique bacterial community structure. It is understood that most root endophytic bacteria are from the soil, and the rhizospheric microbiome is derived from the soil (Lennon and Jones, 2011; Philippot et al., 2013). Furthermore, Spearman’s correlation was used to test the effect of the soil’s physical and chemical properties on the relative abundances of the top 35 bacterial genera. The results showed that AP was significantly positively correlated with 6 bacteria (Arthrobacter, Ellin6055, Acinetobacter, Solirubrobacter, Rubrobacter, and Aeromonas) and negatively correlated with 4 bacteria (Mesorhizobium, Dongia, Subgroup 10, and Chryseolinea) in the rhizosphere soil samples (Figure 8), suggesting that bacterial communities and composition within the three different compartments were significantly related to the different soil properties. It is well known that soil P availability affects and influences the microbial community in roots and soil (Ren et al., 2021). Some reports have shown that soil AP is significantly relevant to the bacterial community of S. davidii, indicating that AP may be an important factor affecting the diversity and composition of the bacterial community in the rhizosphere (Pande et al., 2017; Zhou et al., 2020). It has been found that soil microorganisms can transform the combined phosphorus in soil into water-soluble phosphorus, thus promoting the absorption of phosphorus by plants (He and Zhu, 1998) and shifting the rhizosphere bacterial microbiome (Lennon and Jones, 2011). Soil nitrogen plays an important role in the biological nitrogen fixation of legumes (Taylor and Menge, 2018). In this study, a significant Spearman relationship was observed in which the TN, EC, pH and NN were also important factors correlating with the bacterial community variation (Figure 8).
All of these findings support the conclusion that the highest bacterial diversity observed in the rhizosphere soil may be attributed to the primary site of interaction between plants and soil, and plants play an important role in determining the community composition and structure of the rhizospheric microbiome by recruiting specific microbiomes (Berendsen et al., 2012).
The co-occurrence patterns of microbes in the root nodule, the rhizosphere soil, and the bulk soil
Previous studies found that microbial networks were more complex within bulk soil than in the rhizosphere (Mendes et al., 2014; Fan et al., 2017). In the present study, network analysis of S. davidii-associated microbial communities within the root nodule, the rhizosphere and the bulk soil in hilly and gully regions of the Loess Plateau of China was conducted. Our results revealed that the microbial network diagram of the co-occurrence pattern showed that more complexes were observed in soil than in the root nodules, namely, the links and nodes decreased from the bulk soil to the rhizosphere soil to the root nodules (Figure 9), which demonstrated that the network complexity of the bacterial communities in the three different compartments is consistent with the above result of diversity. The rhizosphere community is generally considered a subset of the bulk soil community. This further proves that the host plants will actively filter bacteria (Bulgarelli et al., 2013), leaving useful bacterial interactions to promote the growth and development of host plants. These findings indicated that only specific microorganisms could colonize and thrive in the rhizosphere through the selection and filtering effects of roots.
Moreover, the microbial network in the root nodules showed more positive correlations. Thus, the NREs in the root nodules may have a higher nutritional dependence on rhizobia, resulting in a higher proportion of NRE-rhizobial edges in the root nodule networks, further indicating that the coexistence of the rhizobia and certain NRE bacteria in the root nodule was cooperative and had beneficial effects on the host plants. Taken together, these results showed that the networks in root nodules had more rhizobia and certain NRE bacterial links than the networks in both the rhizosphere and the bulk soil.
Conclusion
Our results demonstrated that the microbial community composition was distinctly different in root nodules and the rhizosphere and bulk soil environment compartments of S. davidii in hilly and gully regions of the Loess Plateau of China, indicating that the environmental filtering differed in the assembly of bacterial communities. Soil AP influenced the bacterial communities in the bulk soil, suggesting that there were significantly different biogeographic patterns among the three compartments and that they had their own unique bacterial community structure. The networks in the root nodules had more rhizobia and certain NRE bacterial links than the networks in both the rhizosphere and the bulk soil samples. These findings indicated that only specific microorganisms could colonize and thrive in the rhizosphere through the selection and filtering effects of roots. Overall, our study improves the comprehension of the biogeographic patterns and diversity of bacterial microbiota inhabiting root nodules and can help elucidate the root nodule assemblage process of S. davidii.
Data availability statement
The data presented in the study are deposited in the GenBank Data repository, accession number PRJNA851563.
Author contributions
LJ, JY-Y, ZC-C, ZR-H, and AJ-M performed the material preparation, experimental operation, data collection, and analysis. LJ and DZ-S wrote the manuscript. LJ and LX-D performed the project guidance and critical revision of manuscripts. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Grant No. 32160003) and the Graduate Education Innovation Program of Yan’an University (Grant No. YCX2021082).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.945127/full#supplementary-material
Footnotes
- ^ http://ccb.jhu.edu/software/FLASH/
- ^ http://qiime.org/scripts/split_libraries_fastq.html
- ^ https://www.arb-silva.de/
- ^ http://www.drive5.com/usearch/manual/uchime_algo.html
- ^ http://drive5.com/uparse/
- ^ http://www.drive5.com/muscle/
References
Ampomah, O. Y., and Huss-Danell, K. (2011). Genetic diversity of root nodule bacteria nodulating Lotus corniculatus and Anthyllis vulneraria in Sweden. Syst. Appl. Microbiol. 34, 267–275. doi: 10.1016/j.syapm.2011.01.006
Ardley, J. K., Parker, M. A., De Meyer, S. E., Trengove, R. D., O’Hara, G. W., Reeve, W. G., et al. (2012). Microvirga lupini sp. nov., Microvirga lotononidis sp. nov. and Microvirga zambiensis sp. nov. are alphaproteobacterial root-nodule bacteria that specifically nodulate and fix nitrogen with geographically and taxonomically separate legume hosts. Int. J. Syst. Evol. Microbiol. 62, 2579–2588. doi: 10.1099/ijs.0.035097-0
Aserse, A. A., Rasanen, L. A., Aseffa, F., Hailemariam, A., and Lindstrom, K. (2013). Diversity of sporadic symbionts and nonsymbiotic endophytic bacteria isolated from nodules of woody, shrub, and food legumes in Ethiopia. Appl. Microbiol. Biotechnol. 97, 10117–10134. doi: 10.1007/s00253-013-5248-4
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.
Barnawal, D., Bharti, N., Maji, D., Chanotiya, C. S., and Kalra, A. (2014). ACC deaminase-containing Arthrobacter protophormiae induces NaCl stress tolerance through reduced ACC oxidase activity and ethylene production resulting in improved nodulation and mycorrhization in Pisum sativum. J. Plant Physiol. 171, 884–894. doi: 10.1016/j.jplph.2014.03.007
Bastian, M., Heymann, S., and Jacomy, M. (2009). Gephi: an open source software for exploring and manipulating networks. Int. Conf. Weblogs Soc. Media 8, 361–362.
Beiierinck, M. W., and van Delden, A. (1902). Ueber die Assimilation des freien Stickstoffs durch Bakterien. Zentbl. Bakt. Parasitenk. Infekt 9, 3–43.
Benhizia, Y., Benhizia, H., Benguedoua, A., Muresu, R., Giacomini, A., and Squartini, A. (2004). Gamma proteobacteria can nodulate legumes of the genus Hedysarum. Syst. Appl. Microbiol. 27, 462–468. doi: 10.1078/0723202041438527
Berendsen, R. L., Pieterse, C. M. J., and Bakker, P. A. M. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486.
Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10, 57–59. doi: 10.1038/nmeth.2276
Bulgarelli, D., Rott, M., Schlaeppi, K., Loren, E. V., van Themaat, N., Ahmadinejad, F. A., et al. (2012). Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 488, 91–95. doi: 10.1038/nature11336
Bulgarelli, D., Schlaeppi, K., Spaepen, S., Ver Loren van Themaat, E., and Schulze-Lefert, P. (2013). Structure and functions of the bacterial microbiota of plants. Annu. Rev. Plant Biol. 64, 807–838.
Cao, Y., Tie, D., Zhao, J. L., Wang, X. B., Yi, J. J., Chai, Y. F., et al. (2021). Diversity and distribution of Sophora davidii rhizobia in habitats with different irradiances and soil traits in Loess Plateau area of China. Syst. Appl. Microbiol. 44:126224. doi: 10.1016/j.syapm.2021.126224
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.
Chen, L., Fang, K., Zhou, J., Yang, Z. P., Dong, X. F., Dai, G. H., et al. (2019). Enrichment of soil rare bacteria in root by an invasive plant Ageratina adenophora. Sci. Total Environ. 683, 202–209. doi: 10.1016/j.scitotenv.2019.05.220
Cummings, S. P., Gyaneshwar, P., Vinuesa, P., Farruggia, F. T., Andrews, M., Humphry, D., et al. (2009). Nodulation of Sesbania species by Rhizobium (Agrobacterium) strain IRBG74 and other rhizobia. Environ. Microbiol. 11, 2510–2525. doi: 10.1111/j.1462-2920.2009.01975.x
De Meyer, S. E., De Beuf, K., Vekeman, B., and Willems, A. (2015). A large diversity of non-rhizobial endophytes found in legume root nodules in Flanders (Belgium). Soil Biol. Biochem. 83, 1–11.
Deng, Z. S., Zhao, L. F., Kong, Z. Y., Yang, W. Q., Lindström, K., Wang, E. T., et al. (2011). Diversity of endophytic bacteria within nodules of the Sphaerophysa salsula in different regions of Loess Plateau in China. FEMS Microbiol. Ecol. 76, 463–475. doi: 10.1111/j.1574-6941.2011.01063.x
Dennis, P. G., Miller, A. J., and Hirsch, P. R. (2010). Are root exudates more important than other sources of rhizodeposits in structuring rhizosphere bacterial communities. FEMS Microbiol. Ecol. 72, 313–327.
Dessaux, Y., Grandclement, C., and Faure, D. (2016). Engineering the Rhizosphere. Trends Plant Sci. 21, 266–278.
Diouf, D., Fall, D., Chantreuil, C., Ba, A. T., Dreyfus, B., Neyra, M., et al. (2010). Phylogenetic analyses of symbiotic genes and characterization of functional traits of Mesorhizobium spp. strains associated with the promiscuous species Acacia seyal Del. J. Appl. Microbiol. 108, 818–830. doi: 10.1111/j.1365-2672.2009.04500.x
Diouf, D., Samba-Mbaye, R., Lesueur, D., Ba, A. T., Dreyfus, B., de Lajudie, P., et al. (2007). Genetic diversity of Acacia seyal Del. rhizobial populations indigenous to Senegalese soils in relation to salinity and pH of the sampling sites. Microb. Ecol. 54, 553–566. doi: 10.1007/s00248-007-9243-0
Donn, S., Kirkegaard, J. A., Perera, G., Richardson, A. E., and Watt, M. (2015). Evolution of bacterial communities in the wheat crop rhizosphere. Environ. Microbiol. 17, 610–621.
Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797.
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461.
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Fan, K., Cardona, C., Li, Y., Shi, Y., Xiang, X., Shen, C., et al. (2017). Rhizosphere-associated bacterial network structure and spatial distribution differ significantly from bulk soil in wheat crop fields. Soil Biol. Biochem. 113, 275–284.
Fierer, N., Braddord, M. A., and Jackson, R. B. (2007). Toward an ecologyical classificaition of soil bacteria. Ecology 88, 1354–1364.
Gopalakrishnan, S., Vadlamudi, S., Samineni, S., and Sameer Kumar, C. V. (2016). Plant growth promotion and biofortification of chickpea and pigeonpea through inoculation of biocontrol potential bacteria, isolated from organic soils. Springerplus 5:1882. doi: 10.1186/s40064-016-3590-6
Haas, B. J., Gevers, D., Earl, A. M., Feldgarden, M., Ward, D. V., Giannoukos, G., et al. (2011). Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. doi: 10.1101/gr.112730.110
He, Z. L., and Zhu, J. (1998). Microbial utilization and transformation of phosphate adsorbed by variable charge minerals. Soil Biol. Biochem. 30, 917–923.
Hoque, M. S., Broadhurst, L. M., and Thrall, P. H. (2011). Genetic characterization of root-nodule bacteria associated with Acacia salicina and A. stenophylla (Mimosaceae) across south-eastern Australia. Int. J. Syst. Evol. Microbiol. 61, 299–309. doi: 10.1099/ijs.0.021014-0
Ibáñez, F., Angelini, J., Taurian, T., Laura Tonelli, M., and Fabra, A. (2008). Endophytic occupation of peanut root nodules by opportunistic Gammaproteobacteria. Syst. Appl. Microbiol. 32, 49–55. doi: 10.1016/j.syapm.2008.10.001
Jiao, Y. S., Liu, Y. H., Yan, H., Wang, E. T., Tian, C. F., Chen, W. X., et al. (2015). Rhizobial diversity and nodulation characteristics of the extremely promiscuous legume Sophora flavescens. Mol. Plant-Microbe Interact. 28, 1338–1352. doi: 10.1094/MPMI-06-15-0141-R
Kara, E. L., Hanson, P. C., Hu, Y. H., Winslow, L., and McMahon, K. D. (2013). A decade of seasonal dynamics and co-occurrences within freshwater bacterioplankton communities from eutrophic Lake Mendota. WI, USA. ISME J. 7, 680–684. doi: 10.1038/ismej.2012.118
Lai, W. A., Hameed, A., Lin, S. Y., Hung, M. H., Hsu, Y. H., Liu, Y. C., et al. (2015). Paenibacillus medicaginis sp. nov. a chitinolytic endophyte isolated from the root nodule of alfalfa (Medicago sativa L.). Int. J. Syst. Evol. Microbiol. 65, 3853–3860. doi: 10.1099/ijsem.0.000505
Leite, J., Fischer, D., Rouws, L. F. M., Fernandes-Júnior, P. I., Hofmann, A., Kublik, S., et al. (2017). Cowpea nodules harbor non-rhizobial bacterial communities that are shaped by soil type rather than plant genotype. Front. Plant Sci. 7:2064. doi: 10.3389/fpls.2016.02064
Lennon, J. T., and Jones, S. E. (2011). Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat. Rev. Microbiol. 9, 119–130. doi: 10.1038/nrmicro2504
Li, J. H., Wang, E. T., Chen, W. F., and Chen, W. X. (2007). Genetic diversity and potential for promotion of plant growth detected in nodule endophytic bacteria of soybean grown in Heilongjiang province of China. Soil Biol. Biochem. 40, 238–246.
Li, Y. M., Yang, Y., Wu, T. E., Zhang, H., Wei, G. H., and Li, Z. F. (2021). Rhizosphere bacterial and fungal spatial distribution and network pattern of Astragalus mongholicus in representative planting sites differ the bulk soi. Appl. Soil Ecol. 168:104114.
Lin, Y. B., Ye, Y. M., Hu, Y. M., and Shi, H. K. (2019). The variation in microbial community structure under different heavy metal contamination levels in paddy soils. Ecotoxicol. Environ. Saf. 180, 557–564. doi: 10.1016/j.ecoenv.2019.05.057
Liu, Y., Li, D., Qi, J., Peng, Z., Chen, W., Wei, G., et al. (2021). Stochastic processes shape the biogeographic variations in core bacterial communities between aerial and belowground compartments of common bean. Environ. Microbiol. 23, 949–964. doi: 10.1111/1462-2920.15227
Ma, Y. R., Zhou, T. X., Zhao, P., Choi, H. Y., Hao, J., Huang, H. Q., et al. (2021). New flavonoids from the roots of Sophora davidii (Franch.) Skeels and their glucose transporter 4 translocation activities. Bioorg. Chem. 106:104500. doi: 10.1016/j.bioorg.2020.104500
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Martínez-Hidalgo, P., and Hirsch, A. M. (2017). The nodule microbiome: N2-fixing rhizobia do not live alone. Phytobiomes 1, 70–82.
Mendes, L. W., Kuramae, E. E., Navarrete, A. A., van Veen, J. A., and Tsai, S. M. (2014). Taxonomical and functional microbial community selection in soybean rhizosphere. ISME J. 8, 1577–1587.
Mendes, R., Garbeva, P., and Raaijmakers, J. M. (2013). The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663.
Olenska, E., Imperato, V., Małek, W., Włostowski, T., Wójcik, M., Swiecicka, I., et al. (2020). Trifolium repens-associated bacteria as a potential tool to facilitate phytostabilization of zinc and lead polluted waste heaps. Plants 9:1002. doi: 10.3390/plants9081002
Palaniappan, P., Chauhan, P. S., Saravanan, V. S., Anandham, R., and Sa, T. M. (2010). Isolation and characterization of plant growth promoting endophytic bacterial isolates from root nodule of Lespedeza sp. Biol. Fertil. Soils 46, 807–816.
Pande, A., Pandey, P., Mehra, S., Singh, M., and Kaushik, S. (2017). Phenotypic and genotypic characterization of phosphate solubilizing bacteria and their efficiency on the growth of maize. J. Genet. Eng. Biotechnol. 15, 379–391.
Pandya, M., Kumar, G. N., and Rajkumar, S. (2013). Invasion of rhizobial infection thread by non-rhizobia for colonization of Vigna radiata root nodules. FEMS Microbiol. Lett. 348, 58–65. doi: 10.1111/1574-6968.12245
Philippot, L., Raaijmakers, J. M., Lemanceau, P., and van der Putten, W. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, 590–596. doi: 10.1093/nar/gks1219
Rahal, S., and Chekireb, D. (2021). Diversity of rhizobia and non-rhizobia endophytes isolated from root nodules of Trifolium sp. growing in lead and zinc mine site Guelma, Algeria. Arch. Microbiol. 203, 3839–3849. doi: 10.1007/s00203-021-02362-y
Rajendran, G., Sing, F., Desai, A. J., and Archana, G. (2008). Enhanced growth and nodulation of pigeon pea by co-inoculation of Bacillus strains with Rhizobium spp. Bioresour. Technol. 99, 4544–4550. doi: 10.1016/j.biortech.2007.06.057
Rasolomampianina, R., Bailly, X., Fetiarison, R., Rabevohitra, R., Béna, G., Ramaroson, L., et al. (2005). Nitrogen-fixing nodules from rose wood legume trees (Dalbergia spp.) endemic to Madagascar host seven different genera belonging to α-and β-proteobacteria. Mol. Ecol. 14, 4135–4146. doi: 10.1111/j.1365-294X.2005.02730.x
Raupach, G. S., and Kloepper, J. W. (1998). Mixtures of plant growth-promoting rhizobacteria enhance biological control of multiple cucumber pathogens. Phytopathology 88, 1158–1164.
Reinhold-Hurek, B., Bunger, W., Burbano, C. S., Sabale, M., and Hurek, T. (2015). Roots shaping their microbiome: global hotspots for microbial activity. Annu. Rev. Phytopathol. 53, 403–424. doi: 10.1146/annurev-phyto-082712-102342
Ren, C. J., Zhou, Z. H., Guo, Y. X., Yang, G. H., Zhao, F. Z., Wei, G. H., et al. (2021). Contrasting patterns of microbial community and enzyme activity between rhizosphere and bulk soil along an elevation gradient. CATENA 196:104921.
Sarah, R., and Djamel, C. (2021). Diversity of rhizobia and non-rhizobia endophytes isolated from root nodules of Trifolium sp. growing in lead and zinc mine site Guelma, Algeria. Arch. Microbiol. 203, 3839–3849.
Sawana, A., Adeolu, M., and Gupta, R. S. (2014). Molecular signatures and phylogenomic analysis of the genus Burkholderia: proposal for division of this genus into the emended genus Burkholderia containing pathogenic organisms and a new genus Paraburkholderia gen. nov. harboring environmental species. Front. Genet. 5:429. doi: 10.3389/fgene.2014.00429
Schrey, S. D., and Tarkka, M. T. (2008). Friends and foes: streptomycetes as modulators of plant disease and symbiosis. Antonie van Leeuwenhoek 94, 11–19. doi: 10.1007/s10482-008-9241-3
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:60.
Shi, Y. W., Yang, H. M., Zhang, T., Sun, J., and Lou, K. (2014). Illumina-based analysis of endophytic bacterial diversity and space-time dynamics in sugar beet on the north slope of Tianshan mountain. Appl. Microbiol. Biotechnol. 98, 6375–6385. doi: 10.1007/s00253-014-5720-9
Song, C. J., Ma, K. M., Fu, B. J., Qu, L. Y., Xu, X. L., Liu, Y., et al. (2010). Distribution patterns of shrubby N-fixers and non-N fixers in an arid valley in Southwest China: implications for ecological restoration. Ecol. Res. 25, 553–564.
Sreevidya, M., Gopalakrishnan, S., Kudapa, H., and Varshney, R. K. (2016). Exploring plant growth-promotion actinomycetes from vermicompost and rhizosphere soil for yield enhancement in chickpea. Braz. J. Microbiol. 47, 85–95. doi: 10.1016/j.bjm.2015.11.030
Sy, A., Giraud, E., Jourand, P., Garcia, N., Willems, A., de Lajudie, P., et al. (2001). Methylotrophic Methylobacterium bacteria nodulate and fix nitrogen in symbiosis with legumes. J. Bacteriol. 183, 214–220. doi: 10.1128/JB.183.1.214-220.2001
Tagelsir, H. M. A., and Rabab, M. I. M. (2017). Diversity, computability and nitrogen fixation of wild legmes rhizobia. Int. J. Adv. Res. 5, 594–602.
Tai, Z., Cai, L., Dai, L., Dong, L. H., Wang, M. F., Yang, Y. B., et al. (2011). Antioxidant activity and chemical constituents of edible flower of Sophora viciifolia. Food Chem. 126, 1648–1654. doi: 10.1016/j.foodchem.2010.12.048
Tan, H. W., Heenan, P. B., De Meyer, S. E., Willems, A., and Andrews, M. (2015). Diverse novel mesorhizobia nodulate New Zealand native Sophora species. Syst. Appl. Microbiol. 38, 91–98. doi: 10.1016/j.syapm.2014.11.003
Tan, Y., Cui, Y. S., Li, H. Y., Kuang, A. X., Li, X. R., Wei, Y. L., et al. (2017). Diversity and composition of rhizospheric soil and root endogenous bacteria in Panax notoginseng during continuous cropping practices. J. Basic Microbiol. 57, 337–344.
Taylor, B. N., and Menge, D. N. L. (2018). Light regulates tropical symbiotic nitrogen fixation more strongly than soil nitrogen. Nat. Plants 4, 655–661. doi: 10.1038/s41477-018-0231-9
Tian, X. Y., and Zhang, C. S. (2017). Illumina-based analysis of endophytic and rhizosphere bacterial diversity of the Coastal Halophyte Messerschmidia sibirica. Front. Microbiol. 8:2288. doi: 10.3389/fmicb.2017.02288
United States Department of Agriculture [USDA] (1996). Soil Survey Laboratory Methods Manual; Version 3.0; Soil Survey Investigations. Report No. 42. Lincoln, NE: United States Department of Agriculture, Natural Resources Conservation Service, National Soil Survey Centre.
Valverde, A., Velazquez, E., Gutierrez, C., Cervantes, E., Ventosa, A., and Igual, J. M. (2003). Herbaspirillum lusitanum sp. nov., a novel nitrogen-fixing bacterium associated with root nodules of Phaseolus vulgaris. Int. J. Syst. Evol. Microbiol. 53, 1979–1983. doi: 10.1099/ijs.0.02677-0
Vincent, J. M. (1970). A manual for Practical study of Root Nodule Bacteria. IBP Handbook No. 15. Oxford: Blackwell Scientific Publishers, 164.
Wang, X. P., and Zabowski, D. (1998). Nutrient composition of Douglas-fir rhizosphere and bulk soil solutions. Plant Soil 200, 13–20.
Weiss, V. A., Faoro, H., Tadra-Sfeir, M. Z., Raittz, R. T., de Souza, E. M., and Monteiro, R. A. (2012). Draft genome sequence of Herbaspirillum lusitanum P6-12, an endophyte isolated from root nodules of Phaseolus vulgaris. J. Bacteriol. 194, 4136–4137.
Wu, F., Bao, W., Li, F., and Wu, N. (2007). Effects of drought stress and N supply on the growth, biomass partitioning and water-use efficiency of Sophora davidii seedlings. Environ. Exp. Bot. 63, 248–255.
Xiao, X., Chen, W., Zong, L., Yang, J., Jiao, S., Lin, Y., et al. (2017). Two cultivated legume plants reveal the enrichment process of the microbiome in the rhizocompartments. Mol. Ecol. 26, 1641–1651. doi: 10.1111/mec.14027
Xiong, J. B., Liu, Y. Q., Lin, X. G., Zhang, H. Y., Zeng, J., Hou, J. Z., et al. (2012). Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan Plateau. Environ. Microbiol. 14, 2457–2466. doi: 10.1111/j.1462-2920.2012.02799.x
Xu, L., Zhang, Y., Wang, L., Chen, W., and Wei, G. (2014). Diversity of endophytic bacteria associated with nodules of two indigenous legumes at different altitudes of the Qilian Mountains in China. Syst. Appl. Microbiol. 37, 457–465. doi: 10.1016/j.syapm.2014.05.009
Yang, Y., Wang, N., Guo, X., Zhang, Y., and Ye, B. (2017). Comparative analysis of bacterial community structure in the rhizosphere of maize by high-throughput pyrosequencing. PLoS One 12:e0178425. doi: 10.1371/journal.pone.0178425
Yuan, Z. S., Liu, F., Liu, Z. Y., Huang, Q. L., Zhang, G. F., and Pan, H. (2021). Structural variability and differentiation of niches in the rhizosphere and endosphere bacterial microbiome of moso bamboo (Phyllostachys edulis). Sci. Rep. 11, 1574–1574. doi: 10.1038/s41598-021-80971-9
Zai, X. Y., Luo, W., Bai, W. Q., Li, Y. H., Xiao, X., Gao, X., et al. (2021). Effect of root diameter on the selection and network interactions of root-associated bacterial microbiomes in Robinia pseudoacacia L. Microb. Ecol. 82, 391–402. doi: 10.1007/s00248-020-01678-4
Zakhia, F., Jeder, H., Willems, A., Gillis, M., Dreyfus, B., and de Lajudie, P. (2006). Diverse bacteria associated with root nodules of spontaneous legumes in Tunisia and frst report for nifH-like gene within the genera Microbacterium and Starkeya. Microb. Ecol. 51, 375–393. doi: 10.1007/s00248-006-9025-0
Zgadzaj, R., Garrido-Oter, R., Jensen, D. B., Koprivova, A., SchulzeLefert, P., and Radutoiu, S. (2016). Root nodule symbiosis in Lotus japonicus drives the establishment of distinctive rhizosphere, root, and nodule bacterial communities. Proc. Natl. Acad. Sci. U.S.A. 113, 7996–8005. doi: 10.1073/pnas.1616564113
Zgadzaj, R., James, E. K., Kelly, S., Kawaharada, Y., de Jonge, N., Jensen, D. B., et al. (2015). A legume genetic framework controls infection of nodules by symbiotic and endophytic bacteria. PLoS Genet. 11:e1005280. doi: 10.1371/journal.pgen.1005280
Zhang, C. S., Lin, Y., Tian, X. Y., Xu, Q., Chen, Z. H., and Lin, W. (2017). Tobacco bacterial wilt suppression with biochar soil addition associates to improved soil physiochemical properties and increased rhizosphere bacteria abundance. Appl. Soil Ecol. 112, 90–96.
Zhang, K., Adams, J. M., Shi, Y., Yang, T., Sun, R., He, D., et al. (2017). Environment and geographic distance differ in relative importance for determining fungal community of rhizosphere and bulk soil. Environ. Microbiol. 19, 3649–3659.
Zhang, Z. J., Dong, S. W., Gao, D. D., Du, X. Y., Xie, Y. Q., Xia, X. S., et al. (2021). Unusual matrine-adenine hybrids isolated from Sophora davidii and their inhibitory effects on human cytomegalovirus. Phytochemistry 190, 112842–112842. doi: 10.1016/j.phytochem.2021.112842
Zhao, J., Ni, T., Xun, W., Huang, X., Huang, Q., Ran, W., et al. (2017). Influence of straw incorporation with and without straw decomposer on soil bacterial community structure and function in a rice-wheat cropping system. Appl. Microbiol. Biotechnol. 101, 4761–4773. doi: 10.1007/s00253-017-8170-3
Zhao, L., Deng, Z., Yang, W., Cao, Y., Wang, E., and Wei, G. (2010). Diverse rhizobia associated with Sophora alopecuroides grown in different regions of Loess Plateau in China. Syst. Appl. Microbiol. 33, 468–477. doi: 10.1016/j.syapm.2010.08.004
Keywords: Sophora davidii, root nodules, high-throughput sequencing, diversity, biogeographic patterns, network pattern
Citation: Jing L, Jia-min A, Xiao-dong L, Ying-ying J, Chao-chao Z, Rui-hua Z and Zhen-shan D (2022) Environmental filtering drives the establishment of the distinctive rhizosphere, bulk, and root nodule bacterial communities of Sophora davidii in hilly and gully regions of the Loess Plateau of China. Front. Microbiol. 13:945127. doi: 10.3389/fmicb.2022.945127
Received: 16 May 2022; Accepted: 30 June 2022;
Published: 22 July 2022.
Edited by:
Entao Wang, Instituto Politécnico Nacional (IPN), MexicoCopyright © 2022 Jing, Jia-min, Xiao-dong, Ying-ying, Chao-chao, Rui-hua and Zhen-shan. 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.
*Correspondence: Deng Zhen-shan, emhlbnNoYW5kZW5nMjE0QDE2My5jb20=