Inputs drive cell phenotype variability

Computational modelling of cardiac control following myocardial infarction using an in silico patient cohort

Loss of cardiac physiological function following myocardial infarction (MI) is accompanied by neural adaptations in the baroreflex that are compensatory in the short term, but then become associated with long-term disease progression. One marker of these adaptations is decreased baroreflex sensitivity, a strong predictor of post-MI mortality. The relative contributions of cardiac remodelling and neural adaptation in the sensory, central brainstem and peripheral ganglionic loci to baroreflex sensitivity changes remain underexplored. We used a computational model-based approach that accounts for the short-term dynamics of closed-loop human cardiac control to integrate disparate experimental studies on neural adaptation following MI into a unified quantitative framework. We developed an ensemble of 59 distinct model parameterizations that account for the clinically observed heterogeneity of cardiac control in healthy individuals. We simulated an in silico cohort of 35,400 patients with MI, corresponding to six scenarios of one or more loci of neural adaptation coupled with cardiac remodelling. We evaluated the range of MI-induced shifts in arterial pressure, heart rate and baroreflex curve responses. Our results show that adaptation in any single neural locus coupled with cardiac remodelling is sufficient to account for the MI-induced haemodynamic and autonomic changes observed experimentally. Of the adaptation pathways, we found that individuals with central or peripheral vagal efferent adaptation and preserved baroreceptor gain could maintain high baroreflex sensitivity after ischaemic injury. These results suggest that there are a multitude of adaptive pathways for tuning the baroreflex circuit to shift cardiac control physiology, potentially explaining patient heterogeneity post-MI. KEY POINTS: Baroreflex sensitivity is a strong indicator of post-myocardial ischaemia survival and is variable among individuals. We fine-tuned a computational model ensemble based on physiological observations to develop an in silico patient cohort consistent with the range of baroreflex responses observed experimentally. Simulation and analysis of the in silico cohort show that individuals with a functional afferent pathway and the ability to adapt along the vagal efferent pathway can maintain baroreflex sensitivity post-cardiac ischaemia.

From sampling to simulating: Single-cell multiomics in systems pathophysiological modeling

As single-cell omics data sampling and acquisition methods have accumulated at an unprecedented rate, various data analysis pipelines have been developed for the inference of cell types, cell states and their distribution, state transitions, state trajectories, and state interactions. This presents a new opportunity in which single-cell omics data can be utilized to generate high-resolution, high-fidelity computational models. In this review, we discuss how single-cell omics data can be used to build computational models to simulate biological systems at various scales. We propose that single-cell data can be integrated with physiological information to generate organ-specific models, which can then be assembled to generate multi-organ systems pathophysiological models. Finally, we discuss how generic multi-organ models can be brought to the patient-specific level thus permitting their use in the clinical setting.

Unpacking the multimodal, multi-scale data of the fast and slow lanes of the cardiac vagus through computational modelling

The vagus nerve is a key mediator of brain-heart signaling, and its activity is necessary for cardiovascular health. Vagal outflow stems from the nucleus ambiguus, responsible primarily for fast, beat-to-beat regulation of heart rate and rhythm, and the dorsal motor nucleus of the vagus, responsible primarily for slow regulation of ventricular contractility. Due to the high-dimensional and multimodal nature of the anatomical, molecular and physiological data on neural regulation of cardiac function, data-derived mechanistic insights have proven elusive. Elucidating insights has been complicated further by the broad distribution of the data across heart, brain and peripheral nervous system circuits. Here we lay out an integrative framework based on computational modelling for combining these disparate and multi-scale data on the two vagal control lanes of the cardiovascular system. Newly available molecular-scale data, particularly single-cell transcriptomic analyses, have augmented our understanding of the heterogeneous neuronal states underlying vagally mediated fast and slow regulation of cardiac physiology. Cellular-scale computational models built from these data sets represent building blocks that can be combined using anatomical and neural circuit connectivity, neuronal electrophysiology, and organ/organismal-scale physiology data to create multi-system, multi-scale models that enable in silico exploration of the fast versus slow lane vagal stimulation. The insights from the computational modelling and analyses will guide new experimental questions on the mechanisms regulating the fast and slow lanes of the cardiac vagus toward exploiting targeted vagal neuromodulatory activity to promote cardiovascular health.

Unveiling inter-embryo variability in spindle length over time: Towards quantitative phenotype analysis

How can inter-individual variability be quantified? Measuring many features per experiment raises the question of choosing them to recapitulate high-dimensional data. Tackling this challenge on spindle elongation phenotypes, we showed that only three typical elongation patterns describe spindle elongation in C. elegans one-cell embryo. These archetypes, automatically extracted from the experimental data using principal component analysis (PCA), accounted for more than 95% of inter-individual variability of more than 1600 experiments across more than 100 different conditions. The two first archetypes were related to spindle average length and anaphasic elongation rate. The third archetype, accounting for 6% of the variability, was novel and corresponded to a transient spindle shortening in late metaphase, reminiscent of kinetochore function-defect phenotypes. Importantly, these three archetypes were robust to the choice of the dataset and were found even considering only non-treated conditions. Thus, the inter-individual differences between genetically perturbed embryos have the same underlying nature as natural inter-individual differences between wild-type embryos, independently of the temperatures. We thus propose that beyond the apparent complexity of the spindle, only three independent mechanisms account for spindle elongation, weighted differently in the various conditions. Interestingly, the spindle-length archetypes covered both metaphase and anaphase, suggesting that spindle elongation in late metaphase is sufficient to predict the late anaphase length. We validated this idea using a machine-learning approach. Finally, given amounts of these three archetypes could represent a quantitative phenotype. To take advantage of this, we set out to predict interacting genes from a seed based on the PCA coefficients. We exemplified this firstly on the role of tpxl-1 whose homolog tpx2 is involved in spindle microtubule branching, secondly the mechanism regulating metaphase length, and thirdly the central spindle players which set the length at anaphase. We found novel interactors not in public databases but supported by recent experimental publications.

Neuromodulatory co-expression in cardiac vagal motor neurons of the dorsal motor nucleus of the vagus

Vagal innervation is well known to be crucial to the maintenance of cardiac health, and to protect and recover the heart from injury. Only recently has this role been shown to depend on the activity of the underappreciated dorsal motor nucleus of the vagus (DMV). By combining neural tracing, transcriptomics, and anatomical mapping in male and female Sprague-Dawley rats, we characterize cardiac-specific neuronal phenotypes in the DMV. We find that the DMV cardiac-projecting neurons differentially express pituitary adenylate cyclase-activating polypeptide (PACAP), cocaine- and amphetamine-regulated transcript (CART), and synucleins, as well as evidence that they participate in neuromodulatory co-expression involving catecholamines. The significance of these findings is enhanced by previous knowledge of the role of PACAP at the heart and of the other neuromodulators in peripheral vagal targets.

Dynamic dysregulation of transcriptomic networks in brainstem autonomic nuclei during hypertension development in the female spontaneously hypertensive rat

Neurogenic hypertension stems from an imbalance in autonomic function that shifts the central cardiovascular control circuits toward a state of dysfunction. Using the female spontaneously hypertensive rat and the normotensive Wistar-Kyoto rat model, we compared the transcriptomic changes in three autonomic nuclei in the brainstem, nucleus of the solitary tract (NTS), caudal ventrolateral medulla, and rostral ventrolateral medulla (RVLM) in a time series at 8, 10, 12, 16, and 24 wk of age, spanning the prehypertensive stage through extended chronic hypertension. RNA-sequencing data were analyzed using an unbiased, dynamic pattern-based approach that uncovered dominant and several subtle differential gene regulatory signatures. Our results showed a persistent dysregulation across all three autonomic nuclei regardless of the stage of hypertension development as well as a cascade of transient dysregulation beginning in the RVLM at the prehypertensive stage that shifts toward the NTS at the hypertension onset. Genes that were persistently dysregulated were heavily enriched for immunological processes such as antigen processing and presentation, the adaptive immune response, and the complement system. Genes with transient dysregulation were also largely region-specific and were annotated for processes that influence neuronal excitability such as synaptic vesicle release, neurotransmitter transport, and an array of neuropeptides and ion channels. Our results demonstrate that neurogenic hypertension is characterized by brainstem region-specific transcriptomic changes that are highly dynamic with significant gene regulatory changes occurring at the hypertension onset as a key time window for dysregulation of homeostatic processes across the autonomic control circuits.NEW & NOTEWORTHY Hypertension is a major disease and is the primary risk factor for cardiovascular complications and stroke. The gene expression changes in the central nervous system circuits driving hypertension are understudied. Here, we show that coordinated and region-specific gene expression changes occur in the brainstem autonomic circuits over time during the development of a high blood pressure phenotype in a rat model of human essential hypertension.

Cell anatomy and network input explain differences within but not between leech touch cells at two different locations

Mechanosensory cells in the leech share several common features with mechanoreceptors in the human glabrous skin. Previous studies showed that the six T (touch) cells in each body segment of the leech are highly variable in their responses to somatic current injection and change their excitability over time. Here, we investigate three potential reasons for this variability in excitability by comparing the responses of T cells at two soma locations (T2 and T3): (1) Differential effects of time-dependent changes in excitability, (2) divergent synaptic input from the network, and (3) different anatomical structures. These hypotheses were explored with a combination of electrophysiological double recordings, 3D reconstruction of neurobiotin-filled cells, and compartmental model simulations. Current injection triggered significantly more spikes with shorter latency and larger amplitudes in cells at soma location T2 than at T3. During longer recordings, cells at both locations increased their excitability over time in the same way. T2 and T3 cells received the same amount of synaptic input from the unstimulated network, and the polysynaptic connections between both T cells were mutually symmetric. However, we found a striking anatomical difference: While in our data set all T2 cells innervated two roots connecting the ganglion with the skin, 50% of the T3 cells had only one root process. The sub-sample of T3 cells with one root process was significantly less excitable than the T3 cells with two root processes and the T2 cells. To test if the additional root process causes higher excitability, we simulated the responses of 3D reconstructed cells of both anatomies with detailed multi-compartment models. The anatomical subtypes do not differ in excitability when identical biophysical parameters and a homogeneous channel distribution are assumed. Hence, all three hypotheses may contribute to the highly variable T cell responses, but none of them is the only factor accounting for the observed systematic difference in excitability between cells at T2 vs. T3 soma location. Therefore, future patch clamp and modeling studies are needed to analyze how biophysical properties and spatial distribution of ion channels on the cell surface contribute to the variability and systematic differences of electrophysiological phenotypes.

Closed-loop modeling of central and intrinsic cardiac nervous system circuits underlying cardiovascular control

The baroreflex is a multi-input, multi-output control physiological system that regulates blood pressure by modulating nerve activity between the brainstem and the heart. Existing computational models of the baroreflex do not explictly incorporate the intrinsic cardiac nervous system (ICN), which mediates central control of the heart function. We developed a computational model of closed-loop cardiovascular control by integrating a network representation of the ICN within central control reflex circuits. We examined central and local contributions to the control of heart rate, ventricular functions, and respiratory sinus arrhythmia (RSA). Our simulations match the experimentally observed relationship between RSA and lung tidal volume. Our simulations predicted the relative contributions of the sensory and the motor neuron pathways to the experimentally observed changes in the heart rate. Our closed-loop cardiovascular control model is primed for evaluating bioelectronic interventions to treat heart failure and renormalize cardiovascular physiology.

From Seeing to Simulating: A Survey of Imaging Techniques and Spatially-Resolved Data for Developing Multiscale Computational Models of Liver Regeneration

Liver regeneration, which leads to the re-establishment of organ mass, follows a specifically organized set of biological processes acting on various time and length scales. Computational models of liver regeneration largely focused on incorporating molecular and signaling detail have been developed by multiple research groups in the recent years. These modeling efforts have supported a synthesis of disparate experimental results at the molecular scale. Incorporation of tissue and organ scale data using noninvasive imaging methods can extend these computational models towards a comprehensive accounting of multiscale dynamics of liver regeneration. For instance, microscopy-based imaging methods provide detailed histological information at the tissue and cellular scales. Noninvasive imaging methods such as ultrasound, computed tomography and magnetic resonance imaging provide morphological and physiological features including volumetric measures over time. In this review, we discuss multiple imaging modalities capable of informing computational models of liver regeneration at the organ-, tissue- and cellular level. Additionally, we discuss available software and algorithms, which aid in the analysis and integration of imaging data into computational models. Such models can be generated or tuned for an individual patient with liver disease. Progress towards integrated multiscale models of liver regeneration can aid in prognostic tool development for treating liver disease.

Single Cell Scale Neuronal and Glial Gene Expression and Putative Cell Phenotypes and Networks in the Nucleus Tractus Solitarius in an Alcohol Withdrawal Time Series

Alcohol withdrawal syndrome (AWS) is characterized by neuronal hyperexcitability, autonomic dysregulation, and severe negative emotion. The nucleus tractus solitarius (NTS) likely plays a prominent role in the neurological processes underlying these symptoms as it is the main viscerosensory nucleus in the brain. The NTS receives visceral interoceptive inputs, influences autonomic outputs, and has strong connections to the limbic system and hypothalamic-pituitary-adrenal axis to maintain homeostasis. Our prior analysis of single neuronal gene expression data from the NTS shows that neurons exist in heterogeneous transcriptional states that form distinct functional subphenotypes. Our working model conjectures that the allostasis secondary to alcohol dependence causes peripheral and central biological network decompensation in acute abstinence resulting in neurovisceral feedback to the NTS that substantially contributes to the observed AWS. We collected single noradrenergic and glucagon-like peptide-1 (GLP-1) neurons and microglia from rat NTS and measured a subset of their transcriptome as pooled samples in an alcohol withdrawal time series. Inflammatory subphenotypes predominate at certain time points, and GLP-1 subphenotypes demonstrated hyperexcitability post-withdrawal. We hypothesize such inflammatory and anxiogenic signaling contributes to alcohol dependence via negative reinforcement. Targets to mitigate such dysregulation and treat dependence can be identified from this dataset.

A single cell transcriptomics map of paracrine networks in the intrinsic cardiac nervous system

We developed a spatially-tracked single neuron transcriptomics map of an intrinsic cardiac ganglion, the right atrial ganglionic plexus (RAGP) that is a critical mediator of sinoatrial node (SAN) activity. This 3D representation of RAGP used neuronal tracing to extensively map the spatial distribution of the subset of neurons that project to the SAN. RNA-seq of laser capture microdissected neurons revealed a distinct composition of RAGP neurons compared to the central nervous system and a surprising finding that cholinergic and catecholaminergic markers are coexpressed, suggesting multipotential phenotypes that can drive neuroplasticity within RAGP. High-throughput qPCR of hundreds of laser capture microdissected single neurons confirmed these findings and revealed a high dimensionality of neuromodulatory factors that contribute to dynamic control of the heart. Neuropeptide-receptor coexpression analysis revealed a combinatorial paracrine neuromodulatory network within RAGP informing follow-on studies on the vagal control of RAGP to regulate cardiac function in health and disease.

Input-output signal processing plasticity of vagal motor neurons in response to cardiac ischemic injury

Vagal stimulation is emerging as the next frontier in bioelectronic medicine to modulate peripheral organ health and treat disease. The neuronal molecular phenotypes in the dorsal motor nucleus of the vagus (DMV) remain largely unexplored, limiting the potential for harnessing the DMV plasticity for therapeutic interventions. We developed a mesoscale single-cell transcriptomics data from hundreds of DMV neurons under homeostasis and following physiological perturbations. Our results revealed that homeostatic DMV neuronal states can be organized into distinguishable input-output signal processing units. Remote ischemic preconditioning induced a distinctive shift in the neuronal states toward diminishing the role of inhibitory inputs, with concomitant changes in regulatory microRNAs miR-218a and miR-495. Chronic cardiac ischemic injury resulted in a dramatic shift in DMV neuronal states suggestive of enhanced neurosecretory function. We propose a DMV molecular network mechanism that integrates combinatorial neurotransmitter inputs from multiple brain regions and humoral signals to modulate cardiac health.

Investigating the Effects of Brainstem Neuronal Adaptation on Cardiovascular Homeostasis

Central coordination of cardiovascular function is accomplished, in part, by the baroreceptor reflex, a multi-input multi-output physiological control system that regulates the activity of the parasympathetic and sympathetic nervous systems via interactions among multiple brainstem nuclei. Recent single-cell analyses within the brain revealed that individual neurons within and across brain nuclei exhibit distinct transcriptional states contributing to neuronal function. Such transcriptional heterogeneity complicates the task of understanding how neurons within and across brain nuclei organize and function to process multiple inputs and coordinate cardiovascular functions within the larger context of the baroreceptor reflex. However, prior analysis of brainstem neurons revealed that single-neuron transcriptional heterogeneity reflects an adaptive response to synaptic inputs and that neurons organize into distinct subtypes with respect to synaptic inputs received. Based on these results, we hypothesize that adaptation of neuronal subtypes support robust biological function through graded cellular responses. We test this hypothesis by examining the functional impact of neuronal adaptation on parasympathetic activity within the context of short-term baroreceptor reflex regulation. In this work, we extend existing quantitative closed-loop models of the baroreceptor reflex by incorporating into the model distinct input-driven neuronal subtypes and neuroanatomical groups that modulate parasympathetic activity. We then use this extended model to investigate, via simulation, the functional role of neuronal adaptation under conditions of health and systolic heart failure. Simulation results suggest that parasympathetic activity can be modulated appropriately by the coordination of distinct neuronal subtypes to maintain normal cardiovascular functions under systolic heart failure conditions. Moreover, differing degrees of adaptation of these neuronal subtypes contribute to cardiovascular behaviors corresponding to distinct clinical phenotypes of heart failure, such as exercise intolerance. Further, our results suggest that an imbalance between sympathetic and parasympathetic activity regulating ventricular contractility contributes to exercise intolerance in systolic heart failure patients, and restoring this balance can improve the short-term cardiovascular performance of these patients.

Combining Laser Capture Microdissection and Microfluidic qPCR to Analyze Transcriptional Profiles of Single Cells: A Systems Biology Approach to Opioid Dependence

Profound transcriptional heterogeneity in anatomically adjacent single cells suggests that robust tissue functionality may be achieved by cellular phenotype diversity. Single-cell experiments investigating the network dynamics of biological systems demonstrate cellular and tissue responses to various conditions at biologically meaningful resolution. Herein, we explain our methods for gathering single cells from anatomically specific locations and accurately measuring a subset of their gene expression profiles. We combine laser capture microdissection (LCM) with microfluidic reverse transcription quantitative polymerase chain reactions (RT-qPCR). We also use this microfluidic RT-qPCR platform to measure the microbial abundance of gut contents.

Single-Cell Gene Expression Analysis Identifies Chronic Alcohol-Mediated Shift in Hepatocyte Molecular States After Partial Hepatectomy

The analysis of molecular states of individual cells, as defined by their mRNA expression profiles and protein composition, has gained widespread interest in studying biological phenomena ranging from embryonic development to homeostatic tissue function and genesis and evolution of cancers. Although the molecular content of individual cells in a tissue can vary widely, their molecular states tend to be constrained within a transcriptional landscape partly described by the canonical archetypes of a population of cells. In this study, we sought to characterize the effects of an acute (partial hepatectomy) and chronic (alcohol consumption) perturbation on the molecular states of individual hepatocytes during the onset and progression of liver regeneration. We analyzed the expression of 84 genes across 233 individual hepatocytes acquired using laser capture microdissection. Analysis of the single-cell data revealed that hepatocyte molecular states can be considered as distributed across a set of four states irrespective of perturbation, with the proportions of hepatocytes in these states being dependent on the perturbation. In addition to the quiescent, primed, and replicating hepatocytes, we identified a fourth molecular state lying between the primed and replicating subpopulations. Comparison of the proportions of hepatocytes from each experimental condition in these four molecular states suggested that, in addition to aberrant priming, a slower transition from primed to replication state could contribute toward ethanol-mediated suppression of liver regenerative response to partial hepatectomy.

Cellular network modeling and single cell gene expression analysis reveals novel hepatic stellate cell phenotypes controlling liver regeneration dynamics

Background

Recent results from single cell gene and protein regulation studies are starting to uncover the previously underappreciated fact that individual cells within a population exhibit high variability in the expression of mRNA and proteins (i.e., molecular variability). By combining cellular network modeling, and high-throughput gene expression measurements in single cells, we seek to reconcile the high molecular variability in single cells with the relatively low variability in tissue-scale gene and protein expression and the highly coordinated functional responses of tissues to physiological challenges. In this study, we focus on relating the dynamic changes in distributions of hepatic stellate cell (HSC) functional phenotypes to the tightly regulated physiological response of liver regeneration.

Results

We develop a mathematical model describing contributions of HSC functional phenotype populations to liver regeneration and test model predictions through isolation and transcriptional characterization of single HSCs. We identify and characterize four HSC transcriptional states contributing to liver regeneration, two of which are described for the first time in this work. We show that HSC state populations change in vivo in response to acute challenges (in this case, 70% partial hepatectomy) and chronic challenges (chronic ethanol consumption). Our results indicate that HSCs influence the dynamics of liver regeneration through steady-state tissue preconditioning prior to an acute insult and through dynamic control of cell state balances. Furthermore, our modeling approach provides a framework to understand how balances among cell states influence tissue dynamics.

Conclusions

Taken together, our combined modeling and experimental studies reveal novel HSC transcriptional states and indicate that baseline differences in HSC phenotypes as well as a dynamic balance of transitions between these phenotypes control liver regeneration responses.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0605-7) contains supplementary material, which is available to authorized users.

Neurotransmitter identity and electrophysiological phenotype are genetically coupled in midbrain dopaminergic neurons

Most neuronal types have a well-identified electrical phenotype. It is now admitted that a same phenotype can be produced using multiple biophysical solutions defined by ion channel expression levels. This argues that systems-level approaches are necessary to understand electrical phenotype genesis and stability. Midbrain dopaminergic (DA) neurons, although quite heterogeneous, exhibit a characteristic electrical phenotype. However, the quantitative genetic principles underlying this conserved phenotype remain unknown. Here we investigated the quantitative relationships between ion channels’ gene expression levels in midbrain DA neurons using single-cell microfluidic qPCR. Using multivariate mutual information analysis to decipher high-dimensional statistical dependences, we unravel co-varying gene modules that link neurotransmitter identity and electrical phenotype. We also identify new segregating gene modules underlying the diversity of this neuronal population. We propose that the newly identified genetic coupling between neurotransmitter identity and ion channels may play a homeostatic role in maintaining the electrophysiological phenotype of midbrain DA neurons.

Single-Cell RNA Sequencing of Glioblastoma Cells

Single-cell RNA sequencing (sc-RNASeq) is a recently developed technique used to evaluate the transcriptome of individual cells. As opposed to conventional RNASeq in which entire populations are sequenced in bulk, sc-RNASeq can be beneficial when trying to better understand gene expression patterns in markedly heterogeneous populations of cells or when trying to identify transcriptional signatures of rare cells that may be underrepresented when using conventional bulk RNASeq. In this method, we describe the generation and analysis of cDNA libraries from single patient-derived glioblastoma cells using the C1 Fluidigm system. The protocol details the use of the C1 integrated fluidics circuit (IFC) for capturing, imaging and lysing cells; performing reverse transcription; and generating cDNA libraries that are ready for sequencing and analysis.

Novel Influences of IL-10 on CNS Inflammation Revealed by Integrated Analyses of Cytokine Networks and Microglial Morphology

Coordinated interactions between cytokine signaling and morphological dynamics of microglial cells regulate neuroinflammation in CNS injury and disease. We found that pro-inflammatory cytokine gene expression in vivo showed a pronounced recovery following systemic LPS. We performed a novel multivariate analysis of microglial morphology and identified changes in specific morphological properties of microglia that matched the expression dynamics of pro-inflammatory cytokine TNFα. The adaptive recovery kinetics of TNFα expression and microglial soma size showed comparable profiles and dependence on anti-inflammatory cytokine IL-10 expression. The recovery of cytokine variations and microglial morphology responses to inflammation were negatively regulated by IL-10. Our novel morphological analysis of microglia is able to detect subtle changes and can be used widely. We implemented in silico simulations of cytokine network dynamics which showed—counter-intuitively, but in line with our experimental observations—that negative feedback from IL-10 was sufficient to impede the adaptive recovery of TNFα-mediated inflammation. Our integrative approach is a powerful tool to study changes in specific components of microglial morphology for insights into their functional states, in relation to cytokine network dynamics, during CNS injury and disease.

A data-driven modeling approach to identify disease-specific multi-organ networks driving physiological dysregulation

Multiple physiological systems interact throughout the development of a complex disease. Knowledge of the dynamics and connectivity of interactions across physiological systems could facilitate the prevention or mitigation of organ damage underlying complex diseases, many of which are currently refractory to available therapeutics (e.g., hypertension). We studied the regulatory interactions operating within and across organs throughout disease development by integrating in vivo analysis of gene expression dynamics with a reverse engineering approach to infer data-driven dynamic network models of multi-organ gene regulatory influences. We obtained experimental data on the expression of 22 genes across five organs, over a time span that encompassed the development of autonomic nervous system dysfunction and hypertension. We pursued a unique approach for identification of continuous-time models that jointly described the dynamics and structure of multi-organ networks by estimating a sparse subset of ∼12,000 possible gene regulatory interactions. Our analyses revealed that an autonomic dysfunction-specific multi-organ sequence of gene expression activation patterns was associated with a distinct gene regulatory network. We analyzed the model structures for adaptation motifs, and identified disease-specific network motifs involving genes that exhibited aberrant temporal dynamics. Bioinformatic analyses identified disease-specific single nucleotide variants within or near transcription factor binding sites upstream of key genes implicated in maintaining physiological homeostasis. Our approach illustrates a novel framework for investigating the pathogenesis through model-based analysis of multi-organ system dynamics and network properties. Our results yielded novel candidate molecular targets driving the development of cardiovascular disease, metabolic syndrome, and immune dysfunction.

Complex diseases such as hypertension often involve maladaptive autonomic nervous system control over the cardiovascular, renal, hepatic, immune, and endocrine systems. We studied the pathogenesis of physiological homeostasis by examining the temporal dynamics of gene expression levels from multiple organs in an animal model of autonomic dysfunction characterized by cardiovascular disease, metabolic dysregulation, and immune system aberrations. We employed a data-driven modeling approach to jointly predict continuous gene expression dynamics and gene regulatory interactions across organs in the disease and control phenotypes. We combined our analyses of multi-organ gene regulatory network dynamics and connectivity with bioinformatic analyses of genetic mutations that could regulate gene expression. Our multi-organ modeling approach to investigate the mechanisms of complex disease pathogenesis revealed novel candidates for therapeutic interventions against the development and progression of complex diseases involving autonomic nervous system dysfunction.

Embryonic development of selectively vulnerable neurons in Parkinson's disease

A specific set of brainstem nuclei are susceptible to degeneration in Parkinson's disease. We hypothesise that neuronal vulnerability reflects shared phenotypic characteristics that confer selective vulnerability to degeneration. Neuronal phenotypic specification is mainly the cumulative result of a transcriptional regulatory program that is active during the development. By manual curation of the developmental biology literature, we comprehensively reconstructed an anatomically resolved cellular developmental lineage for the adult neurons in five brainstem regions that are selectively vulnerable to degeneration in prodromal or early Parkinson's disease. We synthesised the literature on transcription factors that are required to be active, or required to be inactive, in the development of each of these five brainstem regions, and at least two differentially vulnerable nuclei within each region. Certain transcription factors, e.g., Ascl1 and Lmx1b, seem to be required for specification of many brainstem regions that are susceptible to degeneration in early Parkinson's disease. Some transcription factors can even distinguish between differentially vulnerable nuclei within the same brain region, e.g., Pitx3 is required for specification of the substantia nigra pars compacta, but not the ventral tegmental area. We do not suggest that Parkinson's disease is a developmental disorder. In contrast, we consider identification of shared developmental trajectories as part of a broader effort to identify the molecular mechanisms that underlie the phenotypic features that are shared by selectively vulnerable neurons. Systematic in vivo assessment of fate determining transcription factors should be completed for all neuronal populations vulnerable to degeneration in early Parkinson's disease.

Modeling cytokine regulatory network dynamics driving neuroinflammation in central nervous system disorders

A central goal of pharmacological efforts to treat central nervous system (CNS) diseases is to develop systemic therapeutics that can restore CNS homeostasis. Achieving this goal requires a fundamental understanding of CNS function within the organismal context so as to leverage the mechanistic insights on the molecular basis of cellular and tissue functions towards novel drug target identification. The immune system constitutes a key link between the periphery and CNS, and many neurological disorders and neurodegenerative diseases are characterized by immune dysfunction. We review the salient opportunities for applying computational models to CNS disease research, and summarize relevant approaches from studies of immune function and neuroinflammation. While the accurate prediction of disease-related phenomena is often considered the central goal of modeling studies, we highlight the utility of computational modeling applications beyond making predictions, particularly for drawing counterintuitive insights from model-based analysis of multi-parametric and time series data sets.

Assessing characteristics of RNA amplification methods for single cell RNA sequencing

BACKGROUND

Recently, measurement of RNA at single cell resolution has yielded surprising insights. Methods for single-cell RNA sequencing (scRNA-seq) have received considerable attention, but the broad reliability of single cell methods and the factors governing their performance are still poorly known.

RESULTS

Here, we conducted a large-scale control experiment to assess the transfer function of three scRNA-seq methods and factors modulating the function. All three methods detected greater than 70% of the expected number of genes and had a 50% probability of detecting genes with abundance greater than 2 to 4 molecules. Despite the small number of molecules, sequencing depth significantly affected gene detection. While biases in detection and quantification were qualitatively similar across methods, the degree of bias differed, consistent with differences in molecular protocol. Measurement reliability increased with expression level for all methods and we conservatively estimate measurements to be quantitative at an expression level greater than ~5-10 molecules.

CONCLUSIONS

Based on these extensive control studies, we propose that RNA-seq of single cells has come of age, yielding quantitative biological information.

Single-Cell Transcriptional Analysis Reveals Novel Neuronal Phenotypes and Interaction Networks Involved in the Central Circadian Clock

Single-cell heterogeneity confounds efforts to understand how a population of cells organizes into cellular networks that underlie tissue-level function. This complexity is prominent in the mammalian suprachiasmatic nucleus (SCN). Here, individual neurons exhibit a remarkable amount of asynchronous behavior and transcriptional heterogeneity. However, SCN neurons are able to generate precisely coordinated synaptic and molecular outputs that synchronize the body to a common circadian cycle by organizing into cellular networks. To understand this emergent cellular network property, it is important to reconcile single-neuron heterogeneity with network organization. In light of recent studies suggesting that transcriptionally heterogeneous cells organize into distinct cellular phenotypes, we characterized the transcriptional, spatial, and functional organization of 352 SCN neurons from mice experiencing phase-shifts in their circadian cycle. Using the community structure detection method and multivariate analytical techniques, we identified previously undescribed neuronal phenotypes that are likely to participate in regulatory networks with known SCN cell types. Based on the newly discovered neuronal phenotypes, we developed a data-driven neuronal network structure in which multiple cell types interact through known synaptic and paracrine signaling mechanisms. These results provide a basis from which to interpret the functional variability of SCN neurons and describe methodologies toward understanding how a population of heterogeneous single cells organizes into cellular networks that underlie tissue-level function.

Cellular Deconstruction: Finding Meaning in Individual Cell Variation

The advent of single cell transcriptome analysis has permitted the discovery of cell-to-cell variation in transcriptome expression of even presumptively identical cells. We hypothesize that this variability reflects a many-to-one relation between transcriptome states and the phenotype of a cell. In this relation, the molecular ratios of the subsets of RNA are determined by the stoichiometric constraints of the cell systems, which underdetermine the transcriptome state. Furthermore, the variability is, in part, induced by the tissue context and is important for system-level function. This theory is analogous to theories of literary deconstruction, where multiple 'signifiers' work in opposition to one another to create meaning. By analogy, transcriptome phenotypes should be defined as subsets of RNAs comprising selected RNA systems where the system-associated RNAs are balanced with each other to produce the associated cellular function. This idea provides a framework for understanding cellular heterogeneity in phenotypic responses to variant conditions, such as disease challenge.

Molecular variability elicits a tunable switch with discrete neuromodulatory response phenotypes

Recent single cell studies show extensive molecular variability underlying cellular responses. We evaluated the impact of molecular variability in the expression of cell signaling components and ion channels on electrophysiological excitability and neuromodulation. We employed a computational approach that integrated neuropeptide receptor-mediated signaling with electrophysiology. We simulated a population of neurons in which expression levels of a neuropeptide receptor and multiple ion channels were simultaneously varied within a physiological range. We analyzed the effects of variation on the electrophysiological response to a neuropeptide stimulus. Our results revealed distinct response patterns associated with low versus high receptor levels. Neurons with low receptor levels showed increased excitability and neurons with high receptor levels showed reduced excitability. These response patterns were separated by a narrow receptor level range forming a separatrix. The position of this separatrix was dependent on the expression levels of multiple ion channels. To assess the relative contributions of receptor and ion channel levels to the response profiles, we categorized the responses into six phenotypes based on response kinetics and magnitude. We applied several multivariate statistical approaches and found that receptor and channel expression levels influence the neuromodulation response phenotype through a complex though systematic mapping. Our analyses extended our understanding of how cellular responses to neuromodulation vary as a function of molecular expression. Our study showed that receptor expression and biophysical state interact with distinct relative contributions to neuronal excitability.

Variation is function: Are single cell differences functionally important?: Testing the hypothesis that single cell variation is required for aggregate function

There is a growing appreciation of the extent of transcriptome variation across individual cells of the same cell type. While expression variation may be a byproduct of, for example, dynamic or homeostatic processes, here we consider whether single-cell molecular variation per se might be crucial for population-level function. Under this hypothesis, molecular variation indicates a diversity of hidden functional capacities within an ensemble of identical cells, and this functional diversity facilitates collective behavior that would be inaccessible to a homogenous population. In reviewing this topic, we explore possible functions that might be carried by a heterogeneous ensemble of cells; however, this question has proven difficult to test, both because methods to manipulate molecular variation are limited and because it is complicated to define, and measure, population-level function. We consider several possible methods to further pursue the hypothesis that variation is function through the use of comparative analysis and novel experimental techniques.

Intracellular Information Processing through Encoding and Decoding of Dynamic Signaling Features

Cell signaling dynamics and transcriptional regulatory activities are variable within specific cell types responding to an identical stimulus. In addition to studying the network interactions, there is much interest in utilizing single cell scale data to elucidate the non-random aspects of the variability involved in cellular decision making. Previous studies have considered the information transfer between the signaling and transcriptional domains based on an instantaneous relationship between the molecular activities. These studies predict a limited binary on/off encoding mechanism which underestimates the complexity of biological information processing, and hence the utility of single cell resolution data. Here we pursue a novel strategy that reformulates the information transfer problem as involving dynamic features of signaling rather than molecular abundances. We pursue a computational approach to test if and how the transcriptional regulatory activity patterns can be informative of the temporal history of signaling. Our analysis reveals (1) the dynamic features of signaling that significantly alter transcriptional regulatory patterns (encoding), and (2) the temporal history of signaling that can be inferred from single cell scale snapshots of transcriptional activity (decoding). Immediate early gene expression patterns were informative of signaling peak retention kinetics, whereas transcription factor activity patterns were informative of activation and deactivation kinetics of signaling. Moreover, the information processing aspects varied across the network, with each component encoding a selective subset of the dynamic signaling features. We developed novel sensitivity and information transfer maps to unravel the dynamic multiplexing of signaling features at each of these network components. Unsupervised clustering of the maps revealed two groups that aligned with network motifs distinguished by transcriptional feedforward vs feedback interactions. Our new computational methodology impacts the single cell scale experiments by identifying downstream snapshot measures required for inferring specific dynamical features of upstream signals involved in the regulation of cellular responses.

Single cell studies have shown that differential patterns in the dynamics of signaling proteins, transcription factor activity, gene expression, etc. produce distinct downstream outcomes. The opposite also holds true where particular cellular outcomes have been found to be associated with the dynamical pattern of one or more signaling molecules. Signaling pathways, therefore, serve as signal processing units to inform specific downstream regulation. However, the functional capabilities of the dynamic aspects of signaling are not well understood. To address this issue, we developed a new approach that evaluates information processing between dynamic features in signaling patterns and transcriptional regulatory activity. Our work demonstrates that the information transfer occur through decoding of temporal history of signals rather than only through instantaneous correlations. Moreover, our results identify regulatory network motifs as the critical components in the information processing and filtering of variability in signaling dynamics to produce distinct patterns of downstream transcriptional responses. Our methodology can be broadly applied to single cell scale data on experimentally accessible downstream measures to infer dynamic aspects of upstream signaling.

In Vivo Zonal Variation and Liver Cell-Type Specific NF-κB Localization after Chronic Adaptation to Ethanol and following Partial Hepatectomy

NF-κB is a major inflammatory response mediator in the liver, playing a key role in the pathogenesis of alcoholic liver injury. We investigated zonal as well as liver cell type-specific distribution of NF-κB activation across the liver acinus following adaptation to chronic ethanol intake and 70% partial hepatectomy (PHx). We employed immunofluorescence staining, digital image analysis and statistical distributional analysis to quantify subcellular localization of NF-κB in hepatocytes and hepatic stellate cells (HSCs). We detected significant spatial heterogeneity of NF-κB expression and cellular localization between cytoplasm and nucleus across liver tissue. Our main aims involved investigating the zonal bias in NF-κB localization and determining to what extent chronic ethanol intake affects this zonal bias with in hepatocytes at baseline and post-PHx. Hepatocytes in the periportal area showed higher NF-κB expression than in the pericentral region in the carbohydrate-fed controls, but not in the ethanol group. However, the distribution of NF-κB nuclear localization in hepatocytes was shifted towards higher levels in pericentral region than in periportal area, across all treatment conditions. Chronic ethanol intake shifted the NF-κB distribution towards higher nuclear fraction in hepatocytes as compared to the pair-fed control group. Ethanol also stimulated higher NF-κB expression in a subpopulation of HSCs. In the control group, PHx elicited a shift towards higher NF-κB nuclear fraction in hepatocytes. However, this distribution remained unchanged in the ethanol group post-PHx. HSCs showed a lower NF-κB expression following PHx in both ethanol and control groups. We conclude that adaptation to chronic ethanol intake attenuates the liver zonal variation in NF-κB expression and limits the PHx-induced NF-κB activation in hepatocytes, but does not alter the NF-κB expression changes in HSCs in response to PHx. Our findings provide new insights as to how ethanol treatment may affect cell-type specific processes regulated by NF-κB activation in liver cells.

Multiscale model of dynamic neuromodulation integrating neuropeptide-induced signaling pathway activity with membrane electrophysiology

We developed a multiscale model to bridge neuropeptide receptor-activated signaling pathway activity with membrane electrophysiology. Typically, the neuromodulation of biochemical signaling and biophysics have been investigated separately in modeling studies. We studied the effects of Angiotensin II (AngII) on neuronal excitability changes mediated by signaling dynamics and downstream phosphorylation of ion channels. Experiments have shown that AngII binding to the AngII receptor type-1 elicits baseline-dependent regulation of cytosolic Ca(2+) signaling. Our model simulations revealed a baseline Ca(2+)-dependent response to AngII receptor type-1 activation by AngII. Consistent with experimental observations, AngII evoked a rise in Ca(2+) when starting at a low baseline Ca(2+) level, and a decrease in Ca(2+) when starting at a higher baseline. Our analysis predicted that the kinetics of Ca(2+) transport into the endoplasmic reticulum play a critical role in shaping the Ca(2+) response. The Ca(2+) baseline also influenced the AngII-induced excitability changes such that lower Ca(2+) levels were associated with a larger firing rate increase. We examined the relative contributions of signaling kinases protein kinase C and Ca(2+)/Calmodulin-dependent protein kinase II to AngII-mediated excitability changes by simulating activity blockade individually and in combination. We found that protein kinase C selectively controlled firing rate adaptation whereas Ca(2+)/Calmodulin-dependent protein kinase II induced a delayed effect on the firing rate increase. We tested whether signaling kinetics were necessary for the dynamic effects of AngII on excitability by simulating three scenarios of AngII-mediated KDR channel phosphorylation: (1), an increased steady state; (2), a step-change increase; and (3), dynamic modulation. Our results revealed that the kinetics emerging from neuromodulatory activation of the signaling network were required to account for the dynamical changes in excitability. In summary, our integrated multiscale model provides, to our knowledge, a new approach for quantitative investigation of neuromodulatory effects on signaling and electrophysiology.

MicroRNA network changes in the brain stem underlie the development of hypertension

Hypertension is a major chronic disease whose molecular mechanisms remain poorly understood. We compared neuroanatomical patterns of microRNAs in the brain stem of the spontaneous hypertensive rat (SHR) to the Wistar Kyoto rat (WKY, control). We quantified 419 well-annotated microRNAs in the nucleus of the solitary tract (NTS) and rostral ventrolateral medulla (RVLM), from SHR and WKY rats, during three main stages of hypertension development. Changes in microRNA expression were stage- and region-dependent, with a majority of SHR vs. WKY differential expression occurring at the hypertension onset stage in NTS versus at the prehypertension stage in RVLM. Our analysis identified 24 microRNAs showing time-dependent differential expression in SHR compared with WKY in at least one brain region. We predicted potential gene regulatory targets corresponding to catecholaminergic processes, neuroinflammation, and neuromodulation using the miRWALK and RNA22 databases, and we tested those bioinformatics predictions using high-throughput quantitative PCR to evaluate correlations of differential expression between the microRNAs and their predicted gene targets. We found a novel regulatory network motif consisting of microRNAs likely downregulating a negative regulator of prohypertensive processes such as angiotensin II signaling and leukotriene-based inflammation. Our results provide new evidence on the dynamics of microRNA expression in the development of hypertension and predictions of microRNA-mediated regulatory networks playing a region-dependent role in potentially altering brain-stem cardiovascular control circuit function leading to the development of hypertension.

Identifying functional gene regulatory network phenotypes underlying single cell transcriptional variability

Recent analysis of single-cell transcriptomic data has revealed a surprising organization of the transcriptional variability pervasive across individual neurons. In response to distinct combinations of synaptic input-type, a new organization of neuronal subtypes emerged based on transcriptional states that were aligned along a gradient of correlated gene expression. Individual neurons traverse across these transcriptional states in response to cellular inputs. However, the regulatory network interactions driving these changes remain unclear. Here we present a novel fuzzy logic-based approach to infer quantitative gene regulatory network models from highly variable single-cell gene expression data. Our approach involves developing an a priori regulatory network that is then trained against in vivo single-cell gene expression data in order to identify causal gene interactions and corresponding quantitative model parameters. Simulations of the inferred gene regulatory network response to experimentally observed stimuli levels mirrored the pattern and quantitative range of gene expression across individual neurons remarkably well. In addition, the network identification results revealed that distinct regulatory interactions, coupled with differences in the regulatory network stimuli, drive the variable gene expression patterns observed across the neuronal subtypes. We also identified a key difference between the neuronal subtype-specific networks with respect to negative feedback regulation, with the catecholaminergic subtype network lacking such interactions. Furthermore, by varying regulatory network stimuli over a wide range, we identified several cases in which divergent neuronal subtypes could be driven towards similar transcriptional states by distinct stimuli operating on subtype-specific regulatory networks. Based on these results, we conclude that heterogenous single-cell gene expression profiles should be interpreted through a regulatory network modeling perspective in order to separate the contributions of network interactions from those of cellular inputs.