Novel Synechococcus Genomes Reconstructed from Freshwater Reservoirs

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Genomic insights into the adaptation of Synechococcus to the coastal environment on Xiamen

Synechococcus are widely distributed in the global ocean, from the pelagic zone to coastal waters. However, little is known about Synechococcus in coastal seawater due to limitations in isolation and culture conditions. In this study, a combination of metagenomic sequencing technology, flow cytometry sorting, and multiple displacement amplification was used to investigate Synechococcus in the coastal seawater of Xiamen Island. The results revealed 18 clades of Synechococcus and diverse metabolic genes that appear to contribute to their adaptation to the coastal environment. Intriguingly, some metabolic genes related to the metabolism of carbohydrates, energy, nucleotides, and amino acids were found in 89 prophage regions that were detected in 16,258 Synechococcus sequences. The detected metabolic genes can enable prophages to contribute to the adaptation of Synechococcus hosts to the coastal marine environment. The detection of prophages also suggests that the cyanophages have infected Synechococcus. On the other hand, the identification of 773 genes associated with antiviral defense systems indicates that Synechococcus in Xiamen have evolved genetic traits in response to cyanophage infection. Studying the community diversity and functional genes of Synechococcus provides insights into their role in environmental adaptation and marine ecosystems.

Ecophysiological analysis reveals distinct environmental preferences in closely related Baltic Sea picocyanobacteria

Cluster 5 picocyanobacteria significantly contribute to primary productivity in aquatic ecosystems. Estuarine populations are highly diverse and consist of many co-occurring strains, but their physiology remains largely understudied. In this study, we characterized 17 novel estuarine picocyanobacterial strains. Phylogenetic analysis of the 16S rRNA and pigment genes (cpcB and cpeBA) uncovered multiple estuarine and freshwater-related clusters and pigment types. Assays with five representative strains (three phycocyanin rich and two phycoerythrin rich) under temperature (10-30°C), light (10-190 μmol photons m-2  s-1 ), and salinity (2-14 PSU) gradients revealed distinct growth optima and tolerance, indicating that genetic variability was accompanied by physiological diversity. Adaptability to environmental conditions was associated with differential pigment content and photosynthetic performance. Amplicon sequence variants at a coastal and an offshore station linked population dynamics with phylogenetic clusters, supporting that strains isolated in this study represent key ecotypes within the Baltic Sea picocyanobacterial community. The functional diversity found within strains with the same pigment type suggests that understanding estuarine picocyanobacterial ecology requires analysis beyond the phycocyanin and phycoerythrin divide. This new knowledge of the environmental preferences in estuarine picocyanobacteria is important for understanding and evaluating productivity in current and future ecosystems.

Distinct community assembly processes and habitat specialization driving the biogeographic patterns of abundant and rare bacterioplankton in a brackish coastal lagoon

The ecological diversity patterns and community assembly processes along spatio-temporal scales are least studied in the bacterioplankton sub-communities of brackish coastal lagoons. We examined the biogeographic patterns and relative influences of different assembly processes in structuring the abundant and rare bacterioplankton sub-communities of Chilika, the largest brackish water coastal lagoon of India. Rare taxa demonstrated significantly higher α- and β-diversity and biogeochemical functions than abundant taxa in the high-throughput 16S rRNA gene sequence dataset. The majority of the abundant taxa (91.4 %) were habitat generalists with a wider niche breadth (niche breadth index, B = 11.5), whereas most of the rare taxa (95.2 %) were habitat specialists with a narrow niche breadth (B = 8.9). Abundant taxa exhibited a stronger distance-decay relationship and higher spatial turnover rate than rare taxa. β-diversity partitioning revealed that the contribution of species turnover (72.2-97.8 %) was greater than nestedness (2.2-27.8 %) in causing the spatial variation in both abundant and rare taxa. Null model analyses revealed that the distribution of abundant taxa was mostly structured by stochastic processes (62.8 %), whereas deterministic processes (54.1 %) played a greater role in the rare taxa. However, the balance of these two processes varied across spatio-temporal scales in the lagoon. Salinity was the key deterministic factor controlling the variation of both abundant and rare taxa. Potential interaction networks showed a higher influence of negative interactions, indicating that species exclusion and top-down processes played a greater role in the community assembly. Notably, abundant taxa emerged as keystone taxa across spatio-temporal scales, suggesting their greater influences on other bacterial co-occurrences and network stability. Overall, this study provided detailed mechanistic insights into biogeographic patterns and underlying community assembly processes of the abundant and rare bacterioplankton over spatio-temporal scales in a brackish lagoon.

Adaptive genetic traits in pelagic freshwater microbes

Pelagic microbes have adopted distinct strategies to inhabit the pelagial of lakes and oceans and can be broadly categorized in two groups: free-living, specialized oligotrophs and patch-associated generalists or copiotrophs. In this review, we aim to identify genomic traits that enable pelagic freshwater microbes to thrive in their habitat. To do so, we discuss the main genetic differences of pelagic marine and freshwater microbes that are both dominated by specialized oligotrophs and the difference to freshwater sediment microbes, where copiotrophs are more prevalent. We phylogenomically analysed a collection of >7700 metagenome-assembled genomes, classified habitat preferences on different taxonomic levels, and compared the metabolic traits of pelagic freshwater, marine, and freshwater sediment microbes. Metabolic differences are mainly associated with transport functions, environmental information processing, components of the electron transport chain, osmoregulation and the isoelectric point of proteins. Several lineages with known habitat transitions (Nitrososphaeria, SAR11, Methylophilaceae, Synechococcales, Flavobacteriaceae, Planctomycetota) and the underlying mechanisms in this process are discussed in this review. Additionally, the distribution, ecology and genomic make-up of the most abundant freshwater prokaryotes are described in details in separate chapters for Actinobacteriota, Bacteroidota, Burkholderiales, Verrucomicrobiota, Chloroflexota, and 'Ca. Patescibacteria'.

Diversity dynamics of aerobic anoxygenic phototrophic bacteria in a freshwater lake

Aerobic anoxygenic photoheterotrophic (AAP) bacteria represent a functional group of prokaryotic organisms that harvests light energy using bacteriochlorophyll-containing photosynthetic reaction centers. They represent an active and rapidly growing component of freshwater bacterioplankton, with the highest numbers observed usually in summer. Species diversity of freshwater AAP bacteria has been studied before in lakes, but its seasonal dynamics remain unknown. In this report, we analysed temporal changes in the composition of the phototrophic community in an oligo-mesotrophic freshwater lake using amplicon sequencing of the pufM marker gene. The AAP community was dominated by phototrophic Gammaproteobacteria and Alphaproteobacteria, with smaller contribution of phototrophic Chloroflexota and Gemmatimonadota. Phototrophic Eremiobacteriota or members of Myxococcota were not detected. Interestingly, some AAP taxa, such as Limnohabitans, Rhodoferax, Rhodobacterales or Rhizobiales, were permanently present over the sampling period, while others, such as Sphingomonadales, Rhodospirillales or Caulobacterales appeared only transiently. The environmental factors that best explain the seasonal changes in AAP community were temperature, concentrations of oxygen and dissolved organic matter.

Global Phylogeography of Marine Synechococcus in Coastal Areas Reveals Strong Community Shifts

Marine Synechococcus comprise a numerically and ecologically prominent phytoplankton group, playing a major role in both carbon cycling and trophic networks in all oceanic regions except in the polar oceans. Despite their high abundance in coastal areas, our knowledge of Synechococcus communities in these environments is based on only a few local studies. Here, we use the global metagenome data set of the Ocean Sampling Day (June 21^st, 2014) to get a snapshot of the taxonomic composition of coastal Synechococcus communities worldwide, by recruitment on a reference database of 141 picocyanobacterial genomes, representative of the whole Prochlorococcus, Synechococcus, and Cyanobium diversity. This allowed us to unravel drastic community shifts over small to medium scale gradients of environmental factors, in particular along European coasts. The combined analysis of the phylogeography of natural populations and the thermophysiological characterization of eight strains, representative of the four major Synechococcus lineages (clades I to IV), also brought novel insights about the differential niche partitioning of clades I and IV, which most often co-dominate the Synechococcus community in cold and temperate coastal areas. Altogether, this study reveals several important characteristics and specificities of the coastal communities of Synechococcus worldwide. IMPORTANCE Synechococcus is the second most abundant phytoplanktonic organism on Earth, and its wide genetic diversity allowed it to colonize all the oceans except for polar waters, with different clades colonizing distinct oceanic niches. In recent years, the use of global metagenomics data sets has greatly improved our knowledge of "who is where" by describing the distribution of Synechococcus clades or ecotypes in the open ocean. However, little is known about the global distribution of Synechococcus ecotypes in coastal areas, where Synechococcus is often the dominant phytoplanktonic organism. Here, we leverage the global Ocean Sampling Day metagenomics data set to describe Synechococcus community composition in coastal areas worldwide, revealing striking community shifts, in particular along the coasts of Europe. As temperature appears as an important driver of the community composition, we also characterize the thermal preferenda of 8 Synechococcus strains, bringing new insights into the adaptation to temperature of the dominant Synechococcus clades.

The Evolutionary Kaleidoscope of Rhodopsins

Rhodopsins are widely distributed across all domains of life where they perform a plethora of functions through the conversion of electromagnetic radiation into physicochemical signals. As a result of an extensive survey of available genomic and metagenomic sequencing data, we reported the existence of novel clades and exotic sequence motifs scattered throughout the evolutionary radiations of both Type-1 and Type-3 rhodopsins that will likely enlarge the optogenetics toolbox. We expanded the typical rhodopsin blueprint by showing that a highly conserved and functionally important arginine residue (i.e., Arg82) was substituted multiple times during evolution by an extensive amino acid spectrum. We proposed the umbrella term Alt-rhodopsins (AltRs) for all such proteins that departed Arg82 orthodoxy. Some AltRs formed novel clades in the rhodopsin phylogeny and were found in giant viruses. Some newly uncovered AltRs were phylogenetically close to heliorhodopsins, which allowed a closer examination of the phylogenetic border between Type-1 rhodopsins and heliorhodopsins. Comprehensive phylogenetic trees and ancestral sequence reconstructions allowed us to advance the hypothesis that proto-heliorhodopsins were a eukaryotic innovation before their subsequent diversification into the extant Type-3 rhodopsins. IMPORTANCE The rhodopsin scaffold is remarkably versatile and widespread, coupling light availability to energy production and other light-dependent cellular responses with minor alterations to critical residues. We described an unprecedented spectrum of substitutions at one of the most conserved amino acids in the rhodopsin fold, Arg82. We denoted such phylogenetically diverse rhodopsins with the umbrella name Alt-rhodopsins (AltR) and described a distinct branch of AltRs in giant viruses. Intriguingly, some AltRs were the closest phylogenetic neighbors to Heliorhodopsins (HeRs) whose origins have remained enigmatic. Our analyses of HeR origins in the light of AltRs led us to posit a most unusual evolutionary trajectory that suggested a eukaryotic origin for HeRs before their diversification in prokaryotes.

α-cyanobacteria possessing form IA RuBisCO globally dominate aquatic habitats

RuBisCO (ribulose 1,5-bisphosphate carboxylase/oxygenase) is one the most abundant enzymes on Earth. Virtually all food webs depend on its activity to supply fixed carbon. In aerobic environments, RuBisCO struggles to distinguish efficiently between CO2 and O2. To compensate, organisms have evolved convergent solutions to concentrate CO2 around the active site. The genetic engineering of such inorganic carbon concentrating mechanisms (CCMs) into plants could help facilitate future global food security for humankind. In bacteria, the carboxysome represents one such CCM component, of which two independent forms exist: α and β. Cyanobacteria are important players in the planet's carbon cycle and the vast majority of the phylum possess a β-carboxysome, including most cyanobacteria used as laboratory models. The exceptions are the exclusively marine Prochlorococcus and Synechococcus that numerically dominate open ocean systems. However, the reason why marine systems favor an α-form is currently unknown. Here, we report the genomes of 58 cyanobacteria, closely related to marine Synechococcus that were isolated from freshwater lakes across the globe. We find all these isolates possess α-carboxysomes accompanied by a form 1A RuBisCO. Moreover, we demonstrate α-cyanobacteria dominate freshwater lakes worldwide. Hence, the paradigm of a separation in carboxysome type across the salinity divide does not hold true, and instead the α-form dominates all aquatic systems. We thus question the relevance of β-cyanobacteria as models for aquatic systems at large and pose a hypothesis for the reason for the success of the α-form in nature.

The combined impact of low temperatures and shifting phosphorus availability on the competitive ability of cyanobacteria

In freshwater systems, cyanobacteria are strong competitors under enhanced temperature and eutrophic conditions. Understanding their adaptive and evolutionary potential to multiple environmental states allows us to accurately predict their response to future conditions. To better understand if the combined impacts of temperature and nutrient limitation could suppress the cyanobacterial blooms, a single strain of Microcystis aeruginosa was inoculated into natural phytoplankton communities with different nutrient conditions: oligotrophic, eutrophic and eutrophic with the addition of bentophos. We found that the use of the bentophos treatment causes significant differences in prokaryotic and eukaryotic communities. This resulted in reduced biodiversity among the eukaryotes and a decline in cyanobacterial abundance suggesting phosphorus limitation had a strong impact on the community structure. The low temperature during the experiment lead to the disappearance of M. aeruginosa in all treatments and gave other phytoplankton groups a competitive advantage leading to the dominance of the eukaryotic families that have diverse morphologies and nutritional modes. These results show cyanobacteria have a reduced competitive advantage under certain temperature and nutrient limiting conditions and therefore, controlling phosphorus concentrations could be a possible mitigation strategy for managing harmful cyanobacterial blooms in a future warmer climate.

Elucidating the picocyanobacteria salinity divide through ecogenomics of new freshwater isolates

BACKGROUND

Cyanobacteria are the major prokaryotic primary producers occupying a range of aquatic habitats worldwide that differ in levels of salinity, making them a group of interest to study one of the major unresolved conundrums in aquatic microbiology which is what distinguishes a marine microbe from a freshwater one? We address this question using ecogenomics of a group of picocyanobacteria (cluster 5) that have recently evolved to inhabit geographically disparate salinity niches. Our analysis is made possible by the sequencing of 58 new genomes from freshwater representatives of this group that are presented here, representing a 6-fold increase in the available genomic data.

RESULTS

Overall, freshwater strains had larger genomes (≈2.9 Mb) and %GC content (≈64%) compared to brackish (2.69 Mb and 64%) and marine (2.5 Mb and 58.5%) isolates. Genomic novelties/differences across the salinity divide highlighted acidic proteomes and specific salt adaptation pathways in marine isolates (e.g., osmolytes/compatible solutes - glycine betaine/ggp/gpg/gmg clusters and glycerolipids glpK/glpA), while freshwater strains possessed distinct ion/potassium channels, permeases (aquaporin Z), fatty acid desaturases, and more neutral/basic proteomes. Sulfur, nitrogen, phosphorus, carbon (photosynthesis), or stress tolerance metabolism while showing distinct genomic footprints between habitats, e.g., different types of transporters, did not obviously translate into major functionality differences between environments. Brackish microbes show a mixture of marine (salt adaptation pathways) and freshwater features, highlighting their transitional nature.

CONCLUSIONS

The plethora of freshwater isolates provided here, in terms of trophic status preference and genetic diversity, exemplifies their ability to colonize ecologically diverse waters across the globe. Moreover, a trend towards larger and more flexible/adaptive genomes in freshwater picocyanobacteria may hint at a wider number of ecological niches in this environment compared to the relatively homogeneous marine system.

Seasonal and Spatial Variations in Synechococcus Abundance and Diversity Throughout the Gullmar Fjord, Swedish Skagerrak

The picophytoplankton Synechococcus is a globally abundant autotroph that contributes significantly to primary production in the oceans and coastal areas. These cyanobacteria constitute a diverse genus of organisms that have developed independent niche spaces throughout aquatic environments. Here, we use the 16S V3–V4 rRNA gene region and flow cytometry to explore the diversity of Synechococcus within the picophytoplankton community in the Gullmar Fjord, on the west coast of Sweden. We conducted a station-based 1-year time series and two transect studies of the fjord. Our analysis revealed that within the large number of Synechococcus amplicon sequence variants (ASVs; 239 in total), prevalent ASVs phylogenetically clustered with clade representatives in both marine subcluster 5.1 and 5.2. The near-surface composition of ASVs shifted from spring to summer, when a 5.1 subcluster dominated community developed along with elevated Synechococcus abundances up to 9.3 × 104 cells ml–1. This seasonal dominance by subcluster 5.1 was observed over the length of the fjord (25 km), where shifts in community composition were associated with increasing depth. Unexpectedly, the community shift was not associated with changes in salinity. Synechococcus abundance dynamics also differed from that of the photosynthetic picoeukaryote community. These results highlight how seasonal variations in environmental conditions influence the dynamics of Synechococcus clades in a high latitude threshold fjord.

Diversity and Evolution of Pigment Types in Marine Synechococcus Cyanobacteria

Synechococcus cyanobacteria are ubiquitous and abundant in the marine environment and contribute to an estimated 16% of the ocean net primary productivity. Their light-harvesting complexes, called phycobilisomes (PBS), are composed of a conserved allophycocyanin core, from which radiates six to eight rods with variable phycobiliprotein and chromophore content. This variability allows Synechococcus cells to optimally exploit the wide variety of spectral niches existing in marine ecosystems. Seven distinct pigment types or subtypes have been identified so far in this taxon based on the phycobiliprotein composition and/or the proportion of the different chromophores in PBS rods. Most genes involved in their biosynthesis and regulation are located in a dedicated genomic region called the PBS rod region. Here, we examine the variability of gene content and organization of this genomic region in a large set of sequenced isolates and natural populations of Synechococcus representative of all known pigment types. All regions start with a tRNA-PheGAA and some possess mobile elements for DNA integration and site-specific recombination, suggesting that their genomic variability relies in part on a "tycheposon"-like mechanism. Comparison of the phylogenies obtained for PBS and core genes revealed that the evolutionary history of PBS rod genes differs from the core genome and is characterized by the co-existence of different alleles and frequent allelic exchange. We propose a scenario for the evolution of the different pigment types and highlight the importance of incomplete lineage sorting in maintaining a wide diversity of pigment types in different Synechococcus lineages despite multiple speciation events.

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.

Genome Sequence of Synechococcus sp. Strain LA31, Isolated from a Temperate Estuary

Cluster 5 Synechococcus species are widely acknowledged for their broad distribution and biogeochemical importance. In particular, subcluster 5.2 strains inhabit freshwater, estuarine, and marine environments but are understudied, compared to other subclusters. Here, we present the genome for Synechococcus sp. strain LA31, a strain that was recently isolated from Narragansett Bay, Rhode Island, USA.

Microbial population genomes from the Amazon River reveal possible modulation of the organic matter degradation process in tropical freshwaters

Rivers connect the carbon cycle in land with that in aquatic ecosystems by transporting and transforming terrestrial organic matter (TeOM). The Amazon River receives huge loads of TeOM from the surrounding rainforest, promoting a substantial microbial heterotrophic activity and consequently, CO₂ outgassing. In the Amazon River, microbes degrade up to 55% of the lignin present in the TeOM. Yet, the main microbial genomes involved in TeOM degradation were unknown. Here, we characterize 51 population genomes (PGs) representing some of the most abundant microbes in the Amazon River deriving from 106 metagenomes. The 51 reconstructed PGs are among the most abundant microbes in the Amazon River, and 53% of them are not able to degrade TeOM. Among the PGs capable of degrading TeOM, 20% were exclusively cellulolytic, while the others could also oxidize lignin. The transport and consumption of lignin oxidation byproducts seemed to be decoupled from the oxidation process, being apparently performed by different groups of microorganisms. By connecting the genomic features of abundant microbes in the Amazon River with the degradation machinery of TeOM, we suggest that a complex microbial consortium could explain the quick turnover of TeOM previously observed in this ecosystem.

Effect of hydrogen peroxide on natural phytoplankton and bacterioplankton in a drinking water reservoir: Mesocosm-scale study

Cyanobacterial blooms are increasingly reported worldwide, presenting a challenge to water treatment plants and concerning risks to human health and aquatic ecosystems. Advanced oxidative processes comprise efficient and safe methods for water treatment. Hydrogen peroxide (H₂O₂) has been proposed as a sustainable solution to mitigate bloom-forming cyanobacteria since this group presents a higher sensitivity compared to other phytoplankton, with no major risks to the environment at low concentrations. Here, we evaluated the effects of a single H₂O₂ addition (10 mg L^-1) over 120 h in mesocosms introduced in a reservoir located in a semi-arid region presenting a Planktothrix-dominated cyanobacterial bloom. We followed changes in physical and chemical parameters and in the bacterioplankton composition. H₂O₂ efficiently suppressed cyanobacteria, green algae, and diatoms over 72 h, leading to an increase in transparency and dissolved organic carbon, and a decrease in dissolved oxygen and pH, while nutrient concentrations were not affected. After 120 h, cyanobacterial abundance remained low and green algae became dominant. 16S rRNA sequencing revealed that the original cyanobacterial bloom was composed by Planktothrix, Cyanobium and Microcystis. Only Cyanobium increased in relative abundance at 120 h, suggesting regrowth. A prominent change in the composition of heterotrophic bacteria was observed with Exiguobacterium, Paracoccus and Deinococcus becoming the most abundant genera after the H₂O₂ treatment. Our results indicate that this approach is efficient in suppressing cyanobacterial blooms and improving water quality in tropical environments. Monitoring changes in abiotic parameters and the relative abundance of specific bacterial taxa could be used to anticipate the regrowth of cyanobacteria after H₂O₂ degradation and to indicate where in the reservoir H₂O₂ should be applied so the effects are still felt in the water treatment plant intake.

Seasonal Distribution of Cyanobacteria in Three Urban Eutrophic Lakes Results from an Epidemic-like Response to Environmental Conditions

Cyanobacterial communities of three co-located eutrophic sandpit lakes were surveyed during 2016 and 2017 over season and depth using high-throughput DNA sequencing of the 16S rRNA gene. All three lakes were stratified except during April 2017 when the lakes were recovering from a strong mixing event. 16S rRNA gene V4 sequences were parsed into operational taxonomic units (OTUs) at 99% sequence identity. After rarefaction of 139 samples to 25,000 sequences per sample, a combined total of 921,529 partial 16S rRNA gene sequences were identified as cyanobacteria. They were binned into 19,588 unique cyanobacterial OTUs. Of these OTUs, 11,303 were Cyanobium. Filamentous Planktothrix contributed 1537 and colonial Microcystis contributed 265. The remaining 6482 OTUs were considered unclassified. For Planktothrix and Microcystis one OTU accounted for greater than 95% of the total sequences for each genus. However, in both cases the non-dominant OTUs clustered with the dominant OTUs by date, lake, and depth. All Planktothrix OTUs and a single Cyanobium OTU were detected below the oxycline. All other Cyanobium and Microcystis OTUs were detected above the oxycline. The distribution of Cyanobium OTUs between lakes and seasons can be explained by an epidemic-like response where individual OTUs clonally rise from a diverse hypolimnion population when conditions are appropriate. The importance of using 99% identity over the more commonly used 97% is discussed with respect to cyanobacterial community structure. The approach described here can provide another valuable tool for assessing cyanobacterial populations and provide greater insight into the controls of cyanobacterial blooms.

Responses of abundant and rare bacterioplankton to temporal change in a subtropical urban reservoir

Investigation of bacterial community dynamics across different time scales is important for understanding how environmental conditions drive community change over time. Bacterioplankton from the surface waters of a subtropical urban reservoir in southeast China were analyzed through high-frequency sampling over 13 months to compare patterns and ecological processes between short (0‒8 weeks), medium (9‒24 weeks) and long (25‒53 weeks) time intervals. We classified the bacterial community into different subcommunities: abundant taxa (AT); conditionally rare taxa (CRT); rare taxa (RT). CRT contributed > 65% of the alpha-diversity, and temporal change of beta-diversities was more pronounced for AT and CRT than RT. The bacterial community exhibited a directional change in the short- and medium-time intervals and a convergent dynamic during the long-time interval due to a seasonal cycle. Cyanobacteria exhibited a strong succession pattern than other phyla. CRT accounted for > 76% of the network nodes in three stations. The bacteria-environment relationship and deterministic processes were stronger for large sample size at station G (n = 116) than small sample size at stations C (n = 12) and L (n = 22). These findings suggest that a high-frequency sampling approach can provide a better understanding on the time scales at which bacterioplankton can change fast between being abundant or rare, thus providing the facts about environmental factors driving microbial community dynamics. Patterns and processes in alpha- and beta-diversities and community assembly of bacterioplankton differ among different time intervals (short-, medium- and long-time intervals) and different subcommunities (abundant, conditionally rare and rare taxa) in a subtropical urban reservoir, demonstrating the importance of temporal scale and high-frequency sampling in microbial community ecology.

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.

Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes

Cyanorak v2.1 (http://www.sb-roscoff.fr/cyanorak) is an information system dedicated to visualizing, comparing and curating the genomes of Prochlorococcus, Synechococcus and Cyanobium, the most abundant photosynthetic microorganisms on Earth. The database encompasses sequences from 97 genomes, covering most of the wide genetic diversity known so far within these groups, and which were split into 25,834 clusters of likely orthologous groups (CLOGs). The user interface gives access to genomic characteristics, accession numbers as well as an interactive map showing strain isolation sites. The main entry to the database is through search for a term (gene name, product, etc.), resulting in a list of CLOGs and individual genes. Each CLOG benefits from a rich functional annotation including EggNOG, EC/K numbers, GO terms, TIGR Roles, custom-designed Cyanorak Roles as well as several protein motif predictions. Cyanorak also displays a phyletic profile, indicating the genotype and pigment type for each CLOG, and a genome viewer (Jbrowse) to visualize additional data on each genome such as predicted operons, genomic islands or transcriptomic data, when available. This information system also includes a BLAST search tool, comparative genomic context as well as various data export options. Altogether, Cyanorak v2.1 constitutes an invaluable, scalable tool for comparative genomics of ecologically relevant marine microorganisms.

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.

Picophytoplankton dynamics in a large temperate estuary and impacts of extreme storm events

Picophytoplankton (PicoP) are increasingly recognized as significant contributors to primary productivity and phytoplankton biomass in coastal and estuarine systems. Remarkably though, PicoP composition is unknown or not well-resolved in several large estuaries including the semi-lagoonal Neuse River Estuary (NRE), a tributary of the second largest estuary-system in the lower USA, the Pamlico-Albemarle Sound. The NRE is impacted by extreme weather events, including recent increases in precipitation and flooding associated with tropical cyclones. Here we examined the impacts of moderate to extreme (Hurricane Florence, September 2018) precipitation events on NRE PicoP abundances and composition using flow cytometry, over a 1.5 year period. Phycocyanin-rich Synechococcus-like cells were the most dominant PicoP, reaching ~ 106 cells mL-1, which highlights their importance as key primary producers in this relatively long residence-time estuary. Ephemeral "blooms" of picoeukaryotic phytoplankton (PEUK) during spring and after spikes in river flow were also detected, making PEUK periodically major contributors to PicoP biomass (up to ~ 80%). About half of the variation in PicoP abundance was explained by measured environmental variables. Temperature explained the most variation (24.5%). Change in total dissolved nitrogen concentration, an indication of increased river discharge, explained the second-most variation in PicoP abundance (15.9%). The short-term impacts of extreme river discharge from Hurricane Florence were particularly evident as PicoP biomass was reduced by ~ 100-fold for more than 2 weeks. We conclude that precipitation is a highly influential factor on estuarine PicoP biomass and composition, and show how 'wetter' future climate conditions will have ecosystem impacts down to the smallest of phytoplankton.

Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

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.

Microorganisms of Lake Baikal-the deepest and most ancient lake on Earth

Lake Baikal (Russia) is the largest (by volume) and deepest lake on Earth. The lake remains relatively pristine due to the low population density around its basin. Being very distant from any marine water body but having a remarkable number of similarities to oceans (depth, oxygen content, oligotrophy) provides a unique model of pelagic microbiota that is submitted to marine-like conditions minus the salt content of the water. It is also a model of lakes located at high latitudes and submitted to yearly ice cover (from January to April). The analysis by different approaches has indeed provided a view of the microbiota of this lake. It contains novel microbes that are closely related to marine groups not known to be present in freshwater like Chloroflexi or Pelagibacter. The deep water mass contains large communities of chemolithotrophs that use ammonia generated in the photic zone or methane from the sediments. KEY POINTS: • The chemical composition and limnic features of the deepest lake on Earth determine the vital activity of microorganisms. • The diversity, ecology, and role of individual taxa of microorganisms were studied using cultivation and molecular methods. • Data of large metagenomic datasets in the epipelagic and bathypelagic layers of the water column in southern Baikal were discussed.

Interaction of Cyanobacteria with Nanometer and Micron Sized Polystyrene Particles in Marine and Fresh Water

Microplastics and nanoplastics are emerging pollutants, widespread both in marine and in freshwater environments. Cyanobacteria are also ubiquitous in water and play a vital role in natural ecosystems, using photosynthesis to produce oxygen. Using photography, fluorescence microscopy and cryogenic and scanning electron microscopy (cryo-SEM, SEM) we investigated the physicochemical response of one of the most predominant seawater cyanobacteria (Synechococcus elongatus, PCC 7002) and freshwater cyanobacteria (S. elongatus Nageli PCC 7942) when exposed to 10 μm diameter polystyrene (microPS) and 100 nm diameter polystyrene (nanoPS) particles. Marine and freshwater cyanobacteria formed aggregates with the nanoPS, bound together by extracellular polymeric substances (EPS), and these aggregates sedimented. The aggregates were larger, and the sedimentation was more rapid for the marine system. Aggregate morphologies were qualitatively different for the microPS samples, with the bacteria linking up a small number of particles, all held together by EPS. There was no sedimentation in these samples. The cyanobacteria remained alive after exposure to the particles. The particle size- and salt concentration-dependent response of cyanobacteria to these anthropogenic stressors is an important factor to consider for a proper understanding of the fate of the particles as well as the bacteria.

Seasonal and spatial dynamics of bacterioplankton communities in a brackish water coastal lagoon

Coastal ecosystems, one of the most productive ecosystems, are subjected to natural and anthropogenic stresses. Coastal bacterioplankton communities are highly dynamic due to spatiotemporal heterogeneity in the environmental parameters. We investigated the seasonal and spatial variation in bacterioplankton communities, their abundances and environmental drivers during one year period in Chilika, a brackish water coastal lagoon of India. High-throughput sequencing of 16S rRNA genes of bacterioplankton communities showed that they were dominated by heterotrophs namely α-Proteobacteria SAR11 and their sub-clades (SAR11_Ib, Chesapeake-Delaware_Bay, Candidatus_Pelagibacter, and SAR11_Surface_1), actinobacterial lineages (hgcI, CL500-29, and Candidatus_Aquiluna), β-Proteobacteria MWH-UniP1, β-Proteobacteria OM43, and verrucomicrobial clade Spartobacteria 'LD29'. Synechococcus was the dominant member within autotrophic cyanobacterial community. Response ratio derived from comparisons of taxon-specific absolute abundances and indicator analyses showed that SAR11_Surface_1 sub-clade occupied high-salinity environment especially during summer and winter and emerged as a strong indicator for mesohaline-polyhaline salinity regime. In contrast, Spartobacteria 'LD29', Actinobacteria hgcI, and CL500-29 preferred low-salinity freshwater environment and were strong indicators for oligohaline-mesohaline regimes. Spatiotemporal patterns were governed by 'distance-decay' and 'similarity-time' relationships. Bacterioplankton communities were mostly determined by salinity, dissolved oxygen, phosphate, and pH which resulted 'species sorting' leading to biogeographical patterns in the bacterioplankton communities. Modeling analysis revealed the characteristic shift in the indicator bacterioplankton taxa along with estuarine salinity gradient. This study has provided baseline information on the bacterioplankton communities and their environmental drivers within an anthropogenically impacted cyclone prone coastal lagoon which would be useful in assessing the impact of multiple stressors on this vulnerable ecosystem.

Marine-freshwater prokaryotic transitions require extensive changes in the predicted proteome

BACKGROUND

The adaptation of a marine prokaryote to live in freshwater environments or vice versa is generally believed to be an unusual and evolutionary demanding process. However, the reasons are not obvious given the similarity of both kinds of habitats.

RESULTS

We have found major differences at the level of the predicted metaproteomes of marine and freshwater habitats with more acidic values of the isoelectric points (pI) in marine microbes. Furthermore, by comparing genomes of marine-freshwater phylogenetic relatives, we have found higher pI values (basic shift) in the freshwater ones. This difference was sharper in secreted > cytoplasmic > membrane proteins. The changes are concentrated on the surface of soluble proteins. It is also detectable at the level of total amino acid composition and involves similarly core and flexible genome- encoded proteins.

CONCLUSIONS

The marked changes at the level of protein amino acid composition and pI provide a tool to predict the preferred habitat of a culture or a metagenome-assembled genome (MAG). The exact physiological explanation for such variations in the pIs and electrostatic surface potentials is not known yet. However, these changes might reflect differences in membrane bioenergetics derived from the absence of significant Na⁺ concentrations in most freshwater habitats. In any case, the changes in amino acid composition in most proteins imply that a long evolutionary time is required to adapt from one type of habitat to the other.

High-throughput DNA sequencing reveals the dominance of pico- and other filamentous cyanobacteria in an urban freshwater Lake

The current study presents findings related to algal blooms in a fresh water lake, which has been experiencing severe cyanobacterial blooms (CyanoHABs). Primarily, picocyanobacteria belonging to the genus Synechococcus and filamentous cyanobacterial group belonging to Aphanizomenon and Dolichospermum dominated top water column during non-bloom and bloom periods respectively. The dominance of Synechococcus in early summer informs that blooming in Utah Lake starts in early summer and then later is taken over by other bloom-forming cyanobacteria, such as species belonging to the genus Aphanizomenon. A strong negative correlation (r = -0.9, p < 0.001) was found between the occurrence of Aphanizomenon and Synechococcus which correlates very well with the fact that the blooms of these two different cyanobacteria never coexisted. The predominance of cyanobacteria in 2017 was attributed more to temperature (r = 0.18, p < 0.001). The Actinobacteria was negatively correlated with primary production and high chlorophyll a concentration. Flavobacterium and Limnohabitans were the main phytoplankton colonizers and predators detected that could secrete extracellular enzymes to degrade algal exudates (such as proteins and polysaccharides). Additionally, cyanotoxins producers Microcystis aeruginosa and Planktothrix accounted for up to 12.43% and 7.04% of total cyanobacteria abundance during blooms. The relative abundance of chloroplast reads was overall lower than the cyanobacteria reads, except for the May 5th sampling in 2017. There was inter-annual variability in the bloom-associated heterotrophic bacterial populations, but these populations were consistent with bloom-associated bacterial populations found in other lakes. Community diversity analysis for both Shannon and Simpson indices indicated lower community diversity during the bloom period. The beta diversity conducted by PCoA and UPGMA trees suggested the significant temporal rather than spatial impacts on shaping the phytoplankton community structures during the summer season.

Insights Into the Evolution of Picocyanobacteria and Phycoerythrin Genes (mpeBA and cpeBA)

Marine picocyanobacteria, Prochlorococcus and Synechococcus, substantially contribute to marine primary production and have been the subject of extensive ecological and genomic studies. Little is known about their close relatives from freshwater and non-marine environments. Phylogenomic analyses (using 136 proteins) provide strong support for the monophyly of a clade of non-marine picocyanobacteria consisting of Cyanobium, Synechococcus and marine Sub-cluster 5.2; this clade itself is sister to marine Synechococcus and Prochlorococcus. The most basal lineage within the Syn/Pro clade, Sub-Cluster 5.3, includes marine and freshwater strains. Relaxed molecular clock (SSU, LSU) analyses show that while ancestors of the Syn/Pro clade date as far back as the end of the Pre-Cambrian, modern crown groups evolved during the Carboniferous and Triassic. Comparative genomic analyses reveal novel gene cluster arrangements involved in phycobilisome (PBS) metabolism in freshwater strains. Whilst PBS genes in marine Synechococcus are mostly found in one type of phycoerythrin (PE) rich gene cluster (Type III), strains from non-marine habitats, so far, appear to be more diverse both in terms of pigment content and gene arrangement, likely reflecting a wider range of habitats. Our phylogenetic analyses show that the PE genes (mpeBA) evolved via a duplication of the cpeBA genes in an ancestor of the marine and non-marine picocyanobacteria and of the symbiotic strains Synechococcus spongiarum. A 'primitive' Type III-like ancestor containing cpeBA and mpeBA had thus evolved prior to the divergence of the Syn/Pro clade and S. spongiarum. During the diversification of Synechococcus lineages, losses of mpeBA genes may explain the emergence of pigment cluster Types I, II, IIB, and III in both marine and non-marine habitats, with few lateral gene transfer events in specific taxa.

Complete Genome Sequence and Comparative Analysis of Synechococcus sp. CS-601 (SynAce01), a Cold-Adapted Cyanobacterium from an Oligotrophic Antarctic Habitat

Marine picocyanobacteria belonging to Synechococcus are major contributors to the global carbon cycle, however the genomic information of its cold-adapted members has been lacking to date. To fill this void the genome of a cold-adapted planktonic cyanobacterium Synechococcus sp. CS-601 (SynAce01) has been sequenced. The genome of the strain contains a single chromosome of approximately 2.75 MBp and GC content of 63.92%. Gene prediction yielded 2984 protein coding sequences and 44 tRNA genes. The genome contained evidence of horizontal gene transfer events during its evolution. CS-601 appears as a transport generalist with some specific adaptation to an oligotrophic marine environment. It has a broad repertoire of transporters of both inorganic and organic nutrients to survive in inhospitable environments. The cold adaptation of the strain exhibited characteristics of a psychrotroph rather than psychrophile. Its salt adaptation strategy is likely to rely on the uptake and synthesis of osmolytes, like glycerol or glycine betaine. Overall, the genome reveals two distinct patterns of adaptation to the inhospitable environment of Antarctica. Adaptation to an oligotrophic marine environment is likely due to an abundance of genes, probably acquired horizontally, that are associated with increased transport of nutrients, osmolytes, and light harvesting. On the other hand, adaptations to low temperatures are likely due to prolonged evolutionary changes.

Ecological and genomic features of two widespread freshwater picocyanobacteria

We present two genomes of widespread freshwater picocyanobacteria isolated by extinction dilution from a Spanish oligotrophic reservoir. Based on microscopy and genomic properties, both picocyanobacteria were tentatively designated Synechococcus lacustris Tous, formerly described as a metagenome assembled genome (MAG) from the same habitat, and Cyanobium usitatum Tous, described here for the first time. Both strains were purified in unicyanobacterial cultures, and their genomes were sequenced. They are broadly distributed in freshwater systems; the first seems to be a specialist on temperate reservoirs (Tous, Amadorio, Dexter, Lake Lanier, Sparkling), and the second appears to also be abundant in cold environments including ice-covered lakes such as Lake Baikal, Lake Erie or the brackish Baltic Sea. Having complete genomes provided access to the flexible genome that does not assemble in MAGs. We found several genomic islands in both genomes, within which there were genes for nitrogen acquisition, transporters for a wide set of compounds and biosynthesis of phycobilisomes in both strains. Some of these regions of low coverage in metagenomes also included antimicrobial compounds, transposases and phage defence systems, including a novel type III CRISPR-Cas phage defence system that was only detected in Synechococcus lacustris Tous.

A New Freshwater Cyanosiphovirus Harboring Integrase

Pelagic cyanobacteria are key players in the functioning of aquatic ecosystems, and their viruses (cyanophages) potentially affect the abundance and composition of cyanobacterial communities. Yet, there are few well-described freshwater cyanophages relative to their marine counterparts, and in general, few cyanosiphoviruses (family Siphoviridae) have been characterized, limiting our understanding of the biology and the ecology of this prominent group of viruses. Here, we characterize S-LBS1, a freshwater siphovirus lytic to a phycoerythrin-rich Synechococcus isolate (Strain TCC793). S-LBS1 has a narrow host range, a burst size of ∼400 and a relatively long infecting step before cell lysis occurs. It has a dsDNA 34,641 bp genome with putative genes for structure, DNA packing, lysis, replication, host interactions, DNA repair and metabolism. S-LBS1 is similar in genome size, genome architecture, and gene content, to previously described marine siphoviruses also infecting PE-rich Synechococcus, e.g., S-CBS1 and S-CBS3. However, unlike other Synechococcus phages, S-LBS1 encodes an integrase, suggesting its ability to establish lysogenic relationships with its host. Sequence recruitment from viral metagenomic data showed that S-LBS1-like viruses are diversely present in a wide range of aquatic environments, emphasizing their potential importance in controlling and structuring Synechococcus populations. A comparative analysis with 16 available sequenced cyanosiphoviruses reveals the absence of core genes within the genomes, suggesting high degree of genetic variability in siphoviruses infecting cyanobacteria. It is likely that cyanosiphoviruses have evolved as distinct evolutionary lineages and that adaptive co-evolution occurred between these viruses and their hosts (i.e., Synechococcus, Prochlorococcus, Nodularia, and Acaryochloris), constituting an important driving force for such phage diversification.

Genome analysis of the freshwater planktonic Vulcanococcus limneticus sp. nov. reveals horizontal transfer of nitrogenase operon and alternative pathways of nitrogen utilization

BACKGROUND

Many cyanobacteria are capable of fixing atmospheric nitrogen, playing a crucial role in biogeochemical cycling. Little is known about freshwater unicellular cyanobacteria Synechococcus spp. at the genomic level, despite being recognised of considerable ecological importance in aquatic ecosystems. So far, it has not been shown whether these unicellular picocyanobacteria have the potential for nitrogen fixation. Here, we present the draft-genome of the new pink-pigmented Synechococcus-like strain Vulcanococcus limneticus. sp. nov., isolated from the volcanic Lake Albano (Central Italy).

RESULTS

The novel species Vulcanococcus limneticus sp. nov. falls inside the sub-cluster 5.2, close to the estuarine/marine strains in a maximum-likelihood phylogenetic tree generated with 259 marker genes with representatives from marine, brackish, euryhaline and freshwater habitats. V.limneticus sp. nov. possesses a complete nitrogenase and nif operon. In an experimental setup under nitrogen limiting and non-limiting conditions, growth was observed in both cases. However, the nitrogenase genes (nifHDK) were not transcribed, i.e., V.limneticus sp. nov. did not fix nitrogen, but instead degraded the phycobilisomes to produce sufficient amounts of ammonia. Moreover, the strain encoded many other pathways to incorporate ammonia, nitrate and sulphate, which are energetically less expensive for the cell than fixing nitrogen. The association of the nif operon to a genomic island, the relatively high amount of mobile genetic elements (52 transposases) and the lower observed GC content of V.limneticus sp. nov. nif operon (60.54%) compared to the average of the strain (68.35%) support the theory that this planktonic strain may have obtained, at some point of its evolution, the nif operon by horizontal gene transfer (HGT) from a filamentous or heterocystous cyanobacterium.

CONCLUSIONS

In this study, we describe the novel species Vulcanococcus limneticus sp. nov., which possesses a complete nif operon for nitrogen fixation. The finding that in our experimental conditions V.limneticus sp. nov. did not express the nifHDK genes led us to reconsider the actual ecological meaning of these accessory genes located in genomic island that have possibly been acquired via HGT.

The future of genomics in polar and alpine cyanobacteria

In recent years, genomic analyses have arisen as an exciting way of investigating the functional capacity and environmental adaptations of numerous micro-organisms of global relevance, including cyanobacteria. In the extreme cold of Arctic, Antarctic and alpine environments, cyanobacteria are of fundamental ecological importance as primary producers and ecosystem engineers. While their role in biogeochemical cycles is well appreciated, little is known about the genomic makeup of polar and alpine cyanobacteria. In this article, we present ways that genomic techniques might be used to further our understanding of cyanobacteria in cold environments in terms of their evolution and ecology. Existing examples from other environments (e.g. marine/hot springs) are used to discuss how methods developed there might be used to investigate specific questions in the cryosphere. Phylogenomics, comparative genomics and population genomics are identified as methods for understanding the evolution and biogeography of polar and alpine cyanobacteria. Transcriptomics will allow us to investigate gene expression under extreme environmental conditions, and metagenomics can be used to complement tradition amplicon-based methods of community profiling. Finally, new techniques such as single cell genomics and metagenome assembled genomes will also help to expand our understanding of polar and alpine cyanobacteria that cannot readily be cultured.

Light color acclimation is a key process in the global ocean distribution of Synechococcus cyanobacteria

Marine Synechococcus cyanobacteria are major contributors to global oceanic primary production and exhibit a unique diversity of photosynthetic pigments, allowing them to exploit a wide range of light niches. However, the relationship between pigment content and niche partitioning has remained largely undetermined due to the lack of a single-genetic marker resolving all pigment types (PTs). Here, we developed and employed a robust method based on three distinct marker genes (cpcBA, mpeBA, and mpeW) to estimate the relative abundance of all known Synechococcus PTs from metagenomes. Analysis of the Tara Oceans dataset allowed us to reveal the global distribution of Synechococcus PTs and to define their environmental niches. Green-light specialists (PT 3a) dominated in warm, green equatorial waters, whereas blue-light specialists (PT 3c) were particularly abundant in oligotrophic areas. Type IV chromatic acclimaters (CA4-A/B), which are able to dynamically modify their light absorption properties to maximally absorb green or blue light, were unexpectedly the most abundant PT in our dataset and predominated at depth and high latitudes. We also identified populations in which CA4 might be nonfunctional due to the lack of specific CA4 genes, notably in warm high-nutrient low-chlorophyll areas. Major ecotypes within clades I-IV and CRD1 were preferentially associated with a particular PT, while others exhibited a wide range of PTs. Altogether, this study provides important insights into the ecology of Synechococcus and highlights the complex interactions between vertical phylogeny, pigmentation, and environmental parameters that shape Synechococcus community structure and evolution.

Genomes of Novel Microbial Lineages Assembled from the Sub-Ice Waters of Lake Baikal

We present a metagenomic study of Lake Baikal (East Siberia). Two samples obtained from the water column under the ice cover (5 and 20 m deep) in March 2016 have been deep sequenced and the reads assembled to generate metagenome-assembled genomes (MAGs) that are representative of the microbes living in this special environment. Compared with freshwater bodies studied around the world, Lake Baikal had an unusually high fraction of Verrucomicrobia Other groups, such as Actinobacteria and Proteobacteria, were in proportions similar to those found in other lakes. The genomes (and probably cells) tended to be small, presumably reflecting the extremely oligotrophic and cold prevalent conditions. Baikal microbes are novel lineages recruiting very little from other water bodies and are distantly related to other freshwater microbes. Despite their novelty, they showed the closest relationship to genomes discovered by similar approaches from other freshwater lakes and reservoirs. Some of them were particularly similar to MAGs from the Baltic Sea, which, although it is brackish, connected to the ocean, and much more eutrophic, has similar climatological conditions. Many of the microbes contained rhodopsin genes, indicating that, in spite of the decreased light penetration allowed by the thick ice/snow cover, photoheterotrophy could be widespread in the water column, either because enough light penetrates or because the microbes are already adapted to the summer ice-less conditions. We have found a freshwater SAR11 subtype I/II representative showing striking synteny with Pelagibacterubique strains, as well as a phage infecting the widespread freshwater bacterium PolynucleobacterIMPORTANCE Despite the increasing number of metagenomic studies on different freshwater bodies, there is still a missing component in oligotrophic cold lakes suffering from long seasonal frozen cycles. Here, we describe microbial genomes from metagenomic assemblies that appear in the upper water column of Lake Baikal, the largest and deepest freshwater body on Earth. This lake is frozen from January to May, which generates conditions that include an inverted temperature gradient (colder up), decrease in light penetration due to ice, and, especially, snow cover, and oligotrophic conditions more similar to the open-ocean and high-altitude lakes than to other freshwater or brackish systems. As could be expected, most reconstructed genomes are novel lineages distantly related to others in cold environments, like the Baltic Sea and other freshwater lakes. Among them, there was a broad set of streamlined microbes with small genomes/intergenic spacers, including a new nonmarine Pelagibacter-like (subtype I/II) genome.

Reconstruction of Diverse Verrucomicrobial Genomes from Metagenome Datasets of Freshwater Reservoirs

The phylum Verrucomicrobia contains freshwater representatives which remain poorly studied at the genomic, taxonomic, and ecological levels. In this work we present eighteen new reconstructed verrucomicrobial genomes from two freshwater reservoirs located close to each other (Tous and Amadorio, Spain). These metagenome-assembled genomes (MAGs) display a remarkable taxonomic diversity inside the phylum and comprise wide ranges of estimated genome sizes (from 1.8 to 6 Mb). Among all Verrucomicrobia studied we found some of the smallest genomes of the Spartobacteria and Opitutae classes described so far. Some of the Opitutae family MAGs were small, cosmopolitan, with a general heterotrophic metabolism with preference for carbohydrates, and capable of xylan, chitin, or cellulose degradation. Besides, we assembled large copiotroph genomes, which contain a higher number of transporters, polysaccharide degrading pathways and in general more strategies for the uptake of nutrients and carbohydrate-based metabolic pathways in comparison with the representatives with the smaller genomes. The diverse genomes revealed interesting features like green-light absorbing rhodopsins and a complete set of genes involved in nitrogen fixation. The large diversity in genome sizes and physiological properties emphasize the diversity of this clade in freshwaters enlarging even further the already broad eco-physiological range of these microbes.