Genome-wide protein localization prediction strategies for gram negative bacteria
BMC Genomics volume 12, Article number: S1 (2011)
Genome-wide prediction of protein subcellular localization is an important type of evidence used for inferring protein function. While a variety of computational tools have been developed for this purpose, errors in the gene models and use of protein sorting signals that are not recognized by the more commonly accepted tools can diminish the accuracy of their output.
As part of an effort to manually curate the annotations of 19 strains of Shewanella, numerous insights were gained regarding the use of computational tools and proteomics data to predict protein localization. Identification of the suite of secretion systems present in each strain at the start of the process made it possible to tailor-fit the subsequent localization prediction strategies to each strain for improved accuracy. Comparisons of the computational predictions among orthologous proteins revealed inconsistencies in the computational outputs, which could often be resolved by adjusting the gene models or ortholog group memberships. While proteomic data was useful for verifying start site predictions and post-translational proteolytic cleavage, care was needed to distinguish cellular versus sample processing-mediated cleavage events. Searches for lipoprotein signal peptides revealed that neither TatP nor LipoP are designed for identification of lipoprotein substrates of the twin arginine translocation system and that the +2 rule for lipoprotein sorting does not apply to this Genus. Analysis of the relationships between domain occurrence and protein localization prediction enabled identification of numerous location-informative domains which could then be used to refine or increase confidence in location predictions. This collective knowledge was used to develop a general strategy for predicting protein localization that could be adapted to other organisms.
Improved localization prediction accuracy is not simply a matter of developing better computational algorithms. It also entails gathering key knowledge regarding the host architecture and translocation machinery and associated substrate recognition via experimentation and integration of diverse computational analyses from many proteins and, where possible, that are derived from different species within the same genus.
Knowledge of the subcellular localization of proteins can provide important insights into protein function and thus is particularly useful in the annotation of genomes and the identification of candidate proteins having functions of interest. For example, microbial proteins that are secreted outside the cell are expected to perform functions associated with cell-cell communication and competition, hydrolysis of membrane impermeable polymers, or creating extracellular structures that enable cell motility, attachment to surfaces, or passage of materials between cells. The discovery of novel surface-localized proteins is useful for the development of drug targets, identification of microbial biomarkers and factors contributing to host invasion, and discovery of more efficient enzymes for use in bioprocesses associated with the breakdown of membrane-impermeable polymers, such as those released during the processing of plant materials for alternative fuel production. In some instances, unexpected localization of proteins belonging to a well studied functional class can lead to exciting new discoveries of cellular function. For example, the discovery that c-type cytochromes associated with Mn(IV) and Fe(III) reduction were localized to the cell surface of Shewanella oneidensis MR-1  rather than the inner membrane or periplasm where respiratory proteins are typically found, initiated a whole new field of research in extracellular respiratory metabolism.
A wide variety of computational tools have been developed as a rapid, inexpensive means to predict protein localization using only amino acid sequence information. New tools continue to be developed with improved accuracy or specificity making it difficult to decide which one(s) to use for genome-wide prediction of protein locations. The primary improvements to predictive accuracy center on the identification of the substrates of the Sec inner membrane export system, which is responsible for translocation of the majority of extra-cytoplasmic proteins across the cytoplasmic membrane in bacteria  and the Tat inner membrane export system which translocates a smaller number of proteins in a pre-folded state . However, bacteria with dual membranes also encode additional machinery for export of proteins from the cytoplasm, for inserting them in the outer membrane, or secreting them beyond the outer membrane. The protein substrates of these systems carry N- or C-terminal signal peptides that are distinct from those recognized by Sec and Tat, or lack them all together, thus requiring the application of alternative computational tools or approaches to identify them. Consequently, prediction of protein localization at the genome scale requires combining multiple tools/methods to account for substrates of both the common export systems, such as Sec, and the less frequently used export or secretion systems.
In this report, we present lessons learned while curating protein localization predictions in 19 strains belonging to the gram negative Genus Shewanella and a generalized workflow (Figure 1) for conducting these analyses that incorporates computational predictions of signal peptide occurrence, subcellular localization, protein domain content, and function with experimental data. The combined genomes of these 19 sequenced strains encode an exceptionally diverse set of secretion systems, including all of named types except the type 4a secretion system (T4aSS), and thus this workflow serves as a useful model for developing strain-specific protein localization prediction workflows in other organisms.
Results and discussion
Assessment of the cell architecture
Prior to applying available bioinformatics tools to predict protein localization, it is important to first establish what types of subcellular compartments are present in the organism of interest. The information can then be used to develop a strain-specific strategy for predicting protein localization at the genome scale. Electron microscopy and the genome annotation are useful resources for determining the compartmental organization of the host, but are limited to detection of structures present under the conditions used to generate the sample. However, when supplemented by information garnered from genome annotations this limitation can be overcome. In the sequenced shewanellae, manual curation of the genome annotation suggested that 1) most of the strains harbor at least one bacteriophage within their genomes, some of which have been observed as distinct entities in stressed cultures cells [4, 5] and 2) under selected growth conditions S. benthica and S. putrefaciens strains CN-32 and W3-18-1  will produce cytoplasmic microcompartments that house specific enzymes and associated reactions that benefit from the resulting secluded environment . These observations and sequence-based predictions should be taken into account when predicting protein localization. Bacteriophages encode viral structural proteins that are not components of the cell and, therefore, not appropriate targets for predicting subcellular localization. The genes that encode these structural proteins are frequently co-localized in operons  and can often be identified through blast analysis against domains/proteins stored in the Aclame database . Proteins likely to be encapsulated in microcompartments, on the other hand, can be identified by searching for proteins that exist only in organisms encoding microcompartment structural proteins (identified by hits to pfam00936) and frequently are encoded in the same neighborhood with them.
Identification of protein export and secretion systems
Once the sites that proteins are expected to localize to have been identified, one should proceed with identifying the suite of protein translocation systems that are encoded in the strain(s) of interest. In gram negative organisms, the export of proteins from the cytoplasm to the periplasm is mediated by the inner membrane Sec  or Tat  translocases while their subsequent insertion in the outer membrane requires the Bam  and sometimes Lol  systems. Since these systems are broadly conserved their key components can readily be identified by searching for orthologs of their respective universally conserved protein components (Table 1). For strains harboring double-stranded DNA or RNA phage, an additional route across the inner membrane is expected for export of the endolysin that initiates cell lysis. This translocase is encoded by the phage genome, usually adjacent to the endolysin gene, and comprised of a single protein (holin) which is a small inner membrane protein having a C-terminus enriched in basic amino acids .
In addition to these export and sorting systems, gram negative bacteria may also encode protein secretion systems, named T1SS-T8SS, that translocate proteins to sites beyond the outer membrane . Secretion systems are often poorly annotated by automated pipelines due to the fact that certain components of different classes of secretion systems (e.g. T2SS and T4SS components) have significant sequence similarity to one another while others, that belong to the same class and that are functionally equivalent, have little similarity to one another (e.g. pilin proteins). In addition, many secretion systems have not yet been characterized and/or informative domains that detect their signature components have yet to been defined and deposited in public databases. Fortunately, the genes encoding the key components of these systems are typically co-localized on the genome and thus one can often use genome context analysis to readily identify their constituents and assign them to appropriate secretion classes.
In Shewanella, identification of the outer membrane channel-forming component of these systems (also called secretins or ushers) via domain analysis proved to be an excellent starting point for finding genomic loci that encode secretion systems. Using the 18 domains listed in Table 2, it was possible to identify the genomic loci that encoded 176 putative extracellular secretion systems in 19 sequenced Shewanella. Only 95 of these putative secretion systems belonged to the core secretion systems present in every strain, demonstrating the importance of considering these strain-specific differences when predicting protein localization. With the exception of TolC, all of the secretins were encoded in the same genomic loci as the other secretion system components, making it easier to identify and annotate the other genes associated with these apparatus. The lone exception was a subset of the T1SS apparatus, which were predicted to depend on a common TolC secretin that is also responsible for efflux of non-protein substrates, such as drugs and heavy metals. These T1SS were identified by searching for proteins with domains such as TIGR01843 (HlyD family) or TIGR01842, TIGR01846, TIGR03375, and pfam0341 (PrtD, HlyB, and LssB; bacteriocin exporter families) which identify the membrane fusion and permease components of T1SS systems, respectively. In Shewanella, each of the loci identified also encoded candidate T1SS substrates, which are typically large proteins lacking a signal peptide with an overall amino acid composition typical of extracellular proteins .
While suitable for detecting many of the secretion systems, the domains listed in Table 2 were not able to detect all of the predicted outer membrane protein translocases in the sequenced shewanellae, requiring that other approaches are taken to identify them. For example, protein localization predictions (described below) and comparative genome context analysis can be used to identify commonly occurring genomic loci that encode putative extracellular proteins along with putative outer membrane or lipoproteins. Other types of functional evidence (e.g. domain content, sequence similarity, and literature searches for experimental data on similar proteins) can then be gathered and reviewed for further clues that are indicative of protein secretion machinery. This approach led to the discovery of a conserved five gene locus in two Shewanella that includes proteins (previously annotated as hypothetical) with similarity to the recently identified components of the Fap amyloid fiber .
T5SS systems, in which the secretin and extracellular function are encoded in the same protein were particularly difficult to confidently identify since the channel forming domain of these systems are highly variable in sequence and currently only detectable by two domains, PF03797 and PF03895 . A review of the literature revealed a new T5SS subclass (T5dSS) that is present in all of the sequenced Shewanella and lacks these domains , instead having C-terminal domains (PF07244 and PF01103) that are characteristic of the BamA component of the outer membrane protein assembly complex and an N-terminal patatin domain, which is frequently found in extracellular proteins. The orthologous Shewanella proteins were all predicted to have a Sec signal peptide by SignalP and to reside in the outer membrane by Bomp, but predicted to localize extracellularly by PsortB and Subloc or to a mixture of outer membrane and extracellular environment by Cello and SosuiGramN. Phobius also detected a signal peptide, but suggested that a single transmembrane span remains at the C-terminus. This region matches the TIGR03501 gamma proteobacterial enzyme C-terminal transmembrane domain, an extracellular location informative domain that is predicted to be proteolytically removed prior to protein secretion (Dan Haft, personal communication). These observations suggest that additional novel T5SS can potentially be identified by searching for proteins with similar mixed evidence of location. Another feature to look out for is the occurrence of exceptionally long Sec leaders known to occur in some T5SS proteins [19, 20]. Since its length may preclude its detection by computational tools designed to detect signal peptides (see below), manual inspection of candidate dual domain T5SS translocases for Sec leaders may be necessary.
Detection of signal peptidases and signal peptides
Once the suite of export and secretion systems present in the strain of interest are identified, a review of the recent literature is warranted to determine whether conserved sequence features are expected in their substrates. Most protein localization systems recognize conserved motifs encoded at either the N- or C-terminus of their substrates (Figure 2). Substrates of Sec, the predominant pathway for inner membrane protein translocation, are received in their unfolded conformation and have a characteristic N-terminal signal peptide  that is removed after export by either signal peptidase I (LepB)  or, in the case of lipoprotein substrates, by signal peptidase II (LspA) . Popular localization predictors, such as SignalP [24–26] and PsortB , are designed to detect only LepB processed Sec substrates. Therefore, tools such as LipoP  or Lipo  must be used to identify lipoprotein substrates of the Sec translocator. TatP  was developed to identify signal peptides present in substrates of the Tat exporter, which is responsible for translocation of folded proteins, many of which bind redox cofactors . However, this tool is unable to detect lipoprotein substrates, as it has only recently been recognized they could be substrates of this system. Shewanella sp. are known to use Tat to translocate the molybdopterin-binding lipoprotein subunits of the DMSO  and arsenate  reductases across the inner membrane and hence we expected to find characteristic Tat signal peptides with LspA cleavage sites when analyzing their genomes. Since LipoP is unable to detect Tat signal peptides, Tatfind  was used to identify proteins having them and then manually searched for an adjacent LspA cleavage site. A recent genomic survey using an algorithm based on the DOLOP database of lipoproteins  and TatP rules, suggested that lipoproteins are Tat substrates in numerous other organisms as well . It should also be noted that proteins which form complexes with Tat substrates can be exported by Tat even though they lack a Tat signal peptide. This phenomenon has been demonstrated for translocation of multi-subunit enzymes such as hydrogenase . Since hitchhikers are not detected by TatP or Tatfind, they need to be manually identified through searches for proteins that are encoded in the neighborhood of TAT substrates and having amino acid composition characteristic of extra-cytoplasmic proteins and/or functional annotations suggesting participation in multi-subunit enzymes.
The occurrence of alternative signal peptidases are expected in organisms that possess type II or IV secretion systems since they are necessary for maturation of the pilin/pseudopilin components of these systems. T2SS and T4SS peptidases can be detected by searching for proteins belonging to Merops  families A24A and S26, respectively (Table 3). In Shewanella a single peptidase, PilD, processes the pilin components of all three subclasses of T2SS and most of its substrates are identifiable by matches to pfam07963 and/or TIGR02532, while others (e.g., GspK, PilX, PilW) that have an imperfect match to the PilD cleavage site could be detected only by similarity at the N-terminus of the mature protein to other PilD substrates or matches to pfam03934 (GspK). Pili associated with IncJ and IncP conjugative systems were detected by matches to TIGR02758 and pfam04956, respectively, while IncH and IncF pili were recognizable only by homology to previously characterized pili associated with these systems. Class I and II bacteriocins are processed at the N-terminus by a C39 family peptidase whose activity is encoded in the N-terminus of the permease component of the T1SS system responsible for bacteriocin secretion. These small proteins are often missed during automated annotation, but can typically be found in genomic loci encoding this characteristic T1SS by searching for nearby small open reading frames that encode proteins with a characteristic twin glycine signal peptide .
Gene model improvement
Because the detection of signal peptides is an important step in localization prediction, errors in prediction of the 5’ end of a gene can displace or truncate N-terminal signal peptides and thus impact the accuracy of localization predictions. Significant improvements have been made in the ORF calling algorithms since the advent of whole genome sequencing and, therefore, the gene models for genomes produced with the earlier generation ORF calling algorithms can be readily improved by comparing the output of the newer algorithms to those used in the original Genbank deposit, or simply using the newer gene model predictions. The output of several of these newer algorithms (Glimmer v. 3, Prodigal v. 2, GeneMarkHMM-2.6r, and GeneMark-2.5m) are pre-computed and available to the research community via FTP from NCBI Refseq (ftp://ftp.ncbi.nih.gov/genomes/Bacteria).
Another means to improve the gene model is to map the termini of transposons, insertion sequences, and other mobile elements in the genome as we reported previously for S. oneidensis MR-1 . This task is not routinely part of the automated genome annotation process and results can reveal that seemingly intact genes are truncated at their 5’ end or interrupted and hence localization predictions can be erroneous. Identification of mobile elements is facilitated by the use of resources like ISfinder  and ACLAME  that provide information regarding the sites targeted by and characteristics of the termini of insertion elements and prophage, respectively. Programmed recoding of genes, whereby genes are translated by non-standard rules (e.g., programmed ribosomal frameshifting, translational bypassing, and utilization of alternative tRNAs to decode stop codons as an amino acid) can also be missed during automated annotations, sometimes even resulting in their erroneous annotation as pseudogenes. The Recode database (http://recode.ucc.ie/)  has compiled numerous examples of recoded genes and thus provides a useful resource for identifying genes likely to be subject to recoding.
Comparative analysis of the protein size, domain content, and localization predictions among orthologous proteins can also prove useful for identifying errors in gene models. Inconsistencies in these values among orthologous Shewanella proteins could often be eliminated by adjusting gene start/stop positions or membership within a predicted orthologous group. In some cases, inconsistencies suggested that one or more members of the group possessed longer signal peptides than detectable by programs such as SignalP or LipoP or that a proposed signal peptide was more likely an uncleaved N-terminal transmembrane domain. As mentioned earlier, unusually long leaders would be expected in some T5SS autotransporters and the secreted component of T5SS two partner secretion systems since some members of this class have signal peptides that are preceded by an additional charged (n-region) and hydrophobic domain (h-region) .
Proteomic data can prove especially useful for improving the gene model, but there are several caveats to their use in validation of genes models that one should be aware of. Trypsin, which specifically cleaves proteins C-terminal to arginine (R) or lysine (K) residues, is the most common enzyme used to digest proteins into fragments of suitable number, size, and charge for subsequent sequence identification by gel-free mass spectrometric-based methods for global characterization of proteins. The C-terminus of each peptide generated is expected to be an R or K and the N-terminus should map adjacent to an R or K in the parent protein. In theory, the only peptides with ends that do not match these criteria, should result from host-mediated proteolytical processing (e.g. by LepB) of the parent protein prior to its tryptic digestion and thus detection of partially tryptic peptides should be indicative of host-mediated post-translational processing of proteins or incorrect assignment of a start codon. However, in practice, partially tryptic peptides can also result from the harsh conditions associated with sample processing, sample fragmentation during ionization, or erroneous peptide identification [43–45]. Therefore, when using proteome data for identifying the N-terminus of mature proteins it is prudent to consider only partially tryptic peptides that, among all peptides detected, are the ones mapping most closely to the N-terminus of the parent protein. Furthermore, the N-termini of these peptides should map to a site that is consistent with predicted protease cleavage sites. In Shewanella, the most frequently encountered proteolytic processing event detected in Shewanella was due to cleavage by AmpP or Map (both present in all the Shewanella genomes), which remove the N-terminal methionine when it is adjacent to Pro or a small amino acid (Ala, Ser, Gly, Cys, Thr, Pro, or Val), respectively [46, 47]. In most cases where a partially tryptic peptide did not map to position 2 of the parent protein (AmpP or Map processed) the detected partially tryptic peptides mapped to signal cleavage sites predicted by SignalP or TatP. A notable exception was the long signal peptide (68 amino acids) found in the small subunit of the NiFe hydrogenase, an expected TAT substrate whose cleavage was not recognized by TATP (except in 1 out of 17 strains having this protein) but for which validating partially tryptic peptides were detected in 4 different strains of Shewanella (see Additional file 1) (M. Romine, unpublished results).
Global analyses of cellular proteomes by mass spectrometry uses the protein sequences deduced from the genomic sequence for peptide matching and thus peptides that map outside of the defined gene termini go undetected. Therefore, searches of MS-MS spectra against protein sequences derived from translations between all stop codons (stop-to-stop databases) or between each stop codon and the furthest upstream start codon (start-to-stop databases) have also been used to increase the number of identifiable peptides in hopes of validating earlier start sites or missed open reading frames . However, non-standard start codons, such as GTG and TTG, are frequently used in bacteria and archaea, but would not be translated as methionine in stop-to-stop in-silico translations. Therefore, N-terminal peptides produced from proteins whose translation is initiated at alternative start codons would still go undetected and consequently the returns from such an effort are diminished. Furthermore, since these databases are significantly larger, the chance of erroneous peptide matching is significantly increased and thus warrants manually evaluating each peptide mapping outside pre-defined open reading frame, especially when the peptide is infrequently detected in samples analyzed.
Predicting protein localization
A variety of different computational tools are available for predicting subcellular location, but only a few enable batch analysis via a web interface and each has certain limitations. A comparison of subcellular localization or signal peptide predictions produced with popular computational tools having batch analysis available (Table 4) revealed frequent inconsistencies in location prediction or signal peptide detection among members of the same ortholog group, even after adjusting gene models or ortholog group membership. Disagreements in predictions for a single protein were also common among the predictions generated by different tools designed for the same purpose. The extent of the problem is shown in Tables 5 and 6, which compare results of different analyses among each set of 19 proteins belonging to one of the 1990 core ortholog groups in Shewanella. At best, only 70% of the groups had consistent subcellular localization prediction suggested for all its members. A comparison of the predictions produced by PsortB, Cell, and SosuiGramN for all 81,619 predicted proteins revealed that just under half of them (39,538) were consistent in localization prediction. The disagreements generally reflected that some tools are better suited for certain types of predictions (e.g. 75 of the predicted SignalP false positives were due to incorrect classification of proteins having signal peptides cleaved by PilD or LspA) while inconsistencies in predictions among orthologs simply revealed the uncertainly of these predictions.
To address these issues, a decision tree (Figure 3) was developed as a guide for using predictions of the occurrence of sorting signals or location informative domains to support or refute global subcellular location predictions or proteomics data from subcellular fractions. While majority voting could often be used to predict location, significant manual curation of gene models (4,208 proteins so far), ortholog grouping, and location assignments was necessary to resolve conflicts in location evidence gathered. It should be noted that the tools listed in the decision tree simply reflect those that were used in analysis of Shewanella. Those chosen were limited to ones that allow batch analysis on-line and that are more broadly used by other researchers, but are not necessarily the most accurate ones currently available. The tools listed can be replaced or supplemented by other tools (recently reviewed in [49, 50]) that are better suited to the organism of interest or yield improved accuracy, recognizing that some may require local installation of software to make genome-scale analyses feasible.
The prediction schema is initiated with the curation of secretion systems, whose components often have distinct signal peptides that are not recognized by predictors listed or that are secreted during assembly of the machinery. In addition, the structural components of bacteriophage are identified at this stage as they would otherwise often be erroneously predicted to localize to the cell envelope. Next automated searches for signal peptides are conducted, working first on the less common signal peptides associated with lipoproteins and Tat substrates and then followed by searching for transmembrane spans and Sec signal peptides. A comparison of the latter two results assisted in distinguishing signal peptides from transmembrane spans, but the availability of additional information (e.g., expected location of a protein based on annotation, detection of peptides that map at or near the N-terminus) was generally needed for deciding whether the N-terminus was removed versus being retained for anchoring a protein in the membrane. Domain content and functional annotations were used through-out this decision tree to increase the confidence and accuracy of the predictions. Location informative domains were identified by searching for Pfam and TIGRfam domains that consistently occurred only in proteins predicted to localize to the same site and/or had a known association with proteins found in specific subcellular or extracellular compartments. In addition, results of searches for a C-terminal outer membrane localization signature were used to enhance outer membrane location predictions, recognizing that those having large periplasmic domains (e.g. TolC family proteins) are expected to lack these signatures or contain them at internal sites instead. This species-specific C-terminal signature consists of alternating hydrophobic residues at positions 5, 7, and 9 from the C-terminus and a Phe or Tyr at the terminus [51, 52]. Since shewanellae have numerous TonB receptors (620 in 19 genomes) we used their C-termini to develop a Shewanella-specific signature that could be used to search for additional substrates of this system.
Characterization of the protein content of subcellular fractions by mass spectrometry is also a useful type of evidence for assessing protein localization. This information is particularly useful for identifying proteins that are tethered to the membrane via protein-protein or protein-lipid interactions or for condition-specific changes in protein localization which cannot be revealed by analyses of protein sequence content alone. However, results must be interpreted with caution as there can be significant cross-contamination between subcellular fractions which may vary depending on the protocol used to fractionate and analyze the protein content or the cell type being studied. In reviewing the data from LC/MS-MS analysis of S. oneidensis MR-1 subcellular fractions prepared with a sarkosyl-based method, we found that fractions with the greatest abundance of peptides partitioned were usually consistent with the predicted locations of the parent protein with the notable exception that many more lipoproteins partitioned to the inner membrane than expected . Sarkosyl was chosen over other detergents because of it compatibility with high through-put MS-based proteomic analysis and reduced time and labor required to conduct the cellular fractionation. While this detergent has been shown to preferentially solubilize inner membrane proteins  thus allowing efficient separation of inner and outer membranes, it is possible that it also solubilizes the loosely associated outer membrane lipoproteins.
Alternatively, the predicted localization of these proteins is incorrect. The rules for predicting lipoprotein sorting are based on extensive research on Escherichia coli lipoproteins and suggest that lipoproteins with an aspartic acid (D) at position +2 (D+2) of the mature protein are retained in the inner membrane while the remainder are attached to the outer membrane by Lol . However, numerous exceptions have been found in other organisms [56–58] suggesting that these rules likely only apply to enterobacteria. Indeed, our analysis of over 3000 predicted lipoproteins in this Genus revealed a lack of consistency in occurrence of D+2 in orthologs and that only 5 out of 112 efflux pump membrane fusion lipoproteins, which are expected to be anchored to the inner membrane, have D+2. Furthermore, like selected other bacteria [59–62], Shewanella can also localize lipoproteins to the outer face of the outer membrane and thus must use alternative sorting signals. While it is known that the T2aSS machinery is responsible for their surface translocation in Shewanella[32, 63, 64], the characteristics of the sorting signals used are currently unknown. The large number of putative lipoproteins identified in this genus and combined knowledge available regarding their localization (experimentally validated as well as predicted based on function or domain content), however, provided a more sensitive means to search for conserved sequences that are characteristic of surface lipoproteins. In Shewanella such analyses suggest that enrichment in glycine and serine residues coincides with predicted surface localization (Romine, unpublished results). These same amino acids have recently been reported to be enriched in extracellular proteins  and are commonly found in other sorting signals used for secretion of proteins [66, 67].
While the methodological process described here was derived from studies of a Genus that shares many structural and functional features with organisms from which much of our current understanding of translocation models have been developed, the overall strategy described for predicting protein localization should prove useful for studying other microbes as well. Knowledge gathered regarding distinctive architectural features or unusual translocation machinery content (e.g. missing components, duplications) prior to applying automated sequence analysis methods can significantly impact the choice of computational tools to use and subsequent interpretation of the results. Proteomic analyses can be especially useful for confirming predictions or discovering novel sorting signals, while less costly computational localization predictions, conducted at the genome scale, can reveal novel characteristics of an organism that might not be readily derived from functional annotations derived solely from sequence similarity.
Subcellular localization and ortholog grouping predictions (Additional file 2) and associated protein sequences (Additional file 3) that were used to for making calculations provided in tables 5 and 6 are provided in the supplementary material so that interested parties can use them for evaluating their own prediction strategies to those used by the author. However, it should be noted that updates to the gene models and ortholog membership is an ongoing process, with the most current versions available at http://shewanella-knowledgebase.org:8080/Shewanella/. Updated localization predictions are available through the author.
Myers CR, Myers JM: Localization of cytochromes to the outer membrane of anaerobically grown Shewanella putrefaciens MR-1. J Bacteriol. 1992, 174 (11): 3429-3438.
Driessen AJ, Nouwen N: Protein translocation across the bacterial cytoplasmic membrane. Annu Rev Biochem. 2008, 77: 643-667. 10.1146/annurev.biochem.77.061606.160747.
Sargent F: The twin-arginine transport system: moving folded proteins across membranes. Biochem Soc Trans. 2007, 35 (Pt 5): 835-847.
Wang F, Li Q, Xiao X: A novel filamentous phage from the deep-sea bacterium Shewanella piezotolerans WP3 is induced at low temperature. J Bacteriol. 2007, 189 (19): 7151-7153. 10.1128/JB.00569-07.
Qiu X, Sundin GW, Wu L, Zhou J, Tiedje JM: Comparative analysis of differentially expressed genes in Shewanella oneidensis MR-1 following exposure to UVC, UVB, and UVA radiation. J Bacteriol. 2005, 187 (10): 3556-3564. 10.1128/JB.187.10.3556-3564.2005.
Karpinets TV, Obraztsova AY, Wang Y, Schmoyer DD, Kora GH, Park BH, Serres MH, Romine MF, Land ML, Kothe TB, et al: Conserved synteny at the protein family level reveals genes underlying Shewanella species' cold tolerance and predicts their novel phenotypes. Funct Integr Genomics. 2010, 10 (1): 97-110. 10.1007/s10142-009-0142-y.
Yeates TO, Crowley CS, Tanaka S: Bacterial microcompartment organelles: protein shell structure and evolution. Annu Rev Biophys. 2010, 39: 185-205. 10.1146/annurev.biophys.093008.131418.
Hatfull GF: Bacteriophage genomics. Curr Opin Microbiol. 2008, 11 (5): 447-453. 10.1016/j.mib.2008.09.004.
Leplae R, Lima-Mendez G, Toussaint A: ACLAME: a CLAssification of Mobile genetic Elements, update 2010. Nucleic Acids Res. 2010, 38 (Database issue): D57-61.
du Plessis DJ, Nouwen N, Driessen AJ: The Sec translocase. Biochim Biophys Acta. 2010
Knowles TJ, Scott-Tucker A, Overduin M, Henderson IR: Membrane protein architects: the role of the BAM complex in outer membrane protein assembly. Nat Rev Microbiol. 2009, 7 (3): 206-214. 10.1038/nrmicro2069.
Narita S, Tokuda H: An ABC transporter mediating the membrane detachment of bacterial lipoproteins depending on their sorting signals. FEBS Lett. 2006, 580 (4): 1164-1170. 10.1016/j.febslet.2005.10.038.
Young I, Wang I, Roof WD: Phages will out: strategies of host cell lysis. Trends Microbiol. 2000, 8 (3): 120-128. 10.1016/S0966-842X(00)01705-4.
Desvaux M, Hebraud M, Talon R, Henderson IR: Secretion and subcellular localizations of bacterial proteins: a semantic awareness issue. Trends Microbiol. 2009, 17 (4): 139-145. 10.1016/j.tim.2009.01.004.
Nakashima H, Nishikawa K: Discrimination of intracellular and extracellular proteins using amino acid composition and residue-pair frequencies. J Mol Biol. 1994, 238 (1): 54-61. 10.1006/jmbi.1994.1267.
Dueholm MS, Petersen SV, Sonderkaer M, Larsen P, Christiansen G, Hein KL, Enghild JJ, Nielsen JL, Nielsen KL, Nielsen PH, et al: Functional amyloid in Pseudomonas. Mol Microbiol. 2010
Dautin N, Bernstein HD: Protein secretion in gram-negative bacteria via the autotransporter pathway. Annu Rev Microbiol. 2007, 61: 89-112. 10.1146/annurev.micro.61.080706.093233.
Salacha R, Kovacic F, Brochier-Armanet C, Wilhelm S, Tommassen J, Filloux A, Voulhoux R, Bleves S: The Pseudomonas aeruginosa patatin-like protein PlpD is the archetype of a novel Type V secretion system. Environ Microbiol. 2010, 12 (6): 1498-1512.
Hiss JA, Schneider G: Domain organization of long autotransporter signal sequences. Bioinform Biol Insights. 2009, 3: 189-204.
Desvaux M, Cooper LM, Filenko NA, Scott-Tucker A, Turner SM, Cole JA, Henderson IR: The unusual extended signal peptide region of the type V secretion system is phylogenetically restricted. FEMS Microbiol Lett. 2006, 264 (1): 22-30. 10.1111/j.1574-6968.2006.00425.x.
von Heijne G: The signal peptide. J Membr Biol. 1990, 115 (3): 195-201. 10.1007/BF01868635.
Tuteja R: Type I signal peptidase: an overview. Arch Biochem Biophys. 2005, 441 (2): 107-111. 10.1016/j.abb.2005.07.013.
Hayashi S, Wu HC: Lipoproteins in bacteria. J Bioenerg Biomembr. 1990, 22 (3): 451-471. 10.1007/BF00763177.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.
Nielsen H, Krogh A: Prediction of signal peptides and signal anchors by a hidden Markov model. Proc Int Conf Intell Syst Mol Biol. 1998, 6: 122-130.
Nielsen H, Engelbrecht J, Brunak S, von Heijne G: Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. 1997, 10 (1): 1-6. 10.1093/protein/10.1.1.
Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, Dao P, Sahinalp SC, Ester M, Foster LJ, et al: PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010, 26 (13): 1608-1615. 10.1093/bioinformatics/btq249.
Juncker AS, Willenbrock H, Von Heijne G, Brunak S, Nielsen H, Krogh A: Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci. 2003, 12 (8): 1652-1662. 10.1110/ps.0303703.
Berven FS, Karlsen OA, Straume AH, Flikka K, Murrell JC, Fjellbirkeland A, Lillehaug JR, Eidhammer I, Jensen HB: Analysing the outer membrane subproteome of Methylococcus capsulatus (Bath) using proteomics and novel biocomputing tools. Arch Microbiol. 2006, 184 (6): 362-377. 10.1007/s00203-005-0055-7.
Bendtsen JD, Nielsen H, Widdick D, Palmer T, Brunak S: Prediction of twin-arginine signal peptides. BMC Bioinformatics. 2005, 6: 167-10.1186/1471-2105-6-167.
Palmer T, Sargent F, Berks BC: Export of complex cofactor-containing proteins by the bacterial Tat pathway. Trends Microbiol. 2005, 13 (4): 175-180. 10.1016/j.tim.2005.02.002.
Gralnick JA, Vali H, Lies DP, Newman DK: Extracellular respiration of dimethyl sulfoxide by Shewanella oneidensis strain MR-1. Proc Natl Acad Sci U S A. 2006, 103 (12): 4669-4674. 10.1073/pnas.0505959103.
Saltikov CW, Newman DK: Genetic identification of a respiratory arsenate reductase. Proc Natl Acad Sci U S A. 2003, 100 (19): 10983-10988. 10.1073/pnas.1834303100.
Rose RW, Bruser T, Kissinger JC, Pohlschroder M: Adaptation of protein secretion to extremely high-salt conditions by extensive use of the twin-arginine translocation pathway. Mol Microbiol. 2002, 45 (4): 943-950. 10.1046/j.1365-2958.2002.03090.x.
Babu MM, Priya ML, Selvan AT, Madera M, Gough J, Aravind L, Sankaran K: A database of bacterial lipoproteins (DOLOP) with functional assignments to predicted lipoproteins. J Bacteriol. 2006, 188 (8): 2761-2773. 10.1128/JB.188.8.2761-2773.2006.
Shruthi H, Babu MM, Sankaran K: TAT-pathway-dependent lipoproteins as a niche-based adaptation in prokaryotes. J Mol Evol. 2010, 70 (4): 359-370. 10.1007/s00239-010-9334-2.
Rodrigue A, Chanal A, Beck K, Muller M, Wu LF: Co-translocation of a periplasmic enzyme complex by a hitchhiker mechanism through the bacterial tat pathway. J Biol Chem. 1999, 274 (19): 13223-13228. 10.1074/jbc.274.19.13223.
Rawlings ND: A large and accurate collection of peptidase cleavages in the MEROPS database. Database (Oxford). 2009, bap015-
Dirix G, Monsieurs P, Dombrecht B, Daniels R, Marchal K, Vanderleyden J, Michiels J: Peptide signal molecules and bacteriocins in Gram-negative bacteria: a genome-wide in silico screening for peptides containing a double-glycine leader sequence and their cognate transporters. Peptides. 2004, 25 (9): 1425-1440. 10.1016/j.peptides.2003.10.028.
Romine MF, Carlson TS, Norbeck AD, McCue LA, Lipton MS: Identification of mobile elements and pseudogenes in the Shewanella oneidensis MR-1 genome. Appl Environ Microbiol. 2008, 74 (10): 3257-3265. 10.1128/AEM.02720-07.
Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M: ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34 (Database issue): D32-36.
Bekaert M, Firth AE, Zhang Y, Gladyshev VN, Atkins JF, Baranov PV: Recode-2: new design, new search tools, and many more genes. Nucleic Acids Res. 2010, 38 (Database issue): D69-74.
Duncan MW, Aebersold R, Caprioli RM: The pros and cons of peptide-centric proteomics. Nat Biotechnol. 2010, 28 (7): 659-664. 10.1038/nbt0710-659.
Olsen JV, Ong SE, Mann M: Trypsin cleaves exclusively C-terminal to arginine and lysine residues. Mol Cell Proteomics. 2004, 3 (6): 608-614. 10.1074/mcp.T400003-MCP200.
Gupta N, Tanner S, Jaitly N, Adkins JN, Lipton M, Edwards R, Romine M, Osterman A, Bafna V, Smith RD, et al: Whole proteome analysis of post-translational modifications: applications of mass-spectrometry for proteogenomic annotation. Genome Res. 2007, 17 (9): 1362-1377. 10.1101/gr.6427907.
Frottin F, Martinez A, Peynot P, Mitra S, Holz RC, Giglione C, Meinnel T: The proteomics of N-terminal methionine cleavage. Mol Cell Proteomics. 2006, 5 (12): 2336-2349. 10.1074/mcp.M600225-MCP200.
Yoshimoto T, Orawski AT, Simmons WH: Substrate specificity of aminopeptidase P from Escherichia coli: comparison with membrane-bound forms from rat and bovine lung. Arch Biochem Biophys. 1994, 311 (1): 28-34. 10.1006/abbi.1994.1204.
Kim W, Silby MW, Purvine SO, Nicoll JS, Hixson KK, Monroe M, Nicora CD, Lipton MS, Levy SB: Proteomic detection of non-annotated protein-coding genes in Pseudomonas fluorescens Pf0-1. PLoS One. 2009, 4 (12): e8455-10.1371/journal.pone.0008455.
Choo KH, Tan TW, Ranganathan S: A comprehensive assessment of N-terminal signal peptides prediction methods. BMC Bioinformatics. 2009, 10 Suppl 15: S2-
Gardy JL, Brinkman FS: Methods for predicting bacterial protein subcellular localization. Nat Rev Microbiol. 2006, 4 (10): 741-751. 10.1038/nrmicro1494.
Robert V, Volokhina EB, Senf F, Bos MP, Van Gelder P, Tommassen J: Assembly factor Omp85 recognizes its outer membrane protein substrates by a species-specific C-terminal motif. PLoS Biol. 2006, 4 (11): e377-10.1371/journal.pbio.0040377.
Struyve M, Moons M, Tommassen J: Carboxy-terminal phenylalanine is essential for the correct assembly of a bacterial outer membrane protein. J Mol Biol. 1991, 218 (1): 141-148. 10.1016/0022-2836(91)90880-F.
Brown RN, Romine MF, Schepmoes AA, Smith RD, Lipton MS: Mapping the subcellular proteome of Shewanella oneidensis MR-1 using sarkosyl-based fractionation and LC-MS/MS protein identification. J Proteome Res. 2010, 9 (9): 4454-4463. 10.1021/pr100215h.
Filip C, Fletcher G, Wulff JL, Earhart CF: Solubilization of the cytoplasmic membrane of Escherichia coli by the ionic detergent sodium-lauryl sarcosinate. J Bacteriol. 1973, 115 (3): 717-722.
Gennity JM, Inouye M: The protein sequence responsible for lipoprotein membrane localization in Escherichia coli exhibits remarkable specificity. J Biol Chem. 1991, 266 (25): 16458-16464.
Narita S, Tokuda H: Amino acids at positions 3 and 4 determine the membrane specificity of Pseudomonas aeruginosa lipoproteins. J Biol Chem. 2007, 282 (18): 13372-13378. 10.1074/jbc.M611839200.
Schulze RJ, Zuckert WR: Borrelia burgdorferi lipoproteins are secreted to the outer surface by default. Mol Microbiol. 2006, 59 (5): 1473-1484. 10.1111/j.1365-2958.2006.05039.x.
Wiker HG, Wilson MA, Schoolnik GK: Extracytoplasmic proteins of Mycobacterium tuberculosis - mature secreted proteins often start with aspartic acid and proline. Microbiology. 2000, 146 (Pt 7): 1525-1533.
Sauvonnet N, Pugsley AP: Identification of two regions of Klebsiella oxytoca pullulanase that together are capable of promoting beta-lactamase secretion by the general secretory pathway. Mol Microbiol. 1996, 22 (1): 1-7. 10.1111/j.1365-2958.1996.tb02650.x.
Bryksin AV, Tomova A, Godfrey HP, Cabello FC: BmpA is a surface-exposed outer-membrane protein of Borrelia burgdorferi. FEMS Microbiol Lett. 2010, 309 (1): 77-83.
Roussel-Jazede V, Jongerius I, Bos MP, Tommassen J, van Ulsen P: NalP-mediated proteolytic release of lactoferrin-binding protein B from the meningococcal cell surface. Infect Immun. 2010, 78 (7): 3083-3089. 10.1128/IAI.01193-09.
Haake DA, Matsunaga J: Leptospira: a spirochaete with a hybrid outer membrane. Mol Microbiol. 2010
Shi L, Deng S, Marshall MJ, Wang Z, Kennedy DW, Dohnalkova AC, Mottaz HM, Hill EA, Gorby YA, Beliaev AS, et al: Direct involvement of type II secretion system in extracellular translocation of Shewanella oneidensis outer membrane cytochromes MtrC and OmcA. J Bacteriol. 2008, 190 (15): 5512-5516. 10.1128/JB.00514-08.
DiChristina TJ, Moore CM, Haller CA: Dissimilatory Fe(III) and Mn(IV) reduction by Shewanella putrefaciens requires ferE, a homolog of the pulE (gspE) type II protein secretion gene. J Bacteriol. 2002, 184 (1): 142-151. 10.1128/JB.184.1.142-151.2002.
Smith D, Chapman M: Economical evolution: microbes reduce the synthetic cost of extracellular proteins. MBio. 2010, 1 (2): e00131-e001310.
Delepelaire P: Type I secretion in gram-negative bacteria. Biochim Biophys Acta. 2004, 1694 (1-3): 149-161. 10.1016/j.bbamcr.2004.05.001.
Samudrala R, Heffron F, McDermott J: Accurate prediction of secreted substrates and identification of a conserved putative secretion signal for type III secretion systems. PLoS Pathog. 2009, 5 (4): e1000375-10.1371/journal.ppat.1000375.
Konstantinidis KT, Serres MH, Romine MF, Rodrigues JL, Auchtung J, McCue LA, Lipton MS, Obraztsova A, Giometti CS, Nealson KH, et al: Comparative systems biology across an evolutionary gradient within the Shewanella genus. Proc Natl Acad Sci U S A. 2009, 106 (37): 15909-15914. 10.1073/pnas.0902000106.
Lipton MS, Romine MF, Monroe ME, Elias DA, Pasa-Tolic L, Anderson GA, Anderson DJ, Fredrickson J, Hixson KK, Masselon C, et al: AMT tag approach to proteomic characterization of Deinococcus radiodurans and Shewanella oneidensis. Methods Biochem Anal. 2006, 49: 113-134.
This research was supported by the U.S. Department of Energy (DOE), Office of Biological and Environmental Research (BER), as part of BER’s Genomic Science Program (GSP). This contribution originates from the GSP Foundational Scientific Focus Area (FSFA) at the Pacific Northwest National Laboratory (PNNL). The Pacific Northwest National Laboratory is operated for the DOE by Battelle Memorial Institute under Contract DE-AC05-76RLO 1830. I would like to thank Tatiana Karpinets, Guru Kora, Denise Schmoyer, and Michael Lueze for developing the ortholog and genome editors that I use for curating gene models and ortholog groups and also Mustafa Syed who conducted some of the automated localization predictions. In addition, I would like to thank Margrethe Serres for conducting domain analyses and assisting in curating the functional predictions.
This article has been published as part of BMC Genomics Volume 12 Supplement 1, 2011: Validation methods for functional genome annotation. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/12?issue=S1.
The author declares that they have no competing interests.
Electronic supplementary material
Additional file 1: Alignment of N-termini of the periplasmic [Ni-Fe] hydrogenase large subunit, HyaB. Representative N-terminal amino acid sequences from 14 different Shewanella species are shown adjacent to their corresponding locus tag. Residues found in the conserved TAT motif are shown in bold. The predicted N-termini of the mature proteins are underscored. The sequences corresponding the most N-terminal peptide identified in four of these microbes  using the AMT approach  is shown in bold and underscored. (DOCX 13 KB)
Additional file 2: Curated ortholog grouping and location predictions for 81169 proteins predicted to be encoded in the genomes of 19 Shewanella. Pseudogenes are denoted with an asterisk in the locus tag column. (XLSX 2 MB)
Additional file 3: FastA file of proteins predicted to be encoded in the genomes of 19 Shewanella. This file includes translations of pseudogenes, with internal stop codons assigned the value ‘X’. (ZIP 17 MB)
About this article
Cite this article
Romine, M.F. Genome-wide protein localization prediction strategies for gram negative bacteria. BMC Genomics 12 (Suppl 1), S1 (2011). https://doi.org/10.1186/1471-2164-12-S1-S1