Comparative study of bacterial community dynamics in different soils following application of the herbicide atrazine.

Microbial communities in cultivated soils control the fate of pollutants associated with agricultural practice. The present study was designed to explore the response of bacterial communities to the application of the widely-used herbicide atrazine in three different crop fields that differ significantly in their physicochemical structure and nutritional content: the nutrient-rich (with relatively high carbon and nitrogen content) Newe Yaar (NY) and Ha-Ogen (HO) soils and the nutrient-poor, sandy Sde-Eliyahu (SE) soil. The 16 S rRNA gene amplicon sequencing revealed the nutrient poor HO soil differs in its response to atrazine in comparison to the two nutrient-rich soils both in the shortest persistence of atrazine and its effect on community structure and composition. Potential reported bacterial degraders of atrazine such as Pseudomonas, Clostridium and Bacillus were more abundant in contaminated sandy/poor soils (HO) whereas bacteria known for nitrogen cycling such as Azospirillum, Sinorhizobium, Nitrospira and Azohydromonas were significantly more abundant in the nutrient rich contaminated SE soils. No significant increase of potential indigenous degrader Arthrobacter was detected in SE and NY soils whereas a significant increase was recorded with HO soils. An overall shift in bacterial community composition following atrazine application was observed only in the nutrient poor soil. Understanding atrazine persistence and microbiome response to its application of in dependence with soil types serve the design of precision application strategies.

Diet-induced changes of redox potential underlie compositional shifts in the rumen archaeal community.

Dietary changes are known to affect gut community structure, but questions remain about the mechanisms by which diet induces shifts in microbiome membership. Here, we addressed these questions in the rumen microbiome ecosystem - a complex microbial community that resides in the upper digestive tract of ruminant animals and is responsible for the degradation of the ingested plant material. Our dietary intervention experiments revealed that diet affects the most abundant taxa within the microbiome and that a specific group of methanogenic archaea of the order Methanomicrobiales is highly sensitive to its changes. Using metabolomic analyses together with in vitro microbiology approaches and whole-genome sequencing of Methanomicrobium mobile, a key species within this group, we identified that redox potential changes with diet and is the main factor that causes these dietary induced alternations in this taxa's abundance. Our genomic analysis suggests that the redox potential effect stems from a reduced number of anti-reactive oxygen species proteins coded in this taxon's genome. Our study highlights redox potential as a pivotal factor that could serve as a sculpturing force of community assembly within anaerobic gut microbial communities.

Systems-Wide Prediction of Enzyme Promiscuity Reveals a New Underground Alternative Route for Pyridoxal 5'-Phosphate Production in E. coli.

Recent insights suggest that non-specific and/or promiscuous enzymes are common and active across life. Understanding the role of such enzymes is an important open question in biology. Here we develop a genome-wide method, PROPER, that uses a permissive PSI-BLAST approach to predict promiscuous activities of metabolic genes. Enzyme promiscuity is typically studied experimentally using multicopy suppression, in which over-expression of a promiscuous 'replacer' gene rescues lethality caused by inactivation of a 'target' gene. We use PROPER to predict multicopy suppression in Escherichia coli, achieving highly significant overlap with published cases (hypergeometric p = 4.4e-13). We then validate three novel predicted target-replacer gene pairs in new multicopy suppression experiments. We next go beyond PROPER and develop a network-based approach, GEM-PROPER, that integrates PROPER with genome-scale metabolic modeling to predict promiscuous replacements via alternative metabolic pathways. GEM-PROPER predicts a new indirect replacer (thiG) for an essential enzyme (pdxB) in production of pyridoxal 5'-phosphate (the active form of Vitamin B6), which we validate experimentally via multicopy suppression. We perform a structural analysis of thiG to determine its potential promiscuous active site, which we validate experimentally by inactivating the pertaining residues and showing a loss of replacer activity. Thus, this study is a successful example where a computational investigation leads to a network-based identification of an indirect promiscuous replacement of a key metabolic enzyme, which would have been extremely difficult to identify directly.

Computational evaluation of cellular metabolic costs successfully predicts genes whose expression is deleterious.

Gene suppression and overexpression are both fundamental tools in linking genotype to phenotype in model organisms. Computational methods have proven invaluable in studying and predicting the deleterious effects of gene deletions, and yet parallel computational methods for overexpression are still lacking. Here, we present Expression-Dependent Gene Effects (EDGE), an in silico method that can predict the deleterious effects resulting from overexpression of either native or foreign metabolic genes. We first test and validate EDGE's predictive power in bacteria through a combination of small-scale growth experiments that we performed and analysis of extant large-scale datasets. Second, a broad cross-species analysis, ranging from microorganisms to multiple plant and human tissues, shows that genes that EDGE predicts to be deleterious when overexpressed are indeed typically down-regulated. This reflects a universal selection force keeping the expression of potentially deleterious genes in check. Third, EDGE-based analysis shows that cancer genetic reprogramming specifically suppresses genes whose overexpression impedes proliferation. The magnitude of this suppression is large enough to enable an almost perfect distinction between normal and cancerous tissues based solely on EDGE results. We expect EDGE to advance our understanding of human pathologies associated with up-regulation of particular transcripts and to facilitate the utilization of gene overexpression in metabolic engineering.

A novel nutritional predictor links microbial fastidiousness with lowered ubiquity, growth rate, and cooperativeness.

Understanding microbial nutritional requirements is a key challenge in microbiology. Here we leverage the recent availability of thousands of automatically generated genome-scale metabolic models to develop a predictor of microbial minimal medium requirements, which we apply to thousands of species to study the relationship between their nutritional requirements and their ecological and genomic traits. We first show that nutritional requirements are more similar among species that co-habit many ecological niches. We then reveal three fundamental characteristics of microbial fastidiousness (i.e., complex and specific nutritional requirements): (1) more fastidious microorganisms tend to be more ecologically limited; (2) fastidiousness is positively associated with smaller genomes and smaller metabolic networks; and (3) more fastidious species grow more slowly and have less ability to cooperate with other species than more metabolically versatile organisms. These associations reflect the adaptation of fastidious microorganisms to unique niches with few cohabitating species. They also explain how non-fastidious species inhabit many ecological niches with high abundance rates. Taken together, these results advance our understanding microbial nutrition on a large scale, by presenting new nutrition-related associations that govern the distribution of microorganisms in nature.

The MISOPRED score: Development and validation of a clinical scoring system to predict the effectiveness of Misoprostol treatment for early pregnancy loss.

BACKGROUND: Misoprostol treatment for early pregnancy loss has varied success demonstrated in previous studies. Incorporating predictors in a single clinical scoring system would be highly beneficial in clinical practice. OBJECTIVE: To develop and evaluate the accuracy of a scoring system to predict misoprostol treatment outcomes for managing early pregnancy loss. STUDY DESIGN: Retrospective cohort and validation study. METHODS: Patients discharged from the gynecologic emergency department from 2013 to 2016, diagnosed with early pregnancy loss, who were treated with 800 mcg misoprostol, administrated vaginally were included. All were sonographically reevaluated within 48-72 hours. Patients in whom the gestational sac was not expelled or with endometrial lining >30 mm were offered a repeat dose and returned for reevaluation after seven days. A successful response was defined as complete expulsion. Clinical data were reviewed to identify predictors for successful responses. The scoring system was then retrospectively evaluated on a second cohort to evaluate its accuracy. Multivariate logistic regression was performed to identify factors most predictive of treatment response. RESULTS: The development cohort included 126 patients. Six factors were found to be most predictive of misoprostol treatment effectiveness: nulliparity, prior complete spontaneous abortion, gestational age, vaginal bleeding, abdominal pain, and mean sac diameter, yielding a score of 0-8 (the MISOPRED score), where 8 represents the highest-likelihood of success. The score was validated retrospectively with 119 participants. Successful response in the group with the lowest likelihood score (score 0-3) was 9%, compared with 82% in the highest likelihood score group (score 7-8). Using the MISOPRED score, approximately 15% of patients previously planned to receive misoprostol treatment can be referred for surgical management. CONCLUSIONS: MISOPRED score can be utilized as an adjunct tool for clinical decision-making in cases of Early pregnancy loss. To our knowledge, this is the first scoring system suggested to predict the success rate in these cases.

Modeling-Guided Amendments Lead to Enhanced Biodegradation in Soil.

Extensive use of agrochemicals is emerging as a serious environmental issue coming at the cost of the pollution of soil and water resources. Bioremediation techniques such as biostimulation are promising strategies used to remove pollutants from agricultural soils by supporting the indigenous microbial degraders. Though considered cost-effective and eco-friendly, the success rate of these strategies typically varies, and consequently, they are rarely integrated into commercial agricultural practices. In the current study, we applied metabolic-based community-modeling approaches for promoting realistic in terra solutions by simulation-based prioritization of alternative supplements as potential biostimulants, considering a collection of indigenous bacteria. Efficacy of biostimulants as enhancers of the indigenous degrader Paenarthrobacter was ranked through simulation and validated in pot experiments. A two-dimensional simulation matrix predicting the effect of different biostimulants on additional potential indigenous degraders (Pseudomonas, Clostridium, and Geobacter) was crossed with experimental observations. The overall ability of the models to predict the compounds that act as taxa-selective stimulants indicates that computational algorithms can guide the manipulation of the soil microbiome in situ and provides an additional step toward the educated design of biostimulation strategies. IMPORTANCE Providing the food requirements of a growing population comes at the cost of intensive use of agrochemicals, including pesticides. Native microbial soil communities are considered key players in the degradation of such exogenous substances. Manipulating microbial activity toward an optimized outcome in efficient biodegradation processes conveys a promise of maintaining intensive yet sustainable agriculture. Efficient strategies for harnessing the native microbiome require the development of approaches for processing big genomic data. Here, we pursued metabolic modeling for promoting realistic in terra solutions by simulation-based prioritization of alternative supplements as potential biostimulants, considering a collection of indigenous bacteria. Our genomic-based predictions point at strategies for optimizing biodegradation by the native community. Developing a systematic, data-guided understanding of metabolite-driven targeted enhancement of selected microorganisms lays the foundation for the design of ecologically sound methods for optimizing microbiome functioning.

Strategies for Enhancing in vitro Degradation of Linuron by Variovorax sp. Strain SRS 16 Under the Guidance of Metabolic Modeling.

Phenyl urea herbicides are being extensively used for weed control in both agricultural and non-agricultural applications. Linuron is one of the key herbicides in this family and is in wide use. Like other phenyl urea herbicides, it is known to have toxic effects as a result of its persistence in the environment. The natural removal of linuron from the environment is mainly carried through microbial biodegradation. Some microorganisms have been reported to mineralize linuron completely and utilize it as a carbon and nitrogen source. Variovorax sp. strain SRS 16 is one of the known efficient degraders with a recently sequenced genome. The genomic data provide an opportunity to use a genome-scale model for improving biodegradation. The aim of our study is the construction of a genome-scale metabolic model following automatic and manual protocols and its application for improving its metabolic potential through iterative simulations. Applying flux balance analysis (FBA), growth and degradation performances of SRS 16 in different media considering the influence of selected supplements (potential carbon and nitrogen sources) were simulated. Outcomes are predictions for the suitable media modification, allowing faster degradation of linuron by SRS 16. Seven metabolites were selected for in vitro validation of the predictions through laboratory experiments confirming the degradation-promoting effect of specific amino acids (glutamine and asparagine) on linuron degradation and SRS 16 growth. Overall, simulations are shown to be efficient in predicting the degradation potential of SRS 16 in the presence of specific supplements. The generated information contributes to the understanding of the biochemistry of linuron degradation and can be further utilized for the development of new cleanup solutions without any genetic manipulation.

Genome-scale reconstruction of Paenarthrobacter aurescens TC1 metabolic model towards the study of atrazine bioremediation.

Atrazine is an herbicide and a pollutant of great environmental concern that is naturally biodegraded by microbial communities. Paenarthrobacter aurescens TC1 is one of the most studied degraders of this herbicide. Here, we developed a genome scale metabolic model for P. aurescens TC1, iRZ1179, to study the atrazine degradation process at organism level. Constraint based flux balance analysis and time dependent simulations were used to explore the organism's phenotypic landscape. Simulations aimed at designing media optimized for supporting growth and enhancing degradation, by passing the need in strain design via genetic modifications. Growth and degradation simulations were carried with more than 100 compounds consumed by P. aurescens TC1. In vitro validation confirmed the predicted classification of different compounds as efficient, moderate or poor stimulators of growth. Simulations successfully captured previous reports on the use of glucose and phosphate as bio-stimulators of atrazine degradation, supported by in vitro validation. Model predictions can go beyond supplementing the medium with a single compound and can predict the growth outcomes for higher complexity combinations. Hence, the analysis demonstrates that the exhaustive power of the genome scale metabolic reconstruction allows capturing complexities that are beyond common biochemical expertise and knowledge and further support the importance of computational platforms for the educated design of complex media. The model presented here can potentially serve as a predictive tool towards achieving optimal biodegradation efficiencies and for the development of ecologically friendly solutions for pollutant degradation.

Prebiotic-Induced Anti-tumor Immunity Attenuates Tumor Growth.

Growing evidence supports the importance of gut microbiota in the control of tumor growth and response to therapy. Here, we select prebiotics that can enrich bacterial taxa that promote anti-tumor immunity. Addition of the prebiotics inulin or mucin to the diet of C57BL/6 mice induces anti-tumor immune responses and inhibition of BRAF mutant melanoma growth in a subcutaneously implanted syngeneic mouse model. Mucin fails to inhibit tumor growth in germ-free mice, indicating that the gut microbiota is required for the activation of the anti-tumor immune response. Inulin and mucin drive distinct changes in the microbiota, as inulin, but not mucin, limits tumor growth in syngeneic mouse models of colon cancer and NRAS mutant melanoma and enhances the efficacy of a MEK inhibitor against melanoma while delaying the emergence of drug resistance. We highlight the importance of gut microbiota in anti-tumor immunity and the potential therapeutic role for prebiotics in this process.

Modeling microbial communities from atrazine contaminated soils promotes the development of biostimulation solutions.

Microbial communities play a vital role in biogeochemical cycles, allowing the biodegradation of a wide range of pollutants. The composition of the community and the interactions between its members affect degradation rate and determine the identity of the final products. Here, we demonstrate the application of sequencing technologies and metabolic modeling approaches towards enhancing biodegradation of atrazine-a herbicide causing environmental pollution. Treatment of agriculture soil with atrazine is shown to induce significant changes in community structure and functional performances. Genome-scale metabolic models were constructed for Arthrobacter, the atrazine degrader, and four other non-atrazine degrading species whose relative abundance in soil was changed following exposure to the herbicide. By modeling community function we show that consortia including the direct degrader and non-degrader differentially abundant species perform better than Arthrobacter alone. Simulations predict that growth/degradation enhancement is derived by metabolic exchanges between community members. Based on simulations we designed endogenous consortia optimized for enhanced degradation whose performances were validated in vitro and biostimulation strategies that were tested in pot experiments. Overall, our analysis demonstrates that understanding community function in its wider context, beyond the single direct degrader perspective, promotes the design of biostimulation strategies.

Gut microbiota dependent anti-tumor immunity restricts melanoma growth in Rnf5-/- mice.

Accumulating evidence points to an important role for the gut microbiome in anti-tumor immunity. Here, we show that altered intestinal microbiota contributes to anti-tumor immunity, limiting tumor expansion. Mice lacking the ubiquitin ligase RNF5 exhibit attenuated activation of the unfolded protein response (UPR) components, which coincides with increased expression of inflammasome components, recruitment and activation of dendritic cells and reduced expression of antimicrobial peptides in intestinal epithelial cells. Reduced UPR expression is also seen in murine and human melanoma tumor specimens that responded to immune checkpoint therapy. Co-housing of Rnf5-/- and WT mice abolishes the anti-tumor immunity and tumor inhibition phenotype, whereas transfer of 11 bacterial strains, including B. rodentium, enriched in Rnf5-/- mice, establishes anti-tumor immunity and restricts melanoma growth in germ-free WT mice. Altered UPR signaling, exemplified in Rnf5-/- mice, coincides with altered gut microbiota composition and anti-tumor immunity to control melanoma growth.

Harnessing the landscape of microbial culture media to predict new organism-media pairings.

Culturing microorganisms is a critical step in understanding and utilizing microbial life. Here we map the landscape of existing culture media by extracting natural-language media recipes into a Known Media Database (KOMODO), which includes >18,000 strain-media combinations, >3300 media variants and compound concentrations (the entire collection of the Leibniz Institute DSMZ repository). Using KOMODO, we show that although media are usually tuned for individual strains using biologically common salts, trace metals and vitamins/cofactors are the most differentiating components between defined media of strains within a genus. We leverage KOMODO to predict new organism-media pairings using a transitivity property (74% growth in new in vitro experiments) and a phylogeny-based collaborative filtering tool (83% growth in new in vitro experiments and stronger growth on predicted well-scored versus poorly scored media). These resources are integrated into a web-based platform that predicts media given an organism's 16S rDNA sequence, facilitating future cultivation efforts.

Improved evidence-based genome-scale metabolic models for maize leaf, embryo, and endosperm.

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

Maximal sum of metabolic exchange fluxes outperforms biomass yield as a predictor of growth rate of microorganisms.

Growth rate has long been considered one of the most valuable phenotypes that can be measured in cells. Aside from being highly accessible and informative in laboratory cultures, maximal growth rate is often a prime determinant of cellular fitness, and predicting phenotypes that underlie fitness is key to both understanding and manipulating life. Despite this, current methods for predicting microbial fitness typically focus on yields [e.g., predictions of biomass yield using GEnome-scale metabolic Models (GEMs)] or notably require many empirical kinetic constants or substrate uptake rates, which render these methods ineffective in cases where fitness derives most directly from growth rate. Here we present a new method for predicting cellular growth rate, termed SUMEX, which does not require any empirical variables apart from a metabolic network (i.e., a GEM) and the growth medium. SUMEX is calculated by maximizing the SUM of molar EXchange fluxes (hence SUMEX) in a genome-scale metabolic model. SUMEX successfully predicts relative microbial growth rates across species, environments, and genetic conditions, outperforming traditional cellular objectives (most notably, the convention assuming biomass maximization). The success of SUMEX suggests that the ability of a cell to catabolize substrates and produce a strong proton gradient enables fast cell growth. Easily applicable heuristics for predicting growth rate, such as what we demonstrate with SUMEX, may contribute to numerous medical and biotechnological goals, ranging from the engineering of faster-growing industrial strains, modeling of mixed ecological communities, and the inhibition of cancer growth.

Glycan degradation (GlyDeR) analysis predicts mammalian gut microbiota abundance and host diet-specific adaptations.

UNLABELLED: Glycans form the primary nutritional source for microbes in the human gut, and understanding their metabolism is a critical yet understudied aspect of microbiome research. Here, we present a novel computational pipeline for modeling glycan degradation (GlyDeR) which predicts the glycan degradation potency of 10,000 reference glycans based on either genomic or metagenomic data. We first validated GlyDeR by comparing degradation profiles for genomes in the Human Microbiome Project against KEGG reaction annotations. Next, we applied GlyDeR to the analysis of human and mammalian gut microbial communities, which revealed that the glycan degradation potential of a community is strongly linked to host diet and can be used to predict diet with higher accuracy than sequence data alone. Finally, we show that a microbe's glycan degradation potential is significantly correlated (R = 0.46) with its abundance, with even higher correlations for potential pathogens such as the class Clostridia (R = 0.76). GlyDeR therefore represents an important tool for advancing our understanding of bacterial metabolism in the gut and for the future development of more effective prebiotics for microbial community manipulation. IMPORTANCE: The increased availability of high-throughput sequencing data has positioned the gut microbiota as a major new focal point for biomedical research. However, despite the expenditure of huge efforts and resources, sequencing-based analysis of the microbiome has uncovered mostly associative relationships between human health and diet, rather than a causal, mechanistic one. In order to utilize the full potential of systems biology approaches, one must first characterize the metabolic requirements of gut bacteria, specifically, the degradation of glycans, which are their primary nutritional source. We developed a computational framework called GlyDeR for integrating expert knowledge along with high-throughput data to uncover important new relationships within glycan metabolism. GlyDeR analyzes particular bacterial (meta)genomes and predicts the potency by which they degrade a variety of different glycans. Based on GlyDeR, we found a clear connection between microbial glycan degradation and human diet, and we suggest a method for the rational design of novel prebiotics.

Competitive and cooperative metabolic interactions in bacterial communities.

Revealing the ecological principles that shape communities is a major challenge of the post-genomic era. To date, a systematic approach for describing inter-species interactions has been lacking. Here we independently predict the competitive and cooperative potential between 6,903 bacterial pairs derived from a collection of 118 species' metabolic models. We chart an intricate association between competition and cooperation indicating that the cooperative potential is maximized at moderate levels of resource overlap. Utilizing ecological data from 2,801 samples, we explore the associations between bacterial interactions and coexistence patterns. The high level of competition observed between species with mutual-exclusive distribution patterns supports the role of competition in community assembly. Cooperative interactions are typically unidirectional with no obvious benefit to the giver. However, within their natural communities, bacteria typically form close cooperative loops resulting in indirect benefit to all species involved. These findings are important for the future design of consortia optimized towards bioremediation and bio-production applications.