Genomic encyclopedia of sugar utilization pathways in the Shewanella genus

Rapid quantification of mutant fitness in diverse bacteria by sequencing randomly bar-coded transposons

Transposon mutagenesis with next-generation sequencing (TnSeq) is a powerful approach to annotate gene function in bacteria, but existing protocols for TnSeq require laborious preparation of every sample before sequencing. Thus, the existing protocols are not amenable to the throughput necessary to identify phenotypes and functions for the majority of genes in diverse bacteria. Here, we present a method, random bar code transposon-site sequencing (RB-TnSeq), which increases the throughput of mutant fitness profiling by incorporating random DNA bar codes into Tn5 and mariner transposons and by using bar code sequencing (BarSeq) to assay mutant fitness. RB-TnSeq can be used with any transposon, and TnSeq is performed once per organism instead of once per sample. Each BarSeq assay requires only a simple PCR, and 48 to 96 samples can be sequenced on one lane of an Illumina HiSeq system. We demonstrate the reproducibility and biological significance of RB-TnSeq with Escherichia coli, Phaeobacter inhibens, Pseudomonas stutzeri, Shewanella amazonensis, and Shewanella oneidensis. To demonstrate the increased throughput of RB-TnSeq, we performed 387 successful genome-wide mutant fitness assays representing 130 different bacterium-carbon source combinations and identified 5,196 genes with significant phenotypes across the five bacteria. In P. inhibens, we used our mutant fitness data to identify genes important for the utilization of diverse carbon substrates, including a putative d-mannose isomerase that is required for mannitol catabolism. RB-TnSeq will enable the cost-effective functional annotation of diverse bacteria using mutant fitness profiling.

A large challenge in microbiology is the functional assessment of the millions of uncharacterized genes identified by genome sequencing. Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach to assign phenotypes and functions to genes. However, the current strategies for TnSeq are too laborious to be applied to hundreds of experimental conditions across multiple bacteria. Here, we describe an approach, random bar code transposon-site sequencing (RB-TnSeq), which greatly simplifies the measurement of gene fitness by using bar code sequencing (BarSeq) to monitor the abundance of mutants. We performed 387 genome-wide fitness assays across five bacteria and identified phenotypes for over 5,000 genes. RB-TnSeq can be applied to diverse bacteria and is a powerful tool to annotate uncharacterized genes using phenotype data.

A new regulatory mechanism for bacterial lipoic acid synthesis

Lipoic acid, an essential enzyme cofactor, is required in three domains of life. In the past 60 years since its discovery, most of the pathway for lipoic acid synthesis and metabolism has been elucidated. However, genetic control of lipoic acid synthesis remains unclear. Here, we report integrative evidence that bacterial cAMP-dependent signaling is linked to lipoic acid synthesis in Shewanella species, the certain of unique marine-borne bacteria with special ability of metal reduction. Physiological requirement of protein lipoylation in γ-proteobacteria including Shewanella oneidensis was detected using Western blotting with rabbit anti-lipoyl protein primary antibody. The two genes (lipB and lipA) encoding lipoic acid synthesis pathway were proved to be organized into an operon lipBA in Shewanella, and the promoter was mapped. Electrophoretic mobility shift assays confirmed that the putative CRP-recognizable site (AAGTGTGATCTATCTTACATTT) binds to cAMP-CRP protein with origins of both Escherichia coli and Shewanella. The native lipBA promoter of Shewanella was fused to a LacZ reporter gene to create a chromosome lipBA-lacZ transcriptional fusion in E. coli and S. oneidensis, allowing us to directly assay its expression level by β-galactosidase activity. As anticipated, the removal of E. coli crp gene gave above fourfold increment of lipBA promoter-driven β-gal expression. The similar scenario was confirmed by both the real-time quantitative PCR and the LacZ transcriptional fusion in the crp mutant of Shewanella. Furthermore, the glucose effect on the lipBA expression of Shewanella was evaluated in the alternative microorganism E. coli. As anticipated, an addition of glucose into media effectively induces the transcriptional level of Shewanella lipBA in that the lowered cAMP level relieves the repression of lipBA by cAMP-CRP complex. Therefore, our finding might represent a first paradigm mechanism for genetic control of bacterial lipoic acid synthesis.

Using comparative genomics to drive new discoveries in microbiology

Bioinformatics looks to many microbiologists like a service industry. In this view, annotation starts with what is known from experiments in the lab, makes reasonable inferences of which genes match other genes in function, builds databases to make all that we know accessible, but creates nothing truly new. Experiments lead, then biocuration and computational biology follow. But the astounding success of genome sequencing is changing the annotation paradigm. Every genome sequenced is an intercepted coded message from the microbial world, and as all cryptographers know, it is easier to decode a thousand messages than a single message. Some biology is best discovered not by phenomenology, but by decoding genome content, forming hypotheses, and doing the first few rounds of validation computationally. Through such reasoning, a role and function may be assigned to a protein with no sequence similarity to any protein yet studied. Experimentation can follow after the discovery to cement and to extend the findings. Unfortunately, this approach remains so unfamiliar to most bench scientists that lab work and comparative genomics typically segregate to different teams working on unconnected projects. This review will discuss several themes in comparative genomics as a discovery method, including highly derived data, use of patterns of design to reason by analogy, and in silico testing of computationally generated hypotheses.

Improving mediated electron transport in anodic bioelectrocatalysis

A microbial fuel cell loaded with bio-cocatalyst beads immobilized with recombinant riboflavin-secreting Escherichia coli shows significantly enhanced performance.

Chlorate reduction in Shewanella algae ACDC is a recently acquired metabolism characterized by gene loss, suboptimal regulation and oxidative stress

Previous work on respiratory chlorate reduction has biochemically identified the terminal reductase ClrABC and the chlorite detoxifying enzyme Cld. In Shewanella algae ACDC, genes encoding these enzymes reside on composite transposons whose core we refer to as the chlorate reduction composite transposon interior (CRI). To better understand this metabolism in ACDC, we used RNA-seq and proteomics to predict carbon and electron flow during chlorate reduction and posit that formate is an important electron carrier with lactate as the electron donor, but that NADH predominates on acetate. Chlorate-specific transcription of electron transport chain components or the CRI was not observed, but clr and cld transcription was attenuated by oxygen. The major chlorate-specific response related to oxidative stress and was indicative of reactive chlorine species production. A genetic system based on rpsL-streptomycin counter selection was developed to further dissect the metabolism, but ACDC readily lost the CRI via homologous recombination of the composite transposon's flanking insertion sequences. An engineered strain containing a single chromosomal CRI did not grow on chlorate, but overexpression of cld and its neighbouring cytochrome c restored growth. We postulate that the recently acquired CRI underwent copy-number expansion to circumvent insufficient expression of key genes in the pathway.

Inference of interactions in cyanobacterial-heterotrophic co-cultures via transcriptome sequencing

We used deep sequencing technology to identify transcriptional adaptation of the euryhaline unicellular cyanobacterium Synechococcus sp. PCC 7002 and the marine facultative aerobe Shewanella putrefaciens W3-18-1 to growth in a co-culture and infer the effect of carbon flux distributions on photoautotroph-heterotroph interactions. The overall transcriptome response of both organisms to co-cultivation was shaped by their respective physiologies and growth constraints. Carbon limitation resulted in the expansion of metabolic capacities, which was manifested through the transcriptional upregulation of transport and catabolic pathways. Although growth coupling occurred via lactate oxidation or secretion of photosynthetically fixed carbon, there was evidence of specific metabolic interactions between the two organisms. These hypothesized interactions were inferred from the excretion of specific amino acids (for example, alanine and methionine) by the cyanobacterium, which correlated with the downregulation of the corresponding biosynthetic machinery in Shewanella W3-18-1. In addition, the broad and consistent decrease of mRNA levels for many Fe-regulated Synechococcus 7002 genes during co-cultivation may indicate increased Fe availability as well as more facile and energy-efficient mechanisms for Fe acquisition by the cyanobacterium. Furthermore, evidence pointed at potentially novel interactions between oxygenic photoautotrophs and heterotrophs related to the oxidative stress response as transcriptional patterns suggested that Synechococcus 7002 rather than Shewanella W3-18-1 provided scavenging functions for reactive oxygen species under co-culture conditions. This study provides an initial insight into the complexity of photoautotrophic-heterotrophic interactions and brings new perspectives of their role in the robustness and stability of the association.

Comparisons of Shewanella strains based on genome annotations, modeling, and experiments

BACKGROUND

Shewanella is a genus of facultatively anaerobic, Gram-negative bacteria that have highly adaptable metabolism which allows them to thrive in diverse environments. This quality makes them an attractive bacterial target for research in bioremediation and microbial fuel cell applications. Constraint-based modeling is a useful tool for helping researchers gain insights into the metabolic capabilities of these bacteria. However, Shewanella oneidensis MR-1 is the only strain with a genome-scale metabolic model constructed out of 21 sequenced Shewanella strains.

RESULTS

In this work, we updated the model for Shewanella oneidensis MR-1 and constructed metabolic models for three other strains, namely Shewanella sp. MR-4, Shewanella sp. W3-18-1, and Shewanella denitrificans OS217 which span the genus based on the number of genes lost in comparison to MR-1. We also constructed a Shewanella core model that contains the genes shared by all 21 sequenced strains and a few non-conserved genes associated with essential reactions. Model comparisons between the five constructed models were done at two levels - for wildtype strains under different growth conditions and for knockout mutants under the same growth condition. In the first level, growth/no-growth phenotypes were predicted by the models on various carbon sources and electron acceptors. Cluster analysis of these results revealed that the MR-1 model is most similar to the W3-18-1 model, followed by the MR-4 and OS217 models when considering predicted growth phenotypes. However, a cluster analysis done based on metabolic gene content revealed that the MR-4 and W3-18-1 models are the most similar, with the MR-1 and OS217 models being more distinct from these latter two strains. As a second level of comparison, we identified differences in reaction and gene content which give rise to different functional predictions of single and double gene knockout mutants using Comparison of Networks by Gene Alignment (CONGA). Here, we showed how CONGA can be used to find biomass, metabolic, and genetic differences between models.

CONCLUSIONS

We developed four strain-specific models and a general core model that can be used to do various in silico studies of Shewanella metabolism. The developed models provide a platform for a systematic investigation of Shewanella metabolism to aid researchers using Shewanella in various biotechnology applications.

Metabolically engineered glucose-utilizing Shewanella strains under anaerobic conditions

Comparative genome analysis of Shewanella strains predicted that the strains metabolize preferably two- and three-carbon carbohydrates as carbon/electron source because many Shewanella genomes are deficient of the key enzymes in glycolysis (e.g., glucokinase). In addition, all Shewanella genomes are known to have only one set of genes associated with the phosphotransferase system required to uptake sugars. To engineer Shewanella strains that can utilize five- and six-carbon carbohydrates, we constructed glucose-utilizing Shewanella oneidensis MR-1 by introducing the glucose facilitator (glf; ZMO0366) and glucokinase (glk; ZMO0369) genes of Zymomonas mobilis. The engineered MR-1 strain was able to grow on glucose as a sole carbon/electron source under anaerobic conditions. The glucose affinity (Ks) and glucokinase activity in the engineered MR-1 strain were 299.46 mM and 0.259 ± 0.034 U/g proteins. The engineered strain was successfully applied to a microbial fuel cell system and exhibited current generation using glucose as the electron source.

Comparative genomics and functional analysis of rhamnose catabolic pathways and regulons in bacteria

L-rhamnose (L-Rha) is a deoxy-hexose sugar commonly found in nature. L-Rha catabolic pathways were previously characterized in various bacteria including Escherichia coli. Nevertheless, homology searches failed to recognize all the genes for the complete L-Rha utilization pathways in diverse microbial species involved in biomass decomposition. Moreover, the regulatory mechanisms of L-Rha catabolism have remained unclear in most species. A comparative genomics approach was used to reconstruct the L-Rha catabolic pathways and transcriptional regulons in the phyla Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, and Thermotogae. The reconstructed pathways include multiple novel enzymes and transporters involved in the utilization of L-Rha and L-Rha-containing polymers. Large-scale regulon inference using bioinformatics revealed remarkable variations in transcriptional regulators for L-Rha utilization genes among bacteria. A novel bifunctional enzyme, L-rhamnulose-phosphate aldolase (RhaE) fused to L-lactaldehyde dehydrogenase (RhaW), which is not homologous to previously characterized L-Rha catabolic enzymes, was identified in diverse bacteria including Chloroflexi, Bacilli, and Alphaproteobacteria. By using in vitro biochemical assays we validated both enzymatic activities of the purified recombinant RhaEW proteins from Chloroflexus aurantiacus and Bacillus subtilis. Another novel enzyme of the L-Rha catabolism, L-lactaldehyde reductase (RhaZ), was identified in Gammaproteobacteria and experimentally validated by in vitro enzymatic assays using the recombinant protein from Salmonella typhimurium. C. aurantiacus induced transcription of the predicted L-Rha utilization genes when L-Rha was present in the growth medium and consumed L-Rha from the medium. This study provided comprehensive insights to L-Rha catabolism and its regulation in diverse Bacteria.

Transcriptional regulation of the carbohydrate utilization network in Thermotoga maritima

Hyperthermophilic bacteria from the Thermotogales lineage can produce hydrogen by fermenting a wide range of carbohydrates. Previous experimental studies identified a large fraction of genes committed to carbohydrate degradation and utilization in the model bacterium Thermotoga maritima. Knowledge of these genes enabled comprehensive reconstruction of biochemical pathways comprising the carbohydrate utilization network. However, transcriptional factors (TFs) and regulatory mechanisms driving this network remained largely unknown. Here, we used an integrated approach based on comparative analysis of genomic and transcriptomic data for the reconstruction of the carbohydrate utilization regulatory networks in 11 Thermotogales genomes. We identified DNA-binding motifs and regulons for 19 orthologous TFs in the Thermotogales. The inferred regulatory network in T. maritima contains 181 genes encoding TFs, sugar catabolic enzymes and ABC-family transporters. In contrast to many previously described bacteria, a transcriptional regulation strategy of Thermotoga does not employ global regulatory factors. The reconstructed regulatory network in T. maritima was validated by gene expression profiling on a panel of mono- and disaccharides and by in vitro DNA-binding assays. The observed upregulation of genes involved in catabolism of pectin, trehalose, cellobiose, arabinose, rhamnose, xylose, glucose, galactose, and ribose showed a strong correlation with the UxaR, TreR, BglR, CelR, AraR, RhaR, XylR, GluR, GalR, and RbsR regulons. Ultimately, this study elucidated the transcriptional regulatory network and mechanisms controlling expression of carbohydrate utilization genes in T. maritima. In addition to improving the functional annotations of associated transporters and catabolic enzymes, this research provides novel insights into the evolution of regulatory networks in Thermotogales.

Cloning and characterization of a novel cold-active glycoside hydrolase family 1 enzyme with β-glucosidase, β-fucosidase and β-galactosidase activities

Background

Cold-active enzymes, sourced from cold-adapted organisms, are characterized by high catalytic efficiencies at low temperatures compared with their mesophilic counterparts, which have poor activity. This property makes them advantageous for biotechnology applications as it: (i) saves energy costs, (ii) shortens the times for processes operated at low temperatures, (iii) protects thermosensitive substrates or products of the enzymatic reaction, (iv) prevents undesired chemical transformations, and (v) prevents the loss of volatile compounds.

Results

A bglMKg gene that encodes a monomeric cold-active glycoside hydrolase family 1 enzyme with an apparent molecular mass of 50 kDa was isolated by the functional screening of a marine metagenomic library. The BglMKg enzyme was expressed in E. coli, purified by FPLC and characterized. The recombinant BglMKg could effectively hydrolyze various chromogenic substrates and β-linked oligosaccharides, and had remarkably high β-galactosidase, β-glucosidase and β-fucosidase activities. Because of the lack of information about the usefulness of β-fucosidases in industry, further characterization of the enzymatic properties of BglMKg was only carried out with substrates specific for β-glucosidase or β-galactosidase. The BglMKg had maximal β-galactosidase and β-glucosidase activities at approximately 40°C and 45°C, respectively. The optimum pH for β-galactosidase activity was 6.5, whereas the optimum pH for β-glucosidase activity was 7.5. In general, the enzyme was stable below 30°C and from pHs 6.0 to 8.0. The results of the kinetic studies revealed that BglMKg more efficiently hydrolyzed β-glucosidase substrates than β-galactosidase ones.

Conclusions

BglMKg is a small, monomeric, cold-active β-glucosidase with additional enzymatic activities. It was efficiently expressed in E. coli indicating that BglMKg might be a candidate for industrial applications.

Functional diversification of ROK-family transcriptional regulators of sugar catabolism in the Thermotogae phylum

Large and functionally heterogeneous families of transcription factors have complex evolutionary histories. What shapes specificities toward effectors and DNA sites in paralogous regulators is a fundamental question in biology. Bacteria from the deep-branching lineage Thermotogae possess multiple paralogs of the repressor, open reading frame, kinase (ROK) family regulators that are characterized by carbohydrate-sensing domains shared with sugar kinases. We applied an integrated genomic approach to study functions and specificities of regulators from this family. A comparative analysis of 11 Thermotogae genomes revealed novel mechanisms of transcriptional regulation of the sugar utilization networks, DNA-binding motifs and specific functions. Reconstructed regulons for seven groups of ROK regulators were validated by DNA-binding assays using purified recombinant proteins from the model bacterium Thermotoga maritima. All tested regulators demonstrated specific binding to their predicted cognate DNA sites, and this binding was inhibited by specific effectors, mono- or disaccharides from their respective sugar catabolic pathways. By comparing ligand-binding domains of regulators with structurally characterized kinases from the ROK family, we elucidated signature amino acid residues determining sugar-ligand regulator specificity. Observed correlations between signature residues and the sugar-ligand specificities provide the framework for structure functional classification of the entire ROK family.

Tagaturonate-fructuronate epimerase UxaE, a novel enzyme in the hexuronate catabolic network in Thermotoga maritima

Thermotoga maritima is a marine hyperthermophilic microorganism that degrades a wide range of simple and complex carbohydrates including pectin and produces fermentative hydrogen at high yield. Galacturonate and glucuronate, two abundant hexuronic acids in pectin and xylan, respectively, are catabolized via committed metabolic pathways to supply carbon and energy for a variety of microorganisms. By a combination of bioinformatics and experimental techniques we identified a novel enzyme family (named UxaE) catalysing a previously unknown reaction in the hexuronic acid catabolic pathway, epimerization of tagaturonate to fructuronate. The enzymatic activity of the purified recombinant tagaturonate epimerase from T. maritima was directly confirmed and kinetically characterized. Its function was also confirmed by genetic complementation of the growth of the Escherichia coli uxaB knockout mutant strain on galacturonate. An inferred novel galacturonate to mannonate catabolic pathway in T. maritima was reconstituted in vitro using a mixture of recombinant purified enzymes UxaE, UxaC and UxuB. Members of the newly identified UxaE family were identified in ~50 phylogenetically diverse heterotrophic bacteria from aquatic and soil environments. The genomic context of respective genes and reconstruction of associated pathways suggest that UxaE enzymatic and biological function remains conserved in all of these species.

Microbial metabolomics

Microbial metabolomics constitutes an integrated component of systems biology. By studying the complete set of metabolites within a microorganism and monitoring the global outcome of interactions between its development processes and the environment, metabolomics can potentially provide a more accurate snap shot of the actual physiological state of the cell. Recent advancement of technologies and post-genomic developments enable the study and analysis of metabolome. This unique contribution resulted in many scientific disciplines incorporating metabolomics as one of their “omics” platforms. This review focuses on metabolomics in microorganisms and utilizes selected topics to illustrate its impact on the understanding of systems microbiology.

Diversity and versatility of the Thermotoga maritima sugar kinome

Sugar phosphorylation is an indispensable committed step in a large variety of sugar catabolic pathways, which are major suppliers of carbon and energy in heterotrophic species. Specialized sugar kinases that are indispensable for most of these pathways can be utilized as signature enzymes for the reconstruction of carbohydrate utilization machinery from microbial genomic and metagenomic data. Sugar kinases occur in several structurally distinct families with various partially overlapping as well as yet unknown substrate specificities that often cannot be accurately assigned by homology-based techniques. A subsystems-based metabolic reconstruction combined with the analysis of genome context and followed by experimental testing of predicted gene functions is a powerful approach of functional gene annotation. Here we applied this integrated approach for functional mapping of all sugar kinases constituting an extensive and diverse sugar kinome in the thermophilic bacterium Thermotoga maritima. Substrate preferences of 14 kinases mainly from the FGGY and PfkB families were inferred by bioinformatics analysis and biochemically characterized by screening with a panel of 45 different carbohydrates. Most of the analyzed enzymes displayed narrow substrate preferences corresponding to their predicted physiological roles in their respective catabolic pathways. The observed consistency supports the choice of kinases as signature enzymes for genomics-based identification and reconstruction of sugar utilization pathways. Use of the integrated genomic and experimental approach greatly speeds up the identification of the biochemical function of unknown proteins and improves the quality of reconstructed pathways.

N-acetylgalactosamine utilization pathway and regulon in proteobacteria: genomic reconstruction and experimental characterization in Shewanella

We used a comparative genomics approach to reconstruct the N-acetyl-d-galactosamine (GalNAc) and galactosamine (GalN) utilization pathways and transcriptional regulons in Proteobacteria. The reconstructed GalNAc/GalN utilization pathways include multiple novel genes with specific functional roles. Most of the pathway variations were attributed to the amino sugar transport, phosphorylation, and deacetylation steps, whereas the downstream catabolic enzymes in the pathway were largely conserved. The predicted GalNAc kinase AgaK, the novel variant of GalNAc-6-phosphate deacetylase AgaA(II) and the GalN-6-phosphate deaminase AgaS from Shewanella sp. ANA-3 were validated in vitro using individual enzymatic assays and reconstitution of the three-step pathway. By using genetic techniques, we confirmed that AgaS but not AgaI functions as the main GalN-6-P deaminase in the GalNAc/GalN utilization pathway in Escherichia coli. Regulons controlled by AgaR repressors were reconstructed by bioinformatics in most proteobacterial genomes encoding GalNAc pathways. Candidate AgaR-binding motifs share a common sequence with consensus CTTTC that was found in multiple copies and arrangements in regulatory regions of aga genes. This study provides comprehensive insights into the common and distinctive features of the GalNAc/GalN catabolism and its regulation in diverse Proteobacteria.

The FGGY carbohydrate kinase family: insights into the evolution of functional specificities

Function diversification in large protein families is a major mechanism driving expansion of cellular networks, providing organisms with new metabolic capabilities and thus adding to their evolutionary success. However, our understanding of the evolutionary mechanisms of functional diversity in such families is very limited, which, among many other reasons, is due to the lack of functionally well-characterized sets of proteins. Here, using the FGGY carbohydrate kinase family as an example, we built a confidently annotated reference set (CARS) of proteins by propagating experimentally verified functional assignments to a limited number of homologous proteins that are supported by their genomic and functional contexts. Then, we analyzed, on both the phylogenetic and the molecular levels, the evolution of different functional specificities in this family. The results show that the different functions (substrate specificities) encoded by FGGY kinases have emerged only once in the evolutionary history following an apparently simple divergent evolutionary model. At the same time, on the molecular level, one isofunctional group (L-ribulokinase, AraB) evolved at least two independent solutions that employed distinct specificity-determining residues for the recognition of a same substrate (L-ribulose). Our analysis provides a detailed model of the evolution of the FGGY kinase family. It also shows that only combined molecular and phylogenetic approaches can help reconstruct a full picture of functional diversifications in such diverse families.

The protein universe is under constant expansion and is reshaping through multiple duplication, gene losses, lateral gene transfers, and speciation events. Large and functionally heterogeneous protein families that evolve through these processes contain conserved motifs and structural scaffolds, yet their individual members often perform diverse functions. For this reason, the exact functional annotation for their individual members is difficult without detailed analysis of the family. In our study, we performed such a detailed analysis of a particularly heterogeneous FGGY kinase family through the integration of several computational approaches. The combination of phylogenetic and molecular approaches allowed us to precisely assign function to hundreds of proteins, thus reconstructing carbohydrate utilization pathways in almost 200 bacterial species. This analysis also showed that different molecular mechanisms could evolve within a group of isofunctional proteins. Moreover, based on our experience with this specific protein family of FGGY kinases, we believe that our approach can be generally adapted for the analyses of other protein families and that the accumulation of evolutionary models for various families would lead to a better understanding of the protein universe.

Identification of nicotinamide mononucleotide deamidase of the bacterial pyridine nucleotide cycle reveals a novel broadly conserved amidohydrolase family

The pyridine nucleotide cycle is a network of salvage and recycling routes maintaining homeostasis of NAD(P) cofactor pool in the cell. Nicotinamide mononucleotide (NMN) deamidase (EC 3.5.1.42), one of the key enzymes of the bacterial pyridine nucleotide cycle, was originally described in Enterobacteria, but the corresponding gene eluded identification for over 30 years. A genomics-based reconstruction of NAD metabolism across hundreds of bacterial species suggested that NMN deamidase reaction is the only possible way of nicotinamide salvage in the marine bacterium Shewanella oneidensis. This prediction was verified via purification of native NMN deamidase from S. oneidensis followed by the identification of the respective gene, termed pncC. Enzymatic characterization of the PncC protein, as well as phenotype analysis of deletion mutants, confirmed its proposed biochemical and physiological function in S. oneidensis. Of the three PncC homologs present in Escherichia coli, NMN deamidase activity was confirmed only for the recombinant purified product of the ygaD gene. A comparative analysis at the level of sequence and three-dimensional structure, which is available for one of the PncC family member, shows no homology with any previously described amidohydrolases. Multiple alignment analysis of functional and nonfunctional PncC homologs, together with NMN docking experiments, allowed us to tentatively identify the active site area and conserved residues therein. An observed broad phylogenomic distribution of predicted functional PncCs in the bacterial kingdom is consistent with a possible role in detoxification of NMN, resulting from NAD utilization by DNA ligase.

Control of proteobacterial central carbon metabolism by the HexR transcriptional regulator: a case study in Shewanella oneidensis

Bacteria exploit multiple mechanisms for controlling central carbon metabolism (CCM). Thus, a bioinformatic analysis together with some experimental data implicated the HexR transcriptional factor as a global CCM regulator in some lineages of Gammaproteobacteria operating as a functional replacement of the Cra regulator characteristic of Enterobacteriales. In this study, we combined a large scale comparative genomic reconstruction of HexR-controlled regulons in 87 species of Proteobacteria with the detailed experimental analysis of the HexR regulatory network in the Shewanella oneidensis model system. Although nearly all of the HexR-controlled genes are associated with CCM, remarkable variations were revealed in the scale (from 1 to 2 target operons in Enterobacteriales up to 20 operons in Aeromonadales) and gene content of HexR regulons between 11 compared lineages. A predicted 17-bp pseudo-palindrome with a consensus tTGTAATwwwATTACa was confirmed as a HexR-binding motif for 15 target operons (comprising 30 genes) by in vitro binding assays. The negative effect of the key CCM intermediate, 2-keto-3-deoxy-6-phosphogluconate, on the DNA-regulator complex formation was verified. A dual mode of HexR action on various target promoters, repression of genes involved in catabolic pathways and activation of gluconeogenic genes, was for the first time predicted by the bioinformatic analysis and experimentally verified by changed gene expression pattern in S. oneidensis ΔhexR mutant. Phenotypic profiling revealed the inability of this mutant to grow on lactate or pyruvate as a single carbon source. A comparative metabolic flux analysis of wild-type and mutant strains of S. oneidensis using [(13)C]lactate labeling and GC-MS analysis confirmed the hypothesized HexR role as a master regulator of gluconeogenic flux from pyruvate via the transcriptional activation of phosphoenolpyruvate synthase (PpsA).

Comparative genomic reconstruction of transcriptional networks controlling central metabolism in the Shewanella genus

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in bacteria is one of the critical tasks of modern genomics. The Shewanella genus is comprised of metabolically versatile gamma-proteobacteria, whose lifestyles and natural environments are substantially different from Escherichia coli and other model bacterial species. The comparative genomics approaches and computational identification of regulatory sites are useful for the in silico reconstruction of transcriptional regulatory networks in bacteria.

Results

To explore conservation and variations in the Shewanella transcriptional networks we analyzed the repertoire of transcription factors and performed genomics-based reconstruction and comparative analysis of regulons in 16 Shewanella genomes. The inferred regulatory network includes 82 transcription factors and their DNA binding sites, 8 riboswitches and 6 translational attenuators. Forty five regulons were newly inferred from the genome context analysis, whereas others were propagated from previously characterized regulons in the Enterobacteria and Pseudomonas spp.. Multiple variations in regulatory strategies between the Shewanella spp. and E. coli include regulon contraction and expansion (as in the case of PdhR, HexR, FadR), numerous cases of recruiting non-orthologous regulators to control equivalent pathways (e.g. PsrA for fatty acid degradation) and, conversely, orthologous regulators to control distinct pathways (e.g. TyrR, ArgR, Crp).

Conclusions

We tentatively defined the first reference collection of ~100 transcriptional regulons in 16 Shewanella genomes. The resulting regulatory network contains ~600 regulated genes per genome that are mostly involved in metabolism of carbohydrates, amino acids, fatty acids, vitamins, metals, and stress responses. Several reconstructed regulons including NagR for N-acetylglucosamine catabolism were experimentally validated in S. oneidensis MR-1. Analysis of correlations in gene expression patterns helps to interpret the reconstructed regulatory network. The inferred regulatory interactions will provide an additional regulatory constrains for an integrated model of metabolism and regulation in S. oneidensis MR-1.

Large-scale comparative phenotypic and genomic analyses reveal ecological preferences of shewanella species and identify metabolic pathways conserved at the genus level

The use of comparative genomics for the study of different microbiological species has increased substantially as sequence technologies become more affordable. However, efforts to fully link a genotype to its phenotype remain limited to the development of one mutant at a time. In this study, we provided a high-throughput alternative to this limiting step by coupling comparative genomics to the use of phenotype arrays for five sequenced Shewanella strains. Positive phenotypes were obtained for 441 nutrients (C, N, P, and S sources), with N-based compounds being the most utilized for all strains. Many genes and pathways predicted by genome analyses were confirmed with the comparative phenotype assay, and three degradation pathways believed to be missing in Shewanella were confirmed as missing. A number of previously unknown gene products were predicted to be parts of pathways or to have a function, expanding the number of gene targets for future genetic analyses. Ecologically, the comparative high-throughput phenotype analysis provided insights into niche specialization among the five different strains. For example, Shewanella amazonensis strain SB2B, isolated from the Amazon River delta, was capable of utilizing 60 C compounds, whereas Shewanella sp. strain W3-18-1, isolated from deep marine sediment, utilized only 25 of them. In spite of the large number of nutrient sources yielding positive results, our study indicated that except for the N sources, they were not sufficiently informative to predict growth phenotypes from increasing evolutionary distances. Our results indicate the importance of phenotypic evaluation for confirming genome predictions. This strategy will accelerate the functional discovery of genes and provide an ecological framework for microbial genome sequencing projects.

Inference of the transcriptional regulatory network in Staphylococcus aureus by integration of experimental and genomics-based evidence

Transcriptional regulatory networks are fine-tuned systems that help microorganisms respond to changes in the environment and cell physiological state. We applied the comparative genomics approach implemented in the RegPredict Web server combined with SEED subsystem analysis and available information on known regulatory interactions for regulatory network reconstruction for the human pathogen Staphylococcus aureus and six related species from the family Staphylococcaceae. The resulting reference set of 46 transcription factor regulons contains more than 1,900 binding sites and 2,800 target genes involved in the central metabolism of carbohydrates, amino acids, and fatty acids; respiration; the stress response; metal homeostasis; drug and metal resistance; and virulence. The inferred regulatory network in S. aureus includes ∼320 regulatory interactions between 46 transcription factors and ∼550 candidate target genes comprising 20% of its genome. We predicted ∼170 novel interactions and 24 novel regulons for the control of the central metabolic pathways in S. aureus. The reconstructed regulons are largely variable in the Staphylococcaceae: only 20% of S. aureus regulatory interactions are conserved across all studied genomes. We used a large-scale gene expression data set for S. aureus to assess relationships between the inferred regulons and gene expression patterns. The predicted reference set of regulons is captured within the Staphylococcus collection in the RegPrecise database (http://regprecise.lbl.gov).