Genome sequences of lower Great Lakes Microcystis sp. reveal strain-specific genes that are present and expressed in western Lake Erie blooms

PLoS One. 2017 Oct 11;12(10):e0183859. doi: 10.1371/journal.pone.0183859. eCollection 2017.

Abstract

Blooms of the potentially toxic cyanobacterium Microcystis are increasing worldwide. In the Laurentian Great Lakes they pose major socioeconomic, ecological, and human health threats, particularly in western Lake Erie. However, the interpretation of "omics" data is constrained by the highly variable genome of Microcystis and the small number of reference genome sequences from strains isolated from the Great Lakes. To address this, we sequenced two Microcystis isolates from Lake Erie (Microcystis aeruginosa LE3 and M. wesenbergii LE013-01) and one from upstream Lake St. Clair (M. cf aeruginosa LSC13-02), and compared these data to the genomes of seventeen Microcystis spp. from across the globe as well as one metagenome and seven metatranscriptomes from a 2014 Lake Erie Microcystis bloom. For the publically available strains analyzed, the core genome is ~1900 genes, representing ~11% of total genes in the pan-genome and ~45% of each strain's genome. The flexible genome content was related to Microcystis subclades defined by phylogenetic analysis of both housekeeping genes and total core genes. To our knowledge this is the first evidence that the flexible genome is linked to the core genome of the Microcystis species complex. The majority of strain-specific genes were present and expressed in bloom communities in Lake Erie. Roughly 8% of these genes from the lower Great Lakes are involved in genome plasticity (rapid gain, loss, or rearrangement of genes) and resistance to foreign genetic elements (such as CRISPR-Cas systems). Intriguingly, strain-specific genes from Microcystis cultured from around the world were also present and expressed in the Lake Erie blooms, suggesting that the Microcystis pangenome is truly global. The presence and expression of flexible genes, including strain-specific genes, suggests that strain-level genomic diversity may be important in maintaining Microcystis abundance during bloom events.

MeSH terms

  • Base Sequence
  • CRISPR-Cas Systems / genetics
  • Eutrophication*
  • Gene Expression Regulation, Bacterial*
  • Genes, Bacterial*
  • Great Lakes Region
  • Metagenome
  • Microcystis / genetics*
  • Microcystis / isolation & purification
  • Phylogeny
  • Species Specificity

Grants and funding

This work was supported by a grant from the Erb Family Foundation made through the University of Michigan Water Center (Grant N017871) and by an Environmental Protection Agency Great Lakes Restoration Initiative Grant (2015-062a). Funding was awarded to the Cooperative Institute for Limnology and Ecosystems Research through the National Oceanic and Atmospheric Administration Cooperative Agreement with the University of Michigan (NA12OAR4320071). This project was supported by grants from the University of Michigan Office for Research MCubed program and the Erb Family Foundation made through the University of Michigan Water Center. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.