Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

The Trait Repertoire Enabling Cyanobacteria to Bloom Assessed through Comparative Genomic Complexity and Metatranscriptomics

Huansheng Cao, Yohei Shimura, Morgan M. Steffen, Zhou Yang, Jingrang Lu, Allen Joel, Landon Jenkins, Masanobu Kawachi, Yanbin Yin, Ferran Garcia-Pichel
Stephen J. Giovannoni, Editor
Huansheng Cao
aDepartment of Biological Sciences, Northern Illinois University, DeKalb, Illinois, USA
bBiodesign Center for Fundamental and Applied Microbiomics, Arizona State University, Tempe, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Huansheng Cao
Yohei Shimura
cNational Institute for Environmental Studies, Tsukuba, Ibaraki, Japan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Morgan M. Steffen
dBiology Department, James Madison University, Harrisonburg, Harrisonburg, Virginia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhou Yang
eJiangsu Key Laboratory for Biodiversity and Biotechnology, School of Biological Sciences, Nanjing Normal University, Nanjing, Jiangsu, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jingrang Lu
fU.S. Environmental Protection Agency Office of Research and Development, Cincinnati, Ohio, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Allen Joel
fU.S. Environmental Protection Agency Office of Research and Development, Cincinnati, Ohio, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Landon Jenkins
bBiodesign Center for Fundamental and Applied Microbiomics, Arizona State University, Tempe, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Masanobu Kawachi
cNational Institute for Environmental Studies, Tsukuba, Ibaraki, Japan
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yanbin Yin
aDepartment of Biological Sciences, Northern Illinois University, DeKalb, Illinois, USA
gNebraska Food for Health Center, Department of Food Science and Technology, University of Nebraska-Lincoln, Lincoln, Nebraska, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Yanbin Yin
Ferran Garcia-Pichel
bBiodesign Center for Fundamental and Applied Microbiomics, Arizona State University, Tempe, Arizona, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen J. Giovannoni
Oregon State University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.01155-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Water bloom development due to eutrophication constitutes a case of niche specialization among planktonic cyanobacteria, but the genomic repertoire allowing bloom formation in only some species has not been fully characterized. We posited that the habitat relevance of a trait begets its underlying genomic complexity, so that traits within the repertoire would be differentially more complex in species successfully thriving in that habitat than in close species that cannot. To test this for the case of bloom-forming cyanobacteria, we curated 17 potentially relevant query metabolic pathways and five core pathways selected according to existing ecophysiological literature. The available 113 genomes were split into those of blooming (45) or nonblooming (68) strains, and an index of genomic complexity for each strain’s version of each pathway was derived. We show that strain versions of all query pathways were significantly more complex in bloomers, with complexity in fact correlating positively with strain blooming incidence in 14 of those pathways. Five core pathways, relevant everywhere, showed no differential complexity or correlations. Gas vesicle, toxin and fatty acid synthesis, amino acid uptake, and C, N, and S acquisition systems were most strikingly relevant in the blooming repertoire. Further, we validated our findings using metagenomic gene expression analyses of blooming and nonblooming cyanobacteria in natural settings, where pathways in the repertoire were differentially overexpressed according to their relative complexity in bloomers, but not in nonbloomers. We expect that this approach may find applications to other habitats and organismal groups.

IMPORTANCE We pragmatically delineate the trait repertoire that enables organismal niche specialization. We based our approach on the tenet, derived from evolutionary and complex-system considerations, that genomic units that can significantly contribute to fitness in a certain habitat will be comparatively more complex in organisms specialized to that habitat than their genomic homologs found in organisms from other habitats. We tested this in cyanobacteria forming harmful water blooms, for which decades-long efforts in ecological physiology and genomics exist. Our results essentially confirm that genomics and ecology can be linked through comparative complexity analyses, providing a tool that should be of general applicability for any group of organisms and any habitat, and enabling the posing of grounded hypotheses regarding the ecogenomic basis for diversification.

INTRODUCTION

Functional niche specialization is a ubiquitous trait in microorganisms, a fact that is reflected in their trait and genomic repertoires (1). The presence of a given trait results from positive selection by environmental challenges typical of a particular habitat. In contrast, selection will tend to be negative for traits that do not contribute much to the fitness in the habitat under consideration, since they will be costly to encode and express for no or minimal benefit (2). Two factors determine through selection if a trait will be a part of the specialized repertoire addressing a given environmental challenge: the fitness contribution it provides during challenge and the cost of encoding and expressing the genomic units on which it is based (3). For example, the ability to synthesize highly effective nitrogenous compatible solutes in the adaptation of cyanobacteria to hypersaline media only occurs in extremely halophilic forms (4), whereas moderately halotolerant species settle for the use of disaccharides. However, this dual black-and-white distinction is unlikely to be always patent in nature, because habitats are variable in space and time, and because traits can be of relevance in more than one environment. To follow on our example, the use and expression of disaccharide-compatible solutes also confer significant fitness to soil cyanobacteria (5) even though most soils become hyperosmotic only when severely dry.

Hence, genomic systems may provide variable fitness depending on the frequency and preponderance of an environmental challenge in different environments. Under this framework, it is also logical to hypothesize that pathways or genomic systems that provide specific adaptive solutions to an environmental challenge will be comparatively more complex or sophisticated in species that are often challenged than in those living in habitats where the system brings only marginal or occasional fitness. This is because the increased cost of system maintenance of a complex system will be easily offset by the fitness benefits of a fine-tuned, effectively functioning configuration—if it is used extensively. A high-maintenance, complex system may not be affordable when it is only necessary occasionally, such that organisms settle for less effective or nimble, but simpler, alternatives through adaptive gene loss, or even outsourcing parts of a pathway to cooccurring organisms (6). This hypothesis is consistent with theoretical considerations in systems biology regarding the robustness versus failure of complex evolvable systems (7), the distinction between explicit versus “underground” metabolism (8), and the tendency to add complexity to existing nodes during system growth in metabolic network theory (9). One prominent example is the cellulosome in specialized cellulolytic bacteria (e.g., Clostridium thermocellum), a multienzyme complex with as many as 160 enzymes (10), where each enzyme also has multiple paralog genes (11). Other organisms (yeast, for example) that can grow on cellulose albeit inefficiently, make do with an unorganized subset of these genes (12). Similarly, uptake systems become demonstrably more complex under nutrient limitation (13), as do antibiotic resistance genes in the continued presence of antibiotics (14). Such relationships between complexity and functional relevance should be principally quantifiable.

To test this, we hypothesized that planktonic cyanobacteria that form blooms have and express specific pathways useful with the constraints of life in a bloom, and that these pathways will be more complex than their cognates in other nonblooming forms. We chose this example to study the functional niche specialization of water bloom-forming cyanobacteria at the systems level, because their niche is well defined among the many occupied by cyanobacteria, and only some species are able to form these blooms, even among a much larger diversity of extant planktonic forms (15). Also, cyanobacterial blooms (often referred to as CyanoHABs, for cyano-harmful algal blooms) have been well studied at the ecological, physiological, and genomic levels, so there is a considerable database of information that can help determine correlations of the complexity of certain genomic traits with their known or suspected differential relevance for bloom formation. These blooms are primarily driven by elevated nutrient loading from human activities (16, 17), typically occurring within a short period of time and particularly in waters of long residence time (18), which leaves no time for within-bloom evolution of new functions given the low mutation rate of cyanobacteria (19). Their consistent dominance in eutrophic oceans and inland waters has now expanded to a global scale of increasingly eutrophic waters (20, 21), suggesting that their genomes must be preequipped with effective functional capabilities to take advantage of the elevated resources that were not available in the preeutrophication era.

Some of the important traits for bloom formers have been inferred either through field studies or through nutrient manipulation in laboratory settings. Functions held to be important include those related to resource utilization such as buoyancy generation encoded by a gvp gene cluster (22), utilization of macronutrients (carbon, nitrogen, phosphate, and sulfur), trace elements (17, 23), and small-molecule organics, such as amino acid and simple sugars (24). Stress-resistance pathways are sometimes deemed important, including biosynthesis of ultraviolet sunscreens (25), toxin production (16), heavy metals (26, 27), or antibiotic resistance (28). However, further evaluation is needed because, first, the general occurrence of these pathways in bloomers has not been systematically evaluated, but rather inferred from one or few species (17). Second, metabolomic context and contrast with nonblooming controls are typically lacking (23). Despite these challenges, establishing the consensus of physiology for bloom formation remains a priority for effective bloom control (29).

Here, we apply a comparative genomic perspective derived from our “significance begets complexity” hypothesis, to address this goal using correlational analyses. We categorized 113 cyanobacterial strain genomes according to their tendency to form blooms. We then reconstructed pathways of interest, characterized their complexity pragmatically in each genome using coding lengths and number of identifiable genetic components, and probed the relationships between the two parameters statistically. We then used meta-transcriptomics of natural cyanobacterial planktonic populations to validate our results.

RESULTS

Phylogeny, genomic features, and morphology across the data set.The 113 genomes were initially classified as bloomers (45 genomes) and nonbloomers (68 genomes) (Table S1 in the supplemental material). As shown in Fig. 1 and Fig. S1, both blooming and nonblooming genomes appeared recurrently within the same phylogenetic cluster, such as clusters 1 and 2 of Fig. 1. With respect to morphology, the set of bloomers included unicellular, colonial, and filamentous cyanobacteria. Specifically, bloomers and nonbloomers from our set had representatives in groups I, III, and IV, although the two representatives of group II were both bloomers (Fig. 1). Besides similarities in phylogeny and morphology, the sets of bloomers and nonbloomers also had similar genomic features, such as genome size, gene content or GC content (two-sample Wilcoxon test, P > 0.05), and there were no apparent correlations between their ability to form blooms (Blooming Incident Index, see below) and general genomic features, although bloomers had a narrower range (Fig. S2). All of this supports the notion that our comparative data set did not significantly bias according to phylogeny, morphology, or rough genomic traits.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Phylogeny of the 43 cyanobacterial strains used based on 16S rRNA gene sequences. Blooming species are in blue and nonblooming in red. Morphological types are indicated according to the traditional Rippka groups.

FIG S1

Phylogeny of the 43 cyanobacterial strains used based on 16S rRNA sequences. Download FIG S1, 1.0 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S2

The relationship between bloom incidence index (BII) and cell shape and genomic features in 113 cyanobacterial strains. Download FIG S2, PDF file, 0.2 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S1

The 113 strains used in this study. Download Table S1, DOC file, 0.4 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Overview of the curated ecophysiological pathways in cyanobacteria.We reconstructed 5 core metabolic pathways and 17 query pathways derived from existing ecological tests in the literature (Table S2). Core pathways were used as controls because, if our hypothesis holds, they should not show differential complexity according to habitat. Core pathways included glycolysis, TCA cycle (TCA), pentose phosphate pathway (PPP), the electron transport chain (ETC), and the Calvin cycle (Calvin). These are typically of high similarity and low variability among microbial strains (30). Two modules of photosynthesis—PBS (phycobilisome) and CCM (CO2-concentrating mechanism)—were involved in utilization of external resources (light and CO2), and thus were formally included as query pathways in the comparison between bloomers and nonbloomers, even though one could also argue they are core pathways for phototrophs. The other query pathways selected have been shown to promote cyanobacterial growth and bloom formation (see references in Text S1 the supplemental material; Table S2). We illustrate the query pathways for the case of Aphanizomenon flos-aquae strain NIES-81 (Fig. 2). More detailed illustrations of all core and query pathways are provided in Fig. S3. A. flos-aquae NIES-81 has pathways for the assimilation of carbon (CCM), nitrogen (N), phosphorus (P), sulfur (S), free amino acids (AAT), and sugars, as well as trace elements and vitamins (TEVit). It has capabilities against stressors like heavy metals (MetalR), antibiotics (DrugR), oxidative stress (OSR), UV radiation (MAA), osmotic pressure, and possibly predation (through toxin production). Additionally, buoyancy can be regulated through the biosynthesis of gas vesicles. These pathways had different composition (number of components, e.g., transporters or enzymes) and complexity (components consisting of multiple subunits) among different strains.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

The core and query pathways in Aphanizomenon flos-aquae NIES-81. Each component is detailed in the corresponding tables (Tables S2 to S24). The pathways labeled with dashed borders are not complete (due to the absence of required components) in this strain but may be complete in others.

TEXT S1

Supplemental methods and references. Download Text S1, DOCX file, 0.03 MB.
Copyright © 2020 Cao et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S3

Core and query pathway in this study. Download FIG S3, PDF file, 1.7 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

The classification of central and query pathways and their protein reference sources. Download Table S2, DOC file, 0.3 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Complexity difference in pathways between bloomers and nonbloomers.On average, bloomers had a total of 461 (±57) proteins in the 19 query pathways, accounting for about 14.5% of their entire genomes, whereas nonbloomers encoded only 274 (±42) proteins, accounting for 10.3% of the genes, suggesting there is indeed a tendency for more complex query pathways among bloomers. As a control, core pathways had on average similar number of proteins between bloomers and nonbloomers, 169 versus 160, accounting for 5.3% and 6% of the total genes, respectively. The specific components of the core and query pathways for each strain in our analyses can be found as Tables S3 to S16 at https://github.com/hshcao/CHABpathways.

More specifically, we present the ratio between bloomers and nonbloomers of each of our eight genomic complexity metrics in Table 1. The metrics considered the number of components in each pathway, component complexity (e.g., multisubunits of enzyme complexes or ABC transporters), and the proportion of encoding sequence of the entire genome for these components (see the Materials and Methods). Again here, query pathways (except those related to photosynthesis) were significantly (Wilcoxon sign-rank test, P < 0.05) more complex among bloomers regardless of metric. However, the magnitude of the difference varied considerably with pathway, from almost 7-fold in some (vesicle) to marginal values in others (ORS). The ratios were quite consistent among metrics. As expected, the core metabolic pathways did not differ significantly between bloomers and nonbloomers for any complexity metric. PSII and PSI were only marginally different, however, perhaps because they should be considered core pathways among phototrophs. Table 1 also includes the ratios of a general genomic complexity index (GCI, computed as an average of the eight individual complexity metrics). Clearly, the vesicle pathway had the highest GCI ratio (6.6), followed by toxin, PUFA (polyunsaturated fatty acids), osmosis, AAT, PBS, MAA, nitrogen fixation (NF), and CCM. It is noteworthy that none of the nonbloomers had a complete set of genes for toxin production, so some of the metrics could not be calculated (with zeros as denominator), indicating that our comparison in this case is an underestimate. In Fig. 3 we show the actual distribution of GCI for bloomers and nonbloomers. According to this integrated comparison, none of the core pathways (including PSI and PSII) were different. Among the query pathways, all were significantly more complex among bloomers.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Ratio of complexity metrics of core and query pathways between blooming and nonblooming strains

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Average relative complexity (GCI) in blooming and nonblooming strains for each of the 24 pathways. B, blooming; NB, nonblooming. The background color of each plot represents the P values of the comparison between blooming and nonblooming strains, based on Wilcoxon sign-rank tests.

Because some species were overrepresented with multiple strain genomes in our data set, such as Microcystis aeruginosa, Planktothrix agardhii, Cyanothece sp., Synechococcus sp., Crocosphaera watsonii, Prochlorococcus marinus and Prochlorococcus sp., we carried out additionally analyses with data sets that used only one random genome for each species. A similar pattern in complexity differences was still observed (Fig. S4).

FIG S4

Comparison of the 24 pathways in terms of eight metrics between blooming and nonblooming strains without overrepresented species. Download FIG S4, PDF file, 0.2 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The pattern shown in Fig. 3 may potentially be artifactual due to the presence of evolutionary/phylogenetic legacies, given that the set of strains was not randomly distributed in the cyanobacteria phylogeny, but rather tended to aggregate in certain clades. To test for such an effect, we selected two clusters characterized by different morphology and blooming types (clusters 1 and 2 in Fig. 1) and performed the comparison between blooming and nonblooming species within each. In these analyses, blooming species still had higher metrics (Wilcoxon sign-rank test, all P values < 10−2) even when determined among closely related phylogenetic strains. We conclude that any effects of evolutionary/phylogenetic legacies were insignificant.

We then probed in a more quantitative way the relevance of the query pathways to bloom formation by correlation, because some of the strains are not always exclusively bloomers or nonbloomers and may be referenced in the literature under either class. For each pathway we ran linear correlations between the GCI and a bloom incidence index (BII; derived from the literature, see the Materials and Methods). Expectedly, we found no significant correlation (P > 0.05) for any of the core pathways or PSI and PSII (Fig. 4; see detailed plot in Fig. S5). Among the query pathways, correlations were significant for all pathways except for OSR, DrugR, and MetalR. However, only in 7 pathways (vesicle, toxin, PUFA and AAT, CCM, sulfur, and osmosis) were the correlation coefficients large (R2 >0.4), and only 2 additional pathways had R2 >0.25 (NF and N). Interestingly, besides the significance of the correlations, the slopes of the correlation also showed a similar pattern among the query pathways (Fig. 4B). Vesicle and Calvin were shown as examples of significant correlation (Fig. 4C) and of no correlation (Fig. 4D), respectively. We also calculated the correlation between BII and PCI of all genomes, i.e., a test of the effect of genome redundancy, and a similar result was found to that without redundancy (Fig. S6).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Bubble plots of Pearson’s correlations between GCI and BII in each pathway. The circle size represents either the coefficient of determination (R2) (A) or the slope of the correlation (B) with the significance (P) of the correlation indicated with colors. Two correlations are also shown as examples: Vesicle (C) and Calvin cycle (D).

FIG S5

Pearson’s correlations between GCI and BII in each pathway. Download FIG S5, PDF file, 2.0 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S6

Correlation between the query pathway PCI and BII of all 113 genomes. Download FIG S6, PDF file, 0.5 MB.
Copyright © 2020 Cao et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To illustrate in detail how the query pathways identified above were indeed enriched with more components, we visualized the complete set of components in the significant vesicle pathway as opposed to the nonsignificant OSR, as shown in Fig. 5. For the vesicle pathway genes, all genes were more enriched in bloomers than in nonbloomers (Wilcoxon sign-rank test, P < 0.05). For the OSR genes, only the gene gpx was significantly higher in nonbloomers than bloomers (Wilcoxon sign-rank test, all P values < 10−2), while others were not different (Wilcoxon sign-rank test, P > 0.05).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Gas vesicle (Vesicle) (A) and oxidative stress resistance (OSR) (B) genes present in the 113 genomes studied, sorted by blooming capacity.

Correlation between GCI and the expression level of the query pathways in eutrophic and oligotrophic waters.Finally, to test if the bloom-enabling genetic repertoire determined above was indeed relevant during bloom formation, we tested the prediction that the encoded genetic elements or pathways in the repertoire must be used in proportion to their relevance (and hence complexity). To that end, ecologically relevant metatranscriptomes of our own from eutrophic freshwater blooms in Lake Erie and Harsha Lake were studied in contrast to those of oligotrophic freshwater and oligotrophic ocean plankton which were publicly available. RNAseq reads were mapped to the genomes of the species closest to those dominating each habitat (M. aeruginosa strain NIES-843 for the bloom, Oscillatoria nigro-viridis strain PCC7112 for the oligotrophic freshwater plankton, and Synechococcus sp. strain WH8102 for the oligotrophic marine plankton). Significant correlations between the expression levels of pathways and the pathway GCIs were only found in blooming species in two eutrophic lakes (R2 = 0.47 and 0.38; both P < 0.01) (Fig. 6A to D), but not in nonbloomers in oligotrophic waters (Fig. 6E to H). Particularly, vesicle, PBS, CCM, and N were among those pathways highly expressed in blooming species, while nonsignificant pathways, such OSR, metalR, and even MAA, were expressed at the lowest levels. Interestingly, core pathways were expressed at moderate levels relative to the entire range of query pathway levels.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Correlations between the expression levels and GCI of the pathways studied in three different aquatic habitats. Metatranscriptomes were from eutrophic freshwater (Lake Erie, A and B; Harsha Lake, C and D), oligotrophic freshwater (Sparkling Lake, E and F), oligotrophic ocean water (western Arctic Ocean, Canada, G and H). Reference genomes representing the dominant species in these habitats were used for mapping: M. aeruginosa NIES-843 (Lake Erie and Harsha Lake), Oscillatoria nigro-viridis PCC7112 (Sparkling Lake), and Synechococcus sp. WH8102 (western Arctic Ocean, Canada). Pearson correlations were performed with (A, C, E, and G) or without (B, D, F, and H) core pathways. The coefficient of determination (R2) and the significance (P) of the correlations are on top left corner of each panel. GCI values were subject to angular transformation.

DISCUSSION

A key understudied question in CyanoHAB research is whether there exists a common ecophysiology underlying the broad array of blooming species. We attempted to address this question from the comparative functional genomics perspective, testing the hypothesis that functional niche specialization will result in preferential sophistication of genomic systems that are differentially important for bloom formers. Such differential importance must be manifest in both encoding and expressing these functions. An initial comparison of pathway complexity among bloomers and nonbloomers found that pathway complexity (i.e., more coding investment) was statistically higher among bloomers for most potentially relevant pathways queried (and not for core pathways or core photosynthetic pathways), but also that the effect was not commensurate among pathways. The fact that none of the core metabolism was significantly different agrees with the general pattern of cyanobacterial genome evolution, i.e., the core set of genes is well conserved while accessory pathways are variable (30). The potential relevance of pathways for bloomers was mirrored in the level of expression of the query pathways in eutrophic waters but not in nonblooming species in oligotrophic waters. These results set the stage for direct experimental interrogation of the specific roles and mechanisms of each pathway.

Pathways most enriched among blooming cyanobacteria.We identified vesicle formation for buoyancy as the most enriched and highly expressed pathway, consistent with the ecological literature. Vesicle formation has been well established as part of the common ecophysiological repertoire among various types of bloomers, conferring buoyancy to adjust their position in the water column for optimal resource utilization (17, 23, 31). The toxin production pathway was also strongly enriched. While the exact role or benefit that bloomers derive from toxin production is still a matter of debate (32, 33), our data confirm their differential relevance for CyanoHABs. Their low expression may verify their potential role as a cell-cell signal molecule which is not required in large quantities (34). The relevance of these vesicle and toxin pathways found ample confirmation in how the complexity of these pathways best correlated with BII, and their level of expression in nature was commensurate with its high GCI.

Two pathways that saw significant enrichment among bloomers and a high correlation with BII, were the amino acid uptake systems, and the synthesis of polyunsaturated fatty acids. Cyanobacteria have been known to utilize organic nutrients such as glucose, acetyl-glucosamine, and various amino acids (24, 35), but the ecological roles of these organic nutrients have been little studied in bloom formation. As the level of organic matter, particularly free amino acids, increases in eutrophic waters (36–38), this pathway apparently becomes more important in driving or sustaining blooms, particularly in those cyanobacteria capable of mixotrophy (39). Recovering leaked amino acids from the exometabolome (40) becomes more feasible when cyanobacterial populations attain large density, as they do in blooms, since diffusional losses are counteracted by the high concentration of sources, unlike what would happen with dilute planktonic populations. The significance of PUFA, typically related to an ability to modify membrane fluidity in adaptation to temperature (41, 42), was surprising. Clearly temperature and its fluctuations are an important factor in lakes, but it is not clear to us how they are differentially more relevant than in other habitats. One possibility is that unsaturated fatty acids may help protect photosynthetic machinery in the presence of salt and other stresses (43, 44). It seems that increased attention should be paid to this issue in future research. The prominence of osmosis among the most responsive pathways was also rather surprising, given than salinity fluctuations are not among the most prominent traits in lakes, although they may be important for coastal or estuarine blooms.

Role of macronutrients.The role of macronutrient (C, N, P, or S) acquisition for bloom formation also found confirmation in our study. Particularly for C, N, and S, all correlated rather well with BII and were expressed at levels commensurate with their GCIs. For C, a diffusion limitation effect associated with high population density may enhance the value of sustaining complex carbon concentrating mechanisms (CCM) among bloomers, because heavy drawdown of CO2 can easily result in C limitation. In agreement with our results, rising CO2 level has been found to intensify phytoplankton blooms in eutrophic waters (45), and transporters and regulators of CCM seem very responsive to the changes in dissolved CO2 levels during a dense CyanoHAB (46). Interestingly, the correlation of P with BII was weak and its complexity differences between bloomers and nonbloomers were low (while both significant). According to our results, N and S acquisition would much more strongly define the bloom niche. Strikingly, field meta-analyses of over 1,000 lakes identified total nitrogen (TN), not total phosphorus, as the top nutrient for bloom formation (47). One possible reason for this is that bloom-forming and N-fixing cyanobacteria have a higher N requirement than other cyanobacteria due to their need for protein-based antenna pigments (phycobilisomes), and thus a higher N:P ratio (48) than our set of nonbloomers, in which some strains lack them (Prochlorococcus sp. and Prochlorothrix). But an ad hoc comparison excluding the latter did not support this notion, in that the significance remained below 0.01. In any event, this is not to say that P is not important. The even distribution of genes and particularly the fact that some blooming species have two sets of phosphate ABC transporters (pstABCD) suggests that P may simply not be more significant in blooms that in other habitats. In agreement with this view, external P additions do not necessarily lead to higher Microcystis abundance in Lake Erie where phosphorus level is already high (49), and the expression levels of P are similar to core pathways in all three habitats tested.

Pathways of minor differential relevance.Drug and metal resistance pathways had low correlations with BII and only minor increases in complexity among bloomers, and also showed lowest expression among the query pathways.

The case of MAA sunscreens.With respect to the water-soluble mycosporine pathway, which is common among cyanobacteria of many origins (50), we would like to note that increasing pathway length in this system is associated with the synthesis of mycosporines absorbing at longer wavelengths (from the UVB into the UVA) (51). It is possible that this fact may have contributed spuriously to enhance complexity in bloomers to reach a significance with no environmental meaning, because UVB radiation does not penetrate much in natural waters. This would tend to favor preferentially the use of (genomically complex) bi-substituted iminomycosporines among bloomers. Hence, the significance of MAA differences in complexity, and its correlation with BII, should be taken with caution despite that they were not among the most prominent differences.

A tool of general applicability?More generally, the data presented here provide support for the “relevance begets complexity” in comparative genomics. It delineates a tool to bring together comparative genomics and ecology. While this approach was here applied to the case of bloom-forming cyanobacteria, there is in principle no reason why this approach may not find application in a variety of other microbial habitats and organismal groups. On a broader scale, our approach can be combined with metagenome-assembled genomes to infer the link between genomic complexity and ecological niche, where together they could provide a comprehensive conceptual framework for predictive ecology from genomic data (52).

MATERIALS AND METHODS

Cyanobacterial genome sequences and strain information.Whole-genome sequences of 113 cyanobacteria (Table S1) were downloaded from NCBI in May 2016. These cyanobacteria belong to 45 species. We initially classified them into blooming and nonblooming species based on a curated classification (16, 17, 53). The species not listed in the original curation were classified by querying their species/genus names in the Web of Science; if publications reporting blooms formation were found, they were classified as blooming and otherwise as nonblooming (Table S1).

A more precise ranking followed, in which we examined the blooming incidence index (BII) of all strains in the Web of Science. We used as query either the strain denomination (first), species epithet (if the strain denomination yielded insufficient hits), or, in a few cases, the genus name. The records returned were then further separated into genomic/physiological or environmental studies using appropriate key words. Finally, the records were classified into those that reported a bloom and those that did not. The BII for each strain was simply the proportion of environmental citations that mentioned blooming capacity, and it varied between 0 and 1 (Table S1). We discarded strains in genera that had <2 environmental citations. For convenience, and to check for phylogeny-derived patterns that may confound our search, we also classified the strains into four morpho-typical groups, according to the traditional divisions (54), as well as phylogenetically on the basis of 16S rRNA gene sequences (see details in the Methods in the supplemental material).

Reference proteins of central and query pathways.For core bacterial metabolic pathways, we included glycolysis, the pentose phosphate pathway, the TCA cycle, the Calvin cycle, and the electron transport chain. Additionally, we included four photosynthetic modules: photosystems I and II, phycobilisome, and CO2-concentrating mechanism. The reference protein sequences for the central metabolism were extracted from the SEED subsystems (55) and the KEGG database (56). The rest of the pathways (query pathways), which were presumably involved in differential resource utilization and stress resistance among bloomers, were manually curated from the literature and databases such as the SEED subsystems, KEGG, and the transporter database TCDB (57) (see the pathways in Table S2 and relevant citations in the Methods in the supplemental material).

Homolog search in the cyanobacterial genomes.With the reference proteins as the query, we searched against the 113 cyanobacterial genomes using BLASTP (58). We performed the search first with cyanobacteria-derived protein sequences if available; if not, heterotrophic bacteria-derived proteins were used. To ensure significant hits from BLASTP output, we used different E value and coverage thresholds for central pathways (E value of <10−12 and length coverage of >80%) and query pathways (E value of <10−7 and length coverage of >65%), since proteins in central pathways are more conserved than those in the query pathways (30). Among the 113 genomes, Aphanizomenon flos-aquae NIES-81 (59) was selected to illustrate central and query pathways.

Assessing relevance of query pathways to CyanoHAB formation.To quantify the comparative differences in complexity of any given pathway, we devised eight metrics that considered both the pathway composition and completeness (Table 1). The eight were as follows: NP (total number of proteins in each pathway and strain), NC (total number of protein complexes (including multiprotein complexes and single proteins per pathway and strain), NPC (total number of proteins in the complete multiprotein complexes per pathway and strain), LP (total base-pair length of nucleotide sequences encoding all the proteins in a pathway per strain), FLP (fraction of total coding length of proteins or the ratio of LP to strain’s genome size), LC (total base-pair length of proteins in the complete multiprotein complexes per strain and pathway), FLC (fraction of total coding length of proteins in the complete multiprotein complexes or the ratio of LC to the strain’s genome size), and PC (the ratio of NC to the total number of complexes in the reference protein set in each pathway per strain). Each metric was normalized so the lowest value found was set to 0 and the maximal value found in the set was set to 1. Once normalized, the 8 metrics of each strain and pathway were averaged to attain an overall genomic complexity index (GCI) for a pathway and strain that varied between 0 and 1. For statistical treatments, since GCI and BII are proportions, they were all subject to angular transformation.

Correlation between GCI and the expression level of the query pathways.We tested the correlation with metatranscrioptome of cyanobacteria collected from three types of water: (i) our own metagenomes from a CyanoHAB in eutrophic Lake Erie; (ii) oligotrophic freshwater Sparkling Lake; and (iii) the oligotrophic Arctic Ocean. The reads from each water body were mapped to reference genomes closest to their own dominant cyanobacterial species, M. aeruginosa NIES-843, Oscillatoria nigro-viridis PCC7112, and Synechococcus sp WH8102, respectively (see details in the Methods in the supplementary material). Next, the mean expression levels of each pathway were calculated as the average level of all the genes in the pathways. For each water body, the correlations between mean pathway (log transformed) expression levels and (angular transformed) BIIs were carried out for all (query and core) or only query pathways, using R.

ACKNOWLEDGMENTS

This project was partly supported by the National Natural Science Foundation of China (31730105), by a startup grant from Northern Illinois University, and by a grant from the U.S. National Science Foundation (DBI-1933521).

FOOTNOTES

    • Received 12 May 2020
    • Accepted 26 May 2020
    • Published 30 June 2020
  • Copyright © 2020 Cao et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Sagane Y,
    2. Zech K,
    3. Bouquet J-M,
    4. Schmid M,
    5. Bal U,
    6. Thompson EM
    . 2010. Functional specialization of cellulose synthase genes of prokaryotic origin in chordate larvaceans. Development 137:1483–1492. doi:10.1242/dev.044503.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Vamosi JC,
    2. Poisot T
    . 2016. Evolutionary ecology of specialisation. https://onlinelibrary.wiley.com/doi/full/10.1002/9780470015902.a0026281.
  3. 3.↵
    1. Roop JI,
    2. Chang KC,
    3. Brem RB
    . 2016. Polygenic evolution of a sugar specialization trade-off in yeast. Nature 530:336–339. doi:10.1038/nature16938.
    OpenUrlCrossRef
  4. 4.↵
    1. Abed RM,
    2. Garcia-Pichel F,
    3. Hernández-Mariné M
    . 2002. Polyphasic characterization of benthic, moderately halophilic, moderately thermophilic cyanobacteria with very thin trichomes and the proposal of Halomicronema excentricum gen. nov., sp. nov. Arch Microbiol 177:361–370. doi:10.1007/s00203-001-0390-2.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Baran R,
    2. Lau R,
    3. Bowen BP,
    4. Diamond S,
    5. Jose N,
    6. Garcia-Pichel F,
    7. Northen TR
    . 2017. Extensive turnover of compatible solutes in cyanobacteria revealed by deuterium oxide (D2O) stable isotope probing. ACS Chem Biol 12:674–681. doi:10.1021/acschembio.6b00890.
    OpenUrlCrossRef
  6. 6.↵
    1. Morris JJ,
    2. Lenski RE,
    3. Zinser ER
    . 2012. The black queen hypothesis: evolution of dependencies through adaptive gene loss. mBio 3 doi:10.1128/mBio.00036-12.
    OpenUrlCrossRef
  7. 7.↵
    1. Kitano H
    . 2004. Biological robustness. Nat Rev Genet 5:826–837. doi:10.1038/nrg1471.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Notebaart RA,
    2. Szappanos B,
    3. Kintses B,
    4. Pál F,
    5. Györkei Á,
    6. Bogos B,
    7. Lázár V,
    8. Spohn R,
    9. Csörgő B,
    10. Wagner A,
    11. Ruppin E,
    12. Pál C,
    13. Papp B
    . 2014. Network-level architecture and the evolutionary potential of underground metabolism. Proc Natl Acad Sci U S A 111:11762–11767. doi:10.1073/pnas.1406102111.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Barabási A-L,
    2. Oltvai ZN
    . 2004. Network biology: understanding the cell’s functional organization. Nat Rev Genet 5:101–113. doi:10.1038/nrg1272.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Artzi L,
    2. Bayer EA,
    3. Moraïs S
    . 2017. Cellulosomes: bacterial nanomachines for dismantling plant polysaccharides. Nat Rev Microbiol 15:83–95. doi:10.1038/nrmicro.2016.164.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Tamaru Y,
    2. Miyake H,
    3. Kuroda K,
    4. Nakanishi A,
    5. Matsushima C,
    6. Doi RH,
    7. Ueda M
    . 2011. Comparison of the mesophilic cellulosome-producing Clostridium cellulovorans genome with other cellulosome-related clostridial genomes. Microb Biotechnol 4:64–73. doi:10.1111/j.1751-7915.2010.00210.x.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Zhao Z,
    2. Liu H,
    3. Wang C,
    4. Xu J-R
    . 2013. Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics 14:274. doi:10.1186/1471-2164-14-274.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Gresham D,
    2. Desai MM,
    3. Tucker CM,
    4. Jenq HT,
    5. Pai DA,
    6. Ward A,
    7. DeSevo CG,
    8. Botstein D,
    9. Dunham MJ
    . 2008. The repertoire and dynamics of evolutionary adaptations to controlled nutrient-limited environments in yeast. PLoS Genet 4:e1000303. doi:10.1371/journal.pgen.1000303.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Edlund T,
    2. Normark S
    . 1981. Recombination between short DNA homologies causes tandem duplication. Nature 292:269–271. doi:10.1038/292269a0.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Garcia-Pichel F
    . 2000. Cyanobacteria, p 907–927. In Lederberg L (ed), Encyclopedia of microbiology. Academic Press, San Diego, CA.
  16. 16.↵
    1. Huisman J,
    2. Codd GA,
    3. Paerl HW,
    4. Ibelings BW,
    5. Verspagen JMH,
    6. Visser PM
    . 2018. Cyanobacterial blooms. Nat Rev Microbiol 16:471–483. doi:10.1038/s41579-018-0040-1.
    OpenUrlCrossRef
  17. 17.↵
    1. Oliver R,
    2. Ganf G
    . 2002. Freshwater blooms, p. In Whitton B, Potts M (ed), The ecology of cyanobacteria, Springer, The Netherlands.
  18. 18.↵
    1. Steffen MM,
    2. Belisle BS,
    3. Watson SB,
    4. Boyer GL,
    5. Wilhelm SW
    . 2014. Status, causes and controls of cyanobacterial blooms in Lake Erie. J Great Lakes Res 40:215–225. doi:10.1016/j.jglr.2013.12.012.
    OpenUrlCrossRef
  19. 19.↵
    1. Osburne MS,
    2. Holmbeck BM,
    3. Coe A,
    4. Chisholm SW
    . 2011. The spontaneous mutation frequencies of Prochlorococcus strains are commensurate with those of other bacteria. Environ Microbiol Rep 3:744–749. doi:10.1111/j.1758-2229.2011.00293.x.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Paerl HW
    . 1996. A comparison of cyanobacterial bloom dynamics in freshwater, estuarine and marine environments. Phycologia 35:25–35. doi:10.2216/i0031-8884-35-6S-25.1.
    OpenUrlCrossRef
  21. 21.↵
    1. Paerl HW,
    2. Fulton RS
    . 2006. Ecology of harmful cyanobacteria, p 95–109. In Granéli E, Turner JT (ed), Ecology of harmful algae. Springer Berlin Heidelberg, Berlin, Heidelberg.
  22. 22.↵
    1. Mlouka A,
    2. Comte K,
    3. Castets A-M,
    4. Bouchier C,
    5. Tandeau de Marsac N
    . 2004. The gas vesicle gene cluster from Microcystis aeruginosa and DNA rearrangements that lead to loss of cell buoyancy. J Bacteriol 186:2355–2365. doi:10.1128/jb.186.8.2355-2365.2004.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Carey CC,
    2. Ibelings BW,
    3. Hoffmann EP,
    4. Hamilton DP,
    5. Brookes JD
    . 2012. Eco-physiological adaptations that favour freshwater cyanobacteria in a changing climate. Water Res 46:1394–1407. doi:10.1016/j.watres.2011.12.016.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Stanier RY
    . 1973. Autotrophy and heterotrophy in unicellular blue-green algae, p 501–508. In Carr NG, Whitton BA (ed), The biology of blue-green algae. University of California Press, Berkeley, CA.
  25. 25.↵
    1. Liu Z,
    2. Häder DP,
    3. Sommaruga R
    . 2004. Occurrence of mycosporine-like amino acids (MAAs) in the bloom-forming cyanobacterium Microcystis aeruginosa. J Plankton Res 26:963–966. doi:10.1093/plankt/fbh083.
    OpenUrlCrossRefWeb of Science
  26. 26.↵
    1. Jia Y,
    2. Chen W,
    3. Zuo Y,
    4. Lin L,
    5. Song L
    . 2018. Heavy metal migration and risk transference associated with cyanobacterial blooms in eutrophic freshwater. Sci Total Environ 613-614:1324–1330. doi:10.1016/j.scitotenv.2017.09.180.
    OpenUrlCrossRef
  27. 27.↵
    1. Hoostal MJ,
    2. Bidart-Bouzat MG,
    3. Bouzat JL
    . 2008. Local adaptation of microbial communities to heavy metal stress in polluted sediments of Lake Erie. FEMS Microbiol Ecol 65:156–168. doi:10.1111/j.1574-6941.2008.00522.x.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Guo Y,
    2. Liu M,
    3. Liu L,
    4. Liu X,
    5. Chen H,
    6. Yang J
    . 2018. The antibiotic resistome of free-living and particle-attached bacteria under a reservoir cyanobacterial bloom. Environ Int 117:107–115. doi:10.1016/j.envint.2018.04.045.
    OpenUrlCrossRef
  29. 29.↵
    1. Hudnell HK,
    2. Dortch Q
    . 2008. A synopsis of research needs identified at the interagency, International Symposium on Cyanobacterial Harmful Algal Blooms (ISOC-HAB), p 17–43. In Hudnell HK (ed), Cyanobacterial harmful algal blooms: state of the science and research needs. Springer New York, New York, NY.
  30. 30.↵
    1. Shi T,
    2. Falkowski PG
    . 2008. Genome evolution in cyanobacteria: the stable core and the variable shell. Proc Natl Acad Sci U S A 105:2510–2515. doi:10.1073/pnas.0711165105.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Walsby AE,
    2. Hayes PK,
    3. Boje R,
    4. Stal LJ
    . 1997. The selective advantage of buoyancy provided by gas vesicles for planktonic cyanobacteria in the Baltic Sea. New Phytol 136:407–417. doi:10.1046/j.1469-8137.1997.00754.x.
    OpenUrlCrossRefWeb of Science
  32. 32.↵
    1. Holland A,
    2. Kinnear S
    . 2013. Interpreting the possible ecological role(s) of cyanotoxins: compounds for competitive advantage and/or physiological aide? Mar Drugs 11:2239–2258. doi:10.3390/md11072239.
    OpenUrlCrossRef
  33. 33.↵
    1. Gan N,
    2. Xiao Y,
    3. Zhu L,
    4. Wu Z,
    5. Liu J,
    6. Hu C,
    7. Song L
    . 2012. The role of microcystins in maintaining colonies of bloom-forming Microcystis spp. Environ Microbiol 14:730–742. doi:10.1111/j.1462-2920.2011.02624.x.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Zilliges Y,
    2. Kehr J-C,
    3. Mikkat S,
    4. Bouchier C,
    5. de Marsac NT,
    6. Börner T,
    7. Dittmann E
    . 2008. An extracellular glycoprotein is implicated in cell-cell contacts in the toxic cyanobacterium Microcystis aeruginosa pcc 7806. J Bacteriol 190:2871–2879. doi:10.1128/JB.01867-07.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Dai R,
    2. Liu H,
    3. Qu J,
    4. Zhao X,
    5. Hou Y
    . 2009. Effects of amino acids on microcystin production of the Microcystis aeruginosa. J Hazard Mater 161:730–736. doi:10.1016/j.jhazmat.2008.04.015.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Ylla I,
    2. Sanpera-Calbet I,
    3. Muñoz I,
    4. Romaní AM,
    5. Sabater S
    . 2011. Organic matter characteristics in a Mediterranean stream through amino acid composition: changes driven by intermittency. Aquat Sci 73:523–535. doi:10.1007/s00027-011-0211-x.
    OpenUrlCrossRef
  37. 37.↵
    1. Thomas JD
    . 1997. The role of dissolved organic matter, particularly free amino acids and humic substances, in freshwater ecosystems. Freshw Biol 38:1–36. doi:10.1046/j.1365-2427.1997.00206.x.
    OpenUrlCrossRef
  38. 38.↵
    1. Rosenstock B,
    2. Simon M
    . 2001. Sources and sinks of dissolved free amino acids and protein in a large and deep mesotrophic lake. Limnol Oceanogr 46:644–654. doi:10.4319/lo.2001.46.3.0644.
    OpenUrlCrossRef
  39. 39.↵
    1. Burkholder JM,
    2. Glibert PM,
    3. Skelton HM
    . 2008. Mixotrophy, a major mode of nutrition for harmful algal species in eutrophic waters. Harmful Algae 8:77–93. doi:10.1016/j.hal.2008.08.010.
    OpenUrlCrossRefWeb of Science
  40. 40.↵
    1. Baran R,
    2. Ivanova NN,
    3. Jose N,
    4. Garcia-Pichel F,
    5. Kyrpides NC,
    6. Gugger M,
    7. Northen TR
    . 2013. Functional genomics of novel secondary metabolites from diverse cyanobacteria using untargeted metabolomics. Mar Drugs 11:3617–3631. doi:10.3390/md11103617.
    OpenUrlCrossRef
  41. 41.↵
    1. Sakamoto T,
    2. Bryant DA,
    3. Wada H,
    4. Murata N,
    5. Higashi S
    . 1997. Low-temperature-induced desaturation of fatty acids and expression of desaturase genes in the cyanobacterium Synechococcus sp. PCC 7002. FEMS Microbiol Lett 152:313–320. doi:10.1111/j.1574-6968.1997.tb10445.x.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Gombos Z,
    2. Wada H,
    3. Murata N
    . 1992. Unsaturation of fatty acids in membrane lipids enhances tolerance of the cyanobacterium Synechocystis PCC6803 to low-temperature photoinhibition. Proc Natl Acad Sci U S A 89:9959–9963. doi:10.1073/pnas.89.20.9959.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Gombos Z,
    2. Wada H,
    3. Hideg E,
    4. Murata N
    . 1994. The unsaturation of membrane lipids stabilizes photosynthesis against heat stress. Plant Physiol 104:563–567. doi:10.1104/pp.104.2.563.
    OpenUrlAbstract
  44. 44.↵
    1. Allakhverdiev SI,
    2. Kinoshita M,
    3. Inaba M,
    4. Suzuki I,
    5. Murata N
    . 2001. Unsaturated fatty acids in membrane lipids protect the photosynthetic machinery against salt-induced damage in Synechococcus. Plant Physiol 125:1842–1853. doi:10.1104/pp.125.4.1842.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Verspagen JMH,
    2. Van de Waal DB,
    3. Finke JF,
    4. Visser PM,
    5. Van Donk E,
    6. Huisman J
    . 2014. Rising CO2 levels will intensify phytoplankton blooms in eutrophic and hypertrophic lakes. PLoS One 9:e104325. doi:10.1371/journal.pone.0104325.
    OpenUrlCrossRef
  46. 46.↵
    1. Sandrini G,
    2. Tann RP,
    3. Schuurmans JM,
    4. van Beusekom SAM,
    5. Matthijs HCP,
    6. Huisman J
    . 2016. Diel variation in gene expression of the CO2-concentrating mechanism during a harmful cyanobacterial bloom. Front Microbiol 7:551. doi:10.3389/fmicb.2016.00551.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Beaulieu M,
    2. Pick F,
    3. Gregory-Eaves I
    . 2013. Nutrients and water temperature are significant predictors of cyanobacterial biomass in a 1147 lakes data set. Limnol Oceanogr 58:1736–1746. doi:10.4319/lo.2013.58.5.1736.
    OpenUrlCrossRef
  48. 48.↵
    1. Klausmeier CA,
    2. Litchman E,
    3. Daufresne T,
    4. Levin SA
    . 2004. Optimal nitrogen-to-phosphorus stoichiometry of phytoplankton. Nature 429:171–174. doi:10.1038/nature02454.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Harke MJ,
    2. Davis TW,
    3. Watson SB,
    4. Gobler CJ
    . 2016. Nutrient-controlled niche differentiation of western Lake Erie cyanobacterial populations revealed via metatranscriptomic surveys. Environ Sci Technol 50:604–615. doi:10.1021/acs.est.5b03931.
    OpenUrlCrossRef
  50. 50.↵
    1. Gao Q,
    2. Garcia-Pichel F
    . 2011. Microbial ultraviolet sunscreens. Nat Rev Microbiol 9:791–802. doi:10.1038/nrmicro2649.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Gao Q,
    2. Garcia-Pichel F
    . 2011. An ATP-grasp ligase involved in the last biosynthetic step of the iminomycosporine shinorine in Nostoc punctiforme atcc 29133. J Bacteriol 193:5923–5928. doi:10.1128/JB.05730-11.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Alneberg J,
    2. Bennke C,
    3. Beier S,
    4. Bunse C,
    5. Quince C,
    6. Ininbergs K,
    7. Riemann L,
    8. Ekman M,
    9. Jürgens K,
    10. Labrenz M,
    11. Pinhassi J,
    12. Andersson AF
    . 2020. Ecosystem-wide metagenomic binning enables prediction of ecological niches from genomes. Commun Biol 3:119. doi:10.1038/s42003-020-0856-x.
    OpenUrlCrossRef
  53. 53.↵
    1. Stockner J,
    2. Callieri C,
    3. Cronberg G
    . 2002. Picoplankton and other non-bloom-forming cyanobacteria in lakes, p 195–231. In Whitton BA, Potts M (ed), The ecology of cyanobacteria: their diversity in time and space. Springer Netherlands, Dordrecht, The Netherlands.
  54. 54.↵
    1. Rippka R,
    2. Deruelles J,
    3. Waterbury JB,
    4. Herdman M,
    5. Stanier RY
    . 1979. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology 111:1–61. doi:10.1099/00221287-111-1-1.
    OpenUrlCrossRefWeb of Science
  55. 55.↵
    1. Overbeek R,
    2. Begley T,
    3. Butler RM,
    4. Choudhuri JV,
    5. Chuang H-Y,
    6. Cohoon M,
    7. de Crécy-Lagard V,
    8. Diaz N,
    9. Disz T,
    10. Edwards R,
    11. Fonstein M,
    12. Frank ED,
    13. Gerdes S,
    14. Glass EM,
    15. Goesmann A,
    16. Hanson A,
    17. Iwata-Reuyl D,
    18. Jensen R,
    19. Jamshidi N,
    20. Krause L,
    21. Kubal M,
    22. Larsen N,
    23. Linke B,
    24. McHardy AC,
    25. Meyer F,
    26. Neuweger H,
    27. Olsen G,
    28. Olson R,
    29. Osterman A,
    30. Portnoy V,
    31. Pusch GD,
    32. Rodionov DA,
    33. Rückert C,
    34. Steiner J,
    35. Stevens R,
    36. Thiele I,
    37. Vassieva O,
    38. Ye Y,
    39. Zagnitko O,
    40. Vonstein V
    . 2005. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res 33:5691–5702. doi:10.1093/nar/gki866.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Kanehisa M,
    2. Goto S,
    3. Sato Y,
    4. Kawashima M,
    5. Furumichi M,
    6. Tanabe M
    . 2014. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res 42:D199–D205. doi:10.1093/nar/gkt1076.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    1. Saier MH,
    2. Reddy VS,
    3. Tamang DG,
    4. Västermark Å
    . 2014. The transporter classification database. Nucleic Acids Res 42:D251–D258. doi:10.1093/nar/gkt1097.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Camacho C,
    2. Coulouris G,
    3. Avagyan V,
    4. Ma N,
    5. Papadopoulos J,
    6. Bealer K,
    7. Madden TL
    . 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi:10.1186/1471-2105-10-421.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Cao H,
    2. Shimura Y,
    3. Masanobu K,
    4. Yin Y
    . 2014. Draft genome sequence of the toxic bloom-forming cyanobacterium Aphanizomenon flos-aquae NIES-81. Genome Announc 2:e00044-14. doi:10.1128/genomeA.00044-14.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Download PDF
Citation Tools
The Trait Repertoire Enabling Cyanobacteria to Bloom Assessed through Comparative Genomic Complexity and Metatranscriptomics
Huansheng Cao, Yohei Shimura, Morgan M. Steffen, Zhou Yang, Jingrang Lu, Allen Joel, Landon Jenkins, Masanobu Kawachi, Yanbin Yin, Ferran Garcia-Pichel
mBio Jun 2020, 11 (3) e01155-20; DOI: 10.1128/mBio.01155-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
The Trait Repertoire Enabling Cyanobacteria to Bloom Assessed through Comparative Genomic Complexity and Metatranscriptomics
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
The Trait Repertoire Enabling Cyanobacteria to Bloom Assessed through Comparative Genomic Complexity and Metatranscriptomics
Huansheng Cao, Yohei Shimura, Morgan M. Steffen, Zhou Yang, Jingrang Lu, Allen Joel, Landon Jenkins, Masanobu Kawachi, Yanbin Yin, Ferran Garcia-Pichel
mBio Jun 2020, 11 (3) e01155-20; DOI: 10.1128/mBio.01155-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

cyanobacterial bloom
comparative genomics
ecophysiology
genomic complexity
metatranscriptome
Microcystis aeruginosa
adaptation
cyanobacteria
ecogenomics
water blooms

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511