Front. Microbiol.Frontiers in MicrobiologyFront. Microbiol.1664-302XFrontiers Media S.A.10.3389/fmicb.2022.891245MicrobiologyOriginal ResearchVolatile Organic Compounds of Streptomyces sp. TOR3209 Stimulated Tobacco Growth by Up-Regulating the Expression of Genes Related to Plant Growth and DevelopmentHeYuxi1GuoWenyu12PengJieli1GuoJinying2MaJia1WangXu1ZhangCuimian1JiaNan1WangEntao3HuDong1*WangZhanwu1*1Institute of Agro-Resources and Environment/Hebei Fertilizer Technology Innovation Center, Hebei Academy of Agriculture and Forestry Sciences, Shijiazhuang, China2School of Landscape and Ecological Engineering, Hebei University of Engineering, Handan, China3Departamento de Microbiología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Mexico City, Mexico
Edited by: Reiner Rincón Rosales, Tuxtla Gutierrez Institute of Technology, Mexico
Reviewed by: Laith Khalil Tawfeeq Al-Ani, Universiti Sains Malaysia, Malaysia; Mohsen Mohamed Elsharkawy, Kafrelsheikh University, Egypt
*Correspondence: Dong Hu donghu1983@163.comZhanwu Wang zhanwuw@126.com
This article was submitted to Microbe and Virus Interactions with Plants, a section of the journal Frontiers in Microbiology
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.
To investigate the mechanism underlying the plant growth-promoting (PGP) effects of strain Streptomyces sp. TOR3209, PGP traits responsible for indoleacetic acid production, siderophore production, and phosphate solubilization were tested by culturing the strain TOR3209 in the corresponding media. The effects of volatile organic compounds (VOCs) produced by the strain TOR3209 on plant growth were observed by co-culturing this strain with tobacco seedlings in I-plates. Meanwhile, the effects of VOCs on tobacco gene expression were estimated by performing a transcriptome analysis, and VOCs were identified by the solid-phase micro-extraction (SPME) method. The results showed positive reactions for the three tested PGP traits in the culture of strain TOR3209, while the tobacco seedlings co-cultured with strain TOR3209 revealed an increase in the fresh weight by up to 100% when compared to that of the control plants, demonstrating that the production VOCs was also a PGP trait. In transcriptome analysis, plants co-cultured with strain TOR3209 presented the highest up-regulated expression of the genes involved in plant growth and development processes, implying that the bacterial VOCs played a role as a regulator of plant gene expression. Among the VOCs produced by the strain TOR3209, two antifungal molecules, 2,4-bis(1,1-dimethylethyl)-phenol and hexanedioic acid dibutyl ester, were found as the main compounds. Conclusively, up-regulation in the expression of growth- and development-related genes via VOCs production is an important PGP mechanism in strain TOR3209. Further efforts to explore the effective VOCs and investigate the effects of the two main VOCs in the future are recommended.
RNA sequencingSPMEStreptomycesplant growth-promoting rhizobacteria (PGPR)volatile organic compounds (VOCs)National Key Research and Development Program of China10.13039/501100012166Natural Science Foundation of Hebei Province10.13039/501100003787Introduction
In the context of population growth and rapid economic development, chemical fertilizers/pesticides have been used in agriculture to increase production, while the excessive use of these chemicals has caused serious negative impacts on the agricultural ecosystems (Vejan et al., 2016). To reduce the application of chemical fertilizers/pesticides while maintaining high production, environment-friendly biofertilizers and biopesticides have been explored (Parewa et al., 2021). In the development of biofertilizers and biopesticides, microorganisms and their metabolites have received a lot of attention, since their interactions with the plants could improve the growth and production (Saharan and Nehra, 2011; Bhattacharyya and Jha, 2012) by various mechanisms, as mentioned in the following sections. In general, all the microbes with plant growth-promoting (PGP) abilities are called PGP microorganisms; among them, the bacteria colonized in the rhizosphere, so-called PGP rhizobacteria (PGPR), have been studied extensively and applied widely (Basu et al., 2021; Mohanty et al., 2021).
Many studies on PGPR have been focused on the capacities of the bacteria to promote plant health, such as growth-enhancing potential, stress tolerance, and the production of useful biologically active secondary metabolites (Spaepen and Vanderleyden, 2011). PGPR can stimulate plant growth by various mechanisms, such as nitrogen fixation (Nassar et al., 2005), phosphate dissolution, siderophore production, antibiotic secretion, and production of plant beneficial volatile organic compounds (VOCs) (Oleńska et al., 2020). They can also promote nutrient availability, nutrient uptake (Pailan et al., 2015), and produce plant hormones, such as auxins (Lee et al., 2012), cytokinins (Lugtenberg and Kamilova, 2009), abscisic acid (Sgroy et al., 2009), and gibberellins (Khan et al., 2014). These abilities can significantly help and improve the growth of their host plants. Therefore, the PGP microorganisms, such as Rhizobium, Bacillus, Pseudomonas, Streptomyces, and Trichoderma, and their metabolites have attracted increasing attention in the production of biological fertilizers and/or biopesticides (Atieno et al., 2020).
The members of the genus Streptomyces are Gram-positive, spore-forming bacteria that are widely distributed in the soil and other habitats, and they are well known for their ability to produce a large number of bioactive secondary metabolites (antibiotics) (Quinn et al., 2020), but rarely act as phytopathogens (Li et al., 2019). Recently, Streptomyces strains have been utilized as PGPR for several crops, including wheat (Akbari et al., 2020), tomato (Hu et al., 2020), chickpea (Gopalakrishnan et al., 2015), sorghum (Gopalakrishnan et al., 2013), and rice (Suárez-Moreno et al., 2019). The PGP effects of these Streptomyces strains have been attributed to their production of siderophores and indole-3-acetic acid, nitrogen fixation, and phosphate solubilization, which make them potential biofertilizers (Romano-Armada et al., 2020). In addition, some of the Streptomyces strains could improve plant growth and production by increasing the plant resistance to phytopathogens and to the abiotic stresses, which is related to their production of chitinolytic enzymes (inhibiting insects and fungi), induction of non-specific plant resistance, production of bioactive compounds (VOCs), and antibiotic/siderophore production (Romano-Armada et al., 2020). From Streptomyces species, more than 120 volatile substances belonging to alkenes, esters, alcohols, and ketones have been detected (Schöller et al., 2002), and various functions of these VOCs have been reported in agriculture, including PGP effects and phytopathogen control (Olanrewaju and Babalola, 2019). However, the mechanism by which the VOCs produced by Streptomyces promote plant growth and development is unclear.
Volatile organic compounds are secondary metabolites produced by many fungi and bacteria (Korpi et al., 2009), which are believed to act as signal molecules in interspecific and intraspecific interactions and in intercellular communication (Schmidt et al., 2015). Therefore, the in-depth study on the effects of VOCs on plant gene expression might provide some insight into the interactions between microbes and plants, as well as into the exploration of biofertilizers/biopesticides. Previously, Streptomyces sp. TOR3209 was isolated from the rhizosphere of tomato as a PGP bacterium (Hu et al., 2012), and the inoculation of this strain could significantly promote the growth and enhance the yield of tomato and other crops by changing the structure of the microbial community in the rhizosphere (Hu et al., 2020). However, its capacity to produce VOCs and whether the VOCs play a role in inducing the PGP effects have not been revealed. To check if the VOC production is one of the mechanisms of PGP effects in the strain TOR3209 and to estimate the effects of plant metabolism, we evaluated the PGP ability of VOCs produced by Streptomyces sp. TOR3209 and investigated the molecular composition of these VOCs. We also identify the important molecular events involved in the rhizosphere bacteria–plant interaction by transcriptome analysis.
Materials and MethodsIncubation of Bacteria and Inoculant Preparation
The strain Streptomyces sp. TOR3209 originally isolated from the rhizosphere of tomato plants (Hu et al., 2012) was used in this study, considering its role in the promotion of growth in tomatoes and other plants (Hu et al., 2020). To prepare the inoculant, the strain TOR3209 was massively inoculated on plates containing the Gause's No.1 medium (soluble starch 20 g, KNO3 1.0 g, K2PO4 0.5 g, MgSO4·7H2O 0.5 g, NaCl 0.5 g, FeSO4·7H2O 0.01 g, Agar 20 g, in 1 L of distilled water, pH 7.2–7.4) and incubated at 30°C for 24 h. Then the bacterial biomass was scraped off the plate, suspended in sterilized 10 mM MgCl2, and adjusted to 108–109 CFU/ml based on both optical density (OD600 = 1.0) and plate count. To test the effects of VOCs, strain TOR3209 was cultured overnight in Gause's No.1 liquid medium at 30°C with shaking at 180 rpm.
Plant Growth-Promoting Test
To estimate the PGP effects of strain TOR3209 or its VOCs, tobacco (Nicotiana benthamiana Domin) was used as a model plant, which showed a positive response to inoculation of strain TOR3209 in our previous study. The seeds of tobacco were purchased from Huayueyang Biotechnology (Beijing) Co., Ltd.
Growth Promotion of Plants by Inoculation
The seeds were surface sterilized with 95% ethanol for 30 s, 2% (v/v) sodium hypochlorite for 5 min, and then washed with sterile distilled water (SDW) three times as described previously (Hu et al., 2020). The surface-sterilized seeds were germinated on 0.5 × Murashige and Skoog (MS) medium (PhytoTechnology, United States) supplemented with 1.5% (w/v) sucrose and 1.8% agar (w/v). The MS agar plates with seeds were kept at 4°C for 12 h for vernalization and then kept in a plant growth chamber for germination at 25°C under a long-day photoperiod (day/night: 16/8 h) with a light intensity of 100 μmol m−2 s−1. After 5 days of incubation, the seedlings were transferred to a 24-well plate containing the MS agar and incubated in the plant growth chamber with 16 h light/8 h dark conditions. Five days later, tobacco plants in the plate were inoculated with 20 μl of inoculant (108−9 CFU/ml) around the root, and the negative control plants were inoculated with 20 μl of sterilized 10 mM MgCl2 solution. The tobacco seedlings were harvested 5 days after the inoculation, and total fresh weight was measured. The experiment was repeated three times, and 10 replications (seedlings) for each treatment were included each time.
Growth Promotion of Plants by VOCs of Strain TOR3209
Previously, the VOCs of bacteria have been reported as novel agents in promoting plant growth (Audrain et al., 2015). To evaluate the growth-promoting effect of VOCs on the tobacco plants, specialized plastic Petri dishes (I-plates with 9 cm in diameter) were used, which were divided by a partition into two equal chambers. In these I-plates, Streptomyces sp. TOR3209 and tobacco plants could be co-cultivated without direct contact. One chamber of the I-plate contained half-strength MS agar, and the other chamber contained Gause's No.1 medium with a sterile filter paper disk (8 mm in diameter) on its surface. After 5 days of germination as mentioned earlier, five seedlings were transferred to the chamber containing MS agar. After 3 days of plant growth, 20 μl of overnight incubated strain TOR3209 (see section Incubation of Bacteria and Inoculant Preparation) or sterilized distilled water (SDW) for treatment and control (CK), respectively, were applied drop-wise onto the filter paper disk in the chamber containing Gause's No.1 medium. Then, the I-plates were sealed with parafilm and cultured in the growth chamber under the day/night cycle of 16/8 h, with a light intensity of 100 μmol m−2 s−1 at 25°C. The fresh weight of each plant was measured after 1 week of growth. The experiment was repeated independently for three times, and 10 replications (seedlings) per treatment were included each time.
Characterization of Plant Growth-Promoting Traits
To estimate the possible PGP mechanism of strain TOR3209, the PGP traits mentioned by Romano-Armada et al. (2020), such as the formation of indole-3-acetic acid (IAA), production of siderophores, and solubilization of phosphate, were measured for this strain.
Auxin Production Assay
The production of IAA was quantified by the colorimetric method using the Salkowski reagent (50 ml of 35% HClO4 and 1 ml of 0.5 M FeCl3) (Glickmann and Dessaux, 1995). The strain TOR3209 was activated in LB (Difco) broth for 24 h, and then a full loop of the culture was inoculated in 5 ml of the LB broth for culturing at 30°C with shaking at 180 rpm for 16 h. Then, 20 μl of bacterial suspension was transferred to 5 ml of LB broth supplemented with 0.1% L-tryptophan. After 16 h of incubation, 1 ml of the bacterial culture was centrifuged at 8,000 rpm for 10 min to remove the cells, and 300 μl of the supernatant was mixed with 600 μl of Salkowski reagent. The mixture was left in the dark for 30 min to develop the color (Phi et al., 2010). After the reaction, the absorbance of the mixture was measured at 540 nm using a UV–visible spectrophotometer (MAPADA, China). The concentration of IAA in the culture was calculated (Sigma, USA) based on the standard curve of IAA (serial dilutions from 0 to 50 μg/ml). The experiment was repeated three times, and 10 replications were included for each treatment each time.
Siderophore Production
In this analysis, siderophore production was checked as described previously (Neilands, 1987). Briefly, the strain TOR3209 was cultured on Chrome Azurol S (CAS) agar (Pérez-Miranda et al., 2007) and incubated at 30°C for 7 days. The formation of a yellow ring around the colony indicates the production of siderophore. To determine the siderophore type, a FeCl3 test was performed (Chowdappa et al., 2020). Briefly, the strain was cultured in 5 ml of CAS broth at 30°C with rotation for 5 days, and 1 ml of supernatant of the culture obtained by centrifugation (12,000 rpm for 1 min) was mixed with 2 ml of FeCl3 solution (2%, w/v). Then the absorption was scanned from 300 to 600 nm with a spectrophotometer (MAPADA, China). The occurrence of a peak between 420 and 450 nm and a peak at 495 nm demonstrate the presence of hydroxamate and catecholate types of siderophores, respectively. The CAS broth without inoculation was used as a blank to detect the siderophore type. The experiment was repeated three times, and 10 replications were included for each treatment each time.
Phosphate Solubilization
A single colony of strain TOR3209 grown on Gause's No. 1 medium at 30°C for 3 days was inoculated in an organic phosphorus medium (Jiwen et al., 2011) and an inorganic phosphorus (PKO) medium (Nautiyal, 1999), respectively, to measure the phosphorus solubilization. The inoculated plates were cultured at 30°C for 12 days, and the formation of a transparent ring surrounding the colony demonstrated the dissolution of phosphorus. Since the organophosphorus component in the medium is prepared using egg yolk (lecithin), it may undergo denaturation during long-term culture. Therefore, a control plate was set up for the organophosphorus solubilization to avoid false-positive results. The experiment was repeated three times, and 10 replications were included for each treatment each time.
Expression of Plant Genes Responding to the VOCs
After 7 days of exposure to VOCs (co-cultured with strain TOR3209 in I-plates), as mentioned in Section Growth Promotion of Plants by VOCs of Strain TOR3209, the aerial parts (leaves and stems) of tobacco seedlings were collected from three I-plates of treatment or control in triplicate (15 plants in total for treatment or control). Experiments for treatments and controls were repeated three times per group. Total RNA was extracted from the plant samples using the TIANGEN Magnetic Tissue/Cell/Blood Total RNA Kit (DP762- T1A), according to the manufacturer's guidelines. The concentration of RNA in the extracts was determined using a NanoDrop spectrophotometer (Thermo Scientific), and the quality and integrity of RNA in the extracts were tested by agarose gel electrophoresis (Laila et al., 2017). For each extract, 3 μg of RNA was used as input material to construct the sequencing library using the NEBNext Ultra II RNA Library Prep Kit for Illumina (#E7775). The procedure described by Laila et al. (2017) was followed. Briefly, mRNAs were purified from total RNAs using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in an Illumina proprietary fragmentation buffer. The first strand of cDNA was synthesized using random oligonucleotides and SuperScript II. Synthesis of the second strand of cDNA was subsequently performed using DNA polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities, and the enzymes were later removed. After adenylation of the 3′ ends of the DNA fragments, Illumina PE adapter oligonucleotides were ligated to prepare for hybridization. To select the cDNA fragments of the preferred size (200 bp), the library fragments were purified using the AMPure XP system (Beckman Coulter, Beverly, CA, USA). DNA fragments with ligated adaptor molecules on both ends were selectively enriched using Illumina PCR Primer Cocktail in a 15-cycle PCR reaction. Products were purified (AMPure XP system) and quantified by performing the Agilent High Sensitivity DNA assay on a Bioanalyzer 2100 system (Agilent). The sequencing library was then sequenced commercially on a Novaseq platform (Illumina) by Shanghai Personal Biotechnology Cp. Ltd. Filtration of high-quality reads by eliminating low-quality reads (Q <20), removal of adapter, assembling of the transcriptome, and quantification of transcript abundances were performed as described by Howlader et al. (2020).
The gene expression levels between the samples subjected to VOC treatment and the control were estimated by using fragments per kilobase of transcript per million fragments mapped (FPKM) values, as described by Howlader et al. (2020). In the transcription group, we generally took into account FPKM > 1 when the gene was expressed. Pearson correlation coefficients were used to determine the level of correlation between the expression of genes in the samples to test the reliability of the experiments and whether sample selection was reasonable (van Ruissen et al., 2005). The log2 fold change was calculated by dividing the FPKM values of the treatment groups by the values of the control group. High-quality reads were processed using the Perl script, and the differentially expressed genes (DEGs) were identified using the edgeR package (Boscari et al., 2013), taking a fold change ≥ 2 as the threshold value.
The cluster analysis of FPKM-normalized expression values was carried out, as mentioned by Howlader et al. (2020), to differentiate the up- and down-regulated transcripts between the VOC treatment and the control samples, by using the cluster software version 3.0 and the Java TreeView to visualize the dendrogram. For GSEA, the KEGG pathway collection (Homo sapiens (human) Release 53.0) (Kanehisa and Goto, 2000) was used to define the gene sets using the GSEABase package in Bioconductor/R version 2.9.2 (Gentleman, 2004). Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were used to characterize the differentially expressed genes, and GO functional annotations were obtained from the non-redundant annotation results using the Blast2GO software (Conesa et al., 2005).
Extraction and Identification of VOCs
The strain TOR3209 was activated by culturing in Luria-Bertani (LB) broth (NaCl 10 g, tryptone 10 g, yeast extract 5 g, distilled water 1 L, pH 6.8–7.2) for 16 h at 30°C with shaking at 180 rpm. Then 40 μl of the culture was dropped on the LB agar in bottles (18 × 18 × 3 cm3) and cultured for 3 days at 30°C. The VOCs produced in the bottles were extracted and detected commercially by the solid-phase micro-extraction (SPME) method (Agilent, 6890N-5975B MSD, USA) (Gu et al., 2007). Briefly, SPME fibers were introduced into the atmosphere of the bottle to absorb the VOCs. Then, the SPME fibers with VOCs were incubated at 70°C for 20 min and were desorbed at 250°C for 5 min in the injection port using a DB-WAX column (30 m × 0.25 mm × 0.25 μm - polar, J&W Scientific, Folsom, CA). The initial oven temperature was 40°C (held for 2 min), ramped at 2°C/min to 220°C, and finally ramped at 20°C/min to 240°C (held for 10 min). GC–MS runs were performed for a total of 103 min. Compounds were identified by comparing their retention indices (RI) determined relatively to a homologous series of pure n-alkane standards or values reported in the literature (Jennings and Shibamoto, 1980). Fragmentation patterns in the mass spectra were also compared with those available in the Library of the National Institute of Standards and Technology (NIST05). The effective compounds were identified and estimated according to the methods described in the literature (Li et al., 2018; Chowdhary and Sharma, 2020). Experiments for treatments and controls were repeated three times per group.
Statistical Analysis
The data obtained from the plant growth promotion test were statistically analyzed by ANOVA, and the mean values of the treated samples were separated by the least significant difference (LSD) test at p < 0.05 using SPSS version 19.0 (SPSS Inc., USA). Experiments were done in 10 replicates, and each experiment was repeated for three times. The results of repeated trials of each experiment were similar. Therefore, one representative trial for each experiment is reported in the Section Results.
ResultsEffects of Strain TOR3209 Inoculation and VOCs on the Growth of Plant
In the growth tests (Supplementary Table 1), a tendency toward the promotion of growth was observed in the tobacco seedlings. The seedlings grown in MS medium inoculated with strain TOR3209 showed a 20% increase in their total fresh weight after 5 days of inoculation, in comparison with the controls. However, the increase in growth was not significant according to the LSD test. In I-plates, the total fresh weight of tobacco seedlings co-incubated with strain TOR3209 was significantly elevated up to 100% (Figure 1).
Comparison of effects of strain TOR3209 VOCs on tobacco growth. (A) Comparison of fresh weight between the CK seedlings and treatment seedlings; (B) Photographs showing the growth promotion on growth of tomato seedlings by strain TOR3209 co-culture (bellow) in comparison with CK control (upper). * The same letter followed the values in the same column means no significant difference based on the least significant difference (LSD) test (p < 0.05).
PGP Traits of <italic>Streptomyces</italic> sp. TOR3209
The analysis of IAA production showed that strain TOR3209 produced 12.06 μg/ml of IAA in the culture. In the P solubilization test, the strain TOR3209 grew slowly and weakly on the PKO medium. After 2–7 days of incubation, the colonies gradually expanded, and the phosphate solubilizing ring gradually appeared, but the ring was not clear and remained unchanged after 7 days. On organophosphorus medium, strain TOR3209 presented a more obvious dissolving effect. After 5 days of incubation, the solubilization ring was very obvious (SI = 1.6), and it continuously expanded until the entire plate became almost transparent. On CAS medium, colonies of strain TOR3209 began to grow after 1 week of incubation with the appearance of yellow circles, indicating the production of siderophores. Based on the result of the FeCl3 test, the siderophores produced by the strain TOR3209 were hydroxamate type. Detailed results for PGP traits of Streptomyces sp. TOR3209 are available in Supplementary Figure 1.
Transcriptomal Analysis of Tobacco Plants Responding to VOCs
In transcriptomal analysis, 38,254 genes were detected in the plants co-cultured with strain TOR3209, and 37,636 genes were detected in the control plants. Analysis of RNAseq data revealed that 3,191 differentially expressed transcripts (DEGs) were induced by the VOCs of strain TOR3209, including 2,155 up-regulated and 1,036 down-regulated transcripts (Figure 2). The most important (Log2 fold change >5 or < -4) up-regulated (52) and down-regulated (23) transcripts are listed in Table 1. Among them, UDP-glycosyltransferase, glutamate receptor, RING/U-box protein with C6HC-type zinc finger, and leucine-rich repeat receptor-like protein kinase family protein were the highest up-regulated DEGs in the VOC-treated tobacco seedlings. Chaperone protein ClpB, MYB-like transcription factor family protein, and 17.3 kDa class II heat shock proteins were significantly repressed by the VOC treatment.
Volcanic diagrams of differential expression genes. The horizontal coordinates are log2 (fold change) and the ordinates are -log10 (p-value). The two vertical dashed lines are 2 times the difference threshold, and the horizontal dashed lines are p-value 0.05 thresholds. A red dot indicates the upregulated genes group, a blue dot indicates the downregulated genes group, and a gray dot represents a non-significant difference in the expression of the genes group. The scale was fixed considering the minimum and maximum log2 values (between −7.74 and 12.33) included in the Supplementary Tables 1, 2.
Most important differentially expressed genes in tobacco plants inoculated with TOR3209.
Gene ID
Log2 fold change
Gene description
Up-regulated genes
Niben101Scf00107g03008
12.33288866
Cysteine-rich venom protein
Niben101Scf02349g03003
11.07906488
Peroxidase 4
Niben101Scf00661g00002
8.319303503
UDP-glycosyltransferase 74 F1 LENGTH = 449
Niben101Scf08670g00020
8.139389486
Glutamate receptor 2.9
Niben101Scf07123g01015
7.759232431
Leucine-rich repeat receptor-like protein kinase family protein LENGTH = 1,045
Niben101Scf01463g08009
7.667222902
Calmodulin binding protein-like LENGTH = 494
Niben101Scf01942g04001
7.576358542
WRKY transcription factor 44
Niben101Scf02188g00002
7.509948729
WRKY transcription factor 1
Niben101Scf08526g02006
7.503570319
WRKY transcription factor 1
Niben101Scf12084g00003
7.382594655
WRKY transcription factor 1
Niben101Scf01400g00014
7.377702267
Thaumatin-like protein
Niben101Scf02115g04004
7.269910522
Eukaryotic aspartyl protease family protein LENGTH = 474
Niben101Scf12266g09001
7.039494994
Tetratricopeptide repeat (TPR)-like superfamily protein LENGTH=305
Niben101Scf02406g00002
7.031384344
GDSL esterase/lipase
Niben101Scf07589g01006
7.006395689
conserved hypothetical protein (Ricinus communis) gb|EEF44100.1| conserved hypothetical protein (Ricinus communis)
myb-like transcription factor family protein LENGTH = 233
Niben101Scf04040g08015
−5.793913087
17.3 kDa class II heat shock protein
Niben101Scf06977g00014
−5.659566414
Chaperone protein ClpB
Niben101Scf04040g09011
−5.597721532
17.3 kDa class II heat shock protein
Niben101Scf10384g00015
−5.479676566
17.6 kDa class I heat shock protein
Niben101Scf08675g00002
−5.47869793
BAG family molecular chaperone regulator 6
Niben101Scf03114g03011
−5.377817272
Chaperone protein HtpG
Niben101Scf00298g00001
−5.316649482
Ethylene-responsive transcription factor 1
Niben101Scf05378g03025
−5.236085551
ATP-dependent zinc metalloprotease FtsH
Niben101Scf02175g06018
−5.185219051
BAG family molecular chaperone regulator 6
Niben101Scf04040g09016
−5.028860367
17.3 kDa class II heat shock protein
Niben101Scf10306g00005
−4.958344447
17.3 kDa class II heat shock protein
Niben101Scf02552g00013
−4.659113231
Subtilisin-like serine protease 2 LENGTH = 764
Niben101Scf10384g00016
−4.602974453
17.6 kDa class I heat shock protein
Niben101Scf02124g01003
−4.593158829
Heat shock 70 kDa protein
Niben101Scf07790g00006
−4.426069226
COBRA-like protein 10
Niben101Scf04410g09016
−4.419070922
Chaperone protein ClpB
Niben101Scf02267g00013
−4.344794817
17.3 kDa class II heat shock protein
Niben101Scf01412g00015
−4.328983756
17.7 kDa class II heat shock protein
Niben101Scf02913g00028
−4.256421593
Alpha-galactosidase
Niben101Scf09552g01024
−4.24630503
Protein kinase superfamily protein with octicosapeptide/Phox/Bem1p domain LENGTH = 1,054
Niben101Scf00117g02017
−4.1759724
Octic, putative isoform 1 (Theobroma cacao)
In the functional annotation, DEGs with the largest increase in expression were those related to cell recognition, pollination, pollen–pistil interaction, multi-multicellular organism process, and recognition of pollen based upon the GO enrichment analyses (Figure 3; Table 1). Based on the KEGG pathway, the most significant up-regulated transcripts by VOCs were categorized as those related to glutathione metabolism, taurine and hypotaurine metabolism, and glycosphingolipid biosynthesis pathways (Figure 3). The expression profile of TF genes displays the upcoming transcription actions regulated by these genes. The largest TF families detected in our study were basic helix-loop-helix (bHLH) DNA-binding superfamily protein, WRKY, NAC domain-containing protein (NAC), ethylene-responsive factor (ERF), and zinc finger protein genes (C2H2) (Figure 4).
Gene ontology (GO) (A) and KEGG pathway-rich analysis (B). Based on the GO-rich and KEGG results, the extent of richness is measured by the rich factor, false discovery rate (FDR) values, and the number of genes that were collected into this GO term or KEGG pathway. Rich factor refers to the ratio of the number of different genes collected in the GO term or KEGG pathway to the number of genes annotated.
Transcription factor family. The horizontal coordinates are different transcription factor families, and the ordinates are the number of genes that fall into that transcription factor family.
Analysis of VOCs Produced by <italic>Streptomyces</italic> sp. TOR3209
In GC–MS analysis, 48 compounds were detected, among which eight (accounting for 15.22% of the total VOCs) were reported as PGP effective VOCs, and 14 (accounting for 37.36% of the total VOCs) had the potential to increase the resistance of the plant to phytopathogens (Supplementary Table 2). According to their relative abundances, seven compounds could be identified as the main components, each contributing to more than 5% of the total: the pathogen-inhibiting compounds like trans-1,10-dimethyl-trans-9-decalol (9.43%), 1,4,5,6,7,8,9,9a-octahydro-1,1,7-trimethyl-[3aR-(3a.alpha.,7.alpha.,9a.beta.)]- 3a,7-methano-3aH-cyclopentacyclooctene (6.29%), and 2,4-bis(1,1-dimethylethyl)- phenol (13.61%); the growth-promoting compounds like n-hexadecanoic acid (6.53%) (Chowdhary and Sharma, 2020); and the compounds of an unknown function like (-)-.alpha.-Panasinsen (5.12%), pentadecanoic acid (5.58%), and hexanedioic acid dibutyl ester (25.58%). Among these compounds, 2,4-di-tert-butylphenol [o 2,4-bis(1,1-dimethylethyl)- phenol] and dibutyl adipate (o hexanedioic acid dibutyl ester) were the most abundant VOCs produced by the strain TOR3209 (Figure 5).
Volatile organic compounds of strain TOR3209 using GC-MS. The samples were incubated in bottle containing LB agar for 3 days. 2,4-bis(1,1-dimethylethyl)-phenol and hexanedioic acid dibutyl ester were detected.
Discussion
It is well known that some microorganisms have the potential to promote plant growth through direct and/or indirect effects. The direct effects include the production of molecules that promote plant growth (Goswami et al., 2016). In this study, the growth (total fresh weight) of tobacco seedlings was significantly increased (100%) by co-culturing with the strain Streptomyces sp. TOR3209 (Figure 1; Supplementary Table 1), demonstrating that this strain produced VOCs as PGP agents, consistent with the previous reports for other microbial VOCs (Kanchiswamy et al., 2015). These results suggest an alternative mechanism for the PGP effects of strain TOR3209, other than the previously described regulation of the rhizosphere microbial community (Hu et al., 2020).
As reported previously, the VOCs mainly react with the aerial parts of the plants (Asari et al., 2016). Considering the low amount of VOCs, it could be hypothesized that they cannot function as nutrients, but can serve as signals to regulate the production of plant genes related to growth. Indeed, the transcriptome analysis based upon the KEGG analysis (Kanehisa and Goto, 2000) in the present study revealed that VOCs of strain TOR3209 up-regulated the expression of multiple genes to improve plant growth. The highest increase in the expression was noted for the genes of UDP-glycosyltransferase, glutamate receptor, RING/U-box protein with C6HC-type zinc finger, and leucine-rich repeat receptor-like protein kinase family protein in the VOC treatment samples, which indicated the possible relationships between these genes and the increased growth. The expression of UDP-glucuronosyltransferase (UGT) in meristematic cells is correlated with mitosis and is required for the normal development of plants (Woo et al., 1999). So, the VOCs of strain TOR3209 might increase plant growth by promoting cell division. Also, glycosylation is a universally important mechanism in certain organisms to control the impact of a wide range of hydrophobic compounds that can directly derange cellular metabolism (Nebert, 1991). Therefore, VOCs of strain TOR3209 might also promote plant growth by regulating the UGT-mediated glycosylation in plants. Plant glutamate receptor-like genes (GLRs) coded amino acid receptors related to plant defense and innate immune response (Forde and Roberts, 2014). Recently, the role of GLRs in controlling the plant height has also been established (Kaur, 2020), which further explained the PGP effects of VOCs produced by strain TOR3209. U-box proteins play important myriad functions in plants, such as in growth and development, cell cycle regulation, morphogenesis, stress resistance, and regulating the innate immune response of plants to biotic stress factors. Therefore, strain TOR3209 was capable of promoting growth through VOCs that induce the expression of genes coding U-box proteins and LRR-RLKs (Liu et al., 2017), which further regulate the expression of other genes related to the plant growth and development, including the synthesis of proteins and enzymes.
Transcription factors (TFs) play essential roles in plants, as well as in all other living organisms, by controlling the expression of genes involved in various cellular processes (Han et al., 2014). Numerous processes, such as brassinosteroid signaling, ABA signaling, and axillary meristem formation, are regulated by the members of the bHLH family of TFs during the seedling development (Zhao et al., 2012). A total of 143 transcripts were identified as bHLH TF proteins in the present study, which have been identified as phytochrome-interacting factors (Toledo-Ortiz et al., 2003) and constitute the largest TF family proteins that are produced during the interaction of tobacco plants with the VOCs of strain TOR3209. These results further indicated that the VOCs produced by strain TOR3209 were capable of promoting plant growth by regulating a serial synthesis of TFs. Another large NAC TF family was involved in inducing diverse developmental events in Arabidopsis and soybean (Shamimuzzaman and Vodkin, 2013). This NAC TF family has the capacity to permit crosstalk between different pathways (Olsen et al., 2005). In our present study, 25 putative NAC genes were detected during plant growth. Furthermore, genes encoding C2H2 zinc finger protein play a significant role in the development and differentiation (Razin et al., 2012). During the interaction, 15 putative C2H2 zinc finger genes were identified. In addition, ethylene-responsive factor (ERF) TFs belonging to the Apetala2/ERF TF family were also found to be involved in plant growth. The WRKY genes coding TFs were involved in a broad range of biological processes, including diverse biotic/abiotic stress responses, and developmental and physiological processes (Birkenbihl et al., 2017). A total of 18 putative ERF and 32 putative WRKY genes were detected during plant development in the present study. So, it is clear that the VOCs of Streptomyces sp. TOR3209 could regulate diverse metabolism pathways, which further improved the growth of tobacco plants.
The significantly down-regulated genes in the samples subjected to VOC treatment were chaperone protein ClpB, MYB-like transcription factor family protein, and 17.3 kDa class II heat shock proteins. These three genes and most of the other down-regulated genes are genes that confer adaptation to environmental stress. ClpB (caseinlytic proteinase) is related to the disaggregation of protein complexes and cell adaptation to heat stress (HS) (Mishra and Grover, 2016). A MYB-like transcription factor is involved in plant abiotic stress (like salt) tolerance and phytohormone signal transduction (Zhang et al., 2018). Furthermore, the 17.3 kDa class II heat shock proteins are also stress tolerance proteins related to heat shock (Kang et al., 2021). The reason why the expression of these genes was elevated by VOC treatment is not very clear, but it is possible that the VOCs shuttled down the expression of genes related to some unnecessary cell functions, while stimulating the expression of other genes essential for plant growth.
Previously, more than 120 VOCs have been detected in Streptomyces, which have been classified into different chemical groups, such as alkanes, alkenes, esters, alcohols, and ketones (Schöller et al., 2002). According to the literature, 15.22% (eight) of the VOC molecules detected in the present study presented PGP effects, which might be the potential compounds directly responsible for the growth promotion of tobacco seedlings in the I-plates co-cultured with strain TOR3209. Meanwhile, more VOCs (37.36%, 14 molecules) were related to the resistance to phytopathogens; however, it is not clear whether these VOCs can also promote the growth of plants. As the most abundant compound (25.58%) in the VOCs of strain TOR3209, the effects of hexanedioic acid dibutyl ester on the plant has not been reported and its function in regulating plant growth deserves further study.
In the PGP trait analysis, 12.06 μg/ml of IAA production by strain TOR3209 was recorded, which is relatively a high level since values around 5 to 10 μg/ml are commonly reported. However, some actinobacterial strains capable of producing more than 40 μg/ml have been isolated (Vasconcellos et al., 2010). In addition, the IAA production efficiency of strain TOR3209 might be increased by optimizing the culture conditions, because the production of IAA by a Streptomyces strain was increased from 20.46 to 82.36 μg/ml by modifying the medium components (Myo et al., 2019). Therefore, strain TOR3209 could be a good stimulator of plant growth via the production of IAA, as described for other Streptomyces strains and other bacteria (Hayat et al., 2010). Another PGP trait of strain TOR3209 is to dissolve phosphate, particularly the organic phosphate, which could subsequently promote plant growth (Spaepen et al., 2007). In this study, we also proved the ability of strain TOR3209 to produce siderophore, a low-molecular-weight iron-binding compound that could help the plant chelate iron in the Fe3+ restricted soil (Bottini et al., 2004) and inhibit plant pathogens due to their competitive role (Otieno et al., 2015). Currently, many siderophore molecules, with diverse chemical structures and different bioactivities against a wide spectrum of pathogenic bacteria, have been isolated from Streptomyces species (Terra et al., 2021). Different from the previous reports that the Strteptomyces tendae Tü 901/8c and Streptomyces sp. Tü 6125 isolated from soils produced the catecholate siderophore (Fiedler et al., 2001), the strain TOR3209 produced the hydroxamate siderophore, indicating that different Streptomyces species might have distinct pathways for the synthesis of siderophore. Although the PGP features of strain TOR3209 were detected in the present study, they might not contribute to the growth of tobacco seedlings in the inoculation test in the present study, since they are not necessary for the plants grown in the sterilized medium, which is different from the situation in soils (Hu et al., 2020). However, the capacity of Streptomyces species to produce IAA and siderophore deserves further study to explore their biofertilizer and biopesticide properties.
In our present study, the promotion of plant growth was not statistically significant following the direct inoculation of strain TOR3209, but it was significant in VOC treatments (Supplementary Table 1). This was inconsistent with the previous results that inoculation with strain TOR3209 significantly promoted the growth of tomato plants (Hu et al., 2020). This discrepancy might be related to the differences in the culture conditions because the inoculated plants were cultured in soil (Hu et al., 2020), while the experiments in our study were carried out in medium or artificial substrate. In the medium, the rhizosphere microbiome is absent and the nutrients are not limited, which might also explain why the direct inoculation with the strain could not significantly improve the growth of tobacco seedlings in the present study. Furthermore, plant genotypes (Schwachtje et al., 2012) and environmental conditions (Martínez et al., 2011) also affect the inoculation effects of PGPR, and even the PGP effects of microbial siderophore could be different for distinct tobacco varieties (Robin et al., 2006).
Conclusion
The VOCs of strain TOR3209 could improve the growth of the tested plants and induce the up-regulated expression of a key group of genes associated with plant growth and development and the down-regulation of some adaptation genes. Remarkably, our results suggest that the genes related to UDP-glycosyltransferase, glutamate receptor, RING/U-box protein with C6HC-type zinc finger, and leucine-rich repeat receptor-like protein kinase family protein presented greatly up-regulated expression in the tobacco seedlings in response to the VOCs of strain TOR3209, and all these genes were involved in the plant growth and development processes. In addition, the important VOCs of strain TOR3209 were found to be 2,4-bis(1,1-dimethylethyl)-phenol and hexanedioic acid dibutyl ester (with two compounds exhibiting antifungal functions), and their other PGP effects need to be further investigated. In the long term, these new findings would facilitate the development of biofertilizers and biopesticides.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: National Center for Biotechnology Information (NCBI) BioProject database under accession number PRJNA814410.
Author Contributions
YH conceived the study and wrote the manuscript. YH and ZW designed the experiments. YH, WG, JP, JG, JM, XW, CZ, NJ, EW, DH, and ZW did the experiments. YH and WG assisted in the data analysis. YH, EW, and JG edited the article. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key Research and Development Program of China (2021YFD1901004-4), the HAAFS Science and Technology Innovation Special Project (2022KJCXZX-ZHS-3), and the Project of Natural Science Foundation of Hebei Province (C2020301047). EW was supported by the Project of Sabbatical Year and SIP20200726 authorized by IPN, Mexico.
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.
Sequencing services were provided by Personal Biotechnology Co., Ltd. Shanghai, China. The data were analyzed by using the free online platform Personalbio GenesCloud (http://www.genescloud.cn).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.891245/full#supplementary-material
ReferencesAkbariA.GharanjikS.KoobazP.SadeghiA. (2020). Plant growth promoting Streptomyces strains are selectively interacting with the wheat cultivars especially in saline conditions. Heliyon6, e03445. 10.1016/j.heliyon.2020.e0344532382669AsariS.MatzénS.PetersenM. A.BejaiS.MeijerJ. (2016). Multiple effects of Bacillus amyloliquefaciens volatile compounds: plant growth promotion and growth inhibition of phytopathogens. FEMS Microbiol. Ecol. 92, fiw070. 10.1093/femsec/fiw07027053756AtienoM.HerrmannL.NguyenH. T.PhanH. T.NguyenN. K.SreanP.. (2020). Assessment of biofertilizer use for sustainable agriculture in the Great Mekong Region. J. Environ. Manage.275, 111300. 10.1016/j.jenvman.2020.11130032871522AudrainB.FaragM. A.RyuC.-M.GhigoJ.-M. (2015). Role of bacterial volatile compounds in bacterial biology. FEMS Microbiol. Rev.39, 222–233. 10.1093/femsre/fuu01325725014BasuA.PrasadP.DasS. N.KalamS.SayyedR. Z.ReddyM. S.. (2021). Plant growth promoting rhizobacteria (PGPR) as green bioinoculants: recent developments, constraints, and prospects. Sustainability13, 1140. 10.3390/su13031140BhattacharyyaP. N.JhaD. K. (2012). Plant growth-promoting rhizobacteria (PGPR): emergence in agriculture. World J. Microbiol. Biotechnol.28, 1327–1350. 10.1007/s11274-011-0979-922805914BirkenbihlR. P.KracherB.RoccaroM.SomssichI. E. (2017). Induced genome-wide binding of three Arabidopsis WRKY transcription factors during early MAMP-triggered immunity. Plant Cell29, 20–38. 10.1105/tpc.16.0068128011690BoscariA.FerrariniA.GiudiceJ. dVenturiniL.DelledoneM.. (2013). Mapping and analysis of illumina reads for transcriptome of medicago truncatula during the early organogenesis of the nodule. Bioprotocol3, e966. 10.21769/BioProtoc.966BottiniR.CassánF.PiccoliP. (2004). Gibberellin production by bacteria and its involvement in plant growth promotion and yield increase. Appl. Microbiol. Biotechnol.65, 497–503. 10.1007/s00253-004-1696-115378292ChowdappaS.JagannathS.KonappaN.UdayashankarA. C.JogaiahS. (2020). Detection and characterization of antibacterial siderophores secreted by endophytic fungi from Cymbidium aloifolium. Biomolecules10, 1412. 10.3390/biom1010141233036284ChowdharyK.SharmaS. (2020). Plant growth promotion and biocontrol potential of fungal endophytes in the inflorescence of Aloe vera L. Proc. Natl. Acad. Sci. India Section B Biol. Sci.90, 1045–1055. 10.1007/s40011-020-01173-3ConesaA.GötzS.García-GómezJ. M.TerolJ.TalónM.RoblesM. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics21, 3674–3676. 10.1093/bioinformatics/bti61016081474FiedlerH.-P.KrastelP.MüllerJ.GebhardtK.ZeeckA. (2001). Enterobactin: the characteristic catecholate siderophore of Enterobacteriaceae is produced by Streptomyces species. FEMS Microbiol. Lett.196, 147–151. 10.1111/j.1574-6968.2001.tb10556.x11267771FordeB. G.RobertsM. R. (2014). Glutamate receptor-like channels in plants: a role as amino acid sensors in plant defence?F1000Prime Rep.6, 37–37. 10.12703/P6-3724991414GentlemanR. C.CareyV. J.BatesD. M.BolstadB.DettlingM.DudoitS.. (2004). Bioconductor: open software development for computational biology and bioinformatics. Genome Biol.5, R80. 10.1186/gb-2004-5-10-r8015461798GlickmannE.DessauxY. (1995). A critical examination of the specificity of the salkowski reagent for indolic compounds produced by phytopathogenic bacteria. Appl. Environ. Microbiol.61, 793–796. 10.1128/aem.61.2.793-796.199516534942GopalakrishnanS.SrinivasV.AlekhyaG.PrakashB.KudapaH.RathoreA.. (2015). The extent of grain yield and plant growth enhancement by plant growth-promoting broad-spectrum Streptomyces sp. in chickpea. SpringerPlus4, 1–10. 10.1186/s40064-015-0811-325646153GopalakrishnanS.SrinivasV.Sree VidyaM.RathoreA. (2013). Plant growth-promoting activities of Streptomyces spp. in sorghum and rice. SpringerPlus2, 1–8. 10.1186/2193-1801-2-57424255867GoswamiD.ThakkerJ. N.DhandhukiaP. C. (2016). Portraying mechanics of plant growth promoting rhizobacteria (PGPR): a review. Cogent Food Agric.2, 1127500. 10.1080/23311932.2015.1127500GuY.-Q.MoM.-H.ZhouJ.-P.ZouC.-S.ZhangK.-Q. (2007). Evaluation and identification of potential organic nematicidal volatiles from soil bacteria. Soil Biol. Biochem.39, 2567–2575. 10.1016/j.soilbio.2007.05.011HanX.KumarD.ChenH.WuS.KimJ.-Y. (2014). Transcription factor-mediated cell-to-cell signalling in plants. J. Exp. Bot.65, 1737–1749. 10.1093/jxb/ert42224347464HayatR.AliS.AmaraU.KhalidR.AhmedI. (2010). Soil beneficial bacteria and their role in plant growth promotion: a review. Ann. Microbiol.60, 579–598. 10.1007/s13213-010-0117-1HowladerJ.RobinA. H. K.NatarajanS.BiswasM. K.SumiK. R.SongC. Y.. (2020). Transcriptome analysis by rna–seq reveals genes related to plant height in two sets of parent-hybrid combinations in easter lily (Lilium longiflorum). Sci. Rep.10, 1–15. 10.1038/s41598-020-65909-x32494055HuD.LiS.LiY.PengJ.WeiX.MaJ.. (2020). Streptomyces sp. strain TOR3209: a rhizosphere bacterium promoting growth of tomato by affecting the rhizosphere microbial community. Sci. Rep.10, 20132. 10.1038/s41598-020-76887-533208762HuD.LiX.ChangY.HeH.ZhangC.JiaN.. (2012). Genome sequence of Streptomyces sp. strain TOR3209, a rhizosphere microecology regulator isolated from tomato rhizosphere. J. Bacteriol.194, 1627–1627. 10.1128/JB.06684-1122374957JenningsW.ShibamotoT. (1980). “Part A - analytical considerations,” in Qualitative Analysis of Flavor and Fragrance Volatiles by Glass Capillary Gas Chromatography, eds W. Jennings and T. Shibamoto (New York, NY: Academic Press), 1–27.JiwenW.JingZ.BaoenX.YingyingL.GuanjieL. (2011). Isolation and screening of organic phosphorus-solubilizing bacterium and its initial identification. Henan Sci. 29, 31–34. 10.13537/j.issn.1004-3918.2011.01.016KanchiswamyC. N.MalnoyM.MaffeiM. E. (2015). Bioprospecting bacterial and fungal volatiles for sustainable agriculture. Trends Plant Sci.20, 206–211. 10.1016/j.tplants.2015.01.00425659880KanehisaM.GotoS. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res.28, 27–30. 10.1093/nar/28.1.2710592173KangY.JangS.-W.LeeH. J.BarchengerD. W.JangS. (2021). Expression profiling of heat shock protein genes as putative early heat-responsive members in lettuce. Horticulturae7, 312. 10.3390/horticulturae7090312KaurA.. (2020). A Novel Maize Dwarf Resulting From a Gain-of-Function Mutation in a Glutamate Receptor Gene. West Lafayette: Purdue University Graduate School.KhanA. L.WaqasM.KangS.-M.Al-HarrasiA.HussainJ.Al-RawahiA.. (2014). Bacterial endophyte Sphingomonas sp. LK11 produces gibberellins and IAA and promotes tomato plant growth. J. Microbiol.52, 689-695. 10.1007/s12275-014-4002-724994010KorpiA.JärnbergJ.PasanenA.-L. (2009). Microbial volatile organic compounds. Crit. Rev. Toxicol.39, 139–193. 10.1080/1040844080229149719204852LailaR.RobinA. H. K.YangK.ParkJ.-I.SuhM. C.KimJ.. (2017). Developmental and genotypic variation in leaf wax content and composition, and in expression of wax biosynthetic genes in Brassica oleracea var. capitata. Front. Plant Sci.7, 1972. 10.3389/fpls.2016.0197228119701LeeS.KaJ.-O.SongH.-G. (2012). Growth promotion of Xanthium italicum by application of rhizobacterial isolates of Bacillus aryabhattai in microcosm soil. J. Microbiol.50, 45–49. 10.1007/s12275-012-1415-z22367936LiC.ZhangJ.ZhaoC.JiangH.RenX.SuW.. (2018). Effects of 3-methyl-1-butanol on seed germination and seedling growth of maize and wheat. Bull. Botanical Res.38, 785–789. 10.7525/j.issn.1673-5102.2018.05.019LiY.LiuJ.Díaz-CruzG.ChengZ.BignellD. R. (2019). Virulence mechanisms of plant-pathogenic Streptomyces species: an updated review. Microbiology165, 1025–1040. 10.1099/mic.0.00081831162023LiuP.-L.DuL.HuangY.GaoS.-M.YuM. (2017). Origin and diversification of leucine-rich repeat receptor-like protein kinase (LRR-RLK) genes in plants. BMC Evol. Biol.17, 47. 10.1186/s12862-017-0891-528173747LugtenbergB.KamilovaF. (2009). Plant-growth-promoting rhizobacteria. Annu. Rev. Microbiol.63, 541–556. 10.1146/annurev.micro.62.081307.16291819575558MartínezO. A.JorqueraM. A.CrowleyD. E.de la Luz MoraM. (2011). Influence of nitrogen fertilisation on pasture culturable rhizobacteria occurrence and the role of environmental factors on their potential PGPR activities. Biol. Fertil. Soils47, 875–885. 10.1007/s00374-011-0593-xMishraR. C.GroverA. (2016). ClpB/Hsp100 proteins and heat stress tolerance in plants. Crit. Rev. Biotechnol.36, 862–874. 10.3109/07388551.2015.105194226121931MohantyP.SinghP. K.ChakrabortyD.MishraS.PattnaikR. (2021). Insight into the role of PGPR in sustainable agriculture and environment. Front. Sustain. Food Syst5, 667150. 10.3389/fsufs.2021.667150MyoE. M.GeB.MaJ.CuiH.LiuB.ShiL.. (2019). Indole-3-acetic acid production by Streptomyces fradiae NKZ-259 and its formulation to enhance plant growth. BMC Microbiol.19, 155. 10.1186/s12866-019-1528-131286877NassarA. H.El-TarabilyK. A.SivasithamparamK. (2005). Promotion of plant growth by an auxin-producing isolate of the yeast Williopsis saturnus endophytic in maize (Zea mays L.) roots. Biol. Fertil. Soils42, 97–108. 10.1007/s00374-005-0008-yNautiyalC. S.. (1999). An efficient microbiological growth medium for screening phosphate solubilizing microorganisms. FEMS Microbiol. Lett.170, 265–270. 10.1111/j.1574-6968.1999.tb13383.x9919677NebertD. W.. (1991). MINIREVIEW: proposed role of drug-metabolizing enzymes: regulation of steady state levels of the ligands that effect growth, homeostasis, differentiation, and neuroendocrine functions. Mol. Endocrinol.5, 1203–1214. 10.1210/mend-5-9-12031663211NeilandsJ.. (1987). “Comparative biochemistry of microbial iron assimilation,” in Iron Transport in Microbes, Plants and Animals (Weinheim: Federal Republic of Germany), 3–34.OlanrewajuO. S.BabalolaO. O. (2019). Streptomyces: implications and interactions in plant growth promotion. Appl. Microbiol. Biotechnol.103, 1179–1188. 10.1007/s00253-018-09577-y30594952OleńskaE.MałekW.WójcikM.SwiecickaI.ThijsS.VangronsveldJ. (2020). Beneficial features of plant growth-promoting rhizobacteria for improving plant growth and health in challenging conditions: a methodical review. Sci. Total Environ.743, 140682. 10.1016/j.scitotenv.2020.14068232758827OlsenA. N.ErnstH. A.LeggioL. L.SkriverK. (2005). NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci.10, 79–87. 10.1016/j.tplants.2004.12.01015708345OtienoN.LallyR.KiwanukaS.LloydA.RyanD.GermaineK.. (2015). Plant growth promotion induced by phosphate solubilizing endophytic Pseudomonas isolates. Front. Microbiol.6, 745. 10.3389/fmicb.2015.0074526257721PailanS.GuptaD.ApteS.KrishnamurthiS.SahaP. (2015). Degradation of organophosphate insecticide by a novel Bacillus aryabhattai strain SanPS1, isolated from soil of agricultural field in Burdwan, West Bengal, India. Int. Biodeterior. Biodegradation103, 191–195. 10.1016/j.ibiod.2015.05.006ParewaH. P.JoshiN.MeenaV. S.JoshiS.ChoudharyA.RamM.. (2021). “Chapter 9 - Role of biofertilizers and biopesticides in organic farming,” in Advances in Organic Farming, eds V. S. Meena, S. K. Meena, A. Rakshit, J. Stanley, and C. Srinivasarao (Woodhead Publishing), 133–159.Pérez-MirandaS.CabirolN.George-TéllezR.Zamudio-RiveraL. S.FernándezF. J. (2007). O-CAS, a fast and universal method for siderophore detection. J. Microbiol. Methods70, 127–131. 10.1016/j.mimet.2007.03.02317507108PhiQ.-T.ParkY.-M.SeulK.-J.RyuC.-M.ParkS.-H.KimJ.-G.. (2010). Assessment of root-associated Paenibacillus polymyxa groups on growth promotion and induced systemic resistance in pepper. J. Microbiol. Biotechnol.20, 1605–1613. 10.4014/jmb.1007.0701421193814QuinnG. A.BanatA. M.AbdelhameedA. M.BanatI. M. (2020). Streptomyces from traditional medicine: sources of new innovations in antibiotic discovery. J. Med. Microbiol.69, 1040. 10.1099/jmm.0.00123232692643RazinS. V.BorunovaV. V.MaksimenkoO. G.KantidzeO. L. (2012). Cys2His2 zinc finger protein family: classification, functions, and major members. Biochemistry77, 217–226. 10.1134/S000629791203001722803940RobinA.MougelC.SiblotS.VansuytG.MazurierS.LemanceauP. (2006). Effect of ferritin overexpression in tobacco on the structure of bacterial and pseudomonad communities associated with the roots. FEMS Microbiol. Ecol.58, 492–502. 10.1111/j.1574-6941.2006.00174.x17117991Romano-ArmadaN.Yañez-YazlleM. F.IrazustaV. P.RajalV. B.MoragaN. B. (2020). Potential of bioremediation and PGP traits in Streptomyces as strategies for bio-reclamation of salt-affected soils for agriculture. Pathogens9, 117. 10.3390/pathogens902011732069867SaharanB.NehraV. (2011). Plant growth promoting rhizobacteria: a critical review. Life Sci. Med. Res. 21, 30. Available online at: https://www.researchgate.net/publication/284340739SchmidtR.CordovezV.de BoerW.RaaijmakersJ.GarbevaP. (2015). Volatile affairs in microbial interactions. ISME J.9, 2329–2335. 10.1038/ismej.2015.4226023873SchöllerC. E. G.GürtlerH.PedersenR.MolinS.WilkinsK. (2002). Volatile metabolites from Actinomycetes. J. Agric. Food Chem.50, 2615–2621. 10.1021/jf011675411958631SchwachtjeJ.KarojetS.KunzS.BrouwerS.van DongenJ. T. (2012). Plant-growth promoting effect of newly isolated rhizobacteria varies between two Arabidopsis ecotypes. Plant Signal. Behav.7, 623–627. 10.4161/psb.2017622580689SgroyV.CassánF.MasciarelliO.Del PapaM. F.LagaresA.LunaV. (2009). Isolation and characterization of endophytic plant growth-promoting (PGPB) or stress homeostasis-regulating (PSHB) bacteria associated to the halophyte Prosopis strombulifera. Appl. Microbiol. Biotechnol.85, 371–381. 10.1007/s00253-009-2116-319655138ShamimuzzamanM.VodkinL. (2013). Genome-wide identification of binding sites for NAC and YABBY transcription factors and co-regulated genes during soybean seedling development by ChIP-Seq and RNA-Seq. BMC Genomics14, 477. 10.1186/1471-2164-14-47723865409SpaepenS.VanderleydenJ. (2011). Auxin and plant-microbe interactions. Cold Spring Harb. Perspect. Biol.3, a001438. 10.1101/cshperspect.a00143821084388SpaepenS.VanderleydenJ.RemansR. (2007). Indole-3-acetic acid in microbial and microorganism-plant signaling. FEMS Microbiol. Rev.31, 425–448. 10.1111/j.1574-6976.2007.00072.x17509086Suárez-MorenoZ. R.Vinchira-VillarragaD. M.Vergara-MoralesD. I.CastellanosL.RamosF. A.GuarnacciaC.. (2019). Plant-growth promotion and biocontrol properties of three Streptomyces spp. isolates to control bacterial rice pathogens. Front. Microbiol.10, 290. 10.3389/fmicb.2019.0029030858835TerraL.RatcliffeN.CastroH. C.VicenteA. C. P.DysonP. (2021). Biotechnological potential of Streptomyces siderophores as new antibiotics. Curr. Med. Chem.28, 1407–1421. 10.2174/092986732766620051023551232389112Toledo-OrtizG.HuqE.QuailP. H. (2003). The arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell15, 1749–1770. 10.1105/tpc.01383916896230van RuissenF.RuijterJ. M.SchaafG. J.AsgharnegadL.ZwijnenburgD. A.KoolM.. (2005). Evaluation of the similarity of gene expression data estimated with SAGE and Affymetrix GeneChips. BMC Genomics6, 91. 10.1186/1471-2164-6-9115955238VasconcellosR. L. F. dSilvaM. C. P. dRibeiroC. M.CardosoE. J. B. N. (2010). Isolation and screening for plant growth-promoting (PGP) actinobacteria from Araucaria angustifolia rhizosphere soil. Sci. Agric.67, 743–746. 10.1590/S0103-90162010000600019VejanP.AbdullahR.KhadiranT.IsmailS.Nasrulhaq BoyceA. (2016). Role of plant growth promoting rhizobacteria in agricultural sustainability—a review. Molecules21, 573. 10.3390/molecules2105057327136521WooH.-H.OrbachM. J.HirschA. M.HawesM. C. (1999). Meristem-localized inducible expression of a UDP-glycosyltransferase gene is essential for growth and development in pea and alfalfa. Plant Cell11, 2303–2315. 10.1105/tpc.11.12.230310590160ZhangT.ZhaoY.WangY.LiuZ.GaoC. (2018). Comprehensive analysis of MYB gene family and their expressions under abiotic stresses and hormone treatments in Tamarix hispida. Front. Plant Sci.9, 1303. 10.3389/fpls.2018.0130330283465ZhaoH.LiX.MaL. (2012). Basic helix-loop-helix transcription factors and epidermal cell fate determination in Arabidopsis. Plant Signal. Behav.7, 1556–1560. 10.4161/psb.2240423073001