A global analysis of adaptive evolution of operons in cyanobacteria

Expanding the synthetic biology repertoire of a fast-growing cyanobacterium Synechococcus elongatus PCC 11801

Synechococcus elongatus PCC 11801 is a fast-growing cyanobacterium, exhibiting high tolerance to environmental stresses. We have earlier characterized its genome and analysed its transcriptome and proteome. However, to deploy it as a potential cell factory, it is necessary to expand its synthetic biology toolbox, including promoter elements and ribosome binding sites (RBSs). Here, based on the global transcriptome analysis, 48 native promoters of the genes with high transcript count were characterized using a fluorescent reporter system. The promoters PcpcB, PpsbA1, and P11770 exhibited consistently high fluorescence under all the cultivation conditions. Similarly, from the genome data and proteome analysis, 534 operons were identified. Fifteen intergenic regions exhibiting higher protein expression from the downstream gene were systematically characterized for identifying RBSs, using an operon construct comprising fluorescent protein genes eyfp and mTurq under PcpcB (PcpcB:eyfp:RBS:mTurq:TrrnB). Overall, the work presents promoter and RBS sequence libraries, with varying strengths, to expedite bioengineering of PCC 11801.

A putative bifunctional CPD/ (6-4) photolyase from the cyanobacteria Synechococcus sp. PCC 7335 is encoded by a UV-B inducible operon: New insights into the evolution of photolyases

Photosynthetic organisms are continuously exposed to solar ultraviolet radiation-B (UV-B) because of their autotrophic lifestyle. UV-B provokes DNA damage, such as cyclobutane pyrimidine dimers (CPD) or pyrimidine (6-4) pyrimidone photoproducts (6-4 PPs). The cryptochrome/photolyase family (CPF) comprises flavoproteins that can bind damaged or undamaged DNA. Photolyases (PHRs) are enzymes that repair either CPDs or 6-4 PPs. A natural bifunctional CPD/(6-4)- PHR (PhrSph98) was recently isolated from the UV-resistant bacteria Sphingomonas sp. UV9. In this work, phylogenetic studies of bifunctional CPD/(6-4)- photolyases and their evolutionary relationship with other CPF members were performed. Amino acids involved in electron transfer and binding to FAD cofactor and DNA lesions were conserved in proteins from proteobacteria, planctomycete, bacteroidete, acidobacteria and cyanobacteria clades. Genome analysis revealed that the cyanobacteria Synechococcus sp. PCC 7335 encodes a two-gene assembly operon coding for a PHR and a bifunctional CPD/(6-4) PHR- like. Operon structure was validated by RT-qPCR analysis and the polycistronic transcript accumulated after 15 min of UV-B irradiation. Conservation of structure and evolution is discussed. This study provides evidence for a UV-B inducible PHR operon that encodes a CPD/(6-4)- photolyase homolog with a putative bifunctional role in the repair of CPDs and 6-4 PPs damages in oxygenic photosynthetic organisms.

The Conserved Family of the Pyridoxal Phosphate-Binding Protein (PLPBP) and Its Cyanobacterial Paradigm PipY

The PLPBP family of pyridoxal phosphate-binding proteins has a high degree of sequence conservation and is represented in all three domains of life. PLPBP members, of which a few representatives have been studied in different contexts, are single-domain proteins with no known enzymatic activity that exhibit the fold type III of PLP-holoenzymes, consisting in an α/β barrel (TIM-barrel), where the PLP cofactor is solvent-exposed. Despite the constant presence of cofactor PLP (a key catalytic element in PLP enzymes), PLPBP family members appear to have purely regulatory functions affecting the homeostasis of vitamin B6 vitamers and amino/keto acids. Perturbation of these metabolites and pleiotropic phenotypes have been reported in bacteria and zebrafish after PLPBP gene inactivation as well as in patients with vitamin B6-dependent epilepsy that results from loss-of-function mutations at the PLPBP. Here, we review information gathered from diverse studies and biological systems, emphasizing the structural and functional conservation of the PLPBP members and discussing the informative nature of model systems and experimental approaches. In this context, the relatively high level of structural and functional characterization of PipY from Synechococcus elongatus PCC 7942 provides a unique opportunity to investigate the PLPBP roles in the context of a signaling pathway conserved in cyanobacteria.

OUP accepted manuscript

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion–deletion–excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.

The Exploration of Novel Regulatory Relationships Drives Haloarchaeal Operon-Like Structural Dynamics over Short Evolutionary Distances

Operons are a dominant feature of bacterial and archaeal genome organization. Numerous investigations have related aspects of operon structure to operon function, making operons exemplars for studies aimed at deciphering Nature's design principles for genomic organization at a local scale. We consider this understanding to be both fundamentally important and ultimately useful in the de novo design of increasingly complex synthetic circuits. Here we analyze the evolution of the genomic context of operon-like structures in a set of 76 sequenced and annotated species of halophilic archaea. The phylogenetic depth and breadth of this dataset allows insight into changes in operon-like structures over shorter evolutionary time scales than have been studied in previous cross-species analysis of operon evolution. Our analysis, implemented in the updated software package JContextExplorer finds that operon-like context as measured by changes in structure frequently differs from a sequence divergence model of whole-species phylogeny and that changes seem to be dominated by the exploration of novel regulatory relationships.

RNA helicase-regulated processing of the Synechocystis rimO-crhR operon results in differential cistron expression and accumulation of two sRNAs

The arrangement of functionally-related genes in operons is a fundamental element of how genetic information is organized in prokaryotes. This organization ensures coordinated gene expression by co-transcription. Often, however, alternative genetic responses to specific stress conditions demand the discoordination of operon expression. During cold temperature stress, accumulation of the gene encoding the sole Asp-Glu-Ala-Asp (DEAD)-box RNA helicase in Synechocystis sp. PCC 6803, crhR (slr0083), increases 15-fold. Here, we show that crhR is expressed from a dicistronic operon with the methylthiotransferase rimO/miaB (slr0082) gene, followed by rapid processing of the operon transcript into two monocistronic mRNAs. This cleavage event is required for and results in destabilization of the rimO transcript. Results from secondary structure modeling and analysis of RNase E cleavage of the rimO-crhR transcript in vitro suggested that CrhR plays a role in enhancing the rate of the processing in an auto-regulatory manner. Moreover, two putative small RNAs are generated from additional processing, degradation, or both of the rimO transcript. These results suggest a role for the bacterial RNA helicase CrhR in RNase E-dependent mRNA processing in Synechocystis and expand the known range of organisms possessing small RNAs derived from processing of mRNA transcripts.

Recent advances in synthetic biology of cyanobacteria

Cyanobacteria are attractive hosts that can be engineered for the photosynthetic production of fuels, fine chemicals, and proteins from CO₂. Moreover, the responsiveness of these photoautotrophs towards different environmental signals, such as light, CO₂, diurnal cycle, and metals make them potential hosts for the development of biosensors. However, engineering these hosts proves to be a challenging and lengthy process. Synthetic biology can make the process of biological engineering more predictable through the use of standardized biological parts that are well characterized and tools to assemble them. While significant progress has been made with model heterotrophic organisms, many of the parts and tools are not portable in cyanobacteria. Therefore, efforts are underway to develop and characterize parts derived from cyanobacteria. In this review, we discuss the reported parts and tools with the objective to develop cyanobacteria as cell factories or biosensors. We also discuss the issues related to characterization, tunability, portability, and the need to develop enabling technologies to engineer this "green" chassis.

Genome Sequence and Composition of a Tolyporphin-Producing Cyanobacterium-Microbial Community

The cyanobacterial culture HT-58-2 was originally described as a strain of Tolypothrix nodosa with the ability to produce tolyporphins, which comprise a family of distinct tetrapyrrole macrocycles with reported efflux pump inhibition properties. Upon reviving the culture from what was thought to be a nonextant collection, studies of culture conditions, strain characterization, phylogeny, and genomics have been undertaken. Here, HT-58-2 was shown by 16S rRNA analysis to closely align with Brasilonema strains and not with Tolypothrix isolates. Light, fluorescence, and scanning electron microscopy revealed cyanobacterium filaments that are decorated with attached bacteria and associated with free bacteria. Metagenomic surveys of HT-58-2 cultures revealed a diversity of bacteria dominated by Erythrobacteraceae, 97% of which are Porphyrobacter species. A dimethyl sulfoxide washing procedure was found to yield enriched cyanobacterial DNA (presumably by removing community bacteria) and sequence data sufficient for genome assembly. The finished, closed HT-58-2Cyano genome consists of 7.85 Mbp (42.6% G+C) and contains 6,581 genes. All genes for biosynthesis of tetrapyrroles (e.g., heme, chlorophyll a, and phycocyanobilin) and almost all for cobalamin were identified dispersed throughout the chromosome. Among the 6,177 protein-encoding genes, coding sequences (CDSs) for all but two of the eight enzymes for conversion of glutamic acid to protoporphyrinogen IX also were found within one major gene cluster. The cluster also includes 10 putative genes (and one hypothetical gene) encoding proteins with domains for a glycosyltransferase, two cytochrome P450 enzymes, and a flavin adenine dinucleotide (FAD)-binding protein. The composition of the gene cluster suggests a possible role in tolyporphin biosynthesis.IMPORTANCE A worldwide search more than 25 years ago for cyanobacterial natural products with anticancer activity identified a culture (HT-58-2) from Micronesia that produces tolyporphins. Tolyporphins are tetrapyrroles, like chlorophylls, but have several profound structural differences that reside outside the bounds of known biosynthetic pathways. To begin probing the biosynthetic origin and biological function of tolyporphins, our research has focused on studying the cyanobacterial strain, about which almost nothing has been previously reported. We find that the HT-58-2 culture is composed of the cyanobacterium and a community of associated bacteria, complicating the question of which organisms make tolyporphins. Elucidation of the cyanobacterial genome revealed an intriguing gene cluster that contains tetrapyrrole biosynthesis genes and a collection of unknown genes, suggesting that the cluster may be responsible for tolyporphin production. Knowledge of the genome and the gene cluster sharply focuses research to identify related cyanobacterial producers of tolyporphins and delineate the tolyporphin biosynthetic pathway.

Minimal tool set for a prokaryotic circadian clock

BACKGROUND

Circadian clocks are found in organisms of almost all domains including photosynthetic Cyanobacteria, whereby large diversity exists within the protein components involved. In the model cyanobacterium Synechococcus elongatus PCC 7942 circadian rhythms are driven by a unique KaiABC protein clock, which is embedded in a network of input and output factors. Homologous proteins to the KaiABC clock have been observed in Bacteria and Archaea, where evidence for circadian behavior in these domains is accumulating. However, interaction and function of non-cyanobacterial Kai-proteins as well as homologous input and output components remain mainly unclear.

RESULTS

Using a universal BLAST analyses, we identified putative KaiC-based timing systems in organisms outside as well as variations within Cyanobacteria. A systematic analyses of publicly available microarray data elucidated interesting variations in circadian gene expression between different cyanobacterial strains, which might be correlated to the diversity of genome encoded clock components. Based on statistical analyses of co-occurrences of the clock components homologous to Synechococcus elongatus PCC 7942, we propose putative networks of reduced and fully functional clock systems. Further, we studied KaiC sequence conservation to determine functionally important regions of diverged KaiC homologs. Biochemical characterization of exemplary cyanobacterial KaiC proteins as well as homologs from two thermophilic Archaea demonstrated that kinase activity is always present. However, a KaiA-mediated phosphorylation is only detectable in KaiC1 orthologs.

CONCLUSION

Our analysis of 11,264 genomes clearly demonstrates that components of the Synechococcus elongatus PCC 7942 circadian clock are present in Bacteria and Archaea. However, all components are less abundant in other organisms than Cyanobacteria and KaiA, Pex, LdpA, and CdpA are only present in the latter. Thus, only reduced KaiBC-based or even simpler, solely KaiC-based timing systems might exist outside of the cyanobacterial phylum, which might be capable of driving diurnal oscillations.

PipY, a Member of the Conserved COG0325 Family of PLP-Binding Proteins, Expands the Cyanobacterial Nitrogen Regulatory Network

Synechococcus elongatus PCC 7942 is a paradigmatic model organism for nitrogen regulation in cyanobacteria. Expression of genes involved in nitrogen assimilation is positively regulated by the 2-oxoglutarate receptor and global transcriptional regulator NtcA. Maximal activation requires the subsequent binding of the co-activator PipX. PII, a protein found in all three domains of life as an integrator of signals of the nitrogen and carbon balance, binds to PipX to counteract NtcA activity at low 2-oxoglutarate levels. PII-PipX complexes can also bind to the transcriptional regulator PlmA, whose regulon remains unknown. Here we expand the nitrogen regulatory network to PipY, encoded by the bicistronic operon pipXY in S. elongatus. Work with PipY, the cyanobacterial member of the widespread family of COG0325 proteins, confirms the conserved roles in vitamin B6 and amino/keto acid homeostasis and reveals new PLP-related phenotypes, including sensitivity to antibiotics targeting essential PLP-holoenzymes or synthetic lethality with cysK. In addition, the related phenotypes of pipY and pipX mutants are consistent with genetic interactions in the contexts of survival to PLP-targeting antibiotics and transcriptional regulation. We also showed that PipY overexpression increased the length of S. elongatus cells. Taken together, our results support a universal regulatory role for COG0325 proteins, paving the way to a better understanding of these proteins and of their connections with other biological processes.

Gene essentiality, conservation index and co-evolution of genes in cyanobacteria

Cyanobacteria, a group of photosynthetic prokaryotes, dominate the earth with ~ 1015 g wet biomass. Despite diversity in habitats and an ancient origin, cyanobacterial phylum has retained a significant core genome. Cyanobacteria are being explored for direct conversion of solar energy and carbon dioxide into biofuels. For this, efficient cyanobacterial strains will need to be designed via metabolic engineering. This will require identification of target knockouts to channelize the flow of carbon toward the product of interest while minimizing deletions of essential genes. We propose "Gene Conservation Index" (GCI) as a quick measure to predict gene essentiality in cyanobacteria. GCI is based on phylogenetic profile of a gene constructed with a reduced dataset of cyanobacterial genomes. GCI is the percentage of organism clusters in which the query gene is present in the reduced dataset. Of the 750 genes deemed to be essential in the experimental study on S. elongatus PCC 7942, we found 494 to be conserved across the phylum which largely comprise of the essential metabolic pathways. On the contrary, the conserved but non-essential genes broadly comprise of genes required under stress conditions. Exceptions to this rule include genes such as the glycogen synthesis and degradation enzymes, deoxyribose-phosphate aldolase (DERA), glucose-6-phosphate 1-dehydrogenase (zwf) and fructose-1,6-bisphosphatase class1, which are conserved but non-essential. While the essential genes are to be avoided during gene knockout studies as potentially lethal deletions, the non-essential but conserved set of genes could be interesting targets for metabolic engineering. Further, we identify clusters of co-evolving genes (CCG), which provide insights that may be useful in annotation. Principal component analysis (PCA) plots of the CCGs are demonstrated as data visualization tools that are complementary to the conventional heatmaps. Our dataset consists of phylogenetic profiles for 23,643 non-redundant cyanobacterial genes. We believe that the data and the analysis presented here will be a great resource to the scientific community interested in cyanobacteria.

A Comprehensive Study of Cyanobacterial Morphological and Ecological Evolutionary Dynamics through Deep Geologic Time

Cyanobacteria have exerted a profound influence on the progressive oxygenation of Earth. As a complementary approach to examining the geologic record—phylogenomic and trait evolutionary analyses of extant species can lead to new insights. We constructed new phylogenomic trees and analyzed phenotypic trait data using novel phylogenetic comparative methods. We elucidated the dynamics of trait evolution in Cyanobacteria over billion-year timescales, and provide evidence that major geologic events in early Earth’s history have shaped—and been shaped by—evolution in Cyanobacteria. We identify a robust core cyanobacterial phylogeny and a smaller set of taxa that exhibit long-branch attraction artifacts. We estimated the age of nodes and reconstruct the ancestral character states of 43 phenotypic characters. We find high levels of phylogenetic signal for nearly all traits, indicating the phylogeny carries substantial predictive power. The earliest cyanobacterial lineages likely lived in freshwater habitats, had small cell diameters, were benthic or sessile, and possibly epilithic/endolithic with a sheath. We jointly analyzed a subset of 25 binary traits to determine whether rates of trait evolution have shifted over time in conjunction with major geologic events. Phylogenetic comparative analysis reveal an overriding signal of decreasing rates of trait evolution through time. Furthermore, the data suggest two major rate shifts in trait evolution associated with bursts of evolutionary innovation. The first rate shift occurs in the aftermath of the Great Oxidation Event and “Snowball Earth” glaciations and is associated with decrease in the evolutionary rates around 1.8–1.6 Ga. This rate shift seems to indicate the end of a major diversification of cyanobacterial phenotypes–particularly related to traits associated with filamentous morphology, heterocysts and motility in freshwater ecosystems. Another burst appears around the time of the Neoproterozoic Oxidation Event in the Neoproterozoic, and is associated with the acquisition of traits involved in planktonic growth in marine habitats. Our results demonstrate how uniting genomic and phenotypic datasets in extant bacterial species can shed light on billion-year old events in Earth’s history.

Minimal Tool Set for a Prokaryotic Circadian Clock.. [DOI: 10.1101/075291]

Background

Circadian clocks are found in organisms of almost all domains including photosynthetic Cyanobacteria, whereby large diversity exists within the protein components involved. In the model cyanobacteriumSynechococcus elongatusPCC 7942 circadian rhythms are driven by a unique KaiABC protein clock, which is embedded in a network of input and output factors. Homologous proteins to the KaiABC clock have been observed in Bacteria and Archaea, where evidence for circadian behavior in these domains is accumulating. However, interaction and function of non-cyanobacterial Kai-proteins as well as homologous input and output components remain mainly unclear.

Results

Using a universal BLAST analyses, we identified putative KaiC-based timing systems in organisms outside as well as variations within Cyanobacteria. A systematic analyses of publicly available microarray data elucidated interesting variations in circadian gene expression between different cyanobacterial strains, which might be correlated to the diversity of genome encoded clock components. Based on statistical analyses of co-occurrences of the clock components homologous toSynechococcus elongatusPCC 7942, we propose putative networks of reduced and fully functional clock systems. Further, we studied KaiC sequence conservation to determine functionally important regions of diverged KaiC homologs. Biochemical characterization of exemplary cyanobacterial KaiC proteins as well as homologs from two thermophilic Archaea demonstrated that kinase activity is always present. However, a KaiA-mediated phosphorylation is only detectable in KaiC1 orthologs.

Conclusion

Our analysis of 11,264 genomes clearly demonstrates that components of theSynechococcus elongatusPCC 7942 circadian clock are present in Bacteria and Archaea. However, all components are less abundant in other organisms than Cyanobacteria and KaiA, Pex, LdpA, and CdpA are only present in the latter. Thus, only reduced KaiBC-based or even simpler, solely KaiC-based timing systems might exist outside of the cyanobacterial phylum, which might be capable of driving diurnal oscillations.

Investigating Evolutionary Dynamics of RHA1 Operons

Grouping genes as operons is an important genomic feature of prokaryotic organisms. The comprehensive understanding of the operon organizations would be helpful to decipher transcriptional mechanisms, cellular pathways, and the evolutionary landscape of prokaryotic genomes. Although thousands of prokaryotes have been sequenced, genome-wide investigation of the evolutionary dynamics (division and recombination) of operons among these genomes remains unexplored. Here, we systematically analyzed the operon dynamics of Rhodococcus jostii RHA1 (RHA1), an oleaginous bacterium with high potential applications in biofuel, by comparing 340 prokaryotic genomes that were carefully selected from different genera. Interestingly, 99% of RHA1 operons were observed to exhibit evolutionary events of division and recombination among the 340 compared genomes. An operon that encodes all enzymes related to histidine biosynthesis in RHA1 (His-operon) was found to be segmented into smaller gene groups (sub-operons) in diverse genomes. These sub-operons were further reorganized with different functional genes as novel operons that are related to different biochemical processes. Comparatively, the operons involved in the functional categories of lipid transport and metabolism are relatively conserved among the 340 compared genomes. At the pathway level, RHA1 operons found to be significantly conserved were involved in ribosome synthesis, oxidative phosphorylation, and fatty acid synthesis. These analyses provide evolutionary insights of operon organization and the dynamic associations of various biochemical pathways in different prokaryotes.

The Chthonomonas calidirosea Genome Is Highly Conserved across Geographic Locations and Distinct Chemical and Microbial Environments in New Zealand's Taupō Volcanic Zone

Chthonomonas calidirosea T49^T is a low-abundance, carbohydrate-scavenging, and thermophilic soil bacterium with a seemingly disorganized genome. We hypothesized that the C. calidirosea genome would be highly responsive to local selection pressure, resulting in the divergence of its genomic content, genome organization, and carbohydrate utilization phenotype across environments. We tested this hypothesis by sequencing the genomes of four C. calidirosea isolates obtained from four separate geothermal fields in the Taupō Volcanic Zone, New Zealand. For each isolation site, we measured physicochemical attributes and defined the associated microbial community by 16S rRNA gene sequencing. Despite their ecological and geographical isolation, the genome sequences showed low divergence (maximum, 1.17%). Isolate-specific variations included single-nucleotide polymorphisms (SNPs), restriction-modification systems, and mobile elements but few major deletions and no major rearrangements. The 50-fold variation in C. calidirosea relative abundance among the four sites correlated with site environmental characteristics but not with differences in genomic content. Conversely, the carbohydrate utilization profiles of the C. calidirosea isolates corresponded to the inferred isolate phylogenies, which only partially paralleled the geographical relationships among the sample sites. Genomic sequence conservation does not entirely parallel geographic distance, suggesting that stochastic dispersal and localized extinction, which allow for rapid population homogenization with little restriction by geographical barriers, are possible mechanisms of C. calidirosea distribution. This dispersal and extinction mechanism is likely not limited to C. calidirosea but may shape the populations and genomes of many other low-abundance free-living taxa.

IMPORTANCE This study compares the genomic sequence variations and metabolisms of four strains of Chthonomonas calidirosea, a rare thermophilic bacterium from the phylum Armatimonadetes. It additionally compares the microbial communities and chemistry of each of the geographically distinct sites from which the four C. calidirosea strains were isolated. C. calidirosea was previously reported to possess a highly disorganized genome, but it was unclear whether this reflected rapid evolution. Here, we show that each isolation site has a distinct chemistry and microbial community, but despite this, the C. calidirosea genome is highly conserved across all isolation sites. Furthermore, genomic sequence differences only partially paralleled geographic distance, suggesting that C. calidirosea genotypes are not primarily determined by adaptive evolution. Instead, the presence of C. calidirosea may be driven by stochastic dispersal and localized extinction. This ecological mechanism may apply to many other low-abundance taxa.

Coupling of Cellular Processes and Their Coordinated Oscillations under Continuous Light in Cyanothece sp. ATCC 51142, a Diazotrophic Unicellular Cyanobacterium

Unicellular diazotrophic cyanobacteria such as Cyanothece sp. ATCC 51142 (henceforth Cyanothece), temporally separate the oxygen sensitive nitrogen fixation from oxygen evolving photosynthesis not only under diurnal cycles (LD) but also in continuous light (LL). However, recent reports demonstrate that the oscillations in LL occur with a shorter cycle time of ~11 h. We find that indeed, majority of the genes oscillate in LL with this cycle time. Genes that are upregulated at a particular time of day under diurnal cycle also get upregulated at an equivalent metabolic phase under LL suggesting tight coupling of various cellular events with each other and with the cell's metabolic status. A number of metabolic processes get upregulated in a coordinated fashion during the respiratory phase under LL including glycogen degradation, glycolysis, oxidative pentose phosphate pathway, and tricarboxylic acid cycle. These precede nitrogen fixation apparently to ensure sufficient energy and anoxic environment needed for the nitrogenase enzyme. Photosynthetic phase sees upregulation of photosystem II, carbonate transport, carbon concentrating mechanism, RuBisCO, glycogen synthesis and light harvesting antenna pigment biosynthesis. In Synechococcus elongates PCC 7942, a non-nitrogen fixing cyanobacteria, expression of a relatively smaller fraction of genes oscillates under LL condition with the major periodicity being 24 h. In contrast, the entire cellular machinery of Cyanothece orchestrates coordinated oscillation in anticipation of the ensuing metabolic phase in both LD and LL. These results may have important implications in understanding the timing of various cellular events and in engineering cyanobacteria for biofuel production.

The structural code of cyanobacterial genomes

A periodic bias in nucleotide frequency with a period of about 11 bp is characteristic for bacterial genomes. This signal is commonly interpreted to relate to the helical pitch of negatively supercoiled DNA. Functions in supercoiling-dependent RNA transcription or as a 'structural code' for DNA packaging have been suggested. Cyanobacterial genomes showed especially strong periodic signals and, on the other hand, DNA supercoiling and supercoiling-dependent transcription are highly dynamic and underlie circadian rhythms of these phototrophic bacteria. Focusing on this phylum and dinucleotides, we find that a minimal motif of AT-tracts (AT2) yields the strongest signal. Strong genome-wide periodicity is ancestral to a clade of unicellular and polyploid species but lost upon morphological transitions into two baeocyte-forming and a symbiotic species. The signal is intermediate in heterocystous species and weak in monoploid picocyanobacteria. A pronounced 'structural code' may support efficient nucleoid condensation and segregation in polyploid cells. The major source of the AT2 signal are protein-coding regions, where it is encoded preferentially in the first and third codon positions. The signal shows only few relations to supercoiling-dependent and diurnal RNA transcription in Synechocystis sp. PCC 6803. Strong and specific signals in two distinct transposons suggest roles in transposase transcription and transpososome formation.

The transcriptome landscape of Prochlorococcus MED4 and the factors for stabilizing the core genome

BACKGROUND

Gene gain and loss frequently occurs in the cyanobacterium Prochlorococcus, a phototroph that numerically dominates tropical and subtropical open oceans. However, little is known about the stabilization of its core genome, which contains approximately 1250 genes, in the context of genome streamlining. Using Prochlorococcus MED4 as a model organism, we investigated the constraints on core genome stabilization using transcriptome profiling.

RESULTS

RNA-Seq technique was used to obtain the transcriptome map of Prochlorococcus MED4, including operons, untranslated regions, non-coding RNAs, and novel genes. Genome-wide expression profiles revealed that three factors contribute to core genome stabilization. First, a negative correlation between gene expression levels and protein evolutionary rates was observed. Highly expressed genes were overrepresented in the core genome but not in the flexible genome. Gene necessity was determined as a second powerful constraint on genome evolution through functional enrichment analysis. Third, quick mRNA turnover may increase corresponding proteins' fidelity among genes that were abundantly expressed. Together, these factors influence core genome stabilization during MED4 genome evolution.

CONCLUSIONS

Gene expression, gene necessity, and mRNA turnover contribute to core genome maintenance during cyanobacterium Prochlorococcus genus evolution.

Rhythmic and sustained oscillations in metabolism and gene expression of Cyanothece sp. ATCC 51142 under constant light

Cyanobacteria, a group of photosynthetic prokaryotes, oscillate between day and night time metabolisms with concomitant oscillations in gene expression in response to light/dark cycles (LD). The oscillations in gene expression have been shown to sustain in constant light (LL) with a free running period of 24 h in a model cyanobacterium Synechococcus elongatus PCC 7942. However, equivalent oscillations in metabolism are not reported under LL in this non-nitrogen fixing cyanobacterium. Here we focus on Cyanothece sp. ATCC 51142, a unicellular, nitrogen-fixing cyanobacterium known to temporally separate the processes of oxygenic photosynthesis and oxygen-sensitive nitrogen fixation. In a recent report, metabolism of Cyanothece 51142 has been shown to oscillate between photosynthetic and respiratory phases under LL with free running periods that are temperature dependent but significantly shorter than the circadian period. Further, the oscillations shift to circadian pattern at moderate cell densities that are concomitant with slower growth rates. Here we take this understanding forward and demonstrate that the ultradian rhythm under LL sustains at much higher cell densities when grown under turbulent regimes that simulate flashing light effect. Our results suggest that the ultradian rhythm in metabolism may be needed to support higher carbon and nitrogen requirements of rapidly growing cells under LL. With a comprehensive Real time PCR based gene expression analysis we account for key regulatory interactions and demonstrate the interplay between clock genes and the genes of key metabolic pathways. Further, we observe that several genes that peak at dusk in Synechococcus peak at dawn in Cyanothece and vice versa. The circadian rhythm of this organism appears to be more robust with peaking of genes in anticipation of the ensuing photosynthetic and respiratory metabolic phases.

SHARP: genome-scale identification of gene-protein-reaction associations in cyanobacteria

Genome scale metabolic model provides an overview of an organism's metabolic capability. These genome-specific metabolic reconstructions are based on identification of gene to protein to reaction (GPR) associations and, in turn, on homology with annotated genes from other organisms. Cyanobacteria are photosynthetic prokaryotes which have diverged appreciably from their nonphotosynthetic counterparts. They also show significant evolutionary divergence from plants, which are well studied for their photosynthetic apparatus. We argue that context-specific sequence and domain similarity can add to the repertoire of the GPR associations and significantly expand our view of the metabolic capability of cyanobacteria. We took an approach that combines the results of context-specific sequence-to-sequence similarity search with those of sequence-to-profile searches. We employ PSI-BLAST for the former, and CDD, Pfam, and COG for the latter. An optimization algorithm was devised to arrive at a weighting scheme to combine the different evidences with KEGG-annotated GPRs as training data. We present the algorithm in the form of software "Systematic, Homology-based Automated Re-annotation for Prokaryotes (SHARP)." We predicted 3,781 new GPR associations for the 10 prokaryotes considered of which eight are cyanobacteria species. These new GPR associations fall in several metabolic pathways and were used to annotate 7,718 gaps in the metabolic network. These new annotations led to discovery of several pathways that may be active and thereby providing new directions for metabolic engineering of these species for production of useful products. Metabolic model developed on such a reconstructed network is likely to give better phenotypic predictions.

Highlights from the Indo-US workshop "Cyanobacteria: molecular networks to biofuels" held at Lonavala, India during December 16-20, 2012

An Indo-US workshop on "Cyanobacteria: molecular networks to biofuels" was held December 16-20, 2012 at Lagoona Resort, Lonavala, India. The workshop was jointly organized by two of the authors, PPW, a chemical engineer and LAS, a biologist, thereby ensuring a broad and cross-disciplinary participation. The main objective of the workshop was to bring researchers from academia and industry of the two countries together with common interests in cyanobacteria or microalgae and derived biofuels. An exchange of ideas resulted from a series of oral and poster presentations and, importantly, through one-on-one discussions during tea breaks and meals. Another key objective was to introduce young researchers of India to the exciting field of cyanobacterial physiology, modeling, and biofuels. PhD students and early stage researchers were especially encouraged to participate and about half of the 75 participants belonged to this category. The rest were comprised of senior researchers, including 13-15 invited speakers from each country. Overall, twenty-four institutes from 12 states of India were represented. The deliberations, which are being compiled in the present special issue, revolved mainly around molecular aspects of cyanobacterial biofuels including metabolic engineering, networks, genetic regulation, circadian rhythms, and stress responses. Representatives of some key funding agencies and industry provided a perspective and opportunities in the field and for bilateral collaboration. This article summarizes deliberations that took place at the meeting and provides a bird's eye view of the ongoing research in the field in the two countries.

Metabolic flux analysis of Cyanothece sp. ATCC 51142 under mixotrophic conditions

Cyanobacteria are a group of photosynthetic prokaryotes capable of utilizing solar energy to fix atmospheric carbon dioxide to biomass. Despite several "proof of principle" studies, low product yield is an impediment in commercialization of cyanobacteria-derived biofuels. Estimation of intracellular reaction rates by (13)C metabolic flux analysis ((13)C-MFA) would be a step toward enhancing biofuel yield via metabolic engineering. We report (13)C-MFA for Cyanothece sp. ATCC 51142, a unicellular nitrogen-fixing cyanobacterium, known for enhanced hydrogen yield under mixotrophic conditions. Rates of reactions in the central carbon metabolism under nitrogen-fixing and -non-fixing conditions were estimated by monitoring the competitive incorporation of (12)C and (13)C from unlabeled CO2 and uniformly labeled glycerol, respectively, into terminal metabolites such as amino acids. The observed labeling patterns suggest mixotrophic growth under both the conditions, with a larger fraction of unlabeled carbon in nitrate-sufficient cultures asserting a greater contribution of carbon fixation by photosynthesis and an anaplerotic pathway. Indeed, flux analysis complements the higher growth observed under nitrate-sufficient conditions. On the other hand, the flux through the oxidative pentose phosphate pathway and tricarboxylic acid cycle was greater in nitrate-deficient conditions, possibly to supply the precursors and reducing equivalents needed for nitrogen fixation. In addition, an enhanced flux through fructose-6-phosphate phosphoketolase possibly suggests the organism's preferred mode under nitrogen-fixing conditions. The (13)C-MFA results complement the reported predictions by flux balance analysis and provide quantitative insight into the organism's distinct metabolic features under nitrogen-fixing and -non-fixing conditions.

Metabolic modeling for multi-objective optimization of ethanol production in a Synechocystis mutant

Cyanobacteria have potential to produce drop-in bio-fuels such as ethanol via photoautotrophic metabolism. Although model cyanobacterial strains have been engineered to produce such products, systematic metabolic engineering studies to identify optimal strains for the same have not been performed. In this work, we identify optimal ethanol producing mutants corresponding to appropriate gene deletions that result in a suitable redirection in the carbon flux. In particular, we systematically simulate exhaustive single and double gene deletions considering a genome scale metabolic model of a mutant strain of the unicellular cyanobacterium Synechocystis species strain PCC 6803. Various optimization based metabolic modeling techniques, such as flux balance analysis (FBA), method of minimization of metabolic adjustment (MOMA) and regulatory on/off minimization (ROOM) were used for this analysis. For single gene deletion MOMA simulations, the Pareto front with biomass and ethanol fluxes as the two objectives to be maximized was obtained and analyzed. Points on the Pareto front represent maximal utilization of resources constrained by substrate uptake thereby representing an optimal trade-off between the two fluxes. Pareto analysis was also performed for double gene deletion MOMA and single and double gene deletion ROOM simulations. Based on these analyses, two mutants, with combined gene deletions in ethanol and purine metabolism pathways, were identified as promising candidates for ethanol production. The relevant genes were adk, pta and ackA. An ethanol productivity of approximately 0.15 mmol/(gDW h) was predicted for these mutants which appears to be reasonable based on experimentally reported values in literature for other strains.

Rhythm of carbon and nitrogen fixation in unicellular cyanobacteria under turbulent and highly aerobic conditions

Nitrogen fixing cyanobacteria are being increasingly explored for nitrogenase-dependent hydrogen production. Commercial success however will depend on the ability to grow these cultures at high cell densities. Photo-limitation at high cell densities leads to hindered photoautotrophic growth while turbulent conditions, which simulate flashing light effect, can lead to oxygen toxicity to the nitrogenase enzyme. Cyanothece sp. strain ATCC 51142, a known hydrogen producer, is reported to grow and fix nitrogen under moderately oxic conditions in shake flasks. In this study, we explore the growth and nitrogen fixing potential of this organism under turbulent conditions with volumetric oxygen mass transfer coefficient (KL a) values that are up to 20-times greater than in shake flasks. In a stirred vessel, the organism grows well in turbulent regime possibly due to a simulated flashing light effect with optimal growth at Reynolds number of approximately 35,000. A respiratory burst lasting for about 4 h creates anoxic conditions intracellularly with near saturating levels of dissolved oxygen in the extracellular medium. This is concomitant with complete exhaustion of intracellular glycogen storage and upregulation of nifH and nifX, the genes encoding proteins of the nitrogenase complex. Further, the rhythmic oscillations in exhaust gas CO2 and O2 profiles synchronize faithfully with those in biochemical parameters and gene expression thereby serving as an effective online monitoring tool. These results will have important implications in potential commercial success of nitrogenase-dependent hydrogen production by cyanobacteria.