Deciphering nutritional stress responses via knowledge-enriched transcriptomics for microbial engineering.

Understanding diverse bacterial nutritional requirements and responses is foundational in microbial research and biotechnology. In this study, we employed knowledge-enriched transcriptomic analytics to decipher complex stress responses of Vibrio natriegens to supplied nutrients, aiming to enhance microbial engineering efforts. We computed 64 independently modulated gene sets that comprise a quantitative basis for transcriptome dynamics across a comprehensive transcriptomics dataset containing a broad array of nutrient conditions. Our approach led to the i) identification of novel transporter systems for diverse substrates, ii) a detailed understanding of how trace elements affect metabolism and growth, and iii) extensive characterization of nutrient-induced stress responses, including osmotic stress, low glycolytic flux, proteostasis, and altered protein expression. By clarifying the relationship between the acetate-associated regulon and glycolytic flux status of various nutrients, we have showcased its vital role in directing optimal carbon source selection. Our findings offer deep insights into the transcriptional landscape of bacterial nutrition and underscore its significance in tailoring strain engineering strategies, thereby facilitating the development of more efficient and robust microbial systems for biotechnological applications.

Agent-based model demonstrates the impact of nonlinear, complex interactions between cytokinces on muscle regeneration.

Muscle regeneration is a complex process due to dynamic and multiscale biochemical and cellular interactions, making it difficult to identify microenvironmental conditions that are beneficial to muscle recovery from injury using experimental approaches alone. To understand the degree to which individual cellular behaviors impact endogenous mechanisms of muscle recovery, we developed an agent-based model (ABM) using the Cellular-Potts framework to simulate the dynamic microenvironment of a cross-section of murine skeletal muscle tissue. We referenced more than 100 published studies to define over 100 parameters and rules that dictate the behavior of muscle fibers, satellite stem cells (SSCs), fibroblasts, neutrophils, macrophages, microvessels, and lymphatic vessels, as well as their interactions with each other and the microenvironment. We utilized parameter density estimation to calibrate the model to temporal biological datasets describing cross-sectional area (CSA) recovery, SSC, and fibroblast cell counts at multiple timepoints following injury. The calibrated model was validated by comparison of other model outputs (macrophage, neutrophil, and capillaries counts) to experimental observations. Predictions for eight model perturbations that varied cell or cytokine input conditions were compared to published experimental studies to validate model predictive capabilities. We used Latin hypercube sampling and partial rank correlation coefficient to identify in silico perturbations of cytokine diffusion coefficients and decay rates to enhance CSA recovery. This analysis suggests that combined alterations of specific cytokine decay and diffusion parameters result in greater fibroblast and SSC proliferation compared to individual perturbations with a 13% increase in CSA recovery compared to unaltered regeneration at 28 days. These results enable guided development of therapeutic strategies that similarly alter muscle physiology (i.e. converting extracellular matrix [ECM]-bound cytokines into freely diffusible forms as studied in cancer therapeutics or delivery of exogenous cytokines) during regeneration to enhance muscle recovery after injury.

PERCEPT: Replacing binary p-value thresholding with scaling for more nuanced identification of sample differences.

Key to a biologists' capacity to understand data is the ability to make meaningful conclusions about differences in experimental observations. Typically, data are noisy, and conventional methods rely on replicates to average out noise and enable univariate statistical tests to assign p-values. Yet thresholding p-values to determine significance is controversial and often misleading, especially for omics datasets with few replicates. This study introduces PERCEPT, an alternative that transforms data using an ad-hoc scaling factor derived from p-values. By applying this method, low confidence effects are suppressed compared to high confidence ones, enabling clearer patterns to emerge from noisy datasets. The effectiveness of PERCEPT scaling is demonstrated using simulated datasets and published omics studies. The approach reduces the exclusion of datapoints, enhances accuracy, and enables nuanced interpretation of data. PERCEPT is easy to apply for the non-expert in statistics and provides researchers a straightforward way to improve data-driven analyses.

Statistical and computational methods for integrating microbiome, host genomics, and metabolomics data.

Large-scale microbiome studies are progressively utilizing multiomics designs, which include the collection of microbiome samples together with host genomics and metabolomics data. Despite the increasing number of data sources, there remains a bottleneck in understanding the relationships between different data modalities due to the limited number of statistical and computational methods for analyzing such data. Furthermore, little is known about the portability of general methods to the metagenomic setting and few specialized techniques have been developed. In this review, we summarize and implement some of the commonly used methods. We apply these methods to real data sets where shotgun metagenomic sequencing and metabolomics data are available for microbiome multiomics data integration analysis. We compare results across methods, highlight strengths and limitations of each, and discuss areas where statistical and computational innovation is needed.

Integrated systems biology analysis of acute lymphoblastic leukemia: unveiling molecular signatures and drug repurposing opportunities.

Acute lymphoblastic leukemia (ALL) is a hematological malignancy characterized by aberrant proliferation and accumulation of lymphoid precursor cells within the bone marrow. The tyrosine kinase inhibitor (TKI), imatinib mesylate, has played a significant role in the treatment of Philadelphia chromosome-positive ALL (Ph + ALL). However, the achievement of durable and sustained therapeutic success remains a challenge due to the development of TKI resistance during the clinical course.The primary objective of this investigation is to propose a novel and efficacious treatment approach through drug repositioning, targeting ALL and its Ph + subtype by identifying and addressing differentially expressed genes (DEGs). This study involves a comprehensive analysis of transcriptome datasets pertaining to ALL and Ph + ALL in order to identify DEGs associated with the progression of these diseases to identify possible repurposable drugs that target identified hub proteins.The outcomes of this research have unveiled 698 disease-related DEGs for ALL and 100 for Ph + ALL. Furthermore, a subset of drugs, specifically glipizide for Ph + ALL, and maytansine and isoprenaline for ALL, have been identified as potential candidates for therapeutic intervention. Subsequently, cytotoxicity assessments were performed to confirm the in vitro cytotoxic effects of these selected drugs on both ALL and Ph + ALL cell lines.In conclusion, this study offers a promising avenue for the management of ALL and Ph + ALL through drug repurposed drugs. Further investigations are necessary to elucidate the mechanisms underlying cell death, and clinical trials are recommended to validate the promising results obtained through drug repositioning strategies.

Global epistasis and the emergence of function in microbial consortia.

The many functions of microbial communities emerge from a complex web of interactions between organisms and their environment. This poses a significant obstacle to engineering microbial consortia, hindering our ability to harness the potential of microorganisms for biotechnological applications. In this study, we demonstrate that the collective effect of ecological interactions between microbes in a community can be captured by simple statistical models that predict how adding a new species to a community will affect its function. These predictive models mirror the patterns of global epistasis reported in genetics, and they can be quantitatively interpreted in terms of pairwise interactions between community members. Our results illuminate an unexplored path to quantitatively predicting the function of microbial consortia from their composition, paving the way to optimizing desirable community properties and bringing the tasks of predicting biological function at the genetic, organismal, and ecological scales under the same quantitative formalism.

Pan-cancer analysis of the interplay between mutational signatures and cellular signaling.

Cancer is a multi-faceted disease with intricate relationships between mutagenic processes, alterations in cellular signaling, and the tissue microenvironment. To date, these processes have been largely studied in isolation. A systematic understanding of how they interact and influence each other is lacking. Here, we present a framework for systematically characterizing the interaction between pairs of mutational signatures and between signatures and signaling pathway alterations. We applied this framework to large-scale data from TCGA and PCAWG and identified multiple positive and negative interactions, both cross֊tissue and tissue֊specific, that provide new insights into the molecular routes observed in tumorigenesis and their respective drivers. This framework allows for a more fine-grained dissection of common and distinct etiology of mutational signatures. We further identified several interactions with both positive and negative impacts on patient survival, demonstrating their clinical relevance and potential for improving personalized cancer care.

Spotless, a reproducible pipeline for benchmarking cell type deconvolution in spatial transcriptomics.

Spatial transcriptomics (ST) technologies allow the profiling of the transcriptome of cells while keeping their spatial context. Since most commercial untargeted ST technologies do not yet operate at single-cell resolution, computational methods such as deconvolution are often used to infer the cell type composition of each sequenced spot. We benchmarked 11 deconvolution methods using 63 silver standards, 3 gold standards, and 2 case studies on liver and melanoma tissues. We developed a simulation engine called synthspot to generate silver standards from single-cell RNA-sequencing data, while gold standards are generated by pooling single cells from targeted ST data. We evaluated methods based on their performance, stability across different reference datasets, and scalability. We found that cell2location and RCTD are the top-performing methods, but surprisingly, a simple regression model outperforms almost half of the dedicated spatial deconvolution methods. Furthermore, we observe that the performance of all methods significantly decreased in datasets with highly abundant or rare cell types. Our results are reproducible in a Nextflow pipeline, which also allows users to generate synthetic data, run deconvolution methods and optionally benchmark them on their dataset (https://github.com/saeyslab/spotless-benchmark).

A unified framework for prediction of liver steatosis dynamics in response to different diet and drug interventions.

BACKGROUND & AIMS: Non-alcoholic fatty liver disease (NAFLD) is a common metabolic disorder, characterized by the accumulation of excess fat in the liver, and is a driving factor for various severe liver diseases. These multi-factorial and multi-timescale changes are observed in different clinical studies, but these studies have not been integrated into a unified framework. In this study, we aim to present such a unified framework in the form of a dynamic mathematical model. METHODS: For model training and validation, we collected data for dietary or drug-induced interventions aimed at reducing or increasing liver fat. The model was formulated using ordinary differential equations (ODEs) and the mathematical analysis, model simulation, model formulation and the model parameter estimation were all performed in MATLAB. RESULTS: Our mathematical model describes accumulation of fat in the liver and predicts changes in lipid fluxes induced by both dietary and drug interventions. The model is validated using data from a wide range of drug and dietary intervention studies and can predict both short-term (days) and long-term (weeks) changes in liver fat. Importantly, the model computes the contribution of each individual lipid flux to the total liver fat dynamics. Furthermore, the model can be combined with an established bodyweight model, to simulate even longer scenarios (years), also including the effects of insulin resistance and body weight. To help prepare for corresponding eHealth applications, we also present a way to visualize the simulated changes, using dynamically changing lipid droplets, seen in images of liver biopsies. CONCLUSION: In conclusion, we believe that the minimal model presented herein might be a useful tool for future applications, and to further integrate and understand data regarding changes in dietary and drug induced changes in ectopic TAG in the liver. With further development and validation, the minimal model could be used as a disease progression model for steatosis.

The Wnt-dependent master regulator NKX1-2 controls mouse pre-implantation development.

Embryo size, specification, and homeostasis are regulated by a complex gene regulatory and signaling network. Here we used gene expression signatures of Wnt-activated mouse embryonic stem cell (mESC) clones to reverse engineer an mESC regulatory network. We identify NKX1-2 as a novel master regulator of preimplantation embryo development. We find that Nkx1-2 inhibition reduces nascent RNA synthesis, downregulates genes controlling ribosome biogenesis, RNA translation, and transport, and induces severe alteration of nucleolus structure, resulting in the exclusion of RNA polymerase I from nucleoli. In turn, NKX1-2 loss of function leads to chromosome missegregation in the 2- to 4-cell embryo stages, severe decrease in blastomere numbers, alterations of tight junctions (TJs), and impairment of microlumen coarsening. Overall, these changes impair the blastocoel expansion-collapse cycle and embryo cavitation, leading to altered lineage specification and developmental arrest.

Predicting metabolic modules in incomplete bacterial genomes with MetaPathPredict.

The reconstruction of complete microbial metabolic pathways using 'omics data from environmental samples remains challenging. Computational pipelines for pathway reconstruction that utilize machine learning methods to predict the presence or absence of KEGG modules in incomplete genomes are lacking. Here, we present MetaPathPredict, a software tool that incorporates machine learning models to predict the presence of complete KEGG modules within bacterial genomic datasets. Using gene annotation data and information from the KEGG module database, MetaPathPredict employs deep learning models to predict the presence of KEGG modules in a genome. MetaPathPredict can be used as a command line tool or as a Python module, and both options are designed to be run locally or on a compute cluster. Benchmarks show that MetaPathPredict makes robust predictions of KEGG module presence within highly incomplete genomes.

Experimental evolution for the recovery of growth loss due to genome reduction.

As the genome encodes the information crucial for cell growth, a sizeable genomic deficiency often causes a significant decrease in growth fitness. Whether and how the decreased growth fitness caused by genome reduction could be compensated by evolution was investigated here. Experimental evolution with an Escherichia coli strain carrying a reduced genome was conducted in multiple lineages for approximately 1000 generations. The growth rate, which largely declined due to genome reduction, was considerably recovered, associated with the improved carrying capacity. Genome mutations accumulated during evolution were significantly varied across the evolutionary lineages and were randomly localized on the reduced genome. Transcriptome reorganization showed a common evolutionary direction and conserved the chromosomal periodicity, regardless of highly diversified gene categories, regulons, and pathways enriched in the differentially expressed genes. Genome mutations and transcriptome reorganization caused by evolution, which were found to be dissimilar to those caused by genome reduction, must have followed divergent mechanisms in individual evolutionary lineages. Gene network reconstruction successfully identified three gene modules functionally differentiated, which were responsible for the evolutionary changes of the reduced genome in growth fitness, genome mutation, and gene expression, respectively. The diversity in evolutionary approaches improved the growth fitness associated with the homeostatic transcriptome architecture as if the evolutionary compensation for genome reduction was like all roads leading to Rome.

Model informed precision medicine of Chinese herbal medicines formulas-A multi-scale mechanistic intelligent model.

Recent trends suggest that Chinese herbal medicine formulas (CHM formulas) are promising treatments for complex diseases. To characterize the precise syndromes, precise diseases and precise targets of the precise targets between complex diseases and CHM formulas, we developed an artificial intelligence-based quantitative predictive algorithm (DeepTCM). DeepTCM has gone through multilevel model calibration and validation against a comprehensive set of herb and disease data so that it accurately captures the complex cellular signaling, molecular and theoretical levels of traditional Chinese medicine (TCM). As an example, our model simulated the optimal CHM formulas for the treatment of coronary heart disease (CHD) with depression, and through model sensitivity analysis, we calculated the balanced scoring of the formulas. Furthermore, we constructed a biological knowledge graph representing interactions by associating herb-target and gene-disease interactions. Finally, we experimentally confirmed the therapeutic effect and pharmacological mechanism of a novel model-predicted intervention in humans and mice. This novel multiscale model opened up a new avenue to combine "disease syndrome" and "macro micro" system modeling to facilitate translational research in CHM formulas.

Biomarkers of Alzheimer's Disease Associated with Programmed Cell Death Reveal Four Repurposed Drugs.

Alzheimer's disease (AD) is a neurodegenerative disorder and the most common cause of dementia. Programmed cell death (PCD) is mainly characterized by unique morphological features and energy-dependent biochemical processes. The predominant pathway leading to cell death in AD has not been thoroughly analyzed, although there is evidence of neuron loss in AD and numerous pathways of PCD have been associated with this process. A better understanding of the systems biology underlying the relationship between AD and PCD could lead to the development of new therapeutic approaches. To this end, publicly available transcriptome data were examined using bioinformatic methods such as differential gene expression and weighted gene coexpression network analysis (WGCNA) to find PCD-related AD biomarkers. The diagnostic significance of these biomarkers was evaluated using a logistic regression-based predictive model. Using these biomarkers, a multifactorial regulatory network was developed. Last, a drug repositioning study was conducted to propose new drugs for the treatment of AD targeting PCD. The development of 3PM (predictive, preventive, and personalized) drugs for the treatment of AD would be enabled by additional research on the effects of these drugs on this disease.

Positive and negative feedback regulation of the TGF-ß1 explains two equilibrium states in skin aging.

During aging, skin homeostasis is essential for maintaining appearance, as well as biological defense of the human body. In this study, we identified thrombospondin-1 (THBS1) and fibromodulin (FMOD) as positive and negative regulators, respectively, of the TGF-ß1-SMAD4 axis in human skin aging, based on in vitro and in vivo omics analyses and mathematical modeling. Using transcriptomic and epigenetic analyses of senescent dermal fibroblasts, TGF-ß1 was identified as the key upstream regulator. Bifurcation analysis revealed a binary high-/low-TGF-ß1 switch, with THBS1 as the main controller. Computational simulation of the TGF-ß1 signaling pathway indicated that THBS1 expression was sensitively regulated, whereas FMOD was regulated robustly. Results of sensitivity analysis and validation showed that inhibition of SMAD4 complex formation was a promising method to control THBS1 production and senescence. Therefore, this study demonstrated the potential of combining data-driven target discovery with mathematical approaches to determine the mechanisms underlying skin aging.

Systems Approaches to Cell Culture-Derived Extracellular Vesicles for Acute Kidney Injury Therapy: Prospects and Challenges.

Acute kidney injury (AKI) is a heterogeneous syndrome, comprising diverse etiologies of kidney insults that result in high mortality and morbidity if not well managed. Although great efforts have been made to investigate underlying pathogenic mechanisms of AKI, there are limited therapeutic strategies available. Extracellular vesicles (EV) are membrane-bound vesicles secreted by various cell types, which can serve as cell-free therapy through transfer of bioactive molecules. In this review, we first overview the AKI syndrome and EV biology, with a particular focus on the technical aspects and therapeutic application of cell culture-derived EVs. Second, we illustrate how multi-omic approaches to EV miRNA, protein, and genomic cargo analysis can yield new insights into their mechanisms of action and address unresolved questions in the field. We then summarize major experimental evidence regarding the therapeutic potential of EVs in AKI, which we subdivide into stem cell and non-stem cell-derived EVs. Finally, we highlight the challenges and opportunities related to the clinical translation of animal studies into human patients.

[Prevention and treatment of COVID-19 based on post-genomic pharmacological analysis: Systematic computer analysis of 290,000 scientific articles on COVID-19].

The COVID-19 pandemic has highlighted pressing challenges in biomedical research methodology. It has become obvious that the rapid and effective development of treatments for "new" viral infections is impossible without the coordination of interdisciplinary research and in-depth analysis of data obtained within the framework of the post-genomic paradigm. Presents the results of a systematic computer analysis of 290,000 scientific articles on COVID-19, with an emphasis on the results of post-genomic studies of SARS-CoV-2. The futility of the overly simplified approach, which considers only one "most important receptor protein", only one "key virus gene", etc., is shown. It is shown how post-genomic technologies will make it possible to find informative biomarkers of severe coronavirus infection, including those based on complex immune disorders associated with COVID-19.

ALK inhibitors suppress HCC and synergize with anti-PD-1 therapy and ABT-263 in preclinical models.

Hepatocellular carcinoma (HCC) currently lacks effective therapies, leaving a critical need for new treatment options. A previous study identified the anaplastic lymphoma kinase (ALK) amplification in HCC patients, raising the question of whether ALK inhibitors could be a viable treatment. Here, we showed that both ALK inhibitors and ALK knockout effectively halted HCC growth in cell cultures. Lorlatinib, a potent ALK inhibitor, suppressed HCC tumor growth and metastasis across various mouse models. Additionally, in an advanced immunocompetent humanized mouse model, when combined with an anti-PD-1 antibody, lorlatinib more potently suppressed HCC tumor growth, surpassing individual drug efficacy. Lorlatinib induced apoptosis and senescence in HCC cells, and the senolytic agent ABT-263 enhanced the efficacy of lorlatinib. Additional studies identified that the apoptosis-inducing effect of lorlatinib was mediated via GGN and NRG4. These findings establish ALK inhibitors as promising HCC treatments, either alone or in combination with immunotherapies or senolytic agents.

Single-cell ionic current phenotyping elucidates non-canonical features and predictive potential of cardiomyocytes during automated drug experiments.

All new drugs must go through preclinical screening tests to determine their proarrhythmic potential. While these assays effectively filter out dangerous drugs, they are too conservative, often misclassifying safe compounds as proarrhythmic. In this study, we attempt to address this shortcoming with a novel, medium-throughput drug-screening approach: we use an automated patch-clamp system to acquire optimized voltage clamp (VC) and action potential (AP) data from human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) at several drug concentrations (baseline, 3×, 10× and 20× the effective free plasma concentrations). With our novel method, we show correlations between INa block and upstroke slowing after treatment with flecainide or quinine. Additionally, after quinine treatment, we identify significant reductions in current during voltage steps designed to isolate If and IKs. However, we do not detect any IKr block by either drug, and upon further investigation, do not see any IKr present in the iPSC-CMs when prepared for automated patch experiments (i.e. in suspension) - this is in contrast to similar experiments we have conducted with these cells using the manual patch setup. In this study, we: (1) present a proof-of-concept demonstration of a single-cell medium-throughput drug study, and (2) characterize the non-canonical electrophysiology of iPSC-CMs when prepared for experiments in a medium-throughput setting. KEY POINTS: Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) offer potential as an in vitro model to study the proarrhythmic potential of drugs, but insights from these cells are often limited by the low throughput of manual patch-clamp. In this study, we use a medium-throughput automated patch-clamp system to acquire action potential (AP) and complex voltage clamp (VC) data from single iPSC-CMs at multiple drug concentrations. A correlation between AP upstroke and INa transients was identified and drug-induced changes in ionic currents found. We also characterize the substantially altered physiology of iPSC-CMs when patched in an automated system, suggesting the need to investigate differences between manual and automated patch experiments.