Genetic analysis of housekeeping genes of members of the genus Acholeplasma: phylogeny and complementary molecular markers to the 16S rRNA gene

The partial nucleotide sequences of the rpoB and gyrB genes as well as the complete sequence of the 16S-23S rRNA intergenic transcribed spacer (ITS) were determined for all known Acholeplasma species. The same genes of Mesoplasma and Entomoplasma species were also sequenced and used to infer phylogenetic relationships among the species within the orders Entomoplasmatales and Acholeplasmatales. The comparison of the ITS, rpoB, and gyrB phylogenetic trees with the 16S rRNA phylogenetic tree revealed a similar branch topology suggesting that the ITS, rpoB, and gyrB could be useful complementary phylogenetic markers for investigation of evolutionary relationships among Acholeplasma species. Thus, the multilocus phylogenetic analysis of Acholeplasma multilocale sequence data (ATCC 49900 (T) = PN525 (NCTC 11723)) strongly indicated that this organism is most closely related to the genera Mesoplasma and Entomoplasma (family Entomoplasmataceae) and form the branch with Mesoplasma seiffertii, Mesoplasma syrphidae, and Mesoplasma photuris. The closest genetic relatedness of this species to the order Entomoplasmatales was additionally supported by the finding that A. multilocale uses UGA as the tryptophan codon in its gyrB and gyrA sequences. Use of the UGA codon for encoding tryptophan was previously reported as a unique genetic feature of Entomoplasmatales and Mycoplasmatales but not of Acholeplasmatales. These data, as well as previously published data on metabolic features of A. multilocale, leads to the proposal to reclassify A. multilocale as a member of the family Entomoplasmataceae.     

Use of rpoB gene analysis for identification of nitrogen-fixing Paenibacillus species as an alternative to the 16S rRNA gene

AIM: To avoid the limitations of 16S rRNA-based phylogenetic analysis for Paenibacillus species, the usefulness of the RNA polymerase beta-subunit encoding gene (rpoB) was investigated as an alternative to the 16S rRNA gene for taxonomic studies. METHODS AND RESULTS: Partial rpoB sequences were generated for the type strains of eight nitrogen-fixing Paenibacillus species. The presence of only one copy of rpoB in the genome of P. graminis strain RSA19(T) was demonstrated by denaturing gradient gel electrophoresis and hybridization assays. A comparative analysis of the sequences of the 16S rRNA and rpoB genes was performed and the eight species showed between 91.6-99.1% (16S rRNA) and 77.9-97.3% (rpoB) similarity, allowing a more accurate discrimination between the different species using the rpoB gene. Finally, 24 isolates from the rhizosphere of different cultivars of maize previously identified as Paenibacillus spp. were assigned correctly to one of the nitrogen-fixing species. CONCLUSIONS, SIGNIFICANCE AND IMPACT OF THE STUDY: The data obtained in this study indicate that rpoB is a powerful identification tool, which can be used for the correct discrimination of the nitrogen-fixing species of agricultural and industrial importance within the genus Paenibacillus.     

16S rRNA gene and 16S-23S rRNA gene internal transcribed spacer sequences analysis of the genus Myxococcus.

Phylogenetic relationships of the species belonging to the genus Myxococcus were elucidated based on the sequences of 16S rRNA genes and 16S-23S rRNA gene internal transcribed spacer (ITS) regions. The Myxococcus species were consequently classified into four distinct groups. The type strain of Myxococcus coralloides occupied an independent position (Group 1); it has been recently reclassified as Corallococcus coralloides. Group 2 comprised the type strains of both Myxococcus virescens and Myxococcus xanthus, and some strains assigned to Myxococcus flavescens. The type strain of M. flavescens was contained in Group 3 along with the strains of Myxococcus fulvus. Group 4 included the strains belonging to C. coralloides, M. fulvus, and M. stipitatus. The type strain of M. fulvus that was allocated outside Group 4 in the 16S rRNA gene tree belonged to Group 3 in the ITS tree. These results strongly suggest that the morphological characteristics of Myxococcus species are not consistent with the phylogenetic relationships. The Myxococcus species must therefore be redefined according to the phylogenetic relationships revealed in this study.     

Intrageneric structure of the genus Gluconobacter analyzed by the 16S rRNA gene and 16S-23S rRNA gene internal transcribed spacer sequences

Forty-nine strains belonging to the genus Gluconobacter were re-examined with respect to their species identification based on the sequences of the 16S rDNA and 16S-23S rDNA internal transcribed spacer regions (ITS). A phylogenetic tree constructed from the 16S rDNA sequences indicated the presence of five clusters corresponding, respectively, to the major five species of the genus Gluconobacter, namely G. albidus, G. cerinus, G. frateurii, G. oxydans (type species), and G. thailandicus. The type strain of G. asaii, NBRC 3276T (T=type strain) was included in the G. cerinus cluster, which is consistent with the report that G. asaii is a junior subjective synonym of G. cerinus. Existence of the G. albidus, G. cerinus, G. frateurii, G. oxydans, and G. thailandicus clusters was also recognized by the ITS sequence analysis. Both sequence analyses revealed that the G. cerinus and G. frateurii clusters were heterogeneous. The G. cerinus cluster comprised three strains of G. cerinus and one strain of G. frateurii, while the G. frateurii cluster included ten strains of G. frateurii, three of G. cerinus, and eleven of G. oxydans. These results suggest that phenotypic differences among Gluconobacter species are ambiguous and the species definition must be re-evaluated. The 16S rDNA and ITS sequences determined in this study are valuable for the identification and phylogenetic analysis of Gluconobacter species.     

Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such "heterogeneity hot spots" occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Variations in the 16S-23S rRNA internal transcribed spacer of fibrolytic Butyrivibrio isolates from the reindeer rumen

Strains of Butyrivibrio are principal cellulytic bacteria in the rumen of the High Arctic Svalbard reindeer ( Rangifer tarandus platyrhynchus ). According to phylogenetic analysis based on 16S rRNA gene sequencing, Butyrivibrio can be divided into three subgroups within the Clostridia class of the phylum Firmicutes, but the current phenotypic and genotypic differentiation within the family Lachnospiraceae is insufficient. This current study describes the sequence diversity of the 16S-23S rRNA intergenic transcribed spacer (ITS) region of Butyrivibrio isolates from reindeer. A total of 17 different ITS sequences with sizes between 449 and 784 nt were obtained. Genes encoding tRNA(Ile) and tRNA(Ala) were identified in four of the sequences. Phylogenetic neighbor-joining trees were constructed based on the ITS sequence and compared with a phylogenetic neighbor-joining tree based on 16S rRNA gene sequences previously obtained for the same isolates. These comparisons indicated a better differentiation between strains in the ITS sequence than the 16S rRNA gene based tree. Through this study, a better means for identifying and tracking fibrolytic and potentially probiotic Butyrivibrio strains in reindeer and other ruminants has been provided.     

拟诺卡氏菌16S rRNA,gyrB,sod和rpoB基因的系统发育分析

为了更好地了解拟诺卡氏菌属(Nocardiopsis)各物种间的系统发育关系,该属现有有效描述种的gyrB,sod和rpoB基因的部分序列被测定,结合16S rRNA基因,对拟诺卡氏菌属进行了系统发育重建。研究发现拟诺卡氏菌属gyrB,sod和rpoB基因的平均相似性分别为87.7%、87.3%和94.1%,而16S rRNA基因的平均相似性则达到96.65%,3个看家基因均比16S rRNA具有更高的分歧度。比较基于不同基因的系统树发现,由gyrB基因得到的系统树拓扑结构与16S rRNA得到的结构在亚群上基本一致。因此,gyrB基因在拟诺卡氏菌属的系统分类上比16S rRNA基因更具优越性。     

Comparison of the Biolog OmniLog Identification System and 16S ribosomal RNA gene sequencing for accuracy in identification of atypical bacteria of clinical origin

The Biolog OmniLog Identification System (Biolog) and the 16S ribosomal RNA (rRNA) gene sequencing methods were compared to conventional microbiological methods and evaluated for accuracy of bacterial identification. These methods were evaluated using 159 clinical isolates. Each isolate was initially identified by conventional biochemical tests and morphological characteristics and subsequently placed into one of seven categories: aerobic Actinomycetes, Bacillus, Coryneforms, fastidious Gram-negative rods (GNR), non-fermenting GNR, miscellaneous Gram-positive rods (GPR), and Vibrio/Aeromonas. After comparison to the conventional identification, the Biolog system and 16S rRNA gene sequence identifications were classified as follows: a) correct to the genus and species levels; b) correct to the genus level only; or c) neither (unacceptable) identification. Overall, 16S rRNA gene sequencing had the highest percent accuracy with 90.6% correct identifications, while the Biolog system identified 68.3% of the isolates correctly. For each category, 16S rRNA gene sequencing had a substantially higher percent accuracy compared to the conventional methods. It was determined that the Biolog system is deficient when identifying organisms in the fastidious GNR category (20.0%). The observed data suggest that 16S rRNA gene sequencing provides a more accurate identification of atypical bacteria than the Biolog system.     

rpoB sequence analysis as a novel basis for bacterial identification

Comparison of the sequences of conserved genes, most commonly those encoding 16S rRNA, is used for bacterial genotypic identification. Among some taxa, such as the Enterobacteriaceae, variation within this gene does not allow confident species identification. We investigated the usefulness of RNA polymerase beta-subunit encoding gene ( rpoB  ) sequences as an alternative tool for universal bacterial genotypic identification. We generated a database of partial rpoB for 14 Enterobacteriaceae species and then assessed the intra- and interspecies divergence between the rpoB and the 16S rRNA genes by pairwise comparisons. We found that levels of divergence between the rpoB sequences of different strains were markedly higher than those between their 16S rRNA genes. This higher discriminatory power was further confirmed by assigning 20 blindly selected clinical isolates to the correct enteric species on the basis of rpoB sequence comparison. Comparison of rpoB sequences from Enterobacteriaceae was also used as the basis for their phylogenetic analysis and demonstrated the genus Klebsiella to be polyphyletic. The trees obtained with rpoB were more compatible with the currently accepted classification of Enterobacteriaceae than those obtained with 16S rRNA. These data indicate that rpoB is a powerful identification tool, which may be useful for universal bacterial identification.     

rpoB-based microbial community analysis avoids limitations inherent in 16S rRNA gene intraspecies heterogeneity

Contemporary microbial community analysis frequently involves PCR-amplified sequences of the 16S rRNA gene (rDNA). However, this technology carries the inherent problem of heterogeneity between copies of the 16S rDNA in many species. As an alternative to 16S rDNA sequences in community analysis, we employed the gene for the RNA polymerase beta subunit (rpoB), which appears to exist in one copy only in bacteria. In the present study, the frequency of 16S rDNA heterogeneity in bacteria isolated from the marine environment was assessed using bacterial isolates from the red alga Delisea pulchra and from the surface of a marine rock. Ten strains commonly used in our laboratory were also assessed for the degree of heterogeneity between the copies of 16S rDNA and were used to illustrate the effect of this heterogeneity on microbial community pattern analysis. The rock isolates and the laboratory strains were also used to confirm nonheterogeneity of rpoB, as well as to investigate the versatility of the primers. In addition, a comparison between 16S rDNA and rpoB PCR-DGGE (denaturing gradient gel electrophoresis)-based community analyses was performed using a DNA mixture of nine isolates from D. pulchra. Eight out of 14 isolates from D. pulchra, all rock isolates, and 6 of 10 laboratory strains displayed multiple bands for 16S rDNA when analyzed by DGGE. There was no indication of heterogeneity for either the rock isolates or the laboratory strains when rpoB was used for PCR-DGGE analysis. Microbial community pattern analysis using 16S rDNA PCR-DGGE showed an overestimation of the number of laboratory strains in the sample, while some strains were not represented. Therefore, the 16S rDNA PCR-DGGE-based community analysis was proven to be severely limited by 16S rDNA heterogeneity. The mixture of isolates from D. pulchra proved to be more accurately described using rpoB, compared to the 16S rDNA-based PCR-DGGE.     

The phylogeny of the genus Nitrobacter based on comparative rep-PCR, 16S rRNA and nitrite oxidoreductase gene sequence analysis

Strains of Nitrobacter mediate the second step in the nitrification process by oxidizing nitrite to nitrate. The phylogenetic diversity of the genus is currently not well investigated. In this study, a rep-PCR profile and the nearly complete 16S rRNA gene sequence of 30 strains, comprising a wide physiological as well as ecological diversity and encompassing representatives of the four species, were determined. The sequence diversity of the 16S rRNA gene between different species was low, indicating the need for additional phylogenetic markers. Therefore, primers were developed for amplifying the complete nxrX gene and a 380bp fragment of the nxrB1 gene, which are both genes involved in the nitrite oxidation process. These genes confirmed the division into phylogenetic groups revealed by the 16S rRNA gene but showed a better discriminatory power. They can be a valuable additional tool for phylogenetic analysis within the genus Nitrobacter and can assist in the identification of new Nitrobacter isolates.     

基于16S rRNA和HSP60基因序列分析粘细菌孢囊杆菌亚目形态分类的种属之间的亲缘关系

[目的]利用16S rRNA和HSP60基因分子标记分析鉴定形态分类特征不稳定的粘细菌种属.[方法]利用粘细菌的传统分离纯化方法从土壤中分离粘细菌,根据菌株的形态特征进行分类,PCR方法扩增菌株的16S rRNA和HSP60基因序列并进行系统发育关系分析.[结果]根据形态特征,分离得到的15株粘细菌菌株归入孢囊杆菌亚目(Cystobacterineae)的2个科3个属.其中11株粘细菌具有典型的所在种属的子实体结构,而菌株0085-4、0121-3、NM03和Myx9736的子实体结构发生了不同程度退化.15株粘细菌的16S rRNA基因序列的相似性在95.4%到99.5%之间.而HSP60基因序列差异较大.[结论]在属水平上,粘细菌形态分类特征和16S rRNA基因系统进化关系具有很好的一致性;在揭示粘细菌种间系统发育关系中,HSP60基因序列更为适用.     

Discordant 16S and 23S rRNA gene phylogenies for the genus Helicobacter: implications for phylogenetic inference and systematics

Analysis of 16S rRNA gene sequences has become the primary method for determining prokaryotic phylogeny. Phylogeny is currently the basis for prokaryotic systematics. Therefore, the validity of 16S rRNA gene-based phylogenetic analyses is of fundamental importance for prokaryotic systematics. Discrepancies between 16S rRNA gene analyses and DNA-DNA hybridization and phenotypic analyses have been noted in the genus Helicobacter. To clarify these discrepancies, we sequenced the 23S rRNA genes for 55 helicobacter strains representing 41 taxa (>2,700 bases per sequence). Phylogenetic-tree construction using neighbor-joining, parsimony, and maximum likelihood methods for 23S rRNA gene sequence data yielded stable trees which were consistent with other phenotypic and genotypic methods. The 16S rRNA gene sequence-derived trees were discordant with the 23S rRNA gene trees and other data. Discrepant 16S rRNA gene sequence data for the helicobacters are consistent with the horizontal transfer of 16S rRNA gene fragments and the creation of mosaic molecules with loss of phylogenetic information. These results suggest that taxonomic decisions must be supported by other phylogenetically informative macromolecules, such as the 23S rRNA gene, when 16S rRNA gene-derived phylogeny is discordant with other credible phenotypic and genotypic methods. This study found Wolinella succinogenes to branch with the unsheathed-flagellum cluster of helicobacters by 23S rRNA gene analyses and whole-genome comparisons. This study also found intervening sequences (IVSs) in the 23S rRNA genes of strains of 12 Helicobacter species. IVSs were found in helices 10, 25, and 45, as well as between helices 31' and 27'. Simultaneous insertion of IVSs at three sites was found in H. mesocricetorum.     

基于12S,16S和28S rRNA基因序列的中国小毛瓢虫属系统发育分析

【目的】小毛瓢虫属Scymnus Kugelann昆虫主要捕食蚜虫、蚧虫等害虫,是一类经济上重要的天敌昆虫。目前针对小毛瓢虫属的系统发育研究尚属空白,亚属之间的系统演化关系尚不明确,为了建立合理的分类系统,亟需对小毛瓢虫属的亲缘关系进行研究和探讨。【方法】以华南农业大学馆藏的小毛瓢虫属5亚属共44种为研究对象,采用PCR技术对12S, 16S和28S rRNA基因的部分序列进行扩增;运用MEGA 7.0分析了小毛瓢虫属内12S, 16S和28S rRNA基因的碱基组成,基于K2P模型计算了小毛瓢虫属44种的种间遗传距离;采用最大似然法(maximum-likelihood, ML)和贝叶斯推断法(Bayesian-inference, BI)构建该属的系统发育树。【结果】扩增获得小毛瓢虫属44种的12S rRNA基因序列平均长度为356 bp, 16S rRNA基因序列平均长度为351 bp, 28S rRNA基因序列平均长度为315 bp;序列分析表明,12S rRNA基因的A, T, G和C平均含量分别为38.8%, 43.5%, 11.9%和5.8%, 16S rRNA基因的A, T, G和C平均含量分别为37.6%, 40.3%, 14.4%和7.7%, 28S rRNA基因的A, T, G和C平均含量分别为26.7%, 18.3%, 31.4%和23.5%;基于联合序列分析的种间遗传距离为0.004～0.276,平均遗传距离为0.115。系统发育分析结果表明,小毛瓢虫属为单系起源,而小毛瓢虫亚属Scymnus(Scymnus) Kugelann、毛瓢虫亚属Scymnus(Neopullus) Sasaji、小瓢虫亚属Scymnus(Pullus) Mulsant和拟小瓢虫亚属Scymnus(Parapullus) Yang均为并系起源。【结论】基于12S, 16S和28S rRNA基因序列的小毛瓢虫属系统发育分析显示传统的形态学分类体系与基于分子数据分析的结果部分不一致,这表明应该对该属内各亚属的鉴别特征进行全面检视,筛选并确立各亚属的形态指标,同时也表明该属内的亚属分类单元需重新厘定。     

Phylogenetic relationship of 16 Oedipodidae species （Insecta： Orthoptera） based on the 16S rRNA gene sequences

The sequences of the mitochondrial 16S rRNA gene of 16 Oedipodidae species were amplified and sequenced. All sequences were aligned and analyzed and the phylogenetic relationships were inferred. The properties of 16S gene in Oedipodidae showed typical patterns of many insects such as a high A＋T content and variable distance-dependent transition/transversion ratios. The 0.2 weight for sites of loops may be advisable for phylogeny reconstruction using the maximum parsimony method. The phylogenetic analysis results do not support the current subfamily classification systems of Oedipodidae. Bryodemellinae and Bryodeminae are closely related and should be merged as one subfamily. The status of Oedipodinae and Locustinae is also problematic.     

基于16S rRNA和rpoB基因的分子系统发育分析在铜绿假单胞菌鉴定中的应用

为对比16S rRNA和rpo B基因分子系统发育分析与传统表型分类法对铜绿假单胞菌的鉴定,评估16S rRNA和rpo B基因序列分析在铜绿假单胞菌鉴定中的应用,用表型分类方法对临床自动微生物鉴定系统鉴定为铜绿假单胞菌的23株分离株进行再鉴定,PCR扩增23株分离株16S rRNA和rpo B基因片段,并测序进行系统发育分析。结果表明,表型再鉴定结果与自动微生物鉴定系统鉴定结果一致。基于两个基因的系统发育分析均显示分离株p22与不动杆菌属序列聚为一枝,其余22株分离株与铜绿假单胞菌序列聚为一枝。因此p22应鉴定为不动杆菌,16S rRNA和rpo B基因序列分析均能准确鉴定铜绿假单胞菌并能较好建立假单胞菌属内种间关系。     

Leptospira spp. strain identification by MALDI TOF MS is an equivalent tool to 16S rRNA gene sequencing and multi locus sequence typing (MLST)

ABSTRACT: BACKGROUND: In this study mass spectrometry was used for evaluating extracted leptospiral protein samples and results were compared with molecular typing methods. For this, an extraction protocol for Leptospira spp. was independently established in two separate laboratories. Reference spectra were created with 28 leptospiral strains, including pathogenic, non-pathogenic and intermediate strains. This set of spectra was then evaluated on the basis of measurements with well-defined, cultured leptospiral strains and with 16 field isolates of veterinary or human origin. To verify discriminating peaks for the applied pathogenic strains, statistical analysis of the protein spectra was performed using the software tool ClinProTools. In addition, a dendrogram of the reference spectra was compared with phylogenetic trees of the 16S rRNA gene sequences and multi locus sequence typing (MLST) analysis. RESULTS: Defined and reproducible protein spectra using MALDI-TOF MS were obtained for all leptospiral strains. Evaluation of the newly-built reference spectra database allowed reproducible identification at the species level for the defined leptospiral strains and the field isolates. Statistical analysis of three pathogenic genomospecies revealed peak differences at the species level and for certain serovars analyzed in this study. Specific peak patterns were reproducibly detected for the serovars Tarassovi, Saxkoebing, Pomona, Copenhageni, Australis, Icterohaemorrhagiae and Grippotyphosa. Analysis of the dendrograms of the MLST data, the 16S rRNA sequencing, and the MALDI-TOF MS reference spectra showed comparable clustering. CONCLUSIONS: MALDI-TOF MS analysis is a fast and reliable method for species identification, although Leptospira organisms need to be produced in a time-consuming culture process. All leptospiral strains were identified, at least at the species level, using our described extraction protocol. Statistical analysis of the three genomospecies L. borgpetersenii, L. interrogans and L. kirschneri revealed distinctive, reproducible differentiating peaks for seven leptospiral strains which represent seven serovars. Results obtained by MALDI-TOF MS were confirmed by MLST and 16S rRNA gene sequencing.     

Identification of waterborne bacteria by the analysis of 16S-23S rRNA intergenic spacer region

AIM: In this study, we evaluated, the use of universal primers, specific for the 16S-23S rRNA intergenic region, to detect and identify nine species that are of high interest for the microbiological control of water. METHODS AND RESULTS: The analysis of the fragments was carried out using a High Resolution acrylamide/bisacrylamide gels in a fluorescent automated DNA sequencer. The results showed specific profiles for each of the nine species but this technique failed to detect simultaneously micro-organisms in samples containing a mixed population. CONCLUSION: Nevertheless, the electrophoretic profiles obtained provided a very useful tool for the rapid and specific identification of water isolates. SIGNIFICANCE AND IMPACT OF THE STUDY: A possible new methodology for a rapid identification of pathogens in water.     

Yersinia spp. Identification Using Copy Diversity in the Chromosomal 16S rRNA Gene Sequence

API 20E strip test, the standard for Enterobacteriaceae identification, is not sufficient to discriminate some Yersinia species for some unstable biochemical reactions and the same biochemical profile presented in some species, e.g. Yersinia ferderiksenii and Yersinia intermedia, which need a variety of molecular biology methods as auxiliaries for identification. The 16S rRNA gene is considered a valuable tool for assigning bacterial strains to species. However, the resolution of the 16S rRNA gene may be insufficient for discrimination because of the high similarity of sequences between some species and heterogeneity within copies at the intra-genomic level. In this study, for each strain we randomly selected five 16S rRNA gene clones from 768 Yersinia strains, and collected 3,840 sequences of the 16S rRNA gene from 10 species, which were divided into 439 patterns. The similarity among the five clones of 16S rRNA gene is over 99% for most strains. Identical sequences were found in strains of different species. A phylogenetic tree was constructed using the five 16S rRNA gene sequences for each strain where the phylogenetic classifications are consistent with biochemical tests; and species that are difficult to identify by biochemical phenotype can be differentiated. Most Yersinia strains form distinct groups within each species. However Yersinia kristensenii, a heterogeneous species, clusters with some Yersinia enterocolitica and Yersinia ferderiksenii/intermedia strains, while not affecting the overall efficiency of this species classification. In conclusion, through analysis derived from integrated information from multiple 16S rRNA gene sequences, the discrimination ability of Yersinia species is improved using our method.