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

PatU3 plays a central role in coordinating cell division and differentiation in pattern formation of filamentous cyanobacterium Nostoc sp. PCC 7120

Spatial periodic signal for cell differentiation in some multicellular organisms is generated according to Turing's principle for pattern formation. How a dividing cell responds to the signal of differentiation is addressed with the filamentous cyanobacterium Nostoc sp. PCC 7120, which forms the patterned distribution of heterocysts. We show that differentiation of a dividing cell was delayed until its division was completed and only one daughter cell became heterocyst. A mutant of patU3, which encodes an inhibitor of heterocyst formation, showed no such delay and formed heterocyst pairs from the daughter cells of cell division or dumbbell-shaped heterocysts from the cells undergoing cytokinesis. The patA mutant, which forms heterocysts only at the filament ends, restored intercalary heterocysts by a single nucleotide mutation of patU3, and double mutants of patU3/patA and patU3/hetF had the phenotypes of the patU3 mutant. We provide evidence that HetF, which can degrade PatU3, is recruited to cell divisome through its C-terminal domain. A HetF mutant with its N-terminal peptidase domain but lacking the C-terminal domain could not prevent the formation of heterocyst pairs, suggesting that the divisome recruitment of HetF is needed to sequester HetF for the delay of differentiation in dividing cells. Our study demonstrates that PatU3 plays a key role in cell-division coupled control of differentiation.

HetF defines a transition point from commitment to morphogenesis during heterocyst differentiation in the cyanobacterium Anabaena sp. PCC 7120

The filamentous cyanobacterium Anabaena sp. PCC 7120 is able to form heterocysts for nitrogen fixation. Heterocyst differentiation is initiated by combined-nitrogen deprivation, followed by the commitment step during which the developmental process becomes irreversible. Mature heterocysts are terminally differentiated cells unable to divide, and cell division is required for heterocyst differentiation. Previously, we have shown that the HetF protease regulates cell division and heterocyst differentiation by cleaving PatU3, which is an inhibitor for both events. When hetF is required during the developmental program remains unknown. Here, by controlling the timing of hetF expression during heterocyst differentiation, we provide evidence that hetF is required just before the beginning of heterocyst morphogenesis. Consistent with this finding, transcriptome data show that most of the genes known to be involved in the early step (such as hetR and ntcA) or the commitment step (such as hetP and hetZ) of heterocyst development could be expressed in the ΔhetF mutant. In contrast, most of the genes involved in heterocyst morphogenesis and nitrogen fixation remain repressed in the mutant. These results indicated that in the absence of hetF, heterocyst differentiation is able to be initiated and proceeds to the stage just before heterocyst envelope formation.

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.

Interaction network among factors involved in heterocyst-patterning in cyanobacteria

The genetically regulated pattern of heterocyst formation in multicellular cyanobacteria represents the simplest model to address how patterns emerge and are established, the signals that control them, and the regulatory pathways that act downstream. Although numerous factors involved in this process have been identified, the mechanisms of action of many of them remain largely unknown. The aim of this study was to identify specific relationships between 14 factors required for cell differentiation and pattern formation by exploring their putative physical interactions in the cyanobacterium model Nostoc sp. PCC 7120 and by probing their evolutionary conservation and distribution across the cyanobacterial phylum. A bacterial two-hybrid assay indicated that 10 of the 14 factors studied here are engaged in more than one protein-protein interaction. The transcriptional regulator PatB was central in this network as it showed the highest number of binary interactions. A phylum-wide genomic survey of the distribution of these factors in cyanobacteria showed that they are all highly conserved in the genomes of heterocyst-forming strains, with the PatN protein being almost restricted to this clade. Interestingly, eight of the factors that were shown to be capable of protein interactions were identified as key elements in the evolutionary genomics analysis. These data suggest that a network of 12 proteins may play a crucial role in heterocyst development and patterning. Unraveling the physical and functional interactions between these factors during heterocyst development will certainly shed light on the mechanisms underlying pattern establishment in cyanobacteria.

Heterocyst Development and Diazotrophic Growth of Anabaena variabilis under Different Nitrogen Availability

Nitrogen is globally limiting primary production in the ocean, but some species of cyanobacteria can carry out nitrogen (N) fixation using specialized cells known as heterocysts. However, the effect of N sources and their availability on heterocyst development is not yet fully understood. This study aimed to evaluate the effect of various inorganic N sources on the heterocyst development and cellular growth in an N-fixing cyanobacterium, Anabaena variabilis. Growth rate, heterocyst development, and cellular N content of the cyanobacteria were examined under varying nitrate and ammonium concentrations. A. variabilis exhibited high growth rate both in the presence and absence of N sources regardless of their concentration. Ammonium was the primary source of N in A. variabilis. Even the highest concentrations of both nitrate (1.5 g L-1 as NaNO3) and ammonium (0.006 g L-1 as Fe-NH4-citrate) did not exhibit an inhibitory effect on heterocyst development. Heterocyst production positively correlated with the cell N quota and negatively correlated with vegetative cell growth, indicating that both of the processes were interdependent. Taken together, N deprivation triggers heterocyst production for N fixation. This study outlines the difference in heterocyst development and growth in A. variabilis under different N sources.

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.

Expression from DIF1-motif promoters of hetR and patS is dependent on HetZ and modulated by PatU3 during heterocyst differentiation

HetR and PatS/PatX-derived peptides are the activator and diffusible inhibitor for cell differentiation and patterning in heterocyst-forming cyanobacteria. HetR regulates target genes via HetR-recognition sites. However, some genes (such as patS/patX) upregulated at the early stage of heterocyst differentiation possess DIF1 (or DIF⁺) motif (TCCGGA) promoters rather than HetR-recognition sites; hetR possesses both predicted regulatory elements. How HetR controls heterocyst-specific expression from DIF1 motif promoters remains to be answered. This study presents evidence that the expression from DIF1 motif promoters of hetR, patS and patX is more directly dependent on hetZ, a gene regulated by HetR via a HetR-recognition site. The HetR-binding site upstream of hetR is not required for the autoregulation of hetR. PatU3 (3′ portion of PatU) that interacts with HetZ may modulate the expression of hetR, hetZ and patS. These findings contribute to understanding of the mutual regulation of hetR, hetZ-patU and patS/patX in a large group of multicellular cyanobacteria.

Stress-adaptive gene discovery by exploiting collective decision-making of decentralized plant response systems

Despite having a network of cytoplasmic interconnections (plasmodesmata) facilitating rapid exchange of metabolites and signal molecules, plant cells use the extracellular matrix as an alternative route for cell-cell communication. The need for extracellular signalling in plasmodesmata-networked tissues is baffling. A hypothesis is proposed that this phenomenon defines the plant extracellular matrix as a 'democratic space' for collective decision-making in a decentralized system, similar to quorum-sensing in bacteria. Extracellular communication enables signal integration and coordination across several cell layers through ligand-activated plasma membrane receptors. Recent results from drought stress-adaptive responses and light-mediated signalling in cell death activation show operational utility of this decision-making process. Opportunities are discussed for new innovations in drought gene discovery using platforms targeting the extracellular matrix.

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.

Identification of a gene regulated by HetR, a master regulator of heterocyst differentiation, in the non-heterocyst-forming filamentous cyanobacterium Arthrospira platensis NIES-39

Cyanobacteria are a morphologically and physiologically diverse group of bacteria, which contains unicellular and multicellular filamentous strains. Some filamentous cyanobacteria, such as Anabaena sp. strain PCC 7120, form a differentiated cell called a heterocyst. The heterocyst is a specialized cell for nitrogen fixation and is differentiated from a vegetative cell in response to depletion of combined nitrogen in the medium. In Anabaena PCC 7120, it has been demonstrated that hetR, which encodes a transcriptional regulator, is necessary and sufficient for heterocyst differentiation. However, comprehensive genomic analysis of cyanobacteria has shown that hetR is present in non-heterocyst-forming cyanobacteria. Almost all filamentous cyanobacteria have hetR, but unicellular cyanobacteria do not. In this study, we conducted genetic and biochemical analyses of hetR (NIES39_C03480) of the non-heterocyst-forming cyanobacterium Arthrospira platensis NIES-39. HetR of A. platensis was able to complement the hetR mutation in Anabena PCC 7120 and recognized the same DNA sequence as Anabaena HetR. A search of the A. platensis genome revealed the HetR-recognition sequence within the promoter region of NIES39_O04230, which encodes a protein of unknown function. Expression from the NIES39_O04230 promoter could be suppressed by HetR in Anabaena PCC 7120. These data support the conclusion that NIES39_O04230 is regulated by HetR in A. platensis NIES-39.

patD, a Gene Regulated by NtcA, Is Involved in the Optimization of Heterocyst Frequency in the Cyanobacterium Anabaena sp. Strain PCC 7120

In the filamentous multicellular cyanobacterium Anabaena sp. strain PCC 7120, 5 to 10% of the cells differentiate into heterocysts, which are specialized in N₂ fixation. Heterocysts and vegetative cells are mutually dependent for filament growth through nutrient exchange. Thus, the heterocyst frequency should be optimized to maintain the cellular carbon and nitrogen (C/N) balance for filament fitness in the environment. Here, we report the identification of patD, whose expression is directly activated in developing cells by the transcription factor NtcA. The inactivation of patD increases heterocyst frequency and promotes the upregulation of the positive regulator of heterocyst development hetR, whereas its overexpression decreases the heterocyst frequency. The change in heterocyst frequency resulting from the inactivation of patD leads to the reduction in competitiveness of the filaments under combined-nitrogen-depleted conditions. These results indicate that patD regulates heterocyst frequency in Anabaena sp. PCC 7120, ensuring its optimal filament growth.IMPORTANCE Microorganisms have evolved various strategies in order to adapt to the environment and compete with other organisms. Heterocyst differentiation is a prokaryotic model for studying complex cellular regulation. The NtcA-regulated gene patD controls the ratio of heterocysts relative to vegetative cells on the filaments of Anabaena sp. strain PCC 7120. Such a regulation provides a mechanism through which carbon fixation by vegetative cells and nitrogen fixation by heterocysts are properly balanced to ensure optimal growth and keep a competitive edge for long-term survival.

Preventing Accidental Heterocyst Development in Cyanobacteria

The filamentous cyanobacterium Anabaena can form heterocysts specialized in N₂ fixation, mostly through a cascade of transcriptional activation in response to the nitrogen starvation signal 2-oxoglutarate. It is reported now that a transcription repressor, CalA, acts as a safety device to prevent heterocyst development under certain conditions where the 2-oxoglutarate level may touch the threshold to trigger unnecessary initiation of heterocyst development. Such a control may increase the fitness of Anabaena under a constantly changing environment.

Manipulation of Pattern of Cell Differentiation in a hetR Mutant of Anabaena sp. PCC 7120 by Overexpressing hetZ Alone or with hetP

In the filamentous cyanobacterium, Anabaena sp. PCC 7120, single heterocysts differentiate at semi-regular intervals in response to nitrogen stepdown. HetR is a principal regulator of heterocyst differentiation, and hetP and hetZ are two genes that are regulated directly by HetR. In a hetR mutant generated from the IHB (Institute of Hydrobiology) substrain of PCC 7120, heterocyst formation can be restored by moderate expression of hetZ and hetP. The resulting heterocysts are located at terminal positions. We used a tandem promoter, P_rbcLP_petE, to express hetZ and hetP strongly in the hetR mutant. Co-expression of hetZ and hetP enabled the hetR mutant to form multiple contiguous heterocysts at both terminal and intercalary positions. Expression of hetZ, alone resulted in terminally located heterocysts, whereas expression of hetP, alone produced enlarged cells in strings. In the absence of HetR, formation of heterocysts was insensitive to the peptide inhibitor, RGSGR.

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.

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.

Self-Assembled N-Heterocyclic Carbene-Based Carboxymethylated Dextran Monolayers on Gold as a Tunable Platform for Designing Affinity-Capture Biosensor Surfaces

Sensor surfaces play a predominant role in the development of optical biosensor technologies for the analysis of biomolecular interactions. Thiol-based self-assembled monolayers (SAMs) on gold have been widely used as linker layers for sensor surfaces. However, the degradation of the thiol-gold bond can limit the performance and durability of such surfaces, directly impacting their performance and cost-effectiveness. To this end, a new family of materials based on N-heterocyclic carbenes (NHCs) has emerged as an alternative for surface modification, capable of self-assembling onto a gold surface with higher affinity and superior stability as compared to the thiol-based systems. Here we demonstrate three applications of NHC SAMs supporting a dextran layer as a tunable platform for developing various affinity-capture biosensor surfaces. We describe the development and testing of NHC-based dextran biosensor surfaces modified with each of streptavidin, nitrilotriacetic acid, and recombinant Protein A. These affinity-capture sensor surfaces enable oriented binding of ligands for optimal performance in biomolecular assays. Together, the intrinsic high stability and flexible design of the NHC biosensing platforms show great promise and open up exciting possibilities for future biosensing applications.

Role of PatS and cell type on the heterocyst spacing pattern in a filamentous branching cyanobacterium

Cell differentiation is one of the marks of multicellular organisms. Terminally specialised nitrogen-fixing cells, termed heterocysts, evolved in filamentous cyanobacteria more than 2 Gya. The development of their spacing pattern has been thoroughly investigated in model organisms such as Anabaena sp. PCC 7120. This paper focuses on the more complex, branching cyanobacterium Mastigocladus laminosus (Stigonematales). Contrary to what has been previously published, a heterocyst spacing pattern is present in M. laminosus but it varies with the age of the culture and the morphology of the cells. Heterocysts in young, narrow trichomes were more widely spaced (∼14.8 cells) than those in old, wide trichomes (∼9.4 cells). Biochemical and transgenic experiments reveal that the heterocyst spacing pattern is affected by the heterocyst inhibitor PatS. Addition of the pentapeptide RGSGR (PatS-5) to the growth medium and overexpression of patS from Anabaena sp. PCC 7120 in M. laminosus resulted in the loss of heterocyst differentiation under nitrogen deprivation. Bioinformatics investigations indicated that putative PatS sequences within cyanobacteria are highly diverse, and fall into two main clades. Both are present in most branching cyanobacteria. Despite its more complex, branching phenotype, M. laminosus appears to use a PatS-based pathway for heterocyst differentiation, a property shared by Anabaena/Nostoc.

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.

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.