Polysaccharide utilization loci in Bacteroides determine population fitness and community-level interactions.

Polysaccharide utilization loci (PULs) are co-regulated bacterial genes that sense nutrients and enable glycan digestion. Human gut microbiome members, notably Bacteroides, contain numerous PULs that enable glycan utilization and shape ecological dynamics. To investigate the role of PULs on fitness and inter-species interactions, we develop a CRISPR-based genome editing tool to study 23 PULs in Bacteroides uniformis (BU). BU PULs show distinct glycan-degrading functions and transcriptional coordination that enables the population to adapt upon loss of other PULs. Exploiting a BU mutant barcoding strategy, we demonstrate that in vitro fitness and BU colonization in the murine gut are enhanced by deletion of specific PULs and modulated by glycan availability. PULs mediate glycan-dependent interactions with butyrate producers that depend on the degradation mechanism and glycan utilization ability of the butyrate producer. Thus, PULs determine community dynamics and butyrate production and provide a selective advantage or disadvantage depending on the nutritional landscape.

Design of synthetic human gut microbiome assembly and butyrate production.

The capability to design microbiomes with predictable functions would enable new technologies for applications in health, agriculture, and bioprocessing. Towards this goal, we develop a model-guided approach to design synthetic human gut microbiomes for production of the health-relevant metabolite butyrate. Our data-driven model quantifies microbial interactions impacting growth and butyrate production separately, providing key insights into ecological mechanisms driving butyrate production. We use our model to explore a vast community design space using a design-test-learn cycle to identify high butyrate-producing communities. Our model can accurately predict community assembly and butyrate production across a wide range of species richness. Guided by the model, we identify constraints on butyrate production by high species richness and key molecular factors driving butyrate production, including hydrogen sulfide, environmental pH, and resource competition. In sum, our model-guided approach provides a flexible and generalizable framework for understanding and accurately predicting community assembly and metabolic functions.

Diagnosing and Predicting Mixed-Culture Fermentations with Unicellular and Guild-Based Metabolic Models.

Multispecies microbial communities determine the fate of materials in the environment and can be harnessed to produce beneficial products from renewable resources. In a recent example, fermentations by microbial communities have produced medium-chain fatty acids (MCFAs). Tools to predict, assess, and improve the performance of these communities, however, are limited. To provide such tools, we constructed two metabolic models of MCFA-producing microbial communities based on available genomic, transcriptomic, and metabolomic data. The first model is a unicellular model (iFermCell215), while the second model (iFermGuilds789) separates fermentation activities into functional guilds. Ethanol and lactate are fermentation products known to serve as substrates for MCFA production, while acetate is another common cometabolite during MCFA production. Simulations with iFermCell215 predict that low molar ratios of acetate to ethanol favor MCFA production, whereas the products of lactate and acetate coutilization are less dependent on the acetate-to-lactate ratio. In simulations of an MCFA-producing community fed a complex organic mixture derived from lignocellulose, iFermGuilds789 predicted that lactate was an extracellular cometabolite that served as a substrate for butyrate (C4) production. Extracellular hexanoic (C6) and octanoic (C8) acids were predicted by iFermGuilds789 to be from community members that directly metabolize sugars. Modeling results provide several hypotheses that can improve our understanding of microbial roles in an MCFA-producing microbiome and inform strategies to increase MCFA production. Further, these models represent novel tools for exploring the role of mixed microbial communities in carbon recycling in the environment, as well as in beneficial reuse of organic residues.IMPORTANCE Microbiomes are vital to human health, agriculture, and protecting the environment. Predicting behavior of self-assembled or synthetic microbiomes, however, remains a challenge. In this work, we used unicellular and guild-based metabolic models to investigate production of medium-chain fatty acids by a mixed microbial community that is fed multiple organic substrates. Modeling results provided insights into metabolic pathways of three medium-chain fatty acid-producing guilds and identified potential strategies to increase production of medium-chain fatty acids. This work demonstrates the role of metabolic models in augmenting multi-omic studies to gain greater insights into microbiome behavior.

Understanding and Engineering Distributed Biochemical Pathways in Microbial Communities.

Microbiomes impact nearly every environment on Earth by modulating the molecular composition of the environment. Temporally changing environmental stimuli and spatial organization are major variables shaping the structure and function of microbiomes. The web of interactions among members of these communities and between the organisms and the environment dictates microbiome functions. Microbial interactions are major drivers of microbiomes and are modulated by spatiotemporal parameters. A mechanistic and quantitative understanding of ecological, molecular, and environmental forces shaping microbiomes could inform strategies to control microbiome dynamics and functions. Major challenges for harnessing the potential of microbiomes for diverse applications include the development of predictive modeling frameworks and tools for precise manipulation of microbiome behaviors.

Metatranscriptomic and Thermodynamic Insights into Medium-Chain Fatty Acid Production Using an Anaerobic Microbiome.

Biomanufacturing from renewable feedstocks can offset fossil fuel-based chemical production. One potential biomanufacturing strategy is production of medium-chain fatty acids (MCFA) from organic feedstocks using either pure cultures or microbiomes. While the set of microbes in a microbiome can often metabolize organic materials of greater diversity than a single species can and while the role of specific species may be known, knowledge of the carbon and energy flow within and between organisms in MCFA-producing microbiomes is only now starting to emerge. Here, we integrated metagenomic, metatranscriptomic, and thermodynamic analyses to predict and characterize the metabolic network of an anaerobic microbiome producing MCFA from organic matter derived from lignocellulosic ethanol fermentation conversion residue. A total of 37 high-quality (>80% complete, <10% contamination) metagenome-assembled genomes (MAGs) were recovered from the microbiome, and metabolic reconstruction of the 10 most abundant MAGs was performed. Metabolic reconstruction combined with metatranscriptomic analysis predicted that organisms affiliated with Lactobacillus and Coriobacteriaceae would degrade carbohydrates and ferment sugars to lactate and acetate. Lachnospiraceae- and Eubacteriaceae-affiliated organisms were predicted to transform these fermentation products to MCFA. Thermodynamic analyses identified conditions under which H2 is expected to be either produced or consumed, suggesting a potential role of H2 partial pressure in MCFA production. From an integrated systems analysis perspective, we propose that MCFA production could be improved if microbiomes were engineered to use homofermentative instead of heterofermentative Lactobacillus and if MCFA-producing organisms were engineered to preferentially use a thioesterase instead of a coenzyme A (CoA) transferase as the terminal enzyme in reverse ß-oxidation. IMPORTANCE Mixed communities of microbes play important roles in health, the environment, agriculture, and biotechnology. While tapping the combined activities of organisms within microbiomes may allow the utilization of a wider range of substrates in preference to the use of pure cultures for biomanufacturing, harnessing the metabolism of these mixed cultures remains a major challenge. Here, we predicted metabolic functions of bacteria in a microbiome that produces medium-chain fatty acids from a renewable feedstock. Our findings lay the foundation for efforts to begin addressing how to engineer and control microbiomes for improved biomanufacturing, how to build synthetic mixtures of microbes that produce valuable chemicals from renewable resources, and how to better understand the microbial communities that contribute to health, agriculture, and the environment.

acI Actinobacteria Assemble a Functional Actinorhodopsin with Natively Synthesized Retinal.

Freshwater lakes harbor complex microbial communities, but these ecosystems are often dominated by acI Actinobacteria Members of this cosmopolitan lineage are proposed to bolster heterotrophic growth using phototrophy because their genomes encode actino-opsins (actR). This model has been difficult to validate experimentally because acI Actinobacteria are not consistently culturable. Based primarily on genomes from single cells and metagenomes, we provide a detailed biosynthetic route for members of acI clades A and B to synthesize retinal and its carotenoid precursors. Consequently, acI cells should be able to natively assemble light-driven actinorhodopsins (holo-ActR) to pump protons, unlike many bacteria that encode opsins but may need to exogenously obtain retinal because they lack retinal machinery. Moreover, we show that all acI clades contain genes for a secondary branch of the carotenoid pathway, implying synthesis of a complex carotenoid. Transcription analysis of acI Actinobacteria in a eutrophic lake shows that all retinal and carotenoid pathway operons are transcribed and that actR is among the most highly transcribed of all acI genes. Furthermore, heterologous expression of acI retinal pathway genes showed that lycopene, retinal, and ActR can be made using the genes encoded in these organisms. Model cells producing ActR and the key acI retinal-producing ß-carotene oxygenase formed holo-ActR and acidified solution during illumination. Taken together, our results prove that acI Actinobacteria containing both ActR and acI retinal production machinery have the capacity to natively synthesize a green light-dependent outward proton-pumping rhodopsin.IMPORTANCE Microbes play critical roles in determining the quality of freshwater ecosystems, which are vital to human civilization. Because acI Actinobacteria are ubiquitous and abundant in freshwater lakes, clarifying their ecophysiology is a major step in determining the contributions that they make to nitrogen and carbon cycling. Without accurate knowledge of these cycles, freshwater systems cannot be incorporated into climate change models, ecosystem imbalances cannot be predicted, and policy for service disruption cannot be planned. Our work fills major gaps in microbial light utilization, secondary metabolite production, and energy cycling in freshwater habitats.

TaxAss: Leveraging a Custom Freshwater Database Achieves Fine-Scale Taxonomic Resolution.

Taxonomy assignment of freshwater microbial communities is limited by the minimally curated phylogenies used for large taxonomy databases. Here we introduce TaxAss, a taxonomy assignment workflow that classifies 16S rRNA gene amplicon data using two taxonomy reference databases: a large comprehensive database and a small ecosystem-specific database rigorously curated by scientists within a field. We applied TaxAss to five different freshwater data sets using the comprehensive SILVA database and the freshwater-specific FreshTrain database. TaxAss increased the percentage of the data set classified compared to using only SILVA, especially at fine-resolution family to species taxon levels, while across the freshwater test data sets classifications increased by as much as 11 to 40% of total reads. A similar increase in classifications was not observed in a control mouse gut data set, which was not expected to contain freshwater bacteria. TaxAss also maintained taxonomic richness compared to using only the FreshTrain across all taxon levels from phylum to species. Without TaxAss, most organisms not represented in the FreshTrain were unclassified, but at fine taxon levels, incorrect classifications became significant. We validated TaxAss using simulated amplicon data derived from full-length clone libraries and found that 96 to 99% of test sequences were correctly classified at fine resolution. TaxAss splits a data set's sequences into two groups based on their percent identity to reference sequences in the ecosystem-specific database. Sequences with high similarity to sequences in the ecosystem-specific database are classified using that database, and the others are classified using the comprehensive database. TaxAss is free and open source and is available at https://www.github.com/McMahonLab/TaxAssIMPORTANCE Microbial communities drive ecosystem processes, but microbial community composition analyses using 16S rRNA gene amplicon data sets are limited by the lack of fine-resolution taxonomy classifications. Coarse taxonomic groupings at the phylum, class, and order levels lump ecologically distinct organisms together. To avoid this, many researchers define operational taxonomic units (OTUs) based on clustered sequences, sequence variants, or unique sequences. These fine-resolution groupings are more ecologically relevant, but OTU definitions are data set dependent and cannot be compared between data sets. Microbial ecologists studying freshwater have curated a small, ecosystem-specific taxonomy database to provide consistent and up-to-date terminology. We created TaxAss, a workflow that leverages this database to assign taxonomy. We found that TaxAss improves fine-resolution taxonomic classifications (family, genus, and species). Fine taxonomic groupings are more ecologically relevant, so they provide an alternative to OTU-based analyses that is consistent and comparable between data sets.

Model Communities Hint at Promiscuous Metabolic Linkages between Ubiquitous Free-Living Freshwater Bacteria.

Genome streamlining is frequently observed in free-living aquatic microorganisms and results in physiological dependencies between microorganisms. However, we know little about the specificity of these microbial associations. In order to examine the specificity and extent of these associations, we established mixed cultures from three different freshwater environments and analyzed the cooccurrence of organisms using a metagenomic time series. Free-living microorganisms with streamlined genomes lacking multiple biosynthetic pathways showed no clear recurring pattern in their interaction partners. Free-living freshwater bacteria form promiscuous cooperative associations. This notion contrasts with the well-documented high specificities of interaction partners in host-associated bacteria. Considering all data together, we suggest that highly abundant free-living bacterial lineages are functionally versatile in their interactions despite their distinct streamlining tendencies at the single-cell level. This metabolic versatility facilitates interactions with a variable set of community members.

Metabolic Network Analysis and Metatranscriptomics Reveal Auxotrophies and Nutrient Sources of the Cosmopolitan Freshwater Microbial Lineage acI.

An explosion in the number of available genome sequences obtained through metagenomics and single-cell genomics has enabled a new view of the diversity of microbial life, yet we know surprisingly little about how microbes interact with each other or their environment. In fact, the majority of microbial species remain uncultivated, while our perception of their ecological niches is based on reconstruction of their metabolic potential. In this work, we demonstrate how the "seed set framework," which computes the set of compounds that an organism must acquire from its environment (E. Borenstein, M. Kupiec, M. W. Feldman, and E. Ruppin, Proc Natl Acad Sci U S A 105:14482-14487, 2008, https://doi.org/10.1073/pnas.0806162105), enables computational analysis of metabolic reconstructions while providing new insights into a microbe's metabolic capabilities, such as nutrient use and auxotrophies. We apply this framework to members of the ubiquitous freshwater actinobacterial lineage acI, confirming and extending previous experimental and genomic observations implying that acI bacteria are heterotrophs reliant on peptides and saccharides. We also present the first metatranscriptomic study of the acI lineage, revealing high expression of transport proteins and the light-harvesting protein actinorhodopsin. Putative transport proteins complement predictions of nutrients and essential metabolites while providing additional support of the hypothesis that members of the acI are photoheterotrophs. IMPORTANCE The metabolic activity of uncultivated microorganisms contributes to numerous ecosystem processes, ranging from nutrient cycling in the environment to influencing human health and disease. Advances in sequencing technology have enabled the assembly of genomes for these microorganisms, but our ability to generate reference genomes far outstrips our ability to analyze them. Common approaches to analyzing microbial metabolism require reconstructing the entirety of an organism's metabolic pathways or performing targeted searches for genes involved in a specific process. This paper presents a third approach, in which draft metabolic reconstructions are used to identify compounds through which an organism may interact with its environment. These compounds can then guide more-intensive metabolic reconstruction efforts and can also provide new hypotheses about the specific contributions that microbes make to ecosystem-scale metabolic processes.

Metabolic network analysis reveals microbial community interactions in anammox granules.

Microbial communities mediating anaerobic ammonium oxidation (anammox) represent one of the most energy-efficient environmental biotechnologies for nitrogen removal from wastewater. However, little is known about the functional role heterotrophic bacteria play in anammox granules. Here, we use genome-centric metagenomics to recover 17 draft genomes of anammox and heterotrophic bacteria from a laboratory-scale anammox bioreactor. We combine metabolic network reconstruction with metatranscriptomics to examine the gene expression of anammox and heterotrophic bacteria and to identify their potential interactions. We find that Chlorobi-affiliated bacteria may be highly active protein degraders, catabolizing extracellular peptides while recycling nitrate to nitrite. Other heterotrophs may also contribute to scavenging of detritus and peptides produced by anammox bacteria, and potentially use alternative electron donors, such as H2, acetate and formate. Our findings improve the understanding of metabolic activities and interactions between anammox and heterotrophic bacteria and offer the first transcriptional insights on ecosystem function in anammox granules.

Thermodynamics and H2 Transfer in a Methanogenic, Syntrophic Community.

Microorganisms in nature do not exist in isolation but rather interact with other species in their environment. Some microbes interact via syntrophic associations, in which the metabolic by-products of one species serve as nutrients for another. These associations sustain a variety of natural communities, including those involved in methanogenesis. In anaerobic syntrophic communities, energy is transferred from one species to another, either through direct contact and exchange of electrons, or through small molecule diffusion. Thermodynamics plays an important role in governing these interactions, as the oxidation reactions carried out by the first community member are only possible because degradation products are consumed by the second community member. This work presents the development and analysis of genome-scale network reconstructions of the bacterium Syntrophobacter fumaroxidans and the methanogenic archaeon Methanospirillum hungatei. The models were used to verify proposed mechanisms of ATP production within each species. We then identified additional constraints and the cellular objective function required to match experimental observations. The thermodynamic S. fumaroxidans model could not explain why S. fumaroxidans does not produce H2 in monoculture, indicating that current methods might not adequately estimate the thermodynamics, or that other cellular processes (e.g., regulation) play a role. We also developed a thermodynamic coculture model of the association between the organisms. The coculture model correctly predicted the exchange of both H2 and formate between the two species and suggested conditions under which H2 and formate produced by S. fumaroxidans would be fully consumed by M. hungatei.

Genome-scale reconstruction of metabolic networks of Lactobacillus casei ATCC 334 and 12A.

Lactobacillus casei strains are widely used in industry and the utility of this organism in these industrial applications is strain dependent. Hence, tools capable of predicting strain specific phenotypes would have utility in the selection of strains for specific industrial processes. Genome-scale metabolic models can be utilized to better understand genotype-phenotype relationships and to compare different organisms. To assist in the selection and development of strains with enhanced industrial utility, genome-scale models for L. casei ATCC 334, a well characterized strain, and strain 12A, a corn silage isolate, were constructed. Draft models were generated from RAST genome annotations using the Model SEED database and refined by evaluating ATP generating cycles, mass-and-charge-balances of reactions, and growth phenotypes. After the validation process was finished, we compared the metabolic networks of these two strains to identify metabolic, genetic and ortholog differences that may lead to different phenotypic behaviors. We conclude that the metabolic capabilities of the two networks are highly similar. The L. casei ATCC 334 model accounts for 1,040 reactions, 959 metabolites and 548 genes, while the L. casei 12A model accounts for 1,076 reactions, 979 metabolites and 640 genes. The developed L. casei ATCC 334 and 12A metabolic models will enable better understanding of the physiology of these organisms and be valuable tools in the development and selection of strains with enhanced utility in a variety of industrial applications.

Software platforms to facilitate reconstructing genome-scale metabolic networks.

System-level analyses of microbial metabolism are facilitated by genome-scale reconstructions of microbial biochemical networks. A reconstruction provides a structured representation of the biochemical transformations occurring within an organism, as well as the genes necessary to carry out these transformations, as determined by the annotated genome sequence and experimental data. Network reconstructions also serve as platforms for constraint-based computational techniques, which facilitate biological studies in a variety of applications, including evaluation of network properties, metabolic engineering and drug discovery. Bottom-up metabolic network reconstructions have been developed for dozens of organisms, but until recently, the pace of reconstruction has failed to keep up with advances in genome sequencing. To address this problem, a number of software platforms have been developed to automate parts of the reconstruction process, thereby alleviating much of the manual effort previously required. Here, we review four such platforms in the context of established guidelines for network reconstruction. While many steps of the reconstruction process have been successfully automated, some manual evaluation of the results is still required to ensure a high-quality reconstruction. Widespread adoption of these platforms by the scientific community is underway and will be further enabled by exchangeable formats across platforms.

Quantitative assessment of thermodynamic constraints on the solution space of genome-scale metabolic models.

Constraint-based methods provide powerful computational techniques to allow understanding and prediction of cellular behavior. These methods rely on physiochemical constraints to eliminate infeasible behaviors from the space of available behaviors. One such constraint is thermodynamic feasibility, the requirement that intracellular flux distributions obey the laws of thermodynamics. The past decade has seen several constraint-based methods that interpret this constraint in different ways, including those that are limited to small networks, rely on predefined reaction directions, and/or neglect the relationship between reaction free energies and metabolite concentrations. In this work, we utilize one such approach, thermodynamics-based metabolic flux analysis (TMFA), to make genome-scale, quantitative predictions about metabolite concentrations and reaction free energies in the absence of prior knowledge of reaction directions, while accounting for uncertainties in thermodynamic estimates. We applied TMFA to a genome-scale network reconstruction of Escherichia coli and examined the effect of thermodynamic constraints on the flux space. We also assessed the predictive performance of TMFA against gene essentiality and quantitative metabolomics data, under both aerobic and anaerobic, and optimal and suboptimal growth conditions. Based on these results, we propose that TMFA is a useful tool for validating phenotypes and generating hypotheses, and that additional types of data and constraints can improve predictions of metabolite concentrations.

Identification of functional differences in metabolic networks using comparative genomics and constraint-based models.

Genome-scale network reconstructions are useful tools for understanding cellular metabolism, and comparisons of such reconstructions can provide insight into metabolic differences between organisms. Recent efforts toward comparing genome-scale models have focused primarily on aligning metabolic networks at the reaction level and then looking at differences and similarities in reaction and gene content. However, these reaction comparison approaches are time-consuming and do not identify the effect network differences have on the functional states of the network. We have developed a bilevel mixed-integer programming approach, CONGA, to identify functional differences between metabolic networks by comparing network reconstructions aligned at the gene level. We first identify orthologous genes across two reconstructions and then use CONGA to identify conditions under which differences in gene content give rise to differences in metabolic capabilities. By seeking genes whose deletion in one or both models disproportionately changes flux through a selected reaction (e.g., growth or by-product secretion) in one model over another, we are able to identify structural metabolic network differences enabling unique metabolic capabilities. Using CONGA, we explore functional differences between two metabolic reconstructions of Escherichia coli and identify a set of reactions responsible for chemical production differences between the two models. We also use this approach to aid in the development of a genome-scale model of Synechococcus sp. PCC 7002. Finally, we propose potential antimicrobial targets in Mycobacterium tuberculosis and Staphylococcus aureus based on differences in their metabolic capabilities. Through these examples, we demonstrate that a gene-centric approach to comparing metabolic networks allows for a rapid comparison of metabolic models at a functional level. Using CONGA, we can identify differences in reaction and gene content which give rise to different functional predictions. Because CONGA provides a general framework, it can be applied to find functional differences across models and biological systems beyond those presented here.