Complete genome sequence of Eggerthella lenta type strain (IPP VPI 0255T)

Eggerthella lenta (Eggerth 1935) Wade et al. 1999, emended Würdemann et al. 2009 is the type species of the genus Eggerthella, which belongs to the actinobacterial family Coriobacteriaceae. E. lenta is a Gram-positive, non-motile, non-sporulating pathogenic bacterium that can cause severe bacteremia. The strain described in this study has been isolated from a rectal tumor in 1935. Here we describe the features of this organism, together with the complete genome sequence, and annotation. This is the first complete genome sequence of the genus Eggerthella, and the 3,632,260 bp long single replicon genome with its 3123 protein-coding and 58 RNA genes is part of the Genomic Encyclopedia of Bacteria and Archaea project.


Introduction
Strain VPI 0255 T (= DSM 2243 = ATCC 25559 = JCM 9979) is the type strain of the species Eggerthella lenta, which was first described in 1935 by Eggerth as 'Bacteroides lentus' [1], later in 1938 renamed by Prévot in 'Eubacterium lentum' [2], and was also known under the synonym 'Pseudobacterium lentum' Krasil'nikov 1949 [3]. The strain has been described in detail by Moore et al. in 1971 [4]. Based on 16S rRNA sequence divergence and the presence of unique phenotypic characters the strain was then transferred to the new genus Eggerthella as E. lenta (Kageyama et al. 1999, Wade et al. 1999 [5,6] In 2004 two novel Eggerthella species, E. hongkongensis and E. sinensis were characterized and described in addition [7].
Recently, E. hongkongensis was reclassified as Paraeggerthella hongkongensis [8]. Although the two Eggerthella species and P. hongkongensis are part of the human gut flora, they can be the agent of severe bacteremia. So far the pathogenic potential of the genera are poorly analyzed [7]. Here we present a summary classification and a set of features for E. lenta VPI 0255 T , together with the description of the complete genomic sequencing and annotation.

Classification and features
Members of the species E. lenta have been isolated from several abscesses, from appendix tissues, peritoneal fluid and intestinal tumors. The organ-ism is often involved in mixed infections with less fastidious bacteria. Difficulties in cultivation and identification are probably the reason why bacteremia caused by Eggerthella is rarely reported. Half of the cases of Eggerthella bacteremia are induced by the two novel species: E. sinensis and P. hongkongensis [7]. Stinear et al. described an isolate (AF304434) from human feces resembling E. lenta (98% identity) that carries an enterococcal vanB resistance locus probably received via lateral gene transfer or as a result of genetic mutations [9]. Clavel et al. investigated the occurrence and activity of dietary lignans activating bacterial communities in human feces and identified an E. lenta strain (AY937380) with 98.2% sequence similarity to strain VPI 0255 T [10]. Lignans are a class of phytoestrogen which can be metabolized to the biologically active enterolignans, enterodiol and enterolactone. The human intestinal microbiota is essential for the conversation of the dietary lignans e.g. secoisolariciresinol diglucoside via secoisolariciresinol (SECO) to the enterolignans. Clavel and co-workers also reported that the dehydroxylation of SECO is catalyzed by Eggerthella lenta [11]. Based on 16S rRNA gene sequence analyses another five uncultured clones with 99% identity to E. lenta were reported at the NCBI BLAST server (status June 2009). These clones were derived from the analyses of feces samples from humans. e.g. associated with obesity [12,13], but also from marine metagenomes [14] Figure 1 shows the phylogenetic neighborhood of E. lenta strain VPI 0255 T in a 16S rRNA based tree. The sequences of the three identical copies of the 16S rRNA gene in the genome differ by three nucleotides from the previously published 16S rRNA sequence generated from ATCC 25559 (AF292375). The slight difference between the genome data and the reported 16S rRNA gene sequence is most likely due to sequencing errors in the previously reported sequence data.

Figure 1.
Phylogenetic tree of E. lenta strain VPI 0255 T and all type strains of the genus Eggerthella as well as the type strains from all other genera of the family Coriobacteriaceae inferred from 1,373 aligned characters [15,16] of the 16S rRNA gene under the maximum likelihood criterion [17]. The branches are scaled in terms of the expected number of substitutions per site. Numbers above branches are support values from 1,000 bootstrap replicates if larger than 60%. Lineages with type strain genome sequencing projects registered in GOLD [18] are shown in blue, published genomes in bold, including two of which are reported in this issue of SIGS [19,20] E. lenta strain VPI 0255 T was originally isolated from a rectal tumor and described as Grampositive, non-motile and non-sporulating ( Table 1) [1]. Cells are rod shaped and occur singly or in long chains up to 20 elements ( Figure 2). The cell size and morphology varies depending on the substrate and the age of the culture. Surface colonies were described as circular to slightly scalloped, convex, shiny, gray and translucent. E. lenta is ob-ligately anaerobic and its optimal growth temperature is 37° C [4]. Growth is stimulated by arginine. The existence of the arginine dihydrolase pathway as an important energy source was described by Sperry and Wilkens in 1976 [26]. E. lenta is asaccharolytic [4,26,29], Gelatin is not liquefied, aesculin is not hydrolyzed and nitrate is reduced [29]. E. lenta is bile-resistant and primarily found in human feces [6]. Standards in Genomic Sciences

Chemotaxonomy
The cell wall of E. lenta strain VPI 0255 T contains A1γ-type peptidoglycan glutamic acid occurred in D-form and diaminopimelic acid in meso configuration. Mycolic acids and teichonic acids were not reported. Strain VPI 0255 T contains menaquinone MK-6 as the major respiratory lipoquinone (63.7%) and a lower amount of the methylmenaquinone MMK-6 (36.3%) [8,29,31]. As the predominant fatty acids the unbranched saturated 16:0 DMA (29.4%) and the monounsaturated fatty acid 18:1w9c (22.0%) were identified [5,6]. Polar lipids consist of two phospholipids, phosphatidylglycerol and diphosphatidylglycerol, and four glycolipids GL1-GL4 [8]. Altitude not reported 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 the Gene Ontology project [30]. If the evidence code is IDA, then the property was directly observed for a living isolate by one of the authors, or an expert or reputable institution mentioned in the acknowledgements.

Genome project history
This organism was selected for sequencing on the basis of each phylogenetic position, and is part of the Genomic Encyclopedia of Bacteria and Archaea project. The genome project is deposited in the Genome OnLine Database [18] and the complete genome sequence 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
E. lenta strain VPI 0255 T , DSM 2243, was grown anaerobically in DSMZ medium 209 (Eubacterium lentum Medium [32]) at 37°C. DNA was isolated from 1-1.5 g of cell paste using Qiagen Genomic 500 DNA Kit (Qiagen, Hilden, Germany) following the manufacturer's protocol without modifications.

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. 454 Pyrosequencing reads were assembled using the Newbler assembler version 1.1.02.15 (Roche). Large Newbler contigs were broken into 4,901 overlapping fragments of 1,000 bp and entered into the assembly as pseudo-reads. The sequences were assigned quality scores based on Newbler consensus q-scores with modifications to account for overlap redundancy and to adjust inflated qscores. A hybrid 454/Sanger assembly was made using the parallel phrap assembler (High Performance Software, LLC). Possible mis-assemblies were corrected with Dupfinisher or transposon bombing of bridging clones [33]. Gaps between contigs were closed by editing in Consed, custom primer walk or PCR amplification. A total of 358 Sanger finishing reads were produced to close gaps, to resolve repetitive regions, and to raise the quality of the finished sequence. The final assembly consists of 39,464 Sanger and 471,609 pyrosequence (454) reads. Together all sequence types provided 35.5x coverage of the genome. The error rate of the completed genome sequence is less than 1 in 100,000.

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, UniProt, TIGRFam, Pfam, PRIAM, KEGG, COG, and InterPro 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 is 3,632,260 bp long and comprises one main circular chromosome with a 64.2% GC content ( Table 3 and Figure 3). Of the 3,181 genes predicted, 3,123 were protein coding genes, and 58 RNAs. 53 pseudogenes were also identified. A majority of the genes (70.9%) were assigned with a putative function while the remaining genes were annotated as hypothetical proteins. The properties and the statistics of the genome are summarized in Table 3. The distribution of genes into COGs functional categories is presented in Table 4.