rRNA二级结构序列用于真菌系统学研究的方法初探

本文首次利用核酸二级结构特征代替核酸碱基作为探讨类群之间亲缘关系的信号,构建了基于结构特征的子囊菌部分类群的系统进化树。该方法以S(规范的碱基对),Q(不规范的碱基对),I(单链),B(侧环),M(多分枝环)和H(发卡结构)为代码将二级结构特征区分为6种不同的亚结构类型,然后将二级结构特征转换为结构序列,并进行结构序列分析。该方法使rRNA不只局限于碱基比较,拓展了其应用范围,为揭示分子的功能与进化的关系提供了线索。结果表明,结构序列分析可用于子囊菌的系统学研究;相对于核酸序列分析,结构分析的结果似乎更加清晰地体现子囊果的演化过程。     

Nonuniformity of nucleotide substitution rates in molecular evolution: Computer simulation and analysis of 5S ribosomal RNA sequences

Summary The effects of temporal (among different branches of a phylogeny) and spatial (among different nucleotide sites within a gene) nonuniformities of nucleotide substitution rates on the construction of phylogenetic trees from nucleotide sequences are addressed. Spatial nonuniformity may be estimated by using Shannon's (1948) entropy formula to measure the Relative Nucleotide Variability (RNV) at each nucleotide site in an aligned set of sequences; this is demonstrated by a comparative analysis of 5S rRNAs. New methods of constructing phylogenetic trees are proposed that augment the Unweighted Pair-Group Using Arithmetic Averages (UPGMA) algorithm by estimating and compensating for both spatial and temporal nonuniformity in substitution rates. These methods are evaluated by computer simulations of 5S rRNA evolution that include both kinds of nonuniformities. It was found that the proposed Reference Ratio Method improved both the ability to reconstruct the correct topology of a tree and also the estimation of branch lengths as compared to UPGMA. A previous method (Farris et al. 1970; Klotz et al. 1979; Li 1981) was found to be less successful in reconstructing topologies when there is high probability of multiple mutations at some sites. Phylogenetic analyses of 5S rRNA sequences support the endosymbiotic origins of both chloroplasts and mitochondria, even though the latter exhibit an accelerated rate of nucleotide substitution. Phylogenetic trees also reveal an adaptive radiation within the eubacteria and another within the eukaryotes for the origins of most major phyla within each group during the Precambrian era.     

Compositional Bias May Affect Both DNA-Based and Protein-Based Phylogenetic Reconstructions

It is now well-established that compositional bias in DNA sequences can adversely affect phylogenetic analysis based on those sequences. Phylogenetic analyses based on protein sequences are generally considered to be more reliable than those derived from the corresponding DNA sequences because it is believed that the use of encoded protein sequences circumvents the problems caused by nucleotide compositional biases in the DNA sequences. There exists, however, a correlation between AT/GC bias at the nucleotide level and content of AT- and GC-rich codons and their corresponding amino acids. Consequently, protein sequences can also be affected secondarily by nucleotide compositional bias. Here, we report that DNA bias not only may affect phylogenetic analysis based on DNA sequences, but also drives a protein bias which may affect analyses based on protein sequences. We present a striking example where common phylogenetic tools fail to recover the correct tree from complete animal mitochondrial protein-coding sequences. The data set is very extensive, containing several thousand sites per sequence, and the incorrect phylogenetic trees are statistically very well supported. Additionally, neither the use of the LogDet/paralinear transform nor removal of positions in the protein alignment with AT- or GC-rich codons allowed recovery of the correct tree. Two taxa with a large compositional bias continually group together in these analyses, despite a lack of close biological relatedness. We conclude that even protein-based phylogenetic trees may be misleading, and we advise caution in phylogenetic reconstruction using protein sequences, especially those that are compositionally biased. Received: 19 February 1998 / Accepted: 28 August 1998     

Phylogenetical relationship based on groE genes among phenotypically related Enterobacter,Pantoea, Klebsiella,Serratia and Erwinia species

In an attempt to define the phylogenetical relationship among 17 phenotypically related species of genera Enterobacter, Pantoea, Serratia, Klebsiella and Erwinia, we determined almost all of their groE operon sequences using the polymerase chain reaction direct sequencing method. The number of nucleotide substitutions per site was 0.12+/-0.030. The value was 3.6-fold higher than that of 16S rDNA. As a result, we were successful in constructing molecular phylogenetic trees which had a finer resolution than that based on the 16S rDNA sequences. The phylogenetic trees based on the nucleotide sequences and deduced amino acid sequences of groE operons indicated that the members of genera Enterobacter, Pantoea and Klebsiella were closely related to each other, while Serratia and Erwinia species except Erwinia carotovora, made distinct clades. The close relationship between Enterobacter aerogenes and Klebsiella pneumoniae, that had been suggested by biochemical tests and DNA hybridization, was also supported by our molecular phylogenetic trees.     

Conflict between Amino Acid and Nucleotide Characters

Slowly evolving characters, such as amino acids and replacement substitutions, have generally been favored over faster evolving characters for inferring phylogenetic relationships. However, amino acids constitute composite characters and, because of the degenerate genetic code, are subject to convergence. Based on an analysis of atpB and rbcL in 567 seed plants, we show that silent substitutions may be more phylogenetically informative than replacement substitutions and that artifacts caused by composite characters and/or convergence cause clades on amino acid trees to conflict with nucleotide trees and independent evidence. These findings indicate that coding nucleotide sequences only as amino acid characters for phylogenetic analysis provides little benefit and may yield misleading results.     

TrExML: a maximum-likelihood approach for extensive tree-space exploration

MOTIVATION: Maximum-likelihood analysis of nucleotide and amino acid sequences is a powerful approach for inferring phylogenetic relationships and for comparing evolutionary hypotheses. Because it is a computationally demanding and time-consuming process, most algorithms explore only a minute portion of tree-space, with the emphasis on finding the most likely tree while ignoring the less likely, but not significantly worse, trees. However, when such trees exist, it is equally important to identify them to give due consideration to the phylogenetic uncertainty. Consequently, it is necessary to change the focus of these algorithms such that near optimal trees are also identified. RESULTS: This paper presents the Advanced Stepwise Addition Algorithm for exploring tree-space and two algorithms for generating all binary trees on a set of sequences. The Advanced Stepwise Addition Algorithm has been implemented in TrExML, a phylogenetic program for maximum-likelihood analysis of nucleotide sequences. TrExML is shown to be more effective at finding near optimal trees than a similar program, fastDNAml, implying that TrExML offers a better approach to account for phylogenetic uncertainty than has previously been possible. A program, TreeGen, is also described; it generates binary trees on a set of sequences allowing for extensive exploration of tree-space using other programs. AVAILABILITY: TreeGen, TrExML, and the sequence data used to test the programs are available from the following two WWW sites: http://whitetail.bemidji.msus. edu/trexml/and http://jcsmr.anu.edu.au/dmm/humgen.+ ++html.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

Is 16S rDNA a Reliable Phylogenetic Marker to Characterize Relationships Below the Family Level in the Enterobacteriaceae?

The phylogenetic relationships of multiple enterobacterial species were reconstructed based on 16S rDNA gene sequences to evaluate the robustness of this housekeeping gene in the taxonomic placement of the enteric plant pathogens Erwinia, Brenneria, Pectobacterium, and Pantoea. Four data sets were compiled, two of which consisted of previously published data. The data sets were designed in order to evaluate how 16S rDNA gene phylogenies are affected by the use of different plant pathogen accessions and varying numbers of animal pathogen and outgroup sequences. DNA data matrices were analyzed using maximum likelihood (ML) algorithms, and character support was determined by ML bootstrap and Bayesian analyses. As additional animal pathogen sequences were added to the phylogenetic analyses, taxon placement changed. Further, the phylogenies varied in their placement of the plant pathogen species, and only the genus Pantoea was monophyletic in all four trees. Finally, bootstrap and Bayesian support values were low for most of the nodes, and all nonterminal branches collapsed in strict consensus trees. Inspection of 16S rDNA nucleotide alignments revealed several highly variable blocks punctuated by regions of conserved sequence. These data suggest that 16S rDNA, while effective for both species-level and family-level phylogenetic reconstruction, may underperform for genus-level phylogenetic analyses in the Enterobacteriaceae.     

Linguistic analysis in phylogeny estimation: a case study of mtDNAs of Bovidae

Molecular phylogenetic studies are executed by the alignment of protein or nucleotide sequences, followed by the construction of trees according either to distance, parsimony or maximum likelihood methods. Linguistic analysis was investigated here as an alternative method to aligning sequences. In an empirical study, we inferred trees for a variable number of Bovidae and sister taxa based on three different mitochondrial orthologous sequences. Comparison of our results with existing phylogenies indicated that the method, except for some still disputable points, was able to establish sensible systematic relationships, similar to patterns of radiation of the family found in recent studies.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

Molecular phylogenies of plants and Frankia support multiple origins of actinorhizal symbioses

Molecular phylogenetic trees were reconstructed from nucleotide sequences of nifH and 16S rDNA for Frankia and of rbcL for actinorhizal plants. Comparison of Frankia phylogenetic trees reconstructed using nifH and 16S rDNA sequences indicated that subgroupings of both trees correspond with each other in terms of plant origins of Frankia strains. The results suggested that 16S rDNAs can be utilized for coevolution analysis of actinorhizal symbioses. Frankia and plant phylogenetic trees reconstructed using 16S rDNA and rbcL sequences were compared. The comparison by tree matching and likelihood ratio tests indicated that although branching orders of both trees do not strictly correspond with each other, subgroupings of Frankia and their host plants correspond with each other in terms of symbiotic partnership. Estimated divergence times among Frankia and plant clades indicated that Frankia clades diverged more recently than plant clades. Taken together, actinorhizal symbioses originated more than three times after the four plant clades diverged.     

Phylogenetic analysis of gamma-proteobacteria inferred from nucleotide sequence comparisons of the house-keeping genes adk, aroE and gdh: comparisons with phylogeny inferred from 16S rRNA gene sequences

Nucleotide sequence comparisons of three house-keeping genes, adenylate kinase (adk), shikimate dehydrogenase (aroE), and glucose-6-phosphate dehydrogenase (gdh), were used to infer the phylogeny of 33 gamma-proteobacteria. Phylogenetic trees inferred from each gene, and from the concatenated sequences of all three genes, are, in general, similar to a 16S rRNA gene-inferred tree. Similar grouping of bacteria are revealed at the family, genus, species and strain levels in all five trees. The house-keeping genes, however, show a higher rate of nucleotide sequence substitutions. Consequently, they can possibly probe deeper branches of a phylogenetic tree than the 16S rRNA gene. However, because their nucleotide sequences are not as highly conserved among gamma-proteobacteria, family- or genus-specific primers would need to be designed for the amplification of any of these three house-keeping genes. Since these genes are used in multilocus sequence typing, it is expected that the number of sequences publicly available for many taxa will increase over time proving them very useful either at complementing 16S rRNA-inferred phylogenies or for specific, targeted, phylogenetic analysis.     

线蛱蝶亚科蝶类部分类群线粒体COI基因的分子系统发生分析

以线粒体细胞色素氧化酶I（COI）基因作分子标记,对线蛱蝶亚科蝴蝶进行序列测定．序列分析的结果表明．经比对和处理后的序列总长度是645bp,其中有199个变异位点,147个简约信息位点;所编码的氨基酸序列中有18个变异位点,7个信息位点．A＋T平均含量为69．6％,G＋C平均含量为30．4％,碱基组成出现AT偏斜．以蛱蝶亚科及秀蛱蝶亚科物种为外类群,用NJ、MP及贝叶斯法重建了该亚科的系统发生树,探讨了它们主要类群间的系统发生关系．分子系统树显示,线蛱蝶亚科由以下3大支系：环蛱蝶族＋翠蛱蝶族、线蛱蝶族、丽蛱蝶族构成;其中,环蛱蝶族为单系群（NJ树也支持线蛱蝶族的单系性）;翠蛱蝶族与环蛱蝶族亲缘关系较近：丽蛱蝶族可能是该亚科较早分化出的一支．     

Evaluation of atpB nucleotide sequences for phylogenetic studies of ferns and other pteridophytes

Inferring basal relationships among vascular plants poses a major challenge to plant systematists. The divergence events that describe these relationships occurred long ago and considerable homoplasy has since accrued for both molecular and morphological characters. A potential solution is to examine phylogenetic analyses from multiple data sets. Here I present a new source of phylogenetic data for ferns and other pteridophytes. I sequenced the chloroplast gene atpB from 23 pteridophyte taxa and used maximum parsimony to infer relationships. A 588-bp region of the gene appeared to contain a statistically significant amount of phylogenetic signal and the resulting trees were largely congruent with similar analyses of nucleotide sequences from rbcL. However, a combined analysis of atpB plus rbcL produced a better resolved tree than did either data set alone. In the shortest trees, leptosporangiate ferns formed a monophyletic group. Also, I detected a well-supported clade of Psilotaceae (Psilotum and Tmesipteris) plus Ophioglossaceae (Ophioglossum and Botrychium). The demonstrated utility of atpB suggests that sequences from this gene should play a role in phylogenetic analyses that incorporate data from chloroplast genes, nuclear genes, morphology, and fossil data.     

Multiple local maxima for likelihoods of phylogenetic trees: a simulation study.

Steel demonstrated that the maximum-likelihood function for a phylogenetic tree may have multiple local maxima. If this phenomenon were general, it would compromise the applicability of maximum likelihood as an optimality criterion for phylogenetic trees. In several simulation studies reported on in this paper, the true tree, and other trees of very high likelihood, rarely had multiple maxima. Our results thus provide reassurance that the value of maximum likelihood as a tree selection criterion is not compromised by the presence of multiple local maxima--the best estimates of the true tree are not likely to have them. This result holds true even when an incorrect nucleotide substitution model is used for tree selection.     

The phylogeny of Schistidium (Bryophyta, Grimmiaceae) based on primary and secondary structure study of nuclear rDNA internal transcribed spacers

The phylogeny of Schistidium (Bryophyta, Grimmiaceae) was studied on the basis of nucleotide sequences of internal transcribed spacers ITS1-2 of nuclear DNA and trnT-trnD region of chloroplast DNA. The consistency of phylogenetic trees constructed from nuclear and chloroplast sequences was shown. A basal grade and two large clades were resolved on the phylogenetic trees. Morphological characteristics specific for these clades were described. ITS1 and ITS2 secondary structures of Schistidium species were modeled using thermodynamic criteria. Four different structures of the longest ITS1 hairpin were identified. Possible paths of Schistidium evolution were considered based on the four types of ITS1 secondary structure and phylogenetic trees.     

Molecular evolution of serine/arginine splicing factors family (SR) by positive selection

The serine-rich (SR) protein family is involved in the pre-mRNA splicing process and the DNA sequences of the corresponding genes are highly conserved in the metazoan organisms. The mammalian SR proteins consist of one or two characteristic RNA binding domains (RBD), containing the signature sequences RDAEDA and SWQDLKD and a RS (arginine/serine-rich) domain. We used the amino acid and nucleotide sequences deposited in GenBank and Swiss-Prot databases to perform a phylogenetic analysis using bioinformatics tools. The results of the phylogenetic trees suggest that this family has evolved by several gene duplication events as a result of a positive selection mechanism.     

促葡萄糖转运蛋白基因家族的分子进化研究

葡萄糖转运蛋白是一个在结构上相似功能上不同的多基因家族（ＧＬＵＴ１－ＧＬＵＴ５）。由于这一组蛋白和体内的葡萄糖利用有关,因此被认为是糖尿病胰岛素抵抗（抗性）的一个候选基因。本文比较了不同种生物这一基因家族的氨基酸和核苷酸顺序;推测了亲水性和疏水性分布;计算了蛋白质和核苷酸的进化距离,并在此基础上构建了分子进化树。研究表明：这一基因家族具有高度的同源性、极为相似的亲水性和疏水性分布以及结构的对称性。提示这一基因家族起源于一个共同的祖先并可能通过基因的重复而形成。这一进化机制可能有利于氨基酸结构的稳定及抵抗突变的作用。由于邻元法构建的进化树其分支长度存在差异,提示在这一基因家族的进化过程中,各分支上的进化速率并不相同。蛋白质进化距离和核苷酸进化距离所构建进化树的差异提示了在基因组中可能存在隐匿替换。两种方法构建的进化树都提示了ＧＬＵＴ１、３、４在结构和功能上要更为保守。     

Maximum likelihood estimation of phylogenetic trees is consistent when substitution rates vary according to the invariable sites plus gamma distribution

Maximum likelihood estimation of phylogenetic trees from nucleotide sequences is completely consistent when nucleotide substitution is governed by the general time reversible (GTR) model with rates that vary over sites according to the invariable sites plus gamma (I + gamma) distribution.