Skip to main content
  • Research Article
  • Open access
  • Published:

Parrotfish predation drives distinct microbial communities in reef-building corals

Abstract

Background

Coral-associated microbial communities are sensitive to multiple environmental and biotic stressors that can lead to dysbiosis and mortality. Although the processes contributing to these microbial shifts remain inadequately understood, a number of potential mechanisms have been identified. For example, predation by various corallivore species, including ecologically-important taxa such as parrotfishes, may disrupt coral microbiomes via bite-induced transmission and/or enrichment of potentially opportunistic bacteria. Here, we used a combination of mesocosm experiments and field-based observations to investigate whether parrotfish corallivory can alter coral microbial assemblages directly and to identify the potentially relevant pathways (e.g. direct transmission) that may contribute to these changes.

Results

Our mesocosm experiment demonstrated that predation by the parrotfish Chlorurus spilurus on Porites lobata corals resulted in a 2-4x increase in bacterial alpha diversity of the coral microbiome and a shift in bacterial community composition after 48 h. These changes corresponded with greater abundance of both potentially beneficial (i.e. Oceanospirillum) and opportunistic bacteria (i.e. Flammeovirgaceae, Rhodobacteraceae) in predated compared to mechanically wounded corals. Importantly, many of these taxa were detectable in C. spilurus mouths, but not in corals prior to predation. When we sampled bitten and unbitten corals in the field, corals bitten by parrotfishes exhibited 3x greater microbial richness and a shift in community composition towards greater abundance of both potential beneficial symbionts (i.e. Ruegeria) and bacterial opportunists (i.e. Rhodospiralles, Glaciecola). Moreover, we observed 4x greater community variability in naturally bitten vs. unbitten corals, a potential indicator of dysbiosis. Interestingly, some of the microbial taxa detected in naturally bitten corals, but not unbitten colonies, were also detected in parrotfish mouths.

Conclusions

Our findings suggest that parrotfish corallivory may represent an unrecognized route of bacterial transmission and/or enrichment of rare and distinct bacterial taxa, both of which could impact coral microbiomes and health. More broadly, we highlight how underappreciated pathways, such as corallivory, may contribute to dysbiosis within reef corals, which will be critical for understanding and predicting coral disease dynamics as reefs further degrade.

Background

Reef-building corals host a wide range of microorganisms including endosymbiotic dinoflagellates (Symbiodinaceae), viruses, archaea, and bacteria that collectively comprise the coral holobiont [1]. The relationship between corals and these microbial associates allows reef corals to thrive in nutrient-poor waters and to support high biodiversity [2]. While the coral-dinoflagellate symbiosis is well documented, less is known about the roles of coral-associated bacterial communities in coral health and resilience [3,4,5]. However, an array of mutualistic benefits are suggested, encompassing vital functions such as coral nutrition and immunity that may further impact fundamental ecological processes within coral reefs [4, 6, 7].

Despite their beneficial role in host fitness, coral-associated bacteria are sensitive to numerous environmental and biotic stressors that may lead to microbial dysbiosis (i.e., a shift in either the mean composition or variability of the microbiome, including the loss of beneficial symbionts and/or increase of opportunists) [8, 9]. However, the processes that contribute to dysbiosis remain inadequately understood [8, 10,11,12]. For example, trophic interactions that are common within reef ecosystems, such predation on corals (i.e. corallivory), may favor the disruption of coral microbiomes and subsequently alter coral health [6, 8, 9, 13].

A variety of corallivores, including invertebrates such as fireworms [14], echinoderms [15], and gastropods [16, 17], may serve as reservoirs and/or vectors of opportunistic bacteria (i.e. typically non-pathogenic microorganisms that take advantage of their host under certain circumstances) or pathogens (i.e. microorganisms that cause infection) to corals. Corallivorous, polyp-feeding butterflyfishes may also spread microbes and parasites via their mouthparts [14,15,16,17,18]. However, the fact that butterflyfishes remove limited coral tissue without exposing the underlying skeleton may make them unlikely candidates to transmit microbes among individual corals [16]. In contrast, many parrotfishes scrape or excavate both live coral tissue and skeleton, while also ingesting detritus and turf algae from rocky surfaces [19]. This invasive feeding method and more varied diet make them likely candidates for disrupting coral microbiomes. Parrotfishes play a critical role in structuring the benthic communities of coral reefs and are generally considered to have a net positive effect on ecosystem functions, promoting coral dominance by removing competing macroalgae [19] and acting as important agents of reef bioerosion [20]. That said, parrotfish predation is a chronic stress that may cause significant harm to corals, especially when coupled with other environmental and/or biotic stressors. For example, recent findings suggest that corals exposed to both parrotfish predation and nutrient pollution experienced significantly greater mortality than when exposed to either stressor alone, which was attributed to increased bacterial opportunism [13]. Identifying the mechanisms and conditions in which parrotfishes can disrupt coral microbiomes will improve our ability to predict the microbial impacts associated with corallivory and their potential implications for coral health.

Here, we conducted a series of experiments in Mo’orea, French Polynesia, to assess the effects of parrotfish corallivory on coral microbiomes. Specifically, we tested whether parrotfish can facilitate the enrichment and/or transmission of microbes to corals. We focused on Chlorurus spilurus, a common parrotfish species known to prey on large colonies of Porites on Pacific reefs [21]. We first performed a controlled mesocosm experiment comparing the microbiomes of Porites lobata colonies that were either mechanically wounded or predated by C. spilurus parrotfish. Comparisons were based on samples collected immediately following predation or wounding (Ti) and at 48 h (Tf). To examine microbiome patterns in situ, we collected microbial samples from corals in the field that had either been naturally bitten by parrotfishes or appeared bite free (hereafter “unbitten”). Coral microbiomes from both the mesocosm experiment and field survey were also compared to microbiomes from parrotfish mouthparts to assess potential predation-mediated pathways (e.g. transmission) that may contribute to coral microbiome change. We hypothesized that predation by C. spilurus facilitates the enrichment and/or transmission of microbes to corals, resulting in the following impacts on the coral microbiome: (1) increases in alpha diversity, (2) differences in community composition (3) increases in community variability, and (4) increased abundance of microbial taxa typically found in the mouths of parrotfish but absent in healthy corals.

Results

Experimental overview

The impacts of parrotfish corallivory on coral microbiomes were assessed using a combination of manipulative experiments and field surveys on the north shore of Mo’orea, French Polynesia. First, to assess the ability of C. spilurus to feed on live Porites lobata corals, we conducted a survey in two back reef areas to quantify the number of C. spilurus bites found on live versus dead corals. A manipulative experiment was then performed at the UC Gump Marine station with seven C. spilurus initial phase individuals and ten colonies of Porites lobata (about 20 cm) that were previously collected on a nearby back reef. When at the station, sterile culture swabs were used to collect microbial samples of the mouths of each C. spilurus individual. Each C. spilurus was then coaxed into biting a single Porites lobata colony at two separate locations. The remaining three P. lobata colonies were artificially wounded at two locations using a sterilized bone cutter. Coral samples (mucus, tissue and part of the coral skeleton) were collected: i) immediately following parrotfish biting/mechanical wounding (Ti) and ii) at the end of the experiment (48 h, Tf). Coral samples and fish swabs were stored at − 80 °C prior to laboratory analyses.

To compare our experimental results with conditions in the field, we haphazardly collected P. lobata samples (mucus, tissue and part of the skeleton) that were either unbitten or naturally bitten (n = 10 per status) from a reef on Mo’orea’s north shore. Sterile swabs were also used to collect microbial samples from ten C. spilurus collected haphazardly from the same reef. Four 1 L-water samples were also collected and directly filtered on 0.2 um filters. When on the boat, coral tissues, swabs and filters were placed on ice and stored at − 80 °C immediately upon arrival at the marine station.

DNA extractions on all samples were performed using DNeasy PowerSoil Kit (Qiagen) and AccuStart II PCR ToughMix (Quanta BioSciences, Gaithersburg, Maryland, USA) was used to perform two-step Polymerase Chain Reaction (PCR) on the V4 hypervariable region of the 16S rRNA gene. Data processing and analyses of microbiome diversity, composition and stability metrics of corals, fish mouths and water were then performed using the Delbur workflow [22], QIIME2 [23] pipelines and R [24] for statistical analyses.

Parrotfish feed on live corals in the field

Among the 23 individual fish that were followed, we recorded a total of 5451 bites on either live or dead corals (including rubble and pavement). Of these, 5400 (99%) were taken from dead corals and 51 (~ 1%) were taken from live corals, with the latter comprising mostly massive Porites colonies (49 bites, 96% of the bites on live coral were taken on massive Porites).

Bacterial assemblages differed between parrotfish mouth, coral, and water samples

Mesocosm experiment

Parrotfish mouths showed distinct bacterial communities compared to mechanically wounded corals both at Ti (Additional file 1: Tables S1 and S2; pairwise Adonis; p = 0.02) and Tf (pairwise Adonis; p = 0.012) in the mesocosm experiment. Phylum level assignments in parrotfish mouths showed the dominance of by Proteobacteria and Bacteroidetes (Additional file 2: Figure S1). Among the 83 families characterized in parrotfish mouths, few were present at proportions greater than 1%. These included Flavobacteriaceae (23.1%), Alteromonadaceae (15.3%), Rhodobacteraceae (8.3%), Pseudoalteromonadaceae (5.8%), and Vibrionaceae (5.1%) (Fig. 1; Additional file 1: Table S3). Microbiomes of mechanically wounded corals were populated by Proteobacteria and Bacteroidetes (Additional file 1: Tables S4 and S5), while predated corals were mainly dominated by Proteobacteria (Additional file 1: Tables S6 and S7). More specifically, the common coral symbiont Hahellaceae dominated bacterial communities in mechanically wounded corals at Ti (83.9%) and Tf (59.9%) (Fig. 2; Additional file 1: Tables S4 and S5). Distinct community assemblages were also observed between parrotfish mouths and predated corals at Ti (Fig. 2; Additional file 1: Tables S1 and S2; pairwise Adonis; p = 0.003) and at Tf (pairwise Adonis, p = 0.012). Among the 66 (Ti) and 49 (Tf) bacterial families identified in bitten corals at Ti and Tf, respectively, only a few were represented at relative proportions greater than 1%. These included Hahellaceae (21.6%), Amoebophilaceae (17.5%), and Rivulariaceae (9.7%) (Additional file 1: Table S6) at Ti, and Rhodobacteraceae (13.7%), Pseudoalteromonadaceae (11.0%), Verrucomicrobiaceae (5.9%), Alteromonadaceae (5.5%), Flavobacteriaceae (3.9%), Vibrionaceae (3.9%), Oceanospirillaceae (3.3%), Colwelliaceae (2.5%), Lentisphaeraceae (2.4%), Francisellaceae (1.6%), Paenibacillaceae (1.4%), and Hahellaceae (12.2%) at Tf (Additional file 1: Table S7).

Fig. 1
figure 1

a Bullethead parrotfish Chlorurus spilurus (photo credit: Katrina Munsterman). b Unbitten colonies of Porites lobata (photo credit: Cody Clements) and c naturally bitten colonies by parrotfish as found in our study site in the back reef area of Mo’orea, French Polynesia (photo credit: Mallory Rice)

Fig. 2
figure 2

Heat maps displaying the relative abundance (expressed as proportion) of the 25 most abundant taxa grouped at the family level or to the closest taxonomic rank according to the sample type for (a) the mesocosm experiment (at Tf) and (b) the field survey. P-values represent significant differences, based on pairwise comparisons using the pairwise.adonis function, in microbial community composition between a predated and mechanically wounded corals in the mesocosm experiment or b bitten and unbitten corals in the field

Field survey

Similar to our mesocosm experiment, microbiomes of parrotfish mouths were dominated by Proteobacteria and Bacteroidetes (Additional file 2: Figure S2) and showed distinct bacterial communities compared to naturally unbitten corals in the field (Additional file 1: Tables S8 and S9; pairwise Adonis; p = 0.002). Among the 99 bacterial families identified in fish mouths, only 7 were represented at moderate abundance (1–20%), including Moraxellaceae (16%), Alteromonadaceae (8.5%), Rhodobacteraceae (6.5%), Vibrionaceae (4.6%), Flavobacteriaceae (4.3%), Rhodospirillaceae (4.3%) and Paenibacillaceae (3.3%) (Fig. 2; Additional file 1: Table S10). In naturally unbitten corals, 51 families were identified (Additional file 1: Table S11), and as with the mechanically wounded corals in our mesocosm experiment, microbiomes were mainly populated by Hahellaceae (89.9%) (Additional file 1: Table S11). Differences in bacterial community composition were also observed between parrotfish mouths and naturally bitten corals (Fig. 2; Additional file 1: Tables S8 and S9; pairwise Adonis; p = 0.0015. In naturally bitten corals, families with relative abundances exceeding 1% included Hahellaceae (63.4%), Rhodobacteraceae (7.3%), Paenibacillaceae (3.1%), Flavobacteriaceae (1.9%), Rhodospirillaceae (1.6%), Moraxellaceae (1.4%), Alteromonadaceae (1.4%), and Flammeovirgaceae (1.3%) (Additional file 1: Table S12). Finally, we found that bacterial assemblages of water samples significantly differed from parrotfish mouths (Additional file 1: Table S9; pairwise Adonis; p = 0.002), as well as unbitten (pairwise Adonis; p = 0.003) and bitten corals (pairwise Adonis; p = 0.003), suggesting that changes in the coral’s microbial composition are not solely driven by microbial communities in the surrounding reef environment.

Parrotfish predation increased alpha diversity of P. lobata microbiomes

Parrotfish predation induced significant changes in the alpha diversity of P. lobata microbiomes, both in the mesocosm and in the field. Overall, parrotfish mouth microbiomes showed greater bacterial richness compared to mechanically wounded and predated corals in the mesocosm experiment at both time points (Additional file 1: Tables S13 and S14; pTi = 0.002 and pTf = 0.002) and to unbitten corals in the field (Additional file 1: Tables S15 and S16, p = 0.003). However, Shannon-Wiener bacterial diversity did not significantly differ between parrotfish mouths and both predated corals in the mesocosm (Additional file 1: Table S14, pTi = 0.17, pTf = 0.18) and naturally bitten corals in the field (Additional file 1: Table S16, p = 0.1). While patterns of alpha diversity in the mesocosm were similar between mechanically wounded and predated corals at Ti (Fig. 3a, b; Additional file 1: Tables S13 and S14; Richness – p = 0.15; Shannon – p = 0.13), predated corals at Tf exhibited 2x greater microbial richness (73.4 ± 11) and 4x greater Shannon-Wiener diversity (3.1 ± 0.2) compared to mechanically wounded corals (Richness: 32.6 ± 8.4 and Shannon: 0.72 ± 0.13) (Fig. 3a, b; Additional file 1: Table S14; Richness – p = 0.049; Shannon – p < 0.001). In the field, naturally bitten corals exhibited 3x greater microbial richness (62.1 ± 26.9) and diversity (1.8 ± 0.5) compared to unbitten corals, although only differences in richness were significant (Richness: 19.25 ± 1.8 and Shannon: 0.58 ± 0.1) (Fig. 3c, d; Additional file 1: Tables S15 and S16; Richness – p = 0.04; Shannon – p = 0.08).

Fig. 3
figure 3

Alpha diversity metrics. Observed richness and Shannon-Wiener diversity indices for the mesocosm experiment (a, b) and field survey (c, d) for each sample type and timepoint. P-values represent pairwise comparisons, using Tukey’s or Dunn’s test, at each time point between a, b predated and mechanically wounded corals, c, d bitten and unbitten corals

Corallivory generated a distinct microbiome community structure in P. lobata

Evidence for changes in microbial community composition following parrotfish predation was present in both our mesocosm experiment and field survey (Fig. 2, Additional file 2: Figures S3 and S4). In the former, predated and mechanically wounded corals exhibited similar patterns in their overall microbial community assemblages at Ti (Additional file 1: Table S2, pairwise Adonis; p = 0.07). However, five sub-operational taxonomic units (sOTUs) had greater abundance in predated corals when compared to mechanically wounded corals at Ti. These included members of the Rivulariaceae (genus Rivularia; sOTU_15), Phormidiaceae (sOTU_12) and Amoebophilaceae (clade SGUS912; sOTU_195) families, as well as two taxa from the orders Nostocales (sOTU_18) and Rhizobiales (sOTU_697) (Additional file 1: Table S17; log2 fold-change 7.05 to 22.9). Among these, three taxa were identified only in predated corals (sOTU_15, sOTU_12, sOTU_697; Additional file 1: Table S18). One sequence (sOTU_195) was found both in mechanically wounded and predated corals (Additional file 1: Table S18). Moreover, the specific sOTU_18 was found both in predated corals and in low abundance (0.12%) in fish mouths, (Additional file 1: Table S18) but not in mechanically wounded corals.

At 48 h, mechanically wounded and predated corals exhibited significantly different bacterial community composition (Fig. 2 and Additional file 2: Figure S3; Additional file 1: Table S2; pairwise Adonis; p = 0.018). This coincided with a greater abundance of four sOTUs (Fig. 4a; Additional file 1: Table S19; log2 fold-change 4.6 to 7.6), including members of the families Rhodobacteraceae (genus Phaeobacter; sOTU_771), Oceanospirillaceae (genus Oceanospirillum; sOTU_467), and Lentisphaeraceae (sOTU_39), and the order Rhodospirillales sOTU_480). Of these taxa, two were absent from mechanically wounded corals, but present in relatively low abundance in predated corals (sOTU_771–1.67%; sOTU_467–0.9%), as well as parrotfish mouths (sOTU_771–0.5%; sOTU_467–0.012%) (Additional file 1: Table S18). One taxon (sOTU_480) was found in low abundance in mechanically wounded corals (0.2%) and at moderate levels in predated corals (13%; Additional file 1: Table S18). Finally, sOTU_39 was present at low abundance in mechanically wounded corals (0.04%), but was more abundant in predated corals (2.3%) and parrotfish mouths (2.8%; Additional file 1: Table S18).

Fig. 4
figure 4

Differential abundance analysis (DESeq2) illustrating the sOTUs labeled as genera and families or order that differed significantly between (a) predated and mechanically wounded corals for the mesocosm experiment at Tf (48 h) and (b) naturally unbitten and bitten corals in the field survey. Fish illustration indicates the presence of the corresponding taxa in parrotfish mouths (image credit: Katrina Munsterman)

In the field, naturally bitten and unbitten corals showed distinct patterns in bacterial community composition (Fig. 2 and Additional file 2: Figure S4, Additional file 1: Table S9; pairwise Adonis; p = 0.01). This coincided with a greater abundance of six sOTUs in naturally bitten compared to unbitten corals (Fig. 4b; Additional file 1: Table S20; log2 fold-change 4.9 to 22.7). Among these, three taxa were undetectable in unbitten corals but present in relatively low abundance in parrotfish mouths, including sequences from the families Flammeovirgaceae (genus JTB248 – sOTU_1051–0.44%; Additional file 1: Table S21), Rhodobacteraceae (genus Glaciecola – sOTU_2451–1.7%) and Alteromonadaceae (sOTU_3182–0.27%). Two members from the Rhodobacteraceae family (sOTUs_3439, 3450) were identified in naturally bitten and unbitten corals, as well as in fish mouths (sOTU_3439–1.01%; sOTU_3450–2.41%; Additional file 1: Table S21), while a sequence from the family Alteromonadaceae (sOTU_2321) was found only in naturally bitten corals (0.3%; Additional file 1: Table S21). Finally, we observed 4x greater community variability in the microbiome of naturally bitten corals (0.32 ± 0.04) compared to unbitten corals (0.07 ± 0.005; Fig. 5; Additional file 1: Tables S22 and S23; Tukey HSD - p < 0.001).

Fig. 5
figure 5

Boxplots illustrating the level of community variability among sample types for the field experiment. P-values represent pairwise comparisons, using Tukey’s test, between unbitten and bitten corals

Discussion

A number of corallivores are suspected to facilitate the enrichment and/or transmission of microbes within reef-building corals [14, 16], including consumers such as parrotfishes that play key roles in regulating reef ecosystem processes [13]. Using a combination of mesocosm- and field-based approaches, we demonstrated that corallivory by the parrotfish species Chlorurus spilurus leads to significant changes in bacterial community composition of Porites lobata. In particular, these changes included greater abundances of potential beneficial bacterial taxa and opportunists, some of which were naturally occurring in parrotfish mouths. Our findings indicate that parrotfishes may play an important role in driving the structure of coral microbial communities, either by acting as vectors and/or by facilitating the enrichment of bacteria in reef corals via corallivory.

Parrotfish-induced P. lobata microbiome changes in mesocosm

Patterns of alpha and beta diversity in our mesocosm experiment were similar at Ti for mechanically wounded and predated corals. However, five taxa were already observed in greater abundance in predated corals compared to those that were mechanically wounded. All were present at relatively low abundances (< 10%) in predated corals, but may have the potential to affect coral microbiomes and health. For instance, cyanobacteria from the Nostocales order (sOTU_18) are often found in fish guts [25] and were associated with diseased corals [26]. Members of the clade SGUS912 (sOTU_195) are commonly present in corals exposed to sewage and wastewater outfalls [27]. Taxa from the orders Oscillatoriales (sOTU_12) and Rhizobiales (sOTU_697), and filamentous Cyanobacteria from the genus Rivularia (sOTU_15), were associated with stressed and diseased corals and sponges [9, 28,29,30]. Whether and how these changes affect coral health and fitness, especially when coupled with other stressors, should be investigated further.

At the end of the experiment (Tf), we observed greater bacterial richness and diversity in predated corals compared to mechanically wounded corals. Patterns of increased alpha diversity are often associated with numerous physical and biotic stressors including water pollution [31, 32], elevated temperature [33, 34], ocean acidification [35], algal competition [36, 37], mechanical wounding, and snail corallivory [38, 39]. However, other studies demonstrated no changes or a significant decrease in microbial diversity and/or richness following mechanical injury [38, 40]. These differences among studies may indicate that responses of coral microbiomes differ due to biological vs. mechanical wounding, or that stressor-induced impacts may be variable depending on coral species or genotypes, local environmental conditions, and/or exposure time. In the present study, increases in bacterial richness and diversity coincided with a compositional shift in bacterial assemblages in predated corals compared to mechanically wounded ones. In addition, microbiomes of predated corals were characterized by moderate abundance (59.9% at Tf) of the putative beneficial symbiont Hahellaceae at 48 h when compared to Ti (83.9%). Lower abundance of Hahellaceae bacterial taxa is a pattern previously reported in stressed, mechanical injured, and predated corals [9, 38, 39, 41].

In addition, bacterial communities of corals exposed to predation were dominated by members of the families Rhodobacteraceae, Pseudoalteromonadaceae, Alteromonadaceae, Verrucomicrobiaceae and Flavobacteriaceae – taxa that are often associated with both stressed and healthy coral colonies [32, 42], and were also found in relatively high abundance in parrotfish mouths. Four sOTUs were present in greater abundance among predated corals compared to mechanically wounded ones, including taxa from the genera Phaeobacter (sOTU_771) and Oceanospirillum (sOTU_467), as well as sequences from the Lentisphaerae (sOTU_39) and Rhodospirillales (sOTU_480) orders. Their potential influences on corals may be diverse – ranging from beneficial to opportunistic. Members of the genus Phaeobacter were previously found in corals and jellyfish [43,44,45] and were linked to the production of antibacterial compounds in fishes [46, 47]. Bacteria from the genus Oceanospirillum are frequently observed in healthy coral colonies [48, 49], while members of the phylum Lentisphaerae are common in the fish gut [50] and healthy corals [48]. Sequences from the order Rhodospirillales are commonly found in high abundance in stressed and diseased coral colonies [9, 51,52,53], indicating an opportunistic character. Given that our experiment lasted for 48 h, the persistence of potential beneficial symbionts and opportunistic bacterial taxa and their consequences on coral microbiomes and health will have to be further investigated over longer time period.

Microbiomes of naturally bitten vs. unbitten P. lobata in the field

Microbiomes of P. lobata corals found in the field reinforced findings from our mesocosm, as naturally bitten corals exhibited greater bacterial richness compared to unbitten corals. It is worth noting that corals of all treatments, from both the mesocosm experiment and field survey, exhibited relatively low bacterial richness compared to previous work [54]. However, lower richness has consistently been observed among corals inhabiting reefs in Mo’orea [39, 55] and may be related to greater community dominance by members of the Hahellaceae family.

In our study, increased bacterial richness in naturally bitten was associated with a shift in bacterial community composition compared to unbitten corals. Bitten corals were mainly populated by potential opportunistic bacterial taxa, including sequences from the families Rhodobacteraceae, Paenibacillaceae, Flavobacteriaceae, Rhodospirillaceae, Moraxellaceae, Alteromonadaceae, and Flammeovirgaceae [9, 56,57,58], all of which are taxa that were also present in relatively high abundance in parrotfish mouths. Significant changes in community composition were associated with greater abundance of six taxa in naturally bitten vs. unbitten corals. Among them, three represented the Flammeovirgaceae (genus JTB248; sOTU_1051), and Alteromonadaceae (sOTU_2321; sOTU_3182;) families, that were previously associated with stressed, aged, and/or diseased corals [59,60,61]. Three other taxa were assigned to the Rhodobacteraceae family, taxa commonly associated with both healthy and stressed corals [56, 62]. In particular, a strain from the genus Ruegeria was found in lesioned and diseased corals [41, 56] and is known to inhibit growth of the coral pathogen Vibrio coralliilyticus [63].

As with our mesocosm experiment, our field survey identified taxa with potential beneficial and/or deleterious implications for coral microbiomes, health, and fitness. Further studies are needed to understand the functional roles of these microbes and their interplay with coral stressors. Finally, we observed greater bacterial compositional variability in naturally bitten compared to unbitten corals. Increased microbiome variability is consistent with previous studies showing that numerous animals, including corals, exhibit elevated community variability when exposed to stressors (i.e. the Anna Karenina Principle [64];), such as mechanical wounding [39]. This further indicates the potential for parrotfish to increase dysbiosis susceptibility in corals via corallivory.

Potential parrotfish-mediated bacterial transmission and enrichment pathways in P. lobata

Parrotfish predation could alter the microbiomes of P. lobata via several pathways, including i) direct transmission of bacteria from fish mouths to the coral mucus/tissue layer, ii) indirectly facilitating bacterial invasion from the surrounding environment following wounding, iii) indirectly facilitating growth of bacterial taxa already present within the coral microbiomes or from the surrounding environment, and iv) a combination of these three pathways. We observed evidence for each of these possible pathways in our experiments. For example, evidence that parrotfish may directly transmit bacteria to P. lobata was observed both in the mesocosm experiment and field survey. In the former case, a taxon from the order Nostocales (sOTU_18) at Ti, as well as two taxa from the genera Phaeobacter (sOTU_771) and Oceanospirillum (sOTU_467) at Tf, were both found in predated corals and fish mouths, but not in mechanically wounded corals. This indicates that mechanical wounding was insufficient to introduce these taxa and that they were likely vectored via parrotfish predation. Similar patterns were observed in the field, with sequences from the families Flammeovirgaceae (sOTU_1051; genus JTB248), Rhodobacteraceae (sOTU_2451; genus Glaciecola), and Alteromonadaceae (sOTU_3182) present only in bitten corals and fish mouths – not unbitten corals. We also observed evidence that predation may facilitate the invasion of bacterial taxa from the surrounding environment. At Ti in our mesocosm experiment, three potential opportunistic bacterial coral taxa were found in predated corals, but not in mechanically wounded corals or parrotfish mouths (sOTU_12, sOTU_15 and sOTU_697). Similarly, in the field, sequences from the family Alteromonadaceae (sOTU_2321) were only found in naturally bitten corals, indicating enrichment from the surrounding environment.

We also observed potential enrichment from microbes preexisting on corals and/or from the external environment, such as members from the clade SGUS912 (sOTU_195) and the Rhodospirillales order (sOTU_480), which were identified in predated and mechanically wounded corals – but not fish mouths – at Ti and Tf, respectively. Finally, evidence from both experiments suggested a combination of different pathways including transmission and/or enrichment. In the manipulative experiment, taxon sOTU_39 from the Lentisphaerae order was present in moderate abundance in predated corals, as well as in low abundance in mechanically wounded corals and parrotfish mouths at Tf. In the field, two taxa from the Rhodobacteraceae family (sOTU_3439, sOTU_3450) were present in fish mouths and bitten corals, as well as in relatively low abundance in unbitten corals.

Collectively, our findings suggest that parrotfish corallivory may be an important driver structuring coral-associated bacterial communities. Evidence that parrotfish vector and/or facilitate the enrichment of bacteria within corals, both in our mesocosm experiment and field surveys, was surprisingly consistent – especially given that sampling of corals and parrotfish mouths was conducted haphazardly in the back reef during our field surveys. This suggests that parrotfish mouths may harbor a consistent microbial signature in the studied reef area that allows C. spilurus to vector rare taxa via corallivory. Our findings add to growing body of evidence demonstrating the potential for corallivores, such as snails (Drupella spp., Coralliophila spp.), crown-of-thorn sea stars (Acanthaster spp.) and worms (Hermodice caniculata), to vector and/or facilitate the enrichment of microbes in corals [38, 65, 66]. Our study is the first to document such potential in parrotfishes, adding to their key roles as corallivores, bioeroders, and herbivores on coral reefs. Previous work suggests that other candidate species, such as butterflyfishes, are unlikely to vector microbes [16, 67] – potentially due to their distinct “browser” feeding mode (but see [17]). In contrast, the “scraper” and “excavator” feeding modes of many parrotfishes may make them ideal candidates to transmit microbes to corals. The interplay between these abilities and the other critical roles of parrotfishes on coral reefs will be of considerable interest for reefs of the future.

Conclusion

Our findings provide evidence that parrotfish corallivory can have important effects on coral microbiomes, with the potential to impact coral health. C. spilurus predation both in the laboratory and field induced an increase in alpha diversity and a compositional shift in the microbial assemblages of P. lobata corals, which coincided with a greater abundance of potential beneficial bacteria (i.e. Ruegeria, Phaeobacter) as well as opportunistic taxa (i.e Flammeovirgaceae, Rhodospirillaleceae, Glaciecola). Importantly, several taxa were undetectable on mechanically wounded and naturally unbitten corals but present in predated, naturally bitten corals and in parrotfish mouths, suggesting parrotfish vector new bacteria to corals during predation. However, the ability of C. spilurus to vector and/or facilitate the enrichment of microbial opportunists, as well as increase microbiome variability, in naturally bitten P. lobata corals is consistent with recent findings linking nutrient pollution and parrotfish predation to coral mortality [13]. This suggests that common trophic interactions may increase coral susceptibility to dysbiosis, especially when corals are already stressed from other factors such as nutrient pollution, temperature, or sedimentation. Together, our results shed light on underappreciated pathways linking parrotfishes to microbial enrichment and dysbiosis within reef corals. Future work should investigate the interactive effects of parrotfish corallivory and abiotic stressors (e.g. nutrient pollution and ocean warming) to evaluate their consequences for coral microbiomes and fitness.

Material and methods

Assessing Chlorurus spilurus diet

We quantified feeding behavior of C. spilurus in situ to examine how frequently this species preys on live coral. A diver conducted 20-min timed follows at two backreef locations in the lagoon of Mo’orea (French Polynesia) during the months of July – August in 2017 and 2018. We followed 23 focal individuals during peak grazing hours to control for temporal variation in foraging behavior (1000–1600 [68];). Targeted substrates were binned into bites either on live and dead corals, including rubble and pavement (turf algae, farmers turf and crustose coralline algae). We focused on C. spilurus > 150 mm as these individuals are most likely to bite live corals.

Sample collections and experimental design

Experiments were performed in July 2017 in Mo’orea, French Polynesia at the University of California Gump Research Station. For the mesocosm experiment, we collected 7 individual bullethead parrotfish (Chlorurus spilurus; Fig. 1a) (~ 200 mm total length) at ~ 3 m depth, in the back reef area along the north shore of Mo’orea (17°28′50.6″S 149°48′59.4″W) using hand and barrier nets. We also collected 10 apparently healthy Porites lobata colonies (~ 20 cm diameter) at the same depth and location. Fish and corals were immediately transferred to the Gump Research Marine Station. Once at the research station, parrotfish and coral colonies were placed into two independent mesocosms of 1155 L volume (Pentair AES Polyethylene Tank) that were supplied with flow-through seawater originating from the reef adjacent to the station. Shade cloth was applied uniformly on each mesocosm to moderate light intensity and promote acclimation of corals and fish. The day following collections, we sampled microbes from the mouth of each C. spilurus individual by carefully swabbing the inner side of the beak with sterile culture swabs (BD CultureSwab, BD). After sampling the mouth microbiome of all 7 fish, each fish was assigned to a single P. lobata colony and was coaxed into biting the colony two times at separate locations to create two independent feeding wounds. For microbial analyses of predated coral colonies, a sample of coral tissue was collected from one of the two wounds immediately after parrotfish had bitten the colony (Ti). Each colony’s other bite wound was sampled 48 h later (Tf). Samples were collected using a sterile bone cutter to remove a portion of the coral tissue (tissue, mucus and a small part of the skeleton) approximately 1 × 1 cm at the bite location. Rather than investigating how artificial wounding would affect coral microbiomes [39], we were specifically interested in understanding how parrotfish corallivory changed coral microbiomes. Therefore, instead of comparing the microbiome of corals with parrotfish predation to corals without parrotfish bites, we used mechanically wounded corals as controls. To compare how parrotfish vs. mechanical wounds affected coral microbiomes, the 3 remaining colonies of P. lobata (hereafter “mechanically wounded”) were wounded in two separate locations using sterile bone cutters to mimic the wounds caused by parrotfish. These mechanically wounded corals were then placed in a separate tank and assessed in the same manner as described above. The resulting lesions from both parrotfish and artificial wounding were ~ 2 mm deep and 1 cm long and removed coral tissue as well as part of the skeleton. Due to logistical limitations at the station, all of the parrotfish wounded corals were placed in one mesocosm, while all of the mechanically wounded corals were placed in a second mesocosm. Both mesocosms were fed flow though seawater via a common seawater source system and were physically adjacent to each other assuring similar light levels.

For our field study, we selected a 500 m stretch of shallow back reef area on the north shore of Mo’orea (17°28′35.2″S 149°47′34.9″W). Ten C. spilurus individuals (~ 20 cm total length) were collected using the same methods as described above. On the boat, individual parrotfish were transferred into large coolers equipped with bubblers and filled with seawater from the reef area. Microbial samples from each parrotfish mouth were directly collected using sterile swabs (BD CultureSwab, BD), as described above. Parrotfish were then released back onto the reef. The same day, we haphazardly selected Porites lobata colonies (20 cm diameter) in the back reef area with significant evidence of recent predation by parrotfishes (hereafter “bitten”), as well as colonies with no evidence of predation (hereafter “unbitten”; n = 10 colonies per wound status; Fig. 1b, c). We sampled a segment of coral tissue (1 × 1 cm) from the surface of each colony in situ using a sterile bone cutter. Bitten corals were sampled at a bite location chosen haphazardly on the colony, while unbitten corals were sampled at a haphazard location on the colony. Tissue collection was performed in situ across the designated reef an approximately 500 m stretch of lagoon. One liter water samples (n = 4) were collected haphazardly across the reef and filtered on a 0.2 μm Millipore filter. Sterile swabs and both coral and water samples were placed on ice in coolers until reaching the station where they were frozen at − 80 °C prior to microbial analyses.

DNA extraction and 16S rRNA gene amplification

High-throughput sequencing of the 16S rRNA gene was used to compare microbiome diversity, composition, and stability metrics. DNA extraction was performed using the DNeasy PowerSoil Kit (Qiagen) according to manufacturer instructions. AccuStart II PCR ToughMix (Quanta BioSciences, Gaithersburg, Maryland, USA) was used to perform two-step Polymerase Chain Reaction (PCR) on the V4 hypervariable region of the 16S rRNA gene using the primer pair 515FY (5′-GTGYCAGCMGCCGCGGTAA-3′) [69] and 806RB (5′-GGACTACNVGGGTWTCTAAT-3′) [70] targeting bacterial and archaeal communities. For each reaction, 6.25 μl AccuStart II ToughMix (2X), 1.25 μl forward primer (10 μM), 1.25 μl reverse primer (10 μM), 0.5 μl sample DNA, and 3.25 μl PCR-grade water was used. PCR amplification consisted of a 3 min denaturation at 94 °C followed by 35 cycles of 45 s at 94 °C, 60 s at 50 °C, and 90 s at 72 °C, and ending with 10 min extension step at 72 °C. A 1.5% agarose gel was run with amplified products which were manually excised to purify the 16S target band using Wizard® SV Gel and PCR Clean-Up System (Promega). The resulting products were then custom barcoded in a second PCR reaction with 12.5 μl ToughMix (2X), 9.5 μl water, and 1 μl of gel-purified sample DNA. The 12-cycle barcoding reaction consisted of a 5 min denaturation at 95 °C, 30 s melting at 95 °C, 3 min annealing at 63 °C, 30 s extension at 72 °C, ending with a 10 min hold at 72 °C. Barcoded amplicons were pooled in equivolume ratios and purified using Agencourt® AMPure XP beads. Prepared library pools were sequenced at the Center for Genome Research and Biocomputing (CGRB) at Oregon State University (OSU) on the Illumina MiSeq platform using MiSeq reagent kit v.3 (2 × 300 bp paired-end reads).

Data processing of mesocosm experiment samples

In association with the mesocosm experiment, a total of 50 samples were run through the data processing pipeline, 7 fish samples, 40 coral samples, 2 negative samples, and a positive control. Using VSEARCH v2.8.1 [71], sequences were truncated at the first position having a quality score less than or equal to 10, paired-end reads were merged, and merged reads with a total expected error > 1 per base or with > 1 N were discarded. This resulted in a total of 502,502 reads. Next, the Deblur workflow was used to trim quality-controlled reads to 250 base pairs, to identify exact sequences with single-nucleotide resolution, and to filter de novo chimeras [22]. This process resulted in 42 samples with 179,293 reads after 8 samples were lost in the Deblur workflow.

Next, the QIIME2 pipeline (https://qiime2.org [72];) was then used to process the OTU table resulting from the Deblur workflow. Taxonomy was assigned against the GreenGenes database [73], which is commonly used in microbial analyses [74], using classify-sklearn algorithm in QIIME2. Unassigned OTUs, singletons, and mitochondria or chloroplast sequences were removed from the OTU table. This removed a total of 7149 reads from the dataset. The number of sequences per sample type following filtering varied from 829 to 10,284 for coral tissue and from 3440 to 14,020 for fish mouthpart samples. Samples were then rarefied to a depth of 829 reads which resulted in the loss of 12 samples with insufficient read depth. The pre-filtered unprocessed sOTU table, metadata and associated negative control taxonomy table can be found in the Additional file 1: Tables S24 and S25. Rarefaction was performed using the function rarefy_even depth in the package phyloseq (v.1.26.1) in R. Four more samples were removed from the dataset as they were collected during a sampling timepoint that was ultimately excluded from these analyses due to low replication following the processing described above.

Data processing of field survey

A total of 139 samples were collected in the field and run through the data processing pipeline, including coral tissue, fish mouthparts, fish feces, sediment, water, and negative controls. Only a subset of these samples, 34, were relevant to the questions of this study and included in the analyses.

The following steps represent read counts for the 34 samples included in this analysis. Raw sequences were first demultiplexed then trimmed of primers and adapters resulting in 1,323,828 reads across the 34 samples. Using VSEARCH v2.8.1 [71], sequences were truncated at the first position having a quality score less than or equal to 10, paired-end reads were merged, and merged reads with a total expected error > 1 per base or with > 1 N were discarded. This resulted in a total of 526,544 reads. Next, the Deblur workflow was used to trim quality-controlled reads to 250 base pairs, to identify exact sequences with single-nucleotide resolution, and to filter de novo chimeras [22]. This process resulted in 33 samples with 164,793 after one sample was lost in the Deblur workflow.

Next, the QIIME2 pipeline (https://qiime2.org [72];) was then used to process the OTU table resulting from the Deblur workflow. Taxonomy was assigned against the GreenGenes database [73] using classify-sklearn algorithm in QIIME2. Unassigned OTUs, singletons, and mitochondria or chloroplast sequences were removed from the OTU table. This removed a total of 10,257 reads from the dataset. The number of sequences per sample type following filtering varied from 1551 to 7050 for coral tissue and from 2319 to 10,360 for fish mouthpart samples and 6021 to 8890 reads for water samples. Samples were then rarefied to a depth of 1551 reads which resulting in the loss of 4 samples with insufficient read depth. The pre-filtered unprocessed sub-operational taxonomic unit (sOTU) table, metadata and associated negative control taxonomy table can be found in the Additional file 1: Tables S26 and S27. Rarefaction was performed using the function rarefy_even depth in the package phyloseq (v.1.26.1) in R.

Data analyses and statistics

Following rarefaction, two alpha diversity metrics were computed – the observed species richness and the Shannon-Wiener index. The effects of sample type for i) mesocosm experiment (mechanically wounded, predated coral, fish mouth) within time periods and ii) the field (naturally unbitten, bitten coral, and fish mouth) experiment on diversity metrics were assessed using analysis of variance (ANOVA) with the function aov in the R package stats (v.3.5.3). When significant, pairwise comparisons among groups were performed using Tukey’s Honest significant differences (Tukey HSD). The assumptions of normality and homoscedasticity of residuals were tested using Shapiro-Wilk and Levene tests, respectively. If not fulfilled, nonparametric tests were performed using Kruskal-Wallis and Dunn’s tests using the function dunnTest within the R package FSA (v.0.8.22).

To illustrate the average relative abundance of the 25 most abundant taxa represented in each sample type at Tf (48 h) for both the manipulative and field experiments, we drew two heatmaps and at the family level by agglomerating the 25 taxa using the function tax_glom (including the command NArm = F) within the R package phyloseq (v1.26.1).

Furthermore, to display changes in microbial community composition among samples in either the mesocosm at Tf or the field experiments, two distinct non-metric multidimensional scaling (NMDS) on the Bray-Curtis dissimilarity matrices [75] were performed using the function metaMDS in the R package vegan (v.2.5–4). To test for differences in beta diversity among sample types for the manipulative and field experiments, we computed two permutational analysis of variance (PERMANOVA) based on the Bray-Curtis dissimilarity matrices and 999 permutations using the function Adonis in the R package vegan [76]. Subsequent pairwise differences were tested using the function pairwise.adonis in the R package vegan [76]. P-values were adjusted according to the false discovery rate, accounting for multiple comparisons.

We used the R package DESeq2 (v1.22.2) [77] on a pre-filtered unrarefied sOTU table to identify which sOTUs exhibited significant abundance among sample types both in the mesocosm experiment (at Ti and Tf) and in the field survey. From the pre-filtered unrarefied sOTU table, we used the function tax_glom within the phyloseq package (v1.26.1) to agglomerate taxa at the genus level (including the parameter NArm = F). DESeq2 incorporates a model based on the negative binomial distribution and includes a Wald posthoc test. P-values were adjusted for multiple comparisons using the Benjamini-Hochberg method [78].

Finally, to assess the variability of microbial composition between sampletype we computed an analysis of multivariate homogeneity of group dispersions [79]. This analysis tested whether community variability among samples, measured based on the Bray-Curtis dissimilarity metric, significantly differed between the three sample types, and was performed using the function betadisper in the R package vegan. Due to sample size limitation in the mesocosm experiment, the test was computed on the field survey coral samples only. When community variability significantly differed across sample type, pairwise tests were performed between groups using Tukey HSD. For clarity, findings presented in the result section are described as mean ± SE and p-values were considered significant for p < 0.05.

Availability of data and materials

Raw sequences reads associated to these datasets have been deposited in NCBI Short Read Archive (SRA) under the bioproject PRJNA573999.

The source code associated to this study is hosted in the Git repository: https://github.com/laylaeb/parrotfish.

References

  1. Rohwer F, Seguritan V, Azam F, Knowlton N. Diversity and distribution of coral-associated bacteria. Mar Ecol Prog Ser. 2002;243:1–10.

    Article  Google Scholar 

  2. Muscatine L, Porter JW. Reef corals: mutualistic symbioses adapted to nutrient-poor environments. Bioscience. 1977;27:454–60 American Institute of Biological Sciences Circulation, AIBS, 1313 Dolley Madison Blvd., Suite 402, McLean, VA 22101. USA.

    Article  Google Scholar 

  3. Ainsworth TD, Thurber RV, Gates RD. The future of coral reefs: a microbial perspective. Trends Ecol Evol. 2010;25:233–40 Elsevier.

    Article  PubMed  Google Scholar 

  4. Zhang Y-Y, Ling J, Yang Q-S, Wang Y-S, Sun C-C, Sun H-Y, et al. The diversity of coral associated bacteria and the environmental factors affect their community variation. Ecotoxicology. 2015;24:1467–77. https://doi.org/10.1007/s10646-015-1454-4.

    Article  CAS  PubMed  Google Scholar 

  5. van Oppen MJH, Blackall LL. Coral microbiome dynamics, functions and design in a changing world. Nat Rev Microbiol. 2019;17(9):557–67.

    Article  CAS  PubMed  Google Scholar 

  6. Bourne DG, Morrow KM, Webster NS. Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annu Rev Microbiol. 2016;70:317–40 Annual Reviews.

    Article  CAS  PubMed  Google Scholar 

  7. Rosenberg E, Koren O, Reshef L, Efrony R, Zilber-Rosenberg I. The role of microorganisms in coral health, disease and evolution. Nat Rev Microbiol. 2007;5:355. https://doi.org/10.1038/nrmicro1635 Nature Publishing Group.

    Article  CAS  PubMed  Google Scholar 

  8. Maynard J, Van Hooidonk R, Eakin CM, Puotinen M, Garren M, Williams G, et al. Projections of climate conditions that increase coral disease susceptibility and pathogen abundance and virulence. Nat Clim Chang. 2015;5:688 Nature Publishing Group.

    Article  Google Scholar 

  9. McDevitt-Irwin JM, Baum JK, Garren M, Vega Thurber RL. Responses of coral-associated bacterial communities to local and global stressors. Front Mar Sci. 2017;4:262 Frontiers.

    Article  Google Scholar 

  10. Harvell D, Jordán-Dahlgren E, Merkel S, Rosenberg E, Raymundo L, Smith G, et al. Coral disease, environmental drivers, and the balance between coral and microbial associates. Oceanography. 2007;20:172–95 Oceanography Society.

    Article  Google Scholar 

  11. Richardson LL. Coral diseases: what is really known? Trends Ecol Evol. 1998;13:438–43 Elsevier.

    Article  CAS  PubMed  Google Scholar 

  12. Shore A, Caldwell JM. Modes of coral disease transmission: how do diseases spread between individuals and among populations? Mar Biol. 2019;166:45. https://doi.org/10.1007/s00227-019-3490-8.

    Article  Google Scholar 

  13. Zaneveld JR, Burkepile DE, Shantz AA, Pritchard CE, McMinds R, Payet JP, et al. Overfishing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales. Nat Commun. 2016;7:11833 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Aeby GS, Santavy DL. Factors affecting susceptibility of the coral Montastraea faveolata to black-band disease. Mar Ecol Prog Ser. 2006;318:103–10.

    Article  Google Scholar 

  15. Chong-Seng KM, Cole AJ, Pratchett MS, Willis BL. Selective feeding by coral reef fishes on coral lesions associated with brown band and black band disease. Coral Reefs. 2011;30:473–81 Springer.

    Article  Google Scholar 

  16. Nicolet KJ, Chong-Seng KM, Pratchett MS, Willis BL, Hoogenboom MO. Predation scars may influence host susceptibility to pathogens: evaluating the role of corallivores as vectors of coral disease. Sci Rep. 2018;8:5258 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Raymundo LJ, Halford AR, Maypa AP, Kerr AM. Functionally diverse reef-fish communities ameliorate coral disease. Proc Natl Acad Sci. 2009. https://doi.org/10.1073/pnas.0900365106 National Acad Sciences.

    Article  CAS  Google Scholar 

  18. Martin SB, Sasal P, Cutmore SC, Ward S, Aeby GS, Cribb TH. Intermediate host switches drive diversification among the largest trematode family: evidence from the Polypipapiliotrematinae n. subf.(Opecoelidae), parasites transmitted to butterflyfishes via predation of coral polyps. Int J Parasitol. 2018;48:1107–26 Elsevier.

    Article  PubMed  Google Scholar 

  19. Hoey AS, Bonaldo RM. Biology of parrotfishes. Boca Raton: CRC Press; 2018.

  20. Bellwood DR. Direct estimate of bioerosion by two parrotfish species, Chlorurus gibbus and C. sordidus, on the Great Barrier Reef, Australia. Mar Biol. 1995;121:419–29 Springer.

    Article  Google Scholar 

  21. Bonaldo RM, Bellwood DR. Parrotfish predation on massive Porites on the Great Barrier Reef. Coral Reefs. 2011;30:259–69. https://doi.org/10.1007/s00338-010-0669-3.

    Article  Google Scholar 

  22. Amir A, McDonald D, Navas-Molina JA, Kopylova E, Morton JT, Zech Xu Z, et al. Deblur rapidly resolves single-nucleotide community sequence patterns. In: Gilbert JA, editor. mSystems, vol. 2; 2017. Available from: http://msystems.asm.org/content/2/2/e00191-16.abstract.

    Google Scholar 

  23. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet C, Al-Ghalith GA, et al. QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science. PeerJ. 2018; Preprints.

  24. Team RC. R: a language and environment for statistical computing; 2015; 2018.

    Google Scholar 

  25. Herder F, Freyhof J. Resource partitioning in a tropical stream fish assemblage. J Fish Biol. 2006;69:571–89 Wiley Online Library.

    Article  Google Scholar 

  26. Meirelles PM, Soares AC, Oliveira L, Leomil L, Appolinario LR, Francini-Filho RB, et al. Metagenomics of coral reefs under phase shift and high hydrodynamics. Front Microbiol. 2018:2203 Available from: https://www.frontiersin.org/article/10.3389/fmicb.2018.02203.

  27. Sunagawa S, Woodley CM, Medina M. Threatened corals provide underexplored microbial habitats. PLoS One. 2010;5:e9554. https://doi.org/10.1371/journal.pone.0009554 Public Library of Science.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Sneed JM, Ritson-Williams R, Paul VJ. Crustose coralline algal species host distinct bacterial assemblages on their surfaces. ISME J. 2015;9:2527 Nature Publishing Group.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Sekar R, Kaczmarsky LT, Richardson LL. Microbial community composition of black band disease on the coral host Siderastrea siderea from three regions of the wider Caribbean. Mar Ecol Prog Ser. 2008;362:85–98.

    Article  CAS  Google Scholar 

  30. Ramsby BD, Hoogenboom MO, Whalan S, Webster NS. Elevated seawater temperature disrupts the microbiome of an ecologically important bioeroding sponge. Mol Ecol. 2018;27:2124–37 Wiley Online Library.

    Article  CAS  PubMed  Google Scholar 

  31. Röthig T, Ochsenkühn MA, Roik A, van der Merwe R, Voolstra CR. Long-term salinity tolerance is accompanied by major restructuring of the coral bacterial microbiome. Mol Ecol. 2016;25:1308–23. https://doi.org/10.1111/mec.13567 John Wiley & Sons, Ltd (10.1111).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Ziegler M, Roik A, Porter A, Zubier K, Mudarris MS, Ormond R, et al. Coral microbial community dynamics in response to anthropogenic impacts near a major city in the central Red Sea. Mar Pollut Bull. 2016;105:629–40 Elsevier.

    Article  CAS  PubMed  Google Scholar 

  33. Lee S, Davy SK, Tang S-L, Kench PS. Mucus sugar content shapes the bacterial community structure in thermally stressed Acropora muricata. Front Microbiol. 2016;7:371 Frontiers.

    PubMed  PubMed Central  Google Scholar 

  34. Tout J, Siboni N, Messer LF, Garren M, Stocker R, Webster NS, et al. Increased seawater temperature increases the abundance and alters the structure of natural Vibrio populations associated with the coral Pocillopora damicornis. Front Microbiol. 2015;6:432 Frontiers.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Meron D, Atias E, Kruh LI, Elifantz H, Minz D, Fine M, et al. The impact of reduced pH on the microbial community of the coral Acropora eurystoma. ISME J. 2011;5:51 Nature Publishing Group.

    Article  PubMed  Google Scholar 

  36. Barott KL, Rodriguez-Mueller B, Youle M, Marhaver KL, Vermeij MJA, Smith JE, et al. Microbial to reef scale interactions between the reef-building coral Montastraea annularis and benthic algae. Proc R Soc B Biol Sci. 2011;279:1655–64 The Royal Society.

    Article  Google Scholar 

  37. Thurber RV, Burkepile DE, Correa AMS, Thurber AR, Shantz AA, Welsh R, et al. Macroalgae decrease growth and alter microbial community structure of the reef-building coral, Porites astreoides. PLoS One. 2012;7:e44246 Public Library of Science.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Bettarel Y, Halary S, Auguet J-C, Mai TC, Bui N, Bouvier T, et al. Corallivory and the microbial debacle in two branching scleractinians. ISME J. 2018;12:1109 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Maher RL, Rice MM, McMinds R, Burkepile DE, Vega Thurber R. Multiple stressors interact primarily through antagonism to drive changes in the coral microbiome. Sci Rep. 2019;9:6834. https://doi.org/10.1038/s41598-019-43274-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. van de Water JAJM, Leggat W, Bourne DG, van Oppen MJH, Willis BL, Ainsworth TD. Elevated seawater temperatures have a limited impact on the coral immune response following physical damage. Hydrobiologia. 2015;759:201–14 Springer.

    Article  CAS  Google Scholar 

  41. Apprill A, Hughen K, Mincer T. Major similarities in the bacterial communities associated with lesioned and healthy F ungiidae corals. Environ Microbiol. 2013;15:2063–72 Wiley Online Library.

    Article  CAS  PubMed  Google Scholar 

  42. Meron D, Atias E, Iasur Kruh L, Elifantz H, Minz D, Fine M, et al. The impact of reduced pH on the microbial community of the coral Acropora eurystoma. ISME J. 2010;5:51. https://doi.org/10.1038/ismej.2010.102 International Society for Microbial Ecology.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Huggett MJ, Apprill A. Coral microbiome database: integration of sequences reveals high diversity and relatedness of coral-associated microbes. Environ Microbiol Rep. 2019;11:372–85. https://doi.org/10.1111/1758-2229.12686 John Wiley & Sons, Ltd (10.1111).

    Article  PubMed  Google Scholar 

  44. Raina J-B, Tapiolas D, Motti CA, Foret S, Seemann T, Tebben J, et al. Isolation of an antimicrobial compound produced by bacteria associated with reef-building corals. Rodriguez-Lanetty M, editor. PeerJ. 2016;4:e2275. https://doi.org/10.7717/peerj.2275.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Tinta T, Kogovšek T, Klun K, Malej A, Herndl JG, Turk V. Jellyfish-associated microbiome in the marine environment: exploring its biotechnological potential. Mar Drugs. 2019;17(2):94.

    Article  CAS  PubMed Central  Google Scholar 

  46. D’Alvise PW, Lillebø S, Wergeland HI, Gram L, Bergh Ø. Protection of cod larvae from vibriosis by Phaeobacter spp.: a comparison of strains and introduction times. Aquaculture. 2013;384:82–6 Elsevier.

    Article  Google Scholar 

  47. Robert BBJGLB. Production of antibacterial compounds and biofilm formation by Roseobacter species are influenced by culture conditions. Appl Environ Microbiol. 2007.

  48. de Castro AP, Araújo SD, Reis AMM, Moura RL, Francini-Filho RB, Pappas G, et al. Bacterial community associated with healthy and diseased reef coral Mussismilia hispida from eastern Brazil. Microb Ecol. 2010;59:658–67. https://doi.org/10.1007/s00248-010-9646-1.

    Article  PubMed  Google Scholar 

  49. Mouchka ME, Hewson I, Harvell CD. Coral-associated bacterial assemblages: current knowledge and the potential for climate-driven impacts. Integr Comp Biol. 2010;50:662–74 Oxford University Press.

    Article  PubMed  Google Scholar 

  50. Tyagi A, Singh B. Microbial diversity in Rohu fish gut and inland saline aquaculture sediment and variations associated with next-generation sequencing of 16S rRNA gene. J Fish Life Sci. 2017;2:1–8.

    Google Scholar 

  51. Certner RH, Dwyer AM, Patterson MR, Vollmer SV. Zooplankton as a potential vector for white band disease transmission in the endangered coral, Acropora cervicornis. PeerJ. 2017;5:e3502 PeerJ Inc.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Pootakham W, Mhuantong W, Putchim L, Yoocha T, Sonthirod C, Kongkachana W, et al. Dynamics of coral-associated microbiomes during a thermal bleaching event. Microbiologyopen. 2018;7:e00604 Wiley Online Library.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Ziegler M, Seneca FO, Yum LK, Palumbi SR, Voolstra CR. Bacterial community dynamics are linked to patterns of coral heat tolerance. Nat Commun. 2017;8:14213. https://doi.org/10.1038/ncomms14213 The Author(s).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Pollock FJ, McMinds R, Smith S, Bourne DG, Willis BL, Medina M, et al. Coral-associated bacteria demonstrate phylosymbiosis and cophylogeny. Nat Commun. 2018;9:4921 Nature Publishing Group.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Rice MM, Maher RL, Thurber RV, Burkepile DE. Different nitrogen sources speed recovery from corallivory and uniquely alter the microbiome of a reef-building coral. PeerJ. 2019;7:e8056 PeerJ Inc.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Sunagawa S, DeSantis TZ, Piceno YM, Brodie EL, DeSalvo MK, Voolstra CR, et al. Bacterial diversity and white plague disease-associated community changes in the Caribbean coral Montastraea faveolata. ISME J. 2009;3:512 Nature Publishing Group.

    Article  CAS  PubMed  Google Scholar 

  57. Roder C, Arif C, Bayer T, Aranda M, Daniels C, Shibl A, et al. Bacterial profiling of white plague disease in a comparative coral species framework. ISME J. 2014;8:31 Nature Publishing Group.

    Article  CAS  PubMed  Google Scholar 

  58. Sweet M, Bythell J. Ciliate and bacterial communities associated with white syndrome and brown band disease in reef-building corals. Environ Microbiol. 2012;14:2184–99 Wiley Online Library.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Ng JCY, Chan Y, Tun HM, Leung FCC, Shin PKS, Chiu JMY. 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. 2015:1142 Available from: https://www.frontiersin.org/article/10.3389/fmicb.2015.01142.

  60. Séré M, Wilkinson DA, Schleyer MH, Chabanet P, Quod J-P, Tortosa P. Characterisation of an atypical manifestation of black band disease on Porites lutea in the Western Indian Ocean. PeerJ. 2016;4:e2073 PeerJ Inc.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Sweet M, Bythell J. White syndrome in Acropora muricata: nonspecific bacterial infection and ciliate histophagy. Mol Ecol. 2015;24:1150–9 Available from: https://doi.org/10.1111/mec.13097. John Wiley & Sons, Ltd (10.1111).

    Article  PubMed  PubMed Central  Google Scholar 

  62. Gignoux-Wolfsohn SA, Marks CJ, Vollmer SV. White band disease transmission in the threatened coral, Acropora cervicornis. Sci Rep. 2012;2:804 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Miura N, Motone K, Takagi T, Aburaya S, Watanabe S, Aoki W, et al. Ruegeria sp. strains isolated from the reef-building coral Galaxea fascicularis inhibit growth of the temperature-dependent pathogen Vibrio coralliilyticus. Mar Biotechnol. 2019;21:1–8. https://doi.org/10.1007/s10126-018-9853-1.

    Article  CAS  Google Scholar 

  64. Zaneveld JR, McMinds R, Thurber RV. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol. 2017;2:17121 Nature Publishing Group.

    Article  CAS  PubMed  Google Scholar 

  65. Nicolet KJ, Hoogenboom MO, Gardiner NM, Pratchett MS, Willis BL. The corallivorous invertebrate Drupella aids in transmission of brown band disease on the Great Barrier Reef. Coral Reefs. 2013;32:585–95 Springer.

    Article  Google Scholar 

  66. Sussman M, Loya Y, Fine M, Rosenberg E. The marine fireworm Hermodice carunculata is a winter reservoir and spring-summer vector for the coral-bleaching pathogen Vibrio shiloi. Environ Microbiol. 2003;5:250–5 Wiley Online Library.

    Article  PubMed  Google Scholar 

  67. Nicolet K, Hoogenboom M, Pratchett M, Willis B. Selective feeding by corallivorous fishes neither promotes nor reduces progression rates of black band disease. Mar Ecol Prog Ser. 2018;594:95–106 Available from: https://www.int-res.com/abstracts/meps/v594/p95-106/.

    Article  CAS  Google Scholar 

  68. Bruggemann JH, Van Oppen MJH, Breeman AM. Foraging by the stoplight parrotfish Sparisoma viride. I. Food selection in different, socially determined habitats. Mar Ecol Ser. 1994;106:41 INTER RESEARCH.

    Article  Google Scholar 

  69. Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. 2016;18:1403–14 Wiley Online Library.

    Article  CAS  PubMed  Google Scholar 

  70. Apprill A, McNally S, Parsons R, Weber L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat Microb Ecol. 2015;75:129–37.

    Article  Google Scholar 

  71. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. Hrbek T, editor. PeerJ. 2016;4:e2584. https://doi.org/10.7717/peerj.2584.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610 Nature Publishing Group.

    Article  CAS  PubMed  Google Scholar 

  74. Knight R, Vrbanac A, Taylor BC, Aksenov A, Callewaert C, Debelius J, et al. Best practices for analysing microbiomes. Nat Rev Microbiol. 2018;16:410–22. https://doi.org/10.1038/s41579-018-0029-9.

    Article  CAS  PubMed  Google Scholar 

  75. Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27:325–49 Wiley Online Library.

    Article  Google Scholar 

  76. Martinez Arbizu P. pairwiseAdonis: pairwise multilevel comparison using adonis. R Packag version 00; 2017. p. 1.

    Google Scholar 

  77. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550 BioMed Central.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  78. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57(1):289–300.

    Google Scholar 

  79. Anderson MJ. Distance-based tests for homogeneity of multivariate dispersions. Biometrics. 2006;62:245–53 Wiley Online Library.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

Many thanks to Andrew Shantz, Mark Ladd and Kelly Speare for help with fish collections, as well as Andrew Burns and Deanna Beatty for useful advice in the manuscript conception.

Funding

This work was supported by the Swiss National Science Foundation (Fellowships #P2SKP3_171764 - #P400PB_183867) to LE, a National Sciences Foundation (NSF) CAREER grant (OCE-1547952) to DEB, Dimensions of Biodiversity NSF grant (#1442306) to RVT, and a NSF Graduate Fellowship to RLM (#1314109-DGE). CSC was supported by the Teasley Endowment to the Georgia Institute of Technology and TL by the NASA (Grant #NNX14AR62A), the BOEM (MC15AC00006) and by NOAA.

Author information

Authors and Affiliations

Authors

Contributions

Conceived the experiment: LE, DEB, RVT; Collected the samples: KSM, KML, LE; Performed experiments: LE, KSM, KML; Performed laboratory work: LE, ERS; Performed data analysis: LE, TL, RLM, CSC; Wrote the paper: LE, CSC, RVT, DEB; Revised the manuscript: LE, CSC, RVT, DEB. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Leïla Ezzat.

Ethics declarations

Ethics approval and consent to participate

Research was completed under permits issued by the French Polynesian Government (Délégation à la Recherche) and the Haut-commissariat de la République en Polynésie Francaise (DTRT) (Protocole d’Accueil 2005–2018). This work represents a contribution of the Mo’orea Coral Reef (MCR) LTER Site (NSF OCE 16–37396) and was performed according to the UC Animal Care and Use Committee (IACUC #196, 2016–2019).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Results of the permutational ANOVA on the bacterial assemblages according to the sample type (control, predated coral and fish mouth) assessed at Ti and Tf for the mesocosm experiment. Table S2. Results of pair-wise tests on the effect of the sample type on the bacterial assemblages for the mesocosm experiment at Ti and Tf. Table S3. Average relative abundance of the families present in the fish mouths for the mesocosm experiment. Table S4. Average relative abundance of the families present in mechanically wounded corals at Ti for the mesocosm experiment. Table S5. Average relative abundance of the families present in mechanically wounded corals at Tf. Table S6. Average relative abundance of the families present in the predated corals at Ti for the mesocosm experiment. Table S7. Average relative abundance of families present in predated corals at Tf. Table S8. Results of the permutational ANOVA on the bacterial assemblages according to the sample type assessed for field experiment including or not water samples. Table S9. Results of pair-wise tests on the effect of the sample type on the bacterial assemblages for the field experiment. Table S10. Average relative abundance of the families present in the fish mouths for the field experiment. Table S11. Average relative abundance of the families present in naturally unbitten corals in the field. Table S12. Average relative abundance of families present in bitten corals for the field experiment. Table S13. Results of ANOVA and non-parametric tests of the effect of the sample type on alpha diversity metrics (Observed Richness and Shannon-Wiener Index) for the mesocosm experiment at Ti and Tf. Table S14. Results of posthoc tests assessing the effect of the sample type on alpha diversity metrics (Observed Richness and Shannon-Wiener index) for the mesocosm experiment at Ti and Tf. Table S15. Results of ANOVA and non-parametric tests on the effect of the type of sample on alpha diversity metrics (Observed Richness and Shannon-Wiener index) for the field experiment. Table S16. Results of posthoc tests on the effect of the sample type on alpha diversity metrics (Observed Richness and Shannon-Wiener index) for the field experiment. Table S17. Results from differential abundance analyses (DESeq2) on the effect of the sample type at Ti for the mesocosm experiment. Table S18. Average relative abundance of taxa present in greater differential abundance in predated corals compared to mechanically wounded corals for the mesocosm experiment at Ti and Tf. Table S19. Results from differential abundance analyses (DESeq2) on the effect of the sample type at Tf for the mesocosm experiment. Table S20. Differential abundance analysis for the field experiment according to the sampletype. Table S21. Average relative abundance of taxa present in greater differential abundance in naturally bitten corals compared to controls for the field experiment. Table S22. Results of Permutation test for homogeneity of multivariate dispersions (betadisper) on the effect of the sample type in the field survey. Table S23. Results of Permutation test for homogeneity of multivariate dispersion (betadisper) on the effect of the sample type in the field survey. Table S24. filtered unprocessed sOTU table for the mesocosm experiment. Table S25. Taxa table for the negative control of the mesocosm experiment. Table S26. filtered unprocessed sOTU table for the field survey. Table S27. Taxa table for the negative control in the field survey.

Additional file 2: Figure S1.

Relative abundance (Phyla) of taxa present in mechanically wounded, predated corals and fish mouths for the manipulative experiment at T48h. Figure S2. Relative abundance (Phyla) of taxa present in unbitten, bitten corals and fish mouths in the field. Figure S3. NMDS displaying the microbial assemblages according to the sample type at T48h for the mesocosm experiment. Figure S4. NMDS displaying the microbial assemblages according to the sample type for the field experiment.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ezzat, L., Lamy, T., Maher, R.L. et al. Parrotfish predation drives distinct microbial communities in reef-building corals. anim microbiome 2, 5 (2020). https://doi.org/10.1186/s42523-020-0024-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42523-020-0024-0

Keywords