Bacterial endosymbionts of insects play a central role in upgrading the diet of their hosts. In certain cases, such as aphids and tsetse flies, endosymbionts complement the metabolic capacity of hosts living on nutrient-deficient diets, while the bacteria harbored by omnivorous carpenter ants are involved in nitrogen recycling. In this study, we describe the genome sequence and inferred metabolism of Blattabacterium strain Bge, the primary Flavobacteria endosymbiont of the omnivorous German cockroach Blattella germanica. Through comparative genomics with other insect endosymbionts and free-living Flavobacteria we reveal that Blattabacterium strain Bge shares the same distribution of functional gene categories only with Blochmannia strains, the primary Gamma-Proteobacteria endosymbiont of carpenter ants. This is a remarkable example of evolutionary convergence during the symbiotic process, involving very distant phylogenetic bacterial taxa within hosts feeding on similar diets. Despite this similarity, different nitrogen economy strategies have emerged in each case. Both bacterial endosymbionts code for urease but display different metabolic functions: Blochmannia strains produce ammonia from dietary urea and then use it as a source of nitrogen, whereas Blattabacterium strain Bge codes for the complete urea cycle that, in combination with urease, produces ammonia as an end product. Not only does the cockroach endosymbiont play an essential role in nutrient supply to the host, but also in the catabolic use of amino acids and nitrogen excretion, as strongly suggested by the stoichiometric analysis of the inferred metabolic network. Here, we explain the metabolic reasons underlying the enigmatic return of cockroaches to the ancestral ammonotelic state.

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.

Supplementary figure legends.Supplementary Fig. 1 Number of antibiotic resistance genes (ARG) according to the phylogenetic group among E. coli ST410 isolates. The horizontal lines in the boxes represent the median number of ARG associated to the three clades: Basal (n=22); FQR clade (n=92) and OXA-181 subclade (n=41). The box boundaries represent the first and third quartiles of the distribution and box-plot whiskers span 1.5 times the interquartile range of the distribution. Outliers are denoted as black dots outside whiskers. Statistical significances were tested with a one-sided Wilcoxon rank-sum test. ***, P < 0.001. Supplementary Fig. 2. Recombination events in isolates of the E. coli ST410 FQR clade. Phylogeny of the E. coli ST410 FQR clade (left). On the right side, the recombined blocks identified by Gubbins 1 are highlighted in red when they occurred in an internal branch and affect more than one isolate and in blue when they occurred in a terminal branch affecting a single isolate. The core genome of E. coli ST410 FQR strains is represented above the figure and annotated genes correspond to the limits of major recombination events. Supplementary Fig. 3. Phylogeny and mutations in non-redundant CP-Ec isolates of ST167, ST101 and ST156. ML phylogenies were estimated as for Fig. 1, using 5,843, 21,810, and 18,746 non-recombinant SNPs for a. Ec ST167, rooted with the ST10 strain MG1655 (NC_00913) b. Ec ST101, and c. Ec ST156 respectively, both rooted with the ST1128 strain IAI1 (NC_011741). 14 distantly related non-CP-Ec ST156 isolates have been removed from the phylogenetic analysis. Branch tips indicates the presence and the type of carbapenemase according to the figure key on the left. On the right side of the tree are represented, from left to right: blaCTX-M ESBL, gyrA and parC QRDR mutations, mutations in the ftsI gene, and genetic events affecting ompF and ompC according to the figure key at the bottom. SNPs in the dcw region are represented by small vertical red bars. Genes from the dcw locus are indicated by arrows and the ftsI gene in red. Black arrowheads indicate isolates used as reference for SNPs mapping in the dcw gene cluster.

BackgroundCockroaches 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.ResultsWe 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.ConclusionsThe 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.

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.

No abstract

Blattabacteria are intracellular endosymbionts of cockroaches and primitive termites that belong to the class Flavobacteria and live only in specialized cells in the abdominal fat body of their hosts. In the present study we determined genome sizes as well as genome copy numbers for the endosymbionts of three cockroach species, Blattella germanica, Periplaneta americana and Blatta orientalis. The sole presence of blattabacteria in the fat body was demonstrated by rRNA-targeting techniques. The genome sizes of the three blattabacteria were determined by pulsed field gel electrophoresis. The resulting total genome sizes for the three symbionts were all approximately 650 +/- 15 kb. Comparison of the genome sizes with those of free-living Bacteroidetes shows extended reduction, as occurs in other obligatory insect endosymbionts. Genome copy numbers were determined based on cell counts and determination of DNA amounts via quantitative PCR. Values between 10.2 and 18.3 and between 323 and 353 were found for the symbionts of P. americana and B. orientalis respectively. Polyploidy in intracellular bacteria may play a significant role in the genome reduction process.

The health and environmental risks associated with antibiotic use in aquaculture have promoted bacterial probiotics as an alternative approach to control fish infections in vulnerable larval and juvenile stages. However, evidence-based identification of probiotics is often hindered by the complexity of bacteria-host interactions and host variability in microbiologically uncontrolled conditions. While these difficulties can be partially resolved using gnotobiotic models harboring no or reduced microbiota, most host-microbe interaction studies are carried out in animal models with little relevance for fish farming. Here we studied host-microbiota-pathogen interactions in a germ-free and gnotobiotic model of rainbow trout (Oncorhynchus mykiss), one of the most widely cultured salmonids. We demonstrated that germ-free larvae raised in sterile conditions displayed no significant difference in growth after 35 days compared to conventionally-raised larvae, but were extremely sensitive to infection by Flavobacterium columnare, a common freshwater fish pathogen causing major economic losses worldwide. Furthermore, re-conventionalization with 11 culturable species from the conventional trout microbiota conferred resistance to F. columnare infection. Using mono-re-conventionalized germ-free trout, we identified that this protection is determined by a commensal Flavobacterium strain displaying antibacterial activity against F. columnare. Finally, we demonstrated that use of gnotobiotic trout is a suitable approach for the identification of both endogenous and exogenous probiotic bacterial strains protecting teleostean hosts against F. columnare. This study therefore establishes an ecologically-relevant gnotobiotic model for the study of host-pathogen interactions and colonization resistance in farmed fish.