The evolution of stramenopiles and alveolates as derived by “substitution rate calibration» of small ribosomal subunit RNA

The substitution rate of the individual positions in an alignment of 750 eukaryotic small ribosomal subunit RNA sequences was estimated. From the resulting rate distribution, an equation was derived that gives a more precise relationship between sequence dissimilarity and evolutionary distance than hitherto available. Trees constructed on the basis of evolutionary distances computed by this new equation for small ribosomal subunit RNA sequences from ciliates, apicomplexans, dinoflagellates, oomycetes, hyphochytriomycetes, bicosoecids, labyrinthuloids, and heterokont algae show a more consistent tree topology than trees constructed in the absence of substitution rate calibration. In particular, they do not suffer from anomalies caused by the presence of extremely long branches.     

滑膜肉瘤动物模型的分子细胞遗传学研究

用1例日本滑膜肉瘤（SS）患者的瘤组织标本接种裸鼠,获得了肿瘤生长。亲本肿瘤与裸鼠肿瘤在病理组织形态上存在一定差异,但两者的免疫组织化学特征相同。染色体分析表明,SS细胞株存在染色体数目和结构异常,患者外周血淋巴细胞染色体核型为正常女性,46,XX。 ＳＳ细胞株巴氏小体检出率低于对照,其正常X染色体比易位X染色体晚复制17小时。DNA印迹实验表明,SS DNA存在D2S3座位等位片段丢失,D1S57,D17S5和D13S30座位基因的部分缺失,但无DXS7,DXS14,和D2S44座位基因改变。     

Estimation of evolutionary distances between nucleotide sequences

A formal mathematical analysis of the substitution process in nucleotide sequence evolution was done in terms of the Markov process. By using matrix algebra theory, the theoretical foundation of Barry and Hartigan's (Stat. Sci. 2:191–210, 1987) and Lanave et al.'s (J. Mol. Evol. 20:86–93, 1984) methods was provided. Extensive computer simulation was used to compare the accuracy and effectiveness of various methods for estimating the evolutionary distance between two nucleotide sequences. It was shown that the multiparameter methods of Lanave et al.'s (J. Mol. Evol. 20:86–93, 1984), Gojobori et al.'s (J. Mol. Evol. 18:414–422, 1982), and Barry and Hartigan's (Stat. Sci. 2:191–210, 1987) are preferable to others for the purpose of phylogenetic analysis when the sequences are long. However, when sequences are short and the evolutionary distance is large, Tajima and Nei's (Mol. Biol. Evol. 1:269–285, 1984) method is superior to others.     

蛋白质序列间进化距离的估计

本文考察了目前采用的估计同源蛋白质序列间进化距离的方法缺陷,并提出了几个新的计算公式,它们考虑了氨基酸位点间显然存在的替代速率的差异。另外,提出了一种考虑氨基酸间不同替代概率的最大似然估计方法。文中对这些公式进行了计算比较,并对它在实际中的运用提出了建议。     

A Novel Approach to Phylogeny Reconstruction from Protein Sequences

The reliable reconstruction of tree topology from a set of homologous sequences is one of the main goals in the study of molecular evolution. If consistent estimators of distances from a multiple sequence alignment are known, the distance method is attractive because the tree reconstruction is consistent. To obtain a distance estimate d, the observed proportion of differences p (p-distance) is usually ``corrected' for multiple and back substitutions by means of a functional relationship d=f(p). In this paper the conditions under which this correction of p-distances will not alter the selection of the tree topology are specified. When these conditions are not fulfilled the selection of the tree topology may depend on the correction function applied. A novel method which includes estimates of distances not only between sequence pairs, but between triplets, quadruplets, etc., is proposed to strengthen the proper selection of correction function and tree topology. A ``super' tree that includes all tree topologies as special cases is introduced. Received: 17 February 1998 / Accepted: 20 July 1998     

A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences

Summary Some simple formulae were obtained which enable us to estimate evolutionary distances in terms of the number of nucleotide substitutions (and, also, the evolutionary rates when the divergence times are known). In comparing a pair of nucleotide sequences, we distinguish two types of differences; if homologous sites are occupied by different nucleotide bases but both are purines or both pyrimidines, the difference is called type I (or transition type), while, if one of the two is a purine and the other is a pyrimidine, the difference is called type II (or transversion type). Letting P and Q be respectively the fractions of nucleotide sites showing type I and type II differences between two sequences compared, then the evolutionary distance per site is K = — (1/2) ln {(1 — 2P — Q) }. The evolutionary rate per year is then given by k = K/(2T), where T is the time since the divergence of the two sequences. If only the third codon positions are compared, the synonymous component of the evolutionary base substitutions per site is estimated by K'_S = — (1/2) ln (1 — 2P — Q). Also, formulae for standard errors were obtained. Some examples were worked out using reported globin sequences to show that synonymous substitutions occur at much higher rates than amino acid-altering substitutions in evolution.Contribution No. 1330 from the National Institute of Genetics, Mishima, 411 Japan     

Simulation studies on the evolution of amino acid sequences

Summary A model of molecular evolution in which the parameter (intrinsic rate of amino acid substitution) fluctuates from time to time was investigated by simulating the process. It was found that the usual method of estimation such as Poisson fitting underestimates this variation of the parameter when remote comparisons are made. At the same time, four distance measures (minimum base difference, Poisson fitting, random nucleotide substitutions and negative binomial fitting) were tested for their accuracy. When the substitution rate is not uniform among the amino acid sites, the negative binomial fitting gives most satisfactory results, however, one needs to know the parameter beforehand in order to use this method. It was pointed out that the fluctuation of the evolutionary rate is expected if the nearly neutral but very slightly deleterious mutations play an important role on molecular evolution.Contribution No. 1087 from the National Institute of Genetics, Mishima, Shizuoka-ken, 411 Japan.     

脊髓灰质炎病毒的分子进化

本文根据脊髓灰质炎病毒3个型别参考毒株的基因组核苷酸序列和蛋白质氨基酸序列资料,首次试用Kimura的分子进化理论和计算方法,推算出脊髓灰质炎病毒型间的进化距离、分歧进化时间及病毒蛋白质氨基酸的替换率。结果表明:(1)型间毒株相互进化距离大致相等;(2)三个型病毒是由一个共同祖先病毒在距今约1—2千年以前几乎同时分歧进化而来;(3)型间毒株蛋白质氨基酸替换率也大致相等。     

Towards optimal distance functions for stochastic substitution models

Distance based reconstruction methods of phylogenetic trees consist of two independent parts: first, inter-species distances are inferred assuming some stochastic model of sequence evolution; then the inferred distances are used to construct a tree. In this paper we concentrate on the task of inter-species distance estimation. Specifically, we characterize the family of valid distance functions for the assumed substitution model and show that deliberate selection of distance function significantly improves the accuracy of distance estimates and, consequently, also improves the accuracy of the reconstructed tree.Our contribution consists of three parts: first, we present a general framework for constructing families of additive distance functions for stochastic evolutionary models. Then, we present a method for selecting (near) optimal distance functions, and we conclude by presenting simulation results which support our theoretical analysis.     

Tunicate tails, stolons, and the origin of the vertebrate trunk

Tunicates are primitive chordates that develop a transient 'tail' in the larval stage that is generally interpreted as a rudimentary version of the vertebrate trunk. Not all tunicates have tails, however. The groups that lack them, salps and pyrosomes, instead have a trunk-like reproductive stolon located approximately where the tail would otherwise be. In salps, files of blastozooids are formed along the sides of the stolon. The tail and caudal trunk in more advanced chordates could have evolved from a stolon of this type, an idea referred to here as the 'stolon hypothesis'. This means the vertebrate body could be a composite structure, since there is the potential for each somite to incorporate elements originally derived from a complete functional zooid. If indeed this has occurred, it should be reflected in some fashion in gene expression patterns in the vertebrate trunk. Selected morphological and molecular data are reviewed to show that they provide some circumstantial support for the stolon hypothesis. The case would be stronger if it could be demonstrated that salps and/or pyrosomes are ancestral to other tunicates. The molecular phylogenies so far available generally support the idea of a pelagic ancestor, but offer only limited guidance as to which of the surviving pelagic groups most closely resembles it. The principal testable prediction of the stolon hypothesis is that head structures (or their homologues) should be duplicated in series in the trunk in advanced chordates, and vice versa, i.e. trunk structures should occur in the head. The distribution of both rhabdomeric photoreceptors and nephridia in amphioxus conform with this prediction. Equally striking is the involvement of the Pax2 gene in the development of both the inner ear and nephric ducts in vertebrates. The stolon hypothesis would explain this as a consequence of the common origin of otic capsules and excretory ducts from atrial rudiments: from the paired rudiments of the parent oozooid in the case of the otic capsule (these express Pax2 according to recent ascidian data), and from tubular rudiments in the stolon in the case of the excretory ducts.     

羊轮状病毒NT株VP1基因的测序和分子进化分析

摘要：【目的】为了研究羊轮状病毒NT株VP1基因的遗传进化规律,【方法】根据GenBank中相关VP1基因的保守序列,设计合成引物,扩增NT株VP1基因并进行克隆测序和序列分析。【结果】 氨基酸序列比较表明NT株与其他毒株VP1基因的相似性为77.3%～98.4%,且氨基酸突变多发生在VP1蛋白的非功能区。VP1蛋白进化树表明NT株与牛轮状病毒处于同一进化分支,有较近的亲缘关系。结合26株具有代表性的轮状病毒,计算毒株间VP1基因的核苷酸和氨基酸进化距离,并对核苷酸的同义突变率（dS）和非同义突变率（dN）进行研究,发现dN/dS的比值小于1,说明同义替代是VP1基因在进化过程中的主要变异。【结论】本文首次对羊轮状病毒NT株进行了VP1基因的测序,并对VP1基因的进化距离和进化规律进行深入探讨。     

Early Evolution of Protists Introductory Remarks of Chairman

SYNOPSIS. Phylogenetic relations are ultimately determined by homologous relationships, including both structural and functional data. Following the establishment of those relationships, the direction of evolutionary change must be determined using paleontologic, developmental, and especially morphocline, data. From that perspective the direction of subsequent development becomes clearer and the problems of origins become more explicit. Using the foregoing methodology it has been possible to identify plausibly monophyletic groups of animal protists or protozoa. Allowing for attendant difficulties, there nevertheless emerges certain fairly convincing generalities: (a) the predominantly pseudopodial forms, with a few minor exceptions, have direct origins from apochlorotic algae; (b) the predominantly kinetidal forms (zooflagellates and ciliates), though also derived originally from apochlorotic algae, give evidence of extended evolutionary development with the especially noteworthy emergence of a permanent ingestatory structure; (c) both groups have increased size, a tendency towards multinuclearity and polyploidy, cytoplasmic differentiations of various sorts, and complex life cycles. In terms of further evolution, namely the emergence of multicellular animals, the pseudopodial forms are almost certainly a dead end while the kinetidal forms are arguably the ancestors of at least 2 metazoan phyla.     

Evolutionary Relationships Among the Eukaryotic Crown Taxa Taking into Account Site-to-Site Rate Variation in 18S rRNA

In this study we constructed a bootstrapped distance tree of 500 small subunit ribosomal RNA sequences from organisms belonging to the so-called crown of eukaryote evolution. Taking into account the substitution rate of the individual nucleotides of the rRNA sequence alignment, our results suggest that (1) animals, true fungi, and choanoflagellates share a common origin: The branch joining these taxa is highly supported by bootstrap analysis (bootstrap support [BS] > 90%), (2) stramenopiles and alveolates are sister groups (BS = 75%), (3) within the alveolates, dinoflagellates and apicomplexans share a common ancestor BS > 95%), while in turn they both share a common origin with the ciliates (BS > 80%), and (4) within the stramenopiles, heterokont algae, hyphochytriomycetes, and oomycetes form a monophyletic grouping well supported by bootstrap analysis (BS > 85%), preceded by the well-supported successive divergence of labyrinthulomycetes and bicosoecids. On the other hand, many evolutionary relationships between crown taxa are still obscure on the basis of 18S rRNA. The branching order between the animal-fungal-choanoflagellates clade and the chlorobionts, the alveolates and stramenopiles, red algae, and several smaller groups of organisms remains largely unresolved. When among-site rate variation is not considered, the inferred tree topologies are inferior to those where the substitution rate spectrum for the 18S rRNA is taken into account. This is primarily indicated by the erroneous branching of fast-evolving sequences. Moreover, when different substitution rates among sites are not considered, the animals no longer appear as a monophyletic grouping in most distance trees. Received: 11 June 1997 / Accepted: 21 July 1997     

Synonymous nucleotide substitution rates of β-tubulin and histone genes conform to high overall genomic rates in rodents but not in sea urchins

Summary Sea urchin and rodent genomes have been posited to evolve rapidly as indicated by divergences in single copy nuclear DNA sequences. We have examined whether the synonymous substitution rates of three highly conserved genes, -tubulin, histone H4, and histone H3, adhere to these high genomic substitution rates by comparing sequences between two sea urchins,Strongylocentrotus purpuratus andLytechinus pictus, and between rodents and humans. Whereas the rate of change between the 3 untranslated regions of the -tubulin cDNA ofS. purpuratus (Sp-1), sequenced in this study, and ofL. pictus (Lp-3) was consistent with the overall rate of change estimated from previous DNA hybridization results between these species, the synonymous substitution rates for the carboxyl domains of these -tubulins, as well as for the late histones H4 and H3, were significantly depressed. In contrast, synonymous nucleotide substitution rates between rodents and between rodent and human for the carboxyl domain proper of identical -tubulin isotypes and for histone H4 and H3.1 did not differ from the overall rate of change for the rodent genomes. Moreover, an analysis of paralogous human and mouse -tubulin sequences supported the conclusion that the synonymous substitution rates in the mouse were higher than those in the human. Differences in constraint on evolutionary change were not evident strictly from the conserved amino acid sequences and base compositions of these genes. Other constraining influences seemed more relevant to the departure of the synonymous substitution rates of the sea urchin -tubulin and histone coding regions from the average genomic rate.     

Correlation Between Ka/Ks and Ks is Related to Substitution Model and Evolutionary Lineage

In 2005, Wyckoff and coworkers described a surprisingly strong correlation between Ka/Ks and Ks in several data sets using the LPB93 algorithm. This finding indicated the possibility of a paradigm shift in the way selection strength can be measured using the Ka/Ks ratio. We carried out a calculation of Ka and Ks using six different algorithms on three cross-species orthologous data sets and found a highly variable correlation among the algorithms and lineages. Algorithms based on the GY-HKY substitution model exhibit a weaker positive correlation or a stronger negative correlation than those based on the K2P and JC69 substitution model. Even if one algorithm shows a positive correlation between Ka/Ks and Ks in a warm-blooded lineage, it may show no correlation in a cold-blooded lineage. This algorithm-related and evolutionary lineage-related correlation indicates the need for great caution in drawing conclusions when using only one Ka and Ks algorithm in a genomewide analysis of selection strength. Our results indicated that currently used algorithms for Ka and Ks calculations are flawed and need improvements.     

Phylogenetic Relationships of Maples (Acer L.; Aceraceae) Implied by Nuclear Ribosomal ITS Sequences

Acer. ITS 1 sequences in twenty-eight species of Acer and a species of Dipteronia in the family Aceraceae ranged from 220 to 242 bp and ITS 2 sequences from 215 to 251 bp. The size of the 5.8S coding region was 164 bp for all species examined in the family. Phylogenetic analysis of ITS sequences placed a very robust clade of section Palmata at the base of the tree. Three species of section Parviflora sensu de Jong (1994), A. spicatum, A. distylum and A. nipponicum, did not form a monophyletic clade. Acer spicatum was separated from the robust clade of A. distylum and A. nipponicum. Molecular tree strongly supports the close relationship among section Platanoidea, Glabra series Arguta, and section Macrantha. The close relationship between sections Pentaphylla and Trifoliata was also strongly suggested in ITS tree. Sections Rubra and Hyptiocarpa appeared to be closely allied with each other. The average rate of nucleotide substitution was estimated as (8.0±1.9)×10⁻¹¹ substitutions per site per year for ITS 1 and (9.0±1.6) ×10⁻¹¹ for ITS 2. Received 17 December 1999/ Accepted in revised form 10 March 2000     

Phylogenetic Inference with Weighted Codon Evolutionary Distances

We develop a new approach to estimate a matrix of pairwise evolutionary distances from a codon-based alignment based on a codon evolutionary model. The method first computes a standard distance matrix for each of the three codon positions. Then these three distance matrices are weighted according to an estimate of the global evolutionary rate of each codon position and averaged into a unique distance matrix. Using a large set of both real and simulated codon-based alignments of nucleotide sequences, we show that this approach leads to distance matrices that have a significantly better treelikeness compared to those obtained by standard nucleotide evolutionary distances. We also propose an alternative weighting to eliminate the part of the noise often associated with some codon positions, particularly the third position, which is known to induce a fast evolutionary rate. Simulation results show that fast distance-based tree reconstruction algorithms on distance matrices based on this codon position weighting can lead to phylogenetic trees that are at least as accurate as, if not better, than those inferred by maximum likelihood. Finally, a well-known multigene dataset composed of eight yeast species and 106 codon-based alignments is reanalyzed and shows that our codon evolutionary distances allow building a phylogenetic tree which is similar to those obtained by non-distance-based methods (e.g., maximum parsimony and maximum likelihood) and also significantly improved compared to standard nucleotide evolutionary distance estimates.     

Molecular evolutionary rates of oncogenes

Summary Using nine sets of viral and cellular oncogenes, the rates of nucleotide substitutions were computed by using Gojobori and Yokoyama's (1985) method. The results obtained confirmed our previous conclusion that the rates of nucleotide substitution for the viral oncogenes are about a million times higher than those for their cellular counterparts. For cellular oncogenes and most viral oncogenes, however, the rate of synonymous substitution is higher than that of nonsynonymous substitution. Moreover, the pattern of nucleotide substitutions for viral oncogenes is more similar to that for functional genes (such as cellular oncongenes) than for pseudogenes. This implies that nucleotide substitutions in viral oncogenes may be functionally constrained. Thus, our observation supports that nucleotide substitutions for the oncogenes in those DNA and RNA genomes are consistent with Kimura's neutral theory of molecular evolution (Kimura 1968, 1983).     

Codon usage divergence of homologous vertebrate genes and codon usage clock

Summary This paper is concerned with the divergence of synonymous codon usage and its bias in three homologous genes within vertebrate species. Genetic distances among species are described in terms of synonymous codon usage divergence and the correlation is found between the genetic distances and taxonomic distances among species under study. A codon usage clock is reported in alphaglobin and beta-globin. A method is developed to define the synonymous codon preference bias and it is observed that the bias changes considerably among species.     

A new molecular evolution model for limited insertion independent of substitution

We recently introduced a new molecular evolution model called the IDIS model for Insertion Deletion Independent of Substitution  and . In the IDIS model, the three independent processes of substitution, insertion and deletion of residues have constant rates. In order to control the genome expansion during evolution, we generalize here the IDIS   model by introducing an insertion rate which decreases when the sequence grows and tends to 0 for a maximum sequence length _nmax$n_{\max }$.