Heterocyst patterns without patterning proteins in cyanobacterial filaments

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.

Quantitative models of nitrogen-fixing organisms

Nitrogen-fixing organisms are of importance to the environment, providing bioavailable nitrogen to the biosphere. Quantitative models have been used to complement the laboratory experiments and in situ measurements, where such evaluations are difficult or costly. Here, we review the current state of the quantitative modeling of nitrogen-fixing organisms and ways to enhance the bridge between theoretical and empirical studies.

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.

An integrative approach for modeling and simulation of heterocyst pattern formation in cyanobacteria filaments

Heterocyst differentiation in cyanobacteria filaments is one of the simplest examples of cellular differentiation and pattern formation in multicellular organisms. Despite of the many experimental studies addressing the evolution and sustainment of heterocyst patterns and the knowledge of the genetic circuit underlying the behavior of single cyanobacterium under nitrogen deprivation, there is still a theoretical gap connecting these two macroscopic and microscopic processes. As an attempt to shed light on this issue, here we explore heterocyst differentiation under the paradigm of systems biology. This framework allows us to formulate the essential dynamical ingredients of the genetic circuit of a single cyanobacterium into a set of differential equations describing the time evolution of the concentrations of the relevant molecular products. As a result, we are able to study the behavior of a single cyanobacterium under different external conditions, emulating nitrogen deprivation, and simulate the dynamics of cyanobacteria filaments by coupling their respective genetic circuits via molecular diffusion. These two ingredients allow us to understand the principles by which heterocyst patterns can be generated and sustained. In particular, our results point out that, by including both diffusion and noisy external conditions in the computational model, it is possible to reproduce the main features of the formation and sustainment of heterocyst patterns in cyanobacteria filaments as observed experimentally. Finally, we discuss the validity and possible improvements of the model.

Maintenance of heterocyst patterning in a filamentous cyanobacterium

In the absence of sufficient combined nitrogen, some filamentous cyanobacteria differentiate nitrogen-fixing heterocysts at approximately every 10th cell position. As cells between heterocysts grow and divide, this initial pattern is maintained by the differentiation of a single cell approximately midway between existing heterocysts. This paper introduces a mathematical model for the maintenance of the periodic pattern of heterocysts differentiated by Anabaena sp. strain PCC 7120 based on the current experimental knowledge of the system. The model equations describe a non-diffusing activator (HetR) and two inhibitors (PatS and HetN) that undergo diffusion in a growing one-dimensional domain. The inhibitors in this model have distinct diffusion rates and temporal expression patterns. These unique aspects of the model reflect recent experimental findings regarding the molecular interactions that regulate patterning in Anabaena. Output from the model is in good agreement with both the temporal and spatial characteristics of the pattern maintenance process observed experimentally.

Cyanobacterial cell lineage analysis of the spatiotemporal hetR expression profile during heterocyst pattern formation in Anabaena sp. PCC 7120

Diazotrophic heterocyst formation in the filamentous cyanobacterium, Anabaena sp. PCC 7120, is one of the simplest pattern formations known to occur in cell differentiation. Most previous studies on heterocyst patterning were based on statistical analysis using cells collected or observed at different times from a liquid culture, which would mask stochastic fluctuations affecting the process of pattern formation dynamics in a single bacterial filament. In order to analyze the spatiotemporal dynamics of heterocyst formation at the single filament level, here we developed a culture system to monitor simultaneously bacterial development, gene expression, and phycobilisome fluorescence. We also developed micro-liquid chamber arrays to analyze multiple Anabaena filaments at the same time. Cell lineage analyses demonstrated that the initial distributions of hetR::gfp and phycobilisome fluorescence signals at nitrogen step-down were not correlated with the resulting distribution of developed heterocysts. Time-lapse observations also revealed a dynamic hetR expression profile at the single-filament level, including transient upregulation accompanying cell division, which did not always lead to heterocyst development. In addition, some cells differentiated into heterocysts without cell division after nitrogen step-down, suggesting that cell division in the mother cells is not an essential requirement for heterocyst differentiation.

Modeling heterocyst pattern formation in cyanobacteria

Background

To allow the survival of the population in the absence of nitrogen, some cyanobacteria strains have developed the capability of differentiating into nitrogen fixing cells, forming a characteristic pattern. In this paper, the process by which cyanobacteria differentiates from vegetative cells into heterocysts in the absence of nitrogen and the elements of the gene network involved that allow the formation of such a pattern are investigated.

Methods

A simple gene network model, which represents the complexity of the differentiation process, and the role of all variables involved in this cellular process is proposed. Specific characteristics and details of the system's behavior such as transcript profiles for ntcA, hetR and patS between consecutive heterocysts were studied.

Results

The proposed model is able to capture one of the most distinctive features of this system: a characteristic distance of 10 cells between two heterocysts, with a small standard deviation according to experimental variability. The system's response to knock-out and over-expression of patS and hetR was simulated in order to validate the proposed model against experimental observations. In all cases, simulations show good agreement with reported experimental results.

Conclusion

A simple evolution mathematical model based on the gene network involved in heterocyst differentiation was proposed. The behavior of the biological system naturally emerges from the network and the model is able to capture the spacing pattern observed in heterocyst differentiation, as well as the effect of external perturbations such as nitrogen deprivation, gene knock-out and over-expression without specific parameter fitting.