Metabolic networks of Sodalis glossinidius: a systems biology approach to reductive evolution

Growing Ungrowable Bacteria: Overview and Perspectives on Insect Symbiont Culturability

Insects are often involved in endosymbiosis, that is, the housing of symbiotic microbes within their tissues or within their cells. Endosymbionts are a major driving force in insects' evolution, because they dramatically affect their host physiology and allow them to adapt to new niches, for example, by complementing their diet or by protecting them against pathogens. Endosymbiotic bacteria are, however, fastidious and therefore difficult to manipulate outside of their hosts, especially intracellular species. The coevolution between hosts and endosymbionts leads to alterations in the genomes of endosymbionts, limiting their ability to cope with changing environments. Consequently, few insect endosymbionts are culturable in vitro and genetically tractable, making functional genetics studies impracticable on most endosymbiotic bacteria. However, recently, major progress has been made in manipulating several intracellular endosymbiont species in vitro, leading to astonishing discoveries on their physiology and the way they interact with their host. This review establishes a comprehensive picture of the in vitro tractability of insect endosymbiotic bacteria and addresses the reason why most species are not culturable. By compiling and discussing the latest developments in the design of custom media and genetic manipulation protocols, it aims at providing new leads to expand the range of tractable endosymbionts and foster genetic research on these models.

Simulating the evolutionary trajectories of metabolic pathways for insect symbionts in the genus Sodalis

Insect-bacterial symbioses are ubiquitous, but there is still much to uncover about how these relationships establish, persist and evolve. The tsetse endosymbiont Sodalis glossinidius displays intriguing metabolic adaptations to its microenvironment, but the process by which this relationship evolved remains to be elucidated. The recent chance discovery of the free-living species of the genus Sodalis, Sodalis praecaptivus, provides a serendipitous starting point from which to investigate the evolution of this symbiosis. Here, we present a flux balance model for S. praecaptivus and empirically verify its predictions. Metabolic modelling is used in combination with a multi-objective evolutionary algorithm to explore the trajectories that S. glossinidius may have undertaken from this starting point after becoming internalized. The order in which key genes are lost is shown to influence the evolved populations, providing possible targets for future in vitro genetic manipulation. This method provides a detailed perspective on possible evolutionary trajectories for S. glossinidius in this fundamental process of evolutionary and ecological change.

Insect-Symbiont Gene Expression in the Midgut Bacteriocytes of a Blood-Sucking Parasite

Animals interact with a diverse array of both beneficial and detrimental microorganisms. In insects, these symbioses in many cases allow feeding on nutritionally unbalanced diets. It is, however, still not clear how are obligate symbioses maintained at the cellular level for up to several hundred million years. Exact mechanisms driving host-symbiont interactions are only understood for a handful of model species and data on blood-feeding hosts with intracellular bacteria are particularly scarce. Here, we analyzed interactions between an obligately blood-sucking parasite of sheep, the louse fly Melophagus ovinus, and its obligate endosymbiont, Arsenophonus melophagi. We assembled a reference transcriptome for the insect host and used dual RNA-Seq with five biological replicates to compare expression in the midgut cells specialized for housing symbiotic bacteria (bacteriocytes) to the rest of the gut (foregut-hindgut). We found strong evidence for the importance of zinc in the system likely caused by symbionts using zinc-dependent proteases when acquiring amino acids, and for different immunity mechanisms controlling the symbionts than in closely related tsetse flies. Our results show that cellular and nutritional interactions between this blood-sucking insect and its symbionts are less intimate than what was previously found in most plant-sap sucking insects. This finding is likely interconnected to several features observed in symbionts in blood-sucking arthropods, particularly their midgut intracellular localization, intracytoplasmic presence, less severe genome reduction, and relatively recent associations caused by frequent evolutionary losses and replacements.

Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

The Tsetse Fly Displays an Attenuated Immune Response to Its Secondary Symbiont, Sodalis glossinidius

Sodalis glossinidius, a vertically transmitted facultative symbiont of the tsetse fly, is a bacterium in the early/intermediate state of its transition toward symbiosis, representing an important model for investigating how the insect host immune defense response is regulated to allow endosymbionts to establish a chronic infection within their hosts without being eliminated. In this study, we report on the establishment of a tsetse fly line devoid of S. glossinidius only, allowing us to experimentally investigate (i) the complex immunological interactions between a single bacterial species and its host, (ii) how the symbiont population is kept under control, and (iii) the impact of the symbiont on the vector competence of the tsetse fly to transmit the sleeping sickness parasite. Comparative transcriptome analysis showed no difference in the expression of genes involved in innate immune processes between symbiont-harboring (Gmm^Sod+) and S. glossinidius-free (Gmm^Sod–) flies. Re-exposure of (Gmm^Sod–) flies to the endosymbiotic bacterium resulted in a moderate immune response, whereas exposure to pathogenic E. coli or to a close non-insect associated relative of S. glossinidius, i.e., S. praecaptivus, resulted in full immune activation. We also showed that S. glossinidius densities are not affected by experimental activation or suppression of the host immune system, indicating that S. glossinidius is resistant to mounted immune attacks and that the host immune system does not play a major role in controlling S. glossinidius proliferation. Finally, we demonstrate that the absence or presence of S. glossinidius in the tsetse fly does not alter its capacity to mount an immune response to pathogens nor does it affect the fly’s susceptibility toward trypanosome infection.

A Tale of Three Species: Adaptation of Sodalis glossinidius to Tsetse Biology, Wigglesworthia Metabolism, and Host Diet

The tsetse fly is the insect vector for the Trypanosoma brucei parasite, the causative agent of human African trypanosomiasis. The colonization and spread of the trypanosome correlate positively with the presence of a secondary symbiotic bacterium, Sodalis glossinidius The metabolic requirements and interactions of the bacterium with its host are poorly understood, and herein we describe a metabolic model of S. glossinidius metabolism. The model enabled the design and experimental verification of a defined medium that supports S. glossinidius growth ex vivo This has been used subsequently to analyze in vitro aspects of S. glossinidius metabolism, revealing multiple unique adaptations of the symbiont to its environment. Continued dependence on a sugar, and the importance of the chitin monomer N-acetyl-d-glucosamine as a carbon and energy source, suggests adaptation to host-derived molecules. Adaptation to the amino acid-rich blood diet is revealed by a strong dependence on l-glutamate as a source of carbon and nitrogen and by the ability to rescue a predicted l-arginine auxotrophy. Finally, the selective loss of thiamine biosynthesis, a vitamin provided to the host by the primary symbiont Wigglesworthia glossinidia, reveals an intersymbiont dependence. The reductive evolution of S. glossinidius to exploit environmentally derived metabolites has resulted in multiple weaknesses in the metabolic network. These weaknesses may become targets for reagents that inhibit S. glossinidius growth and aid the reduction of trypanosomal transmission.IMPORTANCE Human African trypanosomiasis is caused by the Trypanosoma brucei parasite. The tsetse fly vector is of interest for its potential to prevent disease spread, as it is essential for T. brucei life cycle progression and transmission. The tsetse's mutualistic endosymbiont Sodalis glossinidius has a link to trypanosome establishment, providing a disease control target. Here, we describe a new, experimentally verified model of S. glossinidius metabolism. This model has enabled the development of a defined growth medium that was used successfully to test aspects of S. glossinidius metabolism. We present S. glossinidius as uniquely adapted to life in the tsetse, through its reliance on the blood diet and host-derived sugars. Additionally, S. glossinidius has adapted to the tsetse's obligate symbiont Wigglesworthia glossinidia by scavenging a vitamin it produces for the insect. This work highlights the use of metabolic modeling to design defined growth media for symbiotic bacteria and may provide novel inhibitory targets to block trypanosome transmission.

The Cost of Metabolic Interactions in Symbioses between Insects and Bacteria with Reduced Genomes

Current understanding of many animal-microbial symbioses involving unculturable bacterial symbionts with much-reduced genomes derives almost entirely from nonquantitative inferences from genome data. To overcome this limitation, we reconstructed multipartner metabolic models that quantify both the metabolic fluxes within and between three xylem-feeding insects and their bacterial symbionts. This revealed near-complete metabolic segregation between cooccurring bacterial symbionts, despite extensive metabolite exchange between each symbiont and the host, suggestive of strict host controls over the metabolism of its symbionts. We extended the model analysis to investigate metabolic costs. The positive relationship between symbiont genome size and the metabolic cost incurred by the host points to fitness benefits to the host of bearing symbionts with small genomes. The multicompartment metabolic models developed here can be applied to other symbioses that are not readily tractable to experimental approaches.

Various intracellular bacterial symbionts that provide their host with essential nutrients have much-reduced genomes, attributed largely to genomic decay and relaxed selection. To obtain quantitative estimates of the metabolic function of these bacteria, we reconstructed genome- and transcriptome-informed metabolic models of three xylem-feeding insects that bear two bacterial symbionts with complementary metabolic functions: a primary symbiont, Sulcia, that has codiversified with the insects, and a coprimary symbiont of distinct taxonomic origin and with different degrees of genome reduction in each insect species (Hodgkinia in a cicada, Baumannia in a sharpshooter, and Sodalis in a spittlebug). Our simulations reveal extensive bidirectional flux of multiple metabolites between each symbiont and the host, but near-complete metabolic segregation (i.e., near absence of metabolic cross-feeding) between the two symbionts, a likely mode of host control over symbiont metabolism. Genome reduction of the symbionts is associated with an increased number of host metabolic inputs to the symbiont and also reduced metabolic cost to the host. In particular, Sulcia and Hodgkinia with genomes of ≤0.3 Mb are calculated to recycle ∼30 to 80% of host-derived nitrogen to essential amino acids returned to the host, while Baumannia and Sodalis with genomes of ≥0.6 Mb recycle 10 to 15% of host nitrogen. We hypothesize that genome reduction of symbionts may be driven by selection for increased host control and reduced host costs, as well as by the stochastic process of genomic decay and relaxed selection.

Genome-wide characterization of Phytophthora infestans metabolism: a systems biology approach

Genome-scale metabolic models (GEMs) provide a functional view of the complex network of biochemical reactions in the living cell. Initially mainly applied to reconstruct the metabolism of model organisms, the availability of increasingly sophisticated reconstruction methods and more extensive biochemical databases now make it possible to reconstruct GEMs for less well-characterized organisms, and have the potential to unravel the metabolism in pathogen-host systems. Here, we present a GEM for the oomycete plant pathogen Phytophthora infestans as a first step towards an integrative model with its host. We predict the biochemical reactions in different cellular compartments and investigate the gene-protein-reaction associations in this model to obtain an impression of the biochemical capabilities of P. infestans. Furthermore, we generate life stage-specific models to place the transcriptomic changes of the genes encoding metabolic enzymes into a functional context. In sporangia and zoospores, there is an overall down-regulation, most strikingly reflected in the fatty acid biosynthesis pathway. To investigate the robustness of the GEM, we simulate gene deletions to predict which enzymes are essential for in vitro growth. This model is an essential first step towards an understanding of P. infestans and its interactions with plants as a system, which will help to formulate new hypotheses on infection mechanisms and disease prevention.

Determinism and Contingency Shape Metabolic Complementation in an Endosymbiotic Consortium

Bacterial endosymbionts and their insect hosts establish an intimate metabolic relationship. Bacteria offer a variety of essential nutrients to their hosts, whereas insect cells provide the necessary sources of matter and energy to their tiny metabolic allies. These nutritional complementations sustain themselves on a diversity of metabolite exchanges between the cell host and the reduced yet highly specialized bacterial metabolism—which, for instance, overproduces a small set of essential amino acids and vitamins. A well-known case of metabolic complementation is provided by the cedar aphid Cinara cedri that harbors two co-primary endosymbionts, Buchnera aphidicola BCc and Ca. Serratia symbiotica SCc, and in which some metabolic pathways are partitioned between different partners. Here we present a genome-scale metabolic network (GEM) for the bacterial consortium from the cedar aphid iBSCc. The analysis of this GEM allows us the confirmation of cases of metabolic complementation previously described by genome analysis (i.e., tryptophan and biotin biosynthesis) and the redefinition of an event of metabolic pathway sharing between the two endosymbionts, namely the biosynthesis of tetrahydrofolate. In silico knock-out experiments with iBSCc showed that the consortium metabolism is a highly integrated yet fragile network. We also have explored the evolutionary pathways leading to the emergence of metabolic complementation between reduced metabolisms starting from individual, complete networks. Our results suggest that, during the establishment of metabolic complementation in endosymbionts, adaptive evolution is significant in the case of tryptophan biosynthesis, whereas vitamin production pathways seem to adopt suboptimal solutions.

Natural selection drove metabolic specialization of the chromatophore in Paulinella chromatophora

Background

Genome degradation of host-restricted mutualistic endosymbionts has been attributed to inactivating mutations and genetic drift while genes coding for host-relevant functions are conserved by purifying selection. Unlike their free-living relatives, the metabolism of mutualistic endosymbionts and endosymbiont-originated organelles is specialized in the production of metabolites which are released to the host. This specialization suggests that natural selection crafted these metabolic adaptations. In this work, we analyzed the evolution of the metabolism of the chromatophore of Paulinella chromatophora by in silico modeling. We asked whether genome reduction is driven by metabolic engineering strategies resulted from the interaction with the host. As its widely known, the loss of enzyme coding genes leads to metabolic network restructuring sometimes improving the production rates. In this case, the production rate of reduced-carbon in the metabolism of the chromatophore.

Results

We reconstructed the metabolic networks of the chromatophore of P. chromatophora CCAC 0185 and a close free-living relative, the cyanobacterium Synechococcus sp. WH 5701. We found that the evolution of free-living to host-restricted lifestyle rendered a fragile metabolic network where >80% of genes in the chromatophore are essential for metabolic functionality. Despite the lack of experimental information, the metabolic reconstruction of the chromatophore suggests that the host provides several metabolites to the endosymbiont. By using these metabolites as intracellular conditions, in silico simulations of genome evolution by gene lose recover with 77% accuracy the actual metabolic gene content of the chromatophore. Also, the metabolic model of the chromatophore allowed us to predict by flux balance analysis a maximum rate of reduced-carbon released by the endosymbiont to the host. By inspecting the central metabolism of the chromatophore and the free-living cyanobacteria we found that by improvements in the gluconeogenic pathway the metabolism of the endosymbiont uses more efficiently the carbon source for reduced-carbon production. In addition, our in silico simulations of the evolutionary process leading to the reduced metabolic network of the chromatophore showed that the predicted rate of released reduced-carbon is obtained in less than 5% of the times under a process guided by random gene deletion and genetic drift. We interpret previous findings as evidence that natural selection at holobiont level shaped the rate at which reduced-carbon is exported to the host. Finally, our model also predicts that the ABC phosphate transporter (pstSACB) which is conserved in the genome of the chromatophore of P. chromatophora strain CCAC 0185 is a necessary component to release reduced-carbon molecules to the host.

Conclusion

Our evolutionary analysis suggests that in the case of Paulinella chromatophora natural selection at the holobiont level played a prominent role in shaping the metabolic specialization of the chromatophore. We propose that natural selection acted as a “metabolic engineer” by favoring metabolic restructurings that led to an increased release of reduced-carbon to the host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0947-6) contains supplementary material, which is available to authorized users.

Nature lessons: The whitefly bacterial endosymbiont is a minimal amino acid factory with unusual energetics

Reductive genome evolution is a universal phenomenon observed in endosymbiotic bacteria in insects. As the genome reduces its size and irreversibly losses coding genes, the functionalities of the cell system, including the energetics processes, are more restricted. Several energetic pathways can also be lost. How do these reduced metabolic networks sustain the energy needs of the system? Among the bacteria with reduced genomes Candidatus Portiera aleyrodidarum, obligate endosymbiont of whiteflies, represents an extreme case since lacks several key mechanisms for ATP generation. Thus, to analyze the cell energetics in this system, a genome-scale metabolic model of this endosymbiont was constructed, and its energy production capabilities dissected using stoichiometric analysis. Our results suggest that energy generation is coupled to the synthesis of essential amino acids and carotenoids, crucial metabolites in the symbiotic association. A deeper insight showed that ATP production via carotenoid synthesis is also connected with amino acid production. This unusual association of energy production with anabolism suggests that, although minimized, endosymbiont metabolic networks maintain a remarkable biosynthetic potential.

Taking the pseudo out of pseudogenes

Pseudogenes are defined as fragments of once-functional genes that have been silenced by one or more nonsense, frameshift or missense mutations. Despite continuing increases in the speed of sequencing and annotating bacterial genomes, the identification and categorisation of pseudogenes remains problematic. Even when identified, pseudogenes are considered to be rare and tend to be ignored. On the contrary, pseudogenes are surprisingly prevalent and can persist for long evolutionary time periods, representing a record of once-functional genetic characteristics. Most importantly, pseudogenes provide an insight into prokaryotic evolutionary history as a record of phenotypic traits that have been lost. Focusing on the intracellular and symbiotic bacteria in which pseudogenes predominate, this review discusses the importance of identifying pseudogenes to fully understand the abilities of bacteria, and to understand prokaryotes within their evolutionary context.

Solving gap metabolites and blocked reactions in genome-scale models: application to the metabolic network of Blattabacterium cuenoti

Background

Metabolic reconstruction is the computational-based process that aims to elucidate the network of metabolites interconnected through reactions catalyzed by activities assigned to one or more genes. Reconstructed models may contain inconsistencies that appear as gap metabolites and blocked reactions. Although automatic methods for solving this problem have been previously developed, there are many situations where manual curation is still needed.

Results

We introduce a general definition of gap metabolite that allows its detection in a straightforward manner. Moreover, a method for the detection of Unconnected Modules, defined as isolated sets of blocked reactions connected through gap metabolites, is proposed. The method has been successfully applied to the curation of iCG238, the genome-scale metabolic model for the bacterium Blattabacterium cuenoti, obligate endosymbiont of cockroaches.

Conclusion

We found the proposed approach to be a valuable tool for the curation of genome-scale metabolic models. The outcome of its application to the genome-scale model B. cuenoti iCG238 is a more accurate model version named as B. cuenoti iMP240.

A systems biology approach to studying the role of microbes in human health

Host-microbe interactions play a crucial role in human health and disease. Of the various systems biology approaches, reconstruction of genome-scale metabolic networks combined with constraint-based modeling has been particularly successful at in silico predicting the phenotypic characteristics of single organisms. Here, we summarize recent studies, which have applied this approach to investigate microbe-microbe and host-microbe metabolic interactions. This approach can be also expanded to investigate the properties of an entire microbial community, as well as single organisms within the community. We illustrate that the constraint-based modeling approach is suitable to model host-microbe interactions at molecular resolution and will enable systematic investigation of metabolic links between the human host and its microbes. Such host-microbe models, combined with experimental data, will ultimately further our understanding of how microbes influence human health.

Contribution of phage-derived genomic islands to the virulence of facultative bacterial pathogens

Facultative pathogens have extremely dynamic pan-genomes, to a large extent derived from bacteriophages and other mobile elements. We developed a simple approach to identify phage-derived genomic islands and apply it to show that pathogens from diverse bacterial genera are significantly enriched in clustered phage-derived genes compared with related benign strains. These findings show that genome expansion by integration of prophages containing virulence factors is a major route of evolution of facultative bacterial pathogens.

Out of fuzzy chemistry: from prebiotic chemistry to metabolic networks

The origin of life on Earth was a chemical affair. So how did primitive biochemical systems originate from geochemical and cosmochemical processes on the young planet? Contemporary research into the origins of life subscribes to the Darwinian principle of material causes operating in an evolutionary context, as advocated by A. I. Oparin and J. B. S. Haldane in the 1920s. In its simplest form (e.g., a bacterial cell) extant biological complexity relies on the functional integration of metabolic networks and replicative genomes inside a lipid boundary. Different research programmes have explored the prebiotic plausibility of each of these autocatalytic subsystems and combinations thereof: self-maintained networks of small molecules, template chemistry, and self-reproductive vesicles. This tutorial review focuses on the debates surrounding the origin of metabolism and offers a brief overview of current studies on the evolution of metabolic networks. I suggest that a leitmotif in the origin and evolution of metabolism is the role played by catalysers' substrate ambiguity and multifunctionality.