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

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.

Cyanometabolites: molecules with immense antiviral potential

Cyanometabolites are active compounds derived from cyanobacteria that include small low molecular weight peptides, oligosaccharides, lectins, phenols, fatty acids, and alkaloids. Some of these compounds may pose a threat to human and environment. However, majority of them are known to have various health benefits with antiviral properties against pathogenic viruses including Human immunodeficiency virus (HIV), Ebola virus (EBOV), Herpes simplex virus (HSV), Influenza A virus (IAV) etc. Cyanometabolites classified as lectins include scytovirin (SVN), Oscillatoria agardhii agglutinin (OAAH), cyanovirin-N (CV-N), Microcystis viridis lectin (MVL), and microvirin (MVN) also possess a potent antiviral activity against viral diseases with unique properties to recognize different viral epitopes. Studies showed that a small linear peptide, microginin FR1, isolated from a water bloom of Microcystis species, inhibits angiotensin-converting enzyme (ACE), making it useful for the treatment of coronavirus disease 2019 (COVID-19). Our review provides an overview of the antiviral properties of cyanobacteria from the late 90s till now and emphasizes the significance of their metabolites in combating viral diseases, particularly severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has received limited attention in previous publications. The enormous medicinal potential of cyanobacteria is also emphasized in this review, which justifies their use as a dietary supplement to fend off pandemics in future.

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.

Mojavia aguilerae and M. dolomitestris - two new Nostocaceae (Cyanobacteria) species from the Americas

While nostocacean cyanobacteria are ubiquitous and play critical roles in terrestrial ecosystems, their taxonomy and biogeography still entail mysteries. We isolated two Nostoc-like cyanobacteria from biological soil crusts of the Atacama (Chile) and Mojave (USA) Deserts. An initial 16S rRNA gene phylogeny placed both in monophyly with Mojavia pulchra. Here, we describe two new species of the previously monotypic Mojavia using a polyphasic approach including morphology, 16S rRNA phylogenies, secondary structure, and percent similarity of the 16S-23S ITS region. Like M. pulchra, both new species produce compact microcolonies, arthrospore-like akinetes, and monocytes, traits characteristic of the genus. Mojavia aguilerae sp. nov. is morphologically distinct from both other species in producing bluntly conical end cells, abundant enlarged akinetes in multiseriate filaments, and gold-colored cells during senescence. Mojavia dolomitestris sp. nov. exhibited distinctly firm, light-colored, compartmentalized mucilage. M. dolomitestris is somewhat cryptic with M. pulchra, but has more densely packed microcolonies, rarity and later onset of brownish sheath pigmentation, and an origin from soils derived from dolomite. The two new species strengthened the position of Mojavia as a robust genus sister to Nostoc. Although 16S rRNA gene data could not separate the Mojavia species from each other, the three species showed distinct dissimilarities in secondary ITS structure and differed greatly from Nostoc sensu stricto. The high dissimilarities between their 16S-23S ITS regions suggest a long evolutionary history of the three species as separate lineages. Mojavia is an evolutionary and ecologically unique nostocacean genus, and its rarity and restricted habitat point to an urgent need for recognition and protection.

Molecular Phylogenetics, Speciation, and Long Distance Dispersal in Tardigrade Evolution: A case study of the genusMilnesium

Microorganisms (sensu lato, i.e., including micrometazoans) are thought to have cosmopolitan geographic distributions due to their theoretically unlimited dispersal capabilities, a consequence of their tiny size, population dynamics, and resistant forms. However, several molecular studies of microorganisms have identified biogeographic patterns indicating cryptic speciation and/or weak species definitions. Using a multi-locus approach with the genus Milnesium (Tardigrada), we aimed to determine the genetic structure of populations worldwide and the effects of long distance dispersal (LDD) on genetic connectivity and relationships across the six continents. Our results on this micrometazoan's genetic structure and LDD at global and micro-local scales indicate contrasting patterns not easily explained by a unique or simple phenomenon. Overall, we report three key findings: (i) confirmation of long distance dispersal for tardigrades, (ii) populations with globally-shared or endemic micro-local haplotypes, and (iii) a supported genetic structure instead of the homogeneous genetic distribution hypothesized for microorganisms with LDD capabilities. Moreover, incongruences between our morphological and molecular results suggest that species delimitation within the genus Milnesium could be problematic due to homoplasy. Duality found for Milnesium populations at the global scale, namely, a molecular phylogenetic structure mixed with widely distributed haplotypes (but without any apparent biogeographic structure), is similar to patterns observed for some unicellular, prokaryotic and eukaryotic, microorganisms. Factors influencing these patterns are discussed within an evolutionary framework.

The rates of global bacterial and archaeal dispersal

The phylogenetic resolution at which microorganisms display geographic endemism, the rates at which they disperse at global scales, and the role of humans on global microbial dispersal are largely unknown. Answering these questions is necessary for interpreting microbial biogeography, ecology, and macroevolution and for predicting the spread of emerging pathogenic strains. To resolve these questions, I analyzed the geographic and evolutionary relationships between 36,795 bacterial and archaeal ("prokaryotic") genomes from ∼7000 locations around the world. I find clear signs of continental-scale endemism, including strong correlations between phylogenetic divergence and geographic distance. However, the phylogenetic scale at which endemism generally occurs is extremely small, and most "species" (defined by an average nucleotide identity ≥ 95%) and even closely related strains (average nucleotide identity ≥ 99.9%) are globally distributed. Human-associated lineages display faster dispersal rates than other terrestrial lineages; the average net distance between any two human-associated cell lineages diverging 50 years ago is roughly 580 km. These results suggest that many previously reported global-scale microbial biogeographical patterns are likely the result of recent or current environmental filtering rather than geographic endemism. For human-associated lineages, estimated transition rates between Europe and North America are particularly high, and much higher than for non-human associated terrestrial lineages, highlighting the role that human movement plays in global microbial dispersal. Dispersal was slowest for hot spring- and terrestrial subsurface-associated lineages, indicating that these environments may act as "isolated islands" of microbial evolution.

Algal and Cyanobacterial Lectins and Their Antimicrobial Properties

Lectins are proteins with a remarkably high affinity and specificity for carbohydrates. Many organisms naturally produce them, including animals, plants, fungi, protists, bacteria, archaea, and viruses. The present report focuses on lectins produced by marine or freshwater organisms, in particular algae and cyanobacteria. We explore their structure, function, classification, and antimicrobial properties. Furthermore, we look at the expression of lectins in heterologous systems and the current research on the preclinical and clinical evaluation of these fascinating molecules. The further development of these molecules might positively impact human health, particularly the prevention or treatment of diseases caused by pathogens such as human immunodeficiency virus, influenza, and severe acute respiratory coronaviruses, among others.

Genotyping and multivariate regression trees reveal ecological diversification within the Microcystis aeruginosa complex along a wide environmental gradient

Addressing the ecological and evolutionary processes underlying biodiversity patterns is essential to identify the mechanisms shaping community structure and function. In bacteria, the formation of new ecologically distinct populations (ecotypes) is proposed as one of the main drivers of diversification. New ecotypes arise when mutations in key functional genes or acquisition of new metabolic pathways by horizontal gene transfer allow the population to exploit new resources, permitting their coexistence with the parental population. We previously reported the presence of microcystin-producing organisms of the Microcystis aeruginosa complex (toxic MAC) through an 800 km environmental gradient ranging from freshwater to estuarine-marine waters in South America. We hypothesize that the success of toxic MAC in such a gradient is due to the existence of very closely related populations that are ecologically distinct (ecotypes), each specialized to a specific arrangement of environmental variables. Here, we analyzed toxic MAC genetic diversity through qPCR and high-resolution melting analysis (HRMA) of a functional gene (mcyJ, microcystin synthetase cluster). We explored the variability of the mcyJ gene along the environmental gradient by multivariate classification and regression trees (mCART). Six groups of mcyJ genotypes were distinguished and associated with different combinations of water temperature, conductivity and turbidity. We propose that each mcyJ variant associated to a defined environmental condition is an ecotype (or species) whose relative abundances vary according to their fitness in the local environment. This mechanism would explain the success of toxic MAC in such a wide array of environmental conditions. Importance Organisms of the Microcystis aeruginosa Complex form harmful algal blooms (HABs) in nutrient-rich water bodies worldwide. MAC HABs are difficult to manage owing to the production of potent toxins (microcystins) that resist water treatment. Besides, the role of microcystins in the ecology of MAC organisms is still elusive, meaning that the environmental conditions driving the toxicity of the bloom are not clear. Furthermore, the lack of coherence between morphology-based and genomic-based species classification makes it difficult to draw sound conclusions about when and where each member species of the MAC will dominate the bloom. Here, we propose that the diversification process and success of toxic MAC in a wide range of waterbodies involves the generation of ecotypes, each specialized in a particular niche, whose relative abundance varies according to its fitness in the local environment. This knowledge can improve the generation of accurate prediction models of MAC growth and toxicity, helping to prevent human and animal intoxication.

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.

Microdiversity and phylogeographic diversification of bacterioplankton in pelagic freshwater systems revealed through long-read amplicon sequencing

BACKGROUND

Freshwater ecosystems are inhabited by members of cosmopolitan bacterioplankton lineages despite the disconnected nature of these habitats. The lineages are delineated based on > 97% 16S rRNA gene sequence similarity, but their intra-lineage microdiversity and phylogeography, which are key to understanding the eco-evolutional processes behind their ubiquity, remain unresolved. Here, we applied long-read amplicon sequencing targeting nearly full-length 16S rRNA genes and the adjacent ribosomal internal transcribed spacer sequences to reveal the intra-lineage diversities of pelagic bacterioplankton assemblages in 11 deep freshwater lakes in Japan and Europe.

RESULTS

Our single nucleotide-resolved analysis, which was validated using shotgun metagenomic sequencing, uncovered 7-101 amplicon sequence variants for each of the 11 predominant bacterial lineages and demonstrated sympatric, allopatric, and temporal microdiversities that could not be resolved through conventional approaches. Clusters of samples with similar intra-lineage population compositions were identified, which consistently supported genetic isolation between Japan and Europe. At a regional scale (up to hundreds of kilometers), dispersal between lakes was unlikely to be a limiting factor, and environmental factors or genetic drift were potential determinants of population composition. The extent of microdiversification varied among lineages, suggesting that highly diversified lineages (e.g., Iluma-A2 and acI-A1) achieve their ubiquity by containing a consortium of genotypes specific to each habitat, while less diversified lineages (e.g., CL500-11) may be ubiquitous due to a small number of widespread genotypes. The lowest extent of intra-lineage diversification was observed among the dominant hypolimnion-specific lineage (CL500-11), suggesting that their dispersal among lakes is not limited despite the hypolimnion being a more isolated habitat than the epilimnion.

CONCLUSIONS

Our novel approach complemented the limited resolution of short-read amplicon sequencing and limited sensitivity of the metagenome assembly-based approach, and highlighted the complex ecological processes underlying the ubiquity of freshwater bacterioplankton lineages. To fully exploit the performance of the method, its relatively low read throughput is the major bottleneck to be overcome in the future. Video abstract.

Genetic diversity of Calliphora vicina (Diptera: Calliphoridae) in the Iberian Peninsula based on cox1, 16S and ITS2 sequences

The study of Diptera at the scene of a crime can provide essential information for the interpretation of evidence. Phylogeographic reconstruction could help differentiate haplotypes of a dipteran species in a geographical area, clarifying, for example, the details of a possible relocation of a corpse. In addition, inferring the ancestral areas of distribution helps to understand the current status of the species and its biogeographic history. One of the most important species in forensic entomology is Calliphora vicina Rovineau-Desvoidy, 1830 (Diptera: Calliphoridae). The aim of this work is to increase our knowledge of this species in the Iberian Peninsula using 464 specimens from Spain and Portugal. These samples were identified using morphological keys and by molecular methods using fragments of the cox1, 16S and ITS2 genes. The phylogeographic history of these populations was inferred from haplotype networks and the reconstruction of ancestral areas of distribution. The molecular results corroborated the morphological identifications of the samples. Phylogeographic networks showed no geographical structure, as haplotypes are shared among almost all populations. reconstruct ancestral state in phylogenies analyses showed a high rate of movement among populations, possibly related to human activity. These results suggest that this species had a very rapid and recent spatial and demographic expansion throughout the Iberian Peninsula.

Effect of Zinc on Microcystis aeruginosa UTEX LB 2385 and Its Toxin Production

Cyanobacteria harmful algal blooms (CHABs) are primarily caused by man-made eutrophication and increasing climate-change conditions. The presence of heavy metal runoff in affected water systems may result in CHABs alteration to their ecological interactions. Certain CHABs produce by-products, such as microcystin (MC) cyanotoxins, that have detrimentally affected humans through contact via recreation activities within implicated water bodies, directly drinking contaminated water, ingesting biomagnified cyanotoxins in seafood, and/or contact through miscellaneous water treatment. Metallothionein (MT) is a small, metal-sequestration cysteine rich protein often upregulated within the stress response mechanism. This study focused on zinc metal resistance and stress response in a toxigenic cyanobacterium, Microcystis aeruginosa UTEX LB 2385, by monitoring cells with (0, 0.1, 0.25, and 0.5 mg/L) ZnCl₂ treatment. Flow cytometry and phase contrast microscopy were used to evaluate physiological responses in cultures. Molecular assays and an immunosorbent assay were used to characterize the expression of MT and MC under zinc stress. The results showed that the half maximal inhibitory concentration (IC₅₀) was 0.25 mg/L ZnCl₂. Flow cytometry and phase contrast microscopy showed morphological changes occurred in cultures exposed to 0.25 and 0.5 mg/L ZnCl₂. Quantitative PCR (qPCR) analysis of selected cDNA samples showed significant upregulation of Mmt through all time points, significant upregulation of mcyC at a later time point. ELISA MC-LR analysis showed extracellular MC-LR (µg/L) and intracellular MC-LR (µg/cell) quota measurements persisted through 15 days, although 0.25 mg/L ZnCl₂ treatment produced half the normal cell biomass and 0.5 mg/L treatment largely inhibited growth. The 0.25 and 0.5 mg/L ZnCl₂ treated cells demonstrated a ~40% and 33% increase of extracellular MC-LR(µg/L) equivalents, respectively, as early as Day 5 compared to control cells. The 0.5 mg/L ZnCl₂ treated cells showed higher total MC-LR (µg/cell) quota yield by Day 8 than both 0 mg/L ZnCl₂ control cells and 0.1 mg/L ZnCl₂ treated cells, indicating release of MCs upon cell lysis. This study showed this Microcystis aeruginosa strain is able to survive in 0.25 mg/L ZnCl₂ concentration. Certain morphological zinc stress responses and the upregulation of mt and mcy genes, as well as periodical increased extracellular MC-LR concentration with ZnCl₂ treatment were observed.

Using distribution models to estimate blooms of phytosanitary cyanobacteria in Brazil

Abstract: The multiple uses of aquatic ecosystems by humankind and the continuous interference of their activities have contributed to the emergence of potentially toxic cyanobacteria blooms. Here, we firstly created a database of occurrences of cyanobacteria blooms in Brazil through a systematic review of the scientific literature available in online platforms (e.g. Web of Science, Capes Thesis Catalogue). Secondly, we carried out ecological niche models with occurrence data obtained from these studies to predict climatically suitable areas for blooms. We select 21 bioclimatic variables input environmental data. We used five modeling methods for the current climate scenario: (1) Maxent; (2) Support Vector Machines; (3) Random Forest; (4) Maximum Likelihood e (5) Gaussian. We found that the number of publications about bloom events was higher in 2009 with a decline in the years 2012, 2013 and 2017. Furthermore, the years with the higher records of blooms in freshwater environments were 2005, 2011 e 2014. These events occurring mainly in public supply reservoirs and are mostly of the genera Microcystis Lemmermann, 1907, Dolichospermum (Ralfs ex Bornet & Flahault) P.Wacklin, L.Hoffmann & J.Komárek, 2009 and Raphidiopsis F.E.Fritsch & F.Rich, 1929. Modeling the potential distribution of blooms, we found sampling gaps that should be targeting for future researches, especially in the Amazon biome. Overall, the models did not predict highly suitable areas in the /north of Brazil, while other regions were relatively well distributed with a higher number of occurrence records in the Southeast region.

Latitudinal gradient of cyanobacterial diversity in tidal flats

Latitudinal diversity gradients are well-known for plants and animals, but only recently similar patterns have been described for some specific microbial communities in distinct habitats. Although microbial diversity is well-investigated worldwide, most of the studies are spatially too restricted to allow general statements about global diversity patterns. Additionally, methodological differences make it hard and often impossible to compare several studies. This study investigated the cyanobacterial diversity in tidal flats along geographical and ecological gradients based on high-throughput sequencing of 16S rRNA gene fragments (Illumina MiSeq) and environmental data on a large spatial scale from the subtropics to the Arctic. Latitude and strongly correlated environmental parameters (e.g. temperature) were identified as important drivers of cyanobacterial diversity on global scale resulting in a latitudinal diversity gradient similar to that known from plants and animals. Other non-correlated parameters (e.g. grain size) were shown to be more important on local scales, although no consistent pattern occurred across different locations. Among a total number of 989 operational taxonomic units (OTUs) only one cosmopolitan (classified as Coleofasciculus chthonoplastes), but many location-specific and putative endemic ones (78%) were detected. High proportions of rare members of the community (up to 86%) were found in all samples. Phylogenetic beta diversity was shown to be influenced by the developmental stage of the mat community becoming increasingly similar with increasing stabilization.

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.

Intraspecific variability in response to phosphorus depleted conditions in the cyanobacteria Microcystis aeruginosa and Raphidiopsis raciborskii

Phosphorus loading plays an important role in the occurrence of cyanobacterial blooms and understanding how this nutrient affects the physiology of cyanobacteria is imperative to manage these phenomena. Microcystis aeruginosa and Raphidiopsis raciborskii are cyanobacterial species that form potentially toxic blooms in freshwater ecosystems worldwide. Blooms comprise numerous strains with high trait variability, which can contribute to the widespread distribution of these species. Here, we explored the intraspecific variability in response to phosphorus depleted conditions (P-) testing five strains of each species. Strains could be differentiated by cell volume or genetic profiles except for those of the same species, sampling location and date, though these presented differences in their response to (P-). Although differently affected by (P-) over 10 days, all strains were able to grow and maintain photosynthetic activity. For most M. aeruginosa and R. raciborskii strains growth rates were not significantly different comparing (P+) and (P-) conditions. After ten days in (P-), only one M. aeruginosa strain and two R. raciborskii strains showed reduction in biovolume yield as compared to (P+) but in most strains chlorophyll-a concentrations were lower in (P-) than in (P+). Reduced photosystem II efficiency was found for only one R. raciborskii strain while all M. aeruginosa strains were affected. Only two M. aeruginosa and one R. raciborskii strain increased alkaline phosphatase activity under (P-) as compared to (P+). Variation in P-uptake was also observed but comparison among strains yielded homogeneous groups comprised of representatives of both species. Comparing the response of each species as a whole, the (P-) condition affected growth rate, biovolume yield and chlorophyll yield. However, these parameters revealed variation among strains of the same species to the extent that differences between M. aeruginosa and R. raciborskii were not significant. Taken together, these results do not support the idea that R. raciborskii, as a species, can withstand phosphorus limitation better than M. aeruginosa and also point that the level of intraspecific variation may preclude generalizations based on studies that use only one or few strains.

Global phylogeography of toxic cyanobacteria Moorea producens reveals distinct genetic partitioning influenced by Proterozoic glacial cycles

Lyngbya majuscula is a marine filamentous cyanobacteria belonging to the family Oscillatoriaceae. Recent phylogenetic analyses of L. majuscula have reclassified a subset of this species into various genera such as Moorea, Okeania and Dapis. From the genus Moorea, Moorea producens is a toxic invasive cyanobacterium that produces bioactive secondary metabolites that can cause severe inflammation and blistering. Despite the global distribution of M. producens, little information is available on their origin, patterns of dispersal and population structure. In this study, the spatial population structure of M. producens was investigated using near-complete 16S rRNA sequences. Analysis of the global population of M. producens by Isolation by Distance and STRUCTURE revealed two significantly distinct cosmopolitan populations that were separated by a genetic break. Lineage-specific divergence estimates of 147 cyanobacterial taxa, based on a relaxed molecular clock indicated the first global emergence of M. producens during the Mesoarchean and a subsequent global recolonization during the Mesoproterozoic period. We conclude that the genetic discontinuity between both cosmopolitan populations is attributed to refugia associated with Proterozoic glacial cycles.

Ecological meta-analysis of bloom-forming planktonic Cyanobacteria in Argentina

The aim of the present research was to summarize the main reasons that explain the distribution of harmful blooms of cyanobacteria in Argentina. It is a large territory with climates ranging from humid tropical to cold temperate. We performed a meta-analysis of the published data and information in technical reports published from 1945 to 2015, and included additional data from personal non-published studies. A total of 122 water bodies affected by planktonic cyanobacterial blooms were recorded and geo-referenced. The analysis showed that blooms, defined as events exceeding 5000 cells/mL, occurred in different types of water bodies, including shallow lakes, rivers, streams, reservoirs, estuaries and storage facilities. Maximum bloom abundance and species and ecological strategies (dispersive, scum-forming, nitrogen fixer) responsible for each event were related to the geographic and climatologic characteristics and type and origin of water bodies. The Puna and the Andean Patagonia eco-regions were mostly free of blooms. The most impaired aquatic systems were shallow lakes and reservoirs (46.7 and 24.6%, respectively). Deep lakes had no reports of blooms and rivers were mainly affected at the regulated reaches, with intensities generally decreasing downstream the dams. Besides, 74.3% of the blooms reported in Argentina exceeded WHO Alert Level 2 for drinking and bathing waters (100,000 cells/mL). Thirty-nine species, identified by Komárek's polyphasic approach to taxonomy, were responsible for the blooms. Microcystis aeruginosa, Dolichospermum spiroides, Dolichospermum circinale, Raphidiopsis mediterranea and Cylindrospermopsis raciborskii were frequently found participating in either mixed or single species blooms. The species distribution was associated with the eco-region and aquatic system typologies and affected by seasonality and climatological and geographic variables. The eco-strategies of cyanobacterial species showed stronger associations with the qualitative and quantitative indicators used in the meta-analysis, and appeared as useful tools for management measures.

Response of bloom-forming cyanobacterium Microcystis aeruginosa to 17β-estradiol at different nitrogen levels

Co-existence of cyanobacterial harmful algal blooms (CyanoHABs) and steroid estrogens (SEs) has been an increasing concern in eutrophic waters. The cellular responses and biodegradation of 17β-estradiol (E2) in cyanobacterium Microcystis aeruginosa were investigated at different nitrogen levels. During the 10-d experiment, the growth of M. aeruginosa was stimulated by 10-100 μg L^-1 of E2 at the lowest nitrogen level of 0.5 mg L^-1, whereas the presence of E2 inhibited the cyanobacterial growth at 5 mg L^-1 of nitrogen. With nitrogen concentration increased to 50 mg L^-1, the impact of E2 on levels of growth rate and chlorophyll a (Chla) alleviated. Exposure to E2 also promoted the superoxide dismutase activity of M. aeruginosa, coupled with cellular oxidative damage as indicated by the increasing malondialdehyde content. A sufficient nitrogen supply mitigated the oxidative stress of E2 through enhancing the synthesis of detoxification-related enzymes. Simultaneously, the secretion of tryptophan-like substances in loosely- and tightly-bound extracellular polymeric substances was triggered for adapting to an E2 addition in the short term. Moreover, significant biodegradation of E2 was observed, and the process followed a first-order kinetic reaction. The obtained half-lives decreased with nitrogen levels and ranged from 2.47 to 2.81 and 3.39-5.04 d, respectively, at 10 and 100 μg L^-1 of E2. These results provide a better understanding of the potential effects of SEs on CyanoHABs formation, as well as the important role of CyanoHABs on SEs removal in aquatic ecosystems, which should be fully considered in the control of combined pollution.

High diversity of thermophilic cyanobacteria in Rupite hot spring identified by microscopy, cultivation, single-cell PCR and amplicon sequencing

Genotypic and morphological diversity of cyanobacteria in the Rupite hot spring (Bulgaria) was investigated by means of optical microscopy, cultivation, single-cell PCR, and 16S rRNA gene amplicon sequencing. Altogether, 34 sites were investigated along the 71-39 °C temperature gradient. Analysis of samples from eight representative sites shown that Illumina, optical microscopy, and Roche 454 identified 72, 45 and 19% respective occurrences of all cumulatively present taxa. Optical microscopy failed to detect species of minor occurrence; whereas, amplicon sequencing technologies suffered from failed primer annealing and the presence of species with extensive extracellular polysaccharides production. Amplicon sequencing of the 16S rRNA gene V5-V6 region performed by Illumina identified the cyanobacteria most reliably to the generic level. Nevertheless, only the combined use of optical microscopy, cultivation and sequencing methods allowed for reliable estimate of the cyanobacterial diversity. Here, we show that Rupite hot-spring system hosts one of the richest cyanobacterial flora reported from a single site above 50 °C. Chlorogloeopsis sp. was the most abundant at the highest temperature (68 °C), followed by Leptolyngbya boryana, Thermoleptolyngbya albertanoae, Synechococcus bigranulatus, Oculatella sp., and Desertifilum sp. thriving above 60 °C, while Leptolyngbya geysericola, Geitlerinema splendidum, and Cyanobacterium aponinum were found above 50 °C.

Environmental factors influencing the quantitative distribution of microcystin and common potentially toxigenic cyanobacteria in U.S. lakes and reservoirs

Many species of cyanobacteria are capable of producing toxins and causing nuisance blooms, however response to environmental conditions is likely taxon-specific. Environmental factors influencing cyanobacterial composition and toxin production in lakes have been examined in many studies; yet are often confined to individual water bodies, or to a small number of systems within the same region. Here, data from the 2012 USEPA National Lakes Assessment are used to examine relationships between biovolume of common potentially-toxigenic cyanobacteria (Aphanizomenon spp., Cylindrospermopsis spp., Dolichospermum spp., Microcystis spp. and Planktothrix spp.) and environmental variables across the entire conterminous United States, and results are compared across nine distinct ecoregions. Total phosphorus and water clarity were identified as the most influential environmental factors correlated with phytoplankton community composition. The Northern, Southern and Temperate Plains ecoregions displayed the highest biovolumes of potentially toxigenic taxa on average, as well as highest mean concentrations of microcystin. In those three ecoregions, samples with microcystin concentrations greater than 1 ppb were primarily dominated by Planktothrix spp. while in all other ecoregions Dolichospermum spp. was the dominant genus. Canonical Correlation Analysis revealed a strong association between high microcystin concentrations and high nutrient concentrations (total nitrogen and total phosphorus), and between high microcystin concentrations and low percentage of watershed forest cover. Results from this study indicate that the likely occurrence of potentially toxigenic taxa in lakes and reservoirs is predictable on a biogeographical basis, depending on morphological and water quality characteristics. Data from this study may be useful to regional managers attempting to prevent or mitigate nuisance cyanobacterial blooms.

The potential for dispersal of microalgal resting cysts by migratory birds

Most microalgal species are geographically widespread, but little is known about how they are dispersed. One potential mechanism for long-distance dispersal is through birds, which may transport cells internally (endozoochory) and deposit them during, or in-between, their migratory stopovers. We hypothesize that dinoflagellates, in particular resting stages, can tolerate bird digestion; that bird temperature, acidity, and retention time negatively affect dinoflagellate viability; and that recovered cysts can germinate after passage through the birds' gut, contributing to species-specific dispersal of the dinoflagellates across scales. Tolerance of two dinoflagellate species (Peridiniopsis borgei, a warm-water species and Apocalathium malmogiense, a cold-water species) to Mallard gut passage was investigated using in vitro experiments simulating the gizzard and caeca conditions. The effect of in vitro digestion and retention time on cell integrity, cell viability, and germination capacity of the dinoflagellate species was examined targeting both their vegetative and resting stages. Resting stages (cysts) of both species were able to survive simulated bird gut passage, even if their survival rate and germination were negatively affected by exposure to acidic condition and bird internal temperature. Cysts of A. malmogiense were more sensitive than P. borgei to treatments and to the presence of digestive enzymes. Vegetative cells did not survive conditions of bird internal temperature and formed pellicle cysts when exposed to gizzard-like acid conditions. We show that dinoflagellate resting cysts serve as dispersal propagules through migratory birds. Assuming a retention time of viable cysts of 2-12 h to duck stomach conditions, cysts could be dispersed 150-800 km and beyond.

Influence of cyanobacteria, mixotrophic flagellates, and virioplankton size fraction on transcription of microcystin synthesis genes in the toxic cyanobacterium Microcystis aeruginosa

Toxic cyanobacteria such as Microcystis aeruginosa are a worldwide concern in freshwater reservoirs. Problems associated with their mass occurrence are predicted to increase in the future due to global warming. The hepatotoxic secondary metabolite microcystin is of particular concern in this context. This study aimed to determine whether co-occurring microorganisms influence the expression of microcystin biosynthesis genes. To this end, we performed cocultivation experiments and measured mcyB and mcyD transcripts in M. aeruginosa using RT-qPCR. We utilized representatives from three different plankton groups: the picocyanobacterium Synechococcus elongatus, the unicellular flagellate grazer Ochromonas danica, and virioplankton from two different lakes. The presence of S. elongatus significantly increased mcyB and mcyD transcription in M. aeruginosa. Cocultivation with the mixotrophic chrysophyte O. danica did not increase the transcription of mcyB and mcyD; in fact, mcyD transcripts decreased significantly. The virioplankton size fraction of environmental water samples induced a significant increase in mcyB and mcyD transcription when obtained from lakes with cyanobacterial blooms. Our results show that co-occurring microorganisms influence the expression of microcystin biosynthesis genes in M. aeruginosa.

A multilevel trait-based approach to the ecological performance of Microcystis aeruginosa complex from headwaters to the ocean

The Microcystis aeruginosa complex (MAC) clusters cosmopolitan and conspicuous harmful bloom-forming cyanobacteria able to produce cyanotoxins. It is hypothesized that low temperatures and brackish salinities are the main barriers to MAC proliferation. Here, patterns at multiple levels of organization irrespective of taxonomic identity (i.e. a trait-based approach) were analyzed. MAC responses from the intracellular (e.g. respiratory activity) to the ecosystem level (e.g. blooms) were evaluated in wide environmental gradients. Experimental results on buoyancy and respiratory activity in response to increased salinity (0-35) and a literature review of maximum growth rates under different temperatures and salinities were combined with field sampling from headwaters (800km upstream) to the marine end of the Rio de la Plata estuary (Uruguay-South America). Salinity and temperature were the major variables affecting MAC responses. Experimentally, freshwater MAC cells remained active for 24h in brackish waters (salinity=15) while colonies increased their flotation velocity. At the population level, maximum growth rate decreased with salinity and presented a unimodal exponential response with temperature, showing an optimum at 27.5°C and a rapid decrease thereafter. At the community and ecosystem levels, MAC occurred from fresh to marine waters (salinity 30) with a sustained relative increase of large mucilaginous colonies biovolume with respect to individual cells. Similarly, total biomass and, specific and morphological richness decreased with salinity while blooms were only detected in freshwater both at high (33°C) and low (11°C) temperatures. In brackish waters, large mucilaginous colonies presented advantages under osmotic restrictive conditions. These traits values have also been associated with higher toxicity potential. This suggest salinity or low temperatures would not represent effective barriers for the survival and transport of potentially toxic MAC under likely near future scenarios of increasing human impacts (i.e. eutrophication, dam construction and climate change).

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.

Biogeography of bloom-forming microcystin producing and non-toxigenic populations of Dolichospermum lemmermannii (Cyanobacteria)

In the last decades, the cyanobacterium Dolichospermum lemmermannii showed an increasing spread to Southern Europe, raising serious concerns due to its ability to produce cyanotoxins. The widening of its geographic distribution and the observation of strains showing high optimum temperature underline its ecological heterogeneity, suggesting the existence of different ecotypes. To investigate its biogeography, new isolates from different European water bodies, together with strains maintained by the Norwegian Institute for Water Research Culture Collection of Algae, were genetically characterised for the 16S rRNA gene and compared with strains obtained from public repositories. Geographic distance highly influenced the differentiation of genotypes, further suggesting the concurrent role of geographic isolation, physical barriers and environmental factors in promoting the establishment of phylogenetic lineages adapted to specific habitats. Differences among populations were also examined by morphological analysis and evaluating the toxic potential of single strains, which revealed the exclusive ability of North European strains to produce microcystins, whereas the populations in Southern Europe tested negative for a wide range of cyanotoxins. The high dispersion ability and the existence of toxic genotypes indicate the possible spread of harmful blooms in other temperate regions.

Impacts of the 2014 severe drought on the Microcystis bloom in San Francisco Estuary

The increased frequency and intensity of drought with climate change may cause an increase in the magnitude and toxicity of freshwater cyanobacteria harmful algal blooms (CHABs), including Microcystis blooms, in San Francisco Estuary, California. As the fourth driest year on record in San Francisco Estuary, the 2014 drought provided an opportunity to directly test the impact of severe drought on cyanobacteria blooms in SFE. A field sampling program was conducted between July and December 2014 to sample a suite of physical, chemical, and biological variables at 10 stations in the freshwater and brackish reaches of the estuary. The 2014 Microcystis bloom had the highest biomass and toxin concentration, earliest initiation, and the longest duration, since the blooms began in 1999. Median chlorophyll a concentration increased by 9 and 12 times over previous dry and wet years, respectively. Total microcystin concentration also exceeded that in previous dry and wet years by a factor of 11 and 65, respectively. Cell abundance determined by quantitative PCR indicated the bloom contained multiple potentially toxic cyanobacteria species, toxic Microcystis and relatively high total cyanobacteria abundance. The bloom was associated with extreme nutrient concentrations, including a 20-year high in soluble reactive phosphorus concentration and low to below detection levels of ammonium. Stable isotope analysis suggested the bloom varied with both inorganic and organic nutrient concentration, and used ammonium as the primary nitrogen source. Water temperature was a primary controlling factor for the bloom and was positively correlated with the increase in both total and toxic Microcystis abundance. In addition, the early initiation and persistence of warm water temperature coincided with the increased intensity and duration of the Microcystis bloom from the usual 3 to 4 months to 8 months. Long residence time was also a primary factor controlling the magnitude and persistence of the bloom, and was created by a 66% to 85% reduction in both the water inflow and diversion of water for agriculture during the summer. We concluded that severe drought conditions can lead to a significant increase in the abundance of Microcystis and other cyanobacteria, as well as their associated toxins.

Influence of temperature, mixing, and addition of microcystin-LR on microcystin gene expression in Microcystis aeruginosa

Cyanobacteria, such as the toxin producer Microcystis aeruginosa, are predicted to be favored by global warming both directly, through elevated water temperatures, and indirectly, through factors such as prolonged stratification of waterbodies. M. aeruginosa is able to produce the hepatotoxin microcystin, which causes great concern in freshwater management worldwide. However, little is known about the expression of microcystin synthesis genes in response to climate change-related factors. In this study, a new RT-qPCR assay employing four reference genes (GAPDH, gltA, rpoC1, and rpoD) was developed to assess the expression of two target genes (the microcystin synthesis genes mcyB and mcyD). This assay was used to investigate changes in mcyB and mcyD expression in response to selected environmental factors associated with global warming. A 10°C rise in temperature significantly increased mcyB expression, but not mcyD expression. Neither mixing nor the addition of microcystin-LR (10 μg L-1 or 60 μg L-1 ) significantly altered mcyB and mcyD expression. The expression levels of mcyB and mcyD were correlated but not identical.

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.

Biogeography of sulfur-oxidizing Acidithiobacillus populations in extremely acidic cave biofilms

Extremely acidic (pH 0–1.5) Acidithiobacillus-dominated biofilms known as snottites are found in sulfide-rich caves around the world. Given the extreme geochemistry and subsurface location of the biofilms, we hypothesized that snottite Acidithiobacillus populations would be genetically isolated. We therefore investigated biogeographic relationships among snottite Acidithiobacillus spp. separated by geographic distances ranging from meters to 1000s of kilometers. We determined genetic relationships among the populations using techniques with three levels of resolution: (i) 16S rRNA gene sequencing, (ii) 16S–23S intergenic transcribed spacer (ITS) region sequencing and (iii) multi-locus sequencing typing (MLST). We also used metagenomics to compare functional gene characteristics of select populations. Based on 16S rRNA genes, snottites in Italy and Mexico are dominated by different sulfur-oxidizing Acidithiobacillus spp. Based on ITS sequences, Acidithiobacillus thiooxidans strains from different cave systems in Italy are genetically distinct. Based on MLST of isolates from Italy, genetic distance is positively correlated with geographic distance both among and within caves. However, metagenomics revealed that At. thiooxidans populations from different cave systems in Italy have different sulfur oxidation pathways and potentially other significant differences in metabolic capabilities. In light of those genomic differences, we argue that the observed correlation between genetic and geographic distance among snottite Acidithiobacillus populations is partially explained by an evolutionary model in which separate cave systems were stochastically colonized by different ancestral surface populations, which then continued to diverge and adapt in situ.

The common bloom-forming cyanobacterium Microcystis is prone to a wide array of microbial antagonists

Many degraded waterbodies around the world are subject to strong proliferations of cyanobacteria - notorious for their toxicity, high biomass build-up and negative impacts on aquatic food webs - the presence of which puts serious limits on the human use of affected water bodies. Cyanobacterial blooms are largely regarded as trophic dead ends since they are a relatively poor food source for zooplankton. As a consequence, their population dynamics are generally attributed to changes in abiotic conditions (bottom-up control). Blooms however generally contain a vast and diverse community of micro-organisms of which some have shown devastating effects on cyanobacterial biomass. For Microcystis, one of the most common bloom-forming cyanobacteria worldwide, a high number of micro-organisms (about 120 taxa) including viruses, bacteria, microfungi, different groups of heterotrophic protists, other cyanobacteria and several eukaryotic microalgal groups are currently known to negatively affect its growth by infection and predation or by the production of allelopathic compounds. Although many of these specifically target Microcystis, sharp declines of Microcystis biomass in nature are only rarely assigned to these antagonistic microbiota. The commonly found strain specificity of their interactions may largely preclude strong antagonistic effects on Microcystis population levels but may however induce compositional shifts that can change ecological properties such as bloom toxicity. These highly specific interactions may form the basis of a continuous arms race (co-evolution) between Microcystis and its antagonists which potentially limits the possibilities for (micro)biological bloom control.

Distribution and Habitat Specificity of Potentially-Toxic Microcystis across Climate, Land, and Water Use Gradients

Toxic cyanobacterial blooms are a growing threat to freshwater bodies worldwide. In order for a toxic bloom to occur, a population of cells with the genetic capacity to produce toxins must be present together with the appropriate environmental conditions. In this study, we investigated the distribution patterns and phylogeny of potentially-toxic Microcystis (indicated by the presence and/or phylogeny of the mcyD and mcyA genes). Samples were collected from the water column of almost 60 water bodies across widely differing gradients of environmental conditions and land use in Israel. Potentially, toxic populations were common but not ubiquitous, detected in ~65% of the studied sites. Local environmental factors, including phosphorus and ammonia concentrations and pH, as well as regional conditions such as the distance from built areas and nature reserves, were correlated with the distribution of the mcyD gene. A specific phylogenetic clade of Microcystis, defined using the sequence of the mcyA gene, was preferentially associated with aquaculture facilities but not irrigation reservoirs. Our results reveal important environmental, geospatial, and land use parameters affecting the geographic distribution of toxinogenic Microcystis, suggesting non-random dispersal of these globally abundant toxic cyanobacteria.

Orthocaspase and toxin-antitoxin loci rubbing shoulders in the genome of Microcystis aeruginosa PCC 7806

Programmed cell death in multicellular organisms is a coordinated and precisely regulated process. On the other hand, in bacteria we have little clue about the network of interacting molecules that result in the death of a single cell within a population or the death of almost complete population, such as often observed in cyanobacterial blooms. With the recent discovery that orthocaspase MaOC1 of the cyanobacterium Microcystis aeruginosa is an active proteolytic enzyme, we have gained a possible hint about at least one step in the process, but the picture is far from complete. Interestingly, the genomic context of MaOC1 revealed the presence of multiple copies of genes that belong to toxin-antitoxin modules. It has been speculated that these also play a role in bacterial programmed cell death. The discovery of two components linked to cell death within the same genomic region could open new ways to deciphering the underlying mechanisms of cyanobacterial cell death.

Palaeoecology of a billion-year-old non-marine cyanobacterium from the Torridon Group and Nonesuch Formation

A new chroococcalean cyanobacterium is described from approximately 1-billion-year-old non-marine deposits of the Torridonian Group of Scotland and the Nonesuch Formation of Michigan, USA. Individual cells of the new microfossil, Eohalothece lacustrina gen. et sp. nov., are associated with benthic microbial biofilms, but the majority of samples are recovered in palynological preparations in the form of large, apparently planktonic colonies, similar to extant species of Microcystis. In the Torridonian, Eohalothece is associated with phosphatic nodules, and we have developed a novel hypothesis linking Eohalothece to phosphate deposition in ancient freshwater settings. Extant cyanobacteria can be prolific producers of extracellular microcystins, which are non-ribosomal polypeptide phosphatase inhibitors. Microcystins may have promoted the retention and concentration of sedimentary organic phosphate prior to mineralization of francolite and nodule formation. This has a further implication that the Torridonian lakes were nitrogen limited as the release of microcystins is enhanced under such conditions today. The abundance and wide distribution of Eohalothece lacustrina attests to the importance of cyanobacteria as oxygen-producing photoautotrophs in lacustrine ecosystems at the time of the Mesoproterozoic-Neoproterozoic transition.

Molecular Analysis of the Cyanobacterial Community in Gastric Contents of Egrets with Symptoms of Steatitis

Many deaths of wild birds that have drunk water contaminated with hepatotoxic microcystin-producing cyanobacteria have been reported. A mass death of egrets and herons with steatitis were found at the agricultural reservoir occurring cyanobacterial waterblooms. This study aimed to verify a hypothesis that the egrets and herons which died in the reservoir drink microcystin-producing cyanobacteria and microcystin involves in the cause of death as well as the symptoms of steatitis. The cyanobacterial community in gastric contents of egrets and herons that died from steatitis was assessed using cyanobacterial 16S rRNA-based terminal-restriction fragment length polymorphism (T-RFLP) profiling and a cyanobacterial 16S rRNA-based clone library analysis. In addition, PCR amplification of the mcyB-C region and the mcyG gene, involved in microcystin biosynthesis, was examined. The cyanobacterial community in the gastric contents of two birds showed a simplistic composition. A comparison of cyanobacterial T-RFLP profiling and cloned sequences suggested that the genus Microcystis predominated in both samples of egrets died. Although we confirmed that two egrets which died in the reservoir have taken in cyanobacterial waterblooms containing the genus Microcystis, no mcy gene was detected in both samples according to the mcy gene-based PCR analysis. This study is the first to show the profiling and traceability of a cyanobacterial community in the gastric contents of wild birds by molecular analysis. Additionally, we consider causing symptoms of steatitis in the dead egrets.

Phylogeography of Yersinia ruckeri reveals effects of past evolutionary events on the current strain distribution and explains variations in the global transmission of enteric redmouth (ERM) disease

Phylogeographic patterns and population genetic structure of Yersinia ruckeri, the pathological agent of enteric redmouth disease (ERM) in salmonids, were investigated on the basis of concatenated multiloci sequences from isolates of different phenotypes obtained between 1965 and 2009 from diverse areas and hosts. Sequence analyses revealed genetic differentiation among subpopulations with the largest genetic distance occurring between subpopulations of Europe and Canada and/or South America. Bayesian analysis indicated the presence of three ancestral population clusters. Mismatch distribution displayed signatures characteristic of changes in size due to demographic and spatial expansions in the overall Y. ruckeri population, and also in the geographically separate subpopulations. Furthermore, a weak signal of isolation by distance was determined. A significant positive correlation between genetic and geographical distances was observed. These results revealed that the population of Y. ruckeri has undergone both ancient and recent population changes that were probably induced by biogeography forces in the past and, much more recently, by adaptive processes forced by aquaculture expansion. These findings have important implications for future studies on Y. ruckeri population dynamics, on the potential role of genetic structure to explain variations in ERM transmission, and on the effect of past evolutionary events on current estimations of gene flow.

Genotypes of ITS region of rRNA in Microcystis (Cyanobacteria) populations in Erhai Lake (China) and their correlation with eutrophication level

Previous studies on spatiotemporal changes of Microcystis genotypes have shown that the existence and succession of dominant genotypes always occur in eutrophicated freshwater bodies. However, few studies have focused on the correlation between genotype composition and eutrophication level. In the present study, clone libraries of the internal transcribed spacer (ITS) of rrn operon were sequenced from Microcystis populations in Erhai Lake, a subtropical plateau lake in the preliminary eutrophication stage. The genotype composition of the Microcystis populations was highly variable at spatiotemporal scales, and 473 ITS genotypes were identified from the 800 ITS sequences obtained. However, no significantly dominant ITS genotypes existed in the lake. Comparison of Erhai Lake with four major lakes in China, namely, Taihu, Chaohu, Gucheng, and Shijiu Lakes, showed that the Microcystis ITS genotypes and genetic diversity were negatively correlated with eutrophication level. Extensive comparison of the Microcystis ITS genotypes from waters worldwide revealed that 440 ITS genotypes were unique to Erhai Lake, and no obvious phylogenetic correlations can be detected among the dominant genotypes from different water bodies. The high genetic diversity of the Microcystis populations in Erhai Lake may have resulted from the effect of the early stage of eutrophication.

Presence of potential toxin-producing cyanobacteria in an oligo-mesotrophic lake in Baltic Lake District, Germany: an ecological, genetic and toxicological survey

Massive developments of potentially toxic cyanobacteria in Lake Stechlin, an oligo-mesotrophic lake in the Baltic Lake District of Germany raised concerns about toxic contamination of these important ecosystems. Field samples in the phase of mass developments of cyanobacteria were used for genetic and toxicological analyses. Microcystins and microcystin genes were detected in field samples of the lake for the first time. However, the toxins were not produced by the dominant taxa (Dolichospermum circinale and Aphanizomenon flos-aquae) but by taxa, which were present only in low biomass in the samples (Microcystis cf. aeruginosa and Planktothrix rubescens). The phytoplankton successions during the study period revealed an increase of cyanobacterial populations. The findings contribute to the changes that have been investigated in Lake Stechlin since the mid-1990s. The possible reasons behind these developments may be climate change, special weather conditions and an increased nutrient pool.

Vectorial capacity and genetic diversity of Anopheles annularis (Diptera: Culicidae) mosquitoes in Odisha, India from 2009 to 2011

Anopheles annularis is one of the major vectors of malaria in Odisha, India. The present study was undertaken to determine the vectorial capacity and assess the genetic diversity of An. annularis collected from different endemic regions of Odisha. Mosquitoes were collected from thirteen endemic districts using standard entomological collection methods from 2009 to 2011. Sibling species of An. annularis were identified by PCR-RFLP and sequencing of D3 region of 28S ribosomal DNA (rDNA) region. Plasmodium falciparum (Pf) sporozoite rate and human blood fed percentage (HBF) were estimated by multiplex PCR using Pf and human specific primers. Genetic diversity of An. annularis was estimated by ISSR markers. Out of 1647 An. annularis collected, 1353 (82.15%) were collected by mechanical aspirators and 294 (17.85%) by light trap. 49 (2.97%) were positive for human blood and 18 (1.09%) were positive for Pf sporozoite. PCR-RFLP and sequencing analyses detected only An annularis A in the study areas. Overall genetic differentiation among An. annularis populations was moderate (FST=0.048) and showed significant correlation between genetic distance and geographic distance (r=0.882; P<0.05). Angul population proved to be genetically unique and was highly divergent FST>0.110) from other populations, suggesting low gene flow between them. The study indicated that only An. annularis A was found in Odisha with potential vectorial capacity that can play a major role in malaria transmission. ISSR markers proved to be useful molecular tools to evaluate genetic variability in An. annularis populations.

Sporadic distribution and distinctive variations of cylindrospermopsin genes in cyanobacterial strains and environmental samples from Chinese freshwater bodies

Increasing reports of cylindrospermopsins (CYNs) in freshwater ecosystems have promoted the demand for identifying all of the potential CYN-producing cyanobacterial species. The present study explored the phylogenetic distribution and evolution of cyr genes in cyanobacterial strains and water samples from China. Four Cylindrospermopsis strains and two Raphidiopsis strains were confirmed to produce CYNs. Mutant cyrI and cyrK genes were observed in these strains. Cloned cyr gene sequences from eight water bodies were clustered with cyr genes from Cylindrospermopsis and Raphidiopsis (C/R group) in the phylogenetic trees with high similarities (99%). Four cyrI sequence types and three cyrJ sequence types were observed to have different sequence insertions and repeats. Phylogenetic analysis of the rpoC1 sequences of the C/R group revealed four conserved clades, namely, clade I, clade II, clade III, and clade V. High sequence similarities (>97%) in each clade and a divergent clade IV were observed. Therefore, CYN producers were sporadically distributed in congeneric and paraphyletic C/R group species in Chinese freshwater ecosystems. In the evolution of cyr genes, intragenomic translocations and intergenomic transfer between local Cylindrospermopsis and Raphidiopsis were emphasized and probably mediated by transposases. This research confirms the existence of CYN-producing Cylindrospermopsis in China and reveals the distinctive variations of cyr genes.

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.