Metagenome-assembled genomes uncover a global brackish microbiome.

BACKGROUND: Microbes are main drivers of biogeochemical cycles in oceans and lakes. Although the genome is a foundation for understanding the metabolism, ecology and evolution of an organism, few bacterioplankton genomes have been sequenced, partly due to difficulties in cultivating them. RESULTS: We use automatic binning to reconstruct a large number of bacterioplankton genomes from a metagenomic time-series from the Baltic Sea, one of world's largest brackish water bodies. These genomes represent novel species within typical freshwater and marine clades, including clades not previously sequenced. The genomes' seasonal dynamics follow phylogenetic patterns, but with fine-grained lineage-specific variations, reflected in gene-content. Signs of streamlining are evident in most genomes, and estimated genome sizes correlate with abundance variation across filter size fractions. Comparing the genomes with globally distributed metagenomes reveals significant fragment recruitment at high sequence identity from brackish waters in North America, but little from lakes or oceans. This suggests the existence of a global brackish metacommunity whose populations diverged from freshwater and marine relatives over 100,000 years ago, long before the Baltic Sea was formed (8000 years ago). This markedly contrasts to most Baltic Sea multicellular organisms, which are locally adapted populations of freshwater or marine counterparts. CONCLUSIONS: We describe the gene content, temporal dynamics and biogeography of a large set of new bacterioplankton genomes assembled from metagenomes. We propose that brackish environments exert such strong selection that lineages adapted to them flourish globally with limited influence from surrounding aquatic communities.

Microbial metagenomics in the Baltic Sea: Recent advancements and prospects for environmental monitoring.

Metagenomics refers to the analysis of DNA from a whole community. Metagenomic sequencing of environmental DNA has greatly improved our knowledge of the identity and function of microorganisms in aquatic, terrestrial, and human biomes. Although open oceans have been the primary focus of studies on aquatic microbes, coastal and brackish ecosystems are now being surveyed. Here, we review so far published studies on microbes in the Baltic Sea, one of the world's largest brackish water bodies, using high throughput sequencing of environmental DNA and RNA. Collectively the data illustrate that Baltic Sea microbes are unique and highly diverse, and well adapted to this brackish-water ecosystem, findings that represent a novel base-line knowledge necessary for monitoring purposes and a sustainable management. More specifically, the data relate to environmental drivers for microbial community composition and function, assessments of the microbial biodiversity, adaptations and role of microbes in the nitrogen cycle, and microbial genome assembly from metagenomic sequences. With these discoveries as background, prospects of using metagenomics for Baltic Sea environmental monitoring are discussed.

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis.

BACKGROUND: Insertion sequences (ISs) are approximately 1 kbp long "jumping" genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I ("cut-and-paste") Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. RESULTS: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location ("local hopping") and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. CONCLUSIONS: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium 'Nostoc azollae' 0708 and in other obligate prokaryotic symbionts.

Functional tradeoffs underpin salinity-driven divergence in microbial community composition.

Bacterial community composition and functional potential change subtly across gradients in the surface ocean. In contrast, while there are significant phylogenetic divergences between communities from freshwater and marine habitats, the underlying mechanisms to this phylogenetic structuring yet remain unknown. We hypothesized that the functional potential of natural bacterial communities is linked to this striking divide between microbiomes. To test this hypothesis, metagenomic sequencing of microbial communities along a 1,800 km transect in the Baltic Sea area, encompassing a continuous natural salinity gradient from limnic to fully marine conditions, was explored. Multivariate statistical analyses showed that salinity is the main determinant of dramatic changes in microbial community composition, but also of large scale changes in core metabolic functions of bacteria. Strikingly, genetically and metabolically different pathways for key metabolic processes, such as respiration, biosynthesis of quinones and isoprenoids, glycolysis and osmolyte transport, were differentially abundant at high and low salinities. These shifts in functional capacities were observed at multiple taxonomic levels and within dominant bacterial phyla, while bacteria, such as SAR11, were able to adapt to the entire salinity gradient. We propose that the large differences in central metabolism required at high and low salinities dictate the striking divide between freshwater and marine microbiomes, and that the ability to inhabit different salinity regimes evolved early during bacterial phylogenetic differentiation. These findings significantly advance our understanding of microbial distributions and stress the need to incorporate salinity in future climate change models that predict increased levels of precipitation and a reduction in salinity.

Picocyanobacteria containing a novel pigment gene cluster dominate the brackish water Baltic Sea.

Photoautotrophic picocyanobacteria harvest light via phycobilisomes (PBS) consisting of the pigments phycocyanin (PC) and phycoerythrin (PE), encoded by genes in conserved gene clusters. The presence and arrangement of these gene clusters give picocyanobacteria characteristic light absorption properties and allow the colonization of specific ecological niches. To date, a full understanding of the evolution and distribution of the PBS gene cluster in picocyanobacteria has been hampered by the scarcity of genome sequences from fresh- and brackish water-adapted strains. To remediate this, we analysed genomes assembled from metagenomic samples collected along a natural salinity gradient, and over the course of a growth season, in the Baltic Sea. We found that while PBS gene clusters in picocyanobacteria sampled in marine habitats were highly similar to known references, brackish-adapted genotypes harboured a novel type not seen in previously sequenced genomes. Phylogenetic analyses showed that the novel gene cluster belonged to a clade of uncultivated picocyanobacteria that dominate the brackish Baltic Sea throughout the summer season, but are uncommon in other examined aquatic ecosystems. Further, our data suggest that the PE genes were lost in the ancestor of PC-containing coastal picocyanobacteria and that multiple horizontal gene transfer events have re-introduced PE genes into brackish-adapted strains, including the novel clade discovered here.

Multiple Modes of Cell Death Discovered in a Prokaryotic (Cyanobacterial) Endosymbiont.

Programmed cell death (PCD) is a genetically-based cell death mechanism with vital roles in eukaryotes. Although there is limited consensus on similar death mode programs in prokaryotes, emerging evidence suggest that PCD events are operative. Here we present cell death events in a cyanobacterium living endophytically in the fern Azolla microphylla, suggestive of PCD. This symbiosis is characterized by some unique traits such as a synchronized development, a vertical transfer of the cyanobacterium between plant generations, and a highly eroding cyanobacterial genome. A combination of methods was used to identify cell death modes in the cyanobacterium. Light- and electron microscopy analyses showed that the proportion of cells undergoing cell death peaked at 53.6% (average 20%) of the total cell population, depending on the cell type and host developmental stage. Biochemical markers used for early and late programmed cell death events related to apoptosis (Annexin V-EGFP and TUNEL staining assays), together with visualization of cytoskeleton alterations (FITC-phalloidin staining), showed that all cyanobacterial cell categories were affected by cell death. Transmission electron microscopy revealed four modes of cell death: apoptotic-like, autophagic-like, necrotic-like and autolytic-like. Abiotic stresses further enhanced cell death in a dose and time dependent manner. The data also suggest that dynamic changes in the peptidoglycan cell wall layer and in the cytoskeleton distribution patterns may act as markers for the various cell death modes. The presence of a metacaspase homolog (domain p20) further suggests that the death modes are genetically programmed. It is therefore concluded that multiple, likely genetically programmed, cell death modes exist in cyanobacteria, a finding that may be connected with the evolution of cell death in the plant kingdom.

Dinitrogen fixation is restricted to the terminal heterocysts in the invasive cyanobacterium Cylindrospermopsis raciborskii CS-505.

The toxin producing nitrogen-fixing heterocystous freshwater cyanobacterium Cylindrospermopsis raciborskii recently radiated from its endemic tropical environment into sub-tropical and temperate regions, a radiation likely to be favored by its ability to fix dinitrogen (diazotrophy). Although most heterocystous cyanobacteria differentiate regularly spaced intercalary heterocysts along their trichomes when combined nitrogen sources are depleted, C. raciborskii differentiates only two terminal heterocysts (one at each trichome end) that can reach >100 vegetative cells each. Here we investigated whether these terminal heterocysts are the exclusive sites for dinitrogen fixation in C. raciborskii. The highest nitrogenase activity and NifH biosynthesis (western-blot) were restricted to the light phase of a 12/12 light/dark cycle. Separation of heterocysts and vegetative cells (sonication and two-phase aqueous polymer partitioning) demonstrated that the terminal heterocysts are the sole sites for nifH expression (RT-PCR) and NifH biosynthesis. The latter finding was verified by the exclusive localization of nitrogenase in the terminal heterocysts of intact trichomes (immunogold-transmission electron microscopy and in situ immunofluorescence-light microscopy). These results suggest that the terminal heterocysts provide the combined nitrogen required by the often long trichomes (>100 vegetative cells). Our data also suggests that the terminal-heterocyst phenotype in C. raciborskii may be explained by the lack of a patL ortholog. These data help identify mechanisms by which C. raciborskii and other terminal heterocyst-forming cyanobacteria successfully inhabit environments depleted in combined nitrogen.

Trichodesmium--a widespread marine cyanobacterium with unusual nitrogen fixation properties.

The last several decades have witnessed dramatic advances in unfolding the diversity and commonality of oceanic diazotrophs and their N2 -fixing potential. More recently, substantial progress in diazotrophic cell biology has provided a wealth of information on processes and mechanisms involved. The substantial contribution by the diazotrophic cyanobacterial genus Trichodesmium to the nitrogen influx of the global marine ecosystem is by now undisputable and of paramount ecological importance, while the underlying cellular and molecular regulatory physiology has only recently started to unfold. Here, we explore and summarize current knowledge, related to the optimization of its diazotrophic capacity, from genomics to ecophysiological processes, via, for example, cellular differentiation (diazocytes) and temporal regulations, and suggest cellular research avenues that now ought to be explored.

Prokaryotic caspase homologs: phylogenetic patterns and functional characteristics reveal considerable diversity.

Caspases accomplish initiation and execution of apoptosis, a programmed cell death process specific to metazoans. The existence of prokaryotic caspase homologs, termed metacaspases, has been known for slightly more than a decade. Despite their potential connection to the evolution of programmed cell death in eukaryotes, the phylogenetic distribution and functions of these prokaryotic metacaspase sequences are largely uncharted, while a few experiments imply involvement in programmed cell death. Aiming at providing a more detailed picture of prokaryotic caspase homologs, we applied a computational approach based on Hidden Markov Model search profiles to identify and functionally characterize putative metacaspases in bacterial and archaeal genomes. Out of the total of 1463 analyzed genomes, merely 267 (18%) were identified to contain putative metacaspases, but their taxonomic distribution included most prokaryotic phyla and a few archaea (Euryarchaeota). Metacaspases were particularly abundant in Alphaproteobacteria, Deltaproteobacteria and Cyanobacteria, which harbor many morphologically and developmentally complex organisms, and a distinct correlation was found between abundance and phenotypic complexity in Cyanobacteria. Notably, Bacillus subtilis and Escherichia coli, known to undergo genetically regulated autolysis, lacked metacaspases. Pfam domain architecture analysis combined with operon identification revealed rich and varied configurations among the metacaspase sequences. These imply roles in programmed cell death, but also e.g. in signaling, various enzymatic activities and protein modification. Together our data show a wide and scattered distribution of caspase homologs in prokaryotes with structurally and functionally diverse sub-groups, and with a potentially intriguing evolutionary role. These features will help delineate future characterizations of death pathways in prokaryotes.

Genome fluctuations in cyanobacteria reflect evolutionary, developmental and adaptive traits.

BACKGROUND: Cyanobacteria belong to an ancient group of photosynthetic prokaryotes with pronounced variations in their cellular differentiation strategies, physiological capacities and choice of habitat. Sequencing efforts have shown that genomes within this phylum are equally diverse in terms of size and protein-coding capacity. To increase our understanding of genomic changes in the lineage, the genomes of 58 contemporary cyanobacteria were analysed for shared and unique orthologs. RESULTS: A total of 404 protein families, present in all cyanobacterial genomes, were identified. Two of these are unique to the phylum, corresponding to an AbrB family transcriptional regulator and a gene that escapes functional annotation although its genomic neighbourhood is conserved among the organisms examined. The evolution of cyanobacterial genome sizes involves a mix of gains and losses in the clade encompassing complex cyanobacteria, while a single event of reduction is evident in a clade dominated by unicellular cyanobacteria. Genome sizes and gene family copy numbers evolve at a higher rate in the former clade, and multi-copy genes were predominant in large genomes. Orthologs unique to cyanobacteria exhibiting specific characteristics, such as filament formation, heterocyst differentiation, diazotrophy and symbiotic competence, were also identified. An ancestral character reconstruction suggests that the most recent common ancestor of cyanobacteria had a genome size of approx. 4.5 Mbp and 1678 to 3291 protein-coding genes, 4%-6% of which are unique to cyanobacteria today. CONCLUSIONS: The different rates of genome-size evolution and multi-copy gene abundance suggest two routes of genome development in the history of cyanobacteria. The expansion strategy is driven by gene-family enlargment and generates a broad adaptive potential; while the genome streamlining strategy imposes adaptations to highly specific niches, also reflected in their different functional capacities. A few genomes display extreme proliferation of non-coding nucleotides which is likely to be the result of initial expansion of genomes/gene copy number to gain adaptive potential, followed by a shift to a life-style in a highly specific niche (e.g. symbiosis). This transition results in redundancy of genes and gene families, leading to an increase in junk DNA and eventually to gene loss. A few orthologs can be correlated with specific phenotypes in cyanobacteria, such as filament formation and symbiotic competence; these constitute exciting exploratory targets.

Genome erosion in a nitrogen-fixing vertically transmitted endosymbiotic multicellular cyanobacterium.

BACKGROUND: An ancient cyanobacterial incorporation into a eukaryotic organism led to the evolution of plastids (chloroplasts) and subsequently to the origin of the plant kingdom. The underlying mechanism and the identities of the partners in this monophyletic event remain elusive. METHODOLOGY/PRINCIPAL FINDINGS: To shed light on this evolutionary process, we sequenced the genome of a cyanobacterium residing extracellularly in an endosymbiosis with a plant, the water-fern Azolla filiculoides Lam. This symbiosis was selected as it has characters which make it unique among extant cyanobacterial plant symbioses: the cyanobacterium lacks autonomous growth and is vertically transmitted between plant generations. Our results reveal features of evolutionary significance. The genome is in an eroding state, evidenced by a large proportion of pseudogenes (31.2%) and a high frequency of transposable elements (approximately 600) scattered throughout the genome. Pseudogenization is found in genes such as the replication initiator dnaA and DNA repair genes, considered essential to free-living cyanobacteria. For some functional categories of genes pseudogenes are more prevalent than functional genes. Loss of function is apparent even within the 'core' gene categories of bacteria, such as genes involved in glycolysis and nutrient uptake. In contrast, serving as a critical source of nitrogen for the host, genes related to metabolic processes such as cell differentiation and nitrogen-fixation are well preserved. CONCLUSIONS/SIGNIFICANCE: This is the first finding of genome degradation in a plant symbiont and phenotypically complex cyanobacterium and one of only a few extracellular endosymbionts described showing signs of reductive genome evolution. Our findings suggest an ongoing selective streamlining of this cyanobacterial genome which has resulted in an organism devoted to nitrogen fixation and devoid of autonomous growth. The cyanobacterial symbiont of Azolla can thus be considered at the initial phase of a transition from free-living organism to a nitrogen-fixing plant entity, a transition process which may mimic what drove the evolution of chloroplasts from a cyanobacterial ancestor.