In Press
Close J, Yao Z, Levi B, Miller J, Bakken T, Menon V, Ting J, Wall A, Krostag A, Thomsen E, et al. Single-cell profiling of an in vitro model of human interneuron development reveals temporal dynamics of cell type production and maturation. Neuron. In Press.
Jang S, Furchtgott L, Choubey S, Zou L-N, Doyle A, Menon V, Loew E, Krostag A-R, Martinez RA, Madisen L, et al. Probabilistic model of gene networks controlling embryonic stem cell differentiation inferred from single-cell transcriptomics. eLife. In Press.
Furchtgott L, Melton S, Menon V, Ramanathan S. Discovering sparse transcription factor codes for cell states and state transitions during development. eLife. In Press.
Yao Z, Menon V, Mich JK, Ku S, Krostag A, Martinez RA, Furchtgott L, Mulholland H, Bort S, Fuqua MA, et al. Neural Lineages of Human Embryonic Stem Cell derived Cell Types revealed by Single Cell RNA Seq. Cell Stem Cell. In Press;20 (1) :120-134.
Thomsen ER, Mich JK, Yao Z, Hodge RD, Doyle AM, Jang S, Shehata SI, Nelson AM, Shapovalova NV, Levi BP, et al. Fixed single-cell transcriptomic characterization of human radial glial diversity. Nat Methods. 2016;13 (1) :87-93.Abstract
The diverse progenitors that give rise to the human neocortex have been difficult to characterize because progenitors, particularly radial glia (RG), are rare and are defined by a combination of intracellular markers, position and morphology. To circumvent these problems, we developed Fixed and Recovered Intact Single-cell RNA (FRISCR), a method for profiling the transcriptomes of individual fixed, stained and sorted cells. Using FRISCR, we profiled primary human RG that constitute only 1% of the midgestation cortex and classified them as ventricular zone-enriched RG (vRG) that express ANXA1 and CRYAB, and outer subventricular zone-localized RG (oRG) that express HOPX. Our study identified vRG and oRG markers and molecular profiles, an essential step for understanding human neocortical progenitor development. FRISCR allows targeted single-cell profiling of any tissues that lack live-cell markers.
Kocabas A, Shen C-H, Guo ZV, Ramanathan S. Controlling interneuron activity in Caenorhabditis elegans to evoke chemotactic behaviour. Nature. 2012;490 (7419) :273-7.Abstract
Animals locate and track chemoattractive gradients in the environment to find food. With its small nervous system, Caenorhabditis elegans is a good model system in which to understand how the dynamics of neural activity control this search behaviour. Extensive work on the nematode has identified the neurons that are necessary for the different locomotory behaviours underlying chemotaxis through the use of laser ablation, activity recording in immobilized animals and the study of mutants. However, we do not know the neural activity patterns in C. elegans that are sufficient to control its complex chemotactic behaviour. To understand how the activity in its interneurons coordinate different motor programs to lead the animal to food, here we used optogenetics and new optical tools to manipulate neural activity directly in freely moving animals to evoke chemotactic behaviour. By deducing the classes of activity patterns triggered during chemotaxis and exciting individual neurons with these patterns, we identified interneurons that control the essential locomotory programs for this behaviour. Notably, we discovered that controlling the dynamics of activity in just one interneuron pair (AIY) was sufficient to force the animal to locate, turn towards and track virtual light gradients. Two distinct activity patterns triggered in AIY as the animal moved through the gradient controlled reversals and gradual turns to drive chemotactic behaviour. Because AIY neurons are post-synaptic to most chemosensory and thermosensory neurons, it is probable that these activity patterns in AIY have an important role in controlling and coordinating different taxis behaviours of the animal.
Müller P, Rogers KW, Jordan BM, Lee JS, Robson D, Ramanathan S, Schier AF. Differential diffusivity of Nodal and Lefty underlies a reaction-diffusion patterning system. Science. 2012;336 (6082) :721-4.Abstract
Biological systems involving short-range activators and long-range inhibitors can generate complex patterns. Reaction-diffusion models postulate that differences in signaling range are caused by differential diffusivity of inhibitor and activator. Other models suggest that differential clearance underlies different signaling ranges. To test these models, we measured the biophysical properties of the Nodal/Lefty activator/inhibitor system during zebrafish embryogenesis. Analysis of Nodal and Lefty gradients revealed that Nodals have a shorter range than Lefty proteins. Pulse-labeling analysis indicated that Nodals and Leftys have similar clearance kinetics, whereas fluorescence recovery assays revealed that Leftys have a higher effective diffusion coefficient than Nodals. These results indicate that differential diffusivity is the major determinant of the differences in Nodal/Lefty range and provide biophysical support for reaction-diffusion models of activator/inhibitor-mediated patterning.
McClean MN, Hersen P, Ramanathan S. Measuring in vivo signaling kinetics in a mitogen-activated kinase pathway using dynamic input stimulation. Methods Mol Biol. 2011;734 :101-19.Abstract
Determining the in vivo kinetics of a signaling pathway is a challenging task. We can measure a property we termed pathway bandwidth to put in vivo bounds on the kinetics of the mitogen-activated protein kinase (MAPk) signaling cascade in Saccharomyces cerevisiae that responds to hyperosmotic stress [the High Osmolarity Glycerol (HOG) pathway]. Our method requires stimulating cells with square waves of oscillatory hyperosmotic input (1 M sorbitol) over a range of frequencies and measuring the activity of the HOG pathway in response to this oscillatory input. The input frequency at which the pathway's steady-state activity drops precipitously because the stimulus is changing too rapidly for the pathway to respond faithfully is defined as the pathway bandwidth. In this chapter, we provide details of the techniques required to measure pathway bandwidth in the HOG pathway. These methods are generally useful and can be applied to signaling pathways in S. cerevisiae and other organisms whenever a rapid reporter of pathway activity is available.
Thomson M, Liu SJ, Zou L-N, Smith Z, Meissner A, Ramanathan S. Pluripotency factors in embryonic stem cells regulate differentiation into germ layers. Cell. 2011;145 (6) :875-89.Abstract
Cell fate decisions are fundamental for development, but we do not know how transcriptional networks reorganize during the transition from a pluripotent to a differentiated cell state. Here, we asked how mouse embryonic stem cells (ESCs) leave the pluripotent state and choose between germ layer fates. By analyzing the dynamics of the transcriptional circuit that maintains pluripotency, we found that Oct4 and Sox2, proteins that maintain ESC identity, also orchestrate germ layer fate selection. Oct4 suppresses neural ectodermal differentiation and promotes mesendodermal differentiation; Sox2 inhibits mesendodermal differentiation and promotes neural ectodermal differentiation. Differentiation signals continuously and asymmetrically modulate Oct4 and Sox2 protein levels, altering their binding pattern in the genome, and leading to cell fate choice. The same factors that maintain pluripotency thus also integrate external signals and control lineage selection. Our study provides a framework for understanding how complex transcription factor networks control cell fate decisions in progenitor cells.
McClean MN, Hersen P, Ramanathan S. In vivo measurement of signaling cascade dynamics. Cell Cycle. 2009;8 (3) :373-6.Abstract
Genetic and biochemical studies yield information about the component proteins and interactions involved in a cellular signaling pathway. However this parts inventory often does not immediately reveal the in vivo signal processing capabilities and function of the pathway. Signaling pathways are complex systems with dynamic behavior and a systems level approach is needed to understand the physiological roles they play within the cell. We recently used such an approach to measure the signal processing behavior of the budding yeast HOG MAP kinase pathway in response to precisely varied temporal stimuli controlled with a microfluidic device. Despite being a well-studied pathway with well-known components the signaling dynamics and biochemical parameters of this pathway were not known. Our approach allowed us to characterize the pathway's in vivo signal processing and put bounds on all of the in vivo reaction rates. The experimental and theoretical techniques used in our study are general and can be applied to understanding other signaling pathways in a range of biological systems.
Mody A, Weiner J, Ramanathan S. Modularity of MAP kinases allows deformation of their signalling pathways. Nat Cell Biol. 2009;11 (4) :484-91.Abstract
Eukaryotic protein kinase pathways have both grown in number and changed their network architecture during evolution. We wondered if there are pivotal proteins in these pathways that have been repeatedly responsible for forming new connections through evolution, thus changing the topology of the network; and if so, whether the underlying properties of these proteins could be exploited to re-engineer and rewire these pathways. We addressed these questions in the context of the mitogen-activated protein kinase (MAPK) pathways. MAPK proteins were found to have repeatedly acquired new specificities and interaction partners during evolution, suggesting that these proteins are pivotal in the kinase network. Using the MAPKs Fus3 and Hog1 of the Saccharomyces cerevisiae mating and hyper-osmolar pathways, respectively, we show that these pivotal proteins can be re-designed to achieve a wide variety of changes in the input-output properties of the MAPK network. Through an analysis of our experimental results and of the sequence and structure of these proteins, we show that rewiring of the network is possible due to the underlying modular design of the MAPKs. We discuss the implications of our findings on the radiation of MAPKs through evolution and on how these proteins achieve their specificity.
Guo ZV, Hart AC, Ramanathan S. Optical interrogation of neural circuits in Caenorhabditis elegans. Nat Methods. 2009;6 (12) :891-6.Abstract
The nematode Caenorhabditis elegans has a compact nervous system with only 302 neurons. Whereas most of the synaptic connections between these neurons have been identified by electron microscopy serial reconstructions, functional connections have been inferred between only a few neurons through combinations of electrophysiology, cell ablation, in vivo calcium imaging and genetic analysis. To map functional connections between neurons, we combined in vivo optical stimulation with simultaneous calcium imaging. We analyzed the connections from the ASH sensory neurons and RIM interneurons to the command interneurons AVA and AVD. Stimulation of ASH or RIM neurons using channelrhodopsin-2 (ChR2) resulted in activation of AVA neurons, evoking an avoidance behavior. Our results demonstrate that we can excite specific neurons expressing ChR2 while simultaneously monitoring G-CaMP fluorescence in several other neurons, making it possible to rapidly decipher functional connections in C. elegans neural circuits.
Nachman I, Ramanathan S. HIV-1 positive feedback and lytic fate. Nat Genet. 2008;40 (4) :382-3.
Hersen P, McClean M, Ramanathan S. [Dynamic analysis of cell signalling: the example of the Saccharomyces cerevisiae osmotic response]. Med Sci (Paris). 2008;24 (12) :1020-2.
Hersen P, McClean MN, Mahadevan L, Ramanathan S. Signal processing by the HOG MAP kinase pathway. Proc Natl Acad Sci U S A. 2008;105 (20) :7165-70.Abstract
Signaling pathways relay information about changes in the external environment so that cells can respond appropriately. How much information a pathway can carry depends on its bandwidth. We designed a microfluidic device to reliably change the environment of single cells over a range of frequencies. Using this device, we measured the bandwidth of the Saccharomyces cerevisiae signaling pathway that responds to high osmolarity. This prototypical pathway, the HOG pathway, is shown to act as a low-pass filter, integrating the signal when it changes rapidly and following it faithfully when it changes more slowly. We study the dependence of the pathway's bandwidth on its architecture. We measure previously unknown bounds on all of the in vivo reaction rates acting in this pathway. We find that the two-component Ssk1 branch of this pathway is capable of fast signal integration, whereas the kinase Ste11 branch is not. Our experimental techniques can be applied to other signaling pathways, allowing the measurement of their in vivo kinetics and the quantification of their information capacity.
McClean MN, Mody A, Broach JR, Ramanathan S. Cross-talk and decision making in MAP kinase pathways. Nat Genet. 2007;39 (3) :409-14.Abstract
Cells must respond specifically to different environmental stimuli in order to survive. The signal transduction pathways involved in sensing these stimuli often share the same or homologous proteins. Despite potential cross-wiring, cells show specificity of response. We show, through modeling, that the physiological response of such pathways exposed to simultaneous and temporally ordered inputs can demonstrate system-level mechanisms by which pathways achieve specificity. We apply these results to the hyperosmolar and pheromone mitogen-activated protein (MAP) kinase pathways in the yeast Saccharomyces cerevisiae. These two pathways specifically sense osmolar and pheromone signals, despite sharing a MAPKKK, Ste11, and having homologous MAPKs (Fus3 and Hog1). We show that in a single cell, the pathways are bistable over a range of inputs, and the cell responds to only one stimulus even when exposed to both. Our results imply that these pathways achieve specificity by filtering out spurious cross-talk through mutual inhibition. The variability between cells allows for heterogeneity of the decisions.
Nachman I, Regev A, Ramanathan S. Dissecting timing variability in yeast meiosis. Cell. 2007;131 (3) :544-56.Abstract
Cell-to-cell variability in the timing of cell-fate changes can be advantageous for a population of single-celled organisms growing in a fluctuating environment. We study timing variability during meiosis in Saccharomyces cerevisiae, initiated upon nutritional starvation. We use time-lapse fluorescence microscopy to measure the timing of meiotic events in single cells and find that the duration of meiosis is highly variable between cells. This variability is concentrated between the beginning of starvation and the onset of early meiosis genes. Cell-cycle variability and nutritional history have little effect on this timing variability. Rather, variation in the production rate of the meiotic master regulator Ime1 and its gradual increase over time govern this variability, and cell size effects are channeled through Ime1. These results tie phenotypic variability with expression dynamics of a transcriptional regulator and provide a general framework for the study of temporal developmental processes.
Hallatschek O, Hersen P, Ramanathan S, Nelson DR. Genetic drift at expanding frontiers promotes gene segregation. Proc Natl Acad Sci U S A. 2007;104 (50) :19926-30.Abstract
Competition between random genetic drift and natural selection play a central role in evolution: Whereas nonbeneficial mutations often prevail in small populations by chance, mutations that sweep through large populations typically confer a selective advantage. Here, however, we observe chance effects during range expansions that dramatically alter the gene pool even in large microbial populations. Initially well mixed populations of two fluorescently labeled strains of Escherichia coli develop well defined, sector-like regions with fractal boundaries in expanding colonies. The formation of these regions is driven by random fluctuations that originate in a thin band of pioneers at the expanding frontier. A comparison of bacterial and yeast colonies (Saccharomyces cerevisiae) suggests that this large-scale genetic sectoring is a generic phenomenon that may provide a detectable footprint of past range expansions.
Ramanathan S, Detwiler PB, Sengupta AM, Shraiman BI. G-protein-coupled enzyme cascades have intrinsic properties that improve signal localization and fidelity. Biophys J. 2005;88 (5) :3063-71.Abstract
G-protein-coupled enzyme cascades are used by eukaryotic cells to detect external signals and transduce them into intracellular messages that contain biological information relevant to the cell's function. Since G-protein-coupled receptors that are designed to detect different kinds of external signals can generate the same kind of intracellular response, effective signaling requires that there are mechanisms to increase signal specificity and fidelity. Here we examine the kinetic equations for the initial three stages in a generic G-protein-coupled cascade and show that the physical properties of the transduction pathway result in two intrinsic features that benefit signaling. 1), The response to a single activated receptor is naturally confined to a localized spatial domain, which could improve signal specificity by reducing cross talk. 2), The peak of the response generated by such a signaling domain is limited. This saturation effect reduces trial-to-trial variability and increases signaling fidelity by limiting the response to receptors that remain active for longer than average. We suggest that this mechanism for reducing response fluctuations may be a contributing factor in making the single photon responses of vertebrate retinal rods so remarkably reproducible.
Ramanathan S, Swain PS. Tracing the sources of cellular variation. Dev Cell. 2005;9 (5) :576-8.Abstract
As the adage says, variety is the spice of life, and despite our best attempts, cells, even those with the same genome, never seem to behave the same. By combining mathematical and experimental analyses, Colman-Lerner and colleagues propose, in a recent issue of Nature, a method to delicately unravel the sources of this variation. Applying their technique to the pheromone response in budding yeast, they show that much of the observed variation originates from cell cycle effects and is dependent on levels of pathway input.