High quality draft genome sequence of Segniliparus rugosus CDC 945T= (ATCC BAA-974T)

Segniliparus rugosus represents one of two species in the genus Segniliparus, the sole genus in the family Segniliparaceae. A unique and interesting feature of this family is the presence of extremely long carbon-chain length mycolic acids bound in the cell wall. S. rugosus is also a medically important species because it is an opportunistic pathogen associated with mammalian lung disease. This report represents the second species in the genus to have its genome sequenced. The 3,567,567 bp long genome with 3,516 protein-coding and 49 RNA genes is part of the NIH Roadmap for Medical Research, Human Microbiome Project.


Introduction
Strain CDC 945 T (= ATCC BAA-974 T = CIP 10838 T = DSM 45345 = CCUG 50838 T = JCM 13579 T ) is the type strain of the species Segniliparus rugosus in the Segniliparaceae family [1]. The genus name was created to acknowledge the presence of novel long carbon-chain fatty acids (mycolic acids) detected using the Mycobacterium species identification method with high performance liquid chromatography (HPLC) [2]. The name was formed from the Latin adjective 'segnis', meaning 'slow' and combined with the Greek adjective 'liparos' for 'fatty', to indicate the 'one with slow fats'. The name relates to the late elution of the apolar, alpha-mycolic acids (fatty acids) during HPLC analysis [1]. The specific epithet for the taxon name is from the Latin adjective 'rugosus', referring to the formation of wrinkled, rough colony morphology [1]. The type strain of S. rugosus, CDC 945 T , was isolated from a human sputum specimen collected in Alabama, USA [1]. S. rugosus has been isolated from multiple patients with cystic fibrosis in the U.S. and Australia and appears to be a respiratory opportunistic pathogen [3,4]. A recent isolation from a ~1 year old sea lion showing third-stage malnutrition with a 30% loss of body weight, moderate bradycardia and severe hypothermia, suggests a possible aquatic or marine niche for the species [5]. The only other validly named species of the genus is Segniliparus rotundus (CDC 1076 T ), which is the type strain of this species. S. rotundus. CDC 1076 shares 98.9% 16S rRNA sequence identity with S. rugosus CDC 945 T , although the DNA-DNA hybridization is less than 28% [1]. The complete genome of S. rotundus was recently reported and has 3,157,527 bp with 3,081 protein-coding and 52 RNA genes [6]. Here we present a summary classification and a set of features for S. rugosus CDC 945 T , together with the description of the high quality draft genomic sequencing and annotation.

Classification and features
The cells of CDC 945 T are irregular rods ranging in length and width from 0.55-0.90 µm by 1.9-4.5 µm (Table 1 and Figure 1). Colonies are wrinkled, rough and form in less than 7 days on 7H10 and 7H11 agar at an optimal temperature of 33 o C [1]. CDC 945 T is aerobic, non-motile, asporogenous, and stains bright red with acid alcohol stain [1]. It is mesophilic and demonstrates a temperature range for growth between 22 and 42 o C [1]. Colonies grown for < 4 weeks are non-pigmented, nonphotochromogenic and do not produce a diagnostic odor [1]. They do not produce aerial mycelium, spores or demonstrate true branching. Young colonies are creamy and smear easily when disturbed. Cell growth at ~4 weeks on Löwenstein-Jensen (LJ) medium produces a diffusible pink color in the agar at the leading edges of mature growth. Aged colonies on LJ develop a light buff pigment and demonstrate 'greening' from uptake of malachite green [1]. CDC 945 T is weakly positive for arylsulfatase at 7 days but is strongly positive in 14 days. No growth on MacConkey agar not containing crystal violet. CDC 945 T grows in the presence of 5% sodium chloride at 7 days, in lysozyme at 21 days. Positive for iron uptake, nitrate reduction, tellurite reduction and tween opacity. Negative for tween hydrolysis [1]. Results with the API CORYNE test kit shows CDC 945 T is positive for β-glucosidase, and pyrazinamidase activities and negative for alkaline phosphatase, β-galactosidase, β-glucuronidase, αglucosidase, N-acetyl-β-glucosaminidase and pyrrolidonyl arylamidase activity at 33 o C [1]. It is susceptible to imipenem 4ug/ml, moxifloxacin 0.5 μg/ml, and trimethoprim-sulfamethoxazole < 4.8 ug/ml, intermediate to cefoxitan 64 ug/ml and resistant to amikacin >128 ug/ml, clarithromycin 32 ug/ml, ciprofloxacin 16 ug/ml, ethambutol >16 ug/ml and tobramycin >64 ug/ml. [1,3]. Strain CDC 945 T uses D-glucose, glycerol, maltose, mannitol, D-sorbitol and trehalose as sole carbon sources with the production of acid. No growth on adonitol, L-arabinose, cellobiose, citrate, dulcitol, i-erythriol, galactose, i-myo-inositol, lactose, mannose, melibiose, raffinose, L-rhamnose, salicin or sodium citrate [1]. The strain hydrolyzes urea but not acetamide adenine, casein, aesculin, hypoxanthine, tyrosine or xanthine [1].
The phylogenetic association of Segniliparus species is shown in Figure 2 in a 16S rRNA based tree. The genus forms a distinct lineage relative to the other mycolic acid containing Actinobacteria. The positioning of Segniliparus in this phylogenetic analysis is consistent with its positioning in the "All-Species Living Tree Project" LTP release 106, August 2011, which is similarly based on 16S rRNA. [11].

Genome sequencing and annotation Genome project history
This organism was selected for sequencing on the basis of its association with respiratory lung disease and is part of the NIH Roadmap for Medical Research, Human Microbiome Project (HMP) [12]. The HMP presents reference genomes in Genome Online Database (GOLD) [10], the Human Microbiome Project Data Analysis and Coordination Center Project Catalog [13] and the complete high quality draft genome sequence is deposited in GenBank [14]. The Broad Institute performed the sequencing and annotation of this high quality draft genome [15]. A summary of the project is given in Table 2.

Growth conditions and DNA isolation
Strain CDC 945 T was grown statically in Middlebrook 7H9 medium at 33 o C until late log. DNA was isolated from whole cells after a chloroform/methanol wash with a disruption solution of guanidine thiocyamate, sarkosyl and mercaptoethanol as described in Mve-Obiang et al. [16]. The purity of DNA was assessed by The Broad Institute using the Quant-iT™ dsDNA Assay High Sensitivity Kit (Invitrogen, Carlsbad, CA) and according to the manufacturer's protocol.   [7] based on 16S rRNA sequences highlighting the position of S. rugosus CDC 945 T relative to the other type strains of mycolic acid containing genera in the suborder Corynebacterineae. GenBank accession numbers are listed after the name. The tree was inferred from 1,468 bp positions aligned using Clustal W [8] in MEGA v4 [9]. Numbers at the branch nodes are support values from 1,000 bootstrap replicates if equal to or greater than 70%. The scale bar indicates substitutions per site. The tree was rooted with Streptomyces coelicolor. Lineages with type strain genome sequencing projects registered in GOLD [10] are shown in blue, published genomes in bold. Standards in Genomic Sciences Table1. Classification and general features of S. rugosus CDC 945 T according to the MIGS recommendations [25].

MIGS ID Property Term Evidence code
Domain Bacteria TAS [26] Phylum Actinobacteria TAS [27] Class Actinobacteria TAS [28] Subclass Actinobacteridae TAS [28] Current classification Order Actinomycetales TAS [28,29] Suborder Corynebacterineae TAS [28,29] Family Segniliparaceae TAS [1,29] Genus 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 the Gene Ontology project [31].

Genome sequencing and assembly
The genome of Segniliparus rugosus ATCC BAA-974 was sequenced using 454 pyrosequence fragment and jump libraries [17]. We assembled the 454 data, consisting of 135,510 fragment reads and 112,271 jump reads, using Newbler Assembler version 2.3 PostRelease-11/19/2009. The assembly is considered High-Quality Draft and consists of 262 contigs arranged in 30 scaffolds with a total size of 3,567,567 bases. The error rate of this draft genome sequence is less than 1 in 10,000 (accuracy of ~ Q40). Average sequence coverage is 13×. Assessment of coverage, GC content, contig BLAST and 16S contig classification were consistent with the species Segniliparus.

Genome annotation
Protein-coding genes were predicted using four ORF-finding tools: GeneMark [18], Glimmer3 [19], Metagene [20], and findBlastOrfs (unpublished). This latter tool builds genes by extending wholegenome blast alignments, in-frame, to include start and stop codons. The final set of non-overlapping ORFs was selected from the output of these tools using an in-house gene-caller, which uses dynamic programming to score candidate gene models based on strength of similarity to entries in UniRef90, then selects non-overlapping genes that, combined, have the highest overall score. In cases where predictions overlapped non-coding RNA features (see below), the genes were manually inspected and removed when necessary. Finally, the gene set was reviewed using both the NCBI discrepancy report and the internal Broad annotation metrics. Ribosomal RNAs (rRNAs) were identified with RNAmmer [21]. The tRNA features were identified using tRNAScan [22]. Other non-coding features were identified with RFAM [23]. The gene product names were assigned based on Hmmer equivalogs from TIGRfam and Pfam, and blast hits to KEGG and SwissProt protein sequence databases. This was done using the naming tool "Pidgin" [24].

Genome properties
This 3,567,567 bp draft genome has high G+C content (Table 3 and Figure 3) and is predicted to encode 3,571 genes, 98% of which are protein coding. Nearly 70% of predicted proteins have a functional prediction and COG functional categories have been assigned to 53% of predicted proteins (Table 4). Standards in Genomic Sciences