Microbiome dynamics in the tissue and mucus of acroporid corals differ in relation to host and environmental parameters

Giulia M. Marchioro; Bettina Glasl; Aschwin H. Engelen; Ester A. Serrão; David G. Bourne; Nicole S. Webster; Pedro R. Frade

doi:10.7717/peerj.9644

Microbiome dynamics in the tissue and mucus of acroporid corals differ in relation to host and environmental parameters

Giulia M. Marchioro^1,2,3, Bettina Glasl ^3,4,5, Aschwin H. Engelen², Ester A. Serrão², David G. Bourne^3,4,5, Nicole S. Webster^3,4,6, Pedro R. Frade²

1University of Algarve, Faro, Portugal

2CCMAR - Centre of Marine Sciences, University of Algarve, Faro, Portugal

3AIMS@JCU, Townsville, Queensland, Australia

4Australian Institute of Marine Science, Townsville, Queensland, Australia

5College of Science and Engineering, James Cook University, Queensland, Townsville, Australia

6Australian Centre for Ecogenomics, University of Queensland, Brisbane, Queensland, Australia

DOI: 10.7717/peerj.9644

Published: 2020-08-17
Accepted: 2020-07-10
Received: 2020-02-07

Academic Editor: Joseph Pawlik

Subject Areas: Ecology, Marine Biology, Microbiology, Molecular Biology
Keywords: Microbial ecology, Coral reefs, Coral symbionts, Coral compartments, 16S rRNA gene sequencing, Great Barrier Reef

Copyright: © 2020 Marchioro et al.
Licence: This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.

Cite this article: Marchioro GM, Glasl B, Engelen AH, Serrão EA, Bourne DG, Webster NS, Frade PR. 2020. Microbiome dynamics in the tissue and mucus of acroporid corals differ in relation to host and environmental parameters. PeerJ 8:e9644 https://doi.org/10.7717/peerj.9644

The authors have chosen to make the review history of this article public.

Abstract

Corals are associated with diverse microbial assemblages; however, the spatial-temporal dynamics of intra-species microbial interactions are poorly understood. The coral-associated microbial community varies substantially between tissue and mucus microhabitats; however, the factors controlling the occurrence, abundance, and distribution of microbial taxa over time have rarely been explored for different coral compartments simultaneously. Here, we test (1) differentiation in microbiome diversity and composition between coral compartments (surface mucus and tissue) of two Acropora hosts (A. tenuis and A. millepora) common along inshore reefs of the Great Barrier Reef, as well as (2) the potential linkage between shifts in individual coral microbiome families and underlying host and environmental parameters. Amplicon based 16S ribosomal RNA gene sequencing of 136 samples collected over 14 months, revealed significant differences in bacterial richness, diversity and community structure among mucus, tissue and the surrounding seawater. Seawater samples were dominated by members of the Synechococcaceae and Pelagibacteraceae bacterial families. The mucus microbiome of Acropora spp. was dominated by members of Flavobacteriaceae, Synechococcaceae and Rhodobacteraceae and the tissue was dominated by Endozoicimonaceae. Mucus microbiome in both Acropora species was primarily correlated with seawater parameters including levels of chlorophyll a, ammonium, particulate organic carbon and the sum of nitrate and nitrite. In contrast, the correlation of the tissue microbiome to the measured environmental (i.e., seawater parameters) and host health physiological factors differed between host species, suggesting host-specific modulation of the tissue-associated microbiome to intrinsic and extrinsic factors. Furthermore, the correlation between individual coral microbiome members and environmental factors provides novel insights into coral microbiome-by-environment dynamics and hence has potential implications for current reef restoration and management efforts (e.g. microbial monitoring and observatory programs).

Introduction

Coral microbiomes include the well-characterized endosymbiotic dinoflagellates of the family Symbiodiniaceae, and a vast diversity of bacteria and archaea (Bourne, Morrow & Webster, 2016; Frade et al., 2016a; Rohwer et al., 2002). The microbiome has a fundamental role in the health and stability of the coral holobiont; it recycles nutrients, removes waste products and defends against pathogens (Lema, Willis & Bourne, 2012; Morris et al., 2011; Rädecker et al., 2015; Rosado et al., 2019). The coral microbiome is influenced by a variety of intrinsic and extrinsic factors. Coral microbiomes are host species-specific and were thought to remain relatively stable over space and time (Frias-Lopez et al., 2002; Rohwer et al., 2002). However, recent studies have proposed that spatial–temporal factors such as environmental parameters (Chen et al., 2011), depth (Glasl et al., 2017), geography (Hong et al., 2009; Littman et al., 2009), seasonality (Ceh, Van Keulen & Bourne, 2011; Chen et al., 2011; Hong et al., 2009; Koren & Rosenberg, 2006), coastal pollution (Klaus et al., 2007), and the physiological status of the host (Grottoli et al., 2018; Littman, Willis & Bourne, 2009) can also influence the occurrence and relative abundance of microbial taxa. For instance, Li et al. (2015) reported a dynamic relationship between the community structure of coral-associated bacteria and the seasonal variation in environmental parameters such as dissolved oxygen and rainfall. Glasl et al. (2019a) showed that although host-associated microbiomes were five-times less responsive to the environment compared to the seawater microbiome, they were still affected by environmental factors (e.g., temperature, turbidity, and nutrient concentration).

The coral provides different microhabitats for its microbial associates, including the surface mucus layer, coral tissue, skeleton and gastrovascular cavity, each differing in microbial richness, diversity and community structure, often assessed through alpha- and beta-diversity metrics (Agostini et al., 2012; Engelen et al., 2018; Pollock et al., 2018; Sweet, Croquer & Bythell, 2011). Each microhabitat has a unique set of biochemical features and harbors a specific microbial community (Engelen et al., 2018; Pollock et al., 2018; Sweet, Croquer & Bythell, 2011). Hence, revealing microhabitat-specific host-microbiome associations and their specific sensitivities to environmental fluctuations is crucial to our understanding of coral holobionts. For example, the coral surface mucus layer is a polysaccharide-protein-lipid complex that provides an interface between the coral epithelium and the surrounding seawater (Brown & Bythell, 2005). Here microbes take advantage of a nutrient-rich medium and particular microbiome members found in the coral mucus overlap with both the tissue and the seawater microbial communities (Bourne & Munn, 2005; Brown & Bythell, 2005; Glasl, Herndl & Frade, 2016; Sweet, Croquer & Bythell, 2011). In contrast to the extracellular polymeric nature of the surface mucus layer, the coral tissue consists of two distinct layers (epidermis and gastrodermis) and a connective-tissue layer, the mesoglea (Muller-Parker, D’Elia & Cook, 2015). The coral tissue harbors photosymbiotic dinoflagellates (family Symbiodiniaceae), that can provide up to 100% of energy required by their coral host (Muller-Parker, D’Elia & Cook, 2015). The Symbiodiniaceae community has been shown to vary in tandem with the bacterial community in early life stages of corals (Quigley et al., 2019) and this may be caused by the release of complex organic molecules such as the organosulfur compound dimethylsulfoniopropionate (DMSP; Bourne et al., 2013; Frade et al., 2016b). The coral tissue microbiome is mostly represented by bacteria belonging to the phyla Proteobacteria and Actinobacteria. For example, the gammaproteobacterial Endozoicomonas are abundant in the coral’s endodermal tissue and are often considered ‘true’ coral symbionts (Bayer et al., 2013; Glasl et al., 2019b; Neave et al., 2016; Neave et al., 2017). When compared to the surface mucus layer, the microbial community in the tissue is significantly less dense and diverse (Bourne & Munn, 2005; Koren & Rosenberg, 2006), likely attributed to the more spatially stable and host controlled environment (Bourne & Munn, 2005), although divergent evidence suggests the mucus is less diverse than the tissue (Pollock et al., 2018). Furthermore, tissue-associated bacterial communities form aggregations within the coral cell layers, also referred to as coral-associated microbial aggregates (CAMAs), and are often co-localized near algal symbiont cells, highlighting potential metabolic interactions between symbionts (Wada et al., 2019).

In this study, we test the hypotheses that different coral compartments (surface mucus layer and tissue) of Acropora spp. harbor distinct microbial communities and that different intrinsic and extrinsic factors explain microbiome dynamics within these compartments. Furthermore, we aim to identify significant correlations of individual bacterial families associated with coral tissue and mucus with host-physiological and seawater parameters.

Materials & Methods

Sample collection

Samples of Acropora millepora, Acropora tenuis and seawater were collected monthly, at Geoffrey Bay (Magnetic Island) in the Great Barrier Reef (Fig. S1), between February 2016 and March 2017, for amplicon based 16S ribosomal RNA (rRNA) gene sequencing along with environmental metadata. All samples were collected under the permit G16/38348.1 issued by the Great Barrier Reef Marine Park Authority.

Samples (n = 3 per sample type and per sampling event) for molecular analysis were collected as part of the Australian Microbiome Initiative and the sample procedure has previously been outlined by Glasl et al. (2019a). In brief, coral nubbins (approximately 5 cm tall) of both Acropora species were collected, rinsed with 0.2 µm filter-sterilized seawater and placed into cryogenic vials. Coral mucus from the same specimens was collected with sterile cotton swabs as previously described by Glasl, Herndl & Frade (2016). Seawater samples for molecular analysis were collected in sterile collapsible bags, pre-filtered through a 50 µm filter mesh to remove large particles, and subsequently filtered onto a 0.2 µm Sterivex filter (Millipore). Coral nubbins, mucus swabs and Sterivex filters were immediately snap frozen in liquid nitrogen after collection and stored at −80 °C until further processing. To acquire environmental information, water and sediment samples were collected in duplicate for each sampling event as described in Glasl et al. (2019a) and further analyzed according to the standard procedures of the Australian Institute of Marine Science (AIMS; Devlin & Lourey, 2000). The environmental information processed includes common reef water quality measures such as salinity, particulate organic carbon, total suspended solids, concentrations of chlorophyll a, ammonium, the sum of nitrite and nitrate, particulate nitrogen, nitrite, total nitrogen, non-purgeable organic carbon, non-purgeable inorganic carbon, phosphate and silica as well as total organic carbon in the sediment, total organic nitrogen in the sediment and grainsize percentage of sediments <0.63 µm, between 0.63 µm and 2 mm, and >2 mm. Seawater temperatures and daylight hours were obtained from AIMS long-term monitoring temperature records (http://eatlas.org.au).

Sample preparation and genetic assays

Frozen coral tissue was airbrushed into a ziploc bag with phosphate-buffered saline (PBS) solution added until all tissue was removed from the skeletal fragment (total PBS volume was recorded). The resulting tissue slurry was homogenized for 1 min at 12,500 rpm using a hand-held tissue homogenizer (Heidolph Silent Crusher M), pelleted (10 min at 16,000 rcf) and snap frozen in liquid nitrogen. DNA from the tissue and mucus samples was extracted using the DNeasy PowerBiofilm kit (QIAGEN). DNA extracts were sent on dry ice to the Ramaciotti Centre for Genomics (Sydney, Australia) for sequencing. The bacterial 16S rRNA gene was sequenced using the 27F (Lane, 1991) and 519R (Turner et al., 1999) primers on the Illumina MiSeq platform using a dual indexed 2 × 300 bp paired-end approach. Primer pairs were selected to warrant comparability across datasets of the Australian Microbiome Initiative (https://www.australianmicrobiome.com).

Sequence analysis

Sequencing data were analyzed as single nucleotide variants following the standardized platform of the Australian Microbiome Initiative (Brown et al., 2018). In brief, paired-end reads were merged using FLASH software (Magoc & Salzberg, 2011) and FASTA formatted sequences were extracted from FASTQ files. Sequences <400 bp in length, and / or containing one or more N’s, or homopolymer runs of >8 bp were removed with MOTHUR (v1.34.1; Schloss et al., 2009). Sequences were de-replicated and ordered by abundance using USEARCH (64 bit v10.0.240; Edgar, 2010). Sequences with less than 4 representatives and Chimeras were removed, and the quality-filtered sequences were mapped to chimera-free zero-radius operational taxonomic units (zOTUs). A table containing the samples and their read abundances was created and the zOTUs were taxonomically classified with SILVA v132 database (Yilmaz et al., 2014) using MOTHUR’s implementation of the Wang classifier (Wang et al., 2007) and a 60% Bayesian probability cut-off. This sequencing dataset has already been used in a previous contribution by the research group (Glasl et al., 2019a), but in the current study it is analyzed from a different perspective aiming at comparing temporal microbiome dynamics between two distinct coral compartments.

Chloroplasts and mitochondria derived reads were removed from the dataset and remaining data was rarefied to a sequencing depth of 3,500 reads per sample in R (R Core Team, 2015) using subset_taxa function in the phyloseq package (McMurdie & Holmes, 2013). Read counts per sample were transformed into relative abundances.

Coral holobiont photopigment quantification

Photopigment (chlorophyll a) concentrations in the tissue of corals were quantified using a spectrophotometric approach (Glasl et al., 2019b). Tissue pellets were thawed on ice to avoid sample degradation and resuspended in 1 ml of 90% ethanol. Samples were sonicated for 1 min and centrifuged for 5 min at 10,000 rcf. Subsequently, 700 µl of the supernatant was removed and transferred to a new tube. The resuspension, sonication and centrifugation were repeated on the remainder of the pellet. The supernatant was recovered again, combined with the previous extraction and mixed by inversion. Sample extract and 90% ethanol (blank read) were loaded in triplicate (200 µl each) to a 96-well plate and the absorbance was recorded at 470, 632, 649, 665, 696 and 750 nm in a Cytation 3 multi-mode microplate reader (BioTek, Winooski, USA) and analyzed using the software Gen5 (BioTek, Winooski, USA). Blank corrected absorbance measures were used to calculate chlorophyll a concentrations (Equation S1).

Coral protein quantification

Soluble protein concentrations of coral tissue samples were quantified using a colorimetric protein assay kit (Pierce BCA Protein Assay Kit; Glasl et al., 2019b). Tissue pellets were thawed on ice and resuspended in 1 ml PBS. The resuspension (25 µl) was added to 200 µl of working reagent from the kit in a 96-well plate. The plate was mixed thoroughly on a plate shaker for 30 s and then incubated at 37 °C for 30 min. The plate was cooled down at room temperature. The absorbance was measured at 563 nm in a Cytation 3 multi-mode microplate reader (BioTek, Winooski, USA) and analyzed using the software Gen5 (BioTek, Winooski, USA). Measurements of the standards and samples were blank corrected to remove background absorbance. For each plate, a protein standard curve was obtained using bovine serum albumin (BSA) solution at concentrations between 25 and 2,000 µg ml⁻¹.

Symbiodiniaceae cell counting

To determine cell numbers of Symbiodiniaceae in the coral tissue, the tissue pellet was thawed on ice, resuspended in 1 ml of 0.2 µl filtered seawater and fixed in 2% formaldehyde (final concentration) to preserve the symbiont cells. The solution was passed through a syringe needle to reduce cell agglomeration and diminish the bias from cell clumps. Samples were then mixed for 1 min and 10 µl of the homogenate was loaded onto a Neubauer haemocytometer (0.100 mm depth). Symbiodiniaceae cells were counted under 40× magnification with an Olympus CX31 light microscope. In total, six independent haemocytometer loadings (24 squares each with 0.1 µl volume) were used per sample to ensure robustness of density determinations.

Statistical analyses

Statistical analyses were performed using RStudio (v1.1.463). Analyses of microbial communities were performed on rarefied relative abundance data at zOTU level. zOTU richness and Shannon-Weaver diversity were compared across host compartments, host species and reference seawater samples using non-parametric Analysis of Variance (Kruskal-Wallis test using function kruskal.test), followed by Dunn’s test for multiple comparisons (function dunn.test). All p-values were adjusted using the Benjamini–Hochberg multiple comparison correction method to decrease the false discovery rate (Benjamini & Hochberg, 1995). A Venn diagram was constructed to describe the shared and unique zOTUs among mucus, tissue and seawater microbiomes using VennDiagram package (Chen & Boutros, 2011) and visualized using eulerr package (Larsson, 2020). Non-Metric Multidimensional Scaling (NMDS) was used to illustrate the microbial community structure among host species and host compartments based on Bray-Curtis dissimilarities (phyloseq package McMurdie & Holmes, 2013). Permutational Multivariate Analysis of Variance (PERMANOVA, 999 permutations) was used to test for differences in microbial structure between host species and host compartments using the adonis2 function of the vegan package (Oksanen et al., 2013).

Physiological variables were normalized (i.e., chlorophyll a normalized to protein content, chlorophyll a normalized to Symbiodiniaceae numbers, Symbiodiniaceae density normalized to protein content) following common procedures in coral physiology studies (Frade et al., 2008; Iglesias-Prieto & Trench, 1997). Due to fragmentation of the collected coral branches, coral surface area could not be measured. Environmental and physiological variables were standardized and checked for collinearity using the Pearson correlation coefficient. Redundant variables based on Pearson’s correlation (>0.7 or <-0.7; Dormann et al., 2013) were removed from the analysis. Non-correlated variables were then used in a Bray-Curtis distance-based Redundancy Analysis (db-RDA), which quantifies the impact of the explanatory variables on the microbiome (dis)similarities (Legendre & Anderson, 1999). zOTU relative abundance, environmental and physiological metadata were used for db-RDA using the phyloseq package (McMurdie & Holmes, 2013). The analysis tests the statistical relationship between microbial community composition and the environmental/physiological variables for each coral compartment and host species combination. A model selection tool (ordiR2step function in the vegan package, sensu Blanchet, Legendre & Borcard, 2008) was performed to select the best db-RDA model (i.e., the best explanatory variables) for variation in microbiome composition of each coral compartment (mucus and tissue) in each host species (Johnson & Omland, 2004). The significance of each explanatory variable was confirmed with an ANOVA-like permutational test (function permutest) for dbRDA. The explanatory value (in %) of significant explanatory variables (e.g., environmental and physiological parameters) on each microbiome was assessed with Variation Partitioning Analysis of the vegan package (Oksanen et al., 2013). A correlation matrix (based on the default Pearson correlation) between the relative abundance of the 20 most abundant microbial families and significant environmental variables was generated using the R package MicrobiomSeq (Ssekagiri, Sloan & Ijaz, 2017), for which p-values were adjusted using the Benjamini–Hochberg multiple comparison correction (Benjamini & Hochberg, 1995).

Results

Composition of coral tissue and mucus microbiomes

The bacterial 16S rRNA genes derived from 136 samples, including coral tissue (n = 24 for A. millepora; n = 30 for A. tenuis), coral mucus layer (n = 24 for A. millepora; n = 28 for A. tenuis) and seawater (n = 30; used as reference samples) were sequenced and 12,051 zOTUs identified as single nucleotide variants.

zOTU richness differed significantly among mucus, tissue and seawater microbiomes (Kruskal-Wallis Chi $_{(2,133)}^{2} = 57.74$ , p = 2.89 × 10⁻¹³), but not between seasons (see Table S1). Coral zOTU richness differed between species (A. millepora vs A. tenuis; Kruskal-Wallis Chi $_{(1,134)}^{2} = 12.23$ , p = 0.00047). Seawater harbored the richest microbial community (558 zOTU ± 54.6), followed by the mucus (A. millepora, 220 zOTU ± 188; A. tenuis 511 zOTU ± 234) and tissue (A. millepora, 125 zOTU ± 31.6; A. tenuis, 173 zOTU ± 146; Table S1). Alpha diversity based on Shannon Index also differed significantly among microbiomes from mucus, tissue and seawater (Kruskal-Wallis Chi $_{(2,133)}^{2} = 53.37$ , p = 2.57 × 10⁻¹²), but not between seasons (see Table S1). Coral zOTU Shannon differed between species (A. millepora vs A. tenuis; Kruskal-Wallis Chi $_{(1,134)}^{2} = 6.002$ , p = 0.01429). Alpha diversity measures of mucus samples were not significantly different (Shannon Index: A. millepora, 4.18 ± 0.83; A. tenuis, 5.15 ± 0.69) from seawater samples (Shannon Index: 4.40 ± 0.209; Table S1). In contrast, the tissue microbiome was dramatically different from the mucus and seawater microbiomes and harbored the lowest microbial diversity (Shannon Index: A. millepora, 3.35 ± 0.63; A. tenuis, 3.54 ± 0.84).

Sequences affiliated to the phyla Proteobacteria dominated the microbial community of all samples (average relative abundance ± SD; mucus: 44.1 ± 11.5%; tissue: 62.8 ± 2%; seawater: 39.6 ± 3.1%), followed in dominance by Bacteroidetes (mucus: 27.5 ± 13.0%; tissue: 9.6 ± 10.9%; seawater: 12.0 ± 11.4%) and Cyanobacteria (mucus: 14.4 ± 9.0%; tissue: 9.8 ± 11.0%; seawater: 38.5 ± 4.0%). Mucus microbiomes for both Acropora species (Fig. 1) were characterized mostly by members of the family Flavobacteriaceae (average relative abundance ± SD; for A. tenuis: 17.3 ± 9.1%; A. millepora: 17.3 ± 12.7%), Synechococcaceae (A. tenuis: 12.3 ± 7.8%; A. millepora: 13.1 ± 10.2%) and Rhodobacteraceae (A. tenuis: 5.7 ± 3.0%; A. millepora: 6.4 ± 6.4%; Fig. 1). In contrast, the Endozoicimonaceae family dominated the tissue microbiome (A. tenuis: 43.2 ± 31.7%; A. millepora: 20.5 ± 19.7%), with additional representation of Flavobacteriaceae (A. tenuis: 7.9 ± 9.6%; A. millepora: 7.2 ± 9.6%), Synechococcaceae (A. tenuis: 5.5 ± 6.8%; A. millepora: 12.3 ± 14.5%) and Rhodobacteraceae (A. tenuis: 6.5 ± 10.4%; A. millepora: 5.3 ± 8.5%; Fig. 1) families. Seawater samples were mostly characterized by members of Synechococcaceae (36.6 ± 3.9%) and Pelagibacteraceae (18.6 ± 4.9%), but also by Rhodobacteraceae (8.6 ± 4.8%) and Flavobacteriaceae (8.0 ± 2.6%; Fig. 1). Tissue and mucus microbiomes exclusively shared 1,193 zOTUs (9.9%), mucus and seawater microbiomes exclusively shared 1,458 zOTUs (12.1%), whereas the tissue and seawater microbiome shared only 66 zOTUs (0.6%; Fig. 2).

Compartment-specific microbiome composition of Acropora tenuis and Acropora millepora. — Figure 1: Compartment-specific microbiome composition of *Acropora tenuis* and *Acropora millepora*.
Microbial community composition (size indicates mean relative abundance and color represents standard error) resolved for the surface mucus layer and tissue of two *Acropora* coral species (*A. tenuis* and *A. millepora*), and surrounding seawater, based on partial 16S rRNA gene amplicon sequencing. Only the 25 most abundant families across all samples are represented.

Download full-size image

DOI: 10.7717/peerj.9644/fig-1

Figure 2: Alpha diversity measures of coral mucus, coral tissue, and seawater microbiomes.
(A) Venn diagram displaying the number of shared, unique, and ubiquitous zOTUs among mucus, tissue and seawater microbiomes. Two *Acropora* species (*A. tenuis* and *A. millepora*) are pooled for the tissue and mucus microbiomes. (B) zOTU richness and (C) Shannon diversity index of microbiomes associated with tissue and mucus of *A. millepora* and *A. tenuis*, as well as with seawater.

Download full-size image

DOI: 10.7717/peerj.9644/fig-2

Microbial community composition (beta-diversity) significantly differed among mucus, tissue and seawater (Fig. 3; PERMANOVA, pseudo-F_(2,126) = 14.53, p = 0.001), between Acropora species (PERMANOVA, pseudo-F_(1,126) = 4.42, p = 0.001), and between seasons (PERMANOVA, pseudo-F_(1,126) = 1.90, p = 0.011). Interaction between species and compartment was also significant (PERMANOVA, pseudo-F_(1,126) = 3.07, p = 0.002; other interactions were not significant; Table S2).

Figure 3: Compositional variability of microbiomes associated with distinct coral-compartments and the ambient seawater.
Two-dimensional non-metric multidimensional scaling (nMDS) ordination depicting variation in microbial community structure between coral compartments (mucus and tissue) of *Acropora tenuis and Acropora millepora*, and seawater samples. “k” is the number of dimensions.

Download full-size image

DOI: 10.7717/peerj.9644/fig-3

Explanatory variables of coral tissue and mucus microbiomes

Physiological parameters of the tissue (i.e., chlorophyll a normalized to protein content, chlorophyll a normalized to Symbiodiniaceae numbers, Symbiodiniaceae density normalized to protein content) remained stable between host species.

Out of a total of 20 environmental variables measured for seawater and sediment, 6 variables were non-mutually collinear and were thus included in the db-RDA analysis. Selected variables were salinity, concentration of particulate organic carbon (POC), total suspended solids (TSS), chlorophyll a (Chla), ammonium (NH $_{4}^{+}$ ) and the sum of nitrite and nitrate concentrations (i.e., NO $_{2}^{-}$ /NO $_{3}^{-}$ ; Table S3).

Environmental/physiological parameters investigated in this study explained a limited amount of variation in the microbial community of mucus and tissue of the two Acropora species studied (Fig. 4). For example, seawater parameters explained 14% (Chla, NH $_{4}^{+}$ and NO $_{2}^{-}$ /NO $_{3}^{-}$ ) and 10% (POC and NO $_{2}^{-}$ /NO $_{3}^{-}$ ) of the compositional variability for the mucus microbiome in A. tenuis and A. millepora, respectively (ANOVA-like permutational test for dbRDA; Table S4); NO $_{2}^{-}$ /NO $_{3}^{-}$ was the only explanatory environmental variable common to the mucus microbiome of both Acropora species (5% of compositional variability explained in each species). In comparison, for the seawater microbiome, environmental parameters (NO $_{2}^{-}$ /NO $_{3}^{-}$ , TSS, POC, Salinity and Chla) explained 32% of the compositional variability of the microbiome (Fig. S2), suggesting greater environmental sensitivity by the microbial community in the seawater compared to the coral-associated communities.

Environmental and physiological drivers of the Acropora tenuis and Acropora millepora microbiomes. — Figure 4: Environmental and physiological drivers of the *Acropora tenuis* and *Acropora millepora* microbiomes.
Distance-based redundancy analysis (db-RDA) of the relationship between environmental/physiological variables and the relative abundance of bacteria in (A) mucus and (B) tissue microbiome of *A. millepora*, and (C) mucus and (D) tissue microbiome of *A. tenuis*. Arrow length indicates the strength of the correlation between the variables and the samples (note only significant variables are shown). The selected variables explain a total of (A) 14.98%, (B) 16.44%, (C) 10.63% and (D) 10.97% of the observed variance, respectively. Environmental/physiological variables represented are the sum of nitrite and nitrate concentrations (NO2–NO3), particulate organic carbon (POC), total suspended solids (TSS), ammonium concentration (NH₄), chlorophyll a concentration (Chla) in seawater and Symbiodiniaceae density normalized to protein contents (Zoox.Prot) of coral tissue.

Download full-size image

DOI: 10.7717/peerj.9644/fig-4

In contrast, tissue microbiomes of A. millepora and A. tenuis differed substantially in their response to environmental and/or to physiological parameters. While host physiology (i.e., Symbiodiniaceae density normalized to protein contents) and environment (TSS and Chla) explained 6% and 10%, respectively, of the variation of the tissue microbiome in A. tenuis, in A. millepora, the compositional variation was solely explained (10%) by environmental parameters (NO $_{2}^{-}$ /NO $_{3}^{-}$ and TSS; Variation Partitioning Analysis and ANOVA-like permutational test for dbRDA; Table S4). TSS was the only explanatory environmental variable common to the tissue microbiomes of both Acropora species (total of 5% and 4% in A. tenuis and in A. millepora, respectively).

Correlation between bacterial families and environmental/ physiological parameters

The relative abundance of Synechococcaceae derived from tissue samples of both Acropora species and the mucus of A. tenuis was negatively correlated with TSS (p = 0.025 − 0.039; Fig. 5 and Tables S5 and S6). In contrast, Synechococcaceae was positively correlated to total NO $_{2}^{-}$ /NO $_{3}^{-}$ in both species (mucus of A. tenuis, p = 0.002, Table S5; and tissue of A. millepora, p = 0.024, Tables S5). For A. tenuis, Synechococcaceae abundance derived from the tissues correlated negatively with the only significant physiological parameter; Symbiodiniaceae density normalized to protein contents (p = 0.025). In the mucus of A. millepora, the abundance of Pirellulaceae was positively correlated with NO $_{2}^{-}$ /NO $_{3}^{-}$ (p = 0.035) and negatively correlated with TSS (p = 0.019), while OCS155 was positively correlated to NO $_{2}^{-}$ /NO $_{3}^{-}$ (p = 0.015). Proteobacteria from the mucus of A. tenuis, Pelagibacteraceae and Halomonadaceae, were both strongly negative correlated with chlorophyll a in the seawater (Pelagibacteraceae, p = 0.013; Halomonadaceae, p = 0.008). Additionally, Halomonadaceae correlated negatively with NH $_{4}^{+}$ (p = 0.005; Fig. 5 and Figs. S5 and S6).

Figure 5: Bacterial taxa significantly correlated with environmental and physiological variables.
Pearson’s coefficient based correlation matrix between the 20 most abundant bacterial families and environmental/physiological variables having a significant effect on the microbiome associated to tissue and surface mucus of (A) *Acropora millepora* (n_tissue = 24, n_mucus = 24)and (B) *Acropora tenuis* (n_tissue = 30, n_mucus = 28). Significant correlations indicated by asterisks at different levels of significance (* for p < 0.05, ** for p < 0.01, *** for p < 0.001) after correction for multiple comparisons (using Benjamini-Hochberg correction). Environmental/physiological variables represented are the sum of nitrite and nitrate concentrations (NO2–NO3), particulate organic carbon (POC), total suspended solids (TSS), ammonium concentration (NH₄), chlorophyll a concentration (Chla) in seawater and Symbiodiniaceae density normalized to protein contents(Zoox.Prot) of coral tissue.

Download full-size image

DOI: 10.7717/peerj.9644/fig-5

Tissue-associated Endozoicimonaceae showed a strong significant positive correlation with Symbiodiniaceae density normalized to protein content in A. tenuis (p = 0.0003). In contrast, in the tissue of A. millepora, Endozoicimonaceae were negatively correlated with NO $_{2}^{-}$ /NO $_{3}^{-}$ (p = 0.020), whereas the abundance of Cryomorphaceae family was negatively correlated with TSS (p = 0.020; Fig. 5, Table S5).

Discussion

Microbial communities associated with corals are continually exposed to fluctuations in the surrounding environment and the physiology of their host. Previous studies have demonstrated changes in the coral microbiome in response to thermal stress (Ainsworth & Hoegh-Guldberg, 2009; Grottoli et al., 2018; Lee et al., 2015; Thurber et al., 2009), ocean acidification (Grottoli et al., 2018; Thurber et al., 2009), organic matter enrichment (Garren & Azam, 2012), bleaching events (Bourne et al., 2008) and other environmental and physiological factors (Glasl et al., 2019a; Guppy & Bythell, 2006; Kelly et al., 2014; Li et al., 2015; Pollock et al., 2018). However, the coral microbiome is not homogenous across the animal and an improved understanding of the sensitivity of the microorganisms inhabiting each coral compartment is needed. This study highlights compositional differences in the bacterial communities associated with coral mucus and coral tissue, as well as with the surrounding seawater, findings that are largely consistent with previous studies (Apprill, Weber & Santoro, 2016; Bourne & Munn, 2005; Engelen et al., 2018; Pollock et al., 2018; Sweet, Croquer & Bythell, 2011). Furthermore, the high similarity between mucus and seawater microbiomes (see Tables S1 and S2, Figs. 2 and 3) and the high dissimilarity between tissue and seawater microbiomes suggests that the mucus microbial community is more strongly influenced by the external environment than the tissue community. Similar results have been reported for other coral species in the Caribbean (Orbicella faveolata, Diploria strigosa, Montastraea cavernosa, Porites porites and Porites astreoides), where mucus and seawater shared significantly more microbial taxa than those shared by tissue and seawater microbiomes (Apprill, Weber & Santoro, 2016). Our results also support that mucus microbiomes are richer and more diverse than tissue microbiomes, which is a pattern corroborated by many previous studies (Bourne & Munn, 2005; Koren & Rosenberg, 2006).

Despite the host species-specificity of the coral microbiomes, some bacterial taxa were ubiquitously associated with a particular coral compartment. For example, Flavobacteriaceae and Synechococcaceae dominated the mucus of both species, while Endozoicimonaceae dominated the tissue microbiome of both Acropora species. However, overall microbiome composition also showed some overlap between host compartments, consistent with previous reports of overlap between the mucus and tissue microbiomes of other coral species (Engelen et al., 2018; Sweet, Croquer & Bythell, 2011). This intersection is a natural feature of the coral holobiont as both compartments are within the same host and because the constituents of the surface mucus layer are originally produced inside the tissue (Bythell & Wild, 2011). The sharing of some microbial taxa between compartments may also arise due to methodological challenges associated with retrieving samples that are exclusively mucus or coral tissue (Sweet, Croquer & Bythell, 2011), and hence these methodological limitations can obscure differences between the mucus and seawater microbiomes (Brown & Bythell, 2005).

Explanatory factors of mucus microbiome variation

We hypothesized that the coral mucus microbiome, which is in direct contact with seawater, would be primarily correlated with seawater parameters, whereas the tissue microbiome would be most affected by the physiological state of the coral host. Mucus is highly hydrated: mucocyte cells release their secretions in a condensed form which then undergo a massive swelling upon hydration, forming a visco-elastic gel (Brown & Bythell, 2005). Surface mucus can therefore be influenced by the presence of nutrients dissolved in the surrounding seawater (Tanaka, Ogawa & Miyajima, 2010). As expected, environmental factors (i.e., seawater parameters) were influential in shaping the mucus microbiome of both species (A. millepora and A. tenuis), consistent with recent studies relating changes in the mucus microbiome with environmental perturbations (Li et al., 2015; Pollock et al., 2018). However, the extent of influence from environmental parameters (10% of variation) on the mucus microbiome was much lower than the influence of environment on the seawater microbiome (32% of variation), suggesting that other factors also play a role in modulating the mucus microbiome. For instance, the surrounding environment may interact with host physiology and together they alter the bacterial community structure of the mucus. Mucus is a nutrient-rich medium fueled by the photosynthetic activity of the Symbiodiniaceae (Brown & Bythell, 2005) and therefore it is expected that some degree of variation in its chemical composition is explained by host-Symbiodiniaceae factors. For example, A. millepora and A. tenuis at the sampling site (Geoffrey Bay at Magnetic Island) associate with distinct Symbiodiniaceae (LaJeunesse et al., 2018; Ulstrup & Van Oppen, 2003; Van Oppen et al., 2001). A. millepora colonies were associated with Durusdinium (Van Oppen et al., 2001) whereas A. tenuis harbored Cladocopium spp. (Ulstrup & Van Oppen, 2003). Links between mucus chemical composition and microbiome community structure have been proposed (Tremblay et al., 2011). Physiological factors regulating the dynamics of production and release of the surface mucus layer could also contribute to regulating mucus microbial composition (Glasl, Herndl & Frade, 2016).

Fluctuations of NH $_{4}^{+}$ , NO $_{2}^{-}$ /NO $_{3}^{-}$ , Chla and POC in the surrounding seawater significantly correlated with the mucus microbiome variation in Acropora species. Li et al. (2015) and Chen et al. (2011) suggested that rainfall had a crucial effect on bacterial community variation in the coral microbiome, being mostly associated with an increase in the relative abundance of the Bacilli group (Chen et al., 2011; Li et al., 2015). In the present study, NO $_{2}^{-}$ /NO $_{3}^{-}$ (and its collinear variables daylight, particulate nitrogen and grainsize of sediments; Table S3) had the greatest influence on microbiome structure, being a significant factor for both studied species. The link between rainfall and increasing nutrients (such as NO $_{2}^{-}$ /NO $_{3}^{-}$ ) is well established for inshore reefs (Fabricius, 2005). In the current study, higher amounts of particulate and dissolved nutrients (but a decrease in TSS), corresponded to an increase in mucus-associated Synechococcaceae, Pirellulaceae, OCS155 and Rhodobacteraceae and a decrease in Halomonadaceae. For instance, Synechococcaceae in the mucus was highly positively correlated with NO $_{2}^{-}$ /NO $_{3}^{-}$ and negatively correlated with TSS. These findings corroborate previous work in which the abundance of free-living Synechococcus in shallow coastal waters decreased significantly under lower nutrient (especially nitrate) and higher TSS concentrations (Uysal & Köksalan, 2006).

Dissolved nutrients, such as nitrogen and phosphorus, can affect coral physiology and drive changes in the associated microbial community (Shaver et al., 2017; Thompson et al., 2015). For example, organic-rich nutrients from terrestrial run-off negatively affect the health of corals and promote rapid growth of opportunistic heterotrophic bacteria (e.g., Vibrionales, Flavobacteriales and Rhodobacterales), thus affecting the overall composition of the coral microbiome (McDevitt-Irwin et al., 2017; Weber et al., 2012). In our study, the abundance of Flavobacteriaceae and Rhodobacteraceae in the mucus of A. tenuis correlated with TSS and NH $_{4}^{+}$ , respectively. The coral holobiont, including cyanobacteria related to Synechococcus spp. (Lesser et al., 2004), can also efficiently take up inorganic nitrogen, for example, as nitrogen is required by the photosynthesis production of its Symbiodiniaceae symbionts (Yellowlees, Rees & Leggat, 2008). In fact, NH $_{4}^{+}$ can be assimilated by both coral and its Symbiodiniaceae (Pernice et al., 2012), and recent work has implicated bacteria such as Vibrio and Alteromonas in the incorporation and translocation of NH $_{4}^{+}$ into coral tissues and associated Symbiodiniaceae (Ceh et al., 2013). Nitrifying members of the mucus microbiome, such as ammonium oxidizing bacteria (e.g., Pirelullaceae) and archaea, are fueled by NH $_{4}^{+}$ (Beman et al., 2007; Siboni et al., 2008; Yang et al., 2013), and NO $_{2}^{-}$ /NO $_{3}^{-}$ can be respired by nitrate reducers putatively active in coral microbiomes (Siboni et al., 2008; Yang et al., 2013). Interestingly, Pirellulaceae abundances in the mucus of A. millepora positively correlated with concentrations of environmental NO $_{2}^{-}$ /NO $_{3}^{-}$ , the products of ammonium oxidation. These nitrogen-cycling processes mediated by microbes are highly dependent on oxygen availability, but because oxygen concentration in the mucus shows strong diel fluctuations (Shashar, Cohen & Loya, 1993), it is possible that both aerobic (e.g., nitrification) and anaerobic (e.g., denitrification) processes happen within the mucus layer at different times of the day. Temporal dynamics in the coral mucus microbiome are thus likely influenced by the individual and collective metabolic capabilities of the diverse assemblage of microbes and by nutrient availability in the surrounding waters.

Explanatory factors of tissue microbiome variation

The statistical relation between the coral tissue microbiome and the environmental and physiological parameters differed between coral species. Whereas the tissue microbiome of A. tenuis corresponded to both environment and host physiology, A. millepora correlated only with environmental parameters. This difference may be associated to specific features of each species, through which A. millepora could modulate the internal environment and create more stable intra-tissue conditions than A. tenuis (e.g., via skeletal light modulation, host morphology and tissue thickness, sensu Enriquez, Mendez & Iglesias-Prieto, 2005). A non-mutually exclusive alternative explanation is the influence of the algal symbiont (Symbiodiniaceae) genotype associated to the host. Little (2004) investigated Symbiodiniaceae communities associated with A. millepora and A. tenuis on Magnetic Island demonstrating that the coral-algal endosymbiotic relationship in Acropora spp. is distinct between species, dynamic and flexible (corals associate with different Symbiodiniaceae types at different life stages, for example), and contributes significantly to physiological attributes of the coral holobiont. For example, different algal genotypes can affect the nutrient availability (e.g., carbon and nitrogen) in the coral holobiont (Pernice et al., 2015; Bayliss et al., 2019). Environmental factors such as seawater temperature can also lead to temporal changes in the symbiont community (Cooper et al., 2011; Howells et al., 2012; Rocker, Willis & Bay, 2012). As the microbiome is strongly associated to the coral holobiont, any disturbance in the host-Symbiodiniaceae relationship may have indirect effects on the microbial composition and its response to environmental and physiological factors. Other studies demonstrate the influence of Symbiodiniaceae on the host microbial community and also support the idea that these two components of the coral holobiont are finely tuned (Glasl et al., 2017; Grottoli et al., 2018; Littman, Bourne & Willis, 2010; Littman, Willis & Bourne, 2009; Quigley et al., 2019). In the present study, Endozoicimonaceae were strongly positively correlated with the Symbiodiniaceae density in the tissue of A. tenuis and negatively correlated with NO $_{2}^{-}$ /NO $_{3}^{-}$ in A. millepora (see Fig. 5). These results are to some extent at odds with experimental results showing a stable dominance of Endozoicimonaceae in tissues of Pocillopora verrucosa irrespective of excess dissolved organic nitrogen and despite a bleaching phenomenon concomitant with structural changes in its Symbiodiniaceae community (Pogoreutz et al., 2018).

Besides the diversity of Symbiodiniaceae associated to each coral species, other factors can affect the coral and its response to environmental parameters, such as photochemical efficiency (Fv/Fm) and symbiont density (Cunning & Baker, 2014; Da-Anoy, Cabaitan & Conaco, 2019). For instance, Da-Anoy, Cabaitan & Conaco (2019) demonstrated a greater reduction of Fv/Fm in A. tenuis in response to elevated temperatures compared to A. millepora and the temperature responses of the corals did not directly correlate with their associated Symbiodiniaceae. This suggests that other species-specific physiological factors could modulate the responses of the coral to the environment and, indirectly, influence the tissue-associated microbiome. One such factor is the way coral-associated microbial aggregates (CAMAs) are distributed throughout the tissue, which varies within populations and can vary among coral species (Work & Aeby, 2014; Wada et al., 2019).

Total suspended solids (TSS) was the only environmental parameter measured in the present study that significantly related to the tissue microbiome of both coral species. TSS can impact corals by limiting light availability for photosynthesis and decreasing Symbiodiniaceae densities, which can indirectly affect microbial communities (Fabricius, 2005; Pollock et al., 2014). High levels of suspended solids characterize the environment of inshore reefs such as those found around Magnetic Island. The decrease in TSS is strongly associated with an increase in the abundance of tissue-associated Synechococcaceae and Cryomorphaceae. Cryomorphaceae are typical copiotrophs in the phylum Bacteroidetes and their increase in the tissue of A. millepora could relate to declines in coral holobiont health.

Conclusions

This study highlights that microbiomes inhabiting different physical microniches within the coral holobiont differ in their linkage between host and environmental factors. Microbiomes of Acropora spp. differed significantly among host compartments (surface mucus layer and tissue) and species (A. tenuis and A. millepora). Seawater parameters had the greatest influence on the mucus microbiome in both species whereas the tissue microbiomes showed differential patterns to environmental/host-physiological parameters, suggesting host-specific modulation of the tissue microbiome. While further research is needed to unequivocally define the drivers of coral microbiome variation, by investigating temporal variation in water quality and coral health measures and correlating these with microbial community dynamics across distinct host compartments in closely related species, this study has identified several intrinsic and extrinsic factors that contribute to microbiome composition in corals.

Supplemental Information

Supplementary Material

Supplementary Tables, Figures, and Equations

DOI: 10.7717/peerj.9644/supp-1

Download

Metadata of the 16S rRNA gene sequencing data and includes information such as environmental, host-physiology parameters, sampling date and sampling location

The unique sample identifier can be matched with the zOTU abundance table.

DOI: 10.7717/peerj.9644/supp-2

Download

Correlation matrix A. millepora

Correlation matrix (based on Pearson’s correlation) between the 20 most abundant bacterial families in mucus and tissue of Acropora millepora and environmental/physiological variables.

DOI: 10.7717/peerj.9644/supp-3

Download

zOTU abundances per sample

The unique sample identifier can be matched with the metadata spreadsheet and the unique zOTU identifier can be matched with the taxonomy spreadsheet.

DOI: 10.7717/peerj.9644/supp-4

Download

Correlation matrix A. tenuis

Correlation matrix (based on Pearson’s correlation) between the 20 most abundant bacterial families in mucus and tissue of Acropora tenuis and environmental/physiological variables.

DOI: 10.7717/peerj.9644/supp-5

Download

Taxonomic affiliation identified of zOTUs

The unique zOTU identifier can be matched with the zOTU abundance table.

DOI: 10.7717/peerj.9644/supp-6

Download

[1] Agostini S, Suzuki Y, Higuchi T, Casareto BE, Yoshinaga K, Nakano Y, Fujimura H. 2012. Biological and chemical characteristics of the coral gastric cavity. Coral Reefs 31(1):147-156

[2] Ainsworth TD, Hoegh-Guldberg O. 2009. Bacterial communities closely associated with coral tissues vary under experimental and natural reef conditions and thermal stress. Aquatic Biology 4:289-296

[3] Apprill A, Weber LG, Santoro AE. 2016. Distinguishing between microbial habitats unravels ecological complexity in coral microbiomes. MSystems 1(5):e00143-16

[4] Bayer T, Neave MJ, Alsheikh-Hussain A, Aranda M, Yum LK, Mincer T, Hughen K, Apprill A, Voolstra CR. 2013. The microbiome of the Red Sea Coral stylophora pistillata is dominated by tissue-associated endozoicomonas bacteria. Applied and Environmental Microbiology 79(15):4759-4762

[5] Bayliss SLJ, Scott ZR, Coffroth MA, TerHorst CP. 2019. Genetic variation in Breviolum antillogorgium, a coral reef symbiont, in response to temperature and nutrients. Ecology and Evolution 9:2803-2813

[6] Beman JM, Roberts KJ, Wegley L, Rohwer F, Francis CA. 2007. Distribution and diversity of archaeal ammonia monooxygenase genes associated with corals. Applied and Environmental Microbiology 73(17):5642-5647

[7] Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: Series B (Methodological) 57(1):289-300

[8] Blanchet FG, Legendre P, Borcard D. 2008. Forward selection of explanatory variables. Ecology 89(9):2623-2632

[9] Bourne DG, Dennis PG, Uthicke S, Soo RM, Tyson GW, Webster N. 2013. Coral reef invertebrate microbiomes correlate with the presence of photosymbionts. The ISME Journal 7(7):1452-1458

[10] Bourne D, Iida Y, Uthicke S, Smith-Keune C. 2008. Changes in coral-associated microbial communities during a bleaching event. The ISME Journal 2(4):350-363

[11] Bourne DG, Morrow KM, Webster NS. 2016. Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annual Review of Microbiology 70(1):317-340

[12] Bourne DG, Munn CB. 2005. Diversity of bacteria associated with the coral Pocillopora damicornis from the Great Barrier Reef. Environmental Microbiology 7(8):1162-1174

[13] Brown B, Bythell J. 2005. Perspectives on mucus secretion in reef corals. Marine Ecology Progress Series 296:291-309

[14] Brown MV, Van de Kamp J, Ostrowski M, Seymour JR, Ingleton T, Messer LF, Jeffries T, Siboni N, Laverock B, Bibiloni-Isaksson J, Nelson TM, Coman F, Davies CH, Frampton D, Rayner M, Goossen K, Robert S, Holmes B, Abell GCJ, Craw P, Kahlke T, Sow SLS, McAllister K, Windsor J, Skuza M, Crossing R, Patten N, Malthouse P, Van Ruth PD, Paulsen I, Fuhrman JA, Richardson A, Koval J, Bissett A, Fitzgerald A, Moltmann T, Bodrossy L. 2018. Systematic, continental scale temporal monitoring of marine pelagic microbiota by the Australian Marine Microbial Biodiversity Initiative. Scientific Data 5(1):180130

[15] Bythell JC, Wild C. 2011. Biology and ecology of coral mucus release. Journal of Experimental Marine Biology and Ecology 408(1-2):88-93

[16] Ceh J, Kilburn MR, Cliff JB, Raina J-B, Van Keulen M, Bourne DG. 2013. Nutrient cycling in early coral life stages: Pocillopora damicornis larvae provide their algal symbiont (Symbiodinium) with nitrogen acquired from bacterial associates. Ecology and Evolution 3(8):2393-2400

[17] Ceh J, Van Keulen M, Bourne DG. 2011. Coral-associated bacterial communities on Ningaloo Reef, Western Australia: Coral bacterial communities, Ningaloo Reef. FEMS Microbiology Ecology 75(1):134-144

[18] Chen C-P, Tseng C-H, Chen CA, Tang S-L. 2011. The dynamics of microbial partnerships in the coral Isopora palifera. The ISME Journal 5:728-740

[19] Chen H, Boutros PC. 2011. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics 12(1):35

[20] Cooper TF, Berkelmans R, Ulstrup KE, Weeks S, Radford B, Jones AM, Doyle J, Canto M, O’Leary RA, Van Oppen MJH. 2011. Environmental factors controlling the distribution of symbiodinium harboured by the coral acropora millepora on the great barrier reef. PLOS ONE 6(10):e25536

[21] Cunning R, Baker AC. 2014. Not just who, but how many: the importance of partner abundance in reef coral symbioses. Frontiers in Microbiology 5:e400

[22] Da-Anoy JP, Cabaitan PC, Conaco C. 2019. Species variability in the response to elevated temperature of select corals in north-western Philippines. Journal of the Marine Biological Association of the United Kingdom 99(6):1273-1279

[23] Devlin MJ, Lourey MJ. 2000. Water quality—field and analytical procedures. In: L-tMotGB R, ed. Standard operational procedure. Townsville: Australian Institute of Marine Science. Vol. 6

[24] Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, Marquéz JRG, Gruber B, Lafourcade B, Leitão PJ, Münkemüller T, McClean C, Osborne PE, Reineking B, Schröder B, Skidmore AK, Zurell D, Lautenbach S. 2013. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography 36(1):27-46

[25] Edgar RC. 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26(19):2460-2461

[26] Engelen AH, Aires T, Vermeij MJA, Herndl GJ, Serrão EA, Frade PR. 2018. Host differentiation and compartmentalization of microbial communities in the azooxanthellate cupcorals tubastrea coccinea and rhizopsammia goesi in the caribbean. Frontiers in Marine Science 5:e391

[27] Enriquez S, Mendez ER, Iglesias-Prieto RI. 2005. Multiple scattering on coral skeletons enhances light absorption by symbiotic algae. Limnology and Oceanography 50(4):1025-1032

[28] Fabricius KE. 2005. Effects of terrestrial runoff on the ecology of corals and coral reefs: review and synthesis. Marine Pollution Bulletin 50(2):125-146

[29] Frade PR, Bongaerts P, Winkelhagen AJS, Tonk L, Bak RPM. 2008. In situ photobiology of corals over large depth ranges: a multivariate analysis on the roles of environment, host, and algal symbiont. Limnology and Oceanography 53(6):2711-2723

[30] Frade PR, Roll K, Bergauer K, Herndl GJ. 2016a. Archaeal and bacterial communities associated with the surface mucus of Caribbean corals differ in their degree of host specificity and community turnover over reefs. PLOS ONE 11(1):e0144702

[31] Frade PR, Schwaninger V, Glasl B, Sintes E, Hill RW, Simó R, Herndl GJ. 2016b. Dimethylsulfoniopropionate in corals and its interrelations with bacterial assemblages in coral surface mucus. Environmental Chemistry 13(2):252-265

[32] Frias-Lopez J, Zerkle AL, Bonheyo GT, Fouke BW. 2002. Partitioning of bacterial communities between seawater and healthy, black band diseased, and dead coral surfaces. Applied and Environmental Microbiology 68(5):2214-2228

[33] Garren M, Azam F. 2012. Corals shed bacteria as a potential mechanism of resilience to organic matter enrichment. The ISME Journal 6(6):1159-1165

[34] Glasl B, Bongaerts P, Elisabeth NH, Hoegh-Guldberg O, Herndl GJ, Frade PR. 2017. Microbiome variation in corals with distinct depth distribution ranges across a shallow-mesophotic gradient (15–85 m) Coral Reefs 36(2):447-452

[35] Glasl B, Bourne DG, Frade PR, Thomas T, Schaffelke B, Webster NS. 2019a. Microbial indicators of environmental perturbations in coral reef ecosystems. Microbiome 7(1):e94

[36] Glasl B, Herndl GJ, Frade PR. 2016. The microbiome of coral surface mucus has a key role in mediating holobiont health and survival upon disturbance. The ISME Journal 10(9):2280-2292

[37] Glasl B, Smith CE, Bourne DG, Webster NS. 2019b. Disentangling the effect of host-genotype and environment on the microbiome of the coral Acropora tenuis. PeerJ 7:e6377

[38] Grottoli AG, Dalcin Martins P, Wilkins MJ, Johnston MD, Warner ME, Cai W-J, Melman TF, Hoadley KD, Pettay DT, Levas S, Schoepf V. 2018. Coral physiology and microbiome dynamics under combined warming and ocean acidification. PLOS ONE 13(1):e0191156

[39] Guppy R, Bythell J. 2006. Environmental effects on bacterial diversity in the surface mucus layer of the reef coral Montastraea faveolata. Marine Ecology Progress Series 328:133-142

[40] Hong M-J, Yu Y-T, Chen CA, Chiang P-W, Tang S-L. 2009. Influence of species specificity and other factors on bacteria associated with the coral Stylophora pistillata in Taiwan. Applied and Environmental Microbiology 75(24):7797-7806

[41] Howells EJ, Beltran VH, Larsen NW, Bay LK, Willis BL, Van Oppen MJH. 2012. Coral thermal tolerance shaped by local adaptation of photosymbionts. Nature Climate Change 2(2):116-120

[42] Iglesias-Prieto R, Trench RK. 1997. Acclimation and adaptation to irradiance in symbiotic dinoflagellates. II. Response of chlorophyll-protein complexes to different photon-flux densities. Marine Biology 130(1):23-33

[43] Johnson JB, Omland KS. 2004. Model selection in ecology and evolution. Trends in Ecology & Evolution 19(2):101-108

[44] Kelly LW, Williams GJ, Barott KL, Carlson CA, Dinsdale EA, Edwards RA, Haas AF, Haynes M, Lim YW, McDole T, Nelson CE, Sala E, Sandin SA, Smith JE, Vermeij MJA, Youle M, Rohwer F. 2014. Local genomic adaptation of coral reef-associated microbiomes to gradients of natural variability and anthropogenic stressors. Proceedings of the National Academy of Sciences of the United States of America 111(28):10227-10232

[45] Klaus JS, Janse I, Heikoop JM, Sanford RA, Fouke BW. 2007. Coral microbial communities, zooxanthellae and mucus along gradients of seawater depth and coastal pollution. Environmental Microbiology 9(5):1291-1305

[46] Koren O, Rosenberg E. 2006. Bacteria Associated with Mucus and Tissues of the Coral Oculina patagonica in Summer and Winter. Applied and Environmental Microbiology 72(8):5254-5259

[47] LaJeunesse TC, Parkinson JE, Gabrielson PW, Jeong HJ, Reimer JD, Voolstra CR, Santos SR. 2018. Systematic revision of Symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Current Biology 28(16):2570-2580

[48] Lane DJ. 1991. 16S/23S rRNA sequencing. 115-175

[49] Larsson J. 2020. eulerr: area-proportional euler and venn diagrams with ellipses. R package version 6.1.0 software

[50] Lee STM, Davy SK, Tang S-L, Fan T-Y, Kench PS. 2015. Successive shifts in the microbial community of the surface mucus layer and tissues of the coral Acropora muricata under thermal stress. FEMS Microbiology Ecology 91(12):fiv142

[51] Legendre P, Anderson MJ. 1999. Distance-based redundancy analysis: testing multispecies responses in multifactorial ecological experiments. Ecological Monographs 69(1):1-24

[52] Lema KA, Willis BL, Bourne DG. 2012. Corals form characteristic associations with symbiotic nitrogen-fixing bacteria. Applied and Environmental Microbiology 78(9):3136-3144

[53] Lesser MP, Mazel CH, Gorbunov MY, Falkowski PG. 2004. Discovery of symbiotic nitrogen-fixing cyanobacteria in corals. Science 305:997-1000

[54] Li J, Chen Q, Long L-J, Dong J-D, Yang J, Zhang S. 2015. Bacterial dynamics within the mucus, tissue and skeleton of the coral Porites lutea during different seasons. Scientific Reports 4(1):7320

[55] Little AF. 2004. Flexibility in algal endosymbioses shapes growth in reef corals. Science 304(5676):1492-1494

[56] Littman RA, Bourne DG, Willis BL. 2010. Responses of coral-associated bacterial communities to heat stress differ with Symbiodinium type on the same coral host. Molecular Ecology 19(9):1978-1990

[57] Littman R, Willis B, Bourne DG. 2009. Bacterial communities of juvenile corals infected with different Symbiodinium (dinoflagellate) clades. Marine Ecology Progress Series 389:45-59

[58] Littman RA, Willis BL, Pfeffer C, Bourne DG. 2009. Diversities of coral-associated bacteria differ with location, but not species, for three acroporid corals on the Great Barrier Reef: Diversity of coral-associated bacteria. FEMS Microbiology Ecology 68(2):152-163

[59] Magoc T, Salzberg SL. 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27(21):2957-2963

[60] McDevitt-Irwin JM, Baum JK, Garren M, Vega Thurber RL. 2017. Responses of coral-associated bacterial communities to local and global stressors. Frontiers in Marine Science 4:e262

[61] McMurdie PJ, Holmes S. 2013. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLOS ONE 8(4):e61217

[62] Morris JJ, Johnson ZI, Szul MJ, Keller M, Zinser ER. 2011. Dependence of the cyanobacterium Prochlorococcus on hydrogen peroxide scavenging microbes for growth at the ocean’s surface. PLOS ONE 6(2):e16805

[63] Muller-Parker G, D’Elia CF, Cook CB. 2015. Interactions between corals and their symbiotic algae. In: Birkeland C, ed. Coral reefs in the Anthropocene. Dordrecht: Springer. 99-116

[64] Neave MJ, Apprill A, Ferrier-Pagès C, Voolstra CR. 2016. Diversity and function of prevalent symbiotic marine bacteria in the genus Endozoicomonas. Applied Microbiology and Biotechnology 100(19):8315-8324

[65] Neave MJ, Rachmawati R, Xun L, Michell CT, Bourne DG, Apprill A, Voolstra CR. 2017. Differential specificity between closely related corals and abundant Endozoicomonas endosymbionts across global scales. The ISME Journal 11(1):186-200

[66] Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. 2013. vegan: community ecology package. R packageversion 2.0-9 software

[67] Pernice M, Dunn SR, Tonk L, Dove S, Domart-Coulon I, Hoppe P, Schintlmeister A, Wagner M, Meibom A. 2015. A nanoscale secondary ion mass spectrometry study of dinoflagellate functional diversity in reef-building corals. Environmental Microbiology 17(10):3570-3580

[68] Pernice M, Meibom A, Van den Heuvel A, Kopp C, Domart-Coulon I, Hoegh-Guldberg O, Dove S. 2012. A single-cell view of ammonium assimilation in coral-dinoflagellate symbiosis. The ISME Journal 6(7):1314-1324

[69] Pogoreutz C, Radecker N, Cardenas A, Gardes A, Wild C, Voolstra CR. 2018. Dominance of Endozoicomonas bacteria throughout coral bleaching and mortality suggests structural inflexibility of the Pocillopora verrucosa microbiome. Ecology and Evolution 8:2240-2252

[70] Pollock FJ, Lamb JB, Field SN, Heron SF, Schaffelke B, Shedrawi G, Bourne DG, Willis BL. 2014. Sediment and turbidity associated with offshore dredging increase coral disease prevalence on nearby reefs. PLOS ONE 9(7):e102498

[71] Pollock FJ, McMinds R, Smith S, Bourne DG, Willis BL, Medina M, Thurber RV, Zaneveld JR. 2018. Coral-associated bacteria demonstrate phylosymbiosis and cophylogeny. Nature Communications 9(1):e4921

[72] Quigley KM, Alvarez-Roa C, Torda G, Bourne DG, Willis BL. 2019. Co-dynamics of Symbiodiniaceae and bacterial populations during the first year of symbiosis with Acropora tenuis juveniles. MicrobiologyOpen 9(2):e959

[73] R Core Team. 2015. R: a language and environment for statistical computing. Vienna: R foundation for statistical computing. software

[74] Rädecker N, Pogoreutz C, Voolstra CR, Wiedenmann J, Wild C. 2015. Nitrogen cycling in corals: the key to understanding holobiont functioning? Trends in Microbiology 23(8):490-497

[75] Rocker MM, Willis BL, Bay LK. 2012. Thermal stress-related gene expression in corals with different Symbiodinium types. In: Proceedings of the 12th international coral reef symposium. 9:1-5

[76] Rohwer F, Seguritan V, Azam F, Knowlton N. 2002. Diversity and distribution of coral-associated bacteria. Marine Ecology Progress Series 243:1-10

[77] Rosado PM, Leite DCA, Duarte GAS, Chaloub RM, Jospin G, Nunes da Rocha U, Saraiva JP, Dini-Andreote F, Eisen JA, Bourne DG, Peixoto RS. 2019. Marine probiotics: increasing coral resistance to bleaching through microbiome manipulation. The ISME Journal 13(4):921-936

[78] Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. 2009. Introducing mothur: open-Source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology 75(23):7537-7541

[79] Shashar N, Cohen Y, Loya Y. 1993. Extreme diel fluctuations of oxygen in diffusive boundary layers surrounding stony corals. The Biological Bulletin 185(3):455-461

[80] Shaver EC, Shantz AA, McMinds R, Burkepile DE, Vega Thurber RL, Silliman BR. 2017. Effects of predation and nutrient enrichment on the success and microbiome of a foundational coral. Ecology 98(3):830-839

[81] Siboni N, Ben-Dov E, Sivan A, Kushmaro A. 2008. Global distribution and diversity of coral-associated Archaea and their possible role in the coral holobiont nitrogen cycle. Environmental Microbiology 10(11):2979-2990

[82] Ssekagiri A, Sloan WT, Ijaz UZ. 2017. microbiomeSeq: an R package for analysis of microbial communities in an environmental context. software

[83] Sweet MJ, Croquer A, Bythell JC. 2011. Bacterial assemblages differ between compartments within the coral holobiont. Coral Reefs 30(1):39-52

[84] Tanaka Y, Ogawa H, Miyajima T. 2010. Effects of nutrient enrichment on the release of dissolved organic carbon and nitrogen by the scleractinian coral Montipora digitata. Coral Reefs 29(3):675-682

[85] Thompson JR, Rivera HE, Closek CJ, Medina M. 2015. Microbes in the coral holobiont: partners through evolution, development, and ecological interactions. Frontiers in Cellular and Infection Microbiology 4:e176

[86] Thurber RV, Willner-Hall D, Rodriguez-Mueller B, Desnues C, Edwards RA, Angly F, Dinsdale E, Kelly L, Rohwer F. 2009. Metagenomic analysis of stressed coral holobionts. Environmental Microbiology 11(8):2148-2163

[87] Tremblay P, Weinbauer MG, Rottier C, Guérardel Y, Nozais C, Ferrier-Pagès C. 2011. Mucus composition and bacterial communities associated with the tissue and skeleton of three scleractinian corals maintained under culture conditions. Journal of the Marine Biological Association of the United Kingdom 91(3):649-657

[88] Turner S, Pryer KM, Miao VP, Palmer JD. 1999. Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis 1. Journal of Eukaryotic Microbiology 46(4):327-338

[89] Ulstrup KE, Van Oppen MJH. 2003. Geographic and habitat partitioning of genetically distinct zooxanthellae (Symbiodinium) in Acropora corals on the Great Barrier Reef. Molecular Ecology 12(12):3477-3484

[90] Uysal Z, Köksalan İ. 2006. The annual cycle of Synechococcus (cyanobacteria) in the northern Levantine Basin shelf waters (Eastern Mediterranean) Marine Ecology 27(3):187-197

[91] Van Oppen MJ, Palstra FP, Piquet AMT, Miller DJ. 2001. Patterns of coral-dinoflagellate associations in Acropora: significance of local availability and physiology of Symbiodinium strains and host-symbiont selectivity. Proceedings of the Royal Society of London. Series B: Biological Sciences 268(1478):1759-1767

[92] Wada N, Ishimochi M, Matsui T, Pollock FJ, Tang SL, Ainsworth TD, Willis BL, Mano N, Bourne DG. 2019. Characterization of coral-associated microbial aggregates (CAMAs) within tissues of the coral Acropora hyacinthus. Scientific Reports 9(1):14662

[93] Wang Q, Garrity GM, Tiedje JM, Cole JR. 2007. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and Environmental Microbiology 73(16):5261-5267