The mesopelagic anoxic Black Sea as an unexpected habitat for Synechococcus challenges our understanding of global "deep red fluorescence"

Uncovering the genomic basis of symbiotic interactions and niche adaptations in freshwater picocyanobacteria

BACKGROUND

Picocyanobacteria from the genera Prochlorococcus, Synechococcus, and Cyanobium are the most widespread photosynthetic organisms in aquatic ecosystems. However, their freshwater populations remain poorly explored, due to uneven and insufficient sampling across diverse inland waterbodies.

RESULTS

In this study, we present 170 high-quality genomes of freshwater picocyanobacteria from non-axenic cultures collected across Central Europe. In addition, we recovered 33 genomes of their potential symbiotic partners affiliated with four genera, Pseudomonas, Mesorhizobium, Acidovorax, and Hydrogenophaga. The genomic basis of symbiotic interactions involved heterotrophs benefiting from picocyanobacteria-derived nutrients while providing detoxification of ROS. The global abundance patterns of picocyanobacteria revealed ecologically significant ecotypes, associated with trophic status, temperature, and pH as key environmental factors. The adaptation of picocyanobacteria in (hyper-)eutrophic waterbodies could be attributed to their colonial lifestyles and CRISPR-Cas systems. The prevailing CRISPR-Cas subtypes in picocyanobacteria were I-G and I-E, which appear to have been acquired through horizontal gene transfer from other bacterial phyla.

CONCLUSIONS

Our findings provide novel insights into the population diversity, ecology, and evolutionary strategies of the most widespread photoautotrophs within freshwater ecosystems. Video Abstract.

Mixotrophy in cyanobacteria

Cyanobacteria evolved the oxygenic photosynthesis to generate organic matter from CO2 and sunlight, and they were responsible for the production of oxygen in the Earth's atmosphere. This made them a model for photosynthetic organisms, since they are easier to study than higher plants. Early studies suggested that only a minority among cyanobacteria might assimilate organic compounds, being considered mostly autotrophic for decades. However, compelling evidence from marine and freshwater cyanobacteria, including toxic strains, in the laboratory and in the field, has been obtained in the last decades: by using physiological and omics approaches, mixotrophy has been found to be a more widespread feature than initially believed. Furthermore, dominant clades of marine cyanobacteria can take up organic compounds, and mixotrophy is critical for their survival in deep waters with very low light. Hence, mixotrophy seems to be an essential trait in the metabolism of most cyanobacteria, which can be exploited for biotechnological purposes.

Effects of environmental parameters on spatial and temporal distribution of marine microbial communities in the southern Black Sea

The Black Sea is a unique environment with strong and permanent vertical stratification, with a thin layer of oxic zone above and a permanent anoxic zone below. Few high-throughput genomic surveys have been conducted to examine microbiota in the Black Sea. Yet, there is no study on the seasonal and vertical variation in microbial community compositions, driving forces and mechanisms of community assembly. In this study, seasonal, vertical, and spatial microbial assemblages were studied in terms of diversity, abundance, and community structure using 16S rRNA metabarcoding. 16S rRNA metabarcoding confirmed seasonal changes in microbial communities and the presence of distinct microbial groups among different water layers. Taxa belonging to Cyanobiaceae contributed a large fraction of the total biomass and were the most abundant autotrophic bacteria found across the whole water column, including hydrogen sulfide-containing anoxic zone. Temperature, salinity, water density, conductivity, light, chlorophyll-a, O2, NO3, NH3, PO4, Si, and H2S had a significant influence on the vertical bacterial community assemblages. The copper mine discharge system at 180 m did not affect microbial community structure and composition. Temperature seemed to be a primary factor in the variance between shallow depths. In conclusion, the lack of light, low dissolved oxygen levels, and low temperature do not restrict microbial diversity, as proven by the higher diversity observed in deeper zones. Wastewater in Black Sea region may be discharged into the Black Sea to depth of 180 m or deeper without impacting microbial ecology.

Suboxic DOM is bioavailable to surface prokaryotes in a simulated overturn of an oxygen minimum zone, Devil's Hole, Bermuda

Oxygen minimum zones (OMZs) are expanding due to increased sea surface temperatures, subsequent increased oxygen demand through respiration, reduced oxygen solubility, and thermal stratification driven in part by anthropogenic climate change. Devil's Hole, Bermuda is a model ecosystem to study OMZ microbial biogeochemistry because the formation and subsequent overturn of the suboxic zone occur annually. During thermally driven stratification, suboxic conditions develop, with organic matter and nutrients accumulating at depth. In this study, the bioavailability of the accumulated dissolved organic carbon (DOC) and the microbial community response to reoxygenation of suboxic waters was assessed using a simulated overturn experiment. The surface inoculated prokaryotic community responded to the deep (formerly suboxic) 0.2 μm filtrate with cell densities increasing 2.5-fold over 6 days while removing 5 μmol L-1 of DOC. After 12 days, the surface community began to shift, and DOC quality became less diagenetically altered along with an increase in SAR202, a Chloroflexi that can degrade recalcitrant dissolved organic matter (DOM). Labile DOC production after 12 days coincided with an increase of Nitrosopumilales, a chemoautotrophic ammonia oxidizing archaea (AOA) that converts ammonia to nitrite based on the ammonia monooxygenase (amoA) gene copy number and nutrient data. In comparison, the inoculation of the deep anaerobic prokaryotic community into surface 0.2 μm filtrate demonstrated a die-off of 25.5% of the initial inoculum community followed by a 1.5-fold increase in cell densities over 6 days. Within 2 days, the prokaryotic community shifted from a Chlorobiales dominated assemblage to a surface-like heterotrophic community devoid of Chlorobiales. The DOM quality changed to less diagenetically altered material and coincided with an increase in the ribulose-1,5-bisphosphate carboxylase/oxygenase form I (cbbL) gene number followed by an influx of labile DOM. Upon reoxygenation, the deep DOM that accumulated under suboxic conditions is bioavailable to surface prokaryotes that utilize the accumulated DOC initially before switching to a community that can both produce labile DOM via chemoautotrophy and degrade the more recalcitrant DOM.

Phytoplankton pigment dynamics in marine lake fluctuating between stratified and holomictic euxinic conditions

Biomass dynamics in the marine lake are strongly dependent on seasonal variability in vertical stratification, indicating rapid adaptation of phytoplankton to short-term changes in the water column. A small marine lake (Rogoznica Lake, Croatia), which fluctuates between stably stratified and holomictic euxinic conditions, was used as a model to study the phytoplankton responses to environmental perturbations, in particular the anoxic stress, caused by periodic holomixia. The epilimnion showed significant temporal and vertical variability with a chlorophyll a subsurface maximum with the highest biomass near the chemocline. Fucoxanthin-containing biomass (diatoms) dominated in the epilimnion in colder seasons and was first to recover after holomictic euxinic events. The shift towards the smaller groups prevailed during highly stratified water column conditions in warmer seasons. Results for the hypolimnion were more enigmatic, with high concentrations of alloxanthin, zeaxanthin, and violaxanthin indicating the presence of a viable small-size mixotrophic community under extreme conditions.

Genomic and Transcriptomic Insights into Salinity Tolerance-Based Niche Differentiation of Synechococcus Clades in Estuarine and Coastal Waters

Cluster 5 Synechococcus is one of the most important primary producers on earth. However, ecotypes of this genus exhibit complex geographical distributions, and the genetic basis of niche partitioning is still not fully understood. Here, we report distinct distributions of subcluster 5.1 (SC5.1) and subcluster 5.2 (SC5.2) Synechococcus in estuarine waters, and we reveal that salinity is the main factor determining their distribution. Clade III (belonging to SC5.1) and CB4 (belonging to SC5.2) are dominant clades in the study region, with different ecological distributions. We further conducted physiological, genomic, and transcriptomic studies of Synechococcus strains YX04-3 and HK05, which are affiliated with clade III and CB4, respectively. Laboratory tests showed that HK05 could grow at low salinity (13 ppt), whereas the growth of YX04-3 was suppressed when salinity decreased to 13 ppt. Genomic and transcriptomic analysis suggested that euryhaline clade CB4 is capable of dealing with a sudden drop of salinity by releasing compatible solutes through mechanosensitive channels that are coded by the mscL gene, decreasing biosynthesis of organic osmolytes, and increasing expression of heat shock proteins and high light-inducible proteins to protect photosystem. Furthermore, CB4 strain HK05 exhibited a higher growth rate when growing at low salinity than at high salinity. This is likely achieved by reducing its biosynthesis of organic osmolyte activity and increasing its photosynthetic activity at low salinity, which allowed it to enhance the assimilation of inorganic carbon and nitrogen. Together, these results provide new insights regarding the ecological distribution of SC5.2 and SC5.1 ecotypes and their underlying molecular mechanisms. IMPORTANCE Synechococcus is a group of unicellular Cyanobacteria that are widely distributed in global aquatic ecosystems. Salinity is a factor that affects the distribution of microorganisms in estuarine and coastal environments. In this study, we studied the distribution pattern of Synechococcus community along the salinity gradient in a subtropical estuary. By using omic methods, we unveiled genetic traits that determine the niche partitioning of euryhaline and strictly marine Synechococcus. We also explored the strategies employed by euryhaline Synechococcus to cope with a sudden drop of salinity, and revealed possible mechanisms for the higher growth rate of euryhaline Synechococcus in low salinity conditions. This study provides new insight into the genetic basis of niche partitioning of Synechococcus clades.

Antibiotic resistance genes correlate with metal resistances and accumulate in the deep water layers of the Black Sea

Seas and oceans are a global reservoir of antibiotic resistance genes (ARGs). Only a few studies investigated the dynamics of ARGs along the water column of the Black Sea, a unique environment, with a peculiar geology, biology and history of anthropogenic pollution. In this study, we analyzed metagenomic data from two sampling campaigns (2013 and 2019) collected across three different sites in the Western Black Sea at depths ranging from 5 to 2000 m. The data were processed to annotate ARGs, metal resistance genes (MRGs) and integron integrase genes. The ARG abundance was significantly higher in the deep water layers and depth was the main driver of beta-diversity both for ARGs and MRGs. Moreover, ARG and MRG abundances strongly correlated (r = 0.95). The integron integrase gene abundances and composition were not influenced by the water depth and did not correlate with ARGs. The analysis of the obtained MAGs showed that some of them harbored intI gene together with several ARGs and MRGs, suggesting the presence of multidrug resistant bacteria and that MRGs and integrons could be involved in the selection of ARGs. These results demonstrate that the Black Sea is not only an important reservoir of ARGs, but also that they accumulate in the deep water layers where co-selection with MRGs could be assumed as a relevant mechanism of their persistence.

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.

mebipred: identifying metal binding potential in protein sequence

Motivation

metal-binding proteins have a central role in maintaining life processes. Nearly one-third of known protein structures contain metal ions that are used for a variety of needs, such as catalysis, DNA/RNA binding, protein structure stability, etc. Identifying metal-binding proteins is thus crucial for understanding the mechanisms of cellular activity. However, experimental annotation of protein metal-binding potential is severely lacking, while computational techniques are often imprecise and of limited applicability.

Results

we developed a novel machine learning-based method, mebipred, for identifying metal-binding proteins from sequence-derived features. This method is over 80% accurate in recognizing proteins that bind metal ion-containing ligands; the specific identity of 11 ubiquitously present metal ions can also be annotated. mebipred is reference-free, i.e. no sequence alignments are involved, and is thus faster than alignment-based methods; it is also more accurate than other sequence-based prediction methods. Additionally, mebipred can identify protein metal-binding capabilities from short sequence stretches, e.g. translated sequencing reads, and, thus, may be useful for the annotation of metal requirements of metagenomic samples. We performed an analysis of available microbiome data and found that ocean, hot spring sediments and soil microbiomes use a more diverse set of metals than human host-related ones. For human microbiomes, physiological conditions explain the observed metal preferences. Similarly, subtle changes in ocean sample ion concentration affect the abundance of relevant metal-binding proteins. These results highlight mebipred’s utility in analyzing microbiome metal requirements.

Availability and implementation

mebipred is available as a web server at services.bromberglab.org/mebipred and as a standalone package at https://pypi.org/project/mymetal/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Picocyanobacterial Synechococcus in marine ecosystem: Insights from genetic diversity, global distribution, and potential function

Marine Synechococcus, a main group of picocyanobacteria, has been ubiquitously observed across the global oceans. Synechococcus exhibits high phylogenetical and phenotypical diversity, and horizontal gene transfer makes its genetic evolution much more intricate. With the development of measurement technologies and analysis methods, the genomic information and niche partition of each Synechococcus lineage tend to be precisely described, but the global analysis is still lacking. Therefore, it is necessary to summarize existing studies and integrate published data to gain a comprehensive understanding of Synechococcus on genetic variation, niche division, and potential functions. In this review, the maximum likelihood trees are constructed based on existing sequence data, including both phylogenetic and pigmentary gene markers. The global distribution characteristics of abundance, lineages, and pigment types are concluded through pooled analysis of more than 700 samples obtained from approximately 50 scientific research cruises. The potential functions of Synechococcus are explored in element cycles and biological interactions. Future work on Synechococcus is suggested to focus on not only elucidating the nature of Synechococcus biodiversity but also demonstrating its interactions with the ecosystem by combining bioinformatics and macroscopic isotope-labeled environmental parameters.

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-Phe_GAA 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.

Metabarcoding Reveals Lacustrine Picocyanobacteria Respond to Environmental Change Through Adaptive Community Structuring

Picocyanobacteria (Pcy) are important yet understudied components of lake foodwebs. While phylogenetic studies of isolated strains reveal a high diversity of freshwater genotypes, little is known about abiotic drivers associated with Pcy in different lakes. Due to methodological limitations, most previous studies assess potential drivers using total cell abundances as a response, with often conflicting and inconsistent results. In the present study, we explored how picocyanobacterial communities respond to environmental change using a combination of epifluorescence microscopy and community data determined using 16S rRNA gene metabarcoding. Temporal shifts in picocyanobacterial abundance, diversity and community dynamics were assessed in relation to potential environmental drivers in five contrasting lakes over 1year. Cell abundances alone were not consistently related to environmental variables across lakes. However, the addition of metabarcoding data revealed diverse picocyanobacterial communities that differed significantly between lakes, driven by environmental variables related to trophic state. Within each lake, communities were temporally dynamic and certain amplicon sequence variants (ASVs) were strongly associated with specific environmental drivers. Rapid shifts in community structure and composition were often related to environmental changes, indicating that lacustrine Pcy can persist at high abundances through collective community adaptation. These results demonstrate that a combination of microscopy and metabarcoding enables an in-depth characterisation of picocyanobacterial communities and reveals strain-specific drivers. We recommend that future studies cease referring to picocyanobacterial as one functional group and take strain specific variability into consideration.

Coping with darkness: The adaptive response of marine picocyanobacteria to repeated light energy deprivation

The picocyanobacteria Prochlorococcus and Synechococcus are found throughout the ocean's euphotic zone, where the daily light:dark cycle drives their physiology. Periodic deep mixing events can, however, move cells below this region, depriving them of light for extended periods of time. Here, we demonstrate that members of these genera can adapt to tolerate repeated periods of light energy deprivation. Strains kept in the dark for 3 d and then returned to the light initially required 18-26 d to resume growth, but after multiple rounds of dark exposure they began to regrow after only 1-2 d. This dark-tolerant phenotype was stable and heritable; some cultures retained the trait for over 132 generations even when grown in a standard 13:11 light:dark cycle. We found no genetic differences between the dark-tolerant and parental strains of Prochlorococcus NATL2A, indicating that an epigenetic change is likely responsible for the adaptation. To begin to explore this possibility, we asked whether DNA methylation-one potential mechanism mediating epigenetic inheritance in bacteria-occurs in Prochlorococcus. LC-MS/MS analysis showed that while DNA methylations, including 6 mA and 5 mC, are found in some other Prochlorococcus strains, there were no methylations detected in either the parental or dark-tolerant NATL2A strains. These findings suggest that Prochlorococcus utilizes a yet-to-be-determined epigenetic mechanism to adapt to the stress of extended light energy deprivation, and highlights phenotypic heterogeneity as an additional dimension of Prochlorococcus diversity.

The Thermosynechococcus Genus: Wide Environmental Distribution, but a Highly Conserved Genomic Core

Cyanobacteria thrive in diverse environments. However, questions remain about possible growth limitations in ancient environmental conditions. As a single genus, the Thermosynechococcus are cosmopolitan and live in chemically diverse habitats. To understand the genetic basis for this, we compared the protein coding component of Thermosynechococcus genomes. Supplementing the known genetic diversity of Thermosynechococcus, we report draft metagenome-assembled genomes of two Thermosynechococcus recovered from ferrous carbonate hot springs in Japan. We find that as a genus, Thermosynechococcus is genomically conserved, having a small pan-genome with few accessory genes per individual strain as well as few genes that are unique to the genus. Furthermore, by comparing orthologous protein groups, including an analysis of genes encoding proteins with an iron related function (uptake, storage or utilization), no clear differences in genetic content, or adaptive mechanisms could be detected between genus members, despite the range of environments they inhabit. Overall, our results highlight a seemingly innate ability for Thermosynechococcus to inhabit diverse habitats without having undergone substantial genomic adaptation to accommodate this. The finding of Thermosynechococcus in both hot and high iron environments without adaptation recognizable from the perspective of the proteome has implications for understanding the basis of thermophily within this clade, and also for understanding the possible genetic basis for high iron tolerance in cyanobacteria on early Earth. The conserved core genome may be indicative of an allopatric lifestyle-or reduced genetic complexity of hot spring habitats relative to other environments.

Dynamics and Distribution of Marine Synechococcus Abundance and Genotypes during Seasonal Hypoxia in a Coastal Marine Ranch

Marine Synechococcus are an ecologically important picocyanobacterial group widely distributed in various oceanic environments. Little is known about the dynamics and distribution of Synechococcus abundance and genotypes during seasonal hypoxia in coastal zones. In this study, an investigation was conducted in a coastal marine ranch along two transects in Muping, Yantai, where hypoxic events (defined here as the dissolved oxygen concentration <3 mg L−1) occurred in the summer of 2015. The hypoxia occurred in the bottom waters from late July and persisted until late August. It was confined at nearshore stations of the two transects, one running across a coastal ranch and the other one outside. During this survey, cell abundance of Synechococcus was determined with flow cytometry, showing great variations ranging from 1 × 104 to 3.0 × 105 cells mL−1, and a bloom of Synechococcus occurred when stratification disappeared and hypoxia faded out outside the ranch. Regression analysis indicated that dissolved oxygen, pH, and inorganic nutrients were the most important abiotic factors in explaining the variation in Synechococcus cell abundance. Diverse genotypes (mostly belonged to the sub-clusters 5.1 and 5.2) were detected using clone library sequencing and terminal restriction fragment length polymorphism analysis of the 16S–23S rRNA internal transcribed spacer region. The richness of genotypes was significantly related to salinity, temperature, silicate, and pH, but not dissolved oxygen. Two environmental factors, temperature and salinity, collectively explained 17% of the variation in Synechococcus genotype assemblage. With the changes in population composition in diverse genotypes, the Synechococcus assemblages survived in the coastal hypoxia event and thrived when hypoxia faded out.

Cyanophage Diversity and Community Structure in Dead Zone Sediments

Cyanophages are viruses that target cyanobacteria and directly control their abundance via viral lysis. Cyanobacteria are known to cause large blooms in water bodies, substantially contributing to oxygen depletion in bottom waters resulting in areas called dead zones.

Up to 20% of prokaryotic organisms in the oceans are estimated to die every day due to viral infection and lysis. Viruses can therefore alter microbial diversity, community structure, and biogeochemical processes driven by these organisms. Cyanophages are viruses that infect and lyse cyanobacterial cells, adding bioavailable carbon and nutrients into the environment. Cyanobacteria are photosynthesizing bacteria, with some species capable of N₂ fixation, which are known to form large blooms as well as resistant resting cells known as akinetes. Here, we investigated cyanophage diversity and community structure plus cyanobacteria in dead zone sediments. We sampled surface sediments and sequenced DNA and RNA, along an oxygen gradient—representing oxic, hypoxic, and anoxic conditions—in one of the world’s largest dead zones located in the Baltic Sea. Cyanophages were detected at all stations and, based on partial genome contigs, had a higher alpha diversity and different beta diversity in the hypoxic-anoxic sediments, suggesting that cyanobacteria in dead zone sediments and/or environmental conditions select for specific cyanophages. Some of these cyanophages can infect cyanobacteria with potential consequences for gene expression related to their photosystem and phosphate regulation. Top cyanobacterial genera detected in the anoxic sediment included Dolichospermum/Anabaena, Synechococcus, and Cyanobium. RNA transcripts classified to cyanobacteria were associated with numerous pathways, including anaerobic carbon metabolism and N₂ fixation. Cyanobacterial blooms are known to fuel oxygen-depleted ecosystems with phosphorus (so-called internal loading), and our cyanophage data indicate the potential for viral lysis of cyanobacteria which might explain the high nutrient turnover in these environments.

IMPORTANCE Cyanophages are viruses that target cyanobacteria and directly control their abundance via viral lysis. Cyanobacteria are known to cause large blooms in water bodies, substantially contributing to oxygen depletion in bottom waters resulting in areas called dead zones. Our knowledge of cyanophages in dead zones is very scarce, and so far, no studies have assembled partial cyanophage genomes and investigated their associated cyanobacteria in these dark and anoxic sediments. Here, we present the first study using DNA and RNA sequencing to investigate in situ diversity of cyanophages and cyanobacteria in dead zones. Our study shows that dead zone sediments contain different cyanophages compared to oxic sediments and suggest that these viruses are able to affect cyanobacterial photosystem and phosphate regulation. Furthermore, cyanophage-controlled lysis of cyanobacteria might also increase the turnover of carbon, phosphorus, and nitrogen in these oxygen-free environments at the bottom of the sea.

The microbiome of the Black Sea water column analyzed by shotgun and genome centric metagenomics

BACKGROUND

The Black Sea is the largest brackish water body in the world, although it is connected to the Mediterranean Sea and presents an upper water layer similar to some regions of the former, albeit with lower salinity and temperature. Despite its well-known hydrology and physicochemical features, this enormous water mass remains poorly studied at the microbial genomics level.

RESULTS

We have sampled its different water masses and analyzed the microbiome by shotgun and genome-resolved metagenomics, generating a large number of metagenome-assembled genomes (MAGs) from them. We found various similarities with previously described Black Sea metagenomic datasets, that show remarkable stability in its microbiome. Our datasets are also comparable to other marine anoxic water columns like the Cariaco Basin. The oxic zone resembles to standard marine (e.g. Mediterranean) photic zones, with Cyanobacteria (Synechococcus but a conspicuously absent Prochlorococcus), and photoheterotrophs domination (largely again with marine relatives). The chemocline presents very different characteristics from the oxic surface with many examples of chemolithotrophic metabolism (Thioglobus) and facultatively anaerobic microbes. The euxinic anaerobic zone presents, as expected, features in common with the bottom of meromictic lakes with a massive dominance of sulfate reduction as energy-generating metabolism, a few (but detectable) methanogenesis marker genes, and a large number of "dark matter" streamlined genomes of largely unpredictable ecology.

CONCLUSIONS

The Black Sea oxic zone presents many similarities to the global ocean while the redoxcline and euxinic water masses have similarities to other similar aquatic environments of marine (Cariaco Basin or other Black Sea regions) or freshwater (meromictic monimolimnion strata) origin. The MAG collection represents very well the different types of metabolisms expected in this kind of environment. We are adding critical information about this unique and important ecosystem and its microbiome.

Anoxic chlorophyll maximum enhances local organic matter remineralization and nitrogen loss in Lake Tanganyika

In marine and freshwater oxygen-deficient zones, the remineralization of sinking organic matter from the photic zone is central to driving nitrogen loss. Deep blooms of photosynthetic bacteria, which form the suboxic/anoxic chlorophyll maximum (ACM), widespread in aquatic ecosystems, may also contribute to the local input of organic matter. Yet, the influence of the ACM on nitrogen and carbon cycling remains poorly understood. Using a suite of stable isotope tracer experiments, we examined the transformation of nitrogen and carbon under an ACM (comprising of Chlorobiaceae and Synechococcales) and a non-ACM scenario in the anoxic zone of Lake Tanganyika. We find that the ACM hosts a tight coupling of photo/litho-autotrophic and heterotrophic processes. In particular, the ACM was a hotspot of organic matter remineralization that controlled an important supply of ammonium driving a nitrification-anammox coupling, and thereby played a key role in regulating nitrogen loss in the oxygen-deficient zone.

Enigmatic blooms of phytoplankton in aquatic oxygen-deficient zones could exacerbate depletion of nitrogen. Here the authors perform stable isotope experiments on the oxygen-deficient waters of Lake Tanganyika in Africa, finding that blooms drive down fixed nitrogen and could expand as a result of climate change.

Discovery of Euryhaline Phycoerythrobilin-Containing Synechococcus and Its Mechanisms for Adaptation to Estuarine Environments

Understanding the strategies developed by different microbial groups to adapt to specific niches is critical. Through genome and transcriptome analyses of two newly isolated novel euryhaline Synechococcus strains, this study revealed that cluster 5 phycoerythrobilin-containing Synechococcus, which are thought to be strictly marine strains, could be abundant in low-salinity waters of the Pearl River estuary (salinity <15 ppt) and explained the molecular mechanisms that enabled them to adapt the low and fluctuating salinity in the estuarine environment.

Synechococcus are among the most abundant and widely distributed picocyanobacteria on earth. Cluster 5 phycoerythrobilin-containing (PEB-containing) Synechococcus, the major marine Synechococcus, were considered to prefer high salinity, and they are absent in estuarine ecosystems. However, we have detected PEB-containing Synechococcus in some low-salinity (<15-ppt) areas of the Pearl River estuary at an abundance up to 1.0 × 10⁵ cells ml⁻¹. Two PEB-containing Synechococcus strains (HK01 and LTW-R) were isolated, and tests on them revealed their ability to cope with variations in the salinity (from 14 to 44 ppt). Phylogenetic analysis showed that HK01 belonged to a novel Synechococcus clade (HK1), whereas LTW-R was clustered with S5.2 strains. Whole-genome analysis revealed that a membrane channel protein with glycine zipper motifs is unique to euryhaline Synechococcus. The upregulation of this protein, the osmotic sensors, and the heat shock protein HSP20 and the downregulation of the osmolyte biosynthesis enable euryhaline Synechococcus to well adapt to the low and fluctuating salinity in the estuarine environment. In addition, decreasing the salinity in LTW-R strongly downregulated several important metabolic pathways, including photosynthesis, and the Calvin-Benson cycle, whereas its growth was not significantly affected. Moreover, obtaining PEB genes from horizontal gene transfer expands the light niche significantly for euryhaline Synechococcus. These results provided new insights into the life strategies and ecological function of marine PEB-containing Synechococcus under the unique environmental condition of estuarine waters, particularly in response to salinity variations.

IMPORTANCE Understanding the strategies developed by different microbial groups to adapt to specific niches is critical. Through genome and transcriptome analyses of two newly isolated novel euryhaline Synechococcus strains, this study revealed that cluster 5 phycoerythrobilin-containing Synechococcus, which are thought to be strictly marine strains, could be abundant in low-salinity waters of the Pearl River estuary (salinity <15 ppt) and explained the molecular mechanisms that enabled them to adapt the low and fluctuating salinity in the estuarine environment. This study expands current understanding on mechanisms involved in niche separation of marine Synechococcus lineages.

Spatial distribution of antibiotic and heavy metal resistance genes in the Black Sea

Antibiotic resistance genes (ARGs) and antibiotic resistant bacteria (ARB) are worldwide considered as emerging contaminants of large interest, and a primary threat to human health. It is becoming clear that the environment plays a central role in the transmission, spread, and evolution of antibiotic resistance. Although marine systems have been largely investigated, only a few studies have considered the presence of ARGs in meso- and bathypelagic waters. To date, no molecular based studies have yet been made to investigate the occurrence of ARGs in the Black Sea, the largest meromictic basin in the world, receiving water from a number of important European rivers and their residues of anthropogenic activities in permanently stratified mesopelagic water masses. In this study, we determined the presence and the abundance of five ARGs (bla_CTXM, ermB, qnrS, sul2, tetA) and of the heavy metal resistance gene (HMRG) czcA, in different sampling sites in the eastern and western Black Sea, at several depths (up to 1000 m) and various distances from the shoreline. Three ARGs (bla_CTXM, sul2, and tetA) and czcA were present in at least 43% of the analysed samples, whereas ermB and qnrS were never detected. In particular, sul2 abundances increased significantly in coastal location, whereas tetA increased with sampling depth. These findings point out the Black Sea as a source of ARGs and HMRGs distributed along the whole water column.

Population genetic structure of turbot (Scophthalmus maximus L., 1758) in the Black Sea

Turbot, Scophthalmus maximus, is a commercially important demersal flatfish species distributed throughout the Black Sea. Several studies performed locally with a limited number of specimens using both mitochondrial DNA (mtDNA) and microsatellite markers evidenced notable genetic variation among populations. However, comprehensive population genetic studies are required to help management of the species in the Black Sea. In the present study eight microsatellite loci were used to resolve the population structure of 414 turbot samples collected from 12 sites across the Black Sea. Moreover, two mtDNA genes, COI and Cyt-b, were used for taxonomic identification. Microsatellite markers of Smax-04 and B12-I GT14 were excluded from analysis due to scoring issues. Data analysis was performed with the remaining six loci. Loci were highly polymorphic (average of 17.8 alleles per locus), indicating high genetic variability. Locus 3/20CA17, with high null allele frequency (>30%), significantly deviated from HW equilibrium. Pairwise comparison of the F_ST index showed significant differences between most of the surveyed sampling sites (P < 0.01). Cluster analysis evidenced the presence of three genetic groups among sampling sites. Significant genetic differentiation between Northern (Sea of Azov and Crimea) and Southern (Turkish Black Sea Coast) Black Sea sampling sites were detected. The Mantel test supported an isolation by distance model of population structure. These findings are vital for long-term sustainable management of the species and development of conservation programs. Moreover, generated mtDNA sequences would be useful for the establishment of a database for S. maximus.

Genomic Comparison and Spatial Distribution of Different Synechococcus Phylotypes in the Black Sea

Picocyanobacteria of the genus Synechococcus are major contributors to global primary production and nutrient cycles due to their oxygenic photoautotrophy, their abundance, and the extensive distribution made possible by their wide-ranging biochemical capabilities. The recent recovery and isolation of strains from the deep euxinic waters of the Black Sea encouraged us to expand our analysis of their adaptability also beyond the photic zone of aquatic environments. To this end, we quantified the total abundance and distribution of Synechococcus along the whole vertical profile of the Black Sea by flow cytometry, and analyzed the data obtained in light of key environmental factors. Furthermore, we designed phylotype-specific primers using the genomes of two new epipelagic coastal strains – first described here – and of two previously described mesopelagic strains, analyzed their presence/abundance by qPCR, and tested this parameter also in metagenomes from two stations at different depths. Together, whole genome sequencing, metagenomics and qPCR techniques provide us with a higher resolution of Synechococcus dynamics in the Black Sea. Both phylotypes analyzed are abundant and successful in epipelagic coastal waters; but while the newly described epipelagic strains are specifically adapted to this environment, the strains previously isolated in mesopelagic waters are able, in low numbers, to withstand the aphotic and oxygen depleted conditions of deep layers. This heterogeneity allows different Synechococcus phylotypes to occupy different niches and underscores the importance of a more detailed characterization of the abundance, distribution, and dynamics of individual populations of these picocyanobacteria.

Influence of the polar light cycle on seasonal dynamics of an Antarctic lake microbial community

BACKGROUND

Cold environments dominate the Earth's biosphere and microbial activity drives ecosystem processes thereby contributing greatly to global biogeochemical cycles. Polar environments differ to all other cold environments by experiencing 24-h sunlight in summer and no sunlight in winter. The Vestfold Hills in East Antarctica contains hundreds of lakes that have evolved from a marine origin only 3000-7000 years ago. Ace Lake is a meromictic (stratified) lake from this region that has been intensively studied since the 1970s. Here, a total of 120 metagenomes representing a seasonal cycle and four summers spanning a 10-year period were analyzed to determine the effects of the polar light cycle on microbial-driven nutrient cycles.

RESULTS

The lake system is characterized by complex sulfur and hydrogen cycling, especially in the anoxic layers, with multiple mechanisms for the breakdown of biopolymers present throughout the water column. The two most abundant taxa are phototrophs (green sulfur bacteria and cyanobacteria) that are highly influenced by the seasonal availability of sunlight. The extent of the Chlorobium biomass thriving at the interface in summer was captured in underwater video footage. The Chlorobium abundance dropped from up to 83% in summer to 6% in winter and 1% in spring, before rebounding to high levels. Predicted Chlorobium viruses and cyanophage were also abundant, but their levels did not negatively correlate with their hosts.

CONCLUSION

Over-wintering expeditions in Antarctica are logistically challenging, meaning insight into winter processes has been inferred from limited data. Here, we found that in contrast to chemolithoautotrophic carbon fixation potential of Southern Ocean Thaumarchaeota, this marine-derived lake evolved a reliance on photosynthesis. While viruses associated with phototrophs also have high seasonal abundance, the negative impact of viral infection on host growth appeared to be limited. The microbial community as a whole appears to have developed a capacity to generate biomass and remineralize nutrients, sufficient to sustain itself between two rounds of sunlight-driven summer-activity. In addition, this unique metagenome dataset provides considerable opportunity for future interrogation of eukaryotes and their viruses, abundant uncharacterized taxa (i.e. dark matter), and for testing hypotheses about endemic species in polar aquatic ecosystems. Video Abstract.

Microbial diversity and long-term geochemical trends in the euxinic zone of a marine, meromictic lake

Hypoxic and anoxic niches of meromictic lakes are important sites for studying the microbial ecology of conditions resembling ancient Earth. The expansion and increasing global distribution of such environments also means that information about them serves to understand future phenomena. In this study, a long-term chemical dataset (1996-2015) was explored together with seasonal (in 2015) information on the diversity and abundance of bacterial and archaeal communities residing in the chemocline, monimolimnion and surface sediment of the marine meromictic Rogoznica Lake. The results of quantitative PCR assays, and high-throughput sequencing, targeting 16S rRNA genes and transcripts, revealed a clear vertical structure of the microbial community with Gammaproteobacteria (Halochromatium) and cyanobacteria (Synechococcus spp.) dominating the chemocline, Deltaproteobacteria and Bacteroidetes dominating the monimolimnion, and significantly more abundant archaeal populations in the surface sediment, most of which affiliated to Nanoarchaeota. Seasonal changes in the community structure and abundance were not pronounced. Diversity in Rogoznica Lake was found to be high, presumably as a consequence of stable environmental conditions accompanied by high dissolved carbon and nutrient concentrations. Long-term data indicated that Rogoznica Lake exhibited climate changes that could alter its physico-chemical features and, consequently, induce structural and physiological changes within its microbial community.