Microbial indicators of environmental perturbations in coral reef ecosystems
Microbiome volume 7, Article number: 94 (2019)
Coral reefs are facing unprecedented pressure on local and global scales. Sensitive and rapid markers for ecosystem stress are urgently needed to underpin effective management and restoration strategies. Although the fundamental contribution of microbes to the stability and functioning of coral reefs is widely recognised, it remains unclear how different reef microbiomes respond to environmental perturbations and whether microbiomes are sensitive enough to predict environmental anomalies that can lead to ecosystem stress. However, the lack of coral reef microbial baselines hinders our ability to study the link between shifts in microbiomes and ecosystem stress. In this study, we established a comprehensive microbial reference database for selected Great Barrier Reef sites to assess the diagnostic value of multiple free-living and host-associated reef microbiomes to infer the environmental state of coral reef ecosystems.
A comprehensive microbial reference database, originating from multiple coral reef microbiomes (i.e. seawater, sediment, corals, sponges and macroalgae), was generated by 16S rRNA gene sequencing for 381 samples collected over the course of 16 months. By coupling this database to environmental parameters, we showed that the seawater microbiome has the greatest diagnostic value to infer shifts in the surrounding reef environment. In fact, 56% of the observed compositional variation in the microbiome was explained by environmental parameters, and temporal successions in the seawater microbiome were characterised by uniform community assembly patterns. Host-associated microbiomes, in contrast, were five-times less responsive to the environment and their community assembly patterns were generally less uniform. By applying a suite of indicator value and machine learning approaches, we further showed that seawater microbial community data provide an accurate prediction of temperature and eutrophication state (i.e. chlorophyll concentration and turbidity).
Our results reveal that free-living microbial communities have a high potential to infer environmental parameters due to their environmental sensitivity and predictability. This highlights the diagnostic value of microorganisms and illustrates how long-term coral reef monitoring initiatives could be enhanced by incorporating assessments of microbial communities in seawater. We therefore recommend timely integration of microbial sampling into current coral reef monitoring initiatives.
Coral reef ecosystems are rapidly degrading due to local and global pressures . Overfishing, pollution, declining water quality, disease and outbreaks of coral predating crown-of-thorns starfish are responsible for localised reef degradation  while climate change is impacting reefs on a global scale, including remote reefs with little local anthropogenic pressure . For example, elevated sea surface temperatures caused back-to-back coral mass bleaching events in 2016 and 2017, resulting in a significant loss of shallow-water corals on the Great Barrier Reef (GBR) . Climate conditions predicted for the end of the century will result in even more frequent and severe coral mass bleaching events with dire projections for the future of coral reefs [5, 6]. This global coral reef crisis is driving the development of new management, reef restoration and bioengineering tools to counteract reef loss and ensure the persistence of coral reefs [7, 8]. Early prediction of ecosystem stress is critical for an effective implementation of local management and restoration strategies on threatened reef sites.
Microorganisms have considerable potential as a monitoring tool for coral reef ecosystem health [9,10,11]. Microorganisms are fundamental drivers of biogeochemical cycling on coral reefs [12,13,14]; they form intimate associations with the coral reef benthos [15,16,17], and they contribute significantly to host health and ecosystem homeostasis [18,19,20]. The constant amendment of microbial communities to exploit available resources  can trigger differential abundances of specific microorganisms; hence, shifts in community composition can provide an early indication of environmental change . For example, compositional and functional shifts of coral-associated microbial communities have been described along gradients of anthropogenic impact [23,24,25] and with changes in water quality . However, despite having many of the useful characteristics required of environmental indicators [9, 27], the diagnostic potential of microorganisms for coral reef monitoring is largely conceptual, with only a few studies elaborating on their potential value. For example, the ‘microbialisation score’ measures human impacts on coral reefs based on the ratio of microbial and fish metabolic rates . The main limitations to further develop and apply microbial-based monitoring approaches are the lack of temporal and spatial baselines for coral reef microbiomes [9, 29].
Coral reefs comprise a complex network of free-living and host-associated microbial communities with strong benthic-pelagic exchange [13, 30]. Therefore, holistic assessments that combine different reef hosts and habitats are required to better understand microbial dynamics and sensitivities to environmental perturbations. The diagnostic value of microbial-based monitoring is likely to vary between distinct habitats of a coral reef ecosystem. For example, microbial communities occurring in seawater may be directly affected by the quality of the ambient reef water or climate conditions; however, the high heterogeneity of seawater due to local hot-spots of available resources [31, 32] may diminish the specificity of these communities. In contrast, microbial communities that dwell in corals live in tight association with the most important frame-builders of reefs  and hence may provide crucial information not only on the environmental conditions but also on the effect of the environment on the coral host itself. Sponges, a highly abundant and diverse component of coral reefs , are renowned for their enormous filtration capacity  and form diverse and intimate associations with microbial communities . Hence, sponge microbiomes may provide suitable indicators to monitor water quality. Host-associated biofilms, such as those inhabiting the mucus layer of corals and the surface of macroalgae, provide another potential niche habitat informative for microbial indicators of environmental state. Coral mucus, for example, has been described as a suitable habitat to screen for enterobacteria from sewage contamination due to its ability to trap bacteria .
Given the complexity of microbial life on coral reefs, we sought to identify the most suitable reef microbiomes for a microbial indicator program to pinpoint environmental state. To do this, we quantified the (1) habitat-specificity, (2) determinacy of microbial community successions and (3) sensitivity towards environmental parameters of multiple free-living and host-associated microbiomes. Subsequently, we tested the microbiome’s ability to infer environmental state using indicator value  and machine learning approaches .
Samples were collected during a 16-month period (February 2016–May 2017), at monthly (Magnetic Island—Geoffrey Bay) and periodic (Orpheus Island—Pioneer Bay and Channel) intervals (Additional file 1: Table S1). The bacterial 16S rRNA genes of 381 samples including seawater, sediment, sponge tissue (Coscinoderma matthewsi and Amphimedon queenslandica), coral tissue and mucus (Acropora tenuis and Acropora millepora) and macroalgal surfaces (Sargassum sp.) were sequenced (Fig. 1). In total 231,316 zero-radius operational taxonomic units (zOTUs) were identified based on 100% sequence similarity .
Coral reef microbiomes are habitat-specific
Habitat-specificity of coral reef microbes was assessed by comparing the similarities of microbial communities associated with seawater (n = 48), sediment (n = 48), A. queenslandica (n = 30), C. matthewsi (n = 42), A. tenuis (tissue n = 48, mucus n = 46), A. millepora (tissue n = 42, mucus n = 42) and Sargassum sp. (n = 35). Non-metric multidimensional scaling based on Bray-Curtis dissimilarities revealed a clear separation of the microbial communities from different reef habitats (Fig. 1), and habitat-specificity was further confirmed with permutational multivariate analysis of variance (PERMANOVA, p = 9.999 × 10−5, Additional file 1: Tables S2-S3). Furthermore, alpha diversities (ANOVA, F(8/372) = 142, p = 2 × 10−16) and zOTU richness (ANOVA, F(8/372) = 369, p = 2 × 10−16) varied significantly between reef habitats (Additional file 1: Figure S1 and Tables S4-S6). Sediment harboured by far the most diverse (Shannon Index 7.4 ± 0.2 SD) bacterial community, although microbial diversity was also high in coral surface mucus (Shannon Index 5.1 ± 0.9 SD), macroalgal biofilms (Shannon Index 4.5 ± 1.4 SD), seawater (Shannon Index 4.4 ± 0.2 SD) and in the tissue of the sponge C. matthewsi (Shannon Index 4.4 ± 0.3 SD). Microbial diversity was lowest in coral tissue (Shannon Index 3.3 ± 0.8 SD) and in the sponge A. queenslandica (Shannon Index 2.7 ± 0.8 SD). These results suggest overall high habitat-specificity of free-living and host-associated microbial communities within coral reef ecosystems.
Uniform vs variable community assembly patterns
The uniformity versus variability of microbial community assembly patterns was explored through comparison of compositional similarity (Bray-Curtis Similarity Index, 0 = dissimilar, 1 = identical) in samples collected monthly at Geoffrey Bay (Magnetic Island). The microbial communities of seawater (n = 30, Wilcoxon Rank-Sum test p = 3.1 × 10−7) and sediment (n = 30; Wilcoxon Rank-Sum test p = 3 × 10−5) had significantly higher similarities ‘within’ than ‘between’ sampling events (Fig. 2a). This uniform response of the free-living microbial communities suggests that deterministic rather than stochastic processes drive their community assembly. For host-associated microbiomes, the overall response pattern varied between species. Microbial communities associated with the sponge C. matthewsi (n = 27; Wilcoxon Rank-Sum test, p = 0.0076), the coral A. tenuis (mucus n = 28, tissue n = 30; Wilcoxon Rank-Sum test, p = 0.0041 and p = 0.0096, respectively) and the macroalga Sargassum sp. (n = 30; Wilcoxon Rank-Sum test, p = 0.00013) followed the same trend as the free-living communities, with significantly higher similarities ‘within’ than ‘between’ sampling events (Fig. 2a). In contrast, the microbiome of the sponge A. queenslandica (n = 30; Wilcoxon Rank-Sum test, p = 0.23) and the coral A. millepora (mucus n = 24, tissue n = 24; Wilcoxon Rank-Sum test, p = 0.15 and p = 0.11 respectively) showed no significant difference in similarities ‘within’ and ‘between’ time points (Fig. 2a). Analysis of the compositional similarity of sample replicates within each sampling time point indicated that the seawater microbial communities not only exhibit an overall higher similarity ‘within’ replicates, but the high compositional similarity is conserved across all sampling events (Fig. 2b). In contrast, host-associated microbial communities showed a generally lower compositional similarity and higher variation between sample replicates within each sampling time point (Fig. 2b).
Trends in the temporal community assembly pattern of free-living, host tissue- and biofilm-associated microbial communities were analysed using analysis of similarity (ANOSIM) as a proxy to describe similarity patterns (R = 0 indicates equal similarity ‘within’ and ‘between’ time point replicates and R = 1 indicates higher ‘within’ than ‘between’ sampling time point similarities; Fig. 2b and Additional file 1: Figure S2). Overall, free-living microbiomes had R values closer to 1 (seawater R = 0.9919 and sediment R = 0.7322), whereas host-associated microbiomes had R values closer to 0 (A. queenslandica R = 0.2927, C. matthewsi R = 0.3449, A. tenuis tissue R = 0.4547, A. millepora tissue R = 0.2151, A. tenuis mucus R = 0.4613, A. millepora mucus R = 0.3090 and Sargassum sp. biofilm R = 0.4440; Fig. 2b and Additional file 1: Figure S2). These results suggest that free-living microbiomes (seawater and sediment) exhibit a uniform compositional succession, whereas host-associated microbiomes (coral, sponge and macroalgae) are more stochastic in their temporal community succession. The uniform temporal response of free-living microbiomes suggests a high diagnostic value of these microbial communities; hence seawater and sediment microbiomes should provide an accurate prediction of environmental variables.
Microbiomes in seawater (n = 48) and sediment (n = 48) were further tested for their compositional similarity between all three sampling sites (Geoffrey Bay, Pioneer Bay and Channel). The microbial community composition of sediment samples varied significantly between sampling sites (ANOSIM R = 0.9430, p = 0.001, Additional file 1: Figure S3a). The seawater microbiome, in contrast, showed high temporal variability (ANOSIM R = 0.9934, p = 0.001) and low spatial variability (ANOSIM R = 0.2343, p = 0.002; Additional file 1: Figure S3b). The high spatial variability of sediment microbiomes indicates that habitat characteristics rather than environmental fluctuations are the main drivers structuring community composition.
Environmental sensitivity of the different microbiomes was assessed by comparing how much of the compositional variation was explained by sea surface temperature, light and water quality parameters (Additional file 1: Figures S4 and S5). The compositional variability of the seawater microbiome (n = 30) was significantly explained by sampling date, season (summer versus winter) and water quality parameters, such as average seawater temperature, average hours of daylight, total suspended solids (TSS), particulate organic carbon (POC), chlorophyll a (Chl a) and non-purgeable organic carbon (NPOC) concentration (permutational ANOVA for Bray Curtis distance-based redundancy analysis (dbRDA); Fig. 3a and Additional file 1: Table S7a-b). In total, these environmental parameters explained 56% of the observed compositional variation in seawater (variation partitioning analysis, Fig. 3b, Additional file 1: Table S7). Season (summer versus winter) and sampling date solely explained 6% and 4%, respectively (variation partitioning analysis, Fig. 3b). In comparison, sampling site significantly explained 24% of the variation in sediment microbial communities (n = 48), which overlapped by 12% with the variation explained by sediment characteristics, such as particle size and total organic carbon (TOC) content (permutational ANOVA for dbRDA and variation partitioning analysis; Additional file 1: Tables S7b and S8). Water quality parameters and sea surface temperature explained only 3% of the observed variability in the sediment microbiome (variation partitioning analysis).
Host-associated microbiomes varied substantially in their response to environmental parameters (permutational ANOVA for dbRDA and variation partitioning analysis, Fig. 3b, c, Additional file 1: Tables S7c-i and S8). On average, 11% of the observed community variations in host-associated microbiomes were explained by the environment (variation partitioning analysis), which is five times less than what we found for the seawater-associated microbial community (Additional file 1: Table S8). This suggests that compositional variations of the seawater microbiome are more likely to reflect environmental changes. Host-associated microbiomes are comparatively stable to changes in environmental factors.
Predictability of environmental metadata
Due to the seawater microbiomes uniform temporal pattern and high sensitivity to changing environmental parameters, the ability to infer environmental state based on microbial community data was tested using an indicator value analysis  and a random forest machine learning approach. In total, 110 zOTUs were identified as significant indicators for temperature (indicator value p < 0.01). Microbial zOTU assemblages that were indicative of high, low and average seawater temperatures (classification based on their variation around observed annual averages) were present throughout the sampling period. However, higher relative abundances and lower variation (as calculated by coefficient of variation) were evident at certain time points (Fig. 4a). Furthermore, we were able to identify microbial indicator taxa for high and low Chl a, TSS and POC levels (Additional file 1: Figure S6). Indicators for low and high seawater temperatures were identified in the bacterial phyla Proteobacteria, Bacteroidetes, Cyanobacteria, Actinobacteria and Planctomycetes (Fig. 4b). High temperatures were indicated by an increase of zOTUs belonging to the bacterial family Rhodobacteraceae and the presence of Cryomorphaceae, Synechococcaeae, Vibrio and Flavobacterium (Fig. 4b). In contrast, the occurrence of zOTUS belonging to the family Pelagibacteriaceae and the genus Prochlorococcus were indicative for low seawater temperatures. The phyla Proteobacteria, Bacteroidetes and Cyanobacteria had the greatest number of indicator zOTUs for temperature and other water quality parameters (Additional file 1: Figure S6). Flavobacteriaceae-affiliated zOTUs were significant indicators for temperature, Chl a, TSS and POC. Halomonadaceae significantly associated with high Chl a and TSS and zOTUs belonging to the phylum Verrucomicrobia were significant indicators for high TSS levels.
The diagnostic value of the seawater microbiome (n = 48) was further evaluated by applying a random forest machine learning classification and regression analysis with 1213 zOTUs preselected based on a non-zero abundance threshold in at least 10% of the samples (n = 48). The seawater microbiome enabled the prediction of seawater temperature classes (low, average, high) with 92% accuracy (Kappa = 88%, Fig. 5a, b and Additional file 1: Figure S7). Highest accuracy (lowest out of bag (OOB) estimated error rate) was achieved with mtry = 100 zOTUS. Random forest regression of the seawater microbiome predicted temperature values (R2 = 0.67, RMSE = 0.5) (Fig. 5c, d and Additional file 1: Figure S8) with the highest accuracy (lowest OOB estimated error rate) when mtry = 400 zOTUs. The effectiveness of zOTUs in reducing uncertainty and variance (also referred to as ‘feature importance’) within the machine learning algorithm was measured by the decrease in mean accuracy for classification and mean-squared error (% Inc. MSE) for regression. The most important zOTUs belong to the bacterial taxa Flavobacteriaceae, Pelagibacteraceae, Cyanobacteria, Rhodobacteraceae, Synechococcaceae and Pirrelulaceae. These results demonstrate that the microbial community associated with coral reef seawater allows for the accurate prediction of fluctuations in sea surface temperature and water quality parameters.
Sensitive and rapidly responding markers of coral ecosystem stress are needed to underpin effective management and restoration strategies. In this study, we used a range of statistical tests and machine learning approaches across multiple free-living and host-associated reef microbiomes to assess their diagnostic value as sensitive indicators of environmental state. Our results show that the microbial community in reef seawater has the highest diagnostic value when compared to other free-living (e.g. sediment) and host-associated microbiomes (e.g. coral, sponge and macroalgae). Our conclusion is based on the microbiome’s (1) habitat-specificity, (2) uniformity of its community assembly, (3) sensitivity towards environmental fluctuations and (4) accuracy to predict environmental parameters. This assessment of the diagnostic capacity of various free-living and host-associated coral reef microbiomes to extrapolate environmental variations provides crucial information for ecosystem management initiatives aimed at incorporating microbial monitoring.
In general, high habitat-specificity was observed across free-living and host-associated microbiomes, confirming previous reports on the compositional variability of microbial communities between coral reef habitats , host species [15, 41,42,43] and even between host compartments . High compositional divergence of microbial communities across different reef habitats can be due to the variation of available resources and/or biotic interactions . High habitat-specificity contributes to the overall high diversity and complexity across different microbial communities on coral reefs, highlighting the importance of holistic studies that focus on microbial interactions across the benthic-pelagic realm.
Bacterial community structure associated with water and sediment is thought to be primarily governed by deterministic processes . Our results are consistent with this, showing uniform community assembly patterns within time point replicates. In contrast, host-associated microbiomes displayed little compositional similarity within a sampling time point, suggesting a non-uniform temporal response. Host-associated microbiomes were also only marginally affected by environmental parameters, indicating that their community assembly pattern is variable between conspecific individuals . A higher variability in community assembly can lead to increased community heterogeneity, also referred to as dispersion, which has been described as a common characteristic of host-associated microbiomes [18, 46,47,48]. Furthermore, lower microbial compositional similarities among replicates may be driven by increased niche space (e.g. host compartments)  and host genotype effects (e.g. host genetics) . Collectively, our results show that free-living microbial communities have a higher potential to infer environmental parameters (such as standard measures in environmental monitoring programs) than host-associated microbial communities due to their higher uniformity and environmental sensitivity. Importantly however, previous metaproteomic research on reef sponges has shown that while microbial community composition can appear stable when seawater temperatures increase, disruption to nutritional interdependence and molecular interactions (such as reduced expression of transporters involved in the uptake of sugars, peptides and other substrates) actually occurs prior to detectable changes in community structure . Hence, considering the importance of microbes to reef invertebrate health, more sensitive transcriptomic/proteomic approaches may still be warranted for sensitive detection of microbial responses to environmental perturbations.
The diagnostic potential of microbial communities, especially in combination with machine learning approaches, has gained momentum across multiple research fields, including disease identification by characterisation of the human gut-microbiome , evaluation of the environment and host genetics on the human microbiome , prediction of hydrological functions in riverine ecosystems  and assessment of macroecological patterns in soil samples . This development of microbial-based diagnostics is largely due to availability of high-throughput sequencing of the 16S rRNA gene and streamlined analytical pipelines that facilitate rapid assessment of microbial community composition [54, 55]. In addition to its utility for inferring environmental fluctuations, the seawater microbiome possesses numerous characteristics desirable for environmental monitoring programs: (i) non-destructive collection and simple processing methods facilitate large-scale collections alongside existing programs that sample water quality measurements, (ii) high fractional contribution of abundant microbes minimises the impacts of sequencing biases (Additional file 1: Figure S9) and (iii) sampling is conducive to future automated, high throughput analyses such as in-line flow cytometry on vessels and real-time DNA/RNA sequencing for community characterisation.
Incorporation of seawater microbial community data into coral reef monitoring approaches should enhance our ability to describe environmental conditions and changes more holistically. For example, temperature fluctuations drive structural variations in seawater microbial communities [56, 57], and elevated seawater temperatures on coral reefs are highly correlated with coral bleaching [1, 58]. The inclusion of microbial community data alongside water quality parameters could therefore improve our ability to predict the likelihood of ecosystem stress. For instance, our sample sites, located in the central sector of the GBR, were not affected by the 2016 bleaching that primarily affected the northern sector ; however, they were impacted by the 2017 bleaching event . In the months prior to bleaching (late December 2016 till March 2017), we observed two to four times higher relative abundances of high temperature indicator assemblages than when compared to the equivalent period at the beginning of 2016 (Fig. 4a), where no bleaching was observed. Interestingly, high temperature indicator assemblages included putative coral pathogens (e.g. Vibrio) and opportunistic bacteria (e.g. Rhodobacteraceae, Verrucomicrobia and Flavobacterium). Coral pathogens, such as Vibrio corallilyticus, increase their efficiency and motility behaviours with rising seawater temperatures [61,62,63], and the higher abundance of these microbes may explain the increased prevalence of coral disease post bleaching . Hence, microbial monitoring could help inform managers about impending disease outbreaks.
Our study provides the first holistic microbial baseline spanning multiple free-living and host-associated microbiomes for selected GBR sites. Results suggest that there is realistic scope to enhance long-term reef monitoring initiatives by incorporating seawater microbiome observations for assessments of environmental change over space and time, especially for rapid and sensitive identification of early signs of declining ecosystem health. The establishment of microbial observatories  and DNA biobanks for long-term biomonitoring  will be paramount to successfully inferring ecosystem state and/or perturbations from microbial communities. We therefore recommend timely integration of microbial sampling into current coral reef monitoring initiatives. Further refinement of the sampling and data analysis techniques should focus on selection and validation of additional indicator taxa as well as assessment of ecologically important microbial functions. A further consideration is to explore which monitoring objectives would benefit most from assessments of microbial communities. For example, it is likely that the rapid response time of microbial indicators makes them better suited to early-warning, impact or compliance monitoring programs than to monitoring of slower, long-term changes.
Samples for microbial community characterisation were collected monthly (Magnetic Island) and periodically (Orpheus Island) from seawater, sediment and multiple host organisms (i.e. corals, sponges and macroalgae), along with environmental metadata, between February 2016 and May 2017 at three Great Barrier Reef sites (Fig. 1). Samples were collected under the permit G16/38348.1 issued by the Great Barrier Reef Marine Park Authority.
Samples (n = 3/sample type/sampling event) for molecular analysis and additional environmental metadata were collected following the standard operational procedures of the Australian Marine Microbial Biodiversity Initiative (AMMBI; https://data.bioplatforms.com/organization/pages/australian-microbiome/methods). In brief, seawater for molecular analysis was collected with collapsible sterile bags close to the reef substrate at 2 m depth and pre-filtered (50 μm) to remove large particles and subsequently filtered (2 L) onto 0.2 μm Sterivex-filters (Millepore). The sediment surface layer was sampled with sterile 50 mL tubes at 2 m depth and subsampled immediately into 2 mL cryogenic vials. The sponges Coscinoderma matthewsi and Amphimedon queenslandica were removed from the substrate (at 7 m and 3 m respectively) with sterile scalpel blades, rinsed with 0.2 μm filter-sterilised seawater and subsampled into 2 mL cryogenic vials. The surface mucus layer of the two acroporid coral species, Acropora tenuis and Acropora millepora, was sampled with sterile cotton swabs . Additionally, coral fragments of each sampled coral were collected at 3 m depth. Coral fragments were rinsed with 0.2 μm filtered-sterilised seawater and placed into 5 mL cryogenic vials. The thallus (including stem, floats and blades) of the macroalgae Sargassum sp. was sampled with sterile scalpels at 3 m depth, rinsed with 0.2 μm filtered-sterilised seawater and placed into 2 mL cryogenic vials. All samples were immediately flash frozen in liquid nitrogen after processing and stored at − 80 °C until DNA extraction.
Additional seawater samples were collected with a diver-operated Niskin bottle close to the reef substrate at 2 m depth at each sampling occasion. Water was subsampled in duplicate for analyses of salinity and concentrations of dissolved organic carbon (DOC), dissolved inorganic carbon (DIC), particulate organic carbon (POC), dissolved inorganic nutrients (DIN), total suspended solids (TSS) and chlorophyll a (Chl a) concentration. Samples were further analysed according to the standard procedures of the Australian Institute of Marine Science (AIMS, Townsville, Australia) . Sediment samples were collected with 100 mL glass jars at 2 m depth and characteristics, such as grain size distribution and total organic carbon (TOC) and nitrogen (TON) content, were assessed for each sampling event. Seawater temperatures were obtained from AIMS long-term monitoring temperature records (http://eatlas.org.au/).
Prior to extraction, the macroalgal biofilm was separated from the algal tissue by overnight incubation at 200 rpm in 10 mL 1x PBS at 37 °C. Coral fragments were defrosted on ice and the tissue was stripped from the skeleton with an airgun into 1x PBS solution, homogenised for 1 min at 12.5 rpm with a tissue homogeniser, pelleted (10 min at 16,000 rcf) and snap frozen in liquid nitrogen prior to DNA extraction. DNA from seawater, sediment, sponge and macroalgal biofilms was extracted with the DNeasy PowerSoil kit (Qiagen) and DNA of coral tissue and mucus samples was extracted using the DNeasy PowerBiofilm kit (Qiagen) following the Manufacturer’s instructions. DNA extracts were stored at − 80 °C until being sent for sequencing.
16S rRNA gene sequencing
DNA extracts were sent on dry ice to the Ramaciotti Centre for Genomics (Sydney, Australia) for sequencing. The bacterial 16S rRNA genes were sequenced using the 27F  and 519R  primer pairs on the Illumina MiSeq platform utilising a duel indexed 2 × 300 bp paired end approach. Further documentation outlining the standard operating procedures for generating and sequencing amplicons is available at https://data.bioplatforms.com/dataset/marine-microbes-methods.
Sequencing data were analysed as single nucleotide variants in a standardised platform alongside other Australian microbial biodiversity initiative samples [39, 70]. In brief, forward and reverse reads were merged using FLASH . FASTA formatted sequences were extracted from FASTQ files and those < 400 bp in length or containing N’s or homopolymer runs of > 8 bp were removed using MOTHUR (v1.34.1) . USEARCH (64 bit v10.0.240)  package was used to de-replicate sequences and to order them by abundance. Sequences with < 4 representatives and Chimeras were removed. Quality-filtered sequences were mapped to chimera-free zero-radius operational taxonomic units (zOTUs) and a sample by read abundance table created. zOTUs were taxonomically classified with SILVA v132  database using MOTHUR’s implementation of the Wang classifier  and a 60% Bayesian probability cut-off.
Chloroplast and mitochondria-derived reads as well as singletons were removed from the dataset. Remaining data were rarefied to 3600 reads per sample and transformed to relative abundances using the phyloseq package  in R .
Habitat and host-specificity
Habitat and host-specificity of a microbiome was assessed by calculating the compositional similarities of all 381 samples with the Bray-Curtis Similarity Index and illustrating them in a non-metric multidimensional scaling (NMDS) plot using the phyloseq package . To confirm habitat and host-specificity, permutational multivariate analysis of variance (PERMANOVA) was applied using the adonis() function of the vegan package  with 10,000 permutations.
Uniform response pattern
The microbiome similarity of replicates for sampling time points versus the microbiome similarity between sampling time points was compared by obtaining the Bray-Curtis Similarity for each habitat individually. The variation between the overall within and between time point replicates was tested with a Wilcoxon Rank-sum test in R . The dispersion of the Bray-Curtis Similarities within a sampling time point was calculated as the coefficient of variation (ratio of the standard deviation to the mean expressed as a percentage). The higher the coefficient of variation, the higher the variability in the microbiome composition among replicates of a time point. Analysis of similarity (ANOSIM; anosim() function of the vegan package ) based on Bray-Curtis Similarities was used to further evaluate within and between time point similarities in the microbial communities.
Environmental metadata were z-score standardised  and checked for collinearity using the Pearson correlation coefficient. Collinearity was assumed if correlation was > 0.7 or < − 0.7 . Collinear variables were considered redundant and removed from the analysis.
zOTU relative abundance, environmental metadata (e.g. average seawater temperature, average hours of daylight, Chl a, POC, NPOC and TSS concentration), season (summer versus winter) and sampling date were used for Bray-Curtis distance-based redundancy analysis (dbRDA) using the phyloseq package . The significance of each response variable was confirmed with an analysis of variance (ANOVA) for the dbRDA (anova.cca() function in the vegan package ). Only significant (p value < 0.05) response variables were kept in the model. The explanatory value (in %) of significant response variables (e.g. environmental parameters, season and sampling date) was assessed with a variation partitioning analysis of the vegan package .
Indicator value analysis
Indicator taxa were identified with the indicator value analysis (indicspecies package ) using the following thresholds: 1000 permutations, minimum specificity (At) and minimum sensitivity (Bt) set to 70% and p value ≤ 0.01.
Random forest machine learning
Random forest machine learning was performed with the caret  and random forest package  in R . zOTUs with non-zero abundance values in at least 10% of the samples (n = 48) were preselected and z-score standardised prior to model training. Random forest (with ntrees = 10,000) prediction error was measured with out-of-bag (OOB) error. Highest accuracy (lowest OOB estimated error rate) for classification was achieved with mtry = 100 zOTUS and for regression with mtry = 400 zOTUs. Importance of zOTUs was measured using the decrease in mean accuracy for classification and mean-squared error (% Inc. MSE) for regression.
Hughes TP, Barnes ML, Bellwood DR, Cinner JE, Cumming GS, Jackson JBC, Kleypas J, van de Leemput IA, Lough JM, Morrison TH, et al. Coral reefs in the Anthropocene. Nature. 2017;546:82–90.
De'ath G, Fabricius KE, Sweatman H, Puotinen M. The 27-year decline of coral cover on the Great Barrier Reef and its causes. Proc Natl Acad Sci U S A. 2012;109:17995–9.
Hoegh-Guldberg O, Mumby PJ, Hooten AJ, Steneck RS, Greenfield P, Gomez E, Harvell CD, Sale PF, Edwards AJ, Caldeira K, et al. Coral reefs under rapid climate change and ocean acidification. Science. 2007;318:1737–42.
Hughes TP, Kerry JT, Baird AH, Connolly SR, Dietzel A, Eakin CM, Heron SF, Hoey AS, Hoogenboom MO, Liu G, et al. Global warming transforms coral reef assemblages. Nature. 2018;556:492–6.
van Hooidonk R, Maynard J, Tamelander J, Gove J, Ahmadia G, Raymundo L, Williams G, Heron SF, Planes S. Local-scale projections of coral reef futures and implications of the Paris Agreement. Sci Rep. 2016;6:39666.
Hughes TP, Anderson KD, Connolly SR, Heron SF, Kerry JT, Lough JM, Baird AH, Baum JK, Berumen ML, Bridge TC, et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science. 2018;359:80–3.
Anthony K, Bay LK, Costanza R, Firn J, Gunn J, Harrison P, Heyward A, Lundgren P, Mead D, Moore T, et al. New interventions are needed to save coral reefs. Nat Ecol Evol. 2017;1:1420–2.
Damjanovic K, Blackall LL, Webster NS, van Oppen MJH. The contribution of microbial biotechnology to mitigating coral reef degradation. Microb Biotechnol. 2017;10:1236–43.
Glasl B, Webster NS, Bourne DG. Microbial indicators as a diagnostic tool for assessing water quality and climate stress in coral reef ecosystems. Mar Biol. 2017;164:91.
Glasl B, Bourne DG, Frade PR, Webster NS. Establishing microbial baselines to identify indicators of coral reef health. Microbiol Aust. 2018;39:42–6.
Roitman S, Joseph Pollock F, Medina M. Coral microbiomes as bioindicators of reef health. Cham: Springer International Publishing; 2018. p. 1–19.
Gast GJ, Wiegman S, Wieringa E, FCv D, Bak RPM. Bacteria in coral reef water types: removal of cells, stimulation of growth and mineralization. Mar Ecol Prog Ser. 1998;167:37–45.
Bourne DG, Webster NS. Coral Reef Bacterial Communities. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes. Berlin, Heidelberg: Springer; 2013. p. 163–87
Sorokin YI. Trophical role of bacteria in ecosystem of coral reef. Nature. 1973;242:415–7.
Rohwer F, Seguritan V, Azam F, Knowlton N. Diversity and distribution of coral-associated bacteria. Mar Ecol Prog Ser. 2002;243:1–10.
Webster NS, Luter HM, Soo RM, Botte ES, Simister RL, Abdo D, Whalan S. Same, same but different: symbiotic bacterial associations in GBR sponges. Front Microbiol. 2012;3:444.
Egan S, Harder T, Burke C, Steinberg P, Kjelleberg S, Thomas T. The seaweed holobiont: understanding seaweed–bacteria interactions. FEMS Microbiol Rev. 2013;37:462–76.
Glasl B, Herndl GJ, Frade PR. The microbiome of coral surface mucus has a key role in mediating holobiont health and survival upon disturbance. ISME J. 2016;10(9):2280–92.
Hentschel U, Schmid M, Wagner M, Fieseler L, Gernert C, Hacker J. Isolation and phylogenetic analysis of bacteria with antimicrobial activities from the Mediterranean sponges Aplysina aerophoba and Aplysina cavernicola. FEMS Microbiol Ecol. 2001;35:305–12.
Webster NS, Reusch TBH. Microbial contributions to the persistence of coral reefs. ISME J. 2017;11:2167–74.
Martiny JB, Jones SE, Lennon JT, Martiny AC. Microbiomes in light of traits: a phylogenetic perspective. Science. 2015;350:aac9323.
Garza DR, van Verk MC, Huynen MA, Dutilh BE. Towards predicting the environmental metabolome from metagenomics with a mechanistic model. Nat Microbiol. 2018;3:456–60.
Ziegler M, Roik A, Porter A, Zubier K, Mudarris MS, Ormond R, Voolstra CR. Coral microbial community dynamics in response to anthropogenic impacts near a major city in the Central Red Sea. Mar Pollut Bull. 2016;105(2):629–40.
Kelly LW, Williams GJ, Barott KL, Carlson CA, Dinsdale EA, Edwards RA, Haas AF, Haynes M, Lim YW, McDole T, et al. Local genomic adaptation of coral reef-associated microbiomes to gradients of natural variability and anthropogenic stressors. Proc Natl Acad Sci. 2014;111:10227–32.
Dinsdale EA, Pantos O, Smriga S, Edwards RA, Angly F, Wegley L, Hatay M, Hall D, Brown E, Haynes M, et al. Microbial ecology of four coral atolls in the northern Line Islands. PLoS One. 2008;3:e 1584.
Angly FE, Heath C, Morgan TC, Tonin H, Rich V, Schaffelke B, Bourne DG, Tyson GW. Marine microbial communities of the great barrier reef lagoon are influenced by riverine floodwaters and seasonal weather events. PeerJ. 2016;4:e1511.
Cooper TF, Gilmour JP, Fabricius KE. Bioindicators of changes in water quality on coral reefs: review and recommendations for monitoring programmes. Coral Reefs. 2009;28:589–606.
McDole T, Nulton J, Barott KL, Felts B, Hand C, Hatay M, Lee H, Nadon MO, Nosrat B, Salamon P, et al. Assessing coral reefs on a pacific-wide scale using the microbialization score. PLoS One. 2012;7:e43233.
Bourne DG, Morrow KM, Webster NS. Coral Holobionts: insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Ann Rev Microbiol. 2016;70:317–40.
Lesser MP. Benthic-pelagic coupling on coral reefs: feeding and growth of Caribbean sponges. J Exp Mar Biol Ecol. 2006;328:277–88.
Azam F. Microbial control of oceanic carbon flux: the plot thickens. Science. 1998;280:694–6.
Stocker R. Marine microbes see a sea of gradients. Science. 2012;338:628–33.
Diaz MC, Rützler K. Sponges: an essential component of Caribbean coral reefs. Bull Mar Sci. 2001;69:535–46.
Reiswig HM. In situ pumping activities of tropical Demospongiae. Mar Biol. 1971;9:38–50.
Taylor MW, Radax R, Steger D, Wagner M. Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol Mol Biol Rev. 2007;71:295–347.
Lipp EK, Griffin DW. Analysis of coral mucus as an improved medium for detection of enteric microbes and for determining patterns of sewage contamination in reef environments. EcoHealth. 2004;1:317–23.
De Cáceres M, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90:3566–74.
Knights D, Costello EK, Knight R. Supervised classification of human microbiota. FEMS Microbiol Rev. 2011;35:343–59.
Brown MV, Jvd K, Ostrowski M, Seymour JR, Ingleton T, Messer LF, Jeffries T, Siboni N, Laverock B, Bibiloni-Isaksson J, et al. Systematic, continental scale temporal monitoring of marine pelagic microbiota by the Australian marine microbial biodiversity initiative. Sci Data. 2018;5:180130.
Tout J, Jeffries TC, Webster NS, Stocker R, Ralph PJ, Seymour JR. Variability in microbial community composition and function between different niches within a coral reef. Microb Ecol. 2014;67:540–52.
Carlos C, Torres TT, Ottoboni LMM. Bacterial communities and species-specific associations with the mucus of Brazilian coral species. Sci Rep. 2013;3:1624.
Glasl B, Smith CE, Bourne DG, Webster NS. Exploring the diversity-stability paradigm using sponge microbial communities. Sci Rep. 2018;8:8425.
Webster NS, Thomas T. The sponge hologenome. MBio. 2016;7:e00135.
Sweet MJ, Croquer A, Bythell JC. Bacterial assemblages differ between compartments within the coral holobiont. Coral Reefs. 2011;30:39–52.
Wang J, Shen J, Wu Y, Tu C, Soininen J, Stegen JC, He J, Liu X, Zhang L, Zhang E. Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes. ISME J. 2013;7:1310–21.
Casey JM, Connolly SR, Ainsworth TD. Coral transplantation triggers shift in microbiome and promotion of coral disease associated potential pathogens. Sci Rep. 2015;5:11903.
Zaneveld JR, McMinds R, Thurber RV. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol. 2017;2:17121.
Zaneveld JR, Burkepile DE, Shantz AA, Pritchard CE, McMinds R, Payet JP, RoryWelsh CAMS, Lemoine NP, Rosales S, et al. Overfishing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales. Nat Commun. 2016;7:11833.
Fan L, Liu M, Simister R, Webster NS, Thomas T. Marine microbial symbiosis heats up: the phylogenetic and functional response of a sponge holobiont to thermal stress. ISME J. 2013;7:991–1002.
Duvallet C, Gibbons SM, Gurry T, Irizarry RA, Alm EJ. Meta-analysis of gut microbiome studies identifies disease-specific and shared responses. Nat Commun. 2017;8:1784.
Rothschild D, Weissbrod O, Barkan E, Kurilshikov A, Korem T, Zeevi D, Costea PI, Godneva A, Kalka IN, Bar N, et al. Environment dominates over host genetics in shaping human gut microbiota. Nature. 2018;555:210–5.
Good SP, URycki DR, Crump BC. Predicting hydrologic function with aquatic gene fragments. Water Resour Res. 2018;54:2424–35.
Ramirez KS, Knight CG, de Hollander M, Brearley FQ, Constantinides B, Cotton A, Creer S, Crowther TW, Davison J, Delgado-Baquerizo M, et al. Detecting macroecological patterns in bacterial communities across independent studies of global soils. Nat Microbiol. 2018;3:189–96.
Schuster SC. Next-generation sequencing transforms today's biology. Nat Methods. 2008;5:16–8.
Waldor MK, Tyson G, Borenstein E, Ochman H, Moeller A, Finlay BB, Kong HH, Gordon JI, Nelson KE, Dabbagh K, Smith H. Where next for microbiome research? PLoS Biol. 2015;13:e1002050.
Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, Djahanschiri B, Zeller G, Mende DR, Alberti A, et al. Ocean plankton. Structure and function of the global ocean microbiome. Science. 2015;348:1261359.
Roik A, Rothig T, Roder C, Ziegler M, Kremb SG, Voolstra CR. Year-long monitoring of physico-chemical and biological variables provide a comparative baseline of coral reef functioning in the Central Red Sea. PLoS One. 2016;11:e0163939.
Brown EB. Coral bleaching: causes and consequences. Coral Reefs. 1997;16:S129–38.
Hughes TP, Kerry JT, Alvarez-Noriega M, Alvarez-Romero JG, Anderson KD, Baird AH, Babcock RC, Beger M, Bellwood DR, Berkelmans R, et al. Global warming and recurrent mass bleaching of corals. Nature. 2017;543:373–7.
ARC Centre of Excellence: Two-thirds of great barrier reef hit by back-to-back mass coral bleaching. Media Release 2017. https://www.coralcoe.org.au/media-releases/two-thirds-of-great-barrier-reef-hit-by-back-to-back-mass-coral-bleaching.
Garren M, Son K, Raina J-B, Rusconi R, Menolascina F, Shapiro OH, Tout J, Bourne DG, Seymour JR, Stocker R. A bacterial pathogen uses dimethylsulfoniopropionate as a cue to target heat-stressed corals. ISME J. 2014;8:999–1007.
Garren M, Son K, Tout J, Seymour JR, Stocker R. Temperature-induced behavioral switches in a bacterial coral pathogen. ISME J. 2016;10:1363–72.
Tout J, Jeffries TC, Petrou K, Tyson GW, Webster NS, Garren M, Stocker R, Ralph PJ, Seymour JR. Chemotaxis by natural populations of coral reef bacteria. ISME J. 2015;9:1764–77.
Muller EM, Rogers CS, Spitzack AS, van Woesik R. Bleaching increases likelihood of disease on Acropora palmata (Lamarck) in Hawksnest Bay, St John, US Virgin Islands. Coral Reefs. 2008;27:191–5.
vButtigieg PL, Fadeev E, Bienhold C, Hehemann L, Offre P, Boetius A. Marine microbes in 4D-using time series observation to assess the dynamics of the ocean microbiome and its links to ocean health. Curr Opin Microbiol. 2018;43:169–85.
Jarman SN, Berry O, Bunce M. The value of environmental DNA biobanking for long-term biomonitoring. Nat Ecol Evol. 2018;2:1192–3.
Devlin MJ, Lourey MJ. Water quality - field and analytical procedures. In: L–tMotGB R, editor. Standard Operational Procedure, vol. 6. Townsville: Australian Institute of Marine Science; 2000.
Lane DJ. 16S/23S rRNA sequencing. In: Stackebrandt E, Goodfellow M, editors. Nucleic acid techniques in bacterial systematics. New York: Wiley; 1991. p. 115–75.
Turner S, Pryer KM, Miao VP, Palmer JD. Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis. J Eukaryot Microbiol. 1999;46:327–38.
Bissett A, Fitzgerald A, Meintjes T, Mele PM, Reith F, Dennis PG, Breed MF, Brown B, Brown MV, Brugger J, et al. Introducing BASE: the biomes of Australian soil environments soil microbial diversity database. Gigascience. 2016;5:21.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Yilmaz P, Parfrey LW, Yarza P, Gerken J, Pruesse E, Quast C, Schweer T, Peplies J, Ludwig W, Glockner FO. The SILVA and “all-species living tree project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014;42:D643–8.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.
R Development Core Team: R: a language and environment for statistical computing. R Foundation for Statistical Computing 2008.
Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. vegan: community ecology package. R package version; 2013. p. 20–9.
Clark-Carter D. z Scores. Wiley StatsRef: Statistics Reference Online: John Wiley & Sons, Ltd; 2014.
Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carre G, Marquez JRG, Gruber B, Lafourcade B, Leitao PJ, et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography. 2013;36:27–46.
Kuhn M. Building predictive models in R using the caret package. 2008;28:1–26
Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2:18–22.
We thank Michele Skuza, Neale Johnston and the AIMS water quality team for their help with analysing the water quality samples. We also thank Heidi Luter, Katarina Damjanovic and Joe Gioffre for their assistance in the field and Sara Bell for her expertise in the laboratory.
We would like to acknowledge the contribution of the Marine Microbes (MM) and Biomes of Australian Soil Environments (BASE) projects, through the Australian Microbiome Initiative in the generation of data used in this publication. The Australian Microbiome Initiative is supported by funding from Bioplatforms Australia through the Australian Government National Collaborative Research Infrastructure Strategy (NCRIS). The study was further funded by the Advance Queensland PhD Scholarship, the Great Barrier Reef Marine Park Authority Management Award and a National Environmental Science Program (NESP) grant awarded to BG.
The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Sequencing data, metadata and protocols are available at the Bioplatforms Australia data portal under the Australian Microbiome project (https://www.bioplatforms.com/australian-microbiome/). Full usage requires free registration. To search for the sequencing data, navigate to ‘Processed data’, select ‘Amplicon is 27f519r_bacteria’ and ‘Environment is Marine’. To search for the Great Barrier Reef sampling sites, add an additional contextual filter, select ‘Sampling Site’ from the dropdown menu and search for ‘Geoffrey Bay’, ‘Pionner Bay’ and ‘Channel’.
Ethics approval and consent to participate
Samples were collected under the permit G16/38348.1 issued by the Great Barrier Reef Marine Park Authority.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary figures and tables. Supplementary material contains additional information on the frequency of sampling (Table S1) and detailed statistical outputs (Table S2-S8). Furthermore, additional supplementary figures are illustrating alpha diversity measures of microbial communities associated with the distinct coral reef habitats (Figure S1), within and between time point similarities of microbial community composition (Figure S2), PCoA plots for sediment and seawater microbiomes (Figure S3), environmental variability at Geoffrey Bay (Magnetic Island) (Figure S4), collinearity of environmental metadata collected at Geoffrey Bay (Magnetic Island) (Figure S5), microbial indicator taxa, calculated with the Indicator Value analysis, for high and low temperature, Chla, POC and TSS concentrations (Figure S6), classification of seawater temperature based on Random Forest machine learning (Figure S7), Random Forest machine learning seawater temperature regression (Figure S8) and the relative fraction of stable and transient microbiomes associated with the distinct coral reef habitats (Figure S9). (DOCX 2596 kb)
About this article
Cite this article
Glasl, B., Bourne, D.G., Frade, P.R. et al. Microbial indicators of environmental perturbations in coral reef ecosystems. Microbiome 7, 94 (2019). https://doi.org/10.1186/s40168-019-0705-7