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

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.

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 evolution of mechanisms to produce phenotypic heterogeneity in microorganisms

In bacteria and other microorganisms, the cells within a population often show extreme phenotypic variation. Different species use different mechanisms to determine how distinct phenotypes are allocated between individuals, including coordinated, random, and genetic determination. However, it is not clear if this diversity in mechanisms is adaptive-arising because different mechanisms are favoured in different environments-or is merely the result of non-adaptive artifacts of evolution. We use theoretical models to analyse the relative advantages of the two dominant mechanisms to divide labour between reproductives and helpers in microorganisms. We show that coordinated specialisation is more likely to evolve over random specialisation in well-mixed groups when: (i) social groups are small; (ii) helping is more "essential"; and (iii) there is a low metabolic cost to coordination. We find analogous results when we allow for spatial structure with a more detailed model of cellular filaments. More generally, this work shows how diversity in the mechanisms to produce phenotypic heterogeneity could have arisen as adaptations to different environments.

Biological stoichiometry and growth dynamics of a diazotrophic cyanobacteria in nitrogen sufficient and deficient conditions

The role of nitrogen (N) fixation in determining the frequency, magnitude, and extent of harmful algal blooms (HABs) has not been well studied. Dolichospermum is a common HAB species that is diazotrophic (capable of N fixation) and thus growth is often considered never to be limited by low combined N sources. However, N fixation is energetically expensive and its cost during bloom formation has not been quantified. Additionally, it is unknown how acclimation to differing nutrient ratios affects growth and cellular carbon (C):N stoichiometry. Here, we test the hypotheses that diazotrophic cyanobacteria are homeostatic for N because of their ability to fix atmospheric N₂ and that previous acclimation to low N environments will result in more fixed N and lower C:N stoichiometry. Briefly, cultures that varied in resource N:phosphorus (P) ranging from 0.01 to 100 (atom), were seeded with Dolichospermum which were previously acclimated to low and high N:P conditions and then sampled temporally for growth and C:N stoichiometry. We found that Dolichospermum was not homeostatic for N and displayed classic signs of N limitation and elevated C:N stoichiometry, highlighting the necessary growth trade-off within cells when expending energy to fix N. Acclimation to N limited conditions caused differences in both C:N and fixed N at various time points in the experiment. These results highlight the importance of environmentally available N to a diazotrophic bloom, as well as how previous growth conditions can influence population growth during blooms experiencing variable N:P.

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.

Reconstruction of a regulated two-cell metabolic model to study biohydrogen production in a diazotrophic cyanobacterium Anabaena variabilis ATCC 29413

Anabaena variabilis is a diazotrophic filamentous cyanobacterium that differentiates to heterocysts and produces hydrogen as a byproduct. Study on metabolic interactions of the two differentiated cells provides a better understanding of its metabolism especially for improving hydrogen production. To this end, a genome-scale metabolic model for Anabaena variabilis ATCC 29413, iAM957, was reconstructed and evaluated in this research. Then, the model and transcriptomic data of the vegetative and heterocyst cells were applied to construct a regulated two-cell metabolic model. The regulated model improved prediction for biomass in high radiation levels. The regulated model predicts that heterocysts provide an oxygen-free environment and then, this model was used to find strategies for improving hydrogen production in heterocysts. The predictions indicate that the removal of uptake hydrogenase improves hydrogen production which is consistent with previous empirical research. Furthermore, the regulated model proposed activation of some reactions to provide redox cofactors which are required for improving hydrogen production up to 60% by bidirectional hydrogenase.

Spectral microscopic imaging of heterocysts and vegetative cells in two filamentous cyanobacteria based on spontaneous Raman scattering and photoluminescence by 976 nm excitation

Photosynthetic pigment-protein complexes are highly concentrated in thylakoid membranes of chloroplasts and cyanobacteria that emit strong autofluorescence (mainly 600-800 nm). In Raman scattering microscopy that enables imaging of pigment concentrations of thylakoid membranes, near infrared laser excitation at 1064 nm or visible laser excitation at 488-532 nm has been often employed in order to avoid the autofluorescence. Here we explored a new approach to Raman imaging of thylakoid membranes by using excitation wavelength of 976 nm. Two types of differentiated cells, heterocysts and vegetative cells, in two diazotrophic filamentous cyanobacteria, Anabaena variabilis, and Rivularia M-261, were characterized. Relative Raman scattering intensities of phycobilisomes of the heterocyst in comparison with the nearest vegetative cells of Rivularia remained at a significantly higher level than those of A. variabilis. It was also found that the 976 nm excitation induces photoluminescence around 1017-1175 nm from the two cyanobacteria, green alga (Parachlorella kessleri) and plant (Arabidopsis thaliana). We propose that this photoluminescence can be used as an index of concentration of chlorophyll a that has relatively small Raman scattering cross-sections. The Rivularia heterocysts that we analyzed were clearly classified into at least two subgroups based on the Chla-associated photoluminescence and carotenoid Raman bands, indicating two physiologically distinct states in the development or aging of the terminal heterocyst.

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.