- Research article
- Open Access
Identification of the type II cytochrome cmaturation pathway in anammox bacteria by comparative genomics
BMC Microbiology volume 13, Article number: 265 (2013)
Anaerobic ammonium oxidizing (anammox) bacteria may contribute up to 50% to the global nitrogen production, and are, thus, key players of the global nitrogen cycle. The molecular mechanism of anammox was recently elucidated and is suggested to proceed through a branched respiratory chain. This chain involves an exceptionally high number of c-type cytochrome proteins which are localized within the anammoxosome, a unique subcellular organelle. During transport into the organelle the c-type cytochrome apoproteins need to be post-translationally processed so that heme groups become covalently attached to them, resulting in mature c-type cytochrome proteins.
In this study, a comparative genome analysis was performed to identify the cytochrome c maturation system employed by anammox bacteria. Our results show that all available anammox genome assemblies contain a complete type II cytochrome c maturation system.
Our working model suggests that this machinery is localized at the anammoxosome membrane which is assumed to be the locus of anammox catabolism. These findings will stimulate further studies in dissecting the molecular and cellular basis of cytochrome c biogenesis in anammox bacteria.
One of the most recent additions to the microbial nitrogen cycle is the anaerobic oxidation of ammonium (anammox), which utilizes nitrite as the electron acceptor and forms dinitrogen gas under anaerobic conditions. Anammox bacteria possess intracellular membrane systems, leading to a remarkable cell compartmentalization . Two membranes on the inner side of the protein-rich cell wall form a ribosome-free peripheral compartment, the paryphoplasm . A third and innermost bilayer membrane exhibits a highly curved configuration and further separates the cytoplasm into two distinct regions, namely the riboplasm and the anammoxosome (Figure 1A). Detailed electron microscopic and labeling studies strongly support the hypothesis of the anammoxosome being a separate organelle, where the central anammox catabolism resides [1, 3, 4]. The annotation of more than 200 genes involved in catabolism and respiration in the genome of the anammox bacterium Kuenenia stuttgartiensis, together with the abundance of 61 genes encoding c-type cytochrome proteins, reflects the complexity of the anammox metabolism and implies the presence of a branched and versatile respiratory chain . This complexity is further confirmed by the genome assemblies of two more anammox species that were recently reported (Scalindua profunda ; strain KSU-1 ). Although c-type cytochrome proteins seem to play a key role in the unique anammox metabolism, the maturation pathway of functional c-type cytochrome holoforms has not been explored. Cytochrome c maturation describes the post-translational process by which b-type hemes (Fe-protoporphyrin IX) are covalently attached to the apoproteins resulting in functional c-type cytochromes. After synthesis, apocytochrome c and heme molecules are independently translocated across the energy-transducing membrane into the bacterial periplasm, the mitochondrial intermembrane space or the thylakoid lumen. Ferric iron of heme(s) and cysteine residues of apocytochrome c are reduced and subsequent thioether linkage formation occurs between the heme vinyl groups and the CX2-4CH sulfhydryls of apocytochrome c, leading to the functional holoform . Three distinct cytochrome c maturation pathways (Systems I, II and III) have been described, each comprising system-specific assembly protein complexes; these biogenesis systems occur in a wide variety of organisms with a complex and unpredictable phylogenetic distribution .
Considering the remarkable anammox cell plan together with the high abundance of cytochrome c-type proteins, determination of the cytochrome c maturation system that these bacteria employ is of particular importance. In this study, comparative computational methods were applied to determine the maturation pathway regulating the assembly of functional c-type cytochrome holoforms in four genera of anammox bacteria, using key protein constituents of maturation Systems I-III as biomarkers. Our analysis showed that all anammox genome assemblies contain at least one full set of System II (Ccs) genes.
All anammox bacteria belong to the order Brocadiales that branches deeply into the phylum Planctomycetes and includes five genera (Kuenenia, Scalindua, Brocadia, Jettenia, and Anammoxoglobus) . In this study draft genomes representative of four anammox genera were analyzed. Kuenenia stuttgartiensis [NCBI bioproject: PRJNA16685 ], Scalindua profunda [JGI: 2017108002 and 2022004002 ], and strain KSU-1 (representing Jettenia genus) [NCBI bioprojects: PRJDA163683 and PRJDB68 ] obtained as described elsewhere. Genomic data for Brocadia fulgida were obtained as described here below.
Brocadia fulgidagenomic data
Library preparation and sequencing
All kits used in this section were obtained from Life technologies (Life technologies, Carlsbad, CA, USA). Genomic DNA, isolated using a CTAB phenol/chloroform based method, was sheared for 5 minutes using the Ion Xpress™ Plus Fragment Library Kit following the manufacturer’s instructions. Further library preparation was performed using the Ion Plus Fragment Library Kit following manufacturer’s instructions. Size selection of the library was performed using an E-gel 2% agarose gel. Emulsion PCR was performed using the Onetouch 200 bp kit and sequencing was performed on an IonTorrent PGM using the Ion PGM 200 bp sequencing kit and an Ion 318 chip, resulting in 5.25 million reads with an average length of 179 bp.
Assembly and annotation
The obtained 5.25 million reads were quality trimmed and all reads below 200 bp were discarded. The remaining 2,22 million reads were assembled using the CLC genomics workbench (v6.5.1, CLCbio, Aarhus, Denmark) with word size 35 and bubble size 5000. Brocadia fulgida accounted for 91% of the assembled reads. Contigs were assigned to Brocadia fulgida based on coverage (>30 fold). The obtained 411 contigs were annotated using Prokka 1.7.2 (Prokka: Prokaryotic Genome Annotation System - http://vicbioinformatics.com/). After annotation, a round of manual curation was performed to correct detected frame shifts. Raw reads and assembled data are available under NCBI bioproject PRJEB4876.
Cytochrome cmaturation pathway
Reference protein datasets for each of the three cytochrome c maturation Systems (I-III) were compiled (Additional file 1), each comprising all protein and polypeptide sequences available at UNIPROT, annotated as any of the defining system-specific components (Additional files 2 and 3). A thioredoxin dataset for maturation System II was also constructed comprising UNIPROT entries for CcsX, DsbD, and CcdA. All abovementioned datasets were limited to peer-reviewed entries.
All anammox gene products were compared to the datasets using blastP (as implemented in the CLC genomics workbench, v6.5.1, CLCbio, Aarhus, Denmark) with an E-value cut off of 10-6. Significant hits were further analyzed by HHpred against all available HMM databases with HHBlits as the MSA generation method . The web server implementation of HMMER (default settings) was also utilized . Protein family matches were identified via Pfam (default settings) . For structure- or sequence-specific feature recognition, transmembrane helical domains were predicted using the TMHMM web server  and potential signal peptides were annotated using SignalP 4.1 . Conserved motifs and critical residues were procured from literature (Additional file 2) and probed in each gene product directly. Multiple alignments of CcsA and CcsB anammox homologs were performed using ClustalW (default settings) and phylogenetic trees were constructed based on the Maximum Likelihood algorithm utilizing the JTT matrix-based model (test of phylogeny: bootstrap method; number of replications: 1000; gaps/missing data treatment: use all sites), both as implemented in MEGA 5.0 . BlastP was also utilized to search for related outgroup sequences in GenBank.
Results & discussion
Assignment of cytochrome cmaturation System II in anammox bacteria
In this study, we applied comparative genomics to predict the maturation pathway of c-type cytochrome proteins in four anammox genera, using key protein components of maturation Systems I-III as biomarkers.
Using our approach, none of the marker genes for System I or III could be identified in the anammox draft genomes. On the contrary, our overall results evinced System II to be the dedicated c-type cytochrome biogenesis pathway that anammox bacteria employ.
System II, (cytochrome c synthesis, 'ccs’) comprises three system-specific proteins (CcsABX) together with a thiol-disulfide membrane transporter (DsbD or CcdA). According to the bacterial working model, two transmembrane proteins (CcsAB), forming a channel entry, facilitate the heme transport and the maintenance of it in a reduced state at the p-side of the membrane . A dedicated membrane-anchored thiol-disulfide oxidoreductase (CcsX) reduces the apocytochrome c cysteines while reducing equivalents are transferred from a non-specific cytoplasmic thioredoxin to the thiol-disulfide membrane transporter (DsbD or CcdA) . Eventually, spontaneous ligation for the thioether linkages formation takes place .
Following the experimental approach described above, homologs of CcsA (sometimes referred to as ResC) were successfully identified in all anammox genera; three putative CcsA proteins were found in Kuenenia, strain KSU-1 and Scalindua and two in Brocadia (Additional file 4). For a functional type II cytochrome c maturation system, complexation of CcsA and CcsB is required . CcsB (sometimes called ResB) exhibits weak sequence conservation although structural homology is observed . Our results further support this, since only one isoform for each Kuenenia, Scalindua, and strain KSU-1 was found by reference database search and two for Brocadia (Additional file 4). Nevertheless, when intra- and intergenome examination with the significant CcsB hit of Kuenenia as query was performed, one more CcsB isoform was retrieved for each Kuenenia, Scalindua and strain KSU-1. Results from HHpred and HMMER annotation were strikingly in agreement with those generated by blastP (compare Additional file 4 with Additional file 5). It is surprising that anammox genera contain multiple CcsB homologs; to the best of our knowledge, only one CcsB homolog has been found in any other organism to date.
Functional assignment of CcsA and CcsB is based on sequence homology , a minimum number of transmembrane helices and the presence of conserved motifs and essential residues (see Additional file 2). The combined results indicate that all anammox genera tested herein share a common protein pattern regarding their cytochrome c maturation system, all coding for two distinct CcsA-CcsB complexes (Table 1). All CcsA and CcsB homologs of Kuenenia and Scalindua were also detected in transcriptome and proteome analyses [6, 20]. In detail, in the genomes of Kuenenia, Brocadia, strain KSU-1 and Scalindua a CcsA homolog, possessing the CcsA-specific tryptophan-rich heme-binding motif (WAXX(A/δ)WGX(F/Y)WXWDXKEXX) and 8 transmembrane helices, is found adjacent to a CcsB homolog possessing 2-4 transmembrane helices and a large soluble domain. Notably, the CcsB sequence motif (VNX1-4P) is found in duplicate in the canonical CcsB from strain KSU-1, whereas in Scalindua only a truncated CcsB motif is retrieved (VN) albeit three times. Intriguingly, the second CcsA-CcsB cytochrome c maturation complex encoded by all four anammox genera displays alterations from the canonical complex  regarding a modified CcsA heme-binding motif:
Published: W A X X (A/S) W G X (F/Y) W X W D X K E X X
Modified: W G X X A W G X Y F L W D A K (V/L) (V/L) W
In the latter, the observed amino acid substitutions may suggest a structurally different heme-binding configuration and/or implications for protein functionality. Nonetheless, the identified CcsA and CcsB homologs are coded adjacent to each other in all anammox genomes. Phylogenetic relationships among the anammox CcsA and CcsB homologs are illustrated in Figure 2A and 2B, respectively.
Along with a functional CcsA-CcsB complex, cytochrome c maturation System II further requires an efficient thiol-reduction pathway through which reducing equivalents are shuttled across the energy-transducing membrane towards the p-side, and are subsequently used for reduction of apocytochrome c cysteines . In all four anammox species we studied, DsbD, a thiol-disulfide membrane transporter involved in the aforementioned pathway, is annotated successfully and with high confidence by a similar comparative methodology adopted for CcsA and CcsB (Table 2 and Additional file 6). In detail, two DsbD homologs are identified in Kuenenia whereas a single copy is retrieved for strain KSU-1 and Brocadia. All DsbD homologs share similar structural features, including 8-11 transmembrane helices and conserved cysteine residues . Scalindua contains a homolog of CcdA, related to but shorter than DsbD, possessing only 6 transmembrane helices along with two cysteine residues . DsbD is a housekeeping thiol-disulphide electron shuttle  and as such it is not an indispensable cytochrome c maturation System II component. In contrast, CcsX (sometimes called ResA) that fulfils the essential role of apocytochrome c reduction in this disulfide bond cascade is a dedicated membrane-anchored thiol-disulfide oxidoreductase of maturation System II. Apart from the conserved thioredoxin cytochrome c recognition motif (CXXC), CcsX also possesses additional cysteine residues and a single transmembrane helix through which it is anchored to the membrane. Our comparative computational approach identified multiple potential CcsX homologs for each anammox genus. Particularly, two CcsX-like homologs for Brocadia, three for Kuenenia and six for each, Scalindua and strain KSU-1, were identified with high confidence (Additional file 6). However, homologs possessing no signal peptide sequences were ruled out from our final collective table (Table 2). Although distinction between the dedicated CcsX proteins and other thioredoxins that might possess similar features cannot be made, the presence of that many CcsX-like homologs suffices for a complete c-type cytochrome maturation System II.
Overall, these results indicate that the assembly of cytochrome c holoforms is achieved by the maturation System II in all anammox bacteria tested herein. All genera code for at least one CcsA-CcsB complex, one DsbD (or CcdA), and one CcsX homolog, all being essential components of a functional cytochrome c maturation System II.
Having analyzed the cytochrome c maturation system in anammox bacteria, it would be stimulating to comprehend how such machinery is localized within the intricate anammox cell plan. A hypothetical cellular pathway for cytochrome c biogenesis is illustrated in Figure 1B. According to our view, the CcsA-CcsB complex, forming the heme channel entry, must be tethered within the anammoxosome membrane. Heme is, thus, translocated into the anammoxosome, with the latter representing the p-side of the anammox cell . This translocation is mediated by selective CcsA heme-binding motifs (as specified in Table 1). Concurrently, housekeeping riboplasmic thioredoxins provide DsbD with the necessary reductants that are shuttled towards the dedicated CcsX thiol-disulfide oxidoreductase. Both DsbD and CcsX possess transmembrane helices spanning the anammoxosome membrane, with the CcsX globular domain facing the inside of the anammoxosome, where apocytochrome c cysteine reduction occurs. Eventually, spontaneous formation of the thioether linkages between the apoprotein and its cofactor takes place, leading to functional cytochrome c holoforms inside the anammoxosome .
These findings suggest that anammox bacteria possess at least one complete machinery for type II cytochrome c biogenesis , adapting it to their complicated cell plan; the anammoxosome membrane is proposed to be the main site of cytochrome c maturation. Our results provide a working model that will be used to guide experimental studies, including protein purification and immunogold electron microscopy, in elucidating both the localization and the function of cytochrome c maturation System II in anammox bacteria.
The data sets supporting the results of this article are included within the article and its additional files.
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. Arch Microbiol. 2001, 175: 413-429. 10.1007/s002030100280.
Jetten MSM, Niftrik LV, Strous M, Kartal B, Keltjens JT, Op den Camp HJM: Biochemistry and molecular biology of anammox bacteria. Critic Rev Biochem Mol Biol. 2009, 44: 65-84.
van Niftrik L, van Helden M, Kirchen S, van Donselaar EG, Harhangi HR, Webb RI, Fuerst J, Op den Camp HJM, Jetten MSM, Strous M: Intracellular localization of membrane-bound ATPases in the compartmentalized anammox bacterium 'Candidatus Kuenenia stuttgartiensis’. Mol Microbiol. 2010, 77: 701-715. 10.1111/j.1365-2958.2010.07242.x.
van Niftrik L, Geerts WJC, van Donselaar EG, Humbel BM, Webb RI, Fuerst J, Verkleij AJ, Jetten MSM, Strous M: Linking ultrastructure and function in four genera of anaerobic ammonium-oxidizing bacteria: cell plan, glycogen storage, and localization of cytochrome c proteins. J Bacteriol. 2008, 190: 708-717. 10.1128/JB.01449-07.
Strous M, Pelletier E, Mangenot S, Rattei T, Lehner A, Taylor MW, Horn M, Daims H, Bartol-Mavel D, Wincker P, Barbe V, Fonknechten N, Vallenet D, Segurens B, Schenowitz-Truong C, Médigue C, Collingro A, Snel B, Dutilh BE, Op den Camp HJM, van der Drift C, Cirpus I, van de Pas-Schoonen KT, Harhangi HR, van Niftrik L, Schmid M, Keltjens J, van de Vossenberg J, Kartal B, Meier H, et al: Deciphering the evolution and metabolism of an anammox bacterium from a community genome. Nature. 2006, 440: 790-794. 10.1038/nature04647.
van de Vossenberg J, Woebken D, Maalcke WJ, Wessels HJ, Dutilh BE, Kartal B, Janssen-Megens EM, Roeselers G, Yan J, Speth D, Gloerich J, Geerts W, van der Biezen E, Pluk W, Francoijs KJ, Russ L, Lam P, Malfatti SA, Tringe SG, Haaijer SC, Op den Camp HJ, Stunnenberg HG, Amann R, Kuypers MM, Jetten MS: The metagenome of the marine anammox bacterium 'Candidatus Scalindua profunda’ illustrates the versatility of this globally important nitrogen cycle bacterium. Environ Microbiol. 2013, 15 (5): 1275-1289. 10.1111/j.1462-2920.2012.02774.x.
Hira D, Toh H, Migita CT, Okubo H, Nishiyama T, Hattori M, Furukawa K, Fujii T: Anammox organism KSU-1 expresses a NirK-type copper-containing nitrite reductase instead of a NirS-type with cytochrome cd1. FEBS Lett. 2012, 586: 1658-1663. 10.1016/j.febslet.2012.04.041.
Hamel P, Corvest V, Giege P, Bonnard G: Biochemical requirements for the maturation of mitochondrial c-type cytochromes. Biochim Biophys Acta. 2009, 1793: 125-138. 10.1016/j.bbamcr.2008.06.017.
Allen J, Ginger M, Ferguson S: Complexity and diversity in c-type cytochrome biogenesis systems. Biochem Soc Trans. 2005, 33: 145-146. 10.1042/BST0330145.
Jetten MSM, Op den Camp HJM, Kuenen JG, Strous M: Description of the order Brocadiales. Bergey’s manual of systematic bacteriology. Volume 4. Edited by: Krieg NR, Staley JT, Brown DR, Hedlund BP, Paster BJ, Ward NL, Ludwig W, Whitman WB. 2010, Heidelberg, Germany: Springer, 596-603.
Soding J, Biegert A, Lupas AN: The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33: W244-W248. 10.1093/nar/gki408.
Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39: W29-W37. 10.1093/nar/gkr367.
Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, Eddy SR, Sonnhammer EL, Bateman A: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, 34: D247-D251. 10.1093/nar/gkj149.
Sonnhammer E, Von Heijne G, Krogh A: A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998, 6: 175-182.
Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786. 10.1038/nmeth.1701.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Frawley ER, Kranz RG: CcsBA is a cytochrome c synthetase that also functions in heme transport. Proc Natl Acad Sci U S A. 2009, 106: 10201-10206. 10.1073/pnas.0903132106.
Beckett CS, Loughman JA, Karberg KA, Donato GM, Goldman WE, Kranz RG: Four genes are required for the system II cytochrome c biogenesis pathway in Bordetella pertussis, a unique bacterial model. Mol Microbiol. 2000, 38: 465-481. 10.1046/j.1365-2958.2000.02174.x.
Kranz RG, Richard-Fogal C, Taylor JS, Frawley ER: Cytochrome c biogenesis: mechanisms for covalent modifications and trafficking of heme and for heme-iron redox control. Microbiol Mol Biol Rev. 2009, 73: 510-528. 10.1128/MMBR.00001-09.
Kartal B, Maalcke WJ, de Almeida NM, Cirpus I, Gloerich J, Geerts W, Op den Camp HJ, Harhangi HR, Janssen-Megens EM, Francoijs KJ, Stunnenberg HG, Keltjens JT, Jetten MS, Strous M: Molecular mechanism of anaerobic ammonium oxidation. Nature. 2011, 479: 127-130. 10.1038/nature10453.
Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comp Appl Biosci: CABIOS. 1992, 8: 275-282.
Stewart EJ, Katzen F, Beckwith J: Six conserved cysteines of the membrane protein DsbD are required for the transfer of electrons from the cytoplasm to the periplasm of Escherichia coli. EMBO J. 1999, 18: 5963-5971. 10.1093/emboj/18.21.5963.
Porat A, Cho SH, Beckwith J: The unusual transmembrane electron transporter DsbD and its homologues: a bacterial family of disulfide reductases. Res Microbiol. 2004, 155: 617-622. 10.1016/j.resmic.2004.05.005.
Ito K, Inaba K: The disulfide bond formation (Dsb) system. Curr Opin Struct Biol. 2008, 18: 450-458. 10.1016/j.sbi.2008.02.002.
The authors thank Boran Kartal and Katinka van de Pas-Schoonen for the enrichment cultures of Brocadia fulgida. Daan R. Speth was supported by BE-Basic fp0702, Joachim Reimann was supported by NWO 62001417, James W. A. Allen was supported by BBSRC, and Mike S. M. Jettten and Christina Ferousi were supported by ERC AG 232937 and Spinoza Premium 2012.
The authors declare that they have no competing interests.
CF, JWAA and MSMJ conceived of the study. DRS sequenced and analyzed the genomic data of Brocadia. CF built the datasets and ran homologue searches. DRS, JR, JTMK, and HJMOC assisted in bioinformatics analysis and data interpretation. CF, JWAA, and MSMJ wrote the manuscript with input from all co-authors. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: c maturation Systems (I-III) and thioredoxin dataset for System II.(ZIP 301 KB)
Additional file 2: c maturation System biomarkers. For each cytochrome c maturation System (I-III), essential protein components that can be used as suitable biomarkers for annotation purposes were selected (for details see Additional file 3) and their defining characteristics are listed herein. (XLSX 12 KB)
Additional file 3: c maturation System biomarkers.(PDF 35 KB)
Additional file 4: Homology identification was performed with blastP as implemented in CLC genomics workbench (v6.5.1, CLCbio, Aarhus, Denmark). Whole anammox genomes are used as queries against a reference database that comprises all reviewed entries for CcsA and CcsB available at UNIPROT. An E-value of 10-6 was set as cut off to prevent ambiguity. (XLSX 14 KB)
Additional file 5: Homology identification was performed with blastP as implemented in CLC genomics workbench (v6.5.1, CLCbio, Aarhus, Denmark). Whole anammox genomes are used as queries against a reference database that comprises all reviewed entries for CcsA and CcsB available at UNIPROT. Intra- and intergenome searches with the significant hits from Kuenenia as queries were also performed (Additional file 4). Retrieved results were further analyzed with HHpred and HMMER. An E-value of 10-3 was set as cut off to prevent ambiguity. (XLSX 14 KB)
Additional file 6: Homology identification was performed with blastP as implemented in CLC genomics workbench (v6.5.1, CLCbio, Aarhus, Denmark). Whole anammox genomes are used as queries against a reference database that comprises all reviewed entries for CcsX and DsbD available at UNIPROT. Retrieved results were further analyzed with HHpred and HMMER. (*): E-value cut off set at 10-6; (**): E-value cut off set at 10-3. (XLSX 14 KB)
About this article
Cite this article
Ferousi, C., Speth, D.R., Reimann, J. et al. Identification of the type II cytochrome cmaturation pathway in anammox bacteria by comparative genomics. BMC Microbiol 13, 265 (2013). https://doi.org/10.1186/1471-2180-13-265
- Cytochrome c biogenesis
- Anaerobic ammonium oxidizing bacteria