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

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00689-4.

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.

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.

A laboratory exercise for detecting microbial contaminants

The development and evaluation of a 6-hours laboratory class, based on capillary electrophoresis (CE) and the detection of microbial contaminants, is described. It can be easily scaled up or down, to suit class sizes up to 188 and completed in a shorter time scale. CE uses narrow-bore fused-silica capillaries to separate a complex array of large and small molecules. A laboratory exercise has been devised to illustrate how CE-based genetic analysis system processes DNA fragment analysis to detect three microbial contaminants. The protocol is relatively inexpensive and uses standard molecular biology reagents and equipment. © 2018 by The International Union of Biochemistry and Molecular Biology, 46(3):279-284, 2018.

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.

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.

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.

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.