- 1Laboratory of Protistology and Aquatic Ecology, Department of Biology, Ghent University, Ghent, Belgium
- 2Laboratory of Microbiology, Department of Biochemistry and Microbiology, Ghent University, Ghent, Belgium
- 3SynBioC Research Group, Department of Green Chemistry and Technology, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
- 4Marine Biology, Department of Biology, Ghent University, Ghent, Belgium
Benthic diatoms are dominant primary producers in intertidal mudflats and constitute a major source of organic carbon to consumers and decomposers residing within these ecosystems. They typically form biofilms whose species richness, community composition and productivity can vary in response to environmental drivers and their interactions with other organisms (e.g., grazers). Here, we investigated whether bacteria can affect diatom community composition and vice versa, and how this could influence the biodiversity-productivity relation. Using axenic experimental communities with three common benthic diatoms (Cylindrotheca closterium, Navicula phyllepta, and Seminavis robusta), we observed an increase in algal biomass production in diatom co-cultures in comparison to monocultures. The presence of bacteria decreased the productivity of diatom monocultures while bacteria did not seem to affect the overall productivity of diatoms grown in co-cultures. The effect of bacteria on diatom growth, however, appeared to be species-specific, resulting in compositional shifts when different diatom species were grown together. The effect of the diatoms on the bacteria also proved to be species-specific as each diatom species developed a bacterial community that differed in its composition. Together, our results suggest that interactions between bacteria and diatoms residing in mudflats are a key factor in the structuring of the benthic microbial community composition and the overall functioning of that community.
Introduction
Intertidal mudflats, found predominantly along estuaries and sea inlets, are highly productive ecosystems (Underwood and Kromkamp, 1999). Their productivity is in part due to benthic microalgal biofilms found on the surface sediments (Decho, 2000; Hochard et al., 2010), with benthic diatoms often being one of the dominant primary producers (Admiraal et al., 1984; MacIntyre et al., 1996; Underwood and Kromkamp, 1999; Bolhuis et al., 2013). These biofilms strongly modulate nutrient fluxes (Cook et al., 2009; Hochard et al., 2010) and provide copious amounts of autochthonously fixed carbon to successive trophic levels within the mudflat (Moerdijk-Poortvliet et al., 2018). The productivity of benthic biofilms, dominated by diatoms, depends on the diatom community composition and structure (Brotas and Plante-Cuny, 2003; Underwood et al., 2005) which in turn is dependent on a combination of both abiotic and biotic factors. Although the variable distribution of diatoms along estuarine gradients can partially be explained by the species-specific tolerances of diatom species for environmental conditions such as salinity and nutrient concentrations (Thornton et al., 2002; Forster et al., 2006; Ribeiro et al., 2013; Sawai et al., 2016; Desianti et al., 2017), biotic factors, such as competition and niche differentiation between different diatom species (Vanelslander et al., 2009) as well as the presence of bacteria (D’Costa and Anil, 2011), can also affect the diatom community structure and composition. There is as yet, however, no overall consensus on the exact relation between primary productivity and algal biodiversity of intertidal benthic biofilms. The biodiversity-productivity relation was shown to vary across sites (Forster et al., 2006) and could be either negative or positive (Forster et al., 2006; Vanelslander et al., 2009). Negative biodiversity effects on productivity are often attributed to competition or chemical interference (i.e., allelopathy) while positive biodiversity effects are the result of selection and/or complementarity effects (Loreau and Hector, 2001). A positive selection effect occurs due to the dominance of highly productive species while a complementarity effect increases productivity by enhancing the use of the available resources as a result of positive interactions or niche differentiation. This study investigates the various interactions taking place within a benthic diatom community and whether the presence of bacteria impacts the community regarding species composition and productivity and vice versa.
Due to the confined space within the matrix of marine benthic biofilms, interactions between the diatoms, and other organisms become inevitable (Clarke, 2016). The excretion of exopolymeric substances (EPS) by diatoms, for example, serves as an important carbon source that can be utilized by the residing benthic organisms (Middelburg et al., 2000; Bellinger et al., 2009) including heterotrophic bacteria (Middelburg et al., 2000; Taylor et al., 2013; Durham et al., 2015) as well as other diatom species (Vanelslander et al., 2009). Since the composition of diatom exudates is both species-specific and dependent on environmental conditions (Durham et al., 2015; Bohórquez et al., 2017), EPS has been suggested to play a major role in determining bacterial community composition and diversity (Haynes et al., 2007; Agogue et al., 2014; Wear et al., 2015; Mühlenbruch et al., 2018). Indeed, studies have shown that diatom species can harbor different associated bacteria (Schafer et al., 2002; Grossart et al., 2005; Doiron et al., 2012; Behringer et al., 2018) and that bacterial community composition in intertidal mudflats strongly co-varies with the composition of the microphytobenthos (Bolhuis et al., 2013; Decleyre et al., 2015; Lavergne et al., 2017).
In turn, the associated bacteria have also been shown to influence diatom growth rates as well as other life cycle features such as sexual reproduction (Grossart, 1999; Amin et al., 2012; Cirri et al., 2018). Despite the general view of heterotrophic bacteria as the primary remineralisers of (e.g., diatom-derived) organic matter (Azam, 1998; Taylor et al., 2013), releasing nutrients back into the environment, the bacterial community can also compete with diatoms for limited nutrients and resources (Thingstad et al., 1993; Grossart, 1999; Havskum et al., 2003; Amin et al., 2012). Other, more specific, interactions have also been observed, such as the ability of strains to produce vitamins (i.e., cobalamin) for auxotrophic diatoms (Haines and Guillard, 1974; Croft et al., 2005), as well as negative interactions such as the production of algicidal metabolites and the induction of diatom cell lysis (Furusawa et al., 2003; Paul and Pohnert, 2011). These effects can be highly specific, with some bacteria having the ability to stimulate the growth of one diatom whilst inhibiting the growth of another (Grossart, 1999; Jung et al., 2008; Paul and Pohnert, 2011; Sison-Mangus et al., 2014).
As many of the above-mentioned studies focused on one-to-one interactions between diatoms and bacteria, they do not reflect the ecological complexity found in nature. Bigalke and Pohnert (2019), for example, recently illustrated that a more complex diatom community reacts differently to bacteria than the individual species did. While biotic and abiotic factors could promote the presence of a potentially harmful bacterial community (Marzinelli et al., 2018) the host-associated bacterial community can also play a role in its resistance and tolerance toward a new or stressful environment (Dittami et al., 2016). Studies like these, stress that the oversimplification of experiments can prevent further insight into relevant ecological interactions. Whether a bacterial community can influence community functioning in terms of productivity and the outcome of interspecific interactions amongst co-occurring biofilm diatom species has to date not been tested. Although D’Costa and Anil (2011) have shown that diatom community structure is related to bacterial community composition, it was not determined whether and how diatom-bacteria associations could affect changes in diatom productivity and interactions. Recent work on plants and macroalgae, suggests that feedback mechanisms between hosts and their specific bacterial community can ultimately determine the outcome of competition amongst host species (Dittami et al., 2016; Gribben et al., 2017; Hortal et al., 2017; Lekberg et al., 2018; Marzinelli et al., 2018).
Using different benthic diatom species grown at three diversity levels (monoculture vs. 2- and 3-species co-cultures) in the presence or absence of a natural bacterial inoculum, we investigated whether bacteria can change the structure and productivity of simple diatom communities. We hypothesize these interkingdom interactions are species-specific and that certain bacteria promote the growth of selected diatoms and vice versa. We therefore predict that the interactions between diatoms and bacteria will change the outcome of competition amongst the different diatom species, favoring the diatom species that benefit the most from the presence of bacteria. This, in turn, could further influence diatom productivity. Finally, we tested whether diatoms can also affect bacterial diversity, predicting an increase in bacterial diversity with increasing diatom diversity and productivity.
Materials and Methods
Experimental Setup
The marine benthic diatoms Cylindrotheca closterium, Navicula phyllepta, and Seminavis robusta were obtained from the BCCM-DCG culture collection (Supplementary Table S1). The three strains were originally not isolated from the same location, but these three diatom species have been observed to co-occur, making them ecologically relevant for this study. Diatoms were cultivated in artificial sea water (ASW; Tropic Marin Bio-Actif Salt) that was enriched with 0.08 g/L NaHCO3 and Guillard’s Marine Water Enrichment Solution (F/2; Sigma) according to the manufacturer’s instructions. Throughout the experiment, cultures were grown in 12:12 h light-dark cycles using cool fluorescent white light (20–25 μmol photons/s/m2) and maintained at 18°C.
Diatoms were made axenic (free of bacteria) by subjecting them to repeated antibiotic treatments. Every 3 days, the medium of the diatoms was refreshed and the next antibiotic treatment was given to the cultures. This process was repeated at least three times. Antibiotic mixes consisted of final concentrations of 100 μg/mL gentamicin, 500 μg/mL streptomycin, and 100 μg/mL neomycin for C. closterium; and 500 μg/mL penicillin, 500 μg/mL ampicillin, 100 μg/mL streptomycin, and 50 μg/mL gentamicin (Sigma-Aldrich) for N. phyllepta and S. robusta. DAPI staining (Shishlyannikov et al., 2011) and plating on Difco Marine agar (BD) was used to confirm axenicity of the diatom cultures after the treatments. Although this study did not specifically test whether the antibiotic treatments resulted in a decline in health or a decreased growth rate for the individual diatom species used within this experiment it was done in the framework of another study (Stock et al., unpublished).
Prior to the start of the experiment, antibiotics in the diatom cultures were washed away by repeated (three times) refreshing of medium with new ASW 3 days after the last antibiotic treatment. The refreshed cultures were diluted to approximately 300, 100, and 85 cells/mL for S. robusta, C. closterium, and N. phyllepta, respectively. Dilutions were based on prior experiments (data not shown), which showed that the coexistence between the diatom species required higher starting densities of S. robusta. At equal starting densities, S. robusta was rapidly outcompeted due to its lower growth rate.
The bacterial inoculum was obtained from intertidal surface mud collected at the Paulina polder, Westerschelde, NL (2 March, 2016 at 51°21′032′′N, 3°43′574′′ E). This sample was taken 1 day before setting up the experiment. During this time, it was stored at 4°C. On the day of the experiment, 20 mL of ASW was added to ± 10 mL of sediment, vortexed and then filtered (3 μm filter size) to separate the bacteria from larger eukaryotic organisms and particles. 1 mL of the bacterial suspension was frozen (-20°C) for DNA extraction (see below).
The experiments were run in 24-well plates (Greiner Bio-one) with the addition of either 1 mL of bacterial inoculum (non-axenic) or 1 mL of ASW (axenic). 1 mL of diatom suspension was then added to obtain a total volume of 2 mL per well. Specifically, 1 mL of a single diatom, 2 × 0.5 mL of two diatoms or 3 × 0.333 mL of all three diatoms were added to obtain monocultures and 2- and 3-species co-cultures, respectively. F/2 stock solution was supplemented to all wells so that a final nutrient concentration of F/20 (one tenth of F/2) was reached. All different combinations were run in triplicate (n = 3), randomized across each plate and grown in the conditions stated above.
Daily biomass estimates were obtained using both pulse amplitude modulation (PAM) fluorometry (IMAGING-PAM M-Series Maxi version, Walz) and cell counts (Stock et al., 2019). Minimal fluorescence (F0) was measured (intensity: 6, gain: 3, and frequency: 1) daily, for a period of 7 days. Photographs were taken of each well in triplicate (Nikon Elemental Camera DS-Fi2 with 10x magnification) using an inverted Axiovert 135 Zeiss microscope. Diatoms were manually counted from the pictures (Supplementary Figure S1) and averaged to obtain daily diatom cell counts.
After 1 week, 1 mL of medium was extracted from every well and stored at -20°C for DNA analysis. DNA extractions were conducted according to Muyzer et al. (1993). The library prep and amplification of the 16S rRNA gene, using primers PA_ill and BKL1_III, were done according to D’Hondt et al. (2018). An artificially created mock community was included to benchmark processing variables (Tytgat et al., 2016; Supplementary Table S2).
Data Analysis and Statistics
Diatoms
Growth curves were constructed from diatom cell counts and PAM fluorometry measurements (F0) using Excel (Version 15.32). The maximal growth rate (μmax) was defined as the highest observed growth rate during the exponential growth phase (the first 4 days as determined after inspection of each individual growth curve) of the experiment and calculated from the slope (LINEST function) after a log2 transformation for every three consecutive days using a moving window. This strategy was applied to both F0 and diatom cell counts. A two-way ANOVA (R-version 3.4.1) was used to test for differences in μmax between diatom combinations and in the presence or absence of bacteria. The assumptions of ANOVA regarding homogeneity of variance and normality were, respectively verified through a Levene and Shapiro–Wilk Test. Where appropriate, a Tukey HSD test was applied as post hoc test.
Diatom productivity was defined as the increase in the diatom biovolume and calculated by multiplying the μmax obtained from the cell counts with the cell biovolume calculated according to Hillebrand et al. (1999). Productivity was expressed in function of diatom biodiversity. Biodiversity effects were further partitioned into selection and complementarity effects according to Loreau and Hector (2001) in R (version 3.4.1). The partitioning of biodiversity effects was done separately for both axenic and non-axenic cultures.
Bacteria
The 300 bp pair-end MiSeq reads were joined and quality-filtered using PEAR (Version 0.9.5). After primer removal, an operational taxonomic unit (OTU) table was constructed by clustering at a 3% divergence level (USearch8). Taxonomic assignment of OTUs was done using MOTHER (Version 1.35.1). Chimeras were removed using the internal check in Usearch whilst mitochondrial and chloroplast reads were removed based on taxonomic assignment. The OTU table was further processed in R (version 3.4.1). To remove potential cross-contamination and sequencing errors, read counts <5 and OTUs that were present in less than three samples were set to 0 (based on 16S mock community data). Relative abundances were derived from the constructed OTU table. Rare OTUs (<1% of number of reads/sample) were removed. A PERMANOVA was run using the Adonis function (Vegan Package 2.4-3) between all diatom combinations and treatments. Significance values were based on 1000 permutations. A constrained correspondence analysis (ANOVA-CCA, Vegan Package 2.4-3) was performed with different diatom combinations set as Boolean variables to identify the influence of the different diatom species on the bacterial community composition. Finally, a Simper Analysis (10000 permutations, Vegan Package 2.4-3) was run to identify the discriminating bacterial OTUs between the diatom monocultures. Bacterial diversity, expressed as both the predicted number of OTUs when rarefying all samples to the same depth and the Simpson’s index, was also linearly expressed in function of both diatom productivity and diatom-diversity. This procedure was separately repeated for the most diverse bacterial group, namely the Alphaproteobacteria.
Results
Diatom Growth
In axenic conditions, C. closterium and N. phyllepta displayed comparably high μmax in monocultures while S. robusta grew considerably slower (Figure 1 and Table 1). The presence of other diatoms did not influence the μmax of C. closterium (p > 0.05) and N. phyllepta (p > 0.05) but significantly increased (p < 0.001) the μmax of S. robusta (Figure 1A). The addition of a mixed natural bacterial inoculum to the monocultures significantly lowered the μmax of C. closterium (p = 0.011) and N. phyllepta (p = 0.007), but not of S. robusta (p > 0.05) (Figure 1A). In the non-axenic co-cultures, the negative effect of bacteria on the growth of N. phyllepta and C. closterium seemed slightly alleviated by the presence of other diatoms (albeit not significantly), although the total cell numbers never exceeded those of their axenic counterparts (Figure 1B). The growth rate of S. robusta was again significantly higher (p ≤ 0.003) in mixed diatom cultures and in the presence of the two other diatom species, S. robusta was the fastest grower (Figure 1A). These species-specific changes in growth rates due to the presence of bacteria resulted in compositional shifts in the diatom co-cultures: the relative abundance of S. robusta was higher when bacteria were present, mainly at the expense of C. closterium (Supplementary Figures S3, S4). Similar trends were observed for the μmax derived from daily fluorescence measurements (Supplementary Figure S5 and Supplementary Table S3) although PAM measurements did not allow separation of biomass increase between diatom species grown in co-cultures. Altogether, changes in μmax, from both cell counts (Figure 1A) or PAM measurements (Supplementary Figure S5) as well as the total biovolume (Supplementary Figure S4) indicate that, while the growth rates of C. closterium and N. phyllepta were significantly repressed in the presence of bacteria, S. robusta was left unaffected while its μmax significantly improved in the presence of other diatom species.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/454984/fmicb-10-01255-HTML/image_m/fmicb-10-01255-g001.jpg)
Figure 1. Bacteria cause a shift in diatom growth rates. (A) Maximum growth rates, derived from daily cell counts for three different diatom species: C. closterium ([C], red), N. phyllepta ([N], blue), and S. robusta ([S], yellow) in the presence (non-axenic) or absence (axenic) of bacteria. The white, gray and dark gray background, respectively show diatoms grown either in monocultures, diatom pairs or the three diatom-species combined. (B) The number of diatom cells per well on the fourth and last day are shown using the same color code as above.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/454984/fmicb-10-01255-HTML/image_m/fmicb-10-01255-t001.jpg)
Table 1. Statistical analysis of the maximum growth rates derived from individual cell counts for S. robusta [S], C. closterium [C], and N. phyllepta [N].
Additive Partitioning of Biodiversity Effects
We determined whether increasing diatom species richness had an impact on productivity, i.e., the algal biomass production, and whether this would change in the presence of bacteria (Figure 2). In general, productivity increased with increasing diatom diversity (p < 0.0001) and the presence of bacteria steepened this relationship (p = 0.007; Figure 2A). This effect of bacteria can largely be attributed to a strong decrease in productivity in the diatom monocultures when exposed to the bacterial inoculum (p = 0.045). In contrast, the bacteria have a reduced negative effect on the diatom species grown in the presence of other diatoms.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/454984/fmicb-10-01255-HTML/image_m/fmicb-10-01255-g002.jpg)
Figure 2. The presence of bacteria steepens the algal diversity-productivity relation. (A) Diatom biovolume production in function of species richness and its slope (m) in the presence (non-axenic, black line) or absence of bacteria (axenic, gray line). (B) Algal biovolume production per diatom species and its slope for each diatom: S. robusta (yellow), C. closterium (red), and N. phyllepta (blue). The presence or absence of bacteria are, respectively depicted as dark- and light-colored lines. The contribution of the selection and complementarity effects to diversity-productivity relation are, respectively shown in (C,D) in the presence (black) or absence (gray) of bacteria.
Analogous to the results obtained from the diatom growth rates, the productivity of S. robusta strongly increased in the presence of other diatoms under axenic conditions, while both C. closterium and N. phyllepta had lower productivity in the co-cultures (Figure 2B). Adding bacteria switched the observed negative diversity effect in C. closterium and N. phyllepta to a positive diversity effect, with both species now showing higher productivity in the co-cultures. This effect, however, was not very strong, and was always less pronounced in comparison to S. robusta (Figure 2B).
The effect of diatom diversity on the productivity was further partitioned into selection and complementarity effects (Loreau and Hector, 2001; Figures 2C,D). The observed increase in production (p < 0.0001) can largely be attributed to positive complementarity effects, which further increased in the presence of bacteria (p > 0.05). A positive but very minor selection effect on productivity was also observed in the presence of bacteria (p > 0.05).
Bacteria Community Composition
After removing non-bacterial and potentially chimeric reads from the 16S rDNA high-throughput sequencing, 8724 ± 5505 (mean ± SD) reads per sample remained. These were assigned to 123 OTUs. Based on the rarefaction curves (Supplementary Figure S2) we determined that the sequencing depth was sufficient to detect the most abundant bacteria present in the algal cultures. The composition of the initial bacterial inoculum was dominated by Rhodospirillaceae (41.5% of the inoculum reads) (Figure 3). Following co-cultivation with diatoms over a week, the bacterial communities became dominated by Rhodobacteraceae, Flavobacteriaceae, and Campylobacteraceae both in terms of diversity (32, 18, and 8%, respectively) and relative abundance (68, 6, and 14%, respectively). Multivariate analyses of the data (Supplementary Figure S6) showed a strong differentiation (PERMANOVA: p = 0.006) between the bacterial communities associated with the different diatom combinations. The strongest difference was noticeable between the S. robusta monocultures and the cultures containing C. closterium and N. phyllepta (ANOVA-CCA: p = 0.01), although there were some minor differences between the latter bacterial communities as well (ANOVA-CCA: p > 0.05). The bacterial communities associated with S. robusta monocultures had a much higher diversity (H’ = 2.01 ± 0.15 based on the rarefied data) in comparison with the bacterial communities associated with the other diatoms (H’ = 1.11 ± 0.21 and H’ = 0.87 ± 0.09 for C. closterium and N. phyllepta monocultures, respectively).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/454984/fmicb-10-01255-HTML/image_m/fmicb-10-01255-g003.jpg)
Figure 3. Bacterial community composition follows the diatom community. Relative bacterial abundances present in the original inoculum (Ino) or when grown with either C. closterium [C], N. phyllepta [N], S. robusta [S], or a combination of these diatoms. The different bacterial orders are depicted using different colors.
Bacterial communities formed in the presence of S. robusta typically contained (amongst others) Alphaproteobacteria (Thalassospira sp., Roseobacter sp., and a Kordiimonadaceae sp.) and Bacteroidetes (Mangrovimonas sp. and Owenweeksia sp.). Especially the Mangrovimonas sp. tended to be highly abundant in S. robusta cultures (up to 42%) whilst being almost absent in C. closterium and N. phyllepta monocultures. In contrast, the bacterial communities of the latter monocultures were dominated (up to 80% of relative abundances) by an Octadecabacter species (Alphaproteobacteria, Rhodobacteraceae) and to a lesser degree also by an unidentified Arcobacter sp. (Epsilonproteobacteria). Bacterial diversity in the diatom co-cultures was not higher than in the monocultures (Supplementary Figure S7). Instead, an intermediate community structure was observed which was more similar to bacterial community present in the N. phyllepta and/or C. closterium monocultures. Finally, no linear relation (p > 0.05) was found between bacterial diversity, as a whole or for the Alphaproteobacteria separately, and diatom productivity.
Discussion
Our results show a positive relationship between benthic diatom diversity and productivity in simple experimental diatom communities. Although the biodiversity-productivity relation has been shown to vary for diatoms at varying spatial scales and conditions (Gessner et al., 2004; Giller et al., 2004; Soininen, 2009; Blanco et al., 2012; Baert et al., 2016), our results are in accordance with what was observed previously for marine benthic diatoms (Forster et al., 2006; Vanelslander et al., 2009). According to the biodiversity-productivity relation, the positive effect of diversity on productivity could largely be attributed to complementarity effects (Loreau and Hector, 2001). This indicated that the competition between diatom species was generally smaller than competition within diatom species and was likely the result of the improved usage of the available resources between the diatom species through cross-feeding and facilitation events (Fridley, 2001; Bruno et al., 2003). Although our experimental setup did not allow further analysis of EPS production and utilization, the presence of chemical cross-talk between diatoms, which include changes in the production, release and utilization of organic compounds, seems to be relatively common (Tuchman et al., 2006) and has been shown to stimulate the growth of co-cultured diatoms (Paul et al., 2009; Vanelslander et al., 2009). For example, the spent medium derived from N. phyllepta, resulted in a mixotrophic switch that improved the growth rate of a specific C. closterium strain, by benefitting from the carbon released by N. phyllepta (Vanelslander et al., 2009). The capability of such a mixotrophic lifestyle, however, appears to be strain dependent (Mensens et al., 2015) and was not observed for the C. closterium strain used in this study (Audoor unpublished data). Regarding our results, a similar process, however, may be occurring in S. robusta, whose growth was positively affected by the presence of other diatom species and would require further testing.
Our experimental setup revealed the recruitment of distinct bacterial assemblages from a common bacterial inoculum by individual diatom species, which is in accordance with previous studies (Schafer et al., 2002; Grossart et al., 2005; Jasti et al., 2005; Eigemann et al., 2013; Bagatini et al., 2014; Sison-Mangus et al., 2014). Interestingly, despite an expected increase in substrate and habitat heterogeneity (Kassen et al., 2000; Giller et al., 2004), the combination of different diatom species did not result in a higher bacterial diversity. Neither could a relation be found between primary production and bacterial diversity, either as a whole or within certain taxonomic groups (Horner-Devine et al., 2003; Stein et al., 2014; Raes et al., 2018). Instead, the bacterial community composition of a mixture of diatom species reflected an intermediate combination of the bacteria that were present in the three separate diatom monocultures.
Although S. robusta and N. phyllepta are more closely related to one another, the bacterial community composition of N. phyllepta was more similar to that of C. closterium. S. robusta harbored a diverse bacterial community in comparison to the bacterial communities of N. phyllepta and C. closterium which were both dominated by Octadecabacter sp., a member of the Roseobacter clade, and Arcobacter sp. The dominance of both bacteria over the other bacteria, coinciding with a negative impact on the growth rate of two out of the three diatom species, is a possible indicator of a negative and potentially allelopathic effect from these bacteria on other microorganisms (Mayali et al., 2008; Slightom and Buchan, 2009). A literature search revealed that not much is known about marine Arcobacter species (Collado and Figueras, 2011), but the genus does include several important mammalian pathogens (Ferreira et al., 2015). Several Roseobacter representatives, on the other hand, are known to reduce both algal and bacterial growth (Mayali et al., 2008; Seyedsayamdost et al., 2011; Sharifah and Eguchi, 2011), e.g., through the production of antibiotics (Ruiz-Ponte et al., 1999; Brinkhoff et al., 2004; Wagner-Dobler et al., 2004; Bruhn et al., 2007). Interestingly, the reduced manifestation of Octadecabacter and Arcobacter in the presence of S. robusta, whose growth rate was left unaffected in the presence of bacteria, could further suggest the ability of either S. robusta or its associated microbial community to suppress the growth of these bacteria or their potential antibiotic activity.
Although several studies have highlighted the species-specific effects of bacteria on diatom growth (e.g., Grossart, 1999; Amin et al., 2012; Doiron et al., 2012) and diatom community composition (D’Costa and Anil, 2011), our results indicate how these diatom-bacteria interactions may further impact the outcome of diatom-diatom competition and the subsequent effect on productivity. Indeed, few studies have investigated the impact of bacteria-host interactions on community functioning and results appear to differ (Hubbard et al., 1986; Grime et al., 1987; Horner-Devine et al., 2003; Pugnaire et al., 2004; Kardol et al., 2007; Gribben et al., 2017; Hortal et al., 2017; Laforest-Lapointe et al., 2017) and are dependent on the environmental conditions (Callaway et al., 2004; Dittami et al., 2016; Marzinelli et al., 2018). In the absence of bacteria, experimental diatom communities were dominated by the species that showed the highest growth rate in monoculture, i.e., C. closterium. In mixed diatom cultures, the addition of bacteria negatively impacted the growth of both C. closterium and N. phyllepta, leading to the dominance of S. robusta. The presence of a bacterial community further altered the algal diversity functioning relationship, steepening the positive relationship found between diatom diversity and productivity as a result of an enhanced complementarity effect. Similar to what has been observed in plants (Lekberg et al., 2018), the most competitive diatoms in our study experienced the most negative bacterial effects which in turn resulted in the bacteria indirectly promoting co-existence amongst the diatom species. Even if our simple experimental setup does not permit further extrapolation to field situations, our findings, using three naturally co-occurring diatom species, call for a broader consideration of the role of benthic microbiota in shaping diatom community structure and function. Our findings fit well within the holobiont concept, which considers both hosts and their microbes as a single integrated unit (Zilber-Rosenberg and Rosenberg, 2008). Although it remains to be shown that co-evolution between diatoms and bacteria is commonplace, as is often stated in this theory (Skillings, 2016), the strong interactions between host-diatoms and their associated bacteria have important implications for the overall diatom fitness and thus their established niche (Vandenkoornhuyse et al., 2015; Kopac and Klassen, 2016) in species-rich natural communities.
Data Availability
The 16S rRNA raw sequence data were deposited in the NCBI Sequence Read Archive under the accession number PRJNA521472. The count data, growth rates and OTU table can be found on Github (https://github.com/willem-stock/Koedooder-et-al.-2019).
Author Contributions
CK, KS, and WS designed the study. CK and WS performed the experimental work and analyzed the data. CK and WS wrote the manuscript, which was revised by all authors. KS and WV provided funding for the research. All authors contributed to the discussion of the results.
Funding
The research leading to results presented in this publication was carried out with infrastructure funded by EMBRC Belgium – FWO project GOH3817N. This research was supported under the BIO-Tide project, funded through the 2015–2016 BiodivERsA COFUND call for research proposals, with the national funders BelSPO, FWO, ANR and SNSF. WS was funded by the Fund for Scientific Research – Flanders (FWO-Flanders, Belgium). The BCCM/DCG diatom culture collection was acknowledged for providing us with the appropriate diatom strains for this study. We are very grateful for help of Sofie D’Hondt with the molecular work and Tine Verstraete for her aid with the experimental work. Bjorn Tytgat provided the mock community.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01255/full#supplementary-material
FIGURE S1 | Example of picture-based diatom cell counting. The image shows a mixed diatom culture inoculated with bacteria at day 3. The culture contains S. robusta (pink, nr. 1), C. closterium (cyan, nr. 2), and N. phyllepta (yellow, nr 3). The image size is a fourth of the pictures taken for the cell counts.
FIGURE S2 | Rarefaction curve for every sample (before removal of rare OTUs) is shown. Colors indicate the different treatments as displayed.
FIGURE S3 | The presence of bacteria causes a shift in the relative proportions of diatoms. The relative proportion of C. closterium (red), S. robusta (yellow), and N. phyllepta (blue) grown in pairs or altogether expressed as the number of cells per species to the total number of cells in the well (Left = T1, Center = T4 in axenic conditions, and Right = T4 in non-axenic conditions). Briefly, in axenic conditions an increase in the proportion of C. closterium is observed when grown in the presence of S. robusta, while this observation was not observed in the presence of bacteria. Instead, S. robusta had an increased overall proportion in comparison to C. closterium.
FIGURE S4 | The change in the total biovolume (μmł) for each of the three diatom species [C. closterium (red), N. phyllepta (blue), and S. robusta (yellow)] 1 day after inoculation (green area) and 4 days after inoculation (purple area). The gray line indicates the biovolumes under non-axenic conditions while the black line specifies those under axenic conditions.
FIGURE S5 | PAM measurements depicting the total fluorescence as a proxy for the total maximum growth rates for different combinations of C. closterium ([C], red), N. phyllepta ([N], blue), and S. robusta ([S], yellow) in the presence (non-axenic) or absence of bacteria (axenic). Combinations of multiple diatoms species together have the appropriate colors combined. The white, gray and dark gray background, respectively show diatoms grown in monocultures, diatom pairs or the three diatom-species combined.
FIGURE S6 | PCA on the relative abundance of the bacterial community from the original inoculum (black circle) and when grown in the presence of different combinations of C. closterium (red), N. phyllepta (blue), and S. robusta (yellow). Samples with more diatoms present have circles with the appropriate colors combined. Gray arrows indicate the vectors of the bacterial OTUs with a cumulative loading larger than 30% on the first two axes. The bacterial genus or family is indicated depending on the closest hit. The proportion of the variance explained by each axis is indicated next to that axis.
FIGURE S7 | Diversity levels of bacteria in the presence of different diatom combinations of C. closterium (C, red), N. phyllepta (N, blue), and S. robusta (S, yellow). Bacterial community diversity was calculated as the average number of predicted OTUs after rarefication of the samples to 781 reads.
TABLE S1 | Information on the three different diatom strains used within our experimental setup. More information can be found at the Belgian Co-ordinated Collections of Micro-organisms website using the diatom accession number (http://bccm.belspo.be/catalogues).
TABLE S2 | Composition of the mock community. The raw fastq data of the mock is available at NCBI Sequence Read Archive under the accession number SUB5539877. Briefly, the DNA was extracted and amplified for each strain individually and equimolar pooled afterward. The collection number refers to the identifier of the strains in the BCCM/LMG collection (http://bccm.belspo.be/about-us/bccm-lmg; LMG) or the research collection of the Laboratory of Microbiology Ugent (http://lmg.ugent.be; R).
TABLE S3 | Statistical analysis of the total maximum growth rates derived from F0 PAM values of 24-plate wells harboring different combinations of S. robusta [S], C. closterium [C], and N. phyllepta [N]. A two-way ANOVA was conducted testing the differences between the total growth rate of the diatom combinations (single, paired, or all together) and the presence or absence of bacteria. A post hoc (Tukey Test) was further conducted on each of the different diatom combinations. An asterisk marks results of significant value where ∗<0.05, ∗∗<0.01, and ∗∗∗<0.001.
References
Admiraal, W., Peletier, H., and Brouwer, T. (1984). The seasonal succession patterns of diatom species on an intertidal mudflat: an experimental analysis. Oikos 42:30. doi: 10.2307/3544606
Agogue, H., Mallet, C., Orvain, F., De Crignis, M., and Mornet, F. (2014). Bacterial dynamics in a microphytobenthic biofilm: a tidal mesocosm approach. J. Sea Res. 92, 36–45. doi: 10.1016/J.SEARES.2014.03.003
Amin, S. A., Parker, M. S., and Armbrust, E. V. (2012). Interactions between diatoms and bacteria. Microbiol. Mol. Biol. Rev. 76, 667–684. doi: 10.1128/MMBR.00007-12
Azam, F. (1998). Microbial control of oceanic carbon flux: the plot thickens. Science (80-) 279, 483–484.
Baert, J. M., Janssen, C. R., Sabbe, K., and De Laender, F. (2016). Per capita interactions and stress tolerance drive stress-induced changes in biodiversity effects on ecosystem functions. Nat. Commun. 7:12486. doi: 10.1038/ncomms12486
Bagatini, I. L., Eiler, A., Bertilsson, S., Klaveness, D., Tessarolli, L. P., and Vieira, A. A. H. (2014). Host-specificity and dynamics in bacterial communities associated with bloom-forming freshwater phytoplankton. PLoS One 9:e85950. doi: 10.1371/journal.pone.0085950
Behringer, G., Ochsenkühn, M. A., Fei, C., Fanning, J., Koester, J. A., and Amin, S. A. (2018). Bacterial communities of diatoms display strong conservation across strains and time. Front. Microbiol. 9:659. doi: 10.3389/fmicb.2018.00659
Bellinger, B., Underwood, G., Ziegler, S., and Gretz, M. (2009). Significance of diatom-derived polymers in carbon flow dynamics within estuarine biofilms determined through isotopic enrichment. Aquat. Microb. Ecol. 55, 169–187. doi: 10.3354/ame01287
Bigalke, A., and Pohnert, G. (2019). Algicidal bacteria trigger contrasting responses in model diatom communities of different composition. Microbiologyopen [Epub ahead of print].
Blanco, S., Cejudo-Figueiras, C., Tudesque, L., Bécares, E., Hoffmann, L., and Ector, L. (2012). Are diatom diversity indices reliable monitoring metrics? Hydrobiologia 695, 199–206. doi: 10.1007/s10750-012-1113-1
Bohórquez, J., McGenity, T. J., Papaspyrou, S., García-Robledo, E., Corzo, A., and Underwood, G. J. C. (2017). Different types of diatom-derived extracellular polymeric substances drive changes in heterotrophic bacterial communities from intertidal sediments. Front. Microbiol. 8:245. doi: 10.3389/fmicb.2017.00245
Bolhuis, H., Fillinger, L., and Stal, L. J. (2013). Coastal microbial mat diversity along a natural salinity gradient. PLoS One 8:e63166. doi: 10.1371/journal.pone.0063166
Brinkhoff, T., Bach, G., Heidorn, T., Liang, L., Schlingloff, A., and Simon, M. (2004). 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. doi: 10.1128/aem.70.4.2560-2565.2003
Brotas, V., and Plante-Cuny, M.-R. (2003). The use of HPLC pigment analysis to study microphytobenthos communities. Acta Oecol. 24, S109–S115. doi: 10.1016/S1146-609X(03)00013-4
Bruhn, J. B., Gram, L., and Belas, R. (2007). Production of antibacterial compounds and biofilm formation by Roseobacter species are influenced by culture conditions. Appl. Environ. Microbiol. 73, 442–450. doi: 10.1128/AEM.02238-06
Bruno, J. F., Stachowicz, J. J., and Bertness, M. D. (2003). Inclusion of facilitation into ecological theory. Trends Ecol. Evol. 18, 119–125. doi: 10.1016/s0169-5347(02)00045-9
Callaway, R. M., Thelen, G. C., Rodriguez, A., and Holben, W. E. (2004). Soil biota and exotic plant invasion. Nature 427, 731–733. doi: 10.1038/nature02322
Cirri, E., Vyverman, W., and Pohnert, G. (2018). Biofilm interactions—bacteria modulate sexual reproduction success of the diatom Seminavis robusta. FEMS Microbiol. Ecol. 94:fiy161. doi: 10.1093/femsec/fiy161
Clarke, E. (2016). Levels of selection in biofilms: multispecies biofilms are not evolutionary individuals. Biol. Philos. 31, 191–212. doi: 10.1007/s10539-016-9517-3
Collado, L., and Figueras, M. J. (2011). Taxonomy, epidemiology, and clinical relevance of the genus Arcobacter. Clin. Microbiol. Rev. 24, 174–192. doi: 10.1128/CMR.00034-10
Cook, P., Van Oevelen, D., Soetaert, K., and Middelburg, J. (2009). Carbon and nitrogen cycling on intertidal mudflats of a temperate Australian estuary. IV. Inverse model analysis and synthesis. Mar. Ecol. Prog. Ser. 394, 35–48. doi: 10.3354/meps08270
Croft, M. T., Lawrence, A. D., Raux-Deery, E., Warren, M. J., and Smith, A. G. (2005). Algae acquire vitamin B12 through a symbiotic relationship with bacteria. Nature 438, 90–93. doi: 10.1038/nature04056
D’Costa, P. M., and Anil, A. C. (2011). The effect of bacteria on diatom community structure – the ‘antibiotics’ approach. Res. Microbiol. 162, 292–301. doi: 10.1016/j.resmic.2010.12.005
Decho, A. W. (2000). Microbial biofilms in intertidal systems: an overview. Cont. Shelf Res. 20, 1257–1273. doi: 10.1016/S0278-4343(00)00022-4
Decleyre, H., Heylen, K., Sabbe, K., Tytgat, B., Deforce, D., Van Nieuwerburgh, F., et al. (2015). A doubling of microphytobenthos biomass coincides with a tenfold increase in denitrifier and total bacterial abundances in intertidal sediments of a temperate estuary. PLoS One 10:e0126583. doi: 10.1371/journal.pone.0126583
Desianti, N., Potapova, M., Enache, M., Belton, T. J., Velinsky, D. J., Thomas, R., et al. (2017). Sediment diatoms as environmental indicators in new jersey coastal lagoons. J. Coast. Res. 78, 127–140. doi: 10.2112/SI78-011.1
D’Hondt, A.-S., Stock, W., Blommaert, L., Moens, T., and Sabbe, K. (2018). Nematodes stimulate biomass accumulation in a multispecies diatom biofilm. Mar. Environ. Res. 140, 78–89. doi: 10.1016/j.marenvres.2018.06.005
Dittami, S. M., Duboscq-Bidot, L., Perennou, M., Gobet, A., Corre, E., Boyen, C., et al. (2016). Host-microbe interactions as a driver of acclimation to salinity gradients in brown algal cultures. ISME J. 10, 51–63. doi: 10.1038/ismej.2015.104
Doiron, K., Linossier, I., Fay, F., Yong, J., Abd Wahid, E., Hadjiev, D., et al. (2012). Dynamic approaches of mixed species biofilm formation using modern technologies. Mar. Environ. Res. 78, 40–47. doi: 10.1016/j.marenvres.2012.04.001
Durham, B. P., Sharma, S., Luo, H., Smith, C. B., Amin, S. A., Bender, S. J., et al. (2015). Cryptic carbon and sulfur cycling between surface ocean plankton. Proc. Natl. Acad. Sci. U.S.A. 112, 453–457. doi: 10.1073/pnas.1413137112
Eigemann, F., Hilt, S., Salka, I., and Grossart, H.-P. (2013). Bacterial community composition associated with freshwater algae: species specificity vs. dependency on environmental conditions and source community. FEMS Microbiol. Ecol. 83, 650–663. doi: 10.1111/1574-6941.12022
Ferreira, S., Queiroz, J. A., Oleastro, M., and Domingues, F. C. (2015). Insights in the pathogenesis and resistance of Arcobacter: a review. Crit. Rev. Microbiol. 42, 1–20. doi: 10.3109/1040841X.2014.954523
Forster, R., Créach, V., Sabbe, K., Vyverman, W., and Stal, L. (2006). Biodiversity-ecosystem function relationship in microphytobenthic diatoms of the Westerschelde estuary. Mar. Ecol. Prog. Ser. 311, 191–201. doi: 10.3354/meps311191
Fridley, J. D. (2001). The influence of species diversity on ecosystem productivity: how, where, and why? Oikos 93, 514–526. doi: 10.1034/j.1600-0706.2001.930318.x
Furusawa, G., Yoshikawa, T., Yasuda, A., and Sakata, T. (2003). Algicidal activity and gliding motility of Saprospira sp. SS98-5. Can. J. Microbiol. 49, 92–100. doi: 10.1139/w03-017
Gessner, M. O., Inchausti, P., Persson, L., Raffaelli, D. G., Giller Gessner, P. S., Persson, P., et al. (2004). Biodiversity effects on ecosystem functioning: insights from aquatic systems. Oikos 104, 419–422. doi: 10.1111/j.0030-1299.1999.13252.x
Giller, P. S., Hillebrand, H., Berninger, U.-G., Gessner, M. O., Hawkins, S., Inchausti, P., et al. (2004). Biodiversity effects on ecosystem functioning: emerging issues and their experimental test in aquatic environments. Oikos 104, 423–436. doi: 10.1111/j.0030-1299.2004.13253.x
Gribben, P. E., Nielsen, S., Seymour, J. R., Bradley, D. J., West, M. N., and Thomas, T. (2017). Microbial communities in marine sediments modify success of an invasive macrophyte. Sci. Rep. 7:9845. doi: 10.1038/s41598-017-10231-2
Grime, J. P., Mackey, J. M. L., Hillier, S. H., and Read, D. J. (1987). Floristic diversity in a model system using experimental microcosms. Nature 328, 420–422. doi: 10.1038/328420a0
Grossart, H.-P. (1999). Interactions between marine bacteria and axenic diatoms (Cylindrotheca fusiformis, Nitzschia laevis, and Thalassiosira weissflogii) incubated under various conditions in the lab. Aquat. Microb. Ecol. 19, 1–11. doi: 10.3354/ame019001
Grossart, H.-P., Levold, F., Allgaier, M., Simon, M., and Brinkhoff, T. (2005). Marine diatom species harbour distinct bacterial communities. Environ. Microbiol. 7, 860–873. doi: 10.1111/j.1462-2920.2005.00759.x
Haines, K. C., and Guillard, R. R. L. (1974). Growth of vitamin B12-requiring marine diatoms in mixed laboratory cultures with vitamin B12- producing marine bacteria. J. Phycol. 10, 245–252. doi: 10.1111/j.1529-8817.1974.tb02709.x
Havskum, H., Thingstad, T. F., Peters, F., Berdalet, E., and Sala, M. M. (2003). Silicate and labile DOC interfere in structuring the microbial food web via algal-bacteria competition for mineral nutrients: results of a mesocosm experiment. Limnol. Ocean. 48, 129–140. doi: 10.4319/lo.2003.48.1.0129
Haynes, K., Hofmann, T. A., Smith, C. J., Ball, A. S., Underwood, G. J. C., and Osborn, A. M. (2007). Diatom-derived carbohydrates as factors affectin‘g bacterial community composition in estuarine sediments. Appl. Environ. Microbiol. 73, 6112–6124. doi: 10.1128/AEM.00551-07
Hillebrand, H., Dürselen, C.-D., Kirschtel, D., Pollingher, U., and Zohary, T. (1999). Biovolume calculations for pelagic and benthic microalgae. J. Phycol. 35, 403–424. doi: 10.1046/j.1529-8817.1999.3520403.x
Hochard, S., Pinazo, C., Grenz, C., Evans, J. L. B., and Pringault, O. (2010). Impact of microphytobenthos on the sediment biogeochemical cycles: a modeling approach. Ecol. Modell. 221, 1687–1701. doi: 10.1016/J.ECOLMODEL.2010.04.002
Horner-Devine, C. M., Leibold, M. A., Smith, V. H., and Bohannan, B. J. M. (2003). Bacterial diversity patterns along a gradient of primary productivity. Ecol. Lett. 6, 613–622. doi: 10.1046/j.1461-0248.2003.00472.x
Hortal, S., Lozano, Y. M., Bastida, F., Armas, C., Moreno, J. L., Garcia, C., et al. (2017). Plant-plant competition outcomes are modulated by plant effects on the soil bacterial community. Sci. Rep. 7:17756. doi: 10.1038/s41598-017-18103-5
Hubbard, J. A. M., Lewandowska, K. B., Hughes, M. N., and Poole, R. K. (1986). Effects of iron-limitation of Escherichia coli on growth, the respiratory chains and gallium uptake. Arch. Microbiol. 146, 80–86. doi: 10.1007/BF00690163
Jasti, S., Sieracki, M. E., Poulton, N. J., Giewat, M. W., and Rooney-Varga, J. N. (2005). Phylogenetic diversity and specificity of bacteria closely associated with Alexandrium spp. and other phytoplankton. Appl. Environ. Microbiol. 71, 3483–3494. doi: 10.1128/AEM.71.7.3483-3494.2005
Jung, S. W., Kim, B.-H., Katano, T., Kong, D.-S., and Han, M.-S. (2008). Pseudomonas fluorescens HYK0210-SK09 offers species-specific biological control of winter algal blooms caused by freshwater diatom Stephanodiscus hantzschii. J. Appl. Microbiol. 105, 186–195. doi: 10.1111/j.1365-2672.2008.03733.x
Kardol, P., Cornips, N. J., van Kempen, M. M. L., Bakx-Schotman, J. M. T., and van der Putten, W. H. (2007). MIcrobe-mediated plant-soil feedback causes historical contingency effects in plant community assembly. Ecol. Monogr. 77, 147–162. doi: 10.1890/06-0502
Kassen, R., Buckling, A., Bell, G., and Rainey, P. B. (2000). Diversity peaks at intermediate productivity in a laboratory microcosm. Nature 406, 508–512. doi: 10.1038/35020060
Kopac, S. M., and Klassen, J. L. (2016). Can they make it on their own? Hosts, microbes, and the holobiont niche. Front. Microbiol. 7:1647. doi: 10.3389/fmicb.2016.01647
Laforest-Lapointe, I., Paquette, A., Messier, C., and Kembel, S. W. (2017). Leaf bacterial diversity mediates plant diversity and ecosystem function relationships. Nature 546, 145–147. doi: 10.1038/nature22399
Lavergne, C., Agogue, H., Leynaert, A., Raimonet, M., De Wit, R., Pineau, P., et al. (2017). Factors influencing prokaryotes in an intertidal mudflat and the resulting depth gradients. Estuar. Coast. Shelf Sci. 189, 74–83. doi: 10.1016/j.ecss.2017.03.008
Lekberg, Y., Bever, J. D., Bunn, R. A., Callaway, R. M., Hart, M. M., Kivlin, S. N., et al. (2018). Relative importance of competition and plant-soil feedback, their synergy, context dependency and implications for coexistence. Ecol. Lett. 21, 1268–1281. doi: 10.1111/ele.13093
Loreau, M., and Hector, A. (2001). Partitioning selection and complementarity in biodiversity experiments. Nature 412, 72–76. doi: 10.1038/35083573
MacIntyre, H. L., Geider, R. J., and Miller, D. C. (1996). Microphytobenthos: the ecological role of the “secret garden” of unvegetated, shallow-water marine habitats. I. Distribution, abundance and primary production. Estuaries 19:186. doi: 10.2307/1352224
Marzinelli, E. M., Qiu, Z., Dafforn, K. A., Johnston, E. L., Steinberg, P. D., and Mayer-Pinto, M. (2018). Coastal urbanisation affects microbial communities on a dominant marine holobiont. Biofilms Microbiomes 4:1. doi: 10.1038/s41522-017-0044-z
Mayali, X., Franks, P. J. S., and Azam, F. (2008). Cultivation and ecosystem role of a marine roseobacter clade-affiliated cluster bacterium. Appl. Environ. Microbiol. 74, 2595–2603. doi: 10.1128/AEM.02191-07
Mensens, C., De Laender, F., Janssen, C. R., Sabbe, K., and De Troch, M. (2015). Stressor-induced biodiversity gradients: revisiting biodiversity-ecosystem functioning relationships. Oikos 124, 677–684. doi: 10.1111/oik.01904
Middelburg, J. J., Barranguet, C., Boschker, H. T. S., Herman, P. M. J., Moens, T., and Heip, C. H. R. (2000). The fate of intertidal microphytobenthos carbon: an in situ 13 C-labeling study. Limnol. Oceanogr. 45, 1224–1234. doi: 10.4319/lo.2000.45.6.1224
Moerdijk-Poortvliet, T., Beauchard, O., Stal, L., and Boschker, H. (2018). Production and consumption of extracellular polymeric substances in an intertidal diatom mat. Mar. Ecol. Prog. Ser. 592, 77–95. doi: 10.3354/meps12481
Mühlenbruch, M., Grossart, H.-P., Eigemann, F., and Voss, M. (2018). Mini-review: phytoplankton-derived polysaccharides in the marine environment and their interactions with heterotrophic bacteria. Environ. Microbiol. 20, 2671–2685. doi: 10.1111/1462-2920.14302
Muyzer, G., de Waal, E. C., and Uitterlinden, A. G. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 59, 695–700.
Paul, C., Barofsky, A., Vidoudez, C., and Pohnert, G. (2009). Diatom exudates influence metabolism and cell growth of co-cultured diatom species. Mar. Ecol. Prog. Ser. 389, 61–70. doi: 10.3354/meps08162
Paul, C., and Pohnert, G. (2011). Interactions of the algicidal bacterium kordia algicida with diatoms: regulated protease excretion for specific algal lysis. PLoS One 6:e21032. doi: 10.1371/journal.pone.0021032
Pugnaire, F. I., Armas, C., and Valladares, F. (2004). Soil as a mediator in plant-plant interactions in a semi-arid community. J. Veg. Sci. 15, 85–92. doi: 10.1111/j.1654-1103.2004.tb02240.x
Raes, E. J., Bodrossy, L., van de Kamp, J., Bissett, A., and Waite, A. M. (2018). Marine bacterial richness increases towards higher latitudes in the eastern Indian Ocean. Limnol. Oceanogr. Lett. 3, 10–19. doi: 10.1002/lol2.10058
Ribeiro, L., Brotas, V., Rincé, Y., and Jesus, B. (2013). Structure and diversity of intertidal benthic diatom assemblages in contrasting shores: a case study from the Tagus estuary 1. J. Phycol. 49, 258–270. doi: 10.1111/jpy.12031
Ruiz-Ponte, C., Samain, J. F., Sánchez, J. L., and Nicolas, J. L. (1999). The benefit of a roseobacter species on the survival of scallop larvae. Mar. Biotechnol. 1, 52–59. doi: 10.1007/PL00011751
Sawai, Y., Horton, B. P., Kemp, A. C., Hawkes, A. D., Nagumo, T., and Nelson, A. R. (2016). Relationships between diatoms and tidal environments in Oregon and Washington, USA. Diatom Res. 31, 17–38. doi: 10.1080/0269249X.2015.1126359
Schafer, H., Abbas, B., Witte, H., and Muyzer, G. (2002). Genetic diversity of satellite bacteria present in cultures of marine diatoms. FEMS Microbiol. Ecol. 42, 25–35. doi: 10.1111/j.1574-6941.2002.tb00992.x
Seyedsayamdost, M. R., Case, R. J., Kolter, R., and Clardy, J. (2011). The jekyll-and-hyde chemistry of Phaeobacter gallaeciensis. Nat. Chem. 3, 331–335. doi: 10.1038/nchem.1002
Sharifah, E. N., and Eguchi, M. (2011). The phytoplankton nannochloropsis oculata enhances the ability of roseobacter clade bacteria to inhibit the growth of fish pathogen Vibrio anguillarum. PLoS One 6:e26756. doi: 10.1371/journal.pone.0026756
Shishlyannikov, S., Zakharova, Y., Volokitina, N., Mikhailov, I., Petrova, D., and Likhoshway, Y. (2011). A procedure for establishing an axenic culture of the diatom Synedra acus subsp. radians (Kütz.) Skabibitsch. from Lake Baikal. Limnol. Ocean. Methods 9, 478–484. doi: 10.4319/lom.2011.9.478
Sison-Mangus, M. P., Jiang, S., Tran, K. N., and Kudela, R. M. (2014). Host-specific adaptation governs the interaction of the marine diatom, Pseudo-nitzschia and their microbiota. ISME J. 8, 63–76. doi: 10.1038/ismej.2013.138
Skillings, D. (2016). Holobionts and the ecology of organisms: multi-species communities or integrated individuals? Biol. Philos. 31, 875–892. doi: 10.1007/s10539-016-9544-0
Slightom, R. N., and Buchan, A. (2009). Surface colonization by marine roseobacters: integrating genotype and phenotype. Appl. Environ. Microbiol. 75, 6027–6037. doi: 10.1128/AEM.01508-09
Soininen, J. (2009). Is diatom diversity driven by productivity in boreal streams? Diatom Res. 24, 197–207. doi: 10.1080/0269249X.2009.9705791
Stein, A., Gerstner, K., and Kreft, H. (2014). Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales. Ecol. Lett. 17, 866–880. doi: 10.1111/ele.12277
Stock, W., Blommaert, L., Daveloose, I., Vyverman, W., and Sabbe, K. (2019). Assessing the suitability of Imaging-PAM fluorometry for monitoring growth of benthic diatoms. J. Exp. Mar. Bio. Ecol. 513, 35–41. doi: 10.1016/j.jembe.2019.02.003
Taylor, J. D., McKew, B. A., Kuhl, A., McGenity, T. J., and Underwood, G. J. C. (2013). Microphytobenthic extracellular polymeric substances (EPS) in intertidal sediments fuel both generalist and specialist EPS-degrading bacteria. Limnol. Oceanogr. 58, 1463–1480. doi: 10.4319/lo.2013.58.4.1463
Thingstad, T. F., Skjoldal, E. F., and Bohne, R. A. (1993). Phosphorus cycling and algal-bacterial competition in Sandsfjord, western Norway. Mar. Ecol. Progr. Ser. 99, 239–259. doi: 10.3354/meps099239
Thornton, D., Dong, L., Underwood, G., and Nedwell, D. (2002). Factors affecting microphytobenthic biomass, species composition and production in the Colne Estuary (UK). Aquat. Microb. Ecol. 27, 285–300. doi: 10.3354/ame027285
Tuchman, N. C., Schollett, M. A., Rier, S. T., and Geddes, P. (2006). Differential heterotrophic utilization of organic compounds by diatoms and bacteria under light and dark conditions. Hydrobiologia 561, 167–177. doi: 10.1007/s10750-005-1612-1614
Tytgat, B., Verleyen, E., Sweetlove, M., D’hondt, S., Clercx, P., Van Ranst, E., et al. (2016). Bacterial community composition in relation to bedrock type and macrobiota in soils from the Sør Rondane Mountains, East Antarctica. FEMS Microbiol. Ecol. 92:fiw126. doi: 10.1093/femsec/fiw126
Underwood, G. J. C., and Kromkamp, J. (1999). Primary production by phytoplankton and microphytobenthos in estuaries. Adv. Ecol. Res. 29, 93–153. doi: 10.1016/S0065-2504(08)60192-0
Underwood, G. J. C., Perkins, R. G., Consalvey, M. C., Hanlon, A. R. M., Oxborough, K., Baker, N. R., et al. (2005). Patterns in microphytobenthic primary productivity: species-specific variation in migratory rhythms and photosynthetic efficiency in mixed-species biofilms. Limnol. Oceanogr. 50, 755–767. doi: 10.4319/lo.2005.50.3.0755
Vandenkoornhuyse, P., Quaiser, A., Duhamel, M., Le Van, A., and Dufresne, A. (2015). The importance of the microbiome of the plant holobiont. New Phytol. 206, 1196–1206. doi: 10.1111/nph.13312
Vanelslander, B., De Wever, A., Van Oostende, N., Kaewnuratchadasorn, P., Vanormelingen, P., Hendrickx, F., et al. (2009). Complementarity effects drive positive diversity effects on biomass production in experimental benthic diatom biofilms. J. Ecol. 97, 1075–1082. doi: 10.1111/j.1365-2745.2009.01535.x
Wagner-Dobler, I., Rheims, H., Felske, A., El-Ghezal, A., Flade-Schröder, D., Laatsch, H., et al. (2004). Oceanibulbus indolifex gen. nov., sp. nov., a North Sea alphaproteobacterium that produces bioactive metabolites. Int. J. Syst. Evol. Microbiol. 54, 1177–1184. doi: 10.1099/ijs.0.02850-0
Wear, E. K., Carlson, C. A., James, A. K., Brzezinski, M. A., Windecker, L. A., and Nelson, C. E. (2015). Synchronous shifts in dissolved organic carbon bioavailability and bacterial community responses over the course of an upwelling-driven phytoplankton bloom. Limnol. Oceanogr. 60, 657–677. doi: 10.1002/lno.10042
Keywords: cross-kingdom interactions, diversity–productivity relationship, microphytobenthos, algae-bacteria relationship, biofilm interactions
Citation: Koedooder C, Stock W, Willems A, Mangelinckx S, De Troch M, Vyverman W and Sabbe K (2019) Diatom-Bacteria Interactions Modulate the Composition and Productivity of Benthic Diatom Biofilms. Front. Microbiol. 10:1255. doi: 10.3389/fmicb.2019.01255
Received: 18 February 2019; Accepted: 20 May 2019;
Published: 05 June 2019.
Edited by:
Ezequiel Miguel Marzinelli, The University of Sydney, AustraliaReviewed by:
Simon M. Dittami, UMR8227 Laboratoire de Biologie Intégrative des Modèles Marins, FranceAlexandra Helene Campbell, University of New South Wales, Australia
Copyright © 2019 Koedooder, Stock, Willems, Mangelinckx, De Troch, Vyverman and Sabbe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Willem Stock, Willem.Stock@ugent.be
†Present address: Coco Koedooder, CNRS, UPMC Univ Paris 06, Laboratoire d’Océanographie Microbienne, Observatoire Océanologique, Sorbonne Université, Banyuls-sur-Mer, France
‡These authors have contributed equally to this work