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.

Acyl-CoA thioesterase-2 facilitates β-oxidation in glycolytic skeletal muscle in a lipid supply dependent manner

Acyl-Coenzyme A (acyl-CoA) thioesters are compartmentalized intermediates that participate in in multiple metabolic reactions within the mitochondrial matrix. The limited availability of free CoA (CoASH) in the matrix raises the question of how the local acyl-CoA concentration is regulated to prevent trapping of CoASH from overload of any specific substrate. Acyl-CoA thioesterase-2 (ACOT2) hydrolyzes long-chain acyl-CoAs to their constituent fatty acids and CoASH, and is the only mitochondrial matrix ACOT refractory to inhibition by CoASH. Thus, we reasoned that ACOT2 may constitutively regulate matrix acyl-CoA levels. Acot2 deletion in murine skeletal muscle (SM) resulted in acyl-CoA build-up when lipid supply and energy demands were modest. When energy demand and pyruvate availability were elevated, lack of ACOT2 activity promoted glucose oxidation. This preference for glucose over fatty acid oxidation was recapitulated in C2C12 myotubes with acute depletion of Acot2 , and overt inhibition of β-oxidation was demonstrated in isolated mitochondria from Acot2 -depleted glycolytic SM. In mice fed a high fat diet, ACOT2 enabled the accretion of acyl-CoAs and ceramide derivatives in glycolytic SM, and this was associated with worse glucose homeostasis compared to when ACOT2 was absent. These observations suggest that ACOT2 supports CoASH availability to facilitate β-oxidation in glycolytic SM when lipid supply is modest. However, when lipid supply is high, ACOT2 enables acyl-CoA and lipid accumulation, CoASH sequestration, and poor glucose homeostasis. Thus, ACOT2 regulates matrix acyl-CoA concentration in glycolytic muscle, and its impact depends on lipid supply.

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.

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.

Conditional expression of endorepellin in the tumor vasculature attenuates breast cancer growth, angiogenesis and hyaluronan deposition

The tumor stroma of most solid malignancies is characterized by a pathological accumulation of pro-angiogenic and pro-tumorigenic hyaluronan driving tumorigenesis and metastatic potential. Of all three hyaluronan synthase isoforms, HAS2 is the primary enzyme that promotes the build-up of tumorigenic HA in breast cancer. Previously, we discovered that endorepellin, the angiostatic C-terminal fragment of perlecan, evokes a catabolic mechanism targeting endothelial HAS2 and hyaluronan via autophagic induction. To explore the translational implications of endorepellin in breast cancer, we created a double transgenic, inducible Tie2CreERT2;endorepellin(ER)Ki mouse line that expresses recombinant endorepellin specifically from the endothelium. We investigated the therapeutic effects of recombinant endorepellin overexpression in an orthotopic, syngeneic breast cancer allograft mouse model. First, adenoviral delivery of Cre evoking intratumor expression of endorepellin in ERKi mice suppressed breast cancer growth, peritumor hyaluronan and angiogenesis. Moreover, tamoxifen-induced expression of recombinant endorepellin specifically from the endothelium in Tie2CreERT2;ERKi mice markedly suppressed breast cancer allograft growth, hyaluronan deposition in the tumor proper and perivascular tissues, and tumor angiogenesis. These results provide insight into the tumor suppressing activity of endorepellin at the molecular level and implicate endorepellin as a promising cancer protein therapy that targets hyaluronan in the tumor microenvironment.

Longitudinal ultrasound imaging and network modeling in rats reveal sex-dependent suppression of liver regeneration after resection in alcoholic liver disease

Liver resection is an important surgical technique in the treatment of cancers and transplantation. We used ultrasound imaging to study the dynamics of liver regeneration following two-thirds partial hepatectomy (PHx) in male and female rats fed via Lieber-deCarli liquid diet protocol of ethanol or isocaloric control or chow for 5–7 weeks. Ethanol-fed male rats did not recover liver volume to the pre-surgery levels over the course of 2 weeks after surgery. By contrast, ethanol-fed female rats as well as controls of both sexes showed normal volume recovery. Contrary to expectations, transient increases in both portal and hepatic artery blood flow rates were seen in most animals, with ethanol-fed males showing higher peak portal flow than any other experimental group. A computational model of liver regeneration was used to evaluate the contribution of physiological stimuli and estimate the animal-specific parameter intervals. The results implicate lower metabolic load, over a wide range of cell death sensitivity, in matching the model simulations to experimental data of ethanol-fed male rats. However, in the ethanol-fed female rats and controls of both sexes, metabolic load was higher and in combination with cell death sensitivity matched the observed volume recovery dynamics. We conclude that adaptation to chronic ethanol intake has a sex-dependent impact on liver volume recovery following liver resection, likely mediated by differences in the physiological stimuli or cell death responses that govern the regeneration process. Immunohistochemical analysis of pre- and post-resection liver tissue validated the results of computational modeling by associating lack of sensitivity to cell death with lower rates of cell death in ethanol-fed male rats. Our results illustrate the potential for non-invasive ultrasound imaging to assess liver volume recovery towards supporting development of clinically relevant computational models of liver regeneration.

Elucidating the Mechanisms of Dynamic and Robust Control of the Liver Homeostatic Renewal Process: Cell Network Modeling and Analysis

Recent experimental investigations of liver homeostatic renewal have identified high replication capacity hepatocyte populations as the primary maintainers of liver mass. However, the molecular and cellular processes controlling liver homeostatic renewal remain unknown. To address this problem, we developed and analyzed a mathematical model describing cellular network interactions underlying liver homeostatic renewal. Model simulation results demonstrate that without feedback control, basic homeostatic renewal is not robust to disruptions, leading to tissue loss under persistent/repetitive insults. Consequently, we extended our basic model to incorporate putative regulatory interactions and investigated how such interactions may confer robustness on the homeostatic renewal process. We utilized a Design of Experiments approach to identify the combination of feedback interactions that yields a cell network model of homeostatic renewal capable of maintaining liver mass robustly during persistent/repetitive injury. Simulations of this robust model indicate that repeated injury destabilizes liver homeostasis within several months, which differs from epidemiological observations of a much slower decay of liver function occurring over several years. To address this discrepancy, we extended the model to include feedback control by liver nonparenchymal cells. Simulations and analysis of the final multicellular feedback control network suggest that achieving robust liver homeostatic renewal requires intrinsic stability in a hepatocellular network combined with feedback control by nonparenchymal cells.

Genome-Scale Metabolic Modeling Reveals Sequential Dysregulation of Glutathione Metabolism in Livers from Patients with Alcoholic Hepatitis

Alcoholic hepatitis (AH) is the most severe form of alcoholic liver disease for which there is no efficacious treatment aiding most patients. AH manifests differently in individuals, with some patients showing debilitating symptoms more so than others. Previous studies showed significant metabolic dysregulation associated with AH. Therefore, we sought to analyze how the activity of metabolic pathways differed in the liver of patients with varying degrees of AH severity. We utilized a genome-scale metabolic modeling approach that allowed for integration of a generic human cellular metabolic model with specific RNA-seq data corresponding to healthy and multiple liver disease states to predict the metabolic fluxes within each disease state. Additionally, we performed a systems-level analysis of the transcriptomic data and predicted metabolic flux data to identify the regulatory and functional differences in liver metabolism with increasing severity of AH. Our results provide unique insights into the sequential dysregulation of the solute transport mechanisms underlying the glutathione metabolic pathway with increasing AH disease severity. We propose targeting of the solute transporters in the glutathione pathway to mimic the flux activity of the healthy liver state as a potential therapeutic intervention for AH.

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.

A Spatial Model of Hepatic Calcium Signaling and Glucose Metabolism Under Autonomic Control Reveals Functional Consequences of Varying Liver Innervation Patterns Across Species

Rapid breakdown of hepatic glycogen stores into glucose plays an important role during intense physical exercise to maintain systemic euglycemia. Hepatic glycogenolysis is governed by several different liver-intrinsic and systemic factors such as hepatic zonation, circulating catecholamines, hepatocellular calcium signaling, hepatic neuroanatomy, and the central nervous system (CNS). Of the factors regulating hepatic glycogenolysis, the extent of lobular innervation varies significantly between humans and rodents. While rodents display very few autonomic nerve terminals in the liver, nearly every hepatic layer in the human liver receives neural input. In the present study, we developed a multi-scale, multi-organ model of hepatic metabolism incorporating liver zonation, lobular scale calcium signaling, hepatic innervation, and direct and peripheral organ-mediated communication between the liver and the CNS. We evaluated the effect of each of these governing factors on the total hepatic glucose output and zonal glycogenolytic patterns within liver lobules during simulated physical exercise. Our simulations revealed that direct neuronal stimulation of the liver and an increase in circulating catecholamines increases hepatic glucose output mediated by mobilization of intracellular calcium stores and lobular scale calcium waves. Comparing simulated glycogenolysis between human-like and rodent-like hepatic innervation patterns (extensive vs. minimal) suggested that propagation of calcium transients across liver lobules acts as a compensatory mechanism to improve hepatic glucose output in sparsely innervated livers. Interestingly, our simulations suggested that catecholamine-driven glycogenolysis is reduced under portal hypertension. However, increased innervation coupled with strong intercellular communication can improve the total hepatic glucose output under portal hypertension. In summary, our modeling and simulation study reveals a complex interplay of intercellular and multi-organ interactions that can lead to differing calcium dynamics and spatial distributions of glycogenolysis at the lobular scale in the liver.

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.

Dysregulation of miR-21-associated miRNA regulatory networks by chronic ethanol consumption impairs liver regeneration

Impaired liver regeneration has been considered as a hallmark of progression of alcohol-associated liver disease. Our previous studies demonstrated that in vivo inhibition of the microRNA (miRNA) miR21 can restore regenerative capacity of the liver in chronic ethanol-fed animals. The present study focuses on the role of microRNA regulatory networks that are likely to mediate the miR-21 action. Rats were chronically fed an ethanol-enriched diet along with pair-fed control animals and treated with AM21 (anti-miR-21), a locked nucleic acid antisense to miR-21. Partial hepatectomy (PHx) was performed and miRNA expression profiling over the course of liver regeneration was assessed. Our results showed dynamic expression changes in several miRNAs after PHx, notably with altered miRNA expression profiles between ethanol and control groups. We found that in vivo inhibition of miR-21 led to correlated differential expression of miR-340-5p and anticorrelated expression of miR-365, let-7a, miR-1224, and miR-146a across all sample groups after PHx. Gene set enrichment analysis identified a miRNA signature significantly associated with hepatic stellate cell activation within whole liver tissue data. We hypothesized that at least part of the PHx-induced miRNA network changes responsive to miR-21 inhibition is localized to hepatic stellate cells. We validated this hypothesis using AM21 and TGF-β treatments in LX-2 human hepatic stellate cells in culture and measured expression levels of select miRNAs by quantitative RT-PCR. Based on the in vivo and in vitro results, we propose a hepatic stellate cell miRNA regulatory network as contributing to the restoration of liver regenerative capacity by miR-21 inhibition.

Retinoic Acid Fluctuation Activates an Uneven, Direction-Dependent Network-Wide Robustness Response in Early Embryogenesis

Robustness is a feature of regulatory pathways to ensure signal consistency in light of environmental changes or genetic polymorphisms. The retinoic acid (RA) pathway, is a central developmental and tissue homeostasis regulatory signal, strongly dependent on nutritional sources of retinoids and affected by environmental chemicals. This pathway is characterized by multiple proteins or enzymes capable of performing each step and their integration into a self-regulating network. We studied RA network robustness by transient physiological RA signaling disturbances followed by kinetic transcriptomic analysis of the recovery during embryogenesis. The RA metabolic network was identified as the main regulated module to achieve signaling robustness using an unbiased pattern analysis. We describe the network-wide responses to RA signal manipulation and found the feedback autoregulation to be sensitive to the direction of the RA perturbation: RA knockdown exhibited an upper response limit, whereas RA addition had a minimal feedback-activation threshold. Surprisingly, our robustness response analysis suggests that the RA metabolic network regulation exhibits a multi-objective optimization, known as Pareto optimization, characterized by trade-offs between competing functionalities. We observe that efficient robustness to increasing RA is accompanied by worsening robustness to reduced RA levels and vice versa. This direction-dependent trade-off in the network-wide feedback response, results in an uneven robustness capacity of the RA network during early embryogenesis, likely a significant contributor to the manifestation of developmental defects.

Analyzing the relationships between city opioid deaths and socioeconomic factors

BACKGROUND

The opioid crisis has been increasing in severity over the past several decades. Every year, thousands of people die due to opioid-related issues. The factors that play into the causes of this are complex. Our study aimed to see how commute times, city budget for roads, and city budget for arts and culture contribute to the city's overall opioid death rates. We know that overdose rates are more common in cities rather than rural areas, therefore making one major city per state as part of our design.

METHODS

We collected data from one city per state (n = 50) and ran two ANOVA tests and 11 logistical regression tests. Both types of tests were run on IBM SPSS Statistical software version 25 at its default settings with a confidence interval set to 95 percent. Opioid deaths were the dependent variable, whereas commute time, budget for roads, and budget for arts and culture were the independent variables.

RESULTS

Commute time yielded a significant result in almost all the tests it was included in: Table 1, 0.033; Table 3, 0.000; Figure 4, 0.000; Figure 5, 0.000. Budget for roads also showed significant results in most of its tests as well: Table 1, 0.003; Table 2, 0.047; Figure 3, 0.001. Budget for arts and culture showed significant results but not in a pattern that we could interpret: Table 1, 0.002; Table 2, 0.021; Table 4, 0.013.

CONCLUSIONS

Commute time and the budget for roads are likely to play a role in their city's opioid crisis. Understanding where a city fits in relation to these results may better help them prepare and reduce opioid death rates.

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.

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.

Innervation and Neuronal Control of the Mammalian Sinoatrial Node a Comprehensive Atlas

[Figure: see text].

Integrated Multiomics Reveals Glucose Use Reprogramming and Identifies a Novel Hexokinase in Alcoholic Hepatitis

BACKGROUND & AIMS

We recently showed that alcoholic hepatitis (AH) is characterized by dedifferentiation of hepatocytes and loss of mature functions. Glucose metabolism is tightly regulated in healthy hepatocytes. We hypothesize that AH may lead to metabolic reprogramming of the liver, including dysregulation of glucose metabolism.

METHODS

We performed integrated metabolomic and transcriptomic analyses of liver tissue from patients with AH or alcoholic cirrhosis or normal liver tissue from hepatic resection. Focused analyses of chromatin immunoprecipitation coupled to DNA sequencing was performed. Functional in vitro studies were performed in primary rat and human hepatocytes and HepG2 cells.

RESULTS

Patients with AH exhibited specific changes in the levels of intermediates of glycolysis/gluconeogenesis, the tricarboxylic acid cycle, and monosaccharide and disaccharide metabolism. Integrated analysis of the transcriptome and metabolome showed the used of alternate energetic pathways, metabolite sinks and bottlenecks, and dysregulated glucose storage in patients with AH. Among genes involved in glucose metabolism, hexokinase domain containing 1 (HKDC1) was identified as the most up-regulated kinase in patients with AH. Histone active promoter and enhancer markers were increased in the HKDC1 genomic region. High HKDC1 levels were associated with the development of acute kidney injury and decreased survival. Increased HKDC1 activity contributed to the accumulation of glucose-6-P and glycogen in primary rat hepatocytes.

CONCLUSIONS

Altered metabolite levels and messenger RNA expression of metabolic enzymes suggest the existence of extensive reprogramming of glucose metabolism in AH. Increased HKDC1 expression may contribute to dysregulated glucose metabolism and represents a novel biomarker and therapeutic target for AH.

Credible practice of modeling and simulation in healthcare: ten rules from a multidisciplinary perspective

The complexities of modern biomedicine are rapidly increasing. Thus, modeling and simulation have become increasingly important as a strategy to understand and predict the trajectory of pathophysiology, disease genesis, and disease spread in support of clinical and policy decisions. In such cases, inappropriate or ill-placed trust in the model and simulation outcomes may result in negative outcomes, and hence illustrate the need to formalize the execution and communication of modeling and simulation practices. Although verification and validation have been generally accepted as significant components of a model’s credibility, they cannot be assumed to equate to a holistic credible practice, which includes activities that can impact comprehension and in-depth examination inherent in the development and reuse of the models. For the past several years, the Committee on Credible Practice of Modeling and Simulation in Healthcare, an interdisciplinary group seeded from a U.S. interagency initiative, has worked to codify best practices. Here, we provide Ten Rules for credible practice of modeling and simulation in healthcare developed from a comparative analysis by the Committee’s multidisciplinary membership, followed by a large stakeholder community survey. These rules establish a unified conceptual framework for modeling and simulation design, implementation, evaluation, dissemination and usage across the modeling and simulation life-cycle. While biomedical science and clinical care domains have somewhat different requirements and expectations for credible practice, our study converged on rules that would be useful across a broad swath of model types. In brief, the rules are: (1) Define context clearly. (2) Use contextually appropriate data. (3) Evaluate within context. (4) List limitations explicitly. (5) Use version control. (6) Document appropriately. (7) Disseminate broadly. (8) Get independent reviews. (9) Test competing implementations. (10) Conform to standards. Although some of these are common sense guidelines, we have found that many are often missed or misconstrued, even by seasoned practitioners. Computational models are already widely used in basic science to generate new biomedical knowledge. As they penetrate clinical care and healthcare policy, contributing to personalized and precision medicine, clinical safety will require established guidelines for the credible practice of modeling and simulation in healthcare.

A Comprehensive Integrated Anatomical and Molecular Atlas of Rat Intrinsic Cardiac Nervous System

We have developed and integrated several technologies including whole-organ imaging and software development to support an initial precise 3D neuroanatomical mapping and molecular phenotyping of the intracardiac nervous system (ICN). While qualitative and gross anatomical descriptions of the anatomy of the ICN have each been pursued, we here bring forth a comprehensive atlas of the entire rat ICN at single-cell resolution. Our work precisely integrates anatomical and molecular data in the 3D digitally reconstructed whole heart with resolution at the micron scale. We now display the full extent and the position of neuronal clusters on the base and posterior left atrium of the rat heart, and the distribution of molecular phenotypes that are defined along the base-to-apex axis, which had not been previously described. The development of these approaches needed for this work has produced method pipelines that provide the means for mapping other organs.

•

Comprehensive single-neuron-scale mapping of the intrinsic cardiac nervous system

•

Whole-organ high-throughput imaging and reconstruction at a cellular resolution

•

3D anatomical framework for spatially tracked single-neuron molecular phenotypes

•

Integrated histology, neuron mapping, and molecular profiles for 3D organ reconstruction

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.

Diurnal Patterns of Gene Expression in the Dorsal Vagal Complex and the Central Nucleus of the Amygdala - Non-rhythm-generating Brain Regions

Genes that establish the circadian clock have differential expression with respect to solar time in central and peripheral tissues. Here, we find circadian-time-induced differential expression in a large number of genes not associated with circadian rhythms in two brain regions lacking overt circadian function: the dorsal vagal complex (DVC) and the central nucleus of the amygdala (CeA). These regions primarily engage in autonomic, homeostatic, and emotional regulation. However, we find striking diurnal shifts in gene expression in these regions of male Sprague Dawley rats with no obvious patterns that could be attributed to function or region. These findings have implications for the design of gene expression studies as well as for the potential effects of xenobiotics on these regions that regulate autonomic and emotional states.

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.

Inflammation-associated suppression of metabolic gene networks in acute and chronic liver disease

Inflammation has been recognized as essential for restorative regeneration. Here, we analyzed the sequential processes during onset of liver injury and subsequent regeneration based on time-resolved transcriptional regulatory networks (TRNs) to understand the relationship between inflammation, mature organ function, and regeneration. Genome-wide expression and TRN analysis were performed time dependently in mouse liver after acute injury by CCl₄ (2 h, 8 h, 1, 2, 4, 6, 8, 16 days), as well as lipopolysaccharide (LPS, 24 h) and compared to publicly available data after tunicamycin exposure (mouse, 6 h), hepatocellular carcinoma (HCC, mouse), and human chronic liver disease (non-alcoholic fatty liver, HBV infection and HCC). Spatiotemporal investigation differentiated lobular zones for signaling and transcription factor expression. Acute CCl₄ intoxication induced expression of gene clusters enriched for inflammation and stress signaling that peaked between 2 and 24 h, accompanied by a decrease of mature liver functions, particularly metabolic genes. Metabolism decreased not only in pericentral hepatocytes that underwent CCl₄-induced necrosis, but extended to the surviving periportal hepatocytes. Proliferation and tissue restorative TRNs occurred only later reaching a maximum at 48 h. The same upstream regulators (e.g. inhibited RXR function) were implicated in increased inflammation and suppressed metabolism. The concomitant inflammation/metabolism TRN occurred similarly after acute LPS and tunicamycin challenges, in chronic mouse models and also in human liver diseases. Downregulation of metabolic genes occurs concomitantly to induce inflammation-associated genes as an early response and appears to be initiated by similar upstream regulators in acute and chronic liver diseases in humans and mice. In the acute setting, proliferation and restorative regeneration associated TRNs peak only later when metabolism is already suppressed.

A Novel Mouse Model of Acute-on-Chronic Cholestatic Alcoholic Liver Disease: A Systems Biology Comparison With Human Alcoholic Hepatitis

BACKGROUND

Alcohol-related liver disease is the main cause of liver-related mortality worldwide. The development of novel targeted therapies for patients with advanced forms (i.e., alcoholic hepatitis, AH) is hampered by the lack of suitable animal models. Here, we developed a novel mouse model of acute-on-chronic alcohol liver injury with cholestasis and fibrosis and performed an extensive molecular comparative analysis with human AH.

METHODS

For the mouse model of acute-on-chronic liver injury, we used 3,5-diethoxycarbonyl-1,4-dihydrocollidine (DDC, 0.05% w/w) diet for 8 weeks to establish cholestatic liver fibrosis. After 1-week washout period, male mice were fed intragastrically for 4 weeks with up to 24 g/kg of ethyl alcohol in a high-fat diet. This animal model was phenotyped using histopathology, clinical chemistry, microbiome, and gene expression approaches. Data were compared to the phenotypes of human alcohol-related liver disease, including AH.

RESULTS

Mice with cholestatic liver fibrosis and subsequent alcohol exposure (DDC + EtOH) exhibited exacerbated liver fibrosis with a pericellular pattern, increased neutrophil infiltration, and ductular proliferation, all characteristics of human AH. DDC administration had no effect on urine alcohol concentration or liver steatosis. Importantly, DDC- and alcohol-treated mice showed a transcriptomic signature that resembled that of patients with AH. Finally, we show that mice in the DDC + EtOH group had an increased gut barrier dysfunction, mimicking an important pathophysiological mechanism of human AH.

CONCLUSIONS

We developed a novel mouse model of acute-on-chronic cholestatic alcoholic liver injury that has considerable translational potential and can be used to test novel therapeutic modalities for AH.

Single-Cell Glia and Neuron Gene Expression in the Central Amygdala in Opioid Withdrawal Suggests Inflammation With Correlated Gut Dysbiosis

Drug-seeking in opioid dependence is due in part to the severe negative emotion associated with the withdrawal syndrome. It is well-established that negative emotional states emerge from activity in the amygdala. More recently, gut microflora have been shown to contribute substantially to such emotions. We measured gene expression in single glia and neurons gathered from the amygdala using laser capture microdissection and simultaneously measured gut microflora in morphine-dependent and withdrawn rats to investigate drivers of negative emotion in opioid withdrawal. We found that neuroinflammatory genes, notably Tnf, were upregulated in the withdrawal condition and that astrocytes, in particular, were highly active. We also observe a decreased Firmicutes to Bacteroides ratio in opioid withdrawal indicating gut dysbiosis. We speculate that these inflammatory and gut microflora changes contribute to the negative emotion experienced in opioid withdrawal that motivates dependence.

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.