Host-Microbiome Interaction and Cancer: Potential Application in Precision Medicine

It has been experimentally shown that host-microbial interaction plays a major role in shaping the wellness or disease of the human body. Microorganisms coexisting in human tissues provide a variety of benefits that contribute to proper functional activity in the host through the modulation of fundamental processes such as signal transduction, immunity and metabolism. The unbalance of this microbial profile, or dysbiosis, has been correlated with the genesis and evolution of complex diseases such as cancer. Although this latter disease has been thoroughly studied using different high-throughput technologies, its heterogeneous nature makes its understanding and proper treatment in patients a remaining challenge in clinical settings. Notably, given the outstanding role of host-microbiome interactions, the ecological interactions with microorganisms have become a new significant aspect in the systems that can contribute to the diagnosis and potential treatment of solid cancers. As a part of expanding precision medicine in the area of cancer research, efforts aimed at effective treatments for various kinds of cancer based on the knowledge of genetics, biology of the disease and host-microbiome interactions might improve the prediction of disease risk and implement potential microbiota-directed therapeutics. In this review, we present the state of the art of sequencing and metabolome technologies, computational methods and schemes in systems biology that have addressed recent breakthroughs of uncovering relationships or associations between microorganisms and cancer. Together, microbiome studies extend the horizon of new personalized treatments against cancer from the perspective of precision medicine through a synergistic strategy integrating clinical knowledge, high-throughput data, bioinformatics and systems biology

MICOM: Metagenome-Scale Modeling To Infer Metabolic Interactions in the Gut Microbiota.

Compositional changes in the gut microbiota have been associated with a variety of medical conditions such as obesity, Crohn\u27s disease, and diabetes. However, connecting microbial community composition to ecosystem function remains a challenge. Here, we introduce MICOM, a customizable metabolic model of the human gut microbiome. By using a heuristic optimization approach based on L2 regularization, we were able to obtain a unique set of realistic growth rates that corresponded well with observed replication rates. We integrated adjustable dietary and taxon abundance constraints to generate personalized metabolic models for individual metagenomic samples. We applied MICOM to a balanced cohort of metagenomes from 186 people, including a metabolically healthy population and individuals with type 1 and type 2 diabetes. Model results showed that individual bacterial genera maintained conserved niche structures across humans, while the community-level production of short-chain fatty acids (SCFAs) was heterogeneous and highly individual specific. Model output revealed complex cross-feeding interactions that would be difficult to measur

Boolean modeling reveals that cyclic attractors in macrophage polarization serve as reservoirs of states to balance external perturbations from the tumor microenvironment

Cyclic attractors generated from Boolean models may explain the adaptability of a cell in response to a dynamical complex tumor microenvironment. In contrast to this idea, we postulate that cyclic attractors in certain cases could be a systemic mechanism to face the perturbations coming from the environment. To justify our conjecture, we present a dynamic analysis of a highly curated transcriptional regulatory network of macrophages constrained into a cancer microenvironment. We observed that when M1-associated transcription factors (STAT1 or NF-κB) are perturbed and the microenvironment balances to a hyper-inflammation condition, cycle attractors activate genes whose signals counteract this effect implicated in tissue damage. The same behavior happens when the M2-associated transcription factors are disturbed (STAT3 or STAT6); cycle attractors will prevent a hyper-regulation scenario implicated in providing a suitable environment for tumor growth. Therefore, here we propose that cyclic macrophage phenotypes can serve as a reservoir for balancing the phenotypes when a specific phenotype-based transcription factor is perturbed in the regulatory network of macrophages. We consider that cyclic attractors should not be simply ignored, but it is necessary to carefully evaluate their biological importance. In this work, we suggest one conjecture: the cyclic attractors can serve as a reservoir to balance the inflammatory/regulatory response of the network under external perturbations

Systems biology of bacterial nitrogen fixation: High-throughput technology and its integrative description with constraint-based modeling

<p>Abstract</p> <p>Background</p> <p>Bacterial nitrogen fixation is the biological process by which atmospheric nitrogen is uptaken by bacteroids located in plant root nodules and converted into ammonium through the enzymatic activity of nitrogenase. In practice, this biological process serves as a natural form of fertilization and its optimization has significant implications in sustainable agricultural programs. Currently, the advent of high-throughput technology supplies with valuable data that contribute to understanding the metabolic activity during bacterial nitrogen fixation. This undertaking is not trivial, and the development of computational methods useful in accomplishing an integrative, descriptive and predictive framework is a crucial issue to decoding the principles that regulated the metabolic activity of this biological process.</p> <p>Results</p> <p>In this work we present a systems biology description of the metabolic activity in bacterial nitrogen fixation. This was accomplished by an integrative analysis involving high-throughput data and constraint-based modeling to characterize the metabolic activity in <it>Rhizobium etli </it>bacteroids located at the root nodules of <it>Phaseolus vulgaris (</it>bean plant). Proteome and transcriptome technologies led us to identify 415 proteins and 689 up-regulated genes that orchestrate this biological process. Taking into account these data, we: 1) extended the metabolic reconstruction reported for <it>R. etli</it>; 2) simulated the metabolic activity during symbiotic nitrogen fixation; and 3) evaluated the <it>in silico </it>results in terms of bacteria phenotype. Notably, constraint-based modeling simulated nitrogen fixation activity in such a way that 76.83% of the enzymes and 69.48% of the genes were experimentally justified. Finally, to further assess the predictive scope of the computational model, gene deletion analysis was carried out on nine metabolic enzymes. Our model concluded that an altered metabolic activity on these enzymes induced different effects in nitrogen fixation, all of these in qualitative agreement with observations made in <it>R. etli </it>and other <it>Rhizobiaceas</it>.</p> <p>Conclusions</p> <p>In this work we present a genome scale study of the metabolic activity in bacterial nitrogen fixation. This approach leads us to construct a computational model that serves as a guide for 1) integrating high-throughput data, 2) describing and predicting metabolic activity, and 3) designing experiments to explore the genotype-phenotype relationship in bacterial nitrogen fixation.</p

Metabolic Reconstruction and Modeling of Nitrogen Fixation in Rhizobium etli

Rhizobiaceas are bacteria that fix nitrogen during symbiosis with plants. This symbiotic relationship is crucial for the nitrogen cycle, and understanding symbiotic mechanisms is a scientific challenge with direct applications in agronomy and plant development. Rhizobium etli is a bacteria which provides legumes with ammonia (among other chemical compounds), thereby stimulating plant growth. A genome-scale approach, integrating the biochemical information available for R. etli, constitutes an important step toward understanding the symbiotic relationship and its possible improvement. In this work we present a genome-scale metabolic reconstruction (iOR363) for R. etli CFN42, which includes 387 metabolic and transport reactions across 26 metabolic pathways. This model was used to analyze the physiological capabilities of R. etli during stages of nitrogen fixation. To study the physiological capacities in silico, an objective function was formulated to simulate symbiotic nitrogen fixation. Flux balance analysis (FBA) was performed, and the predicted active metabolic pathways agreed qualitatively with experimental observations. In addition, predictions for the effects of gene deletions during nitrogen fixation in Rhizobia in silico also agreed with reported experimental data. Overall, we present some evidence supporting that FBA of the reconstructed metabolic network for R. etli provides results that are in agreement with physiological observations. Thus, as for other organisms, the reconstructed genome-scale metabolic network provides an important framework which allows us to compare model predictions with experimental measurements and eventually generate hypotheses on ways to improve nitrogen fixation

Proteomic patterns of cervical cancer cell lines, a network perspective

<p>Abstract</p> <p>Background</p> <p>Cervical cancer is a major mortality factor in the female population. This neoplastic is an excellent model for studying the mechanisms involved in cancer maintenance, because the Human Papilloma Virus (HPV) is the etiology factor in most cases. With the purpose of characterizing the effects of malignant transformation in cellular activity, proteomic studies constitute a reliable way to monitor the biological alterations induced by this disease. In this contextual scheme, a systemic description that enables the identification of the common events between cell lines of different origins, is required to distinguish the essence of carcinogenesis.</p> <p>Results</p> <p>With this study, we sought to achieve a systemic perspective of the common proteomic profile of six cervical cancer cell lines, both positive and negative for HPV, and which differ from the profile corresponding to the non-tumourgenic cell line, HaCaT. Our objectives were to identify common cellular events participating in cancer maintenance, as well as the establishment of a pipeline to work with proteomic-derived results. We analyzed by means of 2D SDS-PAGE and MALDI-TOF mass spectrometry the protein extracts of six cervical cancer cell lines, from which we identified a consensus of 66 proteins. We call this group of proteins, the "central core of cervical cancer". Starting from this core set of proteins, we acquired a PPI network that pointed, through topological analysis, to some proteins that may well be playing a central role in the neoplastic process, such as 14-3-3ζ. <it>In silico </it>overrepresentation analysis of transcription factors pointed to the overexpression of c-Myc, Max and E2F1 as key transcription factors involved in orchestrating the neoplastic phenotype.</p> <p>Conclusions</p> <p>Our findings show that there is a "central core of cervical cancer" protein expression pattern, and suggest that 14-3-3ζ is key to determine if the cell proliferates or dies. In addition, our bioinformatics analysis suggests that the neoplastic phenotype is governed by a non-canonical regulatory pathway.</p

Identifying dynamical modules from genetic regulatory systems: applications to the segment polarity network

BACKGROUND It is widely accepted that genetic regulatory systems are 'modular', in that the whole system is made up of smaller 'subsystems' corresponding to specific biological functions. Most attempts to identify modules in genetic regulatory systems have relied on the topology of the underlying network. However, it is the temporal activity (dynamics) of genes and proteins that corresponds to biological functions, and hence it is dynamics that we focus on here for identifying subsystems. RESULTS Using Boolean network models as an exemplar, we present a new technique to identify subsystems, based on their dynamical properties. The main part of the method depends only on the stable dynamics (attractors) of the system, thus requiring no prior knowledge of the underlying network. However, knowledge of the logical relationships between the network components can be used to describe how each subsystem is regulated. To demonstrate its applicability to genetic regulatory systems, we apply the method to a model of the Drosophila segment polarity network, providing a detailed breakdown of the system. CONCLUSION We have designed a technique for decomposing any set of discrete-state, discrete-time attractors into subsystems. Having a suitable mathematical model also allows us to describe how each subsystem is regulated and how robust each subsystem is against perturbations. However, since the subsystems are found directly from the attractors, a mathematical model or underlying network topology is not necessarily required to identify them, potentially allowing the method to be applied directly to experimental expression data

Transcriptional and Microenvironmental Landscape of Macrophage Transition in Cancer: A Boolean Analysis

The balance between pro- and anti-inflammatory immune system responses is crucial to face and counteract complex diseases such as cancer. Macrophages are an essential population that contributes to this balance in collusion with the local tumor microenvironment. Cancer cells evade the attack of macrophages by liberating cytokines and enhancing the transition to the M2 phenotype with pro-tumoral functions. Despite this pernicious effect on immune systems, the M1 phenotype still exists in the environment and can eliminate tumor cells by liberating cytokines that recruit and activate the cytotoxic actions of TH1 effector cells. Here, we used a Boolean modeling approach to understand how the tumor microenvironment shapes macrophage behavior to enhance pro-tumoral functions. Our network reconstruction integrates experimental data and public information that let us study the polarization from monocytes to M1, M2a, M2b, M2c, and M2d subphenotypes. To analyze the dynamics of our model, we modeled macrophage polarization in different conditions and perturbations. Notably, our study identified new hybrid cell populations, undescribed before. Based on the in vivo macrophage behavior, we explained the hybrid macrophages’ role in the tumor microenvironment. The in silico model allowed us to postulate transcriptional factors that maintain the balance between macrophages with anti- and pro-tumoral functions. In our pursuit to maintain the balance of macrophage phenotypes to eliminate malignant tumor cells, we emulated a theoretical genetically modified macrophage by modifying the activation of NFκB and a loss of function in HIF1-α and discussed their phenotype implications. Overall, our theoretical approach is as a guide to design new experiments for unraveling the principles of the dual host-protective or -harmful antagonistic roles of transitional macrophages in tumor immunoediting and cancer cell fate decisions

A network perspective on the ecology of gut microbiota and progression of type 2 diabetes: Linkages to keystone taxa in a Mexican cohort

IntroductionThe human gut microbiota (GM) is a dynamic system which ecological interactions among the community members affect the host metabolism. Understanding the principles that rule the bidirectional communication between GM and its host, is one of the most valuable enterprise for uncovering how bacterial ecology influences the clinical variables in the host.MethodsHere, we used SparCC to infer association networks in 16S rRNA gene amplicon data from the GM of a cohort of Mexican patients with type 2 diabetes (T2D) in different stages: NG (normoglycemic), IFG (impaired fasting glucose), IGT (impaired glucose tolerance), IFG + IGT (impaired fasting glucose plus impaired glucose tolerance), T2D and T2D treated (T2D with a 5-year ongoing treatment).ResultsBy exploring the network topology from the different stages of T2D, we observed that, as the disease progress, the networks lose the association between bacteria. It suggests that the microbial community becomes highly sensitive to perturbations in individuals with T2D. With the purpose to identify those genera that guide this transition, we computationally found keystone taxa (driver nodes) and core genera for a Mexican T2D cohort. Altogether, we suggest a set of genera driving the progress of the T2D in a Mexican cohort, among them Ruminococcaceae NK4A214 group, Ruminococcaceae UCG-010, Ruminococcaceae UCG-002, Ruminococcaceae UCG-005, Alistipes, Anaerostipes, and Terrisporobacter.DiscussionBased on a network approach, this study suggests a set of genera that can serve as a potential biomarker to distinguish the distinct degree of advances in T2D for a Mexican cohort of patients. Beyond limiting our conclusion to one population, we present a computational pipeline to link ecological networks and clinical stages in T2D, and desirable aim to advance in the field of precision medicine

Modeling Core Metabolism in Cancer Cells: Surveying the Topology Underlying the Warburg Effect

BACKGROUND: Alterations on glucose consumption and biosynthetic activity of amino acids, lipids and nucleotides are metabolic changes for sustaining cell proliferation in cancer cells. Irrevocable evidence of this fact is the Warburg effect which establishes that cancer cells prefers glycolysis over oxidative phosphorylation to generate ATP. Regulatory action over metabolic enzymes has opened a new window for designing more effective anti-cancer treatments. This enterprise is not trivial and the development of computational models that contribute to identifying potential enzymes for breaking the robustness of cancer cells is a priority. METHODOLOGY/PRINCIPAL FINDINGS: This work presents a constraint-base modeling of the most experimentally studied metabolic pathways supporting cancer cells: glycolysis, TCA cycle, pentose phosphate, glutaminolysis and oxidative phosphorylation. To evaluate its predictive capacities, a growth kinetics study for Hela cell lines was accomplished and qualitatively compared with in silico predictions. Furthermore, based on pure computational criteria, we concluded that a set of enzymes (such as lactate dehydrogenase and pyruvate dehydrogenase) perform a pivotal role in cancer cell growth, findings supported by an experimental counterpart. CONCLUSIONS/SIGNIFICANCE: Alterations on metabolic activity are crucial to initiate and sustain cancer phenotype. In this work, we analyzed the phenotype capacities emerged from a constructed metabolic network conformed by the most experimentally studied pathways sustaining cancer cell growth. Remarkably, in silico model was able to resemble the physiological conditions in cancer cells and successfully identified some enzymes currently studied by its therapeutic effect. Overall, we supplied evidence that constraint-based modeling constitutes a promising computational platform to: 1) integrate high throughput technology and establish a crosstalk between experimental validation and in silico prediction in cancer cell phenotype; 2) explore the fundamental metabolic mechanism that confers robustness in cancer; and 3) suggest new metabolic targets for anticancer treatments. All these issues being central to explore cancer cell metabolism from a systems biology perspective