Skip to main content
Download PDF
  • Extended genome report
  • Open access
  • Published:

High quality draft genome sequences of Pseudomonas fulva DSM 17717T, Pseudomonas parafulva DSM 17004T and Pseudomonas cremoricolorata DSM 17059T type strains

Standards in Genomic Sciences volume 11, Article number: 55 (2016) Cite this article

Abstract

Pseudomonas has the highest number of species out of any genus of Gram-negative bacteria and is phylogenetically divided into several groups. The Pseudomonas putida phylogenetic branch includes at least 13 species of environmental and industrial interest, plant-associated bacteria, insect pathogens, and even some members that have been found in clinical specimens. In the context of the Genomic Encyclopedia of Bacteria and Archaea project, we present the permanent, high-quality draft genomes of the type strains of 3 taxonomically and ecologically closely related species in the Pseudomonas putida phylogenetic branch: Pseudomonas fulva DSM 17717T, Pseudomonas parafulva DSM 17004T and Pseudomonas cremoricolorata DSM 17059T. All three genomes are comparable in size (4.6–4.9 Mb), with 4,119–4,459 protein-coding genes. Average nucleotide identity based on BLAST comparisons and digital genome-to-genome distance calculations are in good agreement with experimental DNA-DNA hybridization results. The genome sequences presented here will be very helpful in elucidating the taxonomy, phylogeny and evolution of the Pseudomonas putida species complex.

Introduction

During a taxonomic study of Pseudomonas strains isolated from rice, petroleum fields and oil-brine in Japan, Iizuka and Komagata [1] proposed two new species in 1963, Pseudomonas fulva and Pseudomonas straminea (as cited in Uchino et al. [2]). These new species produced a water-insoluble yellow pigment, but not a water-soluble fluorescent pigment. Later, seven P. fulva strains obtained from culture collections were re-characterized and compared with the strains of related species by Uchino and collaborators [2]. Phylogenetic analysis based on 16S rRNA sequences, experimental DNA-DNA hybridization results and phenotypic characteristics led to the proposal of two new species: Pseudomonas parafulva (2 strains) and Pseudomonas cremoricolorata (1 strain). Three of the four remaining strains were maintained in the species P. fulva , and the other strain was identified as P. straminaea [2]. In a multilocus sequence analysis, the type strains of P. fulva , P. parafulva and P. cremoricolorata clustered in the Pseudomonas putida phylogenetic branch and are considered members of the P. putida group in the Pseudomonas fluorescens lineage [3]. The three species are taxonomically and ecologically closely related. Strains from these species have been isolated from rice paddy samples or from Japanese unhulled rice. P. fulva strains have also been studied for their endophytic growth in Scots pines and for their roles in plant growth promotion and protection against plant pathogenic fungi [4, 5]. The antagonistic effect against plant pathogenic bacteria has also been demonstrated in other strains of Pseudomonas putida [6]. Additionally, P. fulva strains have been isolated from water collected from human-made container habitats of mosquitoes [7]. P. fulva was one of the most abundant species found in a survey of pseudomonads in human homes [8], and very recently P. fulva was identified as a member of a polymicrobial ventriculitis in humans [9]. The difficulty in identifying species closely related to P. putida in the clinical laboratory is highlighted by Rebolledo and collaborators [9]. Biosynthesis of medium-chain-length poly(3-hydroxyalkanoates) by a volatile aromatic hydrocarbons-degrading P. fulva has been proposed as candidate for the biotechnological conversion of toxic petrochemical wastes to valuable biopolymers [10].

In the context of the Genomic Encyclopedia of Bacteria and Archaea (GEBA) project [11], the permanent, high quality draft genomes of the type strains of P. fulva , P. parafulva and P. cremoricolorata are presented. The genome sequences have been annotated, and the results are discussed in relation to the taxonomy of members of the P. putida phylogenetic group.

Organisms information

Classification and features

The type strains of the three species, P. fulva DSM 17717T (=JCM 11242T =NRIC 0180T ), P. parafulva DSM 17004T (=AJ 2129T =JCM 11244T =NRIC 0501T ) and P. cremoricolorata DSM 17059T (=JCM 11246T =NRIC 0181T ), were obtained from the DSMZ. All strains were isolated by Iizuka and Komagata [1, 2] from Japanese rice paddies and were initially proposed as members of the new species P. fulva due to the deep yellow color of their colonies. P. fulva was included in the Approved Lists of Bacterial Names [12]. Uchino and collaborators re-characterized several strains obtained as P. fulva from culture collections and proposed two new species: P. parafulva (2 strains) and P. cremoricolorata (1 strain) [2].

All three type strains shared the basic phenotypic traits of the genus Pseudomonas : Gram-negative rods, motility via polar flagella, with strictly respiratory type of metabolism, catalase and oxidase activity and phylogenetic placement in the genus Pseudomonas on the basis of 16S rRNA gene sequencing. None of the three species produced water-soluble fluorescent pigments but produced a characteristic water-insoluble yellow pigment. Colonies appear smooth, round, flat to convex and pale/creamy yellow on nutrient agar. The three species were differentiated from each other by several phenotypic tests: presence of the arginine dihydrolase pathway, growth at 37 °C and assimilation of D-ribose, D-mannose, adonitol, 2-keto-D-gluconate, butyrate, valerate, caprate, isovalerate, itaconate, citraconate, glycerate, levulinate, Tween 80, p-hydroxybenzoate, inosine, glycine, L-ornithine, L-citrulline and nicotinate. An extensive list of phenotypic characteristics can be found in the original publication by Uchino et al. [2]. The classification and general features of P. fulva , P. parafulva and P. cremoricolorata type strains are given in Tables 1, 2 and 3.

Table 1 Classification and general features of P. fulva DSM 17717T [36]
Full size table
Table 2 Classification and general features of P. parafulva DSM 17004T [36]
Full size table
Table 3 Classification and general features of P. cremoricolorata DSM 17059T [36]
Full size table

Chemotaxonomic data

As reported by Uchino and collaborators [2] the DNA GC-content of the three type strains, as determined by chemical analysis, was 60.0 mol % in P. fulva and P. parafulva and 62.1 mol % in P. cremoricolorata . The percentages of G + C bases based on the genome analysis were 61.71 % for P. fulva DSM 17717T , 62.42 % for P. parafulva DSM 17004T and 63.47 % for P. cremoricolorata DSM 17059T . The GC-contents determined by chemical analysis were slightly lower than those inferred from genome sequences. The predominant respiratory quinone was ubiquinone Q-9, but Q-8 and Q-10 were also detected in smaller amounts. The major cellular fatty acids were C16:0, C16:1 and C18:1, and the major 3-hydroxy fatty acids were C10:0 and C12:0 [2].

For protein analysis, cells were cultured in Luria-Bertani broth aerobically, with shaking at 30 °C, harvested in the exponential growth phase and prepared for Whole-cell MALDI-TOF MS analysis using an Autoflex III mass spectrometer (Bruker Daltonik, Germany) as recommended by the manufacturer. Protein profiles clearly distinguished the type strains in the P. putida phylogenetic group [3]. A list of major proteins that met a minimum intensity threshold of 700, a minimum signal to noise threshold of 15, and a mass to charge ratio (m/z) higher than 3,000 and lower than 10,000 is included in Additional file 1.

Extended feature descriptions

Phylogenetic trees were reconstructed using different methods, namely the maximum-likelihood, maximum-parsimony and neighbor-joining algorithms integrated in MEGA version 6 bioinformatics package [13], and also using the FastME 2.0 phylogeny inference program [14]. All phylogenetic trees tested showed similar topologies and the same strain groupings. The derived phylogeny of the species in the P. putida phylogenetic group based on 16S rDNA gene sequencing had low resolution, and the bootstrap values of branches were low (Fig. 1a). Therefore, a phylogenetic tree based on a multilocus sequence analysis with the partial sequences of three housekeeping genes (16S rDNA, gyrB, and rpoD) was constructed as recommended by Mulet et al. [3] (Fig. 1b). Most branches were supported in the MLSA phylogenetic tree by high bootstrap values, and all type strains were clearly separated in the P. putida phylogenetic group. The strain groupings ( P. putida / Pseudomonas monteilii / P. parafulva / P. fulva ; Pseudomonas soli / Pseudomonas mosselii / Pseudomonas entomophila / Pseudomonas plecoglossicida and Pseudomonas donghuensis / Pseudomonas vranovensis / Pseudomonas alkylphenolica ) were maintained in all trees. ‘ Pseudomonas hunanensis ’ NCCB 100446 (proposed as a new species [15], but not yet validated) and Pseudomonas taiwanensis DSM 21245T branches were supported by low bootstrap values, and their positions varied in the trees.

Fig. 1
figure 1

a Phylogenetic tree based on the almost complete sequence of the 16S rRNA gene of the type strains of the species in the P. putida phylogenetic group. Evolutionary distances were computed with MEGA (version 6) using the Jukes-Cantor method [13]. Dendrogram was generated by the Neighbor-Joining method. The bar indicates the number of base substitutions per site (1264 positions in final dataset). Percentage bootstrap values above 50 % (from 1000 replicates) are indicated at the nodes. b Phylogenetic tree based on concatenated 16S rRNA, gyrB and rpoD partial gene sequences of the type strains of species in the P. putida phylogenetic group. Evolutionary distances were computed with MEGA (version 6) using the Jukes-Cantor method [13]. Dendrogram was generated by the Neighbor-Joining method. The bar indicates the number of base substitutions per site (2758 positions in the final dataset). Percentage bootstrap values above 50 % (from 1000 replicates) are indicated at the nodes. The black dots indicate nodes maintained in all phylogenetic trees tested

Full size image

Genome sequencing information

Genome project history

Sequencing of the three type strains is part of Genomic Encyclopedia of Type Strains, Phase I: the one thousand microbial genomes (KMG-I) project [16], a follow-up of the GEBA pilot project [11, 17]. Project information is deposited in the Genomes on Line Database (GOLD) [18], and the high quality draft genome sequence is deposited in GenBank and in the Integrated Microbial Genomes database (IMG) [19]. Draft sequencing, initial gap closure and annotation were performed by the DOE Joint Genome Institute (JGI) using state-of-the-art sequencing technology [20]. A summary of the project information is shown in Table 4. Genbank IDs are as follows: JHYU00000000 for P. fulva DSM 17717T , AUEB00000000 for P. parafulva DSM 17004T and AUEA00000000 for P. cremoricolorata DSM 17059T .

Table 4 Project information for P. fulva DSM 17717T, P. parafulva DSM 17004T, and P. cremoricolorata DSM 17059T
Full size table

Growth conditions and genomic DNA preparation

P. fulva DSM 17717T , P. parafulva DSM 17004T and P. cremoricolorata DSM 17059T were cultured aerobically in Luria-Bertani medium, with shaking at 30°C, to the early stationary phase. Genomic DNA was extracted and purified with a Promega Wizard® Genomic DNA Purification kit, following the manufacturer’s instructions. DNA quality and quantity were determined with a Nanodrop spectrometer (Thermo Scientific, Wilmington, USA).

Genome sequencing and assembly

An Illumina standard shotgun library was constructed and sequenced using the Illumina HiSeq 2000 platform. For P. fulva , 16,075,374 reads were generated, totaling 2,411.3 Mb, of which 1,500 Mb were used in an assembly, resulting in an average coverage of 312.5x. For P. parafulva , 14,015,748 reads were generated, totaling 2,102.4 Mb, of which 607.0 Mb were used in an assembly, resulting in an average coverage of 122x. For P. cremoricolorata , 11,621,460 reads were generated, totaling 1,743.2 Mb, of which 569.4 Mb were used in an assembly, resulting in an average coverage of 122x. Illumina sequencing and library artifacts were removed using Duk filtering (L. Mingkun, A. Copeland, and H. J. Duk, unpublished data). Filtered Illumina reads were assembled using Velvet (version 1.1.04) [21], simulated paired-end reads were created from Velvet contigs using wgsim and simulated read pairs were reassembled using Allpaths-LG (version r42328) [22].

Genome annotation

Protein-coding genes were identified using Prodigal [23], as part of the DOE-JGI genome annotation pipeline [24]. Additional gene prediction analysis and manual functional annotation were performed within the Integrated Microbial Genomes (IMG) platform, which provides tools for analyzing and reviewing the structural and functional annotations of genomes in a comparative context [19]. Genome annotation procedures are detailed in Markowitz et al. [19] and references therein. Briefly, the predicted CDSs were translated and used to search the NCBI nonredundant database, UNIProt, TIGRFam, Pfam, KEGG, COG and InterPro databases. Transfer RNA genes were identified using the tRNAScan-SE tool and other non-coding RNAs were found using INFERNAL. Ribosomal RNA genes were predicted using hmmsearch against the custom models generated for each type of rRNA.

Genome properties

The assembly of the three genomes consisted of 4.7 Mb in 48 scaffolds for P. fulva , 4.9 Mb in 33 scaffolds for P. parafulva and 4.6 Mb in 27 scaffolds for P. cremoricolorata (Table 5). The G + C content for each strain was 61.72, 62.42 and 63.47 %, respectively. The majority of protein-coding genes (78.96, 80.59 and 79.68 %) were assigned a putative function. The properties and statistics of the genomes are summarized in Table 5, and the number of genes associated with general COG functional categories is shown in Table 6.

Table 5 Genome statistics for P. fulva DSM 17717T, P. parafulva DSM 17004T and P. cremoricolorata DSM 17059T
Full size table
Table 6 Number of genes associated with general COG functional categories for P. fulva DSM 17717T, P. parafulva DSM 17004T, and P. cremoricolorata DSM 17059T
Full size table

Insights from the genome sequence

Experimental DNA-DNA hybridizations were performed by Uchino et al. [2], following the fluorometric procedure proposed by Ezaki et al. [25]. Taxonomic genome comparisons were calculated by two different procedures: Average nucleotide identity based on BLAST was calculated with the JSpecies program [26]. Digital DDH similarities among the genomes of the three type strains were calculated using GGDC web server version 2.0 [27] under recommended settings. The results are given in Table 7 and are highly consistent. Experimental and dDDH values were clearly below the 70 species threshold, and ANIb was below the accepted 95–96 % species threshold. The type strains P. fulva NBRC 16637T and P. parafulva NBRC 16636T have been sequenced at the NBRC, and the ANIb values between the genomes of the type strain pair P. fulva DSM 17717 T and P. fulva NBRC 16637T were in 99.98 % agreement; the dDDH values were 100 % identical. The same results were obtained when comparing the type strain pair P. parafulva DSM 17004T and P. parafulva NBRC 16636T (99.95 and 100 % similarity for ANIb and GGDC results, respectively). Using the whole-genome ANI-based MiSI method [28], which is computed for all bacterial genomes in the Integrated Microbial Genomes system, P. fulva DSM 17717T clustered in the same gANI clique with eight plant-associated genome-sequenced Pseudomonas sp. not yet classified at the species level, with an intra-clique ANI of 99.57 %, indicating that the 9 strains belong genomically to the same species, P. fulva . The strain P. fulva NBRC 16637T is the equivalent type strain of the NITE (Biological Resource Center) and was included in the same clique. The GC-content variation within the clique was less than 1 % (61.58 %–61.88 %), which is proof of the value of draft genomes for taxonomy because the GC-content varies no more than 1 % within species [29], and all strains in the clique should be considered strains in the same genomic species [28]. Three additional genomes of strains identified as P. fulva , P. parafulva and P. cremoricolorata , available in the Genbank database on June 17, 2015, were also analyzed. The completely sequenced genome of P. cremoricolorata ND07 (CP009455) showed ANIb and dDDH values of 92 and 50 %, respectively, with P. cremoricolorata DSM 17059T , indicating a close relationship that is below the species threshold. The complete genome of P. parafulva CRS01-1 (CP009747) showed an ANIb value of 81 % with the type strain of P. parafulva , the closest related type strain. Finally, as was previously documented, the complete sequenced genome of strain P. fulva 12-X (CP002727) demonstrated that it is clearly a distinct species, with an ANIb value of 75.24 % [30] with the P. fulva type strain. In all three cases, the genome comparisons did not support a correct species affiliation of the strains.

Table 7 Experimental and digital genome similarities calculated for the P. fulva DSM 17717T, P. parafulva DSM 17007T and P. cremoricolorata DSM 17059T type strains
Full size table

The presence of genes related to carbohydrate and amino acid transport and metabolism is relevant for the fitness of environmental bacteria. These genes represent 12–13 % of the total genes detected in the three strains, and they also have taxonomic consequences. Substrate utilization is an essential criterion for Pseudomonas taxonomy, and several tests routinely used in Pseudomonas identifications have been employed in the present study.

Catalase and superoxide dismutase are relevant enzymes for protecting the cell against reactive oxygen and are characteristic of most Pseudomonas . Catalase activity was detected by Uchino et al. [2] in P. fulva , P. parafulva , and P. cremoricolorata . Accordingly, 3 genes potentially coding for catalase were found in P. fulva and P. parafulva , but only 2 were found in P. cremoricolorata ; three genes coding for superoxide dismutase were detected in P. fulva and P. cremoricolorta genomes, but only two were found in P. parafulva . Testing for the presence of the arginine dihydrolase (or arginine deiminase) pathway, in combination with other biochemical tests, can also be of diagnostic value [31]. The arginine dihydrolase pathway transforms arginine to ornithine with ATP gain and allows limited growth in several Pseudomonas under anaerobic conditions. The arcA gene is present in the P. fulva and P. parafulva genomes but is absent in P. cremoricolorata , in accordance with the experimental data obtained by Uchino et al. [2]. All three strains were negative for nitrate reduction, nitrate respiration and PHB synthesis, and, accordingly, no gene related to these pathways was detected in any of the genomes. Cleavage of aromatic compounds was also tested using protocatechuate as a substrate; a gene coding for the protocatechuate 3,4-dioxygenase (3-oxoadipate pathway) was found in P. fulva and P. parafulva but was absent in P. cremoricolorata , confirming the ortho cleavage of the aromatic ring as reported by Uchino and collaborators [2]. All three strains possessed genes encoding key enzymes involved in glucose catabolism via the glycolysis, pentose-phosphate and 2-keto-3-deoxi-6 phosphogluconate pathways. The three species were recorded as amylase negative, but an alpha-amylase gene (amyA) was detected in all three genomes, indicating the potential ability to grow with starch as a substrate.

Bacterial secretion systems transport proteins across the cell envelope of Gram-negative bacteria to the external milieu and are considered critical for persistence in an ecological niche and for conquering a new one [32]. Type VI secretion system seems to be the most common and appears to be confined to proteobacteria. The TVISS consists of 13 essential conserved genes, many of which contain a number of functionally accessory elements. Several TVISS are often present in a single genome [33]. They have been mainly studied for their pathogenic role in the interaction between bacteria and hosts, but TVISS seems to play a role in mutualistic relationships between bacteria and eukaryotic cells or between bacteria, as well. A set of 15 conserved TVISS genes were found in P. fulva DSM 17717T but were absent in the other two strains. P. fulva DSM 17717T also has three copies of a Rhs element Vgr protein not present in the other strains that can be exported by the TVISS, but its exact function is still not known. The possible role of TVISS genes in the pathogenesis or in the interactions of P. fulva DSM 17717T with the environment remains to be elucidated.

Prophage-like elements in microbial genomes represent one of the main contributors of mobile DNA, also known as the mobilome [34], and are the main reason for bacterial intraspecies variability. The prophage contribution to the bacterial genome is highly variable. It can represent up to 8 % of the total chromosomal DNA [35], but phages may also be absent. The mobilomes of P. fulva , P. parafulva and P. cremoricolorata were predicted to contain 32, 30, and 44 genes, respectively. In addition to transposases, integrases and regulatory elements, clusters of bacteriophage structural genes (6 to 13 genes in a cluster) were found in the 3 strains: 2 clusters in P. fulva DSM 17717T (6 and 12 genes in each cluster), 2 in P. parafulva DSM 17004T (12 and 9 genes) and 3 in P. cremoricolorata DSM 17059T (9, 13 and 10). CRISPR elements were not detected.

Conclusions

Genome comparisons confirmed the distinct species status of the three type strains analyzed, as well as the close relationships between them. The complete genome analysis also revealed important taxonomic results, highlighting the relevance of the correct species assignation of strains and the need for the genome sequences of species type strains to build a phylogenomic taxonomy. No discrepancies were found between the genome insights and the phenotypic traits previously published for the species. However, the gene content revealed potential properties not yet detected, such as the presence of secretion systems, whose relevance remains to be explored. The genome sequences of the three type strains will be very helpful in elucidating the phylogeny and evolution of the P. putida species complex, a relevant coherent group of closely-related species with important ecological and biotechnological implications.

Abbreviations

ANIb:

Average nucleotide identity based on BLAST

dDDH:

Digital DNA-DNA hybridization

DDH:

DNA-DNA hybridization

gANI:

Whole-genome ANI-based MiSI method

GEBA:

Genomic encyclopedia of Bacteria and Archaea

KMG:

One thousand microbial genomes

MIGS:

Minimum information about a genome sequence

MLSA:

Multilocus sequence analysis

NAS:

Non-traceable

TAS:

Traceable.

References

  1. Iizuka H, Komagata K. New species of Pseudomonas belonging to fluorescent group (Studies on the microorganisms of cereal grains. Part V). Nippon Nogeikagaku Kaishi. 1963;37:137–41 (in Japanese).

    Article  Google Scholar 

  2. Uchino M, Shida O, Uchimura T, Komagata K. Recharacterization of Pseudomonas fulva Iizuka and Komagata, and proposals of Pseudomonas parafulva sp.nov. and Pseudomonas cremocricolorata sp. nov. J Gen Appl Microbiol. 2001;47:247–61.

    Article  CAS  PubMed  Google Scholar 

  3. Mulet M, Gomila M, Scotta C, Sánchez D, Lalucat J, García-Valdés E. Concordance between whole-cell matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and multilocus sequence analysis approaches in species discrimination within the genus Pseudomonas. Syst Appl Microbiol. 2012;35:455–64.

    Article  CAS  PubMed  Google Scholar 

  4. Strzelczyk E, Li CY. Bacterial endobionts in the big non-mycorrhizal roots of Scots pine (Pinus sylvestris L.). Microbiol Res. 2000;155:229–32.

    Article  CAS  Google Scholar 

  5. Pokojska-Burdjiez A, Strzelczyk E, Dahm H, Li CY. Effect of endophytic bacterium Pseudomonas fulva on growth of pine seedlings (Pinus sylvestris), formation of mycorrhizae and protection against pathogens. Phytopathol Pol. 2004;32:33–47.

    Google Scholar 

  6. Liu Q, Zhang Y, Yu N, Bi Z, Zhu A, Zhan X, et al. Genome sequence of Pseudomonas parafulva CRS01-1, an antagonistic bacterium isolated from rice field. J Biotech. 2015;20:89–90.

    Article  Google Scholar 

  7. Ponnusamy L, Xu N, Stav G, Wesson DM, Schal D, Apperson CS. Diversity of bacterial communities in container habitats of mosquitoes. Microb Ecol. 2008;56:593–603.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Remold SK, Brown CK, Farris JE, Hundley TC, Perpich JA, Purdy ME. Differential habitat use and niche partitioning by Pseudomonas species in human homes. Microb Ecol. 2011;62:505–17.

    Article  PubMed  Google Scholar 

  9. Rebolledo PA, Vu CCL, Donahue Carlson R, Kraft CS, Anderson EJ, Burd EM. Polymicrobial ventriculitis in volving Pseudomonas fulva. J Clin Microbiol. 2014;52:2239–41.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Ni YY, Kim Do Y, Chung MG, Lee SH, Park HY, Rhee YH. Biosynthesis of medium-chain-length poly(3-hydroxyalkanoates) by volatile aromatic hydrocarbons-degrading Pseudomonas fulva TY16. Bioresour Technol. 2010;101:8485–8.

    Article  CAS  PubMed  Google Scholar 

  11. Kyrpides NC, Hugenholtz P, Eisen JA, Woyke T, Göker M, Parker CT, et al. Genomic encyclopedia of bacteria and archaea: sequencing a myriad of type strains. PLoS Biol. 2014;12:e1001920.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Skerman VBC, McGowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Bacteriol. 1980;30:225–420.

    Article  Google Scholar 

  13. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Lefort V, Desper R, Gascuel O. FastME 2.0: a comprehensive, accurate and fast distance-based phylogeny inference program. Mol Biol Evol. 2015;32:2798–800.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Gao J, Li BY, Wang HH, Liu ZQ. Pseudomonas hunanensis sp. nov., isolated from soil subjected to long-term manganese pollution. Curr Microbiol. 2014;69:19–24.

    Article  CAS  PubMed  Google Scholar 

  16. Kyrpides NC, Woyke T, Eisen JA, Garrity G, Lilburn TG, Beck BJ, et al. Genomic encyclopedia of type strains, phase I: the one thousand microbial genomes (KMG-I) project. Stand Genomic Sci. 2013;9:628–6234.

    Article  CAS  Google Scholar 

  17. Wu D, Hugenholtz P, Mavromatis K, Pukall R, Dalin E, Ivanova NN, et al. A phylogeny-driven genomic encyclopaedia of Bacteria and Archaea. Nature. 2009;462:1056–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Reddy TBK, Thomas A, Stamatis D, Bertsch J, Isbandi M, Jansson J, et al. The Genomes OnLine Database (GOLD) v.5: a metadata management system based on a four level (meta)genome project classification. Nucl Acids Res. 2015;43:D1099.

    Article  CAS  PubMed  Google Scholar 

  19. Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Pillay M, et al. IMG 4 version of the integrated microbial genomes comparative analysis system. Nucleic Acids Res. 2014;42(Database issue):D560–7.

    Article  CAS  PubMed  Google Scholar 

  20. Mavromatis K, Land M, Brettin TS, Quest DJ, Copeland A, Clum A, et al. The fast changing landscape of sequencing technologies and their impact on microbial assemblies and annotation. PLoS One. 2012;7:e48837.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de bruijn graphs. Genome Res. 2008;18:821–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. MacCallum I, Przybylski D, Gnerre S, Burton J, Shlyakhter I, Gnirke A, et al. ALLPATHS 2: small genomes assembled accurately and with high continuity from short paired reads. Genome Biol. 2009;10:R103.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Huntemann M, Ivanova NN, Mavromatis K, Tripp HJ, Paez-Espino D, Palaniappan K, et al. The standard operating procedure of the DOE-JGI Microbial Genome Annotation Pipeline (MGAP v.4). Stand Genomic Sci. 2015;10:86.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Ezaki T, Hashimoto Y, Yabuuchi E. Fluorometric deoxyribonucleic acid-deoxyribonucleic acid hybridization in microdilution wells as an alternative to membrane filter hybridization in which radioisotopes are used to determine genetic relatedness among bacterial strains. Int J Syst Bacteriol. 1989;39:224–9.

    Article  Google Scholar 

  26. Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009;106:19126–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Varghese NJ, Mukherjee S, Ivanova N, Konstantinidis KT, Mavrommatis K, Kyrpides NC, Pati A. Microbial species delineation using whole genome sequences. Nucleic Acids Res. 2015;43(14):6761–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Meier-Kolthoff JP, Klenk HP. Göker M Taxonomic use of DNA G + C content and DNA-DNA hybridization in the genomic age. Int J System Evol Microbiol. 2014;64:352–6.

    Article  CAS  Google Scholar 

  30. Gomila M, Peña A, Mulet M, Lalucat J, Garcia-Valdes E. Phylogenomics and systematic in Pseudomonas. Front Microbiol. 2015;6:214.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Moore ERB, Tindall BJ, Martins Dos Santos VAP, Pieper DH, Ramos J-L, Palleroni NJ. Nonmedical: Pseudomonas. In: Dworkin M, Falkow S, Rosenberg E, Schleifer K-H, Stackebrandt E, editors. The Prokaryotes: a handbook on the biology of bacteria, vol. 6. 3rd ed. New York, NY: Springer; 2006. p. 646–703.

    Google Scholar 

  32. Decoin V, Barbey C, Bergeau D, Latour X, Feuilloley MG, Orange N, Merieau A. A type VI secretion system is involved in Pseudomonas fluorescens bacterial competition. PLoS One. 2014;9(2):e89411.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Schwarz S, Hood RD, Mougous JD. What is type VI secretion doing in all those bugs? Trends Microbiol. 2010;18:531–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Canchaya C, Fournous G, Chibani-Chennoufi S, Dillmann ML, Brussow H. Phage as agents of lateral gene transfer. Curr Opin Microbiol. 2003;6:417–24.

    Article  CAS  PubMed  Google Scholar 

  35. Ventura M, Canchaya C, Tauch A, Chandra G, Fitzgerald GF, Chater KF, van Sinderen D. Genomics of actinobacteria: tracing the evolutionary history of an ancient phylum. Microbiol Mol Biol Rev. 2007;71:495–548.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26:541–7.

    Article  CAS  PubMed Central  Google Scholar 

  37. Woese CR, Kandler O, Weelis ML. Towards a natural system of organisms. Proposal for the domains Archaea and Bacteria. Proc Natl Acad Sci U S A. 1990;87:4576–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Garrity GM, Bell JA, Lilburn T, Phylum X. Proteobacteria phyl nov. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey’s Manual of Sytematic Bacteriology, vol. 2. 2nd ed. New York: Springer; 2005.

    Google Scholar 

  39. Garrity GM, Bell JA, Lilburn T. Class III. Gammaproteobacteria class nov. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey’s Manual of Sytematic Bacteriology, vol. 2. 2nd ed. New York: Springer; 2005.

    Google Scholar 

  40. Garrity GM, Bell JA, Lilburn T. Order IX. Pseudomonadales. In: Brenner DJ, Krieg NR, Stanley JT, Garrity GM, editors. Bergey’s Manual of Sytematic Bacteriology, The Proteobacteria), part B (The Gammaproteobacteria, vol. 2. Secondth ed. New York: Springer; 2005. p. 323.

    Chapter  Google Scholar 

  41. Garrity GM, Bell JA, Lilburn T. Family I. Pseudomonadaceae. In: Brenner DJ, Krieg NR, Stanley JT, Garrity GM, editors. Bergey’s Manual of Sytematic Bacteriology, The Proteobacteria), part B (The Gammaproteobacteria, vol. 2. Secondth ed. New York: Springer; 2005. p. 323.

    Chapter  Google Scholar 

  42. Palleroni NJ. Genus I. Pseudomonas Mingula 1894. In: Krieg NR NR, Holt JG, editors. Bergey’s Manual of Systematic Bacteriology, vol. I. Baltimore: Williams and Wilkins; 1984. p. 141.

    Google Scholar 

Download references

Acknowledgements

Financial support was obtained from the Spanish MINECO through projects CGL2011-24318, CGL2015-70925-P and Consolider CSD2009-00006, from the EuroNanoMedII project NanoDiaBac (AC14/00017), as well as funds for competitive research groups from the Government of the Balearic Islands (with FEDER cofunding). Margarita Gomila is the recipient of a postdoctoral contract from the Conselleria d’Educació, Cultura i Universitats del Govern de les Illes Balears and the European Social Fund. Arantxa Peña was supported by a postdoctoral contract from the University of the Balearic Islands.

Authors’ contributions

AB and MM cultured the strains, isolated, purified and quality controlled the DNA and performed the multilocus sequence analysis of strains described in 2014. RG performed the WC-MALDI-TOF MS analysis. TBKR, MH, AP, NI, VM and NCK performed draft sequencing, initial gap closure and annotation. AP and MG analyzed the annotation and performed the statistical analysis. EGV, JL, NK and HPK conceived the study and participated in its design and coordination. All authors helped to draft the manuscript. All authors read and approved the final manuscript. None of the authors have any competing interests regarding the manuscript.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

  1. Department of Biology-Microbiology, Universitat de les Illes Balears, Campus UIB, Crtra. Valldemossa km 7.5, 07122, Palma de Mallorca, Spain

    Arantxa Peña, Antonio Busquets, Margarita Gomila, Magdalena Mulet, Elena García-Valdés & Jorge Lalucat

  2. Serveis Cientifico-Tècnics, Universitat de les Illes Balears, Palma de Mallorca, Spain

    Rosa M. Gomila

  3. DOE Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, CA, 94598-1698, USA

    T. B. K. Reddy, Marcel Huntemann, Amrita Pati, Natalia Ivanova, Victor Markowitz, Tanja Woyke & Nikos Kyrpides

  4. Institut Mediterrani d’Estudis Avançats (IMEDEA, CSIC-UIB), Palma de Mallorca, Spain

    Elena García-Valdés & Jorge Lalucat

  5. Leibniz Institute DSMZ - German Collection of Microorganisms and Cell Cultures, 38124, Braunschweig, Germany

    Markus Göker

  6. School of Biology, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK

    Hans-Peter Klenk

  7. Department of Biological Sciences, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia

    Nikos Kyrpides

Authors
  1. Arantxa Peña
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Antonio Busquets
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Margarita Gomila
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Magdalena Mulet
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Rosa M. Gomila
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. T. B. K. Reddy
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Marcel Huntemann
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Amrita Pati
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Natalia Ivanova
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Victor Markowitz
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Elena García-Valdés
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Markus Göker
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Tanja Woyke
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Hans-Peter Klenk
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Nikos Kyrpides
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Jorge Lalucat
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Jorge Lalucat.

Additional file

Additional file 1:

Major protein profiles of P. fulva DSM 17717T, P. parafulva DSM 17004T and P. cremoricolorata DSM 17059T type strains. The intensity value was determined as an average from all spectra containing that peak and the relative intensities with respect to the base peak in percentage are indicated within parenthesis. (m/z: mass to charge ratio; + and -: presence or absence of the corresponding protein). (PDF 99 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Peña, A., Busquets, A., Gomila, M. et al. High quality draft genome sequences of Pseudomonas fulva DSM 17717T, Pseudomonas parafulva DSM 17004T and Pseudomonas cremoricolorata DSM 17059T type strains. Stand in Genomic Sci 11, 55 (2016). https://doi.org/10.1186/s40793-016-0178-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40793-016-0178-2

Keywords

Environmental Microbiome

ISSN: 2524-6372