Genome-wide analysis of light sensing in Prochlorococcus

Emergence of metabolic coupling to the heterotroph Alteromonas promotes dark survival in Prochlorococcus

Prochlorococcus is found throughout the euphotic zone in the oligotrophic open ocean. Deep mixing and sinking while attached to particles can, however, transport Prochlorococcus cells below this sunlit zone, depriving them of light for extended periods of time. Previous work has shown that Prochlorococcus by itself cannot survive extended periods of darkness. However, when co-cultured with a heterotrophic microbe and subjected to repeated periods of extended darkness, Prochlorococcus cells develop an epigenetically inherited dark-tolerant phenotype that can survive longer periods of darkness. Here we examine the metabolic and physiological changes underlying this adaptation using co-cultures of dark-tolerant and parental strains of Prochlorococcus, each grown with the heterotroph Alteromonas under diel light:dark conditions. The relative abundance of Alteromonas was higher in dark-tolerant than parental co-cultures, while dark-tolerant Prochlorococcus cells were larger, contained less chlorophyll, and were less synchronized to the light:dark cycle. Meta-transcriptome analysis revealed that dark-tolerant co-cultures undergo a joint change, in which Prochlorococcus undergoes a relative shift from photosynthesis to respiration, while Alteromonas shifts toward using more organic acids instead of sugars. Furthermore, the transcriptome data suggested enhanced biosynthesis of amino acids and purines in dark-tolerant Prochlorococcus and enhanced degradation of these compounds in Alteromonas. Collectively, our results demonstrate that dark adaptation involves a strengthening of the metabolic coupling between Prochlorococcus and Alteromonas, presumably mediated by an enhanced, and compositionally modified, carbon exchange between the two species.

Transcriptional Mechanisms of Thermal Acclimation in Prochlorococcus

Low temperature limits the growth and the distribution of the key oceanic primary producer Prochlorococcus, which does not proliferate above a latitude of ca. 40°. Yet, the molecular basis of thermal acclimation in this cyanobacterium remains unexplored. We analyzed the transcriptional response of the Prochlorococcus marinus strain MIT9301 in long-term acclimations and in natural Prochlorococcus populations along a temperature range enabling its growth (17 to 30°C). MIT9301 upregulated mechanisms of the global stress response at the temperature minimum (17°C) but maintained the expression levels of genes involved in essential metabolic pathways (e.g., ATP synthesis and carbon fixation) along the whole thermal niche. Notably, the declining growth of MIT9301 from the optimum to the minimum temperature was coincident with a transcriptional suppression of the photosynthetic apparatus and a dampening of its circadian expression patterns, indicating a loss in their regulatory capacity under cold conditions. Under warm conditions, the cellular transcript inventory of MIT9301 was strongly streamlined, which may also induce regulatory imbalances due to stochasticity in gene expression. The daytime transcriptional suppression of photosynthetic genes at low temperature was also observed in metatranscriptomic reads mapping to MIT9301 across the global ocean, implying that this molecular mechanism may be associated with the restricted distribution of Prochlorococcus to temperate zones. IMPORTANCE Prochlorococcus is a major marine primary producer with a global impact on atmospheric CO₂ fixation. This cyanobacterium is widely distributed across the temperate ocean, but virtually absent at latitudes above 40° for yet unknown reasons. Temperature has been suggested as a major limiting factor, but the exact mechanisms behind Prochlorococcus thermal growth restriction remain unexplored. This study brings us closer to understanding how Prochlorococcus functions under challenging temperature conditions, by focusing on its transcriptional response after long-term acclimation from its optimum to its thermal thresholds. Our results show that the drop in Prochlorococcus growth rate under cold conditions was paralleled by a transcriptional suppression of the photosynthetic machinery during daytime and a loss in the organism's regulatory capacity to maintain circadian expression patterns. Notably, warm temperature induced a marked shrinkage of the organism's cellular transcript inventory, which may also induce regulatory imbalances in the future functioning of this cyanobacterium.

Basin-scale biogeography of marine phytoplankton reflects cellular-scale optimization of metabolism and physiology

Extensive microdiversity within Prochlorococcus, the most abundant marine cyanobacterium, occurs at scales from a single droplet of seawater to ocean basins. To interpret the structuring role of variations in genetic potential, as well as metabolic and physiological acclimation, we developed a mechanistic constraint-based modeling framework that incorporates the full suite of genes, proteins, metabolic reactions, pigments, and biochemical compositions of 69 sequenced isolates spanning the Prochlorococcus pangenome. Optimizing each strain to the local, observed physical and chemical environment along an Atlantic Ocean transect, we predicted variations in strain-specific patterns of growth rate, metabolic configuration, and physiological state, defining subtle niche subspaces directly attributable to differences in their encoded metabolic potential. Predicted growth rates covaried with observed ecotype abundances, affirming their significance as a measure of fitness and inferring a nonlinear density dependence of mortality. Our study demonstrates the potential to interpret global-scale ecosystem organization in terms of cellular-scale processes.

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.

A framework for the computational prediction and analysis of non-coding RNAs in microbial environmental populations and their experimental validation

Small regulatory RNAs and antisense RNAs play important roles in the regulation of gene expression in bacteria but are underexplored, especially in natural populations. While environmentally relevant microbes often are not amenable to genetic manipulation or cannot be cultivated in the laboratory, extensive metagenomic and metatranscriptomic datasets for these organisms might be available. Hence, dedicated workflows for specific analyses are needed to fully benefit from this information. Here, we identified abundant sRNAs from oceanic environmental populations of the ecologically important primary producer Prochlorococcus starting from a metatranscriptomic differential RNA-Seq (mdRNA-Seq) dataset. We tracked their homologs in laboratory isolates, and we provide a framework for their further detailed characterization. Several of the experimentally validated sRNAs responded to ecologically relevant changes in cultivation conditions. The expression of the here newly discovered sRNA Yfr28 was highly stimulated in low-nitrogen conditions. Its predicted top targets include mRNAs encoding cell division proteins, a sigma factor, and several enzymes and transporters, suggesting a pivotal role of Yfr28 in the coordination of primary metabolism and cell division. A cis-encoded antisense RNA was identified as a possible positive regulator of atpF encoding subunit b' of the ATP synthase complex. The presented workflow will also be useful for other environmentally relevant microorganisms for which experimental validation abilities are frequently limiting although there is wealth of sequence information available.

A minimum set of regulators to thrive in the ocean

Marine cyanobacteria of the genus Prochlorococcus thrive in high cell numbers throughout the euphotic zones of the world's subtropical and tropical oligotrophic oceans, making them some of the most ecologically relevant photosynthetic microorganisms on Earth. The ecological success of these free-living phototrophs suggests that they are equipped with a regulatory system competent to address many different stress situations. However, Prochlorococcus genomes are compact and streamlined, with the majority encoding only five different sigma factors, five to six two-component systems and eight types of other transcriptional regulators. Here, we summarize the existing information about the functions of these protein regulators, about transcriptomic responses to defined stress conditions, and discuss the current knowledge about riboswitches, RNA-based regulation and the roles of certain metabolites as co-regulators. We focus on the best-studied isolate, Prochlorococcus MED4, but extend to other strains and ecotypes when appropriate, and we include some information gained from metagenomic and metatranscriptomic analyses.

Interplay and Targetome of the Two Conserved Cyanobacterial sRNAs Yfr1 and Yfr2 in Prochlorococcus MED4

The sRNA Yfr1 and members of the Yfr2 sRNA family are almost universally present within cyanobacteria. The conserved motifs of these sRNAs are nearly complementary to each other, suggesting their ability to participate in crosstalk. The conserved motif of Yfr1 is shared by members of the Yfr10 sRNA family, members of which are otherwise less conserved in sequence, structure, and synteny compared to Yfr1. The different structural properties enable the discrimination of unique targets of Yfr1 and Yfr10. Unlike most studied regulatory sRNAs, Yfr1 gene expression only slightly changes under the tested stress conditions and is present at high levels at all times. In contrast, cellular levels of Yfr10 increase during the course of acclimation to darkness, and levels of Yfr2 increase when cells are shifted to high light or nitrogen limitation conditions. In this study, we investigated the targetomes of Yfr2, Yfr1, and Yfr10 in Prochlorococcus MED4, establishing CRAFD-Seq as a new method for identifying direct targets of these sRNAs that is applicable to all bacteria, including those that are not amenable to genetic modification. The results suggest that these sRNAs are integrated within a regulatory network of unprecedented complexity in the adjustment of carbon and nitrogen-related primary metabolism.

The GntR family transcriptional regulator PMM1637 regulates the highly conserved cyanobacterial sRNA Yfr2 in marine picocyanobacteria

Prochlorococcus is a marine picocyanobacterium with a streamlined genome that is adapted to different ecological niches in the oligotrophic oceans. There are currently >20 regulatory small RNAs (sRNAs) that have been identified in the model strain Prochlorococcus MED4. While most of these sRNAs are ecotype-specific, sRNA homologs of Yfr1 and of the Yfr2 family are widely found throughout the cyanobacterial phylum. Although they were identified 13 yrs ago, the functions of Yfr1 and Yfr2 have remained unknown. We observed a strong induction of two Yfr2 sRNA homologs of Prochlorococcus MED4 during high light stress and nitrogen starvation. Several Prochlorococcus and marine Synechococcus yfr2 promoter regions contain a conserved motif we named CGRE1 (cyanobacterial GntR family transcriptional regulator responsive element 1). Using the conserved promoter region as bait in a DNA affinity pull-down assay we identified the GntR family transcriptional regulator PMM1637 as a binding partner. Similar to Yfr2, homologs of PMM1637 are universally and exclusively found in cyanobacteria. We suggest that PMM1637 governs the induction of gene expression of Yfr2 homologs containing CGRE1 in their promoters under nitrogen-depleted and high-light stress conditions.

Heterotroph Interactions Alter Prochlorococcus Transcriptome Dynamics during Extended Periods of Darkness

Microbes evolve within complex ecological communities where biotic interactions impact both individual cells and the environment as a whole. Here we examine how cellular regulation in the marine cyanobacterium Prochlorococcus is influenced by a heterotrophic bacterium, Alteromonas macleodii, under different light conditions. We monitored the transcriptome of Prochlorococcus, grown either alone or in coculture, across a diel light:dark cycle and under the stress of extended darkness-a condition that cells would experience when mixed below the ocean's euphotic zone. More Prochlorococcus transcripts exhibited 24-h periodic oscillations in coculture than in pure culture, both over the normal diel cycle and after the shift to extended darkness. This demonstrates that biotic interactions, and not just light, can affect timing mechanisms in Prochlorococcus, which lacks a self-sustaining circadian oscillator. The transcriptomes of replicate pure cultures of Prochlorococcus lost their synchrony within 5 h of extended darkness and reflected changes in stress responses and metabolic functions consistent with growth cessation. In contrast, when grown with Alteromonas, replicate Prochlorococcus transcriptomes tracked each other for at least 13 h in the dark and showed signs of continued biosynthetic and metabolic activity. The transcriptome patterns suggest that the heterotroph may be providing energy or essential biosynthetic substrates to Prochlorococcus in the form of organic compounds, sustaining this autotroph when it is deprived of solar energy. Our findings reveal conditions where mixotrophic metabolism may benefit marine cyanobacteria and highlight new impacts of community interactions on basic Prochlorococcus cellular processes. IMPORTANCEProchlorococcus is the most abundant photosynthetic organism on the planet. These cells play a central role in the physiology of surrounding heterotrophs by supplying them with fixed organic carbon. It is becoming increasingly clear, however, that interactions with heterotrophs can affect autotrophs as well. Here we show that such interactions have a marked impact on the response of Prochlorococcus to the stress of extended periods of darkness, as reflected in transcriptional dynamics. These data suggest that diel transcriptional rhythms within Prochlorococcus, which are generally considered to be strictly under the control of light quantity, quality, and timing, can also be influenced by biotic interactions. Together, these findings provide new insights into the importance of microbial interactions on Prochlorococcus physiology and reveal conditions where heterotroph-derived compounds may support autotrophs-contrary to the canonical autotroph-to-heterotroph trophic paradigm.

Impact of light intensity and quality on chromatophore and nuclear gene expression in Paulinella chromatophora, an amoeba with nascent photosynthetic organelles

Plastid evolution has been attributed to a single primary endosymbiotic event that occurred about 1.6 billion years ago (BYA) in which a cyanobacterium was engulfed and retained by a eukaryotic cell, although early steps in plastid integration are poorly understood. The photosynthetic amoeba Paulinella chromatophora represents a unique model for the study of plastid evolution because it contains cyanobacterium-derived photosynthetic organelles termed 'chromatophores' that originated relatively recently (0.09-0.14 BYA). The chromatophore genome is about a third the size of the genome of closely related cyanobacteria, but 10-fold larger than most plastid genomes. Several genes have been transferred from the chromatophore genome to the host nuclear genome through endosymbiotic gene transfer (EGT). Some EGT-derived proteins could be imported into chromatophores for function. Two photosynthesis-related genes (psaI and csos4A) are encoded by both the nuclear and chromatophore genomes, suggesting that EGT in Paulinella chromatophora is ongoing. Many EGT-derived genes encode proteins that function in photosynthesis and photoprotection, including an expanded family of high-light-inducible (ncHLI) proteins. Cyanobacterial hli genes are high-light induced and required for cell viability under excess light. We examined the impact of light on Paulinella chromatophora and found that this organism is light sensitive and lacks light-induced transcriptional regulation of chromatophore genes and most EGT-derived nuclear genes. However, several ncHLI genes have reestablished light-dependent regulation, which appears analogous to what is observed in cyanobacteria. We postulate that expansion of the ncHLI gene family and its regulation may reflect the light/oxidative stress experienced by Paulinella chromatophora as a consequence of the as yet incomplete integration of host and chromatophore metabolisms.

Transcriptional response of Prochlorococcus to co-culture with a marine Alteromonas: differences between strains and the involvement of putative infochemicals

Interactions between marine microorganisms may determine the dynamics of microbial communities. Here, we show that two strains of the globally abundant marine cyanobacterium Prochlorococcus, MED4 and MIT9313, which belong to two different ecotypes, differ markedly in their response to co-culture with a marine heterotrophic bacterium, Alteromonas macleodii strain HOT1A3. HOT1A3 enhanced the growth of MIT9313 at low cell densities, yet inhibited it at a higher concentration, whereas it had no effect on MED4 growth. The early transcriptomic responses of Prochlorococcus cells after 20 h in co-culture showed no evidence of nutrient starvation, whereas the expression of genes involved in photosynthesis, protein synthesis and stress responses typically decreased in MED4 and increased in MIT313. Differential expression of genes involved in outer membrane modification, efflux transporters and, in MIT9313, lanthipeptides (prochlorosins) suggests that Prochlorococcus mount a specific response to the presence of the heterotroph in the cultures. Intriguingly, many of the differentially-expressed genes encoded short proteins, including two new families of co-culture responsive genes: CCRG-1, which is found across the Prochlorococcus lineage and CCRG-2, which contains a sequence motif involved in the export of prochlorosins and other bacteriocin-like peptides, and are indeed released from the cells into the media.

Gene Expression Patterns during Light and Dark Infection of Prochlorococcus by Cyanophage

Cyanophage infecting the marine cyanobacteria Prochlorococcus and Synechococcus require light and host photosystem activity for optimal reproduction. Many cyanophages encode multiple photosynthetic electron transport (PET) proteins, which are presumed to maintain electron flow and produce ATP and NADPH for nucleotide biosynthesis and phage genome replication. However, evidence suggests phage augment NADPH production via the pentose phosphate pathway (PPP), thus calling into question the need for NADPH production by PET. Genes implicated in cyclic PET have since been identified in cyanophage genomes. It remains an open question which mode of PET, cyclic or linear, predominates in infected cyanobacteria, and thus whether the balance is towards producing ATP or NADPH. We sequenced transcriptomes of a cyanophage (P-HM2) and its host (Prochlorococcus MED4) throughout infection in the light or in the dark, and analyzed these data in the context of phage replication and metabolite measurements. Infection was robust in the light, but phage were not produced in the dark. Host gene transcripts encoding high-light inducible proteins and two terminal oxidases (plastoquinol terminal oxidase and cytochrome c oxidase)-implicated in protecting the photosynthetic membrane from light stress-were the most enriched in light but not dark infection. Among the most diminished transcripts in both light and dark infection was ferredoxin-NADP+ reductase (FNR), which uses the electron acceptor NADP+ to generate NADPH in linear photosynthesis. The phage gene for CP12, which putatively inhibits the Calvin cycle enzyme that receives NADPH from FNR, was highly expressed in light infection. Therefore, both PET production of NADPH and its consumption by carbon fixation are putatively repressed during phage infection in light. Transcriptomic evidence is thus consistent with cyclic photophosphorylation using oxygen as the terminal electron acceptor as the dominant mode of PET under infection, with ATP from PET and NADPH from the PPP producing the energy and reducing equivalents for phage nucleotide biosynthesis and replication.

Torn apart and reunited: impact of a heterotroph on the transcriptome of Prochlorococcus

Microbial interactions, whether direct or indirect, profoundly affect the physiology of individual cells and ultimately have the potential to shape the biogeochemistry of the Earth. For example, the growth of Prochlorococcus, the numerically dominant cyanobacterium in the oceans, can be improved by the activity of co-occurring heterotrophs. This effect has been largely attributed to the role of heterotrophs in detoxifying reactive oxygen species that Prochlorococcus, which lacks catalase, cannot. Here, we explore this phenomenon further by examining how the entire transcriptome of Prochlorococcus NATL2A changes in the presence of a naturally co-occurring heterotroph, Alteromonas macleodii MIT1002, with which it was co-cultured for years, separated and then reunited. Significant changes in the Prochlorococcus transcriptome were evident within 6 h of initiating co-culture, with groups of transcripts changing in different temporal waves. Many transcriptional changes persisted throughout the 48 h experiment, suggesting that the presence of the heterotroph affected a stable shift in Prochlorococcus physiology. These initial transcriptome changes largely corresponded to reduced stress conditions for Prochlorococcus, as inferred from the depletion of transcripts encoding DNA repair enzymes and many members of the 'high light inducible' family of stress-response proteins. Later, notable changes were seen in transcripts encoding components of the photosynthetic apparatus (particularly, an increase in PSI subunits and chlorophyll synthesis enzymes), ribosomal proteins and biosynthetic enzymes, suggesting that the introduction of the heterotroph may have induced increased production of reduced carbon compounds for export. Changes in secretion-related proteins and transporters also highlight the potential for metabolic exchange between the two strains.

Shedding new light on viral photosynthesis

Viruses infecting the environmentally important marine cyanobacteria Prochlorococcus and Synechococcus encode 'auxiliary metabolic genes' (AMGs) involved in the light and dark reactions of photosynthesis. Here, we discuss progress on the inventory of such AMGs in the ever-increasing number of viral genome sequences as well as in metagenomic datasets. We contextualise these gene acquisitions with reference to a hypothesised fitness gain to the phage. We also report new evidence with regard to the sequence and predicted structural properties of viral petE genes encoding the soluble electron carrier plastocyanin. Viral copies of PetE exhibit extensive modifications to the N-terminal signal peptide and possess several novel residues in a region responsible for interaction with redox partners. We also highlight potential knowledge gaps in this field and discuss future opportunities to discover novel phage-host interactions involved in the photosynthetic process.

Response of Prochlorococcus to varying CO2:O2 ratios

Carbon fixation has a central role in determining cellular redox poise, increasingly understood to be a key parameter in cyanobacterial physiology. In the cyanobacterium Prochlorococcus-the most abundant phototroph in the oligotrophic oceans-the carbon-concentrating mechanism is reduced to the bare essentials. Given the ability of Prochlorococcus populations to grow under a wide range of oxygen concentrations in the ocean, we wondered how carbon and oxygen physiology intersect in this minimal phototroph. Thus, we examined how CO2:O2 gas balance influenced growth and chlorophyll fluorescence in Prochlorococcus strain MED4. Under O2 limitation, per-cell chlorophyll fluorescence fell at all CO2 levels, but still permitted substantial growth at moderate and high CO2. Under CO2 limitation, we observed little growth at any O2 level, although per-cell chlorophyll fluorescence fell less sharply when O2 was available. We explored this pattern further by monitoring genome-wide transcription in cells shocked with acute limitation of CO2, O2 or both. O2 limitation produced much smaller transcriptional changes than the broad suppression seen under CO2 limitation and CO2/O2 co-limitation. Strikingly, both CO2 limitation conditions initially evoked a transcriptional response that resembled the pattern previously seen in high-light stress, but at later timepoints we observed O2-dependent recovery of photosynthesis-related transcripts. These results suggest that oxygen has a protective role in Prochlorococcus when carbon fixation is not a sufficient sink for light energy.

Clade-Specific Quantitative Analysis of Photosynthetic Gene Expression in Prochlorococcus

Newly designed primers targeting rbcL (CO2 fixation), psbA (photosystem II) and rnpB (reference) genes were used in qRT-PCR assays to assess the photosynthetic capability of natural communities of Prochlorococcus, the most abundant photosynthetic organism on Earth and a major contributor to primary production in oligotrophic oceans. After optimizing sample collection methodology, we analyzed a total of 62 stations from the Malaspina 2010 circumnavigation (including Atlantic, Pacific and Indian Oceans) at three different depths. Sequence and quantitative analyses of the corresponding amplicons showed the presence of high-light (HL) and low-light (LL) Prochlorococcus clades in essentially all 182 samples, with a largely uniform stratification of LL and HL sequences. Synechococcus cross-amplifications were detected by the taxon-specific melting temperatures of the amplicons. Laboratory exposure of Prochlorococcus MED4 (HL) and MIT9313 (LL) strains to organic pollutants (PAHs and organochlorine compounds) showed a decrease of rbcL transcript abundances, and of the rbcL to psbA ratios for both strains. We propose this technique as a convenient assay to evaluate effects of environmental stressors, including pollution, on the oceanic Prochlorococcus photosynthetic function.

Dealing with light: the widespread and multitasking cryptochrome/photolyase family in photosynthetic organisms

Light is essential for the life of photosynthetic organisms as it is a source of energy and information from the environment. Light excess or limitation can be a cause of stress however. Photosynthetic organisms exhibit sophisticated mechanisms to adjust their physiology and growth to the local environmental light conditions. The cryptochrome/photolyase family (CPF) is composed of flavoproteins with similar structures that display a variety of light-dependent functions. This family encompasses photolyases, blue-light activated enzymes that repair ultraviolet-light induced DNA damage, and cryptochromes, known for their photoreceptor functions in terrestrial plants. For this review, we searched extensively for CPFs in the available genome databases to trace the distribution and evolution of this protein family in photosynthetic organisms. By merging molecular data with current knowledge from the functional characterization of CPFs from terrestrial and aquatic organisms, we discuss their roles in (i) photoperception, (ii) biological rhythm regulation and (iii) light-induced stress responses. We also explore their possible implication in light-related physiological acclimation and their distribution in phototrophs living in different environments. The outcome of this structure-function analysis reconstructs the complex scenarios in which CPFs have evolved, as highlighted by the novel functions and biochemical properties of the most recently described family members in algae.

Dataset for metatranscriptome analysis of Prochlorococcus-rich marine picoplankton communities in the Gulf of Aqaba, Red Sea

Regulatory RNAs play a central role in the regulation of gene expression and can act on several regulatory levels from transcriptional initiation and RNA processing to the control of initiation of translation and RNA stability. One class of these molecules is non-coding (nc)RNAs in bacteria that typically lack protein-coding potential, range in size between 50 and 500nt and originate from intergenic regions. Common methods for the identification of these RNAs are either based on computational predictions, or on transcriptomic analyses of laboratory cultures, whereas very little is known about ncRNAs in environmental microbial populations. Here, we have combined a metatranscriptomics approach with a selective enrichment protocol for ncRNAs. The primary objective of this study was the identification of novel, environmentally relevant ncRNAs focusing on the cyanobacterium Prochlorococcus, which was one of the dominant microorganisms of the marine community of the Gulf of Aqaba when samples were taken.

Comparative transcriptomics of two environmentally relevant cyanobacteria reveals unexpected transcriptome diversity

Prochlorococcus is a genus of abundant and ecologically important marine cyanobacteria. Here, we present a comprehensive comparison of the structure and composition of the transcriptomes of two Prochlorococcus strains, which, despite their similarities, have adapted their gene pool to specific environmental constraints. We present genome-wide maps of transcriptional start sites (TSS) for both organisms, which are representatives of the two most diverse clades within the two major ecotypes adapted to high- and low-light conditions, respectively. Our data suggest antisense transcription for three-quarters of all genes, which is substantially more than that observed in other bacteria. We discovered hundreds of TSS within genes, most notably within 16 of the 29 prochlorosin genes, in strain MIT9313. A direct comparison revealed very little conservation in the location of TSS and the nature of non-coding transcripts between both strains. We detected extremely short 5' untranslated regions with a median length of only 27 and 29 nt for MED4 and MIT9313, respectively, and for 8% of all protein-coding genes the median distance to the start codon is only 10 nt or even shorter. These findings and the absence of an obvious Shine-Dalgarno motif suggest that leaderless translation and ribosomal protein S1-dependent translation constitute alternative mechanisms for translation initiation in Prochlorococcus. We conclude that genome-wide antisense transcription is a major component of the transcriptional output from these relatively small genomes and that a hitherto unrecognized high degree of complexity and variability of gene expression exists in their transcriptional architecture.

Impact of Solar Radiation on Gene Expression in Bacteria

Microorganisms often regulate their gene expression at the level of transcription and/or translation in response to solar radiation. In this review, we present the use of both transcriptomics and proteomics to advance knowledge in the field of bacterial response to damaging radiation. Those studies pertain to diverse application areas such as fundamental microbiology, water treatment, microbial ecology and astrobiology. Even though it has been demonstrated that mRNA abundance is not always consistent with the protein regulation, we present here an exhaustive review on how bacteria regulate their gene expression at both transcription and translation levels to enable biomarkers identification and comparison of gene regulation from one bacterial species to another.

Genetic diversity in cultured and wild marine cyanomyoviruses reveals phosphorus stress as a strong selective agent

Viruses that infect marine cyanobacteria-cyanophages-often carry genes with orthologs in their cyanobacterial hosts, and the frequency of these genes can vary with habitat. To explore habitat-influenced genomic diversity more deeply, we used the genomes of 28 cultured cyanomyoviruses as references to identify phage genes in three ocean habitats. Only about 6-11% of genes were consistently observed in the wild, revealing high gene-content variability in these populations. Numerous shared phage/host genes differed in relative frequency between environments, including genes related to phosphorous acquisition, photorespiration, photosynthesis and the pentose phosphate pathway, possibly reflecting environmental selection for these genes in cyanomyovirus genomes. The strongest emergent signal was related to phosphorous availability; a higher fraction of genomes from relatively low-phosphorus environments-the Sargasso and Mediterranean Sea-contained host-like phosphorus assimilation genes compared with those from the N. Pacific Gyre. These genes are known to be upregulated when the host is phosphorous starved, a response mediated by pho box motifs in phage genomes that bind a host regulatory protein. Eleven cyanomyoviruses have predicted pho boxes upstream of the phosphate-acquisition genes pstS and phoA; eight of these have a conserved cyanophage-specific gene (PhCOG173) between the pho box and pstS. PhCOG173 is also found upstream of other shared phage/host genes, suggesting a unique regulatory role. Pho boxes are found upstream of high light-inducible (hli) genes in cyanomyoviruses, suggesting that this motif may have a broader role than regulating phosphorous-stress responses in infected hosts or that these hlis are involved in the phosphorous-stress response.

A global analysis of adaptive evolution of operons in cyanobacteria

Operons are an important feature of prokaryotic genomes. Evolution of operons is hypothesized to be adaptive and has contributed significantly towards coordinated optimization of functions. Two conflicting theories, based on (i) in situ formation to achieve co-regulation and (ii) horizontal gene transfer of functionally linked gene clusters, are generally considered to explain why and how operons have evolved. Furthermore, effects of operon evolution on genomic traits such as intergenic spacing, operon size and co-regulation are relatively less explored. Based on the conservation level in a set of diverse prokaryotes, we categorize the operonic gene pair associations and in turn the operons as ancient and recently formed. This allowed us to perform a detailed analysis of operonic structure in cyanobacteria, a morphologically and physiologically diverse group of photoautotrophs. Clustering based on operon conservation showed significant similarity with the 16S rRNA-based phylogeny, which groups the cyanobacterial strains into three clades. Clade C, dominated by strains that are believed to have undergone genome reduction, shows a larger fraction of operonic genes that are tightly packed in larger sized operons. Ancient operons are in general larger, more tightly packed, better optimized for co-regulation and part of key cellular processes. A sub-clade within Clade B, which includes Synechocystis sp. PCC 6803, shows a reverse trend in intergenic spacing. Our results suggest that while in situ formation and vertical descent may be a dominant mechanism of operon evolution in cyanobacteria, optimization of intergenic spacing and co-regulation are part of an ongoing process in the life-cycle of operons.

Acclimation to high-light conditions in cyanobacteria: from gene expression to physiological responses

Photosynthetic organisms have evolved various acclimatory responses to high-light (HL) conditions to maintain a balance between energy supply (light harvesting and electron transport) and consumption (cellular metabolism) and to protect the photosynthetic apparatus from photodamage. The molecular mechanism of HL acclimation has been extensively studied in the unicellular cyanobacterium Synechocystis sp. PCC 6803. Whole genome DNA microarray analyses have revealed that the change in gene expression profile under HL is closely correlated with subsequent acclimatory responses such as (1) acceleration in the rate of photosystem II turnover, (2) downregulation of light harvesting capacity, (3) development of a protection mechanism for the photosystems against excess light energy, (4) upregulation of general protection mechanism components, and (5) regulation of carbon and nitrogen assimilation. In this review article, we survey recent progress in the understanding of the molecular mechanisms of these acclimatory responses in Synechocystis sp. PCC 6803. We also briefly describe attempts to understand HL acclimation in various cyanobacterial species in their natural environments.

ProPortal: a resource for integrated systems biology of Prochlorococcus and its phage

ProPortal (http://proportal.mit.edu/) is a database containing genomic, metagenomic, transcriptomic and field data for the marine cyanobacterium Prochlorococcus. Our goal is to provide a source of cross-referenced data across multiple scales of biological organization—from the genome to the ecosystem—embracing the full diversity of ecotypic variation within this microbial taxon, its sister group, Synechococcus and phage that infect them. The site currently contains the genomes of 13 Prochlorococcus strains, 11 Synechococcus strains and 28 cyanophage strains that infect one or both groups. Cyanobacterial and cyanophage genes are clustered into orthologous groups that can be accessed by keyword search or through a genome browser. Users can also identify orthologous gene clusters shared by cyanobacterial and cyanophage genomes. Gene expression data for Prochlorococcus ecotypes MED4 and MIT9313 allow users to identify genes that are up or downregulated in response to environmental stressors. In addition, the transcriptome in synchronized cells grown on a 24-h light–dark cycle reveals the choreography of gene expression in cells in a ‘natural’ state. Metagenomic sequences from the Global Ocean Survey from Prochlorococcus, Synechococcus and phage genomes are archived so users can examine the differences between populations from diverse habitats. Finally, an example of cyanobacterial population data from the field is included.

Transcriptome response of high- and low-light-adapted Prochlorococcus strains to changing iron availability

Prochlorococcus contributes significantly to ocean primary productivity. The link between primary productivity and iron in specific ocean regions is well established and iron limitation of Prochlorococcus cell division rates in these regions has been shown. However, the extent of ecotypic variation in iron metabolism among Prochlorococcus and the molecular basis for differences is not understood. Here, we examine the growth and transcriptional response of Prochlorococcus strains, MED4 and MIT9313, to changing iron concentrations. During steady state, MIT9313 sustains growth at an order-of-magnitude lower iron concentration than MED4. To explore this difference, we measured the whole-genome transcriptional response of each strain to abrupt iron starvation and rescue. Only four of the 1159 orthologs of MED4 and MIT9313 were differentially expressed in response to iron in both strains. However, in each strain, the expression of over a hundred additional genes changed, many of which are in labile genomic regions, suggesting a role for lateral gene transfer in establishing diversity of iron metabolism among Prochlorococcus. Furthermore, we found that MED4 lacks three genes near the iron-deficiency-induced gene (idiA) that are present and induced by iron stress in MIT9313. These genes are interesting targets for studying the adaptation of natural Prochlorococcus assemblages to local iron conditions as they show more diversity than other genomic regions in environmental metagenomic databases.

Antisense RNA protects mRNA from RNase E degradation by RNA-RNA duplex formation during phage infection

The ecologically important cyanobacterium Prochlorococcus possesses the smallest genome among oxyphototrophs, with a reduced suite of protein regulators and a disproportionately high number of regulatory RNAs. Many of these are asRNAs, raising the question whether they modulate gene expression through the protection of mRNA from RNase E degradation. To address this question, we produced recombinant RNase E from Prochlorococcus sp. MED4, which functions optimally at 12 mM Mg²⁺, pH 9 and 35°C. RNase E cleavage assays were performed with this recombinant protein to assess enzyme activity in the presence of single- or double-stranded RNA substrates. We found that extraordinarily long asRNAs of 3.5 and 7 kb protect a set of mRNAs from RNase E degradation that accumulate during phage infection. These asRNA–mRNA duplex formations mask single-stranded recognition sites of RNase E, leading to increased stability of the mRNAs. Such interactions directly modulate RNA stability and provide an explanation for enhanced transcript abundance of certain mRNAs during phage infection. Protection from RNase E-triggered RNA decay may constitute a hitherto unknown regulatory function of bacterial cis-asRNAs, impacting gene expression.

Learning to read the oceans genomics of marine phytoplankton

The phytoplankton are key members of marine ecosystems, generating about half of global primary productivity, supporting valuable fisheries and regulating global biogeochemical cycles. Marine phytoplankton are phylogenetically diverse and are comprised of both prokaryotic and eukaryotic species. In the last decade, new insights have been gained into the ecology and evolution of these important organisms through whole genome sequencing projects and more recently, through both transcriptomics and targeted metagenomics approaches. Sequenced genomes of cyanobacteria are generally small, ranging in size from 1.8 to 9 million base pairs (Mbp). Eukaryotic genomes, in general, have a much larger size range and those that have been sequenced range from 12 to 57 Mbp. Whole genome sequencing projects have revealed key features of the evolutionary history of marine phytoplankton, their varied responses to environmental stress, their ability to scavenge and store nutrients and their unique ability to form elaborate cellular coverings. We have begun to learn how to read the 'language' of marine phytoplankton, as written in their DNA. Here, we review the ecological and evolutionary insights gained from whole genome sequencing projects, illustrate how these genomes are yielding information on marine natural products and informing nanotechnology as well as make suggestions for future directions in the field of marine phytoplankton genomics.

Ultraviolet stress delays chromosome replication in light/dark synchronized cells of the marine cyanobacterium Prochlorococcus marinus PCC9511

BACKGROUND

The marine cyanobacterium Prochlorococcus is very abundant in warm, nutrient-poor oceanic areas. The upper mixed layer of oceans is populated by high light-adapted Prochlorococcus ecotypes, which despite their tiny genome (approximately 1.7 Mb) seem to have developed efficient strategies to cope with stressful levels of photosynthetically active and ultraviolet (UV) radiation. At a molecular level, little is known yet about how such minimalist microorganisms manage to sustain high growth rates and avoid potentially detrimental, UV-induced mutations to their DNA. To address this question, we studied the cell cycle dynamics of P. marinus PCC9511 cells grown under high fluxes of visible light in the presence or absence of UV radiation. Near natural light-dark cycles of both light sources were obtained using a custom-designed illumination system (cyclostat). Expression patterns of key DNA synthesis and repair, cell division, and clock genes were analyzed in order to decipher molecular mechanisms of adaptation to UV radiation.

RESULTS

The cell cycle of P. marinus PCC9511 was strongly synchronized by the day-night cycle. The most conspicuous response of cells to UV radiation was a delay in chromosome replication, with a peak of DNA synthesis shifted about 2 h into the dark period. This delay was seemingly linked to a strong downregulation of genes governing DNA replication (dnaA) and cell division (ftsZ, sepF), whereas most genes involved in DNA repair (such as recA, phrA, uvrA, ruvC, umuC) were already activated under high visible light and their expression levels were only slightly affected by additional UV exposure.

CONCLUSIONS

Prochlorococcus cells modified the timing of the S phase in response to UV exposure, therefore reducing the risk that mutations would occur during this particularly sensitive stage of the cell cycle. We identified several possible explanations for the observed timeshift. Among these, the sharp decrease in transcript levels of the dnaA gene, encoding the DNA replication initiator protein, is sufficient by itself to explain this response, since DNA synthesis starts only when the cellular concentration of DnaA reaches a critical threshold. However, the observed response likely results from a more complex combination of UV-altered biological processes.

Short RNA half-lives in the slow-growing marine cyanobacterium Prochlorococcus

BACKGROUND

RNA turnover plays an important role in the gene regulation of microorganisms and influences their speed of acclimation to environmental changes. We investigated whole-genome RNA stability of Prochlorococcus, a relatively slow-growing marine cyanobacterium doubling approximately once a day, which is extremely abundant in the oceans.

RESULTS

Using a combination of microarrays, quantitative RT-PCR and a new fitting method for determining RNA decay rates, we found a median half-life of 2.4 minutes and a median decay rate of 2.6 minutes for expressed genes - twofold faster than that reported for any organism. The shortest transcript half-life (33 seconds) was for a gene of unknown function, while some of the longest (approximately 18 minutes) were for genes with high transcript levels. Genes organized in operons displayed intriguing mRNA decay patterns, such as increased stability, and delayed onset of decay with greater distance from the transcriptional start site. The same phenomenon was observed on a single probe resolution for genes greater than 2 kb.

CONCLUSIONS

We hypothesize that the fast turnover relative to the slow generation time in Prochlorococcus may enable a swift response to environmental changes through rapid recycling of nucleotides, which could be advantageous in nutrient poor oceans. Our growing understanding of RNA half-lives will help us interpret the growing bank of metatranscriptomic studies of wild populations of Prochlorococcus. The surprisingly complex decay patterns of large transcripts reported here, and the method developed to describe them, will open new avenues for the investigation and understanding of RNA decay for all organisms.

Spatial patterns and light-driven variation of microbial population gene expression in surface waters of the oligotrophic open ocean

Because bacterioplankton production rates do not vary strongly across vast expanses of the ocean, it is unclear how variability in community structure corresponds with functional variability in the open ocean. We surveyed community transcript functional profiles at eight locations in the open ocean, in both the light and in the dark, using the genomic subsystems approach, to understand variability in gene expression patterns in surface waters. Metatranscriptomes from geographically distinct areas and collected during the day and night shared a large proportion of metabolic functional similarity (74%) at the finest metabolic resolution possible. The variability between metatranscriptomes could be explained by phylogenetic differences between libraries (Mantel test, P < 0.0001). Several key gene expression pathways, including Photosystem I, Photosystem II and ammonium uptake, demonstrated the most variability both geographically and between light and dark. Libraries were dominated by transcripts of the cyanobacterium Prochlorocococcus marinus, where most geographical and diel variability between metatranscriptomes reflected between-station differences in cyanobacterial phototrophic metabolism. Our results demonstrate that active genetic machinery in surface waters of the ocean is dominated by photosynthetic microorganisms and their site-to-site variability, while variability in the remainder of assemblages is dependent on local taxonomic composition.

Prochlorococcus: advantages and limits of minimalism

Prochlorococcus is the key phytoplanktonic organism of tropical gyres, large ocean regions that are depleted of the essential macronutrients needed for photosynthesis and cell growth. This cyanobacterium has adapted itself to oligotrophy by minimizing the resources necessary for life through a drastic reduction of cell and genome sizes. This rarely observed strategy in free-living organisms has conferred on Prochlorococcus a considerable advantage over other phototrophs, including its closest relative Synechococcus, for life in this vast yet little variable ecosystem. However, this strategy seems to reach its limits in the upper layer of the S Pacific gyre, the most oligotrophic region of the world ocean. By losing some important genes and/or functions during evolution, Prochlorococcus has seemingly become dependent on co-occurring microorganisms. In this review, we present some of the recent advances in the ecology, biology, and evolution of Prochlorococcus, which because of its ecological importance and tiny genome is rapidly imposing itself as a model organism in environmental microbiology.

Stress responses in Prochlorococcus MIT9313 vs. SS120 involve differential expression of genes encoding proteases ClpP, FtsH and Lon

Prochlorococcus is a marine cyanobacterium responsible for a significant part of global primary production as well as being one of the most abundant organisms on Earth. Protein turnover is an essential and poorly understood aspect of the cyanobacterial response to environmental stresses. In the present work, cultures of the SS120 and MIT9313 strains were subjected to several conditions, and quantitative real time RT-PCR was used to measure changes in the expression of genes encoding three representative ATP-dependent proteases. We found common responses to conditions such as aging. However, the expression pattern under nutrient starvation was strikingly different in the two strains, probably reflecting the different regulatory backgrounds of the two ecotypes here studied.

Ecological genomics of marine picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus numerically dominate the picophytoplankton of the world ocean, making a key contribution to global primary production. Prochlorococcus was isolated around 20 years ago and is probably the most abundant photosynthetic organism on Earth. The genus comprises specific ecotypes which are phylogenetically distinct and differ markedly in their photophysiology, allowing growth over a broad range of light and nutrient conditions within the 45 degrees N to 40 degrees S latitudinal belt that they occupy. Synechococcus and Prochlorococcus are closely related, together forming a discrete picophytoplankton clade, but are distinguishable by their possession of dissimilar light-harvesting apparatuses and differences in cell size and elemental composition. Synechococcus strains have a ubiquitous oceanic distribution compared to that of Prochlorococcus strains and are characterized by phylogenetically discrete lineages with a wide range of pigmentation. In this review, we put our current knowledge of marine picocyanobacterial genomics into an environmental context and present previously unpublished genomic information arising from extensive genomic comparisons in order to provide insights into the adaptations of these marine microbes to their environment and how they are reflected at the genomic level.

Choreography of the transcriptome, photophysiology, and cell cycle of a minimal photoautotroph, prochlorococcus

The marine cyanobacterium Prochlorococcus MED4 has the smallest genome and cell size of all known photosynthetic organisms. Like all phototrophs at temperate latitudes, it experiences predictable daily variation in available light energy which leads to temporal regulation and partitioning of key cellular processes. To better understand the tempo and choreography of this minimal phototroph, we studied the entire transcriptome of the cell over a simulated daily light-dark cycle, and placed it in the context of diagnostic physiological and cell cycle parameters. All cells in the culture progressed through their cell cycles in synchrony, thus ensuring that our measurements reflected the behavior of individual cells. Ninety percent of the annotated genes were expressed, and 80% had cyclic expression over the diel cycle. For most genes, expression peaked near sunrise or sunset, although more subtle phasing of gene expression was also evident. Periodicities of the transcripts of genes involved in physiological processes such as in cell cycle progression, photosynthesis, and phosphorus metabolism tracked the timing of these activities relative to the light-dark cycle. Furthermore, the transitions between photosynthesis during the day and catabolic consumption of energy reserves at night— metabolic processes that share some of the same enzymes — appear to be tightly choreographed at the level of RNA expression. In-depth investigation of these patterns identified potential regulatory proteins involved in balancing these opposing pathways. Finally, while this analysis has not helped resolve how a cell with so little regulatory capacity, and a ‘deficient’ circadian mechanism, aligns its cell cycle and metabolism so tightly to a light-dark cycle, it does provide us with a valuable framework upon which to build when the Prochlorococcus proteome and metabolome become available.

Carrying photosynthesis genes increases ecological fitness of cyanophage in silico

Several viruses infecting marine cyanobacteria carry photosynthesis genes (e.g. psbA, hli) that are expressed, yield proteins (D1, HLIP) and help maintain the cell's photosynthesis apparatus during the latent period. This increases energy and speeds up virus production, allowing for a reduced latent period (a fitness benefit), but it also increases the DNA size, which slows down new virus production and reduces burst size (a fitness cost). How do these genes affect the net ecological fitness of the virus? Here, this question is explored using a combined systems biology and systems ecology ('systems bioecology') approach. A novel agent-based model simulates individual cyanobacteria cells and virus particles, each with their own genes, transcripts, proteins and other properties. The effect of D1 and HLIP proteins is explicitly considered using a mechanistic photosynthesis component. The model is calibrated to the available database for Prochlorococcus ecotype MED4 and podovirus P-SSP7. Laboratory- and field-scale in silico survival, competition and evolution (gene packaging error) experiments with wild type and genetically engineered viruses are performed to develop vertical survival and fitness profiles, and to determine the optimal gene content. The results suggest that photosynthesis genes are nonessential, increase fitness in a manner correlated with irradiance, and that the wild type has an optimal gene content.

Expression of genes involved in nitrogen assimilation and the C/N balance sensing in Prochlorococcus sp. strain SS120

The expression of five genes involved in nitrogen assimilation in cyanobacteria, namely glnA, glsF, icd, ntcA, and glnB, encoding three key enzymes from that pathway (glutamine synthetase, glutamate synthase, isocitrate dehydrogenase) and two regulatory proteins (NtcA and PII), was studied in this work. Their changes under different conditions were analyzed by quantitative real-time RT-PCR. Nutrient limitation induced clear modifications on the expression of most studied genes: lack of nitrogen provoked an initial increase, followed by a marked decrease; in the cases of phosphorus and iron starvation, a general, stronger expression decrease was observed, particularly striking in the case of iron. Darkness and addition of the photosynthethic inhibitors DCMU and DBMIB also had a strong effect on gene expression. Methionine sulfoximine and azaserine, inhibitors of glutamine synthetase and glutamate synthase, respectively, provoked a sharp increase in icd expression. These results, together with previous studies, suggest that 2-oxoglutarate could be the molecule utilized by Prochlorococcus to sense the C/N balance. Besides, our results confirm the different regulation of nitrogen assimilation in Prochlorococcus with regard to other cyanobacteria.

The challenge of regulation in a minimal photoautotroph: non-coding RNAs in Prochlorococcus

Prochlorococcus, an extremely small cyanobacterium that is very abundant in the world's oceans, has a very streamlined genome. On average, these cells have about 2,000 genes and very few regulatory proteins. The limited capability of regulation is thought to be a result of selection imposed by a relatively stable environment in combination with a very small genome. Furthermore, only ten non-coding RNAs (ncRNAs), which play crucial regulatory roles in all forms of life, have been described in Prochlorococcus. Most strains also lack the RNA chaperone Hfq, raising the question of how important this mode of regulation is for these cells. To explore this question, we examined the transcription of intergenic regions of Prochlorococcus MED4 cells subjected to a number of different stress conditions: changes in light qualities and quantities, phage infection, or phosphorus starvation. Analysis of Affymetrix microarray expression data from intergenic regions revealed 276 novel transcriptional units. Among these were 12 new ncRNAs, 24 antisense RNAs (asRNAs), as well as 113 short mRNAs. Two additional ncRNAs were identified by homology, and all 14 new ncRNAs were independently verified by Northern hybridization and 5'RACE. Unlike its reduced suite of regulatory proteins, the number of ncRNAs relative to genome size in Prochlorococcus is comparable to that found in other bacteria, suggesting that RNA regulators likely play a major role in regulation in this group. Moreover, the ncRNAs are concentrated in previously identified genomic islands, which carry genes of significance to the ecology of this organism, many of which are not of cyanobacterial origin. Expression profiles of some of these ncRNAs suggest involvement in light stress adaptation and/or the response to phage infection consistent with their location in the hypervariable genomic islands.

Patterns and implications of gene gain and loss in the evolution of Prochlorococcus

Prochlorococcus is a marine cyanobacterium that numerically dominates the mid-latitude oceans and is the smallest known oxygenic phototroph. Numerous isolates from diverse areas of the world's oceans have been studied and shown to be physiologically and genetically distinct. All isolates described thus far can be assigned to either a tightly clustered high-light (HL)-adapted clade, or a more divergent low-light (LL)-adapted group. The 16S rRNA sequences of the entire Prochlorococcus group differ by at most 3%, and the four initially published genomes revealed patterns of genetic differentiation that help explain physiological differences among the isolates. Here we describe the genomes of eight newly sequenced isolates and combine them with the first four genomes for a comprehensive analysis of the core (shared by all isolates) and flexible genes of the Prochlorococcus group, and the patterns of loss and gain of the flexible genes over the course of evolution. There are 1,273 genes that represent the core shared by all 12 genomes. They are apparently sufficient, according to metabolic reconstruction, to encode a functional cell. We describe a phylogeny for all 12 isolates by subjecting their complete proteomes to three different phylogenetic analyses. For each non-core gene, we used a maximum parsimony method to estimate which ancestor likely first acquired or lost each gene. Many of the genetic differences among isolates, especially for genes involved in outer membrane synthesis and nutrient transport, are found within the same clade. Nevertheless, we identified some genes defining HL and LL ecotypes, and clades within these broad ecotypes, helping to demonstrate the basis of HL and LL adaptations in Prochlorococcus. Furthermore, our estimates of gene gain events allow us to identify highly variable genomic islands that are not apparent through simple pairwise comparisons. These results emphasize the functional roles, especially those connected to outer membrane synthesis and transport that dominate the flexible genome and set it apart from the core. Besides identifying islands and demonstrating their role throughout the history of Prochlorococcus, reconstruction of past gene gains and losses shows that much of the variability exists at the "leaves of the tree," between the most closely related strains. Finally, the identification of core and flexible genes from this 12-genome comparison is largely consistent with the relative frequency of Prochlorococcus genes found in global ocean metagenomic databases, further closing the gap between our understanding of these organisms in the lab and the wild.

A novel two domain-fusion protein in cyanobacteria with similarity to the CAB/ELIP/HLIP superfamily: evolutionary implications and regulation

Vascular plants contain abundant, light-harvesting complexes in the thylakoid membrane that are non-covalently associated with chlorophylls and carotenoids. These light-harvesting chlorophyll a/b binding (LHC) proteins are members of an extended CAB/ELIP/HLIP superfamily of distantly related polypeptides, which have between one and four transmembrane helices (TMH). This superfamily includes the single TMH, high-light-inducible proteins (Hlips), found in cyanobacteria that are induced by various stress conditions, including high light, and are considered ancestral to the LHC proteins. The roles of, and evolutionary relationships between, these superfamily members are of particular interest, since they function in both light harvesting and photoprotection and may have evolved through tandem gene duplication and fusion events. We have investigated the Hlips (hli gene family) in the thermophilic unicellular cyanobacterium Synechococcus OS-B'. The five hli genes present on the genome of Synechococcus OS-B' are relatively similar, but transcript analyses indicate that there are different patterns of transcript accumulation when the cells are exposed to various growth conditions, suggesting that different Hlips may have specific functions. Hlip5 has an additional TMH at the N-terminus as a result of a novel fusion event. This additional TMH is very similar to a conserved hypothetical, single membrane-spanning polypeptide present in most cyanobacteria. The evolutionary significance of these results is discussed.