Non-contiguous finished genome sequence and contextual data of the filamentous soil bacterium Ktedonobacter racemifer type strain (SOSP1-21T)

Ktedonobacter racemifer corrig. Cavaletti et al. 2007 is the type species of the genus Ktedonobacter, which in turn is the type genus of the family Ktedonobacteraceae, the type family of the order Ktedonobacterales within the class Ktedonobacteria in the phylum ‘Chloroflexi’. Although K. racemifer shares some morphological features with the actinobacteria, it is of special interest because it was the first cultivated representative of a deep branching unclassified lineage of otherwise uncultivated environmental phylotypes tentatively located within the phylum ‘Chloroflexi’. The aerobic, filamentous, non-motile, spore-forming Gram-positive heterotroph was isolated from soil in Italy. The 13,661,586 bp long non-contiguous finished genome consists of ten contigs and is the first reported genome sequence from a member of the class Ktedonobacteria. With its 11,453 protein-coding and 87 RNA genes, it is the largest prokaryotic genome reported so far. It comprises a large number of over-represented COGs, particularly genes associated with transposons, causing the genetic redundancy within the genome being considerably larger than expected by chance. This work is a part of the Genomic Encyclopedia of Bacteria and Archaea project.


Introduction
Strain SOSP1-21 T (= DSM 44963 = NRRL B-41538) is the type strain of the species Ktedonobacter racemifer, which is the type species of the monotypic genus Ktedonobacter, the type genus of the family Ktedonobacteraceae [1]. K. racemifer was first described in 2006 [1,2] as an aerobic, non-motile, filamentous, mesophilic, Gram-positive hetero-troph also capable of growing under microaerophilic conditions [1]. The genus name was derived from the Greek word ktedon -onos, fiber, and the Neo-Latin bacter, a rod, meaning a filamentous rod [1]. The species epithet is derived from the Latin adjective racemifer, carrying clusters of grapes [1].
The original spelling, Ktedobacter racemifer was corrected in 2007 on validation according to Rule 61 and Recommendation 6(7) [2]. Strain SOSP1-21 T was originally isolated from a soil sample of a black locust wood in Gerenzano, Northern Italy. Ten phylogenetically (class level) related strains were also isolated from soil samples collected at different locations in Northern Italy [1]. Only recently, a nearest cultivated neighbor, Thermosporothrix hazakensis, was isolated from hot compost in Japan [3]. Here we present a summary classification and a set of features for K. racemifer strain SOSP1-21 T , together with the description of the complete genomic sequencing and annotation.

Classification and features
Using NCBI BLAST [4], a representative genomic 16S rRNA sequence of K. racemifer SOSP1-21 T was compared under default settings (e.g., considering only the high-scoring segment pairs (HSPs) from the best 250 hits) with the most recent release of the Greengenes database [5] and the relative frequencies of taxa and keywords (reduced to their stem [6]) were determined, weighted by BLAST scores. The most frequently occurring genus was 'Ktedobacter' (100.0%) (1 hit in total; this represents the original, incorrect spelling of Ktedonobacter). No hits to sequences with (other) species names were found. (Note that the Greengenes database uses the INSDC (= EMBL/NCBI/DDBJ) annotation, which is not an authoritative source for nomenclature or classification.) The highestscoring environmental sequence was AM180157 ('New lineage filamentous spore-forming soil isolate SOSP1-30SOSP1-30 str. SOSP1-30'), which showed an identity of 99.0% and an HSP coverage of 95.2%. The most frequently occurring keywords within the labels of environmental samples which yielded hits were 'soil' (11.2%), 'prari, tallgrass' (4.9%), 'miner, weather' (1.9%), 'new' (1.8%) and 'filament, lineag, spore-form' (1.6%) (249 hits in total). These keywords reflect some of the ecological properties reported for strain SOSP1-21 T in the original description [1]. Environmental samples which yielded hits of a higher score than the highest scoring species were not found. Figure 1 shows the phylogenetic neighborhood of K. racemifer in a 16S rRNA based tree. The sequences of the eight 16S rRNA genes copies in the genome differ by up to nine nucleotides from each other and by up to five nucleotides from the previously published 16S rRNA sequence (AM180156), which contains two ambiguous base calls.
K. racemifer strain SOSP1-21 T cells are rodshaped, filamentous and grow both vegetative and aerial mycelia on solid medium (Figure 2a). The large aerial hyphae produce spherical spores that cluster together with a grape-like appearance ( Figure 2b). All other K. racemifer strains produced rounded spores, although they were arranged differently on the aerial hyphae [1]. Filamentous growth of strain SOSP1-21 T also occurred in submerged cultures, which contained the branched mycelia known from actinomycetes [1]. SOSP1-21 T stains Gram-positive and is not acid fast [1]. It produces pigments ranging from cream to pinkish orange on all media [1]. Although essentially aerobic, SOSP1-21 T is capable of growing under microaerophilic conditions [1]. The optimal growth temperature is 28-33°C [1]. It grows well at pH values between 4.8 and 6.8 with an optimum at pH 6 [1]. Salinity up to 10 g per liter does not inhibit the growth of the strain [1]. Strain SOSP1-21 T was capable of hydrolyzing starch, casein, gelatin, and (to a lesser extent) keratin but not cellulose, xylan, or chitin [1]. Strain SOSP1-21 T was catalase positive and produced H2S but could not reduce nitrates [1]. It is sensitive to 5 ug/ml novobiocin or ramoplanin and to 20 mg/ml apramycin and the glycopeptide A40926.

Figure 1.
Phylogenetic tree highlighting the position of K. racemifer relative to the other type strains within the phylum 'Chloroflexi'. The tree was inferred from 1,359 aligned characters [7,8] of the 16S rRNA gene sequence under the maximum likelihood (ML) criterion [9]. Rooting was done initially using the midpoint method [10] and then checked for its agreement with the current classification ( Table 1). The branches are scaled in terms of the expected number of substitutions per site. Numbers above branches are support values from 750 ML bootstrap replicates [11] (left) and from 1,000 maximum parsimony bootstrap replicates [12] (right) if larger than 60%. Lineages with type strain genome sequencing projects registered in GOLD [13] are labeled with one asterisk, those also listed as 'Complete and Published' with two asterisks [14][15][16][17] as well as CP001337, CP000804, CP000909, CP002084, and AP012029.   [18] and the NamesforLife database [19]. 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). These evidence codes are from of the Gene Ontology project [24].

Genome project history
This organism was selected for sequencing on the basis of its phylogenetic position [25], and is part of the Genomic Encyclopedia of Bacteria and Archaea project [26]. The genome project is deposited in the Genomes OnLine Database [13] 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. [27] adjusted to pH 6.0, at 28°C. DNA was isolated from 0.5-1 g of cell paste using Qiagen Genomic 500 DNA Kit (Qiagen 10262) following the manufacturer's protocol, with cell lysis protocol st/LALMP as described in Wu et al. [26]. DNA is available through the DNA Bank Network [28].

Genome sequencing and assembly
The genome was sequenced using a combination of Sanger and 454 sequencing platforms. All general aspects of library construction and sequencing can be found at the JGI website [29]. Pyrosequencing reads were assembled using the Newbler assembler (Roche). The initial Newbler contigs were broken into 14,080 overlapping fragments of 1,000 bp and entered as pseudo-reads into the subsequence assembly. The sequences were assigned quality scores based on Newbler consensus q-scores with modifications to account for overlap redundancy and to adjust inflated q-scores. A hybrid 454/Sanger assembly was produced using parallel phrap (High Performance Software, LLC). Possible mis-assemblies were corrected with Dupfinisher [30], or transposon bombing of bridging clones (Epicentre Biotechnologies, Madison, WI) [31]. Some gaps between contigs were closed by editing in Consed [32], custom primer walking or PCR amplification. A total of 3,354 Sanger finishing reads and five shatter libraries were produced to close gaps, to resolve some repetitive regions, and to raise the quality of the finished sequence. Illumina reads were also used to correct potential base errors and increase consensus quality using a software Polisher developed at JGI [33]. The error rate of the completed genome sequence is less than 1 in 100,000. Together, the combination of the Sanger and 454 sequencing platforms provided 34.7 × coverage of the genome. The final assembly contained 165,050 pyrosequence and 2,305,667 Illumina reads.

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) non-redundant database, Uni-Prot, TIGR-Fam, Pfam, PRIAM, KEGG, COG, and In-terPro databases. Additional gene prediction analysis and functional annotation were performed within the Integrated Microbial Genomes -Expert Review (IMG-ER) platform [36].

Genome properties
The non-contiguous finished genome consists of ten contigs ranging in size from 1,579 bp to almost four Mbp, with five contigs being longer than one Mb (1,302,518 bp, 2,713,222 bp, 2,766,182 bp, 2,916,502 bp, and 3,837,106 bp) and a G+C content of 53.8% (Table 3 and Figure 3). Of the 11,540 genes predicted, 11,453 were protein-coding genes, and 87 RNAs; No pseudogenes were identified. The majority of the protein-coding genes (61.2%) 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.    [41]. However, this genome was also one of the most difficult to assemble. Figure 4 shows the unusually high number of identical sequence fragments across the genome, which caused the termination of the project as noncontiguous finished genome without closure of the last ten sequence gaps.

Comparative genomics
Lacking an available genome sequence of the closest relative of K. racemifer, Thermosporothrix hazakensis [3] (Figure 1), the following comparative analyses were done with Sphaerobacter thermophilus [42] and Thermomicrobium roseum [43], the closest organisms phylogenetically for which there are publically available genome sequences [15,16]. K. racemifer stands out because of its enormous genome size of more than 13 Mbp. The genomes of S. thermophilus and T. roseum are significantly smaller, 3.9 Mbp and 2.9 Mbp, respectively. Whereas S. thermophilus and T. roseum have similar G+C-contents of 68% and 64%, respectively, the G+C-content of the K. racemifer genome is significantly lower (54%). The fraction of shared genes in the three genomes is shown in a Venn diagram ( Figure 5). The numbers of pairwise shared genes were calculated with the phylogenetic profiler function of the IMG-ER platform [36]. Homologous genes within the genomes were detected with a maximum E-value of 10 -5 and a minimum identity of 30%. A total of 1,393 genes are shared by the three genomes, referring to the whole genome sizes 39% and 48% of the genes in S. thermophilus and T. roseum have homologs in the three genomes, in the case of K. racemifer only 12% of the genes are shared by the other two genomes. The pairwise comparison of S. thermophilus and T. roseum revealed 2,249 genes which are shared by these two organisms, referring to the whole genomes 64% of the S. thermophilus genes and 79% of the T. roseum have homologous genes in the respective other genome.
The genome of K. racemifer encodes an enormously high number of transposon-associated genes; its annotation revealed 601 genes encoding transposases, 151 genes encoding integrases and 107 genes encoding resolvases. The genes coding these enzymes are spread over the whole genome with some regions having a higher density than others. The extremely high number of transposases is due to several gene copies that are to a greater or lesser extent similar in their sequences. The presence of that many mobile elements may explain the unusually high number of identical sequence fragments across the genome and the resulting difficulties occurring during the genome assembly.
Within the 9,539 unique genes of K. racemifer that have no detectable homologs in the genomes of S. thermophilus and T. roseum (under the sequence similarity thresholds used for the comparison) the 29 genes encoding xylose isomerases appear to be especially noteworthy; for 27 of these isomerase genes no homologous genes were detected in the other two genomes; only one gene was identified in T. roseum, and two in S. thermophilus. The high number of xylose isomerase genes suggests a strong utilization of pentoses by K. racemifer. To date K. racemifer was not tested regarding xylose utilization, but the close relative T. hazakensis is able to use xylose as the only carbon source [3]. Furthermore, a high number of genes encoding proteins responsible for resistance against several antibiotics were predicted: 61 bleomycin resistance proteins and 41 aminoglycoside phosphotransferases. An estimate of the overall similarity between K. racemifer, S. thermophilus and T. roseum, was generated with the GGDC Genome-to-Genome Distance Calculator [44,45]. This system calculates the distances by comparing the genomes to obtain HSPs (high-scoring segment pairs) and interfering distances from a set of formulas (1, HSP length / total length; 2, identities / HSP length; 3, identities / total length). Table 5 shows the results of the pairwise comparison between the three genomes.   The pairwise comparison ( respectively, is, if genome sizes are taken into consideration. In order to quantify the differences in gene redundancy between the three genomes, as well as to determine over-represented genes, we used approaches based on entropy and the broken-stick distribution, respectively, applied to the set of genes from either genome assigned to COGs. Shannon's entropy (see, e.g., pp. 214, 243 in [46]) H can be used as a measure of disorder for discrete distributions; it is maximum (Hmax) if all categories (COGs in our case) are represented by exactly one item (gene) and then equal to the logarithm of the number of items (genes).
Thus, one can measure the evenness (nonredundancy) within such a distribution as H/Hmax and the corresponding redundancy as 1.0 -H/Hmax. The broken-stick distribution reflects the relative abundance of a given number of categories within a random population of items (see, e.g., p. 244 and 410 in [46]). Over-represented items (here: COGs) are those whose real relative frequencies (here: number of genes assigned to this COG relative to the total number of genes assigned to COGs) are larger than the broken-stick value of the corresponding rank within the list of frequencies sorted in decreasing order. Moreover, the entropy Hexp of the broken-stick distribution can be used as an estimate for the expected entropy, yielding 1.0 -H/Hexp as an alternative measure of redundancy (which becomes negative when the evenness is larger than expected by chance).
The 2,022 genes assigned to 1,300 distinct COGs in the genome of T. roseum corresponded to an entropy of 6.912, an expected entropy of 6.748 and, hence, a redundancy of 9.20% if measured using Hmax and of -2.42% using Hexp, whereas S. thermophilus (2,619 genes assigned to 1,383 COGs) yielded an entropy of 6.837 (expected: 6.810) and a redundancy of 13.14% with Hmax and -0.39% with Hexp. In contrast, the 6,654 genes assigned to 1,731 distinct COGs in the genome of K. racemifer yielded an entropy of only 6.455 (expected: 7.034) and a redundancy of 26.67% (using Hmax) and 8.24% (using Hexp). That is, in contrast to the other two genomes the genes within the genome of K. racemifer are distributed less even than expected by chance.  represented COGs in K. racemifer, S. thermophilus and T. roseum revealed differences between the three organisms. Not surprisingly the genes coding transposases (COG0675; by far the most frequent one), integrases (COG3316) and resolvases (COG2452) can be found among the overrepresented COGs in K. racemifer ( Figure 6).
Our analyses also showed that genes belonging to the category COG3344 are over-represented in the genome of K. racemifer. COG3344 represents retron type reverse transcriptases, which are found in group II introns. Group II introns are large catalytic RNA molecules that act as mobile genetic elements [47]. They were first identified in mitochondria and chloroplast genomes, but with the increasing number of bacterial genome sequencing projects, the number of group II intron sequences in the databases also increased. Dai and Zimmerly reported in 2003 that a quarter of the sequenced bacterial genomes contain group II introns [48,49]. By using the IMG-ER platform [36] we calculated that approximately one third of the 2,727 sequenced bacterial genomes contain group II introns. In the genome of K. racemifer, 34 genes coding reverse transcriptases could be identified, all of them having the same domain structure with the reverse transcriptase domain followed by a maturase-specific domain and the C-terminal HNH-endonuclease domain.