Gene's Functional Arrangement as a Measure of thePhylogenetic Relationships of Microorganisms

With the development of genome sequencing more whole genomes of microorganisms were completed, many methods wereintroduced to reconstruct the phylogenetic tree of those microorganismswith the information extracted from the whole genomes through variousways of transforming or mapping the whole genome sequences into otherforms which can describe the evolutionary distance in a new way. We thinkit might be possible that there exists information buried in the wholegenome transferred along lineage, which remains stable and is moreessential than sequence conservation of individual genes or the arrangementof some genes of a selected set. We need to find one measurement that caninvolve as many phylogenetic features as possible that are beyond thegenome sequence itself. We converted each genome sequence of themicroorganisms into another linear sequence to represent the functionalstructure of the sequence, and we used a new information function tocalculate the discrepancy of sequences and to get one distance matrix of thegenomes, and built one phylogenetic tree with a neighbor joining method.The resulting tree shows that the major lineages are consistent with theresult based on their 16srRNA sequences. Our method discovered onephylogenetic feature derived from the genome sequences and the encodedgenes that can rebuild the phylogenetic tree correctly. The mapping of onegenome sequence to its new form representing the relative positions of thefunctional genes provides a new way to measure the phylogeneticrelationships, and with the more specific classification of gene functions theresult could be more sensitive.     

Phylogenetic tree derived from bacterial, cytosol and organelle 5S rRNA sequences.

A phylogenetic tree was constructed by computer analysis of 47 completely determined 5S rRNA sequences. The wheat mitochondrial sequence is significantly more related to prokaryotic than to eukaryotic sequences, and its affinity to that of the thermophilic Gram-negative bacterium Thermus aquaticus is comparable to the affinity between Anacystis nidulans and chloroplastic sequences. This strongly supports the idea of an endosymbiotic origin of plant mitochondria. A comparison of the plant cytosol and chloroplast sub-trees suggests a similar rate of nucleotide substitution in nuclear genes and chloroplastic genes. Other features of the tree are a common precursor of protozoa and metazoa, which appears to be more related to the fungal than to the plant protosequence, and an early divergence of the archebacterial sequence (Halobacterium cutirubrum) from the prokaryotic branch.     

Phylogenetic Relationships Among Mycoplasmas Based on the Whole Genomic Information

With the rapid increase in available bacterial whole-genome information, comparison of bacteria at the whole-genome level has proven to be highly useful in microbial phylogenetic research. Here we constructed a phylogenetic tree based on 15 whole genomes of Mycoplasma and the related bacteria. First, 143 orthologous gene families that are shared by all of the 15 bacteria were selected and 143 multiple alignments were generated. Next, a concatenated multiple alignment inferred from the 143 multiple alignments was generated. A total of 43,370 amino acid sites were considered in the neighbor-joining analysis. The phylogenetic tree based on the whole-genomic information indicated that the 15 bacteria were divided into four major groups with 100% bootstrap support, i.e., the M. hyopneumoniae (Mhy) group, the M. mycoides (Mmy) group, the M. pneumoniae (Mpn) group, and the Bacillus-Phytoplasma (BP) group. In the phylogenetic tree, the Mhy group was more closely related to the Mpn group than the Mmy group. The relationships among the Mhy, Mmy, Mpn, and BP groups were supported with 100% in bootstrap analysis. The phylogenetic tree based on the whole-genome comparison is different from the 16S rRNA tree. Thirty-nine of the 143 phylogenetic trees had the same type of the topology based on the whole-genome comparison. However, we could not identify a gene family contributing or solely belonging to the topology of the 39 proteins. In this study, we showed that some proteins, such as RpoA, RpoB, RpoC, and RpoD, are not suitable for evolutionary studies on relationships among major groups of mycoplasmas. We also showed that glycolysis-related genes of Ureaplasma have a higher substitution rate than the other bacteria. The phylogenetic approaches at the whole-genome level are very important and will be essential for microbial evolutionary studies. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. John Oakeshott     

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.     

Phylogenetic position of Salinibacter ruber based on concatenated protein alignments

A total of 22 genes from the genome of Salinibacter ruber strain M31 were selected in order to study the phylogenetic position of this species based on protein alignments. The selection of the genes was based on their essential function for the organism, dispersion within the genome, and sufficient informative length of the final alignment. For each gene, an individual phylogenetic analysis was performed and compared with the resulting tree based on the concatenation of the 22 genes, which rendered a single alignment of 10,757 homologous positions. In addition to the manually chosen genes, an automatically selected data set of 74 orthologous genes was used to reconstruct a tree based on 17,149 homologous positions. Although single genes supported different topologies, the tree topology of both concatenated data sets was shown to be identical to that previously observed based on small subunit (SSU) rRNA gene analysis, in which S. ruber was placed together with Bacteroidetes. In both concatenated data sets the bootstrap was very high, but an analysis with a gradually lower number of genes indicated that the bootstrap was greatly reduced with less than 12 genes. The results indicate that tree reconstructions based on concatenating large numbers of protein coding genes seem to produce tree topologies with similar resolution to that of the single 16S rRNA gene trees. For classification purposes, 16S rRNA gene analysis may remain as the most pragmatic approach to infer genealogic relationships.     

枣食芽象甲线粒体基因组全序列测定与系统发育分析

【目的】从线粒体基因组水平上探讨枣食芽象甲Scythropus yasumatsui与近缘种的系统发育关系。【方法】利用Illumina MiSeq测序平台对枣食芽象甲线粒体基因组进行测序,对基因组序列进行拼装、注释和特征分析;利用贝叶斯法和最大似然法构建基于象甲科13个物种的线粒体基因组13个蛋白质编码基因核苷酸序列的系统发育树。【结果】结果表明,枣食芽象甲线粒体基因组全长为16 472 bp (GenBank登录号: MF807224),包含13个蛋白质编码基因、22个tRNA基因、2个rRNA基因和2个非编码控制区,37个基因的排列顺序与祖先昆虫的线粒体基因排列顺序一致。13个蛋白质编码基因的起始密码子为ATN,其中除了cob和nad1基因的完全终止密码子为TAG外,其余11个基因的完全终止密码子为TA(A)。22个tRNA基因中除了trnS1缺少DHU臂,反密码子由GCT变为TCT外,其余均能形成典型的三叶草结构。基于13个蛋白质编码基因序列构建的系统发育树结果显示,象甲科8个亚科系统发育关系为：(((隐喙象亚科(Cryptorhynchinae)+(象虫亚科(Curculioninae)+魔喙象亚科(Molytinae)))+长小蠹亚科(Platypodinae))+(粗喙象亚科(Entiminae)+Cyclominae亚科))+隐颏象亚科(Dryophthorinae)+小蠹亚科(Scolytinae))。【结论】在13种象甲科昆虫物种中,同属于粗喙象亚科的枣食芽象甲与南美果树象甲Naupactus xanthographus在系统发育树中聚为同一分支,表明基于线粒体基因组全序列的分子系统发育结果与传统的形态分类结果是一致的。     

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

【背景】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基因的辅助工具用于干酪乳杆菌的近缘种及亚种的快速准确鉴定。     

根瘤菌新类群代表菌株的16SrDNA全序列测定及其系统发育地位

用双脱氧法测定了一个根瘤菌新类群代表菌株SH2672的16S rDNA全序列,将此全序列与根瘤菌各已知种及相关种的16S rDNA全序列进行了比较及聚类分析,得到系统发育树状图。在系统发育树状图中,菌株SH2672与百脉根中慢生根瘤菌(Mesorhizobium loti),华癸中慢生根瘤菌(M. huakuii)、天山中慢生根瘤菌(M. tianshanense)、地中海中慢生根瘤菌(M. mediterraneum)、鹰嘴豆中慢生根瘤菌(M. ciceri)共同构成一个分支,与各已知种的模式菌株16S rDNA相似性分别为:96.3％,96.4％,97.2％,95.1％,95.6％,均在95％以上,它们应归属于同一属。且分支内各种间DNA同源性低于70％,表明它们分别为不同的种,菌株SH2672代表着一个新的根瘤菌种。     

Toward a more robust assessment of intraspecies diversity, using fewer genetic markers

Phylogenetic sequence analysis of single or multiple genes has dominated the study and census of the genetic diversity among closely related bacteria. It remains unclear, however, how the results based on a few genes in the genome correlate with whole-genome-based relatedness and what genes (if any) best reflect whole-genome-level relatedness and hence should be preferentially used to economize on cost and to improve accuracy. We show here that phylogenies of closely related organisms based on the average nucleotide identity (ANI) of their shared genes correspond accurately to phylogenies based on state-of-the-art analysis of their whole-genome sequences. We use ANI to evaluate the phylogenetic robustness of every gene in the genome and show that almost all core genes, regardless of their functions and positions in the genome, offer robust phylogenetic reconstruction among strains that show 80 to 95% ANI (16S rRNA identity, >98.5%). Lack of elapsed time and, to a lesser extent, horizontal transfer and recombination make the selection of genes more critical for applications that target the intraspecies level, i.e., strains that show >95% ANI according to current standards. A much more accurate phylogeny for the Escherichia coli group was obtained based on just three best-performing genes according to our analysis compared to the concatenated alignment of eight genes that are commonly employed for phylogenetic purposes in this group. Our results are reproducible within the Salmonella, Burkholderia, and Shewanella groups and therefore are expected to have general applicability for microevolution studies, including metagenomic surveys.     

蒙古狼线粒体基因组序列分析及其系统发育地位

应用Long-PCR和克隆测序法得到蒙古狼(Canis lupus chanco)线粒体基因组全序列,结合GenBank中现有犬科动物线粒体基因组数据,应用最大简约法(MP)、最大似然法(ML)和Bayesian分析法对蒙古狼的系统发育地位进行了探讨。结果如下:蒙古狼线粒体基因组全长16709bp,包含13个蛋白质编码基因、22个tRNA基因、2个rRNA基因和1个非编码区。序列碱基的组成存在明显的A-T偏好性。tRNA基因中除tRNA-Ser(AGY)缺少双氢尿嘧啶(DHU)臂以外,其余均能折叠成典型的三叶草二级结构。大多数蛋白质编码基因的起始和终止密码子与犬科动物有报道相同,COXⅡ基因的起始密码子为ATA,与其他犬科动物不同。基于12S rRNA+16S rRNA+H链上的12个蛋白质编码基因的联合数据的系统发育分析发现,在已报道的狼亚种数据中,西藏狼(Canis lupus laniger)的分化时间最早,其次为阿拉伯狼(Canis lupus arabs),蒙古狼与欧亚狼(Canis lupuslupus)的系统发育地位最为接近。     

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

【目的】采用多位点序列分析方法,研究印度洋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基因进化上比较特殊的菌株。【结论】多位点序列分析将实验菌株分为很多不同的类型,成功地将所分离的链霉菌进行了更细的分类,同时也找到部分菌株在个别基因上差异较大。此方法可以用于相近种的快速鉴定。     

The threonine-tRNA ligase gene region is applicable in classification,typing, and phylogenetic analysis of bifidobacteria

In the modern era, molecular genetic techniques are crucial in ecological studies, as well as in the classification, typing, and phylogenetic analysis of prokaryotes. These techniques are mainly aimed at whole genome comparisons and PCR-derived experiments, including amplifying the 16S rRNA and other various housekeeping genes used in taxonomy, as well as MLST (multilocus sequence typing) and MLSA (multilocus sequence analysis) of different taxonomic bacterial groups. The gene encoding threonine-tRNA ligase (thrS) is a gene potentially applicable as an identification and phylogenetic marker in bacteria. It is widely distributed in bacterial genomes and is subject to evolutionary selection pressure due to its important function in protein synthesis. In this study, specific primers were used to amplify a thrS gene fragment (~740 bp) in 36 type and 30 wild strains classified under family Bifidobacteriaceae. The full-length gene has not yet been considered as a possible identification, classification, and phylogenetic marker in bifidobacteria. The thrS sequences revealed higher sequence variability (82.7% of pairwise identities) among members of the family than that shown by 16S rRNA gene sequences (96.0%). Although discrepancies were found between the thrS-derived and previously reported whole genome phylogenetic analyses, the main phylogenetic groups of bifidobacteria were properly assigned. Most wild strains of bifidobacteria were better differentiated based on their thrS sequences than on their 16S rRNA gene identities. Phylogenetic confidence of the evaluated gene with respect to other alternative genetic markers widely used in taxonomy of bifidobacteria (fusA, GroELhsp60, pyrG, and rplB genes) was confirmed using the localized incongruence difference - Templeton analysis.     

Host-symbiont stability and fast evolutionary rates in an ant-bacterium association: cospeciation of camponotus species and their endosymbionts, candidatus blochmannia

Bacterial endosymbionts are widespread across several insect orders and are involved in interactions ranging from obligate mutualism to reproductive parasitism. Candidatus Blochmannia gen. nov. (Blochmannia) is an obligate bacterial associate of Camponotus and related ant genera (Hymenoptera: Formicidae). The occurrence of Blochmannia in all Camponotus species sampled from field populations and its maternal transmission to host offspring suggest that this bacterium is engaged in a long-term, stable association with its ant hosts. However, evidence for cospeciation in this system is equivocal because previous phylogenetic studies were based on limited gene sampling, lacked statistical analysis of congruence, and have even suggested host switching. We compared phylogenies of host genes (the nuclear EF-1alphaF2 and mitochondrial COI/II) and Blochmannia genes (16S ribosomal DNA [rDNA], groEL, gidA, and rpsB), totaling more than 7 kilobases for each of 16 Camponotus species. Each data set was analyzed using maximum likelihood and Bayesian phylogenetic reconstruction methods. We found minimal conflict among host and symbiont phylogenies, and the few areas of discordance occurred at deep nodes that were poorly supported by individual data sets. Concatenated protein-coding genes produced a very well-resolved tree that, based on the Shimodaira-Hasegawa test, did not conflict with any host or symbiont data set. Correlated rates of synonymous substitution (d(S)) along corresponding branches of host and symbiont phylogenies further supported the hypothesis of cospeciation. These findings indicate that Blochmannia-Camponotus symbiosis has been evolutionarily stable throughout tens of millions of years. Based on inferred divergence times among the ant hosts, we estimated rates of sequence evolution of Blochmannia to be approximately 0.0024 substitutions per site per million years (s/s/MY) for the 16S rDNA gene and approximately 0.1094 s/s/MY at synonymous positions of the genes sampled. These rates are several-fold higher than those for related bacteria Buchnera aphidicola and Escherichia coli. Phylogenetic congruence among Blochmannia genes indicates genome stability that typifies primary endosymbionts of insects.     

Ribosomal protein gene cluster analysis in eubacterium genomics: homology between Sinorhizobium meliloti strain 1021 and Bacillus subtilis

The first whole genome sequence of a symbiotic soil bacterium, Sinorhizobium meliloti (formely named Rhizobium meliloti) strain 1021, is due in 2001. As an active participant in the European and North American consortium that has completed this work, our group has sequenced a region on the chromosome containing clusters rpoBC, str, S10, spc and alpha corresponding to 30 protein genes. The structural organization and function of these genes were compared with those of orthologs in another 15 complete eubacterial genomes available in databases. This study, involving the DNA and amino acid sequences as well as the organization of the whole region (gene order, cluster order, etc.), has shown that the phylogenetic tree resulting from a comparison of the amino acid sequence is rather similar to that derived from 16S rRNA sequence data. However, the tree achieved by aligning DNA sequences groups the organisms with a high GC content (>60% GC), while that based on a comparison of gene cluster orientation and organization reveals a greater level of correspondence between the α-proteobacteria S.meliloti and the firmicute Bacillus subtilis.     

Phylogenetic Analysis Reveals Multiple Lateral Transfers of Adenosine-5′-Phosphosulfate Reductase Genes among Sulfate-Reducing Microorganisms

Lateral gene transfer affects the evolutionary path of key genes involved in ancient metabolic traits, such as sulfate respiration, even more than previously expected. In this study, the phylogeny of the adenosine-5'-phosphosulfate (APS) reductase was analyzed. APS reductase is a key enzyme in sulfate respiration present in all sulfate-respiring prokaryotes. A newly developed PCR assay was used to amplify and sequence a fragment ( approximately 900 bp) of the APS reductase gene, apsA, from a taxonomically wide range of sulfate-reducing prokaryotes (n = 60). Comparative phylogenetic analysis of all obtained and available ApsA sequences indicated a high degree of sequence conservation in the region analyzed. However, a comparison of ApsA- and 16S rRNA-based phylogenetic trees revealed topological incongruences affecting seven members of the Syntrophobacteraceae and three members of the Nitrospinaceae, which were clearly monophyletic with gram-positive sulfate-reducing bacteria (SRB). In addition, Thermodesulfovibrio islandicus and Thermodesulfobacterium thermophilum, Thermodesulfobacterium commune, and Thermodesulfobacterium hveragerdense clearly branched off between the radiation of the delta-proteobacterial gram-negative SRB and the gram-positive SRB and not close to the root of the tree as expected from 16S rRNA phylogeny. The most parsimonious explanation for these discrepancies in tree topologies is lateral transfer of apsA genes across bacterial divisions. Similar patterns of insertions and deletions in ApsA sequences of donor and recipient lineages provide additional evidence for lateral gene transfer. From a subset of reference strains (n = 25), a fragment of the dissimilatory sulfite reductase genes (dsrAB), which have recently been proposed to have undergone multiple lateral gene transfers (M. Klein et al., J. Bacteriol. 183:6028-6035, 2001), was also amplified and sequenced. Phylogenetic comparison of DsrAB- and ApsA-based trees suggests a frequent involvement of gram-positive and thermophilic SRB in lateral gene transfer events among SRB.     

一株类芽孢杆菌的分离鉴定及其抗革兰氏阴性菌活性测定

【背景】随着碳青霉烯类和多粘菌素类可转移耐药基因的发现及扩散,多重耐药革兰氏阴性细菌感染更加难以治疗。【目的】筛选有效拮抗革兰氏阴性菌的菌株,为新型抗生素的发掘奠定基础。【方法】利用胰蛋白胨大豆琼脂培养基筛选土壤源细菌,通过16S rRNA基因序列鉴定其种属;通过全基因组测序,antiSMASH比对分析菌株产抗生素潜能,双层琼脂平板法验证其抗菌活性;通过甲醇萃取其次级代谢产物,高效液相色谱串联质谱(HPLC-MS/MS)进行次级代谢产物分析。【结果】从北京周边土壤样品中分离到一株类芽孢杆菌CAU136 (Paenibacillus pabuli CAU136),经过生物信息学分析和antiSMASH比对,表明该菌株有较强合成次级代谢产物的潜能,双层琼脂平板法验证其能抑制多株革兰氏阴性菌生长,HPLC-MS/MS检测结果显示其可能分泌多粘菌素E。【结论】类芽孢杆菌CAU136可能分泌多粘菌素E,能有效拮抗革兰氏阴性菌。     

Chemotaxis away from thiocyanic and isothiocyanic esters in Pseudomonas aeruginosa

Abstract A novel mycoplasmal species designated as Mycoplasma penetrans has recently been isolated from patients infected with human immunodeficiency virus. The 16S rRNA gene from this mycoplasma was cloned and its nucleotide sequence determined. This sequence was aligned with previously published homologous sequences from several mycoplasmas and with related Gram-positive bacteria and a phylogenetic tree was constructed. The results indicate that M. penetrans belongs to the evolutionary group Pneumoniae.     

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.