Bacteria play key roles in the function and diversity of aquatic systems, but aside from study of specific bloom systems, little is known about the diversity or biogeography of bacteria associated with harmful cyanobacterial blooms (cyanoHABs). CyanoHAB species are known to shape bacterial community composition and to rely on functions provided by the associated bacteria, leading to the hypothesized cyanoHAB interactome, a coevolved community of synergistic and interacting bacteria species, each necessary for the success of the others. Here, we surveyed the microbiome associated with Microcystis aeruginosa during blooms in 12 lakes spanning four continents as an initial test of the hypothesized Microcystis interactome. We predicted that microbiome composition and functional potential would be similar across blooms globally. Our results, as revealed by 16S rRNA sequence similarity, indicate that M. aeruginosa is cosmopolitan in lakes across a 280 longitudinal and 90 latitudinal gradient. The microbiome communities were represented by a wide range of operational taxonomic units and relative abundances. Highly abundant taxa were more related and shared across most sites and did not vary with geographic distance, thus, like Microcystis, revealing no evidence for dispersal limitation. High phylogenetic relatedness, both within and across lakes, indicates that microbiome bacteria with similar functional potential were associated with all blooms. While Microcystis and the microbiome bacteria shared many genes, whole-community metagenomic analysis revealed a suite of biochemical pathways that could be considered complementary. Our results demonstrate a high degree of similarity across global Microcystis blooms, thereby providing initial support for the hypothesized Microcystis interactome.

No abstract

SummaryMany podoviruses have been isolated which infect marine picocyanobacteria, and they may play a potentially important role in regulating the biomass and population composition of picocyanobacteria. However, little is known about the diversity and population dynamics of autochthonous cyanopodoviruses in marine environments. Using a set of newly designed PCR primers which specifically amplify the DNA pol from cyanopodoviruses, a total of 221 DNA pol sequences were retrieved from eight Chesapeake Bay virioplankton communities collected at different times and locations. All DNA pol sequences clustered with the eight known podoviruses that infect different marine picocyanobacteria, and could be divided into at least 10 different subclusters (I-X). The presence of these cyanopodovirus genotypes based on PCRamplification of DNA pol gene sequences was supported by the existence of similar DNA pol genotypes with metagenome libraries of Chesapeake Bay virioplankton assemblages. The composition of cyanopodoviruses in the Bay also exhibited distinct winter and summer patterns which were likely related to corresponding seasonal changes in the composition of cyanobacterial populations. Our study suggests that diverse and dynamic populations of cyanopodoviruses are present in the estuarine environment. The PCR method developed in this study provides a specific and sensitive tool to explore the abundance, distribution and phylogenetic diversity of cyanopodoviruses in aquatic environments. Linking the dynamics of host and viral populations in the natural environment is critical to broader characterization of the ecological role of virioplankton within microbial communities.