Abstract

Bacillus subtilis a Gram-positive soil-dwelling bacterium known for its wide range of bioactive secondary metabolites. The lipopeptide plipastatin produced by most B. subtilis isolates have been shown to exhibit potent antifungal activity against plant pathogenic fungi. While the effect of these antifungal compounds are well studied in the context of biocontrol, much less is known of their role in the environment, which also harbor nonproducing strains of these compounds. Fusarium species produce multiple antibacterial compounds resulting in dysbiosis of the plant-associated microbiome and inhibition of plant beneficial bacteria like B. subtilis. While plipastatin is expected to be important for survival of B. subtilis, not all isolates carry the biosynthetic gene cluster for plipastatin suggesting that the protective effect of plipastatin might be shared. In this study, we investigated the protective effect of plipastatin against Fusarium oxysporum in a coculture using a producer and a nonproducer isolate of plipastatin. We tested the survival of single and cocultured strains under Fusarium challenge in liquid media and solid agar plates to dissect the influence of spatial structure. Our results highlights that plipastatin protects the nonproducer strain in a density-dependent manner.

Introduction

Microbial interactions are the drivers of host health and disease, from individual microbe–microbe interactions to the whole network that make up the host microbiome. The many interactions that take place in the microbiome around and in a plant is one of the main factors contributing to host fitness (Mendes et al. 2013, Compant et al. 2019). Numerous bacteria have been shown to possess plant beneficial effects, by promotion of growth, nutrient utilization, or protection against plant pathogens either by induction of systemic resistance or direct inhibition of the plant disease-causing agent (Lugtenberg and Kamilova 2009, Raaijmakers et al. 2009, Blake et al. 2021).

Fusarium oxysporum is a filamentous fungus and causative agent of Fusarium wilt in several economically important crops including banana, tobacco, and tomato (Summerell 2019). The disease is cause of significant yield losses due to clogging the vascular vessels of the plant leading to wilting of the leaves and death of the plant. Among fungal plant pathogens, F. oxysporum is placed as the fifth most scientifically and economically important fungal pathogen due to the breadth of its host range and severe impact on several large sectors of the agricultural market (Dean et al. 2012). Historically, chemical fungicides have been employed in combating Fusarium wilt, which has led to widespread ecological consequences and risk of resistance developing in F. oxysporum (Bawa 2016, Hudson et al. 2021, El-Aswad et al. 2023) underlining the need for more environmentally friendly biocontrol agents.

Bacillus subtilis group species are known to produce a suite of bioactive secondary metabolites, including lipopeptides such as surfactin, fengycins, and iturin (Ongena and Jacques 2008), including plipastatin that belongs to fengycins (Honma et al. 2012, Gong et al. 2015), which together with surfactin elicits systemic resistance in plants (Ongena et al. 2007, Farace et al. 2015). Plipastatin variants exhibit species-specific inhibition of Fusarium, while some plipastatin was reported to only affect F. graminearum and not F. oxysporum (Zhao et al. 2014, Gong et al. 2015), others exhibit activity against both (Kiesewalter et al. 2021) making Bacillus a possible candidate for future biocontrol of the fungal pathogen. Fusarium produces several antibacterial compounds, which negatively affect the plant-associated microbial community (Sondergaard et al. 2016, Venkatesh and Keller 2019), leading to dysbiosis and death of plant beneficial bacteria underlining the importance of effective countermeasures. Therefore, it is also expected that production of Fusarium-inhibiting metabolites, like plipastatin, could benefit the bacterial community by combating the invading fungi.

Interestingly, plipastatin production capability is not present in all B. subtilis isolates, and even in the same soil samples some isolates display plipastatin production while others do not (Kiesewalter et al. 2021), which raises the question whether the protective effect of plipastatin is shared with the species and broadly the microbial community. While sharing of plipastatin could protect the community, it might also suffer from cheating by other nonproducing Bacilli, which could potentially out-compete the producer strain due to fitness costs leading to reduced level of producers and a loss of plipastatin protection in the community. The lack of plipastatin and fengycin production has been broadly proposed in wide variety of isolates of B. subtilis and Bacillus velezensis, respectively, as either frameshift mutations or partial deletion have been identified in the respective biosynthetic gene clusters of these secondary metabolites (Kiesewalter et al. 2021, Steinke et al. 2021).

In this study, we evaluate the role of plipastatin in a coculture system of two B. subtilis isolates when challenged by F. oxysporum. The two B. subtilis isolates originate from the same sample site, but one, MB9_B1, is a producer of plipastatin and thus able to inhibit Fusarium while the other isolate, MB9_B6, is unable to produce this secondary metabolite (Kiesewalter et al. 2021). In our coculture setup, we demonstrate that the plipastatin producer can provide protection to the nonproducer in a density-dependent manner.

Materials and methods

Strains, culturing, and genetic modification

Strains used in this study can be found in Table 1. Bacillus subtilis was routinely cultured in lysogeny broth (Lenox, Carl Roth) for overnight cultures while F. oxysporum was cultured in potato dextrose broth (PDB; Carl Roth). Strains were genetically engineered using the protocol below which is based on Kunst and Rapoport (1995).

Table 1.

Strains used in the study.

StrainCharacteristicsReference
B. subtilis MB9_B1Natural isolate, plipastatin producerThérien et al. (2020)
B. subtilis MB9_B6Natural isolate, plipastatin nonproducerKiesewalter et al. (2021)
B. subtilis MB9_B1 mKate2MB9_B1 transformed with gDNA from TB501 amyE::Phyperspank-mKate2 (specR) (Dragoš et al. 2018)This study
B. subtilis MB9_B1 mKate2 ΔppsCMB9_B1 mKate2 transformed with gDNA from DS4114 ΔppsC::tetR (Müller et al. 2014)This study
B. subtilis MB9_B6 GFPMB9_B6 transformed with gDNA from TB500 amyE::Phyperspank-gfp (specR) (Mhatre et al. 2017).This study
F. oxysporum
IBT 40 872
From surface of an equipment in a Danish factory, 2005, by Anne SvendsenIBT Culture Collection of Fungi, DTU
StrainCharacteristicsReference
B. subtilis MB9_B1Natural isolate, plipastatin producerThérien et al. (2020)
B. subtilis MB9_B6Natural isolate, plipastatin nonproducerKiesewalter et al. (2021)
B. subtilis MB9_B1 mKate2MB9_B1 transformed with gDNA from TB501 amyE::Phyperspank-mKate2 (specR) (Dragoš et al. 2018)This study
B. subtilis MB9_B1 mKate2 ΔppsCMB9_B1 mKate2 transformed with gDNA from DS4114 ΔppsC::tetR (Müller et al. 2014)This study
B. subtilis MB9_B6 GFPMB9_B6 transformed with gDNA from TB500 amyE::Phyperspank-gfp (specR) (Mhatre et al. 2017).This study
F. oxysporum
IBT 40 872
From surface of an equipment in a Danish factory, 2005, by Anne SvendsenIBT Culture Collection of Fungi, DTU

gfp: gfpmut2, specR: spectinomycin resistance, and tetR: tetracycline resistance.

Table 1.

Strains used in the study.

StrainCharacteristicsReference
B. subtilis MB9_B1Natural isolate, plipastatin producerThérien et al. (2020)
B. subtilis MB9_B6Natural isolate, plipastatin nonproducerKiesewalter et al. (2021)
B. subtilis MB9_B1 mKate2MB9_B1 transformed with gDNA from TB501 amyE::Phyperspank-mKate2 (specR) (Dragoš et al. 2018)This study
B. subtilis MB9_B1 mKate2 ΔppsCMB9_B1 mKate2 transformed with gDNA from DS4114 ΔppsC::tetR (Müller et al. 2014)This study
B. subtilis MB9_B6 GFPMB9_B6 transformed with gDNA from TB500 amyE::Phyperspank-gfp (specR) (Mhatre et al. 2017).This study
F. oxysporum
IBT 40 872
From surface of an equipment in a Danish factory, 2005, by Anne SvendsenIBT Culture Collection of Fungi, DTU
StrainCharacteristicsReference
B. subtilis MB9_B1Natural isolate, plipastatin producerThérien et al. (2020)
B. subtilis MB9_B6Natural isolate, plipastatin nonproducerKiesewalter et al. (2021)
B. subtilis MB9_B1 mKate2MB9_B1 transformed with gDNA from TB501 amyE::Phyperspank-mKate2 (specR) (Dragoš et al. 2018)This study
B. subtilis MB9_B1 mKate2 ΔppsCMB9_B1 mKate2 transformed with gDNA from DS4114 ΔppsC::tetR (Müller et al. 2014)This study
B. subtilis MB9_B6 GFPMB9_B6 transformed with gDNA from TB500 amyE::Phyperspank-gfp (specR) (Mhatre et al. 2017).This study
F. oxysporum
IBT 40 872
From surface of an equipment in a Danish factory, 2005, by Anne SvendsenIBT Culture Collection of Fungi, DTU

gfp: gfpmut2, specR: spectinomycin resistance, and tetR: tetracycline resistance.

Bacterial inoculants were prepared from 1 ml overnight cultures, which were spun down and resuspended in 100 µl Milli-Q (MQ) water. 10 µl of the inoculum was transferred into 2 ml of competence medium (80 mmol/l K2HPO4, 38.2 mmol/l KH2PO4, 20 g/l glucose, 3 mmol/l Na3–citrate, 45 µmol/l ferric NH4–citrate, 1 g/l casein hydrolysate, 2 g/l K-glutamate, and 0.335 µmol/l MgSO4·7H2O) and incubated at 37°C with shaking for 3.5 h. Donor DNA was prepared from 1 ml of an overnight culture of the donor strain using Bacterial and Yeast Genomic DNA Purification Kit from EURx, yielding a typical DNA concentration in the range of 50–150 ng/µl. After the initial growth phase, 2 µl donor DNA was mixed with 400 µl competent cells and incubated further for 2 h before plating on selective agar plates which were incubated overnight at 37°C.

For interaction experiments on solid media, potato dextrose agar (PDA; Carl Roth) plates were prepared following the supplier’s instruction and plates were dried in a laminar flow cabinet for 1 h to avoid excessive expansion of B. subtilis colonies.

Population dynamics

Overnight cultures of MB9_B1 mKate2, MB9_B6 GFP, and MB9_B1 mKate2 ΔppsC were adjusted to 0.1 at OD600, and mixed in ratios of 100:1, 1:1, and 1:100 (MB9_B1/MB9_B1 ΔppsC and MB9_B6, respectively), and used in the experimental setups outlined below.

Population dynamics in liquid culture: spores of F. oxysporum were harvested from a 5-day grown culture cultivated in 50 ml PDB with shaking at room temperature and filtered through miracloth to remove remaining mycelia fragments before washing in sterile MQ water. The spore solution was diluted 500× followed by 2-fold serial dilution in PDB in a 96-well microtiter plate and 11 µl of the adjusted single strains and mixes were subsequently added into the microtiter plate to a final volume of 111 µl.

Growth of the bacterial populations in the presence or absence of Fusarium spores were followed using a microplate reader (BioTek Synergy HTX Multi-mode Microplate Reader). The plate was incubated at 28°C with shaking and measurements of OD600, green fluorescence for MB9_B6, and red fluorescence for MB9_B1 were taken every 15 min for 70 h (Optics position: Bottom, GFP: ex: 485/20 em: 528/20, Gain: 50, mKate2: Ex: 590/20, Em: 635/20, Gain 50).

Population dynamics on solid media: 5 µl of the single strains and mixed cocultures (100:1, 1:1, and 1:100) were spotted ~3 cm from the center of PDA medium containing plates and allowed to dry before spotting of 5 µl of F. oxysporum spore suspension in the center of the plates. The samples were incubated at room temperature. Growth and inhibition patterns were imaged every hour for 7 days using a ReShape imaging system (ReShape Biotech, Denmark) and annotated for when Fusarium reached the edge and outer center of the bacterial colonies. For fluorescence imaging, long pass emission filters at 530 and 630 nm and excitation at 467/5 nm and 527/5 nm for green and red fluorescence, respectively. Fluorescence images of mono and co-cultures without Fusarium were processed using ImageJ, briefly, image stacks were converted to 8-bit and colony region of interest (ROI) were selected from the normal light image by autothresholding using Otsu’s algorithm. Each fluorescence channel was measured for total area of the colony and the raw integrated density, which is a sum of pixel values in the selected area.

Results

Plipastatin is protective in a well-mixed liquid culture

The inhibitory effect of the B. subtilis and its secondary metabolite plipastatin on Fusarium is well documented using single bacterial cultures on agar medium (Kiesewalter et al. 2021, Kjeldgaard et al. 2022); however, its role within multistrain systems is less well understood. To dissect how plipastatin influences the dynamics between different B. subtilis strains, we established a two-species coculture system, in the same setup displaying inhibition of Fusarium, consisting of a plipastatin producer (MB9_B1) and a coisolated nonproducer isolate (MB9_B6) in the presence of germinating F. oxysporum spores. MB9_B1 and MB9_B6 have earlier been studied and shown to be close relatives (ANI of 99.95 for the four housekeeping genes gyrA, recA, dnaJ, and rpoB) and nonantagonistic toward each other, with MB9_B6 featuring a point-nonsense mutation in ppsB resulting in inability to produce plipastatin (Kiesewalter et al. 2021). As a control, a plipastatin deficient mutant of MB9_B1 (MB9_B1 ΔppsC) was also cocultivated in the presence of MB9_B6 to detect the influence of any other genetic differences between the two isolates.

Our liquid culture approach assayed survival and growth of the mono- and co-cultures of MB9_B1 or MB9_B1 ΔppsC and MB9_B6 against serial 2-fold dilutions of a freshly prepared F. oxysporum spore suspension (see Materials and Methods) ranging from 2-fold dilution to a 128-fold dilution. Coculture populations were further split into ratios of 100:1, 1:1, and 1:100 of MB9_B1 and MB9_B6, respectively, to assay possible density dependent protection.

In liquid media, the Fusarium challenge often resulted in a collapse of the bacterial culture inferred from the loss of fluorescence signal (Fig. 1, gray lines) or sharply downward trajectory before the end-point compared with the control (Fig. 1, orange lines). The isolate in lower concentration was undetectable from the background in 100:1 and 1:100 ratios, and these were excluded from further analysis (Figs S1 and S2).

Survival of monocultures when challenged with different concentrations of F. oxysporum spores in liquid cultures at 28°C over the course of 70 h. MB9_B1 and ΔppsC mutant was followed by red fluorescence while MB9_B6 was followed by green fluorescence. The bar to the right denotes the Fusarium concentration, with no-spores at the top, followed by the lowest spore dilution and increasing dilutions and thus lower concentration of spores toward the bottom. Surviving populations are colored in green, whereas populations that crashed judging by fluorescence signal are gray, and populations with sharp declines in the late stage are orange. Three replicates of each monoculture were followed.
Figure 1.

Survival of monocultures when challenged with different concentrations of F. oxysporum spores in liquid cultures at 28°C over the course of 70 h. MB9_B1 and ΔppsC mutant was followed by red fluorescence while MB9_B6 was followed by green fluorescence. The bar to the right denotes the Fusarium concentration, with no-spores at the top, followed by the lowest spore dilution and increasing dilutions and thus lower concentration of spores toward the bottom. Surviving populations are colored in green, whereas populations that crashed judging by fluorescence signal are gray, and populations with sharp declines in the late stage are orange. Three replicates of each monoculture were followed.

All three monocultures exhibited similar population crashes at high Fusarium loads of 2-, 4-, and 8-fold (2- and 4-fold included in Figs S3–S5), but differentiated at the 16-fold dilution and lower spore concentrations (Fig 1, green lines). MB9_B1 populations survived with one-third populations at 16-fold, two-thirds at 32- and 64-fold, and all three-thirds at 128-fold (Fig. 1, left panel) whereas MB9_B6 only featured two-third survivors at 128-fold (Fig. 1, right panel) and MB9_B1 ΔppsC featured two sharply decreasing populations even at 128-fold dilution of Fusarium spores (Fig. 1, middle panel).

The differences in survival indicated certain degree of plipastatin-mediated protection against germinating F. oxysporum spores in our setup. Next, we determined whether this effect was shared or privatized by the producing strain MB9_B1.

Cocultures were generally less robust against Fusarium invasion, likely due to halving of the producer density. In the 1:1 mixes, MB9_B1 did not survive at the 16-fold dilution or higher spore concentrations, and featured one-sixth, three-sixth, and five-sixth survivors at 32-, 64-, and 128-fold, respectively (Fig. 2, left panel) while MB9_B1 ΔppsC only showed one-third surviving populations at 128-fold spore dilution (Fig. 2, second panel from the left). Plipastatin producers were beneficial towards nonproducers, as MB9_B6 featured increased survival in the 1:1 coculture with MB9_B1 with one-sixth, three-sixth, and five-sixth surviving populations respectively at 32-, 64-, and 128-fold dilutions (Fig. 2, third panel from the left) versus one-third at the 128-fold dilution when MB9_B6 was cocultured 1:1 with MB9_B1 ΔppsC (Fig. 2, right panel).

Survival of cocultures when challenged with different concentrations of F. oxysporum spores in liquid cultures at 28°C over the course of 70 h. Cocultures with wild-type MB9_B1 are marked 1:1, while cocultures using the nonproducer mutant MB9_B1 ΔppsC are marked 1:1 Δ. MB9_B1 and ΔppsC mutant in the coculture was followed by red fluorescence while MB9_B6 was followed by green fluorescence. The bar to the right denotes the Fusarium concentration, with no-spores at the top, followed by the lowest spore dilution and increasing dilutions and thus lower concentration of spores toward the bottom. Surviving populations are colored in green, whereas populations that crashed judging by fluorescence signal are gray, and populations with sharp declines in the late stage are orange. 1:1 culture was performed in six replicates while 1:1 Δ was performed in three replicates.
Figure 2.

Survival of cocultures when challenged with different concentrations of F. oxysporum spores in liquid cultures at 28°C over the course of 70 h. Cocultures with wild-type MB9_B1 are marked 1:1, while cocultures using the nonproducer mutant MB9_B1 ΔppsC are marked 1:1 Δ. MB9_B1 and ΔppsC mutant in the coculture was followed by red fluorescence while MB9_B6 was followed by green fluorescence. The bar to the right denotes the Fusarium concentration, with no-spores at the top, followed by the lowest spore dilution and increasing dilutions and thus lower concentration of spores toward the bottom. Surviving populations are colored in green, whereas populations that crashed judging by fluorescence signal are gray, and populations with sharp declines in the late stage are orange. 1:1 culture was performed in six replicates while 1:1 Δ was performed in three replicates.

As cocultures proved somewhat beneficial to the nonproducer, we further inspected the growth of each species without Fusarium challenge to determine whether the nonproducer had a fitness benefit over the producer. Here, we noted only minor differences in growth characteristics of the monoculture with MB9_B6 featuring a slightly higher carrying capacity followed by MB9_B1 ΔppsC and close after MB9_B1 (Fig. 3A) although these differences were not significant when comparing at 20 h (Fig. S6). Likewise in the 1:1 cocultures, only minor nonsignificant differences at 20 and 70 h (Fig. S7) were observed. However, MB9_B1 wild-type featured slightly higher fluorescence signal compared to the plipastatin nonproducing MB9_B1 ΔppsC in coculture (Fig. 3B, middle panel), while MB9_B6 performed slightly better in coculture with the plipastatin nonproducer (Fig. 3B, right panel). Overall, we were unable to observe a significant fitness burden associated with plipastatin production and, on the contrary, MB9_B1 featured slightly higher growth compared to MB9_B1 ΔppsC in coculture with MB9_B6.

Growth of mono- and cocultures in the absence of Fusarium challenge, incubated at 28°C for 70 h and measured for optical density at 600 nm and fluorescence signal from RFP and GFP. (A) Monocultures of MB9_B1 (green), MB9_B1 ΔppsC (gray), and MB9_B6 (orange), left panel showing OD600, middle panel showing red fluorescence, and right panel showing green fluorescence. Line shows the mean of three replicates with ribbon denoting standard deviation. (B) Cocultures of MB9_B1 and MB9_B6 (green) and MB9_B1 ΔppsC and MB9_B6 (orange). 1:1 culture was performed in six replicates while 1:1 Δ was performed in three replicates, line shows the mean of these replicates and ribbon represents the standard deviation.
Figure 3.

Growth of mono- and cocultures in the absence of Fusarium challenge, incubated at 28°C for 70 h and measured for optical density at 600 nm and fluorescence signal from RFP and GFP. (A) Monocultures of MB9_B1 (green), MB9_B1 ΔppsC (gray), and MB9_B6 (orange), left panel showing OD600, middle panel showing red fluorescence, and right panel showing green fluorescence. Line shows the mean of three replicates with ribbon denoting standard deviation. (B) Cocultures of MB9_B1 and MB9_B6 (green) and MB9_B1 ΔppsC and MB9_B6 (orange). 1:1 culture was performed in six replicates while 1:1 Δ was performed in three replicates, line shows the mean of these replicates and ribbon represents the standard deviation.

Protection by the shared plipastatin depends on a structured environment

In nature, bacteria are often found in spatially structured environments, including biofilms, which lacks the homogeneous milieu of a well-mixed liquid culture, which might favor cheaters due to constant mixing. Therefore, we next tested whether plipastatin was also protective toward a nonproducer on solid agar plates. On solid agar medium, bacterial colonies and Fusarium spores were spaced 3 cm apart and left at room temperature to grow, expand, and interact with each other and timed for when Fusarium hyphae reached the edge of the colony (Fig. 4A) and the colony center (Fig. 4B).

Interactions between B. subtilis monocultures and Fusarium on solid media PDA, grown for 7 days at room temperature. Representative images of when F. oxysporum was determined to meet B. subtilis colonies at the edge (A) and center (B) with white arrows pointing out the less visible Fusarium front when invading into the center. Scale bars indicate 1 cm. (C) Boxplot displaying the timing of F. oxysporum encountering edge (left panel) and center (right panel) of monoculture colonies of MB9_B1, MB9_B1 ΔppsC, and MB9_B6 with three replicates of each.
Figure 4.

Interactions between B. subtilis monocultures and Fusarium on solid media PDA, grown for 7 days at room temperature. Representative images of when F. oxysporum was determined to meet B. subtilis colonies at the edge (A) and center (B) with white arrows pointing out the less visible Fusarium front when invading into the center. Scale bars indicate 1 cm. (C) Boxplot displaying the timing of F. oxysporum encountering edge (left panel) and center (right panel) of monoculture colonies of MB9_B1, MB9_B1 ΔppsC, and MB9_B6 with three replicates of each.

Generally, given enough time, F. oxysporum was able to overgrow B. subtilis colonies regardless of plipastatin production, but with considerable difference in the time points when hyphae reached the edge and center of the colony. For monoculture colonies, a clear inhibitory zone was observed for plipastatin producers with accompanying delay in Fusarium invasion when compared to nonproducers, taking 25 h longer before Fusarium reached the edge and 51 h before reaching the center of the colony (Fig. 4C) at which time nonproducers have already been fully overgrown by Fusarium. Representative images at days 5 and 7 and colony morphology of mono- and cocultures are displayed in Fig. S8.

Plipastatin production was beneficial in cocultures with invasion delay depending on the relative abundance of the producer strain, Fusarium reached the colony edge of the 1:100 culture first at 90 h, followed by 1:1 at 100 and finally the 100:1 colony at a 108 h delay (Fig. 5A). In all cases, the delay observed was smaller than that of the producer monoculture, which was likely due to the lower producer density compared with monoculture producer strain, MB9_B1. The importance of plipastatin was further underlined by cocultures replacing MB9_B1 with its plipastatin nonproducer derivative, MB9_B1 ΔppsC, in which no inhibitory zone or delay was observed for F. oxysporum (Fig. 5B). The time point when F. oxysporum reached the center of the bacterial colonies followed a similar trend as detected for the colony edge.

Interactions between B. subtilis cocultures and Fusarium on solid media PDA, grown for 7 days at room temperature. Boxplot showing the timing of F. oxysporum encountering edge (left top panel) and center (right top panel) of cocultures of MB9_B1 and MB9_B6 (left subpanels) and MB9_B1 ΔppsC and MB9_B6 (right subpanels) at ratios of 1:100, 1:1, and 100:1 of MB9_B1 and MB9_B6, respectively with three replicates of each.
Figure 5.

Interactions between B. subtilis cocultures and Fusarium on solid media PDA, grown for 7 days at room temperature. Boxplot showing the timing of F. oxysporum encountering edge (left top panel) and center (right top panel) of cocultures of MB9_B1 and MB9_B6 (left subpanels) and MB9_B1 ΔppsC and MB9_B6 (right subpanels) at ratios of 1:100, 1:1, and 100:1 of MB9_B1 and MB9_B6, respectively with three replicates of each.

Fluorescence labeling of the strains allowed us to roughly estimate the abundance of producers and nonproducer populations in the coculture colonies. Using the summed pixel values of the colony adjusted by its area for both fluorophores, we could assay the differences between cocultures without Fusarium challenge at day 5 of the experiment as a control. We noted little difference between cocultures featuring MB9_B1 or MB9_B1 ΔppsC for both red and green fluorescence, suggesting that plipastatin production had little effect on the composition of the cocultured colonies alone (Fig. 6).

Summed fluorescence intensity divided by the area of colonies for each coculture without Fusarium challenge at day 5 with MB9_B6 at ratios of 1:100, 1:1, and 100:1 of MB9_B1 and MB9_B6, respectively. Red fluorescence of MB9_B1 shown in left panel and green fluorescence of MB9_B6 shown in right panel. Cocultures with MB9_B1 shown in light green and cocultures with MB9_B1 ΔppsC shown in purple, boxplots are the result of three replicates.
Figure 6.

Summed fluorescence intensity divided by the area of colonies for each coculture without Fusarium challenge at day 5 with MB9_B6 at ratios of 1:100, 1:1, and 100:1 of MB9_B1 and MB9_B6, respectively. Red fluorescence of MB9_B1 shown in left panel and green fluorescence of MB9_B6 shown in right panel. Cocultures with MB9_B1 shown in light green and cocultures with MB9_B1 ΔppsC shown in purple, boxplots are the result of three replicates.

Discussion

The efficacy of biological control agents are highly dependent on the context of the environment as the resident microbiome members can interact with the biocontrol strain and its effector molecules. Here, we focused on the population dynamics of bioactive secondary metabolite sharing. Beneficial extracellular metabolite producing organisms are at risk of cheating, which could compromise the long-term survival of a producer strain. By cheating or exploitation, nonproducers might gain a fitness advantage over producers through saving metabolic resources and benefiting from the secreted compounds (so called public goods). Evolutionarily, the producer strain is out-competed, leading to the loss of the beneficial trait in the population (i.e. tragedy of commons). Therefore, producers must ensure cooperation or direct beneficial behavior toward kins that are more likely to reciprocate. Understanding how producers of active metabolites such as plipastatin function in the context of other strains is immensely important for formulating effective biosolutions against plant pathogens like F. oxysporum.

Here, we describe that protection against F. oxysporum by plipastatin is shared among B. subtilis isolates with variable production of this antifungal secondary metabolite. Our experiments indicate that protection by plipastatin is dependent on the density of the B. subtilis producer strain and also influenced by the level of Fusarium spores. Under planktonic conditions, we observed either normal growth or collapse of the bacterial coculture containing the mixture of plipastatin producer and nonproducer strains. The fate of the untagged Fusarium strain is difficult to determine in the platereader setup. Therefore, it is possible that Fusarium might be present in the surviving B. subtilis cultures, thus additional culturing or inclusion of a tagged Fusarium strain could provide additional insights into the dynamics between Bacillus and Fusarium and the role of plipastatin. In the solid agar media setup, Fusarium and B. subtilis were spatially separated at the beginning and allowed to meet during growth. Here, presence of plipastatin delayed the fungal invasion, but since Fusarium was not completely eradicated, B. subtilis colonies were eventually invaded by the fungus regardless of plipastatin producer strain proportion. The mechanism by which Fusarium invades and inhibits Bacillus has not yet been elucidated, but determination of those could provide further context to the role of plipastatin, which in our setup and experimental conditions would seem more efficient as a preventative measure against Fusarium establishment rather than during direct competition.

In the planktonic cocultures, efficient inhibition of Fusarium protected the cohabiting nonproducing strain, which raises the question of whether the plipastatin producing strains have fitness a disadvantage compared with the cohabiting nonproducers. However, our results show that the plipastatin producer MB9_B1 featured similar growth compared with nonproducer mutant derivative when cocultured with MB9_B6 under the planktonic conditions, thus we conclude that plipastatin production does not carry a fitness cost in these conditions. Although not significant, MB9_B1 displayed a trend of increased growth compared to the nonproducer, which could implicate plipastatin in intraspecific competition or could otherwise be caused by epistatic effects of plipastatin deletion in the B. subtilis isolate. To exclude other factors between MB9_B1 and MB9_B6 that might have played a role in the lack of fitness difference, a coculture using MB9_B1 and MB9_B1 ΔppsC should be tested in future studies with the addition CFU enumeration or flow cytometry to more accurately determine changes in the coulture population. Further, the timeframe of our experiments were all only one iteration lasting up to 3 days in liquid and 7 days on solid media, which might simply be too short period to detect potential fitness differences or changes in the population dynamics. Potentially, a successive experimental setup, in which the coculture is repeatedly transferred and challenged, could potentially allow to determine the evolutionary stability of plipastatin production in the mixture of plipastatin producer and nonproducer strains.

Additionally, other secondary metabolites might influence the stability of microbial consortia, offering another stabilizing mechanism, like division of labor, whereby each partner organism specialize in production of one or a subset of metabolites or matrix components, which has been documented in B. subtilis in earlier literature (Gestel et al. 2015, Liu et al. 2015, Dragoš et al. 2018, Jautzus et al. 2022, Yannarell et al. 2023). Following expression or directly the concentration of plipastatin and other SMs in each isolate during coculture and when being challenged by Fusarium could shed light on potential task distribution, prudent metabolism, and subpopulations that the coculture might develop and determine whether the presence of Fusarium affects production of plipastatin. Regulation of the biosynthetic gene clusters in B. subtilis has been demonstrated to be influenced by intra- and interspecies signals (Dinesen et al. 2025, Lozano-Andrade et al. 2025).

While we have demonstrated that plipastatin protects kin nonproducers in liquid cultures and on solid agar medium, there are still many questions left unanswered regarding the stability of the behavior, plipastatin expression patterns in the producer population, and the role of other secondary metabolites in inhibiting F. oxysporum. However, we hope this study will aid in future understanding of how B. subtilis functions within a consortium.

Author contributions

Rune Overlund Stannius (Conceptualization, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft), and Ákos T. Kovács (Conceptualization, Funding acquisition, Supervision, Writing – original draft)

Conflict of interest

None declared.

Funding

This project was supported by the Danish National Research Foundation (DNRF137) for the Center for Microbial Secondary Metabolites to Á.T.K.

References

Bawa
 
I
.
Management strategies of Fusarium wilt disease of tomato incited by Fusarium oxysporum f. sp. lycopersici (Sacc.): a review
.
Int J Adv Acad Res
.
2016
;
2
:
32
42
.

Blake
 
C
,
Christensen
 
MN
,
Kovács
 
ÁT
.
Molecular aspects of plant growth promotion and protection by Bacillus subtilis
.
MPMI
.
2021
;
34
:
15
25
.

Compant
 
S
,
Samad
 
A
,
Faist
 
H
 et al.  
A review on the plant microbiome: ecology, functions, and emerging trends in microbial application
.
J Adv Res
.
2019
;
19
:
29
37
.

Dean
 
R
,
van Kal
 
JAL
,
Pretorius
 
ZA
 et al.  
The top 10 fungal pathogens in molecular plant pathology
.
Mol Plant Pathol
.
2012
;
13
:
414
30
.

Dinesen
 
C
,
Vertot
 
M
,
Jarmusch
 
SA
 et al.  
Subtilosin a production is influenced by surfactin levels in Bacillus subtilis
.
Microlife
.
2025
;
6
:
uqae029
.

Dragoš
 
A
,
Kiesewalter
 
H
,
Martin
 
M
 et al.  
Division of labor during biofilm matrix production
.
Curr Biol
.
2018
;
28
:
1903
13.e5
.

El-Aswad
 
AF
,
Aly
 
MI
,
Alsahaty
 
SA
 et al.  
Efficacy evaluation of some fumigants against Fusarium oxysporum and enhancement of tomato growth as elicitor-induced defense responses
.
Sci Rep
.
2023
;
13
:
2479
.

Farace
 
G
,
Fernandez
 
O
,
Jacquens
 
L
 et al.  
Cyclic lipopeptides from Bacillus subtilis activate distinct patterns of defence responses in grapevine
.
Mol Plant Pathol
.
2015
;
16
:
177
87
.

Gestel
 
J
,
Vlamakis
 
H
,
Kolter
 
R
.
From cell differentiation to cell collectives: B acillus subtilis uses division of labor to migrate
.
PLOS Biol
.
2015
;
13
:
e1002141
.

Gong
 
A-D
,
Li
 
H-P
,
Yuan
 
Q-S
 et al.  
Antagonistic mechanism of iturin a and plipastatin a from Bacillus amyloliquefaciens S76-3 from wheat spikes against Fusarium graminearum
.
PLoS One
.
2015
;
10
:
e0116871
.

Honma
 
M
,
Tanaka
 
K
,
Konno
 
K
 et al.  
Termination of the structural confusion between plipastatin A1 and fengycin IX
.
Bioorg Med Chem
.
2012
;
20
:
3793
8
.

Hudson
 
O
,
Waliullah
 
S
,
Ji
 
P
 et al.  
Molecular characterization of laboratory mutants of Fusarium oxysporum f. sp. niveum resistant to prothioconazole, a demethylation inhibitor (DMI) fungicide
.
J Fungi
.
2021
;
7
:
704
.

Jautzus
 
T
,
van Gestel
 
J
,
Kovács
 
ÁT
.
Complex extracellular biology drives surface competition during colony expansion in Bacillus subtilis
.
ISME J
.
2022
;
16
:
2320
8
.

Kiesewalter
 
HT
,
Lozano-Andrade
 
CN
,
Wibowo
 
M
 et al.  
Genomic and chemical diversity of Bacillus subtilis secondary metabolites against plant pathogenic fungi
.
mSystems
.
2021
;
6
:
e00770
20
.

Kjeldgaard
 
B
,
Neves
 
AR
,
Fonseca
 
C
 et al.  
Quantitative high-throughput screening methods designed for identification of bacterial biocontrol strains with antifungal properties
.
Microbiol Spectr
.
2022
;
10
:
e0143321
.

Kunst
 
F
,
Rapoport
 
G
.
Salt stress is an environmental signal affecting degradative enzyme synthesis in Bacillus subtilis
.
J Bacteriol
.
1995
;
177
:
2403
7
.

Liu
 
J
,
Prindle
 
A
,
Humphries
 
J
 et al.  
Metabolic co-dependence gives rise to collective oscillations within biofilms
.
Nature
.
2015
;
523
:
550
4
.

Lozano-Andrade
 
CN
,
Dinesen
 
C
,
Wibowo
 
M
 et al.  
Surfactin facilitates establishment of Bacillus subtilis in synthetic communities
.
ISME J
.
2025
;
19
:
wraf013
.

Lugtenberg
 
B
,
Kamilova
 
F
.
Plant-growth-promoting rhizobacteria
.
Annu Rev Microbiol
.
2009
;
63
:
541
56
.

Mendes
 
R
,
Garbeva
 
P
,
Raaijmakers
 
JM
.
The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms
.
FEMS Microbiol Rev
.
2013
;
37
:
634
63
.

Mhatre
 
E
,
Sundaram
 
A
,
Hölscher
 
T
 et al.  
Presence of calcium lowers the expansion of Bacillus subtilis colony biofilms
.
Microorganisms
.
2017
;
5
:
7
.

Müller
 
S
,
Strack
 
SN
,
Hoefler
 
BC
 et al.  
Bacillaene and sporulation protect Bacillus subtilis from predation by myxococcus xanthus
.
Appl Environ Microbiol
.
2014
;
80
:
5603
10
.

Ongena
 
M
,
Jacques
 
P
.
Bacillus lipopeptides: versatile weapons for plant disease biocontrol
.
Trends Microbiol
.
2008
;
16
:
115
25
.

Ongena
 
M
,
Jourdan
 
E
,
Adam
 
A
 et al.  
Surfactin and fengycin lipopeptides of Bacillus subtilis as elicitors of induced systemic resistance in plants
.
Environ Microbiol
.
2007
;
9
:
1084
90
.

Raaijmakers
 
JM
,
Paulitz
 
TC
,
Steinberg
 
C
 et al.  
The rhizosphere: a playground and battlefield for soilborne pathogens and beneficial microorganisms
.
Plant Soil
.
2009
;
321
:
341
61
.

Sondergaard
 
TE
,
Fredborg
 
M
,
Oppenhagen Christensen
 
A-M
 et al.  
Fast screening of antibacterial compounds from Fusaria
.
Toxins
.
2016
;
8
:
355
.

Steinke
 
K
,
Mohite
 
OS
,
Weber
 
T
 et al.  
Phylogenetic distribution of secondary metabolites in the Bacillus subtilis species complex
.
mSystems
.
2021
;
6
:
e00057
21
.

Summerell
 
BA
.
Resolving Fusarium: current status of the genus
.
Annu Rev Phytopathol
.
2019
;
57
:
323
39
.

Thérien
 
M
,
Kiesewalter
 
HT
,
Auria
 
E
 et al.  
Surfactin production is not essential for pellicle and root-associated biofilm development of Bacillus subtilis
.
Biofilm
.
2020
;
2
:
100021
.

Venkatesh
 
N
,
Keller
 
NP
.
Mycotoxins in conversation with bacteria and fungi
.
Front Microbiol
.
2019
;
10
:
403
.

Yannarell
 
SM
,
Beaudoin
 
ES
,
Talley
 
HS
 et al.  
Extensive cellular multi-tasking within Bacillus subtilis biofilms
.
mSystems
.
2023
;
8
:
e0089122
.

Zhao
 
P
,
Quan
 
C
,
Wang
 
Y
 et al.  
Bacillus amyloliquefaciens Q-426 as a potential biocontrol agent against Fusarium oxysporum f. sp. spinaciae
.
J Basic Microbiol
.
2014
;
54
:
448
56
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License (https://creativecommons.org/licenses/by-nc-nd/4.0/), which permits non-commercial reproduction and distribution of the work, in any medium, provided the original work is not altered or transformed in any way, and that the work is properly cited. For commercial re-use, please contact [email protected]

Supplementary data