A comprehensive analysis of the microbial communities of healthy and diseased marine macroalgae and the detection of known and potential bacterial pathogens

Microorganisms are increasingly being recognized as the causative agents in the diseases of marine higher organisms, such as corals, sponges, and macroalgae. Delisea pulchra is a common, temperate red macroalga, which suffers from a bleaching disease. Two bacterial strains, Nautella italica R11 and Phaeobacter gallaeciensis LSS9, have been shown in vitro to cause bleaching symptoms, but previous work has failed to detect them during a natural bleaching event. To provide a link between in vitro observations and natural occurrences of the disease, we employ here deep-sequencing of the 16S rRNA gene to comprehensively analyze the community composition of healthy and diseased D. pulchra samples from two separate locations. We observed operational taxonomic units (OTUs) with 100% identity and coverage to the 16S RNA gene sequence of both in vitro pathogens, but only the OTU with similarity to strain LSS9 showed a statistically significant higher abundance in diseased samples. Our analysis also reveals the existence of other bacterial groups within the families Rhodobacteraceae and Flavobacteriaceae that strongly contribute to difference between diseased and healthy samples and thus these groups potentially contain novel macroalgal pathogens and/or saprophytes. Together our results provide evidence for the ecological relevance of one kind of in vitro pathogen, but also highlight the possibility that multiple opportunistic pathogens are involved in the bleaching disease of D. pulchra.


INTRODUCTION
Marine sessile macroorganisms, such as seaweed, sponges, and corals, are often colonized by a large number and diversity of bacteria, with which they can have either positive, neutral or negative interactions (Ainsworth et al., 2010;Egan et al., 2012;Hollants et al., 2013). In recent years, disease caused by microorganisms, have been increasingly recognized as a major negative interaction that influences the composition and function of benthic community members (Bourne et al., 2009;Egan et al., 2013). Disease of marine invertebrates and macroalgae has also been clearly linked to changes in the marine environment, such as anthropogenic stressors (pollutants, urbanization etc.) and climate change (e.g., Campbell et al., 2011). However, linking bacterial pathogens to particular diseases in the marine environment is often challenging Egan et al., 2013). This can be due to many factors, including the opportunistic nature of pathogens, the existence of multiple pathogens causing the same disease and/or the inability to sensitively detect a particular pathogen in complex environmental samples.
Delisea pulchra is a marine red macroalgae commonly found across the temperate Eastern coast of Australia, but is also found as widespread as Japan and Antarctica (Papenfuss, 1964). D. pulchra suffers from a bleaching disease, which is characterized by a loss of pigments and which occurs more frequently during summer months, when ultra-violet light radiation and temperature are elevated . The bleaching disease has also been shown to have a significant impact on fecundity and survival of the red algae (Campbell et al., 2014). The involvement of bacteria in the bleaching disease is implied by a significant difference in the bacterial community between healthy and diseased D. pulchra individuals Fernandes et al., 2012). Furthermore, two pathogens, Nautella italica R11 and Phaeobacter gallaeciensis LSS9, have been isolated that can cause the bleaching disease in vitro Fernandes et al., 2011). In the laboratory, these two pathogens can invade the tissue of D. pulchra under conditions of elevated temperature and when the alga's chemical defense based on UV-sensitive molecules called furanones is reduced. These in vitro experiments link pathogen function with the relevant ecological stressors or changes. However, previous studies using denaturing gradient gel electrophoresis (DGGE) from multiple bleaching events  and 16S rRNA gene sequencing of clone libraries from a single bleaching event could not detect the two in vitro pathogens (Fernandes et al., 2012). Therefore, a link between the in vitro bleaching observations Fernandes et al., 2011) and environmental disease events (Campbell et al., 2014) remains to be determined.
The previous inability of Fernandes et al. (2012) to detect N. italica R11 and P. gallaeciensis LSS9 in vivo may be due to the limited number of samples analyzed and/ or the relatively shallow sequencing analysis of the microbial community of D. pulchra. To address this issue, we use here deep-sequencing of the 16S rRNA gene to investigate the microbial community on healthy and bleached individuals from two, natural disease events.

DEEP 16S rRNA GENE SEQUENCING OF MICROBIAL COMMUNITIES SHOWS OTUs MATCHING IN VITRO PATHOGENS IN ENVIRONMENTAL SAMPLES
In this study we analyzed replicate samples for two separate bleaching events of natural D. pulchra populations in two locations (Bare Island and Long Bay) that occurred during the austral summer of 2008 (see Table S1) . Using deep-sequencing of the 16S rRNA gene V4 region with the Illumina HiSeq 2000 platform, we obtained 100 bp reads ranging from 64,151, to 142,324 sequences per microbial community sample after stringent quality filtering (see Experimental Procedures and Table S1). These reads were clustered into operational taxonomic units (OTUs) at both 97% and 100% identity. After removal of spurious OTUs and adjusting for the variation in the 16S rRNA gene copy number (see Experimental Procedures) the sequence abundances per sample ranged from 27,815 to 68,051 for the 97% identity OTUs and from 24,995 to 64,048 for 100% identity OTUs. Rarefaction analysis indicated that the sequencing effort started to saturate the diversity of the 16S rRNA gene fragment for all samples (see Figure 1) and this was further supported by Good's coverage estimates of greater than 97% for both OTU definitions (see Table S1). At the level of phylogenetic resolution provided by the 16S rRNA gene fragment sequenced, our results show an almost complete sampling of the bacterial diversity present on the surface of D. pulchra, which was not achieved with previous Sanger-based sequencing efforts (Fernandes et al., 2012, see thick, dashed line in Figure 1).
The extent of the sampling effort should detect OTUs with similarity to the previously described in vitro pathogens N. italica R11 and P. gallaeciensis LSS9, if they are present. In order to match the in vitro pathogens to OTUs with the highest phylogenetic resolution that can be obtained with our data, only here we used the OTUs defined at a 100% identity cutoff. Using a 100% query coverage and identity cut-off, a single OTU could be found to match the 16S rRNA gene sequence for each pathogen. These two OTUs occurred with abundances between 0.004% (which is the lowest possible given the sequence depth achieved) and 0.196% (Figure 2). Statistical analysis of the relative abundances in healthy and bleached samples of the two different bleaching events showed that only the OTU that matched to P. gallaeciensis LSS9 was significantly more abundant in the bleached samples (p = 0.006) overall, and in particular in Long Bay (p = 0.012) (Figure 2, Table S2). The putative N. italica R11 OTU was only detected in one sample (a bleached individual from Long Bay).

DIFFERENCE ON BACTERIAL COMMUNITIES OF HEALTHY AND BLEACHED D. PULCHRA
Given that the two OTUs assigned to the in vitro pathogens only made up a small fraction of the overall community, we FIGURE 1 | Rarefaction analysis of D. pulchra microbial community samples with OTUs defined at a 16S rRNA V4 region sequence identity of 97%. Bleached (B) and healthy (H) samples from Bare Island (BI) and Long Bay (LB) are shown. The thick dashed line, which essentially overlaps with the y-axis, shows the sequencing depth of a previous study based on 16S rRNA gene clone libraries (Fernandes et al., 2011). The thin dashed line shows the minimum number of sequences for any sample analyzed here and to which level all samples were subsampled for subsequent community comparison. investigated, if other changes occur in the community that correlated with the bleaching of D. pulchra. Based on OTUs defined at a 97% identity cut-off, D. pulchra communities had a richness (Chao1) that ranged between 1940 and 3034 and diversity (inverse Simpson index) were between 10.5 and 99.4 (Table S1), but were not significantly different between healthy and bleached samples overall or in either location ( Table 1).
In contrast, comparisons of community composition (presence/absence data) and structure (abundance data) using a Bray-Curtis similarity revealed clear differences between healthy and bleached algae (see Figures 3, S1). Hypothesis testing based on a multivariate generalized linear model (MGLM) analysis of variance (MGLM-ANOVA) showed a significant interaction between condition and location ( magnitude and/or direction of the differences between microbial communities of healthy and bleached samples can be affected by the location of the samples. Therefore, further comparisons were made within each sampling site to appropriately assess whether there are significant differences between healthy and bleached samples. Comparisons within each site showed significant difference in the community composition and structure between healthy and bleached samples from both sampling sites (Table 1), which is consistent with observations made in two previous studies Fernandes et al., 2012).
Having established general differences in the communities, we next asked which OTUs contribute in a statistically significant way to the observed changes. Univariate (i.e., OTU-by-OTU) results of the MGLM-ANOVA test using a conservative adjustment for multiple testing, showed that 31 OTUs are statistically significant different between healthy and bleached samples across the two bleaching events (Figure 4, Table S3). These OTUs, which represent 0.8% of all OTUs, contribute to the top 5% of differences (as measured by the deviance test statistic) between healthy and bleached samples and 27 of these OTUs were more abundant in bleached samples. These abundant OTUs in bleached samples had a wide taxonomic distribution, belonging to three different phyla, six classes, and nine families (Figure 4). At the phylum level, the majority of the OTUs belong to the Proteobacteria (52%) and Bacteroidetes (44%), with one OTU assigned to the Verrucomicrobia (4%). At the lowest taxonomy level, at which still most of the OTUs could be classified (93% at family level), these OTUs are classified as Rhodobacteraceae (41%) and Flavobacteriaceae (26%). The remaining OTUs belong to the Rhodospirillaceae (4%), Bacteriovoracaceae (4%), Saprospiraceae (7%), Cytophagaceae (4%), Flammeovirgaceae (4%), and Verrumicrobiaceae (4%). The four OTUs that were more abundant in healthy samples can be assigned to the family Rhodobacteraceae, the genera Schleiferia (phylum Flavobacteria), Planctomyces or Blastopirellula (as classified with RDP (Cole et al., 2014) or SILVA (Quast et al., 2013), respectively), and the class Gammaproteobacteria.

DISCUSSION THE PRESENCE OF IN VITRO PATHOGENS ON D. PULCHRA AND THEIR CONTRIBUTION TO CHANGES IN THE MICROBIAL COMMUNITY OF DISEASED INDIVIDUALS
Detection of known pathogens in natural disease events is often limited by sensitivity i.e., the lack of appropriate sequencing depth resulting in incomplete description of the microbial diversity present in a sample. A comprehensive analysis of microbial communities can be achieved by deep-sequencing of part of the 16S rRNA gene with current sequencing technologies Werner et al., 2012), however as with any PCR-based technology bias due to primer binding and amplification efficiencies can potentially misrepresent the relative abundance of certain microbial groups (e.g., Engelbrektson et al., 2010).
After the complete sequence processing, which includes OTU clustering at 100% identity, removal of potentially spurious OTUs and 16S rRNA gene copy number correction, we obtained an average of 46,255 sequence reads per sample (or a 60 times increase to the per sample sequence average of previous study by Fernandes et al., 2011), which resulted in an Frontiers in Microbiology | Aquatic Microbiology February 2015 | Volume 6 | Article 146 | 4 estimated coverage of the community diversity of greater than 97% (with the caveats of bias, mentioned above) and a detection limited for any given OTU at a relative abundance of 0.004%. This sequence analysis allowed us to detect an OTU with similarity to N. italica R11 only in one bleached sample and not in any other of the 14 samples analyzed. Considering the high sampling coverage and the low detection limit reached, this would suggest that N. italica R11 is unlikely to play a major role in the bleaching events observed here and that it constitutes a rare member of the microbial community of the D. pulchra individuals analysed here.
An OTU matching to P. gallaeciensis LSS9 was detected in 9 out of 15 samples (six bleached samples and three healthy ones) and this OTU showed a statistically significant higher abundance in bleached samples in general and in Long Bay when assessing the difference within sites (Table S2). The limited phylogenetic resolution of the V4 region used here however did not allow us to unambiguously assign the OTU to strain LSS9 as BLAST analysis also showed 100% identity to 16S rRNA genes from other members of the Rhodobacteraceae family. Nevertheless, the statistically significant increase of the OTU in one bleaching event would indicate that P. gallaeciensis LSS9 or closely-related organisms might contribute to the disease in natural settings.
The inability to consistently detect higher relative abundance for the OTUs matching to either in vitro pathogen on bleached samples could be due to two alternative, but not mutually exclusive scenarios. Firstly, the bleaching disease of D. pulchra could progress through multiple stages, in which different sets of pathogens colonize and proliferate at any one time. In this scenario, the inability to consistently detect N. italica R11 or P. gallaeciensis LSS9 could be simply because they were not abundant in the disease stages in which the samples were taken. A well-understood example for "staged diseases" is given by the Black Band Disease of corals, which involves distinct morphological stages (e.g., cyanobacterial patches and microbial lesions) and defined microbial succession patterns (Sato et al., 2013). Further temporal studies of the bleaching disease of D. pulchra would be required to define potential disease stages, however field observations have so far failed to observe distinct morphological states prior to bleaching.
The second scenario is that there is a pool of pathogens naturally present on healthy D. pulchra living as commensals. When the natural host defenses are compromised due to detrimental environmental factors (e.g., UV stress etc.), then these pathogens proliferate in an opportunistic fashion (Egan et al., 2012;Fernandes et al., 2012). Under this scenario, each time a bleaching event occurs, a different subset of opportunists dominate based on chance, similar to what has been postulated for the colonization and proliferation of planktonic bacteria on the surface of the green alga Ulva australis (Burke et al., 2011). Support for the second scenario is given by the fact that the majority (80%) of OTUs (defined at a 97% identity cut-off) that contribute most to the overall community difference between healthy and bleached samples (top 20% of the deviance) are in fact different in each bleaching event. Moreover, the majority of OTUs that are abundant on diseased samples were also present in healthy samples (71% and 88% of the OTUs that contributed most to the difference in Bare Island and Long Bay, respectively), which would be consistent with their commensal role in healthy macroalgae.

DETECTION OF MULTIPLE POTENTIAL PATHOGENS OF D. PULCHRA
In addition of attempting to detect in vitro pathogens, we also investigated if and which other OTUs (97% identity cut-off) make a significant contribution to the differences between microbial communities of bleached and healthy D. pulchra (Figure 4). Most OTUs with significant difference in relative abundance where enriched in bleached samples and of those the majority were classified as belonging to the families Rhodobacteraceae (41%) and Flavobacteriaceae (26%).
Within the Rhodobacteraceae, all OTUs that could be classified to the genus level (7 out of 11) belong to the Roseobacter clade. The Roseobacter clade is one of nine major marine clades (Giovannoni and Rappe, 2000) and its members have been found in practically every marine niche, including in associations with marine eukaryotes, such as corals, sponges, cephalopods, scallop larvae, seagrasses as well as micro-and macroalgae (Buchan et al., 2005;Wagner-Döbler and Biebl, 2006). Roseobacter clade bacteria are in fact frequently isolated from macroalgae (Brinkhoff et al., 2008) and have the capacity to utilize algal osmolytes, like putrescine, taurine, creatine, sarcosine, and dimethylsulfoniopropionate (DMSP) (Wagner-Döbler and Biebl, 2006;Kalhoefer et al., 2011;Thole et al., 2012). Of particular interest is DMSP, as several marine micro-and macroalage contain it in high concentrations (Yoch, 2002). Therefore, degrading tissue of bleached D. pulchra could possibly be a source of DMSP (and other osmolytes) for saprophytic behavior of members of the Roseobacter clade.
Alternatively, the Roseobacter-clade OTUs enriched on bleached samples could constitute opportunistic pathogens, with similar properties to N. italica R11 and P. gallaeciensis LSS9, which also belong to this group. The Roseobacter clade contains many other known or putative pathogens, such as Roseovarius crassostreae, which causes the Roseovarius Oyster Disease (Maloy et al., 2007); the strains that are consistently detected in corals affected with white plague-like disease and black band disease (Cooney et al., 2002;Pantos et al., 2003); and diseased individuals of the sponge Rhopaloides odorabile (Webster et al., 2002). As for algal disease, the bacterium Ruegeria atlantica has been shown in vitro to synthesize compounds that lyse the dinoflagellate Alexandrium catenella (Amaro et al., 2005) and Roseobacter strains can cause a tumor-like gall disease in the red alga Prionitis lanceolata (Ashen and Goff, 1998).
Within the Flavobacteriaceae, 3 out of 7 OTUs were assigned to the genus Aquimarina, with one of them (OTU1) representing the most abundant OTU in the community of bleached D. pulchra (relative abundance of 9.3%; see in Figure 4). Members of the Aquimarina have previously been reported to be associated with diseases of marine eukaryotes. For example, Aquimarina homaria dominates (jointly with another Flavobacteriaceae species) the microbial communities associated with shell lesions in American lobster Quinn et al., 2012). In addition, Aquimarina agaralytic, which was isolated from a red macroalga, was found to posses a large number of diverse agarases (Lin

www.frontiersin.org
February 2015 | Volume 6 | Article 146 | 5 et al., 2012a,b) that may function to degrade host tissue, and Aquimarina salinaria was shown to produce metabolites with algicidal activity (Chen et al., 2011). Other OTUs that contributed to the difference in bleached D. pulchra were assigned to the order Cytophagales, the family Saprospiraceae, and the genus Saprospira. These taxa along with the Flavobacteriaceae OTUs, belong to the Cytophaga/Flavobacterium/Bacteroidetes (CFB) group, which contains most of the algicidal bacteria isolated from marine and coastal environments (Fandino et al., 2001;Furusawa et al., 2003;Mayali and Azam, 2004;Roth et al., 2008;Chen et al., 2011). The potential role of such algicidal bacteria as macroalgal pathogens is thus worth further investigation. Additionally, members of the CFB group have been show to contribute to secondary infections of red alga Chondrus crispus after initial infection by an endophytic green algae (Correa and McLachlan, 1994;Craigie and Correa, 1996;Goecke et al., 2010). In conclusions, the contribution that specific OTUs of the Rhodobacteraceae and Flavobacteriaceae/CFB make to the community differences between healthy and diseased samples, together with the fact that members of these groups have previously been implicated in marine diseases, suggests that these groups of organisms contain prime candidates for alternative pathogens of D. pulchra. While the environmental observation made here shows a correlation with disease, isolates from these groups need to be obtained in the future in order to clearly demonstrate them as causative agents in either in vitro or in vivo settings. If more strains with pathogenic properties can be demonstrated to exist, then this would further support the model that environmental diseases, such as the bleaching disease of D. pulchra, result from the action and activities of widespread opportunistic pathogens (Egan et al., 2013). How these potential pathogens respond to environmental conditions and interact with other community members is likely important for the development of disease. Complex shifts in the microbial community have also been observed in diseases affecting other marine organisms, such as corals (Rosenberg et al., 2007;Thurber et al., 2009;Ainsworth et al., 2010;Mouchka et al., 2010;Littman et al., 2011) and sponges (Webster et al., 2008;Angermeier et al., 2011;Fan et al., 2013;Olson et al., 2014), and for those diseases a model involving multiple pathogens might also be applicable.

EXPERIMENTAL PROCEDURES SAMPLING AND MICROBIAL COMMUNITY DNA ISOLATION
Replicate samples of healthy (n = 7) and bleached (n = 8) D. pulchra individuals were collected by SCUBA at depths of around 9 m at Bare Island (S33 • 59 30.80 , E151 • 13 53.60 ) and Long Bay (S 33 • 57 59.79 , E151 • 15 26.11 ), off the coast of Sydney, Australia, during the austral summer (12 February 2008) (see Table S1). Each algal sample was enclosed individually in clip-sealed plastic bags in situ and transported to the laboratory, where they were rinsed in filtered seawater three times to remove any loosely associated epibionts. Algae were then gently patted with sterile paper tissue to remove excess seawater and then freeze-dried. DNA was extracted from freeze-dried algal samples (50-100 mg) using a ZR Soil Microbe DNA extraction kit (Zymo) following the manufacturer's protocol as previously described .

16S rRNA GENE SEQUENCING AND PROCESSING
DNA samples were processed through the Earth Microbiome Project (EMP) (Gilbert et al., 2010) to generate amplicons for the variable region 4 (V4) of the 16S rRNA gene using the universal bacterial/archaeal primers 515F/806R. The amplified samples were sequenced on a HiSeq Illumina 2000 platform with 100 bps from the 515F primer. Reads were demultiplexed and quality trimmed using the QIIME software. For a more detailed description of these methods see www.earthmicrobiome.org. The sequencing data have been deposited in the Short Reads Archive under accession SRX824554 (Bare Island) and SRX824555 (Long Bay).
Sequences were further processed in Mothur (Schloss et al., 2009) using the guidelines of the MiSeq standard operational procedure (http://www.mothur.org/wiki/MiSeq_SOP) with the following modifications: sequences were trimmed with parameters qwindowaverage = 30, qwindowsize = 5, maxambig = 0, maxhomop = 8, minlength = 100. The alignment was done using the reference alignment of Silva release 102 cut to the V4 region of the 16S rRNA gene. Aligned sequences were pre-clustered using diffs = 1. Chimeras were removed using chimera.uchime with dereplicate = t and contaminants (i.e., sequences from chloroplasts, mitochondria, eukaryotes or with unknown taxonomical affiliation) were filtered after classifying the sequences with the RDP version 9 reference taxonomy. Sequences were clustered into OTUs at 97% or 100% using the cluster.split command (splitmethod = classify, taxlevel = 4) and the same reference as with the alignment.

OTU-BASED COMMUNITY ANALYSIS
Due to the large number of reads produced by the HiSeq Illumina platform, sequences that have been filtered for high quality can still produce spurious OTUs. However, an additional abundancebased OTU filtering can produce an OTU collection that better reflects the true diversity of the microbial sample (Bokulich et al., 2013). Here, OTUs were removed that had an absolute abundance across samples lower than 15 (or relative abundance of 0.0009% or lower across samples), which is the total number of samples analyzed here.
The number of reads per OTU was further corrected for the known or inferred 16S rRNA gene copy number of the taxon that the OTU was assigned to. For this, the Greengenes database (version of October 2012) (DeSantis et al., 2006) was downloaded and its taxonomy file clustered at 0.99 identity were formatted according to Mothur specifications. This reference dataset was used in Mothur to classify sequences and obtain a majority consensus taxonomy using the classify.otu command with default settings. This taxonomic classification was used to create a QIIME-formatted OTU table, which was then employed in Copyrighter (Angly et al., 2014) using the default trait estimates file (ssu_img40_gg201210.txt) and option -t. The OTU table with corrected absolute abundances was converted back into the Mothur format and used for all subsequent analysis. After this adjustment, the minimum relative abundance of any OTU in any sample was 0.004%.
The number of OTUs, coverage, Chao1 and the Inverse Simpson diversity index were calculated with Mothur for 1000

Frontiers in Microbiology | Aquatic Microbiology
February 2015 | Volume 6 | Article 146 | 6 random subsamples using the smallest 16S rRNA gene abundance in any sample (i.e., sample BI_H1 with an abundance after 16S rRNA gene copy number correction of 27,815). The average of the subsamples were used to perform Analysis of Variance (ANOVA) tests between samples types. Bray-Curtis similarities were calculated for the communities' composition (presence/absence data) and structure (abundance data) and hierarchical agglomerative clustering dendograms of samples were generated using Primer-E v6 (Clarke and Gorley, 2006). To calculate the Bray-Curtis similarity of community structure, the OTU abundances were square-root transformed. To test whether there were statistically significant differences between healthy and bleached samples in both composition and structure a two-factor design (Location with levels "Bare Island" and "Long Bay" and Condition with levels "healthy" and "bleached") was used to adjust the data to a MGLM using the mvabund package (Wang et al., 2012). In this approach, each OTU is treated as a variable that is fitted to a separate generalized linear model (GLM) using a negative binomial distribution for the analysis of community structure and a binomial distribution for composition analysis. For multivariate hypothesis testing, the ANOVA function (which implements an analysis of deviance) was applied to the MGLM using the p.uni argument set to return univariate OTU-by-OTU results adjusted to control the family-wise error rate across OTUs. Additionally, these univariate ANOVA-like tests were ordered by deviance to identify the OTUs that contribute more strongly to the overall difference between healthy and bleached samples. Using the databases Silva (release 119), RDP (PDS version 10) and Greengenes (release of August 2013) in Mothur, three separate majority-based consensus taxonomic classifications (classify.otu command with default parameters) were obtained for the OTUs that had a statistically significant effect by Condition. A consensus of the three classifications was manually built by reporting only the deepest taxonomic assignment, using the highest consensus confidence results observed and showing alternative taxa when different classifications were obtained with the different reference database.

DEFINITION OF OTUs THAT MATCH THE 16S rRNA GENE OF KNOWN, IN VITRO PATHOGENS
To investigate the presence of the strains N. italica R11 and P. gallaeciensis LSS9 their 16S rRNA gene were searched against the 100% OTU clusters using blastn and a cutoff for coverage and identify of 100% for the OTU sequence. Using this criteria one OTU cluster was found for each pathogen. Because the OTU cluster assigned to N. italica R11 was present in only one sample, no further analysis were done on it. For the OTU assigned to P. gallaeciensis LSS9, the mvabund software package (Wang et al., 2012) was used to fit the by-sample standardized OTU abundance to a GLM using again a negative binomial distribution. This GLM was then used to test the significance of the difference between healthy and bleached samples using the ANOVA function (analysis of deviance).

ACKNOWLEDGMENTS
This work was support by the Earth Microbiome Project and a grant from the Australian Research Council. The author also thanks Alex Campbell for providing the algal samples used in the analysis.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb. 2015.00146/abstract