Quantifying 35 transcripts in a single tube: model-based calibration of the GeXP multiplex RT-PCR assay.

BACKGROUND: Quantitative analysis of differential gene expression is of central importance in molecular life sciences. The Gene eXpression Profiling technology (GeXP) relies on multiplex RT-PCR and subsequent capillary electrophoretic separation of the amplification products and allows to quantify the transcripts of at least 35 genes with a single reaction and one dye. RESULTS: We provide a kinetic model of primer binding and PCR product formation as the rational basis for taking and evaluating calibration curves. The calibration procedure and the model predictions were validated with the help of a purposefully designed data processing workflow supported by easy-to-use Perl scripts for calibration, data evaluation, and quality control. We further demonstrate the robustness and linearity of quantification of individual transcripts at variable relative abundance of other co-amplified transcripts in a complex mixture of RNAs isolated from differentiating Physarum polycephalum plasmodial cells. CONCLUSIONS: We conclude that GeXP analysis is a robust, sensitive, and useful method when the transcripts of tens to few hundred genes are to be precisely quantified in a high number of samples.

Disentangling a complex response in cell reprogramming and probing the Waddington landscape by automatic construction of Petri nets.

We analyzed the developmental switch to sporulation of a multinucleate Physarum polycephalum plasmodial cell, a complex response to phytochrome photoreceptor activation. Automatic construction of Petri nets representing finite state machines assembled from trajectories of differential gene expression in single cells revealed alternative, genotype-dependent interconnected developmental routes and identified reversible steps, metastable states, commitment points, and subsequent irreversible steps together with molecular signatures associated with cell fate decision and differentiation. Formation of cyclic transits identified by transition invariants in mutants that are locked in a proliferative state is remarkable considering the view that oncogenic alterations may cause the formation of cancer attractors. We conclude that the Petri net approach is useful to probe the Waddington landscape of cellular reprogramming, to disentangle developmental routes for the reconstruction of the gene regulatory network, and to understand how genetic alterations or physiological conditions reshape the landscape eventually creating new basins of attraction. Unraveling the complexity of pathogenesis, disease progression, drug response or the analysis of attractor landscapes in other complex systems of uncertain structure might be additional fields of application.

Regulatory Dynamics of Cell Differentiation Revealed by True Time Series From Multinucleate Single Cells.

Dynamics of cell fate decisions are commonly investigated by inferring temporal sequences of gene expression states by assembling snapshots of individual cells where each cell is measured once. Ordering cells according to minimal differences in expression patterns and assuming that differentiation occurs by a sequence of irreversible steps, yields unidirectional, eventually branching Markov chains with a single source node. In an alternative approach, we used multi-nucleate cells to follow gene expression taking true time series. Assembling state machines, each made from single-cell trajectories, gives a network of highly structured Markov chains of states with different source and sink nodes including cycles, revealing essential information on the dynamics of regulatory events. We argue that the obtained networks depict aspects of the Waddington landscape of cell differentiation and characterize them as reachability graphs that provide the basis for the reconstruction of the underlying gene regulatory network.

Transcriptome reprogramming during developmental switching in Physarum polycephalum involves extensive remodeling of intracellular signaling networks.

Activation of a phytochrome photoreceptor triggers a program of Physarum polycephalum plasmodial cell differentiation through which a mitotic multinucleate protoplasmic mass synchronously develops into haploid spores formed by meiosis and rearrangement of cellular components. We have performed a transcriptome-wide RNAseq study of cellular reprogramming and developmental switching. RNAseq analysis revealed extensive remodeling of intracellular signaling and regulation in switching the expression of sets of genes encoding transcription factors, kinases, phosphatases, signal transduction proteins, RNA-binding proteins, ubiquitin ligases, regulators of the mitotic and meiotic cell cycle etc. in conjunction with the regulation of genes encoding metabolic enzymes and cytoskeletal proteins. About 15% of the differentially expressed genes shared similarity with members of the evolutionary conserved set of core developmental genes of social amoebae. Differential expression of genes encoding regulators that act at the transcriptional, translational, and post-translational level indicates the establishment of a new state of cellular function and reveals evolutionary deeply conserved molecular changes involved in cellular reprogramming and differentiation in a prototypical eukaryote.

The Physarum polycephalum Genome Reveals Extensive Use of Prokaryotic Two-Component and Metazoan-Type Tyrosine Kinase Signaling.

Physarum polycephalum is a well-studied microbial eukaryote with unique experimental attributes relative to other experimental model organisms. It has a sophisticated life cycle with several distinct stages including amoebal, flagellated, and plasmodial cells. It is unusual in switching between open and closed mitosis according to specific life-cycle stages. Here we present the analysis of the genome of this enigmatic and important model organism and compare it with closely related species. The genome is littered with simple and complex repeats and the coding regions are frequently interrupted by introns with a mean size of 100 bases. Complemented with extensive transcriptome data, we define approximately 31,000 gene loci, providing unexpected insights into early eukaryote evolution. We describe extensive use of histidine kinase-based two-component systems and tyrosine kinase signaling, the presence of bacterial and plant type photoreceptors (phytochromes, cryptochrome, and phototropin) and of plant-type pentatricopeptide repeat proteins, as well as metabolic pathways, and a cell cycle control system typically found in more complex eukaryotes. Our analysis characterizes P. polycephalum as a prototypical eukaryote with features attributed to the last common ancestor of Amorphea, that is, the Amoebozoa and Opisthokonts. Specifically, the presence of tyrosine kinases in Acanthamoeba and Physarum as representatives of two distantly related subdivisions of Amoebozoa argues against the later emergence of tyrosine kinase signaling in the opisthokont lineage and also against the acquisition by horizontal gene transfer.

Gene expression kinetics in individual plasmodial cells reveal alternative programs of differential regulation during commitment and differentiation.

During its life cycle, the amoebozoon Physarum polycephalum forms multinucleate plasmodial cells that can grow to macroscopic size while maintaining a naturally synchronous population of nuclei. Sporulation-competent plasmodia were stimulated through photoactivation of the phytochrome photoreceptor and the expression of sporulation marker genes was analyzed quantitatively by repeatedly taking samples of the same plasmodial cell at successive time points after the stimulus pulse. Principal component analysis of the gene expression data revealed that plasmodial cells take different trajectories leading to cell fate decision and differentiation and suggested that averaging over individual cells is inappropriate. Queries for genes with pairwise correlated expression kinetics revealed qualitatively different patterns of co-regulation, indicating that alternative programs of differential regulation are operational in individual plasmodial cells. At the single cell level, the response to stimulation of a non-sporulating mutant was qualitatively different as compared to the wild type with respect to the differentially regulated genes and their patterns of co-regulation. The observation of individual differences during commitment and differentiation supports the concept of a Waddington-type quasipotential landscape for the regulatory control of cell differentiation. Comparison of wild type and sporulation mutant data further supports the idea that mutations may impact the topology of this landscape.

EXACT2: the semantics of biomedical protocols.

BACKGROUND: The reliability and reproducibility of experimental procedures is a cornerstone of scientific practice. There is a pressing technological need for the better representation of biomedical protocols to enable other agents (human or machine) to better reproduce results. A framework that ensures that all information required for the replication of experimental protocols is essential to achieve reproducibility. To construct EXACT2 we manually inspected hundreds of published and commercial biomedical protocols from several areas of biomedicine. After establishing a clear pattern for extracting the required information we utilized text-mining tools to translate the protocols into a machine amenable format. We have verified the utility of EXACT2 through the successful processing of previously 'unseen' (not used for the construction of EXACT2)protocols. METHODS: We have developed the ontology EXACT2 (EXperimental ACTions) that is designed to capture the full semantics of biomedical protocols required for their reproducibility. RESULTS: The paper reports on a fundamentally new version EXACT2 that supports the semantically-defined representation of biomedical protocols. The ability of EXACT2 to capture the semantics of biomedical procedures was verified through a text mining use case. In this EXACT2 is used as a reference model for text mining tools to identify terms pertinent to experimental actions, and their properties, in biomedical protocols expressed in natural language. An EXACT2-based framework for the translation of biomedical protocols to a machine amenable format is proposed. CONCLUSIONS: The EXACT2 ontology is sufficient to record, in a machine processable form, the essential information about biomedical protocols. EXACT2 defines explicit semantics of experimental actions, and can be used by various computer applications. It can serve as a reference model for for the translation of biomedical protocols in natural language into a semantically-defined format.

Switch-like reprogramming of gene expression after fusion of multinucleate plasmodial cells of two Physarum polycephalum sporulation mutants.

Nonlinear dynamic processes involving the differential regulation of transcription factors are considered to impact the reprogramming of stem cells, germ cells, and somatic cells. Here, we fused two multinucleate plasmodial cells of Physarum polycephalum mutants defective in different sporulation control genes while being in different physiological states. The resulting heterokaryons established one of two significantly different expression patterns of marker genes while the plasmodial halves that were fused to each other synchronized spontaneously. Spontaneous synchronization suggests that switch-like control mechanisms spread over and finally control the entire plasmodium as a result of cytoplasmic mixing. Regulatory molecules due to the large volume of the vigorously streaming cytoplasm will define concentrations in acting on the population of nuclei and in the global setting of switches. Mixing of a large cytoplasmic volume is expected to damp stochasticity when individual nuclei deliver certain RNAs at low copy number into the cytoplasm. We conclude that spontaneous synchronization, the damping of molecular noise in gene expression by the large cytoplasmic volume, and the option to take multiple macroscopic samples from the same plasmodium provide unique options for studying the dynamics of cellular reprogramming at the single cell level.

JAK/STAT signalling--an executable model assembled from molecule-centred modules demonstrating a module-oriented database concept for systems and synthetic biology.

Mathematical models of molecular networks regulating biological processes in cells or organisms are most frequently designed as sets of ordinary differential equations. Various modularisation methods have been applied to reduce the complexity of models, to analyse their structural properties, to separate biological processes, or to reuse model parts. Taking the JAK/STAT signalling pathway with the extensive combinatorial cross-talk of its components as a case study, we make a natural approach to modularisation by creating one module for each biomolecule. Each module consists of a Petri net and associated metadata and is organised in a database publically accessible through a web interface (). The Petri net describes the reaction mechanism of a given biomolecule and its functional interactions with other components including relevant conformational states. The database is designed to support the curation, documentation, version control, and update of individual modules, and to assist the user in automatically composing complex models from modules. Biomolecule centred modules, associated metadata, and database support together allow the automatic creation of models by considering differential gene expression in given cell types or under certain physiological conditions or states of disease. Modularity also facilitates exploring the consequences of alternative molecular mechanisms by comparative simulation of automatically created models even for users without mathematical skills. Models may be selectively executed as an ODE system, stochastic, or qualitative models or hybrid and exported in the SBML format. The fully automated generation of models of redesigned networks by metadata-guided modification of modules representing biomolecules with mutated function or specificity is proposed.

Physarum polycephalum mutants in the photocontrol of sporulation display altered patterns in the correlated expression of developmentally regulated genes.

Physarum polycephalum is a lower eukaryote belonging to the amoebozoa group of organisms that forms macroscopic, multinucleate plasmodial cells during its developmental cycle. Plasmodia can exit proliferative growth and differentiate by forming fruiting bodies containing mononucleate, haploid spores. This process, called sporulation, is controlled by starvation and visible light. To genetically dissect the regulatory control of the commitment to sporulation, we have isolated plasmodial mutants that are altered in the photocontrol of sporulation in a phenotypic screen of N-ethyl-N-nitrosourea (ENU) mutagenized cells. Several non-sporulating mutants were analyzed by measuring the light-induced change in the expression pattern of a set of 35 genes using GeXP multiplex reverse transcription-polymerase chain reaction with RNA isolated from individual plasmodial cells. Mutants showed altered patterns of differentially regulated genes in response to light stimulation. Some genes clearly displayed pairwise correlation in terms of their expression level as measured in individual plasmodial cells. The pattern of pairwise correlation differed in various mutants, suggesting that different upstream regulators were disabled in the different mutants. We propose that patterns of pairwise correlation in gene expression might be useful to infer the underlying gene regulatory network.

Reconstruction of extended Petri nets from time-series data by using logical control functions.

The aim of this work is to extend a previously presented algorithm (Durzinsky et al. 2008b in Computational methods in systems biology, LNCS, vol 5307. Springer, Heidelberg, pp 328­346; Marwan et al. 2008 in Math Methods Oper Res 67:117­132) for the reconstruction of standard place/transition Petri nets from time-series of experimental data sets. This previously reported method finds provably all networks capable to reproduce the experimental observations. In this paper we enhance this approach to generate extended Petri nets involving mechanisms formally corresponding to catalytic or inhibitory dependencies that mediate the involved reactions. The new algorithm delivers the set of all extended Petri nets being consistent with the time-series data used for reconstruction. It is illustrated using the phosphate regulatory network of enterobacteria as a case study.

A next-generation sequencing approach to study the transcriptomic changes during the differentiation of physarum at the single-cell level.

Physarum polycephalum is a unicellular eukaryote belonging to the amoebozoa group of organisms. The complex life cycle involves various cell types that differ in morphology, function, and biochemical composition. Sporulation, one step in the life cycle, is a stimulus-controlled differentiation response of macroscopic plasmodial cells that develop into fruiting bodies. Well-established Mendelian genetics and the occurrence of macroscopic cells with a naturally synchronous population of nuclei as source of homogeneous cell material for biochemical analyses make Physarum an attractive model organism for studying the regulatory control of cell differentiation. Here, we develop an approach using RNA-sequencing (RNA-seq), without needing to rely on a genome sequence as a reference, for studying the transcriptomic changes during stimulus-triggered commitment to sporulation in individual plasmodial cells. The approach is validated through the obtained expression patterns and annotations, and particularly the results from up- and downregulated genes, which correlate well with previous studies.

Futile attempts to differentiate provide molecular evidence for individual differences within a population of cells during cellular reprogramming.

The heterogeneity of cell populations and the influence of stochastic noise might be important issues for the molecular analysis of cellular reprogramming at the system level. Here, we show that in Physarum polycephalum, the expression patterns of marker genes correlate with the fate decision of individual multinucleate plasmodial cells that had been exposed to a differentiation-inducing photostimulus. For several hours after stimulation, the expression kinetics of PI-3-kinase, piwi, and pumilio orthologs and other marker genes were qualitatively similar in all stimulated cells but quantitatively different in those cells that subsequently maintained their proliferative potential and failed to differentiate accordingly. The results suggest that the population of nuclei in an individual plasmodium behaves synchronously in terms of gene regulation to an extent that the plasmodium provides a source for macroscopic amounts of homogeneous single-cell material for analysing the dynamic processes of cellular reprogramming. Based on the experimental findings, we predict that circuits with switch-like behaviour that control the cell fate decision of a multinucleate plasmodium operate through continuous changes in the concentration of cellular regulators because the nuclear population suspended in a large cytoplasmic volume damps stochastic noise.

Petri nets in Snoopy: a unifying framework for the graphical display, computational modelling, and simulation of bacterial regulatory networks.

Using the example of phosphate regulation in enteric bacteria, we demonstrate the particular suitability of stochastic Petri nets to model biochemical phenomena and their simulative exploration by various features of the software tool Snoopy.

Reconstruction of extended Petri nets from time series data and its application to signal transduction and to gene regulatory networks.

BACKGROUND: Network inference methods reconstruct mathematical models of molecular or genetic networks directly from experimental data sets. We have previously reported a mathematical method which is exclusively data-driven, does not involve any heuristic decisions within the reconstruction process, and deliveries all possible alternative minimal networks in terms of simple place/transition Petri nets that are consistent with a given discrete time series data set. RESULTS: We fundamentally extended the previously published algorithm to consider catalysis and inhibition of the reactions that occur in the underlying network. The results of the reconstruction algorithm are encoded in the form of an extended Petri net involving control arcs. This allows the consideration of processes involving mass flow and/or regulatory interactions. As a non-trivial test case, the phosphate regulatory network of enterobacteria was reconstructed using in silico-generated time-series data sets on wild-type and in silico mutants. CONCLUSIONS: The new exact algorithm reconstructs extended Petri nets from time series data sets by finding all alternative minimal networks that are consistent with the data. It suggested alternative molecular mechanisms for certain reactions in the network. The algorithm is useful to combine data from wild-type and mutant cells and may potentially integrate physiological, biochemical, pharmacological, and genetic data in the form of a single model.

A predictive computational model of the kinetic mechanism of stimulus-induced transducer methylation and feedback regulation through CheY in archaeal phototaxis and chemotaxis.

BACKGROUND: Photo- and chemotaxis of the archaeon Halobacterium salinarum is based on the control of flagellar motor switching through stimulus-specific methyl-accepting transducer proteins that relay the sensory input signal to a two-component system. Certain members of the transducer family function as receptor proteins by directly sensing specific chemical or physical stimuli. Others interact with specific receptor proteins like the phototaxis photoreceptors sensory rhodopsin I and II, or require specific binding proteins as for example some chemotaxis transducers. Receptor activation by light or a change in receptor occupancy by chemical stimuli results in reversible methylation of glutamate residues of the transducer proteins. Both, methylation and demethylation reactions are involved in sensory adaptation and are modulated by the response regulator CheY. RESULTS: By mathematical modeling we infer the kinetic mechanisms of stimulus-induced transducer methylation and adaptation. The model (deterministic and in the form of ordinary differential equations) correctly predicts experimentally observed transducer demethylation (as detected by released methanol) in response to attractant and repellent stimuli of wildtype cells, a cheY deletion mutant, and a mutant in which the stimulated transducer species is methylation-deficient. CONCLUSIONS: We provide a kinetic model for signal processing in photo- and chemotaxis in the archaeon H. salinarum suggesting an essential role of receptor cooperativity, antagonistic reversible methylation, and a CheY-dependent feedback on transducer demethylation.

Snoopy--a unifying Petri net framework to investigate biomolecular networks.

SUMMARY: To investigate biomolecular networks, Snoopy provides a unifying Petri net framework comprising a family of related Petri net classes. Models can be hierarchically structured, allowing for the mastering of larger networks. To move easily between the qualitative, stochastic and continuous modelling paradigms, models can be converted into each other. We get models sharing structure, but specialized by their kinetic information. The analysis and iterative reverse engineering of biomolecular networks is supported by the simultaneous use of several Petri net classes, while the graphical user interface adapts dynamically to the active one. Built-in animation and simulation are complemented by exports to various analysis tools. Snoopy facilitates the addition of new Petri net classes thanks to its generic design. AVAILABILITY: Our tool with Petri net samples is available free of charge for non-commercial use at http://www-dssz.informatik.tu-cottbus.de/snoopy.html; supported operating systems: Mac OS X, Windows and Linux (selected distributions).

Transcriptomic changes arising during light-induced sporulation in Physarum polycephalum.

BACKGROUND: Physarum polycephalum is a free-living amoebozoan protist displaying a complex life cycle, including alternation between single- and multinucleate stages through sporulation, a simple form of cell differentiation. Sporulation in Physarum can be experimentally induced by several external factors, and Physarum displays many biochemical features typical for metazoan cells, including metazoan-type signaling pathways, which makes this organism a model to study cell cycle, cell differentiation and cellular reprogramming. RESULTS: In order to identify the genes associated to the light-induced sporulation in Physarum, especially those related to signal transduction, we isolated RNA before and after photoinduction from sporulation- competent cells, and used these RNAs to synthesize cDNAs, which were then analyzed using the 454 sequencing technology. We obtained 16,669 cDNAs that were annotated at every computational level. 13,169 transcripts included hit count data, from which 2,772 displayed significant differential expression (upregulated: 1,623; downregulated: 1,149). Transcripts with valid annotations and significant differential expression were later integrated into putative networks using interaction information from orthologs. CONCLUSIONS: Gene ontology analysis suggested that most significantly downregulated genes are linked to DNA repair, cell division, inhibition of cell migration, and calcium release, while highly upregulated genes were involved in cell death, cell polarization, maintenance of integrity, and differentiation. In addition, cell death- associated transcripts were overrepresented between the upregulated transcripts. These changes are associated to a network of actin-binding proteins encoded by genes that are differentially regulated before and after light induction.

Quantitative analysis of signal transduction in motile and phototactic cells by computerized light stimulation and model based tracking.

To investigate the responses of Halobacterium salinarum to stimulation with light (phototaxis and photokinesis), we designed an experimental setup consisting of optical devices for automatic video image acquisition and computer-controlled light stimulation, and developed algorithms to analyze physiological responses of the cells. Cells are categorized as motile and nonmotile by a classification scheme based on the square displacement of cell positions. Computerized tracking based on a dynamic model of the stochastic cell movement and a Kalman filter-based algorithm allows smoothed estimates of the cell tracks and the detection of physiological responses to complex stimulus patterns. The setup and algorithms were calibrated which allows quantitative measurements and systematic analysis of cellular sensing and response. Overall, the setup is flexible, extensible, and consists mainly of commercially available products. This facilitates modifications of the setup and algorithms for physiological studies of the motility of cells or microorganisms.