Synechococcus: 3 billion years of global dominance

CABO-16S-a Combined Archaea, Bacteria, Organelle 16S rRNA database framework for amplicon analysis of prokaryotes and eukaryotes in environmental samples

Identification of both prokaryotic and eukaryotic microorganisms in environmental samples is currently challenged by the need for additional sequencing to obtain separate 16S and 18S ribosomal RNA (rRNA) amplicons or the constraints imposed by "universal" primers. Organellar 16S rRNA sequences are amplified and sequenced along with prokaryote 16S rRNA and provide an alternative method to identify eukaryotic microorganisms. CABO-16S combines bacterial and archaeal sequences from the SILVA database with 16S rRNA sequences of plastids and other organelles from the PR2 database to enable identification of all 16S rRNA sequences. Comparison of CABO-16S with SILVA 138.2 results in equivalent taxonomic classification of mock communities and increased classification of diverse environmental samples. In particular, identification of phototrophic eukaryotes in shallow seagrass environments, marine waters, and lake waters was increased. The CABO-16S framework allows users to add custom sequences for further classification of underrepresented clades and can be easily updated with future releases of reference databases. Addition of sequences obtained from Sanger sequencing of methane seep sediments and curated sequences of the polyphyletic SEEP-SRB1 clade resulted in differentiation of syntrophic and non-syntrophic SEEP-SRB1 in hydrothermal vent sediments. CABO-16S highlights the benefit of combining and amending existing training sets when studying microorganisms in diverse environments.

RNA-Binding S1 Domain in Bacterial, Archaeal and Eukaryotic Proteins as One of the Evolutionary Markers of Symbiogenesis

The RNA-binding S1 domain is a β-barrel with a highly conserved RNA-binding site on its surface. This domain is an important part of the structures of different bacterial, archaeal, and eukaryotic proteins. A distinctive feature of the S1 domain is multiple presences (structural repeats) in proteins and protein complexes. Here, we have analyzed all available protein sequences in the UniProt database to obtain data on the distribution of bacterial, eukaryotic and archaeal proteins containing the S1 domain. Mainly, the S1 domain is found in bacterial proteins with the number of domains varying from one to eight. Eukaryotic proteins contain from one to fifteen S1 domains, while in archaeal proteins, only one S1 domain is identified. Analysis of eukaryotic proteins containing S1 domains revealed a group of chloroplast S1 ribosomal proteins (ChRpS1) with characteristic properties of bacterial S1 ribosomal proteins (RpS1) from the Cyanobacteria. Also, in a separate group, chloroplast and mitochondrial elongation factor Ts containing two S1 structural domains were assigned. For mitochondrial elongation factor Ts, the features of S1 in comparison with the RpS1 from Cyanobacteria phylum and the Alphaproteobacteria class were revealed. The data obtained allow us to consider the S1 domain as one of the evolutionary markers of the symbiogenesis of bacterial and eukaryotic organisms.

Novel oceanic cyanobacterium isolated from Bangaram island with profound acid neutralizing ability is proposed as Leptolyngbya iicbica sp. nov. strain LK

An acid-neutralizing, filamentous, non-heterocytous, marine cyanobacterium named 'LK' has been isolated from the seashore of Bangaram Island, an atoll of Lakshadweep, India, and is described here as a novel species. LK has been characterized using morphological, ecological, and genomic features. Based on 16S rRNA, whole-genome sequencing, and marker gene-based analysis, LK has been identified as a new species. LK clustered with Leptolyngbya-like strains belonging to the LPP group but diverged from Leptolyngbya sensu stricto, indicating the polyphyletic nature of the Leptolyngbya genus. Leptolyngbya sp. SIOISBB and Halomicronema sp. CCY15110 were identified as LK's two closest phylogenetic neighbors in various phylogenetic studies. The analysis of 16S rRNA, ITS secondary structures, and genome relatedness indices such as AAI, ANI, and gANI strongly support LK as a novel species of the Leptolyngbya genus. The mechanism behind acid neutralization in LK has been delineated, attributing it to a surface phenomenon most likely due to the presence of salts of calcium, magnesium, sodium, and potassium. We name LK as Leptolyngbya iicbica strain LK which is a novel species with prominent acidic pH-neutralizing properties.

Population genomics resolves cryptic species of the ecologically flexible genus Laspinema (Cyanobacteria)

Cyanobacterial taxonomy is entering the genomic era, but only a few taxonomic studies have employed population genomics, which provides a framework and a multitude of tools to understand species boundaries. Phylogenomic and population genomic analyses previously suggested that several cryptic lineages emerged within the genus Laspinema. Here, we apply population genomics to define boundaries between these lineages and propose two new cryptic species, Laspinema olomoucense and L. palackyanum. Moreover, we sampled soil and puddles across Central Europe and sequenced the 16S rRNA gene and 16S-23S ITS region of the isolated Laspinema strains. Together with database mining of 16S rRNA gene sequences, we determined that the genus Laspinema has a cosmopolitan distribution and inhabits a wide variety of habitats, including freshwater, saline water, mangroves, soil crusts, soils, puddles, and the human body.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Nested PCR Approach for petB Gene Metabarcoding of Marine Synechococcus Populations

The molecular diversity of marine picocyanobacterial populations, an important component of phytoplankton communities, is better characterized using high-resolution marker genes than the 16S rRNA gene as they have greater sequence divergence to differentiate between closely related picocyanobacteria groups. Although specific ribosomal primers have been developed, another general disadvantage of bacterial ribosome-based diversity analyses is the variable number of rRNA gene copies. To overcome these issues, the single-copy petB gene, encoding the cytochrome b6 subunit of the cytochrome b6f complex, has been used as a high-resolution marker gene to characterize Synechococcus diversity. We have designed new primers targeting the petB gene and proposed a nested PCR method (termed Ong_2022) for metabarcoding of marine Synechococcus populations obtained by flow cytometry cell sorting. We evaluated the specificity and sensitivity of Ong_2022 against the standard amplification protocol (termed Mazard_2012) using filtered seawater samples. The Ong_2022 approach was also tested on flow cytometry-sorted Synechococcus populations. Samples (filtered and sorted) were obtained in the Southwest Pacific Ocean, from subtropical (ST) and subantarctic (SA) water masses. The two PCR approaches using filtered samples recovered the same dominant subclades, Ia, Ib, IVa, and IVb, with small differences in relative abundance across the distinct samples. For example, subclade IVa was dominant in ST samples with the Mazard_2012 approach, while the same samples processed with Ong_2022 showed similar contributions of subclades IVa and Ib to the total community. The Ong_2022 approach generally captured a higher genetic diversity of Synechococcus subcluster 5.1 than the Mazard_2012 approach while having a lower proportion of incorrectly assigned amplicon sequence variants (ASVs). All flow cytometry-sorted Synechococcus samples could be amplified only by our nested approach. The taxonomic diversity obtained with our primers on both sample types was in agreement with the clade distribution observed by previous studies that applied other marker genes or PCR-free metagenomic approaches under similar environmental conditions. IMPORTANCE The petB gene has been proposed as a high-resolution marker gene to access the diversity of marine Synechococcus populations. A systematic metabarcoding approach based on the petB gene would improve the characterization/assessment of the Synechococcus community structure in marine planktonic ecosystems. We have designed and tested specific primers to be applied in a nested PCR protocol (Ong_2022) for metabarcoding the petB gene. The Ong_2022 protocol can be applied to samples with low DNA content, such as those obtained by flow cytometry cell sorting, allowing the simultaneous assessment of the genetic diversity of Synechococcus populations and cellular properties and activities (e.g., nutrient cell ratios or carbon uptake rates). Our approach will allow future studies using flow cytometry to investigate the link between ecological traits and taxonomic diversity of marine Synechococcus.

Cartusia hunanesis sp. nov. (Oculatellaceae, Oculatellales) from a Stream in China Based on Polyphasic Approach

Cartusia hunanesis sp. nov. was isolated from a stream in China, and two strains (ZJJ02 and ZJJ03) of which were inquired using morphological features, ecological evidence, and molecular data consisting of the 16S rRNA gene and 16S–23S rRNA intergenic transcribed spacer (ITS) region. Cartusia hunanesis varies from the type species Cartusia fontana by having only a single trichome in the sheath and large granules near the cross wall. The investigated strains of C. hunanesis were revealed to be a sister clade of C. fontana, according to the phylogenetic analysis based on the 16S rRNA gene. In addition, the Cartusia was clustered with the family Oculatellaceae members, genera Pegethrix and Elainella. These two strains displayed 97.6% similarity to C. fontana. The ITS region of C. hunanesis was found to be considerably distinct from that of C. fontana in terms of the secondary structure, which demonstrated that C. hunanesis is a novel species owing to the divergences in its morphological and genetic data compared with the related C. fontana.

An updated classification of cyanobacterial orders and families based on phylogenomic and polyphasic analysis

Cyanobacterial taxonomy is facing a period of rapid changes thanks to the ease of 16S rRNA gene sequencing and established workflows for description of new taxa. Since the last comprehensive review of the cyanobacterial system in 2014 until 2021, at least 273 species in 140 genera were newly described. These taxa were mainly placed into previously defined orders and families although several new families were proposed. However, the classification of most taxa still relied on hierarchical relationships inherited from the classical morphological taxonomy. Similarly, the obviously polyphyletic orders such as Synechococcales and Oscillatoriales were left unchanged. In this study, the rising number of genomic sequences of cyanobacteria and well-described reference strains allowed us to reconstruct a robust phylogenomic tree for taxonomic purposes. A less robust but better sampled 16S rRNA gene phylogeny was mapped to the phylogenomic backbone. Based on both these phylogenies, a polyphasic classification throughout the whole phylum of Cyanobacteria was created, with ten new orders and fifteen new families. The proposed system of cyanobacterial orders and families relied on a phylogenomic tree but still employed phenotypic apomorphies where possible to make it useful for professionals in the field. It was, however, confirmed that morphological convergence of phylogenetically distant taxa was a frequent phenomenon in cyanobacteria. Moreover, the limited phylogenetic informativeness of the 16S rRNA gene, resulting in ambiguous phylogenies above the genus level, emphasized the integration of genomic data as a prerequisite for the conclusive taxonomic placement of a vast number of cyanobacterial genera in the future.

Inherent tendency of Synechococcus and heterotrophic bacteria for mutualism on long-term coexistence despite environmental interference

Mutualism between Synechococcus and heterotrophic bacteria has been found to support their prolonged survival in nutrient-depleted conditions. However, environmental interference on the fate of their mutualism is not understood. Here, we show that exogenous nutrients disrupt their established mutualism. Once the exogenous nutrients were exhausted, Synechococcus and heterotrophic bacteria gradually reestablished their metabolic mutualism during 450 days of culture, which revived unhealthy Synechococcus cells. Using metagenomics, metatranscriptomics, and the ¹⁵N tracer method, we reveal that the associated bacterial nitrogen fixation triggered the reestablishment of the mutualism and revival of Synechococcus health. During this process, bacterial community structure and functions underwent tremendous adjustments to achieve the driving effect, and a cogeneration of nitrogen, phosphorus, iron, and vitamin by the heterotrophic bacteria sustained Synechococcus's prolonged healthy growth. Our findings suggest that Synechococcus and heterotrophic bacteria may have an inherent tendency toward mutualism despite environmental interference. This may exhibit their coevolutionary adaptations in nutrient-deficient environments.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

New cyanobacterial genus Argonema is hidding in soil crusts around the world

Cyanobacteria are crucial primary producers in soil and soil crusts. However, their biodiversity in these habitats remains poorly understood, especially in the tropical and polar regions. We employed whole genome sequencing, morphology, and ecology to describe a novel cyanobacterial genus Argonema isolated from Antarctica. Extreme environments are renowned for their relatively high number of endemic species, but whether cyanobacteria are endemic or not is open to much current debate. To determine if a cyanobacterial lineage is endemic is a time consuming, elaborate, and expensive global sampling effort. Thus, we propose an approach that will help to overcome the limits of the sampling effort and better understand the global distribution of cyanobacterial clades. We employed a Sequencing Read Archive, which provides a rich source of data from thousands of environmental samples. We developed a framework for a characterization of the global distribution of any microbial species using Sequencing Read Archive. Using this approach, we found that Argonema is actually cosmopolitan in arid regions. It provides further evidence that endemic microbial taxa are likely to be much rarer than expected.

Reevaluation of Parasynechococcus-like Strains and Genomic Analysis of Their Microsatellites and Compound Microsatellites

Morphologically similar to Synechococcus, a large number of Parasynechococcus strains were misclassified, resulting in extreme underestimation of their genetic diversity. In this study, 80 Synechococcus-like strains were reevaluated using a combination of 16S rRNA phylogeny and genomic approach, identifying 54 strains as Parasynechococcus-like strains and showing considerably intragenus genetic divergence among the subclades identified. Further, bioinformatics analysis disclosed diversified patterns of distribution, abundance, density, and diversity of microsatellites (SSRs) and compound microsatellites (CSSRs) in genomes of these Parasynechococcus-like strains. Variations of SSRs and CSSRs were observed amongst phylotypes and subclades. Both SSRs and CSSRs were in particular unequally distributed among genomes. Dinucleotide SSRs were the most widespread, while the genomes showed two patterns in the second most abundant repeat type (mononucleotide or trinucleotide SSRs). Both SSRs and CSSRs were predominantly observed in coding regions. These two types of microsatellites showed positive correlation with genome size (p < 0.01) but negative correlation with GC content (p < 0.05). Additionally, the motif (A)n, (AG)n and (AGC)n was a major one in the corresponding category. Meanwhile, distinctive motifs of CSSRs were found in 39 genomes. This study characterizes SSRs and CSSRs in genomes of Parasynechococcus-like strains and will be useful as a prerequisite for future studies regarding their distribution, function, and evolution. Moreover, the identified SSRs may facilitate fast acclimation of Parasynechococcus-like strains to fluctuating environments and contribute to the extensive distribution of Parasynechococcus species in global marine environments.

The "Dark Side" of Picocyanobacteria: Life as We Do Not Know It (Yet)

Picocyanobacteria of the genus Synechococcus (together with Cyanobium and Prochlorococcus) have captured the attention of microbial ecologists since their description in the 1970s. These pico-sized microorganisms are ubiquitous in aquatic environments and are known to be some of the most ancient and adaptable primary producers. Yet, it was only recently, and thanks to developments in molecular biology and in the understanding of gene sequences and genomes, that we could shed light on the depth of the connection between their evolution and the history of life on the planet. Here, we briefly review the current understanding of these small prokaryotic cells, from their physiological features to their role and dynamics in different aquatic environments, focussing particularly on the still poorly understood ability of picocyanobacteria to adapt to dark conditions. While the recent discovery of Synechococcus strains able to survive in the deep Black Sea highlights how adaptable picocyanobacteria can be, it also raises more questions-showing how much we still do not know about microbial life. Using available information from brackish Black Sea strains able to perform and survive in dark (anoxic) conditions, we illustrate how adaptation to narrow ecological niches interacts with gene evolution and metabolic capacity.

Photoautotrophic poly(3-hydroxybutyrate) production by a wild-type Synechococcus elongatus isolated from an extreme environment

The photoautotrophic poly(3-hydroxybutyrate) (PHB) production by cyanobacteria is an attractive option as it only requires CO₂ and light. In this work, a new wild-type strain producing PHB, Synechococcus elongatus UAM-C/S03, was identified using a polyphasic approach. The strain was cultured in a photobioreactor operated under N-sufficiency conditions at different pH values (7 to 11) and fed with CO₂ on demand. We also evaluated the production of PHB under N-starving conditions. Highest biomass productivity, 324 mg L^-1 d^-1, and CO₂ capture, 674 mg L^-1 d^-1, were obtained at pH 7 and under N-sufficiency conditions. The strain accumulated 29.42% of PHB in dry cell weight (DCW) under N-starvation conditions without pH control, and highest PHB productivity was 58.10 mg L^-1 d^-1. The highest carbohydrate content registered at pH 8, 50.84% in DCW, along with a release of carbon-based organic compounds, suggested the presence of exopolysaccharides in the culture medium.

Tenebriella gen. nov. - the dark twin of Oscillatoria

Oscillatoria has long been known to be polyphyletic. After recent resequencing of the reference strain for this genus, many Oscillatoria-like groups phylogenetically distant from the type species O. princeps remained unresolved. Here we describe one of these groups as a new genus Tenebriella. Most of the studied strains originate from Central Europe, where they are able to form prominent microbial mats. Despite the overall Oscillatoria-like morphology, Tenebriella can be distinguished by darker trichomes and forms a separate monophyletic clade in phylogenies inferred from the 16S rRNA gene and two additional loci (rpoC1, rbcLX). Within Tenebriella we recognize two new species differing from each other by morphological and ecological characteristics. First species does not fit any known taxon description, and thus is described as a new species T. amphibia. The latter one corresponds with the information available for Oscillatoria curviceps Agardh ex Gomont, and thus new combination T. curviceps is proposed. The phylogenetic analyses of the 16S-23S ITS region together with the comparison of the hypothetical secondary structures confirmed recognition of these two species and additionally revealed presence of a morphologically cryptic species Tenebriella sp. The results corroborate frequent recurrence of convergent morphotypes in the evolution of cyanobacteria and justify further exploration even of the intensively studied European freshwaters using molecular phylogenetics to discover new and ecologically relevant taxa.

Novel cyanotoxin-producing Synechococcus in tropical lakes

Picocyanobacteria are small cyanobacteria, being about 0.8-1.5 µm in size. They are present in freshwater environments all over the world and are known to cause harmful algal blooms, although their effects are not well understood. Algal blooms are important to manage because they threaten freshwater resources, with potentially severe effects on ecological and human health. There is also increased urgency due to urbanization and climate change trends which are expected to exacerbate these bloom dynamics. These changes are expected to especially favour picocyanobacteria groups, emphasizing the need for better characterization of their effects in the environment. In this study, we report the discovery that Synechococcus sp. could produce cylindrospermopsin (CYN) and anatoxin-a (ATX). This ability had never been previously reported for this species. Their toxin genes were also partial compared to other major producers such as Raphidiopsis sp. and Anabaena sp., demonstrating potentially unique synthesis pathways that provides insight into the various mechanisms of genetic variation that drives toxin synthesis. The Synechococcus sp. strains were found to produce about 9.0 × 10^-5-6.8 × 10^-4 fg CYN cell^-1 and 4.7 × 10^-4-1.5 × 10^-2 fg ATX cell^-1. The potential for Synechococcus sp. to be toxic highlights a global concern due to its widespread distribution, and through environmental trends that increasingly favour its productivity within freshwater systems around the world.

Heterozygous, Polyploid, Giant Bacterium, Achromatium, Possesses an Identical Functional Inventory Worldwide across Drastically Different Ecosystems

Achromatium is large, hyperpolyploid and the only known heterozygous bacterium. Single cells contain approximately 300 different chromosomes with allelic diversity far exceeding that typically harbored by single bacteria genera. Surveying all publicly available sediment sequence archives, we show that Achromatium is common worldwide, spanning temperature, salinity, pH, and depth ranges normally resulting in bacterial speciation. Although saline and freshwater Achromatium spp. appear phylogenetically separated, the genus Achromatium contains a globally identical, complete functional inventory regardless of habitat. Achromatium spp. cells from differing ecosystems (e.g., from freshwater to saline) are, unexpectedly, equally functionally equipped but differ in gene expression patterns by transcribing only relevant genes. We suggest that environmental adaptation occurs by increasing the copy number of relevant genes across the cell's hundreds of chromosomes, without losing irrelevant ones, thus maintaining the ability to survive in any ecosystem type. The functional versatility of Achromatium and its genomic features reveal alternative genetic and evolutionary mechanisms, expanding our understanding of the role and evolution of polyploidy in bacteria while challenging the bacterial species concept and drivers of bacterial speciation.

Sedimentary DNA record of eukaryotic algal and cyanobacterial communities in a shallow Lake driven by human activities and climate change

Global freshwater lakes are changing due to human activities and climate change. Unfortunately, sufficient long-term monitoring is lacking for most lakes. However, lake sedimentary archives can extend the instrumental record and reveal historical environmental trends. In particular, sedimentary DNA analysis of lacustrine sediment cores can aid the reconstruction of past trends in eukaryotic algal and cyanobacterial communities, as was conducted in this study for Lake Chaohu in China. The results presented here indicate that the construction of the Chaohu Dam in 1963 is associated with decreased richness of eukaryotic algal and cyanobacterial communities. Several groups, including the eukaryotic algal taxa, Chlorophyceae, and cyanobacterial groups like Dolichospermum, Microcystis, Planktothricoides, Cyanobium, Pseudanabaena, and Synechococcus, increased in abundance following inferred historical nutrient enrichment. Nutrient concentrations and hydrologic conditions were further implicated as the dominant controls on communities based on Random Forest and generalized additive modeling statistical analyses. In particular, significant increases in lake hydraulic residence times after the construction of the Chaohu Dam were significantly associated with altered biological community structures. Further, phosphorus enrichment was positively associated with increased richness and diversity of these communities following the 1980s. In addition, effects from increased atmospheric temperatures on eukaryotic algal and cyanobacterial communities were apparent. Here, high-throughput sequencing analysis of sedimentary DNA allowed the inference of long-term biodiversity dynamics of Lake Chaohu. These results underscore the important impacts of anthropogenic activities and climate change on aquatic ecosystems at the decadal scale.

Bacterial community composition and diversity in the ballast water of container ships arriving at Yangshan Port, Shanghai, China

Ballast water is a major vector of invasion by protozoans and metazoans. Bacterial invasion is less-well understood. We surveyed the bacterial diversity of ballast water from 26 container ships arriving at the Yangshan Deepwater Port, Shanghai, China during 2015-2016. We characterized the ballast microbiome using high-throughput sequencing (HTS) based on V4-V5 region of 16S rRNA genes. We simultaneously monitored physicochemical parameters of the ballast water, including temperature, pH, dissolved oxygen (DO), salinity, turbidity, total suspended solid (TSS), particulate organic carbon (POC), NO₂, NH₄, PO₄. Proteobacteria was the dominant phylum, comprising more than 50% of the OTUs of almost all vessels, followed by Bacteroidetes (12.08%), Actinobacteria (4.86%) Planctomycetes (3.24%) and Cyanobacteria (1.95%). The relative abundance of Cyanobacteria differed among vessels. It was negatively correlated with temperature, NO₃, pH, TSS, PO₄, and turbidity and positively correlated with NH₄, POC. The genus Synechococcus was the most common Cyanobacteria in our results. Escherichia coli were relatively rare; they are indicator-species of D-2 standards published by the IMO. The relative abundance of the genus Vibrio ranged from 0.003% to 24.88% among different vessels. Our results showed that HTS was able to profile the bacterial communities in ballast-waters, even when the approach was restricted by technical and other obstacles.

A new genomic taxonomy system for the Synechococcus collective

Cyanobacteria of the genus Synechococcus are major contributors to global primary productivity and are found in a wide range of aquatic ecosystems. This Synechococcus collective (SC) is metabolically diverse, with some lineages thriving in polar and nutrient-rich locations and others in tropical or riverine waters. Although many studies have discussed the ecology and evolution of the SC, there is a paucity of knowledge on its taxonomic structure. Thus, we present a new taxonomic classification framework for the SC based on recent advances in microbial genomic taxonomy. Phylogenomic analyses of 1085 cyanobacterial genomes demonstrate that organisms classified as Synechococcus are polyphyletic at the order rank. The SC is classified into 15 genera, which are placed into five distinct orders within the phylum Cyanobacteria: (i) Synechococcales (Cyanobium, Inmanicoccus, Lacustricoccus gen. Nov., Parasynechococcus, Pseudosynechococcus, Regnicoccus, Synechospongium gen. nov., Synechococcus and Vulcanococcus); (ii) Cyanobacteriales (Limnothrix); (iii) Leptococcales (Brevicoccus and Leptococcus); (iv) Thermosynechococcales (Stenotopis and Thermosynechococcus) and (v) Neosynechococcales (Neosynechococcus). The newly proposed classification is consistent with habitat distribution patterns (seawater, freshwater, brackish and thermal environments) and reflects the ecological and evolutionary relationships of the SC.

Presence of toxin-antitoxin systems in picocyanobacteria and their ecological implications

Picocyanobacteria (mainly Synechococcus and Prochlorococcus) contribute significantly to ocean's primary production. Toxin-Antitoxin (TA) systems present in bacteria and archaea are known to regulate cell growth in response to environmental stresses. However, little is known about the presence of TA systems in picocyanobacteria. This study investigated complete genomes of Synechococcus and Prochlorococcus to understand the prevalence of TA systems in picocyanobacteria. Using the TAfinder software, Type II TA systems were predicted in 27 of 33 (81%) Synechococcus strains, but none of 38 Prochlorococcus strains contain TA genes. Synechococcus strains with larger genomes tend to contain more putative type II TA systems. The number of TA pairs varies from 0 to 42 in Synechococcus strains isolated from various environments. A linear correlation between the genome size and the number of putative TA systems in both coastal and freshwater Synechococcus was established. In general, open ocean Synechococcus contain no or few TA systems, while coastal and freshwater Synechococcus contain more TA systems. The type II TA systems inhibit microbial translation via ribonucleases and allow cells to enter the "dormant" stage in adverse environments. Inheritance of TA genes in freshwater and coastal Synechococcus could confer a recoverable persister mechanism important to survive in variable environments.

Spatial heterogeneity and hydrological fluctuations drive bacterioplankton community composition in an Amazon floodplain system

Amazonian floodplains form complex hydrological networks that play relevant roles in global biogeochemical cycles, and bacterial degradation of the organic matter in these systems is key for regional carbon budget. The Amazon undergoes extreme seasonal variations in water level, which produces changes in landscape and diversifies sources of organic inputs into floodplain systems. Although these changes should affect bacterioplankton community composition (BCC), little is known about which factors drive spatial and temporal patterns of bacterioplankton in these Amazonian floodplains. We used high-throughput sequencing (Illumina MiSeq) of the V3-V4 region of the 16S rRNA gene to investigate spatial and temporal patterns of BCC of two size fractions, and their correlation with environmental variables in an Amazon floodplain lake (Lago Grande do Curuai). We found a high degree of novelty in bacterioplankton, as more than half of operational taxonomic units (OTUs) could not be classified at genus level. Spatial habitat heterogeneity and the flood pulse were the main factors shaping free-living (FL) BCC. The gradient of organic matter from transition zone-lake-Amazon River was the main driver for particle-attached (PA) BCC. The BCC reflected the complexity of the system, with more variation in space than in time, although both factors were important drivers of the BCC in this Amazon floodplain system.

Taxonomic resolution of the genus Cyanothece (Chroococcales, Cyanobacteria), with a treatment on Gloeothece and three new genera, Crocosphaera, Rippkaea, and Zehria

The systematics of single-celled cyanobacteria represents a major challenge due to morphological convergence and application of various taxonomic concepts. The genus Cyanothece is one of the most problematic cases, as the name has been applied to oval-shaped coccoid cyanobacteria lacking sheaths with little regard to their phylogenetic position and details of morphology and ultrastructure. Hereby we analyze an extensive set of complementary genetic and phenotypic evidence to disentangle the relationships among these cyanobacteria. We provide diagnostic characters to separate the known genera Cyanothece, Gloeothece, and Aphanothece, and provide a valid description for Crocosphaera gen. nov. We describe two new genera, Rippkaea and Zehria, to characterize two distinct phylogenetic lineages outside the previously known genera. We further describe 13 new species in total including Cyanothece svehlovae, Gloeothece aequatorialis, G. aurea, G. bryophila, G. citriformis, G. reniformis, Gloeothece tonkinensis, G. verrucosa, Crocosphaera watsonii, C. subtropica, C. chwakensis, Rippkaea orientalis, and Zehria floridana to recognize the intrageneric diversity as rendered by polyphasic analysis. We discuss the close relationship of free-living cyanobacteria from the Crocosphaera lineage to nitrogen-fixing endosymbionts of marine algae. The current study includes several experimental strains (Crocosphaera and "Cyanothece") important for the study of diazotrophy and the global oceanic nitrogen cycle, and provides evidence suggesting ancestral N₂ -fixing capability in the chroococcalean lineage.

Primer Design for an Accurate View of Picocyanobacterial Community Structure by Using High-Throughput Sequencing

High-throughput sequencing (HTS) of the 16S rRNA gene has been used successfully to describe the structure and dynamics of microbial communities. Picocyanobacteria are important members of bacterioplankton communities, and, so far, they have predominantly been targeted using universal bacterial primers, providing a limited resolution of the picocyanobacterial community structure and dynamics. To increase such resolution, the study of a particular target group is best approached with the use of specific primers. Here, we aimed to design and evaluate specific primers for aquatic picocyanobacterial genera to be used with high-throughput sequencing. Since the various regions of the 16S rRNA gene have different degrees of conservation in different bacterial groups, we therefore first determined which hypervariable region of the 16S rRNA gene provides the highest taxonomic and phylogenetic resolution for the genera Synechococcus, Prochlorococcus, and Cyanobium An in silico analysis showed that the V5, V6, and V7 hypervariable regions appear to be the most informative for this group. We then designed primers flanking these hypervariable regions and tested them in natural marine and freshwater communities. We successfully detected that most (97%) of the obtained reads could be assigned to picocyanobacterial genera. We defined operational taxonomic units as exact sequence variants (zero-radius operational taxonomic units [zOTUs]), which allowed us to detect higher genetic diversity and infer ecologically relevant information about picocyanobacterial community composition and dynamics in different aquatic systems. Our results open the door to future studies investigating picocyanobacterial diversity in aquatic systems.IMPORTANCE The molecular diversity of the aquatic picocyanobacterial community cannot be accurately described using only the available universal 16S rRNA gene primers that target the whole bacterial and archaeal community. We show that the hypervariable regions V5, V6, and V7 of the 16S rRNA gene are better suited to study the diversity, community structure, and dynamics of picocyanobacterial communities at a fine scale using Illumina MiSeq sequencing. Due to its variability, it allows reconstructing phylogenies featuring topologies comparable to those generated when using the complete 16S rRNA gene sequence. Further, we successfully designed a new set of primers flanking the V5 to V7 region whose specificity for picocyanobacterial genera was tested in silico and validated in several freshwater and marine aquatic communities. This work represents a step forward for understanding the diversity and ecology of aquatic picocyanobacteria and sets the path for future studies on picocyanobacterial diversity.

Evolutionary Patterns of Thylakoid Architecture in Cyanobacteria

While photosynthetic processes have become increasingly understood in cyanobacterial model strains, differences in the spatial distribution of thylakoid membranes among various lineages have been largely unexplored. Cyanobacterial cells exhibit an intriguing diversity in thylakoid arrangements, ranging from simple parietal to radial, coiled, parallel, and special types. Although metabolic background of their variability remains unknown, it has been suggested that thylakoid patterns are stable in certain phylogenetic clades. For decades, thylakoid arrangements have been used in cyanobacterial classification as one of the crucial characters for definition of taxa. The last comprehensive study addressing their evolutionary history in cyanobacteria was published 15 years ago. Since then both DNA sequence and electron microscopy data have grown rapidly. In the current study, we map ultrastructural data of >200 strains onto the SSU rRNA gene tree, and the resulting phylogeny is compared to a phylogenomic tree. Changes in thylakoid architecture in general follow the phylogeny of housekeeping loci. Parietal arrangement is resolved as the original thylakoid organization, evolving into complex arrangement in the most derived group of heterocytous cyanobacteria. Cyanobacteria occupying intermediate phylogenetic positions (greater filamentous, coccoid, and baeocytous types) exhibit fascicular, radial, and parallel arrangements, partly tracing the reconstructed course of phylogenetic branching. Contrary to previous studies, taxonomic value of thylakoid morphology seems very limited. Only special cases such as thylakoid absence or the parallel arrangement could be used as taxonomically informative apomorphies. The phylogenetic trees provide evidence of both paraphyly and reversion from more derived architectures in the simple parietal thylakoid pattern. Repeated convergent evolution is suggested for the radial and fascicular architectures. Moreover, thylakoid arrangement is constrained by cell size, excluding the occurrence of complex architectures in cyanobacteria smaller than 2 μm in width. It may further be dependent on unknown (eco)physiological factors as suggested by recurrence of the radial type in unrelated but morphologically similar cyanobacteria, and occurrence of special features throughout the phylogeny. No straightforward phylogenetic congruences have been found between proteins involved in photosynthesis and thylakoid formation, and the thylakoid patterns. Remarkably, several postulated thylakoid biogenesis factors are partly or completely missing in cyanobacteria, challenging their proposed essential roles.

A singular nitric oxide synthase with a globin domain found in Synechococcus PCC 7335 mobilizes N from arginine to nitrate

The enzyme nitric oxide synthase (NOS) oxidizes L-arginine to NO and citrulline. In this work, we characterise the NOS from the cyanobacteria Synechococcus PCC 7335 (SyNOS). SyNOS possesses a canonical mammalian NOS architecture consisting of oxygenase and reductase domains. In addition, SyNOS possesses an unusual globin domain at the N-terminus. Recombinant SyNOS expressed in bacteria is active, and its activity is suppressed by the NOS inhibitor L-NAME. SyNOS allows E. coli to grow in minimum media containing L-arginine as the sole N source, and has a higher growth rate during N deficiency. SyNOS is expressed in Synechococcus PCC 7335 where NO generation is dependent on L-arginine concentration. The growth of Synechococcus is dramatically inhibited by L-NAME, suggesting that SyNOS is essential for this cyanobacterium. Addition of arginine in Synechococcus increases the phycoerythrin content, an N reservoir. The role of the novel globin domain in SyNOS is discussed as an evolutionary advantage, conferring new functional capabilities for N metabolism.

Metagenomic Analysis of Cyanobacteria in an Oligotrophic Tropical Estuary, South Atlantic

This study assessed the species composition, distribution, and functional profiles of cyanobacteria in Camamu Bay, a tropical oligotrophic estuarine system on the northeast coast of Brazil, using shotgun metagenomic sequencing. Surface-water samples were evaluated in two different rainfall periods (rainy and dry seasons), at nine stations in the three hydrodynamic regions of the bay. At a fixed sampling station, on each season, samples were taken over a tidal cycle at 3-h intervals over 12 h. A total of 219 cyanobacterial taxa were identified, demonstrating a diverse community of freshwater, euryhaline, and marine cyanobacteria. The genera of greater relative abundance, Synechococcus and Prochlorococcus, corresponded to the picoplankton fraction. Although Camamu Bay has conspicuous marine characteristics, the contribution of freshwater during the rainy season caused variation in cyanobacteria community, with an increase in species richness. Due the high prevalence of Synechococcus (90% of the sequences), the functional analysis revealed only minor differences in gene content between the dry and rainy seasons. In both rainy and dry seasons, an increase in Prochlorococcus relative abundance occurred during high tide, demonstrating the tidal influence in the bay. The environmental characteristics of the bay provide niche conditions for a wide variety of cyanobacteria, including freshwater, euryhaline, and marine strains.

Seasonal variability of picophytoplankton under contrasting environments in northern Tunisian coasts, southwestern Mediterranean Sea

We investigated at the single cell level during 16months (June 2012 to September 2013) the temporal distribution of picophytoplankton (picoeukaryotes, Synechococcus and Prochlorococcus) communities in two contrasted ecosystems: the Bay of Bizerte characterised by an oligotrophic regime typical of the Mediterranean Sea and the Bizerte Lagoon that exhibits a mesotrophic/eutrophic state. We aimed at depicting seasonal variations and quantifying the relationships between the environmental factors and the structure and abundance of picophytoplankton communities. Results showed that picophytoplankton groups were able to grow under a wide range of environmental conditions varying seasonally, although their abundances and contributions to the total chlorophyll biomass significantly varied and showed importance in the Bay of Bizerte. Synechococcus was the most abundant group reaching 225∗10³cells·cm^-3 in the Bay and 278∗10³cells·cm^-3 in the lagoon. This group was present all over the year in both ecosystems. Structural equation model results pointed out a different configuration regarding the picophytoplankton environmental drivers. The complexity of the configuration, i.e. number of significant links within the system, decreased under enhanced eutrophication conditions. The less exposure to anthropogenic stress, i.e. in the Bay of Bizerte, highlight a larger role of nutrient and hydrological conditions on the seasonal variations of picophytoplankton, whereas a negative effect of eutrophication on picophytoplankton communities was unveiled in the Bizerte Lagoon. We stress that such influence may be exacerbated under expected scenarios of Mediterranean warming conditions and nutrient release in coastal ecosystems.

Host Traits Drive Viral Life Histories across Phytoplankton Viruses

Viruses are integral to ecological and evolutionary processes, but we have a poor understanding of what drives variation in key traits across diverse viruses. For lytic viruses, burst size, latent period, and genome size are primary characteristics controlling host-virus dynamics. Here we synthesize data on these traits for 75 strains of phytoplankton viruses, which play an important role in global biogeochemistry. We find that primary traits of the host (genome size, growth rate) explain 40%-50% of variation in burst size and latent period. Specifically, burst size and latent period both exhibit saturating relationships versus the host∶virus genome size ratio, with both traits increasing at low genome size ratios while showing no relationship at high size ratios. In addition, latent period declines as host growth rate increases. We analyze a model of latent period evolution to explore mechanisms that could cause these patterns. The model predicts that burst size may often be set by the host genomic resources available for viral construction, while latent period evolves to permit this maximal burst size, modulated by host metabolic rate. These results suggest that general mechanisms may underlie the evolution of diverse viruses. Future extensions of this work could help explain viral regulation of host populations, viral influence on community structure and diversity, and viral roles in biogeochemical cycles.

Dating phototrophic microbial lineages with reticulate gene histories

Phototrophic bacteria are among the most biogeochemically significant organisms on Earth and are physiologically related through the use of reaction centers to collect photons for energy metabolism. However, the major phototrophic lineages are not closely related to one another in bacterial phylogeny, and the origins of their respective photosynthetic machinery remain obscured by time and low sequence similarity. To better understand the co-evolution of Cyanobacteria and other ancient anoxygenic phototrophic lineages with respect to geologic time, we designed and implemented a variety of molecular clocks that use horizontal gene transfer (HGT) as additional, relative constraints. These HGT constraints improve the precision of phototroph divergence date estimates and indicate that stem green non-sulfur bacteria are likely the oldest phototrophic lineage. Concurrently, crown Cyanobacteria age estimates ranged from 2.2 Ga to 2.7 Ga, with stem Cyanobacteria diverging ~2.8 Ga. These estimates provide a several hundred Ma window for oxygenic photosynthesis to evolve prior to the Great Oxidation Event (GOE) ~2.3 Ga. In all models, crown green sulfur bacteria diversify after the loss of the banded iron formations from the sedimentary record (~1.8 Ga) and may indicate the expansion of the lineage into a new ecological niche following the GOE. Our date estimates also provide a timeline to investigate the temporal feasibility of different photosystem HGT events between phototrophic lineages. Using this approach, we infer that stem Cyanobacteria are unlikely to be the recipient of an HGT of photosystem I proteins from green sulfur bacteria but could still have been either the HGT donor or the recipient of photosystem II proteins with green non-sulfur bacteria, prior to the GOE. Together, these results indicate that HGT-constrained molecular clocks are useful tools for the evaluation of various geological and evolutionary hypotheses, using the evolutionary histories of both genes and organismal lineages.

Biodiversity of the microbial mat of the Garga hot spring

BACKGROUND

Microbial mats are a good model system for ecological and evolutionary analysis of microbial communities. There are more than 20 alkaline hot springs on the banks of the Barguzin river inflows. Water temperature reaches 75 °C and pH is usually 8.0-9.0. The formation of microbial mats is observed in all hot springs. Microbial communities of hot springs of the Baikal rift zone are poorly studied. Garga is the biggest hot spring in this area.

RESULTS

In this study, we investigated bacterial and archaeal diversity of the Garga hot spring (Baikal rift zone, Russia) using 16S rRNA metagenomic sequencing. We studied two types of microbial communities: (i) small white biofilms on rocks in the points with the highest temperature (75 °C) and (ii) continuous thick phototrophic microbial mats observed at temperatures below 70 °C. Archaea (mainly Crenarchaeota; 19.8% of the total sequences) were detected only in the small biofilms. The high abundance of Archaea in the sample from hot springs of the Baikal rift zone supplemented our knowledge of the distribution of Archaea. Most archaeal sequences had low similarity to known Archaea. In the microbial mats, primary products were formed by cyanobacteria of the genus Leptolyngbya. Heterotrophic microorganisms were mostly represented by Actinobacteria and Proteobacteria in all studied samples of the microbial mats. Planctomycetes, Chloroflexi, and Chlorobi were abundant in the middle layer of the microbial mats, while heterotrophic microorganisms represented mostly by Firmicutes (Clostridia, strict anaerobes) dominated in the bottom part. Besides prokaryotes, we detect some species of Algae with help of detection their chloroplasts 16 s rRNA.

CONCLUSIONS

High abundance of Archaea in samples from hot springs of the Baikal rift zone supplemented our knowledge of the distribution of Archaea. Most archaeal sequences had low similarity to known Archaea. Metagenomic analysis of microbial communities of the microbial mat of Garga hot spring showed that the three studied points sampled at 70 °C, 55 °C, and 45 °C had similar species composition. Cyanobacteria of the genus Leptolyngbya dominated in the upper layer of the microbial mat. Chloroflexi and Chlorobi were less abundant and were mostly observed in the middle part of the microbial mat. We detected domains of heterotrophic organisms in high abundance (Proteobacteria, Firmicutes, Verrucomicrobia, Planctomicetes, Bacteroidetes, Actinobacteria, Thermi), according to metabolic properties of known relatives, which can form complete cycles of carbon, sulphur, and nitrogen in the microbial mat. The studied microbial mats evolved in early stages of biosphere formation. They can live autonomously, providing full cycles of substances and preventing live activity products poisoning.

A penalty on photosynthetic growth in fluctuating light

Fluctuating light is the norm for photosynthetic organisms, with a wide range of frequencies (0.00001 to 10 Hz) owing to diurnal cycles, cloud cover, canopy shifting and mixing; with broad implications for climate change, agriculture and bioproduct production. Photosynthetic growth in fluctuating light is generally considered to improve with increasing fluctuation frequency. Here we demonstrate that the regulation of photosynthesis imposes a penalty on growth in fluctuating light for frequencies in the range of 0.01 to 0.1 Hz (organisms studied: Synechococcus elongatus and Chlamydomonas reinhardtii). We provide a comprehensive sweep of frequencies and duty cycles. In addition, we develop a 2^nd order model that identifies the source of the penalty to be the regulation of the Calvin cycle – present at all frequencies but compensated at high frequencies by slow kinetics of RuBisCO.

Diversity in photosynthetic electron transport under [CO2]-limitation: the cyanobacterium Synechococcus sp. PCC 7002 and green alga Chlamydomonas reinhardtii drive an O2-dependent alternative electron flow and non-photochemical quenching of chlorophyll fluorescence during CO2-limited photosynthesis

Some cyanobacteria, but not all, experience an induction of alternative electron flow (AEF) during CO₂-limited photosynthesis. For example, Synechocystis sp. PCC 6803 (S. 6803) exhibits AEF, but Synechococcus elongatus sp. PCC 7942 does not. This difference is due to the presence of flavodiiron 2 and 4 proteins (FLV2/4) in S. 6803, which catalyze electron donation to O₂. In this study, we observed a low-[CO₂] induced AEF in the marine cyanobacterium Synechococcus sp. PCC 7002 that lacks FLV2/4. The AEF shows high affinity for O₂, compared with AEF mediated by FLV2/4 in S. 6803, and can proceed under extreme low [O₂] (about a few µM O₂). Further, the transition from CO₂-saturated to CO₂-limited photosynthesis leads a preferential excitation of PSI to PSII and increased non-photochemical quenching of chlorophyll fluorescence. We found that the model green alga Chlamydomonas reinhardtii also has an O₂-dependent AEF showing the same affinity for O₂ as that in S. 7002. These data represent the diverse molecular mechanisms to drive AEF in cyanobacteria and green algae. In this paper, we further discuss the diversity, the evolution, and the physiological function of strategy to CO₂-limitation in cyanobacterial and green algal photosynthesis.

Probabilistic inference of lateral gene transfer events

Background

Lateral gene transfer (LGT) is an evolutionary process that has an important role in biology. It challenges the traditional binary tree-like evolution of species and is attracting increasing attention of the molecular biologists due to its involvement in antibiotic resistance. A number of attempts have been made to model LGT in the presence of gene duplication and loss, but reliably placing LGT events in the species tree has remained a challenge.

Results

In this paper, we propose probabilistic methods that samples reconciliations of the gene tree with a dated species tree and computes maximum a posteriori probabilities. The MCMC-based method uses the probabilistic model DLTRS, that integrates LGT, gene duplication, gene loss, and sequence evolution under a relaxed molecular clock for substitution rates. We can estimate posterior distributions on gene trees and, in contrast to previous work, the actual placement of potential LGT, which can be used to, e.g., identify “highways” of LGT.

Conclusions

Based on a simulation study, we conclude that the method is able to infer the true LGT events on gene tree and reconcile it to the correct edges on the species tree in most cases. Applied to two biological datasets, containing gene families from Cyanobacteria and Molicutes, we find potential LGTs highways that corroborate other studies as well as previously undetected examples.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1268-2) contains supplementary material, which is available to authorized users.

Sedimentary DNA Reveals Cyanobacterial Community Diversity over 200 Years in Two Perialpine Lakes

We reconstructed cyanobacterial community structure and phylogeny using DNA that was isolated from layers of stratified sediments spanning 200 years of lake history in the perialpine lakes Greifensee and Lake Zurich (Switzerland). Community analysis based on amplification and sequencing of a 400-nucleotide (nt)-long 16S rRNA fragment specific to Cyanobacteria revealed operational taxonomic units (OTUs) capturing the whole phylum, including representatives of a newly characterized clade termed Melainabacteria, which shares common ancestry with Cyanobacteria and has not been previously described in lakes. The reconstruction of cyanobacterial richness and phylogenetic structure was validated using a data set consisting of 40 years of pelagic microscopic counts from each lake. We identified the OTUs assigned to common taxa known to be present in Greifensee and Lake Zurich and found a strong and significant relationship (adjusted R² = 0.89; P < 0.001) between pelagic species richness in water and OTU richness in the sediments. The water-sediment richness relationship varied between cyanobacterial orders, indicating that the richness of Chroococcales and Synechococcales may be underestimated by microscopy. PCR detection of the microcystin synthetase gene mcyA confirmed the presence of potentially toxic cyanobacterial taxa over recent years in Greifensee and throughout the last century in Lake Zurich. The approach presented in this study demonstrates that it is possible to reconstruct past pelagic cyanobacterial communities in lakes where the integrity of the sedimentary archive is well preserved and to explore changes in phylogenetic and functional diversity over decade-to-century timescales.

IMPORTANCE

Cyanobacterial blooms can produce toxins that affect water quality, especially under eutrophic conditions, which are a consequence of human-induced climate warming and increased nutrient availability. Lakes worldwide have suffered from regular cyanobacterial blooms over the last century. The lack of long-term data limits our understanding of how these blooms form. We successfully reconstructed the past diversity of whole cyanobacterial communities over two hundred years by sequencing genes preserved in the sediments of two perialpine lakes in Switzerland. We identified changes in diversity over time and validated our results using existing data collected in the same two lakes over the past 40 years. This work shows the potential of our approach for addressing important ecological questions about the effects of a changing environment on lake ecology.