A flexible ontology for inference of emergent whole cell function from relationships between subcellular processes

Essential New Complexity-Based Themes for Patient-Centered Diagnosis and Treatment of Dementia and Predementia in Older People: Multimorbidity and Multilevel Phenomenology

Dementia is a highly prevalent condition with devastating clinical and socioeconomic sequela. It is expected to triple in prevalence by 2050. No treatment is currently known to be effective. Symptomatic late-onset dementia and predementia (SLODP) affects 95% of patients with the syndrome. In contrast to trials of pharmacological prevention, no treatment is suggested to remediate or cure these symptomatic patients. SLODP but not young onset dementia is intensely associated with multimorbidity (MUM), including brain-perturbating conditions (BPCs). Recent studies showed that MUM/BPCs have a major role in the pathogenesis of SLODP. Fortunately, most MUM/BPCs are medically treatable, and thus, their treatment may modify and improve SLODP, relieving suffering and reducing its clinical and socioeconomic threats. Regrettably, the complex system features of SLODP impede the diagnosis and treatment of the potentially remediable conditions (PRCs) associated with them, mainly due to failure of pattern recognition and a flawed diagnostic workup. We suggest incorporating two SLODP-specific conceptual themes into the diagnostic workup: MUM/BPC and multilevel phenomenological themes. By doing so, we were able to improve the diagnostic accuracy of SLODP components and optimize detecting and favorably treating PRCs. These revolutionary concepts and their implications for remediability and other parameters are discussed in the paper.

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.

Biological Aging Acceleration Due to Environmental Exposures: An Exciting New Direction in Toxicogenomics Research

Biological clock technologies are designed to assess the acceleration of biological age (B-age) in diverse cell types, offering a distinctive opportunity in toxicogenomic research to explore the impact of environmental stressors, social challenges, and unhealthy lifestyles on health impairment. These clocks also play a role in identifying factors that can hinder aging and promote a healthy lifestyle. Over the past decade, researchers in epigenetics have developed testing methods that predict the chronological and biological age of organisms. These methods rely on assessing DNA methylation (DNAm) levels at specific CpG sites, RNA levels, and various biomolecules across multiple cell types, tissues, and entire organisms. Commonly known as 'biological clocks' (B-clocks), these estimators hold promise for gaining deeper insights into the pathways contributing to the development of age-related disorders. They also provide a foundation for devising biomedical or social interventions to prevent, reverse, or mitigate these disorders. This review article provides a concise overview of various epigenetic clocks and explores their susceptibility to environmental stressors.

Acute irradiation induces a senescence-like chromatin structure in mammalian oocytes

The mechanisms leading to changes in mesoscale chromatin organization during cellular aging are unknown. Here, we used transcriptional activator-like effectors, RNA-seq and superresolution analysis to determine the effects of genotoxic stress on oocyte chromatin structure. Major satellites are organized into tightly packed globular structures that coalesce into chromocenters and dynamically associate with the nucleolus. Acute irradiation significantly enhanced chromocenter mobility in transcriptionally inactive oocytes. In transcriptionally active oocytes, irradiation induced a striking unfolding of satellite chromatin fibers and enhanced the expression of transcripts required for protection from oxidative stress (Fermt1, Smg1), recovery from DNA damage (Tlk2, Rad54l) and regulation of heterochromatin assembly (Zfp296, Ski-oncogene). Non-irradiated, senescent oocytes exhibit not only high chromocenter mobility and satellite distension but also a high frequency of extra chromosomal satellite DNA. Notably, analysis of biological aging using an oocyte-specific RNA clock revealed cellular communication, posttranslational protein modifications, chromatin and histone dynamics as the top cellular processes that are dysregulated in both senescent and irradiated oocytes. Our results indicate that unfolding of heterochromatin fibers following acute genotoxic stress or cellular aging induced the formation of distended satellites and that abnormal chromatin structure together with increased chromocenter mobility leads to chromosome instability in senescent oocytes.

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.

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.

Measuring biological age using a functionally interpretable multi-tissue RNA clock

The quantification of the biological age of cells yields great promises for accelerating the discovery of novel rejuvenation strategies. Here, we present MultiTIMER, the first multi-tissue aging clock that measures the biological, rather than chronological, age of cells from their transcriptional profiles by evaluating key cellular processes. We applied MultiTIMER to more than 70,000 transcriptional profiles and demonstrate that it accurately responds to cellular stressors and known interventions while informing about dysregulated cellular functions.

Gene editing reverses arrhythmia susceptibility in humanized PLN-R14del mice: modelling a European cardiomyopathy with global impact

AIMS

A mutation in the phospholamban (PLN) gene, leading to deletion of Arg14 (R14del), has been associated with malignant arrhythmias and ventricular dilation. Identifying pre-symptomatic carriers with vulnerable myocardium is crucial because arrhythmia can result in sudden cardiac death, especially in young adults with PLN-R14del mutation. This study aimed at assessing the efficiency and efficacy of in vivo genome editing, using CRISPR/Cas9 and a cardiotropic adeno-associated virus-9 (AAV9), in improving cardiac function in young adult mice expressing the human PLN-R14del.

METHODS AND RESULTS

Humanized mice were generated expressing human wild-type (hPLN-WT) or mutant (hPLN-R14del) PLN in the heterozygous state, mimicking human carriers. Cardiac magnetic resonance imaging at 12 weeks of age showed bi-ventricular dilation and increased stroke volume in mutant vs. WT mice, with no deficit in ejection fraction or cardiac output. Challenge of ex vivo hearts with isoproterenol and rapid pacing unmasked higher propensity for sustained ventricular tachycardia (VT) in hPLN-R14del relative to hPLN-WT. Specifically, the VT threshold was significantly reduced (20.3 ± 1.2 Hz in hPLN-R14del vs. 25.7 ± 1.3 Hz in WT, P < 0.01) reflecting higher arrhythmia burden. To inactivate the R14del allele, mice were tail-vein-injected with AAV9.CRISPR/Cas9/gRNA or AAV9 empty capsid (controls). CRISPR-Cas9 efficiency was evaluated by droplet digital polymerase chain reaction and NGS-based amplicon sequencing. In vivo gene editing significantly reduced end-diastolic and stroke volumes in hPLN-R14del CRISPR-treated mice compared to controls. Susceptibility to VT was also reduced, as the VT threshold was significantly increased relative to controls (30.9 ± 2.3 Hz vs. 21.3 ± 1.5 Hz; P < 0.01).

CONCLUSIONS

This study is the first to show that disruption of hPLN-R14del allele by AAV9-CRISPR/Cas9 improves cardiac function and reduces VT susceptibility in humanized PLN-R14del mice, offering preclinical evidence for translatable approaches to therapeutically suppress the arrhythmogenic phenotype in human patients with PLN-R14del disease.

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.

A reference tissue atlas for the human kidney

Kidney Precision Medicine Project (KPMP) is building a spatially specified human kidney tissue atlas in health and disease with single-cell resolution. Here, we describe the construction of an integrated reference map of cells, pathways, and genes using unaffected regions of nephrectomy tissues and undiseased human biopsies from 56 adult subjects. We use single-cell/nucleus transcriptomics, subsegmental laser microdissection transcriptomics and proteomics, near-single-cell proteomics, 3D and CODEX imaging, and spatial metabolomics to hierarchically identify genes, pathways, and cells. Integrated data from these different technologies coherently identify cell types/subtypes within different nephron segments and the interstitium. These profiles describe cell-level functional organization of the kidney following its physiological functions and link cell subtypes to genes, proteins, metabolites, and pathways. They further show that messenger RNA levels along the nephron are congruent with the subsegmental physiological activity. This reference atlas provides a framework for the classification of kidney disease when multiple molecular mechanisms underlie convergent clinical phenotypes.

High dimensionality reduction by matrix factorization for systems pharmacology

The extraction of predictive features from the complex high-dimensional multi-omic data is necessary for decoding and overcoming the therapeutic responses in systems pharmacology. Developing computational methods to reduce high-dimensional space of features in in vitro, in vivo and clinical data is essential to discover the evolution and mechanisms of the drug responses and drug resistance. In this paper, we have utilized the matrix factorization (MF) as a modality for high dimensionality reduction in systems pharmacology. In this respect, we have proposed three novel feature selection methods using the mathematical conception of a basis for features. We have applied these techniques as well as three other MF methods to analyze eight different gene expression datasets to investigate and compare their performance for feature selection. Our results show that these methods are capable of reducing the feature spaces and find predictive features in terms of phenotype determination. The three proposed techniques outperform the other methods used and can extract a 2-gene signature predictive of a tyrosine kinase inhibitor treatment response in the Cancer Cell Line Encyclopedia.

Data-Driven Math Model of FLT3-ITD Acute Myeloid Leukemia Reveals Potential Therapeutic Targets

FLT3-mutant acute myeloid leukemia (AML) is an aggressive form of leukemia with poor prognosis. Treatment with FLT3 inhibitors frequently produces a clinical response, but the disease nevertheless often recurs. Recent studies have revealed system-wide gene expression changes in FLT3-mutant AML cell lines in response to drug treatment. Here we sought a systems-level understanding of how these cells mediate these drug-induced changes. Using RNAseq data from AML cells with an internal tandem duplication FLT3 mutation (FLT3-ITD) under six drug treatment conditions including quizartinib and dexamethasone, we identified seven distinct gene programs representing diverse biological processes involved in AML drug-induced changes. Based on the literature knowledge about genes from these modules, along with public gene regulatory network databases, we constructed a network of FLT3-ITD AML. Applying the BooleaBayes algorithm to this network and the RNAseq data, we created a probabilistic, data-driven dynamical model of acquired resistance to these drugs. Analysis of this model reveals several interventions that may disrupt targeted parts of the system-wide drug response. We anticipate co-targeting these points may result in synergistic treatments that can overcome resistance and prevent eventual recurrence.

Modelling kidney disease using ontology: insights from the Kidney Precision Medicine Project

An important need exists to better understand and stratify kidney disease according to its underlying pathophysiology in order to develop more precise and effective therapeutic agents. National collaborative efforts such as the Kidney Precision Medicine Project are working towards this goal through the collection and integration of large, disparate clinical, biological and imaging data from patients with kidney disease. Ontologies are powerful tools that facilitate these efforts by enabling researchers to organize and make sense of different data elements and the relationships between them. Ontologies are critical to support the types of big data analysis necessary for kidney precision medicine, where heterogeneous clinical, imaging and biopsy data from diverse sources must be combined to define a patient's phenotype. The development of two new ontologies - the Kidney Tissue Atlas Ontology and the Ontology of Precision Medicine and Investigation - will support the creation of the Kidney Tissue Atlas, which aims to provide a comprehensive molecular, cellular and anatomical map of the kidney. These ontologies will improve the annotation of kidney-relevant data, and eventually lead to new definitions of kidney disease in support of precision medicine.

A Systems Biology Roadmap to Decode mTOR Control System in Cancer

Mechanistic target of rapamycin (mTOR) is a critical protein in the regulation of cell fate decision making, especially in cancer cells. mTOR acts as a signal integrator and is one of the main elements of interactions among the pivotal cellular processes such as cell death, autophagy, metabolic reprogramming, cell growth, and cell cycle. The temporal control of these processes is essential for the cellular homeostasis and dysregulation of mTOR signaling pathway results in different phenotypes, including aging, oncogenesis, cell survival, cell growth, senescence, quiescence, and cell death. In this paper, we have proposed a systems biology roadmap to study mTOR control system, which introduces the theoretical and experimental modalities to decode temporal and dynamical characteristics of mTOR signaling in cancer.

Genomic signatures defining responsiveness to allopurinol and combination therapy for lung cancer identified by systems therapeutics analyses

The ability to predict responsiveness to drugs in individual patients is limited. We hypothesized that integrating molecular information from databases would yield predictions that could be experimentally tested to develop transcriptomic signatures for specific drugs. We analyzed lung adenocarcinoma patient data from The Cancer Genome Atlas and identified a subset of patients in which xanthine dehydrogenase (XDH) expression correlated with decreased survival. We tested allopurinol, an FDA‐approved drug that inhibits XDH, on human non‐small‐cell lung cancer (NSCLC) cell lines obtained from the Broad Institute Cancer Cell Line Encyclopedia and identified sensitive and resistant cell lines. We utilized the transcriptomic profiles of these cell lines to identify six‐gene signatures for allopurinol‐sensitive and allopurinol‐resistant cell lines. Transcriptomic networks identified JAK2 as an additional target in allopurinol‐resistant lines. Treatment of resistant cell lines with allopurinol and CEP‐33779 (a JAK2 inhibitor) resulted in cell death. The effectiveness of allopurinol alone or allopurinol and CEP‐33779 was verified in vivo using tumor formation in NCR‐nude mice. We utilized the six‐gene signatures to predict five additional allopurinol‐sensitive NSCLC cell lines and four allopurinol‐resistant cell lines susceptible to combination therapy. We searched the transcriptomic data from a library of patient‐derived NSCLC tumors from the Jackson Laboratory to identify tumors that would be predicted to be sensitive to allopurinol or allopurinol + CEP‐33779 treatment. Patient‐derived tumors showed the predicted drug sensitivity in vivo. These data indicate that we can use integrated molecular information from cancer databases to predict drug responsiveness in individual patients and thus enable precision medicine.

Dynamic balance between vesicle transport and microtubule growth enables neurite outgrowth

Whole cell responses involve multiple subcellular processes (SCPs). To understand how balance between SCPs controls the dynamics of whole cell responses we studied neurite outgrowth in rat primary cortical neurons in culture. We used a combination of dynamical models and experiments to understand the conditions that permitted growth at a specified velocity and when aberrant growth could lead to the formation of dystrophic bulbs. We hypothesized that dystrophic bulb formation is due to quantitative imbalances between SCPs. Simulations predict redundancies between lower level sibling SCPs within each type of high level SCP. In contrast, higher level SCPs, such as vesicle transport and exocytosis or microtubule growth characteristic of each type need to be strictly coordinated with each other and imbalances result in stalling of neurite outgrowth. From these simulations, we predicted the effect of changing the activities of SCPs involved in vesicle exocytosis or microtubule growth could lead to formation of dystrophic bulbs. siRNA ablation experiments verified these predictions. We conclude that whole cell dynamics requires balance between the higher-level SCPs involved and imbalances can terminate whole cell responses such as neurite outgrowth.

Mechanisms that cause a change of state of a cell arise from unique patterns of interactions between multiple subcellular processes (SCPs). Neurite outgrowth (NOG) is such a change of cell state where a neuron puts out a long process that eventually becomes the axon. We used a top-down based approach to mathematically model interactions between SCPs involved in NOG. These include membrane production at the cell body, membrane delivery from the cell body to the neurite tip and microtubule growth within the neurite. Our analyses show how the different SCPs interact with each other to enable NOG at a given velocity under physiological conditions. This approach is different from the commonly used bottom-up approaches that focus on predicting cell functions based on the activity of molecular interaction networks. Our simulations predict that lower-level sibling SCPs (e.g. vesicle tethering at and vesicle fusion with the plasma membrane) within a group can compensate for each other under physiological conditions, while such simple relationships do not exist between higher level SCPs (e.g. vesicle exocytosis and vesicle transport along microtubules). We predicted that imbalances of activities between higher-level SCPs induce dystrophic bulbs (a pathological response) and validated these predictions via siRNA ablation experiments.

OSCI: standardized stem cell ontology representation and use cases for stem cell investigation

Background

Stem cells and stem cell lines are widely used in biomedical research. The Cell Ontology (CL) and Cell Line Ontology (CLO) are two community-based OBO Foundry ontologies in the domains of in vivo cells and in vitro cell line cells, respectively.

Results

To support standardized stem cell investigations, we have developed an Ontology for Stem Cell Investigations (OSCI). OSCI imports stem cell and cell line terms from CL and CLO, and investigation-related terms from existing ontologies. A novel focus of OSCI is its application in representing metadata types associated with various stem cell investigations. We also applied OSCI to systematically categorize experimental variables in an induced pluripotent stem cell line cell study related to bipolar disorder. In addition, we used a semi-automated literature mining approach to identify over 200 stem cell gene markers. The relations between these genes and stem cells are modeled and represented in OSCI.

Conclusions

OSCI standardizes stem cells found in vivo and in vitro and in various stem cell investigation processes and entities. The presented use cases demonstrate the utility of OSCI in iPSC studies and literature mining related to bipolar disorder.

Systems biology primer: the basic methods and approaches

Systems biology is an integrative discipline connecting the molecular components within a single biological scale and also among different scales (e.g. cells, tissues and organ systems) to physiological functions and organismal phenotypes through quantitative reasoning, computational models and high-throughput experimental technologies. Systems biology uses a wide range of quantitative experimental and computational methodologies to decode information flow from genes, proteins and other subcellular components of signaling, regulatory and functional pathways to control cell, tissue, organ and organismal level functions. The computational methods used in systems biology provide systems-level insights to understand interactions and dynamics at various scales, within cells, tissues, organs and organisms. In recent years, the systems biology framework has enabled research in quantitative and systems pharmacology and precision medicine for complex diseases. Here, we present a brief overview of current experimental and computational methods used in systems biology.