SEARCH · PubMed Health

Results for “bootstrap”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41Linked to original sources

PubMed2005

The comparison of the confidence regions in phylogeny.

In this paper, several different procedures for constructing confidence regions for the true evolutionary tree are evaluated both in terms of coverage and size without considering model misspecification. The regions are constructed on the basis of tests of hypothesis using six existing tests: Shimodaira Hasegawa (SH), SOWH, star form of SOWH (SSOWH), approximately unbiased (AU), likelihood weight (LW), generalized least squares, plus two new tests proposed in this paper: single distribution nonparametric bootstrap (SDNB) and single distribution parametric bootstrap (SDPB). The procedures are evaluated on simulated trees both with small and large number of taxa. Overall, the SH, SSOWH, AU, and LW tests led to regions with higher coverage than the nominal level at the price of including large numbers of trees. Under the specified model, the SOWH test gives accurate coverage and relatively small regions. The SDNB and SDPB tests led to the small regions with occasional undercoverage. These two procedures have a substantial computational advantage over the SOWH test. Finally, the cutoff levels for the SDNB test are shown to be more variable than those for the SDPB test.

Classification↗

PubMed2006

An improved statistical method for detecting heterotachy in nucleotide sequences.

The principle of heterotachy states that the substitution rate of sites in a gene can change through time. In this article, we propose a powerful statistical test to detect sites that evolve according to the process of heterotachy. We apply this test to an alignment of 1289 eukaryotic rRNA molecules to 1) determine how widespread the phenomenon of heterotachy is in ribosomal RNA, 2) to test whether these heterotachous sites are nonrandomly distributed, that is, linked to secondary structure features of ribosomal RNA, and 3) to determine the impact of heterotachous sites on the bootstrap support of monophyletic groupings. Our study revealed that with 21 monophyletic taxa, approximately two-thirds of the sites in the considered set of sequences is heterotachous. Although the detected heterotachous sites do not appear bound to specific structural features of the small subunit rRNA, their presence is shown to have a large beneficial influence on the bootstrap support of monophyletic groups. Using extensive testing, we show that this may not be due to heterotachy itself but merely due to the increased substitution rate at the detected heterotachous sites.

Chi-Square Distribution↗

PubMed1993

MUST, a computer package of Management Utilities for Sequences and Trees.

The MUST package is a phylogenetically oriented set of programs for data management and display, allowing one to handle both raw data (sequences) and results (trees, number of steps, bootstrap proportions). It is complementary to the main available software for phylogenetic analysis (PHYLIP, PAUP, HENNING86, CLUSTAL) with which it is fully compatible. The first part of MUST consists of the acquisition of new sequences, their storage, modification, and checking of sequence integrity in files of aligned sequences. In order to improve alignment, an editor function for aligned sequences offers numerous options, such as selection of subsets of sequences, display of consensus sequences, and search for similarities over small sequence fragments. For phylogenetic reconstruction, the choice of species and portions of sequences to be analyzed is easy and very rapid, permitting fast testing of numerous combinations of sequences and taxa. The resulting files can be formatted for most programs of tree construction. An interactive tree-display program recovers the output of all these programs. Finally, various modules allow an in-depth analysis of results, such as comparison of distance matrices, variation of bootstrap proportions with respect to various parameters or comparison of the number of steps per position. All presently available complete sequences of 28S rRNA are furnished aligned in the package. MUST therefore allows the management of all the operations required for phylogenetic reconstructions.

Amino Acid Sequence↗

PubMed2001

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction.

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.

Animals↗

PubMed2001

Statistics of divergence times.

Given the number of nucleotide substitutions between two species (K) and the substitution rate nu, the expectation of the corresponding divergence time is usually calculated as K/(2 nu). This is strictly true only if nu is regarded as a constant because the ratio of two random variables, such as K/(2 nu), has distributional properties different from those of the distribution of K. Therefore, both the mean and any confidence interval for divergence times are unknown in this situation. We model the distribution of K and nu using the Gamma distribution and calculate the mean and 95% confidence interval for the corresponding divergence time. These calculations are compared with results obtained by bootstrapping sequence data from the model plant Arabidopsis thaliana and its relatives. We show that for nonoverlapping pairs of phylogenetic distances, our method approaches the bootstrap results very closely. In contrast, regarding the mutation rate as a constant leads to strong underestimation of the confidence interval. An implementation of our method of computing divergence times is accessible through a web interface at http://www.soft.ice.mpg.de/cite.

Arabidopsis↗

PubMed2001

Fungal Zuotin proteins evolved from MIDA1-like factors by lineage-specific loss of MYB domains.

Proteins are often characterized by the presence of multiple domains, which make specific contributions to their cellular function. While the gain of domains in proteins by duplication and shuffling is well established, domain loss is poorly documented. Here, we provide evidence that domain loss has played an important role in the evolution of protein architecture and function by demonstrating that fungal Zuotin proteins evolved from MIDA1-like proteins, present in animals and plants, by complete loss of the carboxyl-terminal MYB domains. Phylogenetic analyses of the DnaJ motif (the J domain) present in both Zuotin and MIDA1 proteins were complicated by the limited length and profound differences in evolutionary rates exhibited by this domain. To rigorously examine J domain phylogeny, we combined the nonparametric bootstrap with Monte Carlo simulation. This method, which we have designated the resampled parametric bootstrap, allowed us to assess type I and type II error associated with these analyses. These results revealed significant support for domain loss rather than domain gain or gene loss involving paralogs. The absence of sequences related to the MIDA1 MYB domains in Saccharomyces cerevisiae further indicates that the domains have been completely lost, consistent with known functional differences between Zuotin and MIDA1 proteins. These analyses suggest that the description of additional examples of complete domain loss may provide a method to identify orthologous proteins exhibiting functional differences using genomic sequence data.

Amino Acid Sequence↗

PubMed1996

Molecular divergence and phylogeny: rates and patterns of cytochrome b evolution in cranes.

Analyses of complete cytochrome b sequences from all species of cranes (Aves: Gruidae) reveal aspects of sequence evolution in the early stages of divergence. These DNA sequences are > or = 89% identical, but expected departures from random substitution are evident. Silent, third-position pyrimidine transitions are the dominant substitution type, with transversion comprising only a small fraction of sequence differences. Substitution patterns are not clearly manifested until divergence has reached a moderate level (> 3%), as expected for a stochastic process. Variation in the frequency of mismatch types among lineages decreases at larger divergences, but the level of bias does not decay. Divergence varies up to fivefold among gene regions but is not correlated with structural domain. All protein structural domains except extramembrane 4 display < 20% variable residues. Regions corresponding to putative functional domains show the excepted conservation of amino acids, although the C-terminal portion of the Q0 reaction center displays several nonconservative replacements. Phylogenetic analyses incorporating substitution asymmetries produced mixed results. Distances estimated with multiple parameters (transition, codon-position, composition, and pyrimidine-transition biases) yielded identical additive tree topologies with comparable bootstrap values, all consistent with uncontroversial species relationships. Maximum likelihood analysis incorporating these biases, as well as equally weighted parsimony analysis, produced similar results. Static, differential weighting for parsimony did not improve the phylogenetic signal but produced unusual trees with low bootstraps. The overall rate of nucleotide substitution varies slightly but significantly among cranes, and calibration of distances against fossil dates suggests divergence rates of 0.7%-1.7% per million years.

Animals↗

PubMed1996

Noncoding sequences from the slowly evolving chloroplast inverted repeat in addition to rbcL data do not support gnetalean affinities of angiosperms.

We developed PCR primers against highly conserved regions of the rRNA operon located within the inverted repeat of the chloroplast genome and used these to amplify the region spanning from the 3' terminus of the 23S rRNA gene to the 5' terminus of the 5S rRNA gene. The sequence of this roughly 500-bp region, which includes the 4.5S rRNA gene and two chloroplast intergenic transcribed spacer regions (cpITS2 and cpITS3), was determined from 20 angiosperms, 7 gymnosperms, and 16 ferns (21,700 bp). Sequences for the large subunit of ribulose bisphosphate carboxylase/oxygenase (rbcL) from the same or confamilial genera were analyzed in both separate and combined data sets. Due to the low substitution rate in the inverted repeat region, noncoding sequences in the cpITS region are not saturated with substitutions, in contrast to synonymous sites in rbcL, which are shown to evolve roughly six times faster than noncoding cpITS sequences. Several length polymorphisms with very clear phylogenetic distributions were detected in the data set. Results of phylogenetic analyses provide very strong bootstrap support for monophyly of both spermatophytes and angiosperms. No support for a sister group relationship between Gnetales and angiosperms in either cpITS or rbcL data was found. Rather, weak bootstrap support for monophyly of gymnosperms studied and for a basal position for the aquatic angiosperm Nymphaea among angiosperms studied was observed. Noncoding sequences from the inverted repeat region of chloroplast DNA appear suitable for study of land plant evolution.

Base Sequence↗

PubMed1996

Loculoascomycete origins and evolution of filamentous ascomycete morphology based on 18S rRNA gene sequence data.

The fungal subclass Loculoascomycetes is characterized by asci or sexual spore sacks with two separable wall layers. At maturity, the inner wall layer protrudes out beyond the outer wall as in a jack-in-the-box. If the Loculoascomycetes were monophyletic and their jack-in-the-box type asci evolved once, then taxa from diverse loculoascomycete lineages would cluster together in a DNA sequence-based tree. To evaluate the phylogenetic history of the two-walled asci, I sequenced the 18S nuclear rRNA genes of 16 species from seven families in the loculoascomycete orders Pleosporales, Dothideales, and Chaetothyriales. Within the Loculoascomycetes, the Pleosporales form a monophyletic group in 99% of the bootstrapped parsimony trees. The Dothideales usually appear as a monophyletic group but without statistical support. Capronia pilosella (Herpotrichiellaceae, Order Chaetothyriales) clusters with plectomycete members of the subclass Euascomycetes rather than the other Loculoascomycetes in 99% of parsimony and neighbor-joining bootstrap replicates. Although the jack-in-the-box-type ascus is a good marker for large, monophyletic loculoascomycete orders, it must have evolved at least twice or been lost at least once.

Ascomycota↗

PubMed1998

28S and 18S rDNA sequences support the monophyly of lampreys and hagfishes.

Resolving the interrelationships of three major extant lineages of vertebrates (hagfishes, lampreys, and gnathostomes) is a particularly important issue in evolution, because the basal resolution critically influences our understanding of primitive vertebrate characters. A consensus has emerged over the last 20 years that lampreys are the sister group to the gnathostomes and the hagfishes represent an ancient, basal lineage. This hypothesis has essentially displaced the classical hypothesis of monophyly of the cyclostomes (lampreys plus hagfishes). To test these hypotheses, we compared nearly complete ribosomal DNA sequences from each of these major lineages, as well as those from a cephalochordate and a urochordate, which represent a paraphyletic outgroup for assessing the basal vertebrate relationships. For this comparison, 92%-99% complete 28S rDNA sequences were obtained from the lancelet Branchiostoma floridae, the hagfish Eptatretus stouti, the lamprey Petromyzon marinus, and cartilaginous fishes Hydrolagus colliei and Squalus acanthias and were then analyzed with previously reported 28S and 18S rDNA sequences from other chordates. We conducted conventional (nonparametric) bootstrap analyses, under maximum-likelihood, parsimony, and minimum-evolution (using LogDet distances) criteria, of both 28S and 18S rDNA sequences considered separately and combined. All these analyses provide moderate to very strong support for the monophyly of the cyclostomes. Furthermore, the currently accepted hypothesis of a lamprey-gnathostome clade is moderately rejected by the Kishino-Hasegawa test (P = 0.099) and resoundingly rejected by parametric bootstrap tests (P < 0.01) in favor of monophyly of living cyclostomes. Another significant finding is that the hagfish E. stouti has the longest 28S rDNA gene known in any organism (> 5,200 nt).

Animals↗

PubMed1998

Complete mitochondrial DNA sequence of the fat dormouse, Glis glis: further evidence of rodent paraphyly.

The complete mitochondrial genome of the fat dormouse, Glis glis, has been sequenced (16,602 bp). A total of 23 complete mitochondrial mammalian genomes have been taken into account for phylogenetic reconstruction. Phylogenetic analyses were performed with parsimony, distance (stationary Markov model), and maximum-likelihood methods. In all cases, data strongly support the paraphyly of rodents, with dormouse and guinea pig in a different clade from rat and mouse, reaching bootstrap values of 95%. Rodent monophyly and the existence of Glires (Rodentia and Lagomorpha) are weakly supported, with maximum bootstrap values of 11% and 8.6%, respectively. This result agrees with the analyses of isochore patterns in the nuclear genome and the B2 and B2-like retroposons, which show a close relationship between dormice and guinea pigs rather than between dormice and rats and mice.

Animals↗

PubMed1998

Searching for the closest living relative(s) of tetrapods through evolutionary analyses of mitochondrial and nuclear data.

The phylogenetic relationships of the African lungfish (Protopterus dolloi) and the coelacanth (Latimeria chalumnae) with respect to tetrapods were analyzed using complete mitochondrial genome DNA sequences. A lungfish + coelancanth clade was favored by maximum parsimony (although this result is dependent on which transition:transversion weights are applied), and a lungfish + tetrapod clade was supported by neighbor-joining and maximum-likelihood analyses. These two hypotheses received the strongest statistical and bootstrap support to the exclusion of the third alternative, the coelacanth + tetrapod sister group relationship. All mitochondrial protein coding genes combined favor a lungfish + tetrapod grouping. We can confidently reject the hypothesis that the coelacanth is the closest living relative of tetrapods. When the complete mitochondrial sequence data were combined with nuclear 28S rRNA gene data, a lungfish + coelacanth clade was supported by maximum parsimony and maximum likelihood, but a lungfish + tetrapod clade was favored by neighbor-joining. The seeming conflicting results based on different data sets and phylogenetic methods were typically not statistically strongly supported based on Kishino-Hasegawa and Templeton tests, although they were often supported by strong bootstrap values. Differences in rate of evolution of the different mitochondrial genes (slowly evolving genes such as the cytochrome oxidase and tRNA genes favored a lungfish + coelacanth clade, whereas genes of relatively faster substitution rate, such as several NADH dehydrogenase genes, supported a lungfish + tetrapod grouping), as well as the rapid radiation of the lineages back in the Devonian, rather than base compositional biases among taxa seem to be directly responsible for the remaining uncertainty in accepting one of the two alternate hypotheses.

Animals↗

PubMed1999

Locus specificity of polymorphic alleles and evolution by a birth-and-death process in mammalian MHC genes.

We have conducted an extensive phylogenetic analysis of polymorphic alleles from human and mouse major histocompatibility complex (MHC) class I and class II genes. The phylogenetic tree obtained for 212 complete human class I allele sequences (HLA-A, -B, and -C) has shown that all alleles from the same locus form a single cluster, which is highly supported by bootstrap values, except for one HLA-B allele (HLA-B*7301). Mouse MHC class I loci did not show locus-specific clusters of polymorphic alleles. This was considered to be because of either interlocus genetic exchange or the confusing designation of loci in different haplotypes at the present time. The locus specificity of polymorphic alleles was also observed in human and mouse MHC class II loci. It was therefore concluded that interlocus recombination or gene conversion is not very important for generating MHC diversity, with a possible exception of mouse class I loci. According to the phylogenetic trees of complete coding sequences, we classified human MHC class I (HLA-A, -B, and -C) and class II (DRB1) alleles into three to five major allelic lineages (groups), which were monophyletic with high bootstrap values. Most of these allelic groups remained unchanged even in phylogenetic trees based on individual exons, though this does not exclude the possibility of intralocus recombination involving short DNA segments. These results, together with the previous observation that MHC loci are subject to frequent duplication and deletion, as well as to balancing selection, indicate that MHC evolution in mammals is in agreement with the birth-and-death model of evolution, rather than with the model of concerted evolution.

Alleles↗

PubMed1999

An evaluation of elongation factor 1 alpha as a phylogenetic marker for eukaryotes.

Elongation factor 1 alpha (EF-1 alpha) is a highly conserved ubiquitous protein involved in translation that has been suggested to have desirable properties for phylogenetic inference. To examine the utility of EF-1 alpha as a phylogenetic marker for eukaryotes, we studied three properties of EF-1 alpha trees: congruency with other phyogenetic markers, the impact of species sampling, and the degree of substitutional saturation occurring between taxa. Our analyses indicate that the EF-1 alpha tree is congruent with some other molecular phylogenies in identifying both the deepest branches and some recent relationships in the eukaryotic line of descent. However, the topology of the intermediate portion of the EF-1 alpha tree, occupied by most of the protist lineages, differs for different phylogenetic methods, and bootstrap values for branches are low. Most problematic in this region is the failure of all phylogenetic methods to resolve the monophyly of two higher-order protistan taxa, the Ciliophora and the Alveolata. JACKMONO analyses indicated that the impact of species sampling on bootstrap support for most internal nodes of the eukaryotic EF-1 alpha tree is extreme. Furthermore, a comparison of observed versus inferred numbers of substitutions indicates that multiple overlapping substitutions have occurred, especially on the branch separating the Eukaryota from the Archaebacteria, suggesting that the rooting of the eukaryotic tree on the diplomonad lineage should be treated with caution. Overall, these results suggest that the phylogenies obtained from EF-1 alpha are congruent with other molecular phylogenies in recovering the monophyly of groups such as the Metazoa, Fungi, Magnoliophyta, and Euglenozoa. However, the interrelationships between these and other protist lineages are not well resolved. This lack of resolution may result from the combined effects of poor taxonomic sampling, relatively few informative positions, large numbers of overlapping substitutions that obscure phylogenetic signal, and lineage-specific rate increases in the EF-1 alpha data set. It is also consistent with the nearly simultaneous diversification of major eukaryotic lineages implied by the "big-bang" hypothesis of eukaryote evolution.

Algorithms↗

PubMed2000

Bias in phylogenetic reconstruction of vertebrate rhodopsin sequences.

Two spurious nodes were found in phylogenetic analyses of vertebrate rhodopsin sequences in comparison with well-established vertebrate relationships. These spurious reconstructions were well supported in bootstrap analyses and occurred independently of the method of phylogenetic analysis used (parsimony, distance, or likelihood). Use of this data set of vertebrate rhodopsin sequences allowed us to exploit established vertebrate relationships, as well as the considerable amount known about the molecular evolution of this gene, in order to identify important factors contributing to the spurious reconstructions. Simulation studies using parametric bootstrapping indicate that it is unlikely that the spurious nodes in the parsimony analyses are due to long branches or other topological effects. Rather, they appear to be due to base compositional bias at third positions, codon bias, and convergent evolution at nucleotide positions encoding the hydrophobic residues isoleucine, leucine, and valine. LogDet distance methods, as well as maximum-likelihood methods which allow for nonstationary changes in base composition, reduce but do not entirely eliminate support for the spurious resolutions. Inclusion of five additional rhodopsin sequences in the phylogenetic analyses largely corrected one of the spurious reconstructions while leaving the other unaffected. The additional sequences not only were more proximal to the corrected node, but were also found to have intermediate levels of base composition and codon bias as compared with neighboring sequences on the tree. This study shows that the spurious reconstructions can be corrected either by excluding third positions, as well as those encoding the amino acids Ile, Val, and Leu (which may not be ideal, as these sites can contain useful phylogenetic signal for other parts of the tree), or by the addition of sequences that reduce problems associated with convergent evolution.

Animals↗

PubMed2000

Actin gene family evolution and the phylogeny of coleoid cephalopods (Mollusca: Cephalopoda).

Phylogenetic analysis conducted on a 784-bp fragment of 82 actin gene sequences of 44 coleoid cephalopod taxa, along with results obtained from genomic Southern blot analysis, confirmed the presence of at least three distinct actin loci in coleoids. Actin isoforms were characteri zed through phylogenetic analysis of representative cephalopod sequences from each of the three isoforms, along with translated actin cDNA sequences from a diverse array of metazoan taxa downloaded from GenBank. One of the three isoforms found in cephalopods was closely related to actin sequences expressed in the muscular tissues of other molluscs. A second isoform was most similar to cytoplasmic-specific actin amino acid sequences. The muscle type actins of molluscs were found to be distinct from those of arthropods, suggesting at least two independent derivations of muscle actins in the protostome lineage, although statistical support for this conclusion was lacking. Parsimony and maximum-likelihood analyses of two of the isoforms from which >30 orthologous coleoid sequences had been obtained (one of the cytoplasmic actins and the muscle actin) supported the monophyly of several higher-level coleoid taxa. These included the superorders Octopodiformes and Decapodiformes, the order Octopoda, the octopod suborder Incirrata, and the teuthoid suborder Myopsida. The monophyly of several taxonomic groups within the Decapodiformes was not supported, including the orders Teuthoidea and Sepioidea and the teuthoid suborder Oegopsida. Parametric bootstrap analysis conducted on the simulated cytoplasmic actin data set provided statistical support to reject the monophyly of the Sepioidea. Although parametric bootstrap analysis of the muscle actin isoform did not reject sepioid monophyly at the 5% level, the results (rejection at P: = 0.068) were certainly suggestive of sepioid nonmonophyly.

Actins↗

PubMed2000

NJML: a hybrid algorithm for the neighbor-joining and maximum-likelihood methods.

In the reconstruction of a large phylogenetic tree, the most difficult part is usually the problem of how to explore the topology space to find the optimal topology. We have developed a "divide-and-conquer" heuristic algorithm in which an initial neighbor-joining (NJ) tree is divided into subtrees at internal branches having bootstrap values higher than a threshold. The topology search is then conducted by using the maximum-likelihood method to reevaluate all branches with a bootstrap value lower than the threshold while keeping the other branches intact. Extensive simulation showed that our simple method, the neighbor-joining maximum-likelihood (NJML) method, is highly efficient in improving NJ trees. Furthermore, the performance of the NJML method is nearly equal to or better than existing time-consuming heuristic maximum-likelihood methods. Our method is suitable for reconstructing relatively large molecular phylogenetic trees (number of taxa >/= 16).

Algorithms↗

PubMed1993

Phylogenetic analysis of the outer-membrane-protein genes of Chlamydiae, and its implication for vaccine development.

Examination of 18 complete and 6 partial sequences of the major outer-membrane protein from 24 chlamydiae isolates was used to reconstruct their evolutionary relationships. From this analysis, assuming that the clades with 100% bootstrap support are correct, come the following conclusions: (1) The tree of these sequences is not congruent with the phylogeny of the hosts, and thus host switching would seem to have occurred, thereby limiting the extent to which there has been coevolution of parasite and host. (2) The tree is also noncongruent with clustering by type of cell infected, thereby limiting the extent to which there has been coevolution of parasite and the cell type that it infects. (3) The tree is also noncongruent with clustering by the organ infected (eyes or genitalia), thereby limiting the extent to which there has been coevolution of parasite and the organ that it infects. (4) The tree is also noncongruent with genital strains arising from lymphogranuloma venereum strains. (5) The tree is also noncongruent with the geographic site at which the isolates were obtained, thereby limiting the extent of divergence explained by geographic separation. (6) There are estimated to be 185 amino acid positions that are invariable (as opposed to unvaried) in the major outer-membrane protein. There are 10 unvaried positions in the variable domains, of which 9 appear to be invariable, giving some reason to hope that development of a vaccine might be possible. (7) The rate of change of this protein is too small to see increased divergence over the time span of isolation of these genes, giving hope to any vaccine having longevity. Bootstrapping supports those portions of the tree on which the first five conclusions above depend. The picture that these results provide is more one of pathogen versatility than one of coevolutionary constraints. In addition, we examined 10 60-KDa, outer-membrane protein-2 genes, all but one of which were from these same strains. The tree was not, among the trachomatis strains, congruent with the major-outer-membrane protein tree, suggesting that gene exchange could be occurring among strains. Moreover, there is an apparent slowdown in divergence in this gene, among the trachomatis strains.

Amino Acid Sequence↗