A predicted physicochemically distinct sub-proteome associated with the intracellular organelle of the anammox bacterium Kuenenia stuttgartiensis
- Marnix H Medema†1, 8,
- Miaomiao Zhou†2Email author,
- Sacha AFT van Hijum2, 3, 5,
- Jolein Gloerich4,
- Hans JCT Wessels4,
- Roland J Siezen2, 3, 5 and
- Marc Strous1, 6, 7
© Medema et al; licensee BioMed Central Ltd. 2010
Received: 29 January 2010
Accepted: 12 May 2010
Published: 12 May 2010
Anaerobic ammonium-oxidizing (anammox) bacteria perform a key step in global nitrogen cycling. These bacteria make use of an organelle to oxidize ammonia anaerobically to nitrogen (N2) and so contribute ~50% of the nitrogen in the atmosphere. It is currently unknown which proteins constitute the organellar proteome and how anammox bacteria are able to specifically target organellar and cell-envelope proteins to their correct final destinations. Experimental approaches are complicated by the absence of pure cultures and genetic accessibility. However, the genome of the anammox bacterium Candidatus "Kuenenia stuttgartiensis" has recently been sequenced. Here, we make use of these genome data to predict the organellar sub-proteome and address the molecular basis of protein sorting in anammox bacteria.
Two training sets representing organellar (30 proteins) and cell envelope (59 proteins) proteins were constructed based on previous experimental evidence and comparative genomics. Random forest (RF) classifiers trained on these two sets could differentiate between organellar and cell envelope proteins with ~89% accuracy using 400 features consisting of frequencies of two adjacent amino acid combinations. A physicochemically distinct organellar sub-proteome containing 562 proteins was predicted with the best RF classifier. This set included almost all catabolic and respiratory factors encoded in the genome. Apparently, the cytoplasmic membrane performs no catabolic functions. We predict that the Tat-translocation system is located exclusively in the organellar membrane, whereas the Sec-translocation system is located on both the organellar and cytoplasmic membranes. Canonical signal peptides were predicted and validated experimentally, but a specific (N- or C-terminal) signal that could be used for protein targeting to the organelle remained elusive.
A physicochemically distinct organellar sub-proteome was predicted from the genome of the anammox bacterium K. stuttgartiensis. This result provides strong in silico support for the existing experimental evidence for the existence of an organelle in this bacterium, and is an important step forward in unravelling a geochemically relevant case of cytoplasmic differentiation in bacteria. The predicted dual location of the Sec-translocation system and the apparent absence of a specific N- or C-terminal signal in the organellar proteins suggests that additional chaperones may be necessary that act on an as-yet unknown property of the targeted proteins.
Anaerobic ammonium-oxidizing (anammox) bacteria convert ammonium and nitrite into nitrogen and are major players in the biogeochemical nitrogen cycle [1–4]. They comprise a monophyletic taxon within the Planctomycetes phylum. Like other Planctomycetes, they possess an unusual cellular architecture with a diderm cell envelope and a compartmentalized cytoplasm [5, 6]. More specifically, the cells of anammox bacteria contain a single organelle-like intracytoplasmic compartment bounded by a single bilayer membrane. This compartment is known as the anammoxosome, and was proposed to be the site at which the anammox reaction takes place . This reaction is thought to be performed mainly by cytochrome c enzymes . Within anammox cells, such enzymes have been shown to be present exclusively inside the anammoxosome [9, 10].
If indeed the anammoxosome is a separate compartment in which a distinct and substantial part of the proteome is localized, this would present a situation unique to bacteria. In a thorough electron tomographical study it was reported that, unlike for example the magnetosomes of magnetotactic bacteria  and the chlorosomes of green photosynthetic bacteria , the anammoxosome has no detectable membrane links with the cell envelope during its biogenesis . Furthermore, anammoxosomes divide separately from the cell envelope during cell division .
This leads to two questions regarding the cell biology of anammox bacteria: Firstly, which proteins are targeted to the anammoxosome besides the cytochrome c enzymes? Secondly, by what mechanism are these proteins specifically targeted to the anammoxosome?
One possible answer to the second question is that anammoxosomal proteins might contain specific sorting signals such as targeting motifs, domains or signal peptides [15–19]. For example, in Salmonella, several effectors were reported to contain multifunctional motifs or domains that are responsible for translocation and localization of the effector traits . Moreover, some cases have recently been discovered in which modulation of Sec-signal peptide sequences result in different protein localizations [21, 22]. Most strikingly, in cyanobacteria, signal peptides from proteins targeted to the thylakoid differ from signal peptides of proteins targeted to the cell envelope [23–25].
Progress in the experimental investigation of the cell biology of anammox bacteria is slow because these bacteria grow exceptionally slowly (with a doubling time of two weeks), and are not available in pure culture. However, the genome of the anammox bacterium Candidatus "Kuenenia stuttgartiensis" was recently assembled from a community genome .
Using these genome sequence data, it might be possible to answer the first question. Interestingly, it has been shown experimentally that the anammoxosome may be more acidic than both the cytoplasm and the cell envelope . We reasoned that such a physicochemical difference could be reflected in the amino acid composition of the anammoxosomal sub-proteome , and that this difference could be used to predict this sub-proteome in silico. Therefore, a Random forest (RF) classifier was trained on two sets of anammoxosomal (set A) and cell-envelope (set P) protein sequences, constructed based on existing experimental evidence and comparative genomics. The best RF classifier was successful at predicting the targeting of proteins to the anammoxosome. This approach was complemented by the analysis of the encoded protein translocation machinery. Finally, the predicted signal peptides of the two sub-proteomes were analyzed and compared to unravel the molecular basis of protein sorting in anammox bacteria.
Construction of training sets
We constructed two sets of amino acid sequences from soluble proteins with functions that were known to be specific to either the anammoxosome or the cell envelope (Additional file 1: sheet S1). The anammoxosomal set (termed 'set A') contained the sequences of 30 proteins, including the 26 most highly expressed cytochrome c proteins of K. stuttgartiensis[8, 10] and 4 orthologues of these proteins from another anammox bacterium, Candidatus "Scalindua marina" (data kindly provided by M. Jetten and M. Kuypers). The cell-envelope set (termed 'set P') contained the sequences of 59 proteins that were homologous to proteins with an experimentally validated function specific for the periplasm, cell envelope, or extracellular environment. All proteins of both sets contained a predicted N-terminal signal peptide that can putatively be recognized by the Sec translocon. TatFind  and PilFind  predicted no Tat-system or Type IV system-secreted proteins in these two sets.
Training of the Random forest classifier
To detect the overall compositional differences between the anammoxosomal and cell-envelope protein sets, we constructed Random forest (RF) classifiers based on amino acid composition, using set A and set P as the training sets. Among commonly used analytical machine-learning techniques [30–35], the RF algorithm has shown its power in classifying proteins based on noisy amino acid composition [36–40]. To balance class sizes of sets A and P, set P was first randomized into two sets (P1 and P2) to arrive at three equally sized protein sets (see Material and Methods).
Identification of translocated and membrane proteins
Prediction of the anammoxosomal sub-proteome
Composition of the predicted organellar proteome of K. stuttgartiensis.
Examples of the predicted organellar proteins
Cytochrome C proteins
cytochrome c551 peroxidase
cbb3-type cytochrome c oxidase subunit 1 (CcoN)
hepta heme protein
Respiratory complex proteins
cd1 nitrite reductase (NirS)
NAD(P)H:quinone oxidoreductase chain 5
proton-translocating NADH dehydrogenase I chain A (NuoA)
NADH:ubiquinone oxidoreductase subunit M
Ammonium transporter 1
Nitrite Transporter 1 (FocA)
Nitrite/nitrate antiporter (NarK)
Protein translocation system
SecDF-YajC accessory complex (YajC)
SecYEG translocation complex (SecE)
Transmembrane pore (TatC)
Tat signal recognition (TatA/B)
Cytochrome C maturation system
thiol-disulfide oxidoreductase (ResA)
cytochrome c-type biogenesis protein (ResC)
thiol:disulfide interchange protein (DsbD)
N-acetylglucosaminyltransferase (O-GlcNAc transferase)
kinesin light chain KLC (putative)
Hypothetical and other proteins
Importantly, the predicted anammoxosomal sub-proteome formed a functionally consistent and cohesive set. First, the cytochrome c maturation machinery  was predicted to be anammoxosomal as 6 out of the 8 encoding genes were predicted to be targeted to the anammoxosome (Additional file 2: sheet S4), consistent with the recent proteomics study of Karlsson et al. . Secondly, the destination of proteins encoded in the same genetic neighbourhood (i.e. putative operons) was generally consistent. This makes sense because such proteins are usually subunits of a protein complex or otherwise functionally associated (Additional file 2: sheet S5). Thirdly, we found that only 15 of the 562 predicted anammoxosomal proteins (2.7%) had SMART or PFAM functional domains  (e-value < 0.01) which are on functional grounds incompatible with an anammoxosomal location (Additional file 2: sheet S2).
Finally, the predicted anammoxosomal sub-proteome was consistent with the proposed biological role of the anammoxosome [7, 10, 45]. If the main catabolism of anammox bacteria takes place in the anammoxosome, the respiratory complexes should be associated with this compartment. Indeed, all 14 major respiratory complexes encoded in the K. stuttgartiensis genome were predicted to reside in the anammoxosomal membrane (Additional file 2: sheet S6). Moreover, three out of four ammonium transporters, all nitrite transporters, and all nitrite/nitrate antiporters were predicted to be anammoxosomal (Additional file 2: sheet S7). In contrast, importers of essential trace elements and amino acids, as well as multidrug-efflux proteins were predicted to be located on the outside of the cells, as expected. A single putative copper ATP transporter was predicted to be anammoxosomal consistent with the anammoxosomal destination of some enzymes dependent on copper (or other metal cations).
Overall, the RF classifier predicts that in anammox bacteria the cytoplasmic membrane is mainly used for transport and that essentially all catabolic functions (the anammox reaction, respiration and ATP synthesis) are associated with the intracytoplasmic organelle.
Mechanism of protein translocation
The next point we addressed is the molecular basis for protein sorting in anammox bacteria. We reasoned that comparison of predicted protein features to those of reference bacteria could provide the first clues to how such a sorting system could function.
Protein sorting components encoded in the K. stuttgartiensis genome
Candidatus Kuenenia stuttgartiensis homologues of proteins involved in protein sorting
SecYEG translocation complex
SecYEG translocation complex
SecYEG translocation complex
ATPase motor protein
SecDF-YajC accessory complex
SecDF-YajC accessory complex
membrane protein assembly
Tat signal recognition
Signal Peptidase I
Sec signal peptidase
Signal Peptidase II
lipoprotein signal peptidase
Signal Peptidase IV
prepilin signal peptidase
Major SRP subunit
The identified components of the Sec- and Tat-systems appeared canonical, except for the presence of a C-terminal FecR domain (PF04773) in TatC. This signal-transducing domain is absent in any other TatC protein identified so far. Topology predictions of the TatC-FecR protein (by Phobius , TMHMM , and HMMTOP ) unanimously showed that the FecR domain is non-cytoplasmic. Because of the uniqueness of such a domain combination, it is tempting to speculate that the FecR-like domain may somehow have a role in the evolutionary solution that has been found by anammox bacteria for protein sorting to the anammoxosome.
Interestingly, the RF classifier results for Sec components were ambiguous (only 3 out of 6 subunits were predicted to be anammoxosomal), whereas the Tat system was predicted to be completely anammoxosomal (Additional file 2: sheet S8). Consistently, 7 out of the 9 Tat-substrates predicted in the K. stuttgartiensis genome were also predicted to be anammoxosomal. These include a multi-copper oxidase SufI (kuste4301), a putative superoxide dismutase (kustd1303), two Rieske subunits of the bc1 complex (kuste3096 and kuste4569), and a few hypothetical proteins. The NarG nitrate reductase subunit is very probably a false positive, as has been noted earlier .
In conclusion, the encoded protein translocation machinery shows that no duplication of the Sec- or Tat-systems has taken place to facilitate separate translocation routes, and that both systems may be involved in protein sorting towards the anammoxosome.
The role of signal peptides in protein sorting
Regardless of the protein translocation machinery used, the targeting of proteins to specific subcellular locations is often accomplished by modulation of N- or C-terminal signal peptides, in eukaryotes [15–19] and prokaryotes [23–25] alike. Therefore, we compared the N- and C-termini of the two training sets A and P to identify a possible distinctive property or amino acid motif that could be used to differentiate both sets of protein sequences.
An interesting side observation was that the h-regions of the predicted signal peptides from K. stuttgartiensis proteins contained significantly more phenylalanine residues (2.23 on average for both training sets) than found in E. coli TMHs (1.64 on average, Additional file 3: Supplemental Figure S3). We speculate that this difference may be related to the affinity of trans-membrane helices to the unique ladderane membranes of anammox bacteria that have an exceptionally high density to prevent diffusion .
Yet in conclusion, it is unlikely that the N- or C- terminus plays a role in protein targeting of anammox bacteria.
The anammoxosome of anammox bacteria is one of the best documented cases of organellar biogenesis in Bacteria; electron tomography has shown that this bacterial organelle divides separately from the cytoplasmic membrane, and is not connected to this membrane during a complete cell cycle . Cytochrome c proteins were detected exclusively inside the anammoxosome  and an intracytoplasmic pH gradient was shown to exist . Because experimental investigation of anammox bacteria appeared to be difficult and the genome of K. stuttgartiensis has been sequenced recently, an in silico analysis was therefore obviously the next step forwards in unravelling this interesting biological phenomenon.
The anammoxosomal training set consisted only of amino acid sequences of cytochrome c proteins. It is possible that the classification was biased by the uniqueness of the protein types in this set. However, the cytochrome c protein family is only defined by the presence of a conserved CXXCH heme-binding motif (which we excluded from the RF input data) in an alpha-helical domain. Yet the remainder of the amino acid sequences were vastly variable, some polypeptides even containing regions with different folds or domains [63, 64]. In fact, motif searches showed that no conserved sequence patterns could be found in set A except for the heme-binding motif, and the pairwise sequence identities of proteins in this set were all below 80%, with only 11 of them above 50%.
For the prediction of translocated proteins, a combination of 15 existing signal peptide prediction algorithms was used. These 15 signal peptide predictors reported immensely different predictions (Figure 2) on the presence or absence of a signal peptide in K. stuttgartiensis proteins. Considering the fact that the Gram-negative predictors were mainly trained with sequences from Proteobacteria, which are only very distantly related to anammox bacteria, it is highly unlikely that the reported accuracies of these predictors of 91-95% [65–68] can be valid for Kuenenia proteins. Moreover, among all employed algorithms, the SignalP-HMM algorithm trained on eukaryotes showed the highest true-positive rate by predicting signal peptides in 68 out of 69 of the proteins in set A. This indicates that anammox signal peptides are more similar to those of eukaryotes than to those of Gram-negative or Gram-positive model bacteria.
Analysis of the encoded protein translocation machinery provided some clues as to how this bacterium targets translocated proteins to their proper destination. Because this machinery is non-redundant, an additional layer of chaperoning would be required to achieve specificity. Moreover, because no sorting signal was apparent at the N- or C-termini of the proteins, such chaperones could act on the physicochemical characteristics observed in the amino acid sequence. Alternatively, the signal may act at the level of the messenger RNA that could determine the fate of the protein even before translation starts.
In case of the Sec translocase, some subunits (SecE, SecG and YidC) were predicted to be anammoxosomal, while some others (SecY, SecA, SecDF,YajC) were not. According to the presence of the seemingly canonical Sec-signal peptides on both cell-envelope-targeted proteins and anammoxosome-targeted proteins, it is likely that the Sec translocase has a dual localization on both the anammoxosomal and the periplasmic membranes. However, recent studies have also shown that under certain conditions Sec-system exported proteins could as well be translocated by the Tat-system [69, 70], and therefore we cannot exclude the possibility that the Tat system translocates more than just the predicted twin-arginine-motif-carrying substrates.
An alternative hypothesis for anammoxosomal targeting could be primary translocation to the periplasm through the Sec pathway and secondary retro-translocation to the anammoxosome (or vice versa) through vesicles. However, it must be noted that the application of 3D electron microscopy (tomography) did not reveal any such vesicles in growing or dividing cells .
In order to further elucidate the anammox protein targeting problem, more experiments, such as immunolocalization of the Tat- and Sec-translocase subunits and quantitative proteomics approaches [71, 72], comparing protein concentration levels in purified anammox cells and solutions enriched in anammoxosomes (a purification method which has been described earlier ), are required. The present study provides a clear hypothesis to future experiments: in anammox bacteria, catabolism and respiration are strictly organellar, leaving only transport functions for the cytoplasmic membrane.
The anammoxosome of anammox bacteria is one of the best documented cases of organellar biogenesis in bacteria. Experiments have shown that several key enzymes catalyzing the anammox reaction are present exclusively inside the anammoxosome. The present study makes use of physicochemical characteristics of predicted protein sequences to predict a 90% accurate sub-proteome that constitutes this bacterial organelle. Meanwhile, the mechanism of protein sorting remained largely elusive.
The predicted sub-proteome has been deposited into a freely accessible Microsoft-Excel database (Additional file 2).
Genome sequence of Candidatus K. stuttgartiensis
The complete predicted proteome of K. stuttgartiensis was compiled from all annotated protein sequences (4663 ORFs, accessed 11-02-2008) encoded by K. stuttgartiensis genome fragments kustA - kustE (GenBank accession nrs.: CT030148, CT573074, CT573073, CT573072, CT573071). The proteome has been deposited in the peptidome database (accession number PSE111).
Sets of putatively anammoxosome- and cell-envelope-targeted proteins
The set of amino acid sequences of putative anammoxosomal proteins ("set A") was first constructed from the cytochrome c proteins with a peptide coverage of more than 10% in an experimental analysis of the K. stuttgartiensis proteome (Kartal et al., unpublished data). Next, local BlastP searches were performed with the Scalindu marina metagenome using these protein sequences as queries in order to extend the set A. Reciprocal best Blast hits with the K. stuttgartiensis genome that share an identical gene context with the set A proteins were identified as orthologues and added to set A. By similar methods a set of amino acid sequences of putative cell-envelope and/or excreted proteins ("set P") was also constructed. This set consisted of proteins from the K. stuttgartiensis genome with high similarity to proteins with a validated function in the periplasm, cell envelope or extracellular environment [73–89].
Predicted integral transmembrane proteins (predicted using Phobius , combined with manual inspection) were removed from the sets. The translation start sites of the selected proteins were manually checked and corrected when necessary. Pairwise identities of the sequences from both sets were calculated with MatGAT  and redundant protein sequences (with pairwise identity higher than 85%) were removed.
Prediction of signal peptide-carrying proteins in Candidatus K. stuttgartiensis
Signal peptide predictions on the whole K. stuttgartiensis proteome were performed by the Gram-negative, Gram-positive and eukaryote versions of the algorithm PrediSi , SignalP-HMM and SignalP-NN , Gram-negative and Gram-positive versions of Signal-3L  and Signal-CF , and the general versions of RPSP  and Phobius . All positive predictions from all algorithms were combined into a majority vote decision. Trans-membrane helices (TMH) were predicted by Phobius, TMHMM , and HMMTOP . Tat-secreted proteins were identified with TatFind . The same algorithms were also run on all protein sequences of sets A and P (see below). The n-, h-, and c-regions of the predicted signal peptides were determined manually based on expert knowledge with help of the predictions by Phobius and SignalP. No attempt was made to differentiate between signal peptides (which are cleaved by a SPase) and signal anchors (which are not cleaved) for these sets.
Random Forest classification based on mature protein amino acid composition
The Random Forest classifier package (version 4.5-28) [42, 92] from the R environment (version 2.8.1) [42, 93] was used to train RF classifiers for the separation of anammoxosomal (set A) and cell-envelope (set P) proteins.
In order to reduce the bias due to class size imbalance (larger size of set P), this set was randomly separated into two subsets (set P1 and P2) with sizes similar to that of set A. As input for the RF algorithm, features were determined based on the frequencies (occurrences divided by the sequence length) of amino acids, including the frequencies of two-adjacent-amino-acids, from (i) the SPs, (ii) the full-length amino acid sequences and (iii) the mature protein amino acid sequences of sets A, set P1 and set P2 (except cysteines and histidines, which constitute the heme c binding motif). A three-class (P1, P2, and A) RF model was trained with 1000 trees per forest using each set of input data at each round of P set randomization. The randomization training process was repeated 500 times for each set of input data, after which the votes for classes P1 and P2 were pooled into one merged set P and the overall classification or out-of-bag (OOB) errors were calculated.
The best RF model was selected based on overall accuracy and A protein recall. When a tie situation occurred with the overall accuracy, the model with higher A protein recall was preferred.
Identification of Sec signal peptide cleavage sites by mass-spectrometry
Mass-spectrometry experiments were performed to identify the Sec-signal peptide cleavage sites in proteins from sets A and P. A detailed methodology of these experiments is described in the Additional file 3. The resulting mass spectrometric data files were searched against a database containing the K. stuttgartiensis proteins and known contaminants like human keratins and trypsin using Mascot (Matrix Science Inc., USA, version 2.2) . Variant sequences were modified at the N-terminus by deleting amino acid 1 to 50 and were added to the database in order to search for the Sec-signal peptide cleavage site. The resulting peptide hits were validated using an in-house developed script which selects peptides based on peptide score, the number of variable modifications, the expectation value, and the modified delta score. From the list of validated peptides, a non-redundant N-terminal peptide list was manually extracted by the criteria of: (1), whether the protein was identified with ≥ 3 peptides; (2), whether the peptide was the first detectable peptide (based on calculated m/z values of in silico predicted tryptic peptides in relationship with the m/z detection limits of the mass spectrometer) of the protein; (3), whether the peptide was semi-tryptic (with the non-tryptic side at the N-terminus).
Prediction of the anammoxosomal sub-proteome of K. stuttgartiensis
The Candidatus K. stuttgartiensis proteome was separated into the translocated and non-translocated sets by combining several signal-peptide predictors as follows: firstly, a majority vote of the prediction from 15 signal-peptide prediction algorithms (see above) was used to predict sets of proteins containing SPs (511) and proteins without SPs (4152). Secondly, the proteins without SPs were subdivided into sets of soluble cytoplasmic proteins (3724 proteins, set 4) and SP-less transmembrane proteins (427 proteins, set 3), using predictions of TMHMM 2.0 . Finally the SP-containing proteins were considered subdivided into soluble proteins and proteins containing TMHs, by assessing which proteins had TMHs predicted by a Phobius constrained prediction (constraint: N-terminus = signal peptide). This resulted in a predicted set of 344 SP-containing soluble proteins (set 1) and a predicted set of 167 SP-containing transmembrane proteins (set 2) (Additional file 1: sheet S2). The TMHs (including TMH topology) were predicted by either a constrained Phobius search (for sets 1 and 2) or the TMHMM output (for set 3). Protein sequence composition data of subset 1-3 were then used to predict the anammoxosome proteome.
Identification and sequence analysis of the protein translocation system components in the genome of Candidatus Kuenenia stuttgartiensis
Genes encoding translocation-associated proteins SecY, SecE, SecG, SecA, SecDF, YajC, YidC, TatA, TatC, SPI, SPII, SPIV, FtsY, and Ffh were identified in the K. stuttgartiensis genome by finding reciprocal best Blast hits using BlastP  with the well-studied proteins involved in inner membrane translocation in Escherichia coli K12 as queries. Orthologues from S. marina were identified by performing BlastP analysis on a metagenomic database (M. Jetten and M.Kuypers, unpublished results) constructed with a S. marina enrichment culture , using the above-identified K. stuttgartiensis protein sequences as a query.
Homologues of the Sec- and Tat- translocation system components from other bacterial species were found by PSI-Blast  searches using the Escherichia coli K12 proteins as queries on the GenBank http://www.ncbi.nlm.nih.gov/ genomic data from all bacterial phyla. The K. stuttgartiensis protein translocation-associated proteins were aligned to at least 4 non-planctomycete sequences using Muscle .
Analysis of N-terminal signal peptides
Multiple sequence alignments were made by Muscle 3.6 using standard parameters. Motif searches were performed using MEME , first with the criterion of zero or one motif per sequence, then with one motif per sequence. Phylogenetic footprinting  was performed by aligning the signal peptides of K. stuttgartiensis and S. marina orthologues from set A together with their putative orthologues (based on reciprocal best BlastP hits) from the sequences of planctomycete KSU-1 obtained by Shimamura et al. [100, 101]. The protein physicochemical properties, including the prevalence of general amino acid classes (ILV, FWY, AILVMFWYC, AGS, ST, GNP, DE, DN, KR, EQ, DENQ, HKR, DENQHKR, and DENQHKRST), aliphatic index , and GRAVY index (grand average of hydropathy, according to the Kyte/Doolittle scale)  of the proteins were calculated by custom Python scripts. Statistical analysis on these parameters was performed by calculating the average, standard deviation, and average deviation of the data.
List of Abbreviations
hidden Markov model
open reading frame
out-of-bag error estimate.
M Zhou is funded by the BioRange programme of the Netherlands Bioinformatics Centre (NBIC), which is supported by a BSIK grant through the Netherlands Genomics Initiative (NGI). We thank L. van Niftrik, B. Dutilh, H. op den Camp, and A. Pol for stimulating discussions. B. Kartal and K. van de Pas Schoonen are acknowledged for maintaining the anammox enrichment cultures, for performing the protein extraction for proteomics and providing unpublished proteome data. M. Strous is supported by a VIDI grant from NWO. Preliminary genome data of Candidatus "Scalindua marina" were kindly provided by M. Jetten and M. Kuypers.
- Strous M, Jetten MS: Anaerobic oxidation of methane and ammonium. AnnuRevMicrobiol. 2004, 58: 99-117.
- Jetten MS, Niftrik LV, Strous M, Kartal B, Keltjens JT, Op den Camp HJ: Biochemistry and molecular biology of anammox bacteria. Crit RevBiochemMolBiol. 2009, 1-20.
- Kuypers MM, Lavik G, Woebken D, Schmid M, Fuchs BM, Amann R, Jorgensen BB, Jetten MS: Massive nitrogen loss from the Benguela upwelling system through anaerobic ammonium oxidation. ProcNatlAcadSciUSA. 2005, 102 (18): 6478-6483.View Article
- Op den Camp HJ, Kartal B, Guven D, van Niftrik LA, Haaijer SC, Star van der WR, KT vdP-S, Cabezas A, Ying Z, Schmid MC: Global impact and application of the anaerobic ammonium-oxidizing (anammox) bacteria. BiochemSocTrans. 2006, 34 (Pt 1): 174-178.
- Lindsay MR, Webb RI, Strous M, Jetten MS, Butler MK, Forde RJ, Fuerst JA: Cell compartmentalisation in planctomycetes: novel types of structural organisation for the bacterial cell. ArchMicrobiol. 2001, 175 (6): 413-429.
- Fuerst JA: Intracellular compartmentation in planctomycetes. AnnuRevMicrobiol. 2005, 59: 299-328.
- van Niftrik LA, Fuerst JA, Sinninghe Damste JS, Kuenen JG, Jetten MS, Strous M: The anammoxosome: an intracytoplasmic compartment in anammox bacteria. FEMS MicrobiolLett. 2004, 233 (1): 7-13. 10.1016/j.femsle.2004.01.044.View Article
- Strous M, Pelletier E, Mangenot S, Rattei T, Lehner A, Taylor MW, Horn M, Daims H, Bartol-Mavel D, Wincker P: Deciphering the evolution and metabolism of an anammox bacterium from a community genome. Nature. 2006, 440 (7085): 790-794. 10.1038/nature04647.PubMedView Article
- Karlsson R, Karlsson A, Backman O, Johansson BR, Hulth S: Identification of key proteins involved in the anammox reaction. FEMS Microbiol Lett. 2009, 297 (1): 87-94. 10.1111/j.1574-6968.2009.01677.x.PubMedView Article
- van Niftrik LA, Geerts WJ, van Donselaar EG, Humbel BM, Webb RI, Fuerst JA, Verkleij AJ, Jetten MS, Strous M: Linking ultrastructure and function in four genera of anaerobic ammonium-oxidizing bacteria: cell plan, glycogen storage, and localization of cytochrome C proteins. JBacteriol. 2008, 190 (2): 708-717. 10.1128/JB.01449-07.View Article
- Komeili A: Molecular mechanisms of magnetosome formation. AnnuRevBiochem. 2007, 76: 351-366.
- Blankenship RE, Olson JM, Miller M: Antenna Complexes from Green Photosynthetic Bacteria. Edited by: Blankenship RE, Madigan MT, Bauer CE. 1995, Kluwer Academic Publishers, 399-435.
- van Niftrik L, Geerts WJ, van Donselaar EG, Humbel BM, Yakushevska A, Verkleij AJ, Jetten MS, Strous M: Combined structural and chemical analysis of the anammoxosome: a membrane-bounded intracytoplasmic compartment in anammox bacteria. J Struct Biol. 2008, 161 (3): 401-410. 10.1016/j.jsb.2007.05.005.PubMedView Article
- van Niftrik L, Geerts WJ, van Donselaar EG, Humbel BM, Webb RI, Harhangi HR, Camp HJ, Fuerst JA, Verkleij AJ, Jetten MS: Cell division ring, a new cell division protein and vertical inheritance of a bacterial organelle in anammox planctomycetes. Mol Microbiol. 2009, 73 (6): 1009-1019. 10.1111/j.1365-2958.2009.06841.x.PubMedView Article
- Rapoport TA, Jungnickel B, Kutay U: Protein transport across the eukaryotic endoplasmic reticulum and bacterial inner membranes. AnnuRevBiochem. 1996, 65: 271-303.
- Subramani S: Protein import into peroxisomes and biogenesis of the organelle. AnnuRevCell Biol. 1993, 9: 445-478. 10.1146/annurev.cb.09.110193.002305.
- Jans DA, Xiao CY, Lam MH: Nuclear targeting signal recognition: a key control point in nuclear transport?. Bioessays. 2000, 22 (6): 532-544. 10.1002/(SICI)1521-1878(200006)22:6<532::AID-BIES6>3.0.CO;2-O.PubMedView Article
- Haucke V, Schatz G: Import of proteins into mitochondria and chloroplasts. Trends Cell Biol. 1997, 7 (3): 103-106. 10.1016/S0962-8924(96)10052-0.PubMedView Article
- Pfanner N, Geissler A: Versatility of the mitochondrial protein import machinery. NatRevMolCell Biol. 2001, 2 (5): 339-349.
- Brown NF, Szeto J, Jiang X, Coombes BK, Finlay BB, Brumell JH: Mutational analysis of Salmonella translocated effector members SifA and SopD2 reveals domains implicated in translocation, subcellular localization and function. Microbiology. 2006, 152 (Pt 8): 2323-2343. 10.1099/mic.0.28995-0.PubMedView Article
- Hegde RS, Bernstein HD: The surprising complexity of signal sequences. Trends BiochemSci. 2006, 31 (10): 563-571. 10.1016/j.tibs.2006.08.004.View Article
- Carlsson F, Stalhammar-Carlemalm M, Flardh K, Sandin C, Carlemalm E, Lindahl G: Signal sequence directs localized secretion of bacterial surface proteins. Nature. 2006, 442 (7105): 943-946. 10.1038/nature05021.PubMedView Article
- Rajalahti T, Huang F, Klement MR, Pisareva T, Edman M, Sjostrom M, Wieslander A, Norling B: Proteins in different Synechocystis compartments have distinguishing N-terminal features: a combined proteomics and multivariate sequence analysis. J Proteome Res. 2007, 6 (7): 2420-2434. 10.1021/pr0605973.PubMedView Article
- Aldridge C, Spence E, Kirkilionis MA, Frigerio L, Robinson C: Tat-dependent targeting of Rieske iron-sulphur proteins to both the plasma and thylakoid membranes in the cyanobacterium Synechocystis PCC6803. Mol Microbiol. 2008, 70 (1): 140-150. 10.1111/j.1365-2958.2008.06401.x.PubMedView Article
- Yusa F, Steiner JM, Loffelhardt W: Evolutionary conservation of dual Sec translocases in the cyanelles of Cyanophora paradoxa. BMC Evol Biol. 2008, 8: 304-10.1186/1471-2148-8-304.PubMed CentralPubMedView Article
- Star van der WRL, Dijkema C, de Waard P, Picioreanu C, Strous M, van Loosdrecht MCM: An intracellular pH gradient in the anammox bacterium Kuenenia stuttgartiensis as evaluated by 31P NMR. Appl Microbiol Biotechnol. 2010, 86 (1): 311-7. 10.1007/s00253-009-2309-9.PubMed CentralPubMedView Article
- Pascal G, Medigue C, Danchin A: Persistent biases in the amino acid composition of prokaryotic proteins. Bioessays. 2006, 28 (7): 726-738. 10.1002/bies.20431.PubMedView Article
- 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. MolMicrobiol. 2002, 45 (4): 943-950.
- PilFind. [http://signalfind.org/pilfind.html]
- Kaundal R, Raghava GP: RSLpred: an integrative system for predicting subcellular localization of rice proteins combining compositional and evolutionary information. Proteomics. 2009, 9 (9): 2324-2342. 10.1002/pmic.200700597.PubMedView Article
- Shi JY, Zhang SW, Pan Q, Cheng YM, Xie J: Prediction of protein subcellular localization by support vector machines using multi-scale energy and pseudo amino acid composition. AminoAcids. 2007, 33 (1): 69-74.
- Hua S, Sun Z: Support vector machine approach for protein subcellular localization prediction. Bioinformatics. 2001, 17 (8): 721-728. 10.1093/bioinformatics/17.8.721.PubMedView Article
- Heider D, Appelmann J, Bayro T, Dreckmann W, Held A, Winkler J, Barnekow A, Borschbach M: A computational approach for the identification of small GTPases based on preprocessed amino acid sequences. Technol Cancer Res Treat. 2009, 8 (5): 333-341.PubMedView Article
- Arnold R, Brandmaier S, Kleine F, Tischler P, Heinz E, Behrens S, Niinikoski A, Mewes HW, Horn M, Rattei T: Sequence-based prediction of type III secreted proteins. PLoS Pathog. 2009, 5 (4): e1000376-10.1371/journal.ppat.1000376.PubMed CentralPubMedView Article
- Garg P, Sharma V, Chaudhari P, Roy N: SubCellProt: predicting protein subcellular localization using machine learning approaches. In Silico Biol. 2009, 9 (1-2): 35-44.PubMed
- Kumar KK, Pugalenthi G, Suganthan PN: DNA-Prot: identification of DNA binding proteins from protein sequence information using random forest. J Biomol Struct Dyn. 2009, 26 (6): 679-686.PubMedView Article
- Han P, Zhang X, Feng ZP: Predicting disordered regions in proteins using the profiles of amino acid indices. BMC Bioinformatics. 2009, 10 (Suppl 1): S42-10.1186/1471-2105-10-S1-S42.PubMed CentralPubMedView Article
- Krajewski C, King DG: Molecular divergence and phylogeny: rates and patterns of cytochrome b evolution in cranes. Mol Biol Evol. 1996, 13 (1): 21-30.PubMedView Article
- Cornish-Bowden A: The random character of protein evolution and its effects on the reliability of phylogenetic information deduced from amino acid sequences and compositions. Biochem J. 1980, 191 (2): 349-354.PubMed CentralPubMedView Article
- Statnikov A, Wang L, Aliferis C: A comprehensive comparison of random forests and support vector machines for microarray-based cancer classification. BMC Bioinformatics. 2008, 9 (1): 319-10.1186/1471-2105-9-319.PubMed CentralPubMedView Article
- Segal MA: Machine Learning Benchmarks and Random Forest Regression. Technical Report, Center for Bioinformatics & Molecular Biostatistics, University of California, San Francisco. 2004
- Leo B: Random Forests. Mach Learn. 2001, 45 (1): 5-32. 10.1023/A:1010933404324.View Article
- Kranz R, Lill R, Goldman B, Bonnard G, Merchant S: Molecular mechanisms of cytochrome c biogenesis: three distinct systems. MolMicrobiol. 1998, 29 (2): 383-396.
- Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res. 2000, 28 (1): 263-266. 10.1093/nar/28.1.263.PubMed CentralPubMedView Article
- van Niftrik LA, Geerts WJ, van Donselaar EG, Humbel BM, Yakushevska A, Verkleij AJ, Jetten MS, Strous M: Combined structural and chemical analysis of the anammoxosome: A membrane-bounded intracytoplasmic compartment in anammox bacteria. JStructBiol. 2008, 161 (3): 401-10.
- Veenendaal AK, van der Does DC, Driessen AJ: The protein-conducting channel SecYEG. BiochimBiophysActa. 2004, 1694 (1-3): 81-95.
- Stephenson K: Sec-dependent protein translocation across biological membranes: evolutionary conservation of an essential protein transport pathway (review). MolMembrBiol. 2005, 22 (1-2): 17-28.
- Gold VA, Duong F, Collinson I: Structure and function of the bacterial Sec translocon. MolMembrBiol. 2007, 24 (5-6): 387-394.
- Jongbloed JD, van der Ploeg R, van Dijl JM: Bifunctional TatA subunits in minimal Tat protein translocases. Trends Microbiol. 2006, 14 (1): 2-4. 10.1016/j.tim.2005.11.001.PubMedView Article
- Tuteja R: Type I signal peptidase: an overview. ArchBiochemBiophys. 2005, 441 (2): 107-111.
- Narita S, Matsuyama S, Tokuda H: Lipoprotein trafficking in Escherichia coli. ArchMicrobiol. 2004, 182 (1): 1-6.
- LaPointe CF, Taylor RK: The type 4 prepilin peptidases comprise a novel family of aspartic acid proteases. JBiolChem. 2000, 275 (2): 1502-1510.
- Kall L, Krogh A, Sonnhammer EL: Advantages of combined transmembrane topology and signal peptide prediction--the Phobius web server. Nucleic Acids Res. 2007, W429-W432. 10.1093/nar/gkm256. 35 Web Server
- Krogh A, Larsson B, von Heijne HG, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. JMolBiol. 2001, 305 (3): 567-580.
- Tusnady GE, Simon I: The HMMTOP transmembrane topology prediction server. Bioinformatics. 2001, 17 (9): 849-850. 10.1093/bioinformatics/17.9.849.PubMedView Article
- Martinez-Espinosa RM, Dridge EJ, Bonete MJ, Butt JN, Butler CS, Sargent F, Richardson DJ: Look on the positive side! The orientation, identification and bioenergetics of 'Archaeal' membrane-bound nitrate reductases. FEMS Microbiol Lett. 2007, 276 (2): 129-139. 10.1111/j.1574-6968.2007.00887.x.PubMedView Article
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.PubMed CentralPubMedView Article
- Durbin R, Eddy SR, Krogh A, Mitchison G: Biological sequence analysis: probabilistic models of proteins and nucleic acids. 1998, Cambridge University PressView Article
- Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, W369-W373. 10.1093/nar/gkl198. 34 Web Server
- Sinninghe Damste JS, Strous M, Rijpstra WI, Hopmans EC, Geenevasen JA, van Duin AC, van Niftrik LA, Jetten MS: Linearly concatenated cyclobutane lipids form a dense bacterial membrane. Nature. 2002, 419 (6908): 708-712. 10.1038/nature01128.PubMedView Article
- Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157 (1): 105-132. 10.1016/0022-2836(82)90515-0.PubMedView Article
- Ikai A: Thermostability and aliphatic index of globular proteins. J Biochem. 1980, 88 (6): 1895-1898.PubMed
- Allen JW, Barker PD, Ferguson SJ: A cytochrome b562 variant with a c-type cytochrome CXXCH heme-binding motif as a probe of the Escherichia coli cytochrome c maturation system. J Biol Chem. 2003, 278 (52): 52075-52083. 10.1074/jbc.M307196200.PubMedView Article
- Allen JW, Ferguson SJ: Variation of the axial haem ligands and haem-binding motif as a probe of the Escherichia coli c-type cytochrome maturation (Ccm) system. Biochem J. 2003, 375 (Pt 3): 721-728. 10.1042/BJ20030752.PubMed CentralPubMedView Article
- Chou KC, Shen HB: Signal-CF: a subsite-coupled and window-fusing approach for predicting signal peptides. BiochemBiophysResCommun. 2007, 357 (3): 633-640.
- Bendtsen JD, Nielsen H, von Heijne HG, Brunak S: Improved prediction of signal peptides: SignalP 3.0. JMolBiol. 2004, 340 (4): 783-795.
- Hiller K, Grote A, Scheer M, Munch R, Jahn D: PrediSi: prediction of signal peptides and their cleavage positions. Nucleic Acids Res. 2004, W375-W379. 10.1093/nar/gkh378. 32 Web Server
- Shen HB, Chou KC: Signal-3L: A 3-layer approach for predicting signal peptides. BiochemBiophysResCommun. 2007, 363 (2): 297-303.
- Esbelin J, Martin C, Forano E, Mosoni P: Differential translocation of green fluorescent protein fused to signal sequences of Ruminococcus albus cellulases by the Tat and Sec pathways of Escherichia coli. FEMS Microbiol Lett. 2009, 294 (2): 239-244. 10.1111/j.1574-6968.2009.01576.x.PubMedView Article
- Kouwen TR, Ploeg van der R, Antelmann H, Hecker M, Homuth G, Mader U, van Dijl JM: Overflow of a hyper-produced secretory protein from the Bacillus Sec pathway into the Tat pathway for protein secretion as revealed by proteogenomics. Proteomics. 2009, 9 (4): 1018-1032. 10.1002/pmic.200800580.PubMedView Article
- Lilley KS, Dupree P: Methods of quantitative proteomics and their application to plant organelle characterization. JExpBot. 2006, 57 (7): 1493-1499.
- Nakamura T, Oda Y: Mass spectrometry-based quantitative proteomics. BiotechnolGenetEng Rev. 2007, 24: 147-163.
- Henriques AO, Glaser P, Piggot PJ, Moran CP: Control of cell shape and elongation by the rodA gene in Bacillus subtilis. MolMicrobiol. 1998, 28 (2): 235-247.
- Lutkenhaus J, Addinall SG: Bacterial cell division and the Z ring. AnnuRevBiochem. 1997, 66: 93-116.
- Maeda H, Morihara K: Serralysin and related bacterial proteinases. Methods Enzymol. 1995, 248: 395-413. full_text.PubMedView Article
- Held KG, Postle K: ExbB and ExbD do not function independently in TonB-dependent energy transduction. JBacteriol. 2002, 184 (18): 5170-5173. 10.1128/JB.184.18.5170-5173.2002.View Article
- Johnson JM, Church GM: Alignment and structure prediction of divergent protein families: periplasmic and outer membrane proteins of bacterial efflux pumps. JMolBiol. 1999, 287 (3): 695-715.
- Zgurskaya HI, Nikaido H: Multidrug resistance mechanisms: drug efflux across two membranes. MolMicrobiol. 2000, 37 (2): 219-225.
- Nikaido H, Zgurskaya HI: AcrAB and related multidrug efflux pumps of Escherichia coli. JMolMicrobiolBiotechnol. 2001, 3 (2): 215-218.
- Postle K, Kadner RJ: Touch and go: tying TonB to transport. MolMicrobiol. 2003, 49 (4): 869-882.
- Stevenson G, Andrianopoulos K, Hobbs M, Reeves PR: Organization of the Escherichia coli K-12 gene cluster responsible for production of the extracellular polysaccharide colanic acid. JBacteriol. 1996, 178 (16): 4885-4893.
- Wimley WC: The versatile beta-barrel membrane protein. CurrOpinStructBiol. 2003, 13 (4): 404-411.
- Bateman A, Bycroft M: The structure of a LysM domain from E. coli membrane-bound lytic murein transglycosylase D (MltD). JMolBiol. 2000, 299 (4): 1113-1119.
- Tam R, Saier MH: Structural, functional, and evolutionary relationships among extracellular solute-binding receptors of bacteria. MicrobiolRev. 1993, 57 (2): 320-346.
- Zhang X, Studier FW: Multiple roles of T7 RNA polymerase and T7 lysozyme during bacteriophage T7 infection. JMolBiol. 2004, 340 (4): 707-730.
- Goffin C, Ghuysen JM: Multimodular penicillin-binding proteins: an enigmatic family of orthologs and paralogs. MicrobiolMolBiolRev. 1998, 62 (4): 1079-1093.
- Sara M, Sleytr UB: S-Layer proteins. JBacteriol. 2000, 182 (4): 859-868. 10.1128/JB.182.4.859-868.2000.View Article
- Leahy DJ, Hendrickson WA, Aukhil I, Erickson HP: Structure of a fibronectin type III domain from tenascin phased by MAD analysis of the selenomethionyl protein. Science. 1992, 258 (5084): 987-991. 10.1126/science.1279805.PubMedView Article
- Matsunaga J, Barocchi MA, Croda J, Young TA, Sanchez Y, Siqueira I, Bolin CA, Reis MG, Riley LW, Haake DA: Pathogenic Leptospira species express surface-exposed proteins belonging to the bacterial immunoglobulin superfamily. MolMicrobiol. 2003, 49 (4): 929-945.
- Campanella JJ, Bitincka L, Smalley J: MatGAT: an application that generates similarity/identity matrices using protein or DNA sequences. BMCBioinformatics. 2003, 4: 29-10.1186/1471-2105-4-29.
- Plewczynski D, Slabinski L, Tkacz A, Kajan L, Holm L, Ginalski K, Rychlewski L: The RPSP: Web server for prediction of signal peptides. Polymer. 2007, 48 (19): 5493-5496. 10.1016/j.polymer.2007.07.039.View Article
- Liaw A, Wiener M: Classification and Regression by randomForest. R News. 2002, 2 (3): 18-22.
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2009
- Perkins DN, Pappin DJ, Creasy DM, Cottrell JS: Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. 1999, 20 (18): 3551-3567. 10.1002/(SICI)1522-2683(19991201)20:18<3551::AID-ELPS3551>3.0.CO;2-2.PubMedView Article
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. JMolBiol. 1990, 215 (3): 403-410.
- Vossenberg van de J, Rattray JE, Geerts W, Kartal B, van Niftrik NL, van Donselaar EG, Sinninghe Damsté JS, Strous M, Jetten MS: Enrichment and characterization of marine anammox bacteria associated with global nitrogen gas production. EnvironMicrobiol. 2008, 10 (11): 3120-9.
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView Article
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMCBioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.
- Cliften P, Sudarsanam P, Desikan A, Fulton L, Fulton B, Majors J, Waterston R, Cohen BA, Johnston M: Finding functional features in Saccharomyces genomes by phylogenetic footprinting. Science. 2003, 301 (5629): 71-76. 10.1126/science.1084337.PubMedView Article
- Shimamura M, Nishiyama T, Shigetomo H, Toyomoto T, Kawahara Y, Furukawa K, Fujii T: Isolation of a multiheme protein with features of a hydrazine-oxidizing enzyme from an anaerobic ammonium-oxidizing enrichment culture. Appl Environ Microbiol. 2007, 73 (4): 1065-1072. 10.1128/AEM.01978-06.PubMed CentralPubMedView Article
- Shimamura M, Nishiyama T, Shinya K, Kawahara Y, Furukawa K, Fujii T: Another multiheme protein, hydroxylamine oxidoreductase, abundantly produced in an anammox bacterium besides the hydrazine-oxidizing enzyme. J Biosci Bioeng. 2008, 105 (3): 243-248. 10.1263/jbb.105.243.PubMedView Article
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.