The three-sided right-handed ß-helix is a versatile fold for glycan interactions.

Interactions between proteins and glycans are critical to various biological processes. With databases of carbohydrate-interacting proteins and increasing amounts of structural data, the three-sided right-handed ß-helix (RHBH) has emerged as a significant structural fold for glycan interactions. In this review, we provide an overview of the sequence, mechanistic, and structural features that enable the RHBH to interact with glycans. The RHBH is a prevalent fold that exists in eukaryotes, prokaryotes, and viruses associated with adhesin and carbohydrate-active enzyme (CAZyme) functions. An evolutionary trajectory analysis on structurally characterized RHBH-containing proteins shows that they likely evolved from carbohydrate-binding proteins with their carbohydrate-degrading activities evolving later. By examining three polysaccharide lyase and three glycoside hydrolase structures, we provide a detailed view of the modes of glycan binding in RHBH proteins. The 3-dimensional shape of the RHBH creates an electrostatically and spatially favorable glycan binding surface that allows for extensive hydrogen bonding interactions, leading to favorable and stable glycan binding. The RHBH is observed to be an adaptable domain capable of being modified with loop insertions and charge inversions to accommodate heterogeneous and flexible glycans and diverse reaction mechanisms. Understanding this prevalent protein fold can advance our knowledge of glycan binding in biological systems and help guide the efficient design and utilization of RHBH-containing proteins in glycobiology research.

Large-scale genomic survey with deep learning-based method reveals strain-level phage specificity determinants.

BACKGROUND: Phage therapy, reemerging as a promising approach to counter antimicrobial-resistant infections, relies on a comprehensive understanding of the specificity of individual phages. Yet the significant diversity within phage populations presents a considerable challenge. Currently, there is a notable lack of tools designed for large-scale characterization of phage receptor-binding proteins, which are crucial in determining the phage host range. RESULTS: In this study, we present SpikeHunter, a deep learning method based on the ESM-2 protein language model. With SpikeHunter, we identified 231,965 diverse phage-encoded tailspike proteins, a crucial determinant of phage specificity that targets bacterial polysaccharide receptors, across 787,566 bacterial genomes from 5 virulent, antibiotic-resistant pathogens. Notably, 86.60% (143,200) of these proteins exhibited strong associations with specific bacterial polysaccharides. We discovered that phages with identical tailspike proteins can infect different bacterial species with similar polysaccharide receptors, underscoring the pivotal role of tailspike proteins in determining host range. The specificity is mainly attributed to the protein's C-terminal domain, which strictly correlates with host specificity during domain swapping in tailspike proteins. Importantly, our dataset-driven predictions of phage-host specificity closely match the phage-host pairs observed in real-world phage therapy cases we studied. CONCLUSIONS: Our research provides a rich resource, including both the method and a database derived from a large-scale genomics survey. This substantially enhances understanding of phage specificity determinants at the strain level and offers a valuable framework for guiding phage selection in therapeutic applications.

Annotating microbial functions with ProkFunFind.

Analyzing microbial genomes has become an essential part of microbiology research, giving valuable insights into the functions and evolution of microbial species. Identifying genes of interest and assigning putative annotations to those genes is a central task in genome analysis, and a plethora of tools and approaches have been developed for this task. The ProkFunFind tool was developed to bridge the gap between these various annotation approaches, providing a flexible and customizable search approach to annotate microbial functions. ProkFunFind is designed around hierarchical definitions of biological functions, where individual genes can be identified using heterogeneous search terms consisting of sequences, profile hidden Markov models, protein domains, and orthology groups. This flexible and customizable search approach allows for searches to be tailored to specific biological functions, and the search results are output in multiple formats to facilitate downstream analyses. The utility of the ProkFunFind search tool was demonstrated through its application in searching for bacterial flagella, which are complex organelles composed of multiple genes. Overall, ProkFunFind provides an accessible and flexible way to integrate multiple types of annotation and sequence data while annotating biological functions in microbial genomes.IMPORTANCEGenome sequencing and analysis are increasingly important parts of microbiology, providing a way to predict metabolic functions, identify virulence factors, and understand the evolution of microbes. The expanded use of genome sequencing has also brought an abundance of search and annotation methods, but integrating the information from these different methods can be challenging and is often done through ad hoc approaches. To bridge the gap between different types of annotations, we developed ProkFunFind, a flexible and customizable search tool incorporating multiple search approaches and annotation types to annotate microbial functions. We demonstrated the utility of ProkFunFind by searching for gene clusters encoding flagellar genes using a combination of different annotation types and searches. Overall, ProkFunFind provides a reproducible and flexible way to identify gene clusters of interest, facilitating the meaningful analysis of new and existing microbial genomes.

BilR is a gut microbial enzyme that reduces bilirubin to urobilinogen.

Metabolism of haem by-products such as bilirubin by humans and their gut microbiota is essential to human health, as excess serum bilirubin can cause jaundice and even neurological damage. The bacterial enzymes that reduce bilirubin to urobilinogen, a key step in this pathway, have remained unidentified. Here we used biochemical analyses and comparative genomics to identify BilR as a gut-microbiota-derived bilirubin reductase that reduces bilirubin to urobilinogen. We delineated the BilR sequences from similar reductases through the identification of key residues critical for bilirubin reduction and found that BilR is predominantly encoded by Firmicutes species. Analysis of human gut metagenomes revealed that BilR is nearly ubiquitous in healthy adults, but prevalence is decreased in neonates and individuals with inflammatory bowel disease. This discovery sheds light on the role of the gut microbiome in bilirubin metabolism and highlights the significance of the gut-liver axis in maintaining bilirubin homeostasis.

Deciphering Phage-Host Specificity Based on the Association of Phage Depolymerases and Bacterial Surface Glycan with Deep Learning.

Phage tailspike proteins are depolymerases that target diverse bacterial surface glycans with high specificity, determining the host-specificity of numerous phages. To address the challenge of identifying tailspike proteins due to their sequence diversity, we developed SpikeHunter, an approach based on the ESM-2 protein language model. Using SpikeHunter, we successfully identified 231,965 tailspike proteins from a dataset comprising 8,434,494 prophages found within 165,365 genomes of five common pathogens. Among these proteins, 143,035 tailspike proteins displayed strong associations with serotypes. Moreover, we observed highly similar tailspike proteins in species that share closely related serotypes. We found extensive domain swapping in all five species, with the C-terminal domain being significantly associated with host serotype highlighting its role in host range determination. Our study presents a comprehensive cross-species analysis of tailspike protein to serotype associations, providing insights applicable to phage therapy and biotechnology.

Wastewater surveillance uncovers regional diversity and dynamics of SARS-CoV-2 variants across nine states in the USA.

Wastewater-based epidemiology (WBE) is a non-invasive and cost-effective approach for monitoring the spread of a pathogen within a community. WBE has been adopted as one of the methods to monitor the spread and population dynamics of the SARS-CoV-2 virus, but significant challenges remain in the bioinformatic analysis of WBE-derived data. Here, we have developed a new distance metric, CoVdist, and an associated analysis tool that facilitates the application of ordination analysis to WBE data and the identification of viral population changes based on nucleotide variants. We applied these new approaches to a large-scale dataset from 18 cities in nine states of the USA using wastewater collected from July 2021 to June 2022. We found that the trends in the shift between the Delta and Omicron SARS-CoV-2 lineages were largely consistent with what was seen in clinical data, but that wastewater analysis offered the added benefit of revealing significant differences in viral population dynamics at the state, city, and even neighborhood scales. We also were able to observe the early spread of variants of concern and the presence of recombinant lineages during the transitions between variants, both of which are challenging to analyze based on clinically-derived viral genomes. The methods outlined here will be beneficial for future applications of WBE to monitor SARS-CoV-2, particularly as clinical monitoring becomes less prevalent. Additionally, these approaches are generalizable, allowing them to be applied for the monitoring and analysis of future viral outbreaks.

Discovery of the gut microbial enzyme responsible for bilirubin reduction to urobilinogen.

The degradation of heme and the interplay of its catabolic derivative, bilirubin, between humans and their gut microbiota is an essential facet of human health. However, the hypothesized bacterial enzyme that reduces bilirubin to urobilinogen, a key step that produces the excretable waste products of this pathway, has remained unidentified. In this study, we used a combination of biochemical analyses and comparative genomics to identify a novel enzyme, BilR, that can reduce bilirubin to urobilinogen. We delineated the BilR sequences from other members of the Old Yellow Enzyme family through the identification of key residues in the active site that are critical for bilirubin reduction and found that BilR is predominantly encoded by Firmicutes in the gut microbiome. Our analysis of human gut metagenomes showed that BilR is a common feature of a healthy adult human microbiome but has a decreased prevalence in neonates and IBD patients. This discovery sheds new light on the role of the gut microbiome in bilirubin metabolism and highlights the significance of the gut-liver axis in maintaining bilirubin homeostasis.

Reconstruction and Analysis of Thermodynamically Constrained Models Reveal Metabolic Responses of a Deep-Sea Bacterium to Temperature Perturbations.

Microbial acclimation to different temperature conditions can involve broad changes in cell composition and metabolic efficiency. A systems-level view of these metabolic responses in nonmesophilic organisms, however, is currently missing. In this study, thermodynamically constrained genome-scale models were applied to simulate the metabolic responses of a deep-sea psychrophilic bacterium, Shewanella psychrophila WP2, under suboptimal (4°C), optimal (15°C), and supraoptimal (20°C) growth temperatures. The models were calibrated with experimentally determined growth rates of WP2. Gibbs free energy change of reactions (ΔrG'), metabolic fluxes, and metabolite concentrations were predicted using random simulations to characterize temperature-dependent changes in the metabolism. The modeling revealed the highest metabolic efficiency at the optimal temperature, and it suggested distinct patterns of ATP production and consumption that could lead to lower metabolic efficiency under suboptimal or supraoptimal temperatures. The modeling also predicted rearrangement of fluxes through multiple metabolic pathways, including the glycolysis pathway, Entner-Doudoroff pathway, tricarboxylic acid (TCA) cycle, and electron transport system, and these predictions were corroborated through comparisons to WP2 transcriptomes. Furthermore, predictions of metabolite concentrations revealed the potential conservation of reducing equivalents and ATP in the suboptimal temperature, consistent with experimental observations from other psychrophiles. Taken together, the WP2 models provided mechanistic insights into the metabolism of a psychrophile in response to different temperatures. IMPORTANCE Metabolic flexibility is a central component of any organism's ability to survive and adapt to changes in environmental conditions. This study represents the first application of thermodynamically constrained genome-scale models in simulating the metabolic responses of a deep-sea psychrophilic bacterium to various temperatures. The models predicted differences in metabolic efficiency that were attributed to changes in metabolic pathway utilization and metabolite concentration during growth under optimal and nonoptimal temperatures. Experimental growth measurements were used for model calibration, and temperature-dependent transcriptomic changes corroborated the model-predicted rearrangement of metabolic fluxes. Overall, this study highlights the utility of modeling approaches in studying the temperature-driven metabolic responses of an extremophilic organism.

Putative Host-Derived Insertions in the Genomes of Circulating SARS-CoV-2 Variants.

Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2-infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template-switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution. IMPORTANCE Throughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show that host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic material, potentially enhancing its adaptability, virulence, and spread.

Taxonomic distribution and evolutionary analysis of the equol biosynthesis gene cluster.

BACKGROUND: Equol, an isoflavonoid metabolite with possible health benefits in humans, is known to be produced by some human gut bacteria. While the genes encoding the equol production pathway have been characterized in a few bacterial strains, a systematic analysis of the equol production pathway is currently lacking. RESULTS: This study presents an analysis of the taxonomic distribution and evolutionary history of the gene cluster encoding the equol production pathway. A survey for equol gene clusters within the Genome Taxonomy Database bacterial genomes and human gut metagenomes resulted in the identification of a highly conserved gene cluster found in nine bacterial species from the Eggerthellaceae family. The identified gene clusters from human gut metagenomes revealed potential variations in the equol gene cluster organization and gene content within the equol-producing Eggerthellaceae clades. Subsequent analysis showed that in addition to the four genes directly involved in equol production, multiple other genes were consistently found in the equol gene clusters. These genes were predicted to encode a putative electron transport complex and hydrogenase maturase system, suggesting potential roles for them in the equol production pathway. Analysis of the gene clusters and a phylogenetic reconstruction of a putative NAD kinase gene provided evidence of the recent transfer of the equol gene cluster from a basal Eggerthellaceae species to Slackia_A equolifaciens, Enteroscipio sp000270285, and Lactococcus garvieae 20-92. CONCLUSIONS: This analysis demonstrates that the highly conserved equol gene cluster is taxonomically restricted to the Eggerthellaceae family of bacteria and provides evidence of the role of horizontal gene transfer in the evolutionary history of these genes. These results provide a foundation for future studies of equol production in the human gut and future efforts related to bioengineering and the use of equol-producing bacteria as probiotics.

Putative host-derived insertions in the genomes of circulating SARS-CoV-2 variants.

Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2 infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution. IMPORTANCE: Throughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic materials, potentially enhancing its adaptability, virulence, and spread.

Microbiome Analysis Reveals Diversity and Function of Mollicutes Associated with the Eastern Oyster, Crassostrea virginica.

Marine invertebrate microbiomes play important roles in diverse host and ecological processes. However, a mechanistic understanding of host-microbe interactions is currently available for a small number of model organisms. Here, an integrated taxonomic and functional analysis of the microbiome of the eastern oyster, Crassostrea virginica, was performed using 16S rRNA gene-based amplicon profiling, shotgun metagenomics, and genome-scale metabolic reconstruction. Relatively high variability of the microbiome was observed across individual oysters and among different tissue types. Specifically, a significantly higher alpha diversity was observed in the inner shell than in the gut, gill, mantle, and pallial fluid samples, and a distinct microbiome composition was revealed in the gut compared to other tissues examined in this study. Targeted metagenomic sequencing of the gut microbiota led to further characterization of a dominant bacterial taxon, the class Mollicutes, which was captured by the reconstruction of a metagenome-assembled genome (MAG). Genome-scale metabolic reconstruction of the oyster Mollicutes MAG revealed a reduced set of metabolic functions and a high reliance on the uptake of host-derived nutrients. A chitin degradation and an arginine deiminase pathway were unique to the MAG compared to closely related genomes of Mollicutes isolates, indicating distinct mechanisms of carbon and energy acquisition by the oyster-associated Mollicutes A systematic reanalysis of public eastern oyster-derived microbiome data revealed a high prevalence of the Mollicutes among adult oyster guts and a significantly lower relative abundance of the Mollicutes in oyster larvae and adult oyster biodeposits.IMPORTANCE Despite their biological and ecological significance, a mechanistic characterization of microbiome function is frequently missing from many nonmodel marine invertebrates. As an initial step toward filling this gap for the eastern oyster, Crassostrea virginica, this study provides an integrated taxonomic and functional analysis of the oyster microbiome using samples from a coastal salt pond in August 2017. The study identified high variability of the microbiome across tissue types and among individual oysters, with some dominant taxa showing higher relative abundance in specific tissues. A high prevalence of Mollicutes in the adult oyster gut was revealed by comparative analysis of the gut, biodeposit, and larva microbiomes. Phylogenomic analysis and metabolic reconstruction suggested the oyster-associated Mollicutes is closely related but functionally distinct from Mollicutes isolated from other marine invertebrates. To the best of our knowledge, this study represents the first metagenomics-derived functional inference of Mollicutes in the eastern oyster microbiome.

Correction to: The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest.

An amendment to this paper has been published and can be accessed via the original article.

The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest.

BACKGROUND: Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture. RESULTS: Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world. CONCLUSIONS: The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

FindPrimaryPairs: An efficient algorithm for predicting element-transferring reactant/product pairs in metabolic networks.

The metabolism of individual organisms and biological communities can be viewed as a network of metabolites connected to each other through chemical reactions. In metabolic networks, chemical reactions transform reactants into products, thereby transferring elements between these metabolites. Knowledge of how elements are transferred through reactant/product pairs allows for the identification of primary compound connections through a metabolic network. However, such information is not readily available and is often challenging to obtain for large reaction databases or genome-scale metabolic models. In this study, a new algorithm was developed for automatically predicting the element-transferring reactant/product pairs using the limited information available in the standard representation of metabolic networks. The algorithm demonstrated high efficiency in analyzing large datasets and provided accurate predictions when benchmarked with manually curated data. Applying the algorithm to the visualization of metabolic networks highlighted pathways of primary reactant/product connections and provided an organized view of element-transferring biochemical transformations. The algorithm was implemented as a new function in the open source software package PSAMM in the release v0.30 (https://zhanglab.github.io/psamm/).

Using PSAMM for the Curation and Analysis of Genome-Scale Metabolic Models.

PSAMM is an open source software package that supports the iterative curation and analysis of genome-scale models (GEMs). It aims to integrate the annotation and consistency checking of metabolic models with the simulation of metabolic fluxes. The model representation in PSAMM is compatible with version tracking systems like Git, which allows for full documentation of model file changes and enables collaborative curations of large, complex models. This chapter provides a protocol for using PSAMM functions and a detailed description of the various aspects in setting up and using PSAMM for the simulation and analysis of metabolic models. The overall PSAMM workflow outlined in this chapter includes the import and export of model files, the documentation of model modifications using the Git version control system, the application of consistency checking functions for model curations, and the numerical simulation of metabolic models.

A Genome-Scale Model of Shewanella piezotolerans Simulates Mechanisms of Metabolic Diversity and Energy Conservation.

Shewanella piezotolerans strain WP3 belongs to the group 1 branch of the Shewanella genus and is a piezotolerant and psychrotolerant species isolated from the deep sea. In this study, a genome-scale model was constructed for WP3 using a combination of genome annotation, ortholog mapping, and physiological verification. The metabolic reconstruction contained 806 genes, 653 metabolites, and 922 reactions, including central metabolic functions that represented nonhomologous replacements between the group 1 and group 2 Shewanella species. Metabolic simulations with the WP3 model demonstrated consistency with existing knowledge about the physiology of the organism. A comparison of model simulations with experimental measurements verified the predicted growth profiles under increasing concentrations of carbon sources. The WP3 model was applied to study mechanisms of anaerobic respiration through investigating energy conservation, redox balancing, and the generation of proton motive force. Despite being an obligate respiratory organism, WP3 was predicted to use substrate-level phosphorylation as the primary source of energy conservation under anaerobic conditions, a trait previously identified in other Shewanella species. Further investigation of the ATP synthase activity revealed a positive correlation between the availability of reducing equivalents in the cell and the directionality of the ATP synthase reaction flux. Comparison of the WP3 model with an existing model of a group 2 species, Shewanella oneidensis MR-1, revealed that the WP3 model demonstrated greater flexibility in ATP production under the anaerobic conditions. Such flexibility could be advantageous to WP3 for its adaptation to fluctuating availability of organic carbon sources in the deep sea. IMPORTANCE The well-studied nature of the metabolic diversity of Shewanella bacteria makes species from this genus a promising platform for investigating the evolution of carbon metabolism and energy conservation. The Shewanella phylogeny is diverged into two major branches, referred to as group 1 and group 2. While the genotype-phenotype connections of group 2 species have been extensively studied with metabolic modeling, a genome-scale model has been missing for the group 1 species. The metabolic reconstruction of Shewanella piezotolerans strain WP3 represented the first model for Shewanella group 1 and the first model among piezotolerant and psychrotolerant deep-sea bacteria. The model brought insights into the mechanisms of energy conservation in WP3 under anaerobic conditions and highlighted its metabolic flexibility in using diverse carbon sources. Overall, the model opens up new opportunities for investigating energy conservation and metabolic adaptation, and it provides a prototype for systems-level modeling of other deep-sea microorganisms.

PSAMM: A Portable System for the Analysis of Metabolic Models.

The genome-scale models of metabolic networks have been broadly applied in phenotype prediction, evolutionary reconstruction, community functional analysis, and metabolic engineering. Despite the development of tools that support individual steps along the modeling procedure, it is still difficult to associate mathematical simulation results with the annotation and biological interpretation of metabolic models. In order to solve this problem, here we developed a Portable System for the Analysis of Metabolic Models (PSAMM), a new open-source software package that supports the integration of heterogeneous metadata in model annotations and provides a user-friendly interface for the analysis of metabolic models. PSAMM is independent of paid software environments like MATLAB, and all its dependencies are freely available for academic users. Compared to existing tools, PSAMM significantly reduced the running time of constraint-based analysis and enabled flexible settings of simulation parameters using simple one-line commands. The integration of heterogeneous, model-specific annotation information in PSAMM is achieved with a novel format of YAML-based model representation, which has several advantages, such as providing a modular organization of model components and simulation settings, enabling model version tracking, and permitting the integration of multiple simulation problems. PSAMM also includes a number of quality checking procedures to examine stoichiometric balance and to identify blocked reactions. Applying PSAMM to 57 models collected from current literature, we demonstrated how the software can be used for managing and simulating metabolic models. We identified a number of common inconsistencies in existing models and constructed an updated model repository to document the resolution of these inconsistencies.