Lineage-based identification of cellular states and expression programs

We present a method, LineageProgram, that uses the developmental lineage relationship of observed gene expression measurements to improve the learning of developmentally relevant cellular states and expression programs. We find that incorporating lineage information allows us to significantly improve both the predictive power and interpretability of expression programs that are derived from expression measurements from in vitro differentiation experiments. The lineage tree of a differentiation experiment is a tree graph whose nodes describe all of the unique expression states in the input expression measurements, and edges describe the experimental perturbations applied to cells. Our method, LineageProgram, is based on a log-linear model with parameters that reflect changes along the lineage tree. Regularization with L1 that based methods controls the parameters in three distinct ways: the number of genes change between two cellular states, the number of unique cellular states, and the number of underlying factors responsible for changes in cell state. The model is estimated with proximal operators to quickly discover a small number of key cell states and gene sets. Comparisons with existing factorization, techniques, such as singular value decomposition and non-negative matrix factorization show that our method provides higher predictive power in held, out tests while inducing sparse and biologically relevant gene sets.National Institutes of Health (U.S.) (P01-NS055923)National Institutes of Health (U.S.) (1-UL1-RR024920

Exponential Distribution of Locomotion Activity in Cell Cultures

In vitro velocities of several cell types have been measured using computer controlled video microscopy, which allowed to record the cells' trajectories over several days. On the basis of our large data sets we show that the locomotion activity displays a universal exponential distribution. Thus, motion resulting from complex cellular processes can be well described by an unexpected, but very simple distribution function. A simple phenomenological model based on the interaction of various cellular processes and finite ATP production rate is proposed to explain these experimental results.Comment: 4 pages, 3 figure

Saltatory remodeling of Hox chromatin in response to rostrocaudal patterning signals

Hox genes controlling motor neuron subtype identity are expressed in rostrocaudal patterns that are spatially and temporally collinear with their chromosomal organization. Here we demonstrate that Hox chromatin is subdivided into discrete domains that are controlled by rostrocaudal patterning signals that trigger rapid, domain-wide clearance of repressive histone H3 Lys27 trimethylation (H3K27me3) polycomb modifications. Treatment of differentiating mouse neural progenitors with retinoic acid leads to activation and binding of retinoic acid receptors (RARs) to the Hox1–Hox5 chromatin domains, which is followed by a rapid domain-wide removal of H3K27me3 and acquisition of cervical spinal identity. Wnt and fibroblast growth factor (FGF) signals induce expression of the Cdx2 transcription factor that binds and clears H3K27me3 from the Hox1–Hox9 chromatin domains, leading to specification of brachial or thoracic spinal identity. We propose that rapid clearance of repressive modifications in response to transient patterning signals encodes global rostrocaudal neural identity and that maintenance of these chromatin domains ensures the transmission of positional identity to postmitotic motor neurons later in development.Leona M. and Harry B. Helmsley Charitable TrustNational Institutes of Health (U.S.) (Grant P01 NS055923)Smith Family Foundatio

Microfluidic Perfusion for Regulating Diffusible Signaling in Stem Cells

Background Autocrine & paracrine signaling are widespread both in vivo and in vitro, and are particularly important in embryonic stem cell (ESC) pluripotency and lineage commitment. Although autocrine signaling via fibroblast growth factor-4 (FGF4) is known to be required in mouse ESC (mESC) neuroectodermal specification, the question of whether FGF4 autocrine signaling is sufficient, or whether other soluble ligands are also involved in fate specification, is unknown. The spatially confined and closed-loop nature of diffusible signaling makes its experimental control challenging; current experimental approaches typically require prior knowledge of the factor/receptor in order to modulate the loop. A new approach explored in this work is to leverage transport phenomena at cellular resolution to downregulate overall diffusible signaling through the physical removal of cell-secreted ligands. Methodology/Principal Findings We develop a multiplex microfluidic platform to continuously remove cell-secreted (autocrine\paracrine) factors to downregulate diffusible signaling. By comparing cell growth and differentiation in side-by-side chambers with or without added cell-secreted factors, we isolate the effects of diffusible signaling from artifacts such as shear, nutrient depletion, and microsystem effects, and find that cell-secreted growth factor(s) are required during neuroectodermal specification. Then we induce FGF4 signaling in minimal chemically defined medium (N2B27) and inhibit FGF signaling in fully supplemented differentiation medium with cell-secreted factors to determine that the non-FGF cell-secreted factors are required to promote growth of differentiating mESCs. Conclusions/Significance Our results demonstrate for the first time that flow can downregulate autocrine\paracrine signaling and examine sufficiency of extracellular factors. We show that autocrine\paracrine signaling drives neuroectodermal commitment of mESCs through both FGF4-dependent and -independent pathways. Overall, by uncovering autocrine\paracrine processes previously hidden in conventional culture systems, our results establish microfluidic perfusion as a technique to study and manipulate diffusible signaling in cell systems.National Institutes of Health (U.S.) (NIH grant No. EB007278)Swiss National Science FoundationSwiss National Science Foundatio

Neurotrophic requirements of human motor neurons defined using amplified and purified stem-cell derived cultures

Neurotrophic requirements of human motor neurons defined using amplified and purified stem-cell derived culturesHuman motor neurons derived from embryonic and induced pluripotent stem cells (hESCs and hiPSCs) are a potentially important tool for studying motor neuron survival and pathological cell death. However, their basic survival requirements remain poorly characterized. Here, we sought to optimize a robust survival assay and characterize their response to different neurotrophic factors. First, to increase motor neuron yield, we screened a small-molecule collection and found that the Rho-associated kinase (ROCK) inhibitor Y-27632 enhances motor neuron progenitor proliferation up to 4-fold in hESC and hiPSC cultures. Next, we FACS-purified motor neurons expressing the Hb9::GFP reporter from Y-27632-amplified embryoid bodies and cultured them in the presence of mitotic inhibitors to eliminate dividing progenitors. Survival of these purified motor neurons in the absence of any other cell type was strongly dependent on neurotrophic support. GDNF, BDNF and CNTF all showed potent survival effects (EC(50) 1-2 pM). The number of surviving motor neurons was further enhanced in the presence of forskolin and IBMX, agents that increase endogenous cAMP levels. As a demonstration of the ability of the assay to detect novel neurotrophic agents, Y-27632 itself was found to support human motor neuron survival. Thus, purified human stem cell-derived motor neurons show survival requirements similar to those of primary rodent motor neurons and can be used for rigorous cell-based screening.This work was funded by Project A.L.S., P2ALS and NYSTEM grant number CO24415. The work of N.J.L. was supported by the Portuguese Foundation for Science and Technology SFRH/BD/33421/2008 and the Luso-American Development Foundation. B.J.-K. was supported by the National Institute of Neurological Disorders and Stroke (NINDS). L.R. was supported by the Swedish Brain Foundation/Hjarnfonden. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

An Integrated Model of Multiple-Condition ChIP-Seq Data Reveals Predeterminants of Cdx2 Binding

Regulatory proteins can bind to different sets of genomic targets in various cell types or conditions. To reliably characterize such condition-specific regulatory binding we introduce MultiGPS, an integrated machine learning approach for the analysis of multiple related ChIP-seq experiments. MultiGPS is based on a generalized Expectation Maximization framework that shares information across multiple experiments for binding event discovery. We demonstrate that our framework enables the simultaneous modeling of sparse condition-specific binding changes, sequence dependence, and replicate-specific noise sources. MultiGPS encourages consistency in reported binding event locations across multiple-condition ChIP-seq datasets and provides accurate estimation of ChIP enrichment levels at each event. MultiGPS's multi-experiment modeling approach thus provides a reliable platform for detecting differential binding enrichment across experimental conditions. We demonstrate the advantages of MultiGPS with an analysis of Cdx2 binding in three distinct developmental contexts. By accurately characterizing condition-specific Cdx2 binding, MultiGPS enables novel insight into the mechanistic basis of Cdx2 site selectivity. Specifically, the condition-specific Cdx2 sites characterized by MultiGPS are highly associated with pre-existing genomic context, suggesting that such sites are pre-determined by cell-specific regulatory architecture. However, MultiGPS-defined condition-independent sites are not predicted by pre-existing regulatory signals, suggesting that Cdx2 can bind to a subset of locations regardless of genomic environment. A summary of this paper appears in the proceedings of the RECOMB 2014 conference, April 2–5.National Science Foundation (U.S.) (Graduate Research Fellowship under Grant 0645960)National Institutes of Health (U.S.) (grant P01 NS055923)Pennsylvania State University. Center for Eukaryotic Gene Regulatio

Alteration of inhibitory circuits in the somatosensory cortex of Ts65Dn mice, a model for Down's syndrome

Down's syndrome (DS), with an incidence of one in 800 live births, is the most common genetic disorder associated with mental retardation. This trisomy on chromosome 21 induces a variable phenotype in which the only common feature is the presence of mental retardation. The neural mechanisms underlying mental retardation might include defects in the formation of neuronal networks and neural plasticity. DS patients have alterations in the morphology, the density and the distribution of dendritic spines in the pyramidal neurons of the cortex. Our hypothesis is that the deficits in dendritic arborization observed in the principal neurons of DS patients and Ts65Dn mice (a model for DS that mimics most of the structural alterations observed in humans) may be mediated to some extent by changes in their inhibitory inputs. Different types of interneurons control different types of inhibition. Therefore, to understand well the changes in inhibition in DS, it is necessary to study the different types of interneurons separately. We have studied the expression of synaptophysin, Glutamic acid decarboxylase-67 (GAD-67) and calcium-binding protein-expressing cells in the primary somatosensory cortex of 4¿5 month old Ts65Dn mice. We have observed an increment of GAD67 immunoreactivity that is related mainly to an increment of calretinin-immunoreactive cells and among them the ones with bipolar morphology. Since there is a propensity for epilepsy in DS patients, this increase in interneurons might reflect an attempt by the system to block overexcitation rather than an increment in total inhibition and could explain the deficit in interneurons and principal cells observed in elderly DS patients

Short RNAs Are Transcribed from Repressed Polycomb Target Genes and Interact with Polycomb Repressive Complex-2

Polycomb proteins maintain cell identity by repressing the expression of developmental regulators specific for other cell types. Polycomb repressive complex-2 (PRC2) catalyzes trimethylation of histone H3 lysine-27 (H3K27me3). Although repressed, PRC2 targets are generally associated with the transcriptional initiation marker H3K4me3, but the significance of this remains unclear. Here, we identify a class of short RNAs, ~50–200 nucleotides in length, transcribed from the 5′ end of polycomb target genes in primary T cells and embryonic stem cells. Short RNA transcription is associated with RNA polymerase II and H3K4me3, occurs in the absence of mRNA transcription, and is independent of polycomb activity. Short RNAs form stem-loop structures resembling PRC2 binding sites in Xist, interact with PRC2 through SUZ12, cause gene repression in cis, and are depleted from polycomb target genes activated during cell differentiation. We propose that short RNAs play a role in the association of PRC2 with its target genes.National Institutes of Health (U.S.) (Grant HG002668)National Institutes of Health (U.S.) (Grant NS055923

In vivo assembly of the axon initial segment in motor neurons

International audienceThe axon initial segment (AIS) is responsible for both the modulation of action potentials and the maintenance of neuronal polarity. Yet, the molecular mechanisms controlling its assembly are incompletely understood. Our study in single electroporated motor neurons in mouse embryos revealed that AnkyrinG (AnkG), the AIS master organizer, is undetectable in bipolar migrating motor neurons, but is already expressed at the beginning of axonogenesis at E9.5 and initially distributed homogeneously along the entire growing axon. Then, from E11.5, a stage when AnkG is already apposed to the membrane, as observed by electron microscopy, the protein progressively becomes restricted to the proximal axon. Analysis on the global motor neurons population indicated that Neurofascin follows an identical spatio-temporal distribution, whereas sodium channels and beta 4-spectrin only appear along AnkG(+) segments at E11.5. Early patch-clamp recordings of individual motor neurons indicated that at E12.5 these nascent AISs are already able to generate spikes. Using knock-out mice, we demonstrated that neither beta 4-spectrin nor Neurofascin control the distal-to-proximal restriction of AnkG

Discovery of a potent and selective CDKL5/GSK3 chemical probe that is neuroprotective

Despite mediating several essential processes in the brain, including during development, cyclin-dependent kinase-like 5 (CDKL5) remains a poorly characterized human protein kinase. Accordingly, its substrates, functions, and regulatory mechanisms have not been fully described. We realized that availability of a potent and selective small molecule probe targeting CDKL5 could enable illumination of its roles in normal development as well as in diseases where it has become aberrant due to mutation. We prepared analogs of AT-7519, a compound that has advanced to phase II clinical trials and is a known inhibitor of several cyclin-dependent kinases (CDKs) and cyclin-dependent kinase-like kinases (CDKLs). We identified analog 2 as a highly potent and cell-active chemical probe for CDKL5/GSK3 (glycogen synthase kinase 3). Evaluation of its kinome-wide selectivity confirmed that analog 2 demonstrates excellent selectivity and only retains GSK3α/β affinity. We next demonstrated the inhibition of downstream CDKL5 and GSK3α/β signaling and solved a co-crystal structure of analog 2 bound to human CDKL5. A structurally similar analog (4) proved to lack CDKL5 affinity and maintain potent and selective inhibition of GSK3α/β, making it a suitable negative control. Finally, we used our chemical probe pair (2 and 4) to demonstrate that inhibition of CDKL5 and/or GSK3α/β promotes the survival of human motor neurons exposed to endoplasmic reticulum stress. We have demonstrated a neuroprotective phenotype elicited by our chemical probe pair and exemplified the utility of our compounds to characterize the role of CDKL5/GSK3 in neurons and beyond