Phylogeny Derived from Coding Retroviral Genome Organization

When divergence between viral species is large, the analysis and comparison of nucleotide or protein sequences are dependent on mutation biases and multiple substitutions per site leading, among other things, to the underestimation of branch lengths in phylogenetic trees. To avoid the problem of multiply substituted sites, a method not directly based on the nucleic or protein sequences has been applied to retroviruses. It consisted of asking questions about genome structure or organization, and gene function, the series of answers creating coded sequences analyzed by phylogenic software. This method recovered the principal retroviral groups such as the lentiviruses and spumaviruses and highlighted questions and answers characteristic of each group of retroviruses. In general, there was reasonable concordance between the coded genome methodology and that based on conventional phylogeny of the integrase protein sequence, indicating that integrase was fixing mutations slowly enough to marginalize the problem of multiple substitutions at sites. To a first approximation, this suggests that the acquisition of novel genetic features generally parallels the fixation of amino acid substitutions. Received: 18 May 2001 / Accepted: 7 September 2001     

Evidence for the Role of Recombination in the Regulatory Evolution of Saccharomyces cerevisiae Ty Elements

The recent completion of the sequencing of the Saccharomyces cerevisiae genome provides a unique opportunity to analyze the evolutionary relationships existing among the entire complement of retrotransposons residing within a single genome. In this article we report the results of such an analysis of two closely related families of yeast long terminal repeat (LTR) retrotransposons, Ty1 and Ty2. In our study, we analyzed the molecular variation existing among the 32 Ty1 and 13 Ty2 elements present within the S. cerevisiae genome recently sequenced within the context of the yeast genome project. Our results indicate that while the Ty1 family is most likely ancestral to Ty2 elements, both families of elements are relatively recent components of the S. cerevisiae genome. Our results also indicate that both families of elements have been subject to purifying selection within their protein coding regions. Finally, and perhaps most interestingly, our results indicate that a relatively recent recombination event has occurred between Ty2 and a subclass of Ty1 elements involving the LTR regulatory region. We discuss the possible biological significance of these findings and, in particular, how they contribute to a better overall understanding of LTR retrotransposon evolution. Received: 30 September 1997 / Accepted: 3 February 1998     

The Small-Subunit rRNA Gene Sequences of Venerids and the Phylogeny of Bivalvia

The complete nucleotide sequence of the 18S subunit of ribosomal DNA (rDNA) was determined for the venerid clams Callista chione (Pitarinae) and Venus verrucosa (Venerinae). Comparison of the new sequences with the published sequences of 1 annelid, 2 gastropods, 2 polyplacophorans, and 19 bivalves showed that when the annelids are used as outgroup the gastropods diverge from the bivalves, which form a cluster including the polyplacophorans. When the gastropods alone were compared with the bivalves, the latter split in two groups corresponding to the two subclasses of Heterodonta and Pteriomorpha. The former include two taxa that diverged early, Galeomma and Tridacna, while the Veneridae and Mactridae form two sister groups. In contrast to previous reports and in line with morphological data, the Ostreidae are included in the Pteriomorphia and form a monophyletic group. Received: 16 May 1998 / Accepted: 11 August 1998     

Bootstrap Confidence Levels for HIV-1 Recombination

Recombination has been invoked to explain the disparate evolutionary relationships observed for different genes or sequence segments of a single HIV-1 genome. We present a new method of assessing confidence in HIV-1 recombination as an alternative to the segment-by-segment nonparametric bootstrap commonly applied to confirm HIV-1 recombinant data. Our new method uses the bias-corrected accelerated percentile interval (BC_a) bootstrap method as applied to the ``problem of regions' (Efron and Tibshirani 1998). It is an extension of the BC_a method used in the inference of evolutionary relationships (Efron et al. 1996). This method has two advantages over the traditional bootstrap procedure: (1) it gives a single overall confidence measure rather than segment-by-segment results, and (2) it is more accurate. We test our method on 61 sequences, including 16 with ambiguous recombinant status. Received: 9 August 2000 / Accepted: 6 July 2001     

Phylogeny of Japanese Papilionid Butterflies Inferred from Nucleotide Sequences of the Mitochondrial ND5 Gene

Phylogenetic relationships among the Japanese papilionid butterflies were analyzed by comparing 783 nucleotide sequences of the mitochondrial gene encoding NADH dehydrogenase subunit 5 (ND5). Phylogenetic trees of the representative species from each family in the superfamily Papilionoidea revealed that the species of the family Papilionidae and those of all other families formed distinct clusters, with a few species of the family Hesperiidae (Hesperioidea) as an outgroup. In the phylogenetic trees of most Japanese species of the family Papilionidae with Nymphalis xanthomelas (Nymphalidae) as an outgroup, the tribe Parnassiini (Parnassiinae) formed a cluster, and the rest formed the other cluster in which the tribe Zerynthiini (Parnassiinae) and the subfamily Papilioninae formed different subclusters. In the Papilioninae cluster, the tribes Troidini and Graphiini formed a subcluster, and the tribe Papilionini formed the other subcluster. These results generally agree with the traditional classification of the papilionid butterflies based on their morphological characteristics and support the proposed evolutionary genealogy of the butterflies based on their morphology, behavior, and larval host plants, except that the tribes Parnasiini and Zerynthiini (both Parnassiinae) are not in the same cluster. Received: 16 March 1998 / Accepted: 28 April 1998     

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 Probabilistic Treatment of Phylogeny and Sequence Alignment

Carrying out simultaneous tree-building and alignment of sequence data is a difficult computational task, and the methods currently available are either limited to a few sequences or restricted to highly simplified models of alignment and phylogeny. A method is given here for overcoming these limitations by Bayesian sampling of trees and alignments simultaneously. The method uses a standard substitution matrix model for residues together with a hidden Markov model structure that allows affine gap penalties. It escapes the heavy computational burdens of other models by using an approximation called the ``*' rule, which replaces missing data by a sum over all possible values of variables. The behavior of the model is demonstrated on test sets of globins. Received: 25 May 1998 / Accepted: 8 December 1998     

The Evolution of Hexamerins and the Phylogeny of Insects

The evolutionary relationships among arthropod hemocyanins and insect hexamerins were investigated. A multiple sequence alignment of 12 hemocyanin and 31 hexamerin subunits was constructed and used for studying sequence conservation and protein phylogeny. Although hexamerins and hemocyanins belong to a highly divergent protein superfamily and only 18 amino acid positions are identical in all the sequences, the core structures of the three protein domains are well conserved. Under the assumption of maximum parsimony, a phylogenetic tree was obtained that matches perfectly the assumed phylogeny of the insect orders. An interesting common clade of the hymenopteran and coleopteran hexamerins was observed. In most insect orders, several paralogous hexamerin subclasses were identified that diversified after the splitting of the major insect orders. The dipteran arylphorin/LSP-1-like hexamerins were subject to closer examination, demonstrating hexamerin gene amplification and gene loss in the brachyceran Diptera. The hexamerin receptors, which belong to the hexamerin/hemocyanin superfamily, diverged early in insect evolution, before the radiation of the winged insects. After the elimination of some rapidly or slowly evolving sequences, a linearized phylogenetic tree of the hexamerins was constructed under the assumption of a molecular clock. The inferred time scale of hexamerin evolution, which dates back to the Carboniferous, agrees with the available paleontological data and reveals some previously unknown divergence times among and within the insect orders. Received: 4 August 1997 / Accepted: 29 October 1997     

Actin Phylogeny Identifies Mesostigma viride as a Flagellate Ancestor of the Land Plants

Green algae and land plants trace their evolutionary history to a unique common ancestor. This ``green lineage' is phylogenetically subdivided into two distinct assemblages, the Chlorophyta and the Streptophyta. The Chlorophyta includes the Chlorophyceae, Trebouxiophyceae, Ulvophyceae, and Prasinopohyceae, whereas the Streptophyta includes the Charophyceae plus the bryophytes, ferns, and all other multicellular land plants (Embryophyta). The Prasinophyceae is believed to contain the earliest divergences within the green lineage. Phylogenetic analyses using rDNA sequences identify the prasinophytes as a paraphyletic taxon that diverges at the base of the Chlorophyta. rDNA analyses, however, provide ambiguous results regarding the identity of the flagellate ancestor of the Streptophyta. We have sequenced the actin-encoding cDNAs from Scherffelia dubia (Prasinophyceae), Coleochaete scutata, Spirogyra sp. (Charophyceae), and the single-copy actin gene from Mesostigma viride (Prasinophyceae). Phylogenetic analyses show Mesostigma to be the earliest divergence within the Streptophyta and provide direct evidence for a scaly, biflagellate, unicellular ancestor for this lineage. This result is supported by the existence of two conserved actin-coding region introns (positions 20-3, 152-1), and one intron in the 5′-untranslated region of the actin gene shared by Mesostigma and the embryophytes. Received: 10 July 1997 / Accepted: 9 April 1998     

Molecular Phylogeny and Evolution of the Neurotrophins from Monotremes and Marsupials

We have investigated the phylogenetic relationships of monotremes and marsupials using nucleotide sequence data from the neurotrophins; nerve growth factor (NGF), brain-derived neurotrophic factor (BDNF), and neurotrophin-3 (NT-3). The study included species representing monotremes, Australasian marsupials and placentals, as well as species representing birds, reptiles, and fish. PCR was used to amplify fragments encoding parts of the neurotrophin genes from echidna, platypus, and eight marsupials from four different orders. Phylogenetic trees were generated using parsimony analysis, and support for the different tree structures was evaluated by bootstrapping. The analysis was performed with NGF, BDNF, or NT-3 sequence data used individually as well as with the three neurotrophins in a combined matrix, thereby simultaneously considering phylogenetic information from three separate genes. The results showed that the monotreme neurotrophin sequences associate to either therian or bird neurotrophin sequences and suggests that the monotremes are not necessarily related closer to therians than to birds. Furthermore, the results confirmed the present classification of four Australasian marsupial orders based on morphological characters, and suggested a phylogenetic relationship where Dasyuromorphia is related closest to Peramelemorphia followed by Notoryctemorphia and Diprotodontia. These studies show that sequence data from neurotrophins are well suited for phylogenetic analysis of mammals and that neurotrophins can resolve basal relationships in the evolutionary tree. Received: 27 January 1997 / Accepted: 20 March 1997     

Nuclear and Nucleomorph SSU rDNA Phylogeny in the Cryptophyta and the Evolution of Cryptophyte Diversity

The plastid-bearing members of the Cryptophyta contain two functional eukaryotic genomes of different phylogenetic origin, residing in the nucleus and in the nucleomorph, respectively. These widespread and diverse protists thus offer a unique opportunity to study the coevolution of two different eukaryotic genomes within one group of organisms. In this study, the SSU rRNA genes of both genomes were PCR-amplified with specific primers and phylogenetic analyses were performed on different data sets using different evolutionary models. The results show that the composition of the principal clades obtained from the phylogenetic analyses of both genes was largely congruent, but striking differences in evolutionary rates were observed. These affected the topologies of the nuclear and nucleomorph phylogenies differently, resulting in long-branch attraction artifacts when simple evolutionary models were applied. Deletion of long-branch taxa stabilized the internal branching order in both phylogenies and resulted in a completely resolved topology in the nucleomorph phylogeny. A comparison of the tree topologies derived from SSU rDNA sequences with characters previously used in cryptophyte systematics revealed that the biliprotein type was congruent, but the type of inner periplast component incongruent, with the molecular trees. The latter is indicative of a hidden cellular dimorphism (cells with two periplast types present in a single clonal strain) of presumably widespread occurrence throughout cryptophyte diversity, which, in consequence, has far-reaching implications for cryptophyte systematics as it is practiced today.     

A Molecular Phylogeny of Lilium in the Internal Transcribed Spacer Region of Nuclear Ribosomal DNA

Phylogenetic relationships among 55 species of Lilium, Cardiocrinum giganteum, and Nomocharis saluenensis were inferred from nucleotide sequence variations in the internal transcribed spacer (ITS) regions of 18S–25S nuclear ribosomal DNA. The phylogeny derived from ITS sequences estimated using maximum-likelihood methods indicated that (1) most of the species construct their own clade according to the classification based on morphological features at the section level; (2) section Daurolirion is not independent of Sinomartagon, and it is appropriate to integrate two sections as Sinomartagon; (3) it is appropriate that L. henryi and L. bulbiferum are classified into subsection 6a and Sinomartagon–Daurolirion, respectively; (4) subsection 6b is much closer to Sinomartagon than subsection 6a and Archelirion, and it arose directly from Sinomartagon; and (5) Lilium is much closer to Nomocharis than Cardiocrinum. Phylogenetic estimation using sequences of the ITS region is suitable at the levels of genus, section, and most of subsection. Received: 18 December 1998 / Accepted: 14 March 1999     

Mitochondrial DNA Phylogeny of the Family Cichlidae: Monophyly and Fast Molecular Evolution of the Neotropical Assemblage

A mitochondrial DNA (mtDNA) phylogeny of cichlid fish is presented for the most taxonomically inclusive data set compiled to date (64 taxa). 16S rDNA data establish with confidence relationships among major lineages of cichlids, with a general pattern congruent with previous morphological studies and less inclusive molecular phylogenies based on nuclear genes. Cichlids from Madagascar and India are the most basal groups of the family Cichlidae and sister to African–Neotropical cichlids. The cichlid phylogeny suggests drift-vicariance events, consistent with the fragmentation of Gondwana, to explain current biogeographic distributions. Important phylogenetic findings include the placement of the controversial genus Heterochromis basal among African cichlids, the South American genus Retroculus as the most basal taxon of the Neotropical cichlid assemblage, and the close relationship of the Neotropical genera Cichla with Astronotus rather than with the crenicichlines. Based on a large number of South American genera, the Neotropical cichlids are defined as a monophyletic assemblage and shown to harbor significantly higher levels of genetic variation than their African counterparts. Relative rate tests suggest that Neotropical cichlids have experienced accelerated rates of molecular evolution. But these high evolutionary rates were significantly higher among geophagine cichlids. Received: 18 September 1998 / Accepted: 16 December 1998     

Defining the Core of Nontransferable Prokaryotic Genes: The Euryarchaeal Core

If lateral gene transfer (LGT) has affected all genes over the course of prokaryotic evolution, reconstruction of organismal phylogeny is compromised. However, if a core of genes is immune to transfer, then the evolutionary history of that core might be our most reliable guide to the evolution of organisms. Such a core should be preferentially included in the subset of genes shared by all organisms, but where universally conserved genes have been analyzed, there is too little phylogenetic signal to allow determination of whether or not they indeed have the same history (Hansmann and Martin 2000; Teichmann and Mitchison 1999). Here we look at a more restricted set, 521 homologous genes (COGs) simultaneously present in four sequenced euryarchaeal genomes. Although there is overall little robust phylogenetic signal in this data set, there is, among well-supported trees, strong representation of all three possible four-taxon topologies. ``Informational' genes seem no less subject to LGT than are ``operational genes,' within the euryarchaeotes. We conclude that (i) even in this collection of conserved genes there has been extensive LGT (orthologous gene replacement) and (ii) the notion that there is a core of nontransferable genes (the ``core hypothesis') has not been proven and may be unprovable. Received: 7 November 2000 / Accepted: 20 February 2001     

Evolution of 28S Ribosomal DNA in Chaetognaths: Duplicate Genes and Molecular Phylogeny

The chaetognaths are an extraordinarily homogeneous phylum of animals at the morphological level, with a bauplan that can be traced back to the Cambrian. Despite the attention of zoologists for over two centuries, there is little agreement on classification within the phylum. We have used a molecular biological approach to investigate the phylogeny of extant chaetognaths. A rapidly evolving expansion segment toward the 5′ end of 28S ribosomal DNA (rDNA) was amplified using the polymerase chain reaction (PCR), cloned, and sequenced from 26 chaetognath samples representing 18 species. An unusual finding was the presence of two distinct classes of 28S rDNA gene in chaetognaths; our analyses suggest these arose by a gene (or gene cluster) duplication in a common ancestor of extant chaetognaths. The two classes of chaetognath 28S rDNA have been subject to different rates of molecular evolution; we present evidence that both are expressed and functional. In phylogenetic reconstructions, the two classes of 28S rDNA yield trees that root each other; these clearly demonstrate that the Aphragmophora and Phragmophora are natural groups. Within the Aphragmophora, we find good support for the groupings denoted Solidosagitta, Parasagitta, and Pseudosagitta. The relationships between several well-supported groups within the Aphragmophora are uncertain; we suggest this reflects rapid, recent radiation during chaetognath evolution. Received: 19 March 1996 / Accepted: 5 August 1996     

Characterization of the Hydra Lamin and Its Gene: A Molecular Phylogeny of Metazoan Lamins

We report sequences for nuclear lamins from the teleost fish Danio and six invertebrates. These include two cnidarians (Hydra and Tealia), one priapulid, two echinoderms, and the cephalochordate Branchiostoma. Combining these results with earlier data on Drosophila, Caenorhabditis elegans, and various vertebrates, the following conclusions on lamin evolution can be drawn. First, all invertebrate lamins resemble in size the vertebrate B-type lamin. Second, all lamins described previously for amphibia, birds and mammals as well as the first lamin of a fish, characterized here, show a cluster of 7 to 12 acidic residues in the tail domain. Since this acidic cluster is absent from all invertebrate lamins including that of the cephalochordate Branchiostoma, it was acquired with the vertebrate lineage. The larger A-type lamin of differentiated cells must have arisen subsequently by gene duplication and insertion of an extra exon. This extra exon of the vertebrate A-lamins is the only major change in domain organization in metazoan lamin evolution. Third, the three introns of the Hydra and Priapulus genes correspond in position to the last three introns of vertebrate B-type lamin genes. Thus the entirely different gene organization of the C. elegans and Drosophila Dmo genes seems to reflect evolutionary drift, which probably also accounts for the fact that C. elegans has the most diverse lamin sequence. Finally we discuss the possibility that two lamin types, a constitutively expressed one and a developmentally regulated one, arose independently on the arthropod and vertebrate lineages. Received: 4 February 1999 / Accepted: 1 April 1999     

Phylogeny, Genome Evolution, and Host Specificity of Single-Stranded RNA Bacteriophage (Family Leviviridae)

Bacteriophage of the family Leviviridae have played an important role in molecular biology where representative species, such as Qβ and MS2, have been studied as model systems for replication, translation, and the role of secondary structure in gene regulation. Using nucleotide sequences from the coat and replicase genes we present the first statistical estimate of phylogeny for the family Leviviridae using maximum-likelihood and Bayesian estimation. Our analyses reveal that the coliphage species are a monophyletic group consisting of two clades representing the genera Levivirus and Allolevivirus. The Pseudomonas species PP7 diverged from its common ancestor with the coliphage prior to the ancient split between these genera and their subsequent diversification. Differences in genome size, gene composition, and gene expression are shown with a high probability to have changed along the lineage leading to the Allolevivirus through gene expansion. The change in genome size of the Allolevivirus ancestor may have catalyzed subsequent changes that led to their current genome organization and gene expression. Received: 3 March 2000 / Accepted: 17 October 2000     

The Archaea Monophyly Issue: A Phylogeny of Translational Elongation Factor G(2) Sequences Inferred from an Optimized Selection of Alignment Positions

A global alignment of EF-G(2) sequences was corrected by reference to protein structure. The selection of characters eligible for construction of phylogenetic trees was optimized by searching for regions arising from the artifactual matching of sequence segments unique to different phylogenetic domains. The spurious matchings were identified by comparing all sections of the global alignment with a comprehensive inventory of significant binary alignments obtained by BLAST probing of the DNA and protein databases with representative EF-G(2) sequences. In three discrete alignment blocks (one in domain II and two in domain IV), the alignment of the bacterial sequences with those of Archaea–Eucarya was not retrieved by database probing with EF-G(2) sequences, and no EF-G homologue of the EF-2 sequence segments was detected by using partial EF-G(2) sequences as probes in BLAST/FASTA searches. The two domain IV regions (one of which comprises the ADP-ribosylatable site of EF-2) are almost certainly due to the artifactual alignment of insertion segments that are unique to Bacteria and to Archaea–Eucarya. Phylogenetic trees have been constructed from the global alignment after deselecting positions encompassing the unretrieved, spuriously aligned regions, as well as positions arising from misalignment of the G′ and G″ subdomain insertion segments flanking the ``fifth' consensus motif of the G domain (?varsson, 1995). The results show inconsistencies between trees inferred by alternative methods and alternative (DNA and protein) data sets with regard to Archaea being a monophyletic or paraphyletic grouping. Both maximum-likelihood and maximum-parsimony methods do not allow discrimination (by log-likelihood difference and difference in number of inferred substitutions) between the conflicting (monophyletic vs. paraphyletic Archaea) topologies. No specific EF-2 insertions (or terminal accretions) supporting a crenarchaeal–eucaryal clade are detectable in the new EF-G(2) sequence alignment.     

Estimating the Transition/Transversion Ratio from Independent Pairwise Comparisons with an Assumed Phylogeny

A method is presented for estimating the transition/transversion ratio (TI/TV), based on phylogenetically independent comparisons. TI/TV is a parameter of some models used in phylogeny estimation intended to reflect the fact that nucleotide substitutions are not all equally likely. Previous attempts to estimate TI/TV have commonly faced three problems: (1) few taxa; (2) nonindependence among pairwise comparisons; and (3) multiple hits make the apparent TI/TV between two sequences decrease over time since their divergence, giving a misleading impression of relative substitution probabilities. We have made use of the time dependency, modeling how the observed TI/TV changes over time and extrapolating to estimate the ``instantaneous' TI/TV—the relevant parameter for phylogenetic inference. To illustrate our method, TI/TV was estimated for two mammalian mitochondrial genes. For 26 pairs of cytochrome b sequences, the estimate of TI/TV was 5.5; 16 pairs of 12s rRNA yielded an estimate of 9.5. These estimates are higher than those given by the maximum likelihood method and than those obtained by averaging all possible pairwise comparisons (with or without a two-parameter correction for multiple substitutions). We discuss strengths, weaknesses, and further uses of our method. Received: 22 August 1995 / Accepted: 26 July 1996     

Contributions to the Phylogeny of the Cyclophyllidea (Cestoda) Inferred from Mitochondrial 12S rDNA

A 314-bp fragment of the mitochondrial 12S rRNA gene from 21 cestodes species of eight families was synthesized by PCR with specially designed primers. These allowed amplification of parasite DNA without concomitant synthesis of host DNA. Phylogenetic trees were inferred from the sequence data using three methods (maximum parsimony, maximum likelihood, and Fitch–Margoliash). At the major nodes all three trees were similar. For the first time the genus Mesocestoides could be arranged into the Cyclophyllidea and a narrow relationship between the Mesocestoididae, Taeniidae, Hymenolepididae, Anoplocephalidae, and Dipylidiidae was shown. Members of the families Catenotaeniidae and a cluster of two families (Hymenolepididae and Dilepididae) form two monophyletic groups which derive prior to the remaining families of this phylogenetic study. A third and a fourth clear monophyletic group were formed by the Taeniidae and by the Mesocestoididae. A high degree of variation within the examined 304-bp fragment was observed between two isolates of Taenia taeniaeformis, supporting often discussed genetic heterogeneity within this species. In contrast, only one nucleotide exchange was found in 23 isolates of Echinococcus multilocularis of various geographic origin, indicating that this species is genetically homogenous. Received: 1 October 1997 / Accepted: 4 December 1997