Excess labile carbon promotes the expression of virulence factors in coral reef bacterioplankton

Coastal pollution and algal cover are increasing on many coral reefs, resulting in higher dissolved organic carbon (DOC) concentrations. High DOC concentrations strongly affect microbial activity in reef waters and select for copiotrophic, often potentially virulent microbial populations. High DOC concentrations on coral reefs are also hypothesized to be a determinant for switching microbial lifestyles from commensal to pathogenic, thereby contributing to coral reef degradation, but evidence is missing. In this study, we conducted ex situ incubations to assess gene expression of planktonic microbial populations under elevated concentrations of naturally abundant monosaccharides (glucose, galactose, mannose, and xylose) in algal exudates and sewage inflows. We assembled 27 near-complete (>70%) microbial genomes through metagenomic sequencing and determined associated expression patterns through metatranscriptomic sequencing. Differential gene expression analysis revealed a shift in the central carbohydrate metabolism and the induction of metalloproteases, siderophores, and toxins in Alteromonas, Erythrobacter, Oceanicola, and Alcanivorax populations. Sugar-specific induction of virulence factors suggests a mechanistic link for the switch from a commensal to a pathogenic lifestyle, particularly relevant during increased algal cover and human-derived pollution on coral reefs. Although an explicit test remains to be performed, our data support the hypothesis that increased availability of specific sugars changes net microbial community activity in ways that increase the emergence and abundance of opportunistic pathogens, potentially contributing to coral reef degradation.


Introduction
Over the past decades, coral reefs have been rapidly degrading due to a combination of local and global factors that include overfishing (Jackson et al., 2001), pollution (McCook, 1999), and climate change (Hughes et al., 2003;Hoegh-Guldberg et al., 2007), leading to an increase in coral bleaching and disease emergence (Weil et al., 2006;Baker et al., 2008). Although diseases represent an important cause of deterioration in these productive tropical marine communities, they are still poorly understood (Willis et al., 2004;Weil et al., 2006;Muller and Woesik, 2012). Paleontological and ecological monitoring suggest an increase in the prevalence of reef diseases over the last few decades (Harvell et al., 1999;Porter et al., 2001;Cramer et al., 2012), partially as a consequence of nutrient enrichment (Bruno et al., 2003;Voss and Richardson, 2006;Vega Thurber et al., 2014), and particularly because of increases in dissolved organic carbon (DOC) (Kuntz et al., 2005;Kline et al., 2006;Smith et al., 2006). Autochthonous DOC can enter coral reefs by coral and algal exudation (up to 40% of algal fixed carbon produced in photosynthesis) in the form of a heterogeneous mixture of saccharides, proteins, and lipids (Haas and Wild, 2010;Nelson et al., 2013). Although exudates released by scleractinian corals contain high levels of refractory DOC and are therefore more similar in their sugar composition to reef waters and offshore DOC pools, turf and fleshy macroalgalderived DOC is highly enriched in the labile sugars glucose, galactose, mannose, and xylose . Therefore, algae exudates are composed similarly to sewage discharge (Huang et al., 2010;Wear and Vega Thurber, 2015). Owing to these qualitative differences, coral and macroalgal exudates differentially affect the composition and metabolism of coral reef seawater microbial communities Haas et al., 2016;Lee et al., 2016). The detrimental effects of labile DOC on coral reef functioning are summarized in the DDAM (Disease, Dissolved organic carbon, Algae and Microbe) model (Barott and Rohwer, 2012).
The DDAM model proposes a positive feedback loop in which human-derived stressors, such as overfishing and eutrophication, trigger macroalgal growth and promote the release of DOC-enriched exudates . This mostly labile DOC fosters bacterial growth and oxygen removal, selecting for copiotrophs and opportunistic pathogens both in the seawater and the coral holobiont (Vega Thurber et al., 2009;Haas et al., 2011;Nelson et al., 2013;Haas et al., 2016). This selection for potential pathogens in turn completes the loop as it promotes coral bleaching, disease, and eventually mortality, providing space for algae to grow and become dominant on the reef. The etiology of many coral diseases, however, is poorly understood. Although numerous coral diseases have been described, causative agents have been identified only in a few cases and with incomplete satisfaction of Koch's postulates (Richardson et al., 1998;Alker et al., 2001;Sutherland et al., 2011). Efforts towards the characterization of the whole microbial community rather than looking for single agents have made significant advances in the field. For instance, the appearance of signs of coral disease often correlates with shifts in coral-associated bacterial community composition (Sunagawa et al., 2009;Cárdenas et al., 2012;Roder et al., 2014), suggesting a more complex pathogenesis. To date, these diseases are believed to be opportunistic infections triggered by the exposure to environmental stressors (for example, elevated temperature or nutrient enrichment), reducing host resistance and promoting the uncontrolled growth of opportunistic pathogens (Landsberg, 1995;Rosenberg et al., 2007). These opportunistic pathogens can originate either from the coral holobiont itself (symbionts) or the water column (environmentally acquired) and are also associated with healthy reefs, where they fulfill key functions to support the ecosystem, for example, antibiotic production and nitrogen fixation (Ritchie, 2006;Chimetto et al., 2008;Mohamed et al., 2008).
Virulence can be understood as harm or morbidity caused by a microbe to its host. Thus, identifying virulence determinants or factors (VF) can provide insights into the mechanistic basis of such harm (Brown et al., 2012). VFs are molecules produced by pathogens that contribute to pathogenicity and typically include toxins, exoenzymes, adhesins, and secretion systems (Brown et al., 2012). The expression of VFs is commonly correlated with the ability of a pathogen to invade and exploit host tissues. Nelson et al. (2013) showed that algal exudates support the growth of bacterioplankton populations among the families Vibrionaceae and Pseudoalteromonadaceae harboring high numbers of VFs. Temperature-dependent expression of virulence factors has been reported in several Vibrio species (Banin et al., 2003;Kimes et al., 2012). Despite these studies, there has been no compelling evidence correlating other environmental factors with the activation of VFs in opportunistic coral reef pathogens.
Expression of VFs could provide a selective advantage for opportunistic microbes to gain access to nutrients. Therefore, these genes are often linked to general pathways of nutrient uptake (Görke and Stülke, 2008). As an example, so-called moonlighting proteins are involved in major metabolic processes such as the glycolytic pathway and stress responses that have unexpected functions which contribute to microbial virulence (Henderson and Martin, 2011). Carbon catabolism is linked to bacterial virulence by the sensing of sugars through phosphotransferase system transporters, which in turn activate carbon catabolite repression (Deutscher et al., 2006;Görke and Stülke, 2008). When different carbon sources are available, carbon catabolite repression allows bacteria to control the uptake of the preferred carbon substrate by disabling genes involved in the use of secondary carbon sources (Brückner and Titgemeyer, 2002). Carbon catabolite repression is observed in most free-living heterotrophic bacteria and acts as a global regulator, controlling up to 10% of all bacterial genes, including several virulence factors (Yoshida et al., 2001;Blencke et al., 2003;Görke and Stülke, 2008). The combination of carbon catabolite repression and multiple transcriptional regulation networks allows copiotrophic bacteria to succeed at high-carbon concentrations (Cottrell and Kirchman, 2016).
A prediction of our pathogen emergence model is that elevation of specific sugars will change the microbial community composition and microbial gene expression, thereby promoting the risk of reef degradation. Therefore, in the present work, we characterized the response of microbial communities from coral reef seawater samples before and after the addition of the dominant monosaccharides previously reported for algal exudates and sewage discharge (glucose, galactose, mannose, xylose). This was achieved in an ex situ approach by amending natural coral reef seawater samples in single monosaccharide incubations in the dark in 4 L microcosms over 48 h. Microbial communities were characterized via 16 S rRNA gene amplicon sequencing and reef planktonic population genomes were recovered via metagenomic binning. Microbial mRNA reads were subsequently mapped to the recovered population genomes to assess the activity of distinct microbial community members under sugar enrichment. This allowed us to elucidate potential mechanisms that suggest negative interactions between planktonic microbes and corals driven by high sugar concentrations.

Water sampling and incubations
Incubations consisted of three liters of 0.22 μmfiltered seawater inoculated with one liter of unfiltered reef seawater in a four-liter polycarbonate bottle (sulfuric acid-cleaned and seawater-leached). Unfiltered seawater was collected at 10 m depth from Al-Fahal reef located 13 km offshore in the central Red Sea of Saudi Arabia at 9:30 am on 20 October 2014 and 03 November 2014 for the first and second set of experiments, respectively. Every experiment consisted of four triplicate treatments, two sugaramended and two unamended controls. Experiment 1 consisted of (B1) unamended control before incubation, (C1) unamended control after incubation, (Glu) after glucose incubation and (Gal) after galactose incubation. Experiment 2 consisted of (B2) unamended control before incubation, (C2) unamended control after incubation, (Man) after mannose incubation and (Xyl) after xylose incubation. Each treatment was amended with 500 μM of the corresponding sugars. The concentration was chosen to elicit rapid microbial responses. The incubation time was 48 h at 27°C in the dark to avoid phytoplankton growth. Supplementary Figure 1 shows a schematic representation of the different experiments and enrichments.
Cell density, DOC concentrations, and efficiency calculations Changes in cell abundance and DOC concentrations were measured over time during sugar incubations. Bacterial cell abundance was determined by flow cytometry using the protocol established by Nelson et al. (2011). In brief, unfiltered samples were fixed with paraformaldehyde to a final concentration of 0.4% and frozen (−80°C) within 30 min of fixation. Fixed samples were thawed, mixed, and stained with 1 × SYBR Green for 1 h at room temperature in the dark. Cell counts were done using the Guava easy-Cyte flow cytometer (Millipore, Billerica, MA, USA). To measure DOC concentrations, triplicates of 20 ml samples were filtered through 0.45 μm pore size Minisart-GF filters (Sartorius, Gottingen, Germany). The filtrate was collected in 25 ml pre-rinsed polyethylene HDPE bottles. A total of 100 μl of a 35% H 3 PO 4 solution were added to acidify the samples (pH below 2) and subsequently stored at − 20°C in the dark until processing. DOC concentrations were measured in triplicates using the Apollo 9000 total organic carbon and Total Nitrogen Analyzer (Teledyne Tekmar, United States) at the Analytical Core Lab at KAUST, Saudi Arabia. Bacterial carbon change was determined as the difference in cell carbon (μmol C l − 1 ) before and after the incubation using a previously reported factor for coastal bacteria assemblages of 30 fg C per cell (Fukuda et al., 1998). Bacteria growth efficiency was calculated as the ratio of bacterial carbon production (rate of increase in bacterial carbon) to the rate of DOC removal (Eiler et al., 2003).

Nucleic acid extraction and isolation
Triplicate water samples were pre-filtered subsequently through 10 and 3 μm polycarbonate filters and collected onto 0.22 μm polycarbonate filters. Filters were immediately frozen at − 80°C until nucleic acid extraction. Replicate filters were pooled for nucleic acid extractions using the AllPrep DNA/ RNA kit (Qiagen, Hilden, Germany). Purified DNA was used for rRNA gene amplification and preparation of metagenomic sequencing libraries (see below).
About 200 ng of total RNA were amplified using the MessageAmp II-Bacteria kit (Ambion, Austin, TX, USA) as described previously (Frias-Lopez et al., 2008;Shi et al., 2011). Prokaryotic mRNA was purified according to Stewart et al. (2010) with slight modifications (Daniels et al., 2015). In brief, rRNA was subtracted from the total RNA using anti-sense probes created by the in vitro transcription of ribosomal genes (bacterial and archaeal 16S/23S and eukaryotic18S/28S) at a final template-to-probe ratio of 1:2 (mass, per probe). Subsequently, mRNA was subjected to an oligodT-based separation to exclude polyA+ eukaryotic mRNA from the polyAprokaryotic mRNA. Metatranscriptomic TrueSeq libraries were sequenced using the Illumina MiSeq V3 system (Illumina, San Diego, CA, USA) with paired-end reads (2 × 300 bp). Metagenomic TrueSeq libraries were sequenced using the Illumina HiSeq 2000 system (Illumina) with paired-end reads (2 × 150 bp) according to the manufacturer's specifications.

16S-based diversity analysis
Microbial community composition was assessed via 16S rRNA sequencing before and after sugar enrichments. For this, 16S rRNA V5-V6 hypervariable regions were PCR-amplified for 25 cycles using 16S primers 784 F (5′-AGGATTAGATACCCTGGTA-3′) and 1061 R (5′-CRRCACGAGCTGACGAC-3′) designed by Andersson et al. (2008), which work well in marine environments (Bayer et al., 2013). PCR products were sequenced using Illumina MiSeq V3 (2 × 300 bp) according to the manufacturer's specifications at the Bioscience Core Lab (KAUST, Saudi Arabia). Primer trimming, quality control, clustering and taxonomic classification were done in mothur v1.36.1 (Schloss et al., 2009). Quality control consisted of removing sequences shorter than 300 bp, and/or with quality scores lower than 25 and those that occurred only once (singletons) across the entire dataset. Sequences were aligned to SILVA reference, release 119 (Quast et al., 2013), and clustered into operational taxonomic units (OTUs), defined at 97% similarity. Representative sequences from each cluster were taxonomically assigned using Virulence factor induction in coral reef bacteria A Cárdenas et al Greengenes reference taxonomy (McDonald et al., 2012). Chimeric sequences were identified with the de novo implementation of UCHIME in the mothur interface (Edgar et al., 2011) and subsequently removed. Principal Coordinate Analysis using Bray-Curtis distance and diversity analyses (Shannon-Weaver and Chao1 indexes) was performed on operational taxonomic unit abundance data using the Phyloseq (McMurdie and Holmes, 2013) R package, using output tables generated in mothur.

Metagenome binning
Sequence reads were quality checked and trimmed for low-quality regions, adapter sequences, and a minimum length of 75 bp using Trimmomatic v0.36 (Bolger et al., 2014). MEGAHIT v1.0.2 (Li et al., 2015) was used to produce one assembly per experiment using the 'meta-sensitive' option and kmer length of 81. Metagenome binning was done based on differential coverage and tetranucleotide signatures for each assembly separately using the binning tools GroopM v0.3.4 (Imelfort et al., 2014) and Metabat v0.25.4 (Kang et al., 2015). Genome contamination and completeness were assessed using CheckM v1.0.3 (Parks et al., 2015). To increase genome completeness, the 'merge' option from CheckM was used and population genomes with a delta contamination lower than 3 and a merged contamination lower than 5 were merged. CheckM 'unique' and 'join' commands were used to compare and combine genomes generated from the different binning tools. The metagenomic binning procedure resulted in population genomes representing putative microbial taxa. A first taxonomic approximation was done based on universal, single-copy phylogenetic marker genes using specI (Mende et al., 2013), marker lineages from CheckM, and sequence compositionbased classifier (Patil et al., 2012). Whole-genome based average nucleotide identity and amino-acid identities between the query genomes and NCBI reference genomes were used to validate closest related taxa using JSpecies v1.2.1. (Richter and Rosselló-Móra, 2009) and the enveomics tools (Rodriguez-R and Konstantinidis, 2016). Open reading frames were predicted using Prodigal (Hyatt et al., 2010). Predicted protein-encoding genes were functionally annotated using Gene Ontology (GO) and myRAST server (Aziz et al., 2008) according to SEED subsystems (Overbeek et al., 2014). Specific processes were inferred from GO annotation as listed in the supplement (Supplementary Table S2). For instance, microbial interaction genes were identified based on the PAMGO (Plant-Associated Microbe Gene Ontology) Consortium (Torto-Alalibo et al., 2009). Virulence genes were assigned by a BLASTn v2.2.28 search against MvirDB (Zhou et al., 2007) and VFDB (Chen et al., 2005) databases (both downloaded on 25 March 16), supplemented with GO terms related to pathogenesis (Supplementary  Table S2). Sugar transporters were identified using the Carbohydrate Active Enzymes database (http://www.cazy.org/) (Cantarel et al., 2009). Genes involved in antibiotic resistance were deducted from the BLAST search against the Antibiotic Resistance Genes Database (ARDB) (downloaded on 07 July 16) (Liu and Pop, 2009).

RNA mapping and gene expression analysis
Remaining rRNA in metatranscriptomic data were removed in silico using SortMeRNA v2.0 (Kopylova et al., 2012). Cleaned mRNA fragments were mapped to a multifasta file containing all annotated genomes using the Bowtie aligner v1.0.0 using the 'best' option allowing a single hit (Langmead and Salzberg, 2012). Transcript-level estimates were obtained as mRNA read counts using the transcript abundance estimation tool eXpress v1.5.1 (Roberts and Pachter, 2013) and normalized by the size of the RNA library (sequencing depth) and the coverage of the population genome in the metagenome (population genome abundance). Differential expression, expressed as log2 fold change, was calculated for mRNA reads in each treatment compared with the controls (that is, C1 and C2 samples). A twofold fold change cutoff was used to determine up-and downregulated genes. Upregulated genes (that is, induced genes) were used to establish changes in biological processes after sugar addition.

Results and discussion
The response of reef seawater microbial communities to labile DOC enrichment was investigated to understand the complex interactions and mechanisms involved in the emergence of pathogenicity in these communities. The monosaccharide additions to ex situ microcosms elicited differential, sugarspecific responses reflected in rapid microbial community shifts, gene expression, and associated virulence potential.

Microbial diversity and DOC consumption
The main differences in operational taxonomic unit (OTU) composition between both controls (natural reef water, no sugar added, not incubated) were a higher relative abundance of Alteromonadaceae and lower Oceanospirillaceae in experiments 1 compared with experiment 2. These differences were likely due to small changes in inorganic nutrient concentrations driven by seasonality (Roik et al., 2016) as well as by differences in DOC concentrations in the natural reef water before the start of the Virulence factor induction in coral reef bacteria A Cárdenas et al experiments (Table 1). After 48 h incubations, sugarenriched samples deviated from the controls in a sugar-specific manner, generally towards a less-rich and less-diverse bacterial community (Supplementary Figure S2A and B). Microbial community shifts were characterized by increases in Rhodobacteraceae (Alphaproteobacteria) in galactose, mannose, and xylose, Vibrionaceae (Gammaproteobacteria) in glucose, and Alteromonadaceae (Gammaproteobacteria) in both glucose and mannose enriched samples ( Figure 1). A recent study supporting the microbialization concept showed that oligotrophic coral-dominated waters are characterized by high abundance of Alphaproteobacteria, in contrast to enriched algal-dominated waters where Gammaproteobacteria are abundant (Haas et al., 2016). Linking our study to Haas et al. (2016), our results suggest glucose and mannose are major contributors to bacterial shifts in impacted reefs. Furthermore, microbial communities had significantly faster and more efficient growth in galactose and glucose enrichments (Table 1) suggesting that these sugars play a major role in supporting an increase in microbial biomass. However, efficiency estimates can be influenced by the conversion factor for coastal bacteria assemblages of 30 fg C per cell used here. This is only a rough estimate, and human impact and increased algal cover can also significantly influence cell size (McDole et al., 2012).

Metagenomic binning and microbial population genomes
We assembled metagenomic reads from each experiment and calculated the abundance profile and tetranucleotide frequencies of each contig between treatments. This allowed us to assign co-abundant contigs with similar tetranucleotide frequencies into groups that represent the gene content of distinct microbial population genomes representing putative microbial taxa. We recovered 4400 draft population genomes and only used those with total completeness above 70% for further analyses (Table 2). Annotation and gene expression analyses were conducted for the selected 27 population genomes, after determining nucleic acid and amino-acid identities of complete reference genomes to infer closest related organisms (Supplementary Table S3). Population genomes comprised a wide diversity of taxa including different populations of the copiotrophic Alteromonas and several members of the Rhodobacteraceae, such as Oceanicola, Ruegeria, and Labrenzia (Table 2).
To profile transcriptomic patterns in planktonic microbial populations, bacterial mRNA sequences were mapped to the recovered population genomes. Mapping counts were normalized to the genome abundance in the overall population, thereby uncovering gene expression at an individual rather than population level. Furthermore, the results of the present work are focused on differentially upregulated genes (that is, induced genes) after comparing expression levels in the control versus each treatment. Populations with high numbers of differentially expressed genes were identified using RAST annotation (Aziz et al., 2008) of induced biological processes (Figure 2b and Supplementary Figure S3). These population genomes included members of the Gammaproteobacteria genera Alteromonas, Pseudoalteromonas, Alcanivorax, and Hahella, and of the Alphaproteobacteria genera Oceanicola and Erythrobacter. Most of these are comprised of microbes typically found in coral reefs (Frias-Lopez et al., 2002;Bourne and Munn, 2005;Tout et al., 2015), but are often linked to nutrient-enriched and degraded reefs including fish farm effluents (Garren et al., 2008), oil-contaminated reefs (Al-Dahash and Mahmoud, 2013), and coral diseases (Sunagawa et al., 2009;Cárdenas et al., 2012;Roder et al., 2014).

Sugar transporter diversity and microbial metabolism
The recovery of near-complete population genomes allowed us to investigate cellular mechanisms activated by the sensing and uptake of dissolved monosaccharides. The first step in understanding sugar uptake was to evaluate the diversity of sugar transporters. The most common transporters among the population genomes were galactose ABC transporters, mannose phosphotransferase system, and Virulence factor induction in coral reef bacteria A Cárdenas et al glucose/mannose major facilitator superfamily porters. The mannose phosphotransferase system is not only an important uptake system for mannose and glucose, but it has a major role in the transcriptional regulation and the preferential utilization of monosaccharides (Postma et al., 1993;Vu-Khac and Miller, 2009). The phosphotransferase system transporters were highly induced in population genomes with higher numbers of differentially expressed genes, specifically genes involved in symbiotic interactions and virulence factors (discussed below) (Figures 2a and b). Copiotrophs were expected to possess a wide array of low-affinity and highly specific transporters, contrary to oligotrophs that have a smaller number of broad-specificity and sufficiently high-affinity transporters (Lauro et al., 2009). This was confirmed in the most abundant population genomes (Figure 2c), such as GM_23 Virulence factor induction in coral reef bacteria A Cárdenas et al (Rhodobacteraceae), which possessed five different transporters for sugar uptake (Figure 2a). Most of these transporters belonged to the low-affinity ABC transporters and the multi-sugar major facilitator superfamily, whereas high-affinity transporters were rare. The presence of several glucose uptake systems may allow for flexibility under different environmental conditions. For instance, low-affinity systems are expected to be induced in sugar-rich environments. Contrastingly, high-affinity transporters are used when low amounts of glucose are available (Jahreis et al., 2008). Using metagenomics, Haas et al. (2016) found that coral-dominated reef waters were enriched for genes in the energy efficient EMP (Embden-Meyerhof-Parnas) pathway, whereas algal-dominated reefs waters were enriched with genes in the lessefficient ED (Entner-Doudoroff) and PP (pentose phosphate) pathways. To test whether we saw a similar pattern, we evaluated the expression of genes involved in the EMP, ED, and PP pathways in the population genomes (Figure 3). We chose to focus on the pathway with highest expression levels in each of the population genomes. Although the response was sugar-and population-specific, some general patterns could be established. For instance, population genomes MB_8 (Euryarchaeota), merged_809_810 (Alteromonas), and MB_4b (Alphaproteobacteria) had highest expression levels of the EMP genes. However, genes involved in the ED and PP pathways were induced in population genomes GM_189 (Alteromonas), MB_4a (Erythrobacter), GM_554 (Alcanivorax), GM_293 (Alphaproteobacteria) and MB_66 (Erythrobacter). For the latter population genomes, we could corroborate previous findings that support shifts from EMP to ED and PP. These shifts in the metabolic capacities were suggested to represent a less costly, faster, but lessefficient transformation of organic carbon under high-carbon availability (Flamholz et al., 2013;Haas et al., 2016).

Induction of microbial interaction genes
In general, microbial symbioses can be beneficial, neutral, or pathogenic (Hillman and Goodrich-Blair, 2016). Microbial interaction genes comprise several mechanisms by which microbes colonize, grow in, and occasionally cause detrimental effects to hosts and host tissues (Dale and Moran, 2006). Interaction genes were classified into distinct categories including biofilm formation, adhesion, secretion systems, and genes involved in pathogenesis (among others), according to the PAMGO consortium (Torto-Alalibo et al., 2009). Numerous interaction genes identified in the population genomes belonged to antibiotic transport and production (Supplementary Figure S4). Antibiotic-producing bacteria often associate with marine organisms (Kennedy et al., 2009;Nissimov et al., 2009;Wiese et al., 2009), and  Fold changes o2 (that is, not induced) are denoted in black thin lines. Thicker lines represent fold changes42 (that is, induced) after enriching with sugars colored in blue for galactose, magenta for glucose, yellow for mannose, and green for xylose. (b) Heat map depicting the number of induced genes across the population genomes in biological processes using SEED subsystems annotation. Bacterial population genomes response is clustered by Euclidian distance for galactose and mannose. (c) Relative abundance of each microbial population genome in the metagenome, calculated as the proportion of metagenomic reads mapped to the population genome from the total number of reads in the metagenomic library.
Virulence factor induction in coral reef bacteria A Cárdenas et al Virulence factor induction in coral reef bacteria A Cárdenas et al Figure 4 Interaction and virulence factor genes in microbial population genomes. (a) Change in population abundance between control and treatments (metagenome, y axis) in relation to the differentially expressed bacterial interaction genes (metatranscriptome, x axis). Increasing numbers on the y axis represent organisms with higher changes in abundance after sugar addition, whereas increasing numbers on the x axis denote bacterial species with a high number of expressed interaction genes. Population genomes that are more likely to initiate host interactions (that is, bacteria that grow rapidly under sugar enrichment, expressing a high number of interaction genes) are highlighted in a gray area. (b) Black bars represent the total number of interaction genes (left) and virulence factors (right) present in each population genome. Colored bars represent the number of differentially expressed genes under either galactose (blue), glucose (magenta), mannose (yellow), or xylose (green) enrichments. Stars denote population genomes with the largest number of virulence factors and correspond to the selected population genomes from (A).
Virulence factor induction in coral reef bacteria A Cárdenas et al consequently, coral mucus-associated bacteria are considered the first line of defense (Ritchie, 2006;Shnit-Orland and Kushmaro, 2009). Contrastingly, free-living bacteria are irregular antibiotic producers (Nair and Simidu, 1987;Long and Azam, 2001). Therefore, exhibiting this trait could represent an advantageous strategy to colonize and effectively outcompete native benthic bacteria in new habitats. Another important group of interaction genes, 'viral processes', included viral-host interactions (Supplementary Figure S4). The latter suggests that viral infections are also playing major roles driving bacterioplankton community diversity in this study, as in other marine ecosystems (Middelboe et al., 2001;Weinbauer and Rassoulzadegan, 2004;Holmfeldt et al., 2007;Vega Thurber et al., 2017). A large number of interaction genes were categorized as virulence factors owing to their high homology with previously reported virulence factors in other microbial species. Virulence factors Figure 5 Expression of genes involved in pathogenesis in potential opportunistic pathogens (POPs). Each diagram represents a bacterial cell and each side of the cell membrane represents a group of cellular processes. Antibiotic resistance on the upper side of the membrane, environmental sensing on the right side of the membrane, bacterial secretion systems and adhesion on the lower part, and extracellular proteins and toxic compounds on the left part.
Virulence factor induction in coral reef bacteria A Cárdenas et al represent the subset of interaction genes that include all genes involved in negative microbial interactions according to PAMGO (Torto-Alalibo et al., 2009) and MvirDB (Zhou et al., 2007). We classified the population genomes based on their abundance in the metagenome and the highest number of induced interaction genes and virulence factors into the following microbial lifestyles: (1) the larger the increase in abundance after sugar enrichment, the more likely to have a copiotrophic lifestyle (Figure 4a), (2) the larger the number of interaction genes, the more likely this population genome represents an opportunistic symbiont (Figures 4a  and b), (3) the larger the number of virulence factors, the more likely this population genome represents an opportunistic pathogen (Figure 4b). Based on these criteria, we chose the following six opportunistic pathogens with copiotrophic lifestyles and defined them as Potential Opportunistic Pathogens (POPs): MB_1 (Oceanicola sp.), GM_189 (Alteromonas), GM_439 (Alteromonas), GM_554 (Alcanivorax), MB_4a (Erythrobacter), and GM_66 (Erythrobacter).

Virulence potential in opportunistic pathogens
Much has been discussed regarding opportunistic pathogens in the field of coral disease and the responsible microbes, but little is known about the underlying mechanisms of disease emergence. The POPs we describe above represent candidate bacterial taxa, as they are found commonly associated with healthy organisms and ecosystems as well as disturbed ecosystems and diseased marine organisms. For instance, Erythrobacter and Alteromonas are numerously found in association with healthy marine organisms (Shigemori et al., 1992;Yoon et al., 2004;Koren and Rosenberg, 2006;Lampert et al., 2006;Chiu et al., 2007;Kalimutho et al., 2007;Martínez-García et al., 2007;Vandecandelaere et al., 2008;Reis et al., 2009). However, their participation in coral and other marine animal diseases is strongly suggested (Cipriani et al., 1980;Garland et al., 1983;Pantos et al., 2003;Garren et al., 2008;Sunagawa et al., 2009;Cárdenas et al., 2012;Pollock et al., 2016). In this regard, we highlight the importance of not only correlating the presence of certain bacteria with a disease, but rather looking at possible biotic and abiotic factors turning commensal bacteria into pathogens. Furthermore, the observed relative abundances of these POPs never exceeded 4% of the total community (Supplementary Table S1). This highlights the role of rare bacteria in natural communities and the need to rethink our ideas of how abundant and rare microbes contribute to ecosystem processes (Jousset et al., 2017). In Figure 5, we show genes previously linked with pathogenesis in each POP. Based on each pool of virulence factors, in the following we discuss potential strategies used by the POPs to colonize and effectively outcompete native benthic bacteria on coral surfaces.
Expression of interaction genes to promote virulence potential A common feature among many pathogens is their ability to adhere to surfaces, colonize them, and survive adverse conditions in the form of biofilms (Dunne, 2002;Pizarro-Cerdá and Cossart, 2006). As evidence of these processes, genes involved in adhesion, such as the Type IV secretion system and the translocation and assembly module, were induced under sugar enrichment in all population genomes ( Figure 5). Membrane proteins tamA and B help assemble surface structures essential to hostpathogen interactions, including adhesion and host invasion (Selkrig et al., 2012), and deletion of these elements diminishes the virulence of several pathogens (Struve et al., 2003;Burall et al., 2004;Kelly et al., 2006). Furthermore, GM_189 and GM_439 (Alteromonas sp.) induced eps and ams genes ( Figure 5), which code for exopolysaccharide production and export proteins, that are important for biofilm formation under sugar enrichment. Biofilms not only constitute a protected lifestyle for bacteria and promote pathogenesis, but also cause impacts from an ecological perspective in coral reef ecosystems. For instance, marine biofilms containing Alteromonas species can influence benthic community succession by inducing the settlement and metamorphosis of cnidarian larvae and algal zoospores (Leitz, 1997;Graham et al., 2000) while inhibiting the settlement of polychaetes (Dobretsov and Qian, 2002).
Expression of genes related to iron uptake to promote virulence potential Another group of virulence factors induced after sugar enrichment is related to iron uptake. Iron is limited in the ocean and at the same time is an essential micronutrient. Iron transporters and siderophores (iron chelating compounds) are common strategies among bacteria to acquire iron. Extracellular siderophores are rarely produced by pelagic bacteria, but are more common on particle-attached and benthic marine microorganisms (Hutchins et al., 1999;Hopkinson and Barbeau, 2012). Some pathogens can produce siderophores to sequester iron from the host or other living organisms and use it to support their own growth (Amin et al., 2009). The most common genes related to iron uptake processes found in the POPs were the TonB-dependent transporters and the outer protein IrgA, expressed in the population genomes GM_189 (Alteromonas sp.), MB_4a (Erythrobacter), GM_439 (Alteromonas sp.), and GM_554 (Alcanivorax sp.) ( Figure 5). TonBdependent iron acquisition is one of the main mechanisms permitting bacterial growth in a wide range of iron-limited environments, including host colonization in several pathogens (Goldberg et al., 1992;Torres et al., 2001;Leduc et al., 2008). IrgA is the outer membrane receptor of vibriobactin, a siderophore produced by Vibrio species and extracellularly transported by TonB (Mey et al., 2002;Kustusch et al., 2011). The induction of multiple TonB elements and siderophore receptors in GM_189 and GM_439 (Alteromonas sp.) strongly suggests iron scavenging as an important mechanism used by these POPs for successful growth. However, little is known about marine siderophore substrates and mechanisms by which they contribute to virulence in marine pathogens.
Expression of toxins and proteases to promote virulence potential The third group of sugar-induced virulence factors found in the POPs was 'toxins and proteases'. An example of a non-proteinaceous endotoxin was lipopolysaccharide, evidenced in the population genomes GM_189 (Alteromonas sp.) and GM_554 (Alcanivorax sp.) by the induction of the lipopolysaccharide biosynthetic genes galU, algG, algA ( Figure 5). Lipopolysaccharide is involved in the modulation of marine symbiosis (Foster et al., 2000), but also has an important role in the pathogenesis of several bacterial species, including Vibrio spp. (Nesper et al., 2001), by eliciting strong immune responses in host cells. Among induced extracellular toxins, membrane-disrupting and proteolytic toxins were prevalent (Supplementary Table S4). Two extracellular zinc metalloproteases were induced by the POPs GM_189 (Alteromonas sp.) and MB_1 (Oceanicola sp.) ( Figure 5). Similar extracellular zinc metalloproteases are widely distributed among several marine bacteria and largely in Vibrio species (Denkin and Nelson, 2004;Hasegawa et al., 2008;Labreuche et al., 2010). Only a few mechanisms of these metalloproteases are well described, as in the case of the hemagglutinin produced by Vibrio cholerae to degrade intestinal protective mucus in humans (Benitez and Silva, 2016). In addition, we provide evidence for the induction of subtilisin-like serine proteases in GM_554 (Alcanivorax sp.) that have antifouling activity via degradation of exopolysaccharide biofilms (Leroy et al., 2008). The expression of this protease may negatively regulate the bacterium's own biofilm or compete against other bacteria. Among membrane-disrupting toxins, the cytotoxins TlyA and TlyC were induced in the POPs GM_66 (Erythrobacter) and MB_1 (Oceanicola sp.) ( Figure 5). Homologs of these toxins have the ability to lyse erythrocytes, leukocytes, phagosomal vacuoles, and gastric epithelium cells (Martino et al., 2001;Whitworth et al., 2005). Another group of membrane-disrupting toxins induced by sugar enrichment was the calcium-dependent pore-forming cytotoxins (RTX toxins). Three RTX homologs were expressed in the POP MB_1 (Oceanicola sp.) and two elements involved in RTX toxin translocation were expressed in the POPs GM_189 and GM_439 (Alteromonas sp.) ( Figure 5). RTX proteins are potent virulence factors with cytotoxic and hemolytic activities toward a broad range of animal cells (Lally et al., 1999), but there is little evidence to date about the role of RTX toxins in disease. Among marine bacterial species, RTX toxins have been found in several Vibrio species (Kim et al., 2008;Weynberg et al., 2015), Pseudoalteromonas agarivorans (Choudhury et al., 2015) and in members of the Roseobacter clade, including Phaeobacter (Durighello et al., 2014) and Ruegeria pomeroyi (Christie-Oleza et al., 2012). These proteins are characterized by tandemly repeated nonapeptides with the consensus sequence GGXGXDX[L/I/V/W/Y/ F]X (where X is any amino acid) in the C-terminal half of the protein, as observed in RTX homologs expressed by MB_1 (Oceanicola sp.) (Supplementary Figure S5). Expression of toxin-related proteins is consistent with expression of type I to IV secretion systems that are involved in the transport of several virulence factors including effector proteins and toxins (Backert and Meyer, 2006). For instance, one main feature of all RTX toxins is their conserved mechanism of transport mediated by the type 1 secretion system (Hueck, 1998;Linhartová et al., 2010). Elements of the type 1 secretion system were induced in POPs expressing RTX-related homologs, suggesting the export of these protein homologs during sugar incubations.

Conclusions
Anthropogenic pollution and algal cover are increasingly contributing to the degradation of coral reefs worldwide by the increase of DOC concentrations. By assaying bacterial community composition and gene expression under experimental enrichment of sugars abundant in algal exudates and sewage inflows, we show that high DOC concentrations influence microbial diversity and activity, supporting a microbial lifestyle switch from commensal to pathogenic. In particular, bacteria in the genus Alteromonas, Oceanicola, Erythrobacter, and Alcanivorax provide candidate taxa for this switch, based on their copiothrophic behavior and their potential to become harmful to other marine organisms under high sugar concentrations. These POPs displayed a shift in their metabolic capabilities and expression of a wide set of virulence factors, including metalloproteases, siderophores, and toxins, as well as numerous mechanisms of antibiotic resistance. For the first time, this study shows DOCdependent expression of VFs and shifts in metabolic capacities of coral reef microbial populations within the context of coral reef health. Although the present study has its limitations and requires validation by an in situ approach, our findings corroborate the DDAM and reef microbialization hypothesis, providing new mechanistic insight into the link between elevated DOC concentrations and the switch from free-living to attached stages, as well as from commensal to pathogenic lifestyles.