Volume 12 Supplement 1
Metabolic stasis in an ancient symbiosis: genome-scale metabolic networks from two Blattabacterium cuenoti strains, primary endosymbionts of cockroaches
© González-Domenech et al; licensee BioMed Central Ltd. 2012
Published: 18 January 2012
Cockroaches are terrestrial insects that strikingly eliminate waste nitrogen as ammonia instead of uric acid. Blattabacterium cuenoti (Mercier 1906) strains Bge and Pam are the obligate primary endosymbionts of the cockroaches Blattella germanica and Periplaneta americana, respectively. The genomes of both bacterial endosymbionts have recently been sequenced, making possible a genome-scale constraint-based reconstruction of their metabolic networks. The mathematical expression of a metabolic network and the subsequent quantitative studies of phenotypic features by Flux Balance Analysis (FBA) represent an efficient functional approach to these uncultivable bacteria.
We report the metabolic models of Blattabacterium strains Bge (iCG238) and Pam (iCG230), comprising 296 and 289 biochemical reactions, associated with 238 and 230 genes, and 364 and 358 metabolites, respectively. Both models reflect both the striking similarities and the singularities of these microorganisms. FBA was used to analyze the properties, potential and limits of the models, assuming some environmental constraints such as aerobic conditions and the net production of ammonia from these bacterial systems, as has been experimentally observed. In addition, in silico simulations with the iCG238 model have enabled a set of carbon and nitrogen sources to be defined, which would also support a viable phenotype in terms of biomass production in the strain Pam, which lacks the first three steps of the tricarboxylic acid cycle. FBA reveals a metabolic condition that renders these enzymatic steps dispensable, thus offering a possible evolutionary explanation for their elimination. We also confirm, by computational simulations, the fragility of the metabolic networks and their host dependence.
The minimized Blattabacterium metabolic networks are surprisingly similar in strains Bge and Pam, after 140 million years of evolution of these endosymbionts in separate cockroach lineages. FBA performed on the reconstructed networks from the two bacteria helps to refine the functional analysis of the genomes enabling us to postulate how slightly different host metabolic contexts drove their parallel evolution.
Recently, the genomes of two different strains of Blattabacterium cuenoti (Mercier 1906), Bge and Pam, obligate primary endosymbionts of the cockroaches Blattella germanica and Periplaneta americana, respectively, have been sequenced [1, 2]. Blattabacterium constitutes a clade within the class Flavobacteria, the phylum Bacteroidetes, which contains several instances of symbionts of insects, e.g., “Candidatus Sulcia muelleri”, obligate endosymbiont of cicadas, spittlebugs and leafhoppers , “Candidatus Cardinium”, symbiont of the white fly Bemisia tabaci , and “Candidatus Vestibaculum illigatum”, which establishes a symbiosis with the gut flagellate Staurojoenina sp. associated to the termite Neotermes cubanus . All these endosymbiont bacteria are relatively distant from free-living members within the phylum Bacteroidetes . Thus, if we assume that the age of a symbiotic association of a primary endosymbiont corresponds to the oldest fossil record of its host, we estimate the time of divergence between B. cuenoti and its free-living cousins to be 250 Myr , thus being possibly one of the most ancient mutualistic insect symbioses described so far.
Cockroaches, natural hosts of Blattabacterium sp., excrete waste nitrogen as ammonia [8–11] unlike most terrestrial insects, which eliminate it as uric acid . Ammonotely, the ancestral character present in aquatic animals, has traditionally been considered maladaptive in terrestrial animals . The enigmatic return of cockroaches to ammonotely seems to be related to the role of bacterial endosymbiosis in their nitrogen economy. López-Sánchez et al.  showed the presence of urease activity in endosymbiont-enriched extracts of the cockroaches B. germanica and P. americana. Stoichiometric analysis of the core of the reconstructed metabolic networks would suggest that these endosymbiotic bacteria participate in the nitrogen metabolism of the host. Physiological studies ([1, 8] and references therein) suggest that uric acid may represent a form of nitrogen storage in cockroaches and that B. cuenoti may produce ammonia from uric-derived metabolites provided by the host. In fact, the cockroach fat body contains specialized cells storing uric acid (urocytes) that are in close proximity to the cells containing endosymbionts (bacteriocytes) .
A common feature of genomes from bacterial endosymbionts is their strict conservation of gene order and remarkable differential gene losses in the different lineages [14–16]. In the case of the Bge and Pam strains, comparative genomics reveals both a high degree of conservation in their chromosomal architecture and in the gene repertoires (accounting for a total of 627 and 619 genes in Bge and Pam, respectively) despite the low sequence similarity observed (~85% nucleotide sequence identity) . Thus, the metabolic networks of these endosymbionts should be similar, differing only slightly. These differences might be analyzed from a qualitative point of view by comparison between the inferred metabolic maps, but this approach does not allow quantitative evaluation of how these inequalities might affect the functional capabilities of each microorganism. Constraint-based models of metabolic networks represent an efficient framework for a quantitative understanding of microbial physiology . In fact, computational simulations with constraint-based models are approaches that help to predict cellular phenotypes given particular environmental conditions, with a high correspondence between experimental results and predictions [18–20]. It is worth mentioning that they are especially suitable for reconstructed networks from uncultivable microorganism, as it is the case of primary endosymbionts. Thus, Flux Balance Analysis (FBA) is one of these useful techniques for the study of obligate intracellular bacteria, since it reconstructs fluxes through a network requiring neither kinetic parameters nor other detailed information on enzymes . This modeling method is based on the stoichiometric coefficients of each reaction and the assumption of the system at steady-state . FBA calculates metabolites fluxes through the metabolic reactions that optimize an objective function –usually biomass production–, i.e., how much each reaction contributes to the phenotype desired.
In this study, we have reconstructed the metabolic networks of Bge and Pam strains of B. cuenoti, focusing on their metabolic abilities and relating them to the symbiotic interaction with their host. Following the protocol proposed by Thiele and Palsson , we have quantitatively predicted their biochemical potential by FBA, assuming biomass formation as objective function. In addition, in some simulations we have imposed the constraint of ammonia release from both endosymbionts, in coherence with the physiological observations  and as expected by the measured urease activity and the stoichiometric analysis performed by López-Sánchez et al. . We have performed sensitivity and robustness analyses and deduced how these endosymbionts may be related to their cockroach hosts metabolically. We offer an overview of the remarkably stable metabolic relationships in these old symbioses as well as providing an explanation for a possible environmental cause of the loss of genes coding for enzymes in a central pathway, such as the TCA cycle in one of the endosymbionts.
Metabolic models and FBA simulations
Characteristics of metabolic reconstructions from the strains Bge and Pam of B. cuenoti.
Reactions with protein-encoding gene model assignments (GPRs)
The excretion of ammonia from the system, a phenomenon compatible with the physiological and experimental observations (for review see  and  and references therein), was always observed in simulations with both models under minimal conditions. The efflux of ammonia reached maximum levels when L-Glu uptake was simulated by the system. However, the efflux of ammonia stopped and could even be reversed when 2-oxoglutarate or succinate were provided to both metabolic networks. This was due to an increased assimilation of ammonia through displacement of the glutamate dehydrogenase reaction (EC 188.8.131.52) in the assimilative direction.
Sensitivity and robustness analysis
Finally, we also explored the robustness of both metabolic networks by randomly removing genes. The starting point was the aerobic and minimal conditions already described and then we also evaluated the fragility of the network when additional carbon sources were provided to the system. As cutoff value for gene essentiality a >99% decrease in the biomass production after the gene deletion was used, as described by Thomas et al. . For the Bge strain network, a set of essential genes was determined ranging between 76.1 % (minimal medium) and 72.3 % (with added glycerol) of the total genes comprised in the model. With the Pam network we found a genetic essentiality between 79.6 % (minimal medium) and 73.5 % (with added fumarate, L-malate or glycerol).
Uncultivable bacteria can be studied by in silicosimulations
In this paper we describe the genome-scale metabolic networks corresponding to two strains of B. cuenoti, Bge and Pam, the endosymbiotic bacteria of the cockroaches B. germanica and P. americana, respectively. Despite the approximately 140-Myr of parallel evolution, both metabolic networks showed striking conservation and we decided to compare their functionality by means of a stoichiometric approach such as FBA. This computational methodology has already been successfully used in a study of the metabolic network robustness of B. aphidicola, the primary endosymbiont of aphids, in comparison to E. coli  and for the simulation of reductive evolution in endosymbionts [25, 26]. Thus, FBA represents a valid strategy for the functional study of those bacterial species that pose important obstacles to their empirical study, as it is the case of the uncultivable endosymbionts. In this work we used the E. coli model as a reference since to the best of our knowledge there are no empirical data on the biomass function of any members within the phylum Bacteroidetes. In the absence of information related to real biomass composition of the modeled organism, the use of the equations of E. coli is considered a reliable approach and an acceptable starting point [19, 27–29].
The simulations allowed us to identify the minimal environmental components for a functional metabolic network (Fig. 2). For instance, both Blattabacterium networks show a strict dependence on L-Gln supply from the host due to the absence of glutamine synthase in both endosymbionts. This dependence of the functionality on the availability of some chemical species may also suggest a possible regulatory role of the external medium in the metabolic behavior of the bacterium. In other biological systems, like the nitrogen-fixing nodules of Leguminosae, oxygen availability modulation by the host has been suggested as a mechanism of punishment to cheaters in the symbiotic relationship . Our in silico simulations (Fig. 5) suggest that access to L-Gln and/or oxygen is a good candidate for a control mechanism of cockroaches over their endosymbiotic population.
Why Blattabacteriumstrain Pam has lost the first three enzymatic steps of the TCA cycle
One of the few differences between the metabolic networks inferred for the Bge and Pam strains of Blattabacterium is that the latter has an incomplete TCA cycle (compare  and ). Sabree and coworkers  hypothesized that the two missing enzymes, aconitase (EC 184.108.40.206, acnA) and isocitrate dehydrogenase (EC 220.127.116.11, icd), in the Pam strain metabolic network, can be functionally substituted by the enzymes 3-isopropylmalate isomerase (EC 18.104.22.168, leuC) and 3-isopropylmalate dehydrogenase (EC 22.214.171.124, leuB), respectively. However, the first enzymatic step of the TCA cycle (citrate synthase, EC 126.96.36.199, gltA) is also absent and apparently there is no other alternative solution to this absent activity. Although the functional substitution of two out of three missing metabolic steps in the TCA cycle cannot be excluded, here we have shown the dispensability of all three genes to obtain a functional phenotype in terms of biomass production under certain conditions. Thus, the proposal of functional substitutions by homologous enzymes is an unnecessary conjecture in this case. There are two reasons: (i) as shown in Figure 1, the lack of the three afore-mentioned steps does not generate true dead-end metabolites, and (ii) there is an alternative way to keep a fully functional metabolic network without the first three enzymes in the TCA cycle.
Our simulations show that the Pam network behaves like the Bge network if an anaplerotic reaction (i.e. the uptake of L-Glu or 2-oxoglutarate) is provided. Under these circumstances, the metabolic fluxes are redirected around the TCA cycle (Fig. 4) and, as shown in Figure 6, the sensitivity analysis demonstrates that the flux through the first three enzymatic steps of the TCA cycle can be null. This behavior may explain the dispensability of the corresponding gltA, acnA, and icd genes if the host provides the endosymbiont with any of the above-mentioned compounds. In other words, the provision of a non-essential amino acid to the endosymbiont by the host may offer a set of biochemical conditions favoring the loss of central metabolic genes in one particular evolutionary lineage.
The loss of these three enzymatic steps in the Pam strain of Blattabacterium is an example of how the essentiality of genes may change when the environmental conditions change. Studies of flux connectivity (i.e. reactions that always work together)  and synthetic lethality analysis (i.e. searching the effect of multiple gene deletions)  have shown that in free-living bacteria, such as E. coli or Helicobacter pylori, the enzymes coded by the gltA, acnA and icd genes form a subset of essential steps. This enzymatic subset was also determined during our analysis of elementary flux modes in Blattabacterium Bge . Thus, it is conceivable that during the transition to intracellular lifestyle, the ancestor of Blattabacterium strain Pam found a set of chemical conditions in the host cell making those three formerly essential genes dispensable and thus allowing their loss en bloc.
Blattabacteriummetabolic networks are fragile
The percentage of essential genes predicted for both Blattabacterium strains in gene deletion simulations (ranging from 72.3% to 79.6%) is lower to that determined for B. aphidicola, primary endosymbiont of aphids, which showed a fraction of 84% of essential genes in a similar simulation . Those values of genetic essentiality in endosymbiotic metabolic networks are far from the robustness observed in models of free-living bacteria, e.g., around 15% of essential genes coding for metabolic enzymes in E. coli . Thus, endosymbiotic metabolic networks are less redundant than networks from free-living bacteria. In comparison to the extreme fragility of a minimalist metabolic network, theoretically deduced from comparative genomics  and analyzed by Gabaldón et al. , with 98% of essential genes, endosymbiont metabolic networks show an intermediate degree of robustness, and may represent different stages of the reductive evolutionary process associated to intracellular lifestyle.
Blattabacteriumhas a key role in the nitrogen economy of cockroaches
Our working hypothesis is that Blattabacterium played a key role during the transition from uricotely to a use of urates as nitrogen storage in cockroaches. The elementary flux mode analysis and the enzymatic assays performed by López-Sánchez et al.  indicated that the central metabolism of Blattabacterium can use urea (and some other nitrogen compounds, as non-essential amino acids) and excrete ammonia. As shown in this work, under minimal conditions the reconstructed metabolic networks of the Bge and Pam strains produce ammonia when biomass growth is optimized. This metabolic performance is compatible with the classical physiological observations made by Cochran and coworkers . In addition, physiological studies with cockroaches indicate that uric acid is a form of nitrogen storage instead of a major waste product like in most insects . According to our hypothesis, the fat body metabolism would produce urea from uric acid and the endosymbiont urease would transform urea into ammonia to be used again, partially by the endosymbiont (i.e. synthesis of Glu via the displacement of the Glu dehydrogenase reaction) and partially by the host, especially for glutamine biosynthesis by Gln synthase. It is remarkable that this enzymatic reaction is absent in Blattabacterium, although the metabolic networks of both Bge and Pam strains contain 9 Gln-consuming reactions (in addition to the requirement of Gln for protein synthesis represented by the corresponding tRNA for Gln and a gene coding for glutamine tRNA ligase, glnS).
In that context, the retention of a urease in Blattabacterium makes evolutionary sense as a key piece of the metabolic mosaic of the cockroach nitrogen economy, whereas the bacterial dependence on a Gln supply by the host contributes to the obligate character of this symbiotic association. The dependence on host-supplied Gln has also been recently described in Blochmannia vafer, the primary endosymbiont of the ant Camponotus vafer, which contains urease but lacks both Glu dehydrogenase and Gln synthase .
The released ammonia observed by physiologists would correspond to the escape of some ammonia produced by the system when all the ammonia-utilizing reactions are saturated, a side effect of the serial transformation from uric acid to urea to ammonia to glutamate/glutamine. In this metabolic framework, our in silico modeling was performed with the constraint of ammonia release by the endosymbiont. The mathematical expression of the metabolic networks, thus, helps us understand the systemic properties of the host-endosymbiont relationships. Practically speaking, it serves for the better design of an experimental strategy to functionally characterize the pathway from uric acid to glutamine in cockroaches.
One of our aims was to perform a genome-scale constraint-based modeling of the metabolisms of two different strains of B. cuenoti, Bge and Pam, primary endosymbionts of the cockroaches B. germanica and P. americana, respectively, which are the result of a parallel evolution during the last 140 million years. A striking feature of the two bacteria is not only the genome architecture conservation, as observed in other similar systems, but also the few gene losses undergone in the different lineages. Thus, both metabolic networks differ from each other in only seven enzymatic reactions. The FBA approach has allowed us to evaluate the different host influences that might explain the loss or retention of certain genes, which is not easy to elucidate a priori by visual inspection of the respective metabolic maps. In addition, the fragility shown by the metabolic networks is compatible with a constancy of environmental conditions, and it is the expected outcome for minimal metabolisms derived from the streamlining of endosymbiotic bacterial genomes. The model predictions will allow us to address future functional analyses, and formulate new hypotheses on the metabolic interdependence in the ancient symbiosis between B. cuenoti and cockroaches.
Definition of the iCG238 and iCG230 models and FBA simulations
We reconstructed the iCG238 and iCG230 networks using the E. coli K-12 iJR904 model as a starting point . From this model, we proceeded as Thomas et al.  removing all reactions associated with pseudogenes, genes without homologs in those strains or unconnected with the biomass reaction (e.g., gltX, dna, encoding genes of tRNA ligases and DNA glycosylases). We employed the OrthoMCL algorithm  to search for orthologs between E. coli K-12 and the different strains of Blattabacterium sp. as well as between the two Blattabacterium strains in order to obtain a first draft of the metabolic models (inflation thresholds, between 1.2 and 5, choosing in each case the best, normally 1.5 and 3). In addition, the urease reaction, supported by experimental evidence , was also added to both models, and ubiquinone was replaced by menaquinone in each metabolic reconstruction since this is the unique or majority quinone in the Flavobacteriaceae family . To further curate the models, we performed additional BLAST searches  among the corresponding strain of Blattabacterium, other flavobacteria and E. coli K-12 available in GenBank (e-values below e-11), to incorporate reactions either absent in E. coli or undetected due to the divergence among strains. In addition, we identified functional domains by means of the interface SMART (Simple Modular Architecture Research Tool) (http://smart.emblheidelberg.de/help/smart_about.shtml) [41, 42]. Flux balance analysis (FBA) was performed using the COBRA toolbox , a freely available Matlab toolbox and the models were described using the Systems Biology Markup Language (SBML)  (Additional Files 5 and 6).
We used the biomass equation derived from the iJR904 E. coli model  with a few adaptations derived on updated network of such microorganism, i.e. iAF1260 . In particular we added the cofactors thiamine diphosphate and tetrahydrofolate. Additionally, we adjusted the amounts of the four different deoxynucleotide triphosphates in the biomass equation to reflect the GC content of the Blattabacterium strains (Bge, 27 mol%; Pam, 28 mol%). Furthermore, since Blattabacterium strains are unable to completely synthesize cardiolipin, glycogen, lipopolysaccharide, and spermidine, we removed these components from the biomass equation.
The study of network robustness was performed with the function robustnessAnalysis of the COBRA toolbox . In addition, we evaluated the effect of a gene deletion experiment on cellular growth of the resultant mutant using the option singleGeneDeletion of the COBRA toolbox. We set to zero the upper and lower flux bounds for the reaction(s) corresponding to the simulated deleted gene. If a single gene is associated with multiple reactions, the deletion of that gene will result in the removal of all associated reactions. On the contrary, a reaction that can be catalyzed by multiple non-interacting gene products will not be removed in a single gene deletion. The possible results of a single deletion are unchanged maximal growth (non-lethal), reduced maximal growth or no growth (lethal). We simulated growth and subsequent fragility analysis with all the different sources which enhance/support biomass formation.
CMGD: postdoctoral specialist in Microbiology and Systems Biology; EB: postdoctoral specialist in Bioinformatics, Evolutionary Genomics and Systems Biology; RPN: PhD student specialist in Genetics, ‘omics’ Sciences and Bioinformatics; AM: Full Professor of Genetics; JP: Associate Professor of Biochemistry and Molecular Biology; AL: Full Professor of Genetics.
List of abbreviations used
Blattabacterium cuenoti strain of cockroach Blattella germanica
Flux Balance Analysis
Gene to protein to reaction association
Blattabacterium cuenoti strain of cockroach Periplaneta Americana
- TCA cycle:
Tricarboxylic acid cycle.
Financial support was provided by grants BFU2009-12895-C02-01/BMC (Ministerio de Ciencia e Innovación, Spain) to AL and Prometeo Program (Generalitat Valenciana) to AM. Dr. González-Domenech was supported by grant from the University of Granada. R. Patiño-Navarrete was recipient of a fellowship from Ministerio de Educación y Ciencia, Spain. We also thank to Mr. Alejandro Manzano for his assistance with bioinformatic issues, Dr. Alex Neef for helpful discussions as well as two anonymous reviewers for their valuable comments. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
This article has been published as part of BMC Microbiology Volume 11 Supplement 1, 2012: Arthropod symbioses: from fundamental studies to pest and disease mangement. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2180/12?issue=S1.
- López-Sánchez MJ, Neef A, Peretó J, Patiño-Navarrete R, Pignatelli M, Latorre A, Moya A: Evolutionary convergence and nitrogen metabolism in Blattabacterium strain Bge, primary endosymbiont of the cockroach Blattella germanica. PLoS Genet. 2009, 5: e1000721-10.1371/journal.pgen.1000721.PubMed CentralView ArticlePubMedGoogle Scholar
- Sabree ZL, Kambhampati S, Moran NA: Nitrogen recycling and nutritional provisioning by Blattabacterium, the cockroach endosymbiont. Proc Natl Acad Sci USA. 2009, 106: 19521-1956. 10.1073/pnas.0907504106.PubMed CentralView ArticlePubMedGoogle Scholar
- Gruwell ME, Morse GE, Normark BB: Phylogenetic congruence of armored scale insects (Hemiptera: Diaspididae) and their primary endosymbionts from the phylum Bacteroidetes. Mol Phylogenet Evol. 2007, 44: 267-280. 10.1016/j.ympev.2007.01.014.View ArticlePubMedGoogle Scholar
- Gottlieb Y, Ghanim M, Gueguen G, Kontsedalov S, Vavre F, Fleury F, Zchori-Fein E: Inherited intracellular ecosystem: symbiotic bacteria share bacteriocytes in whiteflies. FASEB J. 2008, 22: 2591-2599. 10.1096/fj.07-101162.View ArticlePubMedGoogle Scholar
- Stingl U, Maass A, Radek R, Brune A: Symbionts of the gut flagellate Staurojoenina sp. from Neotermes cubanus represent a novel, termite-associated lineage of Bacteroidales: description of 'Candidatus Vestibaculum illigatum'. Microbiology. 2004, 150: 2229-2235. 10.1099/mic.0.27135-0.View ArticlePubMedGoogle Scholar
- Sabree ZL, Degnan PH, Moran NA: Chromosome stability and gene loss in cockroach endosymbionts. Appl Environ Microbiol. 2010, 76: 4076-4079. 10.1128/AEM.00291-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Grimaldi D, Engel MS: Evolution of Insects. Edited by: Grimaldi D, Engel MS. 2005, New York/Cambridge: Cambridge University PressGoogle Scholar
- Cochran DG: Nitrogen excretion in cockroaches. Annu Rev Entomol. 1985, 30: 29-49. 10.1146/annurev.en.30.010185.000333.View ArticleGoogle Scholar
- Mullins DE, Cochran DG: Nitrogen excretion in cockroaches: uric acid is not a major product. Science. 1972, 177: 699-701. 10.1126/science.177.4050.699.View ArticlePubMedGoogle Scholar
- Mullins DE, Cochran DG: Nitrogen metabolism in the American cockroach: an examination of whole body and fat body regulation of cations in response to nitrogen balance. J Exp Biol. 1974, 61: 557-570.PubMedGoogle Scholar
- O'Donnell M: Insect excretory mechanisms. Advances in Insect Physiology. Edited by: Simpson SJ. 2008, New York: Academic Press, 35: 1-122.Google Scholar
- Needham J: Contributions of chemical physiology to the problem of reversibility in evolution. Biol Rev. 1938, 13: 225-251. 10.1111/j.1469-185X.1938.tb00515.x.View ArticleGoogle Scholar
- Cochran DG, Mullins DE, Mullins KJ: Cytological changes in the fat body of the American cockroach, Periplaneta americana, in relation to dietary nitrogen levels. Ann Entomol Soci Amer. 1979, 72: 197-205.View ArticleGoogle Scholar
- Moya A, Peretó J, Gil R, Latorre A: Learning how to live together: genomic insights into prokaryote-animal symbioses. Nat Rev Genet. 2008, 9: 218-229. 10.1038/nrg2319.View ArticlePubMedGoogle Scholar
- Moran NA, McCutcheon JP, Nakabachi A: Genomics and evolution of heritable bacterial symbionts. Annu Rev Genet. 2008, 42: 165-190. 10.1146/annurev.genet.41.110306.130119.View ArticlePubMedGoogle Scholar
- Lamelas A, Gosalbes MJ, Moya A, Latorre A: New clues about the evolutionary history of metabolic losses in bacterial endosymbionts, provided by the genome of Buchnera aphidicola from the aphid Cinara tujafilina. Appl Environ Microbiol. 2011, 77: 4446-4454. 10.1128/AEM.00141-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards JS, Covert M, Palsson B: Metabolic modelling of microbes: the flux-balance approach. Environ Microbiol. 2002, 4: 133-140. 10.1046/j.1462-2920.2002.00282.x.View ArticlePubMedGoogle Scholar
- Covert MW, Palsson BO: Transcriptional regulation in constraints-based metabolic models of Escherichia coli. J Biol Chem. 2002, 277: 28058-28064. 10.1074/jbc.M201691200.View ArticlePubMedGoogle Scholar
- Puchalka J, Oberhardt MA, Godinho M, Bielecka A, Regenhardt D, Timmis KN, Papin JA, Martins dos Santos V: Genome-scale reconstruction and analysis of the Pseudomonas putida KT2440 metabolic network facilitates applications in biotechnology. PLoS Comput Biol. 2008, 4: e1000210-10.1371/journal.pcbi.1000210.PubMed CentralView ArticlePubMedGoogle Scholar
- Schilling CH, Covert MW, Famili I, Church GM, Edwards JS, Palsson BO: Genome-scale metabolic model of Helicobacter pylori 26695. J Bacteriol. 2002, 184: 4582-4593. 10.1128/JB.184.16.4582-4593.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Orth JD, Thiele I, Palsson BØ: What is flux balance analysis?. Nat Biotechnol. 2010, 28: 245-248. 10.1038/nbt.1614.PubMed CentralView ArticlePubMedGoogle Scholar
- Thiele I, Palsson BØ: A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010, 5: 93-121.PubMed CentralView ArticlePubMedGoogle Scholar
- Locke M: The fat body. Microscopic anatomy of invertebrates. Insecta Mundi. Edited by: Harrison FW, Locke M. 1998, New York: Wiley, 11B: 641-686.Google Scholar
- Thomas GH, Zucker J, Macdonald SJ, Sorokin A, Goryanin I, Douglas AE: A fragile metabolic network adapted for cooperation in the symbiotic bacterium Buchnera aphidicola. BMC Syst Biol. 2009, 3: 24-10.1186/1752-0509-3-24.PubMed CentralView ArticlePubMedGoogle Scholar
- Pál C, Papp B, Lercher MJ, Csermely P, Oliver SG, Hurst LD: Chance and necessity in the evolution of minimal metabolic networks. Nature. 2006, 440: 667-670. 10.1038/nature04568.View ArticlePubMedGoogle Scholar
- Yizhak K, Tuller T, Papp B, Ruppin E: Metabolic modeling of endosymbiont genome reduction on a temporal scale. Mol Syst Biol. 2011, 7: 479-PubMed CentralView ArticlePubMedGoogle Scholar
- Ates O, Toksoy Oner E, Arga KY: Genome-scale reconstruction of metabolic network for a halophilic extremophile, Chromohalobacter salexigens DSM 3043. BMC Syst Biol. 2011, 5: 12-10.1186/1752-0509-5-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Oberhardt MA, Puchalka J, Fryer KE, Martins dos Santos VA, Papin JA: Genome-scale metabolic network analysis of the opportunistic pathogen Pseudomonas aeruginosa PAO1. J Bacteriol. 2008, 190: 2790-2803. 10.1128/JB.01583-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Thiele I, Weekes D, Li Z, Jaroszewski L, Ginalski K, Deacon AM, Wooley J, Lesley SA, Wilson IA, Palsson B, Osterman A, Godzik A: Three-dimensional structural view of the central metabolic network of Thermotoga maritima. Science. 2009, 325: 1544-1549. 10.1126/science.1174671.PubMed CentralView ArticlePubMedGoogle Scholar
- Kiers ET, Rousseau RA, West SA, Denison RF: Host sanctions and the legume-rhizobium mutualism. Nature. 2003, 425: 78-81. 10.1038/nature01931.View ArticlePubMedGoogle Scholar
- Burgard AP, Nikolaev EV, Schilling CH, Maranas CD: Flux coupling analysis of genome-scale metabolic network reconstructions. Genome Res. 2004, 14: 301-312. 10.1101/gr.1926504.PubMed CentralView ArticlePubMedGoogle Scholar
- Suthers PF, Zomorrodi A, Maranas CD: Genome-scale gene/reaction essentiality and synthetic lethality analysis. Mol Syst Biol. 2009, 5: 301-PubMed CentralView ArticlePubMedGoogle Scholar
- Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, Broadbelt LJ, Hatzimanikatis V, Palsson BO: A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol. 2007, 3: 121-PubMed CentralView ArticlePubMedGoogle Scholar
- Gil R, Silva FJ, Peretó J, Moya A: Determination of the core of a minimal bacterial gene set. Microbiol Mol Biol Rev. 2004, 68: 518-537. 10.1128/MMBR.68.3.518-537.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Gabaldon T, Peretó J, Montero F, Gil R, Latorre A, Moya A: Structural analyses of a hypothetical minimal metabolism. Philos Trans R Soc Lond B Biol Sci. 2007, 362: 1751-1762. 10.1098/rstb.2007.2067.PubMed CentralView ArticlePubMedGoogle Scholar
- Williams LE, Wernegreen JJ: Unprecedented loss of ammonia assimilation capability in a urease-encoding bacterial mutualist. BMC Genomics. 2010, 11: 687-PubMed CentralView ArticlePubMedGoogle Scholar
- Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol. 2003, 4: R54-10.1186/gb-2003-4-9-r54.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen F, Mackey AJ, Stoeckert CJ, Roos DS: OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res. 2006, 34: 363-368. 10.1093/nar/gkj123.View ArticleGoogle Scholar
- Bernardet JF, Nakagawa Y: An introduction to the family Flavobacteriaceae. The Prokaryotes: a Handbook on the Biology of Bacteria. Edited by: Dworkin M et al. 2006, New York: Springer-Verlag, 7: 455-480. 3View ArticleGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: Identification of signaling domains. Proc Natl Acad Sci USA. 1998, 95: 5857-5864. 10.1073/pnas.95.11.5857.PubMed CentralView ArticlePubMedGoogle Scholar
- Letunic I, Goodstadt L, Dickens NJ, Doerks T, Schultz J, Mott R, Ciccarelli F, Copley RR, Ponting CP, Bork P: Recent improvements to the SMART domain based sequence annotation resource. Nucleic Acids Res. 2002, 30: 242-244. 10.1093/nar/30.1.242.PubMed CentralView ArticlePubMedGoogle Scholar
- Becker SA, Feist AM, Mo ML, Hannum G, Palsson BØ, Herrgard MJ: Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc. 2007, 2: 727-738. 10.1038/nprot.2007.99.View ArticlePubMedGoogle Scholar
- Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novère N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J, SBML Forum: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics. 2003, 19: 524-531. 10.1093/bioinformatics/btg015.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.