The heterocyst differentiation transcriptional regulator HetR of the filamentous cyanobacteriumAnabaenaforms tetramers and can be regulated by phosphorylation

Expanding the FurC (PerR) regulon in Anabaena (Nostoc) sp. PCC 7120: Genome-wide identification of novel direct targets uncovers FurC participation in central carbon metabolism regulation

FurC (PerR, Peroxide Response Regulator) from Anabaena sp. PCC 7120 (also known as Nostoc sp. PCC 7120) is a master regulator engaged in the modulation of relevant processes including the response to oxidative stress, photosynthesis and nitrogen fixation. Previous differential gene expression analysis of a furC-overexpressing strain (EB2770FurC) allowed the inference of a putative FurC DNA-binding consensus sequence. In the present work, more data concerning the regulon of the FurC protein were obtained through the searching of the putative FurC-box in the whole Anabaena sp. PCC 7120 genome. The total amount of novel FurC-DNA binding sites found in the promoter regions of genes with known function was validated by electrophoretic mobility shift assays (EMSA) identifying 22 new FurC targets. Some of these identified targets display relevant roles in nitrogen fixation (hetR and hgdC) and carbon assimilation processes (cmpR, glgP1 and opcA), suggesting that FurC could be an additional player for the harmonization of carbon and nitrogen metabolisms. Moreover, differential gene expression of a selection of newly identified FurC targets was measured by Real Time RT-PCR in the furC-overexpressing strain (EB2770FurC) comparing to Anabaena sp. PCC 7120 revealing that in most of these cases FurC could act as a transcriptional activator.

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.

Quantitative Proteomics Reveals the Protein Regulatory Network of Anabaena sp. PCC 7120 under Nitrogen Deficiency

Anabaena sp. PCC 7120 (Anabaena 7120) is a photoautotrophic filamentous cyanobacterium capable of fixing atmospheric nitrogen. It is a model organism used for studying cell differentiation and nitrogen fixation. Under nitrogen deficiency, Anabaena 7120 forms specialized heterocysts capable of nitrogen fixation. However, the molecular mechanisms involved in the cyanobacterial adaptation to nitrogen deficiency are not well understood. Here, we employed a label-free quantitative proteomic strategy to systematically investigate the nitrogen deficiency response of Anabaena 7120 at different time points. In total, 363, 603, and 669 proteins showed significant changes in protein abundance under nitrogen deficiency for 3, 12, and 24 h, respectively. With mapping onto metabolic pathways, we revealed proteomic perturbation and regulation of carbon and nitrogen metabolism in response to nitrogen deficiency. Functional analysis confirmed the involvement of nitrogen stress-responsive proteins in biological processes, including nitrogen fixation, photosynthesis, energy and carbon metabolism, and heterocyst development. The expression of 10 proteins at different time points was further validated by using multiple reaction monitoring assays. In particular, many dysregulated proteins were found to be time-specific and involved in heterocyst development, providing new candidates for future functional studies in this model cyanobacterium. These results provide novel insights into the molecular mechanisms of nitrogen stress responses and heterocyst development in Anabaena 7120.

FurC (PerR) from Anabaena sp. PCC7120: a versatile transcriptional regulator engaged in the regulatory network of heterocyst development and nitrogen fixation

FurC (PerR) from Anabaena sp. PCC7120 was previously described as a key transcriptional regulator involved in setting off the oxidative stress response. In the last years, the cross-talk between oxidative stress, iron homeostasis and nitrogen metabolism is becoming more and more evident. In this work, the transcriptome of a furC-overexpressing strain was compared with that of a wild-type strain under both standard and nitrogen-deficiency conditions. The results showed that the overexpression of furC deregulates genes involved in several categories standing out photosynthesis, iron transport and nitrogen metabolism. The novel FurC-direct targets included some regulatory elements that control heterocyst development (hetZ and asr1734), genes directly involved in the heterocyst envelope formation (devBCA and hepC) and genes which participate in the nitrogen fixation process (nifHDK and nifH2, rbrA rubrerythrin and xisHI excisionase). Likewise, furC overexpression notably impacts the mRNA levels of patA encoding a key protein in the heterocyst pattern formation. The relevance of FurC in these processes is bringing out by the fact that the overexpression of furC impairs heterocyst development and cell growth under nitrogen step-down conditions. In summary, this work reveals a new player in the complex regulatory network of heterocyst formation and nitrogen fixation.

Inactivation of Three RG(S/T)GR Pentapeptide-Containing Negative Regulators of HetR Results in Lethal Differentiation of Anabaena PCC 7120

The filamentous cyanobacterium Anabaena sp. PCC 7120 produces, during the differentiation of heterocysts, a short peptide PatS and a protein HetN, both containing an RGSGR pentapeptide essential for activity. Both act on the master regulator HetR to guide heterocyst pattern formation by controlling the binding of HetR to DNA and its turnover. A third small protein, PatX, with an RG(S/T)GR motif is present in all HetR-containing cyanobacteria. In a nitrogen-depleted medium, inactivation of patX does not produce a discernible change in phenotype, but its overexpression blocks heterocyst formation. Mutational analysis revealed that PatX is not required for normal intercellular signaling, but it nonetheless is required when PatS is absent to prevent rapid ectopic differentiation. Deprivation of all three negative regulators—PatS, PatX, and HetN—resulted in synchronous differentiation. However, in a nitrogen-containing medium, such deprivation leads to extensive fragmentation, cell lysis, and aberrant differentiation, while either PatX or PatS as the sole HetR regulator can establish and maintain a semiregular heterocyst pattern. These results suggest that tight control over HetR by PatS and PatX is needed to sustain vegetative growth and regulated development. The mutational analysis has been interpreted in light of the opposing roles of negative regulators of HetR and the positive regulator HetL.

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.

The Pkn22 Kinase of Nostoc PCC 7120 Is Required for Cell Differentiation via the Phosphorylation of HetR on a Residue Highly Conserved in Genomes of Heterocyst-Forming Cyanobacteria

Hanks-type kinases encoding genes are present in most cyanobacterial genomes. Despite their widespread pattern of conservation, little is known so far about their role because their substrates and the conditions triggering their activation are poorly known. Here we report that under diazotrophic conditions, normal heterocyst differentiation and growth of the filamentous cyanobacterium Nostoc PCC 7120 require the presence of the Pkn22 kinase, which is induced under combined nitrogen starvation conditions. By analyzing the phenotype of pkn22 mutant overexpressing genes belonging to the regulatory cascade initiating the development program, an epistatic relationship was found to exist between this kinase and the master regulator of differentiation, HetR. The results obtained using a bacterial two hybrid approach indicated that Pkn22 and HetR interact, and the use of a genetic screen inducing the loss of this interaction showed that residues of HetR which are essential for this interaction to occur are also crucial to HetR activity both in vitro and in vivo. Mass spectrometry showed that HetR co-produced with the Pkn22 kinase in Escherichia coli is phosphorylated on Serine 130 residue. Phosphoablative substitution of this residue impaired the ability of the strain to undergo cell differentiation, while its phosphomimetic substitution increased the number of heterocysts formed. The Serine 130 residue is part of a highly conserved sequence in filamentous cyanobacterial strains differentiating heterocysts. Heterologous complementation assays showed that the presence of this domain is necessary for heterocyst induction. We propose that the phosphorylation of HetR might have been acquired to control heterocyst differentiation.

cyAbrB Transcriptional Regulators as Safety Devices To Inhibit Heterocyst Differentiation in Anabaena sp. Strain PCC 7120

Cyanobacteria are monophyletic organisms that perform oxygenic photosynthesis. While they exhibit great diversity, they have a common set of genes. However, the essentiality of them for viability has hampered the elucidation of their functions. One example of these genes is cyabrB1 (also known as calA in Anabaena sp. strain PCC 7120), encoding a transcriptional regulator. In the present study, we investigated the function of calA/cyabrB1 in the heterocyst-forming cyanobacterium Anabaena sp. PCC 7120 through CRISPR interference, a method that we recently utilized for the photosynthetic production of a useful chemical in this strain. Conditional knockdown of calA/cyabrB1 in the presence of nitrate resulted in the formation of heterocysts. Two genes, hetP and hepA, which are required for heterocyst formation, were upregulated by calA/cyabrB1 knockdown in the presence of combined nitrogen sources. These genes are known to be induced by HetR, a master regulator of heterocyst formation. hetR was not induced by calA/cyabrB1 knockdown. hetP and hepA were repressed by direct binding of CalA/cyAbrB1 to their promoter regions in a HetR-independent manner. In addition, the overexpression of calA/cyabrB1 abolished heterocyst formation upon nitrogen depletion. Also, knockout of calB/cyabrB2 (a paralogue gene of calA/cyabrB1), in addition to knockdown of calA/cyabrB1, enhanced heterocyst formation in the presence of nitrate, suggesting functional redundancy of cyAbrB proteins. We propose that a balance between amounts of HetR and CalA/cyAbrB1 is a key factor influencing heterocyst differentiation during nitrogen stepdown. We concluded that cyAbrB proteins are essential safety devices that inhibit heterocyst differentiation.IMPORTANCE Spore formation in Bacillus subtilis and Streptomyces has been extensively studied as models of prokaryotic nonterminal cell differentiation. In these organisms, many cells/hyphae differentiate simultaneously, which is governed by a network in which one regulator stands at the top. Differentiation of heterocysts in Anabaena sp. strain PCC 7120 is unique because it is terminal, and only 5 to 10% of vegetative cells differentiate into heterocysts. In this study, we identified CalA/cyAbrB1 as a repressor of two genes that are essential for heterocyst formation independently of HetR, a master activator for heterocyst differentiation. This finding is reasonable for unique cell differentiation of Anabaena because CalA/cyAbrB1 could suppress heterocyst differentiation tightly in vegetative cells, while only cells in which HetR is overexpressed could differentiate into heterocysts.

Ancient association of cyanobacterial multicellularity with the regulator HetR and an RGSGR pentapeptide-containing protein (PatX)

One simple model to explain biological pattern postulates the existence of a stationary regulator of differentiation that positively affects its own expression, coupled with a diffusible suppressor of differentiation that inhibits the regulator's expression. The first has been identified in the filamentous, heterocyst-forming cyanobacterium, Anabaena PCC 7120 as the transcriptional regulator, HetR and the second as the small protein, PatS, which contains a critical RGSGR motif that binds to HetR. HetR is present in almost all filamentous cyanobacteria, but only a subset of heterocyst-forming strains carry proteins similar to PatS. We identified a third protein, PatX that also carries the RGSGR motif and is coextensive with HetR. Amino acid sequences of PatX contain two conserved regions: the RGSGR motif and a hydrophobic N-terminus. Within 69 nt upstream from all instances of the gene is a DIF1 motif correlated in Anabaena with promoter induction in developing heterocysts, preceded in heterocyst-forming strains by an apparent NtcA-binding site, associated with regulation by nitrogen-status. Consistent with a role in the simple model, PatX is expressed dependent on HetR and acts to inhibit differentiation. The acquisition of the PatX/HetR pair preceded the appearance of both PatS and heterocysts, dating back to the beginnings of multicellularity.

Genetic responses to carbon and nitrogen availability in Anabaena

Heterocyst-forming cyanobacteria are filamentous organisms that perform oxygenic photosynthesis and CO₂ fixation in vegetative cells and nitrogen fixation in heterocysts, which are formed under deprivation of combined nitrogen. These organisms can acclimate to use different sources of nitrogen and respond to different levels of CO₂ . Following work mainly done with the best studied heterocyst-forming cyanobacterium, Anabaena, here we summarize the mechanisms of assimilation of ammonium, nitrate, urea and N₂ , the latter involving heterocyst differentiation, and describe aspects of CO₂ assimilation that involves a carbon concentration mechanism. These processes are subjected to regulation establishing a hierarchy in the assimilation of nitrogen sources -with preference for the most reduced nitrogen forms- and a dependence on sufficient carbon. This regulation largely takes place at the level of gene expression and is exerted by a variety of transcription factors, including global and pathway-specific transcriptional regulators. NtcA is a CRP-family protein that adjusts global gene expression in response to the C-to-N balance in the cells, and PacR is a LysR-family transcriptional regulator (LTTR) that extensively acclimates the cells to oxygenic phototrophy. A cyanobacterial-specific transcription factor, HetR, is involved in heterocyst differentiation, and other LTTR factors are specifically involved in nitrate and CO₂ assimilation.

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.

Functional Overlap of hetP and hetZ in Regulation of Heterocyst Differentiation in Anabaena sp. Strain PCC 7120

HetR plays a key role in regulation of heterocyst differentiation and patterning in Anabaena It directly regulates genes involved in heterocyst differentiation (such as hetP and hetZ), genes involved in pattern formation (patA), and many others. In this study, we investigated the functional relationship of hetP and hetZ and their role in HetR-dependent cell differentiation. Coexpression of hetP and hetZ from the promoter of ntcA, which encodes the global nitrogen regulator, enabled a hetR mutant to form heterocysts with low aerobic nitrogenase activity. Overexpression of hetZ restored heterocyst differentiation in a hetP mutant and vice versa. Overexpression of hetR restored heterocyst formation in either a hetP or a hetZ mutant but not in a hetZ hetP double mutant. The functional overlap of hetP and hetZ was further confirmed by reverse transcription-quantitative PCR (RT-qPCR) and transcriptomic analyses of their effects on gene expression. In addition, yeast two-hybrid and pulldown assays showed the interaction of HetZ with HetR. HetP and HetZ are proposed as the two major factors that control heterocyst formation in response to upregulation of hetRIMPORTANCE Heterocyst-forming cyanobacteria contribute significantly to N₂ fixation in marine, freshwater, and terrestrial ecosystems. Formation of heterocysts enables this group of cyanobacteria to fix N₂ efficiently under aerobic conditions. HetR, HetP, and HetZ are among the most important factors involved in heterocyst differentiation. We present evidence for the functional overlap of hetP and hetZ and for the key role of the HetR-HetP/HetZ circuit in regulation of heterocyst differentiation. The regulatory mechanism based on HetR, HetP, and HetZ is probably conserved in all heterocyst-forming cyanobacteria.

Processing of PatS, a morphogen precursor, in cell extracts of Anabaena sp. PCC 7120

Upon N-stepdown, Anabaena sp. PCC 7120 differentiates heterocysts along filaments in a semiregular pattern. A 17-amino acid peptide called PatS is a morphogen precursor for pattern formation. The principal PatS derivative involved in heterocyst patterning has been proposed to be the C-terminal peptide PatS-5 (RGSGR), PatS-6 (ERGSGR), or PatS-8 (CDERGSGR). We present the first evidence for processing of PatS in cell extracts of this cyanobacterium. PatS is probably cleaved between the C-terminal 7th and 8th amino acid residues, producing PatS-7 (DERGSGR), then converted into PatS-6 and PatS-5. The processing site could be changed by a substitution at the C-terminal 8th residue.

The heterocyst regulatory protein HetP and its homologs modulate heterocyst commitment in Anabaena sp. strain PCC 7120

The commitment of differentiating cells to a specialized fate is fundamental to the correct assembly of tissues within a multicellular organism. Because commitment is often irreversible, entry into and progression through this phase of development must be tightly regulated. Under nitrogen-limiting conditions, the multicellular cyanobacterium Anabaena sp. strain PCC 7120 terminally commits ∼10% of its cells to become specialized nitrogen-fixing heterocysts. Although commitment is known to occur 9-14 h after the induction of differentiation, the factors that regulate the initiation and duration of this phase have yet to be elucidated. Here, we report the identification of four genes that share a functional domain and modulate heterocyst commitment: hetP (alr2818), asl1930, alr2902, and alr3234 Epistatic relationships between all four genes relating to commitment were revealed by deleting them individually and in combination; asl1930 and alr3234 acted most upstream to delay commitment, alr2902 acted next in the pathway to inhibit development, and hetP acted most downstream to drive commitment forward. Possible protein-protein interactions between HetP, its homologs, and the heterocyst master regulator, HetR, were assessed, and interaction partners were defined. Finally, patterns of gene expression for each homolog, as determined by promoter fusions to gfp and reverse transcription-quantitative PCR, were distinct from that of hetP in both spatiotemporal organization and regulation. We posit that a dynamic succession of protein-protein interactions modulates the timing and efficiency of the commitment phase of development and note that this work highlights the utility of a multicellular cyanobacterium as a model for the study of developmental processes.

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.