Molecular phylogenies of Parabasalia inferred from four protein genes and comparison with rRNA trees

The molecular phylogeny of parabasalids has mainly been inferred from small subunit (SSU) rRNA sequences and has conflicted substantially with systematics based on morphological and ultrastructural characters. This raises the important question, how congruent are protein and SSU rRNA trees? New sequences from seven diverse parabasalids (six trichomonads and one hypermastigid) were added to data sets of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), enolase, alpha-tubulin and beta-tubulin and used to construct phylogenetic trees. The GAPDH tree was well resolved and identical in topology to the SSU rRNA tree. This both validates the rRNA tree and suggests that GAPDH should be a valuable tool in further phylogenetic studies of parabasalids. In particular, the GAPDH tree confirmed the polyphyly of Monocercomonadidae and Trichomonadidae and the basal position of Trichonympha agilis among parabasalids. Moreover, GAPDH strengthened the hypothesis of secondary loss of cytoskeletal structures in Monocercomonadidae such as Monocercomonas and Hypotrichomonas. In contrast to GAPDH, the enolase and both tubulin trees are poorly resolved and rather uninformative about parabasalian phylogeny, although two of these trees also identify T. agilis as representing the basal-most lineage of parabasalids. Although all four protein genes show multiple gene duplications (for 3-6 of the seven taxa examined), most duplications appear to be relatively recent (i.e., species-specific) and not a problem for phylogeny reconstruction. Only for enolase are there more ancient duplications that may confound phylogenetic interpretation.     

Phylogenetic affiliation of Aeromonas culicicola MTCC 3249(T) based on gyrB gene sequence and PCR-amplicon sequence analysis of cytolytic enterotoxin gene

We determined the gyrB gene sequences of all 17 hybridizations groups of Aeromonas. Phylogenetic trees showing the evolutionary relatedness of gyrB and 16S rRNA genes in the type strains of Aeromonas were compared. Using this approach, we determined the phylogenetic position of Aeromonas culicicola MTCC 3249(T), isolated from midgut of Culex quinquefasciatus. In the gyrB based-analysis A. culicicola MTCC 3249(T) grouped with A. veronii whereas, it grouped with A. jandaei in the 16S rRNA based tree. The number of nucleotide differences in 16S rRNA sequences was less than found with the gyrB sequence data. Most of the observed nucleotide differences in the gyrB gene were synonymous. The Cophenetic Correlation Coefficient (CCC) for gyrB sequences was 0.87 indicating this gene to be a better molecular chronometer compared to 16S rRNA for delineation of Aeromonas species. This strain was found to be positive for the cytolytic enterotoxin gene. PCR-Amplicon Sequence Analysis (PCR-ASA) of this gene showed that the isolate is affiliated to type I and is potentially pathogenic. These PCR-ASA results agreed in part with the gyrB sequence results.     

大黄鱼16SrRNA和18SrRNA基因的测定与序列分析

利用多对引物,扩增并测定出大黄鱼16SrRNA基因和18SrRNA基因的部分序列,其长度分别为1202bp和1275bp,16SrRNA基因序列的GC含量为46．12％,18SrRNA基因的Gc含量为53．oo％。将大黄鱼16SrRNA基因序列与GenBank中15种硬骨鱼类的同源序列结合,同时将其18SrRNA基因序列与GenBank中9种脊索动物的同源序列相结合,运用软件获得各自序列间差异百分比,转换和颠换数值等信息。基于这两种基因序列,利用NJ法和BI法,分别构建16种硬骨鱼类和10种脊索动物的分子系统树。18SrRNA构建的系统树包括三大支,一支为哺乳类、鸟类和爬行类共6个物种,一支为两栖类的1个物种,另一支为2种硬骨鱼类。16SrRNA构建的系统树显示大黄鱼所在的石首鱼科与鲈科和盖刺鱼科亲缘关系较近。此外还讨论了这两个基因的序列特征。     

The phylogenetic relationships of cyanobacteria inferred from 16S rRNA,gyrB, rpoC1 and rpoD1 gene sequences

Phylogenetic analysis of cyanobacteria was carried out using the small subunit rRNA (16S rRNA), DNA gyrase subunit B (gyrB), DNA-dependent RNA polymerase gamma subunit (rpoC1) and a principal sigma factor of E. coli sigma(70) type for DNA-dependent RNA polymerase (rpoD1) gene sequences of 24 strains which contained 5 subgroups of cyanobacteria-3 strains of the Chroococcales, 5 strains of the Pluerocapsales, 7 strains of the Oscillatoriales, 7 strains of the Nostocales and 2 strains of the Stigonematales. Degenerated PCR primers of gyrB, rpoC1 and rpoD1 genes were designed using consensus amino acid sequences registered in GenBank. The phylogenetic positions of cyanobacteria were resolved through phylogenetic analysis based on 16S rDNA, gyrB, rpoC1 and rpoD1 gene sequences. Phylogenies of gyrB, rpoC1 and rpoD1 support 16S rRNA-based classification of cyanobacteria. Interestingly, phylogenies from amino acid sequences deduced from gyrB and combined amino acid sequences deduced from rpoC1 and rpoD1 genes strongly support that of 16S rRNA, but the branching pattens of the trees based on 16S rDNA, GyrB, rpoC1, rpoD1 and combined amino acid sequences deduced from rpoC1 and rpoD1 were not congruent. In this study, we showed the correlation among phylogenetic relationships of 16S rDNA, gyrB, rpoC1 and rpoD1 genes. The phylogenetic trees based on the sequences of 16S rDNA, GyrB, rpoC1, rpoD1 and the combined amino acid sequences deduced from rpoC1 and rpoD1 showed that the lateral gene transfer of rRNA might be suspected for Synechocystis sp. PCC 6803.     

Sequence heterogeneity in the two 16S rRNA genes of Phormium yellow leaf phytoplasma.

Phormium yellow leaf (PYL) phytoplasma causes a lethal disease of the monocotyledon, New Zealand flax (Phormium tenax). The 16S rRNA genes of PYL phytoplasma were amplified from infected flax by PCR and cloned, and the nucleotide sequences were determined. DNA sequencing and Southern hybridization analysis of genomic DNA indicated the presence of two copies of the 16S rRNA gene. The two 16S rRNA genes exhibited sequence heterogeneity in 4 nucleotide positions and could be distinguished by the restriction enzymes BpmI and BsrI. This is the first record in which sequence heterogeneity in the 16S rRNA genes of a phytoplasma has been determined by sequence analysis. A phylogenetic tree based on 16S rRNA gene sequences showed that PYL phytoplasma is most closely related to the stolbur and German grapevine yellows phytoplasmas, which form the stolbur subgroup of the aster yellows group. This phylogenetic position of PYL phytoplasma was supported by 16S/23S spacer region sequence data.     

The phylogenetic significance of peptidoglycan types: Molecular analysis of the genera Microbacterium and Aureobacterium based upon sequence comparison of gyrB, rpoB, recA and ppk and 16SrRNA genes

The type strains of 27 species of the genus Microbacterium, family Microbacteriaceae, were analyzed with respect to the phylogeny of the housekeeping genes coding for DNA gyrase subunit B (gyrB), RNA-polymerase subunit B (rpoB), recombinase A (recA) and polyphosphate kinase (ppk). The resulting gene trees were compared to the 16S rRNA gene phylogeny of the same species. The topology of neighbour-joining and maximum parsimony phylogenetic trees based upon nucleic acid sequences and protein sequences of housekeeping genes differed among each other and no gene tree was identical to that of the 16S rRNA gene tree. Only some species showed consistent clustering by all genes analyzed, but the majority of species branched with different neighbours in most gene trees. The failure to phylogenetically cluster type strains into two groups based upon differences in the amino acid composition of peptidoglycan on the basis of 16S rRNA gene sequence similarity, once leading to the union of the genera Microbacterium and Aureobacterium, was also seen in the analysis of recA, rpoB and gyrB gene and protein phylogenies. Analysis of the pkk gene and protein as well as of a concatenate tree, combining sequences of all five genes (total of 3.700 nucleotides), sees members of the former genus Aureobacterium and other type strains with lysine as diagnostic diamino acid to form a coherent cluster that branches within the radiation of Microbacterium species with ornithine in the peptidoglycan.     

Complete sequence of the Japanese quail (Coturnix japonica) mitochondrial genome and its genetic relationship with related species

The Japanese quail (Coturnix japonica; JQ) is one of the domesticated fowl species of Japan. To provide DNA sequence information for examination of its phylogenetic position in the order Galliformes, the complete sequence of the JQ mitochondria was determined. Sequence analysis revealed that the JQ mitochondrial genome is a circular DNA of 16 697 basepairs (bp), which is smaller than the chicken mitochondrial DNA of 16 775 bp, but the genomic structure of JQ mitochondria was the same as that of the chicken. The sequence homologies of all mitochondrial genes including those for 12S and 16S ribosomal RNA (rRNA), between Japanese quail and chicken ranged from 78.0 to 89.9%. Because the sequences of NADH dehydrogenase subunit 2 and cytochrome b genes had been reported in five species [Phasianus colchicus (ring-neck pheasant: RP), Gallus gallus domesticus (chicken: CH), Perdix perdix (grey partridge: GP), Bambusicola thoracia (Chinese bamboo partridge: CP), and Aythya americana (redhead: RH)], the concatenated nucleotide sequences (2184 bp) and amino acid sequences of these two genes were used in a phylogenetic analysis of JQ against these five species using a maximum likelihood (ML) method. Using the first and second bases of the codons, and the third base of the codons indicated a phylogenic tree of [RH, (RP, GP), (JQ, (CH, CP))]. A phylogenic tree of [RH, JQ, (RP, GP), (CH, CP)] was determined using amino acid sequences. Because the local bootstrap values for the JQ branch in these trees are not high, additional sequence is necessary for construction of a reliable tree.     

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.     

旋唇纲纤毛虫系统发育分子树构建的影响因素

以36种旋唇类高等类群纤毛虫的核糖体小亚基核苷酸(Small subunit ribosomal RNA,SS rRNA)基因序列为素材,比较研究了不同条件(包括外类群、内类群的选择,同一基因不同序列长度的组合,不同建树方法和不同分析软件的使用)对纤毛虫分子系统树构建结果的影响。结果表明,上述因素均可不同程度地影响拓扑结构。结果同时提示,在利用有限数据进行相关研究,特别是在对未明类群的系统关系分析中,必须充分考虑因建树条件的不同所带来的影响。作者同时也建议,在当前可用的分子信息欠充分的前提下,对于纤毛虫任何类群的分子系统学探讨而言,慎重形成结论并尽可能地结合和参照形态学、发生学等资讯,仍是需优先考虑的工作路线。     

拟诺卡氏菌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基因更具优越性。     

分离自印度洋深海沉积物的部分链霉菌分离株多位点序列分析

【目的】采用多位点序列分析方法,研究印度洋3 000 m以下深海沉积物中分离得到的16S rRNA基因比对高度相似的链霉菌菌株的种间系统发育关系,同时探讨各管家基因及多基因聚类分析后的种间区分能力。【方法】以分离自印度洋深海沉积物的7株Streptomyces albidoflavus,11株Streptomyces cavourensis,16株Streptomyces pratensis为研究对象,以16S rRNA、atpD、recA和rpoB基因片段为标记,通过PCR扩增、测序,获得序列。同时从NCBI上下载5株S.pratensis上述4个基因的序列,将所有序列在MLST网站进行比对,并构建系统进化树进行比较。【结果】S.pratensis各菌株种内比较发现,16S rRNA基因构建的系统进化树中相同基因型的菌株没有聚在一起,系统进化树不稳定,区分度不高。其余3个构建的系统进化树稳定,菌株的聚类关系与MLST数据库得到的基因型一致。同时,多基因聚类分析后将菌株分为6个类群。在3个种的种间多位点序列比较中,除区分度明显增加、进化树更加稳定以外,还发现rec A基因进化上比较特殊的菌株。【结论】多位点序列分析将实验菌株分为很多不同的类型,成功地将所分离的链霉菌进行了更细的分类,同时也找到部分菌株在个别基因上差异较大。此方法可以用于相近种的快速鉴定。     

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.     

Molecular phylogenetic relationships of pond frogs distributed in the Palearctic region inferred from DNA sequences of mitochondrial 12S ribosomal RNA and cytochrome b genes

The evolutionary relationships of pond frogs distributed in the Far East and Europe were investigated by analyses of nucleotide sequences of mitochondrial 12S ribosomal RNA (12S rRNA) and cytochrome b (cyt b) genes. The nucleotide sequences of a 412-bp segment of the 12S rRNA gene and a 534-bp segment of the cyt b gene were determined by the PCR-direct sequencing method using 19 frogs belonging to six species and one subspecies distributed in the Palearctic region. Phylogenetic trees were constructed by the neighbor-joining and maximum-likelihood methods using Rana catesbeiana or Xenopus laevis as an outgroup. The 412-bp segment of the 12S rRNA gene contained 65 variable sites including gap sites, and the 534-bp segment of the cyt b gene contained 160 variable sites. The nucleotide sequence divergences of the 12S rRNA gene were 0.25-4.83% within the Far Eastern frogs, 0.25-6.22% within the European frogs, and 8.74-11.24% between the Far Eastern and the European frogs, whereas those of the cyt b gene were 3.64-14.73% within the Far Eastern frogs, 0.38-14.42% within the European frogs, and 16.53-23.58% between the Far Eastern and the European frogs. Although most nucleotide substitutions were at the third codon position of the cyt b gene and were silent mutations, 4 amino acid replacements occurred within the Far Eastern frogs, 4 within the European frogs, and 11 between the Far Eastern and the European frogs. The phylogenetic trees constructed from the nucleotide sequence divergences showed slightly different topologies for the 12S rRNA and cyt b genes. R. esculenta from Ukraine was closely related to R. lessonae from Luxembourg in both the 12S rRNA and the cyt b gene sequences.     

九种绢丝昆虫线粒体12S rRNA基因的序列特征和系统发育分析

为探讨鳞翅目中绢丝昆虫之间的系统发育关系和分子进化特征,本研究测定了中国柞蚕Antheraea pernyi野生型和放养型的线粒体12S rRNA基因的部分序列,结合来自GenBank数据库的17条序列,对总共9种绢丝昆虫（2科3属）的12S rRNA基因序列进行了分析。利用软件MEGA 3.1进行碱基组成、变异位点的统计和分子进化分析,分别用类平均聚类法（UPGMA）、邻接法（NJ）、最小进化法（ME）、最大简约法（MP）重建系统发生树。测定的中国柞蚕野生型的12S rRNA基因序列（427 bp）与放养型“豫早1号”的序列完全一致。序列对齐后共鉴定80个变异位点,50个简约信息位点。碱基组成分析显示在科属间具有明显差异,AT含量蚕蛾科高于大蚕蛾科;在A和T碱基的使用上,大蚕蛾科偏好使用T,而蚕蛾科则偏好使用A。与动物中常见的以转换为主的碱基替换模式不同,所分析的9种昆虫中除桑蚕属内部为转换与颠换基本一致外,其余物种间均是颠换多于转换。进化分析支持柞蚕属、樗蚕属和桑蚕属的单系。基于UPGMA法的进化树支持琥珀蚕是柞蚕属的较原始类型,而NJ、ME和MP法则支持印度柞蚕是较原始的类型,因此,柞蚕属种间的进化关系尚需进一步研究。     

Networks and groups within the genus Neisseria: analysis of argF, recA, rho, and 16S rRNA sequences from human Neisseria species.

To understand the pattern of nucleotide sequence variation among bacteria that frequently exchange chromosomal genes, we analyzed sequences of the recA, argF, and rho genes, as well as part of the small-subunit (16S) rRNA gene, from about 50 isolates of human commensal Neisseria species and the pathogenic N. meningitidis and N. gonorrhoeae. Almost all isolates of these species could be assigned to five phylogenetic groups that are found for all genes examined and generally are supported by high bootstrap values. In contrast, the phylogenetic relationships among groups varied according to the gene analyzed with notable incongruences involving N. cinerea and N. lactamica. Further analysis using split decomposition showed that for each gene, including 16S rRNA, the patterns of sequence divergence within N. meningitidis and closely related species were inconsistent with a bifurcating treelike phylogeny and better represented by an interconnected network. These data indicate that the human commensal Neisseria species can be separated into discrete groups of related species but that the relationships both within and among these groups, including those reconstructed using 16S rRNA, have been distorted by interspecies recombination events.     

干酪乳杆菌的近缘种及亚种部分看家基因的系统发育分析

【背景】16S rRNA基因序列分析已广泛应用于细菌的分类鉴定,但是存在一定局限性,而使用看家基因作为分子标记在近缘种及亚种间的系统发育分析中具有其独特的优势。【目的】研究16S rRNA、uvr C (核酸外切酶ABC,C亚基)和mur E (UDP-N-乙酰胞壁酰三肽合酶)基因序列对干酪乳杆菌的近缘种及亚种的区分能力。【方法】采用分离自传统发酵乳中的6株干酪乳杆菌为研究对象,选取uvr C和mur E基因片段,通过PCR扩增、测序,结合已公布的干酪乳杆菌的近缘种或亚种的相应序列计算遗传距离、构建系统发育树,并与16S rRNA基因序列分析技术进行比较。【结果】研究发现Lactobacilluscasei及相近种间的uvr C、mur E和联合基因(uvr C-mur E)构建的系统发育树拓扑结构与16S rRNA基因结果基本一致,区别在于相似性的不同,其分别为79.00%-99.16%、89.08%-99.20%、76.56%-99.69%和99.58%-100%。基于16S rRNA基因不能区分干酪乳杆菌的近缘种及亚种,而看家基因uvr C和mur E基因序列能够很好地区分干酪乳杆菌的近缘种及亚种,并且将uvr C和mur E基因串联使用后,试验菌株与参考菌株的分类关系更加清晰。【结论】联合基因(uvr C-mur E)可作为16SrRNA基因的辅助工具用于干酪乳杆菌的近缘种及亚种的快速准确鉴定。     

从线粒体16S rDNA序列探讨绒螯蟹类的系统发生关系

测定了绒螯蟹类各物种的线粒体16SrDNA部分片段的序列,构建了NJ树、ML树和MP树。序列歧异数据比较和各系统发生树都支持新绒螯蟹属(Neoeriocheir)为一个独立的属。在3种系统发生树中,直额绒螯蟹(Eriocheir recta)都是绒螯蟹属(Eriocheir)所有其它成员的姐妹群,并且广东珠江1只直额绒螯蟹标本的16SrDNA部分序列与台湾产台湾绒螯蟹(Eriocheir formasa)的相应序列相同。这些结果不支持平绒螯蟹属(Platyeriocheir)是一个有效的属,并表明E．formosa是E．recta的同物异名。绒螯蟹属(Eriocheir)所有其它成员聚为一个单系的分支,支持中华绒螯蟹、合浦绒螯蟹与日本绒螯蟹属于同一个物种Eriocheir japonica。16SrDNA部分序列的比对表明,产于台湾的日本绒螯蟹的此段序列与合浦绒螯蟹的相同,产于崇明岛的和产于美国旧金山海湾的中华绒螯蟹的此段序列与中华绒螯蟹单元型B的序列相同。     

Gene Sequence Phylogenies of the Family <Emphasis Type="Italic">Microbacteriaceae</Emphasis>

The type strains of 32 species of 13 genera of the family Microbacteriaceae were analysed with respect to gene-coding phylogeny for DNA gyrase subunit B (gyrB), RNA-polymerase subunit B (rpoB), recombinase A (recA), and polyphosphate kinase (ppk). The resulting gene trees were compared with the 16S rRNA gene phylogeny of the same strains. The topology of neighbour-joining and maximum parsimony phylogenetic trees, based on nucleic-acid sequences and protein sequences of housekeeping genes, differed from one another, and no gene tree was identical to that of the 16S rRNA gene tree. Most genera analysed containing >1 strain formed phylogenetically coherent taxa. The three pathovars of Curtobacterium flaccumfaciens clustered together to the exclusion of the type strains of other Curtobacterium species in all DNA - and protein-based analyses. In no tree did the distribution of a major taxonomic marker, i.e., diaminobutyric acid versus lysine and/or ornithine in the peptidoglycan, or acyl type of peptidoglycan, correlate with the phylogenetic position of the organisms. The changing phylogenetic position of Agrococcus jenensis was unexpected: This strain defined individual lineages in the trees based on 16S rRNA and gyrB and showed identity with Microbacterium saperdae in the other three gene trees.     

Estimation of bacterial species phylogeny through oligonucleotide frequency distances

Classification of bacteria is mainly based on sequence comparisons of certain homologous genes such as 16S rRNA. Recently there are challenges to classify bacteria using oligonucleotide frequency pattern of nonhomologous sequences. However, the evolutionary significance of oligonucleotides longer than tetra-nucleotide is not studied well. We performed phylogenetic analysis by using the Euclidean distances calculated from the di to deca-nucleotide frequencies in bacterial genomes, and compared these oligonucleotide frequency-based tree topologies with those for 16S rRNA gene and concatenated seven genes. When oligonucleotide frequency-based trees were constructed for bacterial species with similar GC content, their topologies at genus and family level were congruent with those based on homologous genes. Our results suggest that oligonucleotide frequency is useful not only for classification of bacteria, but also for estimation of their phylogenetic relationships for closely related species.