Turbulence-driven shifts in holobionts and planktonic microbial assemblages in St. Peter and St. Paul Archipelago, Mid-Atlantic Ridge, Brazil

The aim of this study was to investigate the planktonic and the holobiont Madracis decactis (Scleractinia) microbial diversity along a turbulence-driven upwelling event, in the world's most isolated tropical island, St Peter and St Paul Archipelago (SPSPA, Brazil). Twenty one metagenomes were obtained for seawater (N = 12), healthy and bleached holobionts (N = 9) before, during and after the episode of high seawater turbulence and upwelling. Microbial assemblages differed between low turbulence-low nutrient (LLR) and high-turbulence-high nutrient (HHR) regimes in seawater. During LLR there was a balance between autotrophy and heterotrophy in the bacterioplankton and the ratio cyanobacteria:heterotrophs ~1 (C:H). Prochlorales, unclassified Alphaproteobacteria and Euryarchaeota were the dominant bacteria and archaea, respectively. Basic metabolisms and cyanobacterial phages characterized the LLR. During HHR C:H < < 0.05 and Gammaproteobacteria approximated 50% of the most abundant organisms in seawater. Alteromonadales, Oceanospirillales, and Thaumarchaeota were the dominant bacteria and archaea. Prevailing metabolisms were related to membrane transport, virulence, disease, and defense. Phages targeting heterotrophs and virulence factor genes characterized HHR. Shifts were also observed in coral microbiomes, according to both annotation–indepent and -dependent methods. HHR bleached corals metagenomes were the most dissimilar and could be distinguished by their di- and tetranucleotides frequencies, Iron Acquision metabolism and virulence genes, such as V. cholerae-related virulence factors. The healthy coral holobiont was shown to be less sensitive to transient seawater-related perturbations than the diseased animals. A conceptual model for the turbulence-induced shifts is put forward.


Introduction
Marine microbial communities are recognized as engines of globally important processes, such as the marine carbon, nitrogen and sulfur cycles (Falkowski et al., 2008;Fuhrman, 2009). Only recently with the introduction of molecular techniques have satisfactory descriptions of natural microbial assemblages been generated (Fierer and Jackson, 2006;Rusch et al., 2007;Costello et al., 2009). Nevertheless, most marine ecosystems are understudied. As a result, there is limited information on the diversity of microbial assemblages in changing environments and on the environmental drivers of microbial diversity shifts (Karl, 2002).
Nutrient dynamics in the sea is inextricably linked to variations in physical processes. Either enhanced nutrient delivery from turbulent mixing or upwelling, or enhanced stratification can lead to shifts in microbial assemblages, with significant consequences for nutrient cycling (Cullen et al., 2002). Episodic mixing events must occur in order to balance supply and demand (Hayward, 1987(Hayward, , 1991. Because open-ocean microbial assemblages are metabolically active with a potential for relatively high specific growth, they are poised to respond quickly and effectively to environmental perturbations (Karl, 2002). A model of turbulence-nutrients regimes decoupled characteristics and adaptations of phytoplankton assemblages and how they relate to food web structure: (i) high turbulence-low nutrients: low biomass, slow turnover, adaptations for efficient use of light and nutrients (e.g., iron-limited, high latitude waters); (ii) LLR: smaller cells, high turnover, competition for nutrients, retention by recycling (microbial loop); (iii) low turbulencehigh nutrients: larger cells, higher biomass, slower turnover, selective pressure to sequester nutrients and minimize losses (e.g., noxious toxic blooms); and (iv) HHR: larger cells, higher biomass, transient, and self-limiting selection for rapid growth (e.g., diatoms). According to this model, the microbial loop is present in all regimes but it dominates the biomass in the low turbulence-high nutrients regime (Cullen et al., 2002).
Some physical mechanisms that vertically supply nutrients from below to the oligotrophic oceanic surface layers are: (i) internal waves and tides, (ii) cyclonic mesoscale eddies, (iii) winddriven Ekman pumping, and (iv) atmospheric storms (Karl, 1999). Internal waves are ubiquitous in deep-ocean environments (Garrett and Munk, 1972), and it appears that the high vertical shear of low-frequency internal waves contributes to occasional pulses of vertical mixing (Gregg et al., 1986;Sherman and Pinkel, 1991). Stochastic events, that may be short-lived, are of major ecological significance. They are undoubtedly undersampled by ship-based observation programs (Platt et al., 1989). Even a monthly sampling schedule, as adopted in the Hawaii Ocean Time-series (HOT) program, is considered too infrequent to register important but intermittent nutrient injections (Karl, 1999). Furthermore, most studies are restricted to the Pacific, Caribbean, and North Atlantic.
St. Peter and St. Paul Archipelago (SPSPA) are the smallest and most isolated tropical islands in the world. It comprises the most seaward Brazilian oceanic islands (∼1000 km from mainland) and the unique on the northern hemisphere, lying 100 km off the equator (00 • 55 ′ N; 29 • 22 ′ W). Data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) evinces the chlorophyll pattern characteristic of mesotrophic waters (0.1-0.3 mg.m −3 Chla) (Supplementary Figures 1A,B), according to the classification of Shushkina et al. (1997). SPSPA is located in the biogeochemical Western Tropical Atlantic Province (WTRA) (Longhurst et al., 1995), under the influence of the Intertropical Convergence Zone. Its formed by minute summits of the Mid-Atlantic Ridge (MAR) within the St Paul Fault Zone (FZ), between the South American and the African Plates ( Supplementary Figures 2A,B), where seismicity is frequent (Campos et al., 2010). The islets are devoid of shore and consist entirely of steep scarpments extending to 60-100 m depth, with the most limited area of shallow habitat among oceanic islands (∼200 m) (Robertson, 2001). Within 2 km diameter bathyal depths are reached and within 5 km depths fall within the abyssal range (−3.600 to 5.000; Supplementary Figure 2B). The South Equatorial Current (SEC) flows E-W superficially and the Atlantic Equatorial Undercurrent (EUC) flows W-E at 60-100 m depth (Edwards and Lubbock, 1983). The latter is one of the fastest, varying and least predictable among the Atlantic currents, which reaches 120 cm.seg −1 above the thermocline (Philander, 1986). A permanent thermocline may prevent deep water masses to emerge  nevertheless the eventual presence of deep water Pleurommama spp. and Heterorhabdus spp at the shallow layer during the diurnal period is contradictory to the former assumption (Macedo-Soares et al., 2009). There are no reports on local hydrodynamics but, seemingly, intermittent vertical flow of deep water masses should result from the violent interaction with the geomorphology-MAR-a perpendicular barrier to the currents. In principle, the friction between the SEC and EUC masses, which flow in opposing directions, should promote and intensify episodic extensions of the thermohaline to the photic zone. Moreira et al. (2014) reported a significant enrichment process along an 8 days period during which an ever-growing turbulence with surge pulses was observed in SPSPA. The process occurred along the lunar phase from crescent to full moon. The work performed the first (and unique to date) characterization of the culturable heterotrophic bacterial community of SPSPA. Bacterial counts (colony forming units, CFU) correlated positively with nutrients in seawater, which in turn correlated positively with turbulence-energy and frequency of the surges. In the present work we analyzed the metagenomic composition and diversity of both the planktonic microbial assemblages and in the scleractinian coral M. decactis along the same period in the same locale. The aim was to characterize the microbial diversity during an upwelling-driven nutrient enrichment. We did not expect to find bleached corals in SPSPA. The coral holobionts were targeted in this survey to investigate whether there was a correlation between the metagenomic features and seawater parameters (vibrio counts, nutrients, bacterioplankton composition). Sampling was performed before, during and after a turbulence surge. It was a short-lived event, which is locally recurrent. We analyzed 21 metagenomic samples of seawater (n = 12), healthy and bleached corals (n = 9). The findings are summarized in a model of the physico-chemical-biological dynamics in SPSPA, where a cyclic recurrent pattern with extreme regimes of low turbulence-low nutrients (LLR) and high turbulence-high nutrients (HHR) is hypothesized to contribute to structure the marine ecosystem in that barren archipelago.

Materials and Methods
Field Sampling (Performed by Moreira et al., 2014) In brief, sampling was performed by SCUBA diving at the Sub-caulerpa zone (mesophotic), according to the zonation of (Edwards and Lubbock, 1983;Moreira et al., 2014). The survey took place along the NW side of the archipelago from the inlet (ca 4500 m 2 ) to Belmonte islet's contiguous vertical wall. The satellite view and topograghy of the inlet are shown in Supplementary Figures 1C, 2B, respectively. The turbulence surge occurred along the lunar phase from crescent (14/Sep/2010) to full moon (22/Sep/2010). The peak of the surge overlapped with cloudiness, wind and rain. Rain and swash flushed guano, a possible additional source of nutrients (Gagnon et al., 2013), from the cays into the bay. Samples were obtained from the onset (LLR, t 1 = 14/Sep/2010), while enhancing (t 2 = 15; t 3 = 18/Sep/2010, HHR) and almost to recovery of LLR condition, or recovery for short (t 4 = 22/set/2010). During climax (20-21/Sep/2010) the strong vortex precluded diving (Supplementary Video 1). Henceforth samples will be referred as 14, 15, 18, and 22. In total, 12 colony fragments (10 × 10 cm) of M. decactis (healthy: n = 8; bleached: n = 4) were collected with hammer and chisel. On days 14 and 15 bleached corals were not found. Coral samples processed for metagenomics were: (i) healthy corals (Mad): 14 (n = 1), 15 (n = 2), 18 (n = 1), and 22 (n = 1); and (ii) bleached corals (MadBle): 18 (n = 2) and 22 (n = 2). Seawater was sampled from the water column immediately above the corals (<1 m) (4 samples: 14, 15, 18, and 22; 20 L/sample; 3 sterivex/sample). All samples were taken immediately to the Scientific Staion (SS) Laboratory, 20 m from the pier (Supplementary Figure 1D; view of the pier and the SS from the water at LLR). Seawater was filtered. Filters and coral samples were preserved in liquid nitrogen until DNA extraction (no longer than 3 months after).

Madracis Decactis (Lyman, 1859) (Scleractinia: Pocilloporideae)
Is a colonial zooxanthellate scleractinian coral. It has a variable bathymetric distribution, from 3 to 30 and up to 100 m (Neves and Johnsson, 2009). It is widespread in Brazil (N to SE), in Caribe, Gulf of Mexico and locally found in the Southeastern Atlantic (West Africa) (World Register of marine Species-WORMS: www.marinespecies.org). Free-living colonies of M. decactis display an unique formation off southern Brazil, at Galé Island. Spheroid shape, a.k.a. circumrotatory colonies, form the first corallith site discovered in the subtropical South Atlantic Ocean (SAO), at 6-15 m depth over 3400 m 2 (coralreeefs-2012). In São Paulo (SE Brazil) it's a major contributor of reef structures, where bleaching has been seriously affecting its populations (Migotto, 1997), adding interest to the study of this coral species. In SPSPA, its one of the two scleractinian species locally found, mostly at the mesophotic zone. Healthy, bleached and with scars left by fish predation (Stegastes sanctipauli, Pomacentridae; Halichoeres radiatus, Labridae) in SPSPA are shown in Supplementary Figure 3.
Seawater: Temperature, Nutrients, and Microbial Abundance Environmental parameters were analyzed by standard oceanographic methods with at least three replicates for each parameter and determined by Moreira et al. (2014). Temperature was recorded in situ with a HOBO UA-002-64/date Logger and UEMIS dive computer from 5, 15, 33, 45, and 65 m depth (n = 5 for each depth), during September/2010 and June/2011 (published in Crespo et al., 2014). Environmental data from (Moreira et al., 2014) and temperature from (Crespo et al., 2014) are summarized in Supplementary Figure 4, for aid in data interpretation.

Metagenomic Library Construction
Metagenomes were obtained by pyrosequencing technology using a 454 GS Junior instrument (Roche) (Margulies et al., 2005). Shotgun libraries were generated with 500 ng of whole metagenome samples, sheared into fragments by nebulization. End-repair and adaptor ligation were performed using GS FLX Titanium kit (Roche). Quality control and quantification were performed with Agilent 2100 Bioanalyzer (Agilent Technologies) and TBS 380 Fluorometer (Turner Biosystems), respectively. After the libraries construction, approximately 10 6 molecules/metagenome were denatured and amplified by emulsion PCR.

Metagenomic Data Analysis
Raw sequences were submitted to quality control using PRINSEQ Standalone Lite (version 0.20.4; available at http://sourceforge. net/projects/prinseq/files/). We analyzed 21 metagenomic samples of seawater (n = 12) and corals (n = 9) from 4 days (t 1 = 14, t 2 = 15, t 3 = 18, t 4 = 22, of set/2010) along the enrichment process (t 4 − t 1 = 8 days). Annotation was performed by Meta-Genome Rapid Annotation using Subsystems Technology (MG-RAST) server (Meyer et al., 2008) version 3.0, using (SEED) Subsystems Technology and the GenBank database for functional and organismal classifications, respectively. For this purpose, all BLAST queries were conducted with a maximum cutoff E-value 0.00001, a minimum identity of 60%, and a minimum alignment length of 20 measured in aa for protein and bp for RNA databases.

Metagenomes Comparison Trough Annotation-independent Analysis
Dinucleotide odds ratio and Karlin distances (δ) between the metagenomes, based on the dinucleotide relative abundances differences (according to Karlin et al., 1997), were calculated using Perl scripts as in Willner et al. (2009). The values in the Karlin matrix were multiplied by a 1000 for easier comparison. Tetranucleotide frequencies were calculated using a Perl script as in Albertsen et al. (2013). The divergence between the observed and expected tetranucleotide frequencies was transferred into z-scores and pairwise comparison of the metagenomic sequences was performed by computing the Pearson's correlation coefficients of the z-scores, both through Python scripts, according to (Teeling et al., 2004). Seawater metagenomes  were compared among each other, as well as coral metagenomes (Mad14, Mad15, Mad18, MadBle18, Mad22, MadBle22).

Phage Detection
The 21 metagenomic libraries were searched for phages using the PHAge Search Tool (PHAST), available at http://phast. wishartlab.com (Zhou et al., 2011). Briefly, pyfasta 0.5.2 (available at http://pypi.python.org/pypi/pyfasta/) was used to split metagenomes into smaller subsets without splitting individual fasta entries, after which PHAST was used to phage search. The tool provides an ensemble of ORF prediction and translation (via GLIMMER; Salzberg et al., 1998), protein identification (via BLASTP; Altschul et al., 1997), phage sequence identification (via BLAST matching to a specific database), tRNA identification (using tRNAscan-SE; Lowe and Eddy, 1997), attachment site recognition (with ARAGORN; Laslett and Canback, 2004) and gene clustering density measurements using density-based spatial clustering of applications with noise (DBSCAN; Ester et al., 1996), and sequence annotation text mining. PHAST's database encloses protein sequences from two sources: the NCBI phage database and the prophage database (Srividhya et al., 2006). Specific keywords (e.g., "protease, " "integrase, " and "tail fiber") are used for screening. Matched phage or phage-like sequences with E < 0.0001 are saved as hits and their positions tracked for subsequent evaluation for local phage density by DBSCAN. Phage schemes shown in Supplementary Figures 11-18 were generated with PHAST.

Statistics
Statistical analysis was conducted using R Version 3.1.3 (Team, 2012) with a suite of packages. The comparison of the correlation coefficients of the z-scores obtained for the tetranucleotide frequencies was visualized through heatmaps using rpy2 and gplots (Gautier, 2008). An exploratory analysis aiming to correlate samples with nutrients' concentrations and metabolisms, according to the level 1 SEED classification, was performed by means of a principal component analysis (PCA) using the rda function of Vegan package (Oksanen et al., 2013). Abundance plots were drawn using the ggplot2 and reshape packages (Wickham, 2007;Wickham and Chang, 2009). The cluster analysis was performed with the APE package (Paradis et al., 2004) using Pearson correlation and ward distance.

Sequence Data
The metagenomic data that we generated are available in the MG-RAST v3 server (

Results
In this study 21 metagenomic libraries were obtained from seawater (n = 12) and corals (n = 9) with a total of 446,129 sequences (2.13 × 10 8 bp) ( Table 1). Through annotation-independent analysis, samples Sw14, -15, and -22 were distinguished from -18 ( Table 1) Figure 6). The same pattern was revealed by the analysis that resulted from the tetranucleotide frequencies estimated for seawater and coral metgenomes. Samples Sw18 (HHR) were the most dissimilar amongst seawater metagenomes (Supplementary Figure 7A) and metagenomes from HHR bleached corals (MadBle18) were the most dissimilar amongst the holobionts' metagenomes (Supplementary Figure 7B). To analyze the overall relationships among the most abundant taxa we performed a clustering analysis (Figure 1). The hierarchical clustering of the 21 metagenomes corroborated the binning based in sequence composition (di-and tetranucleotides frequencies). Two major branches split seawater and coral metagenomes. Two seawater branches were defined, with samples 14 (and Sw15-1) representing the LLR. Samples 15 (Sw15-2 and 4), 18 and 22 reflected the turbulence surge (Figure 1). Coral metagenomes were split into healthy and bleached (with only two exceptions: MadBle18-1 grouped into the healthy corals branch and Mad22 into the bleached corals branch). Healthy corals also showed a trend to group according to the enrichment gradient. The healthy corals branch split Mad14 and Mad15-1 from the remainders Mad15-2, Mad18, and Mad22.

Shifts in Planktonic Assemblages
The taxonomic classification of seawater samples revealed a clear difference between LLR and the remainder groups (Sw15,. In all samples the dominant groups were Cyanobacteria (C) and Proteobacteria (P). At LLR (Sw14), these groups were equally abundant in seawater (C:P ∼1). Along the turbulence gradient, Proteobacteria members increased and Cyanobacteria FIGURE 1 | Hierarchical clustering of the metagenomes using Order level of the taxonomical identification by MG-Rast. The Pearson correlation and ward distance were used to create the dendrogram. The metagenomes from seawater, healthy and bleached M. decactis are indicated by the colors blue, brown, and red, respectively; and by the sampling day (of September/2010). Exceptions to the clustering of corals according to health status are detached (MadBle18-1 from bleached coral clustered within the healthy branch and Mad22 from healthy coral clustered within the bleached branch). The relative abundance of the 10 most abundant orders is presented at right.

Shifts in Metabolisms in Seawater Metagenomes
To investigate the metabolic profile of planktonic dwelling microbes and to correlate differences with the enrichment gradient, we performed a PCA analysis (Figures 3A,B).
Clustering of samples Sw14 (LLR) was not explained by nutrients concentrations, contrarily to the remainder samples ( Figure 3A). Basic cell functions such as Carbohydrates and Respiration explained the clustering of samples Sw14 (LLR), whether enrichment samples Sw15 and Sw18 (HHR) clustered in response to Virulence, Disease, and Defence; Membrane Transport and Nitrogen Metabolism subsystems ( Figure 3B). Samples Sw15 were grouped by the concentration of the nutrients organic phosphorous, ammonia, orthophosphate, and nitrite. The former three showed the highest concentration at Sw15 (Figure 3A, Supplementary Figure 4). HHR samples (Sw18) were coupled by the heterotrophic-characteristic metabolisms, Membrane Transport and Virulence, Disease and Defense. Accordingly, their most abundant gene was the TonB-dependent receptor (Virulence, Disease and Defense). Recovery samples (Sw22) were coupled by the concentration of nitrate (the highest, Figure 3A, Supplementary Figure 4) and Nitrogen Metabolism (Sw22-1; Figure 3B), correspondingly.

Shifts in Microbial Assemblages in Madracis Decactis Metagenomes
At HHR bleached corals (MadBle18) showed more sequences affiliated to Proteobacteria than the healthy corals (Mad18) (Supplementary Figure 19B). The presence of Pasteurellales among the ten most abundant orders in five out of nine metagenomes (in Mad15, -18, and, -22 with 11-14% of counts) appeared to be a diagnostic feature (Figure 1). Pasteurellales was found in low relative abundance (<1%) in the remainder corals and in seawater metagenomes. Thaumarchaeota affiliates in all M. decatis samples (Mad14, approximated Sw18-22 abundance levels ( Figure 2C). Cnidaria and Nematoda metagenomic sequences were more abundant in healthy and bleached corals, respectively (Supplementary Figure 8).

Shifts in the Profiles of Seawater-and Madracis Decactis-Dwelling Communities
The shift from autotroph:heterotroph-balance to offset was further investigated using virulence factors (VFs) as indicators of heterotrophy and risk or threat for corals. In total there were 21,230 significant similarities against the VFDB (Table 3). When normalized to library size, virulence genes were overrepresented in samples Sw 15,) and, to a lesser extent in Sw22 (16.6-16.8) (recovery), compared to LLR (Sw14; 12.7-14.2); and in bleached (MadBle; 0.8-8.5) when compared to healthy corals (Mad; 0.6-3.1). We further investigated the iron uptake system, which is a nonspecific virulence system related to competition skills, and thus suitable to reflect the overall heterotrophic community. Iron related virulence genes comprised ∼9% (n = 1817) of the total hits to the VFDB. Sw18 samples (HHR) presented the lower percentage of iron related genes relative to the total virulence hits per metagenome (P < 0.05), suggesting that other than iron uptake genes were most representative of the surplus heterotrophs. The VFDB lacks genes related to iron from vibrios (which was overly represented in healthy corals), but encompasses genes related to iron from Haemophilus (Pasteurellales). Among all virulence genes related to iron, ∼33% (n = 606) fell into this category. V. cholerae related virulence genes, which confer infective skills, comprised ∼6% (n = 1298) of the total virulence hits. The percentage range of these genes relative to total VFs per metagenome in seawater samples Sw14, -15, and -22 was 4.4-6.5% and in Sw18 (HHR) was 6.6-7.7%. Similarly, in healthy corals this range was 0.1-6.9%, and higher in diseased corals: 4.7-8.1% (Table 3). Heterotrophic populations that overgrew in response to turbulence-nutrient pulses were better represented by pathogenic (e.g., V. cholerae-related) than by non-specific (e.g., iron acquisition) VFs.

Discussion
Microbial assemblages during LLR (Sw14) were comparable to those previously described for the surface layers in the western SAO (South Atlantic Gyral -SATL) (Alves Junior et al., 2015) and within the WTRA (Heywood et al., 2006;Schattenhofer et al., 2009), where the dominant groups (= 50%) detected were also Prochlorales and unclassified Alphaproteobacteria or SAR11 and related. Following LLR the microbial assemblages observed increasingly differed from previous studies focusing the surface layers of those most neighboring locations (Heywood et al., 2006;Schattenhofer et al., 2009;Swan et al., 2011). Alteromonadales appeared as the second most abundant group after Prochlorales, prevailing over unclassified Alphaproteobacteria, and Vibrionales emerged as a new group with >5% relative abundance (Sw15). Comparable Alteromonadales relative abundances, combined with lower abundances of unclassified Alphaproteobacteria (SAR11) were previously reported for the sub-superficial chlorophyll maximum (SCM) layer at higher depths (48-82 m) in the SAO, where, instead of Vibrionales, Pseudomonadales, and Mamiellales emerged as differing groups compared to the surface layers (Alves Junior et al., 2015). At HHR (Sw18), other Gammaproteobacteria appeared with > 5% relative abundance, i.e., Pseudomonadales and Oceanospirillales, whereas Thaumarchaota reached Euryarchaeota relative abundance levels. Pseudomonadales, Oceanospirillales, and Thaumarchaeota characterized deep waters (236-1200 m) in the SAO, and in which water masses Prochlorales was not amongst the 10 most abundant orders (AlvesJr-14). (Schattenhofer et al., 2009) reported a Gammaproteobacteria bloom in the North Atlantic Drift Province (NADR), with a maximum relative abundance of >50% of all picoplankton in surface waters, compared to the average values of 2-5% for all the other Atlantic provinces. Only a minor fraction was identified (Alteromonas/Colwellia and Pseudoalteromonas: 2-5%; Vibrio: 1%, and Oceanospirillum: 4%). The Gammaproteobacteria bloom was attributed to the end of the spring phytoplankton bloom, indicated by declining chlorophyll values. Massive growth of Bacteroidetes was concomitant and deep water Archaea presence at surface was not observed. Gammaproteobacteria have the potential to respond to sudden nutrient pulses released from phytoplankton (Cottrell and Kirchman, 2000). Members of Alteromonas, Pseudoalteromonas, and Vibrio are well known to rapidly respond to excess nutrient supply (Bano and Hollibaugh, 2002;Beardsley et al., 2003;Allers et al., 2007Allers et al., , 2008. Thaumarchaeota are typically more abundant at depths of ≥100 m, as oposed to Euryarchaeota, which is known for decreasing abundance below 100 m (Delong, 1992;Zhang et al., 2009;Santoro et al., 2010;Tseng et al., 2015). In SPSPA, concomitant with Thaumarchaeota increase, ammonia levels decreased, possibly due to its ammonia-oxidizing ability (Francis et al., 2005); and Flavobacteriales (Bacteroidetes) relative abundance decreased. Next (Sw22), the five most abundant groups at LLR recovered relative abundances almost to LLR (Sw14) levels (Prochlorococcus, unclassified Alphaproteobacteria-SAR11, Rhodobacterales, Rhizobiales, and Chroococcales), but Alteromonadales, Pseudomonadales, Vibrionales, and Oceanospirillales remained amongst the 10 most abundant groups, as well as Thaumarchaeota remained abundant, which is a distinctive assemblage for those geographical coordinates (Schattenhofer et al., 2009;Swan et al., 2011;Alves Junior et al., 2015). Shifts in planktonic assemblages at the mesophotic zone in SPSPA were possibly driven by the turbulence surge, meaning that microbes from progressively deeper layers could hitchhike with the vertical flux along the surge. The upwelling is also supported by the enrichment and by the wide variation of water temperatures registered for the sampling depth (Supplementary Figure 4).

A Model of Physical-Chemical-Biological Dynamics in SPSPA
Although sequence similarities to genes do not represent levels of gene expression, metagenomes have been shown to be strong predictors of the biogeochemical conditions driving the microbial community (Dinsdale et al., 2008). According to the lines of evidence garnered the microbiome of the mesophotic waters in SPSPA undergoes cyclic transient shifts in relation to turbulence-nutrients regimes. A microbial succession resulting from the interplay between physical and chemical factors is a plausible scenario. Two extreme turbulencenutrient regimes can be clearly distinguished and alternate with intermediate conditions determining microbial assemblages: (i) When turbulence is low (LLR) at least 50% of the microbiome is composed of Prochlorococcus, followed by unclassified Alphaproteobacteria (SAR11 and related), which are small sized cells, highly adapted to oligotrophic conditions and starvation. Rhodobacterlaes, Rhizobiales, and Chroococcales are typical. In this environment phage genes are the most abundant in seawater, mostly from Prochlorococcus and Synechococcus, following the hosts' abundances. The viral shunt is probably less active toward relatively scarce cells. The proportion of unknown genes is the highest; (ii) Episodic surges promote vertical mixing from the immediate lower water mass to the mixed layer. Waves also wash guano from the cays flushing phosphates and ammonium into the inlet. Heterotrophs (Alteromonadales, Vibrionales) respond quickly and surpass autotrophs, motility, and chemotaxis related genes stand out; (iii) Ongoing eddies and intensified high-energy waves promote entrainment of deep water organisms such as Thaumarchaeota and nutrients (nitrite, nitrate) (HHR). Eventual cloudiness, winds and rain cope with turbulence, irradiance is intermittent and turbidity is enhanced. Heterotrophy predominates with dominance of Alteromonadales, Pseudomonadales, and Oceanospirillales. Gammaproteobacterial groups approximate 50% of the microbial assemblage, resembling the end of the spring phytoplankton bloom in higher latitudes (e.g., NADR). The gene pool in surface waters reflects the shift with membrane transport and virulencerelated genes (e.g., TonB-dependent receptor, V. cholerae virulence genes) surpassing cyanobacterial phages and basic metabolisms genes. Phages targeting heterotrophs are active. The proportion of unknown genes is the lowest; (iv) Turbulence alleviates (e.g., after moon changes toward new). Larger cell sized heterotrophs begin to decline as viral lysis and predation by grazers overrides growth, which is constrained by the paucity of limiting nutrients (e.g., phosphorus). The microbial loop is most prominent at this stage. Autotrophs respond to irradiance and retake growth (if rain, wind, and cloudiness mitigate this response is accelerated). Nitrogen metabolism is intensive. A reversal to autotrophy:heterotrophy equilibrium is triggered (Figure 4). The short-lived but recurrent turbulence-nutrient pulses might be responsible for structuring the marine ecosystem in a bottom-up manner in SPSPA. These pulses might be indispensable to warrant the energy and carbon flow to the higher trophic levels concurring to the observed pelagic fishes biomass around the barren islets (Luiz and Edwards, 2011). On a stable LLR the growth of phytoplankton is largely supported by regenerated nutrients, so only a small proportion of primary production is available to higher trophic levels or for export to the deep sea (Cullen et al., 2002;Karl, 2014). Turbulence is physically forcing the co-ocurrence of nutrients and light in SPSPA, on the other hand, the fact that nutrient resupply is shortlived might concur to retain the local mesotrophic condition. Bacterioplankton shifts were shown to be transient, following the cyclic nutrient-turbulence pulses and other physical parameters (rain, cloudiness, winds, turbidity).

The Holobiont Madracis Decactis
Some bacterial taxa prevalent in bleached corals (Rhodobacterales, Rhizobiales, and Clostridiales) have been previously associated with opportunistic diseases (Frias-Lopez et al., 2002;Rosenberg et al., 2007;Sekar et al., 2008;Sunagawa FIGURE 4 | Cartoon model of the main processes in SPSPA Cyclic turbulence-nutrient pulses determine transient shifts in the bacterioplankton. Circular central blue arrow indicates the cyclic nature of the events. Left panel: during LLR calm tides, weak winds, absence of clouds, and clear waters correlate with lower levels (downwards red arrows) of nitrogen (N) and phosphorous (P), and equivalence between autotrophy and heterotropy, represented by the ratio cyanobacteria:heterotrophs (C:H) equalizing 1. Cyanobacterial phages (green poligons) reflect the abundance of the hosts. Right panel: during HHR harsh conditions are determined by violent turbulence resulting in upwelling and enrichment, represented by high levels (upwards red arrows) of N and P, and a shift to heterotrophy dominance, represented by a low ratio cyanobacteria:heterotrophs (C:H << 1). Intense winds, cloudiness and rain, that washes the rocks covered with guano, contribute to enrichment and the shift to heterotrophy. Proteobacterial phages (red poligons) reflect the abundance of the hosts as well as Iron acquisition (red circle) and virulence factor (red crosses) genes. Vibrios abundance (blue elipses) reflect seawater parameters. Episodic surges frequently correlate with full moon. When turbulence pulses mitigate and weather assuages a recovery takes place in seawater, both in terms of nutrients concentrations and microbial assemblages (e.g., after moon changes toward new).
et al., 2009; Mouchka et al., 2010). On the other hand, a study of the corals microbiome, aiming at distinguishing the core, the symbiotic and the whole community microbiome, suggested that Rhodobacterales pertains to the latter. Conversely, Rhizobiales members were suggested to belong to the symbiotic coral microbiome. Actinomycetales and Burkholderiales (both also prevalent in all coral samples) were characterized as part of the coral core microbiome (D Ainsworth et al., 2015). Pasteurellales was one of the most abundant bacteria in the coral metagenomes, contrasting to its dwindling relative abundance in seawater.
Pasteurellales members can cause disease in a wide range of domestic and wild animals (Wilson and Ho, 2013). They are commonly found in fish tissues (Birkbeck et al., 2002). Reef fishes (Chaetodontidae) have been characterized as major vectors of coral diseases (Raymundo et al., 2009) and damselfish (Stegastes spp.) was shown to increase the prevalence of the coral Black Band Disease (BBD) (Casey et al., 2014). It is plausible that Pasteurellales are frequently transmitted to corals through fishes, possibly by fish bites, since M. decactis is frequently predated by S. sanctipauli, H. radiatus and other fishes in SPSPA (Supplementary Figure 3). This hypothesis explains the uneven distribution of Pasteurellales between healthy and diseased corals, as well as the disconnection to seawater parameters. Bleached corals were distinguished by enhanced Iron Acquisition metabolism. Six functions within this subsystem were represented only in bleached corals. HHR bleached corals samples (MadBle18) were the most dissimilar in terms of sequence composition (di-and tetranucletides frequencies), higher counts of Proteobacteria (including Vibrionales), and higher relative abundance of hits to the VFDB, including V. cholerae-related VFs. This dissimilarity, including the HHR healthy coral (Mad18), indicates that the healthy coral holobiont might be less sensitive to transient seawater-related perturbations than the diseased animals. The distinguishing characteristics of HHR bleached corals agree with the bacterioplankton and seawater features during HHR, reported both in the present and former study (Supplementary Figure 4; Moreira et al., 2014): sequence composition, higher relative abundance of motility and chemotaxis, and of membrane transport and virulence genes (e.g., Ton-B dependent receptor of the Iron Acquisition metabolism, V. cholerae-related VFs), higher vibrio counts and nutrients in seawater. Taken together, the datasets suggest coupling between the benthic and pelagic compartments, as previously reported (Chimetto Tonon et al., 2015).

Caveats
Owing to the remote nature of this site, we do not have complete data sets. Resampling will be needed to strengthen the link between turbulence-upwelling and the shifts in microbial assemblages.

Conclusions
This work analyzed shifts in microbial composition related to physical forcings (turbulence-upwelling and storms) in SPSPA. LLR is characterized by the equilibrium between autotrophyheterotrophy and microbial assemblages that resemble those of surface tropical waters previously characterized in the SAO. At HHR microbial communities shift to heterotrophic and deep-sea characteristic organisms (Thaumarchaota). HHR diseased corals are distinguished by sequence composition and enhanced VFs hits, suggesting some level of coupling between planktonic and coral microbial communities.