Bacterial Community Assembly, Succession, and Metabolic Function during Outdoor Cultivation of Microchloropsis salina

ABSTRACT Outdoor cultivation of microalgae has promising potential for renewable bioenergy, but there is a knowledge gap on the structure and function of the algal microbiome that coinhabits these ecosystems. Here, we describe the assembly mechanisms, taxonomic structure, and metabolic potential of bacteria associated with Microchloropsis salina cultivated outdoors. Open mesocosms were inoculated with algal cultures that were either free of bacteria or coincubated with one of two different strains of alga-associated bacteria and were sampled across five time points taken over multiple harvesting rounds of a 40-day experiment. Using quantitative analyses of metagenome-assembled genomes (MAGs), we tracked bacterial community compositional abundance and taxon-specific functional capacity involved in algal-bacterial interactions. One of the inoculated bacteria (Alteromonas sp.) persisted and dispersed across mesocosms, whereas the other inoculated strain (Phaeobacter gallaeciensis) disappeared by day 17 while a taxonomically similar but functionally distinct Phaeobacter strain became established. The inoculated strains were less abundant than 6 numerically dominant newly recruited taxa with functional capacities for mutualistic or saprophytic lifestyles, suggesting a generalist approach to persistence. This includes a highly abundant unclassified Rhodobacteraceae species that fluctuated between 25% and 77% of the total community. Overall, we did not find evidence for priority effects exerted by the distinct inoculum conditions; all mesocosms converged with similar microbial community compositions by the end of the experiment. Instead, we infer that the 15 total populations were retained due to host selection, as they showed high metabolic potential for algal-bacterial interactions such as recycling alga-produced carbon and nitrogen and production of vitamins and secondary metabolites associated with algal growth and senescence, including B vitamins, tropodithietic acid, and roseobacticides. IMPORTANCE Bacteria proliferate in nutrient-rich aquatic environments, including engineered algal biofuel systems, where they remineralize photosynthates, exchange secondary metabolites with algae, and can influence system output of biomass or oil. Despite this, knowledge on the microbial ecology of algal cultivation systems is lacking, and the subject is worthy of investigation. Here, we used metagenomics to characterize the metabolic capacities of the predominant bacteria associated with the biofuel-relevant microalga Microchloropsis salina and to predict testable metabolic interactions between algae and manipulated communities of bacteria. We identified a previously undescribed and uncultivated organism that dominated the community. Collectively, the microbial community may interact with the alga in cultivation via exchange of secondary metabolites which could affect algal success, which we demonstrate as a possible outcome from controlled experiments with metabolically analogous isolates. These findings address the scalability of lab-based algal-bacterial interactions through to cultivation systems and more broadly provide a framework for empirical testing of genome-based metabolic predictions.

established within an algal culture system. We investigated whether distinct inoculations of bacteria may exhibit priority effects that divergently affect community development over time and whether these system dynamics may affect resulting algal success. Based on the molecular characterizations of the cultivation system, we predicted metabolic interactions between the alga and the bacteria, which we subsequently tested under controlled laboratory cocultivation with metabolically analogous strains. This study not only characterizes the composition, metabolic capacity, and ecological roles of bacteria that are recruited within a biotechnologically relevant algal cultivation system but also more broadly assesses the scalability of laboratory-defined algal-bacterial interactions in complex systems and is relevant to future studies seeking to engineer bacterial communities that enhance algal biofuel productivity.

RESULTS
Fifteen numerically dominant bacterial taxa associated with outdoor M. salina cultures. Microchloropsis salina algal mesocosms were initially established with one of three inoculum treatment conditions: (i) coculture with Phaeobacter gallaeciensis BS107, (ii) coculture with Alteromonas sp. strain I10, or (iii) axenic culture (Fig. 1A). M. salina biomass did not vary by bacterial inoculation treatment except for the final time points, when both mesocosms originally inoculated with the axenic culture had lower biomass than the other four mesocosms (Fig. 1C). While we were unable to attribute the differences in biomass to the presence of a specific bacterial taxon, we subsequently leveraged clues from bacterial community structure and metabolic capacity to infer system dynamics.
Retention of primary inoculated strains varied temporally and by treatment. Since four of the six mesocosms were initially inoculated with alga-associated bacterial strains of the genera Alteromonas and Phaeobacter, we first searched the metagenomes for these organisms. We recovered one nearly complete MAG for each of these two genera. The MAG corresponding to Phaeobacter gallaeciensis BS107 (named MSM14), with closest neighbor GTDB lineage as Phaeobacter gallaeciensis, exhibited 89.71% average pairwise nucleotide identity (ANI) compared to the P. gallaeciensis BS107 reference genome. The MSM14 genome was 78.07% complete, with 1.91% contamination and 59.1% GC content. These data suggest that the Phaeobacter strain initially added to the mesocosms had disappeared within the first 17 days and was replaced with a related yet distinct Phaeobacter sp. (MSM14) that was recruited to all 6 mesocosms and persisted at low to moderate abundances (Table S7A). On the other hand, the MAG corresponding to Alteromonas sp. I10 (MSM1), with closest neighbor GTDB lineage as Alteromonas macleodii, exhibited 99.99% ANI compared to the Alteromonas sp. I10 reference genome. The MSM1 genome was 100% complete, with 0.36% contamination and 44.6% GC content (Table S2). The presence of Alteromonas macleodii (MSM1) in its respective treatment mesocosms, and its dispersal to neighboring mesocosms, varied over time ( Fig. 2; Table S7A and Fig. S1). This shows that the initially added Alteromonas remained present throughout the two mesocosms where it was inoculated and had dispersed to the other four mesocosms as well.
Bacterial recruitment was taxon specific, and prevalence varied temporally. In addition to the two bacterial strains initially inoculated within respective treatments (one remaining throughout and the other showing lack of establishment/persistence), 14 other bacterial taxa were recruited to the 6 mesocosms at high enough abundances to recover medium-to high-quality MAGs (Tables S2 and S7A). Between 7 and 13 taxa cooccurred in any given sample. In order to estimate the contribution of each taxon within the collective community structure, we quantified the abundance of each bacterial MAG with a fold coverage metric (Table S7B). Fold coverage was also transformed relative to all MAGs within each sample, referred to hereafter as "relative coverage", to appropriately compare MAG predominance across treatments and over time (Table S7C). One MAG (unclassified Rhodobacteraceae MSM9) showed remarkably high relative coverage across all samples, as high as ;77% in one sample (;223-fold coverage), and on average representing ;50% of the MAG community. An additional 5 MAGs had intermediate abundance, with their genomes covering at least 1.6-fold in any given sample ( Fig. 2G; Fig. S1). The remaining 9 bacterial MAGs, including the inoculated Alteromonas sp., were less predominant systemwide and encompassed ,1% of the bacterial community ( Fig. 2H; Tables S7B and S7C).
We aimed to test whether MAG community taxonomic composition varied by inoculation treatment and temporally throughout the algal cultivation cycles and used relative coverage for this comparison. While mesocosms were initially inoculated with specific bacterial strain starting conditions (Phaeobacter, Alteromonas, axenic control), all taxa were shared among the treatment mesocosms by the midpoint of the experiment, as the overall MAG composition did not significantly vary with inoculation treatment (Bray-Curtis dissimilarity, permutational multivariate analysis of variance [PERMANOVA]: r 2 = 0.083, P = 0.082). The relative coverage (or abundance) of bacteria did however vary over time  . It subsequently increased in abundance through day 36 in the mesocosms that produced higher biomass (Alteromonas inoculated, 33.47% 6 0.79% relative coverage; Phaeobacter inoculated, 22.63% 6 6.65% relative coverage) but remained at lower abundances in the reduced-biomass mesocosms (control, 9.66% 6 0.60% relative coverage) ( Fig. 1C; Fig. S1). This was the only taxon whose abundance appeared to correlate with algal mesocosm biomass. MAG population-level variation via SNP detection. The high coverage of the unclassified Rhodobacteraceae (MSM9) MAG from all 30 metagenomes allowed us to examine genome heterogeneity across samples, both over time and across inoculation treatment. Since lower sequencing coverage may lead to false-positive detection of single nucleotide polymorphisms (SNPs), we did not perform this analysis on the remaining 14 MAGs. MSM9 represented a very high proportion of the total bacterial community and thus had a high fold coverage overall. As such, our analysis could determine if there were multiple recruitment events of this strain into the mesocosms, with either different strains being recruited to different mesocosms or new strains displacing old strains over time. Our analysis detected just 9 potential SNPs out of the entire 3.76-Mb genome (,0.001% variation). Further, these 6 potential SNPs were distributed seemingly randomly across mesocosms, meaning that there was no discernible SNP pattern by treatment or over time. Due to the lack of sampling earlier in the experiment, it is unknown what the exact order of arrival was of this particular bacterial strain.
M. salina-associated heterotrophic aquatic bacteria capable of aerobic C and N remineralization. Dense algal cultures and natural blooms are overly saturated with oxygen, at least during the day, due to photosynthesis, and algal cells are known to release dissolved organic carbon (DOC), so we expected the M. salina-associated bacteria to exhibit metabolic genes for survival under these conditions. We first examined the general biogeochemical capacity based on marker genes for carbon (C) and nitrogen (N) cycling (43). Overall, the M. salina bacterial community harbored pathways for common aerobic carbon cycling processes, such as aerobic C respiration, and processes less expected to be dominant in algal blooms, such as CO oxidation and aerobic C fixation. Examination of these patterns at the taxon level showed that the different bacterial taxa varied in their metabolic capacity for C utilization and remineralization. For example, while most bacteria (n = 13 MAGs) contained all required marker genes for aerobic respiration, Pseudophaeobacter sp. (MSM11) and Methylophaga sp. (MSM3) lacked 50% and 100% of the required genes for this process, respectively. A smaller subset of bacteria (n = 6 MAGs) were capable of CO oxidation, including P. gallaeciensis, Pseudophaeobacter sp. (MSM11), Roseibaca sp., and three unclassified Rhodobacteraceae members (MSM9, MSM12, MSM13). The capability for aerobic C fixation was restricted to just one bacterium, Methylophaga sp., which appears to be an obligate single-carbon utilizer, able to fix C with the phosphoribulokinase (PRK) and ribulose-bisphosphate carboxylase (rbcS) genes. One bacterium, Phaeobacter sp. (MSM14) may be capable of metabolizing methane via methanotrophy (Fig. 3). All bacterial taxa lacked the full pathways for anaerobic C cycling, such as fermentation, anaerobic C fixation, and methanogenesis (Table S3). However, four of the bacteria were predicted to be capable of photoheterotrophy, including Erythrobacter sp. (MSM4), Roseibaca sp. (MSM6), and two unclassified Rhodobacteraceae members (MSM12, MSM13), which contained 76.92%, 94.23%, 84.62%, and 92.31% of the 52 PGC (photosynthetic gene cluster) gene orthologs, respectively. Bacteria which were not predicted to be photoheterotrophs contained between 7.69% and 69.23% of the PGC genes (Table S4).
With respect to the N cycle, 3 of the 7 primary reactions were represented by M. salinaassociated bacteria: mineralization, N assimilation, and denitrification. Mineralization of N was the most common pathway among bacteria, with the exception of Marinobacter sp.  (Fig. 3B). None of the bacterial taxa were predicted to be capable of ammonification, annamox, N fixation, or nitrification (Table S3).
Bacteria demonstrate metabolic capabilities for influencing algal-bacterial interaction outcomes. In addition to describing the general biogeochemical capabilities of M. salina-inoculated and recruited bacterial taxa, we assessed the bacterial capacity for more specialized secondary metabolisms regarded for algal-bacterial interactions, including production of beneficial B vitamins, the antimicrobial compound tropodithietic acid (TDA), and the algicidal compound roseobacticide. Twelve of the 15 bacterial taxa were capable of production of at least one class of B vitamins, with 7 capable of B1 (thiamine) synthesis, 9 capable of B7 (biotin) synthesis, and 6 capable of B12 (cobalamin) synthesis. Of these 12 B vitamin-synthesizing bacteria, just two were predicted to synthesize all three types: Pseudophaeobacter sp. (MSM11) and Methylophaga sp. (MSM3). Neither Algoriphagus sp. (MSM8) nor two of the unclassified Rhodobacteraceae members (MSM15 and the abundant MSM9) were predicted to biosynthesize any of the B vitamins (Table S5). We also examined the biphasic symbiosis mediated by bacterial metabolism of tropones-algal growth-promoting TDA (tropodithietic acid) and algicidal roseobacticides. This complex process is driven by bacterial metabolism of the alga-derived metabolite p-coumaric acid (pCA), which may occur either through tolerance/transport (44) or degradation (45). Bacteria produce other required precursors/intermediates via phenylacetic acid (PAA) catabolism (46) and sulfur metabolism, and products from those aforementioned steps are collectively incorporated into a tropone backbone (47) for both TDA or roseobacticide secondary metabolites. Ten M. salina-associated genomes (n = 6 Alphaproteobacteria, n = 3 Gammaproteobacteria, n = 1 Bacteroidia) contained the genes for processing pCA, either through reactions for partial degradation to 4-hydroxybenzoate (MSM5, MSM15), full degradation through to the tricarboxylic acid (TCA) cycle (MSM9, MSM11, MSM12, MSM13, MSM14), or tolerance to pCA via phospholipid transport (MSM1, MSM3), with Marinobacter sp. (MSM2) capable of both degradation and tolerance. Five of these taxa plus two others (MSM6, MSM7) contained the full pathway for PAA catabolism. All MAGs were capable of metabolizing sulfur intermediates via the cysI or patB genes. Only one taxon, Phaeobacter gallaeciensis (MSM14), was predicted to synthesize either TDA or roseobacticide compounds via TDA, as this MAG contained the suite of genes for TDA resistance (tdaR1 to tdaR5) (48) and TDA backbone biosynthesis (tdaA to tdaE) ( Fig. 4; Table S6).
Based on the genomic characterizations of the M. salina mesocosm bacterial community, we formed predictions of algal-bacterial interactions mediated by the metabolic activity of the rest of the bacterial community and tested this through a subsequent experiment. Specifically, we hypothesized that pCA-degrading bacteria might provide a mechanism to prevent algicidal activity by roseobacticide producers. We tested this under a more controlled and simplified community than the mesocosm setup, i.e., the alga M. salina was cocultured with the predicted TDA/roseobacticide producer P. gallaeciensis, with or without a bacterial consortium of 8 bacterial strains capable of pCA degradation. While pCA is produced by microalgae, here it was added exogenously to ensure availability at standardized concentrations. Given that pCA is required for bacterial production of roseobacticides, we predicted that supplementing exogenous pCA would stimulate roseobacticide production by P. gallaeciensis, resulting in M. salina growth inhibition. On the other hand, we anticipated that the addition of pCA-degrading bacteria would reduce the availability of this precursor in the medium, slowing or preventing P. gallaeciensis roseobacticide production and instead potentially switching bacterial metabolism to TDA synthesis, and M. salina growth would not be inhibited under these conditions. The experimental results were consistent with this hypothesis. Overall, we observed that P. gallaeciensis reduced M. salina biomass in comparison to axenic algal cultures, with an additional significant reduction in biomass when exogenous pCA was available (Fig. 4A). The pCAdegrading community incapable of roseobacticide synthesis also inhibited algal growth (Fig. 4B), through an as-yet-unknown mechanism. However, the inhibitory effect was not observed when both P. gallaeciensis and the pCA-degrading bacterial community were combined, with algal biomass equal to or slightly enhanced over that of axenic cultures ( Fig. 4C;  Fig. S2).

DISCUSSION
Deterministic community assembly. The addition of mutualistic bacteria in engineered algal cultivation systems has been proposed as a potential strategy for increasing biomass output (for examples, see references 22 and 24). It has been shown previously that the taxa applied in this experiment, Phaeobacter gallaeciensis (49,50) and Alteromonas sp. (51,52), can provide a growth-enhancing effect to M. salina when grown in pure coculture under certain conditions. Therefore, we hypothesized that the growth-enhancing effect of both bacterial strains would be conserved in this scaled-up cultivation system, with inoculated treatments exhibiting a higher algal biomass than that of the initially axenic control mesocosms. As the growth-enhancing effects observed in the laboratory occurred on very short time scales (days), we expected to observe that effect in this system by the end of round 1 of algal cultivation, yet we detected no discernible difference in algal growth or biomass accumulation under the different inoculation treatments up through round 4 (Fig. 1C). Rather, a difference in algal biomass across treatments began to appear only at the very final time points, at the end of round 5 (;38 days postinoculation), when the biomass of M. salina in the control mesocosms was ;18% and ;14% lower than that in the Phaeobacter-and Alteromonas-inoculated treatments, respectively. While this result differed from our expectation, it might suggest that primary inoculation with a mutualist may lead to increased resilience of algae over longer-term time scales and is worth future investigation. These findings may be of relevance for applied biofuel systems, as resilience may prevent sudden crashes in algal cultivation systems that decimate biofuel output (26).
Mutualistic bacteria have also been shown to indirectly increase algal output by altering the composition of the surrounding bacterial community, for example, by excluding competitive or antagonistic bacteria (14). Here, by establishing discrete starting coculture conditions where we initially inoculated M. salina with one of two bacterial strains and then subsequently facilitated recruitment of dispersed bacteria within the open, outdoor cultivation system (Fig. 1A), our study aimed to investigate whether these bacteria may additionally affect M. salina biomass via microbe-microbe interactions. We hypothesized that we would observe altered community assembly trajectory patterns both between our bacterially inoculated treatments and in comparison to the axenic control, potentially due to bacteria exerting priority effects and altering the subsequent establishment and assembly trajectory of the bacterial community (31,53), resulting in divergent community composition patterns across treatments. However, our results do not fully support this hypothesis. As a whole, system-wide community composition did not vary with initial coculture inoculation treatment, as all predominant bacterial taxa were shared across mesocosms by day 17 of the experiment (cultivation round 3) (Fig. 2). There were slight differences in the estimated abundances (via relative coverage) of bacterial community members over time; for example, Erythrobacter sp. (MSM4) showed a lower abundance in the two initially axenic mesocosms than in the four started with the bacterial isolate strains at day 36 (see Fig. S1 in the supplemental material). This decline in Erythrobacter sp. was followed by a subsequent decrease in algal productivity 3 days later (Fig. 1C); however, this result is correlative, and additional variables not measured here may have also contributed to decreased algal productivity. Nevertheless, the correlation data suggest that Erythrobacter sp. could be a potential algal mutualist and justifies subsequent investigations of this understudied genus.
We observed cross-mesocosm dispersal of inoculated strains in addition to environmentally sourced bacteria, specifically the Alteromonas sp. which we detected in both control mesocosms midway through the experiment (day 17) (Table S7). In contrast, the Phaeobacter sp. isolate was undetectable at all time points sampled, suggesting that this bacterium failed to establish and persist. It is unclear to what extent environmental conditions and/or competitive biotic interactions influenced the establishment, abundance, and persistence of inoculated isolates and additionally how these patterns may have influenced downstream community assembly dynamics and algal-bacterial interactions. Thus, we advocate for subsequent quantitative experiments that investigate the interplay between abiotic and biotic interactions on algal-bacterial mutualisms. Interestingly, while the Phaeobacter isolate failed to persist, a related yet distinct Phaeobacter sp. (MSM14) was recruited into the mesocosms. Comparative genomics show that these are divergent strains or species, as evidenced by low nucleotide identity (89%) and distinct functional capacities. Features present in the MAG and absent from the isolate genome included protein families for siderophore and iron transport, methane monooxygenase, urea transporters and urease proteins, and vitamin B12 transporters (Table S3). The presence of methane monooxygenase genes might suggest that this organism is capable of methane metabolism and was potentially selected for in this highly reactive environment where M. salina may release intermediates for this process (54). In combination with complete pathways for CO oxidation and respiration (Fig. 3A), this provides compelling evidence that the M. salina-associated bacteria can utilize diverse single-carbon sources, and their metabolic flexibility may contribute to their enhanced survival here. The other identified functions are suggestive of a mutualistic lifestyle where bacteria potentially influence the mobilization of iron and provide vitamins that promote algal growth. Additional phenotypic and metabolic quantification are needed to test these hypotheses in the future, for example, measuring substrate incorporation and transfer from bacteria to algal host via nanoscale stable isotope tracing (NanoSIMS) (3,23).
Persistent bacteria capable of CO oxidation and denitrification. The low-diversity community associated with M. salina was numerically dominated by one taxon-a previously uncultured and undescribed member of the Rhodobacteraceae (MSM9; Fig. 2)-whose abundance equaled that of the other 14 MAGs combined (Table S7C). We postulate that this curated MAG is representative of one single strain, as there was no evidence of genomic variation across mesocosms or over time (9 SNPs identified out of the entire 3.76-Mb genome, or ,0.001% variation). Thus, after this bacterium was recruited within at least one mesocosm, it quickly became established across the system and was retained over time, without being replaced by other closely related strains. We infer that the persistence of the Rhodobacteraceae species may likely be due in part to its metabolic flexibility for energy production. First, this organism, like many other members of the Rhodobacteraceae (55-57) including the other 4 MAGs identified within this family, is a CO oxidizer. This metabolism is hypothesized to be a physiological adaptation to enhance survival under a wide range of DOC availability (58). Second, Rhodobacteraceae MSM9 also has the genomic capability for denitrification, again an energy-generating process, but in this case one that could contribute to system-wide nitrogen loss to the atmosphere and potentially reduce the availability of nitrate for algae under more N-limited conditions. The fact that the most numerically dominant organism in the outdoor cultures maintains the genetic capability for two major energy-generating metabolisms suggests that the alga-associated bacteria may have been carbon limited.  (49). Four additional taxa (MSM10, MSM12, MSM13, MSM15) could not be assigned to a genus-level annotation despite the nearly complete assemblies, suggesting that neither they nor their close relatives have been cultivated previously. They share metabolic features with previously described alga-associated bacteria, which may have contributed to their recruitment here. Similar to the predominant Rhodobacteraceae, they display diverse carbon biogeochemical capacity, including CO oxidization potential (MSM12, MSM13) ( Fig. 3; Table S3A). Two appear to be photoheterotrophs (MSM12, MSM13), which is a function shared with two other taxa (MSM4, MSM6) via conservation of the PGC gene cluster for anoxygenic photosynthesis (Table S4). The taxa which we identified as potential photoheterotrophs are all classified within the Proteobacteria phylum, with three classified as Rhodobacteraceae, an observation that is in agreement with previous studies on the evolution of photoheterotrophy in bacteria (62,63). One of these taxa (MSM12) appears to be a denitrifying bacterium and, along with another Rhodobacteraceae species (MSM9) and Marinobacter sp. (MSM2), likely contributed to the N cycle via conversion of inorganic nitrate to atmospheric nitrogen gas, potentially promoting nitrogen loss in the system (64). We postulate that three of these four previously undescribed bacterial taxa also likely interact with M. salina via production of essential B vitamins, as the metabolic capacity for biosynthesis of one or more classes of B vitamins was identified, specifically, B1 (MSM10), B7 (MSM10), and B12 (MSM12, MSM13) among other taxa (all MAGs excluding MSM8, MSM9, MSM15). Algae require B vitamins and cannot biosynthesize these compounds; hence, they are dependent on an exogenous source. While the cultivation medium used in this system was replete with vitamins to encourage algal growth, these taxa may be able to maintain algal growth if vitamins become limited over time.
Community cross-feeding potential for TDA/roseobacticide synthesis. The genomic evidence from the outdoor cultures and the simple laboratory experiment together suggest that microbe-microbe interactions may influence bacterial metabolite synthesis, particularly the TDA and roseobacticide system that influences algal growth, but likely others that have not been characterized previously. Notably, the data suggest that multiple, divergent taxa could contribute distinct metabolites shown to be precursors for P. gallaeciensis production of algal-stimulating (TDA) or algal-senescing (roseobacticide) compounds (12,44,49,65). We found it intriguing that different organisms were capable of discrete steps in this metabolic process, as toluene tolerance, PAA catabolism, and TDA production capabilities collectively have only primarily been described for two taxa (49,65). Thus, we expected the collective process to either be conserved or lacking completely in a given genome. TDA and roseobacticide production are energetically expensive processes, and we wonder whether bacterial partitioning of the steps involved may improve community fitness via cross-feeding, a phenomenon that recent studies have shown to be prevalent in microbial community dynamics (66)(67)(68). We have proposed a communitybased cooperative conceptual model (Fig. 4) that is based on our genomic observations and supported by our targeted experimental results. With increasing pCA availability, the predicted roseobacticide-producing isolate P. gallaeciensis BS107 reduced the growth of M. salina compared to axenic algae free of bacteria. However, this growth-inhibiting effect can be alleviated if pCA-degrading bacteria are added to the system. Our experiments and conceptual model suggest a cross-feeding hypothesis that pCA degraders may prevent algicidal roseobacticide production and subsequent algal growth inhibition, potentially through competition for pCA. If the roseobacticide precursor compound pCA were to become scarce due to bacterial degradation, this would instead presumably result in production by P. gallaeciensis of algal growth-promoting TDA compounds, since the metabolic pathway is shared but the product varies depending on availability of precursors (49). Although the pCA degrader synthetic community was predicted to be incapable of roseobacticide production from pCA, it is interesting that this consortium inhibited algal growth only when pCA was added (Fig. S2); however, the mechanism underlying the observed phenotype was not discernible based on genomic features. Since the genomebased metabolic predictions and growth assays were used as the basis for our cooperative model, we acknowledge that we lack corroborating characterizations of bacterial activity and quantification of chemical metabolites, which precludes us from assessing whether metabolites were actively being exchanged and influenced algal biomass as a result. Future investigations should combine a multi-omics (metagenomics, transcriptomics, metabolomics) approach to pinpoint the mechanisms of this proposed cooperative symbiosis.
Through the results presented here, we have begun to describe the metabolic capacity of alga-associated bacteria in high-biomass cultivation systems and to infer the community dynamics both among bacteria and with an algal partner. Given that heterotrophic bacteria tightly associate with bioenergy-relevant organisms, we expect this study to provide a framework for subsequent research to consider how ecological interactions mediate bioenergy output. These data suggest that for complex systems, ecological outcomes are difficult to predict using a reductionist approach, for example, by inoculating bioenergy crops with a single mutualist to increase output. Rather, through the experiment described here, we show that competitive interactions and metabolic complementarity must be considered, as a collective community may not behave as the sum of its parts.

MATERIALS AND METHODS
Algal and bacterial cultivation. Prior to establishment of experimental mesocosms, microalgal and bacterial strains were cultivated in the laboratory as follows. A strain of Microchloropsis salina (CCMP 1776) was obtained from the Provasoli-Guillard National Center for Marine Algae and Microbiota (CCMP-NCMA) at Bigelow Laboratory for Ocean Sciences (West Boothbay Harbor, MA) and was maintained in the laboratory as an axenic culture in enriched seawater medium (ESAW; pH 8.2) (69) at 21°C under constant light (100 mmol photons m 22 s 21 ) and aeration. Two strains of alga-associated aquatic bacteria, Alteromonas sp. I10 and Phaeobacter gallaeciensis BS107, were acquired from the laboratory's collection of bacterial isolate stocks. Isolates were cultivated from glycerol stocks by inoculating a scraping of the stock into 5 mL sterile marine broth liquid medium and shaken overnight at 30°C. Algal-bacterial cocultures were initiated by inoculating a 1-L axenic M. salina culture with 10 mL of the overnight bacterial cultures of either P. gallaeciensis BS107 or Alteromonas sp. I10, and these constituted the coculture inoculations for the subsequent experimental mesocosms.
Mesocosm experimental design. Experimental mesocosms were established 1 August 2013, cultivated in a multiround harvest, and sampled at intervals until conclusion on 10 September 2013 as follows. Mesocosms were inoculated with either (i) a coculture of M. salina and the bacterium Phaeobacter gallaeciensis BS107, (ii) a coculture of M. salina and the bacterium Alteromonas sp. I10, or (iii) an axenic culture of M. salina without bacteria (Fig. 1). Experimental treatments were set up in duplicate, for a total of 6 mesocosms. Mesocosms containing 16 L of ESAW medium with full F/2 nutrients (70) (see Text S1 in the supplemental material) were inoculated with M. salina (either in coculture with bacteria or axenic) at an initial cell density of ;200,000 algal cells mL 21 and were maintained in continuous aeration via ambient air in a semicontinuous mode. Mesocosms were maintained outdoors, uncovered, and exposed to elements, and no additional temperature or lighting controls were necessary.
After initial establishment and inoculation, experimental algal-bacterial mesocosms were grown for a total duration of 40 days, with a harvesting of cells and replenishment of the culture with fresh sterile medium approximately every 8 to 9 days, when cultures presumably reached stationary phase, in order to prevent a population bottleneck. A total of 5 harvest events, in which 90% of culture was removed and 10% was retained, occurred throughout the duration of the experiment, including the final harvest at the conclusion (Fig. 1B). During each harvest cycle for all 6 mesocosms, the biomass of M. salina was measured using relative fluorescence units (RFU) of chlorophyll a and samples were collected for downstream metagenomic sequencing, modified from a previous report (33) as follows. In short, 500 mL from each mesocosm was vacuum filtered onto a 0.2-mm Supor polyethersulfone membrane filter (Pall, Port Washington, NY), with collection of algae along with both free-living and attached bacteria. Membranes were then vortexed in 10 mL artificial seawater, washing cells off the membrane and with the intent to dislodge bacteria attached to algal cells, which was repeated twice for each membrane. The resulting 20 mL of suspension was filtered through a 1-mm membrane filter to remove most M. salina cells while retaining most bacterial cells, and the filtrate was transferred to sterile 50-mL conical tubes, pelleted by centrifugation (10,000 Â g for 5 min), and stored at 280°C prior to DNA extraction.
DNA extraction, library preparation, and metagenomic sequencing. Total nucleic acids were extracted and purified with a metagenomic DNA isolation kit for water (Epicentre, WI) in accordance with the manufacturer's protocol. Metagenomic libraries were constructed using the Nextera XT DNA library prep kit (Illumina, San Diego, CA) with 1 ng starting input double-stranded DNA (dsDNA). Resulting libraries with an average 250-bp insert size were paired-end sequenced (2 Â 150) on a HiSeq 2000 (Illumina, San Diego, CA) at Sandia National Laboratories (Livermore, CA).
MAGs. Using bioinformatics, metagenomes were demultiplexed, filtered, and trimmed prior to assembly and annotation. Resulting high-quality sequences were collectively cross-assembled and binned into metagenome-assembled genomes (MAGs), refined, and assessed for completeness (Text S1). Bacterial MAG taxonomic annotation was assigned via phylogeny inference using GTDB-Tk v1.5.0 and the release 202 database (71,72). Gene calls for MAGs and initial functional annotation were done with PATRIC (73). The metabolic capacity of bacterial MAGs was annotated via HMM or blastp searches against marker gene databases curated with genes and pathways pulled from KEGG (74) for C and N biogeochemical cycles (43), photoheterotrophy (62,63), B vitamin production (9,75), and TDA/roseobacticide production (12,44,65,76). As the MAGs were not fully complete genomes, we considered a metabolic pathway to be functionally complete if $75% of the genes were detected. MAG depth and breadth of read coverage were calculated by mapping individual sample fastq files with bbmap v35.85 and perfectmode=t. MAG presence/absence was determined via breadth calculations (Table S7A). The estimated abundance of each MAG was calculated using the mean of the median fold coverage (Table S7B) and was normalized by transformation to a relative proportion measurement within sample to appropriately compare MAG abundances across samples (Table S7C). The relative proportion of the mean of the median fold coverage is referred to as the "relative coverage" throughout the article for clarity. SNP analysis was done individually for each MAG by using each of the 30 quality-processed sequencing read sets using SNIPPY v4.6 (77). Average nucleotide identity between reference isolate genomes for Phaeobacter gallaeciensis BS107 and Alteromonas sp. I10 and MAGs with similar taxonomy assignments as determined by GTDB were calculated with FastANI (78). All statistical analyses were performed in R (79) and RStudio (80) using the packages "vegan" (81) and "phyloseq" (82).
Experiments testing genome-based predictions of algal-bacterial interactions. Genome-based metabolic predictions underlying algal-bacterial interactions were tested under controlled laboratory conditions. In short, we tested the effect of the algicidal roseobacticide-producing bacterium P. gallaeciensis on the growth of the alga M. salina under the addition of the exogenous roseobacticide precursor p-coumaric acid (pCA) (49). To test the potential for bacterial interference of roseobacticide production via pCA degradation, a portion of the cocultures were inoculated with a laboratory-synthesized mixed community of pCA-degrading bacteria. Roseobacticide-producing P. gallaeciensis BS107 was available for experiments from the original mesocosm experiments; however, the mesocosm-associated bacterial genomes characterized as pCA degraders were uncultivated MAGs. As a proxy, 8 related bacterial strains previously isolated from algae (59) containing the identical and complete suite of genes for pCA degradation (44) and lacking genomic capability for roseobacticide production were used in the experiments (Text S1). Axenic M. salina cultures were either (i) left axenic, (ii) inoculated with P. gallaeciensis, (iii) inoculated with the pCA-degrading community, or (iv) inoculated with both P. gallaeciensis and the pCA degraders. Cultures were grown in F/2 medium either without or with pCA (1.7 mM) (Sigma-Aldrich, St. Louis, MO) and maintained at 22°C on a 14-h-light:10-h-dark cycle. Biomass was measured daily using chlorophyll a relative fluorescence units on a fluorometer (Turner Designs, San Jose, CA).

SUPPLEMENTAL MATERIAL
Supplemental material is available online only. TEXT S1, DOCX file, 0.02 MB. FIG S1, TIF   The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.