Abstract
Microbial secondary metabolites play important roles in biotic interactions in microbial communities and yet, we do not understand how these compounds impact the assembly and development of microbial communities. To address the implications of microbial secondary metabolite production on biotic interactions in the assembly of natural seawater microbiomes, we constructed a model system where the assembly of a natural seawater biofilm community was influenced by the addition of the marine biofilm forming Phaeobacter inhibens that can produce the antibiotic secondary metabolite tropodithietic acid (TDA), or a mutant incapable of TDA production. Because of the broad antibiotic activity of TDA, we hypothesized that the potential of P. inhibens to produce TDA would strongly affect both biofilm and planktonic community assembly patterns. We show that 1.9 % of the microbial composition variance across both environments could be attributed to the presence of WT P. inhibens, and especially genera of the Bacteriodetes were increased by the presence of the TDA producer. Moreover, network analysis with inferred putative microbial interactions revealed that P. inhibens mainly displayed strong positive associations with genera of the Flavobacteriaceae and Alteromonadaceae, and that P. inhibens acts as a keystone OTU in the biofilm exclusively due to its potential to produce TDA. Our results demonstrate the potential impact of microbial secondary metabolites on microbial interactions and assembly dynamics of complex microbial communities.
Similar content being viewed by others
Introduction
Complex and dynamic microbial communities colonize nearly every ecosystem in marine environments, and major biogeochemical processes, such as carbon cycling, are largely driven by microorganisms1. Both stochastic and deterministic abiotic and biotic processes can affect the assembly and composition of microbial communities2. Abiotic factors such as temperature and resource availability have a major influence on the taxonomic composition of microbiomes3,4, but inter-kingdom and interspecific interactions between co-occurring organisms may play equally important roles5,6. Competitive as well as cooperative interactions within microenvironments may be mediated by the production of secondary (or ‘specialized’) metabolites produced by a multitude of microorganisms7,8,9. These metabolites have predominantly been exploited because of their therapeutic, and agri- and aqua-cultural benefits. Studies on their potential ecological role(s) have generally described them as mediators of interference competition between microbes10, and most experiments have been conducted in co-cultures, or as traditional antibiotic screening assays11,12. The perception that microbial interactions are primarily negative has however recently been challenged as growth-promoting positive interactions have proven prevalent among soil community members13.
Mechanistic understanding of how microbial secondary metabolite production drives the assembly of natural microbial communities remains limited, largely due to the methodological challenge of (re)creating complex natural environmental conditions in the laboratory. The introduction of cell-free extracts containing secondary metabolites can either reduce or enrich specific taxa in semi-natural mesocosm communities14. While the addition of defined concentrations of pure compounds resulting in a homogenous suspension may resemble the planktonic environmental conditions of free-living bacteria, it is far from conditions found in bacterial biofilm. The surface associated bacteria in biofilms live in spatially structured communities, in which microscale gradients occur due to the consumption of resources and secretion of secondary metabolites15,16. In this study, we constructed a model system allowing for the formation of such spatially structured communities and introduced a bacterial producer of a potent antibiotic secondary metabolite.
The Rhodobacterales are ubiquitous in the marine environment, where in particular the Roseobacter group such as Tritonibacter and Phaeobacter spp. are primary colonizers on freshly submerged surfaces17. While the Tritonibacter genus (formerly Rugeria) have been isolated from the free-living planktonic phase18 the genus Phaeobacter has typically been isolated from marine biofilms19,20. Also, Phaeobacter are capable of invading preformed bacterial biofilms21,22,23. Several Phaeobacter spp. can produce the antibacterial secondary metabolite tropodithietic acid (TDA)24,25, and previous studies have suggested that TDA has multiple additional functions, such as being a weak iron chelator and having hormetic effects on nearby organisms. Thus, TDA may act as a mediator of antibiosis, or interspecific signaling in sub-inhibitory concentrations26,27.
To determine how such a secondary metabolite producer affects biofilm formation and the proximate planktonic microbiome, we introduced a wild type and a mutant incapable of producing this particular metabolite, respectively, into our model systems. Thus, the TDA producing P. inhibens and the TDA deficient mutant were separately allowed to colonize abiotic surfaces that were then submerged in natural seawater. Here, we hypothesized that the wildtype Phaeobacter, capable of TDA production, would cause changes in the microbiome composition during assembly of both environments potentially eliminating specific taxonomic groups. Deciphering such changes will provide insights into how microbial secondary metabolites can mediate ecological microbial community associations.
Results
A semi-natural marine microbial model system was constructed to determine if, and how, the potential of a bacterium to produce a potent secondary metabolite affected biofilm formation. The WT and the dTDA P. inhibens strains were allowed to pre-form biofilms and the subsequent colonization and change in microbial composition in a natural seawater system was studied and compared to changes in a system without pre-colonization. Since metabolic interactions change over time, we sampled over a 10-day period. To investigate the potential production of TDA in the axenic model system with the WT, we constructed a GFP reporter fusion in the P. inhibens DMS17935 strain to track expression from the tdaCDE promoter as a proxy for TDA production. Qualitative GFP fluorescence detection was observed after day 0, 1, 4 an 10 days, indicating continuous transcription of the tdaCDE operon and suggesting potential production of TDA in this system (Supplementary Fig. 1).
Quantification of absolute microbial abundances and Phaeobacter colonization
In the non-treated system, the biofilm bacterial community from natural seawater reached 7.33 ± 0.3 log10 CFU/cm2 after one day. This was significantly lower than the biofilm level of 8.87 ± 0.1 and 8.41 ± 0.5 log10 CFU/cm2 in the WT and the dTDA (LMM & EMM; p-value < 0.05), respectively. After the first day, the cell densities across all treatments gradually decreased over time. Both the WT and the dTDA treated systems had similar and non-significantly different total abundances up until day 4 (LMM & EMM; p-value = 1.00). At day 10, the total bacterial abundances of the WT treated system had significantly decreased 1 log10 fold compared to the dTDA system (LMM & EMM; p-value = 0.03) and was similar to the control treatment corresponding to 5.36 ± 0.7 log10 CFU/cm2 (Table 1; Fig. 1).
The WT and the dTDA P. inhibens strains readily colonized steel surfaces to the same level (LMM & EMM; p-value = 1) of 8.49 ± 0.3 and of 8.49 ± 0.2 log10 CFU/cm, respectively, before being exposed to the natural sea water (Table 1; Fig. 1).
Even though the surfaces were rinsed, patches of biofilm still detached from the steel coupons upon submersion in seawater and levels of approx. 107 CFU/mL of the WT and dTDA were detected in the planktonic suspensions. After one day, the WT and dTDA were reduced to 8.31 ± 0.2 and 7.93 ± 0.4 log10 CFU/cm in the biofilm (LMM & EMM; p-value=0.951) and to 6.2 ± 0.3 and 6.52 ± 0.3 log10 CFU/mL in the seawater (LMM & EMM; p-value = 1). Phaeobacter abundances continued to decrease after four days in all systems; in the biofilm to 6.12 ± 0.4 and 6.25 ± 0.8 log10 CFU/cm (LMM & EMM; p-value = 1) and in the planktonic suspension 4.52 ± 0.3 and 4.78 ± 0.4 log10 CFU/mL (LMM & EMM; p-value = 0.829). Numbers decreased further after ten days of incubation with a significant additional decrease in the numbers of WT as compared to the dTDA and 2.85 ± 1.5 and 5.06 ± 1.1 log10 CFU/cm2 were detected in the biofilm (LMM & EMM; p-value=0.0001) and 3.02 ± 0.7 and 3.67 ± 0.5 log10 CFU/mL in the planktonic suspension (LMM & EMM; p-value = 0.029). No Phaeobacter spp. were detected by the qPCR in the natural seawater (No ct values), however, a short-term Phaeobacter spp. colonization of stainless steel was observed at day one, corresponding to 3.43 ± 0.1 log10 CFU/cm, yet significantly lower as compared to the inoculated Phaeobacter spp. abundances (p-value < 0.001). The remaining numbers in the Control across all days and treatments were < 0.5 log10 CFU/cm. As expected, the relative abundance of Phaoebacter (WT and dTDA) determined from the amplicon sequencing was high in samples pre-inoculated with the strains (Supplementary Fig. 2B).
Factors driving the microbial community composition
We assumed that time (day 1, 4, or 10), treatment (Control, dTDA or WT systems) and environment (biofilm or planktonic suspension) would predict the microbial community assembly. The P. inhibens OTU were removed from the analyses to exclude the variance caused by its high relative abundance. The beta-dispersion was non-significant (p-value > 0.05) and similar across all experimental combinations.
Multivariate analysis using PERMANOVA revealed that 76.1 % of the variation in the microbial composition across the full data set was explained either individually or in combinations of time, environment and treatment (Table 2). Time was the main factor driving diversity in the bacterial communities (R2 = 0. 417, p-value = 0.001), followed by the treatment (R2 = 0. 122, p-value = 0.001) although the variance explained by the treatment was dependent on time (R2 = 0. 091, p-value = 0.001). Differences were also dictated by the environment e.g. the compositional variance between biofilm compared to the planktonic suspension (R2 = 0.059, p-value = 0.001) and this effect was dependent on time (R2 = 0.041, p-value = 0.001) and treatment (R2 = 0.014, p-value = 0.001). Multilevel pairwise PERMANOVA between the WT and dTDA systems furthermore showed that community composition was directly influenced by the potential of P. inhibens to produce TDA (Table 2; R2 = 0. 019, p-value = 0.001), although this effect was significantly mediated by interaction with time (R2 = 0.015, p-value = 0.001). The significant conditions driving community assembly were further supported by the clustering in nMDS ordination plot (Fig. 2).
Phaeobacter pre-colonization reduces the microbial community complexity
Alpha diversity was determined as the total sum of observed unique genera across all samples. The richness decreased significantly in all systems as compared to the natural seawater samples (Dunnest test; p-values < 0.01; Table 1; Supplementary Fig. 2A). Only the Control system on day 4 had a complexity similar to the natural seawater (Dunnest test; p-value = 0.3061). After the first day of incubation, the difference between the observed richness of the WT and dTDA systems in the biofilm was non-significant, comprising 75 ± 47 and 107 ± 54 unique genera, yet both were substantially lower than the complexity of the Control system 253 ± 18 (LMM & EMM; p-values < 0.001). A similar pattern was seen for the planktonic suspension of the WT and dTDA systems, where the observed richness was reduced to 92 ± 36 and 108 ± 35 unique genera compared to the Control system 206 ± 32 (LMM & EMM; p-values < 0.001), respectively. The complexity remained unchanged and non-significant between treatments on day 4 and day 10 in both environments (LMM & EMM; p-values > 0.05).
Identification of genera that solely respond to the potential TDA production by Phaeobacter inhibens
We used the ANCOM-BC model to analyse differential relative abundance between the WT versus the dTDA treated systems stratified by day and environment. Of 1056 different genera found in the total dataset, 219 unique genera were increased in their relative abundances by the WT and 176 genera were reduced relative to the dTDA system over the 10-day period in both environments. Because the largest sum of differentially abundant genera was observed at day 4, we assumed that this was also the day on which the presence of the TDA producing wildtype had the most marked effect (Fig. 3A). Therefore, we re-analyzed the differential abundance in all three combinations at order level to get an overall view on the taxonomic profile from day 4 (Fig. 3B).
Orders that were specifically increased by the WT were the Flavobacteriales (27.4 %), Sphingobacteriales (8.8 %), Cytophagales (1.8 %), Actinomycetales (0.9 %), Rhizobiales (0.6 %), Phycisphaerae (0.4 %), and Methylophilales (0.1 %), whereas only Lentisphaerales (0.2 %), while only an unclassified order of Deltaproteobacteria (0.6 %) were increased by the dTDA strain in the biofilm. In the planktonic suspension, again the Flavobacteriales (11.9 %) were increased in the presence of the WT compared to dTDA strain, together with the Bdellovibrionales (7.5 %), Chloroplast (0.2 %) and Methylophilales (0.1 %). Several Flavobacteria were significantly increased in the WT treated system in both environments compared to both the dTDA and the control system (Supplementary Fig. 3).
The TDA-producing P. inhibens becomes a keystone OTU in the microbial network in surface colonization
Since the pre-colonization of P. inhibens had a marked effect on microbial community composition, we used network analysis to provide insight into the microbial community structure and to disentangle whether the P. inhibens OTU abundances were directly co-occurring with other genera that also increased in relative abundance in response to the presence of the WT. Twelve paired networks were separately constructed for each day and environment with the network comparison module of NetCoMi, since PERMANOVA and ANCOM-BC analysis revealed that the impact of TDA-producing Phaeobacter depends on its spatial organization and time. For each pairwise network comparison, we compared two network properties that have been associated with stability of ecological community studies: (1) the positive edge percentage and (2) the modularity28,29,30. Next, we identified hubs or ‘keystone genera’ among common nodes between the WT and dTDA systems. All network scores for each node are listed in Supplementary Table 1.
The positive edge percentage and the modularity of each network comparison were similar and non-significant between the WT and the dTDA systems across time and each environment, except for the planktonic samples at day 4 and the biofilm at day 10. A weakly significant difference in the modularity of the dTDA system was calculated in the planktonic suspension at day 4 (Permutation tests; p-value = 0.012). In contrast, the opposite was found for the modularity of the WT in the biofilm at day 10 (Permutation tests; p-value = 0.039), suggesting that neither the WT nor the dTDA had any evident effect on the resilience of the microbiome.
At the first day, both comparison networks were dominated by the most influential nodes being genera of the Bacteroidetes, where in particular genera of the Flavobacteriia were among the most central nodes in both treatments, yet no significant difference was detected between the two treatments.
Interestingly, the WT P. inhibens node was identified as a hub, with a significantly higher eigenvector centrality as compared to the dTDA P. inhibens OTU (Permutation tests; p-value = 0.023) in the biofilm. The WT OTU displayed a high degree of closeness and low betweenness, where 33 of the 55 connections were positively associated with Gammaproteobacteria, Flavobacteriia and Cytophagia genera. All positive edge connections of this node were to other nodes showing high eigenvector centralities as compared to the negative correlations to nodes with very low eigenvector centralities (scores < 0.25; Fig. 4). We could further confirm this by the univariate ANCOM-BC analysis, where we saw the 12 nodes positively connected to the WT P. inhibens node were also significantly increased in their relative abundance at day 4 in the biofilm. The WT node also showed a higher eigenvector centrality in the planktonic suspension at day 4 as compared to the dTDA, however this was non-significant. Only few nodes (6 of 135) were connected to the P. inhibens dTDA node and what instead seemed to be the most central node was the Stenotrophomonas, an opportunistic fish pathogen (Permutation tests; p-value = 0.023), which had no connection to the dTDA node. At day 10, no difference was observed between the centrality measures for the WT and dTDA nodes, with the only exception that the dTDA node showed higher closeness centrality to other nodes in the planktonic suspension as compared to the WT (Permutation tests; p-value = 0.008).
Discussion
Microbial secondary metabolites such as antibiotics and siderophores are important mediators of species interactions in many ecological niches, however, whether these secreted metabolites are predominantly weapons of warfare or have other and broader roles being involved in interspecific signaling between species is being discussed8,31,32,33,34,35. Nevertheless, it is unknown how such metabolites may mediate microbial community structure and assembly patterns. In this study, we constructed a semi-natural model system to mimic the structural heterogeneity of natural bacterial biofilms. In order to determine the microscale effect on the microbiome assembly by a sequence-based approach, we boosted the system, allowing a naturally occurring biofilm-forming secondary metabolite producer to pre-colonize the system. Thus, we manipulated the microbial assembly patterns in semi-synthetic communities to assess whether the potential to produce the secondary metabolite, TDA, affected biofilm assembly and development. We demonstrate a minor, yet significant, level of explained variance (1.9 %) from the potential of P. inhibens to produce TDA across all time points but also observed a significant dependence of time on this effect (1.5 %) in agreement with the most pronounced effect seen at day 4 where P. inhibens temporarily became a keystone species solely through its potential to produce TDA. Eventhough TDA was below detection limit in the chemical bulk-extractions, we captured GFP flouresence indicating transcription activity of genes responsible for the TDA biosynthesis in the GFP-reporter fusion strain, suggesting that TDA could be produced when P. inhibens is established in the biofilm and the surface associated community.
The semi-natural model system likely experienced nutrient depletion over time, and this may explain the overall decrease of the surface colonizing and free-living Phaeobacter spp. in the treated systems, and the increase of relative abundances of autotrophic taxa. The reduction of both Phaeobacter strains could further reflect the transition from the initial inoculated abundances going towards a community composition equilibrium as previously described for similar systems36,37.
Majzoub et al. 38 compared the phenotypic traits between a P. inhibens strain with the capability to produce TDA and mutants with mutations in the TDA gene cluster, and showed that biofilm persistence was independent of the capability to produce TDA. Consistent with this study, no difference was observed between the total abundance of the WT and dTDA strain in the beginning. Yet, the significant decrease of the WT compared to the dTDA strain during the last time point at day 10 suggests that the production of TDA and tolerance maintenance may become a metabolic burden for P. inhibens.
Time was the main factor influencing microbial community assembly and this, combined with the presumed decrease in nutrient availability, reflects that stochastic drift and genetic diversification are inherent properties of time-dependent ecological studies2,6. Nonetheless, the potential to produce TDA had a significant impact on the assembly of the microbial community composition. In some Roseobacter group species, TDA biosynthesis is most prominent under static growth conditions39; however, P. inhibens DSM17935 also produces TDA under shaking conditions40, which may explain the significant change in the planktonic suspension under the assumption that TDA is being produced during the course of the experiment.
The present study specifically aimed to test the potential role of the presence of a secondary metabolite producer in community assembly, and given that TDA affect a wide range of microorganisms25,41, it was surprising that most genera were increased by the WT but not the dTDA strain. The largest differential changes were observed on day four, where members of the Bacteroidetes, i.e. the Flavobacteriales, Sphingobacteriales and Cytophagales were all significantly enriched by the WT compared to the dTDA systems. The positive association between the P. inhibens when it is capable of producing TDA and the Flavobacteriales including Winogradskyella spp. has been observed in previous studies describing different marine systems such as non-axenic algal and oyster microbiomes42,43. The fact that several members of the Bacteroidetes are motile by gliding39 and that the TDA-producing Phaeobacter species most frequently have been isolated from harbor surfaces20,44, suggests that they could compete for surface colonization space in nature. However, the role of TDA as a metabolite broadly facilitating interference competition seems unlikely between P. inhibens and genera of the Bacteriodetes considering the positive effect of the potential TDA production.
Whether the TDA production was directly involved in an actual facilitation of the increase in these different orders is difficult to ascertain, and a variety of mechanisms could explain the direct or indirect increase of these. Members of the Bacteroidetes are proficient biopolymer degraders45,46 and the increase of Bacteroidetes including Flavobacteria spp. could be caused by their capability to exploit the biofilm matrix of the WT as opposed to dTDA strain, as mutations in the TDA gene cluster affect extracellular polymeric substance composition38. TDA or the producing bacteria typically have a reducing effect on fast-growing heterotrophic bacteria and we also observed a significant reduction of several members of the Gammaproteobacteria47,48. Because TDA has been proposed to act as an antiporter49, fast-growing bacteria, such as the vibrios50, with high metabolic turnover may be more antagonized as compared to slower growing genera of the Bacteriodetes. Ultimately, the increase may have been a result of the indirect effect of the reduction of natural parasites, predators or other species competing for the same niche that are sensitive to TDA in combination with the fact that several Flavobacteriia genera and other members of the Bacteriodetes exhibit natural tolerance to TDA51,52. Possible hypotheses for the increase of of Flavobacteria spp. are summarized in Fig. 5.
The Rhodobacteraceae comprising the Phaeobacter species and other TDA-producing bacteria have previously been suggested to serve as keystone taxa in algal microbial communities48. Because of the antagonizing property of TDA, we hypothesized that the WT system could act as an individual stressor and cause perturbation to the system. Indications for perturbed microbiomes have been associated with reduced modularity, decreased microbial diversity and higher positive associations between taxa28, however, no difference between the global network properties were observed between the WT and the dTDA network comparisons (Supplementary Table 1).
Node locations varied across the environments and time within both treatments (Fig. 4A; Supplementary Table 1), indicating that the differences in potential microbial interactions were determined by initial pre-colonization of either the WT or the dTDA. Keystone taxa have previously been defined as taxa that exert their influence on microbiome functioning irrespective of abundance as compared to dominant taxa that often affects the functioning of the ecosystem, or a specific process solely by virtue of abundance53. While the relative abundance of both WT and dTDA on day 4 was reduced to < 16 %, the WT was simultaneously identified as a keystone OTU along with genera of Flavobacteriaceae and Alteromonadaceae previously identified as keystone taxa in aquatic ecosystems48,53. The fact that the WT was solely identified as a hub in the biofilm, and that the dTDA exhibited no significant network importance suggest that TDA may play a significant role in microbial community assembly during surface colonization.
The WT displayed 1.5 times more positive than negative connections to genera that were highly associated in the microbiome (Fig. 4B). In a recent study by Kehe et al 13, growth-promoting positive interactions between metabolic dissimilar genera, and especially between strong and slow-growing bacteria were more common than previously expected5,11. Even though the microbial interactions in Kehe et al. 13 were mainly governed by carbon source utilization capabilities, our results indicate that the potential of microbial secondary metabolite production, such as TDA, may be linked to potential positive interactions between microbes. Inferring true ecological interactions from network analyses and keystones based on empirically chosen cut-offs should be interpreted with caution as correlation does not equal causation. Hence, bottom-up culture-based studies involving isolates from the same system combined with metabolomics should be used to investigate whether TDA is produced and directly involves positive interactions.
In conclusion, this study suggests that the presence of Phaeobacter inhibens, its spatial organization, and its potential to produce TDA can have a considerable impact on the microbial assembly pattern of a marine biofilm. Both differential and network co-occurrence analyses indicated that a large part of the potential interactions between P. inhibens and other community members were positive and either indirectly or directly increased the abundances of other species, indicating that TDA may not only be involved in interference competition.
Materials & methods
Bacterial strains and growth conditions
The wildtype, Phaeobacter inhibens DSM1739554 (WT) and the gentamicin-resistant transposon mutant P. inhibens DSM17395 tdaB::GmR devoid of TDA production55 (dTDA) were revived from frozen -80 °C stocks and streaked on Marine Agar (MA, Difco 2216). The GFP-reporter fusion strain, DSM17395-GFP (pBBR1MCS2-PtdaCDE-GFP) and the control strain, DSM17395-NC (pBBR1MCS2) were grown on MA or marine broth (MB) with 100 μg/ml kanamycin (Kan100). Individual colonies were selected and used as biological triplicates of the wildtype and the mutant, respectively, in the following experiments. All cultures were grown at 25 °C in Marine Broth (MB) or on MA.
Construction of the semi-natural marine model system
Stainless steel coupons (10 × 20 × 1 mm) were cleaned and sterilized as previously described56 and clamped vertically in circular stainless steel racks holding up to 20 coupons and placed in 1 L beakers. Single colonies of the WT and dTDA P. inhibens were grown in MB overnight. One milliliter of each culture with a mean ± standard deviation (sd) of 8.6 ± 0.04 log10 CFU/mL and 8.9 ± 0.15 CFU/mL, respectively, were separately transferred to 250 mL MB in 1 L beakers and allowed to form biofilm on the stainless steel coupons. Non-coated sterile coupons were also immersed in 250 mL MB and used as negative control (Control), and all treatments (3 × 1 L beakers per treatment) were left for three days at 25 °C with agitation (magnetic stirrer; 175 rpm). Coupons with adhered bacteria were gently washed by immersing the circular racks together with coupons twice in 300 mL 3 % Sigma sea salt (Sigma-Aldrich) to remove non-adhering bacteria. Racks and coupons were then submerged in natural seawater collected in Hellerup harbor, Denmark (55.7319, 12.5810), on September 24, 2018 at 10.45 a.m. at approximately 0.1-meter depth (13.5 °C, 3.32 mg L-1 dissolved oxygen [DO], 24.90 [PSU]). The total number of bacterial cells in seawater was determined using SYBR Gold staining and counting by epifluorescence microscopy57. The racks with biofilm-coated coupons and the non-treated Control were separately incubated with 250 mL seawater at 25 °C with agitation (magnetic stirrer; 175 rpm) for 10 days. Bacteria in the biofilm were harvested by transferring single coupons to polystyrene tubes containing 2 mL 3 % SSS. Tubes were sonicated for 4 min at 25 °C (28-kHz, 2 × 150 W sonication bath, Delta 220, Deltasonic, Meaux, France) and vortexed at maximum speed for 15 s to detach bacteria from the surfaces58. The bacteria detached from the biofilm and 2 mL of the planktonic suspension were centrifuged at 8,000 × g for 20 min and the pellet was then used for DNA extractions. Samples were collected from all treatments (in three biological replicates) at day 0, 1, 4, and 10 in three technical replicates per beaker.
Construction of a DSM17395 GFP reporter fusion strain
The promoter region of the tdaCDE operon was amplified using the primer pairs PtdaB-Pinh_Fw (5′-AAAACTCGAGCTGGTTTGAGCGCTGCCT-3′) and PtdaC-gfp-Rv (5′-TTCTCCTTTACGCATGGCCTGTTCCATCGTCC-3′). The promoter-less gfpmut3* open reading frame was amplified using the primer pair GFPmut3-Fw (5′-ACGATGGAACAGGCCATGCGTAAAGGAGAAGAACTTTTC-3′) and GFPmut3-Rv (5′-ATCCCGGGTTATTTGTATAGTTCATCCATGC -3′) and the PCR products were fused to the tdaCDE promoter segment via overlap PCR with primers PtdaB-Pinh and GFPmut3-Rv. The Q5 hot start high fidelity 2 × master mix (New England Biolabs, M0494S, BioNordika, Denmark) was used for the reactions. The fused PtdaCDE-gfp segments were cloned between the XhoI and XmaI sites of the pBBR1-MCS2 vector59 via restriction digestion cloning using the Quick Ligation Kit (New England Biolabs, M2200, BioNordika, Denmark), which resulted in plasmid pBBR1-MCS2-PtdaCDE-gfp. The recombined plasmid, as well as the pBBRMCS2 empty vector, were transferred into strain Escherichia coli WM3064, which was subsequently used as donors to conjugate the plasmids into the wild type DSM17395 strain60. This resulted in the strains DSM17395-GFP and DSM17395-NC.
Expression from the tdaCDE-promoter in an axenic model system
Expression from the tdaCDE promoter was determined in the pre-formed biofilm on the stainless steel coupons of either the DSM17395-GFP or the DSM17395-NC strains. Both strains were incubated for 3 days in 250 mL MB supplemented with Kan100 together with clean coupons in parallel beaker systems as described above. Subsequently, racks and coupons were transferred to 1-L beakers with 250 mL 2.4 % instant ocean (IO) and incubated for 10 days. Green fluorescence was detected from the attached P. inhibens cells after 0, 1, 4, and 10 days. Coupons were washed three times with 2.4 % IO to remove planktonic cells and placed on a cover slide. GFP-fluorescence images were captured using a 60× oil objective on a Nikon Inverted Fluorescence Microscope-EclipseTi2 (Nikon, Tokyo, Japan; WIB ex. 457–487 nm, em. 502–537 nm). The NIS-Element software was used for the acquisition of representative fluorescence micrographs with an exposure time of 170 ms. A total of 2 µL of a DAPI stain solution (of 100 ng/µL) was applied to the same specimen and fluorescence of DAPI (WIB ex. 352–402 nm, em. 417–477 nm) was recorded right after the staining. Images were processed using the Fiji package from ImageJ61 and any adjustments were made to the entire image and equally across samples.
DNA extractions
Cell pellets were lysed in 400 µL lysis buffer (400 mM NaCl, 750 mM sucrose, 20 mM EDTA, and 50 mM Tris-HCl [pH 8.3]) with a final concentration of 1 mg/mL lysozyme (Sigma, St. Louis, MO, USA). Cell lysates were incubated at 37 °C for 30 min, and sodium dodecyl sulfate (SDS) and proteinase K (Sigma, St. Louis, MO, USA) were added to final concentrations of 1 % and 0.1 mg/mL, respectively. After overnight incubation at 55 °C with regular agitation, samples were centrifuged at 3000 × g for 5 min and the supernatants were transferred to the lysate chamber in the Maxwell® 16 LEV Blood DNA Kit (Promega, Madison, WI, USA) and processed with the Maxwell™ 16 instrument. Finally, DNA was diluted into 50 µL TE elution buffer provided in the Maxwell® 16 LEV Blood DNA Kit.
Quantification of bacterial absolute abundances and P. inhibens
The absolute bacterial and P. inhibens abundances in the biofilm and the planktonic suspension samples were quantified by qPCR using the 16 S rRNA gene universal primers (338 F/518 R), as described in Bernbom et al 62, and Phaeobacter spp. specific primers (Pi_Fw/Pi_Rev) developed by Dittmann et al. 47. Standard curves relating CT to CFU/mL were based on gDNA extracted from five tenfold serial dilutions of Phaeobacter strain 27-421, Pseudoalteromonas tunicata D263, and Vibrio anguillarum 90-11-28764 for the estimation of total bacterial abundance and P. inhibens DSM17395 for the estimation of Phaeobacter spp. abundances. Each qPCR reaction contained 7.5 µL PowerUp™ SYBR® Green 2× Master Mix (ThermoFischer; A25742), 1.05 µL of forward (338 F; ACT CCT ACG GGA GGC AGC AG, Pi_Fw; 5′GTG TGT TGC GGT CTT TCA CC′3) and reverse (518 R; ATT ACC GCG GCT GCT GG, Pi_Rev; 5′AGG ACC ATG TCC CCT CTA CC′3) primer (10 µM), 4.4 μL DNAse free water and 1 μL DNA template. Standards and samples were run in duplicates on each plate. The MX3000P instrument (Stratagene, La Jolla, CA) was used with the following cycling conditions: 50 °C for 2 min, 95 °C for 2 min [95 °C for 15 s, 60 °C for 1 min] × 40 cycles.
Amplicon sequencing of the 16 S rRNA V3-V4 region
The V3-V4 region was amplified across all samples, including genomic DNA from P. inhibens DSM17395, using PCR-primers (Fw_V3-V4: 5′-CCTACGGGNGGCWGCAG-3′ and Rv_V3-V4: 5′-GACTACHVGGGTATCTAATCC-3′) tagged with octameric barcodes, amplified, cleaned and pooled in equimolar ratios as described in Bech et al. (2020)65. Amplicons were sequenced by Novogene, Cambridge, United Kingdom on the Illumina NovaSeq 6000 250PE platform.
V3-V4 amplicon processing with QIIME2 and DADA2
The forward and reverse sequences were processed in the QIIME266 pipeline, v2020.866. The 8 nt barcodes were removed using the qiime cutadapt demux-paired plugin. Later, paired-end reads were subjected to quality control followed by denoising, merging, and chimera removal using the DADA2 plugin implemented in QIIME 2 (dada2 denoise-paired with the following parameters trunc_len_f:242, trunc_len_r:242). The final amplicon sequence variant (ASV) table and the corresponding denoised ASVs were imported and analyzed in R. Denoised ASVs were individually taxonomically classified using the RDP naive Bayesian classifier67 trained with the SILVA v132 16 S rRNA gene sequence database68. A total of 50,223,017 reads passed quality control, with a mean ± sd of 250,017 ± 201,619 reads per sample, where the rarefaction curves of the ASV diversity was fully saturated across all samples (https://github.com/PKBech/PRJNA795304)). All ASVs were aggregated to genus level and low abundant genera (<100 reads across all samples) were removed from the dataset (0.0028 % of the total reads) leaving a total of 1056 genera in the final genus table (Table 1). Two ASVs were identified from the V3-V4 region amplicon sequencing of the genomic DNA of P. inhibens DSM17395. Since the V3-V4 region of this species is uniquely identifiable69, the sum of relative abundances of the two specific ASVs were added together as a single OTU referred to as the ‘P. inhibens OTU’ in the genus table.
Beta-diversity
Prior to beta-diversity analysis, the P. inhibens OTU was excluded from the genus table, normalized and scaled to a total sum of 100,000 reads per sample. All the analyses were performed using the ‘vegan’ R package v2.5 (see ‘URLs’ section). Bray-Curtis dissimilarities were calculated from the square-root transformed genus table used to determine the variation in the microbial community composition explained by the three categorical variables: time (day 1, 4, and 10), treatment (Control, WT- and dTDA systems) and environment (biofilm and planktonic). First, we verified that the beta-dispersion was similar and non-significantly different (p-value > 0.05) between the groups across all samples by the betadisper function. Next, the variance explained by the three independent variables was tested using a permutational multivariate analysis of variance (PERMANOVA) with the adonis function. As a posthoc test for pairwise differences within each variable, the pairwise.adonis2 function (see ‘URLs’ section) was used to specifically compare differences between the WT and the dTDA groups. The results were visualized by non-metric multidimensional scaling (nMDS) ordination plots using the metaMDS function with k = 3, since three dimensions achieved a stress value < 0.1.
Statistical analyses on alpha-diversity and absolute abundances
The observed richness was calculated by the estimateR function from the normalized genus table, including the P. inhibens OTU and analyzed with linear mixed-effects models (LMM) using the biological replicate as random effects with the lme4 package70. Differences in the marginal means (EMM) of the linear models of the observed richness and the absolute abundances of the total microbial population and Phaeobacter spp. quantified by qPCR at the individual days were estimated by the emmeans function in the ‘emmeans’ package (see ‘URLs’ section). Differences between seawater samples collected at day 0 and all other samples were tested by the Dunnest test. Means and standard deviations were summarized in Table 1.
Differential relative abundance analyses
We used the ANCOM-BC model, which is optimized to cope with and minimize the challenges associated with the nature of microbiome compositional data71, to identify differentially abundant genera and orders that responded to the pre-colonization of the Phaeobacter strains. The sum of differentially abundant genera was identified between the WT and the dTDA system by each day and environment. The non-normalized genus table excluding the P. inhibens OTU for each day and environment (N = 8–9) were analyzed with the ancombc function with default parameters and a cut-off of Holm adjusted p-values < 0.01 with the WT treated samples set as comparing reference group. The sum of differentially abundant genera found by each comparison were summarized in a balloonplot using the ggballoonplot function from the ggpubr R package (see ‘URLs’ section). To get the taxonomic profile of the log-fold differential changes by each comparison, all samples from day 4 were aggregated to order level and analyzed as described above for both biofilm and the planktonic suspension.
Co-occurrence networks
We constructed different networks to infer keystone taxa in the WT and the dTDA systems, using the NetCoMi package72 developed to compare and quantify network differences between co-occurrences of different genera and the P. inhibens OTU. A total of six comparison networks were constructed with the WT and the dTDA systems by day and environment. The SparCC method was used as the correlation method, which is developed for inferring correlations for compositional data and has been shown to be efficient in detecting co-occurrences in communities exhibiting relatively low Simpson diversity73. Genera having reads >100 in more than 5 replicates and a sparsity threshold of 0.6 were used to obtain the correlation matrix. The netAnalyze function was used for network characterization with the cluster_fast_greedy function to identify clusters and modularity. The scores of betweenness, degree, closeness, and eigenvector centrality have previously been used to capture different aspects of node importance and statistically identify keystone genera74,75,76. In this study, we defined hubs—or ‘keystone OTUs’ with the highest eigenvector centrality value above the empirical 95 % quantile72, combined by high degree, high closeness and low betweenness centrality, which have previously been suggested to be indications for keystone OTUs77. The plot.microNetProps function was used to visualize the network and quantitative network comparisons were done with the netCompare function using 1000 permutations.
URLs
Vegan package v2.5: https://github.com/vegandevs/vegan. Pairwise.adonis2 function: https://github.com/pmartinezarbizu/pairwiseAdonis. The emmeans function in the ‘emmeans’ package: https://github.com/rvlenth/emmeans. The ggballoonplot function from the ggpubr R package: https://rpkgs.datanovia.com/ggpubr/.
Data availability
Sequencing data is available in the Sequencing Read Archive (SRA), amplicon data in BioProject PRJNA795304 and R scripts at https://github.com/PKBech/PRJNA795304. The source data underlying Figs. 1–4 and Supplementary Figs. 2-3 are provided as a Source Data file.
References
Falkowski, P. G., Barber, R. T. & Smetacek, V. Biogeochemical controls and feedbacks on ocean primary production. Science (1979) 281, 200–206 (1998).
Nemergut, D. R. et al. Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356 (2013).
Lucas, J. et al. Annual dynamics of North Sea bacterioplankton: seasonal variability superimposes short-term variation. FEMS Microbiol. Ecol. 91, 1–11 (2015).
Bunse, C. & Pinhassi, J. Marine bacterioplankton seasonal succession dynamics. Trends Microbiol. 25, 494–505 (2017).
Romdhane, S. et al. Unraveling negative biotic interactions determining soil microbial community assembly and functioning. ISME J. 16, 296–306 (2021).
Tobias-Hünefeldt, S. P., Wenley, J., Baltar, F. & Morales, S. E. Ecological drivers switch from bottom–up to top–down during model microbial community successions. ISME J. 15, 1085–1097 (2021).
D’Onofrio, A. et al. Siderophores from neighboring organisms promote the growth of uncultured bacteria. Chem. Biol. 17, 254–264 (2010).
Abrudan, M. I. et al. Socially mediated induction and suppression of antibiosis during bacterial coexistence. Proc. Natl Acad. Sci. USA 112, 11054–11059 (2015).
Yamanaka, K. et al. Desferrioxamine E produced by Streptomyces griseus stimulates growth and development of Streptomyces tanashiensis. Microbiology 151, 2899–2905 (2005).
García-Bayona, L. & Comstock, L. E. Bacterial antagonism in host-associated microbial communities. Science (1979) 361, eaat2456 (2018).
Foster, K. R. & Bell, T. Report competition, not cooperation, dominates interactions among culturable microbial species. Curr. Biol. 22, 1845–1850 (2012).
Ghoul, M. & Mitri, S. The ecology and evolution of microbial competition. Trends Microbiol. 24, 833–845 (2016).
Kehe, J. et al. Positive interactions are common among culturable bacteria. Sci. Adv. 7, eabi7159 (2021).
Patin, N. V. et al. Effects of actinomycete secondary metabolites on sediment microbial communities. Appl. Environ. Microbiol. 83, 1–16 (2017).
Flemming, H. C. et al. Biofilms: an emergent form of bacterial life. Nat. Rev. Microbiol. 14, 563–575 (2016).
Stewart, P. S. & Franklin, M. J. Physiological heterogeneity in biofilms. Nat. Rev. Microbiol. 6, 199–210 (2008).
Dang, H., Li, T., Chen, M. & Huang, G. Cross-ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl. Environ. Microbiol. 74, 52–60 (2008).
Sonnenschein, E. C. et al. Global occurrence and heterogeneity of the Roseobacter-clade species Ruegeria mobilis. ISME J. 11, 569–583 (2017).
Gram, L. et al. Phaeobacter inhibens from the Roseobacter clade have a natural environmental niche as surface colonizer in Danish harbors. Syst. Appl. Microbiol. 38, 483–493 (2015).
Breider, S. et al. Phaeobacter porticola sp. Nov., an antibiotic-producing bacterium isolated from a sea harbour. Int J. Syst. Evol. Microbiol. 67, 2153–2159 (2017).
Hjelm, M. et al. Selection and identification of autochthonous potential probiotic bacteria from turbot larvae (Scophthalmus maximus) rearing units. Syst. Appl. Microbiol. 27, 360–371 (2004).
Porsby, C. H., Nielsen, K. F. & Gram, L. Phaeobacter and Ruegeria species of the Roseobacter clade colonize separate niches in a Danish turbot (Scophthalmus maximus)-rearing farm and antagonize Vibrio anguillarum under different growth conditions. Appl Environ. Microbiol 74, 7356–7364 (2008).
Rao, D. et al. Ability of Pseudoalteromonas tunicata to colonize natural biofilms and its effect on microbial community structure. FEMS Microbiol Ecol. 73, 450–457 (2010).
Brinkhoff, T. et al. Antibiotic Production by a Roseobacter Clade-Affiliated Species from the German Wadden Sea and Its Antagonistic Effects on Indigenous Isolates. Appl Environ. Microbiol 70, 2560–2565 (2004).
Bruhn, J. B. et al. Ecology, inhibitory activity, and morphogenesis of a marine antagonistic bacterium belonging to the Roseobacter clade. Appl Environ. Microbiol 71, 7263–7270 (2005).
Beyersmann, P. G. et al. Dual function of tropodithietic acid as antibiotic and signaling molecule in global gene regulation of the probiotic bacterium Phaeobacter inhibens. Sci. Rep. 7, 730 (2017).
D’Alvise, P. W., Phippen, C. B. W., Nielsen, K. F. & Gram, L. Influence of iron on production of the antibacterial compound tropodithietic acid and its noninhibitory analog in Phaeobacter inhibens. Appl Environ. Microbiol 82, 502–509 (2016).
Hernandez, D. J., David, A. S., Menges, E. S., Searcy, C. A. & Afkhami, M. E. Environmental stress destabilizes microbial networks. ISME J. 15, 1722–1734 (2021).
Grilli, J., Rogers, T. & Allesina, S. Modularity and stability in ecological communities. Nat. Commun. 7, 12031 (2016).
Coyte, K. Z., Schluter, J. & Foster, K. R. The ecology of the microbiome: networks, competition, and stability. Science (1979) 350, 663–6 (2015).
Kramer, J., Özkaya, Ö. & Kümmerli, R. Bacterial siderophores in community and host interactions. Nat. Rev. Microbiol 18, 152–163 (2020).
Davies, J. Specialized microbial metabolites: functions and origins. J. Antibiotics 66, 361–364 (2013).
Butaitė, E., Kramer, J. & Kümmerli, R. Local adaptation, geographical distance and phylogenetic relatedness: assessing the drivers of siderophore-mediated social interactions in natural bacterial communities. J. Evol. Biol. 34, 1266–1278 (2021).
Fajardo, A. & Martínez, J. L. Antibiotics as signals that trigger specific bacterial responses. Curr. Opin. Microbiol 11, 161–167 (2008).
Pishchany, G. & Kolter, R. On the possible ecological roles of antimicrobials. Mol. Microbiol. 113, 580–587 (2020).
Friedman, J., Higgins, L. M. & Gore, J. Community structure follows simple assembly rules in microbial microcosms. Nat. Ecol. Evol. 1, 109 (2017).
Goldford, J. E. et al. Emergent simplicity in microbial community assembly. Science (1979) 361, 469–474 (2018).
Majzoub, M. E., McElroy, K., Maczka, M., Thomas, T. & Egan, S. Causes and consequences of a variant strain of Phaeobacter inhibens with reduced competition. Front Microbiol. 9, 2601 (2018).
Mark, B. M. & Zhu, Y. Gliding motility and por secretion system genes are widespread among members of the phylum bacteroidetes. J. Bacteriol. 195, 270–278 (2013).
Berger, M., Neumann, A., Schulz, S., Simon, M. & Brinkhoff, T. Tropodithietic acid production in Phaeobacter gallaeciensis is regulated by N-acyl homoserine lactone-mediated quorum sensing. J. Bacteriol. 193, 6576–6585 (2011).
Porsby, C. H., Webber, M. A., Nielsen, K. F., Piddock, L. J. V. & Gram, L. Resistance and tolerance to tropodithietic acid, an antimicrobial in aquaculture, is hard to select. Antimicrob. Agents Chemother. 55, 1332–1337 (2011).
Henriksen, N. N. S. E. et al. The ability of Phaeobacter inhibens to produce tropodithietic acid influences the community dynamics of a microalgal microbiome. ISME Commun. 2, 109 (2022).
Dittmann, K. K. et al. Changes in the microbiome of mariculture feed organisms after treatment with a potentially probiotic strain of Phaeobacter inhibens. Appl. Environ. Microbiol 86, 1–15 (2020).
Gram, L. et al. Phaeobacter inhibens from the Roseobacter clade has an environmental niche as a surface colonizer in harbors. Syst. Appl. Microbiol. 38, 483–493 (2015).
Cottrell, M. T. & Kirchman, D. L. Natural assemblages of marine proteobacteria and members of the Cytophaga-flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl Environ. Microbiol 66, 692–697 (2000).
Teeling, H. et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science (1979) 336, 608–611 (2012).
Dittmann, K. K., Sonnenschein, E. C., Egan, S., Gram, L. & Bentzon-Tilia, M. Impact of Phaeobacter inhibens on marine eukaryote-associated microbial communities. Environ. Microbiol Rep. 11, 401–413 (2019).
Geng, H., Tran-Gyamfi, M. B., Lane, T. W., Sale, K. L. & Yu, E. T. Changes in the structure of the microbial community associated with Nannochloropsis salina following treatments with antibiotics and bioactive compounds. Front Microbiol 7, 1155 (2016).
Wilson, M. Z., Wang, R., Gitai, Z. & Seyedsayamdost, M. R. Mode of action and resistance studies unveil new roles for tropodithietic acid as an anticancer agent and the γ-glutamyl cycle as a proton sink. Proc. Natl Acad. Sci. USA 113, 1630–1635 (2016).
D’Alvise, P. W. et al. Phaeobacter gallaeciensis reduces vibrio anguillarum in cultures of microalgae and rotifers, and prevents vibriosis in cod larvae. PLoS One 7, e43996 (2012).
Harrington, C. et al. Characterisation of non-autoinducing tropodithietic acid (TDA) production from marine sponge Pseudovibrio species. Mar. Drugs 12, 5960–5978 (2014).
Tesdorpf, J. E., Geers, A. U., Strube, M. L., Gram, L. & Bentzon-Tilia, M. Roseobacter group probiotics exhibit differential killing of fish pathogenic Tenacibaculum Species. Appl. Environ. Microbiol. 88, e0241821 (2022).
Banerjee, S., Schlaeppi, K. & van der Heijden, M. G. A. Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol 16, 567–576 (2018).
Martens, T. et al. Reclassification of Roseobacter gallaeciensis Ruiz-Ponte et al. 1998 as Phaeobacter gallaeciensis gen. nov., comb. nov., description of Phaeobacter inhibens sp. nov., reclassification of Ruegeria algicola (Lafay et al. 1995) Uchino et al. 1999 as Marinovum algicola gen. nov., comb. nov., and emended descriptions of the genera Roseobacter, Ruegeria and Leisingera. Int J. Syst. Evol. Microbiol. 56, 1293–1304 (2006).
Wang, R., Gallant, É. & Seyedsayamdost, M. R. Investigation of the genetics and biochemistry of roseobacticide production in the Roseobacter clade bacterium Phaeobacter inhibens. mBio 7, e02118 (2016).
Bernbom, N. et al. Adhesion of food-borne bacteria to stainless steel is reduced by food conditioning films. J. Appl. Microbiol. 106, 1268–1279 (2009).
Noble, R. T. & Fuhrman, J. A. Use of SYBR Green I for rapid epifluorescence counts of marine viruses and bacteria. Aquat. Microb. Ecol. 14, 113–118 (1998).
Asséré, A., Oulahal, N. & Carpentier, B. Comparative evaluation of methods for counting surviving biofilm cells adhering to a polyvinyl chloride surface exposed to chlorine or drying. J. Appl. Microbiol. 104, 1692–1702 (2008).
Kovach, M. E. et al. Four new derivatives of the broad-host-range cloning vector pBBR1MCS, carrying different antibiotic-resistance cassettes. Gene 166, 175–176 (1995).
Zhang, S. D., Isbrandt, T., Lindqvist, L. L., Larsen, T. O. & Gram, L. Holomycin, an Antibiotic Secondary Metabolite, Is Required for Biofilm Formation by the Native Producer Photobacterium galatheae S2753. Appl Environ. Microbiol. 87, e00169–21 (2021).
Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).
Bernbom, N., Ng, Y. Y., Olsen, S. M. & Gram, L. Pseudoalteromonas spp. Serve as initial bacterial attractants in mesocosms of coastal waters but have subsequent antifouling capacity in mesocosms and when embedded in paint. Appl Environ. Microbiol. 79, 6885–6893 (2013).
Holmström, C., James, S., Neilan, B. A., White, D. C. & Kjelleberg, S. Pseudoalteromonas tunicata sp. nov., a bacterium that produces antifouling agents. Int J. Syst. Bacteriol. 48, 1205–1212 (1998).
Skov, M. N., Pedersen, K. & Larsen, J. L. Comparison of pulsed-field gel electrophoresis, ribotyping, and plasmid profiling for typing of Vibrio anguillarum serovar O1. Appl Environ. Microbiol. 61, 1540–1545 (1995).
Bech, P. K., Lysdal, K. L., Gram, L., Bentzon-Tilia, M. & Strube, M. L. Marine sediments hold an untapped potential for novel taxonomic and bioactive bacterial diversity. mSystems 5, e00782–20 (2020).
Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).
Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ. Microbiol. 73, 5261–5267 (2007).
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).
Strube, M. L. RibDif: can individual species be differentiated by 16S sequencing? Bioinf. Adv. 1, vbab020 (2021).
Bates, D., Mächler, M., Bolker, B. M. & Walker, S. C. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Lin, H. & Peddada, S. D. Analysis of compositions of microbiomes with bias correction. Nat. Commun. 11, 3514 (2020).
Peschel, S., Müller, C. L., von Mutius, E., Boulesteix, A. L. & Depner, M. NetCoMi: network construction and comparison for microbiome data in R. Brief. Bioinf. 22, bbaa290 (2021).
Weiss, S. et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 10, 1669–1681 (2016).
Liang, Y. et al. Long-term oil contamination alters the molecular ecological networks of soil microbial functional genes. Front. Microbiol. 7, 60 (2016).
Li, F., Chen, L., Zhang, J., Yin, J. & Huang, S. Bacterial community structure after long-term organic and inorganic fertilization reveals important associations between soil nutrients and specific taxa involved in nutrient transformations. Front. Microbiol. 8, 187 (2017).
Costello, E. K., Stagaman, K., Dethlefsen, L., Bohannan, B. J. M. & Relman, D. A. The application of ecological theory toward an understanding of the human microbiome. Science (1979) 336, 1255–1262 (2012).
Berry, D. & Widder, S. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front Microbiol. 5, 219 (2014).
Henriksen, N. N. S. E. et al. Role is in the eye of the beholder—the multiple functions of the antibacterial compound tropodithietic acid produced by marine Rhodobacteraceae. FEMS Microbiol Rev. 46, fuac007 (2022).
Acknowledgements
We thank Dr. Mohammed Seyedsayamdost for kindly sharing the TDA-deficient mutant of Phaeobacter inhibens. The study was supported by the Danish National Research Foundation (DNRF137) for the Center for Microbial Secondary Metabolites and the Villum Kann Rasmussen Annual Award 2016 to Lone Gram.
Author information
Authors and Affiliations
Contributions
Conceived and designed the experiments: P.K.B., L.G., M.B.T. Performed the experiments: P.K.B. Analysed the data: P.K.B. Construction of the GFP strains: S.D.Z. Contributed to the analyses: N.S.H. and M.L.S. Wrote the paper: P.K.B., L.G.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bech, P.K., Zhang, SD., Henriksen, N.N.S.E. et al. The potential to produce tropodithietic acid by Phaeobacter inhibens affects the assembly of microbial biofilm communities in natural seawater. npj Biofilms Microbiomes 9, 12 (2023). https://doi.org/10.1038/s41522-023-00379-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41522-023-00379-3
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.