Mechanical effects modulate drug resistance in MCF-7-derived organoids: Insights into the wnt/β-catenin pathway

Addressing drug resistance poses a significant challenge in cancer treatment, as cancer cells develop diverse mechanisms to evade chemotherapy drugs, leading to treatment failure and disease relapse. Three-dimensional (3D) cell culture has emerged as a valuable model for studying drug resistance, although the underlying mechanisms remain elusive. By obtaining a better understanding of drug resistance within the 3D culture environment, we can develop more effective strategies to overcome it and improve the success of cancer treatments. Notably, the physical structure undergoes notable changes in 3D culture, with mechanical effects believed to play a pivotal role in drug resistance. Hence, our study aimed to explore the influence of mechanical effects on drug resistance by analyzing data related to "drug resistance" and "mechanobiology". Through this analysis, we identified β-catenin and JNK1 as potential factors, which were further examined in MCF-7 cells cultivated under both 2D and 3D culture conditions. Our findings demonstrate that β-catenin is activated through canonical and non-canonical pathways and associated with the drug resistance, particularly in organoids obtained under 3D culture.

Metabolic function-based normalization improves transcriptome data-driven reduction of genome-scale metabolic models

Genome-scale metabolic models (GEMs) are extensively used to simulate cell metabolism and predict cell phenotypes. GEMs can also be tailored to generate context-specific GEMs, using omics data integration approaches. To date, many integration approaches have been developed, however, each with specific pros and cons; and none of these algorithms systematically outperforms the others. The key to successful implementation of such integration algorithms lies in the optimal selection of parameters, and thresholding is a crucial component in this process. To improve the predictive accuracy of context-specific models, we introduce a new integration framework that improves the ranking of related genes and homogenizes the expression values of those gene sets using single-sample Gene Set Enrichment Analysis (ssGSEA). In this study, we coupled ssGSEA with GIMME and validated the advantages of the proposed framework to predict the ethanol formation of yeast grown in the glucose-limited chemostats, and to simulate metabolic behaviors of yeast growth in four different carbon sources. This framework enhances the predictive accuracy of GIMME which we demonstrate for predicting the yeast physiology in nutrient-limited cultures.

Construction of a three-component regulatory network of transcribed ultraconserved regions for the identification of prognostic biomarkers in gastric cancer

Altered expression and functional roles of the transcribed ultraconserved regions (T-UCRs), as genomic sequences with 100% conservation between the genomes of human, mouse, and rat, in the pathophysiology of neoplasms has already been investigated. Nevertheless, the relevance of the functions for T-UCRs in gastric cancer (GC) is still the subject of inquiry. In the current study, we first used a genome-wide profiling approach to analyze the expression of T-UCRs in GC patients. Then, we constructed a three-component regulatory network and investigated potential diagnostic and prognostic values of the T-UCRs. The Cancer Genome Atlas Stomach Adenocarcinoma (TCGA-STAD) dataset was used as a resource for the RNA-sequencing data. FeatureCounts was utilized to quantify the number of reads mapped to each T-UCR. Differential expression analysis was then conducted using DESeq2. In the following, interactions between T-UCRs, microRNAs (miRNAs), and messenger RNAs (mRNAs) were combined into a three-component network. Enrichment analyses were performed and a protein-protein interaction (PPI) network was constructed. The R Survival package was utilized to identify survival-related significantly differentially expressed T-UCRs (DET-UCRs). Using an in-house cohort of GC tissues, expression of two DET-UCRs was furthermore experimentally verified. Our results showed that several T-UCRs were dysregulated in TCGA-STAD tumoral samples compared to nontumoral counterparts. The three-component network was constructed which composed of DET-UCRs, miRNAs, and mRNAs nodes. Functional enrichment and PPI network analyses revealed important enriched signaling pathways and gene ontologies such as "pathway in cancer" and regulation of cell proliferation and apoptosis. Five T-UCRs were significantly correlated with the overall survival of GC patients. While no expression of uc.232 was observed in our in-house cohort of GC tissues, uc.343 showed an increased expression, although not statistically significant, in gastric tumoral tissues. The constructed three-component regulatory network of T-UCRs in GC presents a comprehensive understanding of the underlying gene expression regulation processes involved in tumor development and can serve as a basis to investigate potential prognostic biomarkers and therapeutic targets.

Computational Models for Clinical Applications in Personalized Medicine—Guidelines and Recommendations for Data Integration and Model Validation

The future development of personalized medicine depends on a vast exchange of data from different sources, as well as harmonized integrative analysis of large-scale clinical health and sample data. Computational-modelling approaches play a key role in the analysis of the underlying molecular processes and pathways that characterize human biology, but they also lead to a more profound understanding of the mechanisms and factors that drive diseases; hence, they allow personalized treatment strategies that are guided by central clinical questions. However, despite the growing popularity of computational-modelling approaches in different stakeholder communities, there are still many hurdles to overcome for their clinical routine implementation in the future. Especially the integration of heterogeneous data from multiple sources and types are challenging tasks that require clear guidelines that also have to comply with high ethical and legal standards. Here, we discuss the most relevant computational models for personalized medicine in detail that can be considered as best-practice guidelines for application in clinical care. We define specific challenges and provide applicable guidelines and recommendations for study design, data acquisition, and operation as well as for model validation and clinical translation and other research areas.

A Metabolic Model of Intestinal Secretions: The Link between Human Microbiota and Colorectal Cancer Progression

The human gut microbiota plays a dual key role in maintaining human health or inducing disorders, for example, obesity, type 2 diabetes, and cancers such as colorectal cancer (CRC). High-throughput data analysis, such as metagenomics and metabolomics, have shown the diverse effects of alterations in dynamic bacterial populations on the initiation and progression of colorectal cancer. However, it is well established that microbiome and human cells constantly influence each other, so it is not appropriate to study them independently. Genome-scale metabolic modeling is a well-established mathematical framework that describes the dynamic behavior of these two axes at the system level. In this study, we created community microbiome models of three conditions during colorectal cancer progression, including carcinoma, adenoma and health status, and showed how changes in the microbial population influence intestinal secretions. Conclusively, our findings showed that alterations in the gut microbiome might provoke mutations and transform adenomas into carcinomas. These alterations include the secretion of mutagenic metabolites such as H₂S, NO compounds, spermidine and TMA (trimethylamine), as well as the reduction of butyrate. Furthermore, we found that the colorectal cancer microbiome can promote inflammation, cancer progression (e.g., angiogenesis) and cancer prevention (e.g., apoptosis) by increasing and decreasing certain metabolites such as histamine, glutamine and pyruvate. Thus, modulating the gut microbiome could be a promising strategy for the prevention and treatment of CRC.

Towards biomarkers for outcomes after pancreatic ductal adenocarcinoma and ischaemic stroke, with focus on (co)-morbidity and ageing/cellular senescence (SASKit): protocol for a prospective cohort study

INTRODUCTION

Ageing-related processes such as cellular senescence are believed to underlie the accumulation of diseases in time, causing (co)morbidity, including cancer, thromboembolism and stroke. Interfering with these processes may delay, stop or reverse morbidity. The aim of this study is to investigate the link between (co)morbidity and ageing by exploring biomarkers and molecular mechanisms of disease-triggered deterioration in patients with pancreatic ductal adenocarcinoma (PDAC) and (thromboembolic) ischaemic stroke (IS).

METHODS AND ANALYSIS

We will recruit 50 patients with PDAC, 50 patients with (thromboembolic) IS and 50 controls at Rostock University Medical Center, Germany. We will gather routine blood data, clinical performance measurements and patient-reported outcomes at up to seven points in time, alongside in-depth transcriptomics and proteomics at two of the early time points. Aiming for clinically relevant biomarkers, the primary outcome is a composite of probable sarcopenia, clinical performance (described by ECOG Performance Status for patients with PDAC and the Modified Rankin Scale for patients with stroke) and quality of life. Further outcomes cover other aspects of morbidity such as cognitive decline and of comorbidity such as vascular or cancerous events. The data analysis is comprehensive in that it includes biostatistics and machine learning, both following standard role models and additional explorative approaches. Prognostic and predictive biomarkers for interventions addressing senescence may become available if the biomarkers that we find are specifically related to ageing/cellular senescence. Similarly, diagnostic biomarkers will be explored. Our findings will require validation in independent studies, and our dataset shall be useful to validate the findings of other studies. In some of the explorative analyses, we shall include insights from systems biology modelling as well as insights from preclinical animal models. We anticipate that our detailed study protocol and data analysis plan may also guide other biomarker exploration trials.

ETHICS AND DISSEMINATION

The study was approved by the local ethics committee (Ethikkommission an der Medizinischen Fakultät der Universität Rostock, A2019-0174), registered at the German Clinical Trials Register (DRKS00021184), and results will be published following standard guidelines.

GEMtractor: extracting views into genome-scale metabolic models

SUMMARY

Computational metabolic models typically encode for graphs of species, reactions and enzymes. Comparing genome-scale models through topological analysis of multipartite graphs is challenging. However, in many practical cases it is not necessary to compare the full networks. The GEMtractor is a web-based tool to trim models encoded in SBML. It can be used to extract subnetworks, for example focusing on reaction- and enzyme-centric views into the model.

AVAILABILITY AND IMPLEMENTATION

The GEMtractor is licensed under the terms of GPLv3 and developed at github.com/binfalse/GEMtractor-a public version is available at sbi.uni-rostock.de/gemtractor.

Enzyme-Constrained Models and Omics Analysis of Streptomyces coelicolor Reveal Metabolic Changes that Enhance Heterologous Production

Many biosynthetic gene clusters (BGCs) require heterologous expression to realize their genetic potential, including silent and metagenomic BGCs. Although the engineered Streptomyces coelicolor M1152 is a widely used host for heterologous expression of BGCs, a systemic understanding of how its genetic modifications affect the metabolism is lacking and limiting further development. We performed a comparative analysis of M1152 and its ancestor M145, connecting information from proteomics, transcriptomics, and cultivation data into a comprehensive picture of the metabolic differences between these strains. Instrumental to this comparison was the application of an improved consensus genome-scale metabolic model (GEM) of S. coelicolor. Although many metabolic patterns are retained in M1152, we find that this strain suffers from oxidative stress, possibly caused by increased oxidative metabolism. Furthermore, precursor availability is likely not limiting polyketide production, implying that other strategies could be beneficial for further development of S. coelicolor for heterologous production of novel compounds.

Evolution of computational models in BioModels Database and the Physiome Model Repository

BACKGROUND

A useful model is one that is being (re)used. The development of a successful model does not finish with its publication. During reuse, models are being modified, i.e. expanded, corrected, and refined. Even small changes in the encoding of a model can, however, significantly affect its interpretation. Our motivation for the present study is to identify changes in models and make them transparent and traceable.

METHODS

We analysed 13734 models from BioModels Database and the Physiome Model Repository. For each model, we studied the frequencies and types of updates between its first and latest release. To demonstrate the impact of changes, we explored the history of a Repressilator model in BioModels Database.

RESULTS

We observed continuous updates in the majority of models. Surprisingly, even the early models are still being modified. We furthermore detected that many updates target annotations, which improves the information one can gain from models. To support the analysis of changes in model repositories we developed MoSt, an online tool for visualisations of changes in models. The scripts used to generate the data and figures for this study are available from GitHub https://github.com/binfalse/BiVeS-StatsGenerator and as a Docker image at https://hub.docker.com/r/binfalse/bives-statsgenerator/ . The website https://most.bio.informatik.uni-rostock.de/ provides interactive access to model versions and their evolutionary statistics.

CONCLUSION

The reuse of models is still impeded by a lack of trust and documentation. A detailed and transparent documentation of all aspects of the model, including its provenance, will improve this situation. Knowledge about a model's provenance can avoid the repetition of mistakes that others already faced. More insights are gained into how the system evolves from initial findings to a profound understanding. We argue that it is the responsibility of the maintainers of model repositories to offer transparent model provenance to their users.

Unveiling network-based functional features through integration of gene expression into protein networks

Decoding health and disease phenotypes is one of the fundamental objectives in biomedicine. Whereas high-throughput omics approaches are available, it is evident that any single omics approach might not be adequate to capture the complexity of phenotypes. Therefore, integrated multi-omics approaches have been used to unravel genotype-phenotype relationships such as global regulatory mechanisms and complex metabolic networks in different eukaryotic organisms. Some of the progress and challenges associated with integrated omics studies have been reviewed previously in comprehensive studies. In this work, we highlight and review the progress, challenges and advantages associated with emerging approaches, integrating gene expression and protein-protein interaction networks to unravel network-based functional features. This includes identifying disease related genes, gene prioritization, clustering protein interactions, developing the modules, extract active subnetworks and static protein complexes or dynamic/temporal protein complexes. We also discuss how these approaches contribute to our understanding of the biology of complex traits and diseases. This article is part of a Special Issue entitled: Cardiac adaptations to obesity, diabetes and insulin resistance, edited by Professors Jan F.C. Glatz, Jason R.B. Dyck and Christine Des Rosiers.

Comprehensive functional enrichment analysis of male infertility

Spermatogenesis is a multifactorial process that forms differentiated sperm cells in a complex microenvironment. This process involves the genome, epigenome, transcriptome, and proteome to ensure the stability of the spermatogonia and supporting cells. The identification of signaling pathways linked to infertility has been hampered by the inherent complexity and multifactorial aspects of spermatogenesis. Systems biology is a promising approach to unveil underlying signaling pathways and genes and identify putative biomarkers. In this study, we analyzed thirteen microarray libraries of infertile humans and mice, and different classes of male infertility were compared using differentially expressed genes and functional enrichment analysis. We found regulatory processes, immune response, glutathione transferase and muscle tissue development to be among the most common biological processes in up-regulated genes, and genes involved in spermatogenesis were down-regulated in maturation arrest (MArrest) and oligospermia cases. We also observed the overexpression of genes involved in steroid metabolism in post-meiotic and meiotic arrest. Furthermore, we found that the infertile mouse model most similar to human MArrest was the Dazap1 mutant mouse. The results of this study could help elucidate features of infertility etiology and provide the basis for diagnostic markers.

Meta-Analysis of Gene Expression Profiles in Acute Promyelocytic Leukemia Reveals Involved Pathways

Background: Acute promyelocytic leukemia (APL) is a unique subtype of acute leukemia. APL is a curable disease; however, drug resistance, early mortality, disease relapse and treatment-related complications remain challenges in APL patient management. One issue underlying these challenges is that the molecular mechanisms of the disease are not sufficiently understood. Materials and Methods: In this study, we performed a meta-analysis of gene expression profiles derived from microarray experiments and explored the background of disease by functional and pathway analysis. Results: Our analysis revealed a gene signature with 406 genes that are up or down-regulated in APL. The pathway analysis determined that MAPK pathway and its involved elements such as JUN gene and AP-1 play important roles in APL pathogenesis along with insulin-like growth factor-binding protein-7. Conclusion: The results of this meta-analysis could be useful for developing more effective therapy strategies and new targets for diagnosis and drugs.

Cancerome: A hidden informative subnetwork of the diseasome

Neoplastic disorders are a leading cause of mortality and morbidity worldwide. Studying the relationships between different cancers using high throughput-generated data may elucidate undisclosed aspects of cancer etiology, diagnosis, and treatment. Several studies have described relationships between different diseases based on genes, proteins, pathways, gene ontology, comorbidity, symptoms, and other features. In this study, we first constructed an integrated human disease network based on nine different biological aspects, including molecular, functional, and clinical features. Next, we extracted the cancerome as a cancer-related subnetwork. Further investigation of cancerome could reveal hidden mechanisms of cancer and could be useful in developing new diagnostic tests and effective new drugs.

Impairment of DNA Methylation Maintenance Is the Main Cause of Global Demethylation in Naive Embryonic Stem Cells

Global demethylation is part of a conserved program of epigenetic reprogramming to naive pluripotency. The transition from primed hypermethylated embryonic stem cells (ESCs) to naive hypomethylated ones (serum-to-2i) is a valuable model system for epigenetic reprogramming. We present a mathematical model, which accurately predicts global DNA demethylation kinetics. Experimentally, we show that the main drivers of global demethylation are neither active mechanisms (Aicda, Tdg, and Tet1-3) nor the reduction of de novo methylation. UHRF1 protein, the essential targeting factor for DNMT1, is reduced upon transition to 2i, and so is recruitment of the maintenance methylation machinery to replication foci. Concurrently, there is global loss of H3K9me2, which is needed for chromatin binding of UHRF1. These mechanisms synergistically enforce global DNA hypomethylation in a replication-coupled fashion. Our observations establish the molecular mechanism for global demethylation in naive ESCs, which has key parallels with those operating in primordial germ cells and early embryos.

•

Impaired DNA methylation maintenance is the cause of global demethylation in naive ESCs

•

Loss of H3K9me2 and UHRF1 lead to impaired maintenance targeting to replication foci

•

TET enzymes are not required for global demethylation

•

Mathematical model accurately predicts global 5mC and 5hmC during epigenetic resetting

CentiServer: A Comprehensive Resource, Web-Based Application and R Package for Centrality Analysis

Various disciplines are trying to solve one of the most noteworthy queries and broadly used concepts in biology, essentiality. Centrality is a primary index and a promising method for identifying essential nodes, particularly in biological networks. The newly created CentiServer is a comprehensive online resource that provides over 110 definitions of different centrality indices, their computational methods, and algorithms in the form of an encyclopedia. In addition, CentiServer allows users to calculate 55 centralities with the help of an interactive web-based application tool and provides a numerical result as a comma separated value (csv) file format or a mapped graphical format as a graph modeling language (GML) file. The standalone version of this application has been developed in the form of an R package. The web-based application (CentiServer) and R package (centiserve) are freely available at http://www.centiserver.org/

Alterations in cancer cell metabolism: the Warburg effect and metabolic adaptation

The Warburg effect means higher glucose uptake of cancer cells compared to normal tissues, whereas a smaller fraction of this glucose is employed for oxidative phosphorylation. With the advent of high throughput technologies and computational systems biology, cancer cell metabolism has been reinvestigated over the last decades toward identifying various events underlying "how" and "why" a cancer cell employs aerobic glycolysis. Significant progress has been shaped to revise the Warburg effect. In this study, we have integrated the gene expression of 13 different cancer cells with the genome-scale metabolic network of human (Recon1) based on the E-Flux method, and analyzed them based on constraint-based modeling. Results show that regardless of significant up- and down-regulated metabolic genes, the distribution of metabolic changes is similar in different cancer types. These findings support the theory that the Warburg effect is a consequence of metabolic adaptation in cancer cells.

Computational analysis of reciprocal association of metabolism and epigenetics in the budding yeast: a genome-scale metabolic model (GSMM) approach

Metaboloepigenetics is a newly coined term in biological sciences that investigates the crosstalk between epigenetic modifications and metabolism. The reciprocal relation between biochemical transformations and gene expression regulation has been experimentally demonstrated in cancers and metabolic syndromes. In this study, we explored the metabolism-histone modifications crosstalk by topological analysis and constraint-based modeling approaches in the budding yeast. We constructed nine models through the integration of gene expression data of four mutated histone tails into a genome-scale metabolic model of yeast. Accordingly, we defined the centrality indices of the lowly expressed enzymes in the undirected enzyme-centric network of yeast by CytoHubba plug-in in Cytoscape. To determine the global effects of histone modifications on the yeast metabolism, the growth rate and the range of possible flux values of reactions, we used constraint-based modeling approach. Centrality analysis shows that the lowly expressed enzymes could affect and control the yeast metabolic network. Besides, constraint-based modeling results are in a good agreement with the experimental findings, confirming that the mutations in histone tails lead to non-lethal alterations in the yeast, but have diverse effects on the growth rate and reveal the functional redundancy.

HomoTarget: a new algorithm for prediction of microRNA targets in Homo sapiens

MiRNAs play an essential role in the networks of gene regulation by inhibiting the translation of target mRNAs. Several computational approaches have been proposed for the prediction of miRNA target-genes. Reports reveal a large fraction of under-predicted or falsely predicted target genes. Thus, there is an imperative need to develop a computational method by which the target mRNAs of existing miRNAs can be correctly identified. In this study, combined pattern recognition neural network (PRNN) and principle component analysis (PCA) architecture has been proposed in order to model the complicated relationship between miRNAs and their target mRNAs in humans. The results of several types of intelligent classifiers and our proposed model were compared, showing that our algorithm outperformed them with higher sensitivity and specificity. Using the recent release of the mirBase database to find potential targets of miRNAs, this model incorporated twelve structural, thermodynamic and positional features of miRNA:mRNA binding sites to select target candidates.