PSP: rapid identification of orthologous coding genes under positive selection across multiple closely related prokaryotic genomes
© Su et al.; licensee BioMed Central Ltd. 2013
Received: 24 July 2013
Accepted: 26 December 2013
Published: 27 December 2013
With genomic sequences of many closely related bacterial strains made available by deep sequencing, it is now possible to investigate trends in prokaryotic microevolution. Positive selection is a sub-process of microevolution, in which a particular mutation is favored, causing the allele frequency to continuously shift in one direction. Wide scanning of prokaryotic genomes has shown that positive selection at the molecular level is much more frequent than expected. Genes with significant positive selection may play key roles in bacterial adaption to different environmental pressures. However, selection pressure analyses are computationally intensive and awkward to configure.
Here we describe an open access web server, which is designated as PSP (Positive Selection analysis for Prokaryotic genomes) for performing evolutionary analysis on orthologous coding genes, specially designed for rapid comparison of dozens of closely related prokaryotic genomes. Remarkably, PSP facilitates functional exploration at the multiple levels by assignments and enrichments of KO, GO or COG terms. To illustrate this user-friendly tool, we analyzed Escherichia coli and Bacillus cereus genomes and found that several genes, which play key roles in human infection and antibiotic resistance, show significant evidence of positive selection. PSP is freely available to all users without any login requirement at: http://db-mml.sjtu.edu.cn/PSP/.
PSP ultimately allows researchers to do genome-scale analysis for evolutionary selection across multiple prokaryotic genomes rapidly and easily, and identify the genes undergoing positive selection, which may play key roles in the interactions of host-pathogen and/or environmental adaptation.
KeywordsOrthologous genes Positive selection Synonymous and nonsynonymous substitutions Bacterial microevolution Bacillus cereus Escherichia coli
With the next-generation sequencing data “tsunami” in our midst, sets of closely related prokaryotic genomes suitable for comparative evolutionary studies have been available . To date, well established cases of gene selection have been rare . Big data mining of bacterial genomes has shown that positive selection is more widespread at the molecular level than expected under a restrictive interpretation of the neutral theory . Genome-wide molecular selection analyses, designed to assess selection pressure across the entire genomes of different strains, have attempted to address the role of gene selection in the process of microevolution [4, 5], especially in host-pathogen interactions, and metabolic adaptation to a changing environment (stress, antibiotic). Positive selection studies on model bacteria, such as the species Escherichia coli[4, 6] and Listeria, or the genera Streptococcus[8, 9] and Campylobacter have revealed that positive selection is an essential part of natural selection to fix advantageous mutations, and improves the adaptability of bacteria in a wide range of environmental conditions.
A number of methods have been proposed for detecting positive selection in DNA or protein sequences . The most common approach is to integrate evolutionary features into codon-based models, and to use probability-based theory to estimate the ratio (ω) of nonsynonymous (d N ) and synonymous (d S ) substitutions, such as implemented in the PAML  and FitModel . Estimating the ratio ω gives a measure of selective pressure, indicating neutral evolution (ω = 1), purifying selection (ω < 1) and positive selection (ω > 1). In the model of neutral evolution, the likelihood that a nonsynonymous mutation would go to fixation is the same as that for a synonymous mutation. Purifying selection can result in stabilizing selection through the purging of deleterious variations that arise. Positive selection pressure serves to maintain a given set of adaptive traits that aids in survival.
Several nice tools are currently available, such as IDEA , JCoDA  and WSPMaker . However, they are not set up specifically to examine prokaryotic genomes, and they exhibit two major deficiencies: (i) the evolution selection analysis could be difficult to configure on a local computer for most biologists who are not familiar with phylogenetic or evolutionary theory, and (ii) excessively long computing times for analyzing several genomes at once are prohibitive. In this study, we present an open access web server called PSP (Positive Selection analysis for Prokaryotic genomes) to identify orthologous coding genes under positive selection across closely related prokaryotic genomes. It provides several core functions for in-depth analysis of evolutionary selection: retrieving the orthologous groups, generating codon-delimited and un-gapped alignments, removing recombination, building phylogenetic trees, and estimating ω under different models used by PAML/FitModel. Remarkably, PSP is able to facilitate efficient exploration of the identified orthologous genes at the metabolic pathway level by assignments and enrichments of KO (KEGG Orthology), GO (Gene Ontology) or COG (Clusters of Orthologous Groups) terms. Results are presented in a user-friendly web interface, which provides an efficient visualization of positive selection pressure on each orthologous groups.
Rapid identification of orthologous groups across multiple prokaryotic genomes
The identification of orthologs is an important problem in the field of phylogenetic analyses. Basically, there are three types of relationships between orthologous genes, one-to-one, one-to-many and many-to-many . In the study of Adrian et al., OrthoMCL shows a balanced performance, such as the accuracy, the number of genomes analyzed and usability of the web-interface. Therefore, PSP integrates OrthoMCL to quickly identify the many-to-many orthologous relationship. PSP allows users to select or upload complete sequences and annotation details of closely related genomes for comparison. Users can also simultaneously upload thousands of annotated protein-coding genes as Multi-Fasta formatted files. With the default settings, OrthoMCL recognizes co-ortholog relationships with a BLASTp E-value cutoff of 1e-5 and a minimum of 50% coverage. Then PSP performs homolog grouping using the Markov Cluster algorithm with an inflation value of 1.5. Interestingly, PSP runs on a high-performance server and can accept up to thirty comparator bacterial genomes simultaneously. Moreover, one-to-one orthologous relationship of genes could be identified using reciprocal BLAST best hit, which is the first and most widely used method for automatically establishing orthologous relationships .
Optimization of multiple sequence alignment for automated phylogenetic analysis
The protein-coding genes of the individual orthologous groups can be aligned by using MUSCLE or MAFFT. In connection with the automated phylogenetic analysis, PSP improves the coding sequence alignments by using the two following processes: (i) generation of codon-delimited alignments with ad hoc Perl scripts; (ii) maximization of the un-gapped alignment area by using MaxAlign to remove non-homologous sequences .
Removal of gene recombination off orthologous groups
Recombination (or gene conversion) is 10–50 times more likely to cause changes in nucleotide sequence than mutation . To eliminate the influence of horizontal relationships in the positive selection detection, PSP is able to identify recombination signals among the aligned nucleotide sequences of orthologous groups. PSP performs a statistical test to identify recombination breakpoints by using GeneConv . It pre-defines the g-scale parameter as 1, which allows mismatches within a recombining fragment, and the inner fragment P-values with 10,000 random permutations. In addition, PSP also detects recombination with three other statistical tests implemented with the PhiPack package with default settings, including pairwise homoplasy index (PHI), Max χ2 and neighbor similarity score (NSS) . Finally, if a recombination signal is detected to be significant by all four tests (P-value < 0.05 in each method), the alignment would be splitted into two or more fragments.
Detection of orthologous genes under positive selection
Phylogenetic trees are built by the PHYLIP using maximum parsimony or neighbor-joining method. The trees are also able to be generated with Markov codon models by using CodonPhyML . The evolutionary selection is subsequently implemented in the program PAML or FitModel. Only orthologous groups with enough data (at least 4 protein-coding genes) are used as input to detect positive selection, due to the poor quality of Bayes predictions based on small samples . Because the lack of any methods to deal with alignment gaps properly in both programs, a cutoff of the percentages of sequence have data in PSP is used to filter the alignments column by column. In the PSP server, PAML uses three evolutionary models proposed by Yang et al. (Additional file 1: Table S1): site model, strain-specific branch model and strain-specific site-branch model. In the strain-specific analysis, the branches of selected target strains are specified and referred to as “foreground branches” and the rest as “background branches”, which is a powerful tool to detect the selection pressure during the process of environmental adaptation . The in silico detection of evolutionary selection is computationally intensive, particularly using Bayes empirical Bayes to determine posterior probabilities (PP). PSP, which runs on a high-performance server, is able to rapidly calculate the d N /d S ratios and screen orthologous coding genes under positive selection. Similarly, PSP also can apply switching Markov modulated codon models as implemented in the program FitModel to orthologous coding genes to accurately estimate the strength of selection. To most biologists who are not familiar with phylogenetic or evolutionary theory, PSP pre-defines the evolutionary models as described in Additional file 1: Table S1. PSP is also very flexible to set most key parameters in the PAML/FitModel and run strain-specific analysis freely for the evolutionary researchers. For each pair hypothesis, nested models are calculated by comparing the difference in log likelihood values to a χ2 statistic (LRT) for the detection of positive selection. If there is significant evidence for positive selection of any fragment, the orthologous genes, from which the fragment was separated, are suggested to be under positive selection. Notably, PSP also provides a user-friendly visualization tool for performing evolutionary analysis on orthologous coding genes. The embedded Java applet JalView  reports the PP values for all sites, which is helpful for users to determine nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions. Three-dimensional structural models of the protein of interest are predicted and displayed by HHpred  and Jmol . Additionally, Primer3Plus  is integrated to facilitate design of PCR primers to assay orthologous genes based on similar selection among a panel of strains isolated from the same habitat.
Functional investigation at metabolic pathway level
To explore functions of the identified orthologous coding genes under positive selection, PSP performs KO mapping and GO/COG classification. First, to assign the KO terms, PSP uses the level of sequence identity and ratio of matching length to query length cut-off obtained from BLASTp. For each query against the locally installed KEGG gene database, the simple H a -value homology score  is calculated as follows: H a = i × (l m /l q ), where i is the level of identity between protein sequences in the region with the highest Bit score expressed as a ratio between 0 and 1, l m the length of the highest scoring matching sequence (including gaps), and l q the query length. In this study, the H a -value cutoff of 0.7 was used to assign KO. PSP then enriches metabolic pathways with genes under positive selection by tracing back the hierarchical KO levels. The user can select one or more KEGG-archived closely related genomes (up to 20) as references. PSP calculates the P-value of each pathway category by using a hyper-geometric distribution method . In the same way, PSP provides enrichment analysis for COG functional terms based on RPSBLAST searches against the local CDD database, and GO slims analysis based on GOA database.
In the PSP pipeline, a large number of hypotheses are considered which could result in a high rate of Type-1 error even for a relatively stringent P-value cutoff. To reduce Type-1 errors, PSP corrects the obtained P-values using the Q-value  to produce a Q-value at level of FDR < 0.2.
Results and discussion
The PSP tool is applicable to a wide range of prokaryotic species. In this study, we applied PSP to do genome-wide positive selection analyses in two cases (Escherichia coli and Bacillus cereus) both to benchmark PSP and to illustrate its accuracy and usefulness for the exploration of data.
Positive selection analysis of Escherichia coli
Results of positive selection scanning during the analysis pipeline
Petersen et al.
Chen et al.
Num. of orthologous groups
Num. of orthologous (groups size > = 4)
Num. of genes removed by MaxAlign
Method for detecting recombination
GenConv & PhiPack
GeneConv & Reticulate
Num. of recombination
M1a-M2a & M7-M8
M1a-M2a & Branch-Site Model
Num. of genes under positive selection in E. coli K12
List of additional genes under positive selection across the six Escherichia coli genomes using the M1a-M2a test
Bifunctional cobinamide kinase and cobinamide phosphate guanylyltransferase
Reactivating factor for ethanolamine ammonia lyase
Flagellar filament structural protein
Hsp31 molecular chaperone
Glutamyl-Q tRNA synthetase
Conserved inner membrane protein
Predicted diguanylate cyclase
Putative DNA-binding transcriptional regulator
Genes under positive selection in Bacillus cereusgenomes
Genes that show evidences for positive selection in B. anthracis group
Ribosomal large subunit pseudouridylate synthase D
Glycerol uptake facilitator protein
Protease synthase and sporulation negative regulatory protein
Sensor histidine kinase
Iron (III) transport system permease protein
Multidrug resistance protein
Alkyl hydroperoxide reductase, subunit F
Transcriptional regulator, AraC family
GMPP; mannose-1-phosphate guanylyltransferase
Putative prophage LambdaBa04, major capsid protein
RNA pseudouridine synthase family protein
Sugar-binding transcriptional regulator, LacI family
ABC transporter, permease
ABC transporter, permease
Iron complex transport system permease protein
Nucleoside permease NupC
Mass transport system is another effective way for microorganism to resist antibiotic, although it is probable that they may have other natural physiological functions . There are seven genes involved in the mass transport system, which shows the significant enrichment compared to the references (P-value = 0.0395 in the strain Ames Ancestor). The araJ gene was regarded as nonessential membrane protein of unknown function. Recently, it is believed to belong to a large class of multidrug resistance translocators and in particular to the major facilitator superfamily . The orthologous genes of ptr2, which encode a peptide transporter, present in all B. cereus genomes and are regarded as a multidrug resistance protein. In this research, we found it is under strong positive selection and may have important roles to extrude drugs. Iron acquisition genes are also important contributors to B. anthracis, while iron limitation is a component of host defense against infection . The gene fepD, which shows strong evidence for positive selection, is a virulence-associated gene and involved in iron-chelating ABC uptake systems . The gene afuB was reported to differentially express upon treatment with antibiotic, and is also very important for iron transporter . Porins are also important in interaction with the host immune system and could work as receptors for phages or antibiotics . Aquaglyceroporin glpF selectively conducts the passage of small hydrophilic across the inner membrane of B. cereus. The function of remaining two putative ABC transporters is still unknown, but one of them (COG3127) was predicted to involve in lysophospholipase L1 biosynthesis, which may be a good potential target for new antibiotics.
We also identified two transcriptional repressors (purR and araC) and a two-component system (ntrB-ntrC), which show strong evidences for positive selection. Autoregulation of purR controls the expression of many genes involved in purine biosynthetic pathway in B. subtilis. ntrB, a member of the ntrB-ntrC two-component system, encodes the signal-transducing kinase/phosphatase nitrogen regulator, on the regulated phosphatase activity involved in nitrogen regulation. These genes are known to have multiple functions in different ways, which are likely that the positive selection from the interaction of host-pathogen simultaneously improve the adaptation of microorganisms by acting on the versatile proteins, such as metabolic adaptation.
The PSP web server, which integrates a wide variety of useful analytical and functional tools, has been developed to rapidly identify orthologous coding genes under positive selection across up to thirty user-selected or user-supplied prokaryotic genomes. Hosted by a high-performance server and with easy navigation and flexible input options, we present it as a quick and comprehensive genome microevolution tool for biologists. Remarkably, PSP excludes the effect of gene recombination and incorporates functional investigation at the metabolic pathway level. In the future, we plan to improve the computing power of PSP markedly with Paralleled PAML. The upgraded version will be also able to map the positive selection sites to three-dimensional structures of proteins. We propose that a tool such as PSP will support genome-scale analysis for evolutionary selection, aimed at defining genomic biomarkers of evolutionary lineage, phenotype, pathotype, environmental adaptation and/or disease-association of diverse bacterial species.
Availability and requirements
Project name: PSP:Positive Selection analysis for Prokaryotic genomes.
Operating system(s): CentOS.
License: PSP is available free of charge to academic and non-profit institutions.
Any restrictions to use by non-academics: Please contact authors for commercial use.
This work was supported, in part, by the grants from National Basic Research Program of China (2013CB733901) and from Chinese National Natural Science Foundation (31230002, 31170082, and 31121064).
- Novichkov PS, Ratnere I, Wolf YI, Koonin EV, Dubchak I: ATGC: a database of orthologous genes from closely related prokaryotic genomes and a research platform for microevolution of prokaryotes. Nucleic Acids Res. 2009, 37: D448-D454. 10.1093/nar/gkn684.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Z, Bielawski JP: Statistical methods for detecting molecular adaptation. Trends Ecol Evol. 2000, 15: 496-503. 10.1016/S0169-5347(00)01994-7.View ArticlePubMedGoogle Scholar
- Bustamante CD, Fledel-Alon A, Williamson S, Nielsen R, Hubisz MT, Glanowski S, Tanenbaum DM, White TJ, Sninsky JJ, Hernandez RD, et al: Natural selection on protein-coding genes in the human genome. Nature. 2005, 437: 1153-1157. 10.1038/nature04240.View ArticlePubMedGoogle Scholar
- Chen SL, Hung CS, Xu J, Reigstad CS, Magrini V, Sabo A, Blasiar D, Bieri T, Meyer RR, Ozersky P, et al: Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: a comparative genomics approach. Proc Natl Acad Sci U S A. 2006, 103: 5977-5982. 10.1073/pnas.0600938103.PubMed CentralView ArticlePubMedGoogle Scholar
- Kosiol C, Anisimova M: Selection on the protein-coding genome. Methods Mol Biol. 2012, 856: 113-140. 10.1007/978-1-61779-585-5_5.View ArticlePubMedGoogle Scholar
- Petersen L, Bollback JP, Dimmic M, Hubisz M, Nielsen R: Genes under positive selection in Escherichia coli. Genome Res. 2007, 17: 1336-1343. 10.1101/gr.6254707.PubMed CentralView ArticlePubMedGoogle Scholar
- Orsi RH, Sun Q, Wiedmann M: Genome-wide analyses reveal lineage specific contributions of positive selection and recombination to the evolution of Listeria monocytogenes. BMC Evol Biol. 2008, 8: 233-10.1186/1471-2148-8-233.PubMed CentralView ArticlePubMedGoogle Scholar
- Lefebure T, Stanhope MJ: Evolution of the core and pan-genome of Streptococcus: positive selection, recombination, and genome composition. Genome Biol. 2007, 8: R71-10.1186/gb-2007-8-5-r71.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki H, Lefebure T, Hubisz MJ, Pavinski Bitar P, Lang P, Siepel A, Stanhope MJ: Comparative genomic analysis of the Streptococcus dysgalactiae species group: gene content, molecular adaptation, and promoter evolution. Genome Biol Evol. 2011, 3: 168-185. 10.1093/gbe/evr006.PubMed CentralView ArticlePubMedGoogle Scholar
- Lefebure T, Stanhope MJ: Pervasive, genome-wide positive selection leading to functional divergence in the bacterial genus Campylobacter. Genome Res. 2009, 19: 1224-1232. 10.1101/gr.089250.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Z, Nielsen R: Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol. 2000, 17: 32-43. 10.1093/oxfordjournals.molbev.a026236.View ArticlePubMedGoogle Scholar
- Guindon S, Rodrigo AG, Dyer KA, Huelsenbeck JP: Modeling the site-specific variation of selection patterns along lineages. Proc Natl Acad Sci U S A. 2004, 101: 12957-12962. 10.1073/pnas.0402177101.PubMed CentralView ArticlePubMedGoogle Scholar
- Egan A, Mahurkar A, Crabtree J, Badger JH, Carlton JM, Silva JC: IDEA: Interactive Display for Evolutionary Analyses. BMC Bioinformatics. 2008, 9: 524-10.1186/1471-2105-9-524.PubMed CentralView ArticlePubMedGoogle Scholar
- Steinway SN, Dannenfelser R, Laucius CD, Hayes JE, Nayak S: JCoDA: a tool for detecting evolutionary selection. BMC Bioinformatics. 2010, 11: 284-10.1186/1471-2105-11-284.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee YS, Kim TH, Kang TW, Chung WH, Shin GS: WSPMaker: a web tool for calculating selection pressure in proteins and domains using window-sliding. BMC Bioinformatics. 2008, 9 (Suppl 12): S13-10.1186/1471-2105-9-S12-S13.PubMed CentralView ArticlePubMedGoogle Scholar
- Tatusov RL, Galperin MY, Natale DA, Koonin EV: The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000, 28: 33-36. 10.1093/nar/28.1.33.PubMed CentralView ArticlePubMedGoogle Scholar
- Altenhoff AM, Dessimoz C: Phylogenetic and functional assessment of orthologs inference projects and methods. PLoS Comput Biol. 2009, 5: e1000262-10.1371/journal.pcbi.1000262.PubMed CentralView ArticlePubMedGoogle Scholar
- Gabaldon T: Large-scale assignment of orthology: back to phylogenetics?. Genome Biol. 2008, 9: 235-10.1186/gb-2008-9-10-235.PubMed CentralView ArticlePubMedGoogle Scholar
- Gouveia-Oliveira R, Sackett PW, Pedersen AG: MaxAlign: maximizing usable data in an alignment. BMC Bioinformatics. 2007, 8: 312-10.1186/1471-2105-8-312.PubMed CentralView ArticlePubMedGoogle Scholar
- Feil EJ, Maiden MC, Achtman M, Spratt BG: The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis. Mol Biol Evol. 1999, 16: 1496-1502. 10.1093/oxfordjournals.molbev.a026061.View ArticlePubMedGoogle Scholar
- Sawyer SA: GENECONV: A computer package for the statistical detection of gene conversion. Distributed by the author, Department of Mathematics, Washington University in St Louis. 1999, http://www.math.wustl.edu/~sawyer/geneconv/,Google Scholar
- Bruen T, Bruen T: PhiPack: PHI test and other tests of recombination. 2005, Montreal, Quebec: McGill UniversityGoogle Scholar
- Gil M, Zanetti MS, Zoller S, Anisimova M: CodonPhyML: Fast Maximum Likelihood Phylogeny Estimation under Codon Substitution Models. Mol Biol Evol. 2013, 30: 1270-1280. 10.1093/molbev/mst034.PubMed CentralView ArticlePubMedGoogle Scholar
- Anisimova M, Bielawski JP, Yang Z: Accuracy and power of bayes prediction of amino acid sites under positive selection. Mol Biol Evol. 2002, 19: 950-958. 10.1093/oxfordjournals.molbev.a004152.View ArticlePubMedGoogle Scholar
- Yang Z, Wong WS, Nielsen R: Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.View ArticlePubMedGoogle Scholar
- Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ: Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25: 1189-1191. 10.1093/bioinformatics/btp033.PubMed CentralView ArticlePubMedGoogle Scholar
- Soding J: Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005, 21: 951-960. 10.1093/bioinformatics/bti125.View ArticlePubMedGoogle Scholar
- Herraez A: Biomolecules in the computer: Jmol to the rescue. Biochem Mol Biol Educ. 2006, 34: 255-261. 10.1002/bmb.2006.494034042644.View ArticlePubMedGoogle Scholar
- Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JA: Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007, 35: W71-W74. 10.1093/nar/gkm306.PubMed CentralView ArticlePubMedGoogle Scholar
- Shao Y, He X, Harrison EM, Tai C, Ou HY, Rajakumar K, Deng Z: mGenomeSubtractor: a web-based tool for parallel in silico subtractive hybridization analysis of multiple bacterial genomes. Nucleic Acids Res. 2010, 38: W194-W200. 10.1093/nar/gkq326.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang H, Wu CI: A new method for estimating nonsynonymous substitutions and its applications to detecting positive selection. Mol Biol Evol. 2006, 23: 372-379.View ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong WS, Yang Z, Goldman N, Nielsen R: Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites. Genetics. 2004, 168: 1041-1051. 10.1534/genetics.104.031153.PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan IK, Rogozin IB, Wolf YI, Koonin EV: Essential genes are more evolutionarily conserved than are nonessential genes in bacteria. Genome Res. 2002, 12: 962-968.PubMed CentralView ArticlePubMedGoogle Scholar
- Rogers TJ, Paton JC, Wang H, Talbot UM, Paton AW: Reduced virulence of an fliC mutant of Shiga-toxigenic Escherichia coli O113:H21. Infect Immun. 2006, 74: 1962-1966. 10.1128/IAI.74.3.1962-1966.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Pedersen LB, Murray T, Popham DL, Setlow P: Characterization of dacC, which encodes a new low-molecular-weight penicillin-binding protein in Bacillus subtilis. J Bacteriol. 1998, 180: 4967-4973.PubMed CentralPubMedGoogle Scholar
- Garsin DA: Ethanolamine utilization in bacterial pathogens: roles and regulation. Nature Rev Microbiol. 2010, 8: 290-295. 10.1038/nrmicro2334.View ArticleGoogle Scholar
- Paterson ES, Boucher SE, Lambert IB: Regulation of the nfsA Gene in Escherichia coli by SoxS. J Bacteriol. 2002, 184: 51-58. 10.1128/JB.184.1.51-58.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Tagliabue L, Maciąg A, Antoniani D, Landini P: The yddV‒dos operon controls biofilm formation through the regulation of genes encoding curli fibers’ subunits in aerobically growing Escherichia coli. FEMS Immun Med Microbiol. 2010, 59: 477-484.Google Scholar
- Mujacic M, Baneyx F: Regulation of Escherichia coli hchA, a stress-inducible gene encoding molecular chaperone Hsp31. Mol Microbiol. 2006, 60: 1576-1589. 10.1111/j.1365-2958.2006.05207.x.View ArticlePubMedGoogle Scholar
- Cosgrove K, Coutts G, Jonsson IM, Tarkowski A, Kokai-Kun JF, Mond JJ, Foster SJ: Catalase (KatA) and alkyl hydroperoxide reductase (AhpC) have compensatory roles in peroxide stress resistance and are required for survival, persistence, and nasal colonization in Staphylococcus aureus. J Bacteriol. 2007, 189: 1025-1035. 10.1128/JB.01524-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Succi J, Tenover FC, Koehler TM: Beta-lactamase genes of the penicillin-susceptible Bacillus anthracis Sterne strain. J Bacteriol. 2003, 185: 823-830. 10.1128/JB.185.3.823-830.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Allen HK, Donato J, Wang HH, Cloud-Hansen KA, Davies J, Handelsman J: Call of the wild: antibiotic resistance genes in natural environments. Nature Rev Microbiol. 2010, 8: 251-259. 10.1038/nrmicro2312.View ArticleGoogle Scholar
- Edgar R, Bibi E: MdfA, an Escherichia coli multidrug resistance protein with an extraordinarily broad spectrum of drug recognition. J Bacteriol. 1997, 179: 2274-2280.PubMed CentralPubMedGoogle Scholar
- Maresso AW, Chapa TJ, Schneewind O: Surface protein IsdC and Sortase B are required for heme-iron scavenging of Bacillus anthracis. J Bacteriol. 2006, 188: 8145-8152. 10.1128/JB.01011-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Bhagat N, Virdi JS: Distribution of virulence‒associated genes in Yersinia enterocolitica biovar 1A correlates with clonal groups and not the source of isolation. FEMS Microbiol Lett. 2007, 266: 177-183. 10.1111/j.1574-6968.2006.00524.x.View ArticlePubMedGoogle Scholar
- Melnikow E, Schoenfeld C, Spehr V, Warrass R, Gunkel N, Duszenko M, Selzer PM, Ullrich HJ: A compendium of antibiotic-induced transcription profiles reveals broad regulation of Pasteurella multocida virulence genes. Vet Microbiol. 2008, 131: 277-292. 10.1016/j.vetmic.2008.03.007.View ArticlePubMedGoogle Scholar
- Beijer L, Nilsson RP, Holmberg C, Rutberg L: The glpP and glpF genes of the glycerol regulon in Bacillus subtilis. J Gen Microbiol. 1993, 139: 349-359. 10.1099/00221287-139-2-349.View ArticlePubMedGoogle Scholar
- Saxild HH, Brunstedt K, Nielsen KI, Jarmer H, Nygaard P: Definition of the Bacillus subtilis PurR operator using genetic and bioinformatic tools and expansion of the PurR regulon with glyA, guaC, pbuG, xpt-pbuX, yqhZ-folD, and pbuO. J Bacteriol. 2001, 183: 6175-6183. 10.1128/JB.183.21.6175-6183.2001.PubMed CentralView 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.