3.1. PGPR isolation and characterization
After screening PGPRs from cucumber rhizosphere, BC56 was selected as a PGPR strain based on its ability to promote cucumber seedlings growth (Fig. S1). Specifically, BC56 was identified for its growth-promoting ability, NH3 production, siderophore production and phosphate solubility (Table S1). Analysis of the 16S rDNA sequence showed that BC56 has 100% sequence homology with B. cereus (accession NR074540.1) and the sequence of BC56 was submitted to the GenBank database under accession number OQ271230 (Fig. S1).
3.2. Effect of different NaCl concentration on cucumber growth
To determine the appropriate salt tolerance of cucumber seedlings, two-leaf stage cucumber seedlings were grown in 1/4 Hoagland's solution supplemented with different NaCl concentrations. As the NaCl concentration increased, the growth of the cucumber seedlings decreased progressively. At 125 mM NaCl, the growth of cucumber seedlings almost stopped (Fig. S2). Therefore, we use 100 mM NaCl to treat cucumber for future use.
3.3. Effects of B. cereus BC56 on cucumber growth under salt stress
BC56 had no significant effect on shoot length and fresh weight of cucumber seedlings not exposed to salt (Fig. 1a and 1c). The root length of BC56-uninoculated seedlings was approximately 29.9 cm, whereas the root length of BC56-inoculated seedlings was approximately 32.8 cm at 14 DAI (Fig. 1b). These results showed that BC56 could promote the growth of roots rather than shoots of cucumber seedlings.
Under 100 mM NaCl treatment, the shoot length of BC56-inoculated seedlings was 1.37-fold and 1.28-fold longer than that of BC56-uninoculated seedlings at 10 and 14 DAI, respectively (Fig. 1d). The root length of BC56-inoculated seedlings was 1.81-, 1.16- and 1.16-fold longer than that of BC56-uninoculated seedlings at 6, 10, and 14 DAI, respectively (Fig. 1e). BC56-inoculated seedlings were 1.19- and 1.57-fold heavier than BC56-uninoculated seedlings in shoot fresh weight and root fresh weight at 14 DAI, respectively (Fig. 1f). These results indicated that BC56 could counteract the decrease in biomass of cucumber under salt stress.
3.4. Effects of B. cereus BC56 on biochemical index of cucumber under salt stress
To investigate the effect of BC56 on cucumber photosynthesis under salt stress, several chlorophyll fluorescence parameters were determined. The results showed that the values of ‘ABS/CSm’, ‘TR0/CSm’, and ‘ET0/CSm’ of BC56-inoculated seedlings were 1.19-, 1.22- and 1.52-fold higher than those of BC56-uninoculated seedlings, respectively (Fig. 2a), indicating that BC56 improved the absorption and transmission capacity of light energy in cucumber.
The SPAD values of both BC56-uninoculated and BC56-inoculated seedlings showed a decreasing trend under salt stress, and that of BC56-inoculated seedlings decreased more slowly than that of BC56-uninoculated seedlings. The leaf SPAD value of BC56-inoculated seedlings (27.69) was significantly higher than that of BC56-uninoculated seedlings (24.65) at 14 DAI (Fig. 2b). It can be speculated that BC56 was able to improve the photosynthesis of cucumber leaves under salt stress.
To investigate whether BC56 affects osmotic regulation, ion balance and phytohormone regulation in cucumber under salt stress, TSS, Pro, Na+, K+, and ABA levels were measured. Under salt stress, the TSS content in BC56-inoculated seedlings was 1.36-fold as much as that in BC56-uninoculated seedlings while there was no significant difference in proline content between the two groups (Fig. 2c and 2d). The ABA content in BC56-inoculated seedlings was 0.59-fold that in BC56-uninoculated seedlings (Fig. 2e). Although there was no significant difference in Na+ and K+ content between BC56-uninoculated and BC56-inoculated seedlings, the K+/Na+ ratio in BC56-inoculated seedlings was significantly lower than that in BC56-uninoculated seedlings in roots (Fig. 2f, 2g and 2h). These results inferred that the salt stress of cucumber alleviated by BC56 might be involved in sugar and endogenous hormone metabolism as well as K+ and Na+ homeostasis.
3.5. Effects of B. cereus BC56 on activity of antioxidant enzymes
To investigate the effect of BC56 on antioxidant regulation in cucumber under salt stress, the activity of some antioxidant enzymes was determined. Under salt stress, the activity of POD and GR was 1.17-fold and 2.59-fold higher, respectively, in BC56-inoculated seedlings than in BC56-uninoculated seedlings (Fig. 3a and 3b). No significant difference in CAT and SOD activity was observed between BC56-uninoculated seedlings and BC56-inoculated seedlings (Fig. 3c and 3d). Thus, BC56 increased the activity of antioxidant enzymes under salt stress.
3.6. Transcriptomic analysis of cucumber in response to BC56 inoculation and salt stress
After eliminating low-quality reads and performing quality checks, 489 million high-quality, clean reads were obtained. Principal Component Analysis (PCA) showed that the samples were clustered into three major groups. In addition, CK and BC56 + Salt samples were closely clustered, and the gene expression level in three samples was similar to each other (Fig. 4c and S3), indicating that the sample reproducibility is qualified.
Gene expression was normalized using the fragment per kilobase of transcript per million mapped reads (FPKM) method. A total of 2,547 significant differentially expressed genes (DEGs) with log2-fold change ranging from ≥ + 1.0 to ≤ − 1.0 and P < 0.05 were identified from the comparative analysis of the four datasets (Salt vs CK; BC56 + Salt vs CK; BC56 + Salt vs Salt; BC56 vs CK). There were 845 DEGs (593 up-regulated and 863 down-regulated) in Salt vs CK, 1048 DEGs (494 up-regulated and 554 down-regulated) in BC56 + Salt vs CK, 845 DEGs (593 up-regulated and 252 down-regulated) in BC56 + Salt vs Salt, 1456 DEGs (350 up-regulated and 403 down-regulated) in BC56 vs CK (Fig. 4a). The Venn diagram showed that 953 genes were unique DEGs between Salt vs CK and BC56 + Salt vs CK transcriptional responses, suggesting that BC56 activated these regulatory pathways under salt stress. While a significant overlap of 337 DEGs was observed between BC56 + Salt vs CK and BC56 vs CK transcriptional response, suggesting activation of a common regulatory pathway by BC56 in the presence or absence of salt stress (Fig. 4b).
Expression levels of 9 randomly selected DEGs from each sample were detected by qPCR to further verify the validity of the DEGs characterized by transcriptomic analysis. The correlation between the log2 fold change in qPCR and RNA-seq showed a correlation value of 0.9305 for a total of 36 data points, confirming the validity of the transcriptomic data (Fig. 4d).
3.7. GO and KEGG analysis of the identified DEGs in cucumber
The GO annotation contains three main classification methods: biological process, molecular function, and cellular component. BC56 + Salt vs Salt GO enrichment analysis shows that some DEGs were enriched in response to water deprivation, salt stress, oxidative stress, and some terms related to plant hormones (auxin, abscisic acid, ethylene, salicylic acid) appeared in the GO annotation, suggesting that BC56 may mainly affect the plant hormone regulatory network of cucumber seedlings under salt stress, and other regulatory pathways including antioxidant regulation, water stress, salt stress, etc. In addition, a large number of DEGs were enriched in the chloroplasts GO term, indicating that BC56 also mainly affected the photosynthesis-related pathways (Fig. 5a).
KEGG pathway analysis of DEGs in different comparison groups revealed that in the BC56 + salt vs salt comparison group, DEGs were mainly enriched in eight pathways, namely ‘peroxisome’, ‘plant hormone signal transduction’, ‘MAPK signaling pathway–plant’, ‘metabolic pathway’, ‘biosynthesis of secondary metabolites’, ‘phenylpropanoid biosynthesis’, ‘starch and sucrose metabolism’, ‘flavonoid biosynthesis’ and ‘photosynthesis-antenna proteins’ (Fig. 5b and S4). Most of the regulatory genes were up-regulated in these pathways which were activated by BC56 under salt stress (Table S2).
3.8 B. cereus BC56 partially restored the expression of the genes altered by salt stress in cucumber
The FPKM between these four samples was clustered in the heat map (Fig. S3). Notably, one gene cluster was up-regulated by salt stress, and the expression level of this cluster in BC56 + Salt was partially restored to that of CK and the BC56-inoculated samples (Table S3). According to GO analysis, these genes were related to response to IAA, JA, ABA, ET, GA, transmembrane transport, transcriptional regulation, metabolite synthesis and metabolism, and cellular components. These results indicated that BC56 mainly activated the expression of genes involved in maintaining intracellular communication, cell structure, and cellular metabolism to promote cucumber growth, and thus alleviated the effects of salt stress.
3.9. Functional analysis of DEGs
In order to uncover the mechanism underlying BC56-mediated response to salt stress, some genes related to photosynthesis, phytohormones, peroxidase, sugar metabolism and signaling were selected to analyze their expression levels by qRT-PCR. Three genes LOC101203342 (chlorophyll a-b binding protein of LHCII type 1-like), LOC101218675 (chlorophyll a-b binding protein P4, chloroplastic-like) and LOC101211276 (chlorophyll a-b binding protein of LHCII type 1-like) are related to chlorophyll a-b binding protein, which are involved in the absorption and transfer of light energy. The expression levels of these three genes under BC56 + Salt treatment were 3.3-, 3.43- and 3.21-fold higher than those in salt-treated seedlings, respectively (Fig. 6a). This result indicated that BC56 enhances photosynthesis by promoting the absorption and transfer of light energy in cucumber under salt stress, which is consistent with the improvement of the cucumber chlorophyll content and photosynthetic capacity by BC56 (Fig. 2a and 2b).
The genes LOC101212508 (abscisic acid 8'-hydroxylase 1), LOC101212264 (protein REVEILLE 1-like isoform X1), LOC101204103 (ACC oxidase 2) and LOC101222233 (allene oxide cyclase 3, chloroplastic) are associated with ABA degradation, auxin synthesis, ethylene synthesis and jasmonic acid synthesis, respectively. The expression levels of the four genes under BC56 + Salt treatment were 1.8-, 2.18-, 1.75- and 7.7-fold higher than those in salt-treated seedlings respectively (Fig. 6b), suggesting that phytohormone plays a role in the BC56-mediated response to salt stress in cucumber.
The expression levels of LOC101217567 and LOC101218982 under BC56 + Salt treatment were 3.06- and 3.42-fold higher than those in salt-treated seedlings, respectively (Fig. 6c). The functional annotations of the two genes are ‘peroxidase N1-like’ and ‘peroxidase 45-like’ respectively, suggesting that the alleviation of salt stress by BC56 in cucumber was associated with POD activity.
BC56 also affected the expression of genes responsible for sugar metabolism and signaling in cucumber under salt stress. LOC101218588 (alkaline/neutral invertase A, mitochondrial) and LOC101210545 (inactive beta-amylase 4, chloroplastic) are related to sucrose and starch metabolism. LOC101216977 (2-oxoisovalerate dehydrogenase subunit alpha 1, mitochondrial) is involved in the response to sucrose. The expression levels of these three genes under BC56 + Salt treatment were 4.02-, 11.24- and 1.95-fold higher than those in salt-treated seedlings, respectively (Fig. 6d).