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.

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

NEW FINDINGS: What is the topic of this review? The vagus nerve is a crucial regulator of cardiovascular homeostasis, and its activity is linked to heart health. Vagal activity originates from two brainstem nuclei: the nucleus ambiguus (fast lane) and the dorsal motor nucleus of the vagus (slow lane), nicknamed for the time scales that they require to transmit signals. What advances does it highlight? Computational models are powerful tools for organizing multi-scale, multimodal data on the fast and slow lanes in a physiologically meaningful way. A strategy is laid out for how these models can guide experiments aimed at harnessing the cardiovascular health benefits of differential activation of the fast and slow lanes. ABSTRACT: 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.

Systemic Immune Bias Delineates Malignant Astrocytoma Survival Cohorts.

Patients with grade III anaplastic astrocytomas (AA) separate into survival cohorts based on the presence or absence of mutations in isocitrate dehydrogenase (IDH). Progression to glioblastoma (GBM), morphologically distinguishable by elevated microvascular proliferation, necrosis, and cell division in tumor tissues, is considerably more rapid in IDH wild-type tumors such that their diagnosis as AA is relatively rare. More often initially presenting as GBM, these contain higher numbers of tumor-associated macrophages (TAMs) than most AA, and GBM patients also have higher levels of circulating M2 monocytes. TAM and M2 monocytes share functional properties inhibitory for antitumor immunity. Yet, although there is a wealth of data implicating TAM in tumor-immune evasion, there has been limited analysis of the impact of the circulating M2 monocytes. In the current study, immune parameters in sera, circulating cells, and tumor tissues from patients with primary gliomas morphologically diagnosed as AA were assessed. Profound differences in serum cytokines, glioma extracellular vesicle cross-reactive Abs, and gene expression by circulating cells identified two distinct patient cohorts. Evidence of type 2-immune bias was most often seen in patients with IDH wild-type AA, whereas a type 1 bias was common in patients with tumors expressing the IDH1R132H mutation. Nevertheless, a patient's immune profile was better correlated with the extent of tumor vascular enhancement on magnetic resonance imaging than IDH mutational status. Regardless of IDH genotype, AA progression appears to be associated with a switch in systemic immune bias from type 1 to type 2 and the loss of tumor vasculature integrity.

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.

Systemic leukotriene B4 receptor antagonism lowers arterial blood pressure and improves autonomic function in the spontaneously hypertensive rat.

KEY POINTS: Evidence indicates an association between hypertension and chronic systemic inflammation in both human hypertension and experimental animal models. Previous studies in the spontaneously hypertensive rat (SHR) support a role for leukotriene B4 (LTB4 ), a potent chemoattractant involved in the inflammatory response, but its mode of action is poorly understood. In the SHR, we observed an increase in T cells and macrophages in the brainstem; in addition, gene expression profiling data showed that LTB4 production, degradation and downstream signalling in the brainstem of the SHR are dynamically regulated during hypertension. When LTB4 receptor 1 (BLT1) receptors were blocked with CP-105,696, arterial pressure was reduced in the SHR compared to the normotensive control and this reduction was associated with a significant decrease in systolic blood pressure (BP) indicators. These data provide new evidence for the role of LTB4 as an important neuro-immune pathway in the development of hypertension and therefore may serve as a novel therapeutic target for the treatment of neurogenic hypertension. ABSTRACT: Accumulating evidence indicates an association between hypertension and chronic systemic inflammation in both human hypertension and experimental animal models. Previous studies in the spontaneously hypertensive rat (SHR) support a role for leukotriene B4 (LTB4 ), a potent chemoattractant involved in the inflammatory response. However, the mechanism for LTB4 -mediated inflammation in hypertension is poorly understood. Here we report in the SHR, increased brainstem infiltration of T cells and macrophages plus gene expression profiling data showing that LTB4 production, degradation and downstream signalling in the brainstem of the SHR are dynamically regulated during hypertension. Chronic blockade of the LTB4 receptor 1 (BLT1) receptor with CP-105,696, reduced arterial pressure in the SHR compared to the normotensive control and this reduction was associated with a significant decrease in low and high frequency spectra of systolic blood pressure, and an increase in spontaneous baroreceptor reflex gain (sBRG). These data provide new evidence for the role of LTB4 as an important neuro-immune pathway in the development of hypertension and therefore may serve as a novel therapeutic target for the treatment of neurogenic hypertension.

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.

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.

Coordinated dynamic gene expression changes in the central nucleus of the amygdala during alcohol withdrawal.

BACKGROUND: Chronic alcohol use causes widespread changes in the cellular biology of the amygdala's central nucleus (CeA), a GABAergic center that integrates autonomic physiology with the emotional aspects of motivation and learning. While alcohol-induced neurochemical changes play a role in dependence and drinking behavior, little is known about the CeA's dynamic changes during withdrawal, a period of emotional and physiologic disturbance. METHODS: We used a qRT-PCR platform to measure 139 transcripts in 92 rat CeA samples from control (N = 33), chronically alcohol exposed (N = 26), and withdrawn rats (t = 4, 8, 18, 32, and 48 hours; N = 5, 10, 7, 6, 5). This focused transcript set allowed us to identify significant dynamic expression patterns during the first 48 hours of withdrawal and propose potential regulatory mechanisms. RESULTS: Chronic alcohol exposure causes a limited number of small magnitude expression changes. In contrast, withdrawal results in a greater number of large changes within 4 hours of removal of the alcohol diet. Sixty-five of the 139 measured transcripts (47%) showed differential regulation during withdrawal. Over the 48-hour period, dynamic changes in the expression of γ-aminobutyric acid type A (GABA(A) ), ionotropic glutamate and neuropeptide system-related G-protein-coupled receptor subunits, and the Ras/Raf signaling pathway were seen as well as downstream transcription factors (TFs) and epigenetic regulators. Four temporally correlated gene clusters were identified with shared functional roles including NMDA receptors, MAPKKK and chemokine signaling cascades, and mediators of long-term potentiation, among others. Cluster promoter regions shared overrepresented binding sites for multiple TFs including Cebp, Usf-1, Smad3, Ap-2, and c-Ets, suggesting a potential regulatory role. CONCLUSIONS: During alcohol withdrawal, the CeA experiences rapid changes in mRNA expression of these functionally related transcripts that were not predicted by measurement during chronic exposure. This study provides new insight into dynamic expression changes during alcohol withdrawal and suggests novel regulatory relationships that potentially impact the aspects of emotional modulation.

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.

Temporal changes in innate immune signals in a rat model of alcohol withdrawal in emotional and cardiorespiratory homeostatic nuclei.

BACKGROUND: Chronic alcohol use changes the brain's inflammatory state. However, there is little work examining the progression of the cytokine response during alcohol withdrawal, a period of profound autonomic and emotional upset. This study examines the inflammatory response in the central nucleus of the amygdala (CeA) and dorsal vagal complex (DVC), brain regions neuroanatomically associated with affective and cardiorespiratory regulation in an in vivo rat model of withdrawal following a single chronic exposure. METHODS: For qRT-PCR studies, we measured the expression of TNF-α, NOS-2, Ccl2 (MCP-1), MHC II invariant chain CD74, and the TNF receptor Tnfrsf1a in CeA and DVC samples from adult male rats exposed to a liquid alcohol diet for thirty-five days and in similarly treated animals at four hours and forty-eight hours following alcohol withdrawal. ANOVA was used to identify statistically significant treatment effects. Immunohistochemistry (IHC) and confocal microscopy were performed in a second set of animals during chronic alcohol exposure and subsequent 48-hour withdrawal. RESULTS: Following a chronic alcohol exposure, withdrawal resulted in a statistically significant increase in the expression of mRNAs specific for innate immune markers Ccl2, TNF-α, NOS-2, Tnfrsf1a, and CD74. This response was present in both the CeA and DVC and most prominent at 48 hours. Confocal IHC of samples taken 48 hours into withdrawal demonstrate the presence of TNF-α staining surrounding cells expressing the neural marker NeuN and endothelial cells colabeled with ICAM-1 (CD54) and RECA-1, markers associated with an inflammatory response. Again, findings were consistent in both brain regions. CONCLUSIONS: This study demonstrates the rapid induction of Ccl2, TNF-α, NOS-2, Tnfrsf1a and CD74 expression during alcohol withdrawal in both the CeA and DVC. IHC dual labeling showed an increase in TNF-α surrounding neurons and ICAM-1 on vascular endothelial cells 48 hours into withdrawal, confirming the inflammatory response at the protein level. These findings suggest that an abrupt cessation of alcohol intake leads to an acute central nervous system (CNS) inflammatory response in these regions that regulate autonomic and emotional state.

Adaptive transcriptional dynamics of A2 neurons and central cardiovascular control pathways.

In order to provide an example of useful interaction between systems biology and computational neuroscience traditions, here we aim to identify the molecular response process through which elevated blood pressure induces a temporal sequence of gene expression changes. This initial response may then continually evolve as an adaptive response, altering central blood pressure set-point control. Our approach involves using a set of 96 quantitative PCR gene assays, which represent the molecular process associated with neuronal responses to angiotensin II type 1 receptor (AT1R) activation. We use this set as a probe to search for the AT1R signalling-triggered expression programme in individual neurons and groups of neurons involved in homeostatic regulation of cardiorespiratory function. Specifically, we focus on tissue samples from the nucleus tractus solitarii, as well as groups of A2 neurons and individual A2 cells within the nucleus tractus solitarii, and in the ventrolateral medulla and central amygdala. We assay these neural samples at rest and in response to elevated blood pressure. Analysis of the resulting high-dimensional data set reveals a remarkable complexity and heterogeneity of the samples and of their response to changes in blood pressure. These results demonstrate differential expression programmes for each anatomically distinct neuronal group and neuronal type. Single-cell expression analysis shows that A2 cells also are variable, and that the subset that responds to blood pressure with elevated Fos expression differs from other A2 cells. We present models of gene regulatory networks and of signalling cascades related to AT1R and broadly discuss the opportunities for valuable interactions between systems biology and computational neuroscience.

Rapid temporal changes in the expression of a set of neuromodulatory genes during alcohol withdrawal in the dorsal vagal complex: molecular evidence of homeostatic disturbance.

BACKGROUND: Chronic alcohol exposure produces neuroadaptation, which increases the risk of cellular excitotoxicity and autonomic dysfunction during withdrawal. The temporal progression and regulation of the gene expression that contributes to this physiologic and behavioral phenotype is poorly understood early in the withdrawal period. Further, it is unexplored in the dorsal vagal complex (DVC), a brainstem autonomic regulatory structure. METHODS: We use a quantitative polymerase chain reaction platform to precisely and simultaneously measure the expression of 145 neuromodulatory genes in more than 100 rat DVC samples from control, chronically alcohol-exposed, and withdrawn rats. To gain insight into the dynamic progression and regulation of withdrawal, we focus on the expression of a subset of functionally relevant genes during the first 48 hours, when behavioral symptoms are most severe. RESULTS: In the DVC, expression of this gene subset is essentially normal in chronically alcohol-exposed rats. However, withdrawal results in rapid, large-magnitude expression changes in this group. We observed differential regulation in 86 of the 145 genes measured (59%), some as early as 4 hours into withdrawal. Time series measurements (4, 8, 18, 32, and 48 hours after alcohol removal) revealed dynamic expression responses in immediate early genes, γ-aminobutyric acid type A, ionotropic glutamate, and G-protein coupled receptors and the Ras/Raf signaling pathway. Together, these changes elucidate a complex, temporally coordinated response that involves correlated expression of many functionally related groups. In particular, the expression patterns of Gabra1, Grin2a, Grin3a, and Grik3 were tightly correlated. These receptor subunits share overrepresented transcription factor binding sites for Pax-8 and other transcription factors, suggesting a common regulatory mechanism and a role for these transcription factors in the regulation of neurotransmission within the first 48 hours of alcohol withdrawal. CONCLUSIONS: Expression in this gene set is essentially normal in the alcohol-adapted DVC, but withdrawal results in immediate, large-magnitude, and dynamic changes. These data support both increased research focus on the biological ramifications of alcohol withdrawal and enable novel insights into the dynamic withdrawal expression response in this understudied homeostatic control center.

Investigating Drivers of Antireward in Addiction Behavior with Anatomically Specific Single-Cell Gene Expression Methods.

Increasing rates of addiction behavior have motivated mental health researchers and clinicians alike to understand antireward and recovery. This shift away from reward and commencement necessitates novel perspectives, paradigms, and hypotheses along with an expansion of the methods applied to investigate addiction. Here, we provide an example: A systems biology approach to investigate antireward that combines laser capture microdissection (LCM) and high-throughput microfluidic reverse transcription quantitative polymerase chain reactions (RT-qPCR). Gene expression network dynamics were measured and a key driver of neurovisceral dysregulation in alcohol and opioid withdrawal, neuroinflammation, was identified. This combination of technologies provides anatomic and phenotypic specificity at single-cell resolution with high-throughput sensitivity and specific gene expression measures yielding both hypothesis-generating datasets and mechanistic possibilities that generate opportunities for novel insights and treatments.

Similarities in alcohol and opioid withdrawal syndromes suggest common negative reinforcement mechanisms involving the interoceptive antireward pathway.

Alcohol and opioids are two major contributors to so-called deaths of despair. Though the effects of these substances on mammalian systems are distinct, commonalities in their withdrawal syndromes suggest a shared pathophysiology. For example, both are characterized by marked autonomic dysregulation and are treated with alpha-2 agonists. Moreover, alcohol and opioids rapidly induce dependence motivated by withdrawal avoidance. Resemblances observed in withdrawal syndromes and abuse behavior may indicate common addiction mechanisms. We argue that neurovisceral feedback influences autonomic and emotional circuits generating antireward similarly for both substances. Amygdala is central to this hypothesis as it is principally responsible for negative emotion, prominent in addiction and motivated behavior, and processes autonomic inputs while generating autonomic outputs. The solitary nucleus (NTS) has strong bidirectional connections to the amygdala and receives interoceptive inputs communicating visceral states via vagal afferents. These visceral-emotional hubs are strongly influenced by the periphery including gut microbiota. We propose that gut dysbiosis contributes to alcohol and opioid withdrawal syndromes by contributing to peripheral and neuroinflammation that stimulates these antireward pathways and motivates substance dependence.

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.

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.

3D single cell scale anatomical map of sex-dependent variability of the rat intrinsic cardiac nervous system.

We developed and analyzed a single cell scale anatomical map of the rat intrinsic cardiac nervous system (ICNS) across four male and three female hearts. We find the ICNS has a reliable structural organizational plan across individuals that provide the foundation for further analyses of the ICNS in cardiac function and disease. The distribution of the ICNS was evaluated by 3D visualization and data-driven clustering. The pattern, distribution, and clustering of ICNS neurons across all male and female rat hearts is highly conserved, demonstrating a coherent organizational plan where distinct clusters of neurons are consistently localized. Female hearts had fewer neurons, lower packing density, and slightly reduced distribution, but with identical localization. We registered the anatomical data from each heart to a geometric scaffold, normalizing their 3D coordinates for standardization of common anatomical planes and providing a path where multiple experimental results and data types can be integrated and compared.

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.

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.