Emergence of fractal geometries in the evolution of a metabolic enzyme.

Fractals are patterns that are self-similar across multiple length-scales1. Macroscopic fractals are common in nature2-4; however, so far, molecular assembly into fractals is restricted to synthetic systems5-12. Here we report the discovery of a natural protein, citrate synthase from the cyanobacterium Synechococcus elongatus, which self-assembles into Sierpinski triangles. Using cryo-electron microscopy, we reveal how the fractal assembles from a hexameric building block. Although different stimuli modulate the formation of fractal complexes and these complexes can regulate the enzymatic activity of citrate synthase in vitro, the fractal may not serve a physiological function in vivo. We use ancestral sequence reconstruction to retrace how the citrate synthase fractal evolved from non-fractal precursors, and the results suggest it may have emerged as a harmless evolutionary accident. Our findings expand the space of possible protein complexes and demonstrate that intricate and regulatable assemblies can evolve in a single substitution.

Circadian control of global gene expression by the cyanobacterial master regulator RpaA.

The cyanobacterial circadian clock generates genome-wide transcriptional oscillations and regulates cell division, but the underlying mechanisms are not well understood. Here, we show that the response regulator RpaA serves as the master regulator of these clock outputs. Deletion of rpaA abrogates gene expression rhythms globally and arrests cells in a dawn-like expression state. Although rpaA deletion causes core oscillator failure by perturbing clock gene expression, rescuing oscillator function does not restore global expression rhythms. We show that phosphorylated RpaA regulates the expression of not only clock components, generating feedback on the core oscillator, but also a small set of circadian effectors that, in turn, orchestrate genome-wide transcriptional rhythms. Expression of constitutively active RpaA is sufficient to switch cells from a dawn-like to a dusk-like expression state as well as to block cell division. Hence, complex global circadian phenotypes can be generated by controlling the phosphorylation of a single transcription factor.

Biogenesis of a bacterial organelle: the carboxysome assembly pathway.

The carboxysome is a protein-based organelle for carbon fixation in cyanobacteria, keystone organisms in the global carbon cycle. It is composed of thousands of subunits including hexameric and pentameric proteins that form a shell to encapsulate the enzymes ribulose 1,5-bisphosphate carboxylase/oxygenase and carbonic anhydrase. Here, we describe the stages of carboxysome assembly and the requisite gene products necessary for progression through each. Our results demonstrate that, unlike membrane-bound organelles of eukaryotes, in carboxysomes the interior of the compartment forms first, at a distinct site within the cell. Subsequently, shell proteins encapsulate this procarboxysome, inducing budding and distribution of functional organelles within the cell. We propose that the principles of carboxysome assembly that we have uncovered extend to diverse bacterial microcompartments.

Machine learning reveals the transcriptional regulatory network and circadian dynamics of Synechococcus elongatus PCC 7942.

Synechococcus elongatus is an important cyanobacterium that serves as a versatile and robust model for studying circadian biology and photosynthetic metabolism. Its transcriptional regulatory network (TRN) is of fundamental interest, as it orchestrates the cell's adaptation to the environment, including its response to sunlight. Despite the previous characterization of constituent parts of the S. elongatus TRN, a comprehensive layout of its topology remains to be established. Here, we decomposed a compendium of 300 high-quality RNA sequencing datasets of the model strain PCC 7942 using independent component analysis. We obtained 57 independently modulated gene sets, or iModulons, that explain 67% of the variance in the transcriptional response and 1) accurately reflect the activity of known transcriptional regulations, 2) capture functional components of photosynthesis, 3) provide hypotheses for regulon structures and functional annotations of poorly characterized genes, and 4) describe the transcriptional shifts under dynamic light conditions. This transcriptome-wide analysis of S. elongatus provides a quantitative reconstruction of the TRN and presents a knowledge base that can guide future investigations. Our systems-level analysis also provides a global TRN structure for S. elongatus PCC 7942.

A distinct, high-affinity, alkaline phosphatase facilitates occupation of P-depleted environments by marine picocyanobacteria.

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus, the two most abundant phototrophs on Earth, thrive in oligotrophic oceanic regions. While it is well known that specific lineages are exquisitely adapted to prevailing in situ light and temperature regimes, much less is known of the molecular machinery required to facilitate occupancy of these low-nutrient environments. Here, we describe a hitherto unknown alkaline phosphatase, Psip1, that has a substantially higher affinity for phosphomonoesters than other well-known phosphatases like PhoA, PhoX, or PhoD and is restricted to clade III Synechococcus and a subset of high light I-adapted Prochlorococcus strains, suggesting niche specificity. We demonstrate that Psip1 has undergone convergent evolution with PhoX, requiring both iron and calcium for activity and likely possessing identical key residues around the active site, despite generally very low sequence homology. Interrogation of metagenomes and transcriptomes from TARA oceans and an Atlantic Meridional transect shows that psip1 is abundant and highly expressed in picocyanobacterial populations from the Mediterranean Sea and north Atlantic gyre, regions well recognized to be phosphorus (P)-deplete. Together, this identifies psip1 as an important oligotrophy-specific gene for P recycling in these organisms. Furthermore, psip1 is not restricted to picocyanobacteria and is abundant and highly transcribed in some α-proteobacteria and eukaryotic algae, suggesting that such a high-affinity phosphatase is important across the microbial taxonomic world to occupy low-P environments.

Dual thermal ecotypes coexist within a nearly genetically identical population of the unicellular marine cyanobacterium Synechococcus.

The extent and ecological significance of intraspecific functional diversity within marine microbial populations is still poorly understood, and it remains unclear if such strain-level microdiversity will affect fitness and persistence in a rapidly changing ocean environment. In this study, we cultured 11 sympatric strains of the ubiquitous marine picocyanobacterium Synechococcus isolated from a Narragansett Bay (RI) phytoplankton community thermal selection experiment. Thermal performance curves revealed selection at cool and warm temperatures had subdivided the initial population into thermotypes with pronounced differences in maximum growth temperatures. Curiously, the genomes of all 11 isolates were almost identical (average nucleotide identities of >99.99%, with >99% of the genome aligning) and no differences in gene content or single nucleotide variants were associated with either cool or warm temperature phenotypes. Despite a very high level of genomic similarity, sequenced epigenomes for two strains showed differences in methylation on genes associated with photosynthesis. These corresponded to measured differences in photophysiology, suggesting a potential pathway for future mechanistic research into thermal microdiversity. Our study demonstrates that present-day marine microbial populations can harbor cryptic but environmentally relevant thermotypes which may increase their resilience to future rising temperatures.

Coexisting picoplankton experience different relative grazing pressures across an ocean productivity gradient.

Picophytoplankton populations [Prochlorococcus, Synechococcus (SYN), and picoeukaryotes] are dominant primary producers in the open ocean and projected to become more important with climate change. Their fates can vary, however, with microbial food web complexities. In the California Current Ecosystem, picophytoplankton biomass and abundance peak in waters of intermediate productivity and decrease at higher production. Using experimental data from eight cruises crossing the pronounced CCE trophic gradient, we tested the hypothesis that these declines are driven by intensified grazing on heterotrophic bacteria (HBAC) passed to similarly sized picophytoplankton via shared predators. Results confirm previously observed distributions as well as significant increases in bacterial abundance, cell growth, and grazing mortality with primary production. Mortalities of picophytoplankton, however, diverge from the bacterial mortality trend such that relative grazing rates on SYN compared to HBAC decline by 12-fold between low and high productivity waters. The large shifts in mortality rate ratios for coexisting populations are not explained by size variability but rather suggest high selectivity of grazer assemblages or tightly coupled tradeoffs in microbial growth advantages and grazing vulnerabilities. These findings challenge the long-held view that protistan grazing mainly determines overall biomass of microbial communities while viruses uniquely regulate diversity by "killing the winners".

Synchronization of the circadian clock to the environment tracked in real time.

The circadian system of the cyanobacterium Synechococcus elongatus PCC 7942 relies on a three-protein nanomachine (KaiA, KaiB, and KaiC) that undergoes an oscillatory phosphorylation cycle with a period of ~24 h. This core oscillator can be reconstituted in vitro and is used to study the molecular mechanisms of circadian timekeeping and entrainment. Previous studies showed that two key metabolic changes that occur in cells during the transition into darkness, changes in the ATP/ADP ratio and redox status of the quinone pool, are cues that entrain the circadian clock. By changing the ATP/ADP ratio or adding oxidized quinone, one can shift the phase of the phosphorylation cycle of the core oscillator in vitro. However, the in vitro oscillator cannot explain gene expression patterns because the simple mixture lacks the output components that connect the clock to genes. Recently, a high-throughput in vitro system termed the in vitro clock (IVC) that contains both the core oscillator and the output components was developed. Here, we used IVC reactions and performed massively parallel experiments to study entrainment, the synchronization of the clock with the environment, in the presence of output components. Our results indicate that the IVC better explains the in vivo clock-resetting phenotypes of wild-type and mutant strains and that the output components are deeply engaged with the core oscillator, affecting the way input signals entrain the core pacemaker. These findings blur the line between input and output pathways and support our previous demonstration that key output components are fundamental parts of the clock.

Structure of the antenna complex expressed during far-red light photoacclimation in Synechococcus sp. PCC 7335.

Far-red light photoacclimation, or FaRLiP, is a facultative response exhibited by some cyanobacteria that allows them to absorb and utilize lower energy light (700-800 nm) than the wavelengths typically used for oxygenic photosynthesis (400-700 nm). During this process, three essential components of the photosynthetic apparatus are altered: photosystem I, photosystem II, and the phycobilisome. In all three cases, at least some of the chromophores found in these pigment-protein complexes are replaced by chromophores that have red-shifted absorbance relative to the analogous complexes produced in visible light. Recent structural and spectroscopic studies have elucidated important features of the two photosystems when altered to absorb and utilize far-red light, but much less is understood about the modified phycobiliproteins made during FaRLiP. We used single-particle, cryo-EM to determine the molecular structure of a phycobiliprotein core complex comprising allophycocyanin variants that absorb far-red light during FaRLiP in the marine cyanobacterium Synechococcus sp. PCC 7335. The structure reveals the arrangement of the numerous red-shifted allophycocyanin variants and the probable locations of the chromophores that serve as the terminal emitters in this complex. It also suggests how energy is transferred to the photosystem II complexes produced during FaRLiP. The structure additionally allows comparisons with other previously studied allophycocyanins to gain insights into how phycocyanobilin chromophores can be tuned to absorb far-red light. These studies provide new insights into how far-red light is harvested and utilized during FaRLiP, a widespread cyanobacterial photoacclimation mechanism.

A toolbox to engineer the highly productive cyanobacterium Synechococcus sp. PCC 11901.

Synechococcus sp. PCC 11901 (PCC 11901) is a fast-growing marine cyanobacterial strain that has a capacity for sustained biomass accumulation to very high cell densities, comparable to that achieved by commercially relevant heterotrophic organisms. However, genetic tools to engineer PCC 11901 for biotechnology applications are limited. Here we describe a suite of tools based on the CyanoGate MoClo system to unlock the engineering potential of PCC 11901. First, we characterized neutral sites suitable for stable genomic integration that do not affect growth even at high cell densities. Second, we tested a suite of constitutive promoters, terminators, and inducible promoters including a 2,4-diacetylphloroglucinol (DAPG)-inducible PhlF repressor system, which has not previously been demonstrated in cyanobacteria and showed tight regulation and a 228-fold dynamic range of induction. Lastly, we developed a DAPG-inducible dCas9-based CRISPR interference (CRISPRi) system and a modular method to generate markerless mutants using CRISPR-Cas12a. Based on our findings, PCC 11901 is highly responsive to CRISPRi-based repression and showed high efficiencies for single insertion (31% to 81%) and multiplex double insertion (25%) genome editing with Cas12a. We envision that these tools will lay the foundations for the adoption of PCC 11901 as a robust model strain for engineering biology and green biotechnology.

Deciphering the structure, function, and mechanism of lysine acetyltransferase cGNAT2 in cyanobacteria.

Lysine acetylation is a conserved regulatory posttranslational protein modification that is performed by lysine acetyltransferases (KATs). By catalyzing the transfer of acetyl groups to substrate proteins, KATs play critical regulatory roles in all domains of life; however, no KATs have yet been identified in cyanobacteria. Here, we tested all predicted KATs in the cyanobacterium Synechococcus sp. PCC 7002 (Syn7002) and demonstrated that A1596, which we named cyanobacterial Gcn5-related N-acetyltransferase (cGNAT2), can catalyze lysine acetylation in vivo and in vitro. Eight amino acid residues were identified as the key residues in the putative active site of cGNAT2, as indicated by structural simulation and site-directed mutagenesis. The loss of cGNAT2 altered both growth and photosynthetic electron transport in Syn7002. In addition, quantitative analysis of the lysine acetylome identified 548 endogenous substrates of cGNAT2 in Syn7002. We further demonstrated that cGNAT2 can acetylate NAD(P)H dehydrogenase J (NdhJ) in vivo and in vitro, with the inability to acetylate K89 residues, thus decreasing NdhJ activity and affecting both growth and electron transport in Syn7002. In summary, this study identified a KAT in cyanobacteria and revealed that cGNAT2 regulates growth and photosynthesis in Syn7002 through an acetylation-mediated mechanism.

Elevated ATPase activity of KaiC applies a circadian checkpoint on cell division in Synechococcus elongatus.

A circadian clock coordinates physiology and behavior in diverse groups of living organisms. Another major cyclic cellular event, the cell cycle, is regulated by the circadian clock in the few cases where linkage of these cycles has been studied. In the cyanobacterium Synechococcus elongatus, the circadian clock gates cell division by an unknown mechanism. Using timelapse microscopy, we confirm the gating of cell division in the wild-type and demonstrate the regulation of cytokinesis by key clock components. Specifically, a state of the oscillator protein KaiC that is associated with elevated ATPase activity closes the gate by acting through a known clock output pathway to inhibit FtsZ ring formation at the division site. An activity that stimulates KaiC phosphorylation independently of the KaiA protein was also uncovered. We propose a model that separates the functions of KaiC ATPase and phosphorylation in cell division gating and other circadian behaviors.

Tobacco aquaporin NtAQP1 and human aquaporin hAQP1 contribute to single cell photosynthesis in Synechococcus.

BACKGROUND INFORMATION: Aquaporins are H2O-permeable membrane protein pores. However, some aquaporins are also permeable to other substances such as CO2. In higher plants, overexpression of such aquaporins has already led to an enhanced photosynthetic performance due to improved CO2 mesophyll conductance. In this work, we investigated the effects of such aquaporins on unicellular photosynthetically active organisms, specifically cyanobacteria. RESULTS: Overexpression of aquaporins NtAQP1 or hAQP1 that might have a function to improve CO2 membrane permeability lead to increased photosynthesis rates in the cyanobacterium Synechococcus sp. PCC7002 as concluded by the rate of evolved O2. A shift in the Plastoquinone pool state of the cells supports our findings. Water permeable aquaporins without CO2 permeability, such as NtPIP2;1, do not have this effect. CONCLUSIONS AND SIGNIFICANCE: We conclude that also in single cell organisms like cyanobacteria, membrane CO2 conductivity could be rate limiting and CO2-porins reduce the respective membrane resistance. We could show that besides the tobacco aquaporin NtAQP1 also the human hAQP1 most likely functions as CO2 diffusion facilitator in the photosynthesis assay.

Structural basis for enzymatic photocatalysis in chlorophyll biosynthesis.

The enzyme protochlorophyllide oxidoreductase (POR) catalyses a light-dependent step in chlorophyll biosynthesis that is essential to photosynthesis and, ultimately, all life on Earth1-3. POR, which is one of three known light-dependent enzymes4,5, catalyses reduction of the photosensitizer and substrate protochlorophyllide to form the pigment chlorophyllide. Despite its biological importance, the structural basis for POR photocatalysis has remained unknown. Here we report crystal structures of cyanobacterial PORs from Thermosynechococcus elongatus and Synechocystis sp. in their free forms, and in complex with the nicotinamide coenzyme. Our structural models and simulations of the ternary protochlorophyllide-NADPH-POR complex identify multiple interactions in the POR active site that are important for protochlorophyllide binding, photosensitization and photochemical conversion to chlorophyllide. We demonstrate the importance of active-site architecture and protochlorophyllide structure in driving POR photochemistry in experiments using POR variants and protochlorophyllide analogues. These studies reveal how the POR active site facilitates light-driven reduction of protochlorophyllide by localized hydride transfer from NADPH and long-range proton transfer along structurally defined proton-transfer pathways.

Rubisco condensate formation by CcmM in ß-carboxysome biogenesis.

Cells use compartmentalization of enzymes as a strategy to regulate metabolic pathways and increase their efficiency1. The α- and ß-carboxysomes of cyanobacteria contain ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco)-a complex of eight large (RbcL) and eight small (RbcS) subunits-and carbonic anhydrase2-4. As HCO3- can diffuse through the proteinaceous carboxysome shell but CO2 cannot5, carbonic anhydrase generates high concentrations of CO2 for carbon fixation by Rubisco6. The shell also prevents access to reducing agents, generating an oxidizing environment7-9. The formation of ß-carboxysomes involves the aggregation of Rubisco by the protein CcmM10, which exists in two forms: full-length CcmM (M58 in Synechococcus elongatus PCC7942), which contains a carbonic anhydrase-like domain8 followed by three Rubisco small subunit-like (SSUL) modules connected by flexible linkers; and M35, which lacks the carbonic anhydrase-like domain11. It has long been speculated that the SSUL modules interact with Rubisco by replacing RbcS2-4. Here we have reconstituted the Rubisco-CcmM complex and solved its structure. Contrary to expectation, the SSUL modules do not replace RbcS, but bind close to the equatorial region of Rubisco between RbcL dimers, linking Rubisco molecules and inducing phase separation into a liquid-like matrix. Disulfide bond formation in SSUL increases the network flexibility and is required for carboxysome function in vivo. Notably, the formation of the liquid-like condensate of Rubisco is mediated by dynamic interactions with the SSUL domains, rather than by low-complexity sequences, which typically mediate liquid-liquid phase separation in eukaryotes12,13. Indeed, within the pyrenoids of eukaryotic algae, the functional homologues of carboxysomes, Rubisco adopts a liquid-like state by interacting with the intrinsically disordered protein EPYC114. Understanding carboxysome biogenesis will be important for efforts to engineer CO2-concentrating mechanisms in plants15-19.

An Unstable Singularity Underlies Stochastic Phasing of the Circadian Clock in Individual Cyanobacterial Cells.

The endogenous circadian clock synchronizes with environmental time by appropriately resetting its phase in response to external cues. Of note, some resetting stimuli induce attenuated oscillations of clock output, which has been observed at the population-level in several organisms and in studies of individual humans. To investigate what is happening in individual cellular clocks, we studied the unicellular cyanobacterium S. elongatus. By measuring its phase-resetting responses to temperature changes, we found that population-level arrhythmicity occurs when certain perturbations cause stochastic phases of oscillations in individual cells. Combining modeling with experiments, we related stochastic phasing to the dynamical structure of the cyanobacterial clock as an oscillator and explored the physiological relevance of the oscillator structure for accurately timed rhythmicity in changing environmental conditions. Our findings and approach can be applied to other biological oscillators.

KidA, a multi-PAS domain protein, tunes the period of the cyanobacterial circadian oscillator.

The cyanobacterial clock presents a unique opportunity to understand the biochemical basis of circadian rhythms. The core oscillator, composed of the KaiA, KaiB, and KaiC proteins, has been extensively studied, but a complete picture of its connection to the physiology of the cell is lacking. To identify previously unknown components of the clock, we used KaiB locked in its active fold as bait in an immunoprecipitation/mass spectrometry approach. We found that the most abundant interactor, other than KaiC, was a putative diguanylate cyclase protein predicted to contain multiple Per-Arnt-Sim (PAS) domains, which we propose to name KidA. Here we show that KidA directly binds to the fold-switched active form of KaiB through its N-terminal PAS domains. We found that KidA shortens the period of the circadian clock both in vivo and in vitro and alters the ability of the clock to entrain to light-dark cycles. The dose-dependent effect of KidA on the clock period could be quantitatively recapitulated by a mathematical model in which KidA stabilizes the fold-switched form of KaiB, favoring rebinding to KaiC. Put together, our results show that the period and amplitude of the clock can be modulated by regulating the access of KaiB to the fold-switched form.

Identification and Functional Analysis of Lysine 2-Hydroxyisobutyrylation in Cyanobacteria.

Cyanobacteria are the oldest prokaryotic photoautotrophic microorganisms and have evolved complicated post-translational modification (PTM) machinery to respond to environmental stress. Lysine 2-hydroxyisobutyrylation (Khib) is a newly identified PTM that is reported to play important roles in diverse biological processes, however, its distribution and function in cyanobacteria have not been reported. Here, we performed the first systematic studies of Khib in a model cyanobacterium Synechococcus sp. strain PCC 7002 (Syn7002) using peptide prefractionation, pan-Khib antibody enrichment, and high-accuracy mass spectrometry (MS) analysis. A total of 1875 high-confidence Khib sites on 618 proteins were identified, and a large proportion of Khib sites are present on proteins in the cellular metabolism, protein synthesis, and photosynthesis pathways. Using site-directed mutagenesis and functional studies, we showed that Khib of glutaredoxin (Grx) affects the efficiency of the PS II reaction center and H2O2 resistance in Syn7002. Together, this study provides novel insights into the functions of Khib in cyanobacteria and suggests that reversible Khib may influence the stress response and photosynthesis in both cyanobacteria and plants.

Dynamics of the Synechococcus elongatus cytoskeletal GTPase FtsZ yields mechanistic and evolutionary insight into cyanobacterial and chloroplast FtsZs.

The division of cyanobacteria and their chloroplast descendants is orchestrated by filamenting temperature-sensitive Z (FtsZ), a cytoskeletal GTPase that polymerizes into protofilaments that form a "Z ring" at the division site. The Z ring has both a scaffolding function for division-complex assembly and a GTPase-dependent contractile function that drives cell or organelle constriction. A single FtsZ performs these functions in bacteria, whereas in chloroplasts, they are performed by two copolymerizing FtsZs, called AtFtsZ2 and AtFtsZ1 in Arabidopsis thaliana, which promote protofilament stability and dynamics, respectively. To probe the differences between cyanobacterial and chloroplast FtsZs, we used light scattering to characterize the in vitro protofilament dynamics of FtsZ from the cyanobacterium Synechococcus elongatus PCC 7942 (SeFtsZ) and investigate how coassembly of AtFtsZ2 or AtFtsZ1 with SeFtsZ influences overall dynamics. SeFtsZ protofilaments assembled rapidly and began disassembling before GTP depletion, whereas AtFtsZ2 protofilaments were far more stable, persisting beyond GTP depletion. Coassembled SeFtsZ-AtFtsZ2 protofilaments began disassembling before GTP depletion, similar to SeFtsZ. In contrast, AtFtsZ1 did not alter disassembly onset when coassembled with SeFtsZ, but fluorescence recovery after photobleaching analysis showed it increased the turnover of SeFtsZ subunits from SeFtsZ-AtFtsZ1 protofilaments, mirroring its effect upon coassembly with AtFtsZ2. Comparisons of our findings with previous work revealed consistent differences between cyanobacterial and chloroplast FtsZ dynamics and suggest that the scaffolding and dynamics-promoting functions were partially separated during evolution of two chloroplast FtsZs from their cyanobacterial predecessor. They also suggest that chloroplasts may have evolved a mechanism distinct from that in cyanobacteria for promoting FtsZ protofilament dynamics.

DnaK2 Mediates a Negative Feedback Regulation of the Heat Shock Responsive Hik2-Rre1 Two-Component System in the Cyanobacterium Synechococcus Elongatus PCC 7942.

The two-component system (TCS) is a conserved signal transduction module in bacteria. The Hik2-Rre1 system is responsible for transcriptional activation upon high-temperature shift as well as plastoquinone-related redox stress in the cyanobacterium Synechococcus elongatus PCC 7942. As heat-induced de novo protein synthesis was previously shown to be required to quench the heat-activated response, we investigated the underlying mechanism in this study. We found that the heat-inducible transcription activation was alleviated by the overexpression of dnaK2, which is an essential homolog of the highly conserved HSP70 chaperone and whose expression is induced under the control of the Hik2-Rre1 TCS. Phosphorylation of Rre1 correlated with transcription of the regulatory target hspA. The redox stress response was found to be similarly repressed by dnaK2 overexpression. Considered together with the previous information, we propose a negative feedback mechanism of the Hik2-Rre1-dependent stress response that maintains the cellular homeostasis mediated by DnaK2.