Substitution rate calibration of small subunit ribosomal RNA identifies chlorarachniophyte endosymbionts as remnants of green algae.

Chlorarachniophytes are amoeboid algae with chlorophyll a and b containing plastids that are surrounded by four membranes instead of two as in plants and green algae. These extra membranes form important support for the hypothesis that chlorarachniophytes have acquired their plastids by the ingestion of another eukaryotic plastid-containing alga. Chlorarachniophytes also contain a small nucleus-like structure called the nucleomorph situated between the two inner and the two outer membranes surrounding the plastid. This nucleomorph is a remnant of the endosymbiont's nucleus and encodes, among other molecules, small subunit ribosomal RNA. Previous phylogenetic analyses on the basis of this molecule provided unexpected and contradictory evidence for the origin of the chlorarachniophyte endosymbiont. We developed a new method for measuring the substitution rates of the individual nucleotides of small subunit ribosomal RNA. From the resulting substitution rate distribution, we derived an equation that gives a more realistic relationship between sequence dissimilarity and evolutionary distance than equations previously available. Phylogenetic trees constructed on the basis of evolutionary distances computed by this new method clearly situate the chlorarachniophyte nucleomorphs among the green algae. Moreover, this relationship is confirmed by transversion analysis of the Chlorarachnion plastid small subunit ribosomal RNA.

Pea formaldehyde-active class III alcohol dehydrogenase: common derivation of the plant and animal forms but not of the corresponding ethanol-active forms (classes I and P).

A plant class III alcohol dehydrogenase (or glutathione-dependent formaldehyde dehydrogenase) has been characterized. The enzyme is a typical class III member with enzymatic parameters and substrate specificity closely related to those of already established animal forms. Km values with the pea enzyme are 6.5 microM for NAD+, 2 microM for S-hydroxymethylglutathione, and 840 microM for octanol versus 9, 4, and 1200 microM, respectively, with the human enzyme. Structurally, the pea/human class III enzymes are closely related, exhibiting a residue identity of 69% and with only 3 of 23 residues differing among those often considered in substrate and coenzyme binding. In contrast, the corresponding ethanol-active enzymes, the long-known human liver and pea alcohol dehydrogenases, differ more (47% residue identities) and are also in functionally important active site segments, with 12 of the 23 positions exchanged, including no less than 7 at the usually much conserved coenzyme-binding segment. These differences affect functionally important residues that are often class-distinguishing, such as those at positions 48, 51, and 115, where the plant ethanol-active forms resemble class III (Thr, Tyr, and Arg, respectively) rather than the animal ethanol-active class I forms (typically Ser, His, and Asp, respectively). Calculations of phylogenetic trees support the conclusions from functional residues in subgrouping plant ethanol-active dehydrogenases and the animal ethanol-active enzymes (class I) as separate descendants from the class III line. It appears that the classical plant alcohol dehydrogenases (now called class P) have a duplicatory origin separate from that of the animal class I enzymes and therefore a paralogous relationship with functional convergence of their alcohol substrate specificity. Combined, the results establish the conserved nature of class III also in plants, and contribute to the molecular and functional understanding of alcohol dehydrogenases by defining two branches of plant enzymes into the system.

Equine rhinovirus 1 is more closely related to foot-and-mouth disease virus than to other picornaviruses.

Equine rhinovirus 1 (ERhV1) is a respiratory pathogen of horses which has an uncertain taxonomic status. We have determined the nucleotide sequence of the ERhV1 genome except for a small region at the 5' end. The predicted polyprotein was encoded by 6741 nucleotides and possessed a typical picornavirus proteolytic cleavage pattern, including a leader polypeptide. The genomic structure and predicted amino acid sequence of ERhV1 were more similar to those of foot-and-mouth disease viruses (FMDVs), the only members of the aphthovirus genus, than to those of other picornaviruses. Features which were most similar to FMDV included a 16-amino acid 2A protein which was 87.5% identical in sequence of FMDV 2A, a leader (L) protein similar in size to FMDV Lab and the possibility of a truncated L protein similar in size to FMDV Lb, and a 3C protease which recognizes different cleavage sites. However, unlike FMDV, ERhV1 had only one copy of the 3B (VPg) polypeptide. The phylogenetic relationships of the ERhV1 sequence and nucleotide sequences of representative species of the five genera of the family Picornaviridae were examined. Nucleotide sequences coding for the complete polyprotein, the RNA polymerase, and VP1 were analyzed separately. The phylogenetic trees confirmed that ERhV1 was more closely related to FMDV than to other picornaviruses and suggested that ERhV1 may be a member, albeit very distant, of the aphthovirus genus.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

The vertebrate alcohol dehydrogenase system: variable class II type form elucidates separate stages of enzymogenesis.

A mixed-class alcohol dehydrogenase has been characterized from avian liver. Its functional properties resemble the classical class I type enzyme in livers of humans and animals by exhibiting low Km and kcat values with alcohols (Km = 0.7 mM with ethanol) and low Ki values with 4-methylpyrazole (4 microM). These values are markedly different from corresponding parameters of class II and III enzymes. In contrast, the primary structure of this avian liver alcohol dehydrogenase reveals an overall relationship closer to class II and to some extent class III (69 and 65% residue identities, respectively) than to class I or the other classes of the human alcohol dehydrogenases (52-61%), the presence of an insertion (four positions in a segment close to position 120) as in class II but in no other class of the human enzymes, and the presence of several active site residues considered typical of the class II enzyme. Hence, the avian enzyme has mixed-class properties, being functionally similar to class I, yet structurally similar to class II, with which it also clusters in phylogenetic trees of characterized vertebrate alcohol dehydrogenases. Comparisons reveal that the class II enzyme is approximately 25% more variable than the "variable" class I enzyme, which itself is more variable than the "constant" class III enzyme. The overall extreme, and the unusual chromatographic behavior may explain why the class II enzyme has previously not been found outside mammals. The properties define a consistent pattern with apparently repeated generation of novel enzyme activities after separate gene duplications.

Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

Molecular cloning and sequence analysis of expansins--a highly conserved, multigene family of proteins that mediate cell wall extension in plants.

Expansins are unusual proteins discovered by virtue of their ability to mediate cell wall extension in plants. We identified cDNA clones for two cucumber expansins on the basis of peptide sequences of proteins purified from cucumber hypocotyls. The expansin cDNAs encode related proteins with signal peptides predicted to direct protein secretion to the cell wall. Northern blot analysis showed moderate transcript abundance in the growing region of the hypocotyl and no detectable transcripts in the nongrowing region. Rice and Arabidopsis expansin cDNAs were identified from collections of anonymous cDNAs (expressed sequence tags). Sequence comparisons indicate at least four distinct expansin cDNAs in rice and at least six in Arabidopsis. Expansins are highly conserved in size and sequence (60-87% amino acid sequence identity and 75-95% similarity between any pairwise comparison), and phylogenetic trees indicate that this multigene family formed before the evolutionary divergence of monocotyledons and dicotyledons. Sequence and motif analyses show no similarities to known functional domains that might account for expansin action on wall extension. A series of highly conserved tryptophans may function in expansin binding to cellulose or other glycans. The high conservation of this multigene family indicates that the mechanism by which expansins promote wall extensin tolerates little variation in protein structure.

Evidence from molecular systematics for decreased avian diversification in the pleistocene Epoch.

Pleistocene glaciations have been suggested as major events influencing speciation rates in vertebrates. Avian paleontological studies suggest that most extant species evolved in the Pleistocene Epoch and that species' durations decreased through the Pleistocene because of heightened speciation rates. Molecular systematic studies provide another data base for testing these predictions. In particular, rates of diversification can be determined from molecular phylogenetic trees. For example, an increasing rate of speciation (but constant extinction) requires shorter intervals between successive speciation events on a phylogenetic tree. Examination of the cumulative distribution of reconstructed speciation events in mtDNA phylogenies of 11 avian genera, however, reveals longer intervals between successive speciation events as the present time is approached, suggesting a decrease in net diversification rate through the Pleistocene Epoch. Thus, molecular systematic studies do not indicate a pulse of Pleistocene diversification in passerine birds but suggest, instead, that diversification rates were lower in the Pleistocene than for the preceding period. Documented habitat shifts likely led to the decreased rate of diversification, although from molecular evidence we cannot discern whether speciation rates decreased or extinction rates increased.

The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.

Phylogenetic analyses are increasingly used in attempts to clarify transmission patterns of human immunodeficiency virus type 1 (HIV-1), but there is a continuing discussion about their validity because convergent evolution and transmission of minor HIV variants may obscure epidemiological patterns. Here we have studied a unique HIV-1 transmission cluster consisting of nine infected individuals, for whom the time and direction of each virus transmission was exactly known. Most of the transmissions occurred between 1981 and 1983, and a total of 13 blood samples were obtained approximately 2-12 years later. The p17 gag and env V3 regions of the HIV-1 genome were directly sequenced from uncultured lymphocytes. A true phylogenetic tree was constructed based on the knowledge about when the transmissions had occurred and when the samples were obtained. This complex, known HIV-1 transmission history was compared with reconstructed molecular trees, which were calculated from the DNA sequences by several commonly used phylogenetic inference methods [Fitch-Margoliash, neighbor-joining, minimum-evolution, maximum-likelihood, maximum-parsimony, unweighted pair group method using arithmetic averages (UPGMA), and a Fitch-Margoliash method assuming a molecular clock (KITSCH)]. A majority of the reconstructed trees were good estimates of the true phylogeny; 12 of 13 taxa were correctly positioned in the most accurate trees. The choice of gene fragment was found to be more important than the choice of phylogenetic method and substitution model. However, methods that are sensitive to unequal rates of change performed more poorly (such as UPGMA and KITSCH, which assume a constant molecular clock). The rapidly evolving V3 fragment gave better reconstructions than p17, but a combined data set of both p17 and V3 performed best. The accuracy of the phylogenetic methods justifies their use in HIV-1 research and argues against convergent evolution and selective transmission of certain virus variants.

Calculating the probability of multitaxon evolutionary trees: bootstrappers Gambit.

The reconstruction of multitaxon trees from molecular sequences is confounded by the variety of algorithms and criteria used to evaluate trees, making it difficult to compare the results of different analyses. A global method of multitaxon phylogenetic reconstruction described here, Bootstrappers Gambit, can be used with any four-taxon algorithm, including distance, maximum likelihood, and parsimony methods. It incorporates a Bayesian-Jeffreys'-bootstrap analysis to provide a uniform probability-based criterion for comparing the results from diverse algorithms. To examine the usefulness of the method, the origin of the eukaryotes has been investigated by the analysis of ribosomal small subunit RNA sequences. Three common algorithms (paralinear distances, Jukes-Cantor distances, and Kimura distances) support the eocyte topology, whereas one (maximum parsimony) supports the archaebacterial topology, suggesting that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.