The OS and WT under chilling stress
In order to gain insight about the ability of rice plants to eliminate ROS, we measured the proline (Pro) contents in both WT and OS at the low temperature and normal temperature. There was no significant difference in Pro content between the OS (18.85 µg·g−1·FW) and WT (19.39 µg·g−1·FW) at normal temperature (Figure 1). The content of Pro in OS and WT was significantly increased at low temperature than normal temperature, and the Pro content in OS (41.32 µg·g−1·FW) was significant more than WT (29.25 µg·g−1·FW), increased by 41.26% at low temperature. In addition, it also showed that the Pro content increased in OS plants under low temperature stress, which protected the cell membrane system from damage and enhanced the cold resistance of the plants.
The fluxes of Ca2+, K+, Na+ and H+ in WT and OS
To estimate the flow of some important ions involved in cold response, maintaining osmotic balance and defense signaling in roots, we used the NMT technology. There was no significant difference in Ca2+, K+, Na+ and H+ between WT and OS plants at normal temperature, and a small amount were excreted. There were significant differences observed in the fluxes of ions between the WT and OS (Figure 2 and 3). The flow rate of Ca2+ revealed that the roots of the OS plants assimilated greater Ca2+ as compared to WT. Besides, the roots of the WT plants absorbed Ca2+ and released Ca2+ occasionally. Furthermore, the OS plants showed absorption of Ca2+ as a whole, and the rate reached 54.31 pmol/cm/s. However, the WT plants showed efflux of Ca2+ as a whole, and the rate reached 0.92 pmol·cm−2·s−1. In addition, the OS roots showed significant absorption of K+ and the rate reached 12.78 pmol·cm−2·s−1. However, the roots of the WT plants excreted a more quantity of K+, and the rate reached 157.02 pmol·cm−2·s−1. The flow rate of Na+ shows that the roots of the OS plants excreted more Na+, and the rate reached 43.02 pmol·cm−2·s−1 which was significantly different from the efflux in WT, which was 14.04 pmol·cm−2·s−1. Similarly, In Case of H+ ions there distinct difference between OS (reached 10.21 pmol·cm−2·s−1) and WT (0.09 pmol·cm−2·s−1). These results depicts that under low temperature stress, OS plants absorbed Ca2+, K+ and H+, while efflux of Na+ was observed, which maintains the stability of cell membrane function, and enhance the cold resistance of the plants (reference).
Measurement of Calcium deposition and H+-ATPase activity
We found a lot of Ca2+ ions flow in both genotypes, for further confirmation, we measured the Ca2+ deposition in the root tips of WT and OS. And found that the OS contained higher deposition of Ca2+ than the WT under chilling stress, it indicated the involvement of Ca2+ in the cold response (Figure 4A). The results of ion absorption in this study also indicated that OS plants could absorb Ca2+ in a large amount under low temperature stress. The root calcium deposition experiment confirmed that the calcium content of OS roots was significantly higher than that of WT. Moreover, the H+-ATPase activity was also measured in the roots of both genotypes. The OS exhibited higher H+-ATPase activity than WT which increases the energy level for the influx and efflux of ions is OS (Figure 4B).
OTU cluster and species annotation
Soil bacterial community related to WT and OS defense against cold was unraveled by utilizing deep 16s rDNA pyrosequencing. Across all soil samples, a total of 1,457,835 effective tags containing species annotation was identified, with an average of 76,606 effective tags per sample. Among total effective tags, 40,400 singletons, accounting for 2.9% of total tags, were removed before further analysis (Figure S1). Rarefaction curves demonstrated the sequencing depth, the number of observed species numbers as shown in Figure S2. At a 97% sequence similarity cut-off, we obtained a sum of 54,655 OTUs across the 18 samples. The OTUs numbers in WT R, WT RS, WT N, OS R, OS RS and OS N plots were 3,132, 3,023, 2,905, 3,535, 3,411 and 2,212, respectively. On average, we were able to classify about 99.4%, 97.9%, 94.3%, 85.3%, 62.8% and 23.2% of effective sequences at the phylum, classes, orders, families, and genus and species level, respectively.
Alpha Diversity Indices
Alpha diversity was determined to further investigate species diversity within a sample. The Chao 1 estimator and the ACE were calculated to estimate the abundance of each sample. Both the Shannon and Simpson indices were determined to estimate the diversity within each individual sample. Besides, the bacterial community analysis showed that a significantly higher bacterial diversity was observed in OS R and OS RS than WT R and WT RS. Whereas, less bacterial diversity was observed in OS N than WT N (P ≤ 0.05). Also, there was no significant differences in the identified bacterial species, Shannon, Simpson, Chao 1 and ACE diversity indices between WT R and WT RS (P > 0.05). The richness indices including the Shannon, Simpson, Chao 1 and ACE indices in OS N were lower than other treatments. Both the species observed, Shannon, Simpson, Chao 1 and ACE diversity indices significantly reduced from OS R to OS N (P ≤ 0.05). However, only Shannon and Simpson diversity indices significantly reduce from WT R to WT N (Table 1).
Table 1
Calculations of observed species, richness and diversity in different soil samples.
Treatments
|
Observed species
|
Shannon
|
Simpson
|
Chao 1
|
ACE
|
WT R
|
2686.33b
|
8.46b
|
0.992a
|
3058.32b
|
3126.70bc
|
WT RS
|
2604.67b
|
8.45b
|
0.991a
|
2966.30b
|
3038.65b
|
WT N
|
2483.00b
|
7.85c
|
0.982b
|
2827.76b
|
2900.83b
|
OS R
|
3116.33a
|
9.01a
|
0.994a
|
3485.23a
|
3536.36a
|
OS RS
|
2976.33a
|
8.67b
|
0.992a
|
3329.85a
|
3400.02ab
|
OS N
|
1873.67c
|
6.33d
|
0.934c
|
2222.85c
|
2290.70c
|
Note: WT and OS represent the wild type and Lsi1-Overexpression Dular. R、RS and N represent rhizosphere, rhizoplane and endosphere, respectively. Whereas, “a”, “b”, “c” and “d” shows significance level at P ≤ 0.05. Different letters in columns show significant differences determined by Tukey’s test (P ≤ 0.05, n = 3). |
Beta Diversity Indices
Beta diversity indices were determined by the differences in bacterial community composition and structure between samples. In this study, the weighted Unifrac distance between WT R and WT RS was a minimum of 0.074, but 0.151 in the OS R and OS RS, showed a higher similarity of the soil bacterial community between WT R and WT RS than OS R and OS RS. The weighted Unifrac distances between OS N vs OS R (0.322) and OS N vs OS RS (0.290) more than WT N vs WT R (0.268) and WT N vs WT RS (0.239), showed a higher difference of the soil bacterial community from the R to N in OS than WT. The weighted Unifrac distances between OS R vs WT R, OS RS vs OS RS and OS N vs WT N were 0.093, 0.108 and 0.231, respectively. Which being reason for the highest difference of the soil bacterial community observed in N than others (R and RS) as shown in Figure S3.
PCoA analysis
PCoA was performed based on the weighted UniFrac distance. In addition, the PCoA analysis identified distinct differences in soil bacterial community structure between different treatments. Nevertheless, the first two components (PC1 and PC2) of PCoA explained 60.17% and 18.34% of the total bacterial community variations, respectively as shown in Figure 5. Furthermore, the PCoA analysis showed that the bacterial communities of WT N and OS N were far from other groups and OS N was most diverse among all groups.
Analysis of bacterial communities
The cold significantly changed the bacterial communities in rhizosphere compartments of OS and WT. We identified 52 phyla, however, majority of OTUs belonged to 10 phyla, accounting for more than 95% of the total, including Proteobacteria, Firmicutes, Bacteroidetes, Cyanobacteria, Acidobacteria, Chloroflexi, Verrucomicrobia, Actinobacteria, Planctomycetes and Gemmatimonadetes. The relative abundance of these phyla is shown in (Figure S4). The content of proteobacteria in WT was significantly lower in N than R and RS, but there were no significant differences identified between R and RS, whereas, the proteobacteria content in OS was higher in N and was significantly decreased from RS to R. While the content of firmicutes and Cyanobacteria in WT was significantly increased from R to N. Nevertheless, in OS the firmicutes decreased significantly from R to N. While the Cyanobacteria increased from R to N as shown in Table S1.
Venn diagram analysis of identified genera in OS as compared to WT
We not only identified 84 bacterial genera from R, RS and N but also by Venn diagram analysis we detected the exclusive and shared genera between six different soil samples (Figure 6). The exclusive genera were 20, 18 and 30 in R, RS and N respectively. There were 8, 1 and 4 genera in common in R and RS, RS and N, R and N respectively, and 3 genera in common across in R, RS and N.
Analysis of dominant genera
In this study, 84 bacterial genera have been identified in six different treatments were shown in heat map Figure 7. We classified the identified 84 genera into 10 functional categories, including denitrification (7.14%), defense (5.95%), azotification (8.33%), photosynthesis (4.76%), chemotrophy (11.9), energy (1.19%), biodegradation (9.52%), carbon metabolism (4.76%), pathogenicity (26.19%) and unknown (20.24%). The distribution of microbial communities in R, RS and N were significantly different (Figure 8, S2).
The microorganism for denitrification and biodegradation were accumulated in WT N. On the contrary, the microorganism related to azotification and energy was accumulated in OS N. While the microorganism for defense and photosynthesis were accumulated in OS R, RS and N. Moreover, the pathogenicity related microorganisms were accumulated in WT RS and N. These results indicate that rice plants have a significant shaping effect on the composition and structure of microbial communities in the root environment. In contrast, the results of microbial community in the root environment also affect the growth and development of plants. OS plants can select and enrich a large number of beneficial microorganisms (defense, azotification, photosynthesis and energy) and block the invasion of harmful microorganisms under low temperature stress. These results indicate that rice plants have a significant shaping effect on the composition and structure of microbial communities in the root environment. Moreover, the results of microbial community in the root environment also affect the growth and development of plants.
The OS plants absorb Ca2+, K+, H+ and excrete Na+ thus making OS better equipped to maintain the stability of membrane function. Ion transport in the root system of rice not only an important way to regulate plant response to external low-temperature stress, but also an essential carrier material to mediate material exchange and energy flow between root system and soil (Figure 9). In brief, it can be considered as one of key factors that contributes to the microecological characteristics in different root environments of plants, and it could be the basis for the establishment of plant-microorganism interaction also.