Phylogeny and Rates of Molecular Evolution of Planktonic Foraminifera: SSU rDNA Sequences Compared to the Fossil Record

Planktonic foraminifera are marine protists, whose calcareous shells form oceanic sediments and are widely used for stratigraphic and paleoenvironmental analyses. The fossil record of planktonic foraminifera is compared here to their molecular phylogeny inferred from ribosomal DNA sequences. Eighteen partial SSU rDNA sequences from species representing all modern planktonic families (Globigerinidae, Hastigerinidae, Globorotaliidae, Candeinidae) were obtained and compared to seven sequences representing the major groups of benthic foraminifera. The phylogenetic analyses indicate a polyphyletic origin for the planktonic foraminifera. The Candeinidae, the Globorotaliidae, and the clade Globigerinidae + Hastigerinidae seem to have originated independently, at different epochs in the evolution of foraminifera. Inference of their relationships, however, is limited by substitution rates of heterogeneity. Rates of SSU rDNA evolution vary from 4.0 × 10⁻⁹ substitutions/site/year in the Globigerinidae to less than 1.0 × 10⁻⁹ substitutions/site/year in the Globorotaliidae. These variations may be related to different levels of adaptation to the planktonic mode of life. A clock-like evolution is observed among the Globigerinidae, for which molecular and paleontological data are congruent. Phylogeny of the Globorotaliidae is clearly biased by rapid rates of substitution in two species (G. truncatulinoides and G. menardii). Our study reveals differences in absolute rates of evolution at all taxonomic levels in planktonic foraminifera and demonstrates their effect on phylogenetic reconstructions. Received: 21 January 1997 / Accepted: 17 April 1997     

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     

Phylogeny of Nematoda and Cephalorhyncha Derived from 18S rDNA

Phylogenetic relationships of nematodes, nematomorphs, kinorhynchs, priapulids, and some other major groups of invertebrates were studied by 18S rRNA gene sequencing. Kinorhynchs and priapulids form the monophyletic Cephalorhyncha clade that is the closest to the coelomate animals. When phylogenetic trees were generated by different methods, the position of nematomorphs appeared to be unstable. Inclusion of Enoplus brevis, a representative of a slowly evolving nematode lineage, in the set of analyzed species refutes the tree patterns, previously derived from molecular data, where the nematodes appear as a basal bilateral lineage. The nematodes seem to be closer to the coelomate animals than was speculated earlier. According to the results obtained, nematodes, nematomorphs, tardigrades, arthropods, and cephalorhynchs are a paraphyletic association of closely related taxa. Received: 1 December 1997 / Accepted: 9 April 1998     

Dinoflagellate nuclear SSU rRNA phylogeny suggests multiple plastid losses and replacements

Dinoflagellates are a trophically diverse group of protists with photosynthetic and non-photosynthetic members that appears to incorporate and lose endosymbionts relatively easily. To trace the gain and loss of plastids in dinoflagellates, we have sequenced the nuclear small subunit rRNA gene of 28 photosynthetic and four non-photosynthetic species, and produced phylogenetic trees with a total of 81 dinoflagellate sequences. Patterns of plastid gain, loss, and replacement were plotted onto this phylogeny. With the exception of the apparently early-diverging Syndiniales and Noctilucales, all non-photosynthetic dinoflagellates are very likely to have had photosynthetic ancestors with peridinin-containing plastids. The same is true for all dinoflagellates with plastids other than the peridinin-containing plastid: their ancestors have replaced one type of plastid for another, in some cases most likely through a non-photosynthetic intermediate. Eight independent instances of plastid loss and three of replacement can be inferred from existing data, but as more non-photosynthetic lineages are characterized these numbers will surely grow. Received: 25 September 2000 / Accepted: 24 April 2001     

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     

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     

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     

Mitochondrial DNA Phylogeny and the Evolution of Host-Plant Use in Palearctic Chrysolina (Coleoptera, Chrysomelidae) Leaf Beetles

The genus Chrysolina consists of specialized phytophagous leaf-beetles (Coleoptera, Chrysomelidae) with feed on several plant families. There is no explicit phylogenetic hypothesis available for this genus, which includes 65 subgenera and more than 400 species with a wide distribution. We obtained 839-bp sequence data from the 16S rDNA and cytochrome oxidase subunit I (COI) mitochondrial genes. Thirty Chrysolina taxa representing eight host–plant affiliations, two species of the closely related genus Oreina, and two outgroups were sampled. These data sets were used separately and combined to obtain the mitochondrial cladogram of the group using maximum-parsimony and maximum-likelihood criteria. The results were compared to current proposals for Chrysolina systematics that are based on morphological, ecological, and karyological data. The trees obtained were in the most part congruent with the proposed ancestral association of Chrysolina to Lamiaceae based on chromosome number in several lineages. A minimum of five host-plant switches from the ancestral state inferred at the family level and two at the subclass level suggests the absence of parallel evolution of beetles and their host plants. Another switch leading to oligophagy at the family level was deduced to have occurred in the lineage of the subgenus Chrysolina s.str. Received: 22 May 1998 / Accepted: 16 September 1998     

The Phylogeny of Glyceraldehyde-3-Phosphate Dehydrogenase Indicates Lateral Gene Transfer from Cryptomonads to Dinoflagellates

Sequence analysis of two nuclear-encoded glyceraldehyde-3-phosphate dehydrogenase (GAPDH) genes isolated from the dinoflagellate Gonyaulax polyedra distinguishes them as cytosolic and chloroplastic forms of the enzyme. Distance analysis of the cytosolic sequence shows the Gonyaulax gene branching early within the cytosolic clade, consistent with other analyses. However, the plastid sequence forms a monophyletic group with the plastid isoforms of cryptomonads, within an otherwise cytosolic clade, distinct from all other plastid GAPDHs. This is attributed to lateral gene transfer from an ancestral cryptomonad to a dinoflagellate, providing the first example of genetic exchange accompanying symbiotic associations between the two, which are common in present day cells. Received: 29 April 1998 / Accepted: 7 July 1998     

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     

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     

Multiple Substitutions Affect the Phylogenetic Utility of Cytochrome b and 12S rDNA Data: Examining a Rapid Radiation in Leporid (Lagomorpha) Evolution

Partial sequences of two mitochondrial genes, the 12S ribosomal gene (739 bp) and the cytochrome b gene (672 bp), were analyzed in hopes of reconstructing the evolutionary relationships of 11 leporid species, representative of seven genera. However, partial cytochrome b sequences were of little phylogenetic value in this study. A suite of pairwise comparisons between taxa revealed that at the intergeneric level, the cytochrome b gene is saturated at synonymous coding positions due to multiple substitution events. Furthermore, variation at the nonsynonymous positions is limited, rendering the cytochrome b gene of little phylogenetic value for assessing the relationships between leporid genera. If the cytochrome b data are analyzed without accounting for these two classes of nucleotides (i.e., synonymous and nonsynonymous sites), one may incorrectly conclude that signal exists in the cytochrome b data. The mitochondrial 12S rRNA gene, on the other hand, has not experienced excessive saturation at either stem or loop positions. Phylogenies reconstructed from the 12S rDNA data support hypotheses based on fossil evidence that African rock rabbits (Pronolagus) are outside of the main leporid stock and that leporids experienced a rapid radiation. However, the molecular data suggest that this radiation event occurred in the mid-Miocene several millions of years earlier than the Pleistocene dates suggested by paleontological evidence. Received: 23 April 1998 / Accepted: 14 May 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     

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 Effect of Recombination on the Accuracy of Phylogeny Estimation

Phylogenetic studies based on DNA sequences typically ignore the potential occurrence of recombination, which may produce different alignment regions with different evolutionary histories. Traditional phylogenetic methods assume that a single history underlies the data. If recombination is present, can we expect the inferred phylogeny to represent any of the underlying evolutionary histories? We examined this question by applying traditional phylogenetic reconstruction methods to simulated recombinant sequence alignments. The effect of recombination on phylogeny estimation depended on the relatedness of the sequences involved in the recombinational event and on the extent of the different regions with different phylogenetic histories. Given the topologies examined here, when the recombinational event was ancient, or when recombination occurred between closely related taxa, one of the two phylogenies underlying the data was generally inferred. In this scenario, the evolutionary history corresponding to the majority of the positions in the alignment was generally recovered. Very different results were obtained when recombination occurred recently among divergent taxa. In this case, when the recombinational breakpoint divided the alignment in two regions of similar length, a phylogeny that was different from any of the true phylogenies underlying the data was inferred.     

Short Retroposons of the B2 Superfamily: Evolution and Application for the Study of Rodent Phylogeny

Short retroposons can be used as natural phylogenetic markers. By means of hybridization and PCR analysis, we demonstrate that B2 retroposon copies are present only in the three rodent families: Muridae, Cricetidae, and Spalacidae. This observation highlights the close phylogenetic relation between these families. Two novel B2-related retroposon families, named DIP and MEN elements, are described. DIP elements are found only in the genomes of jerboas (family Dipodidae) and birch mice (family Zapodidae), demonstrating the close relationship between these rodents. MEN element copies were isolated from the squirrel, Menetes berdmorei, but were not detected in three other species from the family Sciuridae. The MEN element has an unusual dimeric structure: the left and right monomers are B2- and B1-related sequences, respectively. Comparison of the B2, DIP, MEN, and 4.5S₁ RNA elements revealed an 80-bp core sequence located at the beginning of the B2 superfamily retroposons. This observation suggests that these retroposon families descended from a common progenitor. A likely candidate for this direct progenitor could be the ID retroposon. Received: 20 December 1996 / Accepted: 17 June 1997     

A Phylogenetic Assessment of the Eukaryotic Light-Harvesting Antenna Proteins, with Implications for Plastid Evolution

The light-harvesting complexes (LHCs) are a superfamily of chlorophyll-binding proteins present in all photosynthetic eukaryotes. The Lhc genes are nuclear-encoded, yet the pigment–protein complexes are localized to the thylakoid membrane and provide a marker to follow the evolutionary paths of plastids with different pigmentation. The LHCs are divided into the chlorophyll a/b-binding proteins of the green algae, euglenoids, and higher plants and the chlorophyll a/c-binding proteins of various algal taxa. This work examines the phylogenetic position of the LHCs from three additional taxa: the rhodophytes, the cryptophytes, and the chlorarachniophytes. Phylogenetic analysis of the LHC sequences provides strong statistical support for the clustering of the rhodophyte and cryptomonad LHC sequences within the chlorophyll a/c-binding protein lineage, which includes the fucoxanthin–chlorophyll proteins (FCP) of the heterokonts and the intrinsic peridinin–chlorophyll proteins (iPCP) of the dinoflagellates. These associations suggest that plastids from the heterokonts, haptophytes, cryptomonads, and the dinoflagellate, Amphidinium, evolved from a red algal-like ancestor. The Chlorarachnion LHC is part of the chlorophyll a/b-binding protein assemblage, consistent with pigmentation, providing further evidence that its plastid evolved from a green algal secondary endosymbiosis. The Chlorarachnion LHC sequences cluster with the green algal LHCs that are predominantly associated with photosystem II (LHCII). This suggests that the green algal endosymbiont that evolved into the Chlorarachnion plastid was acquired following the emergence of distinct LHCI and LHCII complexes. Received: 25 February 1998 / Accepted: 13 May 1998     

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