Full Text

Turn on search term navigation

ARTICLE

Received 3 Nov 2015 | Accepted 4 May 2016 | Published 7 Jun 2016

DOI: 10.1038/ncomms11833 OPEN

Overshing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales

Jesse R. Zaneveld1,*, Deron E. Burkepile2,3,*, Andrew A. Shantz2, Catharine E. Pritchard2,4, Ryan McMinds1, Jrme P. Payet1, Rory Welsh1, Adrienne M.S. Correa1,5, Nathan P. Lemoine2, Stephanie Rosales1, Corinne Fuchs6, Jeffrey A. Maynard7,8 & Rebecca Vega Thurber1,*

Losses of corals worldwide emphasize the need to understand what drives reef decline. Stressors such as overshing and nutrient pollution may reduce resilience of coral reefs by increasing coralalgal competition and reducing coral recruitment, growth and survivorship. Such effects may themselves develop via several mechanisms, including disruption of coral microbiomes. Here we report the results of a 3-year eld experiment simulating overshing and nutrient pollution. These stressors increase turf and macroalgal cover, destabilizing microbiomes, elevating putative pathogen loads, increasing disease more than twofold and increasing mortality up to eightfold. Above-average temperatures exacerbate these effects, further disrupting microbiomes of unhealthy corals and concentrating 80% of mortality in the warmest seasons. Surprisingly, nutrients also increase bacterial opportunism and mortality in corals bitten by parrotsh, turning normal trophic interactions deadly for corals. Thus, overshing and nutrient pollution impact reefs down to microbial scales, killing corals by sensitizing them to predation, above-average temperatures and bacterial opportunism.

1 Department of Microbiology, Oregon State University, 226 Nash Hall, Corvallis, Oregon 97331, USA. 2 Department of Biological Sciences, Florida International University, 3000 NE 151st St, North Miami, Florida 33181, USA. 3 Department of Ecology, Evolution and Marine Biology, University of California, Santa Barbara, Santa Barbara, California 93106-9610, USA. 4 Department of Ecosystem Sciences, Penn State University, 235 Forest Resources Building, University Park, Pennsylvania 16802, USA. 5 Department of BioSciences, Rice University, 6100 Main Street, Houston, Texas 77005, USA. 6 Department of Biology, University of Florida, Gainesville, Florida 32611, USA. 7 SymbioSeas and Marine Applied Research Center, Wilmington, North Carolina, 28411, USA.

8 Laboratoire dExcellence )CORAIL* USR 3278 CNRS EPHE, CRIOBE, Papetoai, Polynsie Franaise. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to D.B. (email: mailto:[email protected]

Web End [email protected] ) or to R.V.T. (email: mailto:[email protected]

Web End [email protected] mailto:[email protected]

Web End =).

NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833

Tropical reefs continue to lose coral cover worldwide due to a variety of anthropogenic stressors, including warming ocean temperatures, nutrient pollution, sedimentation and

overshing13. Thermal stress from sustained periods of high temperature can kill corals by driving widespread coral bleaching events and coral disease outbreaks2,4,5. Localized stressors can also greatly impact corals. For example, elevated nutrients and overshing of herbivorous shes increase algal abundance on reefs68. Competition with these algae can compromise coral recruitment, growth and survivorship9,10. The intersecting impacts of these local stressors may ultimately compromise the resilience of reefs to disturbances, such as hurricanes and thermal anomalies1,11,12.

Algal competition and elevated temperatures may also impact reefs by disrupting corals symbiotic microorganisms, which serve functions ranging from nitrogen xation13 to pathogen inhibition14. Decades of coral microbiology research1326 demonstrate that coral microbiomes are sensitive to multiple stressors, including algal competition, elevated temperatures and disease20,2225. On the basis of this research, several qualitative and quantitative models have been proposed for the effects of algal competition or temperature on the coral microbiome. The DDAM hypothesis (dissolved organic carbon, disease, algae and microbes)27 suggests that turf and macroalgae secrete dissolved organic carbon28, which increases growth and oxygen consumption of bacteria29,30 (but see ref. 31), ultimately harming corals through hypoxia17,26,32. Two additional, non-mutually exclusive hypotheses are that algae produce allelochemicals that directly harm corals, with possible downstream effects on the microbiome20 and that interactions with certain algae lead to transfer of algae-associated microbes to corals19,33.

Sea-surface temperatures also strongly impact the coral microbiome15,16,24,34. While mass coral bleaching caused by high levels of thermal stress has received the greatest attention, even modest increases in temperature appear to make corals more vulnerable to opportunistic bacteria. Elevated temperatures increase the release of dimethylsulfoniopropionate from corals. Many opportunistic microorganisms chemotax along gradients of dimethylsulfoniopropionate, potentially allowing them to target thermally stressed corals34. Elevated temperatures also increase expression of virulence genes in some opportunistic bacteria35, alter innate immune gene expression in corals36 and inhibit the protective role of bacteria in the coral mucus37. These observations have prompted mathematical models of competition between coral mutualists and pathogens, suggesting that temperature variation mediates changes between pathogen- and mutualist-dominated stable states25.

Together, this research strongly suggests that human impacts such as overshing and nutrient pollution may interact with sea-surface temperatures to cause changes in coral reef benthic communities and their microbiomes that together contribute to coral mortality17,18,26. Much of the foundational work on how corals respond to global and local stressors has been done in short-term, small-scale lab and eld experiments based at the organismal level (for example, refs 16,17,19,24,26). Therefore, it is difcult to extrapolate these prior results to determine how factors that are thought to contribute to coral decline interact to alter coral benthic communities and their microbiomes over ecologically realistic temporal and spatial scales. Therefore, we designed a long-term eld experiment to examine how altering herbivory and nutrient pollution impacted coralalgalmicrobe interactions within the context of seasonal temperature variation. Specically, we sought to test the following main hypotheses over ecologically relevant spatial and temporal scales: (1) exclusion of herbivores and nutrient pollution both increase the abundance, and diversity of turf and macroalgae, leading to

intensied coralalgal competition and additive increases in coral tissue loss and mortality; (2) coralalgal competition interacts with above-average temperatures to shift coral microbiomes away from their normal conguration and towards distinct, pathogen-dominated stable states; and (3) alterations to the coral microbiome due to increased algal competition or above-average temperatures positively correlate with increased coral disease, and reduced coral growth and survivorship.

To test these hypotheses, we conducted a 3-year eld experiment (June 2009August 2012) that simulated overshing and nutrient pollution on a reef in the Florida Keys, USA. We tracked their impacts on the benthic community, coralmicrobe dynamics and coral survivorship across multiple seasons. To simulate nutrient pollution, four 9-m2 plots of reef benthos were enriched in nitrogen and phosphorus, while four control plots remained at ambient nutrient levels (Methods). Enrichment increased nitrogen and phosphorous concentrations approximately four- eightfold above ambient, similar to reefs impacted by nutrient pollution (Methods). We simulated overshing of herbivorous shes (for example, parrotshes and surgeonshes) by nesting two 1 m2 herbivore exclosures and two 1 m2 open-topped exclosure controls within each plot (Supplementary Fig. 1). This created factorial treatments of: (1) control,(2) herbivore exclusion, (3) enriched nutrients and (4) herbivore exclusion plus enriched nutrients (Supplementary Notes; Supplementary Data 1 sheet b).

ResultsOvershing and nutrients increase algal cover and richness. Herbivore exclusion rapidly increased algal cover up to sixfold, species richness up to threefold (Fig. 1a,b) and altered algal community composition (Supplementary Fig. 2; Supplementary Data 1 sheets ac). Excluding herbivores increased numerous macroalgae (for example, Dictyota, Halimeda and Amphiroa) and tall lamentous turf algae (X0.5 cm; Supplementary Fig. 2;

Supplementary Notes) at the expense of crustose coralline algae and closely cropped turf algae (o0.5 cm), which are benecial or neutral for coral recruitment, settlement and growth10. Nutrient pollution only slightly increased algal richness (Fig. 1b), instead nutrient pollution increased algal cover by driving seasonal blooms of tall lamentous turf algae, Dictyota and Lyngbya cyanobacteria. Overall, removing herbivorous shes, and to a lesser extent nutrient pollution, facilitated growth of algae known to increase coral tissue loss or mortality via shading, abrasion and allelopathy38,39.

Algae or elevated temperature alter coral microbiomes. To identify how algal communities and nutrient pollution affected the coral microbiome, we collected DNA samples from the surface mucus layer of 80 coral colonies (genera Porites, Siderastrea and Agaricia) at approximately monthly intervals. From these samples, 478 were used for 16S ribosomal RNA (rRNA) gene amplicon analyses following quality control (Supplementary Data 2 sheets ad). The phyla Proteobacteria, Cyanobacteria, Bacteroidetes and Actinobacteria were most abundant in all coral genera, with a core set of 13 bacterial orders in Z95% of all samples (Supplementary Data 3 sheet a). Order Synechococcales (phylum Cyanobacteria), a proposed coral mutualist13, was particularly abundant on control corals.

Increasing algal cover or elevated temperature suppressed the typical, Synechococcus-dominated microbiome of healthy corals and facilitated blooms of other microbes, including many putative opportunists or pathogens (Fig. 2; Supplementary Fig. 3). Alteromonadales, and 10 other orders, many of which were Proteobacteria, increased in abundance as upright algal cover

2 NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833 ARTICLE

a b

225

Exclosure & Nutrients

Nutrients

Control

Pretreatment

200

Nutrients P < 0.001Exclosure P < 0.001 Exclosure x nutrients P < 0.001

Upright algal cover (%)

Richness (algal species)

175

150

Exclosure

125

100

Pretreatment

Summer

Fall

Winter

Spring

Summer

Fall

Winter

Spring

Summer

Fall

Winter

Spring

Summer

Pretreatment

Summer

Fall

Winter

Spring

Summer

Fall

Winter

Spring

Summer

Fall

Winter

Spring

Summer

2009 2010 2011 2012

Nutrients P < 0.001Exclosure P < 0.001 Exclosure x nutrients P = 0.004

2009 2010 2011 2012

Figure 1 | Herbivore exclusion and nutrient pollution alter algal communities. (a) Upright algal cover (macroalgae, cyanobacteria and tall turf algae), with herbivore exclusion and/or nutrient pollution over time. Total cover often exceeds 100% due to the three-dimensional algal canopy. (b) Macroalgal species richness over time. P values are from mixed effect models (Supplementary Table Data 1c). Data are meanss.e.m.

(including tall turf algae, cyanobacteria and macroalgae) increased. Vibrionales and Oscillatoriales increased in abundance with increasing temperature (Fig. 2a; Supplementary Data 3 sheets bg).

Herbivore exclusion and ensuing coralalgal contact also increased the relative abundance of many otherwise rare microbial orders (Supplementary Table 1), increasing microbial community richness and evenness (Supplementary Data 3 sheet f). In contrast, nutrient pollution suppressed many taxa, allowed fewer bacterial orders to dominate (Supplementary Table 1), and decreased community evenness (Supplementary Data 3 sheet f).

Principle coordinates analysis (PCoA) identied transitions from Synechococcus- to Proteobacteria dominance, as the most important axis structuring coral microbiomes (First PCoA axis; Fig. 2a,b). Proteobacteria blooms were accompanied by a drop-off in community evenness (Fig. 2b) and represented 91.9% of bacteria in the least even quartile of samples (Supplementary Fig. 3c). For example, during periods of above-average temperatures, Vibrionales represented up to 62% of the microbial community on individual corals, rising from a mean background level of 5% on unstressed corals.

We used random forest analysis, a machine-learning method, to summarize the extent to which the dominant microbe on coral surfaces could be predicted from external conditions such as algal contact and sea-surface temperature. Random forest analysis predicted displacement of Synechococcus by blooms of other microbes with 78.5% accuracy. Temperature and upright algal cover (including both macroalgae and tall turf algae) were the most informative predictor variables (Supplementary Data 3 sheet h). Bacterial opportunists that displaced Synechococcus bloomed at different combinations of temperature and overall upright algal cover (Fig. 2c; see Supplementary Data 3 sheet f for turf algae and cyanobacteria specically). However, conditions favoring different opportunists overlapped, rendering the specic opportunist that displaced Synechococcales partially stochastic and not predictable from random forest analysis of ecological data (Supplementary Data 3 sheet h).

Antibiotic-producing bacteria may play an important role in protecting corals from outbreaks of harmful bacteria14. In our data, opportunistic Proteobacteria reached their highest abundance in corals, where Actinobacteria were below B2.5%

relative abundance (Supplementary Fig. 4), suggesting that antibiotic-secreting Actinobacteria may help to protect corals from opportunists. Actinobacteria, in turn, decreased with increasing algal cover or elevated temperatures, suggesting that

these stressors remove an important biotic barrier to potential pathogens (Supplementary Fig. 3; Supplementary Data 3 sheets c,d).

Increasing algal cover or elevated temperatures also shifted the predicted functional proles of corals towards microbial pathogenesis. We used the PICRUSt software package40 to estimate the functional consequences of changes in microbial community composition (Methods). According to these estimates, micro-biomes subject to algal competition or above-average temperatures became enriched in pathways involved in opportunism (for example, cell motility and secretion systems) and depleted in pathways for antibiotic production that may help healthy coral microbiomes resist invasion (Supplementary Fig. 5a,b; Supplementary Data 4 sheets a,b). Increasing upright algal cover was also correlated with predicted increases in the abundance of genes for utilization of glycans and pentose sugars in the coral microbiome (Supplementary Fig. 5a). Conversely, upright algal cover was associated with decreases in seven categories of metabolic pathways related to amino acid metabolism and four related to lipid and fatty acid metabolism (Supplementary Fig. 5a). These data support the idea that increased algal cover promotes growth of microbes on the coral surface that are capable of rapid microbial metabolism of algal sugars17,26,32,41.

Algae and temperature increase coral microbiome b-diversity. Stressful conditions may shift coral microbiomes from one stable state to another25. We had originally hypothesized that algal competition would produce such a shift in the structure of coral microbiomes. However, we did not nd any evidence for such treatment-induced shifts between alternative stable states (Supplementary Fig. 6a). Rather, we found that algal contact and increased ambient temperatures reduced the overall stability of the microbiome as a whole. This type of microbiome destabilization was manifest as increased sample-to-sample variability and quantied using measures of microbiome b-diversity (for example, ref. 42). While b-diversity was originally used to describe species turnover among habitats43, in microbial ecology host organisms are often treated as the relevant habitat. This allows b-diversity measures to be used to compare inter-individual variation in the microbiome. This notion has also been generalized to measure species turnover within individuals over time (intra-individual variation). However, few studies (notably ref. 42) have emphasized the relationship between the health of animal hosts and the variability of their microbiomes.

NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833

Color reflects the most abundant microbial order:

Size reflects microbial community evenness:most even

Least even

PC2 (9%)

Dominant bacterial ordersnot enriched by temperature or upright algal cover

Dominant bacterial orders enriched bytemperature

Vibrionales (18)

Alteromonadales (51)

Camplyobacterales (2)

NB1-j (1)

Synechococcales (255)

Bacillales (1)

Oscillatoriales (5)

Pirellulales (3)

Rhodobacterales (74)

Burkholderiales (1)

Pseudomonadales (5)

Rickettsiales (5)

Sphingomonadales (1)

Dominant bacterial orders enriched byupright algal cover

Flavobacterales (5)

Lentisphaerales (2)

Oceanospirallales (5)

Verrucomicrobiales (6)

Cytophagales (5)

Unclassified alphaproteobacteria (2)

Rhizobiales (3)

Thiotrichales (1)

Bacteroidales (3)

Chroococcales (2)

Brachyspirales (1)

PC3 (4%)

PC1 (27%)

180

0 25 26 27 28 29 30

1.0

160

Chroococcales

0.8

140

Bacteroidales

Cytophagales

Upright algal cover (%)

120

Oscillatoriales

Synechococcales

0.6

Evenness

Proteobacteria

100

Camplyobacterales

Alteromonadales

Verrucomicrobiales

0.4

Vibrionales

Pirellulales

Synechococcales

0.2

0.0

0.4 0.2 0.0 0.2 0.4

PC1 (weighted UniFrac)

Temperature (C)

Rhizobiales

Taxa enriched by elevated temperature

Taxa enriched by increased upright algal cover

Synechococcales

Figure 2 | Algal competition and temperature combine to alter coral microbiomes by driving bacterial blooms. (a) Principle coordinates analysis plot, summarizing weighted UniFrac distances between coral microbial samples (n 435). The main pattern we seek to show here is a shift from Synechococcus-

dominated communities (cyan) to dominance by a wide variety of other orders as one moves from left to right along PC1. Points are coloured to reect the most dominant (numerically abundant) microbial order in each sample. Sizes reect quartiles of microbial community evenness (measured by equitability). Orders that were signicantly enriched by temperature or upright algal cover (Pearson regression, FDR qo0.05, Supplementary Table Data 3f and 3g) are grouped together in the legend. Parentheses next to each order note the number of samples, in which that order was dominant. Orders that were dominant in Z15 samples are marked in bold. (b) Displacement of Synechococcales by varied Proteobacteria structured differences between stressed coral microbiomes. Points plot the 1st PC axis (a) against the relative abundance of Synechococcales (cyan squares, dot-dashed line), Proteobacteria (orange triangles, dashed line) and overall microbial community evenness (blue circles, dotted line). Lines show local regression (Loess regression, span 0.75),

with grey bars shading extending to twice the s.e. of the regression. (c) Combinations of algal cover (macroalgae, cyanobacteria and tall turf algae) and naturally occurring seawater temperature that allowed algae- or temperature- enriched bacterial orders (a) to dominate coral microbiomes. The position of each point, and its associated error bars, represents mean algal cover and temperature, and their s.e.s, for all samples that were dominated by the labelled bacterial order. Synechococcales, which dominated most (258/435) coral microbiome samples, are included as a reference. Taxa dominating Z15 samples are marked in bold. Coloured polygons enclose suites of taxa whose mean abundance increased with temperature (orange) or upright algal cover (green;

Supplementary Data 3 sheets fh). Vibrionales and Oscillatoriales responded to both temperature and algal cover, but for Vibrionales only temperature was signicant.

In our experiment, algal contact or elevated temperature drove blooms of microbial opportunists that destabilized coral microbiomes, increasing b-diversity between samples taken from corals subject to these stressors (Supplementary Figs 6ad and 7ab). This increased inter-individual and/or

temporal variation manifested as increased dispersion of samples from stressed hosts around a centroid of samples from healthy hosts, rather than discrete clusters of samples representing healthy versus diseased hosts (Supplementary Fig. 6).

4 NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833 ARTICLE

Algae-driven increases in b-diversity were largely consistent, regardless of whether corals were analysed together or separately by genus (although not individually signicant for Porites; Supplementary Notes). b-diversity increased in samples where

Synechococcus was displaced by Proteobacteria (Supplementary Data 3 sheet l), indicating that blooms of Proteobacteria contribute to microbiome instability. Signicant microbiome destabilization was also found upon reanalysis of a previous study19, in which Porites corals were placed in contact with macroalgae (Supplementary Fig. 7c).

Algal contact induces coral disease and mortality. In this long-term experiment, we found that radical changes to the coral microbiome were strongly correlated with coral tissue loss and mortality in the eld (Fig. 3). Although corals in control plots grew, gaining 36.8% of initial tissue area on average, other treatments caused 30.636.1% tissue loss (Supplementary Data 5 sheets ad) and six- ninefold increases in mortality (Fig. 3a). Importantly, microbiome evenness was lower and instability (b-diversity) higher (permutational t-tests, P 0.033 and

P 0.001, respectively) in corals that lost tissue by the end of the

experiment.

Direct algal contact increased coral microbiome instability (Supplementary Fig. 7a), subsequent tissue loss (Fig. 3b,

Supplementary Data 5) and mortality (Supplementary Fig. 8ac). Competition with Dictyota algae caused the most severe microbiome disruptions (Supplementary Fig. 7b) and the greatest overall coral mortality (Supplementary Fig. 8c). Further, as the diversity of algal competitors increased, so did microbiome instability, the prevalence of tissue loss and coral mortality (Fig. 3c). In Siderastrea corals, algal contact also increased by twofold the prevalence of dark spot syndrome (DSS), a coral disease with poorly understood aetiology (Fig. 3c, Supplementary Fig. 9a). DSS-affected Siderastrea had 5-fold greater prevalence of tissue loss and, on average, lost tissue, while Siderastrea without DSS had positive growth (Supplementary Fig. 9c,d).

Random Forest analysis predicted coral tissue loss with76.80.11% accuracy using a mixture of benthic and microbial data (Supplementary Data 3 sheet j). Because this analysis compared the susceptibility of colonies to tissue loss rather than the timing of tissue loss, temperature was not included as a variable. The abundance of macroalgae and tall turf algae, which increased coral morality (Supplementary Fig. 8c), were the most important predictors of tissue loss. The next best predictors were the abundance of Rhodobacterales and Rhodospirallales, two putatively opportunistic microbial orders, both in the Proteo-bacteria and both commonly enriched in diseased corals (for example ref. 44). In this study, Rhodospirallales were enriched by herbivore exclusion (Supplementary Table 1).

a b

75%

P = 0.006

No algal contact

Nutrients P = 0.001Exclosure P = 0.069Nutrients x exclosure P = 0.034

100%

0.8

Algal contact

90%

50%

B B

50%

80%

0.75

70%

Coral mortality

Tissue gain or loss

40%

0.7

25%

P = 0.164 P < 0.001 P = 0.037

60%

30%

Prevalence

50%

0.65

40%

20%

25%

30%

0.6

20%

10%

0.55

50%

10%

0.5

75%

Control

Exclosure

Nutrients

Exclosure and

nutrients

0 1 2 3

Microbial -diversity

Control

Nutrients

Exclosure

Exclosure and

nutrients

Macroalgal types contacted

Tissue loss Total colony mortality DSS (Siderastrea only) -diversity

d e

0.6

P = 0.001

0 21 22 23 24 25 26 27 28 29 30 31 32

Proteobacteria Cyanobacteria

100%

60%

No tissue loss Tissue loss

P = 0.2289

[afii9851] P < 0.001

0.58

0.9

90%

50%

Microbial -diversity

0.56

0.8

80%

Portion of total

coral mortality

Evenness -diversity

0.7

40%

0.54

70%

0.52

0.6

60%

30%

0.5

50%

P = 0.037

0.4

20%

0.48

40%

0.46

0.3

30%

10%

0.44

0.2

20%

0.1

10%

0.42

Winter

Spring

Summer

Fall

0.4

High ( 30 C)

Mid (2430 C)

Low(< 24 C)

Temperature

Evenness Proteobacteria

-diversity Cyanobacteria

Temperature

Figure 3 | Multiple stressors disrupt coral microbial communities and produce coral mortality. (a) Cumulative coral mortality at end of experiment. P values are from mixed effect models, letters over bars show differences in Tukeys post hoc tests. Herbivore removal signicantly increased coral mortality relative to controls (Tukeys post hoc test Po0.05), but not relative to nutrient pollution alone (post hoc test and mixed effects model P40.05). (b) Effects of algal contact on coral tissue area, across treatments. P values from ANOVAs test the effect of algal contact within each treatment. (c) Number of algal taxa contacting corals versus microbiome b-diversity (weighted UniFrac distance), and the prevalence of coral tissue loss, mortality and Siderastrea dark spot syndrome (DSS). (d) Seasonal distribution of coral mortality, coloured by treatment (a). Red line marks null expectation of equal mortality across seasons. P value is from a w2-test. (e) Microbial community b-diversity for corals with or without tissue loss, split by temperature. P values reect non-parametric t-tests of distances. (f) Temperature effects on coral microbial variability, evenness and relative abundance of Proteobacteria or Cyanobacteria.

Evenness and b-diversity data are meanss.e.m. Microbial and coral health data are averaged within each 1 C interval on the x axis. The vertical red line at 30 C indicates the point nearest to the MMM 1 C value for our site (30.26 C); temperatures beyond this result in accumulation of degree heating

weeks of coral thermal stress (Methods).

NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833

Together, these data show that the rise of algae following chronic nutrient enrichment and removal of consumer species disrupts coral microbiomes, facilitates disease, increases tissue loss and causes coral death.

Overshing or nutrients induce seasonal coral mortality. Surprisingly, coral mortality caused by simulated overshing and nutrient pollution was strongly seasonal. Although we observed minimal coral bleaching, B80% of coral mortality was concentrated in summer and fall (Fig. 3d), which include the warmest months of the year. Yet this seasonal mortality occurred only when herbivores were removed or nutrients increased, suggesting that local stressors and temperature interact to kill corals. If seasonal temperature variation and local stressors do interact to drive coral mortality, we might expect years in which corals suffer greater sub-bleaching thermal stress to correspond to greater coral mortality.

To put treatment outcomes in the context of seasonal temperature variation and coral thermal stress, we calculated raw sea-surface temperatures for our study site using both the Pathnder v5.245 and HYCOM Gulf of Mexico resources46 (Methods; Supplementary Fig. 10). From these we calculated several standard derived measures of coral thermal stress47. The maximum monthly mean (MMM) is the average temperature of the warmest month in long-term climate data for our site, excluding the study period (that is, 19822009). In predictions of coral bleaching, temperatures that exceed the MMM by Z1 C (hereafter MMM 1 C) are considered thermal stress47. The

level of thermal stress is then measured in degree heating weeks (DHWs), which accumulate weekly temperatures above the MMM 1 within a 3-month timeframe (Methods). Because

the threshold for increased pathogenesis may be lower than that for coral bleaching15, we also calculated DHWs based on the MMM (DHWMMM) for reference, although the more standard DHWMMM 1 C was used in all calculations.

Although summertime temperatures during our experiment exceeded the MMM 1 C, and positive DHW values

accumulatedindicating some thermal stressin no cases did these values rise to levels that would predict increased bleaching risk (Supplementary Data 6). However, we did observe correlations between sub-bleaching thermal stress and warm-season coral mortality during 2009, 2010, 2011 and 2012, despite limited

replication of seasons across years (that is, n 4). Across all

treatments, the extent of mortality in summer and fall was positively correlated with the extent of sub-bleaching thermal stress in DHWs (r 0.960, P 0.03), an effect strengthened when

control corals were excluded (r 0.993, P 0.007).

Differences in microbiome stability between healthy and unhealthy corals were exaggerated at temperatures above 30 C (Fig. 3e), as was the displacement of likely benecial Cyanobacteria by opportunistic Proteobacteria (Fig. 3f). Regressions using local temperature data (Supplementary Fig. 10b) indicated that both increasing abundance of Proteobacteria relative to Cyano-bacteria and reduced microbial community evenness became noticeable around the MMM (29.26 C) before becoming pronounced around the MMM 1 C (30.26 C). This nding

suggests that the local MMM may be a critical temperature

threshold for the onset of bacterial opportunism.

Temperature variation explained differences in microbial community structure over time better than other measured seasonal parameters (Supplementary Notes). The accumulation of thermal stress due to periods of above-average temperature also appeared to inuence microbial communities. Samples in which sub-bleaching thermal stress had accumulated (that is, DHWs40) showed increased b-diversity (permutational t-test on weighted UniFrac distance matrices; P 0.001) relative to

samples from periods without thermal stress (DHWs 0).

Cumulative thermal stress also had small, but signicant effects on microbial communities after accounting for daily temperature variation. As the difference in thermal stress between microbiome samples increased, the weighted UniFrac distances between microbial communities also increased (Mantel test, P 0.01

and r 0.176; Methods), even after accounting for the tempera

ture on the day the sample was taken (partial Mantel test, P 0.01 and r 0.17). These microbiological changes were

probably not due to confounding effects of coral bleaching, as cumulative thermal stress stayed well below the four DHW threshold for mass coral bleaching5 (Supplementary Fig. 10a; Supplementary Notes), and little visible bleaching was observed.

Nutrients alter the outcome of predation on corals. Coral mortality patterns showed that nutrient pollution changes the impact of important consumers on reefs. Parrotshes preyed on

P = 0.025

80%

P < 0.001

Coral mortality (Poritesonly)

Proteobacteria:Cyanobacteria ratio

(Porites only)

60%

40%

P = 0.744

20%

P = 0.990

0% 0%

62%

Control

No predation

Control

predation

Nutrients

No predation

Nutrients

predation

Control

No predation

Control

predation

Nutrients

No predation

Nutrients

predation

Figure 4 | Effects of nutrient pollution and parrotsh predation on coral mortality and microbiology. (a) Mortality after predation on Porites corals in control or nutrient pollution plots. P values reect Fishers exact test. (b) Effects of predation on the relative abundance of Proteobacteria and Cyanobacteria in Porites corals in control or nutrient pollution plots. Parrotsh predation columns (in red) reect samples taken after the rst evidence of parrotsh predation. Error bars in b reect the 95% CI of the ratio. P values are from non-parametric t-tests.

6 NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833 ARTICLE

63.4% of Porites corals in open plots, with a marginally higher predation prevalence in nutrient-enriched plots (P 0.08;

Supplementary Data 5 sheet f). However, the ultimate outcome of predation for coral health varied greatly under ambient versus nutrient-enriched conditions. Under ambient conditions, the impact of parrotsh predation was negligible, inducing tissue loss in only 7% of bitten Porites and never inducing mortality. However, in the presence of nutrient pollution, 92% of Porites bitten by parrotshes lost tissue and 62% eventually died (Fig. 4a). Because other stressors increased potentially opportunistic Proteobacteria compared with Cyanobacteria (primarily Synechococcales), we compared the ratio of these bacterial phyla in bitten corals. The combination of predation and nutrient pollution, but not predation alone, increased the ratio of Proteobacteria relative to Cyanobacteria (Fig. 4b).

DiscussionWe demonstrate that overshing and nutrient pollution alter benthic communities, causing coral microbiome disruption, blooms of opportunistic coral pathogens and long-term increases in coral disease, tissue loss and mortality. This experimental framework, which combines classical ecological eld experiments with microbial time series, allowed us to test key predictions of how coral reefs respond to human impacts over ecologically important time scales and with much greater microbiological detail than has previously been possible in eld experiments. The results changed our view of how coral microbiomes respond to environmental stressors; uncovered new interactions among corallivory, nutrient pollution and bacterial opportunism; and quantied how local stressors interact with seasonal temperature variation to impact coral microbiomes and, ultimately, coral survivorship. Together, these results show how altering important trophic interactions can fundamentally reorganize coral reefs down to microbial scales, with multiple negative consequences for reef health.

We predicted that simulated overshing or nutrient pollution would drive distinct benthic community changes and that these benthic changes would additively increase coral mortality. More specically, we predicted that exclusion of herbivores would lead to large increases in the abundance and diversity of macroalgae. In contrast, nutrient pollution would only modestly increase macroalgal cover and diversity, instead driving growth of nutrient-limited but unpalatable members of the benthic community such as Cyanobacteria. Similar to past work (for example, refs 68), we show that herbivore removal rapidly leads to over sixfold increases in algal abundance and threefold increases in algal diversity, as well as intense coralalgal competition. Many of the algal genera that increased, such as Sargassum, Dictyota, Amphiroa and Turbinaria, cause coral tissue loss or mortality via shading, abrasion and allelopathy38,39. Nutrient enrichment led to lower overall increases in algal abundance, but facilitated growth of certain taxa such as Cyanobacteria and lamentous turf algae, which can compete intensely with corals10. However, the consequences of these local stressors for coral mortality deviated greatly from our prior expectations. We expected herbivore removal and nutrient pollution to have additive effects on corals, with their combined effects causing much greater coral mortality than either individually. Yet, while all treatments signicantly increased coral tissue loss and coral mortality above controls, levels of mortality and tissue loss were similar under nutrient pollution, herbivore removal or combined treatments. This pattern likely resulted from different mechanisms of coral tissue loss and mortality operating in the different treatments.

In the nutrient pollution treatment, we traced coral decline to an unexpected interaction among corallivory, nutrient pollution,

bacterial opportunism and coral death. Parrotshes are key herbivores on reefs, but also prey on corals as part of their diet48. This corallivory is relatively intense in the Florida Keys, USA where populations of parrotshes are robust49. Here we found that Porites corals bitten by parrotsh lost tissue and died at much higher rates under nutrient-enriched conditions than in ambient conditions. In fact, parrotsh corallivory caused no mortality in ambient nutrient conditions. Microbial community shifts towards Proteobacteria on the surface of bitten and nutrient-enriched corals (but not bitten corals in ambient nutrient conditions) suggest that increased bacterial opportunism following wounding may cause this increased tissue loss and mortality. It is not clear whether the observed increase in Proteobacteria is due to a proliferation of pathogens vectored by parrotsh (as suggested for butteryshes50) or represents a more general nutrient-driven increase in susceptibility to infection following wounding. In either case, this unexpected observation is especially worrisome, as it shows that nutrient pollution turns parrotshes, which are normally thought of as coral allies, into agents of mortality for some corals. This nding is important given that nutrient pollution is a problem on many reefs and that restoration of parrotsh populations is an important goal as part of coral conservation and management51. Our data suggest that restoration of parrotshes without efforts to combat water quality issues could have surprisingly negative consequences.

When herbivorous shes were removed, coralalgal competition and disruption of the coral microbiome appeared to be the main drivers of coral mortality. We showed that algal contact increases coral microbiome richness, facilitates growth of many conditionally rare taxa and increases overall microbiome destabilization. This is broadly consistent with previous reports of algal alteration of coral microbiomes through mechanisms such as provision of algal neutral sugars41, secondary effects of alleotoxins on the coral microbiome20 and/or transfer of microbes19,33. In keeping with the prediction that local algal competition will increase disease27,33 (for example, through dissolved organic carbon27 or other mechanisms33), we found that the presence of algal competition corresponded to an increased prevalence of DSS52 in Siderastrea corals. These algae-induced changes to coral microbiology and disease prevalence correlate with increased long-term coral tissue loss and mortality in the eld. Thus, our work shows that microbial interactions that were predicted to be important drivers of coral mortality in laboratory studies can be induced by overshing and nutrient pollution, with important and ecologically meaningful effects for the long-term health of corals.

Our long-term results also provide strong contrasts with existing models of coral microbiome dynamics. Qualitative and quantitative models of the coral microbiome have generally assumed that stress would shift coral microbiomes towards stable states dominated by pathogens (for example, ref. 25). On this basis, we predicted that coralalgal competition, nutrient enrichment and/or above-average temperatures would shift coral microbiomes to alternative stable states. While we did conrm the general idea that multiple stressors increase the abundance of fast-growing opportunists on the coral surface, our data do not support a model in which overshing and nutrient pollution, ensuing algal competition, or seasonal temperature variation drive coral microbiomes to specic alternative stable states. Instead, we found that above-average temperatures or algal contact destabilized the microbiome, shifting microbial communities from a stable to an unstable conguration. This nding mirrors recent reports in disturbed human and primate microbiomes42,53, suggesting an underexplored pattern that may be common to many hostmicrobe systems.

NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833

In other animal systems, increased microbiome variability during stress is thought to reect decreased ability of the host, or its native microbiota, to regulate microbial community composition42. Here we show that algal competition and periods of above-average temperature intersect to inuence, which bacteria dominate coral microbiomes. Both of these factors promote stochastic blooms of opportunistic bacteria that displace typically dominant members of the coral surface mucus layer. Microbiomes dominated by putatively opportunistic orders, such as Vibrionales, Rhodobacterales, Flavobacterales and Cytophagales were much less stable (as measured by increased b-diversity of microbial communities) than those dominated by

Synechococcales. This is interesting because some coral pathogens, such as Vibrio coralliilyticus have been shown to suppress coral innate immune pathways36. Immune suppression may, in turn, affect the ability of a coral to regulate its microbiome, increasing temporal and/or inter-colony variation in microbiome composition.

Similarly, many reports have quantied the anti-microbial properties of coral mucus14, including the contribution of antibiotic-producing Actinobacteria54. Our long-term data set suggests that Actinobacteria are important for suppressing opportunists, as outbreaks of Proteobacteria opportunists were more common when Actinobacteria were in low abundance on the coral surface. In agreement with this model, corals exposed to above-average temperature or algal contact showed lower predicted abundances of microbial pathways involved in antibiotic production. Algal contact also reduced the abundance of Actinobacteria, suggesting that algal competition following herbivore removal will lower the natural defenses of corals against potential pathogens, including genera such as Vibrio that are especially problematic at high temperatures.

Finally, we nd that local overshing and nutrient pollution interact with seasonal temperature variation to render corals more vulnerable to blooms of harmful bacteria and increased mortality during the warmest months. In contrast, mild, sub-bleaching thermal stress during summer months did not increase coral mortality under control conditions. Connections between local stressors and temperature variation are often discussed in the eld in the context of climate change, and have been incorporated into models of reef vulnerability or resilience55,56, as well as coral disease susceptibility15. However, experimental evidence connecting laboratory studies of microbial dynamics to coral mortality in the eld has been lacking. In our data, bacterial opportunism increased at temperatures around the local MMM, consistent with past predictions of a threshold at which bacterial pathogenesis becomes especially problematic based on laboratory studies15. Also in agreement with previous coral microbiome eld studies16, we observe blooms of Vibrio during periods of above-average temperatures. Motility genes were enriched in these microbiomes, supporting the idea that chemotaxis towards thermally stressed corals, previously shown in microuidic experiments, may play an important role in coral microbiome dynamics at above-average temperatures34. We extend these observations both by documenting how various combinations of algal competition and temperature favour different bacterial opportunists, and linking these blooms to losses of protective symbionts such as Actinobacteria caused by local stressors. Finally, our data connects microbiome destabilization caused by these blooms to coral tissue loss in the experiment overall, and especially in periods of high temperature. Thus, multiple lines of evidence collected in this study support an ecologically relevant role for coral microbiomes in mediating coral mortality driven by the intersection of local stressors and seasonal temperature variation.

Together, our results provide experimental data linking prevailing models of how human impacts alter reef ecology6 with models of how coral microbiomes respond to algal competition and temperature15,25. They show that overshing and nutrient pollution increase the vulnerability of corals to blooms of opportunistic microorganisms and that the impacts of these local stressors are exacerbated by above-average temperatures. Importantly, the coral species that suffered high mortality rates in our experiments are now some of the most abundant on Caribbean reefs57,58. Thus, some coral species that have withstood the recent decline of more vulnerable relatives may nonetheless be susceptible to increasing local stressors. Clearly, sufciently extreme thermal anomalies and mass bleaching events will kill corals regardless of local factors. However, our work suggests that conserving natural trophic interactions by protecting herbivorous shes and reducing nutrient pollution may help stabilize coral microbiomes and shield corals against temperature-driven bacterial opportunism and mortality, at least in the near term15.

Methods

Experimental design. To simulate the effects of overshing, nutrient loading or the combination of these stressors, we conducted a 3-year eld experiment. Four pairs of 9-m2 plots were established. One member of each of these pairs was enriched with nitrogen and phosphorous, while the other remained at ambient nutrient levels (Supplementary Fig. 1). These plots were 410 m from each other in all cases. Each 9-m2 plot was delineated into nine 1-m2 subplots with metal nails driven into the reef at the corners and centre of each plot. The locations of the plots were selected such that initial variation in rugosity and algal cover within each subplot was minimal. Within each plot, two randomly selected subplots were enclosed with herbivore exclosures, while two other random subplots were selected as exclosure controls. Exclosure controls were tted with open-topped exclosures. These controls allowed access by herbivorous shes, but acted as controls for other potential artifacts of the cages.

All exclosures were made of plastic-coated wire mesh with 2.5-cm diameter holes. This diameter mesh generally excludes most shes 410 cm total length.

Smaller or juvenile herbivorous shes are able to enter the exclosures, but these smaller herbivores generally contribute little to overall grazing rates on reefs and have minimal impacts on the algal communities59. In addition, access by smaller herbivores reects patterns seen under intensive shing, in which larger sh species are preferentially harvested while leaving smaller size classes of sh60,61. We scrubbed both exclosures and exclosure controls every 46 weeks to remove fouling organisms.

Nutrient pollution was simulated using slow-release fertilizer diffusers applied to each nutrient enrichment plot. Each diffuser was a 15-cm diameter PVC tube, perforated with six 1.5 cm holes. The open ends of the PVC tube were wrapped in ne plastic mesh to keep fertilizer pellets inside, but allow diffusion of soluble nutrients. A total of 175 g of Osmocote (19612, NPK) slow-release fertilizer was loaded into each diffuser. PVC enrichment tubes were attached to each metal nail within the 9-m2 enrichment plots for a total of 25 enrichment tubes per enrichment plot. Nutrients were replaced every 3040 days to ensure continued delivery of N and P. Previous studies have shown Osmocote delivery using similar methods to be an effective way of enriching water column nutrients in benthic systems (for example, ref. 62).

Conrmation of the efcacy of nutrient enrichment. Nitrogen and phosphorus levels were assessed in the water column above each enrichment and control plot as in ref. 52. In July of 2009, 2010 and 2011, divers used 60-ml syringes to slowly draw water from B3 cm above the benthos in either control or enriched plots. Samples were taken B30 days after nutrient diffuser deployment to ensure that enrichment occurred across the duration of diffuser deployment. Immediately after collection, samples were ltered (GF/F) into acid-washed bottles, placed on ice, returnedto the laboratory and frozen until analysed. Dissolved inorganic nitrogen(DIN ammonium and nitrite nitrate) and soluble reactive phosphorus

concentrations were determined via autoanalyser. We also assessed nutrient enrichment efciency by analysing tissue carbon:nitrogen (C:N) levels in the common alga Dictyota menstrualis. The nutrient content of macroalgae such asD. menstrualis reects ambient nutrient conditions over longer time scales (that is, weeks to months) than ambient water column nutrients (ref. 63 and references therein). We collected D. menstrualis from both enriched and control treatments during the same months as water samples. Tissues were dried at 60 C, ground to a powder and analysed for %C and %N content with a CHN Carlo-Erba elemental analyzer (NA1500). Nutrient data from both water and algal tissue for each replicate were averaged across summers for statistical analysis via analysis of variance (ANOVA).

8 NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833 ARTICLE

Quantication of the herbivorous sh community. Periodically throughout the study, we used 30 2 m belt transects (n 8) to quantify the density of different

herbivorous shes. Divers slowly swam the length of each transect counting individuals of the different herbivorous shes in the genera Sparisoma, Scarus, Acanthurus and Kyphosus. Fishes were identied to species and their lengthwas estimated to the nearest centimetre. We used published length:weight relationships64 to convert sh densities into herbivore biomass. We analysed these data with one-factor ANOVA examining potential differences in herbivore biomass over time.

Quantication of benthic cover. At least once every season (for example, spring, summer, fall, winter at 1214 week intervals), we visually quantied benthic cover within four, 50 50 cm quadrats in each of the 1 m2 treatment areas. These

quadrats were divided into 49 points, and benthic organisms under each point were identied to species or genus. Algae that are challenging to identify taxonomically under eld conditions (for example, crustose coralline algae and lamentous algae) were classied into algal functional groups. Filamentous algae were classied into short algal turf (o0.5 cm in height) or algal turf (40.5 cm in height) given that taller, thicker algal turf can often be deleterious to coral health and growth10.

Benthic cover was quantied in June 2009 1 week before treatments were initiated to provide a baseline from which to assess changes in algal abundance and community structure. No signicant differences among treatments in algal abundance could be detected at the beginning of the experiment (see initial time points in Fig. 1a,b), as expected given random assignment of subplots to treatment conditions. Further, during the summer of each year (20092012) when algal cover was often at its highest, we also surveyed open areas of reef (areas that did not have three-sided exclosure controls) within the 9-m2 plots to assess whether the exclosure controls had any effect on algal abundance or community composition. We did not detect any differences in algal abundance or community composition between the open unmanipulated areas and exclosure controls (Supplementary Data 1).

Coral tissue growth or loss analyses. At the beginning of the experiment, we mapped each coral colony in the experimental plots that were 42 cm in diameter and took close-up photographs of these corals in situ. Subsequently, we photographed each of these corals every B16 weeks throughout the experiment for a photographic record of changes in coral colony health. In each picture a ruler or object of known size was placed next to the coral to provide scale. In total, we tracked the fate of 226 individual corals spread across each of the treatments for over 3 years. The most common corals were Porites porites (41.1 % of corals), Agaricia spp. (17.7 % of corals), Siderastrea siderea (15.5 % of corals) andP. astreoides (11.5 % of corals).

These corals allowed us to evaluate the impact of the different treatments on

coral growth or tissue loss across the time course of the experiment. We scored growth or tissue loss on a 12-point scale, with bins corresponding to amounts of tissue loss that could be readily observed in photographs (for example,

2 1025% tissue loss). We scored the tissue loss or gain of each coral over the

course of the experiment on the following scale: 6 100% tissue loss, 5

7590% loss, 4 5075% loss, 3 2550% loss, 2 1025% loss, 1

010% loss, 0 0% loss/gain, 1 010% gain, 2 1025% gain, 3 2550% gain,

4 5075% gain, 5 7590% gain and 64100% gain. We then converted these

scores to mean loss/gain by averaging the range corresponding to that score. For example, a coral with a 3 score would be converted to a 37% tissue loss value.

Only nine corals grew 4100% (score 6) over the course of the experiment. For

these corals, we estimated the growth for each coral at 100500% at 50% intervals (for example, 100, 150, 200% and so on). Statistical analyses were conducted based on the raw tissue gain/loss scores, but converted to percentages in the presentation for ease of interpretation. Further, at each time point we scored each coral for:(1) algal competition as measured by direct contact with algal competitors (and the identication of that algal competitor), (2) the presence of overlying sediment on the coral, (3) predation scars from parrotshes and invertebrate corallivores (only the former were observed at appreciable levels), and (4) signs of bleaching or disease. The primary coral disease observed was DSS (see ref. 52 for additional discussion).

Statistical analysis of benthic cover and coral health. Algal cover was analysed using mixed effects models to determine whether the response variable differed among enrichment treatment, herbivore treatment and seasons, as well as whether there was an enrichment herbivore season interaction (Supplementary Data 1

sheet c). The nested, split-plot design of the experiment was incorporated into the model by nesting replicates of the exclosures and exclosure controls within ambient or nutrient-enriched plots. We analysed cover for important species or functional groups, as well as for overall upright algal cover, which is a proxy for the competitive environment of corals. Upright algal cover included all macroalgae and tall lamentous turf, but excluded crustose coralline algae and short lamentous turf as these two functional groups are relatively benign for corals10. In addition, we assessed how herbivore exclusion, nutrient pollution and season impacted algal community structure via PCoA of BrayCurtis divergences, as well as permutational MANOVA (PERMANOVA).

Per cent coral mortality per treatment and coral tissue loss were analysed using similar mixed models to algal cover. For growth measures, corals were nested within ambient or enriched plots, but we did not incorporate season as we only analysed change in tissue for corals at the end of the experiment. We calculated tissue loss statistics either excluding or including corals that suffered total colony mortality. Corals that died suffered total colony morality and therefore 100% loss of live tissue area. Including these corals in coral growth analyses resulted in non-normal distributions that could not be corrected via transformations. Therefore, we analysed coral growth both excluding the corals that died, which satised normality requirements for the analyses, and including the corals that died. Both analyses produced relatively consistent results (Supplementary Data 5), with the exception that only the interaction of herbivory nutrient loading was

signicant in Porites corals (rather than each factor also being individually signicant) when total colony mortality was excluded. We used a w2-test to determine if coral mortality was higher or lower than expected across different seasons given the null hypothesis that coral mortality would be distributed evenly across seasons (25% of total mortality per season).

To assess the impact of algal competition and parrotsh predation on coral mortality for Porites spp., (the corals with the highest mortality rates), we used Fishers exact test to test for differences in the proportion alive versus dead corals across enrichment and/or herbivore exclusion treatments. We also used Fishers exact tests to determine how factors such as competition with different types of algae (that is, corals with or without algal contact) led to loss/gain of coral tissue or coral death/survival.

Mixed effects models were run using the nlme package and post hoc comparisons conducted using the multcomp package in R v3.0.0. Fishers exact tests on contingency tables were run using JMP software (SAS). PCoA of BrayCurtis divergences and PERMANOVA of benthic community data were conducted in QIIME 1.8.

Field sampling of coral mucus for microbial analysis. Coral mucus microbial communities were studied in depth for three coral genera common to the study region and most abundant within the plots: Siderastrea (Siderastrea siderea only), Porites and Agaricia. We used 16S rRNA gene surveys to study the microbial changes in the coral mucus across the course of treatment, focusing especially on whether detectable microbial changes accompanied specic routes to mortality or tissue loss in corals. (For metadata and sequencing statistics, see Supplementary Data 2; for OTU table in BIOM format, see Supplementary Data 7). We focused on mucus communities because these are thought to provide a barrier against invasion by opportunistic pathogens, and can be sampled non-destructively from the same individuals over time. We deemed other methods, such as collecting live coral tissue, too harmful and invasive to the coral for our goal of monitoring the coral microbiome over the long term.

Coral-associated bacteria and archaea were collected using sterile syringe removal of the coral surface mucus layer on SCUBA. A sterile syringe was used to rst agitate the coral and then remove 10 ml of mucus using negative pressure. Once on the surface, the sample was placed in a sterile 15 ml conical tube, frozen on dry ice for the return trip and kept frozen until processing. At a number of time points, 15 ml water sample controls were collected from 41 m above the reef and treated identically to mucus samples in downstream processing. In this study, we did not attempt to evaluate changes in Symbiodinium abundance or taxonomy, as it has not yet been established that mucosal abundances of Symbiodinium types reect abundances within tissues, and destructive tissue sampling would preclude time-series analysis of the same coral colonies.

Quantication of sub-bleaching thermal stress. We calculated cumulative thermal stress using the DHW metric of National oceanic and atmospheric administration (NOAAs) Coral Reef Watch47. These calculations rely on a climatological baseline for the mean temperature of the warmest month known as the MMM. We calculated the MMM for our site using the NOAA Pathnder v5.2 data set, the US ofcial climate data record for sea-surface temperature45, and found that the MMM for our site was 29.26 C (Supplementary Data 6). Data used spanned 19822008, and excluded study dates. Temperatures that exceed the MMM by 1 C are generally regarded as constituting thermal stress47. We

calculated an MMM 1 C threshold of 30.26 C for our site. Temperatures above

the 30.26 C threshold for accumulation of DHWs at our site only occurred during the warm seasons (3 months centred on the warmest month, August).

DHWs represent the extent to which temperatures exceeded the MMM 1 C

in a given season. Temperatures were above 30.26 C for 7 weeks (total) during all sampling years, resulting in the accumulation of o1 DHW during each of the study years (Supplementary Fig. 10a). In predictions of coral bleaching, accumulation of four DHWs is often associated with minor to moderate bleaching5. We saw little bleaching within experimental plots, consistent with sub-bleaching levels of thermal stress. Sub-bleaching thermal stress may nonetheless negatively affect reef health if it increases the abundance and/or virulence of bacterial pathogens. For example, one recent model15 proposes that bacterial pathogens may become problematic for corals at temperatures, exceeding the MMM (a lower threshold than the MMM 1 C threshold used in coral

bleaching studies). We therefore related our microbiological data to both the

NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833

MMM 1 C for coral thermal stress during coral bleaching, and the MMM

threshold proposed for bacterial pathogens (Supplementary Fig. 10b).

Coral mucus microbiome sample processing and sequencing. In the laboratory, coral mucus samples were thawed, centrifuged and supernatant decanted. DNA was puried using an organic extraction as previously described26. After DNA extraction, microbial 16S amplicon libraries were generated using the primers 515F and 806R, both with added 454 sequencing adaptors and with Golay barcodes added to the reverse primer. Triplicate 25 ml reactions were conducted using the

GoTaq Flexi system from Promega (Madison, WI, USA) with the following conditions per reaction: 1 clear buffer, 1 mM dNTPs, 5 mM MgCl, 1 mM of each

primer, 1u Taq polymerase and 1 ml of extracted DNA template. Thermocycling was conducted as follows: 1 cycle of 94 C for 3 min; 35 cycles of 94 C for 45 s, 50 C for 60 s and 72 C for 90 s; and 1 cycle of 72 C for 10 min. Amplication success was checked on a 1.5% agarose gel, and successful triplicate reactions were pooled and cleaned using AMPure magnetic beads from Agencourt. Before sequencing, libraries were quantied using a Qubit dsDNA HS kit from Invitrogen and then pooled into equimolar ratios. The pooled library was checked for amplicon length and purity on an Agilent Bioanalyzer 2100 and then sequenced on a 454 Roche pyrosequencer (GSJunior platform) at the Oregon State Universitys Center for Genome Research and Biocomputing Core Laboratories.

Microbial community data quality control. The QIIME (v.1.8) software pipeline65 was used for quality control, selection of operational taxonomic units and analyses of community diversity (Supplementary Data 2). Sequence libraries were demultiplexed, and sequences with quality scores less than a mean of 35 were removed. Error-correcting barcodes were used to detect and recover sequences whose barcode sequence had exactly one sequencing error. Barcode sequences with two or more errors were removed. Sequences were clustered into operational taxonomic units (OTUs), at a 97% 16S rRNA gene identity threshold using USEARCH 6.1.544, and the subsampled open-reference OTU-picking protocol in QIIME v.1.8 (ref. 65, using greengenes 13_8 as the reference66. This OTU-picking protocol clusters all reads, but assigns reference ids to OTUs in greengenes, which can be useful in comparisons across studies. Chimeric sequences were removed with UCHIME. OTUs represented in the overall analysis by only a single count (singletons) account for a large proportion of noisy reads. Because our emphasis was overall community trends (rather than exploration of the rare biosphere of corals), singleton sequences were removed. Representative sequences for each OTU were classied taxonomically according to the greengenes taxonomy version 13_866 using the Ribosomal Database Project (RDP) Classier software v. 2.2. Sequence alignment and phylogenetic inference for the representative sequence of these OTU is described below in the context of b-diversity analysis.

We took additional steps to account for aspects of the data set unique to host-associated samples. Because coral mucus can contain some amounts of sloughed tissue, we tested whether coral mitochondria were present in any mucus samples. Similarly, because Symbiodinium and other photosynthetic microbial eukaryotes frequently inhabit coral mucus, chloroplast sequences are frequently observed in microbial diversity surveys of corals. As our interest was primarily in cellular bacteria and archaea rather than organelles, coral mitochondrial sequences and 16S sequences classied as chloroplasts with at least 70% condence by the RDP classier were removed before analysis.

Because we observed that many mitochondrial sequences were not efciently identied by RDP, we also removed sequences with very high (1 e 50) sequence

similarity and at least 90% sequence similarity to a reference coral mitochondrial sequence. A table of the coral mitochondrial sequences used is available as Supplementary Data 2 sheet b. These thresholds were selected with care to avoid indiscriminate removal of a-proteobacteria sharing evolutionary ancestry and 16S rRNA gene sequence similarity with mitochondria. The e-value for mitochondrial removal was selected by testing several e-values (1 e 10,1 e 30,1 e 50 and

1 e 100). For each e-value the best BLAST (Basic Local Alignment Search Tool)

similarity against NCBIs nr database was examined. We selected the highest (most lenient) e-value that removed mitochondria, but not related bacteria. At a 1 e 50

BLAST e-value threshold, the best BLAST hit of all removed 16S rRNA sequences was coral mitochondria, and this threshold was therefore selected for the screen. The nal OTU table after quality control (QC) used for our analysis described in this manuscript can be found in Supplementary Data 7.

Analysis of microbial b-diversity. Microbial community b-diversity was calculated based on the weighted UniFrac distance matrix67. This is a phylogenetic measure of community similarity that takes into account organismal abundance and phylogeny67. Phylogenetic trees used to calculate this metric were constructed in QIIME 1.8 (ref. 65) through alignment of representative sequences of each OTU with PyNAST against the greengenes core set alignment66, and approximate maximum likelihood phylogenetic inference with FastTree. We consideredthe pool of distances between samples within each metadata category of interest (for example, algal competition or categories of temperature) using QIIMEs make_distance_boxplots.py. Signicance was assessed by non-parametric t-tests, each with 1,000 Monte Carlo permutations (permutation is important in this instance to account for the non-independence of distances). The effect of this

procedure is to ask whether different factors increase the dispersion of communities. PCoA plots of b-diversity were visualized in the Emperor software.

When multiple categories (for example, different algal types) were tested for effects on b-diversity, the false discovery rate (FDR) for multiple comparisons was controlled at a threshold of q 0.05 using the BejaminiHochberg method.

Comparison of microbial and environmental distance matrices. Mantel tests test for correlation between two distance matrices. For example, a matrix of geographic distances for sample sites might be tested for correlation against a matrix of genetic distances. The partial Mantel test is an extension that tests for correlation between two distance matrices after accounting for the effects of a third, confounding, distance matrix. We used permutational Mantel tests to test whether between-sample variation in continuous environmental factors such as thermal stress, temperature or algal cover correlated with differences in the weighted UniFrac distance matrix67 between coral microbial communities. When data on hypothesized confounding factor was available, we used partial Mantel tests to test signicance after accounting for the confounding parameter. For example, seasonal variation in algal cover might potentially confound the effects of temperature on microbial communities- partial Mantel tests were used to test for such effects. All Mantel and partial Mantel tests were performed in QIIME 1.8 (ref. 65) using the script compare_distance_matrices.py

Analysis of microbial community richness and evenness. Microbial community richness, evenness and b-diversity were calculated in QIIME v.1.8. Richness and evenness were calculated using the chao1 and equitability statistics, respectively, in QIIMEs alpha_diversity.py, collate_alpha.py and compare_alpha.ph scripts. In each case, the data were repeatedly (10 ) raried to 500 reads per sample. Values

for chao1 and equitability were calculated for each raried table, and averaged into a single value and compared across categories. Signicance of results was assessed by FDR-corrected non-parametric t-tests with 1000 Monte Carlo permutations, using qo0.05 as the signicance threshold.

Equitability was calculated as dened in the QIIME software package. Specically, the Shannon entropy (using a base 2 logarithm) was divided by the base 2 logarithm of the number of observed OTUs. Thus, for a sample with n OTUs considered in turn (n1, n2, n3, y, ni), the equitability/evenness was calculated as:

Pni

n n

log2 nn

log2 n

The effects of treatment, temperature, coral genus and individual coral head on richness and evenness were analysed using linear models in R (3.1.1). Seawater temperatures were binned into low temperature, moderate temperature and high temperature categories, with thresholds ato24 C for low temperature

and 430 C for high temperature. The high temperature threshold represents the MMM from the NOAA Pathnder v5.2 SST data set45 (see above).

The specic form of these linear models in the R language was: lm(alpha_diversity B treatment * seawater_ temperature_category * coral_genus

individual_coral_head). The signicance of each effect was compared with ANOVA (Supplementary Data 3 sheet f).

Random forest analysis. Random forest analysis is a machine-learning method for supervised classication. We used random forest analysis (Supplementary Data 3 sheets ij) to (i) determine the conditions under which Synechococcales would be displaced as the most abundant order in the coral microbiome, and (ii) generate a model predicting whether specic coral colonies would lose tissue using ecological (for example, long-term contact with algae from photo series) and microbiological data (the abundance of different bacterial orders).

For both analyses, random forest was conducted in QIIME v.1.8 with the script supervised_learning.py. In both models, the random forest was constructed with 1,000 decision trees and validated by 10-fold cross-validation. This analysis resulted in an overall prediction accuracy, a per-category prediction accuracy and a ranking of feature importance, dened as the amount of accuracy gained or lost when a specic feature was not included in the model. It is worth noting that we employed random forest analysis primarily to provide a summary of the overall predictive power of the data; the internal complexity of the algorithm does not allow inferences to be made about the relationships between variables, beyond their empirical inuence on accuracy. For mechanistic insight we rely on other analyses presented here.

In the rst model, we sought to predict the dominant (most abundant) order of bacteria in samples given only non-microbial data about the coral colony from which the sample was taken at that point in time. A total of 127 features were used in the prediction including: herbivore exclusion, nutrient addition, coral species (species were converted to binary variables), per cent cover of different macroalgae, turf algae and cyanobacterial mats from quadrat surveys, singly and grouped by functional type (for example, total macroalgae), plot and subplot number, measures of instantaneous and weekly average temperatures and salinities, season (quarter), and the simultaneous presence of parrotsh bites and nutrient enrichment. Given these features, we asked the model to predict the dominant order of microorganism in each DNA sample, with our main purpose being to determine whether the

10 NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833 ARTICLE

conditions under which Synechococcus was dominant were predictable. No microbiological data were included in the input feature table, so all predictions about the microbiology were based on externally observable features.

In the second model, we sought to test whether data about the competitive environment and microbial community composition of specic coral colonies could predict whether they would lose tissue over the 3-year study period. This model included features representing herbivore exclusion, nutrient pollution and coral genus; long-term competition with any algae or (as separate features) with turf algae, visible Cyanobacteria mats, Dictyota, Sargassum or Halimeda; the presence of bleaching, DSS or signs of parrotsh predation; the initial diameter of the coral colony in centimeters; and microbiological data (averaged across all samples from each individual) on the abundance of the phyla Proteobacteria, Cyanobacteria, Bacteroidetes, Firmicutes, Acidobacteria and Actionobacteria overall, and six specic microbial orders of interest due to their prominence in this data set and/or the literature13,41,44 (Synechococchales, Vibrionales, Rhodobacterales, Alteromonadales, Rickettsiales, Rhodospirallales and Pseudomonadales).

Predicting coral microbiome function. To predict potential functional consequences of observed changes in microbial taxonomy across treatments, functional proles for each microbial sample were predicted using the PICRUSt tool40. This tool uses hidden state prediction, a form of evolutionary modelling closely related to ancestral state reconstruction, to put bounds on genomic copy numbers of each gene family in uncultivated environmental microorganisms, using their position in a reference bacterial phylogeny relative to all bacteria or archaea with sequenced genomes40,68. These per genome estimates are then combined, taking into account 16S rRNA copy number variation, to estimate a functional prole (or virtual metagenome) for a microbial community based on the information available in all sequenced genomes. The resulting metagenomes have been used in recent analyses of human and environmental microbiomes, and have in at least one study been shown to correlate with overall metabolite proles (reviewed in ref. 68). The accuracy of the method depends on several factors40,68, but the availability of reasonably related reference genomes is generally the most important. Availability of such genomes can be summarized through a Nearest Sequenced Taxon Index (NSTI) score. This score is the average branch length between each member of the community (OTU) and its closest sequenced relative, weighted by abundance. In these data, NSTI scores ranged from 0.03 to 0.18, with a mean of 0.11. This was within the range of soil and mammal microbiomes that have previously been predicted with reasonable accuracy (as assessed using paired 16S rRNA/shotgun metagenomes from the same samples)40.

We used PICRUSt to predict functional proles for all 16S rRNA samples, and summarized these proles of predicted gene family abundance into Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and functional categories. We then tested whether functional changes in coral microbiomes correlated with either temperature increase, extremes of temperature or the abundance of upright algal cover within a subplot (Supplementary Fig. 5a,b). To test the effects of extremes of temperature (either hot or cold) on coral microbiomes, the imputed abundance of KEGG functional categories (level 3) across samples was regressed against the squared deviation of temperature from 28 C, which is the average temperature across samples, and also approximates the annual average temperature at the study site over 32 years (calculated from NOAA Pathnder v5.245). The effects of upright algal cover or increasing temperature were tested by Spearman regression against KEGG functional category abundance. The FDR across KEGG categories was controlled in all cases at qo0.05.

Comparing PICRUSt results with previous experiments. To compare whether the functional categories predicted to change in this eld study were broadly consistent with previous laboratory experiments, we compared PICRUSt data with a previous study that exposed corals to several stressors in aquaria and sequenced their metagenomes. KEGG pathways correlated with extremes of temperature (squared deviation from 28 C, see Predicting coral microbiome function, above) in this study were compared with KEGG categories that increased or decreased in a previous experiment, in which Thurber et al.24, exposed the Pacic coral Porites compressa to thermal stress. Metagenomes from that study24 were annotated with KEGG pathways in MG-RAST using default parameters (BLAST e-value o10 5, 60% coverage and 15% alignable). We compared the 25 categories that increased by 41% with temperature stress to the present study.

Macroalgal contact as a driver of microbial b-diversity. We sought to test whether the increase in b-diversity with macroalgal contact that we saw in this data set is a common pattern in coralalgal competition. We reanalysed data from a previously published experiment19 that studied the effects of macroalgal contact on coral microbiomes of P. astreoides (Supplementary Fig. 7c). In that study, macroalgae were placed in direct contact with P. astreoides corals, and the microbial communities of the macroalgae, corals without algal competitors and corals with algal competitors assessed using terminal restriction fragment length polymorphism19. We translated the table of terminal restriction fragments from this study into a QIIME-compatible .biom format OTU table, and tested whether b-diversity (measured by BrayCurtis divergence in terminal restriction fragment

proles) was greater in corals exposed to macroalgae than control corals or control algae using analyses described above (see Analysis of Microbial b-diversity).

Data availability. Processed microbial data are available as Supplementary Data 7 (OTU tables in BIOM format). Metadata for coral microbiome samples are available as Supplementary Data 2 sheets c,d, respectively. Short-read amplicon sequence data (as fasta and qual les) are deposited in the QIITA database, as study id: 10482. Other relevant data are available from the authors.

References

1. Jackson, J., Donovan, M. K., Cramer, K. L. & Lam, V. V. Status and Trends of Caribbean Coral Reefs: 19702012 (Global Coral Reef Monitoring Network, IUCN, Gland, 2014).

2. Bruno, J. F. & Selig, E. R. Regional decline of coral cover in the Indo-Pacic: timing, extent, and subregional comparisons. PLoS One 2, e711 (2007).

3. Wilkinson, C. Status of Coral Reefs of the World: 2008 (Global Coral Reef Monitoring Network and Reef and Rainforest Research Centre (2008).

4. Bruno, J. F. et al. Thermal stress and coral cover as drivers of coral disease outbreaks. PLoS Biol. 5, 12201227 (2007).

5. Eakin, C. M. et al. Caribbean corals in crisis: record thermal stress, bleaching, and mortality in 2005. PLoS One 5, e13969 (2010).

6. Littler, M. M., Littler, D. S. & Brooks, B. L. Herbivory, nutrients, stochastic events, and relative dominances of benthic indicator groups on coral reefs: a review and recommendations. Smithson. Contrib. Mar. Sci 38, 401414 (2009).

7. Smith, J. E., Hunter, C. L. & Smith, C. M. The effects of top-down versus bottom-up control on benthic coral reef community structure. Oecologia 163, 497507 (2010).

8. Burkepile, D. E. & Hay, M. E. Herbivore vs. nutrient control of marine primary producers: context-dependent effects. Ecology 87, 31283139 (2006).

9. Hughes, T. P. et al. Phase shifts, herbivory, and the resilience of coral reefs to climate change. Curr. Biol. 17, 360365 (2007).

10. Birrell, C. L., McCook, L. J., Willis, B. L. & Diaz-Pulido, G. A. Effects of benthic algae on the replenishment of corals and the implications for the resilience of coral reefs. Oceanogr. Mar. Biol. Annu. Rev 46, 2563 (2008).

11. Graham, N. A., Jennings, S., MacNeil, M. A., Mouillot, D. & Wilson, S. K. Predicting climate-driven regime shifts versus rebound potential in coral reefs. Nature 518, 9497 (2015).

12. Hoegh-Guldberg, O. et al. Coral reefs under rapid climate change and ocean acidication. Science 318, 17371742 (2007).

13. Lesser, M. P., Mazel, C. H., Gorbunov, M. Y. & Falkowski, P. G. Discovery of symbiotic nitrogen-xing cyanobacteria in corals. Science 305, 9971000 (2004).

14. Ritchie, K. B. Regulation of microbial populations by coral surface mucus and mucus-associated bacteria. Mar. Ecol. Prog. Ser. 322, 114 (2006).

15. Maynard, J. et al. Projections of climate conditions that increase coral disease susceptibility and pathogen abundance and virulence. Nat. Clim. Change 5, 688694 (2015).

16. Tout, J. et al. Increased seawater temperature increases the abundance and alters the structure of natural Vibrio populations associated with the coral Pocillopora damicornis. Front. Microbiol 6, 432 (2015).

17. Smith, J. et al. Indirect effects of algae on coral: algae-mediated, microbe-induced coral mortality. Ecol. Lett. 9, 835845 (2006).

18. Vermeij, M. J., Smith, J. E., Smith, C. M., Vega Thurber, R. & Sandin, S. A. Survival and settlement success of coral planulae: independent and synergistic effects of macroalgae and microbes. Oecologia 159, 325336 (2009).

19. Vega Thurber, R. et al. Macroalgae decrease growth and alter microbial community structure of the reef-building coral, Porites astreoides. PLoS One 7, e44246 (2012).

20. Morrow, K., Paul, V., Liles, M. & Chadwick, N. Allelochemicals produced by Caribbean macroalgae and cyanobacteria have species-specic effects on reef coral microorganisms. Coral Reefs 30, 309320 (2011).

21. Jessen, C. et al. In-situ effects of eutrophication and overshing on physiology and bacterial diversity of the Red Sea coral Acropora hemprichii. PLoS One 8, e62091 (2013).

22. Rohwer, F., Breitbart, M., Jara, J., Azam, F. & Knowlton, N. Diversity of bacteria associated with the Caribbean coral Montastraea franksi. Coral Reefs 20, 8591 (2001).

23. Frias-Lopez, J., Zerkle, A. L., Bonheyo, G. T. & Fouke, B. W. Partitioning of bacterial communities between seawater and healthy, black band diseased, and dead coral surfaces. Appl. Environ. Microbiol. 68, 22142228 (2002).

24. Thurber, R. V. et al. Metagenomic analysis of stressed coral holobionts. Environ. Microbiol. 11, 21482163 (2009).

25. Mao-Jones, J., Ritchie, K. B., Jones, L. E. & Ellner, S. P. How microbial community composition regulates coral disease development. PLoS Biol. 8, e1000345 (2010).

26. Kline, D. I., Kuntz, N. M., Breitbart, M., Knowlton, N. & Rohwer, F. Role of elevated organic carbon levels and microbial activity in coral mortality. Mar. Ecol. Prog. Ser. 314, 119125 (2006).

NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11833

27. Rohwer, F., Youle, M. & Vosten, D. Coral Reefs In The Microbial Seas (Plaid Press United States, 2010).

28. Haas, A. F. et al. Effects of coral reef benthic primary producers on dissolved organic carbon and microbial activity. PLoS One 6, e27973 (2011).

29. Gregg, A. et al. Biological oxygen demand optode analysis of coral reef-associated microbial communities exposed to algal exudates. PeerJ 1, e107 (2013).

30. Barott, K. et al. Hyperspectral and physiological analyses of coral-algal interactions. PLoS One 4, e8043 (2009).

31. Wangpraseurt, D. et al. In situ oxygen dynamics in coral-algal interactions. PLoS One 7, e31192 (2012).

32. Kuntz, N. M., Kline, D. I., Sandin, S. A. & Rohwer, F. Pathologies and mortality rates caused by organic carbon and nutrient stressors in three Caribbean coral species. Mar. Ecol. Prog. Ser. 294, 173180 (2005).

33. Nugues, M. M., Smith, G. W., Hooidonk, R. J., Seabra, M. I. & Bak, R. P. Algal contact as a trigger for coral disease. Ecol. Lett. 7, 919923 (2004).

34. Garren, M. et al. A bacterial pathogen uses dimethylsulfoniopropionate as a cue to target heat-stressed corals. ISME J. 8, 9991007 (2014).

35. Kimes, N. E. et al. Temperature regulation of virulence factors in the pathogen Vibrio coralliilyticus. ISME J. 6, 835846 (2012).

36. Vidal-Dupiol, J. et al. Thermal stress triggers broad Pocillopora damicornis transcriptomic remodeling, while Vibrio coralliilyticus infection induces a more targeted immuno-suppression response. PLoS One 9, e107672 (2014).

37. Rypien, K. L., Ward, J. R. & Azam, F. Antagonistic interactions among coral-associated bacteria. Environ. Microbiol. 12, 2839 (2010).

38. Box, S. J. & Mumby, P. J. Effect of macroalgal competition on growth and survival of juvenile Caribbean corals. Mar. Ecol. Prog. Ser. 342, 139149 (2007).

39. Rasher, D. B. & Hay, M. E. Chemically rich seaweeds poison corals when not controled by herbivores. Proc. Natl Acad. Sci. USA 107, 96839688 (2010).40. Langille, M. G. et al. Predictive functional proling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814821 (2013).

41. Nelson, C. E. et al. Coral and macroalgal exudates vary in neutral sugar composition and differentially enrich reef bacterioplankton lineages. ISME J. 7, 962979 (2013).

42. Moeller, A. H. et al. SIV-induced instability of the chimpanzee gut microbiome. Cell Host Microbe 14, 340345 (2013).

43. Whittaker, R. H. Evolution and measurement of species diversity. Taxon 21, 213251 (1972).

44. Roder, C., Arif, C., Daniels, C., Weil, E. & Voolstra, C. R. Bacterial proling of White Plague Disease across corals and oceans indicates a conserved and distinct disease microbiome. Mol. Ecol. 23, 965974 (2014).

45. Casey, K., Brandon, T., Cornillon, P. & Evans, R. in Oceanographyfrom Space. (eds Barale, Vittorio, Gower, J. F. R. & Alberotanza, L.) Ch. 16, 273287 (Springer Netherlands, 2010).

46. Chassignet, E. P. et al. The HYCOM (hybrid coordinate ocean model) data assimilative system. J. Mar. Syst 65, 6083 (2007).

47. Liu, G., Strong, A. E. & Skirving, W. Remote sensing of sea surface temperatures during 2002 Barrier Reef coral bleaching. EOS 84, 137141 (2003).

48. Rotjan, R. D. & Lewis, S. M. Impact of coral predators on tropical reefs. Mar. Ecol. Prog. Ser. 367, 7391 (2008).

49. Burkepile, D. E. Context-dependent corallivory by parrotshes in a Caribbean reef ecosystem. Coral Reefs 31, 111120 (2012).

50. Raymundo, L. J., Halford, A. R., Maypa, A. P. & Kerr, A. M. Functionally diverse reef-sh communities ameliorate coral disease. Proc. Natl Acad. Sci. USA 106, 1706717070 (2009).

51. Mumby, P. J. Herbivory versus corallivory: are parrotsh good or bad for Caribbean coral reefs? Coral Reefs 28, 683690 (2009).

52. Vega Thurber, R. et al. Chronic nutrient enrichment increases prevalence and severity of coral disease and bleaching. Glob. Change Biol. 20, 544554 (2014).

53. Charlson, E. S. et al. Disordered microbial communities in the upper respiratory tract of cigarette smokers. PLoS One 5, e15216 (2010).

54. Nithyanand, P., Manju, S. & Karutha Pandian, S. Phylogenetic characterization of culturable actinomycetes associated with the mucus of the coral Acropora digitifera from Gulf of Mannar. FEMS Microbiol. Lett. 314, 112118 (2011).

55. Edwards, H. J. et al. How much time can herbivore protection buy for coral reefs under realistic regimes of hurricanes and coral bleaching? Glob. Change Biol. 17, 20332048 (2011).

56. Anthony, K. et al. Ocean acidication and warming will lower coral reef resilience. Glob. Change Biol. 17, 17981808 (2011).

57. Ruzicka, R. R. et al. Temporal changes in benthic assemblages on Florida Keys reefs 11 years after the 1997/1998 El Nino. Mar. Ecol. Prog. Ser. 489, 125141 (2013).

58. Green, D., Edmunds, P. J. & Carpenter, R. C. Increasing relative abundance of Porites astreoides on Caribbean reefs mediated by an overall decline in coral cover. Mar. Ecol. Prog. Ser. 359, 110 (2008).

59. Lokrantz, J., Nystrm, M., Thyresson, M. & Johansson, C. The non-linear relationship between body size and function in parrotshes. Coral Reefs 27, 967974 (2008).

60. Mumby, P. J. et al. Fishing, trophic cascades, and the process of grazing on coral reefs. Science 311, 98101 (2006).

61. Bellwood, D. R., Hoey, A. S. & Hughes, T. P. Human activity selectively impacts the ecosystem roles of parrotshes on coral reefs. Proc R Soc B 279, rspb20111906 (2011).

62. Littler, M. M., Littler, D. S., Brooks, B. L. & Lapointe, B. E. Nutrient manipulation methods for coral reef studies: a critical review and experimental eld data. J. Exp. Mar. Bio. Ecol. 336, 242253 (2006).

63. Fourqurean, J. W. & Zieman, J. C. Nutrient content of the seagrass Thalassia testudinum reveals regional patterns of relative availability of nitrogen and phosphorus in the Florida Keys USA. Biogeochemistry 61, 229245 (2002).

64. Bohnsack, J. & Harper, D. NOAA Technical Memorandum (NMFSSEFC-215; 31, 1988).

65. Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335336 (2010).

66. McDonald, D. et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610618 (2012).

67. Lozupone, C. A., Hamady, M., Kelley, S. T. & Knight, R. Quantitative and qualitative b diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 73, 15761585 (2007).

68. Zaneveld, J. R. & Thurber, R. L. Hidden state prediction: a modication of classic ancestral state reconstruction algorithms helps unravel complex symbioses. Front. Microbiol. 5, 431 (2014).

Acknowledgements

We thank numerous volunteers for assistance with the experiment and the Florida Keys National Marine Sanctuary for permits (FKNMS-2009-047 and FKNMS-2011-090). This research was supported by National Science Foundation grant OCE-1130786 to D.E.B and R.V.T., Florida International University and Oregon State University. R. van Hooidonk of NOAA Atlantic oceanographic & meteorologic laboratory (AOML) provided the Pathnder SST time series data.

Author contributions

D.E.B. and R.V.T. designed the research; all authors conducted the research; J.R.Z., D.E.B., R.M., N.P.L., J.P.P., J.A.M. and R.V.T. analysed the data; J.R.Z., D.E.B. and R.V.T. wrote the manuscript; all authors contributed to revisions.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Zaneveld, J.R. et al. Overshing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales. Nat. Commun. 7:11833 doi: 10.1038/ncomms11833 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

12 NATURE COMMUNICATIONS | 7:11833 | DOI: 10.1038/ncomms11833 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

Word count: 13484

Show less

Abstract

Translate

Losses of corals worldwide emphasize the need to understand what drives reef decline. Stressors such as overfishing and nutrient pollution may reduce resilience of coral reefs by increasing coral-algal competition and reducing coral recruitment, growth and survivorship. Such effects may themselves develop via several mechanisms, including disruption of coral microbiomes. Here we report the results of a 3-year field experiment simulating overfishing and nutrient pollution. These stressors increase turf and macroalgal cover, destabilizing microbiomes, elevating putative pathogen loads, increasing disease more than twofold and increasing mortality up to eightfold. Above-average temperatures exacerbate these effects, further disrupting microbiomes of unhealthy corals and concentrating 80% of mortality in the warmest seasons. Surprisingly, nutrients also increase bacterial opportunism and mortality in corals bitten by parrotfish, turning normal trophic interactions deadly for corals. Thus, overfishing and nutrient pollution impact reefs down to microbial scales, killing corals by sensitizing them to predation, above-average temperatures and bacterial opportunism.

Details

Title

Overfishing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales

Author

Zaneveld, Jesse R; Burkepile, Deron E; Shantz, Andrew A; Pritchard, Catharine E; Mcminds, Ryan; Payet, Jérôme P; Welsh, Rory; Correa, Adrienne M S; Lemoine, Nathan P; Rosales, Stephanie; Fuchs, Corinne; Maynard, Jeffrey A; Thurber, Rebecca Vega

Pages

11833

Publication year

2016

Publication date

Jun 2016

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms11833

ProQuest document ID

1794161017

Overfishing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales

Jump to:

Full Text

Abstract

Details

Suggested sources