Draft genome sequence of Gluconobacter thailandicus NBRC 3257.

Gluconobacter thailandicus strain NBRC 3257, isolated from downy cherry (Prunus tomentosa), is a strict aerobic rod-shaped Gram-negative bacterium. Here, we report the features of this organism, together with the draft genome sequence and annotation. The draft genome sequence is composed of 107 contigs for 3,446,046 bp with 56.17% G+C content and contains 3,360 protein-coding genes and 54 RNA genes.

Introduction

Acetic acid bacteria (AAB) are strictly aerobic . AAB are well known for their potential to incompletely oxidize a wide variety of sugars and alcohols. The genus oxidizes a wide range of sugars, sugar alcohols, and sugar acids, and can accumulate a large amount of the corresponding oxidized products in the culture medium [1]. Thus, strains are widely used for the industrial production of pharmaceutical intermediates, such as L-sorbose (vitamin C synthesis), 6-amino-L-sorbose (synthesis of the antidiabetic drug miglitol), and dihydroxyacetone (cosmetics for sunless tanning) [1]. Furthermore, the genera and are widely used for the industrial production of vinegar because of their high ethanol oxidation ability [2].

To date, six genome sequences of strains ( 621H, H24, WSH-003, NBRC 3255, NBRC 101659, and NBRC 103465) are available in the public databases [3-8]. These genomic data are useful for the experimental identification of unique proteins or estimation of the phylogenetic relationship among the related strains [9-11].

NBRC 3257 was isolated from downy cherry (Prunus tomentosa) in Japan [12], and identified based on its 16S rRNA sequence [13]. Here, we present a summary of the classification and a set of features of NBRC 3257, together with a description of the draft genome sequencing and annotation.

Classification and features

A representative genomic 16S rRNA sequence of NBRC 3257 was compared to the 16S rRNA sequences of all known species type strains. The 16S rRNA gene sequence identities between NBRC 3257 and all other type strains of genus species were 97.58-99.85%. species (type strains) exhibiting the highest sequence identities to NBRC 3257 were NBRC 3264T and NBRC 3271T. Figure 1 shows the phylogenetic relationships of NBRC 3257 to other species in a 16S rRNA based tree. All the type strains and ten strains of including NBRC 3257 were used for the analysis [13,17]. Based on this tree, genus is divided into two sub-groups. , , , G. japonicas, , and are classified as clade 1. On the other hand, and belong to the clade 2. All ten strains are closely related to each other, and the 16S rRNA sequences have 100% identities.

Figure 1

Phylogenetic tree highlighting the phylogenetic position of NBRC 3257 relative to other type strains within the . To construct the phylogenetic tree, these sequences were collected and nucleotide sequence alignment was carried out using CLUSTALW [14]. We used the MEGA version 5.05 package to generate phylogenetic trees based on 16S rRNA genes with the neighbor-joining (NJ) approach and 1,000 bootstrap replicates [15,16]. NBRC14818 (X74066) was used as the outgroup.

Although ethanol oxidation ability is a typical feature of AAB, it is a critical feature that NBRC 3257 lacks the ability to oxidize ethanol because it is missing the cytochrome subunit of the alcohol dehydrogenase complex that functions as the primary dehydrogenase in the ethanol oxidase respiratory chain [18]. Despite its inability to oxidize ethanol, NBRC 3257 can efficiently oxidize many unique sugars and sugar alcohols, such as pentitols, D-sorbitol, D-mannitol, glycerol, meso-erythritol, and 2,3-butanediol [19]. Thus, NBRC 3257 has unique characteristic features and the potential for the industrial production of many different oxidized products useful as drug intermediates or commodity chemicals.

G. thanilandicus NBRC 3257 is a strictly aerobic, mesophilic (temperature optimum ≈ 30˚C) organism. Differential interference contrast image of NBRC 3257 cells grown on mannitol medium (25 g of D-manntiol, 5 g of yeast extract, and 3 g of polypeptone per liter) are shown in Figure 2 (A). The cells have short-rod shape with 2.6 ± 0.6 (mean ± SD, n = 10) µm in cell length and 1.2 ± 0.1 (mean ± SD, n = 10) µm in cell width. The flagella stained by the modified Ryu method are shown in Figure 2 (B) and Figure 2 (C) [20]. Singly and multiply flagellated cells were observed frequently. The characteristic features are shown in Table 1.

Figure 2

Cell morphology and flagella of NBRC 3257. (A) Differential interference contrast image of NBRC 3257 grown on mannitol medium. Bar, 5 µm. (B and C) Microscopic images of flagella stained by the modified Ryu method. Singly (B) and multiply (C) flagellated cells were observed. Bars, 5 µm.

Table 1

Classification and general features of NBRC 3257 according to the MIGS recommendations [21]

MIGS ID	Property	Term	Evidence codea
	Current classification	Domain Bacteria	TAS [22]
		Phylum Proteobacteria	TAS [23]
		Class Alphaproteobacteria	TAS [24,25]
		Order Rhodospirillales	TAS [26,27]
		Family Acetobacteraceae	TAS [28,29]
		Genus Gluconobacter	TAS [27,30,31]
		Species Gluconobacter thailandicus	TAS [17,32]
		Strain NBRC 3257	TAS [23]
	Gram stain	Negative	IDA
	Cell shape	Rod-shaped	IDA
	Motility	Motile	NAS
	Sporulation	Not report	NAS
	Temperature range	Mesophilic	IDA
	Optimum temperature	30°C	IDA
	Carbon source	Glucose and/or glycerol	IDA
	Energy source	Glucose and/or glycerol	NAS
MIGS-6	Habitat	Free living	NAS
MIGS-6.3	Salinity	Not report	NAS
MIGS-22	Oxygen	Strict aerobes	IDA
MIGS-15	Biotic relationship	Fruits and Flower	IDA
MIGS-14	Pathogenicity	Non-pathogenic	IDA
MIGS-4	Geographic location	Japan	NAS
MIGS-5	Sample collection time	1954	NAS
MIGS-4.1	Latitude	Not report	NAS
MIGS-4.2	Longitude	Not report	NAS
MIGS-4.3	Depth	Not report	NAS
MIGS-4.4	Altitude	Not report	NAS

a Evidence codes - IDA: Inferred from Direct Assay; 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 the Gene Ontology project [33]. If the evidence code is IDA, then the property should have been directly observed, for the purpose of this specific publication, for a live isolate by one of the authors, or an expert or reputable institution mentioned in the acknowledgements.

Genome sequencing information

Genome project history

This genome was selected for sequencing on the basis of its phylogenetic position and 16S rRNA similarity to other members of the genus. This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession BASM00000000. The version described in this paper is the first version, BASM01000000, and the sequence consists of 107 contigs. Table 2 presents the project information and its association with MIGS version 2.0 compliance [34].

Table 2

Project information

MIGS ID	Property	Term
MIGS-31	Finishing quality	Draft
MIGS-28	Libraries used	Illumina Paired-End library
MIGS-29	Sequencing platforms	Illumina Hiseq 2000
MIGS-31.2	Fold coverage	358 ×
MIGS-30	Assemblers	Velvet ver. 1.2.07
MIGS-32	Gene calling method	Glimmer ver. 3.02
	DDBJ ID	BASM00000000
	DDBJ Date of Release	August 08, 2013
	Project relevance	Industrial

Growth conditions and DNA isolation

The culture of strain NBRC 3257 used to prepare genomic DNA for sequencing was a laboratory stock and grown on ∆P medium [35] at 30°C with vigorous shaking. The genomic DNA was isolated as described in [36] with some modifications [35]. Three ml of culture broth was used to isolate DNA, and the final DNA preparation was dissolved in 10 mM Tris-HCl (pH 8.0) and 1 mM ethylendiamine tetraacetic acid solution. The purity, quality, and size of the genomic DNA preparation were analyzed by Hokkaido System Science Co., Ltd. (Japan) using spectrophotometer, agarose gel electrophoresis, and Qubit (Invitrogen, Carlsbad, CA) according to the their guidelines.

Genome sequencing and assembly

The genome of NBRC 3257 was sequenced using the Illumina Hiseq 2000 sequencing platform by the paired-end strategy (2×100 bp). Paired-end genome fragments were annealed to the flow-cell surface in a cluster station (Illumina). A total of 100 cycles of sequencing-by-synthesis were performed and high-quality sequences were retained for further analysis. The final coverage reached 358-fold for an estimated genome size of 3.44 Mb. The sequence data from Illumina HiSeq 2000 were assembled with Velvet ver. 1.2.07 [37]. The final assembly yielded 107 contigs generating a genome size of 3.44 Mb. The contigs were ordered against the complete genome of 621H [3] using Mauve [38-40].

Genome annotation

Protein-coding genes (ORFs) of draft genome assemblies were predicted using Glimmer version 3.02 with a self-training dataset [41,42]. tRNAs and rRNAs were predicted using ARAGORN and RNAmmer, respectively [43,44]. Functional assignments of the predicted ORFs were based on a BLASTP homology search against two genome sequences, NBRC 3255 and 621H, and the NCBI nonredundant (NR) database [45]. Functional assignment was also performed with a BLASTP homology search against Clusters of Orthologous Groups (COG) databases [46].

Genome properties

The genome of NBRC 3257 is 3,446,046 bp long (107 contigs) with a 56.17% G + C content (Table 3). Of the 3,414 predicted genes, 3,360 were protein coding genes, and 54 were RNAs (3 rRNA genes, and 51 tRNA genes). A total of 2,249 genes (66.93%) were assigned a putative function. The remaining genes were annotated as hypothetical genes. The properties and statistics of the genome are summarized in Table 3. The distribution of genes into COG functional categories is presented in Table 4. Of the 3,360 proteins, 2,669 (79%) were assigned to COG functional categories. Of these, 245 proteins were assigned to multiple COG categories. The most abundant COG category was "General function prediction only" (342 proteins) followed by "Amino acid transport and metabolism" (247 proteins), "Function unknown" (232 proteins), "Cell wall/membrane/envelope biogenesis" (220 proteins), "Inorganic ion transport and metabolism" (210 proteins), and "Replication, recombination and repair" (201 proteins). The genome map of NBRC 3257 is illustrated in Figure 3, which demonstrates that the pattern of GC skew shifts from negative to positive along an ordered set of contigs with some exceptions. This suggests that the draft genome sequences were ordered almost exactly.

Table 3

Nucleotide content and gene count levels of the genome

Attribute	Value	% of totala
Genome size (bp)	3,446,046	-
DNA coding region (bp)	3,118,161	90.48
DNA G+C content (bp)	1,935,814	56.17
Total genesb	3,414	100
RNA genes	54	1.58
Protein-coding genes	3,360	98.42
Genes assigned to COGs	2,669	78.17

a) The total is based on either the size of the genome in base pairs or the total number of protein coding genes in the annotated genome.

Table 4

Number of genes associated with the 25 general COG functional categories

Code	Value	%agea	Description
J	157	4.67	Translation, ribosomal structure and biogenesis
A	0	0.00	RNA processing and modification
K	190	5.65	Transcription
L	201	5.98	Replication, recombination and repair
B	0	0.00	Chromatin structure and dynamics
D	28	0.83	Cell cycle control, cell division, chromosome partitioning
Y	0	0.00	Nuclear structure
V	44	1.31	Defense mechanisms
T	92	2.74	Signal transduction mechanisms
M	220	6.55	Cell wall/membrane/envelope biogenesis
N	43	1.28	Cell motility
Z	0	0.00	Cytoskeleton
W	2	0.06	Extracellular structures
U	95	2.83	Intracellular trafficking, secretion, and vesicular transport
O	120	3.57	Posttranslational modification, protein turnover, chaperones
C	170	5.06	Energy production and conversion
G	194	5.77	Carbohydrate transport and metabolism
E	247	7.35	Amino acid transport and metabolism
F	89	2.65	Nucleotide transport and metabolism
H	129	3.84	Coenzyme transport and metabolism
I	91	2.71	Lipid transport and metabolism
P	210	6.25	Inorganic ion transport and metabolism
Q	64	1.90	Secondary metabolites biosynthesis, transport and catabolism
R	342	10.18	General function prediction only
S	232	6.90	Function unknown
-	691	20.57	No COG assignment
-	245	7.29	Multiple COG assignment

a) The total is based on the total number of protein coding genes in the annotated genome.

Figure 3

Graphical circular map of a simulated draft NBRC 3257 genome. The simulated genome is a set of contigs ordered against the complete genome of 621H [3] using Mauve [38-40]. The circular map was generated using CGview [47]. From the outside to the center: genes on forward strand, genes on reverse strand, GC content, GC skew.

Gene repertoire of NBRC 3257 genome

Annotation of the genome indicated that NBRC 3257 has membrane-bound PQQ-dependent alcohol dehydrogenase, adhAB operon (locus_tag NBRC3257_1377 and NBRC3257_1378) and adh subunit III (NBRC3257_1024). A unique orphan gene of adh subunit I was also identified (NBRC3257_3117). The gene repertories of other membrane-bound PQQ dependent proteins were investigated. Homologous proteins of membrane-bound PQQ-dependent dehydrogenase (NBRC3257_0292), membrane-bound glucose dehydrogenase (PQQ) (NBRC3257_0371), PQQ-dependent dehydrogenase 4 (NBRC3257_0662), and PQQ-dependent dehydrogenase 3 (NBRC3257_1743), were identified. In addition, two paralogous copies of the PQQ-glycerol dehydrogenase sldAB operon (NBRC3257_0924 to NBRC3257_0925 and NBRC3257_1134 to NBRC3257_1135) were identified.

It has been thought that the respiratory chains of species play key roles in respiratory energy metabolism [48-51]. Therefore, the gene repertoires of respiratory chains of NBRC 3257 were also investigated. Besides two type II NADH dehydrogenase homologs (NBRC3257_1995 and NBRC3257_2785) [51], a proton-pumping NADH:ubiquinone oxidoreductase operon (type I NADH dehydrogenase complex) (NBRC3257_2617 to NBRC3257_2629), a cytochrome o ubiquinol oxidase cyoBACD operon (NBRC3257_2304 to NBRC3257_2307), and a cyanide-insensitive terminal oxidase cioAB operon (NBRC3257_0388 to NBRC3257_0389) [48,49], were identified.