- Open Access
Draft genome sequence of Dyadobacter tibetensis type strain (Y620-1) isolated from glacial ice
Standards in Genomic Sciences volume 9, pages883–892 (2014)
Dyadobacter tibetensis Y620-1 is the type strain of the species Dyadobacter tibetensis, isolated from ice at a depth of 59 m from a high altitude glacier in China (5670 m above sea level). It is psychrotolerant with growth temperature ranges of 4 to 35°C. Here we describe the features of this organism, together with the draft genome sequence and annotation. The 5,313,963 bp long genome contains 4,828 protein-coding genes and 39 RNA genes. To the best of our knowledge, this is the first Dyadobacter strain that was isolated from glacial ice. This study provides genetic information of this organism to identify the genes linked to its specific mechanisms for adaption to extreme glacial environment.
Strain Y620-1 (=JCM 18589= CGMCC 1.12215T) is the type strain of the species Dyadobacter tibetensis . The genus Dyadobacter currently has 12 species after it was first proposed by Chelius and Triplett on 2000, and the type species is D. fermentans . Those species isolated from diverse environment, i.e. glacial ice, soil from the Arctic, Colorado Plateau, farm and a ginseng field, desert sand, freshwater and sea water, and plant material [1–12]. So far, however, the genome sequences have been determined for only three Dyadobacter strains (D. alkalitolerans DSM 23607 (GCA_000428845), D. fermentans DSM 18053 (GCA_000023125), D. beijingensis DSM 21582 (GCA_000382205)), and only the complete genome sequence of D. fermentans DSM 18053 has been published .
D. tibetensis strain Y620-1 was isolated from 59 m depth section of an 122 m ice core drilled from Yuzhufeng Glacier at 5670 m above sea level, Tibetan Plateau, China . Glacier ice is an extreme environment with low temperature and nutrients, but high UV radiation, and is a huge reservoir of extremophilic microorganism that have accumulated for hundreds of years . Diverse isolates were recovered from glacial ice, but the genomes of bacteria in the extreme environment were limited [15,16]. Here, we present the genome sequence of psychrotolerant D. tibetensis strain Y620-1 isolated from ice core. This is the first genome sequence of a bacterial isolated from a deep high altitude glacier ice.
Classification and features
The phylogenetic position of genus Dyadobacter is in the Cytophagaceae, a very diverse family within the order Sphingobacteriales, the phylum Bacteroidetes . Closest related genera are Persicitalea and Runella . D. tibetensis strain Y620-1 represents a novel species of the genus Dyadobacter based on reported genotypic and phenotypic feature .
Strain Y620-1 was isolated at 4°C incubation from ice core melt water cultivated on R2A medium . It was Gram-stain-negative, non-motile and rod-shaped with 1 µm to 2 µm length  (Figure 1). Colonies are yellow, round, smooth, convex and opaque on R2A after incubation at 30°C for 2–3 days. Growth occurs at 4–35°C on R2A, with an optimum at 30°C. It produces a flexirubin-like pigment, the same as the other species in the genus. The major fatty acids are summed feature 3 (C16: 1ω7c and/or iso-C15: 0,2-OH), iso-C15: 0, C16: 1ω5c and iso-C17: 0 3-OH. The predominant polar lipid is phosphatidylethanolamine .The important characteristics of the strain based on literature descriptions are summarized in Table 1. The strain exhibited a 16S rRNA gene sequence similarity with other members of the genus Dyadobacter ranging from 95.1% with D. ginsengisoli Gsoil 043 to 93.6% with D. beijingensis A54 (Figure 2).
The utilization of carbon compounds by strain Y620-1 was determined using Generation-III microplates on an OmniLog phenotyping device (BIOLOG Inc., Hayward, CA, USA). The microplates were inoculated at 30°C with a cell suspension at a cell density of 95–96% turbidity and dye IF-A. Strain Y620-1 assimilates dextrin, D-maltose, D-trehalose, D-cellobiose, gentiobiose, D-melibiose, D-salicin, N-Acetyl-D-glucosamine, D-mannose, Glycyl-L-proline, L-alanine, L-histidine, L-serine, methyl, pyruvate, L-lactic acid, citric acid, α-Keto-glutaric acid, L-malic acid, propionic acid and acetic acid, but not stachyose, D-raffinose, N-Acetyl-β-dmannosamine, N-Acetyl-dgalactosamine, N-Acetyl-neuraminic acid, D-galactose, 3-methyl glucose, L-fucose, L-rhamnose, inosine, D-sorbitol, D-mannitol, D-arabitol, myo-Inositol, glycerol, D-glucose-6-PO4, D-aspartic acid, gelatin, L-arginine, L-pyroglutamic acid, pectin, L-galactonic acid lactone, D-gluconic acid, glucuronamide, mucic acid, quinic acid, D-saccharic acid, β-hydroxy-phenylacetic acid, D-lactic acid methyl ester, D-malic acid, γ-amino-butryric acid, β-hydroxy-D,L-butyric Acid and formic acid. It is sensitive to minocycline, potassium tellurite, nalidixic acid, lithium chloride, fusidic acid, D-serine and sodium bromate, but not troleandomycin, rifamycin, lincomycin, guanidine HCl, niaproof 4, vancomycin, tetrazolium violet, tetrazolium blue, aztreonam and sodium butyrate.
Genome sequencing information
Genome project history
The organism was selected for sequencing on the basis of it from extreme deep ice core from high altitude glacier. The shotgun genome sequencing project was completed in December 2012 and has been deposited at DDBJ/EMBL/GenBank under the accession number AZQN00000000. The version described here is the first version, AZQN01000000. The genome sequencing was carried out in Shanghai Majorbio Bio-pharm Technology Co., Ltd (Shanghai, China). A summary of the project information is shown in Table 2
Growth conditions and DNA isolation
Cells of strain Y620-1 were harvested from R2A broth following 2 days incubation at 30 °C with shaking at 180 rpm. The genomic DNA of the strain was extracted according to the method previously described by Marmur et al. . Extraction was started with 100 ml of 48 h culture, centrifuged at 4 °C and 10, 000 rpm for 15 min. Then, cells were washed three times with 5 ml sterile water. The washed cells were resuspended in 1,128 µL Tris-HCl buffer (10 mM) containing 1 mM EDTA (pH 8.0) and 20 µg lysozyme and incubated at 37 °C for 2 h. followed by adding of 6 µL proteinase K (20 mg/mL), 4 µL DNase-free RNase (10 mg/mL), 100 µL SDS (20% w/v) and the cell suspension was incubated at 55 °C for 3 h. The cell lysate was extracted twice with phenol/chloroform/isoamyl alcohol (25:24:1) and once with chloroform/isoamyl alcohol (24:1), and the aqueous layer was separated after centrifugation at 12,000 rpm for 15 min. The DNA was precipitated with 1 volumes of frozen anhydrous ethanol. The purity of genomic DNA was assessed by NanoDrop (2000c, Thermo) with OD 260:280 ratio of 1.8–2. The DNA was stock in TE (pH 8.0) for genome sequencing.
Genome sequencing and assembly
The genome of strain Y620-1 was sequenced using an Illumina GAIIx instrument with two paired-end libraries (170 bp and 800 bp insert size). The raw sequencing data was processed to discard reads containing adaptor sequences, a high rate of ambiguity, and removing the sequence reads which were of low quality. A total of 2,041 Mb high-quality of Illumina data were obtained, providing approximately 384–fold coverage. The high-quality reads were assembled in silico using SOAPdenovo v1.05, resulting in 33 contigs (> 200 bp) with an N50 length of 797,100 bp.
The coding sequences (CDS) were predicted using Glimmer 3.02 , while tRNAscan-SE  and RNAmmer  were used to identify tRNA and rRNA, respectively. The genome sequence was also uploaded into the Rapid Annotation using Subsystem Technology (RAST) system  to check the annotated sequences. The functions of predicted protein-coding genes were then annotated through comparisons with the databases of NCBI-NR , COG , and KEGG . The program TMHMM  and SignalP  were used to identify putative transmembrane helices and signal peptides.
The Y620-1 draft genome sequence has a total of 5,313,963 bp with an average GC content of 43.44%. There are 4,867 predicted genes, of which 4,828 are protein-coding genes, and 39 are RNA genes. A total of 2,844 genes (58.91%) are assigned a putative function. The remaining genes were annotated as either hypothetical proteins or proteins of unknown functions. Using COG functional assignment, 70.55% of protein coding genes could be classified into 20 COG categories. The properties and the statistics of the genome are summarized in Tables 3 and 4. According to the subsystem-based annotation generated by RAST, ∼ 32% protein-coding genes of strain Y620-1 could be assigned to 358 metabolic subsystems. The most abundant of the subsystems are related to carbohydrates (n=323, 7.2% of total protein-coding genes), following by amino acids and derivatives (n=272, 6.0%), cofactors, vitamins, prosthetic groups, pigments (n=184, 4.1%), protein metabolism (n=143, 3.2%), membrane transport (132, 2.9%) and respiration (131, 2.9%).
Although there were 12 species assigned to Dyadobacter, only three have a completed genome sequence and are D. alkalitolerans DSM 23607, D. fermentans DSM 18053, and D. beijingensis DSM 21582. Strain Y620-1 has the smallest genome of D. alkalitolerans DSM 23607, D. fermentans DSM 18053 and D. beijingensis DSM 21582 (6.29 Mbp, 6.97 Mbp and 7.37 Mbp, respectively). The GC content of strain Y620-1 is comparable to that of D. alkalitolerans DSM 23607 (45.66%), but lower than those of D. fermentans DSM 18053 (51.54%) and D. beijingensis DSM 21582 (52.09%). In order to estimate the similarity among the sequenced Dyadobacter strains, an average nucleotide identity (ANI) and Genome-to-Genome Distance Calculator (GGDC) were calculated using the software JSpecies v1.2  and GGDC v2.0 , respectively. Table 5 shows the results of ANI and GGDC.
ANI analysis showed that strain Y620-1 shared a low degree of similarity with other Dyadobacter species (< 69% ANIb and < 84% ANIm), whereas relatively higher ANI value were obtained for D. alkalitolerans DSM 23607, D. fermentans DSM 18053 and D. beijingensis DSM 21582 (Table 5). Although the core concept of GGDC was based on ‘genome blast distance phylogeny’, which is different from ANI , GGDC analysis showed similar results. In both analyses, the highest similarity values were observed in the comparisons of D. fermentans DSM 18053 with D. beijingensis DSM 21582. These results were in line with phylogeny analysis based on 16S rRNA gene, which shows that D. fermentans DSM 18053 and D. beijingensis DSM 21582 form a cluster with Dyadobacter soli MJ20. Moreover, the comparison of distribution of COG categories in the genome of four Dyadobacter strains revealed that there were significant correlations between the distribution of COG categories of strain Y620-1 and other strains (r = 0.970–0.979). However, relatively higher correlation coefficients were observed for D. alkalitolerans DSM 23607 and D. fermentans DSM 18053 (0.988), D. alkalitolerans DSM 23607 and D. beijingensis DSM 21582 (0.991), and D. fermentans DSM 18053 and D. beijingensis DSM 21582 (0.991) (Figure 3).
Five cold-shock proteins were found in this genome including CspA, GyrA, RbfA, and NusA. The proteins coded by gene RecA, RecF, RecG, RecN, RecO, RecQ, RadA and RadC, which play a critical role in recombinational repair of damaged DNA, were also found . Single-stranded-DNA-specific exonuclease RecJ, required for many types of recombination events , and CRISPRs Cas1, interacts with components of the DNA repair systems, also were found . Phage shock protein C existed in the genome, which may play a significant role in the competition for survival under nutrient- or energy-limited conditions . When bacteria deposit on the glacier, the low temperature, high UV radiation, and desiccation could induce the cold-shock and recombinational repair of damaged DNA proteins. Additionally, oligotrophic condition of glacial ice may induce the phage shock protein C. The genome sequence of strain Y620-1 provides genetic information to identify the genes linked to its specific mechanisms for adaption to extreme glacial environment.
Shen L, Liu Y, Yao T, Wang N, Xu B, Jiao N, Liu H, Zhou Y, Liu X, Wang Y. Dyadobacter tibetensis sp. nov., isolated from glacial ice core. Int J Syst Evol Microbiol 2013; 63:3636–3639. PubMed http://dx.doi.org/10.1099/ijs.0.050328-0
Chelius MK, Triplett EW. Dyadobacter fermentans gen. nov., sp. nov., a novel gram-negative bacterium isolated from surface-sterilized Zea mays stems. Int J Syst Evol Microbiol 2000; 50:751–758. PubMed http://dx.doi.org/10.1099/00207713-50-2-751
Baik KS, Kim MS, Kim EM, Kim HR, Seong CN. Dyadobacter koreensis sp. nov., isolated from fresh water. Int J Syst Evol Microbiol 2007; 57:1227–1231. PubMed http://dx.doi.org/10.1099/ijs.0.64902-0
Chaturvedi P, Reddy GSN, Shivaji S. Dyadobacter hamtensis sp. nov., from Hamta glacier, located in the Himalayas, India. Int J Syst Evol Microbiol 2005; 55:2113–2117. PubMed http://dx.doi.org/10.1099/ijs.0.63806-0
Chen L, Jiang F, Xiao M, Dai J, Kan W, Fang C, Peng F. Dyadobacter arcticus sp. nov., isolated from high Arctic soil on the Svalbard Archipelago, Norway. Int J Syst Evol Microbiol 2013; 63:1616–1620. PubMed http://dx.doi.org/10.1099/ijs.0.044198-0
Chun J, Kang JY, Joung Y, Kim H, Joh K, Jahng KY. Dyadobacter jejuensis sp. nov., isolated from seawater. Int J Syst Evol Microbiol 2013; 63:1788–1792. PubMed http://dx.doi.org/10.1099/ijs.0.043620-0
Dong Z, Guo X, Zhang X, Qiu F, Sun L, Gong H, Zhang F. Dyadobacter beijingensis sp. nov., isolated from the rhizosphere of turf grasses in China. Int J Syst Evol Microbiol 2007; 57:862–865. PubMed http://dx.doi.org/10.1099/ijs.0.64754-0
Lee M, Woo SG, Park J, Yoo SA. Dyadobacter soli sp. nov., a starch-degrading bacterium isolated from farm soil. Int J Syst Evol Microbiol 2010; 60:2577–2582. PubMed http://dx.doi.org/10.1099/ijs.0.019653-0
Liu QM, Im WT, Lee M, Yang DC, Lee ST. Dyadobacter ginsengisoli sp. nov., isolated from soil of a ginseng field. Int J Syst Evol Microbiol 2006; 56:1939–1944. PubMed http://dx.doi.org/10.1099/ijs.0.64322-0
Reddy GSN, Garcia-Pichel F. Dyadobacter crusticola sp. nov., from biological soil crusts in the Colorado Plateau, USA, and an emended description of the genus Dyadobacter Chelius and Triplett 2000. Int J Syst Evol Microbiol 2005; 55:1295–1299. PubMed http://dx.doi.org/10.1099/ijs.0.63498-0
Tang Y, Dai J, Zhang L, Mo Z, Wang Y, Li Y, Ji S, Fang C, Zheng C. Dyadobacter alkalitolerans sp. nov., isolated from desert sand. Int J Syst Evol Microbiol 2009; 59:60–64. PubMed http://dx.doi.org/10.1099/ijs.0.001404-0
Zhang DC, Liu HC, Xin YH, Zhou YG, Schinner F, Margesin R. Dyadobacter psychrophilus sp. nov., a psychrophilic bacterium isolated from soil. Int J Syst Evol Microbiol 2010; 60:1640–1643. PubMed http://dx.doi.org/10.1099/ijs.0.017236-0
Lang E, Lapidus A, Chertkov O, Brettin T, Detter JC, Han C, Copeland A, Glavina Del Rio T, Nolan M, Chen F and others. Complete genome sequence of Dyadobacter fermentans type strain (NS114 T). Stand Genomic Sci 2009; 29; 1(2):, 133–140.
Priscu JC, Christner BC, Foreman CM, Royston-Bishop G. Biological material in ice cores. In: Elias SA, editor. Encyclopedia of Quaternary Sciences. Volume 2: Elsevier B.V., UK.; 2007. p 1156–1166.
Koh HY, Lee SG, Lee JH, Doyle S, Christner BC, Kim HJ. Draft genome sequence of Paenisporosarcina sp. strain TG-14, a psychrophilic bacterium isolated from sediment-laden stratified basal ice from Taylor Glacier, McMurdo Dry Valleys, Antarctica. J Bacteriol 2012; 194:6656–6657. PubMed http://dx.doi.org/10.1128/JB.01795-12
Raymond JA, Christner BC, Schuster SC. A bacterial ice-binding protein from the Vostok ice core. Extremophiles 2008; 12:713–717. PubMed http://dx.doi.org/10.1007/s00792-008-0178-2
Reasoner DJ, Geldreich EE. A new medium for the enumeration and subculture of bacteria from potable water. Appl Environ Microbiol 1985; 49:1–7. PubMed
Garrity GM, Lilburn TG, Cole JR, Harrison SH, Euzéby J, Tindall BJ. Taxonomic Outline of the Bacteria and Archaea. Release 7.7 March 6, 2007. Michigan State University Board of Trustees 2007; DOI: 10.1601/TOBA7.7 http://www.taxonomicoutline.org.
Skerman VBD, McGrowan V, Sneath PHA. Approved list of bacterial names. Int J Syst Bacteriol: ASM Press; 1980. 225–420 p.
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. Nat Genet 2000; 25:25–29. PubMed http://dx.doi.org/10.1038/75556
Bruno WJ, Socci ND, Halpern AL. Weighted neighbor joining: a likelihood-based approach to distance-based phylogeny reconstruction. Mol Biol Evol 2000; 17:189–197. PubMed http://dx.doi.org/10.1093/oxfordjournals.molbev.a026231
Cole JR, Chai B, Farris RJ, Wang Q, Kulam-Syed-Mohideen AS, McGarrell DM, Bandela AM, Cardenas E, Garrity GM, Tiedje JM. The ribosomal database project (RDP-II): introducing myRDP space and quality controlled public data. Nucleic Acids Res 2007; 35:D169–D172. PubMed http://dx.doi.org/10.1093/nar/gkl889
Marmur J. A procedure for the isolation of deoxyribonucleic acid from microorganisms. J Mol Biol 1961; 3:208–218. http://dx.doi.org/10.1016/S0022-2836(61)80047-8
Delcher AL, Bratke KA, Powers EC, Salzberg SL. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 2007; 23:673–679. PubMed http://dx.doi.org/10.1093/bioinformatics/btm009
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 1997; 25:0955–0964.
Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res 2007; 35:3100–3108. PubMed http://dx.doi.org/10.1093/nar/gkm160
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics 2008; 9:75. PubMed http://dx.doi.org/10.1186/1471-2164-9-75
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL. GenBank. Nucleic Acids Res 2007; 36:D25–D30. PubMed http://dx.doi.org/10.1093/nar/gkm929
Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science 1997; 278:631–637. PubMed http://dx.doi.org/10.1126/science.278.5338.631
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res 2007; 36:D480–D484. PubMed http://dx.doi.org/10.1093/nar/gkm882
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol 2001; 305:567–580. PubMed http://dx.doi.org/10.1006/jmbi.2000.4315
Bendtsen JD, Nielsen H, von Heijne G, Brunak S. Improved prediction of signal peptides: SignalP 3.0. J Mol Biol 2004; 340:783–795. PubMed http://dx.doi.org/10.1016/j.jmb.2004.05.028
Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci USA 2009; 106:19126–19131. PubMed http://dx.doi.org/10.1073/pnas.0906412106
Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J 2008; 50:346–363. PubMed http://dx.doi.org/10.1002/bimj.200810425
Auch AF, von Jan M, Klenk H-P, Göker M. Digital DNA-DNA hybridization for microbial species delineation by means of genome-to-genome sequence comparison. 2010.
Kunst F, Ogasawara N, Moszer I, Albertini AM, Alloni G, Azevedo V, Bertero MG, Bessieres P, Bolotin A, Borchert S, et al. The complete genome sequence of the Gram-positive bacterium Bacillus subtilis. Nature 1997; 390:249–256. PubMed http://dx.doi.org/10.1038/36786
Lovett ST, Kolodner RD. Identification and purification of a single-stranded-DNA-specific exonuclease encoded by the recJ gene of Escherichia coli. Proc Natl Acad Sci USA 1989; 86:2627–2631. PubMed http://dx.doi.org/10.1073/pnas.86.8.2627
Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, et al. The complete genome sequence of Escherichia coli K-12. Science 1997; 277:1453–1462. PubMed http://dx.doi.org/10.1126/science.277.5331.1453
Brissette JL, Weiner L, Ripmaster TL, Model P. Characterization and sequence of the Escherichia coli stress-induced psp operon. J Mol Biol 1991; 220:35–48. PubMed http://dx.doi.org/10.1016/0022-2836(91)90379-K
This study was financially supported by the National Natural Science Foundation of China (Grant Nos. 41371084, 41171050, 41125003, 40930526, 41025002, 40871045), Strategic Priority Research Program (B) of the Chinese Academy of Sciences (XDB03030100), and the National Basic Research Program of China (No: 2010CB951404).