text:
Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Bacillus subtilis biofilm matrix components target seed oil bodies to promote growth and anti-fungal resistance in melon

## Abstract

Beneficial microorganisms are used to stimulate the germination of seeds; however, their growth-promoting mechanisms remain largely unexplored. Bacillus subtilis is commonly found in association with different plant organs, providing protection against pathogens or stimulating plant growth. We report that application of B. subtilis to melon seeds results in genetic and physiological responses in seeds that alter the metabolic and developmental status in 5-d and 1-month-old plants upon germination. We analysed mutants in different components of the extracellular matrix of B. subtilis biofilms in interaction with seeds and found cooperation in bacterial colonization of seed storage tissues and growth promotion. Combining confocal microscopy with fluorogenic probes, we found that two specific components of the extracellular matrix, amyloid protein TasA and fengycin, differentially increased the concentrations of reactive oxygen species inside seeds. Further, using electron and fluorescence microscopy and metabolomics, we showed that both TasA and fengycin targeted the oil bodies in the seed endosperm, resulting in specific changes in lipid metabolism and accumulation of glutathione-related molecules. In turn, this results in two different plant growth developmental programmes: TasA and fengycin stimulate the development of radicles, and fengycin alone stimulate the growth of adult plants and resistance in the phylloplane to the fungus Botrytis cinerea. Understanding mechanisms of bacterial growth promotion will enable the design of bespoke growth promotion strains.

## Main

Seed germination is a complex biological process that is controlled by interconnected hormone-regulated pathways1. Seed inoculation with plant growth-promoting rhizobacteria is commonly used to enable bacterial colonization of plants and beneficial effects of plant growth-promoting rhizobacteria post germination2. A detailed understanding of the mechanisms that underpin growth-promoting host–microbe interactions is important to enable design of innovative biotechnological probiotics.

Bacillus subtilis and closely related species coexist with plants, providing multiple beneficial services3. Production of polyvalent secondary metabolites, sporulation or biofilm formation underpin bacterial fitness and bioactivity towards plants4. Biofilms are formed by bacterial cells that are embedded in a secreted extracellular matrix (ECM)5 that is essential for efficient colonization of plant organs6. The ECM also comprises secondary metabolites that mediate bacterial communication with the plants by triggering physiological responses associated with defence or growth7.

Here we investigated the mechanism by which melon seeds respond to the stimulatory activity of B. subtilis and evaluated the functions of the ECM in this microbe–host interaction.

## Results

### B. subtilis regulates metabolism and growth of colonized plants

Beneficial bacteria can promote two different and genetically controlled stages of seed germination, namely germination itself and growth of the emergent radicle8. Melon seeds treated with B. subtilis NCIB 3610 produced larger radicles compared with those grown from untreated seeds (Fig. 1a). However, the germination rates (initial emergence of the radicle) of treated seeds did not change compared with untreated seeds. We analysed the expression levels of GA20ox1 and CYP07A1, which are two genes involved in the modulation of the determinant ratio of abscisic acid and gibberellins9, and found that they are statistically similar between treated and untreated seeds (Extended Data Fig. 1a,b). Analysis of transcriptomes of three pools of seeds 16 h after treatment with a suspension of B. subtilis did not identify any changes in the expression levels of the genes involved in the germination-related hormone signalling pathway. However, the upregulation of plant genes involved in carbon metabolism and photosynthesis, together with repression of plant heat shock proteins and structural proteins of lipid storage vesicles (oleosin and caleosin), was detected and enabled activation of seed metabolism (Fig. 1b and Extended Data Fig. 1c,d)10,11,12.

Metabolomic profiles of treated and untreated seeds were compared to reveal increased abundance of energetic resources, such as glycerophospholipids and fatty acyls with analogy to triacylglycerides (TAGs) in the initial stages (0 h) after treatment with B. subtilis cells (Fig. 1c and Extended Data Fig. 2a). At 24 h after the treatment, metabolite levels were similar between both conditions, although specific changes in the abundance of metabolites within these lipidic classes remained clear (Extended Data Fig. 2a).

These post-germination stimulatory effects of treatment with B. subtilis were long-lasting, considering that 1-month-old adult plants emerging from treated seeds developed a more vigorous radicular system and canopy than those emerging from untreated seeds (Fig. 1c). We analysed the metabolome of three adult plants that emerged from treated seeds to define putative metabolic changes that were associated with this long-lasting growth-promoting effect. Carboxylic acids, lipids and lipid-like molecules and organooxygen compounds were the main classes of metabolites differentially detected in aerial regions (leaves and stem) of plants grown from treated and untreated seeds (Extended Data Fig. 2b), while metabolomic composition of the roots was not significantly different between treated and untreated groups.

Changes in fatty acyls, carboxylic acids, organooxygen compounds and prenol lipids were the main metabolic signatures of the leaves from plants grown after bacterial treatment of the seeds (Fig. 1e). Our metabolomic analysis also revealed differential accumulation of l-tryptophan and cinnamic acids in leaves of adult plants derived from treated seeds (Extended Data Fig. 2c). Both families of molecules are biomarkers of beneficial plant–bacteria interactions13,14.

### B. subtilis extracellular matrix and growth promotion

We investigated the growth dynamics of the B. subtilis population during the first 5 d after seed treatment in the seed and emergent radicle (Fig. 2a) by colony-forming unit (c.f.u.) plating over seed and radicle extracts. In the radicles, the population remained unchanged and almost every cell sporulated upon radicle emergence (2 d after treatment). In the seeds, however, the bacterial population increased during the first 5 d post treatment and became almost entirely sporulated 4 d after bacterization. These results suggest that B. subtilis cells, mostly spores, may have been passively dragged by the emergent radicle or might have colonized and proliferated inside seeds. We used scanning electron microscopy analysis (SEM) of treated seeds to confirm that B. subtilis colonizes the inner regions of the seeds (Extended Data Fig. 3a). Confocal laser scanning microscopy (CLSM) analyses of transversely sectioned seeds previously treated with fluorescently labelled B. subtilis (CellTracker cm-DII) indicated that bacterial cells accumulate in the storage tissues near the seed micropyle, which is a natural entry point for bacteria into seeds15 (Extended Data Fig. 3b). The cell densities of the wild-type (WT) strain and ∆hag (flagellum), ∆eps or ∆srf mutant strains, which are known to have altered swimming, sliding or swarming motility, respectively16,17,18, were unchanged in two differentiated parts inside the seeds (Fig. 2b,c). Overall, these findings reveal that growth-promoting activity might be triggered by bacterial cells that enter and colonize the seed storage tissues, which is a process that does not appear to rely specifically on any one type of bacterial motility.

Several studies have reported multifaceted contributions of bacterial biofilms to interactions with hosts6,19,20,21. The cell densities of single structural mutants ∆tasA, ∆tapA, ∆eps and ∆bslA, which form altered biofilms, were substantially decreased in seed extracts compared with that of the wild-type bacteria (set at 100%) 5 d after seed treatment (Fig. 2d). However, this result was not associated with the attachment of bacterial cells to the seeds a few hours after treatments (Extended Data Fig. 3c).

A similar bacterial population pattern was monitored in the emergent radicles, which, as described above for WT, may reflect dragging along of bacterial cells (Extended Data Fig. 3d). The ∆srf, ∆pks or ∆pps strains, which do not produce the secondary metabolites surfactin, bacillaene or fengycin, respectively (also present in the ECM), had unchanged population dynamics compared with the WT strain (Fig. 2d and Extended Data Fig. 3d). The decrease in the bacterial cell population size found in the structural mutants suggests that mutations led to a significant decrease in radicle growth promotion compared with treatment with the wild-type strain (Fig. 2e). In contrast to a correlation between bacterial cell density and plant growth promotion, the poorly persisting ∆tasA strain retained promoting activity, while the ∆pks and ∆pps strains failed to promote radicle growth despite considerable persistence on the seeds (Fig. 2d,e).

The absence of the amyloid protein TasA has recently been reported to provoke important cellular and physiological changes in B. subtilis6. To clarify whether the role of TasA in the structure of the ECM is related to TasA-dependent growth-promoting activity, we used B. subtilis strain JC81, which expresses a version of TasA that fails to restore biofilm formation but reverts the physiological status of ∆tasA cells to wild-type levels6. Treatment of seeds with JC81 failed to promote radicle growth, but the bacteria persisted on the seeds at a level comparable to that of the WT, ∆pks or ∆pps strains (Fig. 2d,e).

These results led us to propose three distinct and complementary contributions of ECM components to the promotion of seed radicle growth: first, EPS, BslA and TapA have a role in the persistence of B. subtilis cells in the seeds; second, fengycin (pps) and bacillaene (pks) have roles in the chemical dialogue of Bacillus with seeds; and third amyloid TasA is involved in both of these functions.

Comparison of the metabolic status of radicles 5 d after treatment of seeds with WT or various ECM mutant strains demonstrated a specific metabolic response of plants to chemicals present in the B. subtilis ECM (Extended Data Fig. 4a). The changes associated with bacterial treatment of seeds included a decrease in organooxygen compound analogues of sphinganines and an increase in abundance of the prenol lipid and stilbene molecular families (Extended Data Fig. 4b). The presence of a functional ECM in the WT strain triggered the accumulation of fatty acyls belonging to two different clusters, and treatments with ECM mutants produced a clear decrease in organooxygen compounds and accumulation of specific prenol lipids (Extended Data Fig. 4b).

### Lysophospholipids, glutathione and growth promotion

According to our data, TasA and fengycin are highly relevant molecules of the ECM with regards to growth promotion of B. subtilis in seeds. TasA, in common with other amyloids, is a polymorphic protein that can adopt a variety of structural conformations that have different biochemistries and functions6,22,23. We polymerized purified homogeneous TasA monomers to show that large aggregates stimulate radicle growth (Extended Data Fig. 5a). Moreover, fresh apoplast fluid extracted from melon seeds promoted the polymerization of TasA into aggregates (Extended Data Fig. 5b), indicating that this largely active polymerized form of TasA may predominate inside seeds. Evaluation of stimulatory activity showed that a 3 µM solution of the most active form of TasA or a 10 µM solution of purified fengycin significantly increased the area of the radicles compared with untreated seeds (Fig. 3a). The lack of stimulatory activity of purified BslA or surfactin confirmed our previous findings using B. subtilis mutants (Figs. 2e and 3a).

The metabolome of individual radicles (n = 3) that emerged from seeds treated with fengycin or TasA was analysed to define their contributions to the metabolic signatures associated with the promotion of radicle growth. Results of principal component analysis and heatmap analyses indicated sample clustering in three groups (Extended Data Fig. 6a,b). Statistical analysis by partial least-squares discriminant analysis indicated that glycerophospholipids, fatty acyls, organooxygen compounds or carboxylic acids were discriminating metabolites in the case of TasA or fengycin treatments, and feature-based molecular networking of these features showed integration into molecular families composed of other metabolites that generally followed the same abundance patterns (Extended Data Fig. 6c–e). Further refinement of this analysis identified two molecular families of interest: lysophosphatidylethanolamine (LysoPE), mainly associated with TasA treatment, and an analogue of reduced glutathione (GSH) that accumulated in the radicles grown from seeds treated with fengycin (Fig. 3b). The treatment of seeds with commercially available LysoPE or GSH increased radicle growth as much as treatment with fengycin or TasA, but did not produce the long-term promoting effects observed in plants 1 month after treatment of seeds with B. subtilis (Fig. 3c).

GSH is an essential and polyvalent metabolite with special antioxidant functions, which maintains cellular redox homoeostasis and development, growth or response of plants to a variety of stimuli24. Further, an increase in the total pool of glutathione usually occurs 3 to 4 d after ROS exposure24,25, therefore we hypothesized that one possible outcome of fengycin treatment might be ROS production. The level of ROS in seeds after treatment with fengycin was statistically higher than that in untreated seeds. An intermediate increase in ROS, which was statistically less than that elicited by fengycin, was also observed in seeds treated with TasA (Fig. 3d).

### TasA and fengycin differentially target seed oil bodies

Accumulation of bacterial cells in the micropylar endosperm (Extended Data Fig. 3b), accumulation of TAGs (Fig. 1c) and changes in the lipid composition of the radicles and adult plants after treatment of seeds with bacteria or ECM components (Extended Data Figs. 2b and 6d,e) led us to hypothesize that seed oil bodies (OBs), which are nutrient reservoir organelles mainly composed of TAGs and other neutral lipids, are the target of fengycin and TasA. During germination, OBs are degraded by interaction with glyoxysomes that feed the plant embryo26,27. We applied CLSM to thin sections of seeds treated with a suspension of B. subtilis or untreated seeds using a neutral lipid stain 16 h after seed treatment and found an abundance of OBs in the endosperm of untreated seeds. In agreement with our hypothesis, certain regions of the endosperm accumulated high levels of disaggregated oil bodies that were surrounded by bacterial aggregates (Extended Data Fig. 7a).

Treatment with purified TasA or fengycin promoted in vivo disaggregation of OBs, similar to the effect observed after treatment with Bacillus cells (Extended Data Fig. 8a). Transmission electron microscopy (TEM) analysis of thin sections of negatively stained seeds showed OB disaggregation and localization of OB-surrounding glyoxysomes in samples treated with TasA (Extended Data Fig. 7b,c). In fengycin-treated seeds, the most predominant changes were disaggregation of OBs and the presence of empty areas between the cell wall and the cytoplasm, suggestive of plasmolysis, a cellular response previously observed in the interaction of fengycin with membranes of fungal cells28 (Extended Data Fig. 7c). The CLSM analysis of OBs purified from melon seeds treated with TasA or fengycin and double-stained with Fast Green FCF and Nile red indicated differences in the size and pattern of aggregation of OBs: fengycin induced disaggregation of OBs as observed in vivo, while TasA preferentially induced a higher level of aggregation of individual OBs than that observed in vivo (Fig. 4a). Higher intensity of the green fluorescence signal near OB aggregates in the samples treated with TasA suggested non-specific staining of the amyloid protein and its localization around the vesicles. TEM analysis of purified OBs and immunogold labelling with anti-TasA antibodies indicated the presence of gold particles decorating the perimeter of the OBs and interconnecting TasA fibres (Extended Data Fig. 7d). The results of pulldown assays of whole protein extracts of the seeds indicated co-elution of purified TasA with oleosins, which are structural proteins that modulate the stability and size of OBs, regulate lipid metabolism and have an important role in OB degradation29,30 (Fig. 4b and Supplementary Table 5). Interaction of the two proteins was confirmed by far-western blot analysis of purified TasA and purified oleosin 1 from Arabidopsis thaliana (Fig. 4c).

We propose that TasA–oleosin interaction is responsible for the aggregation of OBs, favouring the accumulation of lysophospholipids after the optimization of lipid catabolism in glyoxysomes. Fengycin has been described to efficiently interact with and disrupt artificial lipid monolayers or bilayer membranes in a concentration-dependent manner31. In addition to micellar concentrations of fengycin, specific disruption of the membranes relies on the lipid composition of the target32. Exact chemical composition of phospholipids of the membranes of OBs in the seeds is not known and varies between species33; however, we propose that high affinity of fengycin for lipid membranes explains disaggregation and a reduction in the size of OBs.

Non-targeted metabolomic analyses provided additional information about the change in relative abundance of metabolites and potential (bio)chemical modifications of seed metabolites from 0 to 24 h after treatment with TasA or fengycin. The use of the chemical proportionality tool (ChemProp34) showed putative changes in the dynamics of TAGs analogues after seed treatments (Fig. 4d). In control seeds, many modifications in this molecular family denoted an important catalytic activity occurring over key storage metabolites in this time period. Treatments with TasA and, more notably, with fengycin led to an increase in the initial abundance of many metabolites (Fig. 4d, floating bar plots) that would reflect the faster metabolic activity of these seeds. A decrease in the Chemprop score indicated that the significance of the changes between 0 h and 24 h would be directly related with the earlier initiation of the catalytic activity over these molecules. Similar changes to TAGs were also remarkable in the dynamics of glycerophospholipids (Extended Data Fig. 8b).

### Fengycin promotes plant growth and resistance to Botrytis cinerea

Plants that emerged from seeds treated with fengycin were larger than plants grown from untreated or TasA-treated seeds (Fig. 5a). Adult plants grown from seeds treated with LysoPE or GSH, however, did not show a significant increase in growth over time compared with plants from untreated seeds (Fig. 3c). These findings indicate that the presence of different signalling events underpin short-term radicle growth mediated by TasA, fengycin or associated LysoPE and GSH molecules, as well as long-term growth of adult plants specifically associated with fengycin. GSH associated with treatments of the seed with fengycin did not produce a sustained long-term growth-promoting effect, suggesting that a constant endogenous trigger is required that can only be achieved by fengycin treatment. Both TasA and fengycin increased the level of ROS inside seeds, but the highest response was elicited by fengycin (Fig. 3d). ROS production is a normal response after the first water uptake during seed imbibition and, in a certain range of concentrations, ROS has a well-recognized role in endosperm weakening, protection against pathogens, mobilization of seed reserves and interaction with plant hormones during the process of seed germination35,36,37. These results suggest that the increase in ROS levels after seed treatment with fengycin, but not TasA, may act as a beneficial stimulus for further plant development mediated by GSH accumulation.

Previous studies revealed that application of fengycin elicits plant defence responses in adult plants7,38. Therefore, we investigated whether treatment of seeds with fengycin immunized adult plants against aboveground pathogens by inoculating the third leaf of adult plants with a spore suspension of Botrytis cinerea. The size of necrotic lesions induced by the fungus in plants grown from fengycin-treated seeds was significantly smaller than that of control plants or plants grown from TasA-treated seeds (Fig. 5b). Transcriptomic analyses of infected leaves revealed a high basal level of transcripts, some of them related to plant defence mechanisms such as Allene oxide cyclase/synthase, directly implied in the biosynthesis of jasmonic acid39 (Supplementary Table 5) in leaves of plants emerged from seeds treated with fengycin. Their expression was maintained or increased 48 h after inoculation with B. cinerea, while an increase was found in adult plants from non-treated seeds as a result of the response to the infection. Second, upregulation of the genes related to glutathione metabolism, specifically including glutathione S-transferases (GSTs), occurred 48 h after the challenge with B. cinerea (Fig. 5c and Extended Data Fig. 9a). Participation of GSTs in antioxidant reactions, together with the important cellular antioxidant GSH, has been proven to mitigate oxidative stress caused by necrotrophic fungus B. cinerea in infected tissues40. Metabolomic data from the radicles of fengycin-treated seeds demonstrated an increase in the levels of certain molecules of the same cluster that included GSH (Fig. 3b), and adult plants grown from bacterized seeds accumulated certain molecules belonging to the GSH cluster, especially in the aerial region (Extended Data Fig. 9b). In addition, metabolomes of aerial parts of seedlings 5 d after treatment with fengycin also showed an increase in oxidized glutathione and glutathione-related metabolites in addition to defence-related compounds such as flavonoids (that is, mangiferin and 2-O-rhamnosilvitexin) (Fig. 5d and Extended Data Fig. 9c).

## Discussion

In this study, we demonstrated that in addition to a structural role of bacterial ECM in biofilm formation, it is also essential for Bacillus colonization of melon seeds and growth promotion, including induction of anti-fungal resistance in the adult plant. Our findings indicate that the beneficial effects of B. subtilis treatment on melon seeds are associated with changes in the initial ROS increase during seed imbibition, and the contents and specific pools of metabolites released from storage tissues, which were mediated by at least two ECM components, fengycin and TasA.

We propose that specific interactions of TasA or fengycin with OBs defines two different physiological responses. First, a short-term effect on radicle growth, mediated by both components and related to changes in the dynamics of TAGs after differential targeting of oil bodies and lipid catabolism optimization. On one hand, the oleosin–TasA interaction determines the OB aggrupations around glyoxysomes and the subsequent accumulation of lysophospholipids acting as signal molecules for growth. On the other hand, the OB disaggregations caused by the action of fengycin over membranes determine a faster mobilization of their contents. Second, a long-term effect of growth promotion that persists in adult plants and is mediated by fengycin, through an initial increase in ROS that causes accumulation of GSH that is maintained over time and makes the plant more tolerant to biotic stresses (Fig. 6).

Glycerophospholipids have both structural and signalling roles in plants, and this bifunctionality is in part due to continuous synthesis and turnover of endogenous pools41. Accumulating evidence suggests a role for lysophospholipid derivatives in signalling processes in plant cells. LysoPE treatment has been reported to delay fruit softening when used postharvest, mitigate defoliation effects of ethephon, and delay leaf and fruit senescence in tomato and potato42,43 and is used commercially as a plant bioregulator to improve plant product quality41,44. On the other hand, glutathione is an essential metabolite with antioxidant properties, and is involved in cellular redox homoeostasis and many other physiological processes such as the development, growth and environmental response of plants24. The accumulation of glutathione-related metabolites in radicles and aerial parts of seedlings, and the higher expression of GSTs in leaves after seed treatment with fengycin, suggest that initial treatment causing an increase in ROS might act as a beneficial stimulus for plant development, and confer enhanced antioxidant capacity to mitigate an imbalance in the redox status in adult plants imposed by infection with B. cinerea.

We propose that only those seeds containing abundant OBs, and characterized by specific morphology in their primordial tissues, respond to the beneficial interaction with B. subtilis, and that this interaction is mediated at least in part by fengycin and TasA. According to this model, wheat or maize, which are monocotyledonous plants whose seeds are composed of starchy endosperm surrounded by a layer of living cells, did not react to the presence of fengycin, whereas cucumber or soybean seeds, which are anatomically different but do contain OBs, responded with a growth-promoting effect (Supplementary Table 6).

We note however that B. subtilis communities (cells and ECM) still retain stimulatory activity in these and other seeds plants, so it cannot be excluded that additional plant pathways or Bacillus-derived structural or chemical components, such as secondary metabolites or phytohormones, are specifically involved in maintaining these intricate interkingdom mutualistic interactions.

## Methods

### Strains, media and culture conditions

The bacterial strains used in the present study are listed in Supplementary Table 1. Bacterial cultures were grown at 37 °C from frozen stocks on lysogeny broth (LB: 1% tryptone (Oxoid), 0.5% yeast extract (Oxoid) and 0.5% NaCl) plates. Isolated bacteria were inoculated in appropriate media. The necrotrophic fungus B. cinerea was grown at 25 °C from a frozen stock in potato dextrose agar (PDA) plates and maintained until inoculum preparation. Escherichia coli DH5α was used for cloning and plasmid replication. E. coli BL21(DE3) and BL21(AI) were used for protein purification. The final antibiotic concentrations for B. subtilis were: MLS (1 μg m l−1 erythromycin and 25 μg ml−1 lincomycin), spectinomycin (100 μg ml−1), tetracycline (10 μg ml−1), chloramphenicol (5 μg ml−1) and kanamycin (10 μg ml−1).

### Quantitative analysis of bacterial growth and promotion of radicle growth

Bacterial strains were grown in liquid LB at 30 °C overnight. The cells in the cultures were washed twice with sterile distilled water, and bacterial cell suspensions were adjusted to the same optical density (OD600). When using a fluorescent inoculum, CellTracker CM-DiI dye (Invitrogen, C7000) was added to the bacterial suspension. Melon seeds (Rochet Panal - Fitó) were surface-sterilized with 0.1% sodium hypochlorite and bacterized with bacterial cell suspensions by bathing for 1 h at room temperature. The seeds were placed in Petri dishes with wet filter paper, covered with foil and maintained in growth chambers at 25 °C for 5 d. For quantitative analysis of the bacteria, the radicles were separated from the seeds. The radicles were placed in 1.5 ml microcentrifuge tubes with 1 ml PBS medium and 1 g of 3-mm-diameter glass beads, and attached bacteria were recovered by vortexing for 1 min. The seeds were ground using a mortar and a pestle in 1 ml PBS. Total c.f.u. were assayed by plating serial dilutions of resulting suspensions on solid LB medium. C.f.u. of the spores were assayed using the same suspension after thermal shock at 80 °C for 5 min. Seeds treated with water were used as controls.

Analysis of radicle growth-promoting effects was based on radicle areas measured 5 d after seed treatment using ImageJ software. The medium radicle area of the control treatment group was assumed as 100% growth and was used to normalize the percentage of growth-promoting effect corresponding to strain treatment groups. Treatments with compounds purified from the extracellular matrix were performed in a similar manner. BslA and TasA proteins were purified as described above while fengycin, surfactin and iturin were commercially acquired (Biosynth-Carbosynth AF35268, Sigma-Aldrich S3523 and I1774, respectively). Unless indicated, TasA was used at 3 µM and fengycin at 10 µM. Iturin was used for seed treatments at 260 µM, and lysophosphoethanolamine 16:0 (Avanti Polar Lipids, 856705) and reduced glutathione (Sigma-Aldrich G4251) were used at 0.2 mg ml−1 and 0.25 mM, respectively. For the evaluation of long-term growth promotion, seedlings from treated seeds were transferred to pots and maintained in a growth chamber at 25 °C.

### Protein cloning and purification

TasA protein was purified as previously described6. BslA protein was purified from freshly transformed E. coli BL21-AI cells expressing the pDEST17 plasmid harbouring the bslA gene, which was initially subcloned into the donor vector PDONR207 using a Gateway recombination system. The primers designed to amplify the bslA gene were as follows: bslAFw: 5’GGGGACAAGTTTGTACAAAAAAGCAGGCTCAGCTGGAATCTACATCAACTAAAGC3’; bslARv: 5’GGGGACCACTTTGTACAAGAAAGCTGGTTTAGTTGCAACCGCAAGGCT3’. The colonies of E. coli BL21-AI were picked and resuspended in 500 ml LB containing 100 µg ml−1 ampicillin and incubated at 37 °C in an orbital shaker to reach an OD600 of 0.7. Then, the expression of recombinant protein was induced by 0.2% l-arabinose. After 4 h, the cells were collected and resuspended in buffer A (20 mM sodium phosphate, 500 mM NaCl and 20 mM imidazole) containing 400 µl 0.1 M phenylmethylsulfonyl fluoride, 40 µl lysozyme and 4 µl Cellytic 10x (Sigma-Aldrich). After a 45 min incubation with shaking at room temperature, the cells were disrupted by sonication (3 pulses of 1 min at an amplitude of 80%). The cell lysate was centrifuged at 12,000 g for 15 min at 20 °C, and the supernatant was filtered through a 0.45-µm-pore filter. The protein was purified using an AKTA Start fast protein liquid chromatography system (GE Healthcare). The soluble fraction was loaded into a 5 ml HisTrap HP column (GE Healthcare) previously equilibrated with buffer A. The protein was eluted from the column with elution buffer (20 mM sodium phosphate, 500 mM NaCl and 500 mM imidazole, pH 8). After affinity chromatography, the purified protein was loaded into a HiPrep 26/10 desalting column (GE Healthcare), and the buffer was exchanged for 20 mM Tris and 50 mM NaCl.

Oleosin protein was purified from inclusion bodies using freshly transformed E. coli Lemo21(DE3) cells expressing the pDEST17 plasmid harbouring the sequence of OLEOSIN1 from Arabidopsis thaliana, which was initially subcloned in the donor vector pENTR223 acquired from Arabidopsis Biological Resource Center using a Gateway recombination system. E. coli Lemo21(DE3) colonies were picked and resuspended in 500 ml LB containing 100 µg ml−1 ampicillin, 5 µg ml−1 chloramphenicol and 100 µM rhamnose, and incubated at 37 °C with orbital shaking to reach an OD600 of 0.4. Then, the expression of the recombinant protein was induced by 400 mM isopropyl β-d-1-thiogalactopyranoside (IPTG). After 3 h, the cells were collected by centrifugation (5,000 g, 15 min, 4 °C), resuspended in buffer A (50 mM Tris and 150 mM NaCl, pH 8) and then centrifuged again. The pellets were stored at −80 °C until purification or processed after 15 min. After thawing, the cells were resuspended in buffer A, sonicated on ice (3 × 45 s, 60% amplitude) and centrifuged (15,000 g, 60 min, 4 °C). The supernatant was discarded because the protein was mainly expressed in inclusion bodies. The pellet was resuspended in buffer A supplemented with 2% Triton X-100, incubated at 37 °C with shaking for 20 min and centrifuged (15,000 g, 10 min, 4 °C). The pellet was extensively washed with buffer A, centrifuged (15,000 g, 10 min, 4 °C), resuspended in denaturing buffer (50 mM Tris, 500 mM NaCl and 6 M GuHCl) and incubated at 60 °C overnight for complete solubilization. The lysates were clarified by sonication on ice (3 × 45 s, 60% amplitude) and centrifugation (15,000 g, 1 h, 16 °C), and passed through a 0.45 µm filter before affinity chromatography. The protein was purified using an AKTA Start fast protein liquid chromatography system (GE Healthcare). Soluble inclusion bodies were loaded into a 5 ml HisTrap HP column (GE Healthcare) previously equilibrated with binding buffer (50 mM Tris, 0.5 M NaCl, 20 mM imidazole and 8 M urea, pH 8). The protein was eluted from the column with elution buffer (50 mM Tris, 0.5 M NaCl, 500 mM imidazole and 8 M urea, pH 8) and maintained under denaturing conditions.

### Apoplastic fluid extraction from the seeds

The apoplastic fluid (AF) of the seeds was collected by a standard technique allowing the recovery of the components present in the intercellular spaces based on vacuum infiltration and centrifugation45. Briefly, imbibed seeds were carefully decoated and immersed in infiltration buffer (50 mM Tris-HCl, pH 7.5, and 0.6% NaCl). Five vacuum pulses of 10 s (separated by 30 s intervals) were applied using a vacuum pump. Infiltrated seeds were recovered, dried on a filter paper, placed in plastic syringe barrels inside centrifuge tubes and centrifuged for 20 min at 400 g at 4 °C. The AF was recovered and stored at 4 °C until use at various concentration ratios (0, 1:2, 1:5, 1:10, 1:20 and 1:50) in TasA polymerization assays.

### TasA polymerization assays

The kinetics of polymerization of TasA were monitored using thioflavin T (ThT) and detected by fluorescence22. The assay was carried out in a 96-well microtitre plate. Briefly, purified TasA protein was diluted in buffer (20 mM Tris and 50 mM NaCl, pH 7) or AF to final concentrations of 0.4 mg ml−1 or 0.2 mg ml−1. ThT solution was added to the wells to a final ThT concentration of 20 μM in a final volume of 100 μl. The fluorescence signal was measured in a fluorescence microplate spectrophotometer reader (438 nm excitation and 495 nm emission, with a 475 nm cut-off) at 30 °C under shaking. Each protein concentration was assayed in triplicate.

### Infection assays in the plants

Assays of B. cinerea infection were carried out in 5–6-week-old plants. Fungal conidia were collected from light-grown culture in sterile distilled water and filtered through a 40 µm cell strainer to remove remaining hyphae. For inoculation, the conidial suspension was adjusted to 105 conidia per ml in half-strength filtered (0.45 μm) grape juice (100% pure organic). Each leaf was inoculated with 5 μl droplets of conidial suspension. The pots were covered with a plastic dome and placed in a growth chamber. The leaves were imaged 72 h after the inoculation and the size of the lesions was determined using ImageJ software. The size of the lesions in experimental plants was normalized to a medium lesion size in the control plants, which was assumed as 100% disease percentage.

### Purification of oil bodies from the seeds

Oil bodies were purified as previously described46, with some modifications according to ref. 47. Briefly, the seeds were soaked in 0.1 M sodium bicarbonate at pH 9.5 (adjusted with 0.1 M NaOH) for 2 h and soaking medium was discarded. The soaked seeds were ground in the same medium with a mortar and a pestle. The slurry was transferred to 2 ml tubes and centrifuged at 10,000 g for 30 min at 4 °C. The upper layer was transferred to a clean tube, dispersed in washing solution (7 M urea, 1:4 w/w) and centrifuged (10,000 g, 30 min). To remove residue of washing solution, the fat pad was isolated and suspended in water (1:4 w/w) and centrifuged (10,000 g, 30 min). After centrifugation, the creamy layer was collected and stored at 4 °C. Analyses were completed within 24 h.

### Pulldown assay

To study possible TasA interactions with the proteins present in the seeds, we carried out a pulldown assay using a nickel-loaded affinity resin (Protino Ni-TED 2000, Macherey-Nagel); purified TasA was used as a bait protein for protein extract obtained from imbibed seeds. Imbibed and decoated seeds were ground to a fine powder using a mortar and a pestle in liquid nitrogen. The powder was transferred to a clean 15 ml tube and resuspended in extraction buffer (50 mM Tris-HCl, pH 7.5, 150 mM NaCl, 10 mM EDTA, 1% Triton X-100, 10% glycerol, 10 mM dithiothreitol (DTT), 1 mM phenylmethylsulfonyl fluoride, 1 mM NaF and 1% protease inhibitor cocktail). The samples were homogenized using a vortex, sonicated by 5 pulses of 20 s at 80% amplitude and centrifuged for 5 min at 9,000 g at 4 °C. The supernatant was transferred to a new tube and used immediately. Purified TasA was added to the seed protein extract to a final concentration of 50 µM and incubated at 4 °C overnight. Then, the mixture was incubated with equilibrated resin at 4 °C overnight. Control resin that did not contain bait protein was added to the seed protein extract. Flow-through, washing and elution fractions were collected and subjected to trichloroacetic acid protein precipitation. Precipitated proteins were resuspended in 1x Laemmli sample buffer (Bio-Rad) and heated at 100 °C for 5 min. The proteins were separated via SDS–PAGE through 12% polyacrylamide gels.

### Mass spectrometry analysis of protein bands

The protein bands of interest were identified by high-performance liquid chromatography/electrospray ionization tandem mass spectrometry (HPLC-ESI-MS/MS). Briefly, the bands were cut out after electrophoresis, washed and destained. Subsequently, the disulfide bridges were reduced with DTT and cysteines were alkylated by iodoacetamide; in-gel trypsin digestion was performed to extract peptides from the protein samples. The entire process was carried out automatically using an automatic digester (DigestPro, Intavis Bioanalytical Instruments). The peptides were concentrated and desalted using a C18 ZORBAX 300SB-C18 capture column (Agilent Technologies; 5 × 0.3 mm, 5 µm particle diameter and 300 Å pore size) in a gradient of 98% H2O:2% acetonitrile (ACN)/0.1% formic acid (FA) at a flow rate of 20 μl min−1 for 6 min. The capture column was connected through a 6-port valve to a ZORBAX 300SB-C18 analytical column (Agilent Technologies; 150 × 0.075 mm, 3.5 µm particle diameter and 300 Å pore size). Elution of the samples from the capture column was performed by a gradient of 0.1% FA in water as a mobile phase A and 0.1% FA in 80% ACN/20% water as a mobile phase B. The LC system was coupled to a nanospray source (CaptiveSpray, Bruker Daltonics) and a three-dimensional (3D) ion trap mass spectrometer (amaZon speed ETD, Bruker Daltonics) operating in the positive mode at a capillary voltage of 1,500 V and a sweep range of m/z 300–1,500. Data-dependent acquisition was carried out in the automatic mode with sequential collection of the full scan MS spectra (m/z 300–1,400), followed by collection of MS spectra in tandem via collision-induced dissociation of the 8 most abundant ions. ProteinScape 3 software (Bruker Daltonics) coupled to Mascot 3.1 search engine (Matrix Science) was used for identification by matching the MS/MS data against the Cucurbit genomes databases48.

### Far-western blotting

The interaction between TasA and oleosin was examined by far-western blotting49. Purified proteins were separated by SDS–PAGE through a 15% acrylamide gel and electrophoretically transferred onto PVDF membranes using a Trans-Blot Turbo transfer system (Bio-Rad) and PVDF transfer packs (Bio-Rad). The proteins were renatured by incubation of the membrane in denaturing and renaturing buffer (100 mM NaCl, 20 mM Tris, pH 7.6, 0.5 mM EDTA, 10% glycerol, 0.1% Tween 20, 2% skim milk powder and 1 mM DTT) containing decreasing concentrations of guanidine HCl (6, 3, 1, 0.1 and 0 M) and the membrane was blocked with blocking buffer (TBS containing 0.1% Tween 20 and 5% skim milk). Then, the membrane was incubated in protein-binding buffer (100 mM NaCl, 20 mM Tris, pH 7.6, 0.5 mM EDTA, 10% glycerol, 0.1% Tween 20, 3% skim milk powder and 1 mM DTT) containing 10 μg of purified TasA at 4 °C overnight. Excess protein was removed, and bound TasA was subjected to immunodetection. The membrane was probed with an anti-TasA antibody (rabbit) at a 1:20,000 dilution in blocking buffer. A secondary anti-rabbit IgG antibody conjugated to horseradish peroxidase (Bio-Rad) was used at a 1:3,000 dilution in the same buffer. The membranes were developed using Pierce ECL western blotting substrate (Thermo Fisher).

### Gene expression analysis by RT–qPCR

Total RNA was extracted from the frozen tissues collected at various timepoints after the treatments by phenol/chloroform extraction followed by precipitation with 1.5 M LiCl50 and sodium acetate51. For quantitative PCR with reverse transcription (RT–qPCR) assays, the RNA concentration was adjusted to 100 ng µl−1. Then, 1 µg of DNA-free total RNA was reverse transcribed into complementary DNA using SuperScript III reverse transcriptase (Invitrogen) and random primers in a final reaction volume of 20 µl according to the manufacturer’s instructions. RT–qPCR was performed using an iCycler-iQ system and an iQ SYBR Green Supermix kit from Bio-Rad according to the manufacturer’s instructions. The RT–qPCR cycle included the following steps: 95 °C for 3 min, followed by PCR amplification using a 40-cycle amplification programme (95 °C for 20 s, 56 °C for 30 s and 72 °C for 30 s), followed by a third step at 95 °C for 30 s. To evaluate the melting curve, 40 additional cycles were performed for 15 s starting at 75 °C, with stepwise temperature increases of 0.5 °C per cycle. The act7 gene was used as a reference for data normalization. The target genes were amplified using the primers listed in Supplementary Table 2. Primer efficiency tests and confirmation of the specificity of the amplification reactions were performed as previously described52. Relative transcript abundance was estimated using the ΔΔ cycle threshold (Ct) method53. Transcriptional data of the target genes were normalized to the act7 gene and are shown as the fold change in the expression levels of the target genes in each experimental treatment compared to those in the control treatment. The relative expression ratios were calculated as the differences between the qPCR threshold cycles (Ct) of the target gene and the Ct of the act7 gene (ΔCt = Ctgene of interest – Ctact7). The fold-change values were calculated as 2−ΔΔCt, assuming that a single PCR cycle represents a two-fold difference in the template abundance54,55. RT–qPCR analyses were performed in technical triplicates using three independently isolated RNA samples (biological triplicates).

### Total transcriptome analysis

For RNA sequencing analysis, 100 bp single-end read libraries were prepared using a TruSeq stranded total RNA kit (Illumina). The libraries were sequenced using a NextSeq550 instrument (Illumina). Raw reads were preprocessed by NextSeq System Suite v.2.2.0. using specific NGS technology configuration parameters56. This preprocessing removes low-quality, ambiguous and low-complexity stretches, linkers, adapters, vector fragments and contaminated sequences, and preserves the longest informative parts of the reads. SeqTrimNext also discarded sequences shorter than 25 bp. Subsequently, clean reads of the BAM files were aligned and annotated by Bowtie257 using the Cucumis melo genome v.4.058 as the reference; these data were then sorted and indexed using SAMtools v.1.48411059. Uniquely localized reads were used to calculate the read number values for each gene by Sam2counts (https://github.com/vsbuffalo/sam2counts). Differentially expressed genes (DEGs) in the treatment samples were analysed by DEgenes Hunter60, which provides a combined P value calculated on the basis of the Fisher method61 using nominal P values provided by edgeR62 and DEseq263. This combined P value was adjusted by the Benjamini-Hochberg procedure48 (false discovery rate approach) and used to rank all obtained DEGs. For each gene, a combined P value < 0.05 and log2 fold change (FC) >1 or <−1 were considered significance thresholds. A heatmap and DEG clustering were generated using ComplexHeatmap in R Studio and Kobas 2.064. Only profiles with a P value < 0.05 were considered in the present study. DEGs annotated using the Cucumis melo genome were processed to identify the Gene Ontology functional categories using sma3s65 and TopGo software66. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and enrichment were estimated using Bioconductor packages (Bioconductor.org) GGplot2, ClusterProfiler, DOSE and EnrichPlot in R Studio. The data were deposited in the GEO database under the reference GSE175611.

### CLSM

Bacteria, ROS and OBs inside the seeds were visualized by CLSM. Seeds were collected at the specified hours after treatment and were transversally cut with a scalpel. A drop of glycerol was applied to the sections, which were placed into 1.5-mm-thick cover glasses (22 × 22 mm). The bacteria were stained during inoculum preparation by CellTracker CM-DiI dye or in situ by Hoechst solution. Oil bodies and ROS were stained with Nile red (1:1,000 v/v) or dihydrorhodamine 123 (1:500 v/v), respectively. The images were acquired using a Leica SP5 confocal microscope with a ×HCX IRAPO L 25.0×0.95 WATER objective. To image the purified OB suspension, a drop of Nile red and Fast Green FCF-stained preparation was applied onto a patch of polymerized 1% agarose on a glass slide, which was covered with a coverslip. The images were acquired using a Leica SP5 confocal microscope with HCX PL APO lambda blue 63.0×1.40 OIL UV objective. Image processing was performed using Leica Application Suite Advance Fluorescence v.2.7.3.9723 (LCS Lite, Leica Microsystems) and FIJI/ImageJ software. Laser settings, scan speed, photomultiplier detector gain and pinhole aperture were constant for all acquired image stacks for each experiment.

### TEM

For TEM analysis, 16 h after treatment, seeds were fixed directly using 2% paraformaldehyde, 2.5% glutaraldehyde and 0.2 M sucrose mixture in 0.1 M phosphate buffer (PB) overnight at 4 °C. After three washes in PB, the portions were excised from the micropylar endosperm region and postfixed in 1% osmium tetroxide solution in PB for 90 min at room temperature, followed by PB washes and 15 min of stepwise dehydration in an ethanol series (30%, 50%, 70%, 90% and 100% twice). Between the 50% and 70% steps, the samples were incubated in-bloc in 2% uranyl acetate solution in 50% ethanol at 4 °C overnight. After dehydration, the samples were gradually embedded in low-viscosity Spurr’s resin (resin:ethanol, 1:1, 4 h; resin:ethanol, 3:1, 4 h; and pure resin overnight). The sample blocks were embedded in capsule moulds containing pure resin for 72 h at 70 °C. To image purified OBs in suspension or TasA protein under various polymerization conditions, carbon-coated copper grids were deposited over the sample drops and incubated at room temperature for 2 h. After incubation, the grids were washed in phosphate-buffered saline for 5 min, negatively stained with uranyl acetate (1%) for 20 s and washed once with water for 30 s.

For immunolabelling assays, carbon-coated copper grids were deposited over the samples of the corresponding OB treatments. After 2 h of incubation, the grids were washed in phosphate-buffered saline for 5 min and blocked with Pierce protein-free (TBS) blocking buffer (Thermo Fisher) for 30 min. An anti-TasA primary antibody was used at a 1:150 dilution in blocking buffer, and grids were deposited over the drops of antibody solution and incubated for 1 h at room temperature. The samples were washed three times with TBS-T (50 mM Tris-HCl, 150 mM NaCl, pH 7.5, and 0.1% Tween 20) for 5 min and then exposed to a 10-nm-diameter immunogold-conjugated secondary antibody (20 nm goat anti-rabbit conjugate, BBI solutions) for 1 h at a 1:50 dilution. The samples were then washed twice with TBS-T and once with water for 5 min each time. Finally, the grids were treated with 2% glutaraldehyde for 10 min, washed in water for 5 min, negatively stained with 1% uranyl acetate for 20 s and washed once with water for 30 s. The samples were left to dry and imaged under an FEI Tecnai G2 20 TWIN transmission electron microscope at an accelerating voltage of 80 kV. The images were acquired using TIA FEI Imaging Software v.4.14.

### Metabolite extraction from plant tissues

Metabolites from the samples of adult melon plants or radicles (MSV000084674, MSV000084278 and MSV000086360) were extracted by adding 1 ml methanol and vigorous vortexing. At least 2 or 3 replicates were done for each experiment, except for adult melon plants where only 1 replicate was done due to the high number of sections analysed for each plant. Due to the large sample number, weighting and normalization by biomass was not feasible. However, the size of plant sections was kept constant. To enable relative quantitative comparison between samples, we normalized peak areas to total ion count, which is a robust and widely applied strategy in non-targeted metabolomics. In addition, sulfamethazine (1 µM) was added to all samples as an internal standard and was used for normalization. After incubation for 2 h, the extracts were centrifuged at 14,000 g for 30 min. The acquired supernatants were stored at −80 °C until use in liquid chromatography–tandem mass spectrometry (LC–MS/MS).

Metabolites from seeds and seedlings (MSV000088139) were extracted after two rounds of methanol incubation with the powder obtained after grinding the tissue in liquid nitrogen and freeze-drying the samples. Three replicates were done for each condition. The samples were transferred to weighted vials and the methanol was evaporated at room temperature. Extracts were then resuspended in 1 ml methanol and centrifuged at 14,000 g for 20 min at 4 °C. Vials were weighed after the extractions to adjust the final concentrations to 1 mg ml−1, thus no internal standard was used.

In both cases, extraction controls were performed by adding 3 samples with pure methanol used for the extraction as a blank. Differences in abundance of metabolites were estimated on the basis of relative comparisons between equally treated samples due to the absence of authentic standards in the non-targeted metabolomics approach.

### LC–MS/MS

Non-targeted LC–MS/MS analysis was performed using a Q-Exactive Quadrupole-Orbitrap mass spectrometer coupled to a Vanquish ultra high-performance liquid chromatography (UHPLC) system (Thermo Fisher) (MSV000086360, MSV000084674 and MSV000084278) or a Shimadzu Nexera X2 UHPLC system, with attached photodiode array detector coupled to a Shimadzu 9030 QTOF mass spectrometer (MSV000088139).

For the Q-Exactive instrument, 5 µl of the samples were injected for separation on a UHPLC with a C18 core-shell column (Kinetex, 50 × 2 mm, 1.8 μm particle size, 100 Å pore size, Phenomenex). For the mobile phase, we used a flow rate of 0.5 ml min−1 (solvent A: H2O + 0.1% FA; solvent B: acetonitrile + 0.1% FA). During the chromatographic separation, we applied a linear gradient from 0–0.5 min 5% B, 0.5–4 min 5–50% B, 4–5 min 50–99% B, followed by a 2 min washout phase at 99% B and a 2 min re-equilibration phase at 5% B. For positive mode MS/MS acquisition, the electrospray ionization was set to a 35 a.u. sheath gas flow, 10 a.u. auxiliary gas flow, 2 a.u. sweep gas flow and 400 °C auxiliary gas temperature. The spray voltage was set to 3.5 kV with an inlet capillary of 250 °C. The S-lens voltage was set to 50 V. MS/MS product ion spectra were acquired in data-dependent acquisition mode. MS1 survey scans (150–1,500 m/z) and up to 5 MS/MS scans per data-dependent acquisition duty cycle were measured with a resolution of 17,500. The C-trap fill time was set to a maximum of 100 ms or until the Automatic Gain Control target of 5 × 105 ions was reached. The quadrupole precursor selection width was set to 1 m/z. Normalized collision energy was applied stepwise at 20, 30 and 40%, with z = 1 as default charge state. MS/MS scans were triggered with apex mode within 2–15 s from their first occurrence in a survey scan. Dynamic precursor exclusion was set to 5 s. Precursor ions with unassigned charge states and isotope peaks were excluded from MS/MS acquisition.

For the QTOF instrument, 2 μl were injected into a Waters Acquity Peptide BEH C18 column (1.7 μm, 300 Å, 2.1 × 100 mm). The column was maintained at 40 °C and run at a flow rate of 0.5 ml min−1, using 0.1% formic acid in H2O as solvent A and 0.1% formic acid in acetonitrile as solvent B. A gradient was employed for chromatographic separation starting at 5% B for 1 min, then 5% to 85% B for 9 min, 85% to 100% B for 1 min, and finally held at 100% B for 4 min. The column was re-equilibrated to 5% B for 3 min before the next run was started. The LC flow was switched to the waste for the first 0.5 min, then to the MS for 13.5 min, then back to the waste to the end of the run. The photodiode array detector acquisition was performed in the range 200–400 nm at 4.2 Hz, with 1.2 nm slit width. The flow cell was maintained at 40 °C. The MS system was tuned using standard NaI solution (Shimadzu). The same solution was used to calibrate the system before starting. System suitability was checked by including a standard sample made of 5 μg ml−1 thiostrepton, which was analysed regularly in between batches of samples.

All samples were analysed in positive polarity, using data-dependent acquisition mode. In this regard, full scan MS spectra (m/z 400–4,000, scan rate 20 Hz) were followed by 3 data-dependent MS/MS spectra (m/z 400–4,000, scan rate 20 Hz) for the 3 most intense ions per scan. The ions were selected when they reached an intensity threshold of 1,000, isolated at the tuning file Q1 resolution, fragmented using collision-induced dissociation with collision energy ramp (CE 10–40 eV) and excluded for 0.05 s (1 MS scan) before being reselected for fragmentation. The parameters used for the ESI source were: interface voltage 4 kV, interface temperature 300 °C, nebulizing gas flow 3 l min−1 and drying gas flow 10 l min−1.

In both cases, a mixture of 10 mg ml−1 each of sulfamethazine, sulfamethizole, sulfachloropyridazine, sulfadimethoxine, amitriptyline and coumarin was run after every 96 injections for quality control.

### Preprocessing for data analysis and MS/MS network analysis

Raw spectra were converted to .mzXML files using MSconvert (ProteoWizard). Feature identification was performed by MZmine open source software67,68 version 2.37 using the settings listed in Supplementary Table 3; this generated the .mgf files and quantification tables for statistical analyses via Metaboanalyst69 and feature-based molecular network in the Global Natural Product Social Molecular Networking (GNPS) environment70 using Ion Identity Network71.

### Feature-based molecular networking and spectral library search

Molecular networking and library search were performed using the GNPS environment with the settings listed in Supplementary Table 4 and visualized using Cytoscape software72. Putative annotation of detected features was obtained using automatic library search through the GNPS environment73, network annotation propagation74, chemical classification using ClassyFire75 and MolNetEnahancer workflows76.

Mirror plots were done using GNPS and https://metabolomics-usi.ucsd.edu/77, comparing mzspec of the selected features and the metabolites recorded in MS/MS databases (Extended Data Fig. 10). Annotations were done according to guidelines in ref. 78 (Supplementary Table 8), and levels are indicated in each figure.

The automatic workflow for the analysis of the radicles from seeds bacterized with Bacillus subtilis subsp. subtilis NCBI 3610 and mutant strains dataset (MSV000084674) can be accessed at: https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=b2e3636335d24c2da72634b7c6c8b63f; feature-based molecular networking: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=811bd911de4f48a1b59510215559afa2; network annotation propagation: https://proteomics2.ucsd.edu/ProteoSAFe/status.jsp?task=5d2297ea755b4d3986cd92d8e297e603; and MolNetEnhancer: https://proteomics2.ucsd.edu/ProteoSAFe/status.jsp?task=45823b975a8547e4b76b82ce8e037142.

Automatic workflow for the analysis of the radicles from the seeds treated with fengycin and TasA (MSV000084278) can be accessed at: https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=fc05e1b986cb4073b56d9f9a044f0a2f; feature-based molecular networking: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=446bcbde8940460b89eda8c0bddfdd27; network annotation propagation: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=eb776d82cf5f438282514e5d438c851e; and MolNetEnhancer: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=1f316080df4749a3ae00f2d644c95114.

Automatic workflow for the analysis of adult plants grown from control seeds or seeds bacterized with Bacillus subtilis subsp. subtilis NCBI 3610 (MSV000086360) can be accessed at the following addresses: https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=1090c82889a04e4ba9be1058a8a5d6db; feature-based molecular networking: https://proteomics2.ucsd.edu/ProteoSAFe/status.jsp?task=63d6d0c95b3e4d8d846d376499bb544e; network annotation propagation: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=914293e2bc3d4492bb9e0de3c7c49b2c; and MolNetEnhancer: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=e72d0b1205444e77880befcc9d23d5a8.

Automatic workflow for the analysis of the seeds and seedlings after seed treatments with bacteria, fengycin or TasA (MSV000088139) at 0, 24 and 120 h after treatment can be accessed at: https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=72bb762b7b8d45da8d3c5c9005a425d8; feature-based molecular networking: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=c6cce7f2ec2048febbc6b53b09e69525;

### Proportionality score

The proportionality score was calculated between two directly connected nodes across the entire molecular network using the following equation:

$${\rm{Proportionality}} = log\left( {\frac{{{\rm{Ns1}}/{\rm{Ms1}}}}{{{\rm{Ns2}}/{\rm{Ms2}}}}} \right),$$

where Ns1 and Ms1 correspond to the peak area of the detected features N and M in sample S1, while Ns2 and Ms2 correspond to the peak area of the detected features N and M in sample S2. A constant (k = 1.0 × 10−10) is added to each value to avoid including any zero values during the calculation. Chemical proportionality scores are available within the GNPS environment and can be accessed at: https://proteomics2.ucsd.edu/ProteoSAFe/status.jsp?task=7e8109b3ece945979f8b26ab7a48d615 (chemical proportionality table calculated between non-treated seeds at 0 and 24 h); https://proteomics2.ucsd.edu/ProteoSAFe/status.jsp?task=488cd1ca8e7540bb8033f79f0712fa3d (chemical proportionality table calculated between seeds treated with fengycin at 0 and 24 h); and https://proteomics2.ucsd.edu/ProteoSAFe/status.jsp?task=fca204a7579448bc9f821d08863fc0f9 (chemical proportionality table calculated between seeds treated with TasA at 0 and 24 h).

### Collection of plant material for 3D mass spectral molecular cartography

Individual parts of melon plants grown from control and bacterized seeds were sampled. Plants were sectioned in the same number of parts to make results comparative. For this, roots, stems and leaves were numbered, and coordinates were taken for later 3D representation and visualization. Approximately 200 mg of fresh individual plant parts was collected in 1.5 ml microcentrifuge tubes and flash frozen under liquid nitrogen. The samples were stored at −80 °C until further analysis.

### 3D modelling and visualization

Melon plants from non-bacterized seeds were used as the model plant for 3D representation of metabolomic features due to non-significant changes in morphology between plants. Approximately 250 images of melon plants taken from different angles and perspectives were converted into high-definition 3D meshes using Autodesk ReCap. These 3D meshes were then exported to the .obj format and edited using Meshmixer. Point coordinates of sampled plant parts were added using Meshlab79 and the 3D models were exported to the .stl format. Representative samples of roots, stems and leaves were obtained. Point coordinates and LC–MS data were combined into a .csv file and mapped on the 3D models using Ili80.

### Statistics and reproducibility

Statistical significance was assessed by appropriate tests (see figure legends). All analyses were performed using GraphPad Prism version 6. P values < 0.05 were considered significant. Asterisks indicate the level of statistical significance: *P < 0.05, **P < 0.01, ***P < 0.001 and ****P < 0.0001. Statistical analysis of metabolomic datasets was performed by Metaboanalyst v.5.0 (ref. 69), after data filtering by interquartile range (IQR). When needed, data were normalized to the feature corresponding to the internal standard of the samples. Experiments were repeated at least three times independently with similar results.

### Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

The RNA-seq data are deposited in the GEO database under the reference GSE175611. Metabolomics data are deposited at https://massive.ucsd.edu/ with the identifiers MSV000084674 (data from radicles from seeds bacterized with Bacillus subtilis subsp. subtilis NCBI 3610 and mutant strains), MSV000084278 (data from radicles from seeds treated with fengycin and TasA), MSV000086360 (data from adult plants grown from control seeds or seeds bacterized with Bacillus subtilis subsp. subtilis NCBI 3610) and MSV000088139 (data from seeds and seedlings after seed treatments with bacteria, fengycin or TasA). Source data are provided with this paper.

## References

1. Finch-Savage, W. E. & Leubner-Metzger, G. Seed dormancy and the control of germination. New Phytol. 171, 501–523 (2006).

2. O’Callaghan, M. Microbial inoculation of seed for improved crop performance: issues and opportunities. Appl. Microbiol. Biotechnol. 100, 5729–5746 (2016).

3. Blake, C., Christensen, M. N. & Kovacs, A. T. Molecular aspects of plant growth promotion and protection by Bacillus subtilis. Mol. Plant Microbe Interact. 34, 15–25 (2021).

4. Kovács, Á. T. Bacillus subtilis. Trends Microbiol. 27, 724–725 (2019).

5. Arnaouteli, S., Bamford, N. C., Stanley-Wall, N. R. & Kovács, Á. T. Bacillus subtilis biofilm formation and social interactions. Nat. Rev. Microbiol. 19, 600–614 (2021).

6. Cámara-Almirón, J. et al. Dual functionality of the amyloid protein TasA in Bacillus physiology and fitness on the phylloplane. Nat. Commun. 11, 1859 (2020).

7. Ongena, M. et al. Surfactin and fengycin lipopeptides of Bacillus subtilis as elicitors of induced systemic resistance in plants. Environ. Microbiol. 9, 1084–1090 (2007).

8. Rajjou, L. et al. Seed germination and vigor. Annu. Rev. Plant Biol. 63, 507–533 (2012).

9. Finkelstein, R., Reeves, W., Ariizumi, T. & Steber, C. Molecular aspects of seed dormancy. Annu. Rev. Plant Biol. 59, 387–415 (2008).

10. Sung, D. Y., Vierling, E. & Guy, C. L. Comprehensive expression profile analysis of the Arabidopsis hsp70 gene family. Plant Physiol. 126, 789–800 (2001).

11. Soeda, Y. et al. Gene expression programs during Brassica oleracea seed maturation, osmopriming, and germination are indicators of progression of the germination process and the stress tolerance level. Plant Physiol. 137, 354–368 (2005).

12. Rosental, L., Nonogaki, H. & Fait, A. Activation and regulation of primary metabolism during seed germination. Seed Sci. Res. 24, 1–15 (2014).

13. Mhlongo, M. I., Piater, L. A., Steenkamp, P. A., Labuschagne, N. & Dubery, I. A. Metabolic profiling of PGPR-treated tomato plants reveal priming-related adaptations of secondary metabolites and aromatic amino acids. Metabolites 10, 210 (2020).

14. Maier, B. A. et al. A general non-self response as part of plant immunity. Nat. Plants 7, 696–705 (2021).

15. Rodríguez, C. E., Mitter, B., Barret, M., Sessitsch, A. & Compant, S. Commentary: seed bacterial inhabitants and their routes of colonization. Plant Soil 422, 129–134 (2018).

16. Kinsinger, R. F., Shirk, M. C. & Fall, R. Rapid surface motility in Bacillus subtilis is dependent on extracellular surfactin and potassium ion. J. Bacteriol. 185, 5627–5631 (2003).

17. Calvio, C. et al. Swarming differentiation and swimming motility in Bacillus subtilis are controlled by swrA, a newly identified dicistronic operon. J. Bacteriol. 187, 5356–5366 (2005).

18. Grau, R. R. et al. A duo of potassium-responsive histidine kinases govern the multicellular destiny of Bacillus subtilis. mBio 6, e00581–15 (2015).

19. Molina-Santiago, C. et al. The extracellular matrix protects Bacillus subtilis colonies from Pseudomonas invasion and modulates plant co-colonization. Nat. Commun. 10, 1919 (2019).

20. Chen, Y. et al. Biocontrol of tomato wilt disease by Bacillus subtilis isolates from natural environments depends on conserved genes mediating biofilm formation. Environ. Microbiol. 15, 848–864 (2013).

21. Zeriouh, H., de Vicente, A., Pérez-García, A. & Romero, D. Surfactin triggers biofilm formation of Bacillus subtilis in melon phylloplane and contributes to the biocontrol activity. Environ. Microbiol. 16, 2196–2211 (2014).

22. Romero, D., Aguilar, C., Losick, R. & Kolter, R. Amyloid fibers provide structural integrity to Bacillus subtilis biofilms. Proc. Natl Acad. Sci. USA 107, 2230–2234 (2010).

23. Steinberg, N. et al. The extracellular matrix protein TasA is a developmental cue that maintains a motile subpopulation within Bacillus subtilis biofilms. Sci. Signal. 13, eaaw8905 (2020).

24. Noctor, G. et al. Glutathione in plants: an integrated overview. Plant Cell Environ. 35, 454–484 (2012).

25. Queval, G. et al. H2O2-activated up-regulation of glutathione in Arabidopsis involves induction of genes encoding enzymes involved in cysteine synthesis in the chloroplast. Mol. Plant 2, 344–356 (2009).

26. Shimada, T. L., Hayashi, M. & Hara-Nishimura, I. Membrane dynamics and multiple functions of oil bodies in seeds and leaves. Plant Physiol. 176, 199–207 (2018).

27. Graham, I. A. Seed storage oil mobilization. Annu. Rev. Plant Biol. 59, 115–142 (2008).

28. Farzand, A. et al. Suppression of Sclerotinia sclerotiorum by the induction of systemic resistance and regulation of antioxidant pathways in tomato using fengycin produced by Bacillus amyloliquefaciens FZB42. Biomolecules 9, 613 (2019).

29. Miquel, M. et al. Specialization of oleosins in oil body dynamics during seed development in Arabidopsis seeds. Plant Physiol. 164, 1866–1878 (2014).

30. Shao, Q., Liu, X., Su, T., Ma, C. & Wang, P. New insights into the role of seed oil body proteins in metabolism and plant development. Front. Plant Sci. 10, 1568 (2019).

31. Patel, H., Tscheka, C., Edwards, K., Karlsson, G. & Heerklotz, H. All-or-none membrane permeabilization by fengycin-type lipopeptides from Bacillus subtilis QST713. Biochim. Biophys. Acta 1808, 2000–2008 (2011).

32. Sur, S., Romo, T. D. & Grossfield, A. Selectivity and mechanism of fengycin, an antimicrobial lipopeptide, from molecular dynamics. J. Phys. Chem. B 122, 2219–2226 (2018).

33. Tzen, J. T. C., Cao, Y. Z., Laurent, P., Ratnayake, C. & Huang, A. H. C. Lipids, proteins, and structure of seed oil bodies from diverse species. Plant Physiol. 101, 267–276 (1993).

34. Petras, D. et al. Chemical proportionality within molecular networks. Anal. Chem. 93, 12833–12839 (2021).

35. Jeevan Kumar, S. P., Rajendra Prasad, S., Banerjee, R. & Thammineni, C. Seed birth to death: dual functions of reactive oxygen species in seed physiology. Ann. Bot. 116, 663–668 (2015).

36. Bailly, C. The signalling role of ROS in the regulation of seed germination and dormancy. Biochem. J. 476, 3019–3032 (2019).

37. Bailly, C., El-Maarouf-Bouteau, H. & Corbineau, F. From intracellular signaling networks to cell death: the dual role of reactive oxygen species in seed physiology. C. R. Biol. 331, 806–814 (2008).

38. Ongena, M. & Jacques, P. Bacillus lipopeptides: versatile weapons for plant disease biocontrol. Trends Microbiol. 16, 115–125 (2008).

39. Stenzel, I. et al. Jasmonate biosynthesis and the allene oxide cyclase family of Arabidopsis thaliana. Plant Mol. Biol. 51, 895–911 (2003).

40. Gullner, G., Komives, T., Király, L. & Schröder, P. Glutathione S-transferase enzymes in plant-pathogen interactions. Front. Plant Sci. 9, 1836 (2018).

41. Cowan, A. K. Phospholipids as plant growth regulators. Plant Growth Regul. 48, 97–109 (2006).

42. Farag, K. M. & Palta, J. P. Use of lysophosphatidylethanolamine, a natural lipid, to retard tomato leaf and fruit senescence. Physiol. Plant. 87, 515–521 (1993).

43. Özgen, M., Park, S. & Palta, J. P. Mitigation of ethylene-promoted leaf senescence by a natural lipid, lysophosphatidylethanolamine. HortScience 40, 1166–1167 (2005).

44. Cowan, A. K. Plant growth promotion by 18:0-lyso-phosphatidylethanolamine involves senescence delay. Plant Signal. Behav. 4, 324–327 (2009).

45. Regente, M., Corti Monzón, G. & de la Canal, L. Phospholipids are present in extracellular fluids of imbibing sunflower seeds and are modulated by hormonal treatments. J. Exp. Bot. 59, 553–562 (2008).

46. Tzen, J. T. C. & Huang, A. H. C. Surface structure and properties of plant seed oil bodies. J. Cell Biol. 117, 327–335 (1992).

47. De Chirico, S., di Bari, V., Foster, T. & Gray, D. Enhancing the recovery of oilseed rape seed oil bodies (oleosomes) using bicarbonate-based soaking and grinding media. Food Chem. 241, 419–426 (2018).

48. Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300 (1995).

49. Wu, Y., Li, Q. & Chen, X. Z. Detecting protein-protein interactions by far western blotting. Nat. Protoc. 2, 3278–3284 (2007).

50. Lagrimini, L. M., Burkhart, W., Moyer, M. & Rothstein, S. Molecular cloning of complementary DNA encoding the lignin-forming peroxidase from tobacco: molecular analysis and tissue-specific expression. Proc. Natl Acad. Sci. USA 84, 7542–7546 (1987).

51. Viennois, E., Chen, F., Laroui, H., Baker, M. T. & Merlin, D. Dextran sodium sulfate inhibits the activities of both polymerase and reverse transcriptase: lithium chloride purification, a rapid and efficient technique to purify RNA. BMC Res. Notes 6, 360 (2013).

52. Vargas, P., Felipe, A., Michán, C. & Gallegos, M. T. Induction of Pseudomonas syringae pv. tomato DC3000 mexAB-oprM multidrug efflux pump by flavonoids is mediated by the repressor PmeR. Mol. Plant Microbe Interact. 24, 1207–1219 (2011).

53. Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 25, 402–408 (2001).

54. Pfaffl, M. W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29, e45 (2001).

55. Rotenberg, D., Thompson, T. S., German, T. L. & Willis, D. K. Methods for effective real-time RT-PCR analysis of virus-induced gene silencing. J. Virol. Methods 138, 49–59 (2006).

56. Falgueras, J. et al. SeqTrim: a high-throughput pipeline for pre-processing any type of sequence read. BMC Bioinformatics 11, 38 (2010).

57. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).

58. Ruggieri, V. et al. An improved assembly and annotation of the melon (Cucumis melo L.) reference genome. Sci. Rep. 8, 1–9 (2018).

59. Li, H. et al. The sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).

60. González Gayte, I., Bautista Moreno, R., Seoane Zonjic, P. & Claros, M. G. DEgenes Hunter - a flexible R pipeline for automated RNA-seq studies in organisms without reference genome. Genom. Comput. Biol. 3, e31 (2017).

61. Fisher, R. A. in Breakthroughs in Statistics: Methodology and Distribution (eds Kotz, S. & Johnson, N. L.) 66–70 (Springer, 1992).

62. Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2009).

63. Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol. 11, R106 (2010).

64. Xie, C. et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 39, W316–W322 (2011).

65. Casimiro-Soriguer, C. S., Muñoz-Mérida, A. & Pérez-Pulido, A. J. Sma3s: a universal tool for easy functional annotation of proteomes and transcriptomes. Proteomics 17, 12 (2017).

66. Alexa, A. & Rahnenführer, J. Gene Set Enrichment Analysis with topGO http://compdiag.molgen.mpg.de/ngfn/docs/2007/sep/topGO_Exercises.pdf (2007).

67. Pluskal, T., Castillo, S., Villar-Briones, A. & Orešič, M. MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics 11, 395 (2010).

68. Olivon, F., Grelier, G., Roussi, F., Litaudon, M. & Touboul, D. MZmine 2 data-preprocessing to enhance molecular networking reliability. Anal. Chem. 89, 7836–7840 (2017).

69. Chong, J., Wishart, D. S. & Xia, J. Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr. Protoc. Bioinformatics 68, e86 (2019).

70. Nothias, L. F. et al. Feature-based molecular networking in the GNPS analysis environment. Nat. Methods 17, 905–908 (2020).

71. Schmid, R. et al. Ion identity molecular networking in the GNPS environment. Nat. Commun. 12, 3832 (2021).

72. Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003).

73. Wang, M. et al. Sharing and community curation of mass spectrometry data with Global Natural Products Social Molecular Networking. Nat. Biotechnol. 34, 828–837 (2016).

74. da Silva, R. R. et al. Propagating annotations of molecular networks using in silico fragmentation. PLoS Comput. Biol. 14, e1006089 (2018).

75. Djoumbou Feunang, Y. et al. ClassyFire: automated chemical classification with a comprehensive, computable taxonomy. J. Cheminform. 8, 61 (2016).

76. Ernst, M. et al. Molnetenhancer: enhanced molecular networks by integrating metabolome mining and annotation tools. Metabolites 9, 144 (2019).

77. Bittremieux, W. et al. Universal MS/MS visualization and retrieval with the metabolomics spectrum resolver web service. Preprint at bioRxiv https://doi.org/10.1101/2020.05.09.086066 (2020).

78. Sumner, L. W. et al. Proposed minimum reporting standards for chemical analysis. Metabolomics 3, 211–221 (2007).

79. Cignoni, P. et al. MeshLab: an open-source mesh processing tool. In Proc. 6th Eurographics Italian Chapter Conference (eds Scarano, V. et al.) 129–139 (Eurographics Association, 2008).

80. Floros, D. J. et al. Mass spectrometry based molecular 3D-cartography of plant metabolites. Front. Plant Sci. 8, 429 (2017).

## Acknowledgements

We thank S. M. Rojas for technical support, J. G. Maldonado from the Ultrasequencing Unit of the SCBI-UMA for RNA sequencing, C. Cárdenas for protein analysis (Proteomic Unit, SCAI-UMA), J. F. L. Téllez for technical support on TEM analysis, J. Pearson for technical support on confocal laser scanning microscopy analysis (Bionand), S. Elsayed (Leiden University) for technical support on metabolomics analysis and T. Ketelaar (Wageningen University and Research) for suggestions on the interpretation of TEM images. This work was supported by grants from ERC Starting Grant (BacBio 637971) Plan Nacional de I+D+i of Ministerio de Economía y Competitividad and Ministerio de Ciencia e Innovación (AGL2016-78662-R and PID2019-107724GB-I00), Proyectos I+D+I en el marco del Programa Operativo FEDER Andalucía (UMA18-FEDERJA-055) and Ayudas a Proyectos I+D+I Programa PAIDI 2020 of Junta de Andalucía (P20_00479). M.V.B.-C. received an FPU contract (FPU17/03874) from the Ministerio de Ciencia, Innovación y Universidades. C.M.-S. was funded by the programme Juan de la Cierva Incorporación (IJC2018-036923-I). D.P. was supported by the German Research Foundation (DFG) with Grant PE 2600/1. A.M.C.-R. and P.C.D. were supported by the National Institutes of Health (NIH) grant DP2GM137413. P.C.D. was supported by the Gordon and Betty Moore Foundation through Grant GBMF7622, and the U.S. National Institutes of Health for the Center (P41 GM103484, R03 CA211211, R01 GM107550).

## Author information

Authors

### Contributions

D.R. conceived and designed the work, drafted and edited the text; M.V.B.-C. collected most of the experimental data, and drafted the manuscript; C.M.-S. designed, collected and analysed MS data, and edited the manuscript; A.M.C.-R. collected and analysed MS data, and edited the text; D.P. collected MS data and edited the text; L.D.-M. informatically analysed data and drafted figures; V.J.C. supervised MS analyses and edited the text; A.d.V. substantially revised and edited the text; A.P.-G. substantially revised and edited the text; P.C.D. substantially revised and edited the text.

### Corresponding author

Correspondence to D. Romero.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review

### Peer review information

Nature Microbiology thanks Ruifu Zhang and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

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

## Extended data

### Extended Data Fig. 1 Seed treatment with B. subtilis cells does not induce an increase in the germination rate but accelerates post germination progression.

a Germination rates of untreated and bacterized seeds represented as the percentage of germinated seeds versus time. b Relative expression levels of the GA20ox1 and CYP707A1 genes in bacterized seeds 0 and 16 hours after the treatment compared with those in the untreated seeds. The average values of two biological replicates with three technical replicates are shown with error bars representing SD. c Heatmap of fold changes of selected DEGs identified by total transcriptomic analysis of bacterized and control seeds (untreated) 16 hours after the treatment. Colour scale indicates the Log2(FC) of each DEG. d KEGG pathway enrichment analysis of all DEGs identified by RNA-seq of bacterized and control seeds 16 hours after the treatment.

### Extended Data Fig. 2 Carboxylic acids and lipids represent the major metabolic changes in bacterized seeds.

a Top 20 features with high median weighted sum of absolute regression coefficient scores determined by PLS-DA and calculated using MetaboAnalyst, selected as features discriminating bacterized seeds from control seeds 0 hours after the treatment (left) or 24 hours after the treatment (right). Feature IDs are accompanied by chemical class according to Classyfire classification. Features in bold are selected as features discriminating any growth-promoting treatment from the control treatment. b Volcano plots showing the metabolites significatively increased or decreased in leaves (left) and stem (right) of plants emerged from bacterized seeds respecting control plants. Colours indicate the chemical class of each metabolite according to Classyfire. c Distribution of features identified by spectral match as L-tryptophan and cinnamic acid in 3D models representing the plants grown from the control or B. subtilis-treated seeds.

### Extended Data Fig. 3 B. subtilis cells enter the seeds and colonize the inner regions.

a Scanning electron microscopy micrographs of bacterized seeds 5 days after the treatment shows B. subtilis cells colonizing the inner side of the seed coat. Scale bars: 100 = µm (left image) and 30 µm (right image); r: radicle, i: inner seed coat, o- outer seed coat. b Confocal microscopy images of transversally cut bacterized seeds 16 hours after the treatment with fluorescently labeled B. subtilis cells. Scale bars: 100 µm (upper panels) or 10 µm (lower panels); e: embryo. c CFU counts of the surface extracts of the seeds bacterized with the WT (3610) strain and ECM mutants 1 hour after the treatment. Average values are shown with error bars representing SD (n = 6). d Radicle colonization of the ECM mutants relative to that of the WT assumed as 100% (discontinued line) in radicle extracts five days after seed treatment. Average values are shown. Error bars represent SEM. Statistical significance was assessed by one-way ANOVA with post hoc Dunnett’s multiple comparisons test (each treatment vs. WT treatment except for ∆srf and ∆pks) (n = 3 in WT, n = 4 in all ECM mutants except for ∆srf and ∆pks, with n = 2; ∆eps *P = 0.0391, ∆bslA ***P = 0.0007, ∆tasA **P = 0.0014, JC81 **P = 0.0091, ∆tapA ***P = 0.0009).

### Extended Data Fig. 4 Components of the extracellular matrix of B. subtilis trigger metabolic reprogramming of the seed radicles related to plant growth stimulation.

a Heatmap of the hierarchical clustering of the top 50 features of impacted molecular families in the radicles from bacterized seeds. The color code inside the heatmap depicts the relative fold change of each metabolite between groups. Color code accompanying feature names indicates their chemical class according to Classyfire. b Network analysis of representative features related to the presence of bacteria or ECM. Normalized abundances in the radicles in seeds subjected to various treatments are represented in features of all groups. Average values are shown with error bars representing SD. Statistical significance was assessed by one-way ANOVA with post hoc Dunnett’s multiple comparisons test (n = 6 except in control, where n = 12; p < 0.0001). Created with BioRender.com.

### Extended Data Fig. 5 Aggregated forms of TasA are preferred under the apoplast conditions and have the highest radicle growth-promoting activity after seed treatments.

a Left: Polimeryzation dynamics of TasA in the fluorescence emission assays using ThT. Discontinued lines indicate time points when fractions were used for seed treatments. Right: Percentage of radicle growth promotion from the seeds treated with TasA at various polymerization states normalized to the radicle area of the control samples (100%, discontinued line) five days after seed treatment. Average values are shown. Error bars represent SD. Statistical significance was assessed by after one-way ANOVA with post hoc Dunnett’s multiple comparisons test (each treatment vs. control treatment) (n = 11, 12, 10, 11 and 9 for control, n = 9, 10, 9, 12 and 9 for TasA 1.5 µM and n = 11, 10, 9, 10 and 10 for TasA 3 µM, referring to 0, 4, 24, 48 and 120 hours respectively). b Left: Fluorescence emission of ThT over time in the TasA polymerization assay in polymerization buffer (buffer) or in apoplastic fluid (AF). Right: Representative transmission electron microscopy micrographs of TasA after 16 hours of polymerization in buffer or AF.

### Extended Data Fig. 6 Treatments of the seeds with TasA and fengycin influence the final metabolic patterns of emerged radicles.

a PCA 3D score plot of the metabolome of the radicles showing clustering of the samples based on seed treatment. The percentage of variation explained by each principal component is indicated on the axes. b Heatmap of hierarchical clustering results of the top 50 features impacted in the radicles grown from TasA- and fengycin-treated seeds. The color code inside the heatmap depicts the relative fold change of each metabolite between groups. c Top 20 features with high median weighted sum of absolute regression coefficient scores determined by PLS-DA and calculated using MetaboAnalyst. Feature IDs are accompanied by chemical class according to Classyfire classification. Features in bold are selected as features discriminating any growth-promoting treatment from the control treatment. d, e Molecular family of features discriminating fengycin-treated seeds (D) or TasA-treated seeds (E) from the control seeds according to PLS-DA. The chemical structures of annotated features and their average mass based on spectral matches to GNPS libraries are also represented for the corresponding molecular families. Pie charts indicate the peak abundance of each metabolite in the corresponding condition. Node shape indicates the level of identification according to Sumner et al., 78.

### Extended Data Fig. 7 The size of OBs in the seeds is altered after seed treatment with B. subtilis cells or with purified TasA and fengycin.

a Representative CLSM images of transversally cut seeds 16 hours after the treatment with B. subtilis. Oil bodies were stained with Nile red, and bacteria were stained with Hoechst. Scale bar: 20 µm. b Transmission electron microscopy images of thin sections of the seeds (top) or purified OB suspensions (bottom) 16 hours after the treatment with water (control), fengycin, or TasA. Arrows indicate oil bodies. c Representative transmission electron microscopy images of thin sections of the seeds 16 hours after the treatment with water (control), fengycin or TasA. G: glyoxysomes; arrows: oil bodies. d Transmission electron microscopy images of purified OB suspension 16 hours after the treatment with TasA. Samples were immunolabeled with a TasA antibody (1:150) and secondary anti-rabbit antibodies conjugated to 10 nm gold particles (1:50).

### Extended Data Fig. 8 TasA and fengycin induce changes in OBs morphology and glycerophospholipids dynamics in seeds.

a Representative CLSM images of transversally cut seeds 16 hours after the treatment with 10 µM fengycin or 3 µM TasA. Oil bodies were stained with Nile red (scale bar: 20 µm). b Chemical proportionality analysis of a glycerophospholipids molecular family in non-treated seeds (control) and seeds treated with fengycin or TasA from 0 h to 24 h after treatment. The size of the nodes is directly related with the abundance of the features at time 0 h. Arrows indicate chemical directionality of the modifications found and the color scale of the arrow represents the ChemProp score, the value obtained after the measurement of the peak area changes of connected nodes in a molecular network across a sequential data frame by comparing their proportions.

### Extended Data Fig. 9 Treatment of the seeds with fengycin increase the expression and accumulation of metabolites related to glutathione metabolism and flavonoids.

a. KEGG pathway enrichment analysis of all DEGs identified by total transcriptomic analysis of the leaves of adult plants grown from control or fengycin-treated seeds 48 hours after inoculation with B. cinerea. b 3D plant models of adult plants grown from control (left) or B. subtilis-treated seeds (right) showing global distribution of the feature with m/z identical to discriminant feature 1785 of the GSH molecular family identified in the metabolome of the radicles grown from the seeds treated with purified TasA or fengycin. c Molecular family with features differentially abundant in the aerial regions of seedlings from fengycin-treated seeds vs. control seedlings 5 days after seed treatment, annotated as flavonoids according to Classyfire. The chemical structures of annotated features and their average mass based on spectral matches to GNPS libraries are also represented for the corresponding molecular families. Pie charts indicate the peak abundance of each metabolite in the corresponding condition. Node shape indicates the level of identification according to Sumner et al., 78.

### Extended Data Fig. 10 Mirror plots comparing spectra from named features along the manuscript to standard spectra deposited in GNPS.

In the upper part of the plot (black lines) is represented the MS spectra of the candidate feature and in the lower part (green lines); the MS spectra of the standard compound. Mirror plots have been generated using https://metabolomics-usi.ucsd.edu/.

## Supplementary information

### Supplementary Information

Supplementary Tables 1–8, references and source data files (list only).

## Source data

### Source Data Fig. 1

Statistical source data.

### Source Data Fig. 2

Statistical source data.

### Source Data Fig. 3

Statistical source data.

### Source Data Fig. 4

Statistical source data.

### Source Data Fig. 4

Unprocessed gels and western blots.

### Source Data Fig. 5

Statistical source data.

### Source Data Extended Data Fig. 1

Statistical source data.

### Source Data Extended Data Fig. 3

Statistical source data.

### Source Data Extended Data Fig. 4

Statistical source data.

### Source Data Extended Data Fig. 5

Statistical source data.

## Rights and permissions

Reprints and Permissions

Berlanga-Clavero, M.V., Molina-Santiago, C., Caraballo-Rodríguez, A.M. et al. Bacillus subtilis biofilm matrix components target seed oil bodies to promote growth and anti-fungal resistance in melon. Nat Microbiol 7, 1001–1015 (2022). https://doi.org/10.1038/s41564-022-01134-8

• Accepted:

• Published:

• Issue Date:

• DOI: https://doi.org/10.1038/s41564-022-01134-8


text:
# Newton’s method with 10 lines of Python

I’m starting a new series of blog posts, called “XY in less than 10 lines of Python“. This first one is about Newton’s method, which is an old numerical approximation technique that could be used to find the roots of complex polynomials and any differentiable function.

Let’s say we have a complicated polynomial:

$f(x)=6x^5-5x^4-4x^3+3x^2$

and we want to find its roots. Unfortunately we know from the Galois theory that there is no formula for a 5th degree polynomial so we’ll have to use numeric methods. Wait a second. What’s going on? We all learned the quadratic formula in school, and there are formulas for cubic and quartic polynomials, but Galois proved that no such “root-finding” formula exist for fifth or higher degree polynomials, that uses only the usual algebraic operations (addition, subtraction, multiplication, division) and application of radicals (square roots, cube roots, etc).

Some nice guys pointed out on reddit that I didn’t quite get the theory right. Sorry about that, I’m no mathematician by any means. It turns out that this polynomial could be factored into $x^2(x-1)(6x^2 + x - 3)$ and solved with traditional cubic formula. Also the theorem I referred to is the Abel-Ruffini Theorem and it only applies to the solution to the general polynomial of degree five or greater. Nonetheless the example is still valid, and demonstrates how would you apply Newton’s method, to any polynomial, so let’s crack on.

So in these cases we have to resort to numeric linear approximation. A simple way would be to use the intermediate value theorem, which states that if $f(x)$ is continuous on $[a,b]$ and $f(a) < y< f(b)$, then there is an $x$ between $a$ and $b$ so that $f(x)=y$.

We could exploit this by looking for an $x_1$ where $f(x_1)>0$ and an $x_2$ where $f(x_2)<0$, and then we could be certain that $f(x)=0$ must be between $x_1$ and $x_2$. So we could  check $x_3=\frac{x_2-x_1}{2}$, and find that $f(x_3)$ will be positive, then continue with$x_4=\frac{x_2-x_3}{2}$ which would be negative and our proposed $x_n$ would become closer and closer to $f(x_n)=0$. This method however can be pretty slow, so Newton devised a better way to speed things up (when it works).

If you look at the figure below, you’ll see the plot of our polynomial. It has three roots at 0, 1, and somewhere in-between. So how do we find these?

In Newton’s method we take a random point $f(x_0)$, then draw a tangent line through $x_0, f(x_0)$, using the derivative $f'(x_0)$. The point $x_1$ where this tangent line crosses the $x$ axis will become the next proposal we check. We calculate the tangent line at $f'(x_1)$ and find $x_2$. We carry on, and as we do $|f(x_n)| \to 0$, or in other words we can make our approximation as close to zero as we want, provided we are willing to continue with the number crunching. Using the fact that the slope of tangent (by the definition of derivatives) at $x_n$ is $f'(x_n)$ we can derive the formula for $x_{n+1}$, i.e. where the tangent crosses the $x$ axis:

$x_{n+1}=x_n-\frac{f(x_n)}{f'(x_n)}$

With all this in mind it’s easy to write an algorithm that approximates $f(x)=0$ with arbitrary error $\varepsilon$. Obviously in this example depending on where we start $x_0$ we might find different roots.

There you go, we’are done in 10 lines (9 without the blank line), even less without the print statements. So in order to use this, we need two functions, $f(x)$ and $f'(x)$. For this polynomial these are:

Now we can simply find the three roots with:

All of the above code, and some additional comparison test with the scipy.optimize.newton method can be found in this Gist. And don’t forget, if you find it too much trouble differentiating your functions, just use SymPy, I wrote about it here.

Newton’s method is pretty powerful but there could be problems with the speed of convergence, and awfully wrong initial guesses might make it not even converge ever, see here. Nonetheless I hope you found this relatively useful.. Let me know in the comments.

Newton’s method with 10 lines of Python
• Analytics2Go

Love it!

• Very cool. Congrats!

• Ericracia

Hey, How do you put code in your blog’s posts like this one? Is it a plug-in?

• Силантий

Why function dx returns the absolute value of the difference between zero and the value of the function? If you can just return the module value of the function. It’s the same thing.

• Filib Rinaldi

Great Job!!

• Kartik Mukkavilli

I would like to show the number of iterations done in this program. How would I do that?

• Bjørn Egil Johannessen

I think it was Nils Henrik Abel who proved that the 5th degree polynomial cannot be solved using the elemental calculus. He is a well-known Norwegian mathematician who unfortunately died at a young age (26 years).

format: plain