Search

The human transcriptome comprises a complex network of coding and non-coding RNAs implicated in a myriad of biological functions. Non-coding RNAs exhibit highly organized spatial and temporal expression patterns and are emerging as critical regulators of differentiation, homeostasis, and pathological states, including in the cardiovascular system. This review defines the current knowledge gaps, unmet methodological needs, and describes the challenges in dissecting and understanding the role and regulation of the non-coding transcriptome in cardiovascular disease. These challenges include poor annotation of the non-coding genome, determination of the cellular distribution of transcripts, assessment of the role of RNA processing and identification of cell-type specific changes in cardiovascular physiology and disease. We highlight similarities and differences in the hurdles associated with the analysis of the non-coding and protein-coding transcriptomes. In addition, we discuss how the lack of consensus and absence of standardized methods affect reproducibility of data. These shortcomings should be defeated in order to make significant scientific progress and foster the development of clinically applicable non-coding RNA-based therapeutic strategies to lessen the burden of cardiovascular disease. [ABSTRACT FROM AUTHOR]

Halomonas elongata thrives in hypersaline environments producing polyhydroxyalkanoates (PHAs) and osmoprotectants such as ectoine. Despite its biotechnological importance, several aspects of the dynamics of its metabolism remain elusive. Here, we construct and validate a genome-scale metabolic network model for H. elongata 153B. Then, we investigate the flux distribution dynamics during optimal growth, ectoine, and PHA biosynthesis using statistical methods, and a pipeline based on shadow prices. Lastly, we use optimization algorithms to uncover novel engineering targets to increase PHA production. The resulting model (iEB1239) includes 1534 metabolites, 2314 reactions, and 1239 genes. iEB1239 can reproduce growth on several carbon sources and predict growth on previously unreported ones. It also reproduces biochemical phenotypes related to Oad and Ppc gene functions in ectoine biosynthesis. A flux distribution analysis during optimal ectoine and PHA biosynthesis shows decreased energy production through oxidative phosphorylation. Furthermore, our analysis unveils a diverse spectrum of metabolic alterations that extend beyond mere flux changes to encompass heightened precursor production for ectoine and PHA synthesis. Crucially, these findings capture other metabolic changes linked to adaptation in hypersaline environments. Bottlenecks in the glycolysis and fatty acid metabolism pathways are identified, in addition to PhaC, which has been shown to increase PHA production when overexpressed. Overall, our pipeline demonstrates the potential of genome-scale metabolic models in combination with statistical approaches to obtain insights into the metabolism of H. elongata. Our platform can be exploited for researching environmental adaptation, and for designing and optimizing metabolic engineering strategies for bioproduct synthesis., (© 2024 The Author(s). Biotechnology Journal published by Wiley‐VCH GmbH.)

Deep-learning tools that extract prognostic factors derived from multi-omics data have recently contributed to individualized predictions of survival outcomes. However, the limited size of integrated omics-imaging-clinical datasets poses challenges. Here, we propose two biologically interpretable and robust deep-learning architectures for survival prediction of non-small cell lung cancer (NSCLC) patients, learning simultaneously from computed tomography (CT) scan images, gene expression data, and clinical information. The proposed models integrate patient-specific clinical, transcriptomic, and imaging data and incorporate Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathway information, adding biological knowledge within the learning process to extract prognostic gene biomarkers and molecular pathways. While both models accurately stratify patients in high- and low-risk groups when trained on a dataset of only 130 patients, introducing a cross-attention mechanism in a sparse autoencoder significantly improves the performance, highlighting tumor regions and NSCLC-related genes as potential biomarkers and thus offering a significant methodological advancement when learning from small imaging-omics-clinical samples., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Gut microbiome dysbiosis is linked to many neurological disorders including Alzheimer's disease (AD). A major risk factor for AD is polymorphism in the apolipoprotein E ( APOE ) gene, which affects gut microbiome composition. To explore the gut-brain axis in AD, long-lived animal models of naturally developing AD-like pathologies are needed. Octodon degus (degu) exhibit spontaneous AD-like symptoms and ApoE mutations, making them suitable for studying the interplay between AD genetic determinants and gut microbiome. We analyzed the association between APOE genotype and gut microbiome in 50 humans and 32 degu using16S rRNA gene amplicon sequencing. Significant associations were found between the degu ApoE mutation and gut microbial changes in degu, notably a depletion of Ruminococcaceae and Akkermansiaceae and an enrichment of Prevotellaceae , mirroring patterns seen in people with AD. The altered taxa were previously suggested to be involved in AD, validating the degu as an unconventional model for studying the AD/microbiome crosstalk., Competing Interests: The authors declare no competing interests., (© 2024 The Authors.)

Alterations in Dp71 expression, the most ubiquitous dystrophin isoform, have been associated with patient survival across tumours. Intriguingly, in certain malignancies, Dp71 acts as a tumour suppressor, while manifesting oncogenic properties in others. This diversity could be explained by the expression of two Dp71 splice variants encoding proteins with distinct C-termini, each with specific properties. Expression of these variants has impeded the exploration of their unique roles. Using CRISPR/Cas9, we ablated the Dp71f variant with the alternative C-terminus in a sarcoma cell line not expressing the canonical C-terminal variant, and conducted molecular (RNAseq) and functional characterisation of the knockout cells. Dp71f ablation induced major transcriptomic alterations, particularly affecting the expression of genes involved in calcium signalling and ECM-receptor interaction pathways. The genome-scale metabolic analysis identified significant downregulation of glucose transport via membrane vesicle reaction (GLCter) and downregulated glycolysis/gluconeogenesis pathway. Functionally, these molecular changes corresponded with, increased calcium responses, cell adhesion, proliferation, survival under serum starvation and chemotherapeutic resistance. Knockout cells showed reduced GLUT1 protein expression, survival without attachment and their migration and invasion in vitro and in vivo were unaltered, despite increased matrix metalloproteinases release. Our findings emphasise the importance of alternative splicing of dystrophin transcripts and underscore the role of the Dp71f variant, which appears to govern distinct cellular processes frequently dysregulated in tumour cells. The loss of this regulatory mechanism promotes sarcoma cell survival and treatment resistance. Thus, Dp71f is a target for future investigations exploring the intricate functions of specific DMD transcripts in physiology and across malignancies., (© The Author(s) 2024. Published by Oxford University Press.)

Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels., Competing Interests: Declaration of interests The authors declare no conflicts of interest., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Artificial intelligence (AI) is widely used for exploiting multimodal biomedical data, with increasingly accurate predictions and model-agnostic interpretations, which are however also agnostic to biological mechanisms. Combining metabolic modelling, 'omics, and imaging data via multimodal AI can generate predictions that can be interpreted mechanistically and transparently, therefore with significantly higher therapeutic potential., Competing Interests: Declaration of interests None declared by authors., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Background: Melanoma, the most lethal skin cancer type, occurs more frequently in Parkinson's disease (PD), and PD is more frequent in melanoma patients, suggesting disease mechanisms overlap. α-synuclein, a protein that accumulates in PD brain, and the oncogene DJ-1, which is associated with PD autosomal recessive forms, are both elevated in melanoma cells. Whether this indicates melanoma progression or constitutes a protective response remains unclear. We hereby investigated the molecular mechanisms through which α-synuclein and DJ-1 interact, suggesting novel biomarkers and targets in melanoma., Methods: The Cancer Genome Atlas (TCGA) expression profiles derived from UCSC Xena were used to obtain α-synuclein and DJ-1 expression and correlated with survival in skin cutaneous melanoma (SKCM). Immunohistochemistry determined the expression in metastatic melanoma lymph nodes. Protein-protein interactions (PPIs) and molecular docking assessed protein binding and affinity with chemotherapeutic drugs. Further validation was performed using in vitro cellular models and ELISA immunoassays., Results: α-synuclein and DJ-1 were upregulated in primary and metastatic SKCM. Aggregated α-synuclein was selectively detected in metastatic melanoma lymph nodes. α-synuclein overexpression in SK-MEL-28 cells induced the expression of DJ-1, supporting PPI and a positive correlation in melanoma patients. Molecular docking revealed a stable protein complex, with differential binding to chemotherapy drugs such as temozolomide, dacarbazine, and doxorubicin. Parallel reduction of both proteins in temozolomide-treated SK-MEL-28 spheroids suggests drug binding may affect protein interaction and/or stability., Conclusion: α-synuclein, together with DJ-1, may play a role in melanoma progression and chemosensitivity, constituting novel targets for therapeutic intervention, and possible biomarkers for melanoma., (© 2024 The Authors. Cancer Medicine published by John Wiley & Sons Ltd.)

Data are the most important elements of bioinformatics: Computational analysis of bioinformatics data, in fact, can help researchers infer new knowledge about biology, chemistry, biophysics, and sometimes even medicine, influencing treatments and therapies for patients. Bioinformatics and high-throughput biological data coming from different sources can even be more helpful, because each of these different data chunks can provide alternative, complementary information about a specific biological phenomenon, similar to multiple photos of the same subject taken from different angles. In this context, the integration of bioinformatics and high-throughput biological data gets a pivotal role in running a successful bioinformatics study. In the last decades, data originating from proteomics, metabolomics, metagenomics, phenomics, transcriptomics, and epigenomics have been labelled -omics data, as a unique name to refer to them, and the integration of these omics data has gained importance in all biological areas. Even if this omics data integration is useful and relevant, due to its heterogeneity, it is not uncommon to make mistakes during the integration phases. We therefore decided to present these ten quick tips to perform an omics data integration correctly, avoiding common mistakes we experienced or noticed in published studies in the past. Even if we designed our ten guidelines for beginners, by using a simple language that (we hope) can be understood by anyone, we believe our ten recommendations should be taken into account by all the bioinformaticians performing omics data integration, including experts., Competing Interests: The authors declare they have no conflict of interest., (Copyright: © 2023 Chicco et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The number of "omics" approaches is continuously growing. Among others, epigenetics has appeared as an attractive area of investigation by the cardiovascular research community, notably considering its association with disease development. Complex diseases such as cardiovascular diseases have to be tackled using methods integrating different omics levels, so called "multi-omics" approaches. These approaches combine and co-analyze different levels of disease regulation. In this review, we present and discuss the role of epigenetic mechanisms in regulating gene expression and provide an integrated view of how these mechanisms are interlinked and regulate the development of cardiac disease, with a particular attention to heart failure. We focus on DNA, histone, and RNA modifications, and discuss the current methods and tools used for data integration and analysis. Enhancing the knowledge of these regulatory mechanisms may lead to novel therapeutic approaches and biomarkers for precision healthcare and improved clinical outcomes., (© 2023. The Author(s).)

Multi-strain probiotics are widely regarded as effective products for improving gut microbiota stability and host health, providing advantages over single-strain probiotics. However, in general, it is unclear to what extent different strains would cooperate or compete for resources, and how the establishment of a common biofilm microenvironment could influence their interactions. In this work, we develop an integrative experimental and computational approach to comprehensively assess the metabolic functionality and interactions of probiotics across growth conditions. Our approach combines co-culture assays with genome-scale modelling of metabolism and multivariate data analysis, thus exploiting complementary data- and knowledge-driven systems biology techniques. To show the advantages of the proposed approach, we apply it to the study of the interactions between two widely used probiotic strains of Lactobacillus reuteri and Saccharomyces boulardii, characterising their production potential for compounds that can be beneficial to human health. Our results show that these strains can establish a mixed cooperative-antagonistic interaction best explained by competition for shared resources, with an increased individual exchange but an often decreased net production of amino acids and short-chain fatty acids. Overall, our work provides a strategy that can be used to explore microbial metabolic fingerprints of biotechnological interest, capable of capturing multifaceted equilibria even in simple microbial consortia., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Gliomas are the most common brain tumours with the high-grade glioblastoma representing the most aggressive and lethal form. Currently, there is a lack of specific glioma biomarkers that would aid tumour subtyping and minimally invasive early diagnosis. Aberrant glycosylation is an important post-translational modification in cancer and is implicated in glioma progression. Raman spectroscopy (RS), a vibrational spectroscopic label-free technique, has already shown promise in cancer diagnostics., Methods: RS was combined with machine learning to discriminate glioma grades. Raman spectral signatures of glycosylation patterns were used in serum samples and fixed tissue biopsy samples, as well as in single cells and spheroids., Results: Glioma grades in fixed tissue patient samples and serum were discriminated with high accuracy. Discrimination between higher malignant glioma grades (III and IV) was achieved with high accuracy in tissue, serum, and cellular models using single cells and spheroids. Biomolecular changes were assigned to alterations in glycosylation corroborated by analysing glycan standards and other changes such as carotenoid antioxidant content., Conclusion: RS combined with machine learning could pave the way for more objective and less invasive grading of glioma patients, serving as a useful tool to facilitate glioma diagnosis and delineate biomolecular glioma progression changes., (© 2023. Crown.)

Multi-omics data integration via mechanistic models of metabolism is a scalable and flexible framework for exploring biological hypotheses in microbial systems. However, although most microorganisms are unculturable, such multi-omics modeling is limited to isolate microbes or simple synthetic communities. Here, we developed an approach for modeling microbial activity and interactions that leverages the reconstruction of metagenome-assembled genomes and associated genome-centric metatranscriptomes. At its core, we designed a method for condition-specific metabolic modeling of microbial communities through the integration of metatranscriptomic data. Using this approach, we explored the behavior of anaerobic digestion consortia driven by hydrogen availability and human gut microbiota dysbiosis associated with Crohn's disease, identifying condition-dependent amino acid requirements in archaeal species and a reduced short-chain fatty acid exchange network associated with disease, respectively. Our approach can be applied to complex microbial communities, allowing a mechanistic contextualization of multi-omics data on a metagenome scale., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

Breast cancer is one of the most common cancers in women worldwide, which causes an enormous number of deaths annually. However, early diagnosis of breast cancer can improve survival outcomes enabling simpler and more cost-effective treatments. The recent increase in data availability provides unprecedented opportunities to apply data-driven and machine learning methods to identify early-detection prognostic factors capable of predicting the expected survival and potential sensitivity to treatment of patients, with the final aim of enhancing clinical outcomes. This tutorial presents a protocol for applying machine learning models in survival analysis for both clinical and transcriptomic data. We show that integrating clinical and mRNA expression data is essential to explain the multiple biological processes driving cancer progression. Our results reveal that machine-learning-based models such as random survival forests, gradient boosted survival model, and survival support vector machine can outperform the traditional statistical methods, i.e., Cox proportional hazard model. The highest C-index among the machine learning models was recorded when using survival support vector machine, with a value 0.688, whereas the C-index recorded using the Cox model was 0.677. Shapley Additive Explanation (SHAP) values were also applied to identify the feature importance of the models and their impact on the prediction outcomes., (© 2023. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Background: Recently, multi-omic machine learning architectures have been proposed for the early detection of cancer. However, for rare cancers and their associated small datasets, it is still unclear how to use the available multi-omics data to achieve a mechanistic prediction of cancer onset and progression, due to the limited data available. Hepatoblastoma is the most frequent liver cancer in infancy and childhood, and whose incidence has been lately increasing in several developed countries. Even though some studies have been conducted to understand the causes of its onset and discover potential biomarkers, the role of metabolic rewiring has not been investigated in depth so far., Methods: Here, we propose and implement an interpretable multi-omics pipeline that combines mechanistic knowledge from genome-scale metabolic models with machine learning algorithms, and we use it to characterise the underlying mechanisms controlling hepatoblastoma., Results and Conclusions: While the obtained machine learning models generally present a high diagnostic classification accuracy, our results show that the type of omics combinations used as input to the machine learning models strongly affects the detection of important genes, reactions and metabolic pathways linked to hepatoblastoma. Our method also suggests that, in the context of computer-aided diagnosis of cancer, optimal diagnostic accuracy can be achieved by adopting a combination of omics that depends on the patient's clinical characteristics., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Glioblastoma is the most aggressive form of brain cancer, presenting poor prognosis despite current advances in treatment. There is therefore an urgent need for novel biomarkers and therapeutic targets. Interactions between mucin 4 (MUC4) and the epidermal growth factor receptor (EGFR) are involved in carcinogenesis, and may lead to matrix metalloproteinase-9 (MMP9) overexpression, exacerbating cancer cell invasiveness. In this study, the role of MUC4, MMP9, and EGFR in the progression and clinical outcome of glioma patients was investigated. Immunohistochemistry (IHC) and immunofluorescence (IF) in fixed tissue samples of glioma patients were used to evaluate the expression and localization of EGFR, MMP9, and MUC4. Kaplan-Meier survival analysis was also performed to test the prognostic utility of the proteins for glioma patients. The protein levels were assessed with enzyme-linked immunosorbent assay (ELISA) in serum of glioma patients, to further investigate their potential as non-invasive serum biomarkers. We demonstrated that MUC4 and MMP9 are both significantly upregulated during glioma progression. Moreover, MUC4 is co-expressed with MMP9 and EGFR in the proliferative microvasculature of glioblastoma, suggesting a potential role for MUC4 in microvascular proliferation and angiogenesis. The combined high expression of MUC4/MMP9, and MUC4/MMP9/EGFR was associated with poor overall survival (OS). Finally, MMP9 mean protein level was significantly higher in the serum of glioblastoma compared with grade III glioma patients, whereas MUC4 mean protein level was minimally elevated in higher glioma grades (III and IV) compared with control. Our results suggest that MUC4, along with MMP9, might account for glioblastoma progression, representing potential therapeutic targets, and suggesting the 'MUC4/MMP9/EGFR axis' may play a vital role in glioblastoma diagnostics., (© 2022. Crown.)

Salt tolerant organisms are increasingly being used for the industrial production of high-value biomolecules due to their better adaptability compared to mesophiles. Chromohalobacter canadensis is one of the early halophiles to show promising biotechnology potential, which has not been explored to date. Advanced high throughput technologies such as whole-genome sequencing allow in-depth insight into the potential of organisms while at the frontiers of systems biology. At the same time, genome-scale metabolic models (GEMs) enable phenotype predictions through a mechanistic representation of metabolism. Here, we sequence and analyze the genome of C. canadensis 85B, and we use it to reconstruct a GEM. We then analyze the GEM using flux balance analysis and validate it against literature data on C. canadensis. We show that C. canadensis 85B is a metabolically versatile organism with many features for stress and osmotic adaptation. Pathways to produce ectoine and polyhydroxybutyrates were also predicted. The GEM reveals the ability to grow on several carbon sources in a minimal medium and reproduce osmoadaptation phenotypes. Overall, this study reveals insights from the genome of C. canadensis 85B, providing genomic data and a draft GEM that will serve as the first steps towards a better understanding of its metabolism, for novel applications in industrial biotechnology., (© 2022 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd.)

Duchenne muscular dystrophy (DMD) affects myofibers and muscle stem cells, causing progressive muscle degeneration and repair defects. It was unknown whether dystrophic myoblasts-the effector cells of muscle growth and regeneration-are affected. Using transcriptomic, genome-scale metabolic modelling and functional analyses, we demonstrate, for the first time, convergent abnormalities in primary mouse and human dystrophic myoblasts. In Dmd mdx myoblasts lacking full-length dystrophin, the expression of 170 genes was significantly altered. Myod1 and key genes controlled by MyoD ( Myog, Mymk, Mymx , epigenetic regulators, ECM interactors, calcium signalling and fibrosis genes) were significantly downregulated. Gene ontology analysis indicated enrichment in genes involved in muscle development and function. Functionally, we found increased myoblast proliferation, reduced chemotaxis and accelerated differentiation, which are all essential for myoregeneration. The defects were caused by the loss of expression of full-length dystrophin, as similar and not exacerbated alterations were observed in dystrophin-null Dmd mdx-βgeo myoblasts. Corresponding abnormalities were identified in human DMD primary myoblasts and a dystrophic mouse muscle cell line, confirming the cross-species and cell-autonomous nature of these defects. The genome-scale metabolic analysis in human DMD myoblasts showed alterations in the rate of glycolysis/gluconeogenesis, leukotriene metabolism, and mitochondrial beta-oxidation of various fatty acids. These results reveal the disease continuum: DMD defects in satellite cells, the myoblast dysfunction affecting muscle regeneration, which is insufficient to counteract muscle loss due to myofiber instability. Contrary to the established belief, our data demonstrate that DMD abnormalities occur in myoblasts, making these cells a novel therapeutic target for the treatment of this lethal disease., Competing Interests: MG, VM, MB, SV, AO, JR, LB, MK, CA, CP, DG No competing interests declared, SR SR received royalties from Cancer Research Technology Ltd through University of Cambridge Enterprise Ltd under agreement A13086 for "Chromatin-bound METTL3 m6A methyltransferase facilitates myeloid leukaemogenesis by modulating ribosome stalling", (© 2022, Gosselin et al.)

The degu ( Octodon degus ) is a diurnal long-lived rodent that can spontaneously develop molecular and behavioral changes that mirror those seen in human aging. With age some degu, but not all individuals, develop cognitive decline and brain pathology like that observed in Alzheimer's disease including neuroinflammation, hyperphosphorylated tau and amyloid plaques, together with other co-morbidities associated with aging such as macular degeneration, cataracts, alterations in circadian rhythm, diabetes and atherosclerosis. Here we report the whole-genome sequencing and analysis of the degu genome, which revealed unique features and molecular adaptations consistent with aging and Alzheimer's disease. We identified single nucleotide polymorphisms in genes associated with Alzheimer's disease including a novel apolipoprotein E ( Apoe ) gene variant that correlated with an increase in amyloid plaques in brain and modified the in silico predicted degu APOE protein structure and functionality. The reported genome of an unconventional long-lived animal model of aging and Alzheimer's disease offers the opportunity for understanding molecular pathways involved in aging and should help advance biomedical research into treatments for Alzheimer's disease., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Hurley, Urra, Garduno, Bruno, Kimbell, Wilkinson, Marino-Buslje, Ezquer, Ezquer, Aburto, Poulin, Vasquez, Deacon, Avila, Altimiras, Whitney Vanderklish, Zampieri, Angione, Constantino, Holmes, Coba, Xu and Cogram.)

The continuous spread and evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and the rapid surge in infection cases in the coronavirus disease 2019 (COVID-19) evoke a dire need for effective therapeutics. In this study, we explored the inhibitory potential of a library of 605 phytocompounds, selected from Indian medicinal plants with reported antiviral and anti-inflammatory activities, against the receptor-binding domain of spike proteins of the SARS-CoV-2 wild-type and the variants of concern, including variants B.1.1.7 (Alpha), B.1.351 (Beta), P.1 (Gamma), B.1.617.2 (Delta), and B.1.1.529 (Omicron). Our approach was based on extensive molecular docking, assessment of drug-likeness, and robust molecular dynamics simulations. We also identified promising inhibitory candidates against the host (human) proteins associated with SARS-CoV-2 spike activation and attachment, namely, ACE2 receptor, proteases TMPRSS2 and CTSL, and the endocytic regulator AAK1. In addition, we screened promising inhibitory compounds against the human proinflammatory cytokines- IL-6, IL-1β, TNF-α, and IFN-γ, that are associated with the adverse cytokine storm in COVID-19 patients. Our analysis returned an encouraging list of promising inhibitory candidates that includes: abietatriene against the spike proteins of the SARS-CoV-2 wild-type and the variants of concern; taraxerol against the human ACE2, CTSL and TNF-α; β-amyrin against the human TMPRSS2; cynaroside against the human AAK1 and IL-1β; and friedelin against the human IL-6 and IFN-γ. Our findings provide substantial evidence for the inhibitory potential of these compounds and encourage further in vitro and in vivo studies to validate their use as safe and effective therapeutics against COVID-19., (© 2022 Wiley Periodicals LLC.)

In this proof-of-concept work, we evaluate the performance of multiple machine-learning methods as surrogate models for use in the analysis of agent-based models (ABMs). Analysing agent-based modelling outputs can be challenging, as the relationships between input parameters can be non-linear or even chaotic even in relatively simple models, and each model run can require significant CPU time. Surrogate modelling, in which a statistical model of the ABM is constructed to facilitate detailed model analyses, has been proposed as an alternative to computationally costly Monte Carlo methods. Here we compare multiple machine-learning methods for ABM surrogate modelling in order to determine the approaches best suited as a surrogate for modelling the complex behaviour of ABMs. Our results suggest that, in most scenarios, artificial neural networks (ANNs) and gradient-boosted trees outperform Gaussian process surrogates, currently the most commonly used method for the surrogate modelling of complex computational models. ANNs produced the most accurate model replications in scenarios with high numbers of model runs, although training times were longer than the other methods. We propose that agent-based modelling would benefit from using machine-learning methods for surrogate modelling, as this can facilitate more robust sensitivity analyses for the models while also reducing CPU time consumption when calibrating and analysing the simulation., Competing Interests: The authors have declared that no competing interests exist.

Motivation: Gene regulation is responsible for controlling numerous physiological functions and dynamically responding to environmental fluctuations. Reconstructing the human network of gene regulatory interactions is thus paramount to understanding the cell functional organization across cell types, as well as to elucidating pathogenic processes and identifying molecular drug targets. Although significant effort has been devoted towards this direction, existing computational methods mainly rely on gene expression levels, possibly ignoring the information conveyed by mechanistic biochemical knowledge. Moreover, except for a few recent attempts, most of the existing approaches only consider the information of the organism under analysis, without exploiting the information of related model organisms., Results: We propose a novel method for the reconstruction of the human gene regulatory network, based on a transfer learning strategy that synergically exploits information from human and mouse, conveyed by gene-related metabolic features generated in silico from gene expression data. Specifically, we learn a predictive model from metabolic activity inferred via tissue-specific metabolic modelling of artificial gene knockouts. Our experiments show that the combination of our transfer learning approach with the constructed metabolic features provides a significant advantage in terms of reconstruction accuracy, as well as additional clues on the contribution of each constructed metabolic feature., Availability and Implementation: The method, the datasets and all the results obtained in this study are available at: https://doi.org/10.6084/m9.figshare.c.5237687., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Complex, distributed, and dynamic sets of clinical biomedical data are collectively referred to as multimodal clinical data. In order to accommodate the volume and heterogeneity of such diverse data types and aid in their interpretation when they are combined with a multi-scale predictive model, machine learning is a useful tool that can be wielded to deconstruct biological complexity and extract relevant outputs. Additionally, genome-scale metabolic models (GSMMs) are one of the main frameworks striving to bridge the gap between genotype and phenotype by incorporating prior biological knowledge into mechanistic models. Consequently, the utilization of GSMMs as a foundation for the integration of multi-omic data originating from different domains is a valuable pursuit towards refining predictions. In this chapter, we show how cancer multi-omic data can be analyzed via multimodal machine learning and metabolic modeling. Firstly, we focus on the merits of adopting an integrative systems biology led approach to biomedical data mining. Following this, we propose how constraint-based metabolic models can provide a stable yet adaptable foundation for the integration of multimodal data with machine learning. Finally, we provide a step-by-step tutorial for the combination of machine learning and GSMMs, which includes: (i) tissue-specific constraint-based modeling; (ii) survival analysis using time-to-event prediction for cancer; and (iii) classification and regression approaches for multimodal machine learning. The code associated with the tutorial can be found at https://github.com/Angione-Lab/Tutorials_Combining_ML_and_GSMM ., (© 2022. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

Computational modelling of metabolic processes has proven to be a useful approach to formulate our knowledge and improve our understanding of core biochemical systems that are crucial to maintaining cellular functions. Towards understanding the broader role of metabolism on cellular decision-making in health and disease conditions, it is important to integrate the study of metabolism with other core regulatory systems and omics within the cell, including gene expression patterns. After quantitatively integrating gene expression profiles with a genome-scale reconstruction of human metabolism, we propose a set of combinatorial methods to reverse engineer gene expression profiles and to find pairs and higher-order combinations of genetic modifications that simultaneously optimize multi-objective cellular goals. This enables us to suggest classes of transcriptomic profiles that are most suitable to achieve given metabolic phenotypes. We demonstrate how our techniques are able to compute beneficial, neutral or "toxic" combinations of gene expression levels. We test our methods on nine tissue-specific cancer models, comparing our outcomes with the corresponding normal cells, identifying genes as targets for potential therapies. Our methods open the way to a broad class of applications that require an understanding of the interplay among genotype, metabolism, and cellular behaviour, at scale.

Motivation: High-throughput biological data, thanks to technological advances, have become cheaper to collect, leading to the availability of vast amounts of omic data of different types. In parallel, the in silico reconstruction and modeling of metabolic systems is now acknowledged as a key tool to complement experimental data on a large scale. The integration of these model- and data-driven information is therefore emerging as a new challenge in systems biology, with no clear guidance on how to better take advantage of the inherent multisource and multiomic nature of these data types while preserving mechanistic interpretation., Results: Here, we investigate different regularization techniques for high-dimensional data derived from the integration of gene expression profiles with metabolic flux data, extracted from strain-specific metabolic models, to improve cellular growth rate predictions. To this end, we propose ad-hoc extensions of previous regularization frameworks including group, view-specific and principal component regularization and experimentally compare them using data from 1143 Saccharomyces cerevisiae strains. We observe a divergence between methods in terms of regression accuracy and integration effectiveness based on the type of regularization employed. In multiomic regression tasks, when learning from experimental and model-generated omic data, our results demonstrate the competitiveness and ease of interpretation of multimodal regularized linear models compared to data-hungry methods based on neural networks., Availability and Implementation: All data, models and code produced in this work are available on GitHub at https://github.com/Angione-Lab/HybridGroupIPFLasso_pc2Lasso., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Combining a computational framework for flux balance analysis with machine learning improves the accuracy of predicting metabolic activity across conditions, while enabling mechanistic interpretation. This protocol presents a guide to condition-specific metabolic modeling that integrates regularized flux balance analysis with machine learning approaches to extract key features from transcriptomic and fluxomic data. We demonstrate the protocol as applied to Synechococcus sp. PCC 7002; we also outline how it can be adapted to any species or community with available multi-omic data. For complete details on the use and execution of this protocol, please refer to Vijayakumar et al. (2020)., Competing Interests: The authors declare no competing interests., (© 2021 The Author(s).)

Cancer is considered a high-risk condition for severe illness resulting from COVID-19. The interaction between severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and human metabolism is key to elucidating the risk posed by COVID-19 for cancer patients and identifying effective treatments, yet it is largely uncharacterised on a mechanistic level. We present a genome-scale map of short-term metabolic alterations triggered by SARS-CoV-2 infection of cancer cells. Through transcriptomic- and proteomic-informed genome-scale metabolic modelling, we characterise the role of RNA and fatty acid biosynthesis in conjunction with a rewiring in energy production pathways and enhanced cytokine secretion. These findings link together complementary aspects of viral invasion of cancer cells, while providing mechanistic insights that can inform the development of treatment strategies., (© 2021 The Authors. FEBS Letters published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Today's most troublesome population health challenges are often driven by social and environmental determinants, which are difficult to model using traditional epidemiological methods. We agree with those who have argued for the wider adoption of agent-based modelling (ABM) in taking on these challenges. However, while ABM has been used occasionally in population health, we argue that for ABM to be most effective in the field it should be used as a means for answering questions normally inaccessible to the traditional epidemiological toolkit. In an effort to clearly illustrate the utility of ABM for population health research, and to clear up persistent misunderstandings regarding the method's conceptual underpinnings, we offer a detailed presentation of the core concepts of complex systems theory, and summarise why simulations are essential to the study of complex systems. We then examine the current state of the art in ABM for population health, and propose they are well-suited for the study of the 'wicked' problems in population health, and could make significant contributions to theory and intervention development in these areas., (© 2021. The Author(s).)

Cancer is a complex disease that deregulates cellular functions at various molecular levels (e.g., DNA, RNA, and proteins). Integrated multi-omics analysis of data from these levels is necessary to understand the aberrant cellular functions accountable for cancer and its development. In recent years, Deep Learning (DL) approaches have become a useful tool in integrated multi-omics analysis of cancer data. However, high dimensional multi-omics data are generally imbalanced with too many molecular features and relatively few patient samples. This imbalance makes a DL based integrated multi-omics analysis difficult. DL-based dimensionality reduction technique, including variational autoencoder (VAE), is a potential solution to balance high dimensional multi-omics data. However, there are few VAE-based integrated multi-omics analyses, and they are limited to pancancer. In this work, we did an integrated multi-omics analysis of ovarian cancer using the compressed features learned through VAE and an improved version of VAE, namely Maximum Mean Discrepancy VAE (MMD-VAE). First, we designed and developed a DL architecture for VAE and MMD-VAE. Then we used the architecture for mono-omics, integrated di-omics and tri-omics data analysis of ovarian cancer through cancer samples identification, molecular subtypes clustering and classification, and survival analysis. The results show that MMD-VAE and VAE-based compressed features can respectively classify the transcriptional subtypes of the TCGA datasets with an accuracy in the range of 93.2-95.5% and 87.1-95.7%. Also, survival analysis results show that VAE and MMD-VAE based compressed representation of omics data can be used in cancer prognosis. Based on the results, we can conclude that (i) VAE and MMD-VAE outperform existing dimensionality reduction techniques, (ii) integrated multi-omics analyses perform better or similar compared to their mono-omics counterparts, and (iii) MMD-VAE performs better than VAE in most omics dataset.

Machine learning has recently emerged as a promising tool for inferring multi-omic relationships in biological systems. At the same time, genome-scale metabolic models (GSMMs) can be integrated with such multi-omic data to refine phenotypic predictions. In this work, we use a multi-omic machine learning pipeline to analyze a GSMM of Synechococcus sp. PCC 7002, a cyanobacterium with large potential to produce renewable biofuels. We use regularized flux balance analysis to observe flux response between conditions across photosynthesis and energy metabolism. We then incorporate principal-component analysis, k -means clustering, and LASSO regularization to reduce dimensionality and extract key cross-omic features. Our results suggest that combining metabolic modeling with machine learning elucidates mechanisms used by cyanobacteria to cope with fluctuations in light intensity and salinity that cannot be detected using transcriptomics alone. Furthermore, GSMMs introduce critical mechanistic details that improve the performance of omic-based machine learning methods., Competing Interests: The authors declare no competing interests., (© 2020 The Author(s).)

Metabolic modeling and machine learning are key components in the emerging next generation of systems and synthetic biology tools, targeting the genotype-phenotype-environment relationship. Rather than being used in isolation, it is becoming clear that their value is maximized when they are combined. However, the potential of integrating these two frameworks for omic data augmentation and integration is largely unexplored. We propose, rigorously assess, and compare machine-learning-based data integration techniques, combining gene expression profiles with computationally generated metabolic flux data to predict yeast cell growth. To this end, we create strain-specific metabolic models for 1,143 Saccharomyces cerevisiae mutants and we test 27 machine-learning methods, incorporating state-of-the-art feature selection and multiview learning approaches. We propose a multiview neural network using fluxomic and transcriptomic data, showing that the former increases the predictive accuracy of the latter and reveals functional patterns that are not directly deducible from gene expression alone. We test the proposed neural network on a further 86 strains generated in a different experiment, therefore verifying its robustness to an additional independent dataset. Finally, we show that introducing mechanistic flux features improves the predictions also for knockout strains whose genes were not modeled in the metabolic reconstruction. Our results thus demonstrate that fusing experimental cues with in silico models, based on known biochemistry, can contribute with disjoint information toward biologically informed and interpretable machine learning. Overall, this study provides tools for understanding and manipulating complex phenotypes, increasing both the prediction accuracy and the extent of discernible mechanistic biological insights., Competing Interests: The authors declare no competing interest.