Metabolic stasis in an ancient symbiosis: genome-scale metabolic networks from two Blattabacterium cuenoti strains, primary endosymbionts of cockroaches

Comparative Transcriptomics of Fat Bodies between Symbiotic and Quasi-Aposymbiotic Adult Females of Blattella germanica with Emphasis on the Metabolic Integration with Its Endosymbiont Blattabacterium and Its Immune System

We explored the metabolic integration of Blattella germanica and its obligate endosymbiont Blattabacterium cuenoti by the transcriptomic analysis of the fat body of quasi-aposymbiotic cockroaches, where the endosymbionts were almost entirely removed with rifampicin. Fat bodies from quasi-aposymbiotic insects displayed large differences in gene expression compared to controls. In quasi-aposymbionts, the metabolism of phenylalanine and tyrosine involved in cuticle sclerotization and pigmentation increased drastically to compensate for the deficiency in the biosynthesis of these amino acids by the endosymbionts. On the other hand, the uricolytic pathway and the biosynthesis of uric acid were severely decreased, probably because the reduced population of endosymbionts was unable to metabolize urea to ammonia. Metabolite transporters that could be involved in the endosymbiosis process were identified. Immune system and antimicrobial peptide (AMP) gene expression was also reduced in quasi-aposymbionts, genes encoding peptidoglycan-recognition proteins, which may provide clues for the maintenance of the symbiotic relationship, as well as three AMP genes whose involvement in the symbiotic relationship will require additional analysis. Finally, a search for AMP-like factors that could be involved in controlling the endosymbiont identified two orphan genes encoding proteins smaller than 200 amino acids underexpressed in quasi-aposymbionts, suggesting a role in the host-endosymbiont relationship.

Genome-Scale Model of Rhizopus microsporus: Metabolic integration of a fungal holobiont with its bacterial and viral endosymbionts

Rhizopus microsporus often lives in association with bacterial and viral symbionts that alter its biology. This fungal model represents an example of the complex interactions established among diverse organisms in functional holobionts. We constructed a Genome-Scale Model (GSM) of the fungal-bacterial-viral holobiont (iHol). We employed a constraint-based method to calculate the metabolic fluxes to decipher the metabolic interactions of the symbionts with their host. Our computational analyses of iHol simulate the holobiont's growth and the production of the toxin rhizoxin. Analyses of the calculated fluxes between R. microsporus in symbiotic (iHol) versus asymbiotic conditions suggest that changes in the lipid and nucleotide metabolism of the host are necessary for the functionality of the holobiont. Glycerol plays a pivotal role in the fungal-bacterial metabolic interaction, as its production does not compromise fungal growth, and Mycetohabitans bacteria can efficiently consume it. Narnavirus RmNV-20S and RmNV-23S affected the nucleotide metabolism without impacting the fungal-bacterial symbiosis. Our analyses highlighted the metabolic stability of Mycetohabitans throughout its co-evolution with the fungal host. We also predicted changes in reactions of the bacterial metabolism required for the active production of rhizoxin. This iHol is the first GSM of a fungal holobiont.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Exploring Gut Microbial Dynamics and Symbiotic Interaction in Blattella germanica Using Rifampicin

Blattella germanica harbours two cohabiting symbiotic systems: an obligate endosymbiont, Blattabacterium, located inside bacteriocytes and vertically transmitted, which is key in nitrogen metabolism, and abundant and complex gut microbiota acquired horizontally (mainly by coprophagy) that must play an important role in host physiology. In this work, we use rifampicin treatment to deepen the knowledge on the relationship between the host and the two systems. First, we analysed changes in microbiota composition in response to the presence and removal of the antibiotic with and without faeces in one generation. We found that, independently of faeces supply, rifampicin-sensitive bacteria are strongly affected at four days of treatment, and most taxa recover after treatment, although some did not reach control levels. Second, we tried to generate an aposymbiotic population, but individuals that reached the second generation were severely affected and no third generation was possible. Finally, we established a mixed population with quasi-aposymbiotic and control nymphs sharing an environment in a blind experiment. The analysis of the two symbiotic systems in each individual after reaching the adult stage revealed that endosymbiont's load does not affect the composition of the hindgut microbiota, suggesting that there is no interaction between the two symbiotic systems in Blattella germanica.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

Identification of the Gene Repertoire of the IMD Pathway and Expression of Antimicrobial Peptide Genes in Several Tissues and Hemolymph of the Cockroach Blattella germanica

Antimicrobial peptide (AMP) genes, triggered by Toll and IMD pathways, are essential components of the innate immune system in the German cockroach Blattella germanica. Besides their role in killing pathogenic bacteria, AMPs could be involved in controlling its symbiotic systems (endosymbiont and microbiota). We found that the IMD pathway was active in the adult female transcriptomes of six tissues (salivary glands, foregut, midgut, hindgut, Malpighian tubules and fat body) and hemolymph. Total expression of AMP genes was high in hemolymph and salivary glands and much lower in the other sample types. The expression of specific AMP genes was very heterogeneous among sample types. Two genes, defensin_g10 and drosomycin_g5, displayed relevant expression in the seven sample types, although higher in hemolymph. Other genes only displayed high expression in one tissue. Almost no expression of attacin-like and blattellicin genes was observed in any sample type, although some of them were among the genes with the highest expression in adult female whole bodies. The expression of AMP genes in salivary glands could help control pathogens ingested with food and even determine gut microbiota composition. The low expression levels in midgut and hindgut are probably related to the presence of beneficial microbiota. Furthermore, a reduction in the expression of AMP genes in fat body could be the way to prevent damage to the population of the endosymbiont Blattabacterium cuenoti within bacteriocytes.

Of Cockroaches and Symbionts: Recent Advances in the Characterization of the Relationship between Blattella germanica and Its Dual Symbiotic System

Mutualistic stable symbioses are widespread in all groups of eukaryotes, especially in insects, where symbionts have played an essential role in their evolution. Many insects live in obligate relationship with different ecto- and endosymbiotic bacteria, which are needed to maintain their hosts’ fitness in their natural environment, to the point of even relying on them for survival. The case of cockroaches (Blattodea) is paradigmatic, as both symbiotic systems coexist in the same organism in two separated compartments: an intracellular endosymbiont (Blattabacterium) inside bacteriocytes located in the fat body, and a rich and complex microbiota in the hindgut. The German cockroach Blattella germanica is a good model for the study of symbiotic interactions, as it can be maintained in the laboratory in controlled populations, allowing the perturbations of the two symbiotic systems in order to study the communication and integration of the tripartite organization of the host–endosymbiont–microbiota, and to evaluate the role of symbiotic antimicrobial peptides (AMPs) in host control over their symbionts. The importance of cockroaches as reservoirs and transmission vectors of antibiotic resistance sequences, and their putative interest to search for AMPs to deal with the problem, is also discussed.

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.

Computing irreversible minimal cut sets in genome-scale metabolic networks via flux cone projection

MOTIVATION

Minimal cut sets (MCSs) for metabolic networks are sets of reactions which, if they are removed from the network, prevent a target reaction from carrying flux. To compute MCSs different methods exist, which may fail to find sufficiently many MCSs for larger genome-scale networks.

RESULTS

Here we introduce irreversible minimal cut sets (iMCSs). These are MCSs that consist of irreversible reactions only. The advantage of iMCSs is that they can be computed by projecting the flux cone of the metabolic network on the set of irreversible reactions, which usually leads to a smaller cone. Using oriented matroid theory, we show how the projected cone can be computed efficiently and how this can be applied to find iMCSs even in large genome-scale networks.

AVAILABILITY AND IMPLEMENTATION

Software is freely available at https://sourceforge.net/projects/irreversibleminimalcutsets/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Finding MEMo: minimum sets of elementary flux modes

Metabolic network reconstructions are widely used in computational systems biology for in silico studies of cellular metabolism. A common approach to analyse these models are elementary flux modes (EFMs), which correspond to minimal functional units in the network. Already for medium-sized networks, it is often impossible to compute the set of all EFMs, due to their huge number. From a practical point of view, this might also not be necessary because a subset of EFMs may already be sufficient to answer relevant biological questions. In this article, we study MEMos or minimum sets of EFMs that can generate all possible steady-state behaviours of a metabolic network. The number of EFMs in a MEMo may be by several orders of magnitude smaller than the total number of EFMs. Using MEMos, we can compute generating sets of EFMs in metabolic networks where the whole set of EFMs is too large to be enumerated.

Flux tope analysis: studying the coordination of reaction directions in metabolic networks

Motivation

Elementary flux mode (EFM) analysis allows an unbiased description of metabolic networks in terms of minimal pathways (involving a minimal set of reactions). To date, the enumeration of EFMs is impracticable in genome-scale metabolic models. In a complementary approach, we introduce the concept of a flux tope (FT), involving a maximal set of reactions (with fixed directions), which allows one to study the coordination of reaction directions in metabolic networks and opens a new way for EFM enumeration.

Results

A FT is a (nontrivial) subset of the flux cone specified by fixing the directions of all reversible reactions. In a consistent metabolic network (without unused reactions), every FT contains a 'maximal pathway', carrying flux in all reactions. This decomposition of the flux cone into FTs allows the enumeration of EFMs (of individual FTs) without increasing the problem dimension by reaction splitting. To develop a mathematical framework for FT analysis, we build on the concepts of sign vectors and hyperplane arrangements. Thereby, we observe that FT analysis can be applied also to flux optimization problems involving additional (inhomogeneous) linear constraints. For the enumeration of FTs, we adapt the reverse search algorithm and provide an efficient implementation. We demonstrate that (biomass-optimal) FTs can be enumerated in genome-scale metabolic models of B.cuenoti and E.coli, and we use FTs to enumerate EFMs in models of M.genitalium and B.cuenoti.

Availability and implementation

The source code is freely available at https://github.com/mpgerstl/FTA.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Genome analysis of new Blattabacterium spp., obligatory endosymbionts of Periplaneta fuliginosa and P. japonica

The successful adaptation of cockroaches is, in part, dependent of the activity of their obligatory endosymbionts, Blattabacterium spp., which are involved in uric acid degradation, nitrogen assimilation and nutrient provisioning. Their strategic localization, within bacteriocytes in the proximities of uric acid storage cells (urocytes), highlights their importance in the recycling of nitrogen from urea and ammonia, end-products not secreted by their host insects. In this study, we present the complete genome sequence of two new Blattabacterium spp. from Periplaneta fuliginosa (BPfu) and P. japonica (BPja), and detailed comparison with other Blattabacterium strains from different cockroach species. The genomes of BPfu and BPja show a high degree of stability as showed with for other Blattabacterium representatives, only presenting a 19-kb fragment inversion between BPja and BPfu. In fact, the phylogenomics showed BPja as an ancestor species of BPfu, BPLAN (P. americana) and BBor (Blatta orientalis), in congruence with their host cockroach phylogeny. Their functional profile is similar and closest to the omnivorous strain BBge (Blattella germanica). Interesting, BPja possesses the complete set of enzymes involved sulfate assimilatory pathway only found in BBge and BMda (Mastotermes darwiniensis). The newly sequenced genomes of BPja and BPfu emphasise the remarkable stability of Blattabacterium genomes supported by their long-term coevolution and obligatory lifestyle in their host insect.

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.

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

An open question in evolutionary biology is how does the selection–drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host–symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shift.

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.

Exact quantification of cellular robustness in genome-scale metabolic networks

Motivation: Robustness, the ability of biological networks to uphold their functionality in spite of perturbations, is a key characteristic of all living systems. Although several theoretical approaches have been developed to formalize robustness, it still eludes an exact quantification. Here, we present a rigorous and quantitative approach for the structural robustness of metabolic networks by measuring their ability to tolerate random reaction (or gene) knockouts.

Results: In analogy to reliability theory, based on an explicit consideration of all possible knockout sets, we exactly quantify the probability of failure for a given network function (e.g. growth). This measure can be computed if the network’s minimal cut sets (MSCs) are known. We show that even in genome-scale metabolic networks the probability of (network) failure can be reliably estimated from MSCs with lowest cardinalities. We demonstrate the applicability of our theory by analyzing the structural robustness of multiple Enterobacteriaceae and Blattibacteriaceae and show a dramatically low structural robustness for the latter. We find that structural robustness develops from the ability to proliferate in multiple growth environments consistent with experimentally found knowledge.

Conclusion: The probability of (network) failure provides thus a reliable and easily computable measure of structural robustness and redundancy in (genome-scale) metabolic networks.

Availability and implementation: Source code is available under the GNU General Public License at https://github.com/mpgerstl/networkRobustnessToolbox.

Contact:juergen.zanghellini@boku.ac.at

Supplementary information:Supplementary data are available at Bioinformatics online.

Physiology of environmental adaptations and resource acquisition in cockroaches

Cockroaches are a group of insects that evolved early in geological time. Because of their antiquity, they for the most part display generalized behavior and physiology and accordingly have frequently been used as model insects to examine physiological and biochemical mechanisms involved with water balance, nutrition, reproduction, genetics, and insecticide resistance. As a result, a considerable amount of information on these topics is available. However, there is much more to be learned by employing new protocols, microchemical analytical techniques, and molecular biology tools to explore many unanswered questions.

Small but powerful, the primary endosymbiont of moss bugs, Candidatus Evansia muelleri, holds a reduced genome with large biosynthetic capabilities

Moss bugs (Coleorrhyncha: Peloridiidae) are members of the order Hemiptera, and like many hemipterans, they have symbiotic associations with intracellular bacteria to fulfill nutritional requirements resulting from their unbalanced diet. The primary endosymbiont of the moss bugs, Candidatus Evansia muelleri, is phylogenetically related to Candidatus Carsonella ruddii and Candidatus Portiera aleyrodidarum, primary endosymbionts of psyllids and whiteflies, respectively. In this work, we report the genome of Candidatus Evansia muelleri Xc1 from Xenophyes cascus, which is the only obligate endosymbiont present in the association. This endosymbiont possesses an extremely reduced genome similar to Carsonella and Portiera. It has crossed the borderline to be considered as an autonomous cell, requiring the support of the insect host for some housekeeping cell functions. Interestingly, in spite of its small genome size, Evansia maintains enriched amino acid (complete or partial pathways for ten essential and six nonessential amino acids) and sulfur metabolisms, probably related to the poor diet of the insect, based on bryophytes, which contains very low levels of nitrogenous and sulfur compounds. Several facts, including the congruence of host (moss bugs, whiteflies, and psyllids) and endosymbiont phylogenies and the retention of the same ribosomal RNA operon during genome reduction in Evansia, Portiera, and Carsonella, suggest the existence of an ancient endosymbiotic Halomonadaceae clade associated with Hemiptera. Three possible scenarios for the origin of these three primary endosymbiont genera are proposed and discussed.

The cockroach Blattella germanica obtains nitrogen from uric acid through a metabolic pathway shared with its bacterial endosymbiont

Uric acid stored in the fat body of cockroaches is a nitrogen reservoir mobilized in times of scarcity. The discovery of urease in Blattabacterium cuenoti, the primary endosymbiont of cockroaches, suggests that the endosymbiont may participate in cockroach nitrogen economy. However, bacterial urease may only be one piece in the entire nitrogen recycling process from insect uric acid. Thus, in addition to the uricolytic pathway to urea, there must be glutamine synthetase assimilating the released ammonia by the urease reaction to enable the stored nitrogen to be metabolically usable. None of the Blattabacterium genomes sequenced to date possess genes encoding for those enzymes. To test the host's contribution to the process, we have sequenced and analysed Blattella germanica transcriptomes from the fat body. We identified transcripts corresponding to all genes necessary for the synthesis of uric acid and its catabolism to urea, as well as for the synthesis of glutamine, asparagine, proline and glycine, i.e. the amino acids required by the endosymbiont. We also explored the changes in gene expression with different dietary protein levels. It appears that the ability to use uric acid as a nitrogen reservoir emerged in cockroaches after its age-old symbiotic association with bacteria.

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.

Comparative genomics of Blattabacterium cuenoti: the frozen legacy of an ancient endosymbiont genome

Many insect species have established long-term symbiotic relationships with intracellular bacteria. Symbiosis with bacteria has provided insects with novel ecological capabilities, which have allowed them colonize previously unexplored niches. Despite its importance to the understanding of the emergence of biological complexity, the evolution of symbiotic relationships remains hitherto a mystery in evolutionary biology. In this study, we contribute to the investigation of the evolutionary leaps enabled by mutualistic symbioses by sequencing the genome of Blattabacterium cuenoti, primary endosymbiont of the omnivorous cockroach Blatta orientalis, and one of the most ancient symbiotic associations. We perform comparative analyses between the Blattabacterium cuenoti genome and that of previously sequenced endosymbionts, namely those from the omnivorous hosts the Blattella germanica (Blattelidae) and Periplaneta americana (Blattidae), and the endosymbionts harbored by two wood-feeding hosts, the subsocial cockroach Cryptocercus punctulatus (Cryptocercidae) and the termite Mastotermes darwiniensis (Termitidae). Our study shows a remarkable evolutionary stasis of this symbiotic system throughout the evolutionary history of cockroaches and the deepest branching termite M. darwiniensis, in terms of not only chromosome architecture but also gene content, as revealed by the striking conservation of the Blattabacterium core genome. Importantly, the architecture of central metabolic network inferred from the endosymbiont genomes was established very early in Blattabacterium evolutionary history and could be an outcome of the essential role played by this endosymbiont in the host's nitrogen economy.

The nutrient supplying capabilities of Uzinura, an endosymbiont of armoured scale insects

An emerging common physiological feature of plant sap-feeding insects is the presence of bacterial endosymbionts capable of providing essential nutrients to their host. These microbial partners are inviable outside of specialized host tissues, and therefore a cultivation-independent approach, namely high-throughput next-generation genome sequencing, can be used to characterize their gene content and metabolic potential. To this end, we sequenced the first complete genome of the obligate endosymbiont, Candidatus 'Uzinura diaspidicola', of armoured scale insects. At 263 431 bp, Uzinura has an extremely reduced genome that is composed largely of genes encoding enzymes involved in translation and amino acid biosynthesis. The tiny size of the Uzinura genome parallels that observed in some other insect endosymbionts. Despite this extreme genome reduction, the absence of a known obligate partner bacterial symbiont suggests that Uzinura alone can supply sufficient nutrients to its host.

Making the most of "omics" for symbiosis research

Omics, including genomics, proteomics, and metabolomics, enable us to explain symbioses in terms of the underlying molecules and their interactions. The central task is to transform molecular catalogs of genes, metabolites, etc., into a dynamic understanding of symbiosis function. We review four exemplars of omics studies that achieve this goal, through defined biological questions relating to metabolic integration and regulation of animal-microbial symbioses, the genetic autonomy of bacterial symbionts, and symbiotic protection of animal hosts from pathogens. As omic datasets become increasingly complex, computationally sophisticated downstream analyses are essential to reveal interactions not evident from visual inspection of the data. We discuss two approaches, phylogenomics and transcriptional clustering, that can divide the primary output of omics studies-long lists of factors-into manageable subsets, and we describe how they have been applied to analyze large datasets and generate testable hypotheses.