By modulating the human gut microbiome, prebiotics and probiotics (combinations of which are called synbiotics) may be used to treat diseases such as colorectal cancer (CRC). Methodological limitations have prevented determining the potential combinatorial mechanisms of action of such regimens. We expanded our HuMiX gut-on-a-chip model to co-culture CRC-derived epithelial cells with a model probiotic under a simulated prebiotic regimen, and we integrated the multi-omic results with in silico metabolic modeling. In contrast to individual prebiotic or probiotic treatments, the synbiotic regimen caused downregulation of genes involved in procarcinogenic pathways and drug resistance, and reduced levels of the oncometabolite lactate. Distinct ratios of organic and short-chain fatty acids were produced during the simulated regimens. Treatment of primary CRC-derived cells with a molecular cocktail reflecting the synbiotic regimen attenuated self-renewal capacity. Our integrated approach demonstrates the potential of modeling for rationally formulating synbiotics-based treatments in the future.

Low oxygen concentrations (hypoxia) are known to affect the cellular metabolism and have been suggested to regulate a subpopulation of cancer cells with tumorigenic properties, the so-called tumor-initiating cells (TICs). To better understand the mechanism of hypoxia-induced TIC activation, we set out to study the role of hypoxia-responsive miRNAs in recently established colon cancer patient-derived TICs. We were able to show that low oxygen concentrations consistently lead to the upregulation of miR-210 in different primary TIC-enriched cultures. Both stable overexpression of miR-210 and knockdown of its target gene ISCU resulted in enhanced TIC self-renewal. We could validate the tumorigenic properties of miR- 210 in in vivo experiments by showing that ectopic expression of miR-210 results in increased tumor incidence. Furthermore, enhanced miR-210 expression correlated with reduced TCA cycle activity and increased lactate levels. Importantly, by blocking lactate production via inhibition of LDHA, we could reverse the promoting effect of miR-210 on self-renewal capacity, thereby emphasizing the regulatory impact of the glycolytic phenotype on colon TIC properties. Finally, by assessing expression levels in patient tissue, we could demonstrate the clinical relevance of the miR-210/ISCU signaling axis for colorectal carcinoma. Taken together, our study highlights the importance of hypoxia-induced miR-210 in the regulation of colon cancer initiation.

The mechanism of each chemical reaction in a metabolic network can be represented as a set of atom mappings, each of which relates an atom in a substrate metabolite to an atom of the same element in a product metabolite. Genome-scale metabolic network reconstructions typically represent biochemistry at the level of reaction stoichiometry. However, a more detailed representation at the underlying level of atom mappings opens the possibility for a broader range of biological, biomedical and biotechnological applications than with stoichiometry alone. Complete manual acquisition of atom mapping data for a genome-scale metabolic network is a laborious process. However, many algorithms exist to predict atom mappings. How do their predictions compare to each other and to manually curated atom mappings? For more than four thousand metabolic reactions in the latest human metabolic reconstruction, Recon 3D, we compared the atom mappings predicted by six atom mapping algorithms. We also compared these predictions to those obtained by manual curation of atom mappings for over five hundred reactions distributed among all top level Enzyme Commission number classes. Five of the evaluated algorithms had similarly high prediction accuracy of over 91% when compared to manually curated atom mapped reactions. On average, the accuracy of the prediction was highest for reactions catalysed by oxidoreductases and lowest for reactions catalysed by ligases. In addition to prediction accuracy, the algorithms were evaluated on their accessibility, their advanced features, such as the ability to identify equivalent atoms, and their ability to map hydrogen atoms. In addition to prediction accuracy, we found that software accessibility and advanced features were fundamental to the selection of an atom mapping algorithm in practice.

Cardiomyopathy is a common complication associated with increased mortality in sepsis, but lacks specific therapy. Here, using genetic and pharmacological approaches, we explored the therapeutic effect of α2A-adrenergic receptor (AR) blockade on septic cardiomyopathy. CLP-induced septic rats were treated with BRL44408 (α2A-AR antagonist), prazosin (α1-AR antagonist) and/or reserpine. CLP-induced cardiomyopathy, indicated by reduced dP/dt and increased cardiac troponin I phosphorylation, was attenuated by BRL44408, this was associated with reduced cardiac TNF-α and endothelial VCAM-1 expression, cardiomyocyte apoptosis and related signal molecule phosphorylation. BRL44408 increased cardiac norepinephrine (NE) concentration in CLP rats. Pretreatment with reserpine that exhausts cardiac NE without affecting the circulating NE concentration or with prazosin partially abolished the cardioprotection of BRL44408 and reversed its inhibitory effects on myocardial TNF-α, apoptosis and related signal molecule phosphorylation, but not on VCAM-1 expression in septic rats. These effects of BRL44408 were confirmed by α2A-AR gene deletion in septic mice. Furthermore, α2-AR agonist not only enhanced LPS-induced TNF-α and VCAM-1 expression in cardiac endothelial cells that express α2A-AR, but also enhanced LPS-induced cardiac dysfunction in isolated rat hearts. Our data indicate that α2A-AR blockade attenuates septic cardiomyopathy by promoting cardiac NE release that activates myocardial α1-AR and suppressing cardiac endothelial activation.

Most melanoma patients with BRAFV600E positive tumors respond well to a combination of BRAF kinase and MEK inhibitors. However, some patients are intrinsically resistant while the majority of patients eventually develop drug resistance to the treatment. For patients insufficiently responding to BRAF and MEK inhibitors, there is an ongoing need for new treatment targets. Cellular metabolism is such a promising new target line: mutant BRAFV600E has been shown to affect the metabolism.

Development of heart failure (HF) after myocardial infarction (MI) is responsible for premature death. Complex cellular and molecular mechanisms are involved in this process. A number of studies have linked the epitranscriptomic RNA modification N6-methyladenosine (m6A) with HF, but it remains unknown how m6A affects the risk of developing HF after MI. We addressed the regulation of m6A and its demethylase fat mass and obesity-associated (FTO) after MI and their association with HF. Using liquid chromatography coupled to mass spectrometry, we observed an increase of m6A content in the infarcted area of rat hearts subjected to coronary ligation and a decrease in blood. FTO expression measured by quantitative PCR was downregulated in the infarcted hearts. In whole blood samples collected at the time of reperfusion in MI patients, m6A content was lower in patients who developed HF as attested by a 4-month ejection fraction (EF) of ≤40% as compared to patients who did not develop HF (EF > 50%). M6A content was higher in females. These results show that m6A measured in blood is associated with HF development after MI and motivate further investigation of the potential role of m6A as a novel epitranscriptomics biomarker and therapeutic target of HF.

Mutations in ATP13A2 (PARK9) are causally linked to the rare neurodegenerative disorders Kufor-Rakeb syndrome, hereditary spastic paraplegia and neuronal ceroid lipofuscinosis. This suggests that ATP13A2, a lysosomal cation-transporting ATPase, plays a crucial role in neuronal cells. The heterogeneity of the clinical spectrum of ATP13A2-associated disorders is not yet well understood and currently, these diseases remain without effective treatment. Interestingly, ATP13A2 is widely conserved among eukaryotes, and the yeast model for ATP13A2 deficiency was the first to indicate a role in heavy metal homeostasis, which was later confirmed in human cells. In this study, we show that the deletion of YPK9 (the yeast orthologue of ATP13A2) in Saccharomyces cerevisiae leads to growth impairment in the presence of Zn2+, Mn2+, Co2+ and Ni2+, with the strongest phenotype being observed in the presence of zinc. Using the ypk9Δ mutant, we developed a high-throughput growth rescue screen based on the Zn2+ sensitivity phenotype. Screening of two libraries of Food and Drug Administration-approved drugs identified 11 compounds that rescued growth. Subsequently, we generated a zebrafish model for ATP13A2 deficiency and found that both partial and complete loss of atp13a2 function led to increased sensitivity to Mn2+. Based on this phenotype, we confirmed two of the drugs found in the yeast screen to also exert a rescue effect in zebrafish-N-acetylcysteine, a potent antioxidant, and furaltadone, a nitrofuran antibiotic. This study further supports that combining the high-throughput screening capacity of yeast with rapid in vivo drug testing in zebrafish can represent an efficient drug repurposing strategy in the context of rare inherited disorders involving conserved genes. This work also deepens the understanding of the role of ATP13A2 in heavy metal detoxification and provides a new in vivo model for investigating ATP13A2 deficiency.

Tumours are composed of various cancer cell populations with different mutation profiles, phenotypes and metabolism that cause them to react to drugs in diverse manners. Increasing the resolution of metabolic models based on single-cell expression data will provide deeper insight into such metabolic differences and improve the predictive power of the models. scFASTCORMICS is a network contextualization algorithm that builds multi-cell population genome-scale models from single-cell RNAseq data. The models contain a subnetwork for each cell population in a tumour, allowing to capture metabolic variations between these clusters. The subnetworks are connected by a union compartment that permits to simulate metabolite exchanges between cell populations in the microenvironment. scFASTCORMICS uses Pareto optimization to simultaneously maximise the compactness, completeness and specificity of the reconstructed metabolic models. scFASTCORMICS is implemented in MATLAB and requires the installation of the COBRA toolbox, rFASTCORMICS and the IBM CPLEX solver.

Parkinson's disease (PD)-specific neurons, grown in standard 2D cultures, typically only display weak endophenotypes. The cultivation of PD patient-specific neurons, derived from induced pluripotent stem cells carrying the LRRK2-G2019S mutation, is optimized in 3D microfluidics. The automated image analysis algorithms are implemented to enable pharmacophenomics in disease-relevant conditions. In contrast to 2D cultures, this 3D approach reveals robust endophenotypes. High-content imaging data show decreased dopaminergic differentiation and branching complexity, altered mitochondrial morphology, and increased cell death in LRRK2-G2019S neurons compared to isogenic lines without using stressor agents. Treatment with the LRRK2 inhibitor 2 (Inh2) rescues LRRK2-G2019S-dependent dopaminergic phenotypes. Strikingly, a holistic analysis of all studied features shows that the genetic background of the PD patients, and not the LRRK2-G2019S mutation, constitutes the strongest contribution to the phenotypes. These data support the use of advanced in vitro models for future patient stratification and personalized drug development.

Curcumin, a natural polyphenol, has many biological properties, such as anti-inflammatory, antioxidant, and anti-carcinogenic properties, yet, its sensitivity to light, oxygen, and heat, and its low solubility in water renders its preservation and bioavailability challenging. To increase its bioaccessibility, we fabricated nanoliposomes and chitosan-coated nanoliposomes encapsulating curcumin, and we evaluated the systems in terms of their physicochemical characteristics and release profiles in simulated gastrointestinal mediums. Chitosan-coating enhanced the stability of nanoliposomes and slowed the release of curcumin in the simulated gastrointestinal (GI) environment. This study demonstrates that nanoliposomes and chitosan-coated nanoliposomes are promising carriers for poorly soluble lipophilic compounds with low oral bioavailability, such as curcumin.

The new European legislation on data protection, namely, the General Data Protection Regulation (GDPR), has introduced comprehensive requirements for the documentation about the processing of personal data as well as informing the data subjects of its use. GDPR's accountability principle requires institutions, projects, and data hubs to document their data processings and demonstrate compliance with the GDPR. In response to this requirement, we see the emergence of commercial data-mapping tools, and institutions creating GDPR data register with such tools. One shortcoming of this approach is the genericity of tools, and their process-based model not capturing the project-based, collaborative nature of data processing in biomedical research.

A multitude of factors contribute to complex diseases and can be measured with 'omics' methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources 'Human metabolism', 'Gut microbiome', 'Disease', 'Nutrition', and 'ReconMaps'. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH's unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

Genome-scale metabolic reconstructions include all known biochemical reactions occurring in a cell. A typical application is the prediction of potential drug targets for cancer treatment. The precision of these predictions relies on the definition of the objective function. Generally, the biomass reaction is used to illustrate the growth capacity of a cancer cell. Today, seven human biomass reactions can be identified in published metabolic models. The impact of these differences on the metabolic model predictions has not been explored in detail. We explored this impact on cancer metabolic model predictions and showed that the metabolite composition and the associated coefficients had a large impact on the growth rate prediction accuracy, whereas gene essentiality predictions were mainly affected by the metabolite composition. Our results demonstrate the importance of defining a consensus biomass reaction compatible with most human models, which would contribute to ensuring the reproducibility and consistency of the results.

Human brain organoid models that recapitulate the physiology and complexity of the human brain have a great potential for in vitro disease modeling, in particular for neurodegenerative diseases, such as Parkinson disease. In the present study, we compare single-cell RNA-sequencing data of human midbrain organoids to the developing human embryonic midbrain. We demonstrate that the in vitro model is comparable to its in vivo equivalents in terms of developmental path and cellular composition. Moreover, we investigate the potential of midbrain organoids for modeling early developmental changes in Parkinson disease. Therefore, we compare the single-cell RNA-sequencing data of healthy-individual-derived midbrain organoids to their isogenic LRRK2-p.Gly2019Ser-mutant counterparts. We show that the LRRK2 p.Gly2019Ser variant alters neurodevelopment, resulting in an untimely and incomplete differentiation with reduced cellular variability. Finally, we present four candidate genes, APP, DNAJC6, GATA3, and PTN, that might contribute to the LRRK2-p.Gly2019Ser-associated transcriptome changes that occur during early neurodevelopment.

S100A9, a Ca2+-binding protein, is tightly associated to neutrophil pro-inflammatory functions when forming a heterodimer with its S100A8 partner. Upon secretion into the extracellular environment, these proteins behave like damage-associated molecular pattern molecules, which actively participate in the amplification of the inflammation process by recruitment and activation of pro-inflammatory cells. Intracellular functions have also been attributed to the S100A8/A9 complex, notably its ability to regulate nicotinamide adenine dinucleotide phosphate (NADPH) oxidase activation. However, the complete functional spectrum of S100A8/A9 at the intracellular level is far from being understood. In this context, we here investigated the possibility that the absence of intracellular S100A8/A9 is involved in cytokine secretion. To overcome the difficulty of genetically modifying neutrophils, we used murine neutrophils derived from wild-type and S100A9-/- Hoxb8 immortalized myeloid progenitors. After confirming that differentiated Hoxb8 neutrophil-like cells are a suitable model to study neutrophil functions, our data show that absence of S100A8/A9 led to a dysregulation of cytokine secretion after lipopolysaccharide (LPS) stimulation. Furthermore, we demonstrate that S100A8/A9-induced cytokine secretion was regulated by the nuclear factor kappa B (NF-κB) pathway. These results were confirmed in human differentiated HL-60 cells, in which S100A9 was inhibited by shRNAs. Finally, our results indicate that the degranulation process could be involved in the regulation of cytokine secretion by S100A8/A9.

Although most cases of Parkinson´s disease (PD) are idiopathic with unknown cause, an increasing number of genes and genetic risk factors have been discovered that play a role in PD pathogenesis. Many of the PD-associated proteins are involved in mitochondrial quality control, e.g., PINK1, Parkin, and LRRK2, which were recently identified as regulators of mitochondrial-endoplasmic reticulum (ER) contact sites (MERCs) linking mitochondrial homeostasis to intracellular calcium handling. In this context, Miro1 is increasingly recognized to play a role in PD pathology. Recently, we identified the first PD patients carrying mutations in RHOT1, the gene coding for Miro1. Here, we describe two novel RHOT1 mutations identified in two PD patients and the characterization of the cellular phenotypes.

Changes in the human gastrointestinal microbiome are associated with several diseases. To infer causality, experiments in representative models are essential, but widely used animal models exhibit limitations. Here we present a modular, microfluidics-based model (HuMiX, human-microbial crosstalk), which allows co-culture of human and microbial cells under conditions representative of the gastrointestinal human-microbe interface. We demonstrate the ability of HuMiX to recapitulate in vivo transcriptional, metabolic and immunological responses in human intestinal epithelial cells following their co-culture with the commensal Lactobacillus rhamnosus GG (LGG) grown under anaerobic conditions. In addition, we show that the co-culture of human epithelial cells with the obligate anaerobe Bacteroides caccae and LGG results in a transcriptional response, which is distinct from that of a co-culture solely comprising LGG. HuMiX facilitates investigations of host-microbe molecular interactions and provides insights into a range of fundamental research questions linking the gastrointestinal microbiome to human health and disease.

Reduced levels of intratumoural oxygen are associated with hypoxia-induced pro-oncogenic events such as invasion, metabolic reprogramming, epithelial-mesenchymal transition, metastasis and resistance to therapy, all favouring cancer progression. Small extracellular vesicles (EV) shuttle various cargos (proteins, miRNAs, DNA and others). Tumour-derived EVs can be taken up by neighbouring or distant cells in the tumour microenvironment, thus facilitating intercellular communication. The quantity of extracellular vesicle secretion and their composition can vary with changing microenvironmental conditions and disease states. Here, we investigated in melanoma cells the influence of hypoxia on the content and number of secreted EVs. Whole miRNome and proteome profiling revealed distinct expression patterns in normoxic or hypoxic growth conditions. Apart from the well-known miR-210, we identified miR-1290 as a novel hypoxia-associated microRNA, which was highly abundant in hypoxic EVs. On the other hand, miR-23a-5p and -23b-5p were consistently downregulated in hypoxic conditions, while the protein levels of the miR-23a/b-5p-predicted target IPO11 were concomitantly upregulated. Furthermore, hypoxic melanoma EVs exhibit a signature consisting of six proteins (AKR7A2, DDX39B, EIF3C, FARSA, PRMT5, VARS), which were significantly associated with a poor prognosis for melanoma patients, indicating that proteins and/or miRNAs secreted by cancer cells may be exploited as biomarkers.

Parkinson's disease (PD) is the second most common neurodegenerative disorder. Accumulating evidences suggest that PD might have a strong neurodevelopmental component. Among the genetic cases, mutations in the leucine-rich repeat kinase 2 (LRRK2) are well known to be disease causing. Although the molecular mechanism of the pathogenic LRRK2 function is not fully clear, inhibition of microRNA (miRNA) activity has been suggested to be among the pathogenic LRRK2 targets. Here, we demonstrate that the miRNA activity inhibition function of pathogenic LRRK2 is directly antagonized by the neuronal cell fate determinant TRIM32. These findings suggest that TRIM32 might be a modifier for PD and could be a novel therapeutic target.

Melanoma is the most aggressive and deadly form of skin cancer with increasing case numbers worldwide. The development of inhibitors targeting mutated BRAF (found in around 60% of melanoma patients) has markedly improved overall survival of patients with late-stage tumors, even more so when combined with MEK inhibitors targeting the same signaling pathway. However, invariably patients become resistant to this targeted therapy resulting in rapid progression with treatment-refractory disease. The purpose of this study was the identification of new kinase inhibitors that do not lead to the development of resistance in combination with BRAF inhibitors (BRAFi), or that could be of clinical benefit as a 2nd line treatment for late-stage melanoma patients that have already developed resistance.