Spatial fluctuations in expression of the heterocyst differentiation regulatory gene hetR in Anabaena filaments

Growth couples temporal and spatial fluctuations of tissue properties during morphogenesis

Living tissues display fluctuations-random spatial and temporal variations of tissue properties around their reference values-at multiple scales. It is believed that such fluctuations may enable tissues to sense their state or their size. Recent theoretical studies developed specific models of fluctuations in growing tissues and predicted that fluctuations of growth show long-range correlations. Here, we elaborated upon these predictions and we tested them using experimental data. We first introduced a minimal model for the fluctuations of any quantity that has some level of temporal persistence or memory, such as concentration of a molecule, local growth rate, or mechanical property. We found that long-range correlations are generic, applying to any such quantity, and that growth couples temporal and spatial fluctuations, through a mechanism that we call "fluctuation stretching"-growth enlarges the length scale of variation of this quantity. We then analyzed growth data from sepals of the model plant Arabidopsis and we quantified spatial and temporal fluctuations of cell growth using the previously developed cellular Fourier transform. Growth appears to have long-range correlations. We compared different genotypes and growth conditions: mutants with lower or higher response to mechanical stress have lower temporal correlations and longer-range spatial correlations than wild-type plants. Finally, we used theoretical predictions to merge experimental data from all conditions and developmental stages into a unifying curve, validating the notion that temporal and spatial fluctuations are coupled by growth. Altogether, our work reveals kinematic constraints on spatiotemporal fluctuations that have an impact on the robustness of morphogenesis.

Spatio-temporal coherence of circadian clocks and temporal control of differentiation in Anabaena filaments

Circadian clock arrays in multicellular filaments of the heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120 display remarkable spatio-temporal coherence under nitrogen-replete conditions. To shed light on the interplay between circadian clocks and the formation of developmental patterns, we followed the expression of a clock-controlled gene under nitrogen deprivation, at the level of individual cells. Our experiments showed that differentiation into heterocysts took place preferentially within a limited interval of the circadian clock cycle, that gene expression in different vegetative intervals along a developed filament was discoordinated, and that the circadian clock was active in individual heterocysts. Furthermore, Anabaena mutants lacking the kaiABC genes encoding the circadian clock core components produced heterocysts but failed in diazotrophy. Therefore, genes related to some aspect of nitrogen fixation, rather than early or mid-heterocyst differentiation genes, are likely affected by the absence of the clock. A bioinformatics analysis supports the notion that RpaA may play a role as master regulator of clock outputs in Anabaena, the temporal control of differentiation by the circadian clock and the involvement of the clock in proper diazotrophic growth. Together, these results suggest that under nitrogen-deficient conditions, the clock coherent unit in Anabaena is reduced from a full filament under nitrogen-rich conditions to the vegetative cell interval between heterocysts.IMPORTANCECircadian clocks, from unicellular organisms to animals, temporally align biological processes to day and night cycles. We study the dynamics of a circadian clock-controlled gene at the individual cell level in the multicellular filamentous cyanobacterium Anabaena, under nitrogen-stress conditions. Under these conditions, some cells along filaments differentiate to carry out atmospheric nitrogen fixation and lose their capability for oxygenic photosynthesis. We found that clock synchronization is limited to organismic units of contiguous photosynthetic cells, contrary to nitrogen-replete conditions in which clocks are synchronized over a whole filament. We provided evidence that the circadian clock regulates the process of differentiation, allowing it to occur preferentially within a limited time window during the circadian clock period. Lastly, we present evidence that the signal from the core clock to clock-regulated genes is conveyed in Anabaena as in unicellular cyanobacteria.

Spatial self-organization of metabolism in microbial systems: A matter of enzymes and chemicals

Most bacteria live in dense, spatially structured communities such as biofilms. The high density allows cells to alter the local microenvironment, whereas the limited mobility can cause species to become spatially organized. Together, these factors can spatially organize metabolic processes within microbial communities so that cells in different locations perform different metabolic reactions. The overall metabolic activity of a community depends both on how metabolic reactions are arranged in space and on how they are coupled, i.e., how cells in different regions exchange metabolites. Here, we review mechanisms that lead to the spatial organization of metabolic processes in microbial systems. We discuss factors that determine the length scales over which metabolic activities are arranged in space and highlight how the spatial organization of metabolic processes affects the ecology and evolution of microbial communities. Finally, we define key open questions that we believe should be the main focus of future research.

Mathematical models of nitrogen-fixing cell patterns in filamentous cyanobacteria

The Anabaena genus is a model organism of filamentous cyanobacteria whose vegetative cells can differentiate under nitrogen-limited conditions into a type of cell called a heterocyst. These heterocysts lose the possibility to divide and are necessary for the filament because they can fix and share environmental nitrogen. In order to distribute the nitrogen efficiently, heterocysts are arranged to form a quasi-regular pattern whose features are maintained as the filament grows. Recent efforts have allowed advances in the understanding of the interactions and genetic mechanisms underlying this dynamic pattern. Here, we present a systematic review of the existing theoretical models of nitrogen-fixing cell differentiation in filamentous cyanobacteria. These filaments constitute one of the simplest forms of multicellular organization, and this allows for several modeling scales of this emergent pattern. The system has been approached at three different levels. From bigger to smaller scale, the system has been considered as follows: at the population level, by defining a mean-field simplified system to study the ratio of heterocysts and vegetative cells; at the filament level, with a continuous simplification as a reaction-diffusion system; and at the cellular level, by studying the genetic regulation that produces the patterning for each cell. In this review, we compare these different approaches noting both the virtues and shortcomings of each one of them.

Terminal heterocyst differentiation in the Anabaena patA mutant as a result of post-transcriptional modifications and molecular leakage

The Anabaena genus is a model organism of filamentous cyanobacteria whose vegetative cells can differentiate under nitrogen-limited conditions into a type of cell called heterocyst. These heterocysts lose the possibility to divide and are necessary for the colony because they can fix and share environmental nitrogen. In order to distribute the nitrogen efficiently, heterocysts are arranged to form a quasi-regular pattern whose features are maintained as the filament grows. Recent efforts have allowed advances in the understanding of the interactions and genetic mechanisms underlying this dynamic pattern. However, the main role of the patA and hetF genes are yet to be clarified; in particular, the patA mutant forms heterocysts almost exclusively in the terminal cells of the filament. In this work, we investigate the function of these genes and provide a theoretical model that explains how they interact within the broader genetic network, reproducing their knock-out phenotypes in several genetic backgrounds, including a nearly uniform concentration of HetR along the filament for the patA mutant. Our results suggest a role of hetF and patA in a post-transcriptional modification of HetR which is essential for its regulatory function. In addition, the existence of molecular leakage out of the filament in its boundary cells is enough to explain the preferential appearance of terminal heterocysts, without any need for a distinct regulatory pathway.

Understanding multicellular pattern formation is key for the study of both natural and synthetic developmental processes. Arguably one of the simplest model systems for this is the filamentous cyanobacterium Anabaena, that in conditions of nitrogen deprivation undergoes a dynamical differentiation process that differentiates roughly one in every ten cells into nitrogen-fixing heterocysts, in a quasi-regular pattern that is maintained as the filament keeps growing. One of the most characteristic mutations affecting this process forms heterocysts mostly constrained to the terminal cells of the filament. We have used experimental observations to propose a mathematical model of heterocyst differentiation able to reproduce this striking phenotype. The model extends our understanding of the regulations in this pattern-forming system and makes several predictions on molecular interactions. Importantly, a key aspect is the boundary condition at the filament’s ends: inhibitors of differentiation should be able to leak out of the filament, or otherwise the terminal cells would not differentiate. This highlights, in a very clear example, the importance of considering physical constraints in developmental processes.

The Making of a Heterocyst in Cyanobacteria

Heterocyst differentiation that occurs in some filamentous cyanobacteria, such as Anabaena sp. PCC 7120, provides a unique model for prokaryotic developmental biology. Heterocyst cells are formed in response to combined-nitrogen deprivation and possess a microoxic environment suitable for nitrogen fixation following extensive morphological and physiological reorganization. A filament of Anabaena is a true multicellular organism, as nitrogen and carbon sources are exchanged among different cells and cell types through septal junctions to ensure filament growth. Because heterocysts are terminally differentiated cells and unable to divide, their activity is an altruistic behavior dedicated to providing fixed nitrogen for neighboring vegetative cells. Heterocyst development is also a process of one-dimensional pattern formation, as heterocysts are semiregularly intercalated among vegetative cells. Morphogens form gradients along the filament and interact with each other in a fashion that fits well into the Turing model, a mathematical framework to explain biological pattern formation. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Robust, coherent, and synchronized circadian clock-controlled oscillations along Anabaena filaments

Circadian clocks display remarkable reliability despite significant stochasticity in biomolecular reactions. We study the dynamics of a circadian clock-controlled gene at the individual cell level in Anabaena sp. PCC 7120, a multicellular filamentous cyanobacterium. We found significant synchronization and spatial coherence along filaments, clock coupling due to cell-cell communication, and gating of the cell cycle. Furthermore, we observed low-amplitude circadian oscillatory transcription of kai genes encoding the post-transcriptional core oscillatory circuit and high-amplitude oscillations of rpaA coding for the master regulator transducing the core clock output. Transcriptional oscillations of rpaA suggest an additional level of regulation. A stochastic one-dimensional toy model of coupled clock cores and their phosphorylation states shows that demographic noise can seed stochastic oscillations outside the region where deterministic limit cycles with circadian periods occur. The model reproduces the observed spatio-temporal coherence along filaments and provides a robust description of coupled circadian clocks in a multicellular organism.

Impaired cell-cell communication in the multicellular cyanobacterium Anabaena affects carbon uptake, photosynthesis, and the cell wall

Cell-cell communication is an essential attribute of multicellular organisms. The effects of perturbed communication were studied in septal protein mutants of the heterocyst-forming filamentous cyanobacterium Anabaena sp. PCC 7120 model organism. Strains bearing sepJ and sepJ/fraC/fraD deletions showed differences in growth, pigment absorption spectra, and spatial patterns of expression of the hetR gene encoding a heterocyst differentiation master regulator. Global changes in gene expression resulting from deletion of those genes were mapped by RNA sequencing analysis of wild-type and mutant strains, both under nitrogen-replete and nitrogen-poor conditions. The effects of sepJ and fraC/fraD deletions were non-additive, and perturbed cell-cell communication led to significant changes in global gene expression. Most significant effects, related to carbon metabolism, included increased expression of genes encoding carbon uptake systems and components of the photosynthetic apparatus, as well as decreased expression of genes encoding cell wall components related to heterocyst differentiation and to polysaccharide export.

HetL, HetR and PatS form a reaction-diffusion system to control pattern formation in the cyanobacterium nostoc PCC 7120

Local activation and long-range inhibition are mechanisms conserved in self-organizing systems leading to biological patterns. A number of them involve the production by the developing cell of an inhibitory morphogen, but how this cell becomes immune to self-inhibition is rather unknown. Under combined nitrogen starvation, the multicellular cyanobacterium Nostoc PCC 7120 develops nitrogen-fixing heterocysts with a pattern of one heterocyst every 10-12 vegetative cells. Cell differentiation is regulated by HetR which activates the synthesis of its own inhibitory morphogens, diffusion of which establishes the differentiation pattern. Here, we show that HetR interacts with HetL at the same interface as PatS, and that this interaction is necessary to suppress inhibition and to differentiate heterocysts. hetL expression is induced under nitrogen-starvation and is activated by HetR, suggesting that HetL provides immunity to the heterocyst. This protective mechanism might be conserved in other differentiating cyanobacteria as HetL homologues are spread across the phylum.

Self-Organization in Pattern Formation

Self-organization is pervasive in development, from symmetry breaking in the early embryo to tissue patterning and morphogenesis. For a few model systems, the underlying molecular and cellular processes are now sufficiently characterized that mathematical models can be confronted with experiments, to explore the dynamics of pattern formation. Here, we review selected systems, ranging from cyanobacteria to mammals, where different forms of cell-cell communication, acting alone or together with positional cues, drive the patterning of cell fates, highlighting the insights that even very simple models can provide as well as the challenges on the path to a predictive understanding of development.

Pentapeptide-repeat, cytoplasmic-membrane protein HglK influences the septal junctions in the heterocystous cyanobacterium Anabaena

N₂ -fixing heterocystous cyanobacteria grow as chains of cells that are connected by proteinaceous septal junctions, which traverse the septal peptidoglycan through nanopores and mediate intercellular molecular transfer. In the model organism Anabaena sp. strain PCC 7120, proteins SepJ, FraC and FraD, which are localized at the cell poles in the intercellular septa, are needed to produce septal junctions. The pentapeptide-repeat, membrane-spanning protein HglK has been described to be involved in the deposition of the heterocyst-specific glycolipid layer, but the hglK mutant also showed intercellular septa broader than in the wild type. Here we found that hglK mutant of Anabaena is impaired in the expression of heterocyst-related genes coxB2A2C2 (cytochrome c oxidase) and nifHDK (nitrogenase), indicating a defect in heterocyst differentiation. HglK was predominantly localized at the intercellular septa and was required to make long filaments, produce a normal number of nanopores and express full intercellular molecular transfer activity. However, the effects of hglK inactivation were not additive to those of the inactivation of sepJ and/or fraC-fraD. We suggest that HglK contributes to the architecture of the intercellular septa with an impact on the function of septal junctions.

Cyanobacterial Septal Junctions: Properties and Regulation

Heterocyst-forming cyanobacteria are multicellular organisms that grow as chains of cells (filaments or trichomes) in which the cells exchange regulators and nutrients. In this article, we review the morphological, physiological and genetic data that have led to our current understanding of intercellular communication in these organisms. Intercellular molecular exchange appears to take place by simple diffusion through proteinaceous structures, known as septal junctions, which connect the adjacent cells in the filament and traverse the septal peptidoglycan through perforations known as nanopores. Proteins that are necessary to produce, and that may be components of, the septal junctions-SepJ, FraC and FraD-have been identified in the heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120 model. Additionally, several proteins that are necessary to produce a normal number of nanopores and functional septal junctions have been identified, including AmiC-type amidases, peptidoglycan-binding proteins and some membrane transporters. Available reports and reevaluation of intercellular molecular transfer data for some mutants of Anabaena suggest that the septal junctions can be regulated, likely by a mechanism of gating.

Noise⁻Seeded Developmental Pattern Formation in Filamentous Cyanobacteria

Under nitrogen-poor conditions, multicellular cyanobacteria such as Anabaena sp. PCC 7120 undergo a process of differentiation, forming nearly regular, developmental patterns of individual nitrogen-fixing cells, called heterocysts, interspersed between intervals of vegetative cells that carry out photosynthesis. Developmental pattern formation is mediated by morphogen species that can act as activators and inhibitors, some of which can diffuse along filaments. We survey the limitations of the classical, deterministic Turing mechanism that has been often invoked to explain pattern formation in these systems, and then, focusing on a simpler system governed by birth-death processes, we illustrate pedagogically a recently proposed paradigm that provides a much more robust description of pattern formation: stochastic Turing patterns. We emphasize the essential role that cell-to-cell differences in molecular numbers—caused by inevitable fluctuations in gene expression—play, the so called demographic noise, in seeding the formation of stochastic Turing patterns over a much larger region of parameter space, compared to their deterministic counterparts.

Specific mutations in the permease domain of septal protein SepJ differentially affect functions related to multicellularity in the filamentous cyanobacterium Anabaena

Filamentous, heterocyst-forming cyanobacteria are multicellular organisms in which growth requires the activity of two interdependent cell types that exchange nutrients and regulators. Vegetative cells provide heterocysts with reduced carbon, and heterocysts provide vegetative cells with fixed nitrogen. Additionally, heterocyst differentiation from vegetative cells is regulated by inhibitors of differentiation produced by prospective heterocysts and heterocysts. Proteinaceous structures known as septal junctions join the cells in the filament. The SepJ protein is involved in formation of septal junctions in the model heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120. SepJ bears extra-membrane and membrane (permease) domains and is located at the cell poles in the intercellular septa of the filament. Here we created Anabaena mutants that produce SepJ proteins altered in the permease domain. Some of these mutant SepJ proteins did not provide functions needed for Anabaena to form long filaments and (in some cases) differentiate heterocysts, identifying amino acids and amino acid stretches that are important for the structure or function of the protein. Some other mutant SepJ proteins fulfilled filamentation and heterocyst differentiation functions but failed to provide normal communication function assessed via the intercellular transfer of the fluorescent marker calcein. These mutant SepJ proteins bore mutations in amino acids located at the cytoplasmic face of the permease, which could affect access of the fluorescent marker to the septal junctions. Overall, the data are consistent with the idea that SepJ carries out multiple roles in the multicellular function of the Anabaena filament.

The Ser/Thr Kinase PknH Is Essential for Maintaining Heterocyst Pattern in the Cyanobacterium Anabaena sp. Strain PCC 7120

In the filamentous cyanobacterium Anabaena sp. strain, PCC 7120, heterocysts (which are nitrogen-fixing cells) are formed in the absence of combined nitrogen in the medium. Heterocysts are separated from one another by 10 to 15 vegetative cells along the filaments, which consist of a few hundred of cells. hetR is necessary for heterocyst differentiation; and patS and hetN, expressed in heterocysts, play important roles in heterocyst pattern formation by laterally inhibiting the expression of hetR in adjacent cells. The results of this study indicated that pknH, which encodes a Ser/Thr kinase, was also involved in heterocyst pattern formation. In the pknH mutant, the heterocyst pattern was normal within 24 h after nitrogen deprivation, but multiple contiguous heterocysts were formed from 24 to 48 h. A time-lapse analysis of reporter strains harboring a fusion between gfp and the hetR promoter indicated that pknH was required to suppress hetR expression in cells adjacent to the preexisting heterocysts. These results indicated that pknH was necessary for the lateral inhibition of heterocyst differentiation to maintain the heterocyst pattern.

Soft X-Ray Imaging of Cellular Carbon and Nitrogen Distributions in Heterocystous Cyanobacteria

Soft x-ray microscopy (SXM) is a minimally invasive technique for single-cell high-resolution imaging as well as the visualization of intracellular distributions of light elements such as carbon, nitrogen, and oxygen. We used SXM to observe photosynthesis and nitrogen fixation in the filamentous cyanobacterium Anabaena sp. PCC 7120, which can form heterocysts during nitrogen starvation. Statistical and spectroscopic analyses from SXM images around the K-absorption edge of nitrogen revealed a significant difference in the carbon-to-nitrogen (C/N) ratio between vegetative cells and heterocysts. Application of this analysis to soft x-ray images of Anabaena sp. PCC 7120 revealed inhomogenous C/N ratios in the cells. Furthermore, soft x-ray tomography of Anabaena sp. PCC 7120 revealed differing cellular C/N ratios, indicating different carbon and nitrogen distributions between vegetative cells and heterocysts in three dimensions.

Robust stochastic Turing patterns in the development of a one-dimensional cyanobacterial organism

Under nitrogen deprivation, the one-dimensional cyanobacterial organism Anabaena sp. PCC 7120 develops patterns of single, nitrogen-fixing cells separated by nearly regular intervals of photosynthetic vegetative cells. We study a minimal, stochastic model of developmental patterns in Anabaena that includes a nondiffusing activator, two diffusing inhibitor morphogens, demographic fluctuations in the number of morphogen molecules, and filament growth. By tracking developing filaments, we provide experimental evidence for different spatiotemporal roles of the two inhibitors during pattern maintenance and for small molecular copy numbers, justifying a stochastic approach. In the deterministic limit, the model yields Turing patterns within a region of parameter space that shrinks markedly as the inhibitor diffusivities become equal. Transient, noise-driven, stochastic Turing patterns are produced outside this region, which can then be fixed by downstream genetic commitment pathways, dramatically enhancing the robustness of pattern formation, also in the biologically relevant situation in which the inhibitors' diffusivities may be comparable.

Overexpression of SepJ alters septal morphology and heterocyst pattern regulated by diffusible signals in Anabaena

Filamentous, N2 -fixing, heterocyst-forming cyanobacteria grow as chains of cells that are connected by septal junctions. In the model organism Anabaena sp. strain PCC 7120, the septal protein SepJ is required for filament integrity, normal intercellular molecular exchange, heterocyst differentiation, and diazotrophic growth. An Anabaena strain overexpressing SepJ made wider septa between vegetative cells than the wild type, which correlated with a more spread location of SepJ in the septa as observed with a SepJ-GFP fusion, and contained an increased number of nanopores, the septal peptidoglycan perforations that likely accommodate septal junctions. The septa between heterocysts and vegetative cells, which are narrow in wild-type Anabaena, were notably enlarged in the SepJ-overexpressing mutant. Intercellular molecular exchange tested with fluorescent tracers was increased for the SepJ-overexpressing strain specifically in the case of calcein transfer between vegetative cells and heterocysts. These results support an association between calcein transfer, SepJ-related septal junctions, and septal peptidoglycan nanopores. Under nitrogen deprivation, the SepJ-overexpressing strain produced an increased number of contiguous heterocysts but a decreased percentage of total heterocysts. These effects were lost or altered in patS and hetN mutant backgrounds, supporting a role of SepJ in the intercellular transfer of regulatory signals for heterocyst differentiation.

Formation and maintenance of nitrogen-fixing cell patterns in filamentous cyanobacteria

Cyanobacteria forming one-dimensional filaments are paradigmatic model organisms of the transition between unicellular and multicellular living forms. Under nitrogen-limiting conditions, in filaments of the genus Anabaena, some cells differentiate into heterocysts, which lose the possibility to divide but are able to fix environmental nitrogen for the colony. These heterocysts form a quasiregular pattern in the filament, representing a prototype of patterning and morphogenesis in prokaryotes. Recent years have seen advances in the identification of the molecular mechanism regulating this pattern. We use these data to build a theory on heterocyst pattern formation, for which both genetic regulation and the effects of cell division and filament growth are key components. The theory is based on the interplay of three generic mechanisms: local autoactivation, early long-range inhibition, and late long-range inhibition. These mechanisms can be identified with the dynamics of hetR, patS, and hetN expression. Our theory reproduces quantitatively the experimental dynamics of pattern formation and maintenance for wild type and mutants. We find that hetN alone is not enough to play the role as the late inhibitory mechanism: a second mechanism, hypothetically the products of nitrogen fixation supplied by heterocysts, must also play a role in late long-range inhibition. The preponderance of even intervals between heterocysts arises naturally as a result of the interplay between the timescales of genetic regulation and cell division. We also find that a purely stochastic initiation of the pattern, without a two-stage process, is enough to reproduce experimental observations.

Structural insights into HetR-PatS interaction involved in cyanobacterial pattern formation

The one-dimensional pattern of heterocyst in the model cyanobacterium Anabaena sp. PCC 7120 is coordinated by the transcription factor HetR and PatS peptide. Here we report the complex structures of HetR binding to DNA, and its hood domain (HetRHood) binding to a PatS-derived hexapeptide (PatS6) at 2.80 and 2.10 Å, respectively. The intertwined HetR dimer possesses a couple of novel HTH motifs, each of which consists of two canonical α-helices in the DNA-binding domain and an auxiliary α-helix from the flap domain of the neighboring subunit. Two PatS6 peptides bind to the lateral clefts of HetRHood, and trigger significant conformational changes of the flap domain, resulting in dissociation of the auxiliary α-helix and eventually release of HetR from the DNA major grove. These findings provide the structural insights into a prokaryotic example of Turing model.