Metabolic network percolation quantifies biosynthetic capabilities across the human oral microbiome

The biosynthetic capabilities of microbes underlie their growth and interactions, playing a prominent role in microbial community structure. For large, diverse microbial communities, prediction of these capabilities is limited by uncertainty about metabolic functions and environmental conditions. To address this challenge, we propose a probabilistic method, inspired by percolation theory, to computationally quantify how robustly a genome-derived metabolic network produces a given set of metabolites under an ensemble of variable environments. We used this method to compile an atlas of predicted biosynthetic capabilities for 97 metabolites across 456 human oral microbes. This atlas captures taxonomically-related trends in biomass composition, and makes it possible to estimate inter-microbial metabolic distances that correlate with microbial co-occurrences. We also found a distinct cluster of fastidious/uncultivated taxa, including several Saccharibacteria (TM7) species, characterized by their abundant metabolic deficiencies. By embracing uncertainty, our approach can be broadly applied to understanding metabolic interactions in complex microbial ecosystems.T32GM008764 - NIGMS NIH HHS; T32 GM008764 - NIGMS NIH HHS; R01 DE024468 - NIDCR NIH HHS; R01 GM121950 - NIGMS NIH HHS; DE-SC0012627 - Biological and Environmental Research; RGP0020/2016 - Human Frontier Science Program; NSFOCE-BSF 1635070 - National Science Foundation; HR0011-15-C-0091 - Defense Advanced Research Projects Agency; R37DE016937 - NIDCR NIH HHS; R37 DE016937 - NIDCR NIH HHS; R01GM121950 - NIGMS NIH HHS; R01DE024468 - NIDCR NIH HHS; 1457695 - National Science FoundationPublished versio

Reconstructing phylogeny from metabolic substrate-product relationships

<p> Abstract</p> <p>Background</p> <p>Many approaches utilize metabolic pathway information to reconstruct the phyletic tree of fully sequenced organisms, but how metabolic networks can add information to original genomic annotations has remained open.</p> <p>Methods</p> <p>We translated enzyme reactions assigned in 1075 organisms into substrate-product relationships to represent the metabolic information at a finer resolution than enzymes and compounds. Each organism was represented as a vector of substrate-product relationships and the phyletic tree was reconstructed by a simple hierarchical method. Obtained results were compared with several other approaches that use genome information and network properties.</p> <p>Results</p> <p>Phyletic trees without consideration of network properties can already extract organisms in anomalous environments. This efficient method can add insights to traditional genome-based phylogenetic reconstruction.</p> <p>Conclusions</p> <p>Structural relationship among metabolites can highlight parasitic or symbiont species such as spirochaete and clamydia. The method assists understanding of species-environment interaction when used in combination with traditional phylogenetic methods.</p

Biochemical Networks Across Planets and Scales

abstract: Biochemical reactions underlie all living processes. Their complex web of interactions is difficult to fully capture and quantify with simple mathematical objects. Applying network science to biology has advanced our understanding of the metabolisms of individual organisms and the organization of ecosystems, but has scarcely been applied to life at a planetary scale. To characterize planetary-scale biochemistry, I constructed biochemical networks using global databases of annotated genomes and metagenomes, and biochemical reactions. I uncover scaling laws governing biochemical diversity and network structure shared across levels of organization from individuals to ecosystems, to the biosphere as a whole. Comparing real biochemical reaction networks to random reaction networks reveals the observed biological scaling is not a product of chemistry alone, but instead emerges due to the particular structure of selected reactions commonly participating in living processes. I perform distinguishability tests across properties of individual and ecosystem-level biochemical networks to determine whether or not they share common structure, indicative of common generative mechanisms across levels. My results indicate there is no sharp transition in the organization of biochemistry across distinct levels of the biological hierarchy—a result that holds across different network projections. Finally, I leverage these large biochemical datasets, in conjunction with planetary observations and computational tools, to provide a methodological foundation for the quantitative assessment of biology’s viability amongst other geospheres. Investigating a case study of alkaliphilic prokaryotes in the context of Enceladus, I find that the chemical compounds observed on Enceladus thus far would be insufficient to allow even these extremophiles to produce the compounds necessary to sustain a viable metabolism. The environmental precursors required by these organisms provides a reference for the compounds which should be prioritized for detection in future planetary exploration missions. The results of this framework have further consequences in the context of planetary protection, and hint that forward contamination may prove infeasible without meticulous intent. Taken together these results point to a deeper level of organization in biochemical networks than what has been understood so far, and suggests the existence of common organizing principles operating across different levels of biology and planetary chemistry.Dissertation/ThesisDoctoral Dissertation Geological Sciences 201

A census analysis of the 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase and EPSP-associated domains

Glyphosate is one of the most used herbicides against weeds that targets the enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). EPSPS is the central enzyme in the shikimate pathway to synthesize 3 essential amino acids in plants, fungi, and prokaryotes. Although this pathway is not found in animals, herbicide may affect the biodiversity of environmental and host-associated microorganisms. In this master thesis I will survey the distribution of the EPSPS enzyme in thousands of microorganisms and I will analyse the evolution of the multi domain structure of the EPSPS enzyme in fungi. Data was gathered from public databases of proteins (e.g., Pfam and COG). The analysis of the distribution of the EPSPS was performed using Excel functions and a bipartite network was analysed with the program Cytoscape. The Count program was used to assess evolutionary scenarios by Dollo’s maximum parsimony, and the phylogenetic trees were visualized with iTOL. The EPSPS enzyme is widely distributed in archaea, bacteria, plants, and fungi. The multi domain structure of the EPSPS in fungi is strongly associated with six other genes of the shikimate pathway. The most common multi domain structure is composed by a group of five enzymes (HQ synthase, EPSPS, SKI, DHquinase I and Shikimate DH), which I call in this thesis the “Major 5”. The EPSPS multi domain structure in fungi ranges between two to eight domains. The evolutionary analysis shows that the ancestral of fungi had a multi domain structure of six domains. Thus, there have been domain gains and losses throughout the evolution of the EPSPS in fungi. Further investigations are needed to determine the effect of the EPSPS-associated domains to glyphosate resistance. A scientific article that includes data from this master thesis is publicly available as a preprint at biorxiv and submitted to a peer-reviewed Nature Methods journal

Quantifying biosynthetic network robustness across the human oral microbiome

Metabolic interactions, such as cross-feeding, play a prominent role in microbial communitystructure. For example, they may underlie the ubiquity of uncultivated microorganisms. We investigated this phenomenon in the human oral microbiome, by analyzing microbial metabolic networks derived from sequenced genomes. Specifically, we devised a probabilistic biosynthetic network robustness metric that describes the chance that an organism could produce a given metabolite, and used it to assemble a comprehensive atlas of biosynthetic capabilities for 88 metabolites across 456 human oral microbiome strains. A cluster of organisms characterized by reduced biosynthetic capabilities stood out within this atlas. This cluster included several uncultivated taxa and three recently co-cultured Saccharibacteria (TM7) phylum species. Comparison across strains also allowed us to systematically identify specific putative metabolic interdependences between organisms. Our method, which provides a new way of converting annotated genomes into metabolic predictions, is easily extendible to other microbial communities and metabolic products.https://www.biorxiv.org/content/10.1101/392621v1First author draf

Adaptive Evolution of Phosphorus Metabolism in Prochlorococcus.

Inorganic phosphorus is scarce in the eastern Mediterranean Sea, where the high-light-adapted ecotype HLI of the marine picocyanobacterium Prochlorococcus marinus thrives. Physiological and regulatory control of phosphorus acquisition and partitioning has been observed in HLI both in culture and in the field; however, the optimization of phosphorus metabolism and associated gains for its phosphorus-limited-growth (PLG) phenotype have not been studied. Here, we reconstructed a genome-scale metabolic network of the HLI axenic strain MED4 (iJC568), consisting of 568 metabolic genes in relation to 794 reactions involving 680 metabolites distributed in 6 subcellular locations. iJC568 was used to quantify metabolic fluxes under PLG conditions, and we observed a close correspondence between experimental and computed fluxes. We found that MED4 has minimized its dependence on intracellular phosphate, not only through drastic depletion of phosphorus-containing biomass components but also through network-wide reductions in phosphate-reaction participation and the loss of a key enzyme, succinate dehydrogenase. These alterations occur despite the stringency of having relatively few pathway redundancies and an extremely high proportion of essential metabolic genes (47%; defined as the percentage of lethal in silico gene knockouts). These strategies are examples of nutrient-controlled adaptive evolution and confer a dramatic growth rate advantage to MED4 in phosphorus-limited regions. IMPORTANCE Microbes are known to employ three basic strategies to compete for limiting elemental resources: (i) cell quotas may be adjusted by alterations to cell physiology or by substitution of a more plentiful resource, (ii) stressed cells may synthesize high-affinity transporters, and (iii) cells may access more costly sources from internal stores, by degradation, or by petitioning other microbes. In the case of phosphorus, a limiting resource in vast oceanic regions, the cosmopolitan cyanobacterium Prochlorococcus marinus thrives by adopting all three strategies and a fourth, previously unknown strategy. By generating a detailed model of its metabolism, we found that strain MED4 has evolved a way to reduce its dependence on phosphate by minimizing the number of enzymes involved in phosphate transformations, despite the stringency of nearly half of its metabolic genes being essential for survival. Relieving phosphorus limitation, both physiologically and throughout intermediate metabolism, substantially improves phosphorus-specific growth rates

Analysis of Generative Chemistries

For the modelling of chemistry we use undirected, labelled graphs as explicit models of molecules and graph transformation rules for modelling generalised chemical reactions. This is used to define artificial chemistries on the level of individual bonds and atoms, where formal graph grammars implicitly represent large spaces of chemical compounds. We use a graph rewriting formalism, rooted in category theory, called the Double Pushout approach, which directly expresses the transition state of chemical reactions. Using concurrency theory for transformation rules, we define algorithms for the composition of rewrite rules in a chemically intuitive manner that enable automatic abstraction of the level of detail in chemical pathways. Based on this rule composition we define an algorithmic framework for generation of vast reaction networks for specific spaces of a given chemistry, while still maintaining the level of detail of the model down to the atomic level. The framework also allows for computation with graphs and graph grammars, which is utilised to model non-trivial chemical systems. The graph generation relies on graph isomorphism testing, and we review the general individualisation-refinement paradigm used in the state-of-the-art algorithms for graph canonicalisation, isomorphism testing, and automorphism discovery. We present a model for chemical pathways based on a generalisation of network flows from ordinary directed graphs to directed hypergraphs. The model allows for reasoning about the flow of individual molecules in general pathways, and the introduction of chemically motivated routing constraints. It further provides the foundation for defining specialised pathway motifs, which is illustrated by defining necessary topological constraints for both catalytic and autocatalytic pathways. We also prove that central types of pathway questions are NP-complete, even for restricted classes of reaction networks. The complete pathway model, including constraints for catalytic and autocatalytic pathways, is implemented using integer linear programming. This implementation is used in a tree search method to enumerate both optimal and near-optimal pathway solutions. The formal methods are applied to multiple chemical systems: the enzyme catalysed beta-lactamase reaction, variations of the glycolysis pathway, and the formose process. In each of these systems we use rule composition to abstract pathways and calculate traces for isotope labelled carbon atoms. The pathway model is used to automatically enumerate alternative non-oxidative glycolysis pathways, and enumerate thousands of candidates for autocatalytic pathways in the formose process

Seeing the forest for the trees : retrieving plant secondary biochemical pathways from metabolome networks

Over the last decade, a giant leap forward has been made in resolving the main bottleneck in metabolomics, i.e., the structural characterization of the many unknowns. This has led to the next challenge in this research field: retrieving biochemical pathway information from the various types of networks that can be constructed from metabolome data. Searching putative biochemical pathways, referred to as biotransformation paths, is complicated because several flaws occur during the construction of metabolome networks. Multiple network analysis tools have been developed to deal with these flaws, while in silico retrosynthesis is appearing as an alternative approach. In this review, the different types of metabolome networks, their flaws, and the various tools to trace these biotransformation paths are discussed

Automation of gene assignments to metabolic pathways using high-throughput expression data

BACKGROUND: Accurate assignment of genes to pathways is essential in order to understand the functional role of genes and to map the existing pathways in a given genome. Existing algorithms predict pathways by extrapolating experimental data in one organism to other organisms for which this data is not available. However, current systems classify all genes that belong to a specific EC family to all the pathways that contain the corresponding enzymatic reaction, and thus introduce ambiguity. RESULTS: Here we describe an algorithm for assignment of genes to cellular pathways that addresses this problem by selectively assigning specific genes to pathways. Our algorithm uses the set of experimentally elucidated metabolic pathways from MetaCyc, together with statistical models of enzyme families and expression data to assign genes to enzyme families and pathways by optimizing correlated co-expression, while minimizing conflicts due to shared assignments among pathways. Our algorithm also identifies alternative ("backup") genes and addresses the multi-domain nature of proteins. We apply our model to assign genes to pathways in the Yeast genome and compare the results for genes that were assigned experimentally. Our assignments are consistent with the experimentally verified assignments and reflect characteristic properties of cellular pathways. CONCLUSION: We present an algorithm for automatic assignment of genes to metabolic pathways. The algorithm utilizes expression data and reduces the ambiguity that characterizes assignments that are based only on EC numbers