Evaluating Gene Expression Dynamics Using Pairwise RNA FISH Data

Recently, a novel approach has been developed to study gene expression in single cells with high time resolution using RNA Fluorescent In Situ Hybridization (FISH). The technique allows individual mRNAs to be counted with high accuracy in wild-type cells, but requires cells to be fixed; thus, each cell provides only a “snapshot” of gene expression. Here we show how and when RNA FISH data on pairs of genes can be used to reconstruct real-time dynamics from a collection of such snapshots. Using maximum-likelihood parameter estimation on synthetically generated, noisy FISH data, we show that dynamical programs of gene expression, such as cycles (e.g., the cell cycle) or switches between discrete states, can be accurately reconstructed. In the limit that mRNAs are produced in short-lived bursts, binary thresholding of the FISH data provides a robust way of reconstructing dynamics. In this regime, prior knowledge of the type of dynamics – cycle versus switch – is generally required and additional constraints, e.g., from triplet FISH measurements, may also be needed to fully constrain all parameters. As a demonstration, we apply the thresholding method to RNA FISH data obtained from single, unsynchronized cells of Saccharomyces cerevisiae. Our results support the existence of metabolic cycles and provide an estimate of global gene-expression noise. The approach to FISH data presented here can be applied in general to reconstruct dynamics from snapshots of pairs of correlated quantities including, for example, protein concentrations obtained from immunofluorescence assays

DNA splicing: computing by observing

Motivated by several techniques for observing molecular processes in real-time we introduce a computing device that stresses the role of the observer in biological computations and that is based on the observed behavior of a splicing system. The basic idea is to introduce a marked DNA strand into a test tube with other DNA strands and restriction enzymes. Under the action of these enzymes the DNA starts to splice. An external observer monitors and registers the evolution of the marked DNA strand. The input marked DNA strand is then accepted if its observed evolution follows a certain expected pattern. We prove that using simple observers (finite automata), applied on finite splicing systems (finite set of rules and finite set of axioms), the class of recursively enumerable languages can be recognized. © Springer Science+Business Media B.V. 2007

An Information Theoretic, Microfluidic-Based Single Cell Analysis Permits Identification of Subpopulations among Putatively Homogeneous Stem Cells

An incomplete understanding of the nature of heterogeneity within stem cell populations remains a major impediment to the development of clinically effective cell-based therapies. Transcriptional events within a single cell are inherently stochastic and can produce tremendous variability, even among genetically identical cells. It remains unclear how mammalian cellular systems overcome this intrinsic noisiness of gene expression to produce consequential variations in function, and what impact this has on the biologic and clinical relevance of highly ‘purified’ cell subgroups. To address these questions, we have developed a novel method combining microfluidic-based single cell analysis and information theory to characterize and predict transcriptional programs across hundreds of individual cells. Using this technique, we demonstrate that multiple subpopulations exist within a well-studied and putatively homogeneous stem cell population, murine long-term hematopoietic stem cells (LT-HSCs). These subgroups are defined by nonrandom patterns that are distinguishable from noise and are consistent with known functional properties of these cells. We anticipate that this analytic framework can also be applied to other cell types to elucidate the relationship between transcriptional and phenotypic variation

Programmable in situ amplification for multiplexed imaging of mRNA expression

In situ hybridization methods enable the mapping of mRNA expression within intact biological samples. With current approaches, it is challenging to simultaneously map multiple target mRNAs within whole-mount vertebrate embryos, representing a significant limitation in attempting to study interacting regulatory elements in systems most relevant to human development and disease. Here, we report a multiplexed fluorescent in situ hybridization method based on orthogonal amplification with hybridization chain reactions (HCR). With this approach, RNA probes complementary to mRNA targets trigger chain reactions in which fluorophore-labeled RNA hairpins self-assemble into tethered fluorescent amplification polymers. The programmability and sequence specificity of these amplification cascades enable multiple HCR amplifiers to operate orthogonally at the same time in the same sample. Robust performance is achieved when imaging five target mRNAs simultaneously in fixed whole-mount and sectioned zebrafish embryos. HCR amplifiers exhibit deep sample penetration, high signal-to-background ratios and sharp signal localization

High Throughput Ratio Imaging to Profile Caspase Activity: Potential Application in Multiparameter High Content Apoptosis Analysis and Drug Screening

Recent advancement in the area of green fluorescent protein techniques coupled with microscopic imaging has significantly contributed in defining and dissecting subcellular changes of apoptosis with high spatio-temporal resolution. Although single cell based studies using EGFP and associated techniques have provided valuable information of initiation and hierarchical changes of apoptosis, they are yet to be exploited for multiparameter cell based real time analysis for possible drug screening or pathway defining in a high throughput manner. Here we have developed multiple cancer cell lines expressing FRET sensors for active caspases and adapted them for high throughput live cell ratio imaging, enabling high content image based multiparameter analysis. Sensitivity of the system to detect live cell caspase activation was substantiated by confocal acceptor bleaching as well as wide field FRET imaging. Multiple caspase-specific activities of DEVDase, IETDase and LEHDase were analysed simultaneously with other decisive events of cell death. Through simultaneous analysis of caspase activation by FRET ratio change coupled with detection of mitochondrial membrane potential loss or superoxide generation, we identified several antitumor agents that induced caspase activation with or without membrane potential loss or superoxide generation. Also, cells that escaped the initial drug-induced caspase activation could be easily followed up for defining long term fate. Employing such a revisit imaging strategy of the same area, we have tracked the caspase surviving fractions with multiple drugs and its subsequent response to retreatment, revealing drug-dependent diverging fate of surviving cells. This thereby indicates towards a complex control of drug induced tumor resistance. The technique described here has wider application in both screening of compound libraries as well as in defining apoptotic pathways by linking multiple signaling to identify non-classical apoptosis inducing agents, the greatest advantage being that the high content information obtained are from individual cells rather than being population based

Microgenomic Analysis in Skeletal Muscle: Expression Signatures of Individual Fast and Slow Myofibers

BACKGROUND: Skeletal muscle is a complex, versatile tissue composed of a variety of functionally diverse fiber types. Although the biochemical, structural and functional properties of myofibers have been the subject of intense investigation for the last decades, understanding molecular processes regulating fiber type diversity is still complicated by the heterogeneity of cell types present in the whole muscle organ. METHODOLOGY/PRINCIPAL FINDINGS: We have produced a first catalogue of genes expressed in mouse slow-oxidative (type 1) and fast-glycolytic (type 2B) fibers through transcriptome analysis at the single fiber level (microgenomics). Individual fibers were obtained from murine soleus and EDL muscles and initially classified by myosin heavy chain isoform content. Gene expression profiling on high density DNA oligonucleotide microarrays showed that both qualitative and quantitative improvements were achieved, compared to results with standard muscle homogenate. First, myofiber profiles were virtually free from non-muscle transcriptional activity. Second, thousands of muscle-specific genes were identified, leading to a better definition of gene signatures in the two fiber types as well as the detection of metabolic and signaling pathways that are differentially activated in specific fiber types. Several regulatory proteins showed preferential expression in slow myofibers. Discriminant analysis revealed novel genes that could be useful for fiber type functional classification. CONCLUSIONS/SIGNIFICANCE: As gene expression analyses at the single fiber level significantly increased the resolution power, this innovative approach would allow a better understanding of the adaptive transcriptomic transitions occurring in myofibers under physiological and pathological condition