Dissecting host-microbiome interaction effects on phytoplankton microbiome composition and diversity.

Phytoplankton and their associated microbiomes of heterotrophic bacteria are foundational to primary production, energy transfer, and biogeochemical cycling in aquatic systems. While it is known that these microbiomes are shaped by host-released dissolved organic matter (DOM), the extent to which dynamic phytoplankton-bacteria interactions shape bacterial community assembly remains to be examined. Here, we investigated the effects of two mechanisms in host-microbiome interactions on phytoplankton bacterial microbiome formation: (i) innate host selection and (ii) host-microbiome feedback. For the former, phytoplankton-produced DOM composition is based solely on the host's properties (species or physiological state); for the latter, the presence of the microbiome modifies host DOM production. The microbiome of Chlorella sorokiniana was extracted and exposed to six ratios of the two effects. We found that microbiome composition changed along with the six host-microbiome feedback versus innate host selection ratios, with the highest compositional distance between communities under the strongest and the weakest ratio of the two effects. This indicates that each mechanism selects for different bacterial species. In addition, our findings showed that when both selective forces were applied, it led to a higher community richness, while host-microbiome feedback alone reduces community evenness due to its strong species-specific selection.

Comparative genomic analysis of Microcystis strain diversity using conserved marker genes.

Microcystis-dominated cyanobacterial harmful algal blooms (cyanoHABs) have a global impact on freshwater environments, affecting both wildlife and human health. Microcystis diversity and function in field samples and laboratory cultures can be determined by sequencing whole genomes of cultured isolates or natural populations, but these methods remain computationally and financially expensive. Amplicon sequencing of marker genes is a lower cost and higher throughput alternative to characterize strain composition and diversity in mixed samples. However, the selection of appropriate marker gene region(s) and primers requires prior understanding of the relationship between single gene genotype, whole genome content, and phenotype. To identify phylogenetic markers of Microcystis strain diversity, we compared phylogenetic trees built from each of 2,351 individual core genes to an established phylogeny and assessed the ability of these core genes to predict whole genome content and bioactive compound genotypes. We identified single-copy core genes better able to resolve Microcystis phylogenies than previously identified marker genes. We developed primers suitable for current Illumina-based amplicon sequencing with near-complete coverage of available Microcystis genomes and demonstrate that they outperform existing options for assessing Microcystis strain composition. Results showed that genetic markers can be used to infer Microcystis gene content and phenotypes such as potential production of bioactive compounds , although marker performance varies by bioactive compound gene and sequence similarity. Finally, we demonstrate that these markers can be used to characterize the Microcystis strain composition of laboratory or field samples like those collected for surveillance and modeling of Microcystis-dominated cyanobacterial harmful algal blooms.

The Western Lake Erie culture collection: A promising resource for evaluating the physiological and genetic diversity of Microcystis and its associated microbiome.

Cyanobacterial harmful algal blooms (cyanoHABs) dominated by Microcystis spp. have significant public health and economic implications in freshwater bodies around the world. These blooms are capable of producing a variety of cyanotoxins, including microcystins, that affect fishing and tourism industries, human and environmental health, and access to drinking water. In this study, we isolated and sequenced the genomes of 21 primarily unialgal Microcystis cultures collected from western Lake Erie between 2017 and 2019. While some cultures isolated in different years have a high degree of genetic similarity (genomic Average Nucleotide Identity >99%), genomic data show that these cultures also represent much of the breadth of known Microcystis diversity in natural populations. Only five isolates contained all the genes required for microcystin biosynthesis while two isolates contained a previously described partial mcy operon. Microcystin production within cultures was also assessed using Enzyme-Linked Immunosorbent Assay (ELISA) and supported genomic results with high concentrations (up to 900 µg L⁻¹) in cultures with complete mcy operons and no or low toxin detected otherwise. These xenic cultures also contained a substantial diversity of bacteria associated with Microcystis, which has become increasingly recognized as an essential component of cyanoHAB community dynamics. These results highlight the genomic diversity among Microcystis strains and associated bacteria in Lake Erie, and their potential impacts on bloom development, toxin production, and toxin degradation. This culture collection significantly increases the availability of environmentally relevant Microcystis strains from temperate North America.

Variation in resource competition traits among Microcystis strains is affected by their microbiomes.

Freshwater harmful algal blooms are often dominated by Microcystis, a phylogenetically cohesive group of cyanobacteria marked by extensive genetic and physiological diversity. We have previously shown that this genetic diversity and the presence of a microbiome of heterotrophic bacteria influences competitive interactions with eukaryotic phytoplankton. In this study, we sought to explain these observations by characterizing Monod equation parameters for resource usage (maximum growth rate µ max, half-saturation value for growth K s, and quota) as a function of N and P levels for four strains (NIES-843, PCC 9701, PCC 7806 [WT], and PCC 7806 ΔmcyB) in presence and absence of a microbiome derived from Microcystis isolated from Lake Erie. Results indicated limited differences in maximum growth rates but more pronounced differences in half-saturation values among Microcystis strains. The largest impact of the microbiome was reducing the minimal nitrogen concentration sustaining growth and reducing half saturation values, with variable results depending on the Microcystis strain. Microcystis strains also differed from each other in their N and P quotas and the extent to which microbiome presence affected them. Our data highlight the importance of the microbiome in altering Microcystis-intrinsic traits, strain competitive hierarchies, and thus bloom dynamics. As quota, µ max, and K s are commonly used in models for harmful algal blooms, our data suggest that model improvement may be possible by incorporating genotype dependencies of resource-use parameters.

Impacts of an invasive filter-feeder on bacterial biodiversity are context dependent.

Bacteria represent most of the biodiversity and play key roles in virtually every ecosystem. In doing so, bacteria act as part of complex communities shaped by interactions across all domains of life. Here, we report on direct interactions between bacteria and dreissenid mussels, a group of invasive filter-feeders threatening global aquatic systems due to high filtration rates. Previous studies showed that dreissenids can impact bacterial community structure by changing trait distributions and abundances of specific taxa. However, studies on bacterial community effects were conducted using water from Lake Michigan (an oligotrophic lake) only, and it is unknown whether similar patterns are observed in systems with differing nutrient regimes. We conducted ten short-term dreissenid grazing experiments in 2019 using water from two eutrophic lake regions-the western basin of Lake Erie and Saginaw Bay in Lake Huron. Predation by dreissenids led to decline in overall bacterial abundance and diversity in both lakes. However, feeding on bacteria was not observed during every experiment. We also found that traits related to feeding resistance are less phylogenetically conserved than previously thought. Our results highlight the role of temporal, spatial, and genomic heterogeneity in bacterial response dynamics to a globally important invasive filter feeder.

Effects of Phycosphere Bacteria on Their Algal Host Are Host Species-Specific and Not Phylogenetically Conserved.

Phytoplankton is fundamental to life on Earth. Their productivity is influenced by the microbial communities residing in the phycosphere surrounding algal cells. Expanding our knowledge on how algal-bacterial interactions affect algal growth to more hosts and bacteria can help elucidate general principles of algal-host interactions. Here, we isolated 368 bacterial strains from phycosphere communities, right after phycosphere recruitment from pond water and after a month of lab cultivation and examined their impacts on growth of five green algal species. We isolated both abundant and rare phycosphere members, representing 18.4% of the source communities. Positive and neutral effects predominated over negative effects on host growth. The proportion of each effect type and whether the day of isolation mattered varied by host species. Bacteria affected algal carrying capacity more than growth rate, suggesting that nutrient remineralization and toxic byproduct metabolism may be a dominant mechanism. Across-host algal fitness assays indicated host-specific growth effects of our isolates. We observed no phylogenetic conservation of the effect on host growth among bacterial isolates. Even isolates with the same ASV had divergent effects on host growth. Our results emphasize highly specific host-bacterial interactions in the phycosphere and raise questions as to which mechanisms mediate these interactions.

The genetic and ecophysiological diversity of Microcystis.

Microcystis is a cyanobacterium that forms toxic blooms in freshwater ecosystems around the world. Biological variation among taxa within the genus is apparent through genetic and phenotypic differences between strains and via the spatial and temporal distribution of strains in the environment, and this fine-scale diversity exerts strong influence over bloom toxicity. Yet we do not know how varying traits of Microcystis strains govern their environmental distribution, the tradeoffs and links between these traits, or how they are encoded at the genomic level. Here we synthesize current knowledge on the importance of diversity within Microcystis and on the genes and traits that likely underpin ecological differentiation of taxa. We briefly review spatial and environmental patterns of Microcystis diversity in the field and genetic evidence for cohesive groups within Microcystis. We then compile data on strain-level diversity regarding growth responses to environmental conditions and explore evidence for variation of community interactions across Microcystis strains. Potential links and tradeoffs between traits are identified and discussed. The resulting picture, while incomplete, highlights key knowledge gaps that need to be filled to enable new models for predicting strain-level dynamics, which influence the development, toxicity and cosmopolitan nature of Microcystis blooms.

Host specificity of microbiome assembly and its fitness effects in phytoplankton.

Insights into symbiosis between eukaryotic hosts and their microbiomes have shifted paradigms on what determines host fitness, ecology, and behavior. Questions remain regarding the roles of host versus environment in shaping microbiomes, and how microbiome composition affects host fitness. Using a model system in ecology, phytoplankton, we tested whether microbiomes are host-specific, confer fitness benefits that are host-specific, and remain conserved in time in their composition and fitness effects. We used an experimental approach in which hosts were cleaned of bacteria and then exposed to bacterial communities from natural environments to permit recruitment of microbiomes. We found that phytoplankton microbiomes consisted of a subset of taxa recruited from these natural environments. Microbiome recruitment was host-specific, with host species explaining more variation in microbiome composition than environment. While microbiome composition shifted and then stabilized over time, host specificity remained for dozens of generations. Microbiomes increased host fitness, but these fitness effects were host-specific for only two of the five species. The shifts in microbiome composition over time amplified fitness benefits to the hosts. Overall, this work solidifies the importance of host factors in shaping microbiomes and elucidates the temporal dynamics of microbiome compositional and fitness effects.

Genotype and host microbiome alter competitive interactions between Microcystis aeruginosa and Chlorella sorokiniana.

Cyanobacterial harmful algal blooms (cyanoHABs) continue to increase in frequency and magnitude, threatening global freshwater ecosystems and services. In north-temperate lakes cyanobacteria appear in early summer, succeeding green algae as the dominant phytoplankton group, a pattern thought to be mediated by changes in temperature and bioavailable nutrients. To understand additional drivers of this successional pattern our study used reciprocal invasion experiments to examine the competitive interaction between Microcystis aeruginosa, a dominant contributor to cyanoHABs, and the green alga Chlorella sorokiniana. We considered two factors that may impact these interactions: (1) strain variation, with a specific emphasis on the presence or absence of the gene for the hepatotoxin microcystin, and (2) host-associated bacteria. We used toxic M. aeruginosa PCC 7806 (microcystin producing strain), a non-toxic mutant of PCC 7806, non-toxic M. aeruginosa PCC 9701 (non-microcystin producing strain), and C. sorokiniana. Each organism was available free of all bacteria (i.e., axenic) and with a re-introduced defined bacterial community to generate their xenic counterparts. Competitive interactions were assessed with reciprocal invasion experiments between paired xenic and paired axenic populations of C. sorokiniana and one of the two Microcystis strains, each assessed separately. Flow cytometry and random forest models were used to rapidly discriminate and quantify phytoplankton population densities with 99% accuracy. We found that M. aeruginosa PCC 7806, but not strain PCC 9701, could proliferate from low abundance in a steady-state population of C. sorokiniana. Further, the presence of bacteria allowed M. aeruginosa PCC 7806 to grow to a higher population density into an established C. sorokiniana population than when grown axenic. Conversely, when M. aeruginosa was dominant, C. sorokiniana was only able to proliferate from low density into the PCC 9701 strain, and only when axenic. The mutant of PCC 7806 lacking the ability to produce microcystin behaved similarly to the toxic wild-type, implying microcystin is not responsible for the difference in competitive abilities observed between the two wild-type strains. Quantification of microcystins (MCs) when PCC 7806 M. aeruginosa was introduced into the C. sorokiniana culture showed two-fold more MCs per cell when host-associated bacteria were absent compared to present in both species cultures. Our results show that the ability of M. aeruginosa to compete with C. sorokiniana is determined by genomic differences beyond genes involved in microcystin toxin generation and indicate an important role of host-associated bacteria in mediating phytoplankton interspecies interactions. These results expand our understanding of the key drivers of phytoplankton succession and the establishment and persistence of freshwater harmful cyanobacterial blooms.

To Dereplicate or Not To Dereplicate?

Metagenome-assembled genomes (MAGs) expand our understanding of microbial diversity, evolution, and ecology. Concerns have been raised on how sequencing, assembly, binning, and quality assessment tools may result in MAGs that do not reflect single populations in nature. Here, we reflect on another issue, i.e., how to handle highly similar MAGs assembled from independent data sets. Obtaining multiple genomic representatives for a species is highly valuable, as it allows for population genomic analyses; however, when retaining genomes of closely related populations, it complicates MAG quality assessment and abundance inferences. We show that (i) published data sets contain a large fraction of MAGs sharing >99% average nucleotide identity, (ii) different software packages and parameters used to resolve this redundancy remove very different numbers of MAGs, and (iii) the removal of closely related genomes leads to losses of population-specific auxiliary genes. Finally, we highlight some approaches that can infer strain-specific dynamics across a sample series without dereplication.

Temperature and Nutrient Levels Correspond with Lineage-Specific Microdiversification in the Ubiquitous and Abundant Freshwater Genus Limnohabitans.

Most freshwater bacterial communities are characterized by a few dominant taxa that are often ubiquitous across freshwater biomes worldwide. Our understanding of the genomic diversity within these taxonomic groups is limited to a subset of taxa. Here, we investigated the genomic diversity that enables Limnohabitans, a freshwater genus key in funneling carbon from primary producers to higher trophic levels, to achieve abundance and ubiquity. We reconstructed eight putative Limnohabitans metagenome-assembled genomes (MAGs) from stations located along broad environmental gradients existing in Lake Michigan, part of Earth's largest surface freshwater system. De novo strain inference analysis resolved a total of 23 strains from these MAGs, which strongly partitioned into two habitat-specific clusters with cooccurring strains from different lineages. The largest number of strains belonged to the abundant LimB lineage, for which robust in situ strain delineation had not previously been achieved. Our data show that temperature and nutrient levels may be important environmental parameters associated with microdiversification within the Limnohabitans genus. In addition, strains predominant in low- and high-phosphorus conditions had larger genomic divergence than strains abundant under different temperatures. Comparative genomics and gene expression analysis yielded evidence for the ability of LimB populations to exhibit cellular motility and chemotaxis, a phenotype not yet associated with available Limnohabitans isolates. Our findings broaden historical marker gene-based surveys of Limnohabitans microdiversification and provide in situ evidence of genome diversity and its functional implications across freshwater gradients.IMPORTANCELimnohabitans is an important bacterial taxonomic group for cycling carbon in freshwater ecosystems worldwide. Here, we examined the genomic diversity of different Limnohabitans lineages. We focused on the LimB lineage of this genus, which is globally distributed and often abundant, and its abundance has shown to be largely invariant to environmental change. Our data show that the LimB lineage is actually comprised of multiple cooccurring populations for which the composition and genomic characteristics are associated with variations in temperature and nutrient levels. The gene expression profiles of this lineage suggest the importance of chemotaxis and motility, traits that had not yet been associated with the Limnohabitans genus, in adapting to environmental conditions.

Microhabitats are associated with diversity-productivity relationships in freshwater bacterial communities.

Eukaryotic communities commonly display a positive relationship between biodiversity and ecosystem function (BEF) but the results have been mixed when assessed in bacterial communities. Habitat heterogeneity, a factor in eukaryotic BEFs, may explain these variable observations but it has not been thoroughly evaluated in bacterial communities. Here, we examined the impact of habitat on the relationship between diversity assessed based on the (phylogenetic) Hill diversity metrics and heterotrophic productivity. We sampled co-occurring free-living (more homogenous) and particle-associated (more heterogeneous) bacterial habitats in a freshwater, estuarine lake over three seasons: spring, summer and fall. There was a strong, positive, linear relationship between particle-associated bacterial richness and heterotrophic productivity that strengthened when considering dominant taxa. There were no observable BEF trends in free-living bacterial communities for any diversity metric. Biodiversity, richness and Inverse Simpson's index, were the best predictors of particle-associated production whereas pH was the best predictor of free-living production. Our findings show that heterotrophic productivity is positively correlated with the effective number of taxa and that BEF relationships are associated with microhabitats. These results add to the understanding of the highly distinct contributions to diversity and functioning contributed by bacteria in free-living and particle-associated habitats.

Microbiomes Reduce Their Host's Sensitivity to Interspecific Interactions.

Bacteria associated with eukaryotic hosts can affect host fitness and trophic interactions between eukaryotes, but the extent to which bacteria influence the eukaryotic species interactions within trophic levels that modulate biodiversity and species coexistence is mostly unknown. Here, we used phytoplankton, which are a classic model for evaluating interactions between species, grown with and without associated bacteria to test whether the bacteria alter the strength and type of species interactions within a trophic level. We demonstrate that host-associated bacteria alter host growth rates and carrying capacity. This did not change the type but frequently changed the strength of host interspecific interactions by facilitating host growth in the presence of an established species. These findings indicate that microbiomes can regulate their host species' interspecific interactions. As between-species interaction strength impacts their ability to coexist, our findings show that microbiomes have the potential to modulate eukaryotic species diversity and community composition.IMPORTANCE Description of the Earth's microbiota has recently undergone a phenomenal expansion that has challenged basic assumptions in many areas of biology, including hominid evolution, human gastrointestinal and neurodevelopmental disorders, and plant adaptation to climate change. By using the classic model system of freshwater phytoplankton that has been drawn upon for numerous foundational theories in ecology, we show that microbiomes, by facilitating their host population, can also influence between-species interactions among their eukaryotic hosts. Between-species interactions, including competition for resources, has been a central tenet in the field of ecology because of its implications for the diversity and composition of communities and how this in turn shapes ecosystem functioning.

Giant virus diversity and host interactions through global metagenomics.

Our current knowledge about nucleocytoplasmic large DNA viruses (NCLDVs) is largely derived from viral isolates that are co-cultivated with protists and algae. Here we reconstructed 2,074 NCLDV genomes from sampling sites across the globe by building on the rapidly increasing amount of publicly available metagenome data. This led to an 11-fold increase in phylogenetic diversity and a parallel 10-fold expansion in functional diversity. Analysis of 58,023 major capsid proteins from large and giant viruses using metagenomic data revealed the global distribution patterns and cosmopolitan nature of these viruses. The discovered viral genomes encoded a wide range of proteins with putative roles in photosynthesis and diverse substrate transport processes, indicating that host reprogramming is probably a common strategy in the NCLDVs. Furthermore, inferences of horizontal gene transfer connected viral lineages to diverse eukaryotic hosts. We anticipate that the global diversity of NCLDVs that we describe here will establish giant viruses-which are associated with most major eukaryotic lineages-as important players in ecosystems across Earth's biomes.

Randomized Lasso Links Microbial Taxa with Aquatic Functional Groups Inferred from Flow Cytometry.

High-nucleic-acid (HNA) and low-nucleic-acid (LNA) bacteria are two operational groups identified by flow cytometry (FCM) in aquatic systems. A number of reports have shown that HNA cell density correlates strongly with heterotrophic production, while LNA cell density does not. However, which taxa are specifically associated with these groups, and by extension, productivity has remained elusive. Here, we addressed this knowledge gap by using a machine learning-based variable selection approach that integrated FCM and 16S rRNA gene sequencing data collected from 14 freshwater lakes spanning a broad range in physicochemical conditions. There was a strong association between bacterial heterotrophic production and HNA absolute cell abundances (R 2 = 0.65), but not with the more abundant LNA cells. This solidifies findings, mainly from marine systems, that HNA and LNA bacteria could be considered separate functional groups, the former contributing a disproportionately large share of carbon cycling. Taxa selected by the models could predict HNA and LNA absolute cell abundances at all taxonomic levels. Selected operational taxonomic units (OTUs) ranged from low to high relative abundance and were mostly lake system specific (89.5% to 99.2%). A subset of selected OTUs was associated with both LNA and HNA groups (12.5% to 33.3%), suggesting either phenotypic plasticity or within-OTU genetic and physiological heterogeneity. These findings may lead to the identification of system-specific putative ecological indicators for heterotrophic productivity. Generally, our approach allows for the association of OTUs with specific functional groups in diverse ecosystems in order to improve our understanding of (microbial) biodiversity-ecosystem functioning relationships.IMPORTANCE A major goal in microbial ecology is to understand how microbial community structure influences ecosystem functioning. Various methods to directly associate bacterial taxa to functional groups in the environment are being developed. In this study, we applied machine learning methods to relate taxonomic data obtained from marker gene surveys to functional groups identified by flow cytometry. This allowed us to identify the taxa that are associated with heterotrophic productivity in freshwater lakes and indicated that the key contributors were highly system specific, regularly rare members of the community, and that some could possibly switch between being low and high contributors. Our approach provides a promising framework to identify taxa that contribute to ecosystem functioning and can be further developed to explore microbial contributions beyond heterotrophic production.

Genome evolution and host-microbiome shifts correspond with intraspecific niche divergence within harmful algal bloom-forming Microcystis aeruginosa.

Intraspecific niche divergence is an important driver of species range, population abundance and impacts on ecosystem functions. Genetic changes are the primary focus when studying intraspecific divergence; however, the role of ecological interactions, particularly host-microbiome symbioses, is receiving increased attention. The relative importance of these evolutionary and ecological mechanisms has seen only limited evaluation. To address this question, we used Microcystis aeruginosa, the globally distributed cyanobacterium that dominates freshwater harmful algal blooms. These blooms have been increasing in occurrence and intensity worldwide, causing major economic and ecological damages. We evaluated 46 isolates of M. aeruginosa and their microbiomes, collected from 14 lakes in Michigan, USA, that vary over 20-fold in phosphorus levels, the primary limiting nutrient in freshwater systems. Genomes of M. aeruginosa diverged along this phosphorus gradient in genomic architecture and protein functions. Fitness in low-phosphorus lakes corresponded with additional shifts within M. aeruginosa including genome-wide reductions in nitrogen use, an expansion of phosphorus assimilation genes and an alternative life history strategy of nonclonal colony formation. In addition to host shifts, despite culturing in common-garden conditions, host-microbiomes diverged along the gradient in taxonomy, but converged in function with evidence of metabolic interdependence between the host and its microbiome. Divergence corresponded with a physiological trade-off between fitness in low-phosphorus environments and growth rate in phosphorus-rich conditions. Co-occurrence of genotypes adapted to different nutrient environments in phosphorus-rich lakes may have critical implications for understanding how M. aeruginosa blooms persist after initial nutrient depletion. Ultimately, we demonstrate that the intertwined effects of genome evolution, host life history strategy and ecological interactions between a host and its microbiome correspond with an intraspecific niche shift with important implications for whole ecosystem function.

Gene Expansion and Positive Selection as Bacterial Adaptations to Oligotrophic Conditions.

We examined the genomic adaptations of prevalent bacterial taxa in a highly nutrient- and ion-depleted freshwater environment located in the secondary cooling water system of a nuclear research reactor. Using genome-centric metagenomics, we found that none of the prevalent bacterial taxa were related to typical freshwater bacterial lineages. We also did not identify strong signatures of genome streamlining, which has been shown to be one of the ecoevolutionary forces shaping the genome characteristics of bacterial taxa in nutrient-depleted environments. Instead, focusing on the dominant taxon, a novel Ramlibacter sp. which we propose to name Ramlibacter aquaticus, we detected extensive positive selection on genes involved in phosphorus and carbon scavenging pathways. These genes were involved in the high-affinity phosphate uptake and storage into polyphosphate granules, metabolism of nitrogen-rich organic matter, and carbon/energy storage into polyhydroxyalkanoate. In parallel, comparative genomics revealed a high number of paralogs and an accessory genome significantly enriched in environmental sensing pathways (i.e., chemotaxis and motility), suggesting extensive gene expansions in R. aquaticus The type strain of R. aquaticus (LMG 30558T) displayed optimal growth kinetics and productivity at low nutrient concentrations, as well as substantial cell size plasticity. Our findings with R. aquaticus LMG 30558T demonstrate that positive selection and gene expansions may represent successful adaptive strategies to oligotrophic environments that preserve high growth rates and cellular productivity.IMPORTANCE By combining a genome-centric metagenomic approach with a culture-based approach, we investigated the genomic adaptations of prevalent populations in an engineered oligotrophic freshwater system. We found evidence for widespread positive selection on genes involved in phosphorus and carbon scavenging pathways and for gene expansions in motility and environmental sensing to be important genomic adaptations of the abundant taxon in this system. In addition, microscopic and flow cytometric analysis of the first freshwater representative of this population (Ramlibacter aquaticus LMG 30558T) demonstrated phenotypic plasticity, possibly due to the metabolic versatility granted by its larger genome, to be a strategy to cope with nutrient limitation. Our study clearly demonstrates the need for the use of a broad set of genomic tools combined with culture-based physiological characterization assays to investigate and validate genomic adaptations.

Correction: Verrucomicrobia are prevalent in north-temperate freshwater lakes and display class-level preferences between lake habitats.

[This corrects the article DOI: 10.1371/journal.pone.0195112.].

Ecological Engineering Helps Maximize Function in Algal Oil Production.

Algal biofuels have the potential to curb the emissions of greenhouse gases from fossil fuels, but current growing methods fail to produce fuels that meet the multiple standards necessary for economical industrial use. For example, algae grown as monocultures for biofuel production have not simultaneously and economically achieved high yields of the high-quality lipid-rich biomass desired for the industrial-scale production of bio-oil. Decades of study in the field of ecology have demonstrated that simultaneous increases in multiple functions, such as the quantity and quality of biomass, can occur in natural ecosystems by increasing biological diversity. Here, we show that species consortia of algae can improve the production of bio-oil, which benefits from both a high biomass yield and a high quality of biomass rich in fatty acids. We explain the underlying causes of increased quantity and quality of algal biomass among species consortia by showing that, relative to monocultures, species consortia can differentially regulate lipid metabolism genes while growing to higher levels of biomass, in part due to a greater utilization of nutrient resources. We identify multiple genes involved in lipid biosynthesis that are frequently upregulated in bicultures and further show that these elevated levels of gene expression are highly predictive of the elevated levels in biculture relative to that in monoculture of multiple quality metrics of algal biomass. These results show that interactions between species can alter the expression of lipid metabolism genes and further demonstrate that our understanding of diversity-function relationships from natural ecosystems can be harnessed to improve the production of bio-oil.IMPORTANCE Algal biofuels are one of the more promising forms of renewable energy. In our study, we investigate whether ecological interactions between species of microalgae regulate two important factors in cultivation-the biomass of the crop produced and the quality of the biomass that is produced. We found that species interactions often improved production yields, especially the fatty acid content of the algal biomass, and that differentially expressed genes involved in fatty acid metabolism are predictive of improved quality metrics of bio-oil. Other studies have found that diversity often improves productivity and stability in agricultural and natural ecosystems. Our results provide further evidence that growing multispecies crops of microalgae may improve the production of high-quality biomass for bio-oil.