Complete genome sequence of Acetohalobium arabaticum type strain (Z-7288T)

Acetohalobium arabaticum Zhilina and Zavarzin 1990 is of special interest because of its physiology and its participation in the anaerobic C1-trophic chain in hypersaline environments. This is the first completed genome sequence of the family Halobacteroidaceae and only the second genome sequence in the order Halanaerobiales. The 2,469,596 bp long genome with its 2,353 protein-coding and 90 RNA genes is a part of the Genomic Encyclopedia of Bacteria and Archaea project.


Introduction
Strain Z-7288 T (= DSM 5501 = ATCC 49924) is the type strain of the species Acetohalobium arabaticum, which is the type species of the genus Acetohalobium [1,2]. The genus name derives from the Latin word 'acetum', meaning vinegar, and the Greek words 'halos' and 'bios', meaning salt and life, respectively, in order to indicate an acetateproducing organism living in salt [3]. The species name derives from Arabat, a peninsula between the Sea of Azov and Sivash [3], since the strain was isolated from lagoons of the Arabat spit (East Crimea) which separates Sivash lake from the Sea of Azov [2]. Currently, this is the only known strain in the genus Acetohalobium. A. arabaticum partici-pates together with other halophilic bacteria and the genera Methanohalophilus and Methanohalobium in the C1-trophic chain in hypersaline environments [2]. Here we present a summary classification and a set of features for A. arabaticum Z-7288 T , together with the description of the complete genomic sequencing and annotation.

Classification and features
The cells of A. arabaticum are bent rods, motile by one to two subterminal flagella (Table 1) [2]. The flagella are stated in the original description [2], though they are not visible in our study ( Figure 1). The cells are single, in pairs or form short chains, being 0.7-1 µm in diameter and 1-5 µm in length [2]. Other typical cell aggregates are palisades and ribbons, which are formed by adhesion of cells having intimate contact ( Figure 1) [2]. The multiplication is by binary fission. The outer membrane is typical of a Gram-negative organism [2]. Growth is completely inhibited by 100 µM/ml streptomycin, benzylpenicillin, bacitracin, erythromycin, gentamycin, kanamycin, vancomycin or tetracyclin [2]. Strain Z-7288 T is obligately anaerobic, tolerating up to 12 mM H2S. Neither O2, S2O3 2-, SO4 2-, nor S 0 can serve as electron acceptors. Strain Z-7288 T requires a salt concentration of 10-25% NaCl, the optimum is 15-18% NaCl [2]. The optimal pH is between 7.6 and 8.0 [2].
A. arabaticum exhibits three modes of nutrition [2]: It is chemolithoautotrophic using H2 together with CO2 or CO; it is methylotrophic using trimethylamine (TMA); and it is organotrophic using betaine, lactate, pyruvate or histidine. Carbohydrates are not utilized. No growth occurs on methanol, monomethylamine (MMA), dimethylamine (DMA), dimethylglycine, choline or sarcosine [2]. When grown on TMA, an equimolar amount of acetate is formed along with lesser amounts of DMA and MMA [2]. Betaine is degraded mainly to acetate and minor amounts of methylamines [2]. The 16S rRNA genes of the other type strains in the family Halobacteroidaceae share between 85.9% (Orenia sivashensis [19]) and 95.1% (Sporohalobacter lortetii [20]) sequence identity with strain Z-7288 T [21]. Uncultured clone sequences from environmental samples and metagenomic surveys do not surpass 84-86% sequence similarity to the 16S rRNA gene sequence of strain Z-7288 T , indicating a lack of further members of the genus Acetohalobium in the habitats screened thus far (status June 2010). Figure 2 shows the phylogenetic neighborhood of A. arabaticum Z-7288 T in a 16S rRNA based tree. The sequences of the five 16S rRNA gene copies in the genome of Acetohalobium arabaticum Z-7288 T differ from each other by up to one nucleotide, and differ by up to three nucleotides from the previously published 16S rRNA sequence generated from DSM 5501 (X89077).

Chemotaxonomy
No chemotaxonomic data are currently available for the genus Acetohalobium.  Phylogenetic tree highlighting the position of A. arabaticum Z-7288 T relative to the type strains of the other genera within the order Halanaerobiales. The trees were inferred from 1,308 aligned characters [22,23] of the 16S rRNA gene sequence under the maximum likelihood criterion [24] and rooted in accordance with the current taxonomy [25]. The branches are scaled in terms of the expected number of substitutions per site. Numbers above branches are support values from 300 bootstrap replicates [26] if larger than 60%. Lineages with type strain genome sequencing projects registered in GOLD [27] are shown in blue, published genomes [28] in bold.

Genome project history
This organism was selected for sequencing on the basis of its phylogenetic position [29], and is part of the Genomic Encyclopedia of Bacteria and Archaea project [30]. The genome project is deposited in the Genome OnLine Database [27] and the complete genome sequence is deposited in GenBank. Sequencing, finishing and annotation were performed by the DOE Joint Genome Institute (JGI). A summary of the project information is shown in Table 2.

Growth conditions and DNA isolation
A. arabaticum Z-7288 T , DSM 5501, was grown anaerobically in DSMZ medium 494 (Acetohalobium medium) [31] at 37°C. DNA was isolated from 0.5-1 g of cell paste using MasterPure Gram Positive DNA Purification Kit (Epicentre MGP04100). Two µl lysozyme and five µl mutanolysin were added to the standard lysis solution for 40min at 37°C followed by 1 hour incubation on ice after the MPC-step.

Genome sequencing and assembly
The genome of A. arabaticum Z-7288 T was sequenced at using a combination of Illumina and 454 technologies. An Illumina GAii shotgun library with reads of 483 Mb a 454 Titanium draft library with average read length of 341 bases, and a paired end 454 library with average insert size of 10 kb were generated for this genome. All general aspects of library construction and sequencing can be found at http://www.jgi.doe.gov/. Illumina sequencing data was assembled with VELVET and the consensus sequences were shredded into 1.5 kb overlapped fake reads and assembled together with the 454 data. Draft assemblies were based on 241 Mb 454 draft data, and 454 paired end data. Newbler parameters are -consed -a 50 -l 350 -g -m -ml 20. The initial assembly contained 72 contigs in one scaffold. The initial 454 assembly was converted into a phrap assembly by making fake reads from the consensus, collecting the read pairs in the 454 paired end library. The Phred/Phrap/Consed software package was used for sequence assembly and quality assessment in the following finishing process [32]. After the shotgun stage, reads were assembled with parallel phrap (High Performance Software, LLC). Possible mis-assemblies were corrected with gapResolution (http://www.jgi.doe.gov/), Dupfinisher [32], or sequencing cloned bridging PCR fragments with subcloning or transposon bombing (Epicentre Biotechnologies, Madison, WI). Gaps between contigs were closed by editing in Consed, by PCR and by Bubble PCR primer walks (J.-F. Cheng, unpublished). A total of 292 additional reactions were necessary to close gaps and to raise the quality of the finished sequence. Illumina reads were used to improve the final consensus quality using an inhouse developed tool (the Polisher [33], ). The completed genome sequences have an error rate of less than 1 in 100,000 bp. Evidence codes -IDA: Inferred from Direct Assay (first time in publication); TAS: Traceable Author Statement (i.e., a direct report exists in the literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from of the Gene Ontology project [15]. If the evidence code is IDA, then the property was directly observed by one of the authors or an expert mentioned in the acknowledgements. Standards in Genomic Sciences

Genome annotation
Genes were identified using Prodigal [34] as part of the Oak Ridge National Laboratory genome annotation pipeline, followed by a round of manual curation using the JGI GenePRIMP pipeline [35]. The predicted CDSs were translated and used to search the National Center for Biotechnology Information (NCBI) nonredundant database, Uni-Prot, TIGRFam, Pfam, PRIAM, KEGG, COG, and In-terPro databases. Additional gene prediction analysis and functional annotation was performed within the Integrated Microbial Genomes -Expert Review (IMG-ER) platform [36].

Genome properties
The genome consists of a 2,469,596 bp long chromosome with a 36.6% GC content (Table 3 and Figure 3). Of the 2,443 genes predicted, 2,353 were protein-coding genes, and 90 RNAs; Seventyone pseudogenes were also identified. The majority of the protein-coding genes (76.4%) were assigned a putative function while the remaining ones were annotated as hypothetical proteins. The distribution of genes into COGs functional categories is presented in Table 4.