Multiscale mapping of transcriptomic signatures for cardiotoxic drugs.

Drug-induced gene expression profiles can identify potential mechanisms of toxicity. We focus on obtaining signatures for cardiotoxicity of FDA-approved tyrosine kinase inhibitors (TKIs) in human induced-pluripotent-stem-cell-derived cardiomyocytes, using bulk transcriptomic profiles. We use singular value decomposition to identify drug-selective patterns across cell lines obtained from multiple healthy human subjects. Cellular pathways affected by cardiotoxic TKIs include energy metabolism, contractile, and extracellular matrix dynamics. Projecting these pathways to published single cell expression profiles indicates that TKI responses can be evoked in both cardiomyocytes and fibroblasts. Integration of transcriptomic outlier analysis with whole genomic sequencing of our six cell lines enables us to correctly reidentify a genomic variant causally linked to anthracycline-induced cardiotoxicity and predict genomic variants potentially associated with TKI-induced cardiotoxicity. We conclude that mRNA expression profiles when integrated with publicly available genomic, pathway, and single cell transcriptomic datasets, provide multiscale signatures for cardiotoxicity that could be used for drug development and patient stratification.

Endogenous adenine is a potential driver of the cardiovascular-kidney-metabolic syndrome.

Mechanisms underlying the cardiovascular-kidney-metabolic (CKM) syndrome are unknown, although key small molecule metabolites may be involved. Bulk and spatial metabolomics identified adenine to be upregulated and specifically enriched in coronary blood vessels in hearts from patients with diabetes and left ventricular hypertrophy. Single nucleus gene expression studies revealed that endothelial methylthioadenosine phosphorylase (MTAP) was increased in human hearts with hypertrophic cardiomyopathy. The urine adenine/creatinine ratio in patients was predictive of incident heart failure with preserved ejection fraction. Heart adenine and MTAP gene expression was increased in a 2-hit mouse model of hypertrophic heart disease and in a model of diastolic dysfunction with diabetes. Inhibition of MTAP blocked adenine accumulation in the heart, restored heart dysfunction in mice with type 2 diabetes and prevented ischemic heart damage in a rat model of myocardial infarction. Mechanistically, adenine-induced impaired mitophagy was reversed by reduction of mTOR. These studies indicate that endogenous adenine is in a causal pathway for heart failure and ischemic heart disease in the context of CKM syndrome.

Pharmacology of boldine: summary of the field and update on recent advances.

Over the past decade, boldine, a naturally occurring alkaloid found in several plant species including the Chilean Boldo tree, has garnered attention for its efficacy in rodent models of human disease. Some of the properties that have been attributed to boldine include antioxidant activities, neuroprotective and analgesic actions, hepatoprotective effects, anti-inflammatory actions, cardioprotective effects and anticancer potential. Compelling data now indicates that boldine blocks connexin (Cx) hemichannels (HCs) and that many if not all of its effects in rodent models of injury and disease are due to CxHC blockade. Here we provide an overview of boldine's pharmacological properties, including its efficacy in rodent models of common human injuries and diseases, and of its absorption, distribution, pharmacokinetics, and metabolism.

From transcriptomics to digital twins of organ function.

Cell level functions underlie tissue and organ physiology. Gene expression patterns offer extensive views of the pathways and processes within and between cells. Single cell transcriptomics provides detailed information on gene expression within cells, cell types, subtypes and their relative proportions in organs. Functional pathways can be scalably connected to physiological functions at the cell and organ levels. Integrating experimentally obtained gene expression patterns with prior knowledge of pathway interactions enables identification of networks underlying whole cell functions such as growth, contractility, and secretion. These pathways can be computationally modeled using differential equations to simulate cell and organ physiological dynamics regulated by gene expression changes. Such computational systems can be thought of as parts of digital twins of organs. Digital twins, at the core, need computational models that represent in detail and simulate how dynamics of pathways and networks give rise to whole cell level physiological functions. Integration of transcriptomic responses and numerical simulations could simulate and predict whole cell functional outputs from transcriptomic data. We developed a computational pipeline that integrates gene expression timelines and systems of coupled differential equations to generate cell-type selective dynamical models. We tested our integrative algorithm on the eicosanoid biosynthesis network in macrophages. Converting transcriptomic changes to a dynamical model allowed us to predict dynamics of prostaglandin and thromboxane synthesis and secretion by macrophages that matched published lipidomics data obtained in the same experiments. Integration of cell-level system biology simulations with genomic and clinical data using a knowledge graph framework will allow us to create explicit predictive models that mechanistically link genomic determinants to organ function. Such integration requires a multi-domain ontological framework to connect genomic determinants to gene expression and cell pathways and functions to organ level phenotypes in healthy and diseased states. These integrated scalable models of tissues and organs as accurate digital twins predict health and disease states for precision medicine.

Ontology-based modeling, integration, and analysis of heterogeneous clinical, pathological, and molecular kidney data for precision medicine.

Many data resources generate, process, store, or provide kidney related molecular, pathological, and clinical data. Reference ontologies offer an opportunity to support knowledge and data integration. The Kidney Precision Medicine Project (KPMP) team contributed to the representation and addition of 329 kidney phenotype terms to the Human Phenotype Ontology (HPO), and identified many subcategories of acute kidney injury (AKI) or chronic kidney disease (CKD). The Kidney Tissue Atlas Ontology (KTAO) imports and integrates kidney-related terms from existing ontologies (e.g., HPO, CL, and Uberon) and represents 259 kidney-related biomarkers. We have also developed a precision medicine metadata ontology (PMMO) to integrate 50 variables from KPMP and CZ CellxGene data resources and applied PMMO for integrative kidney data analysis. The gene expression profiles of kidney gene biomarkers were specifically analyzed under healthy control or AKI/CKD disease states. This work demonstrates how ontology-based approaches support multi-domain data and knowledge integration in precision medicine.

Memory discrimination is promoted by the expression of the transcription repressor WT1 in the dentate gyrus.

The hippocampus is critical for the precise formation of contextual memories. Overlapping inputs coming from the entorhinal cortex are processed by the trisynaptic pathway to form distinct memories. Disruption in any step of the circuit flow can lead to a lack of memory precision, and to memory interference. We have identified the transcriptional repressor Wilm's Tumor 1 (WT1) as an important regulator of synaptic plasticity involved in memory discrimination in the hippocampus. In male mice, using viral and transgenic approaches, we showed that WT1 deletion in granule cells of the dentate gyrus (DG) disrupts memory discrimination. With electrophysiological methods, we then identified changes in granule cells' excitability and DG synaptic transmission indicating that WT1 knockdown in DG granule cells disrupts the inhibitory feedforward input from mossy fibers to CA3 by decreasing mIPSCs and shifting the normal excitatory/inhibitory (E/I) balance in the DG → CA3 circuit in favor of excitation. Finally, using a chemogenetic approach, we established a causal link between granule cell hyperexcitability and memory discrimination impairments. Our results suggest that WT1 enables a circuit-level computation that drives pattern discrimination behavior.

A spatially specified systems pharmacology therapy for axonal recovery after injury.

There are no known drugs or drug combinations that promote substantial central nervous system axonal regeneration after injury. We used systems pharmacology approaches to model pathways underlying axonal growth and identify a four-drug combination that regulates multiple subcellular processes in the cell body and axons using the optic nerve crush model in rats. We intravitreally injected agonists HU-210 (cannabinoid receptor-1) and IL-6 (interleukin 6 receptor) to stimulate retinal ganglion cells for axonal growth. We applied, in gel foam at the site of nerve injury, Taxol to stabilize growing microtubules, and activated protein C to clear the debris field since computational models predicted that this drug combination regulating two subcellular processes at the growth cone produces synergistic growth. Physiologically, drug treatment restored or preserved pattern electroretinograms and some of the animals had detectable visual evoked potentials in the brain and behavioral optokinetic responses. Morphology experiments show that the four-drug combination protects axons or promotes axonal regrowth to the optic chiasm and beyond. We conclude that spatially targeted drug treatment is therapeutically relevant and can restore limited functional recovery.

Spinal cord injury regulates circular RNA expression in axons.

Introduction: Neurons transport mRNA and translational machinery to axons for local translation. After spinal cord injury (SCI), de novo translation is assumed to enable neurorepair. Knowledge of the identity of axonal mRNAs that participate in neurorepair after SCI is limited. We sought to identify and understand how axonal RNAs play a role in axonal regeneration. Methods: We obtained preparations enriched in axonal mRNAs from control and SCI rats by digesting spinal cord tissue with cold-active protease (CAP). The digested samples were then centrifuged to obtain a supernatant that was used to identify mRNA expression. We identified differentially expressed genes (DEGS) after SCI and mapped them to various biological processes. We validated the DEGs by RT-qPCR and RNA-scope. Results: The supernatant fraction was highly enriched for mRNA from axons. Using Gene Ontology, the second most significant pathway for all DEGs was axonogenesis. Among the DEGs was Rims2, which is predominately a circular RNA (circRNA) in the CNS. We show that Rims2 RNA within spinal cord axons is circular. We found an additional 200 putative circRNAs in the axonal-enriched fraction. Knockdown in primary rat cortical neurons of the RNA editing enzyme ADAR1, which inhibits formation of circRNAs, significantly increased axonal outgrowth and increased the expression of circRims2. Using Rims2 as a prototype we used Circular RNA Interactome to predict miRNAs that bind to circRims2 also bind to the 3'UTR of GAP-43, PTEN or CREB1, all known regulators of axonal outgrowth. Axonally-translated GAP-43 supports axonal elongation and we detect GAP-43 mRNA in the rat axons by RNAscope. Discussion: By enriching for axonal RNA, we detect SCI induced DEGs, including circRNA such as Rims2. Ablation of ADAR1, the enzyme that regulates circRNA formation, promotes axonal outgrowth of cortical neurons. We developed a pathway model using Circular RNA Interactome that indicates that Rims2 through miRNAs can regulate the axonal translation GAP-43 to regulate axonal regeneration. We conclude that axonal regulatory pathways will play a role in neurorepair.

Endogenous adenine mediates kidney injury in diabetic models and predicts diabetic kidney disease in patients.

Diabetic kidney disease (DKD) can lead to end-stage kidney disease (ESKD) and mortality; however, few mechanistic biomarkers are available for high-risk patients, especially those without macroalbuminuria. Urine from participants with diabetes from the Chronic Renal Insufficiency Cohort (CRIC) study, the Singapore Study of Macro-angiopathy and Micro-vascular Reactivity in Type 2 Diabetes (SMART2D), and the American Indian Study determined whether urine adenine/creatinine ratio (UAdCR) could be a mechanistic biomarker for ESKD. ESKD and mortality were associated with the highest UAdCR tertile in the CRIC study and SMART2D. ESKD was associated with the highest UAdCR tertile in patients without macroalbuminuria in the CRIC study, SMART2D, and the American Indian study. Empagliflozin lowered UAdCR in nonmacroalbuminuric participants. Spatial metabolomics localized adenine to kidney pathology, and single-cell transcriptomics identified ribonucleoprotein biogenesis as a top pathway in proximal tubules of patients without macroalbuminuria, implicating mTOR. Adenine stimulated matrix in tubular cells via mTOR and stimulated mTOR in mouse kidneys. A specific inhibitor of adenine production was found to reduce kidney hypertrophy and kidney injury in diabetic mice. We propose that endogenous adenine may be a causative factor in DKD.

Muscle-restricted Nox4 knockout partially corrects muscle contractility following spinal cord injury in mice.

Spinal cord injury (SCI) results in severe atrophy of skeletal muscle in paralyzed regions, and a decrease in the force generated by muscle per unit of cross-sectional area. Oxidation of skeletal muscle ryanodine 1 receptors (RyR1) reduces contractile force due to reduced binding of calstabin 1 to RyR1 together with altered gating of RyR1. One cause of RyR1 oxidation is NADPH oxidase 4 (Nox4). We have previously shown that in rats, RyR1 was oxidized and bound less calstabin 1 at 56 days after spinal cord injury (SCI) by transection. Here, we used a conditional knock-out mouse model of Nox4 in muscle to investigate the role of Nox4 in reduced muscle specific force after SCI. Peak twitch force in control mice after SCI was reduced by 42% compared to sham-operated controls but was increased by approximately 43% in SCI Nox4 conditional KO mice compared to SCI controls although it remained less than that for sham-operated controls. Unlike what observed in rats, after SCI the expression of Nox4 was not increased in gastrocnemius muscle and binding of calstabin 1 to RyR1 was not reduced in this muscle. The results suggest a link between Nox4 expression in muscle tissue and reduction in muscle twitch force, however further studies are needed to understand the mechanistic basis for this linkage.

Synaptojanin 1 Modulates Functional Recovery After Incomplete Spinal Cord Injury in Male Apolipoprotein E Epsilon 4 Mice.

Apolipoprotein E epsilon 4 (ApoE4) is the second most common variant of ApoE, being present in ∼14% of the population. Clinical reports identify ApoE4 as a genetic risk factor for poor outcomes after traumatic spinal cord injury (SCI) and spinal cord diseases such as cervical myelopathy. To date, there is no intervention to promote recovery of function after SCI/spinal cord diseases that is specifically targeted at ApoE4-associated impairment. Studies in the human and mouse brain link ApoE4 to elevated levels of synaptojanin 1 (synj1), a lipid phosphatase that degrades phosphoinositol 4,5-bisphosphate (PIP2) into inositol 4-monophosphate. Synj1 regulates rearrangements of the cytoskeleton as well as endocytosis and trafficking of synaptic vesicles. We report here that, as compared to ApoE3 mice, levels of synj1 messenger RNA and protein were elevated in spinal cords of healthy ApoE4 mice associated with lower PIP2 levels. Using a moderate-severity model of contusion SCI in mice, we found that genetic reduction of synj1 improved locomotor function recovery at 14 days after SCI in ApoE4 mice without altering spared white matter. Genetic reduction of synj1 did not alter locomotor recovery of ApoE3 mice after SCI. Bulk RNA sequencing revealed that at 14 days after SCI in ApoE4 mice, genetic reduction of synj1 upregulated genes involved in glutaminergic synaptic transmission just above and below the lesion. Overall, our findings provide evidence for a link between synj1 to poor outcomes after SCI in ApoE4 mice, up to 14 days post-injury, through mechanisms that may involve the function of excitatory glutaminergic neurons.

Machine learning using institution-specific multi-modal electronic health records improves mortality risk prediction for cardiac surgery patients.

Background: The Society of Thoracic Surgeons risk scores are widely used to assess risk of morbidity and mortality in specific cardiac surgeries but may not perform optimally in all patients. In a cohort of patients undergoing cardiac surgery, we developed a data-driven, institution-specific machine learning-based model inferred from multi-modal electronic health records and compared the performance with the Society of Thoracic Surgeons models. Methods: All adult patients undergoing cardiac surgery between 2011 and 2016 were included. Routine electronic health record administrative, demographic, clinical, hemodynamic, laboratory, pharmacological, and procedural data features were extracted. The outcome was postoperative mortality. The database was randomly split into training (development) and test (evaluation) cohorts. Models developed using 4 classification algorithms were compared using 6 evaluation metrics. The performance of the final model was compared with the Society of Thoracic Surgeons models for 7 index surgical procedures. Results: A total of 6392 patients were included and described by 4016 features. Overall mortality was 3.0% (n = 193). The XGBoost algorithm using only features with no missing data (336 features) yielded the best-performing predictor. When applied to the test set, the predictor performed well (F-measure = 0.775; precision = 0.756; recall = 0.795; accuracy = 0.986; area under the receiver operating characteristic curve = 0.978; area under the precision-recall curve = 0.804). eXtreme Gradient Boosting consistently demonstrated improved performance over the Society of Thoracic Surgeons models when evaluated on index procedures within the test set. Conclusions: Machine learning models using institution-specific multi-modal electronic health records may improve performance in predicting mortality for individual patients undergoing cardiac surgery compared with the standard-of-care, population-derived Society of Thoracic Surgeons models. Institution-specific models may provide insights complementary to population-derived risk predictions to aid patient-level decision making.

Boldine modulates glial transcription and functional recovery in a murine model of contusion spinal cord injury.

Membrane channels such as those formed by connexins (Cx) and P2X7 receptors (P2X7R) are permeable to calcium ions and other small molecules such as adenosine triphosphate (ATP) and glutamate. Release of ATP and glutamate through these channels is a key mechanism driving tissue response to traumas such as spinal cord injury (SCI). Boldine, an alkaloid isolated from the Chilean boldo tree, blocks both Cx and Panx1 hemichannels (HCs). To test if boldine could improve function after SCI, boldine or vehicle was administered to treat mice with a moderate severity contusion-induced SCI. Boldine led to greater spared white matter and increased locomotor function as determined by the Basso Mouse Scale and horizontal ladder rung walk tests. Boldine treatment reduced immunostaining for markers of activated microglia (Iba1) and astrocytic (GFAP) markers while increasing that for axon growth and neuroplasticity (GAP-43). Cell culture studies demonstrated that boldine blocked glial HC, specifically Cx26 and Cx30, in cultured astrocytes and blocked calcium entry through activated P2X7R. RT-qPCR studies showed that boldine treatment reduced expression of the chemokine Ccl2, cytokine IL-6 and microglial gene CD68, while increasing expression of the neurotransmission genes Snap25 and Grin2b, and Gap-43. Bulk RNA sequencing revealed that boldine modulated a large number of genes involved in neurotransmission in spinal cord tissue just caudal from the lesion epicenter at 14 days after SCI. Numbers of genes regulated by boldine was much lower at 28 days after injury. These results indicate that boldine treatment ameliorates injury and spares tissue to increase locomotor function.

Role of endogenous adenine in kidney failure and mortality with diabetes.

Diabetic kidney disease (DKD) can lead to end-stage kidney disease (ESKD) and mortality, however, few mechanistic biomarkers are available for high risk patients, especially those without macroalbuminuria. Urine from participants with diabetes from Chronic Renal Insufficiency Cohort (CRIC), Singapore Study of Macro-Angiopathy and Reactivity in Type 2 Diabetes (SMART2D), and the Pima Indian Study determined if urine adenine/creatinine ratio (UAdCR) could be a mechanistic biomarker for ESKD. ESKD and mortality were associated with the highest UAdCR tertile in CRIC (HR 1.57, 1.18, 2.10) and SMART2D (HR 1.77, 1.00, 3.12). ESKD was associated with the highest UAdCR tertile in patients without macroalbuminuria in CRIC (HR 2.36, 1.26, 4.39), SMART2D (HR 2.39, 1.08, 5.29), and Pima Indian study (HR 4.57, CI 1.37-13.34). Empagliflozin lowered UAdCR in non-macroalbuminuric participants. Spatial metabolomics localized adenine to kidney pathology and transcriptomics identified ribonucleoprotein biogenesis as a top pathway in proximal tubules of patients without macroalbuminuria, implicating mammalian target of rapamycin (mTOR). Adenine stimulated matrix in tubular cells via mTOR and stimulated mTOR in mouse kidneys. A specific inhibitor of adenine production was found to reduce kidney hypertrophy and kidney injury in diabetic mice. We propose that endogenous adenine may be a causative factor in DKD.

Predicting individual-specific cardiotoxicity responses induced by tyrosine kinase inhibitors.

Introduction: Tyrosine kinase inhibitor drugs (TKIs) are highly effective cancer drugs, yet many TKIs are associated with various forms of cardiotoxicity. The mechanisms underlying these drug-induced adverse events remain poorly understood. We studied mechanisms of TKI-induced cardiotoxicity by integrating several complementary approaches, including comprehensive transcriptomics, mechanistic mathematical modeling, and physiological assays in cultured human cardiac myocytes. Methods: Induced pluripotent stem cells (iPSCs) from two healthy donors were differentiated into cardiac myocytes (iPSC-CMs), and cells were treated with a panel of 26 FDA-approved TKIs. Drug-induced changes in gene expression were quantified using mRNA-seq, changes in gene expression were integrated into a mechanistic mathematical model of electrophysiology and contraction, and simulation results were used to predict physiological outcomes. Results: Experimental recordings of action potentials, intracellular calcium, and contraction in iPSC-CMs demonstrated that modeling predictions were accurate, with 81% of modeling predictions across the two cell lines confirmed experimentally. Surprisingly, simulations of how TKI-treated iPSC-CMs would respond to an additional arrhythmogenic insult, namely, hypokalemia, predicted dramatic differences between cell lines in how drugs affected arrhythmia susceptibility, and these predictions were confirmed experimentally. Computational analysis revealed that differences between cell lines in the upregulation or downregulation of particular ion channels could explain how TKI-treated cells responded differently to hypokalemia. Discussion: Overall, the study identifies transcriptional mechanisms underlying cardiotoxicity caused by TKIs, and illustrates a novel approach for integrating transcriptomics with mechanistic mathematical models to generate experimentally testable, individual-specific predictions of adverse event risk.

Dissecting aortic aneurysm in Marfan syndrome is associated with losartan-sensitive transcriptomic modulation of aortic cells.

To improve our limited understanding of the pathogenesis of thoracic aortic aneurysm (TAA) that leads to acute aortic dissection, single-cell RNA sequencing (scRNA-seq) was employed to profile disease-relevant transcriptomic changes of aortic cell populations in a well-characterized mouse model of the most commonly diagnosed form of Marfan syndrome (MFS). As result, 2 discrete subpopulations of aortic cells (SMC3 and EC4) were identified only in the aorta of Fbn1mgR/mgR mice. SMC3 cells highly express genes related to extracellular matrix formation and nitric oxide signaling, whereas the EC4 transcriptional profile is enriched in smooth muscle cell (SMC), fibroblast, and immune cell-related genes. Trajectory analysis predicted close phenotypic modulation between SMC3 and EC4, which were therefore analyzed together as a discrete MFS-modulated (MFSmod) subpopulation. In situ hybridization of diagnostic transcripts located MFSmod cells at the intima of Fbn1mgR/mgR aortas. Reference-based data set integration revealed transcriptomic similarity between MFSmod- and SMC-derived cell clusters modulated in human TAA. Consistent with the angiotensin II type I receptor (At1r) contribution to TAA development, MFSmod cells were absent in the aorta of Fbn1mgR/mgR mice treated with the At1r antagonist losartan. Altogether, our findings indicate that a discrete dynamic alteration of aortic cell identity is associated with dissecting TAA in MFS mice and increased risk of aortic dissection in MFS patients.

FN (Fibronectin)-Integrin α5 Signaling Promotes Thoracic Aortic Aneurysm in a Mouse Model of Marfan Syndrome.

BACKGROUND: Marfan syndrome, caused by mutations in the gene for fibrillin-1, leads to thoracic aortic aneurysms (TAAs). Phenotypic modulation of vascular smooth muscle cells (SMCs) and ECM (extracellular matrix) remodeling are characteristic of both nonsyndromic and Marfan aneurysms. The ECM protein FN (fibronectin) is elevated in the tunica media of TAAs and amplifies inflammatory signaling in endothelial and SMCs through its main receptor, integrin α5ß1. We investigated the role of integrin α5-specific signals in Marfan mice in which the cytoplasmic domain of integrin α5 was replaced with that of integrin α2 (denoted α5/2 chimera). METHODS: We crossed α5/2 chimeric mice with Fbn1mgR/mgR mice (mgR model of Marfan syndrome) to evaluate the survival rate and pathogenesis of TAAs among wild-type, α5/2, mgR, and α5/2 mgR mice. Further biochemical and microscopic analysis of porcine and mouse aortic SMCs investigated molecular mechanisms by which FN affects SMCs and subsequent development of TAAs. RESULTS: FN was elevated in the thoracic aortas from Marfan patients, in nonsyndromic aneurysms, and in mgR mice. The α5/2 mutation greatly prolonged survival of Marfan mice, with improved elastic fiber integrity, mechanical properties, SMC density, and SMC contractile gene expression. Furthermore, plating of wild-type SMCs on FN decreased contractile gene expression and activated inflammatory pathways whereas α5/2 SMCs were resistant. These effects correlated with increased NF-kB activation in cultured SMCs and mgR aortas, which was alleviated by the α5/2 mutation or NF-kB inhibition. CONCLUSIONS: FN-integrin α5 signaling is a significant driver of TAA in the mgR mouse model. This pathway thus warrants further investigation as a therapeutic target.

Boldine modulates glial transcription and functional recovery in a murine model of contusion spinal cord injury.

Membrane channels such as connexins (Cx), pannexins (Panx) and P2X 7 receptors (P2X 7 R) are permeable to calcium ions and other small molecules such as ATP and glutamate. Release of ATP and glutamate through these channels is a key mechanism driving tissue response to traumas such as spinal cord injury (SCI). Boldine, an alkaloid isolated from the Chilean boldo tree, blocks both Cx hemichannels (HC) and Panx. To test if boldine could improve function after SCI, boldine or vehicle was administered to treat mice with a moderate severity contusion-induced SCI. Boldine led to greater spared white matter and increased locomotor function as determined by the Basso Mouse Scale and horizontal ladder rung walk tests. Boldine treatment reduced immunostaining for markers of activated microglia (Iba1) and astrocytic (GFAP) markers while increasing that for axon growth and neuroplasticity (GAP-43). Cell culture studies demonstrated that boldine blocked glial HC, specifically Cx26 and Cx30, in cultured astrocytes and blocked calcium entry through activated P2X 7 R. RT-qPCR studies showed that boldine treatment reduced expression of the chemokine Ccl2, cytokine IL-6 and microglial gene CD68, while increasing expression of the neurotransmission genes Snap25 and Grin2b, and Gap-43. Bulk RNA sequencing (of the spinal cord revealed that boldine modulated a large number of genes involved in neurotransmission in in spinal cord tissue just below the lesion epicenter at 14 days after SCI. Numbers of genes regulated by boldine was much lower at 28 days after injury. These results indicate that boldine treatment ameliorates injury and spares tissue to increase locomotor function.

Relating individual cell division events to single-cell ERK and Akt activity time courses.

Biochemical correlates of stochastic single-cell fates have been elusive, even for the well-studied mammalian cell cycle. We monitored single-cell dynamics of the ERK and Akt pathways, critical cell cycle progression hubs and anti-cancer drug targets, and paired them to division events in the same single cells using the non-transformed MCF10A epithelial line. Following growth factor treatment, in cells that divide both ERK and Akt activities are significantly higher within the S-G2 time window (~ 8.5-40 h). Such differences were much smaller in the pre-S-phase, restriction point window which is traditionally associated with ERK and Akt activity dependence, suggesting unappreciated roles for ERK and Akt in S through G2. Simple metrics of central tendency in this time window are associated with subsequent cell division fates. ERK activity was more strongly associated with division fates than Akt activity, suggesting Akt activity dynamics may contribute less to the decision driving cell division in this context. We also find that ERK and Akt activities are less correlated with each other in cells that divide. Network reconstruction experiments demonstrated that this correlation behavior was likely not due to crosstalk, as ERK and Akt do not interact in this context, in contrast to other transformed cell types. Overall, our findings support roles for ERK and Akt activity throughout the cell cycle as opposed to just before the restriction point, and suggest ERK activity dynamics may be more important than Akt activity dynamics for driving cell division in this non-transformed context.

Whole cell response to receptor stimulation involves many deep and distributed subcellular biochemical processes.

Neurite outgrowth is an integrated whole cell response triggered by the cannabinoid-1 receptor. We sought to identify the many different biochemical pathways that contribute to this whole cell response. To understand underlying mechanisms, we identified subcellular processes (SCPs) composed of one or more biochemical pathways and their interactions required for this response. Differentially expressed genes and proteins were obtained from bulk transcriptomics and proteomic analysis of extracts from cells stimulated with a cannabinoid-1 receptor agonist. We used these differentially expressed genes and proteins to build networks of interacting SCPs by combining the expression data with prior pathway knowledge. From these SCP networks, we identified additional genes that when ablated, experimentally validated the SCP involvement in neurite outgrowth. Our experiments and informatics modeling allowed us to identify diverse SCPs such as those involved in pyrimidine metabolism, lipid biosynthesis, and mRNA splicing and stability, along with more predictable SCPs such as membrane vesicle transport and microtubule dynamics. We find that SCPs required for neurite outgrowth are widely distributed among many biochemical pathways required for constitutive cellular functions, several of which are termed 'deep', since they are distal to signaling pathways and the key SCPs directly involved in extension of the neurite. In contrast, 'proximal' SCPs are involved in microtubule growth and membrane vesicle transport dynamics required for neurite outgrowth. From these bioinformatics and dynamical models based on experimental data, we conclude that receptor-mediated regulation of subcellular functions for neurite outgrowth is both distributed, that is, involves many different biochemical pathways, and deep.