Multilocus sequence typing (MLST) reveals high genetic diversity and clonal population structure of the toxic cyanobacterium Microcystis aeruginosa

The occurrence of positive selection on BicA transporter of Microcystis aeruginosa

The rapid growth of cyanobacteria, particularly Microcystis aeruginosa, poses a significant threat to global water security. The proliferation of toxic Microcystis aeruginosa raises concerns due to its potential harm to human health and socioeconomic impacts. Dense blooms contribute to spatiotemporal inorganic carbon depletion, promoting interest in the roles of carbon-concentrating mechanisms (CCMs) for competitive carbon uptake. Despite the importance of HCO3- transporters, genetic evaluations and functional predictions in M. aeruginosa remain insufficient. In this study, we explored the diversity of HCO3- transporters in the genomes of 46 strains of M. aeruginosa, assessing positive selection for each. Intriguingly, although the Microcystis BicA transporter became a partial gene in 23 out of 46 genomic strains, we observed significant positive sites. Structural analyses, including predicted 2D and 3D models, confirmed the structural conservation of the Microcystis BicA transporter. Our findings suggest that the Microcystis BicA transport likely plays a crucial role in competitive carbon uptake, emphasizing its ecological significance. The ecological function of the Microcystis BicA transport in competitive growth during cyanobacterial blooms raises important questions. Future studies require experimental confirmation to better understand the role of the Microcysits BicA transporter in cyanobacterial blooms dynamics.

Trade-off in genome turnover events leading to adaptive evolution of Microcystis aeruginosa species complex

BACKGROUND

Numerous studies in the past have expanded our understanding of the genetic differences of global distributed cyanobacteria that originated around billions of years ago, however, unraveling how gene gain and loss drive the genetic evolution of cyanobacterial species, and the trade-off of these evolutionary forces are still the central but poorly understood issues.

RESULTS

To delineate the contribution of gene flow in mediating the hereditary differentiation and shaping the microbial evolution, a global genome-wide study of bloom-forming cyanobacterium, Microcystis aeruginosa species complex, provided robust evidence for genetic diversity, reflected by enormous variation in gene repertoire among various strains. Mathematical extrapolation showed an 'open' microbial pan-genome of M. aeruginosa species, since novel genes were predicted to be introduced after new genomes were sequenced. Identification of numerous horizontal gene transfer's signatures in genome regions of interest suggested that genome expansion via transformation and phage-mediated transduction across bacterial lineage as an evolutionary route may contribute to the differentiation of Microcystis functions (e.g., carbohydrate metabolism, amino acid metabolism, and energy metabolism). Meanwhile, the selective loss of some dispensable genes at the cost of metabolic versatility is as a mean of adaptive evolution that has the potential to increase the biological fitness.

CONCLUSIONS

Now that the recruitment of novel genes was accompanied by a parallel loss of some other ones, a trade-off in gene content may drive the divergent differentiation of M. aeruginosa genomes. Our study provides a genetic framework for the evolution of M. aeruginosa species and illustrates their possible evolutionary patterns.

Microcystis pangenome reveals cryptic diversity within and across morphospecies

Microcystis, a common harmful algal bloom (HAB) taxon, threatens water supplies and human health, yet species delimitation is contentious in this taxon, leading to challenges in research and management of this threat. Historical and common morphology-based classifications recognize multiple morphospecies, most with variable and diverse ecologies, while DNA sequence-based classifications indicate a single species with multiple ecotypes. To better delimit Microcystis species, we conducted a pangenome analysis of 122 genomes. Core- and non-core gene phylogenetic analyses placed 113 genomes into 23 monophyletic clusters containing at least two genomes. Overall, genome-related indices revealed that Microcystis contains at least 16 putative genospecies. Fifteen genospecies included at least one Microcystis aeruginosa morphospecies, and 10 genospecies included two or more morphospecies. This classification system will enable consistent taxonomic identification of Microcystis and thereby aid in resolving some of the complexities and controversies that have long characterized eco-evolutionary research and management of this important HAB taxon.

Statistical Inferences Applying Non-Parametric Data on Cyanobacterial Investigations: Contributions to Water Quality and New Trends under Global Changes on Portuguese Freshwater Ecosystems

Cyanobacteria are a bloom-forming ancient group of photosynthetic prokaryotes. A rise in temperature is a major contributor to its massive proliferation, namely on freshwater ecosystems, with social and economic impacts. Thus, reliable and cost-effective tools can permit the fast surveillance and assessment of temperature effects on potentially toxic cyanobacteria distribution and impacts. The occurrence of three potentially toxic cyanobacteria species was assessed on seven sampling points across three sampling years. Moreover, the association between the occurrence of those cyanobacteria species with climate change events was addressed. Here, we combined molecular and statistical methods to study the impacts of temperature on the occurrence of three globally occurring cyanotoxin-producing cyanobacteria species—Microcystis aeruginosa (microcystins), Raphidiopsis raciborskii (cylindrospermopsins and saxitoxins) and Planktothrix agardhii (microcystins and saxitoxins). Samples were collected on seven European temperate freshwater systems located on the North and Centre regions of Portugal, across three distinct sampling years with distinct ranges of air temperature. Data support that M. aeruginosa is still a common inhabitant of Portuguese freshwater ecosystems and a new trend was found on R. raciborskii recent invasion and establishment on the colder north ecosystems of Portugal. Additionally, the highest frequency of detection of both cyanobacteria was associated with warmer years. P. agardhii also revealed a new trend, being reported for the first time on North and Centre Regions of Portugal, however with no statistical relation with air temperature, demonstrating a higher ecological fitness. Distinct profiles of the statistical analysis on the three tested cyanobacteria species contribute to deepen the studies on other species as well as of our analyzed species on a global level. This assessment may help to anticipate possible repercussions on water quality and public health due to most probable alterations on cyanotoxins profile given the ecological fitness established among air temperature and PCR detection of potentially toxic cyanobacteria.

Comprehensive and functional analyses reveal the genomic diversity and potential toxicity of Microcystis

Microcystis is a cyanobacteria that is widely distributed across the world. It has attracted great attention because it produces the hepatotoxin microcystin (MC) that can inhibit eukaryotic protein phosphatases and pose a great risk to animal and human health. Due to the high diversity of morphospecies and genomes, it is still difficult to classify Microcystis species. In this study, we investigated the pangenome of 23 Microcystis strains to detect the genetic diversity and evolutionary dynamics. Microcystis revealed an open pangenome containing 22,009 gene families and exhibited different functional constraints. The core-genome phylogenetic analysis accurately differentiated the toxic and nontoxic strains and could be used as a taxonomic standard at the genetic level. We also investigated the functions of HGT events, of which were mostly conferred from cyanobacteria and closely related species. In order to detect the potential toxicity of Microcystis, we searched and characterized MC biosynthetic gene clusters and other secondary metabolite gene clusters. Our work provides insights into the genetic diversity, evolutionary dynamics, and potential toxicity of Microcystis, which could benefit the species classification and development of new methods for drinking water quality control and management of bloom formation in the future.

Genetic detection of freshwater harmful algal blooms: A review focused on the use of environmental DNA (eDNA) in Microcystis aeruginosa and Prymnesium parvum

Recurrence and severity of harmful algal blooms (HABs) are increasing due to a number of factors, including human practices and climate change. Sensitive and robust methods that allow for early and expedited HAB detection across large landscape scales are needed. Among the suite of HAB detection tools available, a powerful option exists in genetics-based approaches utilizing environmental sampling, also termed environmental DNA (eDNA). Here we provide a detailed methodological review of three HAB eDNA approaches (quantitative PCR, high throughput sequencing, and isothermal amplification). We then summarize and synthesize recently published eDNA applications covering a variety of HAB surveillance and research objectives, all with a specific emphasis in the detection of two widely problematic freshwater species, Microcystis aeruginosa and Prymnesium parvum. In our summary and conclusion we build on this literature by discussing ways in which eDNA methods could be advanced to improve HAB detection. We also discuss ways in which eDNA data could be used to potentially provide novel insight into the ecology, mitigation, and prediction of HABs.

Cyanobacterial blooms in China: diversity, distribution, and cyanotoxins

Cyanobacterial blooms, which refer to the massive growth of harmful cyanobacteria, have altered the global freshwater ecosystems during the past decades. China has the largest population in the world, and it is suffering from the harmful effect of water eutrophication and cyanobacterial blooms along with rapid development of the economy and society. Research on cyanobacterial blooms and cyanotoxins in China have been overwhelmingly enhanced and emphasized during the past decades. In the present review, the research on cyanobacterial blooms in China is generally introduced, including the history of cyanobacterial bloom studies, the diversity of the bloom-forming cyanobacteria species (BFCS), and cyanotoxin studies in China. Most studies have focused on Microcystis, its blooms, and microcystins. Newly emerging blooms with the dominance of non-Microcystis BFCS have been gradually expanding to wide regions in China. Understanding the basic features of these non-Microcystis BFCS and their blooms, including their diversity, occurrence, physio-ecology, and harmful metabolites, will provide direction on future studies of cyanobacterial blooms in China.

Microcoleus (Cyanobacteria) form watershed-wide populations without strong gradients in population structure

The relative importance of separation by distance and by environment to population genetic diversity can be conveniently tested in river networks, where these two drivers are often independently distributed over space. To evaluate the importance of dispersal and environmental conditions in shaping microbial population structures, we performed genome‐resolved metagenomic analyses of benthic Microcoleus‐dominated cyanobacterial mats collected in the Eel and Russian River networks (California, USA). The 64 Microcoleus genomes were clustered into three species that shared >96.5% average nucleotide identity (ANI). Most mats were dominated by one strain, but minor alleles within mats were often shared, even over large spatial distances (>300 km). Within the most common Microcoleus species, the ANI between the dominant strains within mats decreased with increasing spatial separation. However, over shorter spatial distances (tens of kilometres), mats from different subwatersheds had lower ANI than mats from the same subwatershed, suggesting that at shorter spatial distances environmental differences between subwatersheds in factors like canopy cover, conductivity, and mean annual temperature decreases ANI. Since mats in smaller creeks had similar levels of nucleotide diversity (π) as mats in larger downstream subwatersheds, within‐mat genetic diversity does not appear to depend on the downstream accumulation of upstream‐derived strains. The four‐gamete test and sequence length bias suggest recombination occurs between almost all strains within each species, even between populations separated by large distances or living in different habitats. Overall, our results show that, despite some isolation by distance and environmental conditions, sufficient gene‐flow occurs among cyanobacterial strains to prevent either driver from producing distinctive population structures across the watershed.

Cyanobacteria and Cyanotoxins in a Changing Environment: Concepts, Controversies, Challenges

Concern is widely being published that the occurrence of toxic cyanobacteria is increasing in consequence of climate change and eutrophication, substantially threatening human health. Here, we review evidence and pertinent publications to explore in which types of waterbodies climate change is likely to exacerbate cyanobacterial blooms; whether controlling blooms and toxin concentrations requires a balanced approach of reducing not only the concentrations of phosphorus (P) but also those of nitrogen (N); how trophic and climatic changes affect health risks caused by toxic cyanobacteria. We propose the following for further discussion: (i) Climate change is likely to promote blooms in some waterbodies—not in those with low concentrations of P or N stringently limiting biomass, and more so in shallow than in stratified waterbodies. Particularly in the latter, it can work both ways—rendering conditions for cyanobacterial proliferation more favourable or less favourable. (ii) While N emissions to the environment need to be reduced for a number of reasons, controlling blooms can definitely be successful by reducing only P, provided concentrations of P can be brought down to levels sufficiently low to stringently limit biomass. Not the N:P ratio, but the absolute concentration of the limiting nutrient determines the maximum possible biomass of phytoplankton and thus of cyanobacteria. The absolute concentrations of N or P show which of the two nutrients is currently limiting biomass. N can be the nutrient of choice to reduce if achieving sufficiently low concentrations has chances of success. (iii) Where trophic and climate change cause longer, stronger and more frequent blooms, they increase risks of exposure, and health risks depend on the amount by which concentrations exceed those of current WHO cyanotoxin guideline values for the respective exposure situation. Where trophic change reduces phytoplankton biomass in the epilimnion, thus increasing transparency, cyanobacterial species composition may shift to those that reside on benthic surfaces or in the metalimnion, changing risks of exposure. We conclude that studying how environmental changes affect the genotype composition of cyanobacterial populations is a relatively new and exciting research field, holding promises for understanding the biological function of the wide range of metabolites found in cyanobacteria, of which only a small fraction is toxic to humans. Overall, management needs case-by-case assessments focusing on the impacts of environmental change on the respective waterbody, rather than generalisations.

Preliminary evidence on the presence of cyanobacteria and cyanotoxins from culture enrichments followed by PCR analysis: new perspectives from Africa (Mali) and South Pacific (Fiji) countries

Cyanobacteria are a group of microorganisms that can be found in a diverse range of biogeographical areas and produce potent and damaging cyanotoxins, which reveal importance for continuous studies and surveillance efforts. In this study, we analyzed worldwide two-month culture-enriched water samples collected from 12 distinct countries (Costa Rica, Cuba, Fiji, France, Indonesia, Mali, Portugal, South Africa, Spain, Thailand, USA, Vietnam) including two undisclosed areas (Fiji and Mali). We performed a PCR-based molecular multi-step scheme that consisted in the detection of the main cyanobacterial species, genera, and cyanotoxins biosynthesis genes. Results from this study indicate that Microcystis aeruginosa followed by Planktothrix agardhii were the most prevalent species of all the 12 countries analyzed. Cylindrospermospis raciborskii was detected in Costa Rica, while P. agardhii was detected in Fiji and South Africa. M. aeruginosa was detected in Fiji and Mali. Regarding the main cyanotoxins biosynthesis genes, a cyrC gene fragment (cylindrospermopsins) was amplified in the African continent (South Africa), while anaC (anatoxin-a) was detected in two distinct locations, Mali and Vietnam. Saxitoxins biosynthesis gene was also detected in Fiji and Vietnam. Microcystins biosynthesis gene (mcyA) was co-detected with anatoxin-a biosynthesis gene in Mali and with saxitoxins biosynthesis gene (sxtI) in Portugal. This study therefore constitutes a major contribution to the global biogeography of cyanobacteria and its cyanotoxins and recommends continuous vigilance of toxic cyanobacteria particularly in the more undisclosed areas of the world. The PCR analysis data obtained in our 2-month culture-enriched water samples supports molecular methods as a preliminary tool in the environmental surveillance of cyanobacteria and cyanotoxins in undisclosed locations, particularly since the several positive amplifications detected may indicate that though samples were collected under non-bloom conditions, if environmental conditions change in the ecosystem, there is a risk that bloom-forming species may arose along with their detected cyanotoxicity.

Draft Genome Sequences of Four Microcystis aeruginosa Strains (NIES-3787, NIES-3804, NIES-3806, and NIES-3807) Isolated from Lake Kasumigaura, Japan

Microcystis aeruginosa is a bloom-forming cyanobacterium found in freshwater environments. The draft genomes of the M. aeruginosa strains NIES-3787, NIES-3804, NIES-3806, and NIES-3807, which were isolated from Lake Kasumigaura, Japan, were sequenced. The genome sizes of NIES-3787, NIES-3804, NIES-3806, and NIES-3807 were 4,524,637, 4,522,701, 4,370,004, and 4,378,226 bp, respectively.

Microcystis aeruginosa is a bloom-forming cyanobacterium found in freshwater environments. The draft genomes of the M. aeruginosa strains NIES-3787, NIES-3804, NIES-3806, and NIES-3807, which were isolated from Lake Kasumigaura, Japan, were sequenced. The genome sizes of NIES-3787, NIES-3804, NIES-3806, and NIES-3807 were 4,524,637, 4,522,701, 4,370,004, and 4,378,226 bp, respectively.

Genomic Characteristics of the Toxic Bloom-Forming Cyanobacterium Microcystis aeruginosa NIES-102

Microcystis aeruginosa, a bloom-forming cyanobacterium distributed mainly in freshwater environments, can be divided into at least 12 groups (A-K and X) based on multi-locus phylogenetic analyses. In this study, we characterized the genome of microcystin-producing M. aeruginosa NIES-102, assigned to group A, isolated from Lake Kasumigaura, Japan. The complete genome sequence of M. aeruginosa NIES-102 comprised a 5.87-Mbp circular chromosome containing 5,330 coding sequences. The genome was the largest among all sequenced genomes for the species. In a comparison with the genome of M. aeruginosa NIES-843, which belongs to the same group, the microcystin biosynthetic gene cluster and CRISPR-Cas locus were highly similar. A synteny analysis revealed small-scale rearrangements between the two genomes. Genes encoding transposases were more abundant in these two genomes than in other Microcystis genomes. Our results improve our understanding of structural genomic changes and adaptation to a changing environment in the species.

First Report of Microcystis Strains Producing MC-FR and -WR Toxins in Japan

Microcystins (MCs) are a group of cyclic heptapeptide hepatotoxins produced by Microcystis and several other genera of cyanobacteria. Many structural variants have been characterized using various methods such as liquid chromatography-mass spectrometry (LC-MS) analysis, enzyme-linked immunosorbent assay (ELISA) and protein phosphatase 2A (PP2A) inhibition assay. The representative MC, MC-LR, and related cyanobacterial toxins strongly inhibit PP2A activity and can therefore be assayed by measuring the extent of PP2A inhibition. However, these methods require reference toxin standards for the quantification and identification of known MCs. To obtain various MC-producing cyanobacterial strains, we surveyed and collected MC-producing cyanobacteria from environmental sources of water in Okinawa, Japan. Using a dual assay (LC-MS analysis and PP2A inhibition assay), we identified and isolated Microcystis strains producing five MC variants (MC-LR, -RR, -LA, -FR and -WR). Approximately 4 mg of MC-WR and -FR toxins were purified from the laboratory culture of the Microcystis isolate NIES-4344. Pure MC-WR and -FR variants were prepared for future use as toxin standards in LC-MS analysis. Phylogenetic analysis based on ftsZ revealed that the NIES-4344 strain belongs to the identified groups in Microcystis aeruginosa. This is the first report of Microcystis strains producing mainly MC-WR and -FR toxins in Japan.

A novel salt-tolerant genotype illuminates the sucrose gene evolution in freshwater bloom-forming cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is a water bloom-forming cyanobacterium found in fresh and brackish water ecosystems worldwide. Previously, we showed that several instances of M. aeruginosa bloom in brackish water can be explained by the proliferation of salt-tolerant M. aeruginosa strains harboring genes for a compatible solute sucrose. However, evolutionary history of sucrose genes in M. aeruginosa remains unclear because salt-tolerant strains have been poorly described. Here, we characterized a novel salt-tolerant strain of M. aeruginosa (NIES-4325) isolated from the brackish water of Lake Abashiri, Japan. A whole-genome analysis of M. aeruginosa NIES-4325 identified genes for sucrose synthesis (sppA, spsA and susA). Quantitative sucrose and gene expression analyses suggested that sucrose is implicated in acclimation to high salt in NIES-4325. Notably, the sucrose genes of M. aeruginosa are monophyletic, yet sucrose genes of NIES-4325 are highly divergent from those of other salt-tolerant M. aeruginosa strains. This suggests an early sucrose gene import into M. aeruginosa from other cyanobacteria, followed by multiple losses during intraspecific diversification. One of a few survivors of salt-tolerant strains is a likely donor of recent horizontal spreads of sucrose genes across M. aeruginosa lineages.

Coherence of Microcystis species revealed through population genomics

Microcystis is a genus of freshwater cyanobacteria, which causes harmful blooms in ecosystems worldwide. Some Microcystis strains produce harmful toxins such as microcystin, impacting drinking water quality. Microcystis colony morphology, rather than genetic similarity, is often used to classify Microcystis into morphospecies. Yet colony morphology is a plastic trait, which can change depending on environmental and laboratory culture conditions, and is thus an inadequate criterion for species delineation. Furthermore, Microcystis populations are thought to disperse globally and constitute a homogeneous gene pool. However, this assertion is based on relatively incomplete characterization of Microcystis genomic diversity. To better understand these issues, we performed a population genomic analysis of 33 newly sequenced genomes mainly from Canada and Brazil. We identified 17 Microcystis clusters of genomic similarity, five of which correspond to monophyletic clades containing at least three newly sequenced genomes. Four out of these five clades match to named morphospecies. Notably, M. aeruginosa is paraphyletic, distributed across 12 genomic clusters, suggesting it is not a coherent species. A few clades of closely related isolates are specific to a unique geographic location, suggesting biogeographic structure over relatively short evolutionary time scales. Higher homologous recombination rates within than between clades further suggest that monophyletic groups might adhere to a Biological Species-like concept, in which barriers to gene flow maintain species distinctness. However, certain genes-including some involved in microcystin and micropeptin biosynthesis-are recombined between monophyletic groups in the same geographic location, suggesting local adaptation.

Polyphasic approach using multilocus analyses supports the establishment of the new aerophytic cyanobacterial genus Pycnacronema (Coleofasciculaceae, Oscillatoriales)

A new Phormidium-like genus was found during an investigation of Oscillatoriales diversity in Brazil. Eight aerophytic populations from south and southeastern regions were isolated in monospecific cultures and submitted to polyphasic evaluation. The populations presented homogeneous morphology with straight trichomes, not attenuated, and apical cell with thickened cell wall. Phylogenetic analyses based on 16S rRNA gene sequences showed that these populations, plus the Brazilian strain Phomidium sp. B-Tom from GenBank, formed a highly supported and distinctive clade, which corresponds to the new genus Pycnacronema, comprising six new species: P. brasiliensis (type species), P. arboriculum, P. conicum, P. marmoreum, P. rubrum, and P. savannensis. These results were confirmed and supported by rpoC1 and rbcL genes evaluated independently and by the concatenated analysis of 16S rRNA, rpoC1 and rbcL genes (for all species but P. savannensis). Secondary structures of the D1-D1', box-B, and V3 regions of the internal transcribed spacer were informative at specific level, being conserved in P. brasiliensis and variable among the other strains, also confirming the phylogenetic analyses. The generic name and specific epithets of the new taxa are proposed under the provisions of the International Code of Nomenclature of algae, fungi, and plants.

Adaptation of the Freshwater Bloom-Forming Cyanobacterium Microcystis aeruginosa to Brackish Water Is Driven by Recent Horizontal Transfer of Sucrose Genes

Microcystis aeruginosa is a bloom-forming cyanobacterium found in eutrophic water bodies worldwide. M. aeruginosa blooms usually occur in freshwater; however, they have also been reported to occur in brackish water. Because M. aeruginosa often produces the cyanotoxin microcystin, they are a major concern to public health and environment. Despite this, the ecology, genomic basis, and evolutionary process underlying the M. aeruginosa bloom invasion from fresh to brackish water have been poorly investigated. Hence, in the present study, we have sequenced and characterized genomes of two newly discovered salt-tolerant M. aeruginosa strains obtained from Japanese brackish water lakes (Lakes Shinji and Tofutsu). Both genomes contain a set of genes for the synthesis of osmolyte sucrose (sppA, spsA, and susA), hitherto identified in only one strain (PCC 7806) of M. aeruginosa. Chemical and gene expression analyses confirmed sucrose accumulation induced by salt. A comprehensive genetic survey of >200 strains indicated that sucrose genes are extremely rare in M. aeruginosa. Most surprisingly, comparative genome analyses of the three strains indicated extremely low genetic diversity in the sucrose genes compared with other core genome genes, suggesting very recent acquisitions via horizontal transfer. Invasion of M. aeruginosa blooms into brackish water may be a recent event triggered by anthropogenic eutrophication of brackish water.

Complete Genome Sequence of Microcystis aeruginosa NIES-2481 and Common Genomic Features of Group G M. aeruginosa

Microcystis aeruginosa is a freshwater bloom-forming cyanobacterium that is distributed worldwide. M. aeruginosa can be divided into at least 8 phylogenetic groups (A-G and X) at the intraspecific level. Here, we report the complete genome sequence of M. aeruginosa NIES-2481, which was isolated from Lake Kasumigaura, Japan, and is assigned to group G. The complete genome sequence of M. aeruginosa NIES-2481 comprises a 4.29-Mbp circular chromosome and a 147,539-bp plasmid; the circular chromosome and the plasmid contain 4,332 and 167 protein-coding genes, respectively. Comparative analysis with the complete genome of M. aeruginosa NIES-2549, which belongs to the same group with NIES-2481, showed that the genome size is the smallest level in previously sequenced M. aeruginosa strains, and the genomes do not contain a microcystin biosynthetic gene cluster in common. Synteny analysis revealed only small-scale rearrangements between the two genomes.

Genomic divergence and cohesion in a species of pelagic freshwater bacteria

Background

In many prokaryotic genera a clustered phylogeny is observed, akin to the occurrence of species in sexually reproducing organisms. For some taxa, homologous recombination has been invoked as the underlying mechanism providing genomic cohesion among conspecific individuals. Whether this mechanism is applicable to prokaryotes in freshwaters with low habitat connectivity – i.e. elevated geographic barriers to gene flow – is unclear. To investigate further we studied genomic trends within the globally abundant PnecC cluster (genus Polynucleobacter, Betaproteobacteria) and analyzed homologous recombination within the affiliated species P. asymbioticus.

Results

Comparisons among 20 PnecC genomes revealed a clearly discontinuous distribution of nucleotide sequence similarities. Among the nine conspecific individuals (P. asymbioticus) all average nucleotide identity (ANI) values were greater than 97%, whereas all other comparisons exhibited ANI values lower than 85%. The reconstruction of recombination and mutation events for the P. asymbioticus core genomes yielded an r/m ratio of 7.4, which is clearly above estimated thresholds for recombination to act as a cohesive force. Hotspots of recombination were found to be located in the flanking regions of genomic islands. Even between geographically separated habitats a high flux of recombination was evident. While a biogeographic population structure was suggested from MLST data targeting rather conserved loci, such a structure was barely visible when whole genome data was considered. However, both MLST and whole genome data showed evidence of differentiation between two lineages of P. asymbioticus. The ratios of non-synonymous to synonymous substitution rates as well as growth rates in transplantation experiments suggested that this divergence was not selectively neutral.

Conclusions

The high extent of homologous recombination among P. asymbioticus bacteria can act as a cohesive force that effectively counteracts genetic divergence. At least on a regional scale, homologous recombination can act across geographically separated ecosystems and therefore plays an important role in the evolution and consistency of bacterial freshwater species. A species model akin to the biological species concept may be applicable for P. asymbioticus. Nonetheless, two genetically distinct lineages have emerged and further research may clarify if their divergence has been initiated by reinforced geographical barriers or has been evolving in sympatry.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4199-z) contains supplementary material, which is available to authorized users.

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

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.

Neutral Evolution and Dispersal Limitation Produce Biogeographic Patterns in Microcystis aeruginosa Populations of Lake Systems

Molecular observations reveal substantial biogeographic patterns of cyanobacteria within systems of connected lakes. An important question is the relative role of environmental selection and neutral processes in the biogeography of these systems. Here, we quantify the effect of genetic drift and dispersal limitation by simulating individual cyanobacteria cells using an agent-based model (ABM). In the model, cells grow (divide), die, and migrate between lakes. Each cell has a full genome that is subject to neutral mutation (i.e., the growth rate is independent of the genome). The model is verified by simulating simplified lake systems, for which theoretical solutions are available. Then, it is used to simulate the biogeography of the cyanobacterium Microcystis aeruginosa in a number of real systems, including the Great Lakes, Klamath River, Yahara River, and Chattahoochee River. Model output is analyzed using standard bioinformatics tools (BLAST, MAFFT). The emergent patterns of nucleotide divergence between lakes are dynamic, including gradual increases due to accumulation of mutations and abrupt changes due to population takeovers by migrant cells (coalescence events). The model predicted nucleotide divergence is heterogeneous within systems, and for weakly connected lakes, it can be substantial. For example, Lakes Superior and Michigan are predicted to have an average genomic nucleotide divergence of 8200 bp or 0.14%. The divergence between more strongly connected lakes is much lower. Our results provide a quantitative baseline for future biogeography studies. They show that dispersal limitation can be an important factor in microbe biogeography, which is contrary to the common belief, and could affect how a system responds to environmental change.

Genetic divergence among toxic and non-toxic cyanobacteria of the dry zone of Sri Lanka

Sri Lanka has rich cyanobacterial diversity, however, only few studies have been conducted to identify the potential toxin producers in water bodies used for human consumption. As the detection of cyanotoxin is vital in water quality management, a study was done by employing 16S rRNA gene to explore the genetic divergence, phylogenetic relationships and potential toxin producing cyanobacteria in reservoirs and well waters in the dry zone of Sri Lanka. Forty five, 16S rRNA gene sequences were assayed and phylogenetic tree was constructed. Among 45 isolates, 20 isolates were classified as unidentified cyanobacteria and considered as novel cyanobacterial genera. Of 25 identified isolates, seven isolates were identified up to species level. With 16S rRNA phylogeny, 20 unidentified cyanobacterial isolates were able to place on their taxonomic positions up to order level. Results revealed that water samples understudy had vast cyanobacterial diversity with potential microcystin (MC) and cylindrospermopsin (CYN) producers and eleven clusters clearly demonstrated five cyanobacterial orders with more than 90% similarity irrespective to their toxicity which showed the suitability of 16S rRNA gene for taxonomic differentiation. Sixteen isolates had the potential to produce MC and two isolates to produce CYN. Findings of the study confirm the rich cyanobacterial diversity and the divergence among the potential cyanotoxin producers in the dry zone water bodies of Sri Lanka.

Genetic characterization of Microcystis aeruginosa isolates from Portuguese freshwater systems

Cyanobacteria are microorganisms that pose a serious threat to the aquatic waterways through the production of dense blooms under eutrophic conditions and the release of toxic secondary metabolites-cyanotoxins. Within cyanobacteria, the colonial planktonic Microcystis aeruginosa is widely distributed in both fresh and brackish aquatic environments throughout the world being frequently observed in the Portuguese water systems. Apart from the well-established distribution of M. aeruginosa in Portugal, knowledge of its genetic diversity and population structure is unknown. Therefore, in this study twenty-seven strains were obtained from the North, Centre and South regions of Portugal and were subjected to extensive phylogenetic analyses using simultaneously four distinct genetic markers (16S rRNA, 16S-23S ITS, DNA gyrase subunit ß and cell division protein (ftsZ)) encompassing in total 2834 bp. With this work we characterized the phylogenetic relationship among the Portuguese strains, with the southern strains showing higher genetic structure relatively to the North and Centre strains. A total of fifteen genotypes were determined for M. aeruginosa in Portuguese water systems revealing a high genetic diversity. This is also the first study to report geographic variation on the population structure of the Portuguese M. aeruginosa.

Rapid Classification and Identification of Microcystis aeruginosa Strains Using MALDI-TOF MS and Polygenetic Analysis

Matrix-assisted laser desorption-ionization-time-of-flight mass spectrometry (MALDI-TOF MS) was used to establish a rapid, simple, and accurate method to differentiate among strains of Microcystis aeruginosa, one of the most prevalent types of bloom-forming cyanobacteria. M. aeruginosa NIES-843, for which a complete genome has been sequenced, was used to characterize ribosomal proteins as biomarkers and to optimize conditions for observing ribosomal proteins as major peaks in a given mass spectrum. Thirty-one of 52 ribosomal subunit proteins were detected and identified along the mass spectrum. Fifty-five strains of M. aeruginosa from different habitats were analyzed using MALDI-TOF MS; among these samples, different ribosomal protein types were observed. A polygenetic analysis was performed using an unweighted pair-group method with arithmetic means and different ribosomal protein types to classify the strains into five major clades. Two clades primarily contained toxic strains, and the other three clades contained exclusively non-toxic strains. This is the first study to differentiate cyanobacterial strains using MALDI-TOF MS.

Multilocus sequence typing of Lactobacillus casei isolates from naturally fermented foods in China and Mongolia

Lactobacillus casei is a lactic acid bacterium used in manufacturing of many fermented food products. To investigate the genetic diversity and population biology of this food-related bacterium, 224 Lb. casei isolates and 5 reference isolates were examined by multilocus sequence typing (MLST). Among them, 224 Lb. casei isolates were isolated from homemade fermented foods, including naturally fermented dairy products, acidic gruel, and Sichuan pickles from 38 different regions in China and Mongolia. The MLST scheme was developed based on the analysis of 10 selected housekeeping genes (carB, clpX, dnaA, groEL, murE, pyrG, pheS, recA, rpoC, and uvrC). All 229 isolates could be allocated to 171 unique sequence types, including 25 clonal complexes and 71 singletons. The high index of association value (1.3524) and standardized index of association value (0.1503) indicate the formation of an underlying clonal population by all the isolates. However, split-decomposition, relative frequency of occurrence of recombination and mutation, and relative effect of recombination and mutation in the diversification values confirm that recombination may have occurred, and were more frequent than mutation during the evolution of Lb. casei. Results from Structure analyses (version 2.3; http://pritch.bsd.uchicago.edu/structure.html) demonstrated that there were 5 lineages in the Lb. casei isolates, and the overall relatedness built by minimum spanning tree showed no clear relationship between the clonal complexes with either the isolation sources or sampling locations of the isolates. Our newly developed MLST scheme of Lb. casei was an easy and valuable tool that, together with the construction of an MLST database, will contribute to further detailed studies on the evolution and population genetics of Lb. casei from various niches.

Assessment of the genetic and phenotypic diversity among rhizogenic Agrobacterium biovar 1 strains infecting solanaceous and cucurbit crops

Rhizogenic Agrobacterium biovar 1 strains have been found to cause extensive root proliferation on hydroponically grown Cucurbitaceae and Solanaceae crops, resulting in substantial economic losses. As these agrobacteria live under similar ecological conditions, infecting a limited number of crops, it may be hypothesized that genetic and phenotypic variation among such strains is relatively low. In this study we assessed the phenotypic diversity as well as the phylogenetic and evolutionary relationships of several rhizogenic Agrobacterium biovar 1 strains from cucurbit and solanaceous crops. A collection of 41 isolates was subjected to a number of phenotypic assays and characterized by MLSA targeting four housekeeping genes (16S rRNA gene, recA, rpoB and trpE) and two loci from the root-inducing Ri-plasmid (part of rolB and virD2). Besides phenotypic variation, remarkable genotypic diversity was observed, especially for some chromosomal loci such as trpE. In contrast, genetic diversity was lower for the plasmid-borne loci, indicating that the studied chromosomal housekeeping genes and Ri-plasmid-borne loci might not exhibit the same evolutionary history. Furthermore, phylogenetic and network analyses and several recombination tests suggested that recombination could be contributing in some extent to the evolutionary dynamics of rhizogenic Agrobacterium populations. Finally, a genomospecies-level identification analysis revealed that at least four genomospecies may occur on cucurbit and tomato crops (G1, G3, G8 and G9). Together, this study gives a first glimpse at the genetic and phenotypic diversity within this economically important plant pathogenic bacterium.

A polyphasic approach leading to the revision of the genus Planktothrix (Cyanobacteria) and its type species, P. agardhii, and proposal for integrating the emended valid botanical taxa, as well as three new species, Planktothrix paucivesiculata sp. nov.ICNP, Planktothrix tepida sp. nov.ICNP, and Planktothrix serta sp. nov.ICNP, as genus and species names with nomenclatural standing under the ICNP

Twenty strains of Planktothrix and five of 'Oscillatoria' were characterized by a polyphasic approach, for clarification of their taxonomic relationships. Emphasis was given to the strains (17) of the Pasteur Culture Collection of Cyanobacteria (PCC). Phenotypic characters analyzed comprised morphology, phycobiliprotein composition, temperature and salinity tolerance. The gvpA gas vesicle gene was detected by PCR in all strains, and transmission electron microscopy confirmed gas vesicle formation in the strains of 'Oscillatoria'. MALDI-TOF mass spectrometry revealed 13 chemotypes, nine of which produce microcystins. A multi-locus sequence typing (MLST) analysis was conducted using individual and concatenated nucleotide sequences of the 16S rDNA, internal transcribed spacer (ITS), gyrB, rpoC1 and rpoB. The results highlighted an unexpected diversity within the genus Planktothrix, showing that the five strains of 'Oscillatoria' need to be included in this taxon. Consequently, the genus consists of seven phylogenetic clusters, three of which represent new species, named Planktothrix paucivesiculata sp. nov.ICNP (type strain: PCC 8926T), Planktothrix tepida sp. nov.ICNP (type strain: PCC 9214T) and Planktothrix serta sp. nov.ICNP (type strain: PCC 8927T). These, together with the emended genus Planktothrix and its type species P. agardhii, valid taxa under the ICN, are described/re-described for gaining nomenclatural standing under the ICNP.

Secondary metabolite gene expression and interplay of bacterial functions in a tropical freshwater cyanobacterial bloom

Cyanobacterial harmful algal blooms (cyanoHABs) appear to be increasing in frequency on a global scale. The Cyanobacteria in blooms can produce toxic secondary metabolites that make freshwater dangerous for drinking and recreation. To characterize microbial activities in a cyanoHAB, transcripts from a eutrophic freshwater reservoir in Singapore were sequenced for six samples collected over one day-night period. Transcripts from the Cyanobacterium Microcystis dominated all samples and were accompanied by at least 533 genera primarily from the Cyanobacteria, Proteobacteria, Bacteroidetes and Actinobacteria. Within the Microcystis population, abundant transcripts were from genes for buoyancy, photosynthesis and synthesis of the toxin microviridin, suggesting that these are necessary for competitive dominance in the Reservoir. During the day, Microcystis transcripts were enriched in photosynthesis and energy metabolism while at night enriched pathways included DNA replication and repair and toxin biosynthesis. Microcystis was the dominant source of transcripts from polyketide and non-ribosomal peptide synthase (PKS and NRPS, respectively) gene clusters. Unexpectedly, expression of all PKS/NRPS gene clusters, including for the toxins microcystin and aeruginosin, occurred throughout the day-night cycle. The most highly expressed PKS/NRPS gene cluster from Microcystis is not associated with any known product. The four most abundant phyla in the reservoir were enriched in different functions, including photosynthesis (Cyanobacteria), breakdown of complex organic molecules (Proteobacteria), glycan metabolism (Bacteroidetes) and breakdown of plant carbohydrates, such as cellobiose (Actinobacteria). These results provide the first estimate of secondary metabolite gene expression, functional partitioning and functional interplay in a freshwater cyanoHAB.

African origin and europe-mediated global dispersal of the cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is a bloom-forming cyanobacteria, which currently has a cosmopolitan distribution. Since M. aeruginosa can produce toxic compounds across all continents that it inhabits, it is of major public health relevance to assess its origin and dispersal. Thus, we conducted a worldwide study using 29 isolates representative of all the main continents, and used a concatenated genetic system for phylogenetic analyses consisting of four genetic markers (spanning ca. 3,485 bp). Our results support an early origin of M. aeruginosa in the African continent, with a subsequent dispersal to establish a second genetic pool in the European continent, from where M. aeruginosa then colonized the remaining continental regions. Our findings indicate that the European population has a cosmopolitan distribution, and is genetically closer to populations from Africa and North America. Our study also highlights the utility of using a concatenated dataset for phylogenetic inferences in cyanobacteria.

Frequent recombination shapes the epidemic population structure of Planktothrix (Cyanoprokaryota) in Italian subalpine lakes

The planktonic genus Planktothrix, as other cyanobacteria, shows signals of both homologous and nonhomologous recombination. However, the frequency of recombination and its effect on Planktothrix population structuring is unknown. We isolated 290 Planktothrix strains from seven neighboring lakes in the subalpine Italian region and analyzed these using multilocus sequence typing. Four of six loci analyzed were polymorphic, resulting in 20 distinct multilocus genotypes. Association indices among alleles at different loci were suggestive of an "epidemic population structure," resulting from an explosive (and temporary) dominance of one genotype against a panmictic background. ClonalFrame analyses supported this view by detecting: (i) three major clades affected by three distinct recombination events, (ii) a recombination rate about equal to the mutation rate, and (iii) the fact that recombination had an impact on introducing molecular diversity more than double the mutation rate. Furthermore, analysis of molecular variance over an annual cycle in three of seven lakes revealed that both local clonal expansion and recombination processes affected among-lake diversity. Our observations suggest that recombination affects microevolution of Planktothrix and that an epidemic structure can emerge in populations of this genus.

Phylogeny and biogeography of cyanobacteria and their produced toxins

Phylogeny is an evolutionary reconstruction of the past relationships of DNA or protein sequences and it can further be used as a tool to assess population structuring, genetic diversity and biogeographic patterns. In the microbial world, the concept that everything is everywhere is widely accepted. However, it is much debated whether microbes are easily dispersed globally or whether they, like many macro-organisms, have historical biogeographies. Biogeography can be defined as the science that documents the spatial and temporal distribution of a given taxa in the environment at local, regional and continental scales. Speciation, extinction and dispersal are proposed to explain the generation of biogeographic patterns. Cyanobacteria are a diverse group of microorganisms that inhabit a wide range of ecological niches and are well known for their toxic secondary metabolite production. Knowledge of the evolution and dispersal of these microorganisms is still limited, and further research to understand such topics is imperative. Here, we provide a compilation of the most relevant information regarding these issues to better understand the present state of the art as a platform for future studies, and we highlight examples of both phylogenetic and biogeographic studies in non-symbiotic cyanobacteria and cyanotoxins.

A tribute to disorder in the genome of the bloom-forming freshwater cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is one of the most common bloom-forming cyanobacteria in freshwater ecosystems worldwide. This species produces numerous secondary metabolites, including microcystins, which are harmful to human health. We sequenced the genomes of ten strains of M. aeruginosa in order to explore the genomic basis of their ability to occupy varied environments and proliferate. Our findings show that M. aeruginosa genomes are characterized by having a large open pangenome, and that each genome contains similar proportions of core and flexible genes. By comparing the GC content of each gene to the mean value of the whole genome, we estimated that in each genome, around 11% of the genes seem to result from recent horizontal gene transfer events. Moreover, several large gene clusters resulting from HGT (up to 19 kb) have been found, illustrating the ability of this species to integrate such large DNA molecules. It appeared also that all M. aeruginosa displays a large genomic plasticity, which is characterized by a high proportion of repeat sequences and by low synteny values between the strains. Finally, we identified 13 secondary metabolite gene clusters, including three new putative clusters. When comparing the genomes of Microcystis and Prochlorococcus, one of the dominant picocyanobacteria living in marine ecosystems, our findings show that they are characterized by having almost opposite evolutionary strategies, both of which have led to ecological success in their respective environments.

Genotypic composition and the relationship between genotypic composition and geographical proximity of the cyanobacterium Microcystis aeruginosa in western Japan

Microcystis aeruginosa is one of the bloom-forming harmful algae in freshwater ecosystems. We genetically characterized Microcystis populations during bloom-forming periods in various reservoirs, lakes, and ponds in Japan during 2009. Using phylogenetic analysis, we evaluated the relationship between current genotype expansions and geographic location within western Japan and intraspecific variation. Microcystis aeruginosa colonies were isolated at 15 sites and were analyzed by sequencing the 16S-23S internal transcribed spacer (ITS) region of the ribosomal operon, and the potential to produce toxins was assessed by PCR-based detection of the microcystin synthetase gene mcyG. In total, 171 colonies were separated into 41 genotypes. The highest genotypic composition was detected in the south basin of Lake Biwa and the lowest in Lagoon Iba. Cluster analysis indicated no obvious association between genotypic composition and geographic distance. Thus, clear genetic differentiation accompanied by geographic origins was not found in western Japan. The resulting neighbor-joining tree revealed 3 clusters, 2 of which contained strains that showed both nonamplification and amplification of the mcyG gene.

Identification of free-living amoebae and amoeba-associated bacteria from reservoirs and water treatment plants by molecular techniques

The occurrence of free-living amoebae (FLA) was investigated in 83 water samples from reservoirs and water treatment plants, with culture positive in 64 of them (77.1%). Polymerase chain reaction (PCR) of partial 18S rRNA gene and ITS region was performed in order to identify amoeba isolates, and the presence of Legionella pneumophila , Mycobacterium spp., Pseudomonas spp., and Microcystis aeruginosa was investigated in 43 isolates of amoebae by multiplex PCR. Of the isolated amoebae, 31 were Acanthamoeba spp., 21 were Hartmannella vermiformis, 13 were Naegleria spp., and one was Vanella spp. T2, T4, and T5 genotypes of Acanthamoeba have been identified, and T4 isolates were grouped into five subgenotypes and graphically represented with a Weblog application. Inside amoebae, L. pneumophila was detected in 13.9% (6/43) of the isolates, and Pseudomonas spp. and Mycobacterium spp. were detected in 32.6% (14/43) and 41.9% (18/43), respectively. No statistical correlation was demonstrated between FLA isolation and seasonality, but the presence of intracellular bacteria was associated with warm water temperatures, and also the intracellular presence of Mycobacterium spp. and Pseudomonas spp. were associated. These results highlight the importance of amoebae in natural waters as reservoirs of potential pathogens and its possible role in the spread of bacterial genera with interest in public and environmental health.

Patterns of nucleotide diversity of the ldpA circadian gene in closely related species of cyanobacteria from extreme cold deserts

In the circadian system of cyanobacteria, the ldpA gene is a component of the input to the clock. We comparatively analyzed nucleotide polymorphism of this gene in populations of two closely related species of cyanobacteria (denoted as Synechococcus species S1 and S2, respectively) from extreme cold deserts in Antarctica, the Canadian Arctic, and Tibet. Although both species manifested similarly high haplotype diversities (0.990 and 0.809, respectively), the nucleotide diversity differed significantly (0.0091 in S1 and 0.0037 in S2). The populations of species S2 were more differentiated (F(ST) = 0.2242) compared to those of species S1 (F(ST) between 0.0296 and 0.1188). An analysis of positive selection with several tests yielded highly significant values (P < 0.01) for both species. On the other hand, these results may be somewhat compromised by fluctuating population sizes of the species. The apparent selection pressure coupled with the pronounced demographic factors, such as population expansion, small effective population size, and genetic drift, may thus result in the observed significant interpopulation differentiation and subsequent speciation of cyanobacteria.

Recombination Signals In The rpoC1 Gene Indicate Gene-Flow Between Planktothrix (Cyanoprokaryota) Species

The delineation of species boundaries in the potentially harmful cyanobacterium Planktothrix Anagnostidis et Komárek 1988 is particularly tangled. Genetic recombination has been invoked to explain the occurrence of overlapping biological traits among recognized species. Although horizontal gene transfer is shown as a driver of diversification in this genus, clear evidence for homologous recombination at the single gene level is still lacking. Several Planktothrix strains (n = 244) sampled in eight fresh water lakes in north Italy were characterized by sequencing the rpoC1 gene, a molecular marker previously proposed to discriminate between species. Six haplotypes were detected, four of which are newly described. A relevant number of rpoC1 sequences (n = 54) showed evidence of homologous recombination. By comparing the sequences produced in the work presented here to those available on GenBank for the genus, multiple recombination events were tracked between haplotypes associated to P. rubescens, P. suspensa and P. agardhii, the latter being a species not found in our survey. Recombination signals were in form of (i) a vast mosaic structure present in the alignment of rpoC1 haplotypes, (ii) multiple and statistically significant paths in the split decomposition network connecting these haplotypes and (iii) many individual crossing-over events detected by means of recombination detection tests. Data suggest that the molecular evolution of the rpoC1 gene in the genus Planktothrix appears as strongly influenced by homologous recombination. In addition, rpoC1 diversity effectively tracks recombinational processes among species in the complex made up by P. rubescens, P. agardhii and P. suspensa, which are not isolated in terms of gene-flow.

Gene flow, recombination, and selection in cyanobacteria: population structure of geographically related Planktothrix freshwater strains

Several Planktothrix strains, each producing a distinct oligopeptide profile, have been shown to coexist within Lake Steinsfjorden (Norway). Using nonribosomal peptide synthetase (NRPS) genes as markers, it has been shown that the Planktothrix community comprises distinct genetic variants displaying differences in bloom dynamics, suggesting a Planktothrix subpopulation structure. Here, we investigate the Planktothrix variants inhabiting four lakes in southeast of Norway utilizing both NRPS and non-NRPS genes. Phylogenetic analyses showed similar topologies for both NRPS and non-NRPS genes, and the lakes appear to have similar structuring of Planktothrix genetic variants. The structure of distinct variants was also supported by very low genetic diversity within variants compared to the between-variant diversity. Incongruent topologies and split decomposition revealed recombination events between Planktothrix variants. In several strains the gene variants seem to be a result of recombination. Both NRPS and non-NRPS genes are dominated by purifying selection; however, sites subjected to positive selection were also detected. The presence of similar and well-separated Planktothrix variants with low internal genetic diversity indicates gene flow within Planktothrix populations. Further, the low genetic diversity found between lakes (similar range as within lakes) indicates gene flow also between Planktothrix populations and suggests recent, or recurrent, dispersals. Our data also indicate that recombination has resulted in new genetic variants. Stability within variants and the development of new variants are likely to be influenced by selection patterns and within-variant homologous recombination.

Recovery of viable cyanophages from the sediments of a eutrophic lake at decadal timescales

Cyanobacteria and their associated viruses (cyanophages) are abundant throughout the world in both marine and freshwater environments. The predator-prey relationship influences population dynamics within these ecosystems and facilitates the co-evolution of both organisms. Evidence of the close-linked interactions between cyanobacteria and viruses has been found extensively throughout marine biomes, but freshwater systems are less well studied. Eutrophic lake sediments potentially allow the preservation of cyanophages. If historic cyanophages could be isolated, they could provide insights into the evolution, biology and population dynamics over defined timescales. To determine whether viable phages are present in this environment, sectioned sediment cores (~ 50 cm in length) were taken from a eutrophic, stratifying lake (Rostherne Mere, Cheshire, UK). They were examined under the transmission electron microscope, and phages were isolated on two Microcystis strains PCC 7820 and BC 84/1. Viable phages were recovered from ~ 33- and ~ 50-year-old sediments. This is the first known study to investigate the viability of freshwater cyanophages recovered from dated lake sediments.

Cyanophage infection in the bloom-forming cyanobacteria Microcystis aeruginosa in surface freshwater

Host-like genes are often found in viral genomes. To date, multiple host-like genes involved in photosynthesis and the pentose phosphate pathway have been found in phages of marine cyanobacteria Synechococcus and Prochlorococcus. These gene products are predicted to redirect host metabolism to deoxynucleotide biosynthesis for phage replication while maintaining photosynthesis. A cyanophage, Ma-LMM01, infecting the toxic cyanobacterium Microcystis aeruginosa, was isolated from a eutrophic freshwater lake and assigned as a member of a new lineage of the Myoviridae family. The genome encodes a host-like NblA. Cyanobacterial NblA is known to be involved in the degradation of the major light harvesting complex, the phycobilisomes. Ma-LMM01 nblA gene showed an early expression pattern and was highly transcribed during phage infection. We speculate that the co-option of nblA into Microcystis phages provides a significant fitness advantage to phages by preventing photoinhibition during infection and possibly represents an important part of the co-evolutionary interactions between cyanobacteria and their phages.

Extreme conservation and non-neutral evolution of the cpmA Circadian locus in a globally distributed Chroococcidiopsis sp. from naturally stressful habitats

Cyanobacteria are among the most ancient organisms known to have circadian rhythms. The cpmA gene is involved in controlling the circadian output signal. We studied polymorphism and divergence of this gene in six populations of a stress-tolerant cyanobacterium, Chroococcidiopsis sp., sampled in extreme habitats across the globe. Despite high haplotype diversity (0.774), nucleotide diversity of cpmA is very low (π = 0.0034): the gene appears to be even more conserved than housekeeping genes. Even though the populations were sampled thousands kilometers apart, they manifested virtually no genetic differentiation at this locus (F(ST) = 0.0228). Using various tests for neutrality, we determined that evolution of cpmA significantly departures from the neutral model and is governed by episodic positive selection.

Intricate interactions between the bloom-forming cyanobacterium Microcystis aeruginosa and foreign genetic elements, revealed by diversified clustered regularly interspaced short palindromic repeat (CRISPR) signatures

Clustered regularly interspaced short palindromic repeats (CRISPR) confer sequence-dependent, adaptive resistance in prokaryotes against viruses and plasmids via incorporation of short sequences, called spacers, derived from foreign genetic elements. CRISPR loci are thus considered to provide records of past infections. To describe the host-parasite (i.e., cyanophages and plasmids) interactions involving the bloom-forming freshwater cyanobacterium Microcystis aeruginosa, we investigated CRISPR in four M. aeruginosa strains and in two previously sequenced genomes. The number of spacers in each locus was larger than the average among prokaryotes. All spacers were strain specific, except for a string of 11 spacers shared in two closely related strains, suggesting diversification of the loci. Using CRISPR repeat-based PCR, 24 CRISPR genotypes were identified in a natural cyanobacterial community. Among 995 unique spacers obtained, only 10 sequences showed similarity to M. aeruginosa phage Ma-LMM01. Of these, six spacers showed only silent or conservative nucleotide mutations compared to Ma-LMM01 sequences, suggesting a strategy by the cyanophage to avert CRISPR immunity dependent on nucleotide identity. These results imply that host-phage interactions can be divided into M. aeruginosa-cyanophage combinations rather than pandemics of population-wide infectious cyanophages. Spacer similarity also showed frequent exposure of M. aeruginosa to small cryptic plasmids that were observed only in a few strains. Thus, the diversification of CRISPR implies that M. aeruginosa has been challenged by diverse communities (almost entirely uncharacterized) of cyanophages and plasmids.

The distribution of a phage-related insertion sequence element in the cyanobacterium, Microcystis aeruginosa

The cyanophage Ma-LMM01, specifically-infecting Microcystis aeruginosa, has an insertion sequence (IS) element that we named IS607-cp showing high nucleotide similarity to a counterpart in the genome of the cyanobacterium Cyanothece sp. We tested 21 strains of M. aeruginosa for the presence of IS607-cp using PCR and detected the element in strains NIES90, NIES112, NIES604, and RM6. Thermal asymmetric interlaced PCR (TAIL-PCR) revealed each of these strains has multiple copies of IS607-cp. Some of the ISs were classified into three types based on their inserted positions; IS607-cp-1 is common in strains NIES90, NIES112 and NIES604, whereas IS607-cp-2 and IS607-cp-3 are specific to strains NIES90 and RM6, respectively. This multiplicity may reflect the replicative transposition of IS607-cp. The sequence of IS607-cp in Ma-LMM01 showed robust affinity to those found in M. aeruginosa and Cyanothece spp. in a phylogenetic tree inferred from counterparts of various bacteria. This suggests the transfer of IS607-cp between the cyanobacterium and its cyanophage. We discuss the potential role of Ma-LMM01-related phages as donors of IS elements that may mediate the transfer of IS607-cp; and thereby partially contribute to the genome plasticity of M. aeruginosa.

Lack of phylogeographic structure in the freshwater cyanobacterium Microcystis aeruginosa suggests global dispersal

BACKGROUND

Free-living microorganisms have long been assumed to have ubiquitous distributions with little biogeographic signature because they typically exhibit high dispersal potential and large population sizes. However, molecular data provide contrasting results and it is far from clear to what extent dispersal limitation determines geographic structuring of microbial populations. We aimed to determine biogeographical patterns of the bloom-forming freshwater cyanobacterium Microcystis aeruginosa. Being widely distributed on a global scale but patchily on a regional scale, this prokaryote is an ideal model organism to study microbial dispersal and biogeography.

METHODOLOGY/PRINCIPAL FINDINGS

The phylogeography of M. aeruginosa was studied based on a dataset of 311 rDNA internal transcribed spacer (ITS) sequences sampled from six continents. Richness of ITS sequences was high (239 ITS types were detected). Genetic divergence among ITS types averaged 4% (maximum pairwise divergence was 13%). Preliminary analyses revealed nearly completely unresolved phylogenetic relationships and a lack of genetic structure among all sequences due to extensive homoplasy at multiple hypervariable sites. After correcting for this, still no clear phylogeographic structure was detected, and no pattern of isolation by distance was found on a global scale. Concomitantly, genetic differentiation among continents was marginal, whereas variation within continents was high and was mostly shared with all other continents. Similarly, no genetic structure across climate zones was detected.

CONCLUSIONS/SIGNIFICANCE

The high overall diversity and wide global distribution of common ITS types in combination with the lack of phylogeographic structure suggest that intercontinental dispersal of M. aeruginosa ITS types is not rare, and that this species might have a truly cosmopolitan distribution.

Local expansion of a panmictic lineage of water bloom-forming cyanobacterium Microcystis aeruginosa

In previous studies, we have demonstrated that the population structure of the bloom-forming cyanobacterium Microcystis aeruginosa is clonal. Expanded multilocus sequence typing analysis of M. aeruginosa using 412 isolates identified five intraspecific lineages suggested to be panmictic while maintaining overall clonal structure probably due to a reduced recombination rate between lineages. Interestingly, since 2005 most strains belonging to one of these panmictic clusters (group G) have been found in a particular locality (Lake Kasumigaura Basin) in Japan. In this locality, multiple, similar but distinct genotypes of this lineage predominated in the bloom, a pattern that is unprecedented for M. aeruginosa. The population structure underlying blooms associated with this lineage is comparable to epidemics of pathogens. Our results may reveal an expansion of the possible adaptive lineage in a localized aquatic environment, providing us with a unique opportunity to investigate its ecological and biogeographical consequences.

Spatial divergence in the proportions of genes encoding toxic peptide synthesis among populations of the cyanobacterium Planktothrix in European lakes

It has been frequently reported that seasonal changes in toxin production by cyanobacteria are due to changes in the proportion of toxic/nontoxic genotypes in parallel to increases or decreases in population density during the seasonal cycle of bloom formation. In order to find out whether there is a relationship between the proportion of genes encoding toxic peptide synthesis and population density of Planktothrix spp. we compared the proportion of three gene regions that are indicative of the synthesis of the toxic heptapeptide microcystin (mcyB), and the bioactive peptides aeruginoside (aerB) and anabaenopeptin (apnC) in samples from 23 lakes of five European countries (n=153). The mcyB, aerB, and apnC genes occurred in 99%, 99%, and 97% of the samples, respectively, and on average comprised 60 ± 3%, 22 ± 2%, and 54 ± 4% of the total population, respectively. Although the populations differed widely in abundance (10⁻³–10³ mm³ L⁻¹) no dependence of the proportion of the mcyB, aerB, and apnC genes on the density of the total population was found. In contrast populations differed significantly in their average mcyB, aerB, and apnC gene proportions, with no change between prebloom and bloom conditions. These results emphasize stable population-specific differences in mcyB, aerB, and apnC proportions that are independent from seasonal influences.

Characterization of epidemiologically unrelated Acinetobacter baumannii isolates from four continents by use of multilocus sequence typing, pulsed-field gel electrophoresis, and sequence-based typing of bla(OXA-51-like) genes

This study used a diverse collection of epidemiologically unrelated Acinetobacter baumannii isolates to compare the robustness of a multilocus sequence typing (MLST) scheme, based on conserved regions of seven housekeeping genes, gltA, gdhB, recA, cpn60, rpoD, gyrB, and gpi, with that of sequence-based typing of bla(OXA-51-like) genes (SBT-bla(OXA-51-like) genes). The data obtained by analysis of MLST and SBT-bla(OXA-51-like) genes were compared to the data generated by pulsed-field gel electrophoresis (PFGE). The topologies of the phylogenetic trees generated for the gyrB and gpi genes showed evidence of recombination and were inconsistent with those of the trees generated for the other five genes. MLST identified 24 sequence types (STs), of which 19 were novel, and 5 novel alleles. Clonality was demonstrated by eBURST analysis and standardized index of association values of >1 (P < 0.001). MLST data revealed that all isolates harboring the major bla(OXA-51-like) alleles OXA-66, OXA-69, and OXA-71 fell within the three major European clonal lineages. However, the MLST data were not always in concordance with the PFGE data, and some isolates containing the same bla(OXA-51-like) allele demonstrated <50% relatedness by PFGE. It was concluded that the gyrB and gpi genes are not good candidates for use in MLST analysis and that a SBT-bla(OXA-51-like) gene scheme produced results comparable to those produced by MLST for the identification of the major epidemic lineages, with the advantage of having a significantly reduced sequencing cost and time. It is proposed that studies of A. baumannii epidemiology could involve initial screening of bla(OXA-51-like) alleles to identify isolates belonging to major epidemic lineages, followed by MLST analysis to categorize isolates from common lineages, with PFGE being reserved for fine-scale typing.

An rpoB signature sequence provides unique resolution for the molecular typing of cyanobacteria

The use of morphological characters for the classification of cyanobacteria has often led to ambiguous strain assignment. In the past two decades, the availability of sequences, such as those of the 16S rRNA, nif, cpc and rpoC1 genes, and the use of metagenomics, has steadily increased and has made the reconstruction of evolutionary relationships of some cyanobacterial groups possible in addition to improving strain assignment. Conserved indels (insertions/deletions) are present in all cyanobacterial RpoB (β subunit of RNA polymerase) sequences presently available in public databases. These indels are located in the Rpb2_6 domain of RpoB, which is involved in DNA binding and DNA-directed RNA polymerase activity. They are variable in length (6-44 aa) and sequence, and form part of what appears to be a longer signature sequence (43-81 aa). Indeed, a number of these sequences turn out to be distinctive among several strains of a given genus and even among strains of a given species. These signature sequences can thus be used to identify cyanobacteria at a subgenus level and can be useful molecular markers to establish the taxonomic positions of cyanobacterial isolates in laboratory cultures, and/or to assess cyanobacterial biodiversity in space and time in natural ecosystems.