Abstract
Host-associated microbiota are critical for eukaryotic host functioning, to the extent that hosts and their associated microbial communities are often considered “holobionts”. Most studies of holobionts have focused on descriptive approaches or have used model systems, usually in the laboratory, to understand host-microbiome interactions. To advance our understanding of host-microbiota interactions and their wider ecological impacts, we need experimental frameworks that can explore causation in non-model hosts, which often have highly diverse microbiota, and in their natural ecological setting (i.e. in the field). We used a dominant habitat-forming seaweed, Hormosira banksii, to explore these issues and to experimentally test host-microbiota interactions in a non-model holobiont. The experimental protocols were aimed at trying to disentangle microbially mediated effects on hosts from direct effects on hosts associated with the methods employed to manipulate host-microbiota. This was done by disrupting the microbiome, either through removal/disruption using a combination of antimicrobial treatments, or additions of specific taxa via inoculations, or a combination of thew two. The experiments were done in mesocosms and in the field. Three different antibiotic treatments were used to disrupt seaweed-associated microbiota to test whether disturbances of microbiota, particularly bacteria, would negatively affect host performance. Responses of bacteria to these disturbances were complex and differed substantially among treatments, with some antibacterial treatments having little discernible effect. However, the temporal sequence of responses antibiotic treatments, changes in bacterial diversity and subsequent decreases in host performance, strongly suggested an effect of the microbiota on host performance in some treatments, as opposed to direct effects of the antibiotics. To further test these effects, we used 16S-rRNA-gene sequencing to identify bacterial taxa that were either correlated, or uncorrelated, with poor host performance following antibiotic treatment. These were then isolated and used in inoculation experiments, independently or in combination with the previously used antibiotic treatments. Negative effects on host performance were strongest where specific microbial antimicrobials treatments were combined with inoculations of strains that were correlated with poor host performance. For these treatments, negative host effects persisted the entire experimental period (12 days), even though treatments were only applied at the beginning of the experiment. Host performance recovered in all other treatments. These experiments provide a framework for exploring causation and disentangling microbially mediated vs. direct effects on hosts for ecologically important, non-model holobionts in the field. This should allow for better predictions of how these systems will respond to, and potentially mitigate, environmental disturbances in their natural context.
Similar content being viewed by others
Introduction
Evidence from systems as diverse as coral reefs and the human gut indicates that host-associated microbial communities are fundamentally important for the normal development and functioning of eukaryotic hosts1,2,3,4. This is leading to a paradigm shift in biology, where eukaryotes are no longer viewed as single entities, but rather as part of a coherent biological association comprised of the host and their microbiome, i.e. as “holobionts”5,6. Most studies of holobiont ecology use observational approaches, that is, descriptions of host or microbial characteristics of the holobiont, often measured in association with variation in environmental parameters or host performance (e.g refs. 7,8). Such approaches typically involve sequencing and ‘omics’ techniques that have allowed detailed characterisations of microbial diversity and, in some cases, predicted function of holobionts9,10,11. Where experiments are used (see refs. 12,13,14), they are often on model organisms (although see refs. 13,14,15,16) in the laboratory17,18,19.
The use of model systems has enabled us to determine direct, microbially mediated effects, and their underlying mechanisms in a number of model eukaryotic hosts, particularly those associated with one or a few microbial symbionts, such as for the bobtail squid Euphrymna scolopes and its bacterial symbiont Aliivibrio fisheri20. Within macroalgal systems, the clearest demonstrations of the role of complex microbial communities on host function (reviewed in ref. 21) are from the model algae Ulva spp22,23,24 and Ectocarpus spp12,25. Beyond model systems, there has been work exploring the role of microbiota in the function of hosts such as mosquitos15,26,27,28 and daphnia16,29. However, in many of these examples of non-model hosts, the experimental manipulations occurred in the laboratory, removed from their natural environment, limiting the extent of our inferences and our capacity to predict responses of holobionts in nature.
Descriptive studies are important in providing knowledge of the diversity (who is there?) and the functional potential (what can/do they do?) of host-associated microorganisms and how microbial diversity and function vary with host performance. These studies, however, do not untangle cause-effect relationships, where the effects of complex microbial communities on hosts (and vice versa) can be distinguished and tested30,31. While description of patterns is a fundamental first step in the scientific framework32, to understand complex systems we need to move beyond descriptions to experimental manipulations that determine causation and from which potential mechanisms can subsequently be explored33,34,35. Causation indicates that one event—change in host performance - is the result of the occurrence of the other event—change in host-associated microbiota, i.e. there is a causal relationship between the two events. Importantly, not understanding the mechanisms underlying microbial effects on hosts does not impede establishing causation36. Relying on inferred causation only from observations (e.g. associations between variables) hinders, for example, our capacity to predict the responses of holobionts to disturbances or changing environmental conditions31,37. The most notable, if rare, example of an experimental system for the manipulation of complex microbiota is the use of “germ-free” mice in biomedical systems, which require intensive, whole-of-life-cycle facilities to raise mice without their normal microbiome17,18. While laboratory systems such as these provide important opportunities to explore direct versus indirect effects of microbial manipulation in controlled settings, their findings often do not represent the responses of natural populations38,39,40. Thus, to understand the natural ecology of, for example, key habitat forming holobionts, we need to do experiments in the field. This would align the holobiont perspective with that of macrobial ecology, when half a century ago a largely descriptive field was transformed into a strong field based experimental discipline (e.g refs. 41,42,43). Unlike in macrobial ecology, the question of how to experimentally approach holobiont ecology remains a bottleneck and has become a major focus of modern microbial ecology31. This is partly due to the complexity of most holobionts which harbour thousands of taxa44,45,46, giving rise to the potential for millions of interactions among these taxa and between the microbiota and the host, but is also due to the challenges of developing appropriate experimental protocols.
While significant advances on these issues have been made in model and/or laboratory systems the process of unravelling the complex associations and interactions between hosts and their microbiota is at a very early stage for most non-model holobionts living in open, “real-world” environmental systems where large, structurally complex, habitat-forming holobionts such as trees, kelps or corals create the biogenic structure of ecosystems. Understanding of host-microbiome interactions for such foundational holobionts is particularly significant because the impacts on the interaction between these hosts and their microbiota can cascade throughout entire ecosystems (e.g. ref. 47).
Experimental frameworks that allow establishing causation of host-microbiota interactions in holobionts are beginning to be developed30,48. In general, to determine microbially mediated effects of microbiota on host functioning (as opposed to direct effects on the host per se), we need experimental approaches that manipulate the microbiota in a controlled manner (where there are no effects of the manipulations other than on the microbiota) and, if possible, in a specific fashion using targeted manipulations of specific microorganisms correlated with host performance. Moreover, ideally for the environmental systems of concern here, manipulations would be conducted in the field (e.g. refs. 49,50,51) to avoid unwanted effects that are linked to the laboratory environment, as has been the case for many other biotic interactions in ecology32. This would have the additional advantage of integrating experiments on holobionts into the established existing framework of field experiments for macrobial ecology.
One particular challenge for such manipulations is that many of the experimental protocols do not disentangle direct effects of treatments on the hosts from effects of manipulations of the microbiota on hosts. Experimental manipulations on holobionts can involve the selective removal/reduction of microorganisms of interest (e.g., those found to correlate either positively/negatively with host performance), the addition of such microorganisms or a combination of both. The removal or reduction/disruption of host-associated microorganisms using antibiotics is a common approach (perhaps the most common approach in the field; reviewed in17) that has been used to study, for example, effects on host immunity, metabolism and disease52,53,54. However, the key question which arises from their use is: how do we decouple the potential direct physical/chemical effects of antimicrobials on host performance (e.g. refs. 55,56) from indirect, microbially-mediated effects? This is a complex and ongoing issue as different antimicrobials cause vastly different effects, both within the host and on its microbiota57,58.
We used a dominant, habitat-forming marine holobiont to develop and test a general experimental framework to determine direct versus indirect effects of the microbiota, i.e. cause-effect relationships between microbiota and host performance, in non-model holobionts in the laboratory and in the field. We combined microbial reduction/disruption, using antibiotics, with inoculations of specific microbial taxa associated with a widespread marine habitat-forming fucoid alga, Hormosira banksii (hereafter Hormosira). Hormosira forms extensive intertidal forests along Australian and New Zealand rocky shores but has declined in many urbanised areas due to multiple anthropogenic stressors (e.g. refs. 59,60). As with other similar seaweeds, such declines may be mediated by disruptions of their associated microbiome8.
We first used different antibiotic treatments that reduced/disrupted different components of Hormosira’s surface-associated microbiota to test the hypothesis that such microbial disruption would negatively affect host performance. To distinguish direct impacts of antibiotics from indirect, microbially-mediated effects, we isolated, and subsequently inoculated, bacterial taxa whose abundance were either correlated (and thus predicted to affect the host when added) or uncorrelated (no effect predicted) with effects on the host. By inoculating such taxa, independently or in combination with microbiota reductions over an ecologically relevant temporal scale, we were able to establish that microbial changes could affect host performance. The benefits and challenges of these approaches, individually and in concert, are explored.
Results
Mesocosm experiment
Antibiotic treatments significantly changed the bacterial community structure on Hormosira in the mesocosm experiment, although, there were significant differences among different antibiotics and over time (PERMANOVA pseudo-F9, 32 = 1.29, p < 0.001; Supplementary Table 7; Fig. 1A). The effects of the antibiotics were not immediate, with no effect of treatments at time 1. As of day 2 Hormosira treated with antibiotic mix 2 had a significantly different community structure than all other treatments and that difference remained until the end of the experiment at day 12 (F9, 32 = 2.71, p < 0.001). Antibiotic mix 3’s microbial community differed significantly from controls at day 2 which lasted till day 12 (F9, 32 = 1.94, p < 0.001; Fig. 1A) and was significantly different from controls and all other treatments at day 12 (F 9, 32 = 3.63, p < 0.001; Fig. 1A).
Community richness was significantly lower in the treatment antibiotic mix 2 than the control (ANOVA, F3, 80 = 6.05, p < 0.001; Supplementary Table 8), Antibiotic mix 3 had significantly lower richness than controls (ANOVA, F3, 80 = 4.84, p < 0.016). Simpson’s diversity index was significantly lower in the antibiotic mix 3 than all other treatments at the end of the experiment (day 12, ANOVA, F3, 80 = 5.49, p < 0.001; Supplementary Table 8).
These changes in the microbiome correlated with changes in host photosynthesis, which occurred ~3 days after treatment with AB2. Hormosira’s photosynthetic yield under antibiotic mix 2 experienced a continuous and significant decrease after 5 days (ANOVA, F9,32 = 8.56, p < 0.001). Algae treated with antibiotic mix 3 also had decreased photosynthetic yield compared to controls, but only towards the end of the experiment, with the most significant decrease at day 12 (ANOVA, F9,32 = 7.16, p < 0.001; Fig. 1B). Algae treated with the antibiotic mix 1 or in the control did not show any change in photosynthetic yield over time (Fig. 1A).
Field experiment
In the field, bacterial relative abundance and community structure differed significantly across time and among treatments (PERMANOVA, pseudo-F20,98 = 1.6505, p < 0.001, Fig. 1C Supplementary Table 9). Antibiotics and iodine both caused significant changes in community structure from day 7 which lasted for the 46 days of the experiment (PERMANOVA, pseudo-F20,98 = 3.35, p < 0.001). Multiple applications of iodine caused greater changes in community structure than a single application alone (PERMANOVA, F5,98 = 3.506, p < 0.001). Antibiotics, when administered once, had a similar effect to multiple applications of iodine (PERMANOVA, t = 1.706, p > 0.05; Supplementary Table 9). The number of bacterial ASVs significantly decreased in algae treated with antibiotic and iodine relative to controls (ANOVA, F20,98 = 4.92, p < 0.001).
Photosynthetic yield of the host differed significantly among treatments, with individuals in the antibiotic treatment AB2 showing a significant decline in photosynthetic efficiency over time (ANOVA, F4,162 = 6.145, p < 0.001; Supplementary Table 6), starting at day 9 with the greatest effect occurring at day 21. This effect remained significant till day 46, when photosynthetic efficiency did not differ from controls (Fig. 1D; Supplementary Table 6). Algae treated multiple times with iodine (IM) also showed lower photosynthetic yield than controls and single application of iodine but did not recover by the end of the experiment. Again, this was after a lag period of 21 days (ANOVA, F20,162 = 3.201, p < 0.001; Fig. 1C; Supplementary Table 6). No other treatments differed from the controls which was coupled with no significant differences in the community composition within these treatments.
Bacterial Isolation
To identify strains correlated with host function we isolated bacteria, leading to over 140 isolates, of which 22 strains were negatively correlated with poor host performance within both the mesocosm and field experiments. The abundance/presence of Vibrio genomosp. F10 was most strongly correlated with poor host performance (r = −0.68, p < 0.0001; Supplementary Fig. 2) with its abundance significantly increasing in all disruption treatments except AB1 (ANOVA, F9,32 = 4.39, p < 0.001). There were also many bacterial taxa whose abundances were not correlated with host performance, e.g., Vibrio chagasii (r = 0.14, p < 0.12; Supplementary Fig. 2).
Inoculation experiment
GLM analyses identified 172 ASVs (~1.8% of a total of 9138 ASVs) whose abundances differed significantly among treatments over the duration of the inoculation experiment. Of those, we found a strong main effect of treatment on 10 ASVs in the classes Gammaproteobacteria, Alphaproteobacteria and Cyanophyceae, with abundances of the latter two being lower in all treatments with microbiome disruption. There was a significant increase in Vibrio species alongside a significant increase of the inoculants V. chagasii and V. genomosp. F10 whose abundances were significantly higher within inoculation treatments (ANOVA, F5,24 = 1.77, p < 0.001; Fig. 2A).
Inoculation with V. genomosp. F10, when combined with a disruption of the microbiome, had the strongest effect on microbial community structure, which differed from all other treatments at day 2 and lasted until the end of the experiment (ANOVA, F5,24 = 4.28, p < 0.001; Fig. 3). The microbiome structure of algae disrupted with antibiotics and inoculated with V. chagasii did not differ from disruption with antibiotics but was different from controls (ANOVA, F5,24 = 1.31, p < 0.01; Fig. 3C) indicating that whilst the strain was present it did not influence community structure. No significant differences were found throughout the experiment among the three control treatments (Control, Procedural control, and inoculation control, Supplementary Table 12). Furthermore, all treatments with antibiotics had significantly lower richness than controls until day 8 where there was no significant difference among treatments (ANOVA, F5,24 = 0.969, p = 0.21).
All antibiotic treatments, including both treatments where the algae were further inoculated with both Vibrio spp., had significantly lower photosynthetic yield than those where microbiomes were left undisturbed and this effect was observed ~2 days after observed microbial changes (ANOVA, F5,24 = 4.25, p < 0.001, Fig. 3A). Inoculation with V. genomosp. F10 had a similar effect on photosynthesis as treatments with antibiotics but without inoculants or inoculated with V. chagasii (ANOVA, F5,24 = 1.58, p < 0.03, Fig. 3A). All treatments, except microbial disruption followed by inoculation with V. genomosp. F10, recovered to the same level of photosynthetic yield that was present at the start of the experiment by day 12.
Discussion
The structure and function of host-associated microbiomes are strongly related with host performance, but as is highlighted in the recent literature61, we still have very limited understanding of cause-effect in these interactions in complex, non-model holobionts within their natural ecological setting. Here, we identified the effect of the host-associated bacterial community, or subsets thereof, on the performance (as measured by photosynthetic efficiency) of a dominant, marine habitat-forming holobiont via a combination of removals and additions of microbial taxa associated with the host’s performance in mesocosms and field experiments. We disrupted the host’s microbiota using a series of different antibiotic mixes with several methods of action, but which are not thought to have direct toxic effects on multicellular algae at the concentrations used (see refs. 62,63). Some of these treatments had no effect on the host even though they affected the hosts microbiota providing further evidence that the observed decreases in host performance were indeed microbially mediated. We showed that the reduction/disruption of different components of the microbiota affected host performance in different ways, which changed in a non-linear fashion over time. Temporal sampling and the sequence of temporal sampling added further evidence for microbially-mediated effect being the primary consequence of these treatments. This is because any direct chemical effect of the antimicrobials on the host would likely be immediate63,64 (although see refs. 65,66,67) whereas microbially-mediated effects would take time following the removal/decline of taxa68. This was the case for our experiments in both the lab and the field. Furthermore, using inoculants which were either correlated or uncorrelated with host performance further confirmed that it was a microbially mediated effect.
In the mesocosm experiment, when exposed to antibiotics, all bacterial communities were significantly altered. These changes were followed within some treatments by decreases in host performance, i.e., there was a lag between the bacterial changes and the effect on host performance. In addition, different antibiotic mixes affected the bacterial community differently despite having similar microbial targets (Supplementary Table 2). For instance, AB1 influenced the bacterial community but not the host, while AB2 and AB3 affected the bacterial community and subsequently the host, but host responses had different lag times (3 vs 5 days after microbial changes, respectively).
These different combinations of treatments were all done to try and disentangle direct vs. Indirect (microbially mediated) effects of the treatments on the host. Because (i) bacterial changes occurred before host changes, (ii) the antibiotics had no discernible immediate effect on the host which would be expected from previous work and (iii) the three antibiotic mixes had similar MOAs but varying targets, all this suggests no (undesired) direct effects of the antibiotics on the host and changes in host performance were therefore likely caused by changes in the microbiome. These experiments and our results also speak to the challenges associated with controlling for undesired, direct effects of antimicrobial treatments on hosts. From our results, we argue that the combination of the use of several antimicrobials with similar targets but varying MOAs and chemical structures, and temporal sampling of the microbiota with host performance in parallel provides a robust approach to understanding microbially-mediated effects on host performance for non-model holobionts. The antibiotics used in our work were selected because previous studies suggested that they do not have chemical effects on the host25.
Other approaches for disentangling the chemical effect of antibiotics include reintroducing the microbiota in its entirety following its removal, i.e. conventionalisation26,29. Whilst this technique of reintroducing the microbiota has been highly successful in a number of hosts ranging from daphnia16,29 and mosquitos15,26, to zebra fish69 and mice17,54, applying this technique in the field can be highly challenging. Whilst there has been some success in hosts such as corals51,53,70 through the homogenisation of the hosts tissue, its applicability to other highly complex microbiota for hosts in the field is yet to be tested. Future studies could attempt this technique in the field, but should also corroborate that re-inoculating entire, complex microbial assemblages re-instates the community at similar levels of relative abundances, spatial distribution, ecological interactions and functioning of all microbial taxa (typically hundreds to thousands of ASVs) on the host as those found on unmanipulated hosts in nature.
The antibiotic mix which had the most significant and lasting effect was antibiotic mix 2 (AB2), whose difference from the other antibiotic mixes was the inclusion of the microbial RNA polymerase inhibiting Rifampicin. Antibiotic mix 3 (AB3) also contained a similar antibiotic target (Chloramphenicol) and was close to AB2’s efficacy in significantly decreasing relevant bacterial taxa from the microbiome. Indeed, the main group negatively affected by these two treatments were nitrogen-fixing cyanobacteria which we were unable to get into a pure culture. The most abundant cyanobacteria ranged from 1–3% (Pleurocapsa sp) of the total abundance. N-fixing bacteria such as these have been implicated in numerous symbioses with terrestrial, and more recently, marine plants71,72,73. N is a limiting nutrient within intertidal environments with many algal species often relying on extracellular N fixation to meet their physiological needs74,75. This may explain the observed decline in host performance.
Microbiota manipulations in the field experiment supported the results from our mesocosm experiment as microbial disruption was followed by negative effects on host performance. As in the lab, the microbial disruption caused by AB2 had a significant negative effect on host function, which eventually recovered over time (45 Days). Antiseptics (Povidone iodine 10%) had a similar, negative effect on host function when applied once (pulse) but caused lasting, negative effects on the host when applied multiple times (press). This approach mimics the natural environment closer where Hormosira is exposed to both press (e.g., increasing ocean temperatures76,77, micrograzers78,79,80) and pulse (e.g., desiccation at low tide81 or extreme rainfall) stressors over varying timescales.
Interestingly, cyanobacteria were the main bacteria which remained at low abundances during the press treatments both in the laboratory and in the field. Whilst our disruptions within both the lab and field show consistent patterns, many others have shown difficulty in relating findings in the lab to real-world settings (see ref. 82), highlighting the need to run these experiments in the field wherever possible. This is particularly important in the context of environmental change because stressors can interact in complex ways, leading to different responses of organisms depending on the environmental context83,84.
Disruption of the microbiota allowed for our second experimental approach—inoculations of relevant taxa – to be done in an ecologically realistic way. From the lab and field experiments we were able to isolate 22 unique strains associated with poor host performance. These isolates were dominated by readily culturable taxa from the orders Altermonadales and Vibrionales, with several ASVs being assigned to taxa that are known pathogens. Although the isolated taxa made up ~1% of total taxonomic diversity within the community, their abundance contributed to ~5% of the total. We used some of these taxa in inoculation experiments combined with microbial disruptions to determine their independent or combined effect on host performance. We were unable to culture bacterial taxa positively correlated with host function in these experiments, though we did culture taxa that were uncorrelated to performance. The inoculation of beneficial bacteria would be a logical next step for the experimental protocol developed here and could be a focus of future work on the positive effects of microbiota on host performance. Alongside this, due to the increasing recognition of the prevalence of marine fungi associated with hosts (see refs. 85,86), future work should also include attempts to isolate and reintroduce members of the community.
Microbial disruption, when combined with inoculations of bacterial taxa associated with poor host performance (Vibrio genomosp. F10), led to a much stronger negative effect on host performance than disruption or inoculation alone, and host performance did not recover during our experiment. This effect may be due to competitive release (see refs. 87,88,89), where inoculants may have been able to establish and affect the host because of the disruption/reduction of other taxa, e.g. competitors. This treatment also led to a shift in the bacterial community structure that did not recover throughout the experiment, suggesting that a threshold may have been reached where the resilience of the microbial community was overcome. Environmental stressors can lead to dysbiosis or an imbalance in the microbiota of holobionts as diverse as corals, kelps and humans, which, in turn, can make them more susceptible to being affected by pathogens77,90. Thus, the combination of microbial disturbance and inoculations with taxa related to poor host performance can help understand and better predict responses of holobionts to environmental stressors and the potential mechanisms underlying such responses.
Bacteria belonging to the genera Vibrio, Pseudoalteromonas and Phaeobacter have been described as antagonistic towards their hosts91, although this is species-specific as certain Phaeobacter species have been shown to provide antifouling and antilarval abilities92. The specificity of the effects of different strains is further demonstrated as inoculation with our control inoculant Vibrio chagasii spp, despite its close taxonomic relation to Vibrio genomosp. F10, did not induce a negative effect on host performance – this was expected as this taxon was chosen as an inoculation control because its abundance did not correlate with host performance in our first set of experiments. Sequencing analyses after inoculation showed a significant increase in the abundance of both Vibrio chagsii spp and Vibrio genomosp. F10 indicating that our inoculation was successful in inducing temporary establishment within the host microbiota. However, this establishment did not last throughout the experiment; abundances of both inoculants returned to control levels by the end of the experimental period. Both strains have been reported as opportunists within other host species and are not found in high abundances within the undisturbed microbiota of Hormosira, which may be the result of them being outcompeted by other taxa in the microbial assemblage. This initial increase in abundance followed by rapid decrease is supported by recent work showing that commensal and mutualistic bacterial taxa switch to antagonistic pathogens when disrupted93.
To enable us to understand how holobionts will respond to environmental change, it is important to study them within their natural ecological settings. This is inherently difficult in natural, open systems, where thousands of taxa may contribute to millions of potential interactions. Here we have developed an experimental protocol that combines microbiota disruption using antibiotics, and the identification and subsequent inoculation of presumed harmful bacteria. This was done in both the laboratory and the field. We believe this combined approach, applied (wherever possible) in real-world settings, provides a powerful and necessary approach towards understanding of these interactions.
We showed that disruption of the host-associated microbiota using different antibiotics and subsequent temporal sampling of microbial changes and host performance provided critical information on the microorganisms potentially associated with host function, a pattern that was consistent in mesocosms and in the field. We then targeted, isolated and cultured relevant microbial taxa associated with changes in host performance (or lack thereof, as controls) and used them in controlled experiments to formally test their impact on the host. While not definitive, this approach increased the strength of our inferences as we can start to disentangle the effects of disturbance on the broader microbial community (e.g., dysbiosis caused by the antimicrobials) from the effect of individual strains of relevant microbial taxa, or both combined.
There are critical aspects that still need to be considered when interpreting outcomes of these experimental approaches, particularly the complex responses of microorganisms to different antimicrobials, the potential direct (not microbially mediated) effects on the host of the methodology employed to manipulate microbiota, and the temporal sequence of events and the environmental context (i.e. laboratory vs field). Future work should also incorporate, where possible, techniques used in laboratory experiments, such as conventionalisation, which provides a further control for manipulative experiments involving the removal/disruption of the a microbiota17,54,69. Another logical progression of this work is to incorporate “natural” disruptions, as those influenced by environmental changes, allowing better predictions of how these ecosystems will respond to, and potentially mitigate, future environmental disturbances.
Methods
Relationships between microbiome and host performance
Mesocosm experiment
Hormosira individuals of similar frond size (mean diameter: 7.8 + /− S.E 2.2 cm, length: 11.2 + /− S.E 0.81 cm; N = 120) were haphazardly collected by carefully detaching their holdfast from intertidal rocks during low tide at Cronulla beach, Sydney, Australia (34°03’22.8” S 151°09’19.7” E) on 4 December 2018. Hormosira were transported in seawater within 1 h to the flow-through aquarium facility at the Sydney Institute of Marine Science (SIMS) where they were thoroughly rinsed in autoclaved filtered (0.2 µm) seawater (AFSW) to remove fouling organisms and loosely associated microorganisms. Hormosira individuals were then placed into 12, 2 litre tanks (10 individuals 2.5 cm apart per tank) by cable-tying their holdfasts to weighted mesh (e.g. ref. 94) and left overnight to acclimatize. Tanks were supplied by filtered (5 µm), UV-treated, running seawater. The following day, algae were rinsed with AFSW and randomly assigned to one of three antibiotic combinations (AB1, AB2 or AB3; see details in Supplementary Table 2) or an AFSW control. Algae were left to soak in the different antibiotic treatments or in AFSW (control) for 24 h, then rinsed with AFSW and placed back into the 12 independent tanks (n = 3 tanks per treatment; 10 individuals per tank) with UV treated, 5 µm filtered seawater.
One Hormosira individual was randomly selected and removed from each tank prior to the addition of any antibiotics (day 0), just after their treatment with antibiotics (day 1) and every second day for 1 week and on day 12 to destructively sample individuals to characterise surface-associated microbiome and host performance (see details in Supplementary Information). Temporal sampling allows for the detection of unwanted direct effects as it is predicted that if antibiotics had any direct (e.g., chemical) effects on the host, such effects would be immediate (see ref. 62). On the other hand, if effects on host are microbially-mediated, it is predicted that there would be a lag between the application of the treatment and the response of the host (i.e. time for the microbiota to respond to the treatment and consequentially elicit a response on the host). Whilst temporal sampling provides evidence of direct effects if present, it is fundamentally important to also use suitable controls, i.e., in this case multiple antibiotics with similar methods of action (MOAs) but varying bacterial targets.
To characterise the associated microbiome, a consistent area of the algal surface at mid-thallus (13 + /− SE 2.7 cm2; older, non-meristematic tissue adjacent to, but independent from the area where photosynthesis was measured, see below) was swabbed for 30 s with sterile cotton swabs, which were immediately placed in sterile cryogenic tubes in liquid nitrogen and then stored at −80 °C until DNA extraction. After bacterial sampling, Host performance was assessed by quantifying the maximum photosynthetic quantum yield of each individual (n = 3, after individuals were dark-adapted for 15 min) using a Pulse Amplitude Modulated (PAM) fluorometer (WALZ, Germany), a metric of host performance widely used for photosynthetic organisms (e.g. plants and macroalgae) particularly under different environmental conditions and in relation with changes to their surface-associated microbiota8,95,96,97,98,99.
Field experiment
Forty-two Hormosira individuals of similar size (~6.8 + /− 2.3 cm in diameter and 8.5 + /− SE 1.3 cm in length) ~1 m apart were tagged at Cape Banks, Sydney, Australia (33°59’55.3” S 151°14’53.6” E) during low tide on 13 May 2018. To apply treatment to the individual alga, a sponge halo was placed over the individual alga before they were rinsed with AFSW for 1 minute and then subjected to one of 6 different treatments: (1) the antibiotic combination AB2 (chosen due to its effect on the microbiome and host performance in the mesocosm experiment above) applied once for 2 h during exposure at low tide, (2) Povidone Iodine 10% w/v (‘iodine’ hereafter) applied once for 2 h (IO), (3) iodine applied every ~2 days for 15 minutes (IM), (4) a procedural control for the treatment application once (PC; applying AFSW once for 2 h), (5) a procedural control for the continuous treatment application (IMC; applying AFSW for 15 mins every ~2 days), (6) control (C, undisturbed individual; see details in Supplementary Table 1).
Half of the surface area of each individual (~12 + /− 2.7 cm2) was swabbed for 30 seconds to quantify the surface-associated microbiome as described for the mesocosm experiment. Swab controls (i.e., used similarly to those for swabbing alga except no algae was swabbed) were collected. The other half was swabbed for bacterial culturing (details below); these swabs were stored in 1.5 ml Eppendorf filled with 1 ml of AFSW and placed in ice until arrival to the laboratory. Additional swab controls were collected to test for potential contamination in cultures due to the placement in AFSW. Samples were collected every 2nd day for a total of 46 days (see Supplementary Table 2 for specific dates).
Bacterial isolation and culture
For both the laboratory and field experiments, swabs in 1.5 ml Eppendorf tubes filled with 1 ml of AFSW were immediately taken to the laboratory for spreading on plates. Briefly, the suspended swab was vortexed and 10 ul of the solution was then plated on standard half strength marine broth agar (MP Biomedicals LLC). Individual morphologically distinct colonies were picked after 24 h. These colonies were suspended in liquid marine broth and left to grow for 48 h at room temperature on a shaker plate. 10 ul of the liquid culture was then plated as above on a new sterile agar plate. Colonies were replated until only 1 morphological form remained which was then extracted and sequenced (details below). Isolates were stored in glycerol at −80 °C until resuspension in marine broth for inoculations (see Supplementary Information for growth conditions).
Isolates were sequenced (details below) and the identity of the taxa were matched with those from the mesocosm and field experiment. Microbial amplicon sequence variant (hereafter, “ASV”) abundance was correlated with host photosynthetic yield to identify isolates for the inoculation experiments. Two strains were used for inoculations: Vibrio genomosp (F10 str. 9ZC157) and Vibrio chagasii (str ECSMB 14107). V. genomosp F10 was chosen as its increase in abundance post-antibiotic treatment was directly correlated with poor host performance (PAM) and had the greatest change in total abundance over the experiment (Supplementary Fig. 1). V. chagasii was chosen as it was phylogenetically similar to V. genomosp. F10 but remained at a constant abundance throughout the experiment and was thus uncorrelated to variation in host performance. Furthermore, the use of a control inoculation (Vibrio chagasii) was important to control for the potential effect of adding in a high density bacterial inoculant. For all of these reasons, Vibrio chagasii was therefore used as a control inoculant (Supplementary Fig. 1). In order to use both strains as inoculants they were grown to an OD of 600 before passing the cultures through a 0.2 µm filter which was being rinsed off with AFSW to the same OD.
Testing for direct effects of components of the microbiome
Inoculation experiment
Hormosira individuals of similar length (7.5 cm +/− SE 1.2; N = 110) were haphazardly collected as described above from the rocky shore at Cape Banks on 25 January 2021 and transported to the SIMS aquarium within an hour, where they were rinsed with AFSW. Hormosira individuals (N = 110) were then placed into 15 ml falcon tubes which had the bottoms removed (1 individual per tube) and attached via their holdfast to a stainless-steel rod at the bottom of the tube (Supplementary Fig. 3). Each tube was individually fed filtered (50 µm 1 L/min) seawater through aquarium airlines from four replicate line splitters. Once all samples were attached in the aquarium, they were left overnight to acclimatise. On the next day, all tubes were rinsed with AFSW and the algae inside was randomly assigned to one of seven treatments (Supplementary Table 1):
1. Microbiome disruption with antibiotic combination AB2 applied once for 25 h (as in the mesocosm experiment; Supplementary Table 1);
2. Microbiome disruption as in (1) followed by inoculation with bacteria whose relative abundances correlated with poor host performance in the experiments above (Vibrio genomosp. F10, see Supplementary Fig. 1; cell density 5 × 108 CFU);
3. Microbiome disruption as in (1) followed by inoculation with bacteria unrelated with host performance in the experiments above (negative control: Vibrio chagsii sp.; see Supplementary Fig. 1; cell density 5 × 108 CFU);
4. Microbiome left undisturbed (AFSW) for 24 h followed by inoculation with V. genomosp. F10 as in (2);
5. Microbiome left undisturbed (AFSW) for 24 h followed by inoculation with the negative control strain as in (3);
6. Procedural control, AFSW for 25 h;
7. Control (undisturbed algae).
Inoculations were done by immersing Hormosira for 1 h within the high-density cell culture which allowed sufficient time for the inoculated bacteria to attach to algal surfaces (See Supplementary Fig. 1). This timing was confirmed by leaving individuals within high cell density washed cultures for varying times before being swabbed as above. DNA was extracted and custom qPCR primers targeting the V3-V4 regions of the bacterial and archaeal 16 S rRNA gene for each of the inoculants were used to determine their abundance over time (see Supplementary Information for further details on primer design)
Samples were collected following the same timepoints as in the previous experiment: 5 individuals were collected at day 0 before treatments were applied and 3 independent individuals were randomly selected from each treatment on days 1, 2, 5, 8 and 12, at the same time of day on each sampling occasion. These time-points and the overall duration of this experiment were based on results from the previous experiments above which showed effects on microbiome and hosts within ~3–5 days (see Results) and our capacity to maintain independent replicate algae for each treatment x time combination. Surface-associated bacteria and algal photosynthetic efficiency was then characterised for each sample as described above (See mesocosm experiment above).
DNA extraction and sequencing
For characterisation of microbial communities in all experiments, microbial DNA was extracted from each swab sample in a randomised order using a PowerSoil DNA Isolation kit (Qiagen) following the manufacturers protocol. DNA extracts were quantified using spectrophotometry (NanoDrop 1000) and stored at -20 °C until sequencing.
The extracted DNA samples were amplified with Polymerase Chain Reaction (PCR) using the 16 S rRNA gene primers 341 (F) (5’- CCTACGGGNGGCWGCAG-3’) and 805(R) – (5’-GACTACHVGGGTATCTAATCC-‘3), covering the V3-V4 regions of the bacterial and archaeal 16 S rRNA gene100. The PCR conditions involved a pre-heating step to 95 °C for 3 min followed by 35 cycles of 95 °C for 15 s, 55 °C for 1 min and 73 °C for 30 s. Both positive (with known DNA sequence) and negative controls (nuclease-free water, control swabs) were used. The negative controls did not amplify DNA, suggesting no contamination on swabs/materials or during extraction and amplification. Agarose gel electrophoresis and Nanodrop 1000 were used to ensure the quantity and quality of the amplicons before they were sent for sequencing via the Illumina MiSeq 2000 platform at the Ramaciotti Centre for Genomics (UNSW, Sydney).
Bioinformatics
Raw sequencing data was quality filtered using Trimmomatic101 with a sliding window trim of 4:15 base pairs (bp) and removal of sequences with <36 bp. Paired-end reads were merged with a minimum length of 400 bp and maximum of 500 bp using USEARCH102. UNOISE was then used to remove chimeras and produce amplicon sequence variants (ASVs), i.e., operational taxonomic units at a unique sequence level (0% distance)102. USEARCH was used to map the original reads back to ASVs, generating a table of 5812, 7637 and 9138 ASVs for the mesocosm, field and inoculation experiments, respectively. ASV sequences were searched with BlastN against the SILVA SSU Ref NR99 database for taxonomic classification to classify and remove chloroplasts, the Genome Taxonomy Database (GTDB) was then used for taxonomic assignment. Singletons and low abundance taxa (<0.01% of reads) were removed from the dataset for statistical analyses, resulting in 4112, 6488 and 7028 ASVs (mesocosm, field and inoculation experiments respectively).
Estimation of absolute bacterial abundance
Total abundance of the 16 S rRNA gene was quantified for each sample by qPCR using the primers 341 F/805 R103. Gene amplification and analysis were performed using the QuantStudio 3 thermocycler (Thermo Fisher with PrimeTime® Gene Expression Master Mix, Integrated DNA Technologies) and associated software. The reaction conditions for amplification of DNA were 56 °C for 2 min, 95 °C for 10 min and 40 cycles of 95 °C for 15 s and 60 °C for 1 min. The final gene copy number per sample was corrected for the total extraction volume, the surface area and the dilution factor and DNA yield per sample (see ref. 104 for further details on the normalisation) and were used to estimate absolute abundances of ASVs and inoculants.
Statistical analyses
To account for uneven sequencing depth among samples, data were normalised using the counts per million reads method in the R package DESeq2105. This resulted in a total of 4,423,807, 10,443,709 and 9,921,043 sequences for the mesocosm, field experiment and inoculation datasets, respectively.
To test for effects of microbiome manipulations on host photosynthetic efficiency, we used a linear model with the factors treatment (fixed) and time (fixed, crossed). For the first mesocosm experiment where multiple Hormosira individuals were sampled from each tank, ‘tank’ was fitted as a random factor nested within treatment using the lme4 R package (v4.0.3). To meet the model’s assumptions, data were square root transformed. p-values were obtained using the anova function within the car package in R with significance being tested using F-tests, or likelihood ratio tests for the mixed model that included tank as a random effect.
Alpha diversity measures of richness (i.e., number of unique sequences) and Simpson’s diversity index were calculated using the ‘vegan’ R package106 and differences between treatment (fixed) and time (fixed, crossed) were examined using a linear model in the R GAD package. Where appropriate, tank was included as a random factor.
To determine differences in the structure of the associated bacterial assemblages, the normalised ASV data were analysed using permutational multivariate analysis of variance (PERMANOVA)107 in the R vegan package106, with the factors treatment, time and tank (mesocosm experiment only) as above. Similarity matrices were calculated using Bray-Curtis measure on square-root transformed data and visualised through non-metric multi-dimensional scaling (nMDS) ordinations. We also calculated and plotted mean Bray-Curtis similarities to the control (undisturbed algae) treatment. To determine which bacterial taxa’s abundance differed among treatments and times, we used multivariate generalised linear models (GLMs) using the R package ‘mvabund’108, assuming a negative-binomial distribution to account for over-dispersion. To determine how the total abundance of the two inoculants changed over the course of the experiment, linear models were fitted with the factors treatment (fixed) and time (fixed, crossed) in R (v4.0.3).
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
All data, sampling locations, amplicon sequence data, host response and timepoints are available on zenodo https://doi.org/10.5281/zenodo.8031031.
Change history
03 May 2024
A Correction to this paper has been published: https://doi.org/10.1038/s41522-024-00515-7
References
Rosenberg, E., Koren, O., Reshef, L., Efrony, R. & Zilber-Rosenberg, I. The role of microorganisms in coral health, disease and evolution. Nat. Rev. Microbiol. 5, 355–362 (2007).
Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K. & Knight, R. Diversity, stability and resilience of the human gut microbiota. Nature 489, 220–230 (2012).
McFall-Ngai, M. et al. Animals in a bacterial world, a new imperative for the life sciences. Proc. Natl Acad. Sci. USA 110, 3229–3236 (2013).
Wilkins, L. G. E., Ettinger, C. L., Jospin, G. & Eisen, J. A. Metagenome-assembled genomes provide new insight into the microbial diversity of two thermal pools in Kamchatka, Russia. Sci. Rep. 9, 3059 (2019).
Rosenberg, E. & Zilber-Rosenberg, I. The hologenome concept of evolution after 10 years. Microbiome 6, 78 (2018).
Dittami, S. M. et al. A community perspective on the concept of marine holobionts: current status, challenges, and future directions. Peer J. 9, e10911 (2021).
Nielsen, S. et al. Disrupted progression of the intestinal microbiota with age in children with cystic fibrosis. Sci. Rep. 6, 24857 (2016).
Qiu, Z. et al. Future climate change is predicted to affect the microbiome and condition of habitat-forming kelp. Proc. R. Soc. B Biol. Sci. 286, 20181887 (2019).
Gilbert, J. A. & Dupont, C. L. Microbial metagenomics: beyond the genome. Annu. Rev. Mar. Sci. 3, 347–371 (2011).
Hug, L. A. et al. A new view of the tree of life. Nat. Microbiol. 1, 16048 (2016).
Lloyd-Price, J. et al. Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases. Nature 569, 655–662 (2019).
KleinJan, H. et al. Changing microbial activities during low salinity acclimation in the brown alga ectocarpus subulatus. bioRxiv https://doi.org/10.1101/2021.04.13.439635 (2021).
Houwenhuyse, S., Stoks, R., Mukherjee, S. & Decaestecker, E. Locally adapted gut microbiomes mediate host stress tolerance. ISME J. 15, 2401–2414 (2021).
Coon, K. L., Vogel, K. J., Brown, M. R. & Strand, M. R. Mosquitoes rely on their gut microbiota for development. Mol. Ecol. 23, 2727–2739 (2014).
Coon, K. L., Brown, M. R. & Strand, M. R. Mosquitoes host communities of bacteria that are essential for development but vary greatly between local habitats. Mol. Ecol. 25, 5806–5826 (2016).
Decaestecker, E., Declerck, S., De Meester, L. & Ebert, D. Ecological implications of parasites in natural daphnia populations. Oecologia 144, 382–390 (2005).
Kennedy, E. A., King, K. Y. & Baldridge, M. T. Mouse microbiota models: comparing germ-free mice and antibiotics treatment as tools for modifying gut bacteria. Front. Physiol. 9, 1534 (2018).
Fontaine, C. A. et al. How free of germs is germ-free? detection of bacterial contamination in a germ free mouse unit. Gut Microbes 6, 225–233 (2015).
Douglas, A. E. Simple animal models for microbiome research. Nat. Rev. Microbiol. 17, 764–775 (2019).
Koch, E. J. & McFall-Ngai, M. Model systems for the study of how symbiotic associations between animals and extracellular bacterial partners are established and maintained. Drug Discov. Today Dis. Models 28, 3–12 (2018).
Egan, S. et al. The seaweed holobiont: understanding seaweed-bacteria interactions. FEMS Microbiol. Rev. 37, 462–476 (2013).
Provasoli, L. & Pintner, I. J. Bacteria Iinduced polymorphism in an axenic laboratory stain of ulva lactuca (Chlorophyceae)1. J. Phycol. 16, 196–201 (1980).
Nakanishi, K., Nishijima, M., Nishimura, M., Kuwano, K. & Saga, N. Bacteria that induce morphogenesis in Ulva pertusa (Chlorophyta) grown under axenic conditions. J. Phycol. 32, 479–482 (1996).
Ghaderiardakani, F., Coates, J. C. & Wichard, T. Bacteria-induced morphogenesis of Ulva intestinalis and Ulva mutabilis (Chlorophyta): a contribution to the lottery theory. FEMS Microbiol. Ecol. 93, fix094 (2017).
Dittami, S. M. et al. Host–microbe interactions as a driver of acclimation to salinity gradients in brown algal cultures. ISME J. 10, 51–63 (2016).
Cansado-Utrilla, C., Zhao, S. Y., McCall, P. J., Coon, K. L. & Hughes, G. L. The microbiome and mosquito vectorial capacity: rich potential for discovery and translation. Microbiome 9, 111 (2021).
Dada, N. et al. Considerations for mosquito microbiome research from the mosquito microbiome consortium. Microbiome 9, 36 (2021).
Rodríguez-Ruano, S. M., Juhaňáková, E., Vávra, J. & Nováková, E. Methodological insight Into mosquito microbiome studies. Front. Cell. Infect. Microbiol. 10, 86 (2020).
Massol, F., Macke, E., Callens, M. & Decaestecker, E. A methodological framework to analyse determinants of host–microbiota networks, with an application to the relationships between Daphnia magna’ s gut microbiota and bacterioplankton. J. Anim. Ecol. 90, 102–119 (2021).
Dittami, S. M., Eveillard, D. & Tonon, T. A metabolic approach to study algal-bacterial interactions in changing environments. Mol. Ecol. 23, 1656–1660 (2014).
Prosser, J. I. & Martiny, J. B. H. Conceptual challenges in microbial community ecology. Philos. Trans. R. Soc. B Biol. Sci. 375, 20190241 (2020).
Underwood, A. J., Chapman, M. G. & Connell, S. D. Observations in ecology: you can’t make progress on processes without understanding the patterns. J. Exp. Mar. Biol. Ecol. 250, 97–115 (2000).
Underwood, A. J., Underwood, A. L. & Wnderwood, A. Experiments In Ecology: Their Logical Design And Interpretation Using Analysis Of Variance, 524 (Cambridge university press, 1997).
Jansson, J. K. & Prosser, J. I. The life beneath our feet. Nature 494, 40–41 (2013).
Prosser, J. I. How and why in microbial ecology: an appeal for scientific aims, questions, hypotheses and theories. Environ. Microbiol. 24, 4973–4980 (2022).
Ross, L. N. Causal concepts in biology: how pathways differ from mechanisms and why it matters. Br. J. Philos. Sci. 72, 131–158 (2021).
Trevathan-Tackett, S. M. et al. A horizon scan of priorities for coastal marine microbiome research. Nat. Ecol. Evol. 3, 1509–1520 (2019).
Morrow, J. L., Frommer, M., Shearman, D. C. A. & Riegler, M. The microbiome of field-caught and laboratory-adapted Australian tephritid fruit fly species with different host plant use and specialisation. Microb. Ecol. 70, 498–508 (2015).
Hegde, S. et al. Microbiome interaction networks and community structure from laboratory-reared and field-collected Aedes aegypti, Aedes albopictus, and Culex quinquefasciatus mosquito vectors. Front. Microbiol. 9, 2160 (2018).
Calisi, R. M. & Bentley, G. E. Lab and field experiments: are they the same animal? Horm. Behav. 56, 1–10 (2009).
Connell, J. H. The influence of interspecific competition and other factors on the distribution of the barnacle chthamalus stellatus. Ecology 42, 710–723 (1961).
Paine, R. T. Food web complexity and species diversity. Am. Nat. 100, 65–75 (1966).
Simberloff, D. S. & Wilson, E. O. Experimental zoogeography of islands: the colonization of empty islands. Ecology 50, 278–296 (1969).
Donaldson, G. P., Lee, S. M. & Mazmanian, S. K. Gut biogeography of the bacterial microbiota. Nat. Rev. Microbiol. 14, 20–32 (2016).
Marzinelli, E. M., Leong, M. R., Campbell, A. H., Steinberg, P. D. & Vergés, A. Does restoration of a habitat-forming seaweed restore associated faunal diversity? Restor. Ecol. 24, 81–90 (2016).
Bourne, D. G., Morrow, K. M. & Webster, N. S. Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annu. Rev. Microbiol. 70, 317–340 (2016).
Laforest-Lapointe, I., Paquette, A., Messier, C. & Kembel, S. W. Leaf bacterial diversity mediates plant diversity and ecosystem function relationships. Nature 546, 145–147 (2017).
Carrier, T. J. & Reitzel, A. M. The hologenome across environments and the implications of a host-associated microbial repertoire. Front. Microbiol. 8, 802 (2017).
Campbell, A. H., Harder, T., Nielsen, S., Kjelleberg, S. & Steinberg, P. D. Climate change and disease: bleaching of a chemically defended seaweed: climate change and seaweed bleaching. Glob. Change Biol. 17, 2958–2970 (2011).
Campbell, A. H., Marzinelli, E. M., Vergés, A., Coleman, M. A. & Steinberg, P. D. Towards restoration of missing underwater forests. PLoS ONE 9, e84106 (2014).
Peixoto, R. S. et al. Coral probiotics: premise, promise, prospects. Annu. Rev. Anim. Biosci. 9, 265–288 (2021).
Glasl, B., Herndl, G. J. & Frade, P. R. The microbiome of coral surface mucus has a key role in mediating holobiont health and survival upon disturbance. ISME J. 10, 2280–2292 (2016).
Peixoto, R. S., et al Beneficial microorganisms for corals (BMC): proposed mechanisms for coral health and resilience. Front. Microbiol. https://doi.org/10.3389/fmicb.2017.00341 (2017).
Blacher, E. et al. Potential roles of gut microbiome and metabolites in modulating ALS in mice. Nature 572, 474–480 (2019).
Rubin, B. K. & Tamaoki, J. Antibiotics As Anti-Inflammatory And Immunomodulatory Agents 5th edn, 276 (Birkhauser Verlag AG, 2005).
Kohanski, M. A., Dwyer, D. J. & Collins, J. J. How antibiotics kill bacteria: from targets to networks. Nat. Rev. Microbiol. 8, 423–435 (2010).
Luo, Y. & Song, Y. Mechanism of antimicrobial peptides: antimicrobial, anti-inflammatory and antibiofilm activities. Int. J. Mol. Sci. 22, 11401 (2021).
Looft, T. et al. In-feed antibiotic effects on the swine intestinal microbiome. Proc. Natl Acad. Sci. USA 109, 1691–1696 (2012).
Bellgrove, A., McKenzie, P. F., Cameron, H. & Pocklington, J. B. Restoring rocky intertidal communities: lessons from a benthic macroalgal ecosystem engineer. Mar. Pollut. Bull. 117, 17–27 (2017).
Doblin, M. A. & Clayton, M. N. Effects of secondarily-treatedsewage effluent on the early life-history stages of two species of brown macroalgae: Hormosirabanks#and durvillaea potatorum. Mar. Biol. 122, 10 (1995).
Marchesi, J. R. et al. The gut microbiota and host health: a new clinical frontier. Gut 65, 330–339 (2016).
Li, S., Show, P. L., Ngo, H. H. & Ho, S.-H. Algae-mediated antibiotic wastewater treatment: a critical review. Environ. Sci. Ecotechnol. 9, 100145 (2022).
Fu, L. et al. Toxicity of 13 different antibiotics towards freshwater green algae pseudokirchneriella subcapitata and their modes of action. Chemosphere 168, 217–222 (2017).
Ebert, I. et al. Toxicity of the fluoroquinolone antibiotics enrofloxacin and ciprofloxacin to photoautotrophic aquatic organisms. Environ. Toxicol. Chem. 30, 2786–2792 (2011).
Leclercq, S. et al. Low-dose penicillin in early life induces long-term changes in murine gut microbiota, brain cytokines and behavior. Nat. Commun. 8, 15062 (2017).
Schokker, D. et al. Long-lasting effects of early-life antibiotic treatment and routine animal handling on gut microbiota composition and immune system in pigs. PLOS ONE 10, e0116523 (2015).
Cheng, G. et al. Antibiotic alternatives: the substitution of antibiotics in animal husbandry? Front. Microbiol. https://doi.org/10.3389/fmicb.2014.00217 (2014).
Schmidt, R. & Saha, M. Infochemicals in terrestrial plants and seaweed holobionts: current and future trends. N. Phytol. 229, 1852–1860 (2021).
Levraud, J.-P., Rawls, J. F. & Clatworthy, A. E. Using zebrafish to understand reciprocal interactions between the nervous and immune systems and the microbial world. J. Neuroinflammation 19, 170 (2022).
Peixoto, R. S. et al. Harnessing the microbiome to prevent global biodiversity loss. Nat. Microbiol. 7, 1726–1735 (2022).
Mohr, W. et al. Terrestrial-type nitrogen-fixing symbiosis between seagrass and a marine bacterium. Nature 600, 105–109 (2021).
Fiore, C. L., Jarett, J. K., Olson, N. D. & Lesser, M. P. Nitrogen fixation and nitrogen transformations in marine symbioses. Trends Microbiol. 18, 455–463 (2010).
Nutman, P. S. Genetics of symbiosis and nitrogen fixation in legumes. Proc. R. Soc. Lond. B Biol. Sci. 172, 417–437 (1969).
Bracken, M. E. S. & Nielsen, K. J. Diversity of intertidal macroalgae increases with nitrogen loading by invertebrates. Ecology 85, 2828–2836 (2004).
Pfister, C. A. Intertidal invertebrates locally enhance primary production. Ecology 88, 1647–1653 (2007).
Smale, D. A. Impacts of ocean warming on kelp forest ecosystems. N. Phytol. 225, 1447–1454 (2020).
Greenspan, S. E. et al. Warming drives ecological community changes linked to host-associated microbiome dysbiosis. Nat. Clim. Change 10, 1057–1061 (2020).
Schiel, D. R. & Foster, M. S. The population biology of large brown seaweeds: ecological consequences of multiphase life histories in dynamic coastal environments. Annu. Rev. Ecol. Evol. Syst. 37, 343–372 (2006).
Underwood, A. J. Grazing and disturbance: an experimental analysis of patchiness in recovery from a severe storm by the intertidal alga Hormosira banksii on rocky shores in New South Wales. J. Exp. Mar. Biol. Ecol. 231, 291–306 (1998).
Wootton, H. & Keough, M. Disturbance type and intensity combine to affect resilience of an intertidal community. Mar. Ecol. Prog. Ser. 560, 121–133 (2016).
Macinnis-Ng, C. M. O., Morrison, D. A. & Ralph, P. J. Temporal and spatial variation in the morphology of the brown macroalga Hormosira banksii (Fucales, Phaeophyta). Bot. Mar. 48, 198–207 (2005).
Stachowicz, J. J., Graham, M., Bracken, M. E. & Szoboszlai, A.I. Diversity enchances cover and stability of seaweed assemblages: the role of heterogeneity and time. Ecology 89, 3008–3019 (2008).
Mayer-Pinto, M., Underwood, A. J. & Marzinelli, E. M. The matrix influences direct and indirect effects of an anthropogenic disturbance on marine organisms. Environ. Res. 136, 15–20 (2015).
Orr, J. A. et al. Towards a unified study of multiple stressors: divisions and common goals across research disciplines. Proc. R. Soc. B Biol. Sci. 287, 20200421 (2020).
Tourneroche, A. et al. Bacterial–fungal Interactions in the kelp endomicrobiota drive autoinducer-2 quorum sensing. Front. Microbiol. 10, 1693 (2019).
Amend, A. S. et al. Fungi in the marine environment: open questions and unsolved problems. mBio 10, e01189–18 (2019).
Varga, J. J. et al. Antibiotics drive expansion of rare pathogens in a chronic infection microbiome model. mSphere 7, e00318–e00322 (2022).
Coyte, K. Z., Schluter, J. & Foster, K. R. The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666 (2015).
Foster, K. R., Schluter, J., Coyte, K. Z. & Rakoff-Nahoum, S. The evolution of the host microbiome as an ecosystem on a leash. Nature 548, 43–51 (2017).
Minich, J. J. et al. Elevated temperature drives kelp microbiome dysbiosis, while elevated carbon dioxide induces water microbiome disruption. PLOS ONE 13, e0192772 (2018).
Li, J. et al. Bacterial controlled mitigation of dysbiosis in a seaweed disease. ISME J. 16, 378–387 (2022).
Rao, D. et al. Low densities of epiphytic bacteria from the marine alga Ulva australis Inhibit settlement of fouling organisms. Appl. Environ. Microbiol. 73, 7844–7852 (2007).
Hoeger, A.-L., Jehmlich, N., Kipping, L., Griehl, C. & Noll, M. Associated bacterial microbiome responds opportunistic once algal host Scenedesmus vacuolatus is attacked by endoparasite amoeboaphelidium protococcarum. Sci. Rep. 12, 13187 (2022).
Marzinelli, E. M., Zagal, C. J., Chapman, M. G. & Underwood, A. J. Do modified habitats have direct or indirect effects on epifauna? Ecology 90, 2948–2955 (2009).
Straub, S. C. et al. Resistance, extinction, and everything in between—the diverse responses of seaweeds to marine heatwaves. Front. Mar. Sci. 6, 763 (2019).
Williams, S. L. & Dethier, M. N. High and dry: variation in net photosynthesis of the intertidal seaweed Fucus gardneri. Ecology 86, 2373–2379 (2005).
Lamote, M., Johnson, L. E. & Lemoine, Y. Interspecific differences in the response of juvenile stages to physical stress: fluorometric responses of fucoid embryos to variation in meteorological conditions 1. J. Phycol. 43, 1164–1176 (2007).
Ferreira, J. G., Arenas, F., Martínez, B., Hawkins, S. J. & Jenkins, S. R. Physiological response of fucoid algae to environmental stress: comparing range centre and southern populations. N. Phytol. 202, 1157–1172 (2014).
Ferrari, J. et al. Molecular analysis of a fungal disease in the habitat‐forming brown macroalga Phyllospora comosa (Fucales) along a latitudinal gradient. J. Phycol. 57, 1504–1516 (2021).
Klindworth, A. et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res 41, e1–e1 (2013).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
Edgar, R. C. Accuracy of microbial community diversity estimated by closed- and open-reference OTUs. PeerJ. 5, e3889 (2017).
Thijs, S. et al. Comparative evaluation of four bacteria-specific primer pairs for 16S rRNA gene surveys. Front. Microbiol. https://doi.org/10.3389/fmicb.2017.00494 (2017).
Nappi, J. et al. Differential priority effects impact taxonomy and functionality of host‐associated microbiomes. Mol. Ecol. https://doi.org/10.1111/mec.16336 (2022).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Oksanen, J. et al. Package ‘vegan’. Community Ecol. Package Version 2, 1–295 (2013).
Anderson, M. J. & Walsh, D. C. I. Permanova, anosim, and the Mantel test in the face of heterogeneous dispersions: what null hypothesis are you testing? Ecol. Monogr. 83, 557–574 (2013).
Wang, Y., Naumann, U., Wright, S. T. & Warton, D. I. mvabund—an R package for model-based analysis of multivariate abundance data: the mvabund R package. Methods Ecol. Evol. 3, 471–474 (2012).
Acknowledgements
We thank the volunteers who contributed to the fieldwork and experimental setup within the mesocosm. We would also like to thank the staff at the Sydney Institute of Marine Science for all of their help and expertise in the technical aspects of this work. This work was funded by the Australian Research Council through a Discovery Project to P.D.S., E.M.M. and S.K. (DP180104041), and the Ecological Society of Australia through a Holsworth Wildlife Research Endowment to A.H.M. and E.M.M. A.H.M. was further supported by a scholarship granted by the University of Sydney through the George William Michael Bequest.
Author information
Authors and Affiliations
Contributions
Study conception and design: A.H.M., K.L., E.M.M., P.D.S., S.E., S.K. Sample collection and processing, isolation and cultivation of bacterial isolates: A.H.M., K.L., G.W. Bioinformatic analyses and data visualisation: A.H.M. and S.V.G. Drafting of manuscript: A.H.M., K.L., G.W., S.E., S.V.G., P.S., E.M. All authors edited and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethics approval and consent to participate
Ethics approval was not required for the study.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
41522_2024_503_MOESM1_ESM.pdf
Supplementary information: Disentangling direct vs indirect effects of microbiome manipulations in a habitat-forming marine holobiont
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
McGrath, A.H., Lema, K., Egan, S. et al. Disentangling direct vs indirect effects of microbiome manipulations in a habitat-forming marine holobiont. npj Biofilms Microbiomes 10, 33 (2024). https://doi.org/10.1038/s41522-024-00503-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41522-024-00503-x