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.

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.

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.

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.

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.

A first glimpse at the transcriptome of Physarum polycephalum.

BACKGROUND: Physarum polycephalum, an acellular plasmodial species belongs to the amoebozoa, a major branch in eukaryote evolution. Its complex life cycle and rich cell biology is reflected in more than 2500 publications on various aspects of its biochemistry, developmental biology, cytoskeleton, and cell motility. It now can be genetically manipulated, opening up the possibility of targeted functional analysis in this organism. METHODS: Here we describe a large fraction of the transcribed genes by sequencing a cDNA library from the plasmodial stage of the developmental cycle. RESULTS: In addition to the genes for the basic metabolism we found an unexpected large number of genes involved in sophisticated signaling networks and identified potential receptors for environmental signals such as light. In accordance with the various developmental options of the plasmodial cell we found that many P. polycephalum genes are alternatively spliced. Using 30 donor and 30 acceptor sites we determined the splicing signatures of this species. Comparisons to various other organisms including Dictyostelium, the closest relative, revealed that roughly half of the transcribed genes have no detectable counterpart, thus potentially defining species specific adaptations. On the other hand, we found highly conserved proteins, which are maintained in the metazoan lineage, but absent in D. discoideum or plants. These genes arose possibly in the last common ancestor of Amoebozoa and Metazoa but were lost in D. discoideum. CONCLUSION: This work provides an analysis of up to half of the protein coding genes of Physarum polycephalum. The definition of splice motifs together with the description of alternatively spliced genes will provide a valuable resource for the ongoing genome project.

Automatic reconstruction of molecular and genetic networks from discrete time series data.

We apply a mathematical algorithm which processes discrete time series data to generate a complete list of Petri net structures containing the minimal number of nodes required to reproduce the data set. The completeness of the list as guaranteed by a mathematical proof allows to define a minimal set of experiments required to discriminate between alternative network structures. This in principle allows to prove all possible minimal network structures by disproving all alternative candidate structures. The dynamic behaviour of the networks in terms of a switching rule for the transitions of the Petri net is part of the result. In addition to network reconstruction, the algorithm can be used to determine how many yet undetected components at least must be involved in a certain process. The algorithm also reveals all alternative structural modifications of a network that are required to generate a predefined behaviour.

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.

Theory of time-resolved somatic complementation and its use to explore the sporulation control network in Physarum polycephalum.

Mutants of Physarum polycephalum can be complemented by fusion of plasmodial cells followed by cytoplasmic mixing. Complementation between strains carrying different mutational defects in the sporulation control network may depend on the signaling state of the network components. We have previously suggested that time-resolved somatic complementation (TRSC) analysis with such mutants may be used to probe network architecture and dynamics. By computer simulation it is now shown how and under which conditions the regulatory hierarchy of genes can be determined experimentally. A kinetic model of the sporulation control network is developed, which is then used to demonstrate how the mechanisms of TRSC can be understood and simulated at the kinetic level. On the basis of theoretical considerations, experimental parameters that determine whether functional complementation of two mutations will occur are identified. It is also shown how gene dosage-effect relationships can be employed for network analysis. The theoretical framework provided may be used to systematically analyze network structure and dynamics through time-resolved somatic complementation studies. The conclusions drawn are of general relevance in that they do not depend on the validity of the model from which they were derived.

Estimating the number of plasmids taken up by a eukaryotic cell during transfection and evidence that antisense RNA abolishes gene expression in Physarum polycephalum.

We have estimated the statistical distribution of the number of plasmids taken up by individual Jurkat lymphoma cells during electroporation in the presence of two plasmids, one encoding for yellow (EYFP) the other for cyan (ECFP) fluorescent protein. The plasmid concentration at which most of the cells take up only one plasmid or several molecules was determined by statistical analysis. We found that cells behaved slightly heterogeneous in plasmid uptake and describe how the homogeneity of a cell population can be quantified by Poisson statistics in order to identify experimental conditions that yield homogeneously transfection-competent cell populations. The experimental procedure worked out with Jurkat cells was applied to assay the effectiveness of antisense RNA in knocking down gene expression in Physarum polycephalum. Double transfection of flagellates with vectors encoding EYFP and antisense-EYFP revealed for the first time that gene expression can be suppressed by co-expression of antisense RNA in Physarum. Quantitative analysis revealed that one copy of antisense expressing gene per EYFP gene was sufficient to completely suppress formation of the EYFP protein in Physarum.

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.

The sequence of regulatory events in the sporulation control network of Physarum polycephalum analysed by time-resolved somatic complementation of mutants.

The developmental decision for sporulation of Physarum polycephalum plasmodia is under sensory control by environmental factors like visible light or heat shock and endogenous signals like glucose starvation. Several hours after perceiving an inductive stimulus, plasmodia become committed to sporulation; thereby, they lose their unlimited replicative potential and execute a developmental program that involves differentiation into various cell types required to form a mature fruiting body. Plasmodia are multinuclear single cells which spontaneously fuse upon physical contact. Fusion of mutant plasmodia and cytoplasmic mixing allows complementation studies to be performed at the functional level. Mutant cells altered in their ability to sporulate in response to phytochrome activation by far-red light were cured by fusion with wild-type or other mutant plasmodia. Phytochrome activation in one plasmodium and subsequent fusion with a non-induced plasmodium revealed that complementation of the two mutations depended on (i) which of two genetically distinct plasmodial cells was stimulated; and (ii) on the delay time elapsed between stimulation and cytoplasmic mixing. Such experiments allow us to determine the kinetics and the causal sequence of the regulatory events tagged by mutation.

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.

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.

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.