The Genome of Deep-Sea Vent Chemolithoautotroph Thiomicrospira crunogena XCL-2

Presented here is the complete genome sequence of Thiomicrospira crunogena XCL-2, representative of ubiquitous chemolithoautotrophic sulfur-oxidizing bacteria isolated from deep-sea hydrothermal vents. This gammaproteobacterium has a single chromosome (2,427,734 base pairs), and its genome illustrates many of the adaptations that have enabled it to thrive at vents globally. It has 14 methyl-accepting chemotaxis protein genes, including four that may assist in positioning it in the redoxcline. A relative abundance of coding sequences (CDSs) encoding regulatory proteins likely control the expression of genes encoding carboxysomes, multiple dissolved inorganic nitrogen and phosphate transporters, as well as a phosphonate operon, which provide this species with a variety of options for acquiring these substrates from the environment. Thiom. crunogena XCL-2 is unusual among obligate sulfur-oxidizing bacteria in relying on the Sox system for the oxidation of reduced sulfur compounds. The genome has characteristics consistent with an obligately chemolithoautotrophic lifestyle, including few transporters predicted to have organic allocrits, and Calvin-Benson-Bassham cycle CDSs scattered throughout the genome

Transporter Gene Frequencies within the Genomes of <i>Thiom. crunogena</i> XCL-2 (Arrow) and Other Proteobacteria

<p><i>Nitrob. winogradskyi (Nitrobacter winogradskyi)</i> is an alphaproteobacterium, <i>Nitros. europaea (Nitrosomonas europaea)</i> is a betaproteobacterium, and <i>Nitrosoc. oceani (Nitrosococcus oceani)</i> and <i>Methylo. capsulatus (Methylococcus capsulatus)</i> are gammaproteobacteria. Bars for intracellular pathogens are lighter red than the other heterotrophic gammaproteobacteria.</p

Circular Map of the <i>Thiom. crunogena</i> XCL-2 Genome

<p>The outer two rings (rings 1 and 2) are protein-encoding genes, which are color-coded according to COG category. Rings 3 and 4 are tRNA and rRNA genes. Ring 5 indicates the location of a prophage (magenta), phosphonate/heavy metal resistance island (cyan), and four insertion sequences (red; two insertion sequences at 2028543 and 2035034 are superimposed on this figure). The black circle indicates the deviation from the average %GC, and the purple and green circle is the GC skew (= [G − C]/[G + C]). Both the %GC and GC skew were calculated using a sliding window of 10,000 bp with a window step of 100.</p

Prophage Genome within the <i>Thiom. crunogena</i> XCL-2 Genome

<p>Lysogenic and lytic genes are delineated, as are predicted gene functions.</p

Numbers of MCP Genes in <i>Thiom. crunogena</i> XCL-2 and Other Proteobacteria

<p><i>Thiom. crunogena</i> is marked with an arrow.</p

Cell Model for <i>Thiom. crunogena</i> XCL-2, with an Emphasis on Ultrastructure, Transport, Energy, Carbon Metabolism, and Chemotaxis

<p>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 are yellow, and abbreviations are as follows: bc₁, <i>bc₁</i> complex; cbb₃, <i>cbb₃</i>-type cytochrome C oxidase; cytC, cytochrome C; NDH, NADH dehydrogenase; Sox, Sox system; UQ, ubiquinone. MCPs are fuchsia, as are MCPs with PAS domains or PAS folds. Influx and efflux transporter families with representatives in this genome are indicated on the figure, with the number of each type of transporter in parentheses. ATP-dependent transporters are red, secondary transporters are sky blue, ion channels are light green, and unclassified transporters are purple. Abbreviations for transporter families are as follows: ABC, ATP-binding cassette superfamily; AGCS, alanine or glycine:cation symporter family; AMT, ammonium transporter family; APC, amino acid-polyamine-organocation family; ATP syn, ATP synthetase; BASS, bile acid:Na⁺ symporter family; BCCT, betaine/carnitine/choline transporter family; CaCA, Ca²⁺:cation antiporter family; CDF, cation diffusion facilitator family; CHR, chromate ion transporter family; CPA, monovalent cation:proton antiporter-1, −2, and −3 families; DAACS, dicarboxylate/amino acid:cation symporter family; DASS, divalent anion:Na⁺ symporter family; DMT, drug/metabolite transporter superfamily; FeoB, ferrous iron uptake family; IRT, iron/lead transporter superfamily; MATE, multidrug/oligosaccharidyl-lipid/polysaccharide (MOP) flippase superfamily, MATE family; McsS, small conductance mechanosensitive ion channel family; MFS, major facilitator superfamily; MgtE, Mg²⁺ transporter-E family; MIT, CorA metal ion transporter family; NCS2, nucleobase:cation symporter-2 family; NRAMP, metal ion transporter family; NSS, neurotransmitter:sodium symporter family; P-ATP, P-type ATPase superfamily; Pit, inorganic phosphate transporter family; PNaS, phosphate:Na⁺ symporter family; PnuC, nicotamide mononucleotide uptake permease family; RhtB, resistance to homoserine/threonine family; RND, resistance-nodulation-cell division superfamily; SSS, solute:sodium symporter family; SulP, sulfate permease family; TRAP, tripartite ATP-independent periplasmic transporter family; TRK, K⁺ transporter family; VIC, voltage-gated ion channel superfamily.</p

<i>Thiom. crunogena</i> XCL-2 Phosphonate Operon

<p>The <i>DBR-1</i> genes are identical to each other, as are the <i>DBR-2</i> genes. Gene abbreviations are: <i>chp</i>, conserved hypothetical protein; <i>DBR-1</i> and -<i>2</i>, DNA breaking-rejoining enzymes; <i>hyp</i>, hypothetical protein; <i>phnFDCEGHIJKLMNP</i>, phosphonate operon. An asterisk (*) marks the location of a region (within and upstream of tRNA-phe) with a high level of similarity to the 5′ ends of the two direct repeat sequences noted in the figure. The transposase and integrase are actually a single CDS separated by a frameshift.</p

Models for Glycolysis, Gluconeogenesis, and the CAC in <i>Thiom. crunogena</i> XCL-2

<p>Models for central carbon metabolism for cells under environmental conditions with (A) sufficient reduced sulfur and oxygen; (B) sulfide scarcity; and (C) oxygen scarcity; green arrows represent the two “non-canonical” CAC enzymes, 2-oxoglutarate oxidoreductase (2-OG OR) and malate: quinone oxidoreductase (MQO).</p

Calvin-Benson-Bassham Cycle Gene Organization in Proteobacteria

<p>RubisCO genes <i>(cbbLS</i> and <i>cbbM)</i> are green, phosphoribulokinase genes <i>(cbbP)</i> are red, other genes encoding Calvin-Benson-Bassham cycle enzymes are black, and carboxysome structural genes are grey. For species in which <i>cbbP</i> is not near <i>cbbLS</i> or <i>cbbM,</i> the distance from the RubisCO gene to <i>cbbP</i> in kbp is indicated in parentheses. <i>Thiob. denitrificans</i> has two <i>cbbP</i> genes, so two distances are indicated for this species. Names of organisms that are unable to grow well as organoheterotrophs are boxed. Abbreviations and accession numbers for the 16S sequences used to construct the cladogram are as follows: <i>A. ehrlichei, Alkalilimnicola ehrlichei,</i> AF406554; Brady. sp., Bradyrhizobium sp., AF338169;<i>B. japonicum, Bradyrhizobium japonicum,</i> D13430; <i>B. xenovorans, Burkholderia xenovorans</i>, U86373; <i>D. aromatica, Dechloromonas aromatica,</i> AY032610; <i>M. magneticum, Magnetospirillum magneticum,</i> D17514; <i>M. capsulatus, Methylococcus capsulatus</i> BATH, AF331869; <i>N. hamburgensis, Nitrobacter hamburgensis,</i> L11663; <i>N. winogradskyi, Nitrobacter winogradskyi,</i> L11661; <i>N. oceani, Nitrosococcus oceani,</i> AF363287; <i>N. europaea, Nitrosomonas europaea,</i> BX321856; <i>N. multiformis, Nitrosospira multiformis,</i> L35509; <i>P. denitrificans, Paracoccus denitrificans,</i> X69159; <i>R. sphaeroides, Rhodobacter sphaeroides,</i> CP000144; <i>R. ferrireducens, Rhodoferax ferrireducens,</i> AF435948; <i>R. palustris, Rhodopseudomonas palustris,</i> NC 005296; <i>R. rubrum, Rhodospirillum rubrum,</i> D30778; <i>R. gelatinosus, Rubrivivax gelatinosus,</i> M60682; <i>S. meliloti, Sinorhizobium meliloti,</i> D14509; <i>T. denitrificans, Thiobacillus denitrificans,</i> AJ43144; <i>T. crunogena, Thiomicrospira crunogena,</i> AF064545. The cladogram was based on an alignment of 1,622 bp of the 16S rRNA genes, and is the most parsimonious tree (length 2,735) resulting from a heuristic search with 100 replicate random step-wise addition and TBR branch swapping (PAUP*4.0b10 [<a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.0040383#pbio-0040383-b113" target="_blank">113</a>]). Sequences were aligned using ClustalW [<a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.0040383#pbio-0040383-b114" target="_blank">114</a>], as implemented in BioEdit. Percent similarities and identities for <i>cbbL, cbbM,</i> and <i>cbbP</i> gene products, as well as gene locus tags, are provided as supporting information (<a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.0040383#pbio-0040383-st004" target="_blank">Table S4</a>).</p