
==== Front
BMC Microbiol
BMC Microbiol
BMC Microbiology
1471-2180
BioMed Central London

3479
10.1186/s12866-024-03479-y
Research
Mechanisms of ROS-mediated interactions between Bacillus aryabhattai LAD and maize roots to promote plant growth
Deng Chao 12
Zeng Nan 1
Li Chunji 3456
Pang Jiahe 7
Zhang Ning zhangning@syau.edu.cn

7
Li Bingxue libingxue@syau.edu.cn

1
1 https://ror.org/01n7x9n08 grid.412557.0 0000 0000 9886 8131 College of Land and Environment, Shenyang Agricultural University, Shenyang, 110866 People’s Republic of China
2 https://ror.org/04s5wwh39 grid.469526.a 0000 0001 2065 4085 School of Chemistry and Life Science, Anshan Normal University, Anshan, 114007 People’s Republic of China
3 https://ror.org/000b7ms85 grid.449900.0 0000 0004 1790 4030 College of Agriculture and Biology, Zhongkai University of Agriculture and Engineering, Guangzhou, 510225 People’s Republic of China
4 https://ror.org/000b7ms85 grid.449900.0 0000 0004 1790 4030 Innovative Institute for Plant Health, Zhongkai University of Agriculture and Engineering, Guangzhou, 510225 People’s Republic of China
5 https://ror.org/05ckt8b96 grid.418524.e 0000 0004 0369 6250 Key Laboratory of Green Prevention and Control on Fruits and Vegetables in South China, Ministry of Agriculture and Rural Affairs, Guangzhou, 510225 People’s Republic of China
6 Guangdong University Key Laboratory for Sustainable Control of Fruit and Vegetable Diseases and Pests, Guangzhou, 510225 China
7 https://ror.org/01n7x9n08 grid.412557.0 0000 0000 9886 8131 College of Bioscience and Biotechnology, Shenyang Agricultural University, Shenyang, 110866 People’s Republic of China
6 9 2024
6 9 2024
2024
24 32725 3 2024
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Background

Plant growth-promoting rhizobacteria (PGPR), as a group of environmentally friendly bacteria growing in the rhizosphere of plants, play an important role in plant growth and development and resistance to environmental stresses. However, their limited understanding has led to the fact that their large-scale use in agriculture is still scarce, and the mechanisms by which beneficial bacteria are selected by plants and how they interact with them are still unclear.

Method

In this study, we investigated the interaction between the auxin-producing strain Bacillus aryabhattai LAD and maize roots, and performed transcriptomic and metabolomic analyses of Bacillus aryabhattai LAD after treatment with maize root secretions(RS).

Results

Our results show that there is a feedback effect between the plant immune system and bacterial auxin. Bacteria activate the immune response of plant roots to produce reactive oxygen species(ROS), which in turn stimulates bacteria to synthesize IAA, and the synthesized IAA further promotes plant growth. Under the condition of co-culture with LAD, the main root length, seedling length, root surface area and root volume of maize increased by 197%, 107%, 89% and 75%, respectively. In addition, the results of transcriptome metabolome analysis showed that LAD was significantly enriched in amino acid metabolism, carbohydrate metabolism and lipid metabolism pathways after RS treatment, including 93 differentially expressed genes and 45 differentially accumulated metabolites.

Conclusion

Our findings not only provide a relevant model for exploring the effects of plant-soil microbial interactions on plant defense functions and thereby promoting plant growth, but also lay a solid foundation for the future large-scale use of PGPR in agriculture for sustainable agricultural development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-024-03479-y.

Keywords

Bacillus aryabhattai
PGPR
Maize
IAA
ROS
Feedback effect
National Natural Science Foundation of China31271818 International Cooperation Project of Universities in Liaoning Province2023-01 Liaoning Province Rural Science and Technology Special Action Project2022030673-JH5/104 Shenyang Science and Technology Project22-319-2-13 China Scholarship CouncilCSC202208850002 issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcBackground

Plants have been in a symbiotic relationship with soil microorganisms during their growth and development. As a key member of the soil ecosystem, plant growth-promoting rhizobacteria (PGPR) is mainly distributed in the roots and rhizosphere of plants and participates in a variety of biological activities in the soil ecosystem. PGPR can improve crop productivity and resistance to pathogens by inducing complex changes in plant growth and development [1, 2]. In addition, PGPR also competitively settles in plant roots and plays a crucial role in promoting plant growth through different mechanisms such as phosphate solubilization, nitrogen fixation, IAA production, siderophores, biofilm synthesis, VOCs production, and induction of plant systemic resistance [1–3]. To date, many PGPR strains have been isolated for extensive and in-depth research, which is of great significance for promoting sustainable agricultural development [4, 5]. However, so far, the mechanism of plant-rhizosphere microbiome interaction is extremely complex, and little is known about how beneficial bacteria interact with the plant immune system.

In general, the growth and stress effects of plants under the action of PGPR could be determined by changes in hormone levels and related signaling pathways. Each plant has its specific microflora, and the structure of the microflora in the rhizosphere of plants is composed of root secretions produced by plant roots and signal molecules produced by plants and microorganisms [6–8]. With the help of these signaling molecules, a variety of close relationships have been established between plant roots and microorganisms. Acturally, plant roots can interact with a large number of bacteria in the surrounding soil. When bacteria migrate from bulk soil to the root surface and colonize further, the bacterial diversity decreases, which seems to indicate that plants exert selective forces on their colonizing bacteria [9]. An early filter used by plants to identify and respond to bacteria is their immune system [10], which recognizes bacterial flagella, peptidoglycans and bacterial elongation factors through MAMPs [11]. When plant roots recognize these molecules, a series of molecular events will occur. The earliest stage is the outflow of calcium ions and the outbreak of reactive oxygen species (ROS) [12, 13]. Under the condition of pathogen infection, the recognition and activation of plant immune system has been widely described, and MAMP receptors can recognize molecules in all bacteria, and it is found that beneficial bacteria can also induce immune responses similar to pathogenic bacteria [14, 15]. Therefore, the effect of plant immune system on healthy root microbiota and how beneficial microorganisms respond to plant immune system is still an active research field, and there are still many issues that need to be explored in depth [16].

Previously, we found that B. aryabhattai LAD, as a good PGPR strain, significantly promotes maize growth [17, 18]. Therefore, in order to understand the interaction between B. aryabhattai LAD and maize roots to further understand the interaction between beneficial bacteria and plant immune system, we co-cultured the strain LAD with maize to explore its IAA synthesis level, intracellular ROS level and the relationship between ROS and IAA synthesis. Subsequently, we further analyzed the changes of intracellular transcriptome and metabolome of B. aryabhattai LAD under the stimulation of maize root secretions. These results provide a relevant model for exploring the effects of plant-soil microbial interactions on plant defense functions and thus promoting plant growth, and also lay a solid foundation for the future application of PGPR to field production for sustainable agricultural development.

Materials and methods

Strain and culture conditions

The strain B. aryabhattai LAD used in this study was isolated from the rhizosphere soil of maize and preserved in the China General Microbiological Culture Collection Center ( CGMCC18653 ), the NCBI accession number is PRJNA716506. B. aryabhattai LAD was inoculated in beef extract peptone medium and cultured to logarithmic phase at 37 °C and 180 rpm. Subsequently, the LAD bacterial suspension was transferred to 50 mL beef extract peptone liquid medium at 1% inoculation amount, and cultured at 37 °C and 180 rpm for 24 h as seed solution.

The plant material used in this study was maize and the seeds were sourced from Dalian Baisite Seed Co., Ltd. The maize seeds were surface sterilized by soaking in 70% ethanol for 2 min and then in 1% sodium hypochlorite solution for 20 min. After soaking, the seeds were rinsed three times in sterile desalinated water and prepared for use. Plants were kept in climate cabinets at 21 °C (180 µmol light m− 2·s− 1 at plant level; 16 h:8 h, light: dark).

Effects of exogenous IAA and B. aryabhattai LAD on the growth of maize seedlings

The B. aryabhattai LAD was inoculated into beef extract peptone medium at 1% inoculation amount, and cultured at 37 °C and 180 rpm for 72 h until OD600 nm was 1.4–1.5. The fermentation broth was taken and the IAA content was detected by Angilent ultra-high performance liquid chromatograph 1290 Infinity II ( Agilent Technologies Inc., California, USA ). The chromatographic column was Agilent Polaris C18-A ( 250 mm × 4.6 mm ), the flow rate was 1 mL / min, the injection volume was 10 µL, the column temperature was 35 °C, the wavelength was 254 nm, the mobile phase A was pure methanol, and the mobile phase B was glacial acetic acid aqueous solution with pH value of 3.2 [ V ( A ) : V ( B ) = 5.5 : 4.5 ].

In addition, the bacterial suspension cultured for 72 h was diluted in a gradient of 10 times, and the B. aryabhattai LAD bacterial suspension diluted to 102 CFU mL− 1 was used to soak the maize seeds for 24 h. The seeds were placed in a hydroponic tank containing B. aryabhattai LAD bacterial suspension and IAA for 14 days. The normal culture medium was used as a blank control to measure the roots of maize seedlings. The root development was measured by WinRHIZO Reg STD4800 Root Analysis System (Beijing Ruiding Environmental Technology Co., Ltd., Beijing, China). The growth promotion indexes of maize seedlings and their roots were counted using percentage values calculated as [(treatment group - control group)/control group]*100%.

Collection of maize root secretions

To test the effect of maize root secretions on IAA synthesis in LAD, root secretions were collected by hydroponics. The seeds soaked with LAD were cultured in a hydroponic tank containing B. aryabhattai LAD bacterial suspension, and the normal culture medium was used as a blank control. After 14 days of culture, the root secretions were collected and filtered with 0.22 μm filter membrane and stored at 4 °C. At the same time, the IAA content was detected by Angilent ultra-high performance liquid chromatograph 1290 Infinity II ( Agilent Technologies Inc., California, USA ).

Effects of maize root secretions and H2O2 onB. aryabhattai LAD

The B. aryabhattai LAD was inoculated into beef extract peptone medium supplemented with 0.1%, 0.2% and 0.4% H2O2 and 1%, 2% and 4% maize root secretions at a 1% inoculation amount, respectively. The blank beef extract peptone medium inoculated with B. aryabhattai LAD bacterial liquid was used as a blank control. One part was placed in a 12-well plate at 37 °C, 180 rpm for 24 h in a microplate reader, and the OD ( 600 nm ) was measured every 1 h. The other part was in a 50 mL conical flask. IAA production was determined by high performance liquid chromatography after 72 h of shaking culture at 180 rpm.

Determination of ROS levels in cells and maize roots

Maize roots were collected at 2,4,6,8,10 and 12 h, respectively. At the same time, bacterial cells were collected by centrifugation at 8000 rpm for 10 min, and then the bacterial cells were fully washed with PBS ( pH 7.0 ) and suspended in 1 mL phosphate buffer.

The relative levels of intracellular reactive oxygen species (ROS) in B. aryabhattai LAD were determined by fluorescence assay, mainly using the oxidant-sensitive probe 2′,7′-dichlorodihydrofluorescein diacetate ( DCFH-DA ) (Beijing Solarbio Science & Technology Co., Ltd, Beijing, China). The fluorescence intensity was detected using a Bio Tek Cytation5 microplate reader (Beijing Derica Biotechnology Co., Ltd, Beijing, China), and the emission wavelength was set at 525 nm and the excitation wavelength was set at 488 nm. The relative level of ROS in maize roots was measured by Plant Reactive Oxygen Species ELISA Kit instruction at 450 nm to reflect the ROS level. Each sample was performed using three independent biological replicates.

Total RNA extraction and transcriptome analysis

B. aryabhattai LAD treated with maize root secretions was selected as the treatment group (G), and B. aryabhattai LAD cultured under normal conditions was used as the control group (CK). 1 mL bacterial suspension was centrifuged at 10,000 r min− 1 for 2 min at 4 °C in a 1.5 mL enzyme-free tube, and the supernatant was discarded to collect the bacteria and quickly frozen in liquid nitrogen for 15 min. The total RNA was extracted using the Trizol kit ( takara, dalian, China ) according to the kit instructions.

High-quality RNA samples were collected by NEB Next Ultra RNA Library Preparation Kit ( NEB, CA ) for the preparation of cDNA library construction. Next-Generation Sequencing ( NGS ) was used to perform transcriptome sequencing analysis and comparison of samples based on the Illumina HiSeq sequencing platform. GO and KEGG enrichment analysis was performed using Blast2GO ( v2.5 ) and KOBAS ( v3.0 ) software. Differentially expressed genes ( DEGs ) were performed using the DESeq ( v1.32.0 ) program, and log2 ( fold change ) ≥ 1 and p ( padj ) ≤ 0.05 were considered to be DEGs. Pathway analysis of DEGs was performed using the KEGG PATHWAY database and Blast2GO.

Metabolite extraction and metabolome analysis

2 mL of pre-cooled 60% methanol water was placed in a 5 mL centrifuge tube, and then 2 mL of bacterial solution was quickly added to it. The mixture was mixed manually for 5 s. The quenched bacterial solution was centrifuged at 3,000 r min− 1, 4 °C for 10 min, and the supernatant was removed and placed in dry ice. The metabolites were detected by Thermo Vanquish ( Thermo Fisher Scientific, USA ) ultra-high performance liquid chromatography system combined with Thermo Orbitrap Exploris 120 mass spectrometry detector ( Thermo Fisher Scientific, USA ).

The identification of relevant metabolites in the metabolome was based on the databases of HMDB (http://www.hmdb.ca), Metlin (http://metlin.scripps.edu), MassBank (http://www.massbank.jp), LipidMaps (http://www.lipidmaps.org) and mzClound (https://www.mzcloud.org) databases were performed. Differently accumulated metabolites (DAMs) were assessed by partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA) models using variable importance in prediction (VIP) scores. Thresholds for the analytical models were set at VIP ≥ 1, fold change ≥ 1.5 or ≤ 0.667, and p ≤ 0.05. Integration analyses of metabolomics and transcriptomics were performed via Pearson correlation coefficients, as previously described [19, 20].

Real-time quantitative polymerase chain reaction (qPCR) analysis

To further elucidate the ROS-mediated interactions between B. aryabhattai LAD and maize rhizosphere to promote growth, we first examined the changes in the expression of genes related to growth hormone synthesis in B. aryabhattai LAD after exogenous additions of root secretions and H2O2. Subsequently, changes in the expression of growth hormone synthesis and antioxidant-related genes in the maize root system after exogenous addition of B. aryabhattai LAD treatment were further examined to determine the relationship between activation of the plant immune system and growth hormone synthesis following B. aryabhattai LAD interactions with the maize root system.

The SYBR Green Premix Pro Taq HS qPCR Kit (AG11718, ACCURATE BIOTECHNOLOGY, HUNAN, Co., Ltd, Changsha, China) was used to perform qPCR in a Quant Studio6 Flex system (Thermo Fisher Scientific, USA). 16 S rDNA and Actin were selected as the internal reference genes for B. aryabhattai LAD and maize roots, respectively. To analyze the qPCR results, the relative expression of each gene was calculated by the comparative crossover (CP) method and expressed as 2−∆∆CT. Each gene expression analysis was performed using three independent biological replicates. Specific primers used for qPCR are listed in Table S1.

Statistical analysis

IBM SPSS 25.0 (SPSS Inc., Chicago, IL, USA) software was used to statistically analyze the experimental data. One-way ANOVA ( One-way ANOVA ) was used to determine significant differences between the following: different treatments on maize root growth, different treatments on LAD growth and its IAA yield, and different treatments on the level of ROS in LAD and in the maize root system. All trials were conducted in triplicate and data were averaged and expressed as mean ± standard error (SE).

The qPCR data of selected gene expression profiles were analyzed using the software package IBM SPSS 25.0 of Windows (SPSS Inc., Chicago, IL, USA). Prior to analysis, data were normalized by transforming as ln (x + 1). All data were expressed as means ± standard error of at least triplicates of independent cultures. The means of different treatments were compared by one-way ANOVA using Turkey’s honestly signifcant difference test at the 5% probability level (p < 0.05). Metabolomics analysis was evaluated using variable importance in the projection (VIP) in the frst principal component of each OPLS-DA model. Thresholds with P < 0.05 and VIP > 1 were considered statistically significant. MetaboAnalyst software (version 3.0) was used for pathway enrichment analysis [21]. Differentially expressed genes were analyzed for differences using DESeq, and the threshold for identifying DEGs was set at |log2 (fold change)|≥1,padj ≤ 0.05.

Results

Effect of exogenous addition of B. aryabhattai LAD and IAA on the growth of maize roots

The results of root development of maize seedlings under different treatment conditions measured by root analyzer are shown in Fig. 1 and Table S2. The main root length of maize seedlings increased by 197% and 167% respectively after exogenous addition of B. aryabhattai LAD and IAA. The main root thickness increased by 22% and 16% respectively, and the seedling length of maize seedlings increased by 107% and 78% respectively. Exogenous addition of B. aryabhattai LAD and IAA increased root surface area by 89% and 61%, root volume by 75% and 34%, and total root length by 97% and 63%, respectively. Exogenous addition of B. aryabhattai LAD and direct exogenous addition of IAA with the same content could promote the growth of maize seedlings, but the growth of maize roots in the exogenous addition of B. aryabhattai LAD group was significantly higher than that in the IAA group with the same concentration.

Effect of root secretion and H2O2 addition on the growth and IAA production of B. aryabhattai LAD

The effects of root secretion and H2O2 addition on the growth of B. aryabhattai LAD and its IAA production are shown in Fig. 2A, B. The growth of B. aryabhattai LAD was strongly inhibited by the addition of H2O2 treatment, and its level of IAA production was significantly lower than that of the control group, and the levels of IAA synthesized after the addition of 0.1%, 0.2%, and 0.4% of H2O2 to the B. aryabhattai LAD culture for 72 h were 15.1, 17.5, and 15.8 µg/mL, respectively, while the growth inhibition of B. aryabhattai LAD was strongly alleviated after the addition of IAA to its culture. Subsequently, the growth of B. aryabhattai LAD was also inhibited by the addition of root secretion, but the level of IAA production was significantly higher than that of the control group and the other treatment groups, and the levels of IAA synthesized after 72 h of incubation of B. aryabhattai LAD with 1%, 2%, and 4% of root secretion were 26.7, 24.2, and 23.4 µg/mL, respectively.

Fig. 1 Effect of exogenous addition of B. aryabhattai LAD and IAA on root growth of maize seedlings. All data are from three independent biological replicate means ± SE (n = 3). Different letters (a-b) indicate the degree of significant difference between treatments. Values with the same letter above the bar graph indicate no significant difference according to the Tukey test

Fig. 2 Effects of root secretions and H2O2 addition on the growth (A) and IAA production (B) of B. aryabhattai LAD and changes in ROS levels in cells after interaction between B. aryabhattai LAD (C) and maize roots (D). All data are from three independent biological replicate means ± SE (n = 3). Different letters (a-b) indicate the degree of significant difference between treatments. Values with the same letter above the bar graph indicate no significant difference according to the Tukey test

Effect of different treatment conditions on ROS levels in B. aryabhattai LAD and maize roots

The ROS levels in the roots of maize at different times are shown in Fig. 2C, and the ROS levels were determined by taking the maize roots on the third day from the start of the submerged seed culture, and both the control group and the treatment group showed a trend of increasing and then decreasing, and the ROS levels of the maize roots were significantly higher than those of the control group after the addition of the B. aryabhattai LAD treatment.

The ROS levels of B. aryabhattai LAD strains at different culture stages were shown in Fig. 2D, which combined with the B. aryabhattai LAD growth curve showed that the intracellular ROS in the culture to 2 h, the concentration of the bacteria was small, and the difference of its intracellular ROS levels was not significant, after 2 h, the bacteria multiplied a lot and entered into the logarithmic growth period, and the ROS levels in the control group began to be significantly higher than that of the treatment group, and by 8 h, the growth of the strains began to slow down, and the ROS levels in the treatment group began to decrease and the control group was still at a high level.

Transcriptome analysis

To identify DEGs in B. aryabhattai LAD under maize root secretion treatment conditions, we constructed 6 cDNA libraries and filtered the sequencing data of transcriptome samples to obtain a total of 17.71 Gb of data. Gene expression was differentially analyzed using DESeq, and the threshold for identifying DEGs was set at |log2 (fold change)| ≥ 1, padj ≤ 0.05. Compared with the control group, there were 93 DEGs in B. aryabhattai LAD under root secretion treatment conditions, of which 57 genes were up-regulated and 36 genes were down-regulated .

Subsequently, we performed GO and KEGG enrichment analysis to further understand the biological functions of these DEGs and the related biological processes involved. As shown in Fig. 3A, in terms of GO enrichment analysis, the DEGs of B. aryabhattai LAD under root secretions treatment were mainly enriched in molecular functions (MF) and biological processes (BP). The top three GO terms significantly enriched in the MF category were phosphoribosylformylglycinamidine synthase activity (GO:0004642), carbon-nitrogen ligase activity, with glutamine as amido-N-donor (GO:0016884), ligase activity, forming carbon-nitrogen bonds (GO:0016879). In the BP category, the top three GOs significantly enriched were phospholipid catabolic process (GO:0009395), urea cycle (GO:000050), and urea metabolic process (GO:0019627). KEGG enrichment analysis showed (Fig. 3B) that the top 20 significantly enriched pathways of DEGs in B. aryabhattai LAD under root secretion treatment conditions were mainly associated with Metabolism, Genetic Information Processing, Environmental Information Processing in KEGG level 1 and Purine metabolism (ko00230), Inositol phosphate metabolism (ko00562), Alanine, aspartate and glutamate metabolism (ko00250), Glycerophospholipid metabolism (ko00564), ABC transporters (ko02010), Tyrosine metabolism (ko00350), Lysine degradation (ko00310), Oxidative phosphorylation (ko00190) Glycolysis / Gluconeogenesis (ko00010), Pyruvate metabolism (ko00620) in KEGG level 2. These results may be importantly related to the promotion of maize growth and resistance to environmental stresses by B. aryabhattai LAD.

Fig. 3 Transcriptome analysis of G and CK in B. aryabhattai LAD. (A) Circos plots depicting GO enrichment analysis of DEGs. (B) Circos plots illustrating KEGG enrichment analysis of DEGs. Starting from the outer circle and moving inward the representation includes the top 20 enriched GO terms or KEGG pathways (ring 1), the number of background genes in the genome (ring 2), the number of upregulated and downregulated genes (ring 3), and the rich factor of DEGs in corresponding GO terms or KEGG pathways (ring 4) are represented

Metabolome analysis

In order to evaluate the changes of intracellular metabolism of B. aryabhattai LAD after root secretions treatment, LC-MS / MS was used to evaluate the metabolic profile. The results showed that a total of 1,572 metabolites were identified based on the MS2 spectral database. The identified metabolites were classified according to chemical classification statistics, including 82 classifications, of which Carboxylic acids and derivatives were up to 361. Secondly, there are 113 Steroids and steroid derivatives, 105 Benzene and substituted derivatives, 105 Fatty Acyls, 105 Prenol lipids, 94 Organooxygen compounds, 46 Indoles and derivatives, 45 Glycerophospholipids and 30 Organonitrogens ( Fig. 4A ). At a threshold setting of VIP ≥ 1, fold change ≥ 1.5 or ≤ 0.667, and p ≤ 0.05, a total of 1,574 metabolites were detected in the treatment and control groups, of which 967 were detected in the positive ion mode and 607 were detected in the negative ion mode. Differential analysis of the detected metabolites showed a total of 45 differential metabolites, of which 8 differential metabolites were up-regulated and 37 differential metabolites were down-regulated (Fig. 4B). Subsequently, we performed KEGG enrichment analysis to better understand the biological functions of these differential metabolites and their involvement in relevant biological processes. The results showed that a total of 87 KEGG pathways were enriched by differential metabolites under root secretion treatment conditions. As shown in the figure, the first 20 significantly enriched pathways were mainly distributed in Metabolism, Environmental Information Processing and Organismal Systems at KEGG level 1, and Amino acid metabolism, Carbohydrate metabolism and Lipid metabolism at KEGG level 2 (Fig. 4C). Therefore, these metabolic pathways are important for B. aryabhattai LAD to respond to root secretions to resist environmental stresses and promote plant growth.

Fig. 4 Metabolome analysis of G and CK in B. aryabhattai LAD. (A) Quantitative weight of the identified metabolites in each chemical classification. (B) Heat map of hierarchical clustering of significantly different metabolites. (C) Circos plots illustrating KEGG enrichment analysis of DAMs. Starting from the outer circle and moving inward the representation includes the top 20 enriched GO terms or KEGG pathways (ring 1), the number of background genes in the genome (ring 2), the number of upregulated and downregulated genes (ring 3), and the rich factor of DEGs in corresponding GO terms or KEGG pathways (ring 4) are represented

Furthermore, the differentially expressed metabolites and all transcripts were simultaneously mapped to the KEGG pathway database, and there were 69 relevant metabolic pathways, including 213 genes. The main metabolic pathways included Tryptophan biosynthesis, Tryptophan metabolism, Zeatin synthesis, ABC transport system and Two-component system. Among them, 9 differential genes and 14 differential metabolites were annotated in the Tryptophan Metabolism and Tryptophan Biosynthesis pathways, which are closely related to IAA synthesis. Among them, K03781 (catalase: katE, CAT, catB, srpA), KO00453 (tryptophan 2,3 dioxygenase: TOD2, kynA), K00817 (histidyl-phosphate aminotransferase: hisC), K00891 (manganic acid kinase: aroK, aroL), K03786 (3-dehydroquinate dehydratase II: aroQ, qutE), KO01426 (amidase: amiE), and KO00128 (aldehyde dehydrogenase: ALDH) genes were significantly up-regulated after root secretion treatment; K16901 (monooxygenase), K01696 (tryptophan synthase β-chain: trpB), and K01817 (phosphoribosyl o-aminobenzoate isoform: trpF) were significantly down-regulated after root secretion treatment, and the results of integration of differential genes and differential metabolites related to IAA synthesis are shown in Fig. 5.

Fig. 5 Analysis of IAA synthesis pathway in B. aryabhattai LAD after root secretions treatment

Changes in the expression of genes related to IAA synthesis in B. aryabhattai LAD and maize roots

The expression changes of IAA synthesis related genes YUC, GH3 and PIN genes in maize roots under the condition of interaction with B. aryabhattai LAD are shown in the Fig. 6A. Among them, YUC gene was up-regulated on the 6d and 9d, and down-regulated on the 12d. The GH3 gene was down-regulated on 6d and day 9d, and the GH3 gene began to be up-regulated on 12d. PIN gene showed an upward trend at 6,9 and 12 d. The changes in the expression of antioxidant genes POD, SDO and CAT genes in maize roots showed that the expression of antioxidant-related genes decreased gradually on the 6d, 9d and 12d after adding B. aryabhattai LAD to maize roots (Fig. 6B).

Fig. 6 Changes in the expression of IAA (A) and antioxidant genes (B) in maize roots and genes related to IAA synthesis (C) in B. aryabhattai LAD

In the B. aryabhattai LAD strain, both root secretion and hydrogen peroxide treatments up-regulated the amiE and ALDH genes (Fig. 6C). The amiE gene was significantly up-regulated by root secretion, and the strain was induced by hydrogen peroxide to synthesize more growth hormone by up-regulating the transcript level of ALDH gene.

Discussion

Our results indicate that there is a feedback loop between the plant immune system and bacterial auxin synthesis (Fig. 7). When bacteria colonize the roots of plants, they will trigger the immune response of plants and the production of ROS, while ROS in turn will induce bacteria to synthesize auxin to reduce the toxicity of ROS to cells. At the same time, it promotes the spread of bacteria on roots and the formation of colonies, induces the expression of immune receptors in plants, and further accelerates the feedback loop [9]. Moreover, auxin, as an important plant hormone, plays a vital role in plant development [22]. Many bacterial species, such as Agrobacterium tumefaciens and Pseudomonas syringae, can control plant growth by synthesizing and secreting auxin [23–25]. However, despite decades of research on the production of bacterial auxin and its effects on plants, little is known about the role of bacterial auxin in bacterial physiology and its interaction with plants.

Fig. 7 Regulatory mechanisms of ROS and IAA in the promotion of maize root development by B. aryabhattai LAD

In this study, we found that when B. aryabhattai LAD responded to maize root secretions, the root secretion produced a stress response to B. aryabhattai LAD resulting in a significant increase in intracellular ROS levels, which acted as a feedback to regulate B. aryabhattai LAD intracellularly to respond to this stress effect, such as the increase in the synthesis of IAA production, while amino acid metabolism, carbohydrate metabolism and lipid metabolism pathways were significantly enriched. Free amino acids, carbohydrates and lipids are important regulators of cellular responses to environmental stresses and promote plant growth [21, 26–28]. The joint analysis showed that tryptophan synthesis, tryptophan metabolism, glutamate metabolism, cysteine metabolism, glycine metabolism and glutathione metabolism were significantly enriched in amino acid metabolic pathways and metabolites such as glutathione were significantly accumulated. Tryptophan synthesis and tryptophan metabolism are the main metabolic pathways for IAA synthesis. When B. aryabhattai LAD feels environmental stress, cells act as a feedback effect to activate IAA synthesis-related pathways to synthesize IAA to help cells resist environmental stress and promote plant growth. In the B. aryabhattai LAD tryptophan metabolic pathway, the tryptophan content was significantly elevated and the amiE gene transcript level was significantly up-regulated. Under the action of amiE, indoleacetamide was converted to indoleacetic acid, which in turn was secreted extracellularly through the ABC transport system and was absorbed and utilized by plants. In addition, the enrichment of glutamate metabolism, cysteine metabolism, and glycine metabolism may be mainly used for the synthesis of glutathione, because glutathione, as a substance possessing an extremely strong antioxidant capacity to scavenge ROS, inhibit the formation of free radicals, and lipid peroxidation, helps the B. aryabhattai LAD to play an important role in regulating the intracellular ROS in the event of environmental stress [29, 30].

In carbohydrate metabolism, glycolysis and gluconeogenesis, fructose and mannose metabolism, galactose metabolism, starch and sucrose metabolism, glyoxylate and dicarboxylic acid metabolism, oxidative phosphorylation and other metabolic pathways were significantly enriched, among which glucose, fructose, UDP-D-G, NAD and so on were significantly accumulated. As energy substances of cells, on the one hand, they provide energy sources for their own growth to accelerate cell metabolism, on the other hand, they can help cells improve their ability to resist osmotic stress and help cells resist environmental stress [31–33]. Moreover, it is worth noting that glyoxalate and dicarboxylic acid metabolism were also significantly enriched in B. aryabhattai LAD, and the antioxidant enzyme CAT was significantly up-regulated. Previous studies have shown that the glyoxylate cycle plays an important role in stress defense and resistance to environmental stress, especially when some bacteria activate the glyoxylate cycle to help cells resist environmental stress and activate related antioxidant enzyme systems to scavenge the toxic effects of ROS on the cells when they are subjected to environmental stress [34, 35]. Furthermore, inositol is also significantly accumulated in B. aryabhattai LAD, which is involved in a variety of physiological and biochemical processes as a precursor of many important metabolites. Cells resist the stress response produced by plant roots to maintain their own growth by up-regulating inositol [36, 37]. In the lipid metabolism pathway, the metabolic pathways related to glycerol synthesis were also significantly enriched, including glycerolipid metabolism and glycerol phospholipid metabolism, and the accumulation of lipids such as glycerol was significantly increased. It can not only help B. aryabhattai LAD resist the stress effect of plant roots, but also provide substrates and energy for other metabolic pathways, and improve the antioxidant capacity of cells in an indirect way.

The cross-talk between auxin and ROS has been increasingly emphasized. On the one hand, stress treatments cause an increase in ROS in the plant, which in turn regulates the activity of the root zone by controlling the accumulation of auxin; on the other hand, auxin can also reduce the redox state in the cell by modulating the activity of scavenging enzymes in vivo, which can alter the growth and development of the plant. In addition, ROS produced by plant roots also activates the immune response of rhizosphere bacteria to synthesize more IAA to further promote plant growth [38, 39]. In this study, it was shown that exogenous addition of B. aryabhattai LAD could significantly promote the growth of maize seedlings. The verification of the expression changes of IAA synthesis-related genes in maize showed that the IAA synthesis-related genes YUC, GH3 and PIN in maize roots were significantly up-regulated. They are closely related to the synthesis of IAA in plants and play an important role in the growth and development of plants [40–42], in which the YUC gene was up-regulated on the 6th and 9th days, indicating that IAA was synthesized in large quantities in maize roots at this time, while on the 12 th day, the YUC gene was down-regulated. It may be that the accumulation of IAA in the plant itself reaches a certain level, and the synthesis rate begins to decline. The GH3 gene was down-regulated on the 6th and 9th days, indicating that the IAA metabolism was relatively slow at the accumulation stage. At 12 days, the GH3 gene began to be up-regulated, and the IAA effect of plant roots may be relatively obvious. After exogenous addition of B. aryabhattai LAD, the PIN gene in maize roots mainly controlled IAA transport, and most of its genes were up-regulated in each period to promote maize root development.

Together, plants promote their own growth through this ROS-mediated interaction with B. aryabhattai LAD, and plants interact with a wide variety of bacteria in nature. The composition of these microbial populations is affected by microbial diversity, immune system activation, bacteria and other bacteria, fungi and phages. The interaction of other organisms and other factors, understanding each of these components can reasonably manipulate the plant microbiome to benefit plants. Therefore, these results may provide some ideas for understanding microbial-plant interactions and provide a theoretical basis for the development of biofertilizers and sustainable agriculture.

Conclusion

Here, we explored the interaction between B. aryabhattai LAD and maize roots. The results showed that the interaction between auxin-secreting bacteria B. aryabhattai LAD and maize roots promoted their connection with plants to further help plant growth, and found that there was a feedback effect between plant immune system and bacterial auxin. Bacteria activate the immune response of plant roots to produce ROS, which in turn stimulates bacteria to synthesize IAA, and the synthesized IAA further promotes plant growth. This result provides a theoretical basis for exploring the endogenous mechanism of maize development, and also provides a relevant model for exploring the effect of plant-soil microbial interaction on plant defense function and promoting plant growth.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Acknowledgements

We extend our sincere gratitude to Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China) for their expert assistance in conducting Next-Generation Sequencing (NGS) and facilitating the bioinformatic analyses, which played an instrumental role in advancing the outcomes of this research.

Author contributions

Chao Deng: Conceptualization, Investigation, Data curation, Writing - Original draft. Nan Zeng: Conceptualization, Software, Data curation, Writing - Original draft, Writing—review and editing. Chunji Li: Writing—review and editing. Jiahe Pang: Data curation. Ning Zhang: Investigation, Supervision. Bingxue Li: Supervision, Resources, Project administration, Visualization, Writing—review and editing.

Funding

This work was supported by the National Natural Science Foundation of China (31271818); International Cooperation Project of Universities in Liaoning Province (2023-01); Liaoning Province Rural Science and Technology Special Action Project (2022030673-JH5/104); Shenyang Science and Technology Project (22-319-2-13), and China Scholarship Council (CSC202208850002).

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. The raw RNA-seq data has been formally deposited with the National Center for Biotechnology Information under project PRJNA1091454.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Chao Deng and Nan Zeng contributed equally to this work.
==== Refs
References

1. Prakash R, Subramani R, Krodi, Anusha, Berde CV, Chandrasekhar T, Prathyusha AMVN, Kariali E, Bramhachari PV. Rhizobacteriome: Plant Growth-Promoting Traits and Its Functional Mechanism in Plant Growth, Development, and Defenses. 2022.
2. Mashabela MD, Piater LA, Dubery IA, Tugizimana F, Mhlongo MI. Rhizosphere tripartite interactions and PGPR-Mediated metabolic reprogramming towards ISR and Plant Priming: a Metabolomics Review. Biology. 2022;11(3).
3. Kapadia C, Patel N, Rana A, Vaidya H, Alfarraj S, Ansari MJ, Gafur A, Poczai P, Sayyed RZ. Evaluation of plant growth-promoting and salinity ameliorating potential of halophilic Bacteria isolated from saline soil. Front Plant Sci. 2022;13.
4. Upadhyay MMSD. PGPR-mediated induction of systemic resistance and physiochemical alterations in plants against the pathogens: current perspectives. J Basic Microbiology: Int J Morphology Physiol Genet Ecol Microorganisms. 2020;60(10).
5. Khan MS, Musarrat J, Zaidi A. Role of phosphate-solubilizing Bacteria in Legume Improvement. Springer Vienna. 2010;Chapter 11273–92. 10.1007/978-3-211-99753-6.
6. Jalmi SK Rhizosphere signaling nurturing phyto-microbiome niche Trop Plant Res 2020 7 2 522 8 10.22271/tpr.2020.v7.i2.064
Jalmi SK. Rhizosphere signaling nurturing phyto-microbiome niche. Trop Plant Res. 2020;7(2):522–8.10.22271/tpr.2020.v7.i2.064
7. Nelson MS, Sadowsky MJ. Secretion systems and signal exchange between nitrogen-fixing rhizobia and legumes. Front Plant Sci. 2015;6.
8. Zhang X Zhang R Gao J Wang X Fan F Ma X Yin H Zhang C Feng K Deng Y Thirty-one years of rice-rice-green manure rotations shape the rhizosphere microbial community and enrich beneficial bacteria Soil Biol Biochem 2017 104 208 17 10.1016/j.soilbio.2016.10.023
Zhang X, Zhang R, Gao J, Wang X, Fan F, Ma X, Yin H, Zhang C, Feng K, Deng Y. Thirty-one years of rice-rice-green manure rotations shape the rhizosphere microbial community and enrich beneficial bacteria. Soil Biol Biochem. 2017;104:208–17.10.1016/j.soilbio.2016.10.023
9. Tzipilevich E Russ D Dangl JL Benfey PN Plant immune system activation is necessary for efficient root colonization by auxin-secreting beneficial bacteria Cell Host Microbe 2021 29 10 1507 20 10.1016/j.chom.2021.09.005 34610294
Tzipilevich E, Russ D, Dangl JL, Benfey PN. Plant immune system activation is necessary for efficient root colonization by auxin-secreting beneficial bacteria. Cell Host Microbe. 2021;29(10):1507–20.34610294 10.1016/j.chom.2021.09.005
10. DanCouo CZ. Regulation of pattern recognition receptor signalling in plants. Nat Rev Immunol. 2016.
11. Zipfel C Plant pattern-recognition receptors Trends Immunol 2014 35 7 345 51 10.1016/j.it.2014.05.004 24946686
Zipfel C. Plant pattern-recognition receptors. Trends Immunol. 2014;35(7):345–51.24946686 10.1016/j.it.2014.05.004
12. Zipfel C Early molecular events in PAMP-triggered immunity Curr Opin Plant Biol 2009 12 4 414 20 10.1016/j.pbi.2009.06.003 19608450
Zipfel C. Early molecular events in PAMP-triggered immunity. Curr Opin Plant Biol. 2009;12(4):414–20.19608450 10.1016/j.pbi.2009.06.003
13. Spoel SH Dong X How do plants achieve immunity? Defence without specialized immune cells Nat Rev Immunol 2012 12 2 89 100 10.1038/nri3141 22273771
Spoel SH, Dong X. How do plants achieve immunity? Defence without specialized immune cells. Nat Rev Immunol. 2012;12(2):89–100.22273771 10.1038/nri3141
14. Xin XF, Kvitko B, He SY. Pseudomonas syringae: what it takes to be a pathogen. Nat Rev Microbiol. 2018.
15. Stringlis, Ioannis A, Proietti, Silvia. Hickman, Richard, Van, Verk, Marcel. Root transcriptional dynamics induced by beneficial rhizobacteria and microbial immune elicitors reveal signatures of adaptation to mutualists. Plant J. 2018.
16. Chen T Nomura K Wang X Sohrabi R Xu J Yao L Paasch BC Ma L Kremer J Cheng Y A plant genetic network for preventing dysbiosis in the phyllosphere Nature 2020 580 7805 653 7 10.1038/s41586-020-2185-0 32350464
Chen T, Nomura K, Wang X, Sohrabi R, Xu J, Yao L, Paasch BC, Ma L, Kremer J, Cheng Y, et al. A plant genetic network for preventing dysbiosis in the phyllosphere. Nature. 2020;580(7805):653–7.32350464 10.1038/s41586-020-2185-0
17. Deng C Zhang N Liang X Huang T Li B Bacillus aryabhattai LAD impacts rhizosphere bacterial community structure and promotes maize plant growth J Sci Food Agric 2022 102 14 6650 7 10.1002/jsfa.12032 35603593
Deng C, Zhang N, Liang X, Huang T, Li B. Bacillus aryabhattai LAD impacts rhizosphere bacterial community structure and promotes maize plant growth. J Sci Food Agric. 2022;102(14):6650–7.35603593 10.1002/jsfa.12032
18. Deng C, Liang X, Zhang N, Li B, Wang X, Zeng N. Molecular mechanisms of plant growth promotion for methylotrophic Bacillus aryabhattai LAD. Front Microbiol. 2022;13.
19. Li C Cheng P Li Z Xu Y Sun Y Qin D Yu G Transcriptomic and metabolomic analyses provide insights into the enhancement of Torulene and Torularhodin Production in Rhodotorula glutinis ZHK under Moderate Salt conditions J Agric Food Chem 2021 69 38 11523 33 10.1021/acs.jafc.1c04028 34545740
Li C, Cheng P, Li Z, Xu Y, Sun Y, Qin D, Yu G. Transcriptomic and metabolomic analyses provide insights into the enhancement of Torulene and Torularhodin Production in Rhodotorula glutinis ZHK under Moderate Salt conditions. J Agric Food Chem. 2021;69(38):11523–33.34545740 10.1021/acs.jafc.1c04028
20. Zeng N, Zhang N, Ma X, Wang Y, Zhang Y, Wang D, Pu F, Li B. Transcriptomics Integrated with metabolomics: assessing the Central Metabolism of different cells after cell differentiation in Aureobasidium pullulans NG. J Fungi (Basel). 2022;8(8).
21. Li C Zhao D Yan J Zhang N Li B Metabolomics integrated with transcriptomics: assessing the central metabolism of marine red yeast Sporobolomyces pararoseus under salinity stress Arch Microbiol 2021 203 3 889 99 10.1007/s00203-020-02082-9 33074377
Li C, Zhao D, Yan J, Zhang N, Li B. Metabolomics integrated with transcriptomics: assessing the central metabolism of marine red yeast Sporobolomyces pararoseus under salinity stress. Arch Microbiol. 2021;203(3):889–99.33074377 10.1007/s00203-020-02082-9
22. Paponov IA, Palme K, Teale WD. Auxin in action: signalling, transport and the control of plant growth and development. 2013.
23. Spaepen V, Remans. Indole-3-acetic acid in microbial and microorganism-plant signaling. FEMS Microbiol Rev. 2007;2007,31(4)(-):425–448.
24. Costacurta A Vanderleyden J Synthesis of Phytohormones by Plant-Associated Bacteria Crc Crit Reviews Microbiol 1994 21 1 1 18 10.3109/10408419509113531
Costacurta A, Vanderleyden J. Synthesis of Phytohormones by Plant-Associated Bacteria. Crc Crit Reviews Microbiol. 1994;21(1):1–18.10.3109/10408419509113531
25. Kunkel BN Harper CP The roles of auxin during interactions between bacterial plant pathogens and their hosts J Exp Bot 2018 69 2 245 54 10.1093/jxb/erx447 29272462
Kunkel BN, Harper CP. The roles of auxin during interactions between bacterial plant pathogens and their hosts. J Exp Bot. 2018;69(2):245–54.29272462 10.1093/jxb/erx447
26. Fan J Zheng L Acclimation to NaCl and light stress of heterotrophic Chlamydomonas reinhardtii for lipid accumulation J Biosci Bioeng 2017 124 3 302 8 10.1016/j.jbiosc.2017.04.009 28483385
Fan J, Zheng L. Acclimation to NaCl and light stress of heterotrophic Chlamydomonas reinhardtii for lipid accumulation. J Biosci Bioeng. 2017;124(3):302–8.28483385 10.1016/j.jbiosc.2017.04.009
27. Liu Z, Zhou Z, Wang L, Li M, Wang W, Yi Q, Huang S, Song L. Dopamine and serotonin modulate free amino acids production and Na+/K + pump activity in Chinese Mitten crab Eriocheir sinensis under Acute Salinity stress. Front Physiol. 2018; 9.
28. Lee M Park JC Kim D Kang S Shin K Park HG Han J Lee J Interrelationship of salinity shift with oxidative stress and lipid metabolism in the monogonont rotifer Brachionus Koreanus Comp Biochem Physiol A: Mol Integr Physiol 2017 214 79 84 10.1016/j.cbpa.2017.09.014 28951139
Lee M, Park JC, Kim D, Kang S, Shin K, Park HG, Han J, Lee J. Interrelationship of salinity shift with oxidative stress and lipid metabolism in the monogonont rotifer Brachionus Koreanus. Comp Biochem Physiol A: Mol Integr Physiol. 2017;214:79–84.28951139 10.1016/j.cbpa.2017.09.014
29. Quintana-Cabrera R Bolaños JP Glutathione and γ-glutamylcysteine in the antioxidant and survival functions of mitochondria Biochem Soc Trans 2013 41 1 106 10 10.1042/BST20120252 23356267
Quintana-Cabrera R, Bolaños JP. Glutathione and γ-glutamylcysteine in the antioxidant and survival functions of mitochondria. Biochem Soc Trans. 2013;41(1):106–10.23356267 10.1042/BST20120252
30. Iskusnykh IY, Zakharova AA, Pathak D. Glutathione in Brain disorders and Aging. Molecules. 2022; 27(1).
31. Ahmadi FI Karimi K Struik PC Effect of exogenous application of methyl jasmonate on physiological and biochemical characteristics of Brassica napus L. Cv. Talaye under salinity stress S Afr J Bot 2018 115 5 11 10.1016/j.sajb.2017.11.018
Ahmadi FI, Karimi K, Struik PC. Effect of exogenous application of methyl jasmonate on physiological and biochemical characteristics of Brassica napus L. Cv. Talaye under salinity stress. S Afr J Bot. 2018;115:5–11.10.1016/j.sajb.2017.11.018
32. Xiao X, Wang Q, Ma X, Lang D, Guo Z, Zhang X. Physiological biochemistry-combined transcriptomic analysis reveals mechanism of Bacillus cereus G2 improved salt-stress tolerance of Glycyrrhiza Uralensis Fisch. Seedlings by balancing Carbohydrate Metabolism. Front Plant Sci. 2022; 12.
33. Zeng N Zhang N Wang D Long J Wang Y Zhang Y Pu F Li Z Baloch FB Li B Regulation of cell differentiation to promote pullulan synthesis in Aureobasidium pullulans NG Appl Microbiol Biotechnol 2023 107 22 6761 73 10.1007/s00253-023-12758-z 37698607
Zeng N, Zhang N, Wang D, Long J, Wang Y, Zhang Y, Pu F, Li Z, Baloch FB, Li B. Regulation of cell differentiation to promote pullulan synthesis in Aureobasidium pullulans NG. Appl Microbiol Biotechnol. 2023;107(22):6761–73.37698607 10.1007/s00253-023-12758-z
34. Nandakumar M Nathan C Rhee KY Isocitrate lyase mediates broad antibiotic tolerance in Mycobacterium tuberculosis Nat Commun 2014 5 1 4306 10.1038/ncomms5306 24978671
Nandakumar M, Nathan C, Rhee KY. Isocitrate lyase mediates broad antibiotic tolerance in Mycobacterium tuberculosis. Nat Commun. 2014;5(1):4306.24978671 10.1038/ncomms5306
35. Lorenz MC, Fink GR, Lorenz MC, Fink GR. The glyoxylate cycle is required for fungal virulence. Nature 412, 8386. Nature. 2001;412(6842):83–86.
36. Li L Wang F Yan P Jing W Zhang C Kudla J Zhang W A phosphoinositide-specific phospholipase C pathway elicits stress-induced Ca2 + signals and confers salt tolerance to rice New Phytol 2017 214 3 1172 87 10.1111/nph.14426 28157263
Li L, Wang F, Yan P, Jing W, Zhang C, Kudla J, Zhang W. A phosphoinositide-specific phospholipase C pathway elicits stress-induced Ca2 + signals and confers salt tolerance to rice. New Phytol. 2017;214(3):1172–87.28157263 10.1111/nph.14426
37. Santos T Marinho HS Cyrne L Regulation of the inositol transporter Itr1p by hydrogen peroxide in Saccharomyces cerevisiae Arch Microbiol 2019 201 1 123 34 10.1007/s00203-018-1584-1 30283989
Santos T, Marinho HS, Cyrne L. Regulation of the inositol transporter Itr1p by hydrogen peroxide in Saccharomyces cerevisiae. Arch Microbiol. 2019;201(1):123–34.30283989 10.1007/s00203-018-1584-1
38. Yang L Zhang J He J Qin Y Hua D Duan Y Chen Z Gong Z ABA-Mediated ROS in Mitochondria regulate Root Meristem Activity by Controlling PLETHORA expression in Arabidopsis PLoS Genet 2014 10 12 e1004791 10.1371/journal.pgen.1004791 25522358
Yang L, Zhang J, He J, Qin Y, Hua D, Duan Y, Chen Z, Gong Z. ABA-Mediated ROS in Mitochondria regulate Root Meristem Activity by Controlling PLETHORA expression in Arabidopsis. PLoS Genet. 2014;10(12):e1004791.25522358 10.1371/journal.pgen.1004791
39. Correa-Aragunde N, Foresi N, Delledonne M, Lamattina L. Auxin induces redox regulation of ascorbate peroxidase 1 activity by S-nitrosylation/denitrosylation balance resulting in changes of root growth pattern in Arabidopsis. J Exp Bot. 2013(11).
40. Jez JM Connecting primary and specialized metabolism: amino acid conjugation of phytohormones by GRETCHEN HAGEN 3 (GH3) acyl acid amido synthetases Curr Opin Plant Biol 2022 66 102194 10.1016/j.pbi.2022.102194 35219141
Jez JM. Connecting primary and specialized metabolism: amino acid conjugation of phytohormones by GRETCHEN HAGEN 3 (GH3) acyl acid amido synthetases. Curr Opin Plant Biol. 2022;66:102194.35219141 10.1016/j.pbi.2022.102194
41. Qin M, Wang J, Zhang T, Hu X, Liu R, Gao T, Zhao S, Yuan Y, Zheng J, Wang Z. Genome-wide identification and analysis on YUCCA Gene Family in Isatis Indigotica Fort. And IiYUCCA6-1 Functional Exploration. Int J Mol Sci. 2020; 21(6).
42. Marhava P Recent developments in the understanding of PIN polarity New Phytol 2022 233 2 624 30 10.1111/nph.17867 34882802
Marhava P. Recent developments in the understanding of PIN polarity. New Phytol. 2022;233(2):624–30.34882802 10.1111/nph.17867
