A computational study of Shewanella oneidensis MR-1: structural prediction and functional inference of hypothetical proteins

The genomes of many organisms have been sequenced in the last 5 years. Typically about 30% of predicted genes from a newly sequenced genome cannot be given functional assignments using sequence comparison methods. In these situations three-dimensional structural predictions combined with a suite of computational tools can suggest possible functions for these hypothetical proteins. Suggesting functions may allow better interpretation of experimental data (e.g., microarray data and mass spectroscopy data) and help experimentalists design new experiments. In this paper, we focus on three hypothetical proteins of Shewanella oneidensis MR-1 that are potentially related to iron transport/metabolism based on microarray experiments. The threading program PROSPECT was used for protein structural predictions and functional annotation, in conjunction with literature search and other computational tools. Computational tools were used to perform transmembrane domain predictions, coiled coil predictions, signal peptide predictions, sub-cellular localization predictions, motif prediction, and operon structure evaluations. Combined computational results from all tools were used to predict roles for the hypothetical proteins. This method, which uses a suite of computational tools that are freely available to academic users, can be used to annotate hypothetical proteins in general.     

Global detection and characterization of hypothetical proteins in Shewanella oneidensis MR-1 using LC-MS based proteomics

The availability of whole genome sequences has enabled the application of powerful tools for assaying global expression patterns in environmentally relevant bacteria such as Shewanella oneidensis MR-1. A large number of genes in prokaryote genomes, including MR-1, have been annotated as hypothetical, indicating that no similar protein has yet been identified in other organisms. Using high-sensitivity MS coupled with accurate mass and time (AMT) tag methodology, 1078 tryptic peptides were collectively detected in MR-1 cultures, 671 of which were unique to their parent protein. Using only these unique tryptic peptides and a minimum of two peptides per protein, we identified, with high confidence, the expression of 258 hypothetical proteins. These proteins ranged from 3.5 to 139 kDa, with 47 being 100 amino acid residues or less. Using a combination of information including detection in cells grown under specific culture conditions, presence within a specific cell fraction, and predictive algorithms such as PSORT and PSORT-B, possible/plausible functions are proposed for some hypothetical proteins. Further, by applying this approach a number of proteins were found not only to be expressed, but only expressed under certain culturing conditions, thereby suggesting function while at the same time isolating several proteins to distinct locales of the cell. These results demonstrate the utility of the AMT tag methodology for comprehensive profiling of the microbial proteome while confirming the expression of a large number of hypothetical genes.     

Exploring the Roles of DNA Methylation in the Metal-Reducing Bacterium Shewanella oneidensis MR-1

We performed whole-genome analyses of DNA methylation in Shewanella oneidensis MR-1 to examine its possible role in regulating gene expression and other cellular processes. Single-molecule real-time (SMRT) sequencing revealed extensive methylation of adenine (N6mA) throughout the genome. These methylated bases were located in five sequence motifs, including three novel targets for type I restriction/modification enzymes. The sequence motifs targeted by putative methyltranferases were determined via SMRT sequencing of gene knockout mutants. In addition, we found that S. oneidensis MR-1 cultures grown under various culture conditions displayed different DNA methylation patterns. However, the small number of differentially methylated sites could not be directly linked to the much larger number of differentially expressed genes under these conditions, suggesting that DNA methylation is not a major regulator of gene expression in S. oneidensis MR-1. The enrichment of methylated GATC motifs in the origin of replication indicates that DNA methylation may regulate genome replication in a manner similar to that seen in Escherichia coli. Furthermore, comparative analyses suggest that many Gammaproteobacteria, including all members of the Shewanellaceae family, may also utilize DNA methylation to regulate genome replication.     

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is ‘open’, with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.     

Anaerobic respiration of Shewanella putrefaciens requires both chromosomal and plasmid-borne genes

Abstract Pleiotropic respiratory mutants, incapable of growth on any electron acceptor other than oxygen, were isolated from two strains of Shewanella putrefaciens (MR-1 and sp200). All anaerobic respiratory functions were restored by complementation of the mutants with specific cloned DNA fragments. Southern hybridization experiments revealed that the fragment that complements the MR-1 mutant was localized on the megaplasmids of both strains, while the fragment that complements the sp200 mutant was chromosomal. Neither of these fragments hybridized with the anaerobic regulatory genes of S. putrefaciens ( etrA ) or E. coli ( fnr ).     

Combined Genomics and Experimental Analyses of Respiratory Characteristics of Shewanella putrefaciens W3-18-1

It has previously been shown that the Shewanella putrefaciens W3-18-1 strain produces remarkably high current in microbial fuel cells (MFCs) and can form magnetite at 0°C. To explore the underlying mechanisms, we developed a genetic manipulation method by deleting the restriction-modification system genes of the SGI1 (Salmonella genome island 1)-like prophage and analyzed the key genes involved in bacterial respiration. W3-18-1 has less respiratory flexibility than the well-characterized S. oneidensis MR-1 strain, as it possesses fewer cytochrome c genes and lacks the ability to oxidize sulfite or reduce dimethyl sulfoxide (DMSO) and timethylamine oxide (TMAO). W3-18-1 lacks the hydrogen-producing Fe-only hydrogenase, and the hydrogen-oxidizing Ni-Fe hydrogenase genes were split into two separate clusters. Two periplasmic nitrate reductases (NapDAGHB and NapDABC) were functionally redundant in anaerobic growth of W3-18-1 with nitrate as the electron acceptor, though napDABC was not regulated by Crp. Moreover, nitrate respiration started earlier in W3-18-1 than in MR-1 (with NapDAGHB only) under microoxic conditions. These results indicate that Shewanella putrefaciens W3-18-1 is well adapted to habitats with higher oxygen levels. Taken together, the results of this study provide valuable insights into bacterial genome evolution.     

Comparative genomics and experimental characterization of N-acetylglucosamine utilization pathway of Shewanella oneidensis

We used a comparative genomics approach implemented in the SEED annotation environment to reconstruct the chitin and GlcNAc utilization subsystem and regulatory network in most proteobacteria, including 11 species of Shewanella with completely sequenced genomes. Comparative analysis of candidate regulatory sites allowed us to characterize three different GlcNAc-specific regulons, NagC, NagR, and NagQ, in various proteobacteria and to tentatively assign a number of novel genes with specific functional roles, in particular new GlcNAc-related transport systems, to this subsystem. Genes SO3506 and SO3507, originally annotated as hypothetical in Shewanella oneidensis MR-1, were suggested to encode novel variants of GlcN-6-P deaminase and GlcNAc kinase, respectively. Reconstitution of the GlcNAc catabolic pathway in vitro using these purified recombinant proteins and GlcNAc-6-P deacetylase (SO3505) validated the entire pathway. Kinetic characterization of GlcN-6-P deaminase demonstrated that it is the subject of allosteric activation by GlcNAc-6-P. Consistent with genomic data, all tested Shewanella strains except S. frigidimarina, which lacked representative genes for the GlcNAc metabolism, were capable of utilizing GlcNAc as the sole source of carbon and energy. This study expands the range of carbon substrates utilized by Shewanella spp., unambiguously identifies several genes involved in chitin metabolism, and describes a novel variant of the classical three-step biochemical conversion of GlcNAc to fructose 6-phosphate first described in Escherichia coli.     

Selection of Shewanella oneidensis MR-1 gene-knockout mutants that adapt to an electrode-respiring condition

Mutants of Shewanella oneidensis MR-1 that adapted to an electrode-respiring condition were selected from a random transposon-insertion mutant library to obtain active current-generating mutants and identify relevant cellular components. The mutants were selected in the presence of an electrode (poised at +0.2 V vs. an Ag/AgCl reference electrode) as the sole electron acceptor, and they were isolated on agar plates. Transposon-insertion sites in the isolated mutants were identified by inverse PCR coupled to sequence analyses. Southern blotting using a transposon probe was also performed to detect mutants that grew abundantly on the electrode. These analyses revealed that in many isolated mutants transposons were inserted in genes relevant to the synthesis of cell-surface structures, including SO_3350 (pilus synthesis), SO_3171 (polysaccharide synthesis), SO_3174 (polysaccharide synthesis), and SO_0165 (general secretion pathway). In microbial fuel cells, some of these (the SO_3350 and SO_4704 mutants) generated higher electrical outputs than wild-type MR-1, while the others generated lower outputs. The results suggest that cell-surface structures have a large influence on microbial current generation.     

DNase-chip: a high-resolution method to identify DNase I hypersensitive sites using tiled microarrays

Mapping DNase I hypersensitive sites is an accurate method of identifying the location of gene regulatory elements, including promoters, enhancers, silencers and locus control regions. Although Southern blots are the traditional method of identifying DNase I hypersensitive sites, the conventional manual method is not readily scalable to studying large chromosomal regions, much less the entire genome. Here we describe DNase-chip, an approach that can rapidly identify DNase I hypersensitive sites for any region of interest, or potentially for the entire genome, by using tiled microarrays. We used DNase-chip to identify DNase I hypersensitive sites accurately from a representative 1% of the human genome in both primary and immortalized cell types. We found that although most DNase I hypersensitive sites were present in both cell types studied, some of them were cell-type specific. This method can be applied globally or in a targeted fashion to any tissue from any species with a sequenced genome.     

金属还原菌希瓦氏菌厌氧呼吸能力及其在环境修复中的研究进展

Shewanella oneidensis MR-1是一种模式金属还原菌,它能够在厌氧条件下,将多种金属化合物和人工合成染料等作为电子受体还原代谢。因此,该菌常常被用于生态修复等研究。厌氧条件下,S.oneidensis MR-1能够将细胞质内或细胞内膜产生的电子通过定位于细胞内膜、细胞膜周质和细胞外膜上的c-血红色素蛋白或还原酶所组成的具有多样性的电子传递系统,最终传递到存在于细菌细胞外环境中的电子受体。通过对多种电子传递过程的介绍,进一步阐明其对污染物修复和纳米材料合成的机理,从而为未来对该类微生物的利用和开发提供更为充分的理论依据。     

产电微生物基因调控网络的构建和特异性通路分析

研究产电微生物胞外电子传递过程和机制,发现与产电效率相关的关键基因、通路和代谢物,是微生物燃料电池研究中的关键技术。为了发现在胞外电子传递过程中起到关键作用的基因以及通路,首先利用比较基因组学的方法,以模式微生物大肠杆菌和同属希瓦氏菌的其他菌株为参考,构建了Shewanella.onedensis MR-1的全基因组基因转录调控网络,大大扩展了目前已知的基因调控关系。然后以此网络为基础,结合基于蛋白质相互作用分析得到的胞外电子传递通路,构建了与胞外电子传递直接传递密切相关的细胞色素C编码基因及其相关调控基因构成的子网络,结合全基因组基因表达数据,研究了特异性条件下胞外电子传递的可能通路和基因调控过程。     

E622, a miniature, virulence-associated mobile element

Miniature inverted terminal repeat elements (MITEs) are nonautonomous mobile elements that have a significant impact on bacterial evolution. Here we characterize E622, a 611-bp virulence-associated MITE from Pseudomonas syringae, which contains no coding region but has almost perfect 168-bp inverted repeats. Using an antibiotic coupling assay, we show that E622 is transposable and can mobilize an antibiotic resistance gene contained between its borders. Its predicted parent element, designated TnE622, has a typical transposon structure with a three-gene operon, consisting of resolvase, integrase, and exeA-like genes, which is bounded by the same terminal inverted repeats as E622. A broader genome level survey of the E622/TnE622 inverted repeats identified homologs in Pseudomonas, Salmonella, Shewanella, Erwinia, Pantoea, and the cyanobacteria Nostoc and Cyanothece, many of which appear to encompass known virulence genes, including genes encoding toxins, enzymes, and type III secreted effectors. Its association with niche-specific genetic determinants, along with its persistence and evolutionary diversification, indicates that this mobile element family has played a prominent role in the evolution of many agriculturally and clinically relevant pathogenic bacteria.     

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.     

High frequency of glucose-utilizing mutants in Shewanella oneidensis MR-1

Shewanella oneidensis MR-1 has conventionally been considered unable to use glucose as a carbon substrate for growth. The genome sequence of S. oneidensis MR-1 however suggests the ability to use glucose. Here, we demonstrate that during initial glucose exposure, S. oneidensis MR-1 quickly and frequently gains the ability to utilize glucose as a sole carbon source, in contrast to wild-type S. oneidensis, which cannot immediately use glucose as a sole carbon substrate. High-performance liquid chromatography and (14)C glucose tracer studies confirm the disappearance in cultures and assimilation and respiration, respectively, of glucose. The relatively short time frame with which S. oneidensis MR-1 gained the ability to use glucose raises interesting ecological implications.     

Functional specificity of extracellular nucleases of Shewanella oneidensis MR-1

Bacterial species such as Shewanella oneidensis MR-1 require extracellular nucleolytic activity for the utilization of extracellular DNA (eDNA) as a source of nutrients and for the turnover of eDNA as a structural matrix component during biofilm formation. We have previously characterized two extracellular nucleases of S. oneidensis MR-1, ExeM and ExeS. Although both are involved in biofilm formation, they are not specifically required for the utilization of eDNA as a nutrient. Here we identified and characterized EndA, a third extracellular nuclease of Shewanella. The heterologously overproduced and purified protein was highly active and rapidly degraded linear and supercoiled DNAs of various origins. Divalent metal ions (Mg(2+) or Mn(2+)) were required for function. endA is cotranscribed with phoA, an extracellular phosphatase, and is not upregulated upon phosphostarvation. Deletion of endA abolished both extracellular degradation of DNA by S. oneidensis MR-1 and the ability to use eDNA as a sole source of phosphorus. PhoA is not strictly required for the exploitation of eDNA as a nutrient. The activity of EndA prevents the formation of large cell aggregates during planktonic growth. However, in contrast to the findings for ExeM, endA deletion had only minor effects on biofilm formation. The findings strongly suggest that the extracellular nucleases of S. oneidensis exert specific functions required under different conditions.