Plant growth-promoting effect of S. rochei S32 on different crops
Compared with the control group, the addition of 400-fold diluted cell-free fermentation broth significantly promoted the growth of wheat plants (Fig. S1). This effect was manifested by significantly increased plant height (24.7%) and root length (17.4%) after 15 d of treatment (Fig. 1a). The addition of cell-free fermentation broth also improved the growth of tomato plants mainly in terms of root length (Fig. S2). When treated with 100-, 200-, and 400-fold diluted cell-free fermentation broth, the root length of tomato increased by 9.7%, 40.9%, and 7.8%, respectively, compared with the control group (Fig. 1b).
Antagonistic activity of S. rochei S32 against soil-borne pathogens
S. rochei S32 had no or weak antagonistic effect against Sclerotinia sclerotiorum, F. graminearum, and F. oxysporum. Its inhibitory effect on Altenaria alternata, Macrophoma kawatsukai, Rhizoctonia cerealis, Colletotrichum orbiculare, and F. incarnatum was striking (Fig. 2). The highest inhibition rate was observed for M. kawatsukai (65.0%), followed by F. incarnatum (49.8%) and A. alternate (48.1%; Table 1).
Table 1
Growth inhibition rate of pathogenic fungi by S. rochei S32 grown on potato dextrose agar plates (28°C for 5 d)
Soil-borne pathogen | Plant disease | Growth inhibition rate (%) |
Altenaria alternata | Tobacco brown spot | 48.1 ± 4.0 |
Macrophoma kawatsukai | Apple ring rot | 65.0 ± 2.3 |
Rhizoctonia cerealis | Wheat sharp eyespot | 33.4 ± 5.1 |
Colletotrichum orbiculare | Cucumber anthracnose | 41.3 ± 2.3 |
Fusarium incarnatum | Muskmelon fruit rot | 49.8 ± 4.8 |
Note: Values are the means ± standard deviation |
Genomic information of S. rochei S32
The whole genome of S. rochei S32 contains a linear chromosome (Fig. 3) and two plasmids. The total length of the chromosome is 8,041,158 bp with a GC ratio of 72.5%. The size and GC ratio of plasmid 1 are 171,061 bp and 69.4%, respectively. The size and GC ratio of plasmid 2 are 95,155 bp and 69.4%, respectively. The general characteristics of the bacterial genome are summarized in Table 2.
A total of 7726 coding genes were predicted in the whole genome, and 7486 (96.8%) of them were annotated with specific functions. The number of genes annotated based on different databases was 338 (4.4%) in VFDB, 32 (0.4%) in ARDB, 433 (5.6%) in CAZy, 6133 (79.4%) in IPR, 2364 (30.6%) in Swiss-Prot, 3 (0.0%) in CARD, 7387 (95.6%) in NR, 1459 (18.9%) in T3SS, 5293 (68.5%) in COG, 4208 (54.5%) in GO, and 3866 (50.0%) in KEGG.
In the GO-based gene function classification, there were 6335, 1155, and 5054 genes involved in biological processes, cellular components, and molecular functions, respectively (Fig. 4). The three functional categories were further divided into 17, 3, and 11 subcategories, respectively. The main subcategories included “cellular processes” (1959) and “metabolic processes” (2179) in biological processes, as well as “binding” (1837) and “catalytic activity” (2495) in molecular functions.
Table 2
Genome characteristics of S. rochei S32
Characteristic | Chromosome | Plasmid 1 | Plasmid 2 |
Genome size (bp) | 8,041,158 | 171,061 | 95,155 |
G + C content (%) | 72.5% | 69.4% | 69.4% |
Ribosomal RNAs | 18 | – | – |
Transfer RNAs | 66 | – | – |
Small non-coding RNAs | 51 | – | – |
Coding sequences | 7437 | 180 | 109 |
Prophage | 8 | 1 | 1 |
CRISPR | 5 | 10 | 1 |
GenBank accession No. | CP133098 | CP133099 | CP133100 |
Prediction of gene functions in plant growth promotion
Phytohormone biosynthesis. KEGG pathway enrichment analysis revealed the biosynthesis of L-TRP and indole (two precursors for IAA biosynthesis) via the TRP biosynthesis pathway in S. rochei S32. Among the common pathways of IAA biosynthesis with TRP as a substrate, the indole-3-acetamide (IAM) pathway, the indole-3-pyruvate (IPA) pathway, the tryptamine (TAM) pathway, and the indole-3-acetaldehyde oxime (IAOx) pathway were all incomplete in the TRP metabolic pathway (Figs. S3, S4). However, only the iaam gene was missing in the IAM pathway (Fig. S3), so IAA was more likely to be synthesized via the IAM pathway at the gene level.
Biosynthesis of fungal cell wall hydrolases. The results of CAZy enzyme activity analysis showed the presence of numerous genes encoding fungal cell wall hydrolases in the chromosome of S. rochei S32. There were 185 genes for glycoside hydrolases, 113 genes for glycosyl transferases, 30 genes for carbohydrate esterases, 132 genes for carbohydrate-binding modules, 9 genes for auxiliary activities, and 12 genes for polysaccharide lyases. Genes encoding the main enzymes that hydrolyze fungal cell walls, including chitinase, cellulase, and β-1,3-glucanase, were predicted (Fig. 5).
Secondary metabolite biosynthesis. Based on the antiSMASH database, a total of 30 biosynthetic gene clusters (BGCs) were predicted in the chromosome of S. rochei S32. Twelve BGCs had > 80% similarity to known clusters with high confidence (Table 3). Among them, the products of six BGCs were identified: hopene, geosmin, albaflavenone, ectoine, isorenieratene, and 7-prenylisatin. Despite the unknown products of the remaining BGCs, they showed high similarities to those encoding borrelidin, lipopeptide, streptothricin, candicidin, and desferrioxamine B/E antibiotics. These results indicate that S. rochei S32 has the potential to produce antibiotics targeting nucleic acids, macrolides, and peptides, as well as siderophores. Other BGCs were not annotated in available databases or had low similarity to known natural BGCs, suggesting that S. rochei S32 is likely to synthesize a wide range of unknown natural products.
Table 3
Twelve predicted gene clusters of secondary metabolites in the genome of S. rochei S32
Region | Cluster category | Length (bp) | Metabolite type | Metabolite name | Cluster blast similarity (%) |
1 | Lanthipeptide-class-iii | 20,781 | Lanthipeptide | SAL-2242 | 88% |
2 | T1PKS/NRPS | 97,198 | Polyketide | Borrelidin | 81% |
3 | Terpene | 25,919 | Terpene | Hopene | 100% |
4 | NRPS | 79,373 | Non-ribosomal peptide | Lipopeptide 8D1-1/8D1-2 | 84% |
5 | Terpene | 19,913 | Terpene | Geosmin | 100% |
6 | Terpene | 20,630 | Terpene | Albaflavenone | 100% |
7 | Siderophore | 11,772 | Other | Desferrioxamine B/E | 83% |
8 | Ectoine | 10,398 | Other | Ectoine | 100% |
9 | NRPS-like | 41,145 | Non-ribosomal peptide | Streptothricin | 95% |
10 | hglE-KS/T1PKS/NRPS-like/NRPS | 208,122 | Polyketide | Candicidin | 95% |
11 | Terpene | 25,581 | Terpene | Isorenieratene | 100% |
12 | Indole | 21,178 | Other | 7-Prenylisatin | 100% |
Note: T1PKS, Type I polyketide synthases; NRPS, Non-ribosomal peptide synthetase; and KS, Heterocyst glycolipid synthase-like PKS. |
Phytohormone and metabolome profiles
Thirty-nine phytohormones were detected in the culture filtrate of S. rochei S32 (Fig. 6). The phytohormones were classified into seven categories: auxins, cytokinins (CKs), jasmonates (JAs), gibberellins (GAs), salicylic acid (SA), strigolactones (SLs), and abscisic acid (ABA). Additionally, 2205 secondary metabolites were detected. Among the top 200 most frequently detected metabolites (Fig. 7), “benzene and substituted derivatives” (28), “heterocyclic compounds” (25), “organic acid and its derivatives” (22), “small peptides” (19), and “sugars” (15) were relatively abundant. These metabolites included biologically active compounds, such as phytosphingosine, fosfomycin calcium, aspulvinone E, haloprogin, acivicin, nordihydrocapsiate, and corynebactin (Table S1).
Mechanisms of plant growth promotion
S. rochei S32 grew normally on the organophosphorus-solubilizing medium, inorganic phosphorus-solubilizing medium, and potassium-solubilizing medium, despite no halos (Fig. 8a). Normal bacterial growth also occurred on the nitrogen-free medium after three subcultures (Fig. 8b). These results indicate that S. rochei S32 has the ability to fix atmospheric nitrogen, but not to solubilize phosphorus or potassium. Poor colony growth was observed on the ADF medium with ACC as the sole nitrogen source (Fig. 8c), indicating no production of ACC deaminase. However, siderophores were produced on CAS plates, as indicated by the halo with a distinct orange color around colonies (Fig. 8d).