Modeling heterocyst pattern formation in cyanobacteria

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.

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.

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.

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.

Intraspecific chemical communication in microalgae

Contents 516 I. 516 II. 518 III. 518 IV. 521 V. 523 VI. 523 VII. 526 526 References 526 SUMMARY: The relevance of infochemicals in the relationships between organisms is emerging as a fundamental aspect of aquatic ecology. Exchanges of chemical cues are likely to occur not only between organisms of different species, but also between conspecific individuals. Especially intriguing is the investigation of chemical communication in microalgae, because of the relevance of these organisms for global primary production and their key role in trophic webs. Intraspecific communication between algae has been investigated mostly in relation to sexuality and mating. The literature also contains information on other types of intraspecific chemical communication that have not always been explicitly tagged as ways to communicate to conspecifics. However, the proposed role of certain compounds as intraspecific infochemicals appears questionable. In this article, we make use of this plethora of information to describe the various instances of intraspecific chemical communication between conspecific microalgae and to identify the common traits and ecological significance of intraspecific communication. We also discuss the evolutionary implications of intraspecific chemical communication and the mechanisms by which it can be inherited. A special focus is the genetic diversity among conspecific algae, including the possibility that genetic diversity is an absolute requirement for intraspecific chemical communication.

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.

Differences in cell division rates drive the evolution of terminal differentiation in microbes

Multicellular differentiated organisms are composed of cells that begin by developing from a single pluripotent germ cell. In many organisms, a proportion of cells differentiate into specialized somatic cells. Whether these cells lose their pluripotency or are able to reverse their differentiated state has important consequences. Reversibly differentiated cells can potentially regenerate parts of an organism and allow reproduction through fragmentation. In many organisms, however, somatic differentiation is terminal, thereby restricting the developmental paths to reproduction. The reason why terminal differentiation is a common developmental strategy remains unexplored. To understand the conditions that affect the evolution of terminal versus reversible differentiation, we developed a computational model inspired by differentiating cyanobacteria. We simulated the evolution of a population of two cell types –nitrogen fixing or photosynthetic– that exchange resources. The traits that control differentiation rates between cell types are allowed to evolve in the model. Although the topology of cell interactions and differentiation costs play a role in the evolution of terminal and reversible differentiation, the most important factor is the difference in division rates between cell types. Faster dividing cells always evolve to become the germ line. Our results explain why most multicellular differentiated cyanobacteria have terminally differentiated cells, while some have reversibly differentiated cells. We further observed that symbioses involving two cooperating lineages can evolve under conditions where aggregate size, connectivity, and differentiation costs are high. This may explain why plants engage in symbiotic interactions with diazotrophic bacteria.

The evolution of multicellularity is one of the most fascinating topics of evolutionary biology. Without multicellularity the incredible diversity of extant life would not be possible. In many multicellular organisms with specialized cells, some cell types become terminally differentiated (somatic cells) and lose the ability to reproduce new organisms while other cells maintain this ability (germline). Little is known about the conditions that favor the evolution of terminal differentiation in multicellular organisms. To understand this problem we have developed a computational model, inspired by multicellular cyanobacteria, in which the cells in an organism composed of two cell types (photosynthetic and nitrogen fixing) are allowed to evolve from germline to soma cells. We find three striking results. First, faster dividing cell types always evolve to become the germline. Second, the conditions under which we find different outcomes from the model are in good agreement with the different forms of development observed in multicellular cyanobacteria. Third, some conditions lead to a symbiotic state in which the two cell types separate into different lineages evolving independently of one another. Remarkably, cyanobacteria are also known to engage in symbiotic relationships with plants, producing fixed nitrogen for the plant in exchange for carbohydrates.

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.

The 20th anniversary of EMBnet: 20 years of bioinformatics for the Life Sciences community

The EMBnet Conference 2008, focusing on 'Leading Applications and Technologies in Bioinformatics', was organized by the European Molecular Biology network (EMBnet) to celebrate its 20th anniversary. Since its foundation in 1988, EMBnet has been working to promote collaborative development of bioinformatics services and tools to serve the European community of molecular biology laboratories. This conference was the first meeting organized by the network that was open to the international scientific community outside EMBnet. The conference covered a broad range of research topics in bioinformatics with a main focus on new achievements and trends in emerging technologies supporting genomics, transcriptomics and proteomics analyses such as high-throughput sequencing and data managing, text and data-mining, ontologies and Grid technologies. Papers selected for publication, in this supplement to BMC Bioinformatics, cover a broad range of the topics treated, providing also an overview of the main bioinformatics research fields that the EMBnet community is involved in.