Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2009, p. 4835-4852, Vol. 75, No. 14
0099-2240/09/$08.00+0 doi:10.1128/AEM.02874-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Timothy J. Williams,2
Gary Xie,3,4
Hajnalka E. Kiss,3,4
Thomas S. Brettin,3,4
Diego Martinez,3,4
Christian A. Ross,1
Dirk Schüler,5
B. Lea Cox,6
Kenneth H. Nealson,6 and
Dennis A. Bazylinski1*
University of Nevada, Las Vegas, School of Life Sciences, Las Vegas, Nevada 89154-4004,1 School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW 2052, Australia,2 Bioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico,3 DOE Joint Genome Institute, Department of Energy, Los Alamos, New Mexico,4 Department of Biology I, Ludwig Maximilians University Munich, 80638 Munich, Germany,5 Department of Earth Sciences, University of Southern California, 3616 Trousdale Parkway, Los Angeles, California 90089-03716
Received 17 December 2008/ Accepted 15 May 2009
|
|
|---|
|
|
|---|
Cells of strain MC-1 are not metabolically versatile, based on previous growth experiments (14, 113) and genomic information presented here. Cells of this species grow chemolithoautotrophically, with thiosulfate or sulfide as an electron donor, and chemoheterotrophically, with acetate as the electron donor and carbon source (14). Unlike some other autotrophic marine MTB, such as vibrio strain MV-1, which represents a new genus of bacteria in the Alphaproteobacteria ("Candidatus Magnetovibrio blakemorei" [D. A. Bazylinski, unpublished data]), cells of strain MC-1 do not utilize the Calvin-Benson-Bassham cycle for CO2 fixation but rely on the reverse tricarboxylic acid (rTCA) cycle (113).
The salient feature of MTB is their ability to biomineralize magnetosomes. These intracellular structures consist of well-ordered crystals of magnetite (Fe3O4) or greigite (Fe3S4) surrounded by a lipid bilayer membrane (5, 42, 50, 76). Different strains of MTB synthesize magnetosome crystals not only of a specific composition but also of a specific morphology and size (8, 13). These consistent, species-specific crystal morphologies and sizes indicate that biomineralization of magnetite is genetically controlled in MTB, and it is likely that a common set of genes is responsible for magnetosome biomineralization, although there are presumably some species-specific genes that might be involved in some aspects of biomineralization, such as the control of crystal shape and/or size.
In the last decade, genomic and genetic studies have confirmed that many cultured MTB contain a common set of genes that appear to be involved in magnetite biomineralization (89). It has been shown that the genes encoding magnetosome (membrane) proteins are localized on a 130-kb magnetosome island (MAI) in Magnetospirillum gryphiswaldense (94, 111). In addition, this genomic island is conserved among other Magnetospirillum strains, including Magnetospirillum magneticum strain AMB-1 and Magnetospirillum magnetotacticum strain MS-1 (48, 61). This may not be surprising since these strains are closely related. Recent studies show that the magnetotactic marine vibrio strain MV-1 contains an MAI as well. Comparisons to other MTB strains revealed that strain MV-1 contains a highly conserved mamAB cluster, but other magnetosome gene clusters are arranged differently or are absent (60). A comparative genomic analysis between these Magnetospirillum species and the distantly related strain MC-1 revealed some significant differences in the set of magnetosome genes. Richter et al. (89) compared the gene contents based on reciprocal best matches and showed that strain MC-1 shares only about one-half the number of magnetosome genes with M. gryphiswaldense and M. magneticum that the Magnetospirillum strains share with each other. This study also revealed a "Magnetospirillum-specific" set of genes as well as some MTB-specific genes that are also present in strain MC-1 (89). There is evidence that the MC-1 genome is mosaic to the highest degree found within the alphaproteobacteria. Only 33% of the genes showed a best BLAST hit to those in a collection of 18 alphaproteobacterial sequences (39). Earlier studies showed that strain MC-1 is one of the more divergent alphaproteobacteria (117).
In this study, we analyzed the complete genome sequence of the marine magnetotactic coccus strain MC-1. This analysis provided us with considerable new insights into the organization of the magnetosome genes in MC-1 and the MAI in which they reside, as well as allowing us to compare the magnetosome genes to those of other MTB strains to determine whether novel magnetosome genes are present in MC-1 and not in other MTB strains. In addition, we sought to examine whether the genome sequence (through the presence and absence of specific genes) supports the known physiological and metabolic aspects of MC-1 as well as its lack of metabolic versatility and if the sequence could reveal or suggest novel metabolic or physiological features of this organism. It has clearly been shown that growth and environmental conditions that affect physiology also greatly influence magnetosome production in some species (21, 97). The results from this study might prove invaluable in the understanding of environmental control of and physiological effects on magnetosome synthesis as well as for the isolation of other MTB, in particular the ubiquitous magnetotactic cocci, in terms of their use of electron donors and acceptors.
|
|
|---|
Genomic DNA from strain MC-1 was isolated as previously described by Kimble et al. (64). The quality of the DNA was determined by agarose gel electrophoresis.
Genome sequencing and assembly.
The MC-1 genome was sequenced using a combination of 3-kb, 8-kb, and 40-kb (fosmid) DNA libraries. The number of clones made for these libraries was enough for about 15x coverage of the genome, using the three libraries. All general aspects of library construction and sequencing were performed at the Joint Genome Institute (JGI), and details can be found at http://www.jgi.doe.gov/.
Draft assemblies were based on 126,624 total reads. The Phred/Phrap/Consed software package was used for sequence assembly and quality assessment (40, 41, 51). After the shotgun stage, the 126,624 reads were assembled with parallel Phrap (High Performance Software, LLC). Possible misassemblies were corrected with Dupfinisher (55) or transposon bombing of bridging clones (Epicentre Biotechnologies, Madison, WI). Gaps between contigs were closed by editing in Consed, custom primer walking, or PCR amplification (Roche Applied Science, Indianapolis, IN). A total of 27,823 additional sequencing reactions were necessary to close gaps and to raise the quality of the finished sequence. The final assembly of the genome of the magnetic coccus strain MC-1 contained 107,421 reads after removal of the majority of the draft reads that landed in duplicated regions during the finishing process. The error rate of the final genome sequence is less than 1 in 100,000 bp.
Automated gene prediction was performed by using the output of Critica (4) complemented with the output of Generation and Glimmer (32). The tRNAScanSE tool (72) was used to find tRNA genes, whereas rRNAs were found by using BLASTN searches of 16S and 23S rRNA databases. Other "standard" structural RNAs (e.g., 5S rRNA, rnpB, transfer messenger RNA, and signal recognition particle RNA) were found by using covariance models with the internal search tool (37). The assignment of product descriptions was made by using search results for the following curated databases, in the indicated order: TIGRFam, PRIAM (e–30 cutoff), Pfam, Smart, COGs, Swissprot/TrEMBL (SPTR), and KEGG. If there was no significant similarity to any protein in another organism, it was described as a "hypothetical protein." The term "conserved hypothetical protein" was used if at least one match was found to a hypothetical protein in another organism. EC numbering was based on searches in PRIAM with a cutoff of e–10; COG and KEGG functional classifications were based on homology searches in the respective databases.
Comparative analysis of strain MC-1 and related organisms was performed by using a set of tools available at the JGI Integrated Microbial Genomes (IMG) website (http://img.jgi.doe.gov/cgi-bin/pub/main.cgi). Unique and orthologous MC-1 genes were identified by using BLASTP (using cutoffs of E scores of <10–2 and 20% identity and reciprocal hits with cutoffs of E scores of <10–5 and 30% identity, respectively). Signal peptides and transmembrane helices were predicted using SignalP 3.0 (16) and TMHMM 2.0 (68), set at default values. Protein localizations were predicted with PSORTb (49), and twin-arginine translocation systems were identified using the TatP program (17). Insertion sequence elements were identified by using the ISFinder database (102). Metabolic pathways were constructed using MetaCyc as a reference data set (25). Other details regarding genome properties and genome annotation can be obtained from the JGI's IMG website (77).
Prophage prediction.
The prediction of potential prophage loci was done with Prophage Finder under default conditions, using an E value of 0.5, five hits per prophage, and 5,500-bp hit spacing (22; http://bioinformatics.uwp.edu/
phage/ProphageFinder.php).
Nucleotide sequence accession number.
The sequence of strain MC-1 can be accessed using GenBank accession number CP000471.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. General genome features of strain MC-1a
|
![]() View larger version (57K): [in a new window] |
FIG. 1. Circular representation of the genome of the magnetotactic coccus strain MC-1. The location of the MAI on the chromosome is indicated in red outside the circles. The outer two rings (ring 1 and ring 2) depict predicted protein-encoding and structural RNA genes on the plus and minus strands, respectively, and are color coded according to COG category. The black circle (ring 3) indicates the deviation from the average %GC, and the purple and green circle (ring 4) represents the GC skew.
|
Magnetosome genes.
The genomes of three strains of MTB, all in the genus Magnetospirillum (M. magnetotacticum strain MS-1, M. gryphiswaldense strain MSR-1, and M. magneticum strain AMB-1), have been sequenced at least partially, revealing an MAI containing the regulatory and structural genes coding for magnetosome (membrane) proteins. The genome of MC-1 reveals a putative MAI of about 102 kb, located on the chromosome between 2.80 and 2.93 Mb (Fig. 1). Strong evidence from tetranucleotide usage patterns in Magnetospirillum species and strain MV-1 showing that the MAIs in these MTB are transferred by horizontal gene transfer was recently described (60). However, the same analysis showed that the putative MAI in strain MC-1 lacks a conspicuous tetranucleotide signal (60). The GC bias of the MAI in strain MC-1 shows the same random distribution (Fig. 1). These observations might result from a very early transfer of the island to MC-1-like organisms or from additional horizontal gene transfer events after the transfer of the MAI. These possibilities might also explain the mosaic structure of the MC-1 genome (39). Nonetheless, the presence of an MAI in MC-1 is based on the very close proximity of the majority of the magnetosome genes in a relatively small genomic region and the presence of two integrase encoding genes, IS4 and IS3, that border this genomic region (Fig. 2). Within the putative MAI, the genes are arranged in four clusters (mamAB, mamCEIH, mamDXZ, and mms), with a fifth magnetosome-related cluster (mtx) situated outside the MAI. In contrast to many hypervariable genomic islands (18, 52, 94), the MC-1 putative MAI contains no transposon genes within the island, and thus the position of the MAI within the genome is likely to be more stable than those in other MTB. However, some characteristic features of genomic islands are in close proximity to the island; for example, two genes with homology to insertion sequence elements of the IS4 family (35) and three tRNA genes that often serve as attachment sites for phages are located upstream of the MAI. In addition, 14 prophages were identified distributed throughout the genome. This is in contrast to the MAIs in Magnetospirillum gryphiswaldense and M. magneticum, which contain numerous integrase and transposase genes, presumably involved in the mobility of the MAI, that make up more than 20% of the coding sequence of the MAI. Spontaneous deletions and mutations resulting from these integrase and transposase genes have been shown to cause nonmagnetic mutants to occur at a rate of 10–2 in M. gryphiswaldense and M. magneticum (48, 111). Nonmagnetic mutants have never been observed in strain MC-1. Another key difference in the MC-1 genome compared to the genomes of the Magnetospirillum strains is the organization of the magnetosome genes within the MAI (Fig. 2).
![]() View larger version (26K): [in a new window] |
FIG. 2. Genomic organization of the magnetosome genes in MC-1. The magnetosome genes are organized into five different clusters localized on a putative 102-kb MAI. Green, genes of the mamAB cluster; purple, genes of the mms cluster; dark blue, genes of the mamDXZ cluster; orange, genes of the mamCEIH cluster; light blue, putative new magnetosome genes; yellow, tRNA genes; red, mobility genes; gray, hypothetical genes; white, housekeeping genes.
|
|
View this table: [in a new window] |
TABLE 2. Features of deduced genes in the 102-kb MAI in strain MC-1
|
(ii) mamCEIH cluster.
The genes mamC, mamD, and mamF are unique to MTB and, along with mamG, form a single cluster in Magnetospirillum strains (89). There is no mamGFDC cluster in MC-1; instead, mamC, mamD, and mamF are scattered throughout the MAI, and mamG is absent from the genome (this study). In MC-1, mamC belongs to a gene cluster 12.4 kb downstream of the mamAB cluster, along with mamE, mamI, and mamH as well as six additional ORFs (indicated in light blue in Fig. 2). Thus, strain MC-1 appears to contain a new magnetosome gene cluster, here termed the mamCEIH cluster (Mmc12263 to Mmc12272), that is transcribed divergently as well. MamH belongs to the major facilitator superfamily, a family of transporters capable of transporting small solutes (89). MamE contains a PDZ domain and shows sequence similarity to trypsin-like (HtrA-like) serine proteases (53). Also within the MC-1 mamCEIH cluster is a gene that encodes a protein with a ferritin-like domain (Mmc12264), a gene (Mmc12263) that encodes a hypothetical protein, and a gene that encodes a Na+/Ala symporter (Mmc12272); none of these are found in Magnetospirillum strains. mamD belongs to the mamDXZ cluster (see below) and also shows weak similarity to Mmc12270 in the mamCEIH cluster. mamF is present as three different copies in the MC-1 genome: a mamF-like gene (Mmc12258) lies within the mamAB cluster, mamF (Mmc12234) is adjacent to the mamDXZ cluster, and mmsF (Mmc12274) is located within the mms6 cluster (see below).
(iii) mamDXZ cluster.
The mamDXZ cluster (Mmc12237 to Mmc12240), a cluster containing four genes, was formerly designated the mamXY cluster, which is conserved in Magnetospirillum strains (89). In strain MC-1, mamY is replaced by mamD, and we could not identify the mamY gene or any homologous genes in the MC-1 genome. This cluster is located 15.4 kb upstream of the mamAB cluster, transcribed in the right orientation. The MamD protein appears to be involved in the control of the size of the magnetosome magnetite crystals in M. gryphiswaldense (93). This cluster also contains another mamD-like gene within 5 kb of mamD, mamX (similar to mamE and mamS), and mamZ (a mamH-like gene) (see above). The genes in the mamXY cluster of Magnetospirillum strains encode an FtsZ-like protein, a MamH-like protein, MamX, and MamY (89). A gene encoding an FtsZ-like protein was not identified in the MAI of MC-1, though one is present elsewhere in the genome (Mmc10747).
(iv) mms genes.
The genomes of both M. gryphiswaldense and M. magneticum contain an mms gene cluster consisting of five genes (48, 111). In M. magneticum, mms6 encodes a magnetosome membrane-associated protein implicated in magnetite crystal nucleation (3). In MC-1, the mms cluster is located 589 bp downstream of the mamCEIH cluster, transcribed in the opposite direction, and contains two genes: one is homologous to mms6 (Mmc12275), and the other is homologous to mmsF (Mmc12274). The remaining three genes of the mms cluster in Magnetospirillum could not be identified in the genome of strain MC-1.
(v) mtx cluster.
A gene cluster referred to as the mtx cluster was found to be conserved among different MTB strains, including MC-1, but is not located within the putative MAI in any MTB (89). In MC-1, there is an mtx cluster (Mmc13695 to Mmc13697) located about 1.7 Mb downstream of the MAI. The gene cluster consists of a TPR protein gene (Mmc13695); the mtxA gene (Mmc13696), which encodes a magnetosome membrane protein conserved among different MTB strains; and a gene encoding an adenylate cyclase family protein presumed to be involved in signal transduction.
(vi) Putative new magnetosome genes in strain MC-1.
Within the MAI in strain MC-1, 14 previously unrecognized putative genes (ORFs) cluster together with the mam and mms genes. Because most genes within the MAI encode magnetosome membrane proteins in strain MC-1 and/or other MTB, it seems likely that these new genes encode such proteins and therefore might be involved in magnetosome biomineralization in strain MC-1. Between the mamAB and mamDXZ cluster is a cluster of five genes (Mmc12241 to Mmc12245) which are transcribed in the opposite direction. The mamCEIH cluster contains a total of 10 genes arranged in the same orientation for transcription. Aside from the four genes which give the cluster its name, there are six other genes, including a gene that encodes a protein with a ferritin-like domain (Mmc12264) and a MamD-like protein gene (Mmc12270) (Table 2). The presence of genes in the MAI of MC-1 that are absent from Magnetospirillum strains and/or are unique to MC-1 suggests that at least some of these genes are specific to the magnetosome biomineralization process in strain MC-1 and, for example, may be involved in synthesis of elongated magnetite crystals (versus the equidimensional crystals synthesized by Magnetospirillum strains) (13). This may be especially true of genes that encode proteins such as TPR proteins and the Mmc12260 protein (see above), which show motifs consistent with a role in multiprotein assembly and/or cytoskeletal structures. Conversely, genes for certain proteins that have been ascribed an important (and even essential) role in magnetosome chain formation and/or are associated with the magnetosome membranes of other MTB could not be identified in the MC-1 genome. These include the gene for MamJ (discussed earlier) and the gene for MamW. The latter was identified as a magnetosome membrane protein in M. gryphiswaldense by Ullrich et al. (111) and is encoded by a gene found within the MAI but outside previously identified gene clusters. The absence of such genes from the MC-1 genome further supports the hypothesis that the processes of magnetosome biomineralization and chain assembly in MC-1 may exhibit marked differences from those in Magnetospirillum strains, possibly resulting in magnetosome crystal size and morphology differences.
Iron metabolism.
Although magnetosome biomineralization was discussed in general in a previous section, we focus here on iron uptake and its oxidation and reduction, processes which must be important in magnetosome biomineralization in strain MC-1 and other MTB. In the cell, iron can be oxidized and reduced either chemically (without the use of enzymes) or enzymatically, through respiration and/or metabolism. Important processes involved in the biomineralization of magnetite by MTB include iron uptake by the cell, transport into the magnetosome vesicles (unless the vesicles are permanent invaginations of the cell membrane, as suggested for M. magneticum [65, 66]), and the controlled precipitation of magnetite.
(i) Iron uptake into the cell.
Iron transport is obviously important in MTB, as it is known that up to 3% of the dry weight of cells of these organisms consists of iron (20), which is several orders of magnitude higher than the level in nonmagnetotactic bacteria (12). Even though MTB have been shown to be especially facile with regard to iron uptake, enabling the synthesis of up to 16.7 mg magnetite liter–1 day–1 (108), currently recognized iron uptake systems in MTB seem not to be unique. A comparison of the genomes of four different MTB strains revealed a group of MTB-specific genes that included some magnetosome membrane protein genes but no unique iron or metal transporter genes (89). Surprisingly, only a few specific metal ion transporter genes were identified in the MC-1 genome. There are four genes encoding CDF transporters, including mamB (Mmc12250), mamM (Mmc12256), Mmc13400, and Mmc13150, two ferrous iron transport-encoding genes (feoA and feoB) (Mmc11730 and Mmc11729), and three genes encoding outer membrane receptors, two for ferrienterochelin (Mmc11596 and Mmc12372) and one for iron (Mmc13560). Two Fe(III) uptake regulator-encoding genes (Mmc10894 and Mmc13182) belonging to the Fur (ferric uptake regulator) family were identified as well. Interestingly, although 52 of the 115 transporter-encoding genes present in the genome of MC-1 are part of putative ATP-binding cassette (ABC) systems, none of them appear to be specific for iron.
The genomes of other MTB contain more iron transporter genes than the genome of strain MC-1. For example, M. magneticum contains 10 iron transport-encoding genes. The putative iron transporters are divided into five CDF transporters, two copies of both feoA and feoB, and an additional Fe(II) transporter. The genome of M. magnetotacticum contains 87 genes involved in iron transport, including 27 specific for Fe(III)-siderophore transport systems. In contrast to the MTB, Nitrosomonas europaea has over 100 genes involved in iron transport, including genes for receptors for Fe(III), but it is apparently incapable of synthesizing siderophores (26). It is not known whether cells of strain MC-1 are able to synthesize siderophores. Because MTB take up roughly 100 times more iron than other bacteria (12), we expected to identify more metal transporter-encoding genes, especially specific iron transporter-encoding genes, in the genomes of MC-1 and other MTB. Our analysis shows the opposite trend, with nonmagnetotactic bacteria exhibiting an equal number of, and in some cases more, iron transporters.
(ii) Iron redox reactions.
We could not identify any Fe(II) oxidase genes in the genome of strain MC-1. The copper-containing, periplasmic Fe(II) oxidase of strain MV-1, which may be partly responsible for Fe(II) uptake and processing in this strain (36), has no similarity to any protein predicted by the genome of MC-1. However, two putative Fe(III) reductase-encoding genes (Mmc12239 and Mmc12484) were identified (Fig. 3). The Mmc12239 protein was earlier referred to as MamZ and may be involved in magnetosome biomineralization, as discussed above. The amino acid sequence of the Mmc12484 protein shows high sequence similarity to a putative Fe(III) reductase of other MTB (e.g., M. magnetotacticum [E value, 8e–33], M. magneticum [E value, 3e–32], and M. gryphiswaldense [E value, 2e–25]). For M. magnetotacticum, it was suggested that a cytoplasmic Fe(III) reductase was involved in the biomineralization process (119). Therefore, putative Fe(III) reductase genes, such as that for MamZ (Mmc12239) and Mmc12484, one of which seems to be conserved in these MTB strains, might be important in magnetite biomineralization. We also identified an Fe-S oxidoreductase-encoding gene (Mmc10050) that may also be involved in magnetite biomineralization, as it shows very high similarity to Fe-S oxidoreductases from both M. magneticum and M. gryphiswaldense (E value, 0.0) and a lower similarity to that of M. magnetotacticum, with an E value of 4e–146. For M. magnetotacticum, there is evidence suggesting that an Fe(II):nitrite oxidoreductase is involved in cytoplasmic Fe(II) oxidation, which appears to be essential for magnetosome biomineralization (81).
![]() View larger version (50K): [in a new window] |
FIG. 3. Cell model of strain MC-1, with an emphasis on ultrastructure, transport, energy, carbon metabolism, and chemotaxis. Genes encoding virtually all of the steps for the synthesis of nucleotides and amino acids by canonical pathways are present in the bacterium but are omitted here for simplicity. Electron transport components include NADH dehydrogenase (NDH), periplasmic nitrate reductase (Nap), nitric oxide reductase (cNor), and the sox system (SOX). MCPs are methyl-accepting chemotaxis proteins. CheABRYW form a chemotaxis protein complex (see text). Influx and efflux transporter families with representatives in this genome are indicated on the figure, with the number of genes encoding each type of transporter shown in parentheses. Abbreviations: ABC, ATP binding cassette superfamily; ATP syn, ATP synthetase; CES, cation efflux system protein; CDF, cation diffusion facilitator; FeoA and -B, ferrous iron uptake family proteins; HME, heavy metal efflux pump CzcA family; MFS, major facilitator superfamily; MIT, CorA metal ion transporter family; OME, outer membrane efflux protein; Pho, regulon for phosphate uptake; RND, resistance-nodulation-cell division superfamily.
|
Because we identified very few of the known iron uptake proteins, iron transporters, and iron reductases and no iron oxidases in MC-1, it is possible that MC-1 contains unknown genes that encode novel proteins involved in these processes in the cell.
Transport.
The genome of strain MC-1 contains 115 genes predicted to encode transporters (Fig. 3). This is a relatively small number compared to those in other bacteria, including other MTB. For example, 259 and 694 transporter-encoding genes were identified in the genomes of M. magneticum and M. magnetotacticum, respectively (see Table S1 in the supplemental material). The genomes of some nonmagnetotactic bacteria that are not coccoid in morphology appear to contain many more of these genes. For example, the genome of Escherichia coli contains 405 transporter-encoding genes, and those of Rhodospirillum rubrum and Roseobacter denitrificans contain 241 and 395, respectively. Nonmagnetotactic cocci contain fewer transporters; only 50 to 200 transporter-encoding genes have been found in the genomes of these organisms (see Table S1 in the supplemental material). The small number of transporter genes in MC-1 is also comparable with those in other (albeit obligate) autotrophic proteobacteria and cyanobacteria as well as intracellular pathogenic bacteria (82, 86). Heterotrophic bacteria are known to have more transporter-encoding genes than autotrophs, and these likely function in the uptake of multiple organic carbon and energy sources (e.g., sugars, amino acids, and organic acids) (86, 100). Alternatively, it is possible that MC-1 has transporters which have broad specificity and transport a number of different solutes. In general, prokaryotes vary widely with regard to the number of transporters (see Table S1 in the supplemental material), and the small number of transporters in MC-1 seems to be in keeping with its very limited heterotrophic abilities and its autotrophic lifestyle.
The number of transporters in bacteria has been reported to correlate with the genome size (86), but in a comparison of the numbers of transporter-encoding genes in the genomes of a number of different bacteria relative to the genome size (per Mb) (see Table S1 in the supplemental material), we found a large amount of variation. Values ranged from 16.3 (Mycobacterium leprae) to 161.4 (M. magnetotacticum), and the average for 30 species was 58.0 per Mb. Again, strain MC-1 was at the low end, having 24.5 transporter-encoding genes per Mb.
Genes for metal transporters other than iron transporters (discussed in an earlier section) in the MC-1 genome include those for magnesium (Mmc13074) and zinc (Mmc12530), as well as ABC transport systems for cobalt and molybdate/molybdenum (discussed in the next paragraph). Also notable is the presence of 10 resistance-nodulation-cell division (RND) superfamily genes (Mmc10309, Mmc10980, Mmc11818, Mmc11894, Mmc13248, Mmc12511, Mmc12524, Mmc13162, Mmc13187, and Mmc12361). In Sulfurimonas denitrificans, these genes are predicted to encode transporters involved in metal efflux (101). In MC-1, the presence of these genes could be important in magnetite synthesis if they transport iron; they might perhaps act as an efflux system involved in transporting iron from the cytoplasm into magnetosomes. The system may be unique to strain MC-1, as only four RND superfamily genes were found in M. magnetotacticum and none were identified in the genomes of the other MTB strains.
Fifty-two of the 115 transporter-encoding genes in the genome of MC-1 are for putative ABC transport systems. ABC systems are ubiquitous in prokaryotes; most of these are transmembrane proteins that function in the transport of a wide variety of substrates, at the cost of ATP hydrolysis (30). The majority of ABC transporter genes in MC-1 appear to be for the export of compounds, with only two components (inner membrane/permease component and ATP-binding cassette [ATPase]). The genome indicates that MC-1 has a relatively small number of ABC transport systems for uptake (inner membrane/permease, ATPase, and periplasmic solute-binding components are all present), including transporters for the uptake of amino acids (two transporters [Mmc10486, Mmc10487, Mmc10488, Mmc10489, and Mmc10490 and Mmc11019, Mmc11020, Mmc11021, Mmc11022, and Mmc11023]), zinc (Mmc13180, Mmc13181, and Mmc13183), cobalt (Mmc13138, Mmc13139, and Mmc13720), molybdate/molybdenum (Mmc11183, Mmc11184, Mmc11732, and Mmc12961), and phosphate (Mmc11539, Mmc11540, and Mmc11541). ABC transporters typically have high affinities for solutes and consequently can target compounds at low concentrations in the environment. There is also a symporter available for the uptake of acetate by secondary transport (Mmc11042), which correlates with the ability of MC-1 to use this compound as a substrate (see above). In strain MC-1, several transporters appear to be specific for phosphate/phosphonate transport. Other predicted transporters include transport proteins for ammonium, phosphate, citrate, and sulfate ions.
We also identified several genes encoding proteins that might play a role in transport across the outer membrane; these include three ExbB-, two ExbD-, three TonB family protein-, four TonB-dependent receptor-, and three TolC-encoding genes. Translocation of proteins across the cytoplasmic membrane is apparently mediated by two different systems that can be identified from the MC-1 genome, namely, (i) the preprotein translocase system (SecABDEFGY system; Mmc13327, Mmc13422, Mmc13212, Mmc10834, Mmc13211, Mmc11884, and Mmc10867) and (ii) the TatABC (twin-arginine translocation) (Mmc11619 to Mmc11617) system. Genes encoding the major components included in type II secretion (Gsp proteins) (Mmc10931, Mmc10932, Mmc10933, and Mmc13531) are also present and presumably involved in protein translocation across the outer membrane.
Signal transduction.
The MC-1 genome contains a variety of genes that encode proteins involved in signal transduction processes, including the chemotaxis protein CheC (of which there are 3), 5 chemotaxis protein histidine kinases, 5 chemotaxis response regulators containing a CheY-like receiver domain, 9 chemotaxis signal transduction proteins, 7 CheY-like receivers, and 32 methyl-accepting chemotaxis proteins. The MC-1 genome appears to contain 34 PAS/PAC domain-encoding genes, which may function as redox sensors (120). The genome also harbors 73 genes encoding response regulator proteins that contain a GGDEF domain, an HD-GYP domain, a helix-turn-helix DNA-binding domain, an AAA-type ATPase and DNA-binding domain, or a winged-helix DNA-binding domain. In addition, 15 genes that encode proteins with both EAL and GGDEF domains are present. Proteins with these two domains have been described to be involved in the synthesis (GGDEF domain) and phosphorylation (EAL domains that function as phosphodiesterases) of cyclic di-GMP, a novel bacterial signaling system (91).
It is noteworthy that the genomes of MTB contain a large number of genes that encode signal transduction histidine kinases responsible for signal detection (1). For example, the genome of strain MC-1 contains 87 of this type of genes, while there are 107 signal transduction histidine kinases in M. magneticum and 212 in M. magnetotacticum. The genomes of nonmagnetotactic bacteria such as S. denitrificans, Thiomicrospira crunogena, and Nitrosococcus oceani generally contain much smaller numbers of signal transduction histidine kinases, with these three having 36, 17, and 18, respectively (101). Overall, the genomes of these three species also contain fewer genes encoding signaling proteins (101).
Regulation of gene transcription.
As previously discussed, the genome of MC-1 appears to contain many genes that encode receptors involved in the sensing of environmental changes, and therefore it is expected that the genome contains an equal or similar number of response regulator genes. In MTB, oxygen and iron concentrations are known to control magnetite biomineralization, whose level can be determined by measuring both magnetism and iron content of cells (21, 57, 97, 98). In M. gryphiswaldense, a high oxygen concentration and iron starvation downregulate mam and mms gene transcription (95). Characterizing the signal transduction pathways in MC-1 and other MTB is necessary to understand which environmental conditions trigger the specific biochemical pathways resulting in magnetosome magnetite biomineralization. The MC-1 genome contains four genes that encode members of the FNR (fumarate and nitrate reductase) family of transcription regulators. These regulators might function as redox sensors but also as DNA-binding transcription factors (103). An important regulator for iron uptake is Fur (ferric uptake regulator) (69). The MC-1 genome contains two genes that encode proteins that are members of the Fur family (Mmc10894 and Mmc13182). Another known iron uptake regulator is AraC, which is coded for by the Mmc10422 and Mmc10423 genes in MC-1. It has been postulated that AraC transcriptional regulators allow bacteria to fine-tune the expression of siderophore biosynthesis and transport genes in such a manner that maximum expression of these genes occurs only in environments where the corresponding ferrisiderophore complex can serve as an efficient iron source (23). It was also reported that regulators of the AraC family are able to activate the expression of multidrug efflux genes (80). An interesting possibility is that this regulator might be involved in the regulation of iron transport into the cell and then efflux into magnetosomes. Besides the major sigma factors,
70 and
54, we also identified the genes for
28 (FliA), which controls the transcription of flagellin genes, along with its antagonist, the anti-
28 factor FlgM.
Motility.
Like all magnetotactic cocci, cells of strain MC-1 always possess two bundles of flagella. A flagellar bundle consists of seven to nine separate flagella, each with a diameter of about 10 nm (44). Each bundle originates from a pit-like structure (D. A. Bazylinski and T. J. Beveridge, unpublished data). Genes encoding the necessary components of a complete flagellar apparatus (fli, flh, and flg genes) form five distinct clusters within the MC-1 genome (see Table S2 and Fig. S1 in the supplemental material), with the largest being 29.9 kb and containing 29 genes. Genes in these clusters encode proteins that make up the flagellar filament, basal body, and hook and two proteins that make up the motor of the flagellum (see Table S2 in the supplemental material).
Chemotaxis.
Like those of other MTB, cells of strain MC-1 passively align and swim along geomagnetic field lines. However, magnetotactic cocci exhibit a polar preference in their swimming direction that is dependent upon the local O2 (electron acceptor) concentration, the concentration of electron donor (e.g., S2–), and possibly redox conditions (44). Cells use magnetotaxis in conjunction with aerotaxis (i.e., magneto-aerotaxis) to more efficiently position themselves at an optimal local O2 concentration (44) where their metabolism also appears to be optimal and there is continual energy generation. The ability of cells to navigate to niches that support maximum energy levels has been termed energy taxis (2). In the case of MC-1, the optimal O2 concentration is very low, and cells form microaerobic bands at the OAI in semisolid O2-gradient cultures (44). Because of the cells' unidirectional swimming pattern and polar preference in swimming direction in magnetic fields, the term polar magneto-aerotaxis has been used to describe the behavior of cells of strain MC-1. MC-1 has been shown to use the magnetic field for direction as well as for an axis for a presumably more efficient aerotaxis response. This is in contrast to the response of the magnetospirilla, which upon repeated cultivation in liquid medium swim in both directions and appear to use the magnetic field only as an axis (i.e., they display axial magneto-aerotaxis) (44). Cells of strain MC-1 display strong aerotaxis under very low concentrations of oxygen and appear to quickly locate and maintain a position at the OAI even when a magnetic field is completely absent (Bazylinski, unpublished data).
In comparison to many chemotactic bacterial species (see Table S3 in the supplemental material), strain MC-1 contains a relatively large number of known chemotaxis genes, i.e., 65 of them, possibly suggesting an important role for these genes in the ecology of MC-1 (locating and maintaining their optimal position at the OAI), and probably other MTB as well. Indeed, the genomes of some Magnetospirillum species, for example, M. magneticum and M. magnetotacticum, contain more than twice the number of recognized chemotaxis genes (162 and 326 genes, respectively). Like strain MC-1, these species also show a strong aerotactic response, but they are not polar magneto-aerotactic. Analysis of other motile bacteria, however (see Table S3 in the supplemental material), reveals that the number of chemotaxis genes in MC-1 is not exceptionally high. When the number of chemotaxis genes is normalized to genome size (number of chemotaxis genes per Mb), strain MC-1 has a slightly larger number of these genes than most of the other motile bacteria examined. In a comparison of 29 motile prokaryotes, the number of chemotaxis genes per Mb ranged from 0 (Pyrococcus furiosus DSM 3638 and Dehalococcoides sp. strain BAV1) to 75.8 (M. magnetotacticum), while the average was 12.2 chemotaxis genes per Mb (see Table S3 in the supplemental material). MC-1 contains 13.8 chemotaxis genes per Mb.
The genome of strain MC-1 contains genes for 18 chemotaxis sensory transducers, 12 CheW domain proteins (8 CheW and 4 CheA), and 7 hemerythrin-like metal-binding proteins. Chemotaxis sensory transducers consist of two domains, a methyl-accepting chemotaxis protein and a histidine kinase (HAMP) domain. This protein is thought to transduce the signal to CheA. CheW is an adaptor protein and forms a complex with CheA. In other bacteria, CheA associates with a receptor (CheY) and the phosphate group from CheA is transferred to CheY. This activates the flagellum, where it modulates the motor reaction (63).
Hemerythrin proteins are oxygen-binding proteins common in invertebrates and a large number of prokaryotes, including members of both the Bacteria and the Archaea (46). Each has a nonheme di-iron site that reversibly binds oxygen. Two iron ions are bound to the protein via a motif comprised of seven conserved amino acid residues (62). Despite the ubiquity of these proteins in prokaryotes, specific functions of these proteins have not been elucidated, although a number of them have been proposed, including binding of oxygen or iron (or other metals) as a storage or detoxification mechanism (46). The putative methyl-accepting protein DcrH of Desulfovibrio vulgaris contains a C-terminal hemerythrin-like domain and has been ascribed a role as an oxygen sensor for aerotaxis in this organism (58, 118). Given the strong aerotactic response by MTB under conditions of low oxygen, it is tempting to speculate on a similar function in the MTB.
Many of the chemotaxis genes are clustered in the genome of strain MC-1 (see Table S2 and Fig. S1 in the supplemental material). There are four clusters of chemotaxis genes in MC-1. Each cluster, except one, contains one CheW, CheA, CheB, CheR, and CheY gene and a chemotaxis sensory transducer-encoding gene, although the genes are not in the same order in the three clusters. The fourth cluster is missing the CheA protein gene (see Fig. S1 in the supplemental material). The clustering of the che genes is conserved in many prokaryotes, and the respective encoded proteins represent a pathway common to most chemotactic bacteria and archaea (63, 115). Thus, it is likely that these proteins in MC-1 have similar functions to the proteins in other bacteria.
Metabolic functions. (i) Nitrogen metabolism.
We previously found that cells of strain MC-1 show nitrogenase activity (i.e., cells reduce acetylene to ethylene under nitrogen-limited conditions) (14) and thus can fix atmospheric nitrogen. The genome of MC-1 contains a cluster of 19 nitrogenase genes (nif) that include nifZVXNEYTKDHABQ (Mmc11189 to Mmc11209) (see Fig. S2 in the supplemental material). In addition to the nif genes, the cluster also contains dinitrogenase and dinitrogenase reductase genes, known to be involved in the nitrogen fixation process. The genes nifR (Mmc13670) and nifU (Mmc10723) are also present in the genome but are not part of this cluster. M. gryphiswaldense, M. magnetotacticum, and M. magneticum display nitrogenase activity as well (10, 11). The organization and presence of the nif genes in these strains are different from those in strain MC-1. A cluster of 20 genes encoding nitrogenases, dinitrogenases, and nitrogenase reductases is present in M. magnetotacticum. This cluster also contains 13 nif genes, including nifWVQXNEKDHHTZ and nifB (14). The genome of M. magneticum contains one cluster of 23 genes encoding enzymes involved in nitrogen fixation, including nifWQXXNHZ and nifB. Seven additional nif genes, such as nifR and nifU, are distributed throughout the genome in M. magneticum. The presence of specific genes within the genome indicates that there are three available pathways for the assimilation of ammonium (either taken up from the environment or synthesized de novo via nitrogen fixation), namely, glutamine synthetase and glutamate synthase (GS-GOGAT cycle); glutamate dehydrogenase; and alanine dehydrogenase.
Analysis of the MC-1 genome revealed two gene clusters encoding enzymes necessary for dissimilatory nitrate reduction and no apparent ORFs that might carry genes involved in assimilatory nitrate reduction. Magnetospirillum strains are capable of assimilatory and dissimilatory nitrate reduction and can grow either anaerobically or microaerobically (9) with nitrate as a terminal electron acceptor. All are known to be capable of complete denitrification (converting nitrate to dinitrogen gas). Cells of strain MC-1 do not grow anaerobically with nitrate as a terminal electron acceptor under chemolithoautotrophic (thiosulfate as an electron donor) or chemoorganoheterotrophic (acetate as an electron donor) conditions (14). However, MC-1 is fastidious with regard to growth in general, and it is possible that the culture conditions were not optimal for anaerobic growth with nitrate. Interestingly, nitrate added to the growth medium did not inhibit acetylene reduction by cells of strain MC-1 as it does in other organisms that fix dinitrogen (e.g., Magnetospirillum species) (10, 14; Bazylinski, unpublished data), further indicating that MC-1 cannot utilize nitrate as a sole source of nitrogen.
The dissimilatory nitrate reductases have been classified into periplasmic (Nap) enzymes and membrane-bound (Nar) enzymes (88, 107). Strain MC-1 appears to have the ability to synthesize a functional periplasmic nitrate reductase, as seven nap genes (napCBHGADF) (Mmc11587 to Mmc11593) are present in the genome of MC-1 (see Fig. S2 in the supplemental material). Prokaryotic nap gene clusters are complex and diverse. M. magnetotacticum and M. magneticum each contain a nap gene cluster of seven genes (napFDAGHBC) similar to that in strain MC-1, but in reversed order. The napFDAGHBC arrangement of nap genes is also present in E. coli (59), in which napA, napB, napC, and napD are essential for enzymatic activity, while napF, napG, and napH are not (83). In contrast, the genome of Bradyrhizobium japonicum contains a different arrangement, napEDABC, and all genes are essential for nitrate respiration (33).
A gene encoding nitrite reductase, the enzyme that defines denitrification and catalyzes the reduction of nitrite to nitric oxide (NO), appears to be missing from the genome. Two types of NO reductases (Nor), catalyzing the reduction of NO to nitrous oxide (N2O), have been classified, including (i) cNor, which receives electrons from cytochrome c; and (ii) qNor, which receives electrons from quinol (15, 121). A gene cluster is present in the genome of MC-1 that includes norCBQ, an Fe-S cluster-binding protein gene, and norD (Mmc10121 to Mmc10117). The norCB genes encode the small and large subunits of cNor, respectively (15). However, even though the cluster contains all genes required for cNor assembly and activation, activity of this protein has not yet been observed. Nitrous oxide reductase, the enzyme catalyzing the last step for complete denitrification (the reduction of N2O to N2), appears to be missing from the genome.
(ii) Sulfur metabolism.
Our analysis shows that the genes encoding the first two enzymes of the assimilatory sulfate reduction pathway are present in the genome of MC-1. These are the small and large subunits of sulfate adenylyltransferase, encoded by the cysN and cysD genes (Mmc11018 and Mmc11017), respectively. This enzyme activates sulfate and catalyzes the synthesis of adenosine-5'-phosphosulfate (APS). APS becomes phosphorylated by the APS kinase, encoded by cysC (Mmc12549), and 3'-phosphoadenosine-5'-phosphosulfate (PAPS) is produced. Genes encoding the PAPS reductase, CysH, which catalyzes the conversion of PAPS to sulfite, as well as the sulfite reductase CysI, which reduces sulfite to sulfide, appear to be missing from the genome, and thus MC-1 appears to be incapable of assimilatory sulfate reduction to H2S. MC-1 should still be able to synthesize cysteine from serine, however, if another source of S2– is present (possibly from thiosulfate in the growth medium) or is derived from sulfate in a pathway different from the recognized pathway of assimilatory sulfate reduction, because the genome contains cysE (gene for serine O-acetyltransferase; Mmc13065) and cysK (gene for cysteine synthase; Mmc13416). Despite MC-1's apparent ability to synthesize cysteine, this amino acid is required for the growth of MC-1 in [O2]-gradient medium, at least as a reducing agent. It may also utilize cysteine as a general source of cell sulfur, which might partially explain the requirement for the amino acid, based on the presumed inability of the organism to utilize sulfate as a sole source of sulfur (based on the absence of specific genes as described above). However, although it was present in the growth medium, cysteine was not originally added to the medium as a source of sulfur but as a reducing agent, and cysteine may have fortuitously served a dual role in the growth medium. Alternatively, perhaps other novel genes are present in MC-1 that encode currently unrecognized enzymes for the reduction of PAPS to sulfide. In contrast, the genomes of M. magneticum and M. magnetotacticum contain all genes necessary for assimilatory sulfate reduction, including the genes encoding PAPS reductase (cysH) and assimilatory sulfite reductase (cysI).
There is no evidence for the presence of any of the genes for dissimilatory sulfate reduction in the MC-1 genome, consistent with the fact that cells of MC-1 do not grow anaerobically with sulfate as a terminal electron acceptor.
Cells of MC-1 are able to use thiosulfate or sulfide as an electron donor when they grow chemolithoautotrophically (14, 113). The MC-1 genome contains 8 of the 15 sox genes, which in Paracoccus pantotrophus encode enzymes for sulfur oxidation (47). These genes include soxXYZAB (Mmc11905 to Mmc11909), soxW (Mmc13719 and Mmc13718), soxE (Mmc12995), and soxF (Mmc12996). In addition, the MC-1 genome contains another copy of the soxA and soxX genes (Mmc13714 and Mmc13717) and two additional copies of the soxZ and soxY genes (Mmc12998, Mmc13715, Mmc12997, and Mmc13716) (see Fig. S3 in the supplemental material). SoxZY are known to interact with SoxAB, and their duplication may facilitate differential regulation of different loci, as proposed for specific homologues in S. denitrificans (101). The sox genes are organized into three clusters located at three different positions on the chromosome (see Fig. S3 in the supplemental material). A similar organization of different clusters also appears in S. denitrificans (101). The soxXYZAB genes encode three periplasmic proteins (38, 85), and soxF encodes a monomeric flavoprotein that has sulfide dehydrogenase activity (7, 84). soxW encodes a periplasmic thioredoxin described as essential for chemotrophic growth with thiosulfate (6). However, the soxCD genes that encode the sulfur dehydrogenase, which are important for electron transfer (47), appear to be missing from the genome. This might indicate that MC-1 uses soxBXAYZ gene products to oxidize thiosulfate to sulfur, which is either stored inside the cell or excreted (56). In the case of MC-1, sulfur globules accumulate intracellularly. Another possibility is that strain MC-1 accomplishes sulfur oxidation in combination with the dsr genes, as described for Chlorobaculum tepidum (47). The MC-1 genome contains nine of the dsr genes, including dsrAB (Mmc12156 and Mmc12157, respectively), dsrC (Mmc10048), and dsrEFHJOP (Mmc11704 to Mmc11709). The complete dsr gene cluster in Allochromatium vinosum is comprised of 15 genes that are essential for sulfur oxidation in that species (28). The M. magneticum chromosome contains 12 dsr genes (dsrABEFHCMKLOOP) organized in a single cluster. In contrast, only three dsr genes (dsrABC) were found in M. magnetotacticum. In addition to these sulfide-oxidizing enzyme-encoding genes, we also identified seven genes that encode sulfite oxidases, including yedY and yedZ (Mmc12485 and Mmc12484), which are also present in M. magneticum.
(iii) Phosphonate metabolism.
Phosphorus is a growth-limiting factor in many ecosystems, particularly some marine habitats (109, 110, 116). Therefore, it is especially important for marine bacteria to degrade organic phosphates and phosphonates for the release and assimilation of phosphorus. In other characterized systems, the uptake of phosphate is controlled by the Pho regulon, a set of genes and operons that are regulated by the extracellular concentration of phosphate (112). The regulation of this response is controlled by the PstSCAB transporter, which "senses" different Pi levels and communicates through PhoU to the two-component system PhoR and PhoB (74, 75). A similar regulation might occur in MC-1, since its genome contains all of the genes for these components. However, the pst cluster is interrupted by a rhodanese domain-containing protein gene, as follows: pstCAB-phoU-RHOD-pstS (Mmc11539 to Mmc11544); phoRB (Mmc11492 and Mmc11491) are further upstream. Another operon that might be regulated by the Pho regulon, as described for other characterized systems, contains the phn genes that regulate phosphonate uptake. The MC-1 genome contains phnC (Mmc10734), phnA (Mmc12533), and four additional phosphonate transporter genes (Mmc11819, Mmc11923, Mmc13404, and Mmc13535). The pstSCAB, phoU, and phoRB transporter system is also present in M. magnetotacticum, M. magneticum, and M. gryphiswaldense.
All three Magnetospirillum strains also contain genes encoding polyphosphate kinase and exopolyphosphatase. These enzymes are involved in the accumulation of phosphate in polyphosphate granules, which have been shown to be present in these species (99; this study). Cells of strain MC-1 also produce polyphosphate inclusions, and we identified two genes that might contribute to the synthesis of these structures, as they are putative exopolyphosphatase-encoding genes (Mmc11700 and Mmc12960).
(iv) Carbon fixation and central carbon metabolism.
Cells of strain MC-1 have been shown to grow chemolithoautotrophically with sulfide or thiosulfate as an electron donor (113). It was initially suspected that strain MC-1 used the Calvin-Benson-Bassham cycle for autotrophy, as shown for nearly all aerobic chemolithoautotrophic bacteria (14). However, genes with significant homology to either form I or form II ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) genes were not identified in the MC-1 genome, and Rubisco activity was not observed during autotrophic growth (113). Recently, it was reported that genes for the enzymes of the rTCA cycle were present in MC-1, and activities of these enzymes were demonstrated (113). Thus, we infer that MC-1 cells use only the rTCA cycle for CO2 fixation. Our genome analysis shows that the genes for all enzymes necessary for the TCA cycle to operate in the reductive mode (101) are present in strain MC-1 (see Fig. S4 in the supplemental material). These include pyruvate:ferredoxin oxidoreductase (Mmc11741 and Mmc11742), 2-oxoglutarate:ferredoxin oxidoreductase (Mmc11749 and Mmc11750), and ATP-dependent citrate lyase (Mmc13638 and Mmc13639).
Acetyl-coenzyme A (acetyl-CoA) and oxaloacetate produced by the rTCA cycle likely feed into central carbon metabolism. Given that acetate can support heterotrophic growth of MC-1, it is probably converted to acetyl-CoA through the action of acetyl-CoA synthetase (Mmc10598). Neither of the genes associated with operation of the glyoxylate bypass (encoding isocitrate lyase and malate synthase) are present in the MC-1 genome. Acetyl-CoA is likely converted to pyruvate via pyruvate:ferredoxin oxidoreductase (see above). Oxaloacetate could be used to form phosphoenolpyruvate via phosphoenolpyruvate carboxykinase (Mmc10369) or fed into the TCA cycle in the forward (oxidative) direction. Acetyl-CoA is also likely used in the synthesis of fatty acids (acetyl-CoA carboxylase; Mmc12330).
Carbon fixed by the rTCA cycle is likely converted to pyruvate that can then be channeled in different directions for the synthesis of lipids, proteins, nucleotides, and cell wall components. All of the enzymes necessary to direct carbon through gluconeogenesis are encoded by the MC-1 genome. A gene for the key enzyme fructose-1,6-bisphosphatase (EC 3.1.3.11) was not recognized by the automated annotation, but one of the two genes annotated as inositol monophosphatase genes (Mmc10078) appears to be an archaeal protein-type fructose-1,6-bisphosphatase gene. All of the genes necessary for glycolysis (albeit only from phosphorylated sugars) are present, even though MC-1, like other MTB strains investigated thus far, is unable to utilize sugars as sole sources of carbon (14, 113). Of special note in the glycolytic pathway of MC-1 is the presence of a 2,3-bisphosphoglycerate-independent phosphoglycerate mutase. The glycolytic pathway is used to generate metabolic intermediates to replenish levels of cellular reductant that can be used by the rTCA cycle or the citric acid cycle.
The genes encoding glucose permease (EC 2.7.1.69) and components of the phosphotransferase system are not present in the MC-1 genome. This suggests that the reason that MC-1 cannot utilize sugars as a source of energy and carbon is that cells simply are not able to take up sugars from the extracellular environment. Genes encoding the enzymes necessary to direct pyruvate to the citric acid cycle are present. The gene encoding citrate synthase (Mmc10897), a key enzyme of the oxidative citric acid cycle, is present, and therefore the citric acid cycle is complete in strain MC-1.
Other cellular processes.
As expected, the MC-1 genome contains all genes necessary for the F1-F0 ATP synthase and for oxidative phosphorylation for energy. We identified all genes necessary for the synthesis of the 20 common amino acids. All 20 tRNAs are present as well. MC-1 also contains the genes necessary to synthesize nucleotides, fatty acids, phospholipids, peptidoglycan, and vitamins and cofactors, such as riboflavin, nicotinic acid, pantothenate, biotin, and heme. However, no genes encoding the biosynthesis of enzymes involved in secondary metabolism or enzymes for the biodegradation and metabolism of xenobiotics were identified in the MC-1 genome.
Oxygen as a terminal electron acceptor and oxygen stress.
Experimental evidence shows that cells of strain MC-1 use only O2 as a terminal electron acceptor, not nitrate, nitrite, nitrous oxide, sulfate, sulfite, thiosulfate, dimethyl sulfoxide, or trimethylamine oxide (Bazylinski, unpublished data). Genomic data appear to also support this (with the exception of nitrate reduction), and the Mmc12353, Mmc12354, and Mmc12355 genes appear to represent the three subunits of a cytochrome c oxidase which could be the terminal oxidase for O2 respiration.
Genes encoding superoxide dismutase and catalase, proteins involved in the detoxification of toxic, reactive by-products of oxygen respiration found in almost all aerobic bacteria, were not present in strain MC-1. The only genes involved in the detoxification of reactive oxygen species that were identified were those that encode cytochrome c peroxidase (Mmc12947, Mmc11488, and Mmc12527) and alkyl hydroperoxide reductase (Mmc11472 and Mmc11679). These enzymes remove hydrogen peroxide from the cell. The small number of genes to detoxify reactive oxygen species might explain the obligately microaerophilic nature of strain MC-1.
Prophages.
Another striking feature of the MC-1 genome is the large number (14) of predicted prophages (22) distributed widely over the entire genome, with lengths of 5,287 to 42,314 bp and G+C contents of 52.14 to 58.34% (see Table S4 in the supplemental material). One example of these, predicted prophage 7, is depicted in Fig. S5 of the supplemental material. Prophage 7 is 42,314 bp in length, contains 54 coding sequences, and is transcribed in opposite orientation. This prophage "genome" starts with a lysogeny module containing an integrase gene (Mmc12087) and a cI-like repressor gene (Mmc12084). This type of lysogeny module is characteristic of lambdoid prophages (73). In most cases, the lysogeny module is located between lysis and DNA replication modules (73). However, no genes encoding lysine and holin, representing the lysis module, could be identified. Downstream of the lysogeny module, a DNA replication module was identified. This cluster includes all genes involved in DNA replication and metabolism, such as peptidase (Mmc12081), helicase (Mmc12074), DNA primase (Mmc12070), DNA-binding protein (Mmc12066), and DNase I (Mmc12065) genes. Some genes encoding regulatory proteins that control gene expression are scattered throughout the prophage genome. These early expressed genes encode nonvirion components and have nonprophage bacterial homologues. Some phage virion assembly proteins encoded by the late operons are more conserved, including the large terminase, portal protein, head maturation protein, and tail tape-measure protein (24). Even though the above proteins appear to be more highly conserved, many of the virion assembly proteins are very variable and difficult to identify from sequence information (24, 90). The phage genomes show striking gene clustering with a well-conserved gene order that we were able to identify in prophage 7 as well. The late operons start with the gene for a terminase (Mmc12050) that cleaves virion-length molecules from replicating DNA and is involved in packaging the phage genome into the phage head. This gene is followed by a cluster of genes that encode head structural components and assembly components, including the portal protein (Mmc12048), which forms the hole through which DNA is packaged into the capsid, and head maturation protease (Mmc12046), which controls proteolytic cleavage of the phage head. A cluster of genes encoding tail structure components and assembly components follows the capsid assembly genes. The most significant of these is the tail tape-measure protein (Mmc12034), which determines the length of the tail shaft (24). The high similarity of the ORFs to those of lambdoid temperate phages, which are generally members of the Siphoviridae, allows us to assume that MC-1 prophage 7 belongs to this family.
In contrast to that of strain MC-1, the genome of M. magneticum AMB-1 contains 12 putative prophages and that of M. gryphiswaldense MSR-1 contains 6, with its plasmid containing 1 (this study). To our knowledge, this is the first time that prophages have been reported within the genomes of MTB. The presence of these prophages may prove important in future studies (e.g., genetics) involving MTB. Prophages might also be another cause of horizontal gene transfer of the MAI. During times of stress, phage genes may become activated and cause phages to enter the lytic cycle. During this process, some phages may encapsulate bacterial genes and pass them on to new microbes (114). In addition, the MC-1 chromosome harbors 135 transposase/integrase-related genes that are distributed throughout the chromosome and not localized in hot spots as in M. magneticum AMB-1 (78).
Conclusions.
The marine magnetotactic coccus strain MC-1 has a number of unique features that clearly differentiate it from other MTB whose genomes have been sequenced. Although the size of the genome of MC-1 is comparable to those of other MTB, the organization of the magnetosome genes is quite different. In addition, there are a number of genes within the MAI of strain MC-1 that may represent magnetosome genes not found in other MTB, and the large number of mobile elements found in and surrounding the MAI in other MTB are not present in MC-1. There are also a large number of predicted prophage sequences in the genome of strain MC-1. Although the lack of mobile elements (e.g., insertion sequences) in and surrounding the MAI might suggest that magnetosome genes in MC-1 did not originate from lateral gene transfer, the presence of these prophages might suggest, in fact, that they did or that they are transferred from MC-1 to other prokaryotes. Lastly, genomic data support the lack of metabolic versatility of strain MC-1 (e.g., in the use of different sources of carbon and electron sources and donors) compared to Magnetospirillum species.
On the other hand, there are some significant similarities between MC-1 and Magnetospirillum species. For example, there are numerous chemotaxis proteins probably used in energy taxis (2). Some of these have redox sensory domains which likely function to efficiently position cells at their optimal position at the OAI, thereby enabling them to obtain the electron donors (S2–) and acceptors (O2) required for energy generation and growth (44). Magnetotaxis and the ability to biomineralize magnetosomes are currently thought to make chemotaxis even more efficient for organisms such as MTB, whose survival seems to rely on locating and maintaining an optimal position in vertical chemical gradients in aquatic habitats (44). In sum, the information derived from the genome of MC-1 helps to emphasize the metabolic and physiological diversity of MTB and highlights the unique adaptations of magnetotactic cocci such as MC-1, for example, its exquisitely sensitive, fine-tuned aerotactic response. In addition, genomic information from MC-1 might provide insights leading to the isolation of other magnetotactic cocci.
Sequencing was performed under the auspices of the U.S. Department of Energy's Office of Science, Biological and Environmental Research Program and was supported by the University of California Lawrence Berkeley National Laboratory (contracts DE-AC02-05CH11231 and DE-AC52-07NA27344) and the Los Alamos National Laboratory (contract DE-AC02-06NA25396). This work was also partially supported by U.S. National Science Foundation grant EAR-0715492 to D.A.B.
Published ahead of print on 22 May 2009. ![]()
Present address: INM, Leibniz Institute for New Materials, 66123 Saarbrücken, Germany. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»