Incident heart failure in chronic kidney disease: proteomics informs biology and risk stratification.

BACKGROUND AND AIMS: Incident heart failure (HF) among individuals with chronic kidney disease (CKD) incurs hospitalizations that burden patients and health care systems. There are few preventative therapies, and the Pooled Cohort equations to Prevent Heart Failure (PCP-HF) perform poorly in the setting of CKD. New drug targets and better risk stratification are urgently needed. METHODS: In this analysis of incident HF, SomaScan V4.0 (4638 proteins) was analysed in 2906 participants of the Chronic Renal Insufficiency Cohort (CRIC) with validation in the Atherosclerosis Risk in Communities (ARIC) study. The primary outcome was 14-year incident HF (390 events); secondary outcomes included 4-year HF (183 events), HF with reduced ejection fraction (137 events), and HF with preserved ejection fraction (165 events). Mendelian randomization and Gene Ontology were applied to examine causality and pathways. The performance of novel multi-protein risk models was compared to the PCP-HF risk score. RESULTS: Over 200 proteins were associated with incident HF after adjustment for estimated glomerular filtration rate at P < 1 × 10-5. After adjustment for covariates including N-terminal pro-B-type natriuretic peptide, 17 proteins remained associated at P < 1 × 10-5. Mendelian randomization associations were found for six proteins, of which four are druggable targets: FCG2B, IGFBP3, CAH6, and ASGR1. For the primary outcome, the C-statistic (95% confidence interval [CI]) for the 48-protein model in CRIC was 0.790 (0.735, 0.844) vs. 0.703 (0.644, 0.762) for the PCP-HF model (P = .001). C-statistic (95% CI) for the protein model in ARIC was 0.747 (0.707, 0.787). CONCLUSIONS: Large-scale proteomics reveal novel circulating protein biomarkers and potential mediators of HF in CKD. Proteomic risk models improve upon the PCP-HF risk score in this population.

Drugst.One - a plug-and-play solution for online systems medicine and network-based drug repurposing.

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.

WikiPathways 2024: next generation pathway database.

WikiPathways (wikipathways.org) is an open-source biological pathway database. Collaboration and open science are pivotal to the success of WikiPathways. Here we highlight the continuing efforts supporting WikiPathways, content growth and collaboration among pathway researchers. As an evolving database, there is a growing need for WikiPathways to address and overcome technical challenges. In this direction, WikiPathways has undergone major restructuring, enabling a renewed approach for sharing and curating pathway knowledge, thus providing stability for the future of community pathway curation. The website has been redesigned to improve and enhance user experience. This next generation of WikiPathways continues to support existing features while improving maintainability of the database and facilitating community input by providing new functionality and leveraging automation.

In vivo protein turnover rates in varying oxygen tensions nominate MYBBP1A as a mediator of the hyperoxia response.

Oxygen deprivation and excess are both toxic. Thus, the body's ability to adapt to varying oxygen tensions is critical for survival. While the hypoxia transcriptional response has been well studied, the post-translational effects of oxygen have been underexplored. In this study, we systematically investigate protein turnover rates in mouse heart, lung, and brain under different inhaled oxygen tensions. We find that the lung proteome is the most responsive to varying oxygen tensions. In particular, several extracellular matrix (ECM) proteins are stabilized in the lung under both hypoxia and hyperoxia. Furthermore, we show that complex 1 of the electron transport chain is destabilized in hyperoxia, in accordance with the exacerbation of associated disease models by hyperoxia and rescue by hypoxia. Moreover, we nominate MYBBP1A as a hyperoxia transcriptional regulator, particularly in the context of rRNA homeostasis. Overall, our study highlights the importance of varying oxygen tensions on protein turnover rates and identifies tissue-specific mediators of oxygen-dependent responses.

Using published pathway figures in enrichment analysis and machine learning.

Pathway Figure OCR (PFOCR) is a novel kind of pathway database approaching the breadth and depth of Gene Ontology while providing rich, mechanistic diagrams and direct literature support. Here, we highlight the utility of PFOCR in disease research in comparison with popular pathway databases through an assessment of disease coverage and analytical applications. In addition to common pathway analysis use cases, we present two advanced case studies demonstrating unique advantages of PFOCR in terms of cancer subtype and grade prediction analyses.

BioThings Explorer: a query engine for a federated knowledge graph of biomedical APIs.

SUMMARY: Knowledge graphs are an increasingly common data structure for representing biomedical information. These knowledge graphs can easily represent heterogeneous types of information, and many algorithms and tools exist for querying and analyzing graphs. Biomedical knowledge graphs have been used in a variety of applications, including drug repurposing, identification of drug targets, prediction of drug side effects, and clinical decision support. Typically, knowledge graphs are constructed by centralization and integration of data from multiple disparate sources. Here, we describe BioThings Explorer, an application that can query a virtual, federated knowledge graph derived from the aggregated information in a network of biomedical web services. BioThings Explorer leverages semantically precise annotations of the inputs and outputs for each resource, and automates the chaining of web service calls to execute multi-step graph queries. Because there is no large, centralized knowledge graph to maintain, BioThings Explorer is distributed as a lightweight application that dynamically retrieves information at query time. AVAILABILITY AND IMPLEMENTATION: More information can be found at https://explorer.biothings.io and code is available at https://github.com/biothings/biothings_explorer.

Drugst.One - A plug-and-play solution for online systems medicine and network-based drug repurposing.

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.

Defining blood-induced microglia functions in neurodegeneration through multiomic profiling.

Blood protein extravasation through a disrupted blood-brain barrier and innate immune activation are hallmarks of neurological diseases and emerging therapeutic targets. However, how blood proteins polarize innate immune cells remains largely unknown. Here, we established an unbiased blood-innate immunity multiomic and genetic loss-of-function pipeline to define the transcriptome and global phosphoproteome of blood-induced innate immune polarization and its role in microglia neurotoxicity. Blood induced widespread microglial transcriptional changes, including changes involving oxidative stress and neurodegenerative genes. Comparative functional multiomics showed that blood proteins induce distinct receptor-mediated transcriptional programs in microglia and macrophages, such as redox, type I interferon and lymphocyte recruitment. Deletion of the blood coagulation factor fibrinogen largely reversed blood-induced microglia neurodegenerative signatures. Genetic elimination of the fibrinogen-binding motif to CD11b in Alzheimer's disease mice reduced microglial lipid metabolism and neurodegenerative signatures that were shared with autoimmune-driven neuroinflammation in multiple sclerosis mice. Our data provide an interactive resource for investigation of the immunology of blood proteins that could support therapeutic targeting of microglia activation by immune and vascular signals.

BioThings Explorer: a query engine for a federated knowledge graph of biomedical APIs.

Knowledge graphs are an increasingly common data structure for representing biomedical information. These knowledge graphs can easily represent heterogeneous types of information, and many algorithms and tools exist for querying and analyzing graphs. Biomedical knowledge graphs have been used in a variety of applications, including drug repurposing, identification of drug targets, prediction of drug side effects, and clinical decision support. Typically, knowledge graphs are constructed by centralization and integration of data from multiple disparate sources. Here, we describe BioThings Explorer, an application that can query a virtual, federated knowledge graph derived from the aggregated information in a network of biomedical web services. BioThings Explorer leverages semantically precise annotations of the inputs and outputs for each resource, and automates the chaining of web service calls to execute multi-step graph queries. Because there is no large, centralized knowledge graph to maintain, BioThing Explorer is distributed as a lightweight application that dynamically retrieves information at query time. More information can be found at https://explorer.biothings.io, and code is available at https://github.com/biothings/biothings_explorer.

The scalable precision medicine open knowledge engine (SPOKE): a massive knowledge graph of biomedical information.

MOTIVATION: Knowledge graphs (KGs) are being adopted in industry, commerce and academia. Biomedical KG presents a challenge due to the complexity, size and heterogeneity of the underlying information. RESULTS: In this work, we present the Scalable Precision Medicine Open Knowledge Engine (SPOKE), a biomedical KG connecting millions of concepts via semantically meaningful relationships. SPOKE contains 27 million nodes of 21 different types and 53 million edges of 55 types downloaded from 41 databases. The graph is built on the framework of 11 ontologies that maintain its structure, enable mappings and facilitate navigation. SPOKE is built weekly by python scripts which download each resource, check for integrity and completeness, and then create a 'parent table' of nodes and edges. Graph queries are translated by a REST API and users can submit searches directly via an API or a graphical user interface. Conclusions/Significance: SPOKE enables the integration of seemingly disparate information to support precision medicine efforts. AVAILABILITY AND IMPLEMENTATION: The SPOKE neighborhood explorer is available at https://spoke.rbvi.ucsf.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

scNetViz: from single cells to networks using Cytoscape.

Single-cell RNA-sequencing (scRNA-seq) has revolutionized molecular biology and medicine by enabling high-throughput studies of cellular heterogeneity in diverse tissues. Applying network biology approaches to scRNA-seq data can provide useful insights into genes driving heterogeneous cell-type compositions of tissues. Here, we present scNetViz-  a Cytoscape app to aid biological interpretation of cell clusters in scRNA-seq data using network analysis. scNetViz calculates the differential expression of each gene across clusters and then creates a cluster-specific gene functional interaction network between the significantly differentially expressed genes for further analysis, such as pathway enrichment analysis. To automate a complete data analysis workflow, scNetViz integrates parts of the Scanpy software, which is a popular Python package for scRNA-seq data analysis, with Cytoscape apps such as stringApp, cyPlot, and enhancedGraphics. We describe our implementation of methods for accessing data from public single cell atlas projects, differential expression analysis, visualization, and automation. scNetViz enables users to analyze data from public atlases or their own experiments, which we illustrate with two use cases. Analysis can be performed via the Cytoscape GUI or CyREST programming interface using R (RCy3) or Python (py4cytoscape).

SARS-CoV-2 spike protein induces abnormal inflammatory blood clots neutralized by fibrin immunotherapy.

Blood clots are a central feature of coronavirus disease-2019 (COVID-19) and can culminate in pulmonary embolism, stroke, and sudden death. However, it is not known how abnormal blood clots form in COVID-19 or why they occur even in asymptomatic and convalescent patients. Here we report that the Spike protein from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) binds to the blood coagulation factor fibrinogen and induces structurally abnormal blood clots with heightened proinflammatory activity. SARS-CoV-2 Spike virions enhanced fibrin-mediated microglia activation and induced fibrinogen-dependent lung pathology. COVID-19 patients had fibrin autoantibodies that persisted long after acute infection. Monoclonal antibody 5B8, targeting the cryptic inflammatory fibrin epitope, inhibited thromboinflammation. Our results reveal a procoagulant role for the SARS-CoV-2 Spike and propose fibrin-targeting interventions as a treatment for thromboinflammation in COVID-19. ONE-SENTENCE SUMMARY: SARS-CoV-2 spike induces structurally abnormal blood clots and thromboinflammation neutralized by a fibrin-targeting antibody.

Proteogenomic insights into the biology and treatment of HPV-negative head and neck squamous cell carcinoma.

We present a proteogenomic study of 108 human papilloma virus (HPV)-negative head and neck squamous cell carcinomas (HNSCCs). Proteomic analysis systematically catalogs HNSCC-associated proteins and phosphosites, prioritizes copy number drivers, and highlights an oncogenic role for RNA processing genes. Proteomic investigation of mutual exclusivity between FAT1 truncating mutations and 11q13.3 amplifications reveals dysregulated actin dynamics as a common functional consequence. Phosphoproteomics characterizes two modes of EGFR activation, suggesting a new strategy to stratify HNSCCs based on EGFR ligand abundance for effective treatment with inhibitory EGFR monoclonal antibodies. Widespread deletion of immune modulatory genes accounts for low immune infiltration in immune-cold tumors, whereas concordant upregulation of multiple immune checkpoint proteins may underlie resistance to anti-programmed cell death protein 1 monotherapy in immune-hot tumors. Multi-omic analysis identifies three molecular subtypes with high potential for treatment with CDK inhibitors, anti-EGFR antibody therapy, and immunotherapy, respectively. Altogether, proteogenomics provides a systematic framework to inform HNSCC biology and treatment.

Introducing R as a smart version of calculators enables beginners to explore it on their own.

Rapid technological advances in the past decades have enabled molecular biologists to generate large-scale and complex data with affordable resource investments, or obtain such data from public repositories. Yet, many graduate students, postdoctoral scholars, and senior researchers in the biosciences find themselves ill-equipped to analyze large-scale data. Global surveys have revealed that active researchers prefer short training workshops to fill their skill gaps. In this article, we focus on the challenge of delivering a short data analysis workshop to absolute beginners in computer programming. We propose that introducing R or other programming languages for data analysis as smart versions of calculators can help lower the communication barrier with absolute beginners. We describe this comparison with a few analogies and hope that other instructors will find them useful. We utilized these in our four-hour long training workshops involving participatory live coding, which we delivered in person and via videoconferencing. Anecdotal evidence suggests that our exposition made R programming seem easy and enabled beginners to explore it on their own.

Mir-30d Regulates Cardiac Remodeling by Intracellular and Paracrine Signaling.

RATIONALE: Previous translational studies implicate plasma extracellular microRNA-30d (miR-30d) as a biomarker in left ventricular remodeling and clinical outcome in heart failure (HF) patients, although precise mechanisms remain obscure. OBJECTIVE: To investigate the mechanism of miR-30d-mediated cardioprotection in HF. METHODS AND RESULTS: In rat and mouse models of ischemic HF, we show that miR-30d gain of function (genetic, lentivirus, or agomiR-mediated) improves cardiac function, decreases myocardial fibrosis, and attenuates cardiomyocyte (CM) apoptosis. Genetic or locked nucleic acid-based knock-down of miR-30d expression potentiates pathological left ventricular remodeling, with increased dysfunction, fibrosis, and cardiomyocyte death. RNA sequencing of in vitro miR-30d gain and loss of function, together with bioinformatic prediction and experimental validation in cardiac myocytes and fibroblasts, were used to identify and validate direct targets of miR-30d. miR-30d expression is selectively enriched in cardiomyocytes, induced by hypoxic stress and is acutely protective, targeting MAP4K4 (mitogen-associate protein kinase 4) to ameliorate apoptosis. Moreover, miR-30d is secreted primarily in extracellular vesicles by cardiomyocytes and inhibits fibroblast proliferation and activation by directly targeting integrin α5 in the acute phase via paracrine signaling to cardiac fibroblasts. In the chronic phase of ischemic remodeling, lower expression of miR-30d in the heart and plasma extracellular vesicles is associated with adverse remodeling in rodent models and human subjects and is linked to whole-blood expression of genes implicated in fibrosis and inflammation, consistent with observations in model systems. CONCLUSIONS: These findings provide the mechanistic underpinning for the cardioprotective association of miR-30d in human HF. More broadly, our findings support an emerging paradigm involving intercellular communication of extracellular vesicle-contained miRNAs (microRNAs) to transregulate distinct signaling pathways across cell types. Functionally validated RNA biomarkers and their signaling networks may warrant further investigation as novel therapeutic targets in HF.

The Dynamic Platelet Transcriptome in Obesity and Weight Loss.

OBJECTIVE: Adiposity is associated with oxidative stress, inflammation, and glucose intolerance. Previous data suggest that platelet gene expression is associated with key cardiometabolic phenotypes, including body mass index but stable in healthy individuals over time. However, modulation of gene expression in platelets in response to metabolic shifts (eg, weight reduction) is unknown and may be important to defining mechanism. Approach and Results: Platelet RNA sequencing and aggregation were performed from 21 individuals with massive weight loss (>45 kg) following bariatric surgery. Based on RNA sequencing data, we measured the expression of 67 genes from isolated platelet RNA using high-throughput quantitative reverse transcription quantitative PCR in 1864 FHS (Framingham Heart Study) participants. Many transcripts not previously studied in platelets were differentially expressed with bariatric surgical weight loss, appeared specific to platelets (eg, not differentially expressed in leukocytes), and were enriched for a nonalcoholic fatty liver disease pathway. Platelet aggregation studies did not detect alteration in platelet function after significant weight loss. Linear regression models demonstrated several platelet genes modestly associated with cross-sectional cardiometabolic phenotypes, including body mass index. There were no associations between studied transcripts and incident diabetes or cardiovascular end points. CONCLUSIONS: In summary, while there is no change in platelet aggregation function after significant weight loss, the human platelet experiences a dramatic transcriptional shift that implicates pathways potentially relevant to improved cardiometabolic risk postweight loss (eg, nonalcoholic fatty liver disease). Further studies are needed to determine the mechanistic importance of these observations.

WikiPathways: connecting communities.

WikiPathways (https://www.wikipathways.org) is a biological pathway database known for its collaborative nature and open science approaches. With the core idea of the scientific community developing and curating biological knowledge in pathway models, WikiPathways lowers all barriers for accessing and using its content. Increasingly more content creators, initiatives, projects and tools have started using WikiPathways. Central in this growth and increased use of WikiPathways are the various communities that focus on particular subsets of molecular pathways such as for rare diseases and lipid metabolism. Knowledge from published pathway figures helps prioritize pathway development, using optical character and named entity recognition. We show the growth of WikiPathways over the last three years, highlight the new communities and collaborations of pathway authors and curators, and describe various technologies to connect to external resources and initiatives. The road toward a sustainable, community-driven pathway database goes through integration with other resources such as Wikidata and allowing more use, curation and redistribution of WikiPathways content.

Pathway information extracted from 25 years of pathway figures.

Thousands of pathway diagrams are published each year as static figures inaccessible to computational queries and analyses. Using a combination of machine learning, optical character recognition, and manual curation, we identified 64,643 pathway figures published between 1995 and 2019 and extracted 1,112,551 instances of human genes, comprising 13,464 unique NCBI genes, participating in a wide variety of biological processes. This collection represents an order of magnitude more genes than found in the text of the same papers, and thousands of genes missing from other pathway databases, thus presenting new opportunities for discovery and research.

Shared Mechanisms Govern HIV Transcriptional Suppression in Circulating CD103+ and Gut CD4+ T Cells.

Latent HIV infection is the main barrier to cure, and most HIV-infected cells reside in the gut, where distinct but unknown mechanisms may promote viral latency. Transforming growth factor ß (TGF-ß), which induces the expression of CD103 on tissue-resident memory T cells, has been implicated in HIV latency. Using CD103 as a surrogate marker to identify cells that have undergone TGF-ß signaling, we compared the HIV RNA/DNA contents and cellular transcriptomes of CD103+ and CD103- CD4 T cells from the blood and rectum of HIV-negative (HIV-) and antiretroviral therapy (ART)-suppressed HIV-positive (HIV+) individuals. Like gut CD4+ T cells, circulating CD103+ cells harbored more HIV DNA than did CD103- cells but transcribed less HIV RNA per provirus. Circulating CD103+ cells also shared a gene expression profile that is closer to that of gut CD4 T cells than to that of circulating CD103- cells, with significantly lower expression levels of ribosomal proteins and transcriptional and translational pathways associated with HIV expression but higher expression levels of a subset of genes implicated in suppressing HIV transcription. These findings suggest that blood CD103+ CD4 T cells can serve as a model to study the molecular mechanisms of HIV latency in the gut and reveal new cellular factors that may contribute to HIV latency.IMPORTANCE The ability of HIV to establish a reversibly silent, "latent" infection is widely regarded as the main barrier to curing HIV. Most HIV-infected cells reside in tissues such as the gut, but it is unclear what mechanisms maintain HIV latency in the blood or gut. We found that circulating CD103+ CD4+ T cells are enriched for HIV-infected cells in a latent-like state. Using RNA sequencing (RNA-seq), we found that CD103+ T cells share a cellular transcriptome that more closely resembles that of CD4+ T cells from the gut, suggesting that they are homing to or from the gut. We also identified the cellular genes whose expression distinguishes gut CD4+ or circulating CD103+ T cells from circulating CD103- T cells, including some genes that have been implicated in HIV expression. These genes may contribute to latent HIV infection in the gut and may serve as new targets for therapies aimed at curing HIV.