# Genome sequence of the pink to light reddish-pigmented Rubellimicrobium mesophilum type strain (DSM 19309T), a representative of the Roseobacter group isolated from soil, and emended description of the species

Thomas Riedel, Stefan Spring, Anne Fiebig, Jörn Petersen, Markus Göker, Hans-Peter Klenk

► To cite this version:

Thomas Riedel, Stefan Spring, Anne Fiebig, Jörn Petersen, Markus Göker, et al.. Genome sequence of the pink to light reddish-pigmented Rubellimicrobium mesophilum type strain (DSM 19309T), a representative of the Roseobacter group isolated from soil, and emended description of the species. Standards in Genomic Sciences, 2014, 9 (3), pp.902-913. 10.4056/sigs.5621012. hal-01230502

HAL Id: hal-01230502

https://hal.sorbonne-universite.fr/hal-01230502

Submitted on 18 Nov 2015

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genome sequence of the pink to light reddish-pigmented *Rubellimicrobium mesophilum* type strain (DSM 19309<sup>T</sup>), a representative of the *Roseobacter* group isolated from soil, and emended description of the species

Thomas Riedel<sup>1,2</sup>, Stefan Spring<sup>3</sup>, Anne Fiebig<sup>3</sup>, Jörn Petersen<sup>3</sup>, Markus Göker<sup>3*</sup>, Hans-Peter Klenk<sup>3</sup>

**Keywords:** Irregular rod-shaped, motile, aerobic, stenohaline, chemoorganotroph, *Roseobacter* group, *Roseobacter* group, *Rhodobacteraceae*, *Alphaproteobacteria*.

*Rubellimicrobium mesophilum* Dastager et al. 2008 is a mesophilic and light reddish-pigmented representative of the *Roseobacter* group within the alphaproteobacterial family *Rhodobacteraceae*. Representatives of the *Roseobacter* group play an important role in the marine biogeochemical cycles and were found in a broad variety of marine environments associated with algal blooms, different kinds of sediments, and surfaces of invertebrates and vertebrates. Roseobacteria were shown to be widely distributed, especially within the total bacterial community found in coastal waters, as well as in mixed water layers of the open ocean. Here we describe the features of *R. mesophilum* strain MSL-20T together with its genome sequence and annotation generated from a culture of DSM 19309T. The 4,927,676 bp genome sequence consists of one chromosome and probably one extrachromosomal element. It contains 5,082 protein-coding genes and 56 RNA genes. As previously reported, the G+C content is significantly different from the actual genome sequence-based G+C content and as the type strain tests positively for oxidase, the species description is emended accordingly. The genome was sequenced as part of the activities of the Transregional Collaborative Research Centre 51 (TRR51) funded by the German Research Foundation (DFG).

## Introduction

Strain MSL-20T (=DSM 19309T = KCTC 22012T) is the type strain of the species *Rubellimicrobium mesophilum* [1], one of four species with validly published names in the genus *Rubellimicrobium* [2,3]; the other three species in the genus are *R. thermophilum* [3], *R. aerolatum* [4] and *R. roseum* [5]. *Rubellimicrobium* belongs to the abundant marine *Roseobacter* group [6]. The species epithet *mesophilum* refers to the Greek adjective *mesos*, middle, as well as from the Neo-Latin adjective ‘philus -a -um’, friend/loving [1], the middle (temperature-) loving. Strain MSL-20T was isolated from soil located at Bigeum Island, Republic of Korea [1], whereas the other type strains within the genus *Rubellimicrobium* were isolated from a

paper mill (*R. thermophilum* [3]), air (*R. aerolatum* [4]) and forest soil (*R. roseum* [5]), which indicates rather diverse habitats for *Rubellimicrobium*. Current PubMed records do not indicate any follow-up research with strain MSL-20T since the initial description of *R. mesophilum* [1]. Here we present a summary classification and a set of features for *R. mesophilum* MSL-20T, together with the description of the complete genomic sequencing and annotation.

## Classification and features
16S rRNA gene analysis

Figure 1 shows the phylogenetic neighborhood of *R. mesophilum* in a 16S rRNA gene sequence-based

<figure><img src="image_2.png" /></figure>

<footnote n="1"><sup>1</sup> Sorbonne Universités, UPMC Univ Paris 06, USR3579, LBBM, Observatoire Océanologique, Banyuls/Mer, France</footnote>

<footnote n="2"><sup>2</sup> CNRS, USR3579, LBBM, Observatoire Océanologique, Banyuls/Mer, France</footnote>

<footnote n="3"><sup>3</sup> Leibniz Institute DSMZ - German Collection of Microorganisms and Cell Cultures, Braunschweig, Germany</footnote>

<footnote>*Correspondence: Markus Göker (mgo08@dsmz.de)</footnote>
tree. The sequence of the single 16S rRNA gene in the DSM 19309T genome does not differ from the previously published 16S rRNA gene sequence (EF547368), which contains four ambiguous base calls.

The genomic 16S rRNA gene sequence of *R. mesophilum* DSM 19309T was compared with the Greengenes database for determining the weighted relative frequencies of taxa and (truncated) keywords as previously described [7]. The most frequently occurring genera were *Paracoccus* (45.3%), *Loktanella* (30.3%), *Rubellimicrobium* (14.0%), *Methylarcula* (8.4%) and *'Pararubellimicrobium'* (2.0%) (58 hits in total). Regarding the five hits to sequences from other members of the genus, the average identity within HSPs was 94.9%, whereas the average coverage by HSPs was 99.3%. Among all other species, the one yielding the highest score was

*'Pararubellimicrobium aerilata'* (EU338486), which corresponded to an identity of 96.3% and a HSP coverage of 98.0%. (Note that the Greengenes database uses the INSDC (=EMBL/NCBI/DDBJ) annotation, which is not an authoritative source for nomenclature or classification). The highest-scoring environmental sequence was JF417792 (Greengenes short name 'microbial structures coalbeds located Eerduosi Basin China coalbed clone QQSB73'), which showed an identity of 98.7% and a HSP coverage of 99.6%. The most frequently occurring keywords within the labels of all environmental samples which yielded hits were 'skin' (10.6%), 'fossa' (5.9%), 'poplit' (4.2%), 'forearm, volar' (3.3%) and 'sea' (2.8%) (192 hits in total). Environmental samples which yielded hits of a higher score than the highest scoring species were not found, indicating that *R. mesophilum* has rarely been detected in the environment.

<figure><img src="image_3.png" /></figure>

**Figure 1.** Phylogenetic tree highlighting the position of *R. mesophilum* relative to the type strains of the other species within the genus *Rubellimicrobium* and the neighboring genera *Citreicella* and *Wenxinia*. The tree was inferred from 1,381 aligned characters of the 16S rRNA gene sequences under the maximum likelihood (ML) criterion as previously described [7]. The branches are scaled in terms of the expected number of substitutions per site. Numbers adjacent to the branches are support values from 1,000 ML bootstrap replicates (left) and from 1,000 maximum-parsimony bootstrap replicates (right) if larger than 60% [7]. Lineages with type strain genome sequencing projects registered in GOLD [8] are labeled with one asterisk [9].
**Morphology and physiology**

Cells of strain MSL-20T stain Gram-negative, are described to be motile (without a flagellum) [1], and ovoid or rod-shaped, 1.6-3.4 µm in length and 0.4-0.7 µm in width (Figure 2 and Table 1). On Reasoner's 2A (R2A) agar they form pink to light red-pigmented colonies. According to [1], cells are negative for oxidase (but see below) and nitrate reduction activities, but show only weak catalase activity. They hydrolyze starch and Tween 80, assimilate cellulose, histidine, leucine and fructose, but do not utilize citrate and propionate. Cells test

positive for leucine arylamidase, naphthol-AS-BI-phosphohydrolase and α-glucosidase. Growth is observed in a temperature range of 20-37°C with an optimum at 28°C. The pH range for growth is between pH 7-11 with an optimum at pH 7.0 ± 0.2. No growth occurs in the presence of NaCl in concentrations of 0.5% and above. Cells of strain MSL-20T do not utilize the carbohydrates cellobiose, D-mannose, salicin, D-xylose, α-melibiose, D-sorbitol, L-malate and D-ribose, which are utilized by its close relative *R. thermophilum* DSM 16684T (all data from [1]).

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Micrograph of *R. mesophilum* DSM 19309T.</figcaption></figure>

**Chemotaxonomy**

The principal cellular fatty acids of strain MSL-20T are C16:0 (36.9%), C18:1ω7c (36.5%), 11-methyl C18:1ω7c (12.4%), C18:0 (3.6%), C10:0 (1.3%), C12:0 (1.3%) and C17:0 (1.2%) and differ significantly from those detected in *R. thermophilum*. The major respiratory lipoquinone is ubiquinone Q-10, which is a common feature of alphaproteobacterial representatives (all data from [1]).

**Genome sequencing and annotation**
**Genome project history**

The genome of strain *R. mesophilum* DSM 19309T was first selected for genome sequencing in phase I of the one thousand microbial genomes (KMG-I) project [20], an extension of the Genomic Ency-

clopaedia of *Bacteria* and *Archaea* (GEBA) [21], but ultimately sequenced within the DFG funded project "Ecology, Physiology and Molecular Biology of the *Roseobacter* clade: Towards a Systems Biology Understanding of a globally Important Clade of Marine Bacteria". The strain was chosen for genome sequencing according to a phylogeny-driven target selection procedure for large scale genome-sequencing (and other) projects as routinely used for the KMG-I project [20,22].

The project information can be found in the Genome OnLine Database [8]. The Whole Genome Shotgun (WGS) sequence is deposited in GenBank and the Integrated Microbial Genomes database (IMG) [23]. A summary of the project information is shown in Table 2.
<table><caption>**Table 1.** Classification and general features of *R. mesophilum* MSL-20T according the MIGS recommendations [10] published by the Genome Standards Consortium [11].</caption><thead><tr><th>MIGS ID</th><th>Property</th><th>Term</th><th>Evidence code</th></tr></thead><tbody><tr><td></td><td></td><td>Domain Bacteria</td><td>TAS [12]</td></tr><tr><td></td><td></td><td>Phylum Proteobacteria</td><td>TAS [13]</td></tr><tr><td></td><td></td><td>Class Alphaproteobacteria</td><td>TAS [14,15]</td></tr><tr><td></td><td></td><td>Order Rhodobacterales</td><td>TAS [15,16]</td></tr><tr><td></td><td></td><td>Family Rhodobacteraceae</td><td>TAS [15,17]</td></tr><tr><td></td><td></td><td>Genus Rubellimicrobium</td><td>TAS [3]</td></tr><tr><td></td><td></td><td>Species Rubellimicrobium mesophilum</td><td>TAS [1]</td></tr><tr><td></td><td>Current classification</td><td>Strain MSL-20T</td><td>TAS [1]</td></tr><tr><td></td><td>Gram stain</td><td>negative</td><td>TAS [1]</td></tr><tr><td></td><td>Cell shape</td><td>irregular rod-shaped</td><td>TAS [1]</td></tr><tr><td></td><td>Motility</td><td>motile</td><td>TAS [1]</td></tr><tr><td></td><td>Sporulation</td><td>non-sporulating</td><td>NAS</td></tr><tr><td></td><td>Temperature range</td><td>20-37°C</td><td>TAS [1]</td></tr><tr><td></td><td>Optimum temperature</td><td>28°C</td><td>TAS [1]</td></tr><tr><td></td><td>Salinity</td><td>stenohaline</td><td>TAS [1]</td></tr><tr><td>MIGS-22</td><td>Oxygen requirement</td><td>aerobic</td><td>TAS [1]</td></tr><tr><td></td><td>Carbon source</td><td>carbohydrates, amino acids</td><td>TAS [1]</td></tr><tr><td></td><td>Energy metabolism</td><td>chemoorganotroph</td><td>NAS</td></tr><tr><td>MIGS-6</td><td>Habitat</td><td>soil</td><td>TAS [1]</td></tr><tr><td>MIGS-15</td><td>Biotic relationship</td><td>free living</td><td>TAS [1]</td></tr><tr><td>MIGS-14</td><td>Pathogenicity</td><td>none</td><td>NAS</td></tr><tr><td></td><td>Biosafety level</td><td>1</td><td>TAS [18]</td></tr><tr><td>MIGS-231</td><td>Isolation</td><td>soil</td><td>TAS [1]</td></tr><tr><td>MIGS-4</td><td>Geographic location</td><td>Bigeum island (Republic of Korea)</td><td>TAS [1]</td></tr><tr><td>MIGS-5</td><td>Sample collection time</td><td>April 2006</td><td>NAS</td></tr><tr><td>MIGS-4.1</td><td>Latitude</td><td>34.739</td><td>NAS</td></tr><tr><td>MIGS-4.2</td><td>Longitude</td><td>125.920</td><td>NAS</td></tr><tr><td>MIGS-4.3</td><td>Depth</td><td>not reported</td><td>NAS</td></tr><tr><td>MIGS-4.4</td><td>Altitude</td><td>not reported</td><td>NAS</td></tr></tbody></table>

Evidence codes - 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). Evidence codes are from of the Gene Ontology project [19]

.
<table><caption>**Table 2.** Genome sequencing project information</caption><thead><tr><th>MIGS ID</th><th>Property</th><th>Term</th></tr></thead><tbody><tr><td>MIGS-31</td><td>Finishing quality</td><td>Non-contiguous finished</td></tr><tr><td>MIGS-28</td><td>Libraries used</td><td>Two genomic libraries: one Illumina PE library (420 bp insert size), one 454 PE library (3 kb insert size)</td></tr><tr><td>MIGS-29</td><td>Sequencing platforms</td><td>Illumina GA IIx, Illumina MiSeq, 454 GS-FLX+Titanium</td></tr><tr><td>MIGS-31.2</td><td>Sequencing coverage</td><td>129x</td></tr><tr><td>MIGS-30</td><td>Assemblers</td><td>Velvet version 1.1.36, Newbler version 2.3, Consed 20.0</td></tr><tr><td>MIGS-32</td><td>Gene calling method</td><td>Prodigal 1.4</td></tr><tr><td></td><td>INSDC ID</td><td>AOSK00000000</td></tr><tr><td></td><td>GenBank Date of Release</td><td>pending publication</td></tr><tr><td></td><td>GOLD ID</td><td>Gi0042374</td></tr><tr><td></td><td>NCBI project ID</td><td>188767</td></tr><tr><td></td><td>Database: IMG</td><td>2523533591</td></tr><tr><td>MIGS-13</td><td>Source material identifier</td><td>DSM 19309<sup>T</sup></td></tr><tr><td></td><td>Project relevance</td><td>Tree of Life, biodiversity</td></tr></tbody></table>

## Growth conditions and DNA isolation

A culture of DSM 19309<sup>T</sup> was grown aerobically in DSMZ medium 830 (R2A medium) [24] at 28°C. Genomic DNA was isolated using Jetflex Genomic DNA Purification Kit (GENOMED 600100) following the standard protocol provided by the manufacturer, but modified by an incubation time of 60 min, an overnight incubation on ice on a shaker, the use of additional 50 µl proteinase K, and the addition of 100 µl protein precipitation buffer. DNA is available from DSMZ through the DNA Bank Network [25].

## Genome sequencing and assembly

The genome was sequenced using a combination of two libraries (Table 2). The paired-end library contained inserts of an average of 420 bp in length. Illumina sequencing was performed on a GA IIx platform with 150 cycles. The first run on the Illumina GA IIx platform delivered 3.6 million reads. In order to increase the sequencing depth, a second Illumina run was performed, providing another 7.0 million reads. Error correction and clipping were performed by fastq-mcf [26] and quake [27]. The data was assembled using Velvet [28]. The first draft assembly from 5,400,234 filtered reads (median read length of 132 nt) resulted in more than 143 unordered contigs. To gain information about the contig arrangement an additional 454 run was performed. The paired-end

jumping library of 3 kb insert size was sequenced on 1/8 of a lane. Pyrosequencing resulted in 102,695 reads with an average read length of 199 bp, assembled with Newbler (Roche Diagnostics). The resulting assembly consisted of 261 scaffolds. Both draft assemblies (Illumina and 454 sequenc- es) were fractionated into artificial Sanger reads of 1,000 nt in length plus 75 bp overlap on each site. These artificial reads served as an input for the phred/phrap/consed package [29]. By manual editing, 138 contigs could be assembled on 127 scaffolds. The combined sequences provided a 129x coverage of the genome.

## Genome annotation

Genes were identified using Prodigal [30] as part of the JGI genome annotation pipeline. The predicted CDSs were translated and used to search the National Center for Biotechnology Information (NCBI) non-redundant database, UniProt, TIGR-Fam, Pfam, PRIAM, KEGG, COG, and InterPro databases. Identification of RNA genes was carried out by using HMMER 3.0rc1 [31] (rRNAs) and tRNAscan-SE 1.23 [32] (tRNAs). Other non-coding genes were predicted using INFERNAL 1.0.2 [33]. Additional gene prediction analysis and functional annotation was performed within the Integrated Microbial Genomes - Expert Review (IMG-ER) platform [34]. CRISPR elements were detected using CRT [35] and PILER-CR [36].
## Genome properties

The genome statistics are provided in Table 3 and Figure 3. The genome of strain DSM 19309T has a total length of 4,927,676 bp and a G+C content of 69.7%. Of the 5,138 genes predicted, 5,082 were identified as protein-coding genes, and 56 as

RNAs. The majority of the protein-coding genes (56.7%) 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.

<table><caption>**Table 3.** Genome statistics</caption><thead><tr><td>Attribute</td><td>Value</td><td>% of Total</td></tr></thead><tbody><tr><td>Genome size (bp)</td><td>4,927,676</td><td>100.00</td></tr><tr><td>DNA coding region (bp)</td><td>4,254,404</td><td>86.34</td></tr><tr><td>DNA G+C content (bp)</td><td>3,431,981</td><td>69.65</td></tr><tr><td>Number of scaffolds MIGS-9</td><td>127</td><td></td></tr><tr><td>Extrachromosomal elements MIGS-10</td><td>1</td><td></td></tr><tr><td>Total genes</td><td>5,138</td><td>100.00</td></tr><tr><td>RNA genes</td><td>56</td><td>1.09</td></tr><tr><td>rRNA operons</td><td>1</td><td></td></tr><tr><td>tRNA genes</td><td>45</td><td>0.88</td></tr><tr><td>Protein-coding genes</td><td>2,915</td><td>56.73</td></tr><tr><td>Genes with function prediction (proteins)</td><td>2,167</td><td>42.18</td></tr><tr><td>Genes in paralog clusters</td><td>4,172</td><td>81.20</td></tr><tr><td>Genes assigned to COGs</td><td>3,818</td><td>74.31</td></tr><tr><td>Genes assigned Pfam domains</td><td>3,977</td><td>77.40</td></tr><tr><td>Genes with signal peptides</td><td>384</td><td>7.47</td></tr><tr><td>Genes with transmembrane helices</td><td>966</td><td>18.80</td></tr><tr><td>CRISPR repeats</td><td>0</td><td></td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Graphical map of the largest, 267,932 bp long scaffold. From bottom to the top: Genes on forward strand (colored by COG categories), Genes on reverse strand (colored by COG categories), RNA genes (tRNAs green), GC content (black), GC skew (purple/olive).</figcaption></figure>
<table><caption>**Table 4.** Number of genes associated with the general COG functional categories</caption><thead><tr><th>Code</th><th>Value</th><th>%age</th><th>Description</th></tr></thead><tbody><tr><td>J</td><td>186</td><td>4.4</td><td>Translation, ribosomal structure and biogenesis</td></tr><tr><td>A</td><td>3</td><td>0.1</td><td>RNA processing and modification</td></tr><tr><td>K</td><td>279</td><td>6.6</td><td>Transcription</td></tr><tr><td>L</td><td>250</td><td>5.9</td><td>Replication, recombination and repair</td></tr><tr><td>B</td><td>4</td><td>0.1</td><td>Chromatin structure and dynamics</td></tr><tr><td>D</td><td>35</td><td>0.8</td><td>Cell cycle control, cell division, chromosome partitioning</td></tr><tr><td>Y</td><td>0</td><td>0.0</td><td>Nuclear structure</td></tr><tr><td>V</td><td>34</td><td>0.8</td><td>Defense mechanisms</td></tr><tr><td>T</td><td>176</td><td>4.2</td><td>Signal transduction mechanisms</td></tr><tr><td>M</td><td>223</td><td>5.3</td><td>Cell wall/membrane/envelope biogenesis</td></tr><tr><td>N</td><td>27*</td><td>0.6</td><td>Cell motility</td></tr><tr><td>Z</td><td>0</td><td>0.0</td><td>Cytoskeleton</td></tr><tr><td>W</td><td>0</td><td>0.0</td><td>Extracellular structures</td></tr><tr><td>U</td><td>51</td><td>1.2</td><td>Intracellular trafficking and secretion, and vesicular transport</td></tr><tr><td>O</td><td>148</td><td>3.5</td><td>Posttranslational modification, protein turnover, chaperones</td></tr><tr><td>C</td><td>251</td><td>6.0</td><td>Energy production and conversion</td></tr><tr><td>G</td><td>453</td><td>10.7</td><td>Carbohydrate transport and metabolism</td></tr><tr><td>E</td><td>478</td><td>11.3</td><td>Amino acid transport and metabolism</td></tr><tr><td>F</td><td>94</td><td>2.2</td><td>Nucleotide transport and metabolism</td></tr><tr><td>H</td><td>152</td><td>3.6</td><td>Coenzyme transport and metabolism</td></tr><tr><td>I</td><td>143</td><td>3.4</td><td>Lipid transport and metabolism</td></tr><tr><td>P</td><td>182</td><td>4.3</td><td>Inorganic ion transport and metabolism</td></tr><tr><td>Q</td><td>124</td><td>2.9</td><td>Secondary metabolites biosynthesis, transport and catabolism</td></tr><tr><td>R</td><td>508</td><td>12.0</td><td>General function prediction only</td></tr><tr><td>S</td><td>421</td><td>10.0</td><td>Function unknown</td></tr><tr><td>-</td><td>1,320</td><td>25.7</td><td>Not in COGs</td></tr></tbody></table>

## Insights into the genome

### Plasmids

The identification of plasmids is difficult because typical replication modules comprising the characteristic replicase and the adjacent *parAB* partitioning operon are missing [36]. However, comprehensive BLASTP searches with plasmid replicases from *Rhodobacterales* revealed the presence of one RepB gene (rumeso_01479), whereas RepA-, RepABC-type and DnaA-like

replicases are absent from the genome. The localization of the chromosomal replication initiator DnaA documents that scaffold 15 is part of the chromosome (Table 5).

The 119 kb RepB type plasmid contains a post-segregational killing system (PSK) consisting of a typical operon with two small genes encoding a stable toxin and an unstable antitoxin (rumesco_01477/78 [37];).

<table><caption>**Table 5.** General genomic location and features of the chromosomal and one extrachromosomal replicon from *R. mesophilum* strain DSM 19309T.</caption><thead><tr><th>Replicon</th><th>Scaffold</th><th>Replicase</th><th>Length (bp)</th><th>GC (%)</th><th>Topology</th><th>No. Genes*</th></tr></thead><tbody><tr><td>Chromosome<fnref n="1" /></td><td>15</td><td>DnaA</td><td>102,082</td><td>71</td><td>linear*</td><td>105</td></tr><tr><td>Plasmid</td><td>9</td><td>RepB</td><td>119,205</td><td>68</td><td>linear*</td><td>141</td></tr></tbody></table>

<footnote>*Only one gene each for flagellar motor and flagellar hook capping, no structural genes for flagella.</footnote>

<footnote>*circularity not experimentally validated</footnote>

<footnote>†deduced from automatic annotation</footnote>

<footnote n="1">¹partial sequence including the replicase *dnaA* (rumeso_02152).</footnote>
## Phages

Phages are widely distributed and abundant in marine and freshwater environments [38-40] and are known to be horizontal gene transfer agents that drive bacterial diversity [40,41]. Temperate phage genomes can be integrated in the host genome as prophages and perform a symbiotic relationship with their hosts [42].

Several phage-associated gene sequences were detected in the genome sequence of strain DSM 19309T, particularly in “genomic islands” (e.g., rumeso_00405, rumeso_00407 rumeso_01586 to rumeso_01600).

## Quorum Sensing

Several Gram-negative bacteria produce and release chemical signal molecules called autoinducers. In correlation to the population density they detect those signal molecules and respond with an alteration of gene expression and therefore with diverse behaviors (e.g., luminescence, virulence, antibiotic resistance, changes in morphology and cell division) [43-46].

Genome analysis of strain DSM 19309T revealed the presence of gene-encoding sequences associated with the mechanism of quorum sensing e.g. N-homoserine-lactone synthetase, rumeso_02218 (LuxI homologue); probably involved in response and transcriptional regulators, rumeso_02217 (luxR homologue).

## Metabolic plasticity

Unlike many representatives of the Roseobacter group [6], R. mesophilum DSM 19309T encodes no genes involved in the harvesting of light and photoheterotrophic growth, which reflect its occurrence in niches within soil that are characterized by the absence of light. Nevertheless, the annotated genome sequence reveals a high metabolic versatility that was not expected by the phenotypic characterization presented in the species description [1].

The genome encodes a large number of diverse ABC transporters facilitating the uptake of various substrates like carbohydrates (e.g., rumeso_04497 to 04500), polyamines (e.g., rumeso_04716 to 04719), peptides (e.g., rumeso_00087 to 00090), amino acids (e.g., rumeso_00231 to 00234) and sulfonates (e.g., rumeso_05058 to 05059). Sulfonates could represent unexpected but common substrates for this species. The organic sulfonates taurine and cysteic acid are widely dis-

tributed in animal tissue and can enter soil by fces. In some soil bacteria, these compounds are used as sole source of carbon, nitrogen and sulfur [47]. Indeed, a complete degradation pathway for taurine was detected in the genome of strain DSM 19309T. Taurine is first converted by a taurine-pyruvate aminotransferase (rumeso_05057) to sulfoacetaldehyde, which in turn is cleaved by the enzyme sulfoacetaldehyde acetyltransferase (rumeso_03970) into sulfite and acetyl-phosphate. Acetyl-phosphate can be either converted to acetyl-CoA by a phosphotransacetylase (rumeso_03968) and funneled into the intermediary metabolism or is used for the generation of ATP by the enzyme acetate kinase (rumeso_03967). The potentially toxic compound sulfite can be oxidized to sulfate by various sulfite oxidases (e.g., rumeso_03951).

In addition, the utilization of electron acceptors seems to be variable and not restricted to oxygen. Genes encoding at least two predicted cytochrome c oxidases, one of the cbb₃-type (rumeso_00470 to 00472) and the other of the aa₃-type (rumeso_02204 to 02206), which terminate the electron transport chain with oxygen, were detected. However, according to the species description strain MSL-20T should be oxidase negative [1], we have found that the oxidase test for this strain is positive, which is in line with the results of the genome analysis.

Under periodic anoxic conditions that frequently occur in wet soils, nitrate could be used as alternative electron acceptor. According to the genome sequence, the denitrification pathway of this strain is probably incomplete and terminates with the greenhouse gas nitrous oxide (N₂O), as has been previously demonstrated for Ottowia thiooxydans [48]. Only genes encoding a respiratory nitrate reductase (rumeso_02471 to 02474), nitrite reductase (rumeso_02669) and nitric oxide reductase (rumeso_00142 to 00145) were detected, whereas no genes for the terminal nitrous oxide reductase were found.

## Comparison of *Rubellimicrobium* genomes

Recently the genome sequence of the type strain for second representative of the genus Rubellimicrobium, R. thermophilum DSM 16684T became available [9]. Lifestyle, habitat and preferred temperature range of R. thermophilum differ significantly from the ones of R. mesophilum [3]. The genome sequences of both strains were compared using the digital DNA-DNA hybridiza-
tion (dDDH) tool GGDC server version 2.0, an online tool provided through the DSMZ web pages [49]. The resulting dDDH value of 19.3 ± 2.3% according to distance formula 2 (as described in [50]), confirmed that both strains belong to independent species.

Figure 4 depicts the fraction of shared genes between the two genome-sequenced *Rubellimicrobium* type strains and the type strain of *Wenxinia marina* [51], another closely related member of the *Roseobacter* group (see Figure 1).

The number of pairwise genes was inferred from the phylogenetic profiler tool of the IMG platform. Homologous genes were detected with an E-value cutoff of 10⁻⁵ and a minimum identity of 30%. Proportions of 56% and 45% of the gene count in *W. marina* and *R. mesophilum*, respectively, are shared between all three genomes. In the case of *R. thermophilum*, a fraction of homologous genes of 70% is present in the other two genomes. Very few genes are shared only between *R. thermophilum* and *W. marina*.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Venn diagram depicting the intersections of proteins sets (total numbers in parentheses) of the two *Rubellimicrobium* species and *W. marina*.</figcaption></figure>

Although both genomes differ significantly in size (3.2 Mbp for *R. mesophilum* and 4.9 Mbp for *R. thermophilum*), the proportions of genes per COG category is very similar (Table 3 and [9]). The IMG Abundance Profile [34] demonstrated some differences, however. Enzymes for transport and utilization of amino acids and polyamines (COG1173, COG0747, COG3842) were present in higher abundance in *R. thermophilum*, which is in agreement with the results from wet-lab substrate tests [1,3]. Huge differences in the abundance of proteins can be found within the class of transposases (COG2801, COG 3436, COG2936, COG0665, COG0404). While *R. thermophilum* codes for two transposases, more than 30 tranposase genes were identified in *R. mesophilum*. Combined with the presence of the site-specific recombinase XerD (involved in the recombination of plasmids [20]) this indicates a high level of genetic recombination within *R. mesophilum*. Furthermore, 23 genes coding for RTX toxins and Ca⁺-binding proteins (COG 2931) were found. These proteins are structurally diverse, playing an important role in the colonization of various habitats and surfaces [50].

Additively, 14 proteins of the xenobiotic-degrading glutathion-S-transferases were present in *R. mesophilum*. The occurrence of these proteins may enable the bacteria to grow in polluted areas.

**Taxonomic note**

The G +C content of the genomic DNA of strain MSL-20T is given in the species description as 72.3 mol% [1], which represents a discrepancy of more than 2% from the value of 69.7 mol% deduced from the genome sequence. In addition to the deviant oxidase test this calls for an emendation of the species description according to the proposal of Meier-Kolthoff et al. [47].

**Emended description of *Rubellimicrobium mesophilum* Dastager et al. 2008**

The description of the species *Rubellimicrobium mesophilum* is the one given by Dastager *et al.* 2008 [1], with the following modifications.

Oxidase test is positive. The G+C content, rounded to zero decimal places, is 70%.
Acknowledgements

The authors gratefully acknowledge Evelyne Brambilla (DSMZ) for DNA extraction and quality control. This work was performed under the aus-

pices of the German Research Foundation (DFG)
Transregio-SFB 51 *Roseobacter* grant.

References

1. Dastager SG, Lee JC, Ju YJ, Park DJ, Kim CJ. *Rubellimicrobium mesophilum* sp. nov., a mesophilic, pigmented bacterium isolated from soil. *Int J Syst Evol Microbiol* 2008; **58**:1797-1800. PubMed http://dx.doi.org/10.1099/ijs.0.65590-0

2. Euzeby JP. List of prokaryotic names with standing in nomenclature. www.bacterio.net/qr/rubellimicrobium.html

3. Denner EB, Kolari M, Hoornstra D, Tsitko I, Kämpfer P, Busse HJ, Salkinoja-Salonen M. *Rubellimicrobium thermophilum* gen. nov., sp. nov., a red-pigmented, moderately thermophilic bacterium isolated from coloured slime deposits in paper machines. *Int J Syst Evol Microbiol* 2006; **56**:1355-1362. PubMed http://dx.doi.org/10.1099/ijs.0.63751-0

4. Weon HY, Son JA, Yoo SH, Hong SB, Jeon YA, Kwon SW, Koo BS. *Rubellimicrobium aerolatum* sp. nov., isolated from an air sample in Korea. *Int J Syst Evol Microbiol* 2009; **59**:406-410. PubMed http://dx.doi.org/10.1099/ijs.0.65856-0

5. Cao YR, Jiang Y, Wang Q, Tang SK, He WX, Xue QH, Xu LH, Jiang CL. *Rubellimicrobium roseum* sp. nov., a Gram-negative bacterium isolated from the forest soil sample. *Antonie van Leeuwenhoek* 2010; **98**:389-394. PubMed http://dx.doi.org/10.1007/s10482-010-9452-2

6. Buchan A, González JM, Moran MA. Overview of the marine *Roseobacter* lineage. *Appl Environ Microbiol* 2005; **71**:5665-5677. PubMed http://dx.doi.org/10.1128/AEM.71.10.5665-5677.2005

7. Göker M, Cleland D, Saunders E, Lapidus A, No- lan M, Lucas S, Hammon N, Deshpande S, Cheng JF, Tapia R, et al. Complete genome sequence of *Isosphaera pallida* type strain (IS1B<sup>T</sup>). *Stand Ge- nomic Sci* 2011; **4**:63-71. PubMed http://dx.doi.org/10.4056/sigs.1533840

8. Pagani I, Liolios K, Jansson J, Chen IM, Smirnova T, Nosrat B, Markowitz VM, Kyrpides NC. The Genomes OnLine Database (GOLD) v.4: status of genomic and metagenomic projects and their as- sociated metadata. *Nucleic Acids Res* 2012; **40**:D571-D579. PubMed http://dx.doi.org/10.1093/nar/gkr1100

9. Fiebig A, Riedel T, Gronow S, Petersen J, Klenk HP, Göker M. Genome sequence of the reddish-pigmented *Rubellimicrobium thermophilum* type strain (DSM 16684<sup>T</sup>), a member *Roseobacter* clade. *Stand Genomic Sci* 2013; **7**:107-119.

10. Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV, et al. The minimum information about a genome sequence (MIGS) specification. *Nat Biotechnol* 2008; **26**:541-547. PubMed http://dx.doi.org/10.1038/nbt1360

11. Field D, Amaral-Zettler L, Cochrane G, Cole JR, Dawyndt P, Garrity GM, Gilbert J, Glöckner FO, Hirschman L, Karsch-Mzrachi I, et al. Clarifying Concepts and Terms in Biodiversity Informatics. *PLoS Biol* 2013; **9**:e1001088. PubMed http://dx.doi.org/10.1371/journal.pbio.1001088

12. Woese CR, Kandler O, Weelis ML. Towards a natural system of organisms. Proposal for the domains *Archaea*, *Bacteria* and *Eucarya*. *Proc Natl Acad Sci USA* 1990; **87**:4576-4579. PubMed http://dx.doi.org/10.1073/pnas.87.12.4576

13. Garrity GM, Bell JA, Lilburn T. Phylum XIV. *Proteobacteria* phyl nov. *In*: Brenner DJ, Krieg NR, Stanley JT, Garrity GM (*eds*), Bergey's Manual of Systematic Bacteriology, second edition. Vol. 2 (The *Proteobacteria*), part B (The *Gammaproteobacteria*), Springer, New York, 2005, p. 1.

14. Garrity GM, Bell JA, Lilburn T. Class I. *Alphaproteobacteria* class. nov. *In*: Brenner DJ, Krieg NR, Stanley JT, Garrity GM (*eds*), Bergey's Manual of Systematic Bacteriology, second edition. Vol. 2 (The *Proteobacteria*), part C (The *Alphaproteobacteria*), Springer, New York, 2005, p. 1.

15. Validation List No. 107. List of new names and new combinations previously effectively, but not validly, published. *Int J Syst Evol Microbiol* 2006; **56**:1-6. PubMed http://dx.doi.org/10.1099/ijs.0.64188-0

16. Garrity GM, Bell JA, Lilburn T. Order III. *Rhodobacterales* ord. nov. *In*: Brenner DJ, Krieg NR, Staley JT, Garrity GM (*eds*), Bergey's Manual of Systematic Bacteriology, second edition. vol. 2 (The *Proteobacteria*), part C (The *Alpha*-, *Beta*-,
*Delta-*, and *Epsilonproteobacteria*), Springer, New York, 2005, p. 161.

17. Garrity GM, Bell JA, Lilburn T. Family I. *Rhodobacteraceae* fam. nov. *In*: Brenner DJ, Krieg NR, Staley JT, Garrity GM (*eds*), Bergey's Manual of Systematic Bacteriology, second edition. vol. 2 (The *Proteobacteria*), part C (The *Alpha-*, *Beta-*, *Delta-*, and *Epsilonproteobacteria*), Springer, New York, 2005, p. 161.

18. BAuA. Classification of *Bacteria* and *Archaea* in risk groups. *TRBA* 2010; **466**:93.

19. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. *Nat Genet* 2000; **25**:25-29. PubMed http://dx.doi.org/10.1038/75556

20. Kyrpides NC, Woyke T, Eisen JA, Garrity GM, Whitman WB, Klenk HP. Genomic encyclopedia of type strains, phase I: the one thousand microbial genomes (KNG-I) project. *Stand Genomic Sci* 2013; **9**:628-634. http://dx.doi.org/10.4056/sigs.5068949

21. Wu D, Hugenholtz P, Mavromatis K, Pukall R, Dalin E, Ivanova NN, Kunin V, Goodwin L, Wu M, Tindall BJ, et al. A phylogeny-driven Genomic Encyclopaedia of *Bacteria* and *Archaea*. *Nature* 2009; **462**:1056-1060. PubMed http://dx.doi.org/10.1038/nature08656

22. Göker M, Klenk HP. Phylogeny-driven target selection for large scale genome-sequencing (and other) projects. *Stand Genomic Sci* 2013; **8**:360-374. PubMed http://dx.doi.org/10.4056/sigs.3446951

23. Markowitz VM, Chen IMA, Palaniappan K, Chu K, Szeto E, Grechkin Y, Ratner A, Jacob B, Huang J, Williams P, et al. IMG: the integrated microbial genomes database and comparative analysis system. *Nucleic Acids Res* 2012; **40**:D115-D122. PubMed http://dx.doi.org/10.1093/nar/gkr1044

24. List of growth media used at the DSMZ: http://www.dmsz.de/catalogues/catalogo-microorganisms/culture-technology/list-of-media-for-microorganisms.html.

25. Gemeinholzer B, Dröge G, Zetzsche H, Haszprunar G, Klenk HP, Güntsch A, Berendsohn WG, Wägele JW. The DNA Bank Network: the start from a German initiative. *Biopreserv Biobank* 2011; **9**:51-55. http://dx.doi.org/10.1089/bio.2010.0029

26. Aronesty E. *ea-utils*: Command-line tools for processing biological sequencing data; 2011, http://code.google.com/p/ea-utils.

27. Kelley DR, Schatz MC, Salzberg SL. Quake: quality-aware detection and correction of sequencing errors. *Genome Biol* 2010; **11**:R116. PubMed http://dx.doi.org/10.1186/gb-2010-11-11-r116

28. Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. *Genome Res* 2008; **18**:821-829. PubMed http://dx.doi.org/10.1101/gr.074492.107

29. Gordon D, Abajian C, Green P. Consed: a graphical tool for sequence finishing. *Genome Res* 1998; **8**:195-202. PubMed http://dx.doi.org/10.1101/gr.8.3.195

30. Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. *BMC Bioinformatics* 2010; **11**:119. PubMed http://dx.doi.org/10.1186/1471-2105-11-119

31. Finn DR, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. *Nucleic Acids Research* 2011, Web Server Issue **39**:W29-W37.

32. Lowe TM, Eddy SR. tRNAscan-SE: A Program for Improved Detection of Transfer RNA Genes in Genomic Sequence. *Nucleic Acids Res* 1997; **25**:955-964. PubMed http://dx.doi.org/10.1093/nar/25.5.0955

33. Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: Inference of RNA alignments. *Bioinformatics* 2009; **25**:1335-1337. PubMed http://dx.doi.org/10.1093/bioinformatics/btp157

34. Markowitz VM, Ivanova NN, Chen IMA, Chu K, Kyrpides NC. IMG ER: a system for microbial genome annotation expert review and curation. *Bioinformatics* 2009; **25**:2271-2278. PubMed http://dx.doi.org/10.1093/bioinformatics/btp393

35. Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, Hugenholtz P. CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. *BMC Bioinformatics* 2007; **8**:209. PubMed http://dx.doi.org/10.1186/1471-2105-8-209

36. Petersen J. Phylogeny and compatibility: plasmid classification in the genomics era. *Arch Microbiol* 2011; **193**:313-321. PubMed

37. Zielenkiewicz U, Ceglowski P. Mechanisms of plasmid stable maintenance with special focus on plasmid addiction systems. *Acta Biochim Pol* 2001; **48**:1003-1023. PubMed
38. Wommack KE, Colwell RR. Viroplankton: viruses in aquatic ecosystems. *Microbiol Mol Biol Rev* 2000; **64**:69-114. PubMed http://dx.doi.org/10.1128/MMBR.64.1.69-114.2000

39. Proctor LM, Fuhrman JA, Ledbetter MC. Marine bacteriophages and bacterial mortality. *Eos* 1988; **69**:1111-1112.

40. Paul JH. Prophages in marine bacteria: dangerous molecular bombs or the key to survival in the seas? *ISME J* 2008; **2**:579-589. PubMed http://dx.doi.org/10.1038/ismej.2008.35

41. Canchaya C, Proux C, Fournous G, Bruttin A. Prophage genomics. *Microbiol Mol Biol Rev* 2003; **67**:238-276. PubMed http://dx.doi.org/10.1128/MMBR.67.2.238-276.2003

42. Ackermann HW, DuBow MS. Viruses of Prokaryotes. Vol I. General Properties of Bacteriophages. CRC Press 1987; Inc.: Boca Raton, FL.

43. Miller MB, Bassler BL. Quorum sensing in bacteria. *Annu Rev Microbiol* 2001; **55**:165-199. PubMed http://dx.doi.org/10.1146/annurev.micro.55.1.165

44. Bassler BL. How bacteria talk to each other: regulation of gene expression by quorum sensing. *Curr Opin Microbiol* 1999; **2**:582-587. PubMed http://dx.doi.org/10.1016/S1369-5274(99)00025-9

45. Waters CM, Bassler BL. Quorum Sensing: Cell-to-Cell Communication in Bacteria. *Annu Rev Cell Dev Biol* 2005; **21**:319-246. PubMed http://dx.doi.org/10.1146/annurev.cellbio.21.012074.131001

46. Patzelt D, Wang H, Buchholz I, Rohde M, Gröbe L, Pradella S, Neumann A, Schulz S, Heyber S, Münch K, et al. You are what you talk: quorum sensing induces individual morphologies and cell division modes in *Dinoroseobacter shibae*. *ISME J* 2013; **7**:2274-2286. PubMed http://dx.doi.org/10.1038/ismej.2013.107

47. Stapley EO, Starkey RL. Decomposition of cysteic acid and taurine by soil microorganisms. *J Gen Microbiol* 1970; **64**:77-84. http://dx.doi.org/10.1099/00221287-64-1-77

48. Spring S, Jäckel U, Wagner M, Kämpfer P. *Ottowia thiooxydans* gen. nov., sp. nov., a novel facultatively anaerobic, N₂O-producing bacterium isolated from activated sludge, and transfer of *Aquaspirillum gracile* to *Hylemonella gracilis* gen. nov., comb. nov. *Int J Syst Evol Microbiol* 2004; **54**:99-106. PubMed http://dx.doi.org/10.1099/ijs.0.02727-0

49. Meier-Kolthoff JP, Auch AF, Klenk HP, Göker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. *BMC Bioinformatics* 2013; **14**:60. PubMed http://dx.doi.org/10.1186/1471-2105-14-60

50. Auch AF, Klenk HP, Göker M. Standard operating procedure for calculating genome-to-genome distances based on high-scoring segment pairs. *Stand Genomic Sci* 2010; **2**:142-148. PubMed http://dx.doi.org/10.4056/sigs.541628

51. Riedel T, Fiebig A, Han J, Huntemann M, Spring S, Petersen J, Ivanova NN, Markowitz V, Göker M, et al. Genome sequence of the *Wenxinia marina* type strain (DSM 24838T), a representative of the *Roseobacter* clade isolated from oilfield sediments. *Stand Genomic Sci* (under review).

52. Edgar RC. PILER-CR Fast and accurate identification of CRISPR repeats. *BMC Bioinformatics* 2007; **8**:18. PubMed http://dx.doi.org/10.1186/1471-2105-8-18

53. Meier-Kolthoff JP, Klenk HP, Göker M. Taxonomic use of the G+C content and DNA:DNA hybridization in the genomic age. *Int J Syst Evol Microbiol* 2014; **64**:352-356. PubMed http://dx.doi.org/10.1099/ijs.0.056994-0