|Home | About | Journals | Submit | Contact Us | Français|
Strain Ola 51T (=LMG 24251T=CGMCC 1.7012T) is the type strain of the species Kosakonia oryzae and was isolated from surface-sterilized roots of the wild rice species Oryza latifolia grown in Guangdong, China. Here we summarize the features of the strain Ola 51T and describe its complete genome sequence. The genome contains one circular chromosome of 5,303,342 nucleotides with 54.01% GC content, 4773 protein-coding genes, 16 rRNA genes, 76 tRNA genes, 13 ncRNA genes, 48 pseudo genes, and 1 CRISPR array.
The online version of this article (doi:10.1186/s40793-017-0240-8) contains supplementary material, which is available to authorized users.
Enterobacter cowanii , E. radicincitans , E. oryzae , E. arachidis , E. sacchari , E. oryziphilus [6, 7], and E. oryzendophyticus [6, 7] have been transferred into the novel genus Kosakonia of the family “Enterobacteriaceae” [8–10]. A novel species “Kosakonia pseudosacchari”  closely related to K. sacchari was recently proposed. With the exception of the type species K. cowanii , which was originally obtained from clinical samples , the other members of the genus Kosakonia are nitrogen-fixing bacteria associated with plants [2–6, 11] and commonly occur in the nitrogen-fixing bacterial community of some non-legume crops, such as rice  and sugarcane . Some nitrogen-fixing Kosakonia strains are able to promote crop growth [12–14].
Strain Ola 51T (=LMG 24251 T=CGMCC 1.7012 T) is the type strain of the species Kosakonia oryzae and was isolated from surface-sterilized roots of the wild rice species Oryza latifolia grown in Guangdong, China . Here we present the summary of the features of the K. oryzae type strain Ola 51T and its complete genome sequence, which provides a reference for resolving the phylogeny and taxonomy of closely related strains and the genetic information to study its plant growth-promoting potential and its plant-associated life style.
K. oryzae strain Ola 51T is a Gram-negative, non-spore-forming, motile rod with peritrichous flagella (Fig. 1). It grows aerobically but reduces N2 to NH3 at a low pO2. It forms circular, convex, smooth colonies with entire margins on nutrient agar [3, 8]. It grows best around 30 °C and pH 7 (Table 1) . K. oryzae Ola 51T has the typical biochemical phenotypes of the genus Kosakonia: positive for acetoin production (Voges-Proskauer test) while negative for indole production; positive for β-galactosidase and arginine dihydrolase while negative for lysine decarboxylase; positive for oxidation of arabinose, cellobiose, citrate, fructose, galactose, gluconate, glucose, glycerol, lactose, malate, maltose, mannitol, mannose, sorbitol, sucrose and trehalose (Table 1) [3, 8].
The 16S rRNA gene sequence of K. oryzae Ola 51T was deposited in GenBank under the accession number EF488759 . A phylogenetic analysis of the 16S rRNA gene sequences from the strains belonging to the genus Kosakonia and Escherichia coli ATCC11775 T (the type strain of the type species of the type genus of the family Enterobacteriaceae) showed that K. oryzae Ola 51T is most closely related to the strains belonging to the species K. radicincitans (Fig. 2) [3, 8–11].
Whole-cell fatty acids were extracted from cells grown aerobically at 28 °C for 24 h on the TSA medium according to the recommendations of the Microbial Identification System (MIDI Inc., Delaware USA). The whole-cell fatty acid composition was determined using a 6890 N gas chromatograph (Agilent Technologies, Santa Clara, USA) and the peaks of the profiles were identified using the TSBA50 identification library version 5.0 (MIDI). K. oryzae Ola 51T shows the typical cell fatty acid profile of the genus Kosakonia . The major fatty acids are C16:0, C18:1 ω7c, C16:1 ω7c/15:0 iso 2OH, C17:0 cyclo and C14:0 3OH/16:1 iso I [8, 11].
K. oryzae Ola 51T was selected for sequencing based on its taxonomic significance. The genome sequence is deposited in GenBank under the accession number CP014007. A summary of the genome sequencing project information and its association with MIGS version 2.0  is shown in Table 2.
K. oryzae Ola 51T was grown aerobically in liquid Luria-Bertani medium at 30 °C until early stationary phase. The genome DNA was extracted from the cells by using a TIANamp bacterial DNA kit (Tiangen Biotech, Beijing, China). DNA quality (OD260/OD280=1.8) and quantity (22 μg) were determined with a Nanodrop spectrometer (Thermo Scientific, Wilmington, USA).
The genomic DNA of K. oryzae Ola 51T was constructed into 8 – 11 kb insert libraries and sequenced using PacBio SMRT sequencing technology  at the Duke University Genome Sequencing & Analysis Core Resource. Sequencing was run on two SMRT cells and resulted in 124,997 high-quality filtered reads with an average length of 8,260 bp. High-quality reads were assembled by the RS_HGAP_Assembly.3 in the SMRT analysis v2.3.0. The final assembly produced 128-fold coverage of the genome.
Automated genome annotation was done using the NCBI Prokaryotic Genome Annotation Pipeline . Functional annotations were done by searching against the KEGG , InterPro , and COG  databases. Genes with signal peptides were predicted using SignalP . Genes with transmembrane helices were predicted using TMHMM .
The genome of K. oryzae Ola 51T contains one circular chromosome (Fig. 3). The chromosome contains 5,303,342 nucleotides with 54.0% G+C content. The genome contains 4,926 predicted genes, 4773 protein-coding genes, 105 RNA genes (16 rRNA genes, 76 tRNA genes, and 13 ncRNA genes), 48 pseudo genes, and 1 CRISPR repeats. Among the 4,773 protein-coding genes, 3,765 genes (78.88%) have been assigned functions, while 1008 genes (21.12%) have been annotated as hypothetical or unknown proteins (Table 3). The distribution of genes into COG functional categories is presented in Table 4 and Fig. 3.
The genome sequences of K. cowanii JCM 10956 T, K. radicincitans DSM 16656 T (=D5/23T) , K. radicincitans UMEnt01/12 , K. radicincitans YD4 , K. sacchari SP1T , “K. pseudosacchari” JM-387T , K. oryzae KO348 , and Enterobacter sp. R4-368  which was close to K. sacchari SP1T  had been deposited in the GenBank database. The genome ANIs (Additional file 1: Table S1) between Ola 51T and the other strains belonging to the genus Kosakonia were calculated using the Orthologous Average Nucleotide Identity tool . The cut-off ANI value for species boundary was set at 95% - 96% . The ANI value (95.85%) between K. oryzae Ola 51T and K. radicincitans DSM 16656 T is in the fuzzy zone 95% - 96%. The digital DDH value between Ola 51T and DSM 16656 T calculated by the Genome-to-Genome Distance Calculator  with the Formula 2 is 66.2%, below the 70% cut-off value for species boundary. Moreover, Ola 51T and DSM 16656 T were differentiated by metabolic phenotypes [3, 11] and ribosomal protein mass profiles . Therefore, K. oryzae and K. radicincitans are closely related sister species.
Strain YD4 was closer to K. radicincitans DSM 16656 T than K. oryzae Ola 51T on the phylogenetic tree based on the 16S rRNA genes (Fig. 2). However, the ANI value and the digital DDH value between YD4 and K. radicincitans DSM 16656 T is 95.56% and 64.4%, respectively, while between YD4 and K. oryzae Ola 51T is 97.04% and 74.3%, respectively. Therefore, the strain YD4 belongs to K. oryzae but not K. radicincitans.
Strain KO348 was grouped with K. sacchari SP1T, Enterobacter sp. R4-368, and “K. pseudosacchari” JM-387T on the phylogenetic tree based on the 16S rRNA genes (Fig. 2). The ANI value between KO348 and K. oryzae Ola 51T is 84.04%. The strain KO348 thus does not belong to K. oryzae. The ANI value between KO348 and Enterobacter sp. R4-368 , K. sacchari SP1T, or “K. pseudosacchari” JM-387T is 98.80%, 94.56%, or 94.05%, respectively. Therefore, KO348 and R4-368 belong to the same species, likely a novel species closely related to K. sacchari and “K. pseudosacchari”.
K. oryzae Ola 51T and YD4, K. radicincitans DSM 16656 T and UMEnt01/12, K. sacchari SP1T, “K. pseudosacchari” JM-387T, and Kosakonia sp. KO348 and R4-368 were all isolated from plants. Their genomes contain genes encoding multiple enzymes degrading plant cell wall polysaccharides and removing reactive oxygen species, likely facilitating endophytic colonization . They all contain genes encoding the regulatory protein (Fha1) and structural proteins (Lip, IcmF, DotU and ClpV) and secreted proteins (VgrG and Hcp) of the type VI secretion system, which may play a role in the plant-associated lifestyle . Except K. radicincitans DSM 16656 T and UMEnt01/12, these strains contain the most structural proteins (YscCJRSTUVN) of the type III secretion system, which is not widespread among the previously studied endophytic bacteria .
These plant-associated Kosakonia strains contain genes contributing to multiple plant growth-promoting activities. They all contain the nif gene cluster (nifJHDKTYENXUSVWZMFLABQ) for the Mo-Fe nitrogenase-dependent nitrogen fixation, the genes encoding indole-3-acetaldehyde dehydrogenase, aspartate aminotransferase, aromatic amino acid aminotransferase and phenylpyruvate decarboxylase for producing the phytohormone auxin, and the budABC genes for producing volatile acetoin and 2,3-butanediol which induce plant systemic resistance to pathogens . In addition, K. oryzae Ola 51T and YD4, and K. radicincitans DSM 16656 T and UMEnt01/12 also contain the anf gene cluster (anfHDGK) for the Fe-Fe nitrogenase-dependent nitrogen fixation. In contrast, the clinical strain K. cowanii JCM 10956 T does not contain the nif gene cluster.
The phylogeny of the members of the genus Kosakonia based on the 16S rRNA gene sequences is roughly in agreement with their overall genome relatedness. The complete genome sequence of K. oryzae Ola 51T provides the reference genome for genomic identification of strains belonging to K. oryzae. Analyses of the overall genome relatedness indices (ANI and digital DDH values), easily and reliably show that K. oryzae and K. radicincitans are closely related sister species and that the strain YD4, which shows close 16S rRNA gene-based phylogeney to K. radicincitans and was classified into K. radicincitans, belongs to K. oryzae. As well as YD4, which is able to promote growth of the yerba mate plants in low-fertility soils , K. oryzae Ola 51T contains both the nif gene cluster and the anf gene cluster for nitrogen fixation and genes contributing to production of auxin and volatile acetoin and 2,3-butanediol. Therefore, K. oryzae Ola 51T may be able to promote plant growth. Genomic analyses also show that K. oryzae Ola 51T and YD4 may have the type III and VI secretion systems and thus motivate us to study the functions of the type III and VI secretion systems in the interactions between beneficial Kosakonia bacteria and plants.
This work was supported by the National Natural Science Foundation of China (31171504, 31370052 and 31471449), Zhejiang Provincial Natural Science Foundation of China (LY14C010002), Science and Technology Planning Project of Guangdong Province (2014A030313459 and 2014A050503058), and the State Key Laboratory of Rice Biology, China.
YL, SL and MC assembled the sequencing data and completed the genome analysis; GP did the microbiological studies and obtained the organism information; ZT and QA designed the study and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.(39K, doc) Average nucleotide identities (ANIs) between genomes of the strains belonging to the genus Kosakonia
. (DOC 38 kb)
Average nucleotide identities (ANIs) between genomes of the strains belonging to the genus Kosakonia . (DOC 38 kb)