Search

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. Here we examine the success rate of methods and data for differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBseq, AME, DREME, HOMER, KMAC, diffTF and DeepAccess) to discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that use gene expression, biological networks and chromatin accessibility data, and comprehensively test parameter and preprocessing of input data to optimize performance. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top ten candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF have higher correlation with the ranked significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and diffTF are optimal methods for transcription factor recovery that will allow for systematic prioritization of transcription factor candidates to aid in the design of new reprogramming protocols.

A key mechanism in cellular regulation is the ability of the transcriptional machinery to physically access DNA. Pioneer transcription factors interact with DNA to open chromatin, which subsequently enables changes to gene expression during development, disease, or as a response to environmental stimuli. However, the regulation of DNA accessibility via the recruitment of transcription factors is difficult to understand in the context of the native genome because every genomic site is distinct in multiple ways. Here we introduce the Multiplexed Integrated Accessibility Assay (MIAA), a multiplexed parallel reporter assay which measures changes to genome accessibility as a result of the integration of synthetic oligonucleotide phrase libraries into a controlled, natively inaccessible genomic context. We apply MIAA to measure the effects of sequence motifs on cell type-specific DNA accessibility between mouse embryonic stem cells and embryonic stem cell-derived definitive endoderm cells, screening a total of 7,905 distinct phrases. MIAA is able to recapitulate differential accessibility patterns of 100-nt sequences derived from natively differential genomic regions, identifying the presence of E-box motifs common to epithelial-mesenchymal transition driver transcription factors in stem cell-specific accessible regions that become repressed during differentiation to endoderm. We further present causal evidence that a single binding motif for a key regulatory transcription factor is sufficient to open chromatin, and classify sets of stem cell-specific, endoderm-specific, and shared pioneer factor motifs. We also demonstrate that over-expression of two definitive endoderm transcription factors, Brachyury and FoxA2, results in changes to accessibility in phrases containing their respective DNA-binding motifs. Finally, we use MIAA results to explore the order of motif interactions and identify preferential motif ordering arrangements that appear to have an effect on accessibility.

Genomic interactions provide important context to our understanding of the state of the genome. One question is whether specific transcription factor interactions give rise to genome organization. We introduce spatzie, an R package and a website that implements statistical tests for significant transcription factor motif cooperativity between enhancer-promoter interactions. We conducted controlled experiments under realistic simulated data from ChIP-seq to confirm spatzie is capable of discovering co-enriched motif interactions even in noisy conditions. We then use spatzie to investigate cell type specific transcription factor cooperativity within recent human ChIA-PET enhancer-promoter interaction data. The method is available online at https://spatzie.mit.edu.

Summary Proper assembly and function of the nervous system requires the generation of a uniquely diverse population of neurons expressing a cell-type-specific combination of effector genes that collectively define neuronal morphology, connectivity, and function. How countless partially overlapping but cell-type-specific patterns of gene expression are controlled at the genomic level remains poorly understood. Here we show that neuronal genes are associated with highly complex gene regulatory systems composed of independent cell-type- and cell-stage-specific regulatory elements that reside in expanded non-coding genomic domains. Mapping enhancer-promoter interactions revealed that motor neuron enhancers are broadly distributed across the large chromatin domains. This distributed regulatory architecture is not a unique property of motor neurons but is employed throughout the nervous system. The number of regulatory elements increased dramatically during the transition from invertebrates to vertebrates, suggesting that acquisition of new enhancers might be a fundamental process underlying the evolutionary increase in cellular complexity.

Calcium is a key signaling agent in animals and plants. Its involvement in the regulation of a wide range of processes has led to the question of how calcium signals can activate stimulus-specific responses. We introduce a computational framework for studying intracellular calcium signaling using elements of information theory. We use mutual information to quantify the differential activation of proteins in response to different calcium signals to provide an operational definition of specificity . Using optimization procedures this framework allows us to explore the biochemical determinants of calcium decoding. We explore simple toy models and general binding kinetics approaches to demonstrate the utility and limitations of the proposed framework. Unravelling signaling specificity is key for understanding information processing within cells and for the future design of synthetic nanodevices for molecular communications.

Motivation Sequence models based on deep neural networks have achieved state-of-the-art performance on regulatory genomics prediction tasks, such as chromatin accessibility and transcription factor binding. But despite their high accuracy, their contributions to a mechanistic understanding of the biology of regulatory elements is often hindered by the complexity of the predictive model and thus poor interpretability of its decision boundaries. To address this, we introduce seqgra, a deep learning pipeline that incorporates the rule-based simulation of biological sequence data and the training and evaluation of models, whose decision boundaries mirror the rules from the simulation process. Results We show that seqgra can be used to (i) generate data under the assumption of a hypothesized model of genome regulation, (ii) identify neural network architectures capable of recovering the rules of said model and (iii) analyze a model’s predictive performance as a function of training set size and the complexity of the rules behind the simulated data. Availability and implementation The source code of the seqgra package is hosted on GitHub (https://github.com/gifford-lab/seqgra). seqgra is a pip-installable Python package. Extensive documentation can be found at https://kkrismer.github.io/seqgra. Supplementary information Supplementary data are available at Bioinformatics online.

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. We examine the success rate of methods and data for directed differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBSeq, AME, DREME, HOMER, KMAC, diffTF, and DeepAccess) to correctly discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that utilize gene expression, biological networks, and chromatin accessibility data to identify eight sets of known reprogramming factors and comprehensively test parameter and pre-processing of input data to optimize performance of these methods. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top 10 candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF are more likely to consistently correctly rank the significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and DeepAccess are optimal methods for transcription factor recovery and ranking which will allow for systematic prioritization of transcription factor candidates to aid in the design of novel reprogramming protocols.

SummaryBuilding a nervous system is a protracted process that starts with the specification of individual neuron types and ends with the formation of mature neural circuits. The molecular mechanisms that regulate the temporal progression of maturation in individual cell types remain poorly understood. In this work, we have mapped the gene expression and chromatin accessibility changes in mouse spinal motor neurons throughout their lifetimes. We found that both motor neuron gene expression and putative regulatory elements are dynamic during the first three weeks of postnatal life, when motor circuits are maturing. Genes that are up-regulated during this time contribute to adult motor neuron diversity and function. Almost all of the chromatin regions that gain accessibility during maturation are motor neuron specific, yet a majority of the transcription factor binding motifs enriched in these regions are shared with other mature neurons. Collectively, these findings suggest that a core transcriptional program operates in a context-dependent manner to access cell-type-specific cis-regulatory systems associated with maturation genes. Discovery of general principles governing neuronal maturation might inform methods for transcriptional reprogramming of neuronal age and for improved modelling of age-related neurodegenerative diseases.

Planaria are complex metazoans that repair damage to their bodies and cease remodeling when a correct anatomy has been achieved. This model system offers a unique opportunity to understand how large‐scale anatomical homeostasis emerges from the activities of individual cells. Much progress has been made on the molecular genetics of stem cell activity in planaria. However, recent data also indicate that the global pattern is regulated by physiological circuits composed of ionic and neurotransmitter signaling. Here, we overview the multi‐scale problem of understanding pattern regulation in planaria, with specific focus on bioelectric signaling via ion channels and gap junctions (electrical synapses), and computational efforts to extract explanatory models from functional and molecular data on regeneration. We present a perspective that interprets results in this fascinating field using concepts from dynamical systems theory and computational neuroscience. Serving as a tractable nexus between genetic, physiological, and computational approaches to pattern regulation, planarian pattern homeostasis harbors many deep insights for regenerative medicine, evolutionary biology, and engineering.

Automated methods for the reverse-engineering of complex regulatory networks are paving the way for the inference of mechanistic comprehensive models directly from experimental data. These novel methods can infer not only the relations and parameters of the known molecules defined in their input datasets, but also unknown components and pathways identified as necessary by the automated algorithms. Identifying the molecular nature of these unknown components is a crucial step for making testable predictions and experimentally validating the models, yet no specific and efficient tools exist to aid in this process. To this end, we present here MoCha (Molecular Characterization), a tool optimized for the search of unknown proteins and their pathways from a given set of known interacting proteins. MoCha uses the comprehensive dataset of protein-protein interactions provided by the STRING database, which currently includes more than a billion interactions from over 2,000 organisms. MoCha is highly optimized, performing typical searches within seconds. We demonstrate the use of MoCha with the characterization of unknown components from reverse-engineered models from the literature. MoCha is useful for working on network models by hand or as a downstream step of a model inference engine workflow and represents a valuable and efficient tool for the characterization of unknown pathways using known data from thousands of organisms. MoCha and its source code are freely available online under the GPLv3 license.

Artificial neural networks are both a well-established tool in machine learning and a mathematical model of distributed information processing. Developmental and regenerative biology is in desperate need of conceptual models to explain how some species retain memories despite drastic reorganization, remodeling, or regeneration of the brain. Here, we formalize a method of artificial neural network perturbation and quantitatively analyze memory persistence during different types of topology change. We introduce this system as a computational model of the complex information processing mechanisms that allow memories to persist during significant cellular and morphological turnover in the brain. We found that perturbations in artificial neural networks have a general negative effect on the preservation of memory, but that the removal of neurons with different firing patterns can effectively minimize this memory loss. The training algorithms employed and the difficulty of the pattern recognition problem tested are key factors determining the impact of perturbations. The results show that certain perturbations, such as neuron splitting and scaling, can achieve memory persistence by functional recovery of lost patterning information. The study of models integrating both growth and reduction, combined with distributed information processing is an essential first step for a computational theory of pattern formation, plasticity, and robustness.

The shape of an animal body plan is constructed from protein components encoded by the genome. However, bioelectric networks composed of many cell types have their own intrinsic dynamics, and can drive distinct morphological outcomes during embryogenesis and regeneration. Planarian flatworms are a popular system for exploring body plan patterning due to their regenerative capacity, but despite considerable molecular information regarding stem cell differentiation and basic axial patterning, very little is known about how distinct head shapes are produced. Here, we show that after decapitation in G. dorotocephala, a transient perturbation of physiological connectivity among cells (using the gap junction blocker octanol) can result in regenerated heads with quite different shapes, stochastically matching other known species of planaria (S. mediterranea, D. japonica, and P. felina). We use morphometric analysis to quantify the ability of physiological network perturbations to induce different species-specific head shapes from the same genome. Moreover, we present a computational agent-based model of cell and physical dynamics during regeneration that quantitatively reproduces the observed shape changes. Morphological alterations induced in a genomically wild-type G. dorotocephala during regeneration include not only the shape of the head but also the morphology of the brain, the characteristic distribution of adult stem cells (neoblasts), and the bioelectric gradients of resting potential within the anterior tissues. Interestingly, the shape change is not permanent, after regeneration is complete, intact animals remodel back to G. dorotocephala-appropriate head shape within several weeks in a secondary phase of remodeling following initial complete regeneration. We present a conceptual model to guide future work to delineate the molecular mechanisms by which bioelectric networks stochastically select among a small set of discrete head morphologies. Taken together, these data and analyses shed light on important physiological modifiers of morphological information in dictating species-specific shape, and reveal them to be a novel instructive input into head patterning in regenerating planaria.