Genome-driven evolutionary game theory helps understand the rise of metabolic interdependencies in microbial communities

Metabolite exchanges in microbial communities give rise to ecological interactions that govern ecosystem diversity and stability. It is unclear, however, how the rise of these interactions varies across metabolites and organisms. Here we address this question by integrating genome-scale models of metabolism with evolutionary game theory. Specifically, we use microbial fitness values estimated by metabolic models to infer evolutionarily stable interactions in multi-species microbial “games”. We first validate our approach using a well-characterized yeast cheater-cooperator system. We next perform over 80,000 in silico experiments to infer how metabolic interdependencies mediated by amino acid leakage in Escherichia coli vary across 189 amino acid pairs. While most pairs display shared patterns of inter-species interactions, multiple deviations are caused by pleiotropy and epistasis in metabolism. Furthermore, simulated invasion experiments reveal possible paths to obligate cross-feeding. Our study provides genomically driven insight into the rise of ecological interactions, with implications for microbiome research and synthetic ecology.We gratefully acknowledge funding from the Defense Advanced Research Projects Agency (Purchase Request No. HR0011515303, Contract No. HR0011-15-C-0091), the U.S. Department of Energy (Grants DE-SC0004962 and DE-SC0012627), the NIH (Grants 5R01DE024468 and R01GM121950), the national Science Foundation (Grants 1457695 and NSFOCE-BSF 1635070), MURI Grant W911NF-12-1-0390, the Human Frontiers Science Program (grant RGP0020/2016), and the Boston University Interdisciplinary Biomedical Research Office ARC grant on Systems Biology Approaches to Microbiome Research. We also thank Dr Kirill Korolev and members of the Segre Lab for their invaluable feedback on this work. (HR0011515303 - Defense Advanced Research Projects Agency; HR0011-15-C-0091 - Defense Advanced Research Projects Agency; DE-SC0004962 - U.S. Department of Energy; DE-SC0012627 - U.S. Department of Energy; 5R01DE024468 - NIH; R01GM121950 - NIH; 1457695 - national Science Foundation; NSFOCE-BSF 1635070 - national Science Foundation; W911NF-12-1-0390 - MURI; RGP0020/2016 - Human Frontiers Science Program; Boston University Interdisciplinary Biomedical Research Office ARC)Published versio

Modulation of functional network properties in major depressive disorder following electroconvulsive therapy (ECT): a resting-state EEG analysis

Electroconvulsive therapy (ECT) is a highly effective neuromodulatory intervention for treatment-resistant major depressive disorder (MDD). Presently, however, understanding of its neurophysiological effects remains incomplete. In the present study, we utilised resting-state electroencephalography (RS-EEG) to explore changes in functional connectivity, network topology, and spectral power elicited by an acute open-label course of ECT in a cohort of 23 patients with treatment-resistant MDD. RS-EEG was recorded prior to commencement of ECT and again within 48 h following each patient’s final treatment session. Our results show that ECT was able to enhance connectivity within lower (delta and theta) frequency bands across subnetworks largely confined to fronto-central channels, while, conversely, more widespread subnetworks of reduced connectivity emerged within faster (alpha and beta) bands following treatment. Graph-based topological analyses revealed changes in measures of functional segregation (clustering coefficient), integration (characteristic path length), and small-world architecture following ECT. Finally, post-treatment enhancement of delta and theta spectral power was observed, which showed a positive association with the number of ECT sessions received. Overall, our findings indicate that RS-EEG can provide a sensitive measure of dynamic neural activity following ECT and highlight network-based analyses as a promising avenue for furthering mechanistic understanding of the effects of convulsive therapies

Delayed auditory feedback and transcranial direct current stimulation treatment for the enhancement of speech fluency in adults who stutter: Protocol for a randomized controlled trial

Background: Stuttering is a complex speech disorder that affects speech fluency. Recently, it has been shown that noninvasive brain stimulation may be useful to enhance the results of fluency interventions in adults who stutter. Delayed auditory feedback (DAF) is a method to enhance speech fluency in individuals who stutter. Adjunctive interventions are warranted to enhance the efficacy of this intervention. Objective: Individuals who stutter have pathological activation patterns in the primary and secondary auditory areas. Consequently, in this study, we hypothesize that stimulation of these areas might be promising as an adjunctive method to fluency training via DAF to enhance speech therapy success in individuals with a stutter. We will systematically test this hypothesis in this study. Methods: This study is designed as a randomized, double-blind, sham-controlled clinical trial. All participants will receive DAF. The intervention group will additionally receive real transcranial direct current stimulation, while the control group will be exposed to sham stimulation. The assignment of the participants to one of these groups will be randomized. Before starting the treatment program, 2 preintervention assessments will be conducted to determine the severity of stuttering. Once these assessments are completed, each subject will participate in 6 intervention sessions. Postintervention assessments will be carried out immediately and 1 week after the last intervention session. Subsequently, to explore the long-term stability of the treatment results, the outcome parameters will be obtained in follow-up assessments 6 weeks after the treatment. The primary outcome measurement�the percentage of stuttered syllables�will be calculated in pre-, post-, and follow-up assessments; the secondary outcomes will be the scores of the following questionnaires: the Stuttering Severity Instrument�Fourth Edition and the Overall Assessment of the Speaker�s Experience of Stuttering. Results: This protocol was funded in 2019 and approved by the Research Ethics Committee of the Iran University of Medical Sciences in June 2019. Data collection started in October 2019. As of February 2020, we have enrolled 30 participants. We expect data analysis to be completed in April 2020, and results will be published in summer 2020. Conclusions: We anticipate that this study will show an adjunctive effect of transcranial direct current stimulation, when combined with DAF, on stuttering. This should include not only a reduction in the percentage of stuttered syllables but also improved physical behavior and quality of life in adults who stutter. © 2020 Journal of Medical Internet Research. All rights reserved

Universality of Level Spacing Distributions in Classical Chaos

We suggest that random matrix theory applied to a classical action matrix can be used in classical physics to distinguish chaotic from non-chaotic behavior. We consider the 2-D stadium billiard system as well as the 2-D anharmonic and harmonic oscillator. By unfolding of the spectrum of such matrix we compute the level spacing distribution, the spectral auto-correlation and spectral rigidity. We observe Poissonian behavior in the integrable case and Wignerian behavior in the chaotic case. We present numerical evidence that the action matrix of the stadium billiard displays GOE behavior and give an explanation for it. The findings present evidence for universality of level fluctuations - known from quantum chaos - also to hold in classical physics

Genome-scale gene/reaction essentiality and synthetic lethality analysis

Synthetic lethals are to pairs of non-essential genes whose simultaneous deletion prohibits growth. One can extend the concept of synthetic lethality by considering gene groups of increasing size where only the simultaneous elimination of all genes is lethal, whereas individual gene deletions are not. We developed optimization-based procedures for the exhaustive and targeted enumeration of multi-gene (and by extension multi-reaction) lethals for genome-scale metabolic models. Specifically, these approaches are applied to iAF1260, the latest model of Escherichia coli, leading to the complete identification of all double and triple gene and reaction synthetic lethals as well as the targeted identification of quadruples and some higher-order ones. Graph representations of these synthetic lethals reveal a variety of motifs ranging from hub-like to highly connected subgraphs providing a birds-eye view of the avenues available for redirecting metabolism and uncovering complex patterns of gene utilization and interdependence. The procedure also enables the use of falsely predicted synthetic lethals for metabolic model curation. By analyzing the functional classifications of the genes involved in synthetic lethals, we reveal surprising connections within and across clusters of orthologous group functional classifications

Impact of stoichiometry representation on simulation of genotype-phenotype relationships in metabolic networks.

<div><p>Genome-scale metabolic networks provide a comprehensive structural framework for modeling genotype-phenotype relationships through flux simulations. The solution space for the metabolic flux state of the cell is typically very large and optimization-based approaches are often necessary for predicting the active metabolic state under specific environmental conditions. The objective function to be used in such optimization algorithms is directly linked with the biological hypothesis underlying the model and therefore it is one of the most relevant parameters for successful modeling. Although linear combination of selected fluxes is widely used for formulating metabolic objective functions, we show that the resulting optimization problem is sensitive towards stoichiometry representation of the metabolic network. This undesirable sensitivity leads to different simulation results when using numerically different but biochemically equivalent stoichiometry representations and thereby makes biological interpretation intrinsically subjective and ambiguous. We hereby propose a new method, Minimization of Metabolites Balance (MiMBl), which decouples the artifacts of stoichiometry representation from the formulation of the desired objective functions, by casting objective functions using metabolite turnovers rather than fluxes. By simulating perturbed metabolic networks, we demonstrate that the use of stoichiometry representation independent algorithms is fundamental for unambiguously linking modeling results with biological interpretation. For example, MiMBl allowed us to expand the scope of metabolic modeling in elucidating the mechanistic basis of several genetic interactions in <em>Saccharomyces cerevisiae</em>.</p> </div

Conserved properties of dendritic trees in four cortical interneuron subtypes

Dendritic trees influence synaptic integration and neuronal excitability, yet appear to develop in rather arbitrary patterns. Using electron microscopy and serial reconstructions, we analyzed the dendritic trees of four morphologically distinct neocortical interneuron subtypes to reveal two underlying organizational principles common to all. First, cross-sectional areas at any given point within a dendrite were proportional to the summed length of all dendritic segments distal to that point. Consistent with this observation, total cross-sectional area was almost perfectly conserved at bifurcation points. Second, dendritic cross-sections became progressively more elliptical at more proximal, larger diameter, dendritic locations. Finally, computer simulations revealed that these conserved morphological features limit distance dependent filtering of somatic EPSPs and facilitate distribution of somatic depolarization into all dendritic compartments. Because these features were shared by all interneurons studied, they may represent common organizational principles underlying the otherwise diverse morphology of dendritic trees

Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data

<p>Abstract</p> <p>Background</p> <p><it>Saccharomyces cerevisiae </it>is the first eukaryotic organism for which a multi-compartment genome-scale metabolic model was constructed. Since then a sequence of improved metabolic reconstructions for yeast has been introduced. These metabolic models have been extensively used to elucidate the organizational principles of yeast metabolism and drive yeast strain engineering strategies for targeted overproductions. They have also served as a starting point and a benchmark for the reconstruction of genome-scale metabolic models for other eukaryotic organisms. In spite of the successive improvements in the details of the described metabolic processes, even the recent yeast model (i.e., <it>i</it>MM904) remains significantly less predictive than the latest <it>E. coli </it>model (i.e., <it>i</it>AF1260). This is manifested by its significantly lower specificity in predicting the outcome of grow/no grow experiments in comparison to the <it>E. coli </it>model.</p> <p>Results</p> <p>In this paper we make use of the automated GrowMatch procedure for restoring consistency with single gene deletion experiments in yeast and extend the procedure to make use of synthetic lethality data using the genome-scale model <it>i</it>MM904 as a basis. We identified and vetted using literature sources 120 distinct model modifications including various regulatory constraints for minimal and YP media. The incorporation of the suggested modifications led to a substantial increase in the fraction of correctly predicted lethal knockouts (i.e., specificity) from 38.84% (87 out of 224) to 53.57% (120 out of 224) for the minimal medium and from 24.73% (45 out of 182) to 40.11% (73 out of 182) for the YP medium. Synthetic lethality predictions improved from 12.03% (16 out of 133) to 23.31% (31 out of 133) for the minimal medium and from 6.96% (8 out of 115) to 13.04% (15 out of 115) for the YP medium.</p> <p>Conclusions</p> <p>Overall, this study provides a roadmap for the computationally driven correction of multi-compartment genome-scale metabolic models and demonstrates the value of synthetic lethals as curation agents.</p