- 1Department of Biology, Hong Kong Baptist University, Hong Kong, Hong Kong
- 2Oral Biosciences, Faculty of Dentistry, The University of Hong Kong, Hong Kong, Hong Kong
- 3School of Applied Sciences, Institute for Applied Ecology New Zealand, Auckland University of Technology, Auckland, New Zealand
- 4School of Biological Sciences, The University of Hong Kong, Hong Kong, Hong Kong
- 5Department of Animal Science, University of Manitoba, Winnipeg, MB, Canada
- 6Department of Biology and Chemistry, City University of Hong Kong, Hong Kong, Hong Kong
- 7State Key Laboratory in Marine Pollution, Hong Kong, Hong Kong
Corals are rapidly declining globally due to coral diseases. Skeletal growth anomalies (SGA) or “coral tumors” are a group of coral diseases that affect coral reefs worldwide, including Hong Kong waters in the Indo-Pacific region. To better understand how bacterial communities may vary in corals with SGA, for the first time, we examined the bacterial composition associated with the apparently healthy and the diseased tissues of SGA-affected Platgyra carnosus using 16S ribosomal rRNA gene pyrosequencing. Taxonomic analysis revealed Proteobacteria, Bacteroidetes, Cyanobacteria, and Actinobacteria as the main phyla in both the apparently healthy and the diseased tissues. A significant difference in the bacterial community composition was observed between the two conditions at the OTU level. Diseased tissues were associated with higher abundances of Acidobacteria and Gemmatimonadetes, and a lower abundance of Spirochaetes. Several OTUs belonging to Rhodobacteraceae, Rhizobiales, Gammaproteobacteria, and Cytophaga-Flavobacterium-Bacteroidetes (CFB) were strongly associated with the diseased tissues. These groups of bacteria may contain potential pathogens involved with the development of SGA or opportunistic secondary or tertiary colonizers that proliferated upon the health-compromised coral host. We suggest that these bacterial groups to be further studied based on inoculation experiments and testing of Koch's postulates in efforts to understand the etiology and progression of SGA.
Introduction
In the past two decades, coral reefs have been rapidly declining due to environmental impacts, such as infectious diseases (Harvell et al., 2007), bleaching (Baker et al., 2008), fishing overexploitation (Willis et al., 2004), predator outbreaks (Kayal et al., 2012), and global climate change (Hoegh-Guldberg et al., 2007). The reduction of reef-building corals contributes to the overall degradation of marine ecosystems, including the decline in abundance and diversity of reef fish (Jones et al., 2004; Paddak et al., 2009). Coral diseases are known as one of the most significant contributors to the decline of corals worldwide (Miller and Richardson, 2011), resulting in the total loss of coral cover in some regions in the Indo-Pacific and Australia (Bruno and Selig, 2007; Sweatman et al., 2011).
Skeletal growth anomalies (SGA) or “coral tumors” are a group of coral diseases affecting the major reefs of the Indo-Pacific (Chiu et al., 2012; Tavakoli-Kolour et al., 2015), Australia (Haapkylä et al., 2010), Hawaii (Stimson, 2011), Costa Rica (Gateño et al., 2003), and the Philippines (Kaczmarsky and Richardson, 2011). SGA are abnormal tissue growths which result in raised and enlarged areas that differ markedly in morphology from the adjacent healthy tissues (Breitbart et al., 2005). Corals with SGA are associated with rapid growth (Gateño et al., 2003), a reduced number of polyps and symbiotic zooxanthellae, finer skeletal structures, and a decline in fecundity and calcification rates (Stimson, 2011). Although the etiology of SGA is unknown, studies have shown that the disease is linked to environmental factors, such as climate change (Rosenberg and Ben-Haim, 2002), sea surface temperature, and human population size (Aeby et al., 2011). For instance, the increase in sea surface temperature is often attributed for inducing physiological stress, which allows for transmission of pathogenic agents in compromised corals (Rosenberg and Ben-Haim, 2002).
Several studies have highlighted the role of bacteria in coral diseases (Sunagawa et al., 2010; Chiu et al., 2012; Godwin et al., 2012). Corals harbor a large diversity of bacteria, many of which have been shown to be beneficial by promoting coral health, defense, and nitrogen fixation (Richie, 2006; Chimetto et al., 2008). It has been suggested that bacteria contribute to the corals' ability to adapt to changes in the environment (Reshef et al., 2006), which may result in disease when under environmental stress (Rosenberg et al., 2007). Coral diseases have been shown to cause changes to the structure, diversity, and abundance of bacterial communities (Harvell et al., 2007). Meta-analysis on coral-associated bacterial communities revealed that corals with diseases generally have more Rhodobacter and Cyanobacteria sequences (Mouchka et al., 2010; Miller and Richardson, 2011). Members of the Vibrio have been extensively studied due to their identification as etiological agents in a number of coral diseases (Luna et al., 2010; Sweet et al., 2014). Other microorganisms targeted as potential pathogens include Thalassomonas loyana for white plague disease (Thompson et al., 2006), Cytophaga sp. for black band disease (Cooney et al., 2002; Miller and Richardson, 2011), and Serratia marcescens for white box disease (Harvell et al., 2007). While no clear evidence of proliferation of bacterial, viral and fungal pathogens has been observed in corals with SGA yet, a faster growth rate of bacterial cells has been observed in corals with SGA (Breitbart et al., 2005) and SGA appeared to be transmissible between colonies (Kaczmarsky and Richardson, 2007).
Very little is known about the bacterial communities of corals with SGA in the Indo-Pacific. SGA have been documented to affect 26 scleractinian coral species in the region (Sutherland et al., 2004), including Platygyra carnosus (Chiu et al., 2012), P. daedalea (Tavakoli-Kolour et al., 2015), P. pini, and P. sinensis (Sutherland et al., 2004). In Southern China, P. carnosus is a dominant scleractinian coral, forming the major structural framework of coral communities (Veron, 2000). A recent survey of Hong Kong reef-building corals indicated that 63.5% of P. carnosus colonies have developed SGA (Chiu et al., 2012), which highlights the importance of deepening our understanding of the bacterial communities associated with SGA.
We recently observed no difference between the bacterial communities associated with remote healthy and SGA-affected P. carnosus corals using the culture-dependent method (Chiu et al., 2012). Except for our previous work conducted in Hong Kong, no other previous microbiological work has pursued to compare healthy and SGA-affected coral colonies. In addition, no studies are available on the diversity of the coral bacterial community for P. carnosus. Thus, the bacterial composition involved with SGA in P. carnosus remains largely unknown. To gain better understanding of how bacterial communities associated with corals may vary within colonies affected by SGA, this paper presents the first culture-independent study to provide detailed characterization of the bacterial community composition associated with the apparently healthy and the diseased tissues of SGA-affected P. carnosus colonies. Bacterial communities were detected using 454 pyrosequencing of the 16S ribosomal rRNA genes, which have been previously employed by various studies to investigate coral-associated communities (Cárdenas et al., 2011; Cróquer et al., 2012; Godwin et al., 2012).
Materials and Methods
Site Description and Sample Collection
The sampling of Platygyra carnosus colonies for this study was authorized by the Agriculture, Fisheries and Conservation Department of Hong Kong (Permit AF GE MPA 01/5/2 pt11).
Samples of Platygyra carnosus exhibiting signs of SGA were collected on May 30, 2011 at Hoi Ha Wan Marine Park (22° 28.896′ N, 114° 19.996′ E) off the northwest side of Mo Chau (Moon Island) in Hong Kong. Samples were collected from the apparently healthy and the diseased polyps of four replicate SGA-affected P. carnosus colonies. Diseased samples were collected from tissues that displayed abnormal raised growths (Figure 1), and apparently healthy samples were collected from tissues with no visible signs of disease. Apparently healthy and diseased samples were at least 20 cm apart from each other. A modified air-driven industrial drill was used to drive coral cores of 12 mm diameter and 30–40 mm thick. All samples were washed twice by submersion in autoclaved and membrane-filtered (pore size 0.22 μm) seawater to remove loosely attached microorganisms. The core holes were subsequently filled with epoxy paste (PC11; Protective Coating Co., Allentown, PA, USA) to prevent infection. All samples were immediately transported from the collection site to the laboratory in extraction buffer (100 mM Tris–HCl, 100 mM Na2 · EDTA, 100 mM NaH2PO4, 1.5 M NaCl, and 1% CTAB) in sterile conical tubes in an icebox.
Figure 1. Picture of Platgyra carnosus skeletal growth anomaly (SGA) as indicated by an arrow and normal tissue on the remaining areas.
Extraction of Genomic DNA, PCR Amplification and Pyrosequencing
Coral samples were crushed using a mortar and pestle and vortexed for 3 min at maximum speed. Bacterial genomic DNA was extracted with phenol-chloroform-isoamyl alcohol (Liu et al., 1997) and purified with ammonium acetate precipitation (Miller et al., 1999). PCR amplification of 16S rRNA genes was performed using the FastStart High Fidelity PCR system (Roche Molecular Diagnostics, Branchburg, NJ, USA) with primers 341F (5′- ACTCCTACGGGAGGCAGCAG-3′) and 926R (5′- CCGTCAATTCMTTTGAGTTT-3′) targeting the hypervariable V3–V5 region. Both forward and reverse primers were fused to Roche Fusing primers A and B, with each Fusion primer A having a unique 10 base DNA barcode specific for each sample. PCR products were verified for quality, size, and quantity by gel electrophoresis and PicoGreen Assay. Equimolar quantities of PCR products were pooled and separated by 1% agarose gel electrophoresis. DNA corresponding to the amplified 16S rRNA gene was excised from the gels and purified using PureLink® Quick Gel Extraction Kit (Invitrogen Life Technologies, Hong Kong, China). Amplicon libraries for all samples were subjected to pyrosequencing using a bench-top 454 GS Junior (454 Life Sciences-a Roche Company, Branford, CT, USA). All sequences obtained were submitted into the National Centre for Biotechnology Information (NCBI) database (accession number SRP028284).
Bioinformatics and Statistical Analyses
Raw pyrosequence data were separated into three files (.fasta, .qual, and .txt) using sff software tools from Roche/454. Subsequent data were processed and analyzed using the Quantitative Insights Into Microbial Ecology (QIIME) pipeline (Caporaso et al., 2010b). Sequences with a mean quality score lower than 25, length of < 200 bp or >900 bp, missing primer sequences, contained ambiguous bases, contained homopolymer runs exceeding eight nucleotides or uncorrectable barcodes were removed from further analysis. Demultiplexing was performed using sample-specific barcode sequences. Denoising of sequences was performed using DENOISER v. 0.9.1 (Reeder and Knight, 2010; Quince et al., 2011), as implemented in the QIIME platform. Chimeric sequences were removed using Chimera Slayer (Hass et al., 2011). The remaining good quality sequences were rarefied and assigned into Operational Taxonomic Units (OTUs) at a threshold of 97% pair-wise nucleotide sequence identity using UCLUST (Edgar, 2010). Sequences were subjected to BLAST searches at NCBI to determine their closest relatives in the GenBank nt database (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Representative sequences at the cluster centroid for each OTU were chosen and classified to genus level by the Ribosomal Database Project (RDP) Classifier at a confidence threshold of 80% (Wang et al., 2007). The OTUs were aligned using PyNAST (Caporaso et al., 2010a) with a minimum alignment length of 150 bp and a minimum identity at 75%. PH LANE mask was performed to screen away hypervariable regions. An approximately-maximum-likelihood phylogenetic tree was constructed using FastTree (Price et al., 2010) with Kimura's 2-parameter model. Alpha diversity analysis was performed by computing rarefaction curves, diversity estimates (Shannon-Weaver index) (Magurran, 1988) and nonparametric richness estimates (Chao1) (Lee and Chao, 1994). Similarity in community composition between the apparently healthy and diseased tissues of SGA-affected colonies was compared with random even subsampling and 1000 Monte Carlo iterations using unweighted and weighted UniFrac distances visualized in Principal Coordinates Analysis (PCoA) (Lozupone et al., 2006). Similarity Percentage (SIMPER) analysis was performed to identify the taxa that contributed to most of the variance among samples (Clarke, 1993). Wilcoxon-Mann-Whitney rank sum test and Metastats test (White et al., 2009) were performed to determine the differences between the apparently healthy and diseased tissue samples for individual taxa. Analysis of similarity (ANOSIM) and Permutational Multivariate Analysis of Variance (PERMANOVA) based on the Bray-Curtis dissimilarity index (Bray and Curtis, 1957) were performed to test for the overall community differences in composition and abundance of all taxa between the apparently healthy and diseased samples with 9999 permutations. All statistical tests, unless specified, were performed in R (version 2.15.2, R Core Team, 2012).
Results
Diversity of Bacteria Associated with the Apparently Healthy and the Diseased Tissues of P. carnosus
A total of 108,015 good quality sequences, subdivided into 45,077 sequences for the apparently healthy tissues, and 62,938 sequences for the diseased tissues, were obtained from SGA-affected colonies. Diversity analysis and number of observed Operational Taxonomic Units at 97% sequence similarity were performed by averaging four random subsamples of 1802 sequences each, corresponding to the lowest number of sequences obtained from one of the coral tissue samples. Average values after even subsampling were 14,416 sequences, including 7208 sequences from both the apparently healthy and diseased tissues (Table 1). Sequences clustered into a total of 1573 OTUs, of which 550 were exclusively associated with the apparently healthy tissues, 788 were exclusively associated with the diseased tissues, and 235 were shared by the two respective groups (Table 1, Figure 2A).
Table 1. Summary of the bacterial diversity found by pyrosequencing of the apparently healthy and the diseased tissues of SGA-affected P. carnosus colonies.
Figure 2. Venn Diagrams illustrating the number of shared OTUs (97% sequence similarity) (A) between the apparently healthy and diseased P. carnosus tissues; and among the replicates of (B) the apparently healthy and (C) the diseased tissues detected in P. carnosus by pyrosequencing. Healthy, apparently healthy tissues of SGA affected P. carnosus colonies; Diseased, diseased tissues of SGA affected colonies.
Taxonomic classification of the 16S rRNA gene sequences using the RDP classifier showed that 99.95% of total sequences belonged to the domain Bacteria, 0.01% to Archaea, and 0.03% to others. Only bacterial sequences were considered in this study. Sequences were classified into 17 phyla, 31 classes, 56 orders, and 83 families for the apparently healthy tissues, and 17 phyla, 38 classes, 67 orders, and 96 families for the diseased tissues (Table 1). A high alpha-diversity (Shannon-Weaver's Index) and richness estimate (Chao1) were observed in both the diseased tissues (mean and standard error of 6.6 ± 0.4 and 579.8 ± 56.3, respectively) and the apparently healthy tissues (5.9 ± 0.4 and 349.4 ± 107, respectively) of SGA-affected colonies (Table 1). The bacterial communities associated with the apparently healthy tissues was dominated by members of the Alphaproteobacteria (25.0%), followed by Bacteroidetes (20.9%), Gammaproteobacteria (13.9%), Firmicutes (12.7%), Cyanobacteria (4.6%), and Actinobacteria (4.1%) (Figure 3, Figure S1). The phylum Proteobacteria represented 42.4% of all analyzed sequences obtained from the apparently healthy tissues and 4.3% were unclassifiable. Diseased tissues exhibited high diversity and was also dominated by members of the Alphaproteobacteria (45.1%), followed by Gammaproteobacteria (17.1%), Bacteroidetes (13.8%), Deltaproteobacteria (4.7%), Actinobacteria (3.1%), Acidobacteria (2.3%), Chloroflexi (2.3%), and Cyanobacteria (1.5%) (Figure 3, Figure S1). The phylum Proteobacteria represented 67.6% of the total number of sequences obtained from the diseased tissues. Firmicutes comprised less than 1% of all the sequences. Of the analyzed sequences from diseased samples, 4.5% were unclassifiable at the phylum level.
Figure 3. Relative abundance (%) of major bacterial phyla (except Proteobacteria into class) in sequenced OTUs from the apparently healthy and the diseased tissue samples. Healthy, apparently healthy tissues of SGA affected P. carnosus colonies; Diseased, diseased tissues of SGA affected colonies. Stacked bars were calculated from the mean of relative sequence abundance belonging to each phylum in each sample condition.
Comparison of Bacterial Communities Associated with the Apparently Healthy and the Diseased Tissues
Alpha diversity estimates and rarefaction curves were performed by subsampling to an even number of sequences. The Chao1, Shannon-Weaver, and Phylogenetic diversity indices were comparable across samples indicating that the bacterial community richness and diversity at the selected number of sequences is not different between the apparently healthy and the diseased tissues of SGA-affected colonies (Table 1).
Multivariate ANOSIM and PERMANOVA analysis based on the Bray-Curtis dissimilarity index revealed clear differences between the two conditions at the OTU (equivalent to genus-to-species) level (one-way ANOSIM, R = 0.4688, p = 0.0286; PERMANOVA, F = 1.879, p = 0.028, permutation 9999). This result is supported by UniFrac measures visualized in Principal Coordinates Analysis, which indicates that the bacterial communities generally divided into two groups (Figure 4). However, we observed variable distances in the bacterial communities between replicates of the two conditions that a natural variation between colonies existed (Figure 4). In particular, a higher degree of variation was observed among the bacterial communities of the apparently healthy tissues than the diseased tissues as indicated by the closer clustering of diseased samples in Principal Coordinates Analysis (Figure 4). These results corroborate the Venn Diagrams constructed on the apparently healthy tissues and the diseased tissues (Figures 2B,C), which indicates that there were a low number of OTUs shared by the coral colonies in each condition.
Figure 4. Principal Coordinates Analysis (PCoA) of the apparently healthy tissues and diseased tissues of tumor affected P. carnosus coral based on unweighted UniFrac matrices (weighted data not shown due to similar results). Analysis showed clustering of diseased tissues from apparently healthy tissues, suggesting that the presence or absence of certain bacterial lineages may be critical for coral's health.
Differentially Abundant Bacterial Taxa in the Apparently Healthy vs. the Diseased Coral Tissues
Wilcoxon rank-sum test and Metastats test identified a number of bacterial taxa that were differentially abundant between the apparently healthy and the diseased tissues (Table 2). SIMPER analysis was performed to determine which taxa contributed most to the differences between the two conditions. Analysis at phylum level indicated that Proteobacteria contributed most of the differences in the relative abundance between the apparently healthy and diseased tissues (37.6%), followed by Firmicutes (16.7%), Bacteroidetes (13.4%), Cyanobacteria (5.1%), and Actinobacteria (4.8%) (Table 2). Comparison of the apparently healthy and the diseased tissues showed a higher abundance of bacteria belonging to the phyla Acidobacteria and Gemmatimonadetes, and a lower abundance of bacteria belonging to Spirochaetes in the diseased samples (p < 0.05, both Wilcoxon rank-sum test and Metastats test, Table 2). SIMPER analysis at family level indicated that Rhodobacteraceae contributed to most of the differences in the relative abundance between the two conditions (14.6%), followed by Flammeovirgaceae (6.4%), Clostridiales (5.9%), Flavobacteriaceae (4.8%), Rhizobiales (4.2%), Rhodospirillaceae (3.7%), and Chromatiales (3.7%) (Table S1). Wilcoxon rank-sum test and Metastat test at family level showed that Rhizobiales, Rhodospirillaceae, and Chromatiales were differentially more abundant in the diseased tissues than the apparently healthy sections of SGA-affected P. carnosus colonies (p < 0.05, Table S1).
Identities of Bacteria Exclusively Associated with the Diseased Tissues of SGA-affected P. carnosus
A total of 788 OTUs were identified by 16S ribosomal RNA gene sequencing as exclusively associated with the diseased tissues and were absent in the apparently healthy tissues of SGA-affected colonies (Figure 2A). Of these 788 OTUs, 42 were consistently detected in all four replicates of the diseased condition (Table 3). BLAST analysis indicated that all 42 OTUs were of marine origin, and there were no evidence of bacteria of human or terrestrial origin (Table 3). Twenty-one of the 42 OTUs (50.0%) were identified as members of the class Alphaproteobacteria. Other bacterial taxa represented included Gammaproteobacteria (21.4% of OTUs), Bacteroidetes (9.5%), Actinobacteria (4.8%), Cyanobacteria (4.8%), Deltaproteobacteria (4.8%), and Acidobacteria (2.4%). The phylum Proteobacteria represented 76.2% of the sequenced OTUs found exclusively associated with the diseased tissues (Table 3). Of the analyzed OTUs, 2.4% had no close relatives in the NCBI database and could only be classified as unknown bacterial clones.
The Alphaproteobacteria exhibited high diversity and represented 21 (50%) OTUs and 75.0% of total sequences (4496 of 5994 sequences) identified as exclusively associated with diseased tissues. The class was dominated by members of the Rhodobacterales (42.9% of OTUs and 56.3% of sequences among Alphaproteobacteria), followed by Rhizobiales (33.3 and 7.9%), Rhodospirillales (9.5 and 3.8%), and Sphingomonadales (4.8 and 0.9%). Unclassified alphaproteobacterial sequences constituted 9.5% of OTUs and 6.2% of all analyzed sequences obtained from the diseased tissues. Twenty-two sequences of OTU#2586 and 448 sequences of OTU#3330 exhibited high levels of sequence similarity to an uncultured bacterium (accession no. EF123384; sequence identity 99% and KF179946; 99%, respectively) in the genus Roseovarius associated with black band disease, and Ruegeria species associated with white patch syndrome (Table 3). Members of Roseovarius (order Rhodobacterales) were the most numerically dominant genus (26.5%) among all sequences detected as exclusively associated with the diseased tissues. Also of interest are 35 sequences of OTU#3086, which were closely related to an uncultured bacterium (JF835684; sequence identity 99%) associated with bleached colonies in the coral Siderastrea stellata. All remaining OTUs belonging to Rhodobacterales were closely related to the uncultured bacterium (DQ200606; sequence identity 97%, GU066430; 100%, and AB694457; 98%, respectively) isolated from the coral Montastraea annularis, marine biofilm, and deep-sea sediment, respectively.
Wilcoxon rank-sum test and Metastats test identified a number of bacterial taxa which were significantly more abundant in the diseased tissues than in the apparently healthy tissues of SGA-affected colonies. The most noteworthy of these were OTUs belonging to the order Rhizobiales and Rhodospirillales in the class Alphaproteobacteria. Of particular interest are 48 sequences of OTU#1986 which exhibited high levels of sequence similarity to Hoeflea species associated with black band disease in the coral Siderastrea siderea. Also of interest are 119 sequences of OTU#1468, and 107 sequences of OTU#1644, which were closely related to Rhodospirillaceae species (GU118834, sequence identity 99% and JQ516560; 99%) isolated from the bacterial community of threatened Caribbean corals in the genus Pelagibius, and the coral Montastraea faveolata, respectively. Two OTUs (OTU#1880 and OTU#2325) with 99% sequence identity within Alphaproteobacteria were designated unclassified by the RDP (Table 3). BLAST analysis indicated that 357 sequences of OTU#2325 were closely related to an uncultured bacterium (FJ202551; sequence identity 99%) associated with white plague disease in the Caribbean coral Montastraea faveolata.
The class Gammaproteobacteria represented 9 (21.4%) OTUs and 18.2% of total sequences (1089 of 5647 sequences) identified as exclusively associated with diseased tissues (Table 3). Six OTUs (OTU#1891, OTU#1917, OTU#3469, OTU#689, OTU#1015, and OTU#1495) with high sequence identity (≥97%) represented an organism that was unrelated to described families within the Gammaproteobacteria and was designated unclassified by the RDP. Other OTUs identified (OTU#588, OTU#1277, and OTU#2513) belonged to the uncertain status of incertae sedis of the Gammaproteobacteria. OTU#2513 was the only OTU assigned to the genus Arenicella, while all other OTUs were taxonomically unclassified at the species level. Both OTU#1277 and OTU#1891 were closely related to species (FJ203434; sequence identity 99% and FJ202524; 99%) found associated with white plague disease in corals.
All the remaining OTUs isolated from the diseased tissues and absent from the apparently healthy tissues belonged to the phyla Bacteroidetes, Actinobacteria, Cyanobacteria, and Acidobacteria (Table 3). The most noteworthy of these was OTU#2093, which belonged to the family Flavobacteriaceae in Bacteroidetes. BLAST analysis indicated that this OTU was closely related to an uncultured bacterium (GQ412799; sequence identity 98%) from the coral-associated bacterial communities of Porites cylindrica. All of these OTUs were present only in samples of the diseased tissues and absent from samples of the apparently healthy tissues of SGA-affected colonies.
Discussion
Bacterial Diversity of the Apparently Healthy and the Diseased Tissues of P. carnosus
Coral-associated bacterial communities form complex interactions that vary with species, geographical location, seasons, and health states (Kvennefors et al., 2010; Ceh et al., 2011; Morrow et al., 2012; Roder et al., 2014; McKew et al., 2015; Zhang et al., 2015). Several investigators have shown that the diversity of the bacterial communities associated with identical coral species varied with geographic location and seasonal changes (Morrow et al., 2012; Carlos et al., 2013; McKew et al., 2015; Pantos et al., 2015; Zhang et al., 2015). For instance, a recent study has demonstrated that the bacterial communities associated with Acropora and Porites corals were more diverse in the Mexican Caribbean than their Indonesian counterparts (McKew et al., 2015), suggesting that coral-associated bacterial communities are highly dynamic and diverse. In this study, we have observed a diverse suite of bacterial species from a number of phyla, including Proteobacteria, Bacteroidetes, Cyanobacteria, and Actinobacteria. Members of these phyla have also been reported to be associated with corals of other species (Cárdenas et al., 2011; Garcia et al., 2013; Séré et al., 2013; Kellogg et al., 2014; Meyer et al., 2014; Fernando et al., 2015). It is noted that the sample size in this study was limited due to the conservation value of the coral (replicate number = 4), and this may have led to underestimating the overall diversity of the bacterial communities.
Alphaproteobacteria was identified as the predominant class, representing almost twice the number of sequences compared to the Gammaproteobacteria. These findings contradict with the results of our previous work, which have revealed dominance of gammaproteobacterial sequences among remote healthy and diseased P. caronsus colonies, in particular Vibrio species (Chiu et al., 2012). Analysis of 16S rRNA gene sequences indicated that the Gammaproteobacteria represented the second most dominant group, with Vibrionaceae representing only 0.02% of the analyzed sequences in the diseased tissues. The difference in dominance is suggested to be related to the use of culture-independent or culture-dependent approaches for identifying bacterial communities. Consistent with our results, studies based on culture-independent methods using direct amplification and sequencing in coral tissues have observed dominance of alphaproteobacterial sequences (Pantos et al., 2003; Godwin et al., 2012). On the contrary, culture-dependent studies have detected gammaproteobacterial sequences as the dominant cultivable group in reef-building corals (Pantos et al., 2003; Chiu et al., 2012; Godwin et al., 2012). The dominance of this group in culturing techniques, particularly members of the order Vibrionales, is suggested to be associated to the groups' resilience to antibiotics and their ability to produce antibiotics that inhibit the proliferation of other bacteria (Long and Azam, 2001).
Diversity and richness of the bacterial communities associated with the apparently healthy and the diseased tissues of SGA-affected colonies were comparable at the selected number of sequences, indicating similar diversity and richness of the bacterial communities. These results are inconsistent with findings which have reported higher bacterial diversity in the diseased tissues than the apparently healthy tissues of corals affected by diseases, such as white plague disease (Pantos et al., 2003) and yellow band disease (Cróquer et al., 2012; Closek et al., 2014). However, comparison of the bacterial community composition between the two conditions using Bray-Curtis similarity revealed clear differences at the OTU level.
Differentially Abundant Bacterial Taxa in Apparently Healthy vs. Diseased Tissues
Our statistical analysis indicated that the largest contributor to the difference in abundance between the apparently healthy and the diseased tissues was Proteobacteria. This result is consistent with previously reported observations of bacterial communities associated with black band disease (Sato et al., 2010; Miller and Richardson, 2011), white plague disease (Sunagawa et al., 2009; Cárdenas et al., 2011), and white patch syndrome (Séré et al., 2013). Members of the Acidobacteria, Rhizobiales, and Rhodospirillaceae were more abundant in the diseased tissues, as have been shown previously in bleached corals (Mouchka et al., 2010) and in corals with white plague disease (Cárdenas et al., 2011; Roder et al., 2014). Furthermore, we identified Rhodobacteraceae as the prominent family contributing to the most difference in abundance between the apparently healthy and the diseased tissues. Our results well relate with studies which have observed Rhodobacteraceae contributing as much as 75% of the difference between healthy and diseased corals (Roder et al., 2014).
Bacterial Groups Strongly Associated with SGA in P. carnosus
There were a number of OTUs that were strongly associated with the diseased tissues and have a high level of 16S rRNA gene sequence similarity to bacteria associated with marine diseases or known pathogens. Most OTUs with high sequence similarity to bacteria previously isolated from marine diseases were classified as belonging to the Rhodobacteraceae, Rhizobiales, Gammaproteobacteria, and Cytophaga-Flavobacterium-Bacteroidetes (CFB).
Within the Rhodobacteraceae, all OTUs that could be classified to the genus level (4 out of 9) belong to the Roseobacter clade, whose members are closely related to bacteria implicated in other marine diseases. Of particular importance is OTU#2586, which exhibited high levels of sequence similarity to an alphaproteobacterium associated with black band disease in the coral Siderastrea siderea (Sekar et al., 2008). OTU#2586 belong to the genus Roseovarius, which contains the known pathogen Roseovarius crassostreae, the causative agent associated with Juvenile Oyster Disease (JOD) in the cultured Eastern oyster Crassostrea virginica (Boettcher et al., 2005; Sekar et al., 2006, 2008). It is worthwhile noting that the relative abundance of Roseovarius makes up 25.6% of all sequences identified in the diseased tissues, but is absent in the apparently healthy tissue. Furthermore, OTU#3330 represented a bacterium assigned to the genus Ruegeria found in association with white patch syndrome in the coral Porites lutea (Séré et al., 2013). Members of the Ruegeria have previously been reported to be associated with yellow band disease in Fungiidae corals (Apprill et al., 2013). Moreover, the bacterium Ruegeria atlantica has been shown to be associated with shellfish poisoning by producing compounds that lyse the toxic dinoflagellate Alexandrium catenella (Amaro et al., 2005). Roseobacter species are known to contain phenotypes for the production of antibacterial compounds (Bruhn et al., 2007), which may have caused disruption of the native commensal coral bacterial community.
Within the Rhizobiales, OTU#1986 was assigned to the genus Hoeflea representing an uncultured bacterium found in association with corals affected by black band disease (Sekar et al., 2008). Members of the Hoeflea have previously been reported to be associated with marine diseases (Palacios et al., 2006; Fiebig et al., 2013). For example, Hoeflea alexandrii, which was isolated from the toxin-producing dinoflagellate Alexandrium minutum, was found to be associated with paralytic shellfish poisoning events (Palacios et al., 2006). In addition, Hoeflea phototrophica have previously been identified to be associated with the toxic dinoflagellate Prorocentrum lima (Fiebig et al., 2013). The Rhizobiales also contains many known pathogens, such as Aurantimonas coralicida, which causes the white plague type II disease in the coral Dichocoenia stokesi (Denner et al., 2003), and Brucella spp., which is responsible for brucellosis in a wide range of marine vertebrates (Carvalho et al., 2010).
Within the Gammaproteobacteria, OTU#1277 and OTU#1891 were found closely related to the uncultured bacterium associated with the coral Montastraea faveolata affected by white plague disease (Sunagawa et al., 2009). It is worthwhile noting that we did not detect any OTUs belonging to Vibionaceae in the bacterial community exclusively associated with the diseased tissues. The pathogenicity of Vibrio spp. in coral diseases is well documented. For instance, V. coralliilyticus and V. harveyi have been attributed for white syndrome and white band disease (Sussman et al., 2008; Luna et al., 2010; Sweet et al., 2014), and V. shiloi has been found in association with bleaching in the coral species Oculina patagonica (Kushmaro et al., 2006). Many studies have detected Vibrio spp. in both healthy and diseased corals (Breitbart et al., 2005; Chiu et al., 2012; Wilson et al., 2012). Yet the role of Vibrionaceae in coral disease remains unclear, as is the role of other species within the Gammaproteobacteria in the development of SGA.
All remaining OTUs were assigned to the Flammeovirgaceae and Flavobacteriaceae. These taxa belong to the CFB group, which have been implicated in marine diseases (Pinhassi et al., 1997; Romero et al., 2010). For example, Cytophaga species have previously been proposed to be the causative agents of black band disease (Miller and Richardson, 2011), and Flavobacteriaceae species have been found in association with marine eukaryote diseases (Chistoserdov et al., 2012; Quinn et al., 2012). The CFB group also contains a number of known bacterial pathogens, including Tenacibaculum maritimum, which is responsible for the fish disease tenacibaculosis in turbot (Romero et al., 2010; Faílde et al., 2013). In many marine ecosystems, CFB members are known to be important for the turnover of organic matter (Kirchman, 2002), which may contribute to success in colonization of CFB bacteria via the production of rich organic habitats. Additionally, the group contains most of the algicidal bacteria isolated from marine environments (Mayali and Azam, 2004; Tian et al., 2012). Members possessing algicial properties have been proposed to play a role in the causation of Yellow Blotch Disease (Cervino et al., 2004), and the rapid progression of tissue whitening in white plague disease in corals (Sunagawa et al., 2009).
These bacterial groups strongly associated with the diseased coral tissues may contain potential pathogens involved in the development of SGA in P. carnosus. Nevertheless, it is possible that the corals' defense capabilities may be compromised at the diseased status, allowing the colonization and proliferation of competitive bacteria (Pantos et al., 2003; Cárdenas et al., 2011; Godwin et al., 2012). Future work should focus on enriching isolates from these groups so that they may be used in infection experiments to attempt to determine the bacterial pathogen(s) of SGA.
Conclusion
We established that SGA in P. carnosus produces significant differences in bacterial community composition between the apparently healthy and the diseased tissues at OTU level. We also report differentially abundant members between the two conditions in Acidobacteria, Rhizobiales and Rhodospirillaceae. There were certain specific OTUs within the Rhodobacteraceae, Rhizobiales, Gammaproteobacteria, and CFB that found to be exclusively associated with the diseased tissues, which may contain potential pathogens or opportunistic colonizers that proliferate during the development of SGA. Nevertheless, our one-time sampling snapshot data makes it impossible to determine whether OTUs in these groups are primary pathogens or bacteria fueled by metabolic by-products of the disease. Since the samples were only obtained from SGA-affected colonies, we were also not able to determine whether distinct differences are present between the bacterial communities of tissues from healthy corals and apparently healthy tissues from diseased corals. The corals may have been infected and already progressed into a diseased state prior to the appearance of any visible signs. Future studies based on experimental inoculation and testing of Koch's postulates are needed in order to derive a conclusion about the role of these bacteria in the pathogenesis and the etiology of the SGA.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The work described in this paper was fully supported by the State Key Laboratory in Marine Pollution, Hong Kong, Hong Kong. We would like to thank for the critical and constructive comments from the Editor and the two Reviewers that helped improve the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2015.01142
References
Aeby, G. S., Williams, G. J., Franklin, E. C., Haapkyla, J., Harvell, C. D., Neale, S., et al. (2011). Growth anomalies on the coral genera Acropora and Porites are strongly associated with host density and human population size across the Indo-Pacific. PLoS ONE 6:e16887. doi: 10.1371/journal.pone.0016887
Amaro, A. M., Fuentes, M. S., Ogalde, S. R., Venegas, J. A., and Suárez-Isla, B. A. (2005). Identification and characterization of potentially alga-lytic marine bacteria strongly associated with the toxic dinoflagellate Alexandium catenella. J. Eukaryot. Microbiol. 52, 191–200. doi: 10.1111/j.1550-7408.2005.00031.x
Apprill, A., Hughen, K., and Mincer, T. (2013). Major similarities in the bacterial communities associated with lesioned and healthy Fungiidae corals. Environ. Microbiol. 15, 2063–2072. doi: 10.1111/1462-2920.12107
Baker, A. C., Glynn, P. W., and Riegl, B. (2008). Climate change and coral reef bleaching: an ecological assessment of long-term impacts, recovery trends and future outlook. Estuar. Coast Shelf Sci. 80, 435–371. doi: 10.1016/j.ecss.2008.09.003
Boettcher, K, J., Geaghan, K. K., Maloy, A. P., and Barber, B. J. (2005). Roseovarius crassostreae sp. nov., a member of the Roseobacter clade and the apparent cause of juvenile oyster disease (JOD) in cultured Eastern oysters. Int. J. Syst. Evol. Microl. 55, 1531–1537. doi: 10.1099/ijs.0.63620-0
Bray, R. J., and Curtis, J. T. (1957). An ordination of the upland forest communities of southern Wisconsin. Ecol. Monogr. 27, 325–249. doi: 10.2307/1942268
Breitbart, M., Bhagooli, R., Griffin, S., Johnston, I., and Rohwer, F. (2005). Microbial communities associated with skeletal tumors on Porites compressa. FEMS Microbiol. Lett. 243, 431–436. doi: 10.1016/j.femsle.2005.01.004
Bruhn, J. B., Gram, L., and Belas, R. (2007). Production of antibacterial compounds and biofilm formation by Roseobacter species are influenced by culture conditions. Appl. Environ. Microbiol. 73, 442–450. doi: 10.1128/AEM.02238-06
Bruno, J. F., and Selig, E. R. (2007). Regional decline of coral cover in the Indo-Pacific: timing, extent, and subregional comparisons. PLoS ONE 2:e711. doi: 10.1371/journal.pone.0000711
Caporaso, J. G., Bittinger, K., Bushman, F. D., DeSantis, T. Z,., Andersen, G. L., and Knight, R. (2010a). PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 2, 266–267. doi: 10.1093/bioinformatics/btp636
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010b). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Cárdenas, A., Rodriguez-R, L. M., Pizarro, V., Cadavid, L. F., and Arévalo-Ferro, C. (2011). Shifts in bacterial communities of two Caribbean reef-building coral species affected by white plague disease. ISME J. 6, 502–512. doi: 10.1038/ismej.2011.123
Carlos, C., Torres, T. T., and Ottoboni, L. M. M. (2013). Bacterial communities and species-specific associations with the mucus of Brazilian coral species. Sci. Rep. 3:1624. doi: 10.1038/srep01624
Carvalho, F. M., Souza, R. C., Barcellos, F. G., Hungria, M., and Vasconcelos, A. T. (2010). Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales. BMC Microbiol. 10:37. doi: 10.1186/1471-2180-10-37
Ceh, J., Van Keulen, M., and Bourne, D. G. (2011). Coral-associated bacterial communities on Ningaloo Reef, Western Australia. FEMS Microbiol. Ecol. 75, 134–144. doi: 10.1111/j.1574-6941.2010.00986.x
Cervino, J. M., Hayes, R. I., Polson, S. W., Polson, S. C., Goreau, T. J., Martinez, R. J., et al. (2004). Relationship of Vibrio species infection and elevated temperature to yellow blotch/band disease in Caribbean corals. Appl. Environ. Microbiol. 70, 6855–6864. doi: 10.1128/AEM.70.11.6855-6864.2004
Chimetto, L. A., Brocchi, M., Thompson, C. C., Martins, R. C. R., Ramos, H. R., and Thompson, F. L. (2008). Vibrios dominate as culturable nitrogen-fixing bacteria of the Brazilian coral Mussismilia hispida. Syst. Appl. Microbiol. 31, 312–319. doi: 10.1016/j.syapm.2008.06.001
Chistoserdov, A. Y., Quinn, R. A., Gubbala, S. L., and Smolowitz, R. (2012). Bacterial communities associated with lesions of shell disease in the American lobster, Homarus americanus Milne-Edwards. J. Shellfish Res. 31, 449–462. doi: 10.2983/035.031.0205
Chiu, J. M. Y., Li, S., Li, A., Po, B., Zhang, R., Shin, P. K. S., et al. (2012). Bacteria associated with skeletal tissue growth anomalies in the coral Platygyra carnosus. FEMS Microbiol. Ecol. 79, 380–391. doi: 10.1111/j.1574-6941.2011.01225.x
Clarke, K. R. (1993). Non-parametric multivariate analysis of changes in community structure. Aust. J. Ecol. 18, 117–143. doi: 10.1111/j.1442-9993.1993.tb00438.x
Closek, C. J., Sunagawa, S., DeSalvo, M. K., Piceno, Y. M., DeSantis, T. Z., Brodie, E. L., et al. (2014). Coral transcriptome and bacterial community profiles reveal distinct Yellow Band Disease states in Orbicella faveolata. ISME J. 8, 2411–2422. doi: 10.1038/ismej.2014.85
Cooney, R. P., Pantos, O., Le Tissier, M. D. A., Barer, M. R., O'Donnell, A. G., and Bythell, J. C. (2002). Characterization of the bacterial consortium associated with black band disease in coral using molecular microbiological techniques. Environ. Microbiol. 4, 401–413. doi: 10.1046/j.1462-2920.2002.00308.x
Cróquer, A., Bastidas, C., Elliott, A., and Sweet, M. (2012). Bacterial assemblages shifts from healthy to yellow band disease states in the dominant reef coral Montastraea faveolata. Environ. Microbiol. Rep. 5, 90–96. doi: 10.1111/j.1758-2229.2012.00397.x
Denner, E. B. M., Smith, G. W., Busse, H. J., Schumann, P., Narzt, T., Polson, S. W., et al. (2003). Aurantimonas coralicida gen. nov., sp. nov., the causative agent of white plague type II on Caribbean scleractinian corals. Int. J. Syst. Evol. Microbiol. 53, 1115–1122. doi: 10.1099/ijs.0.02359-0
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Faílde, L. D., Losada, A. P., Bermúdez, R., Santos, Y., and Quiroga, M. I. (2013). Tenacibaculum maritimum infection: pathology and immunohistochemistry in experimentally challenged turbot (Psetta maxima L.). Microb. Pathog. 65, 82–88. doi: 10.1016/j.micpath.2013.09.003
Fernando, S. C., Wang, J., Sparling, K., Garcia, G. D., Francini-Filho, R. B., de Moura, R. L., et al. (2015). Microbiota of the major South Atlantic reef building coral Mussismillia. Microb. Ecol. 69, 267–280. doi: 10.1007/s00248-014-0474-6
Fiebig, A., Pradella, S., Petersen, J., Michael, V., Päuker, O., Rohde, M., et al. (2013). Genome of the marine alphaproteobacterium Hoeflea phototrophica type strain (DFL-43(T)). Stand. Genomic Sci. 7, 440–448. doi: 10.4056/sigs.3486982
Garcia, G. D., Gregoracci, G. B., Santos Ede, O., Meirelles, P. M., Silva, G. G., Edwards, R., et al. (2013). Metagenomic analysis of healthy and white plague-affected Mussismillia braziliensis corals. Microb. Ecol. 65, 1076–1086. doi: 10.1007/s00248-012-0161-4
Gateño, D., León, A., Barki, Y., Cortés, J., and Rinkevich, B. (2003). Skeletal tumor formations in the massive coral Pavona clavus. Mar. Ecol. Prog. Ser. 258, 97–108. doi: 10.3354/meps258097
Godwin, S., Bent, E., Borneman, J., and Pereg, L. (2012). The role of coral-associated bacterial communities in Australian Subtropical White Syndrome of Turbinaria mesenterina. PLoS ONE 7:e4423. doi: 10.1371/journal.pone.0044243
Haapkylä, J., Melbourne-Thomas, J., Flavell, M., and Willis, B. L. (2010). Spatiotemporal patterns of coral disease prevalence on Heron Island, Great Barrier Reef, Australia. Coral Reefs 29, 1035–1045. doi: 10.1007/s00338-010-0660-z
Harvell, C. D., Jordan-Dahlgren, E., Merkel, S., Rosenberg, E., Raymundo, L., Smith, G., et al. (2007). Coral disease, environmental drivers, and the balance between coral and microbial associates. Oceanography 20, 172–195. doi: 10.5670/oceanog.2007.91
Hass, B. J., Gevers, D., Earl, A. M., Feldgarden, M., Ward, D. V., Giannoukos, G., et al. (2011). Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. doi: 10.1101/gr.112730.110
Hoegh-Guldberg, O., Mumby, P. J., Hooten, A. J., Steneck, R. S., Greenfield, P., Gomez, E., et al. (2007). Coral reefs under rapid climate change and ocean acidification. Science 318, 1737–1742. doi: 10.1126/science.1152509
Jones, G. P., McCornick, M. I., Srinivasan, M., and Eagle, J. V. (2004). Coral decline threatens fish biodiversity in marine reserves. Proc. Natl. Acad. Sci. U.S.A. 101, 8251–8253. doi: 10.1073/pnas.0401277101
Kaczmarsky, L., and Richardson, L. L. (2007). Transmission of growth anomalies between Indo-Pacific Porites corals. J. Invertebr. Pathol. 94, 218–221. doi: 10.1016/j.jip.2006.11.007
Kaczmarsky, L., and Richardson, L. L. (2011). Do elevated nutrients and organic carbon on Philippine reefs increase the prevalence of coral disease? Coral Reefs 30, 253–257. doi: 10.1007/s00338-010-0686-2
Kayal, M., Vercelloni, J., Lison de Loma, T., Bosserelle, P., Chancerelle, Y., Geoffroy, S., et al. (2012). Predator crown-of-thorns starfish (Acanthaster planci) outbreak, mass mortality of corals, and cascading effects on reef fish and benthic communities. PLoS ONE 7:e47363. doi: 10.1371/journal.pone.0047363
Kellogg, C. A., Piceno, Y. M., Tom, L. M., DeSantis, T. Z., Gray, M. A., and Anderson, G. L. (2014). Comparing bacterial community composition of healthy and dark spot-affected Siderastrea sidereal in Florida and the Caribbean. PLoS ONE 9:e108767. doi: 10.1371/journal.pone.0108767
Kirchman, D. L. (2002). The ecology of Cytophaga-Flavobacteria in aquatic environments. FEMS Microbiol. Ecol. 39, 91–100. doi: 10.1016/s0168-6496(01)00206-9
Kushmaro, A., Banin, E., Loya, Y., Stackebrandt, E., and Rosenberg, E. (2006). Vibrio shiloi sp. nov., the causative agent of beaching of the coral Oculina patagonica. Int. J. Syst. Evol. Microbiol. 51, 1383–1388. doi: 10.1099/00207713-51-4-1383
Kvennefors, E. C. E., Sampayo, E., Ridgway, T., Barnes, A. C., Hoegh-Guldberg, O., and Yang, C. H. (2010). Bacterial communities of two ubiquitous Great Barrier Reef corals reveals both site- and species-specificity of common bacterial associates. PLoS ONE 5:e10401. doi: 10.1371/journal.pone.0010401
Lee, S. M., and Chao, A. (1994). Estimating population-size via sample coverage for closed capture-recapture models. Biometrics 50, 88–97. doi: 10.2307/2533199
Liu, W. T., Marsh, T. L., Cheng, H., and Forney, L. J. (1997). Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl. Environ. Microbiol. 63, 4516–4522.
Long, R. A., and Azam, F. (2001). Antagonistic interactions among marine pelagic bacteria. Appl. Environ. Micrbiol. 67, 4975–4983. doi: 10.1128/AEM.67.11.4975-4983.2001
Lozupone, C., Hamady, M., and Knight, R. (2006). UniFrac - an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinformatics 7:371. doi: 10.1186/1471-2105-7-371
Luna, G. M., Bongiorni, L., Biavasco, F., and Danovaro, R. (2010). Vibrio harveyi as a causative agent of the White Syndrome in tropical stony corals. Environ. Microbiol. Rep. 2, 120–127. doi: 10.1111/j.1758-2229.2009.00114.x
Magurran, A. E. (1988). Ecological Diversity and Its Measurement. Princeton, NJ: Princeton University Press.
Mayali, X., and Azam, F. (2004). Algicidal bacteria in the sea and their impact on algal blooms. J. Eukaryot. Microbiol. 51, 139–144. doi: 10.1111/j.1550-7408.2004.tb00538.x
McKew, B. A., Dumbrell, A. J., Daud, S. D., Hepburn, E., Thorpe, L., Mogensen, L., et al. (2015). Characterization of geographically distinct bacterial communities associated with coral mucus produced by Acropora spp. and Porites spp. Appl. Environ. Microbiol. 78, 5229–5237. doi: 10.1128/AEM.07764-11
Meyer, J. L., Paul, V. J., and Teplitski, M. (2014). Community shifts in the surface microbiomes of the coral Porites astreoides with unusual lesions. PLoS ONE 9:e100316. doi: 10.1371/journal.pone.0100316
Miller, A. W., and Richardson, L. L. (2011). A meta-analysis of 16S rRNA gene clone libraries from the polymicrobial black band disease of corals. FEMS Microbiol. Ecol. 75, 231–241. doi: 10.1111/j.1574-6941.2010.00991.x
Miller, D. N., Bryant, J. E., Madsen, E. L., and Ghiorse, W. C. (1999). Evaluation and optimization of DNA extraction and purification procedures for soil and sediment samples. Appl. Environ. Microbiol. 65, 4715–4724.
Morrow, K. M., Moss, A. G., Chadwick, N. E., and Liles, M. R. (2012). Bacterial associates of two Caribbean coral species reveal species-specific distribution and geographic variability. Appl. Environ. Microb. 78, 6438–6449. doi: 10.1128/AEM.01162-12
Mouchka, M. E., Hewson, I., and Harvell, D. (2010). Coral-associated bacterial assemblages: current knowledge and the potential for climate-driven impacts. Integr. Comp. Biol. 50, 662–674. doi: 10.1093/icb/icq061
Paddak, M. J., Reynolds, J. D., Aguilar, C., Appeldoorn, R. S., Beets, J., Burkett, E. W., et al. (2009). Recent region-wide declines in Caribbean reef fish abundance. Curr. Biol. 19, 590–595. doi: 10.1016/j.cub.2009.02.041
Palacios, L., Arahai, D. R., Reguera, B., and Marín, I. (2006). Hoeflea alexandrii sp. nov., isolated from the toxic dinoflagellate Alexandrium minutum AL1V. Int. J. Syst. Evol. Microbiol. 56, 1991–1995. doi: 10.1099/ijs.0.64238-0
Pantos, O., Bongaerts, P., Dennis, P. G., Tyson, G. W., and Hoegh-Guldberg, O. (2015). Habitat-specific environmental conditions primary control the mcrobiomes of the coral Seriatopora hystrix. ISME J. 9, 1916–1927. doi: 10.1038/ismej.2015.3
Pantos, O., Cooney, R. P., Le Tissier, M. D. A., Barer, M. R., O'Donnell, A. G., and Bythell, J. C. (2003). The bacterial ecology of a plague-like disease affecting the Caribbean coral Montastrea annularis. Environ. Microbiol. 5, 370–382. doi: 10.1046/j.1462-2920.2003.00427.x
Pinhassi, J., Zweifei, U. L., and Hagström, Å. (1997). Dominant marine bacterioplankton species found among colony-forming bacteria. Appl. Environ. Microbiol. 63, 2259–3366.
Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS ONE 5:e9490. doi: 10.1371/journal.pone.0009490
Quince, C., Lanzen, A., Davenport, R. J., and Turnbaug, P. J. (2011). Removing noise from pyrosequenced amplicons. BMC Bioinformatics 12:38. doi: 10.1186/1471-2105-12-38
Quinn, R. A., Metzler, A., Smoloqitz, R. M., Tlusty, M., and Chistoserdov, A. Y. (2012). Exposures of homarus americanus shell to three bacteria isolated from naturally occurring epizootic shell disease lesions. J. Shellfish Res. 31, 485–493. doi: 10.2983/035.031.0208
R Core Team (2012). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Reeder, J., and Knight, R. (2010). Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat. Meth. 7, 668–669. doi: 10.1038/nmeth0910-668b
Reshef, L., Koren, O., Loya, Y., Zilber-Rosenberg, I., and Rosenberg, E. (2006). The coral probiotic hypothesis. Environ. Microbiol. 8, 2068–2073. doi: 10.1111/j.1462-2920.2006.01148.x
Richie, K. B. (2006). Regulation of microbial populations by coral surface mucus and mucus-associated bacteria. Mar. Ecol. Prog. Ser. 332, 1–14. doi: 10.3354/meps322001
Roder, C., Arif, C., Daniels, C., Wil, E., and Voolstra, C. R. (2014). Bacterial profiling of white plague disease across corals and oceans indicates a conserved and distinct disease microbiome. Mol. Ecol. 23, 965–975. doi: 10.1111/mec.12638
Romero, M., Avendaño-Herrera, R., Magariños, B., Cámara, M., and Otero, A. (2010). Acylhomoserine lactone production and degradation by the fish pathogen Tenacibaculum maritimum, a member of the Cytophaga-Flavobacterium-Bacteroides (CFB) group. FEMS Microbiol. Lett. 304, 131–139. doi: 10.1111/j.1574-6968.2009.01889.x
Rosenberg, E., and Ben-Haim, Y. (2002). Microbial diseases of corals and global warming. Environ. Microbiol. 4, 318–326. doi: 10.1046/j.1462-2920.2002.00302.x
Rosenberg, E., Koren, O., Reshef, L., Efrony, R., and Zilber-Rosenberg, I. (2007). The role of microorganisms in coral health, disease and evolution. Nat. Rev. Microbiol. 5, 355–362. doi: 10.1038/nrmicro1635
Sato, Y., Willis, B. L., and Bourne, D. G. (2010). Successional changes in bacterial communities during the development of black band disease on the reef coral, Montipora hispida. ISME J. 4, 203–214. doi: 10.1038/ismej.2009.103
Sekar, R., Kaczmarsky, L. T., and Richardson, L. L. (2008). Microbial community composition of black band disease on the coral host Siderastrea siderea from three regions of the wider Caribbean. Mar. Ecol. Prof. Ser. 362, 85–98. doi: 10.3354/meps07496
Sekar, R., Mills, D. K., Remily, E. R., Voss, J. D., and Richardson, L. L. (2006). Microbial communities in the surface mucopolysaccharide layer and the black band microbial mat of black band-diseased Siderastrea siderea. Appl. Environ. Microbiol. 72, 5963–5973. doi: 10.1128/AEM.00843-06
Séré, G. M., Tortosa, P., Chabanet, P., Turquet, J., Quod, J. P., and Schleyer, M. H. (2013). Bacterial communities associated with Porites white patch syndrome (PWPS) on three western indian ocean (WIO) coral reefs. PLoS ONE 8:e83746. doi: 10.1371/journal.pone.0083746
Stimson, J. (2011). Ecological characterization of coral growth anomalies on Porites compressa in Hawaii. Coral Reefs 30, 133–142. doi: 10.1007/s00338-010-0672-8
Sunagawa, S., DeSantis, T. Z., Piceno, Y. M., Brodie, E. L., DeSalvo, M. K., Voolstra, C. R., et al. (2009). Bacterial diversity and White Plague Disease-associated community changes in the Caribbean coral Montastraea faveolata. ISME J. 3, 512–521. doi: 10.1038/ismej.2008.131
Sunagawa, S., Woodley, C. M., and Medina, M. (2010). Threatened corals provide underexplored microbial habitats. PLoS ONE 5:e9554. doi: 10.1371/journal.pone.0009554
Sussman, M., Willis, B. L., Victor, S., and Bourne, D. G. (2008). Coral pathogens identified for White Syndrome (WS) epizootics in the Indo-Pacific. PLoS ONE 3:e2393. doi: 10.1371/journal.pone.0002393
Sutherland, K. P., Porter, J. W., and Torres, C. (2004). Disease and immunity in Caribbean and Indo-Pacific zooxanthellate corals. Mar. Ecol. Prog. Ser. 266, 273–302. doi: 10.3354/meps266273
Sweatman, H., Delean, S., and Syms, C. (2011). Assessing loss of coral cover on Australia's Great Barrier Reef over two decades, with implications for longer-term trends. Coral Reefs 30, 521–531. doi: 10.1007/s00338-010-0715-1
Sweet, M. J., Croquer, A., and Bythell, J. C. (2014). Experimental antibiotic treatment identifies potential pathogens of white band disease in the endangered Caribbean coral Acropora cervicornis. Proc. R. Soc. B. 281:20140094. doi: 10.1098/rspb.2014.0094
Tavakoli-Kolour, P., Kavousi, J., and Rezai, H. (2015). Outbreak of growth anomalies in coral communities of Qeshm Island, Persian Gulf. Int. Aquat. Res. 7, 151–156. doi: 10.1007/s40071-015-0100-3
Thompson, F. L., Barash, Y., Sawabe, T., Sharon, G., Swings, J., and Rosenberg, E. (2006). Thalassomonas loyana sp. nov., a causative agent of the white plague-like disease of corals on the Eilat coral reef. Int. J. Syst. Evol. Microbiol. 56, 365–368. doi: 10.1099/ijs.0.63800-0
Tian, C., Liu, X., Tan, J., Lin, S., Li, D., and Yang, H. (2012). Isolation, identification and characterization of an algicidal bacterium from Lake Taihu and preliminary studies on its algicidal compounds. J. Environ. Sci. (China) 24, 1823–1831. doi: 10.1016/S1001-0742(11)60983-2
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
White, J. R., Nagarajan, N., and Pop, M. (2009). Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput. Biol. 5:e1000352. doi: 10.1371/journal.pcbi.1000352
Willis, B. L., Page, C. A., and Dinsdale, E. A. (2004). “Coral disease on the great barrier reef,” in Coral Health and Disease, eds E. Rosenberg and Y. Loya (Berlin: Springer-Verlag), 69–104.
Wilson, B., Aeby, G. S., Work, T. M., and Bourne, D. G. (2012). Bacterial communities associated with healthy and Acropora white syndrome-affected corals from American Samoa. FEMS Microbiol. Ecol. 80, 509–520. doi: 10.1111/j.1574-6941.2012.01319.x
Keywords: Platygyra carnosus, coral disease, skeletal growth anomalies, microbial communities, metagenomics
Citation: Ng JCY, Chan Y, Tun HM, Leung FCC, Shin PKS and Chiu JMY (2015) Pyrosequencing of the bacteria associated with Platygyra carnosus corals with skeletal growth anomalies reveals differences in bacterial community composition in apparently healthy and diseased tissues. Front. Microbiol. 6:1142. doi: 10.3389/fmicb.2015.01142
Received: 05 June 2015; Accepted: 05 October 2015;
Published: 20 October 2015.
Edited by:
Hongyue Dang, Xiamen University, ChinaReviewed by:
Fabiano Thompson, Federal University of Rio de Janeiro, BrazilMelissa Garren, Pelagic Data Systems and California State University Monterey Bay, USA
Copyright © 2015 Ng, Chan, Tun, Leung, Shin and Chiu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jill M. Y. Chiu, jillchiu@hkbu.edu.hk