Higher-level phylogeny of Foraminifera inferred from the RNA polymerase II (RPB1) gene

Macroevolutionary relations among main lineages of Foraminifera have traditionally been inferred from the small subunit ribosomal genes (SSU rDNA). However, important discrepancies in the rates of SSU rDNA evolution between major lineages led to difficulties in accurate interpretation of SSU-based phylogenetic reconstructions. Recently, actin and beta-tubulin sequences have been used as alternative markers of foraminiferal phylogeny and their analyses globally confirm results obtained with SSU rDNA. In order to test new protein markers, we sequenced a fragment of the largest subunit of the RNA polymerase II (RPB1), a nuclear encoded single copy gene, for 8 foraminiferal species representing major orders of Foraminifera. Analyses of our data robustly confirm previous SSU rDNA and actin phylogenies and show (i) the paraphyly and ancestral position of monothalamid Foraminifera; (ii) the independent origin of miliolids; (iii) the monophyly of rotaliids, including buliminids and globigerinids; and (iv) the polyphyly of planktonic families Globigerinidae and Candeinidae. Additionally, the RPB1 phylogeny suggests Allogromiidae as the most ancestral foraminiferal lineage. In the light of our study, RPB1 appears as a valuable phylogenetic marker, particularly useful for groups of protists showing extreme variations of evolutionary rates in ribosomal genes.     

Extreme differences in rates of molecular evolution of foraminifera revealed by comparison of ribosomal DNA sequences and the fossil record

Foraminifera have one of the best known fossil records among the unicellular eukaryotes. However, the origin and phylogenetic relationships of the extant foraminiferal lineages are poorly understood. To test the current paleontological hypotheses on evolution of foraminifera, we sequenced about 1,000 base pairs from the 3' end of the small subunit rRNA gene (SSU rDNA) in 22 species representing all major taxonomic groups. Phylogenies were derived using neighbor- joining, maximum-parsimony, and maximum-likelihood methods. All analyses confirm the monophyletic origin of foraminifera. Evolutionary relationships within foraminifera inferred from rDNA sequences, however, depend on the method of tree building and on the choice of analyzed sites. In particular, the position of planktonic foraminifera shows important variations. We have shown that these changes result from the extremely high rate of rDNA evolution in this group. By comparing the number of substitutions with the divergence times inferred from the fossil record, we have estimated that the rate of rDNA evolution in planktonic foraminifera is 50 to 100 times faster than in some benthic foraminifera. The use of the maximum-likelihood method and limitation of analyzed sites to the most conserved parts of the SSU rRNA molecule render molecular and paleontological data generally congruent.      

Characterization of Mitochondrial Small-Subunit Ribosomal RNAs from Holoparasitic Plants

Mitochondrial small-subunit (19S) rDNA sequences were obtained from 10 angiosperms to further characterize sequence divergence levels and structural variation in this molecule. These sequences were derived from seven holoparasitic (nonphotosynthetic) angiosperms as well as three photosynthetic plants. 19S rRNA is composed of a conservative core region (ca. 1450 nucleotides) as well as two variable regions (V1 and V7). In pairwise comparisons of photosynthetic angiosperms to Glycine, the core 19S rDNA sequences differed by less than 1.4%, thus supporting the observation that variation in mitochondrial rDNA is 3–4 times lower than seen in protein coding and rDNA genes of other subcellular organelles. Sequences representing four distinct lineages of nonasterid holoparasites showed significantly increased numbers of substitutions in their core 19S rDNA sequences (2.3–7.6%), thus paralleling previous findings that showed accelerated rates in nuclear (18S) and plastid (16S) rDNA from the same plants. Relative rate tests confirmed the accelerated nucleotide substitution rates in the holoparasites whereas rates in nonparasitic plants were not significantly increased. Among comparisons of both parasitic and nonparasitic plants, transversions outnumbered transitions, in many cases more than two to one. The core 19S rRNA is conserved in sequence and structure among all nonparasitic angiosperms whereas 19S rRNA from members of holoparasitic Balanophoraceae have unique extensions to the V5 and V6 variable domains. Substitution and insertion/deletion mutations characterized the V1 and V7 regions of the nonasterid holoparasites. The V7 sequence of one holoparasite (Scybalium) contained repeat motifs. The cause of substitution rate increases in the holoparasites does not appear to be a result of RNA editing, hence the underlying molecular mechanism remains to be fully documented. Received: 18 May 1997 / Accepted: 11 July 1997     

Molecular versus Taxonomic Rates of Evolution in Planktonic Foraminifera

Neogene planktonic foraminifera are among the most widely used microfossils in the study of tempo and mode of evolution. Comparisons of taxonomic rates between the two major clades in this group have shown that the nonspinose globorotaliids have undergone a significantly more rapid evolutionary turnover than the spinose globigerinids (S. M. Stanleyet al., 1988,Paleobiology14, 235–249). In order to test if similar fluctuations are observed in molecular data, we have used different methods to calculate absolute and relative rates of substitutions based on 16 partial SSU rDNA sequences from representatives of both groups. According to our data, rates of substitution are relatively constant within the globigerinids with a mean value of 4.3 subst./site/10⁹years, but vary in the globorotaliid clade with three species having a rate of about 1 subst./site/10⁹years and two species evolving much faster with rates of more than 7 subst./site/10⁹years. Assuming that the fast rates result from recent accelerations, the globorotaliids have basically much slower molecular evolutionary rates than the globigerinids, in opposition to the fossil data.     

Substitution Model of Sequence Evolution for the Human Immunodeficiency Virus Type 1 Subtype B gp120 Gene over the C2-V5 Region

Phylogenetic analyses frequently rely on models of sequence evolution that detail nucleotide substitution rates, nucleotide frequencies, and site-to-site rate heterogeneity. These models can influence hypothesis testing and can affect the accuracy of phylogenetic inferences. Maximum likelihood methods of simultaneously constructing phylogenetic tree topologies and estimating model parameters are computationally intensive, and are not feasible for sample sizes of 25 or greater using personal computers. Techniques that initially construct a tree topology and then use this non-maximized topology to estimate ML substitution rates, however, can quickly arrive at a model of sequence evolution. The accuracy of this two-step estimation technique was tested using simulated data sets with known model parameters. The results showed that for a star-like topology, as is often seen in human immunodeficiency virus type 1 (HIV-1) subtype B sequences, a random starting topology could produce nucleotide substitution rates that were not statistically different than the true rates. Samples were isolated from 100 HIV-1 subtype B infected individuals from the United States and a 620 nt region of the env gene was sequenced for each sample. The sequence data were used to obtain a substitution model of sequence evolution specific for HIV-1 subtype B env by estimating nucleotide substitution rates and the site-to-site heterogeneity in 100 individuals from the United States. The method of estimating the model should provide users of large data sets with a way to quickly compute a model of sequence evolution, while the nucleotide substitution model we identified should prove useful in the phylogenetic analysis of HIV-1 subtype B env sequences. Received: 4 October 2000 / Accepted: 1 March 2001     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

Molecular evolution of nitrate reductase genes

To understand the evolutionary mechanisms and relationships of nitrate reductases (NRs), the nucleotide sequences encoding 19 nitrate reductase (NR) genes from 16 species of fungi, algae, and higher plants were analyzed. The NR genes examined show substantial sequence similarity, particularly within functional domains, and large variations in GC content at the third codon position and intron number. The intron positions were different between the fungi and plants, but conserved within these groups. The overall and nonsynonymous substitution rates among fungi, algae, and higher plants were estimated to be 4.33 × 10⁻¹⁰ and 3.29 × 10⁻¹⁰ substitutions per site per year. The three functional domains of NR genes evolved at about one-third of the rate of the N-terminal and the two hinge regions connecting the functional domains. Relative rate tests suggested that the nonsynonymous substitution rates were constant among different lineages, while the overall nucleotide substitution rates varied between some lineages. The phylogenetic trees based on NR genes correspond well with the phylogeny of the organisms determined from systematics and other molecular studies. Based on the nonsynonymous substitution rate, the divergence time of monocots and dicots was estimated to be about 340 Myr when the fungi–plant or algae–higher plant divergence times were used as reference points and 191 Myr when the rice–barley divergence time was used as a reference point. These two estimates are consistent with other estimates of divergence times based on these reference points. The lack of consistency between these two values appears to be due to the uncertainty of the reference times. Received: 10 April 1995 / Accepted: 10 September 1995     

Variable Substitution Rates of the 18 Domain Sequences in Artemia Hemoglobin

The Artemia hemoglobin is a dimer comprising two nine-domain covalent polymers in quaternary association. Each polymer is encoded by a gene representing nine successive globin domains which have different sequences and are presumed to have been copied originally from a single-domain gene. Two different polymers exist as the result of a complete duplication of the nine-domain gene, allowing the formation of either homodimers or the heterodimer. The total population size of 18 domains comprising nine corresponding pairs, coupled with the probability that they reflect several hundred million years of evolution in the same lineage, provides a unique model in which the process of gene multiplication can be analyzed. The outcome has important implications for the reliability of local molecular clocks. The two polymers differ from each other at 11.7% of amino acid sites; however when corresponding individual domains are compared between polymers, amino acid substitution fluctuates by a factor of 2.7-fold from lowest to highest. This variation is not obvious at the DNA level: Domain pair identity values fluctuate by 1.3-fold. Identity values are, however, uncorrected for multiple substitutions, and both silent and nonsilent changes are pooled. Therefore, to determine the variability in relative substitution rates at the DNA level, we have used the method of Li (1993, J Mol Evol 36:96–99) to determine estimates of nonsynonymous (K _A ) and synonymous (K _S ) substitutions per site for the nine pairs of domains. As expected, the overall level of silent substitutions (K _S of 56.9%) far exceeded nonsilent substitutions (K _A of 6.7%); however, for corresponding domain pairs, K _A fluctuates by 2.3-fold and K _S by 1.7-fold. The large discrepancies reflected in the expressed protein have accrued within a single lineage and the implication is that divergence dates of different genera based on amino acid sequences, even with well-studied proteins of reasonable size, can be wrong by a factor well in excess of 2. Received: 4 June 1997 / Accepted: 17 December 1997     

Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates. Received: 15 December 1997 / Accepted: 24 February 1998     

A revised evolutionary history of hepatitis B virus (HBV)

Previous studies of the evolutionary history of hepatitis B virus (HBV) have been compromised by intergenotype recombination and complex patterns of nucleotide substitution, perhaps caused by differential selection pressures. We examined the phylogenetic distribution of recombination events among human HBV genotypes and found that genotypes A plus D, and genotypes B plus C, had distinct patterns of recombination suggesting differing epidemiological relationships among them. By analyzing the nonoverlapping regions of the viral genome we found strong bootstrap support for some intergenotypic groupings, with evidence of a division between human genotypes A–E from the viruses sampled from apes and human genotype F. However, the earliest events in the divergence of HBV remain uncertain. These uncertainties could not be explained by differential selection pressures, as the ratio of nonsynonymous-to-synonymous substitutions (d _N/d _S) did not vary extensively among lineages and there is no strong evidence for positive selection across the whole tree. Finally, we provide a new estimate of the mean substitution rate in HBV, 4.2 × 10⁻⁵, which suggests that divergence of HBV in humans and apes has occurred only in the last 6000 years.     

Complex Patterns of Plastid 16S rRNA Gene Evolution in Nonphotosynthetic Green Algae

This study provides a phylogenetic/comparative approach to deciphering the processes underlying the evolution of plastid rRNA genes in genomes under relaxed functional constraints. Nonphotosynthetic green algal taxa that belong to two distinct classes, Chlorophyceae (Polytoma) and Trebouxiophyceae (Prototheca), were investigated. Similar to the situation described previously for plastid 16S rRNA genes in nonphotosynthetic land plants, nucleotide substitution levels, extent of structural variations, and percentage AT values are increased in nonphotosynthetic green algae compared to their closest photosynthetic relatives. However, the mutational processes appear to be different in many respects. First, with the increase in AT content, more transversions are noted in Polytoma and holoparasite angiosperms, while more transitions characterize the evolution of the 16S rDNA sequences in Prototheca. Second, although structural variations do accumulate in both Polytoma and Prototheca (as well as holoparasitic plastid 16S rRNAs), insertions as large as 1.6 kb characterize the plastid 16S rRNA genes in the former, whereas significantly smaller indels (not exceeding 24 bp) seem to be more prevalent in the latter group. The differences in evolutionary rates and patterns within and between lineages might be due to mutations in replication/repair-related genes; slipped-strand mispairing is likely the mechanism responsible for the expansion of insertions in Polytoma plastid 16S rRNA genes. Received: 29 December 2000 / Accepted: 18 May 2001     

Molecular evolution of a portion of the mitochondrial 16S ribosomal gene region in scleractinian corals

Relationships among families and suborders of scleractinian corals are poorly understood because of difficulties 1) in making inferences about the evolution of the morphological characters used in coral taxonomy and 2) in interpreting their 240-million-year fossil record. Here we describe patterns of molecular evolution in a segment of the mitochondrial (mt) 16S ribosomal gene from taxa of 14 families of corals and the use of this gene segment in a phylogenetic analysis of relationships within the order. We show that sequences obtained from scleractinians are homologous to other metazoan 16S ribosomal sequences and fall into two distinct clades defined by size of the amplified gene product. Comparisons of sequences from the two clades demonstrate that both sets of sequences are evolving under similar evolutionary constraints: they do not differ in nucleotide composition, numbers of transition and transversion substitutions, spatial patterns of substitutions, or in rates of divergence. The characteristics and patterns observed in these sequences as well as the secondary structures, are similar to those observed in mt 16S ribosomal DNA sequences from other taxa. Phylogenetic analysis of these sequences shows that they are useful for evaluating relationships within the order. The hypothesis generated from this analysis differs from traditional hypotheses for evolutionary relationships among the Scleractinia and suggests that a reevaluation of evolutionary affinities in the order is needed. Received: 4 September 1996 / Accepted: 7 April 1997     

Phylogenetic Analysis of Reptilian Hemoglobins: Trees, Rates, and Divergences

Phylogenetic relationships among reptiles were examined using previously published and newly determined hemoglobin sequences. Trees reconstructed from these sequences using maximum-parsimony, neighbor-joining, and maximum-likelihood algorithms were compared with a phylogenetic tree of Amniota, which was assembled on the basis of published morphological data. All analyses differentiated α chains into α^A and α^D types, which are present in all reptiles except crocodiles, where only α^A chains are expressed. The occurrence of the α^D chain in squamates (lizards and snakes only in this study) appears to be a general characteristic of these species. Lizards and snakes also express two types of β chains (βI and βII), while only one type of β chain is present in birds and crocodiles. Reconstructed hemoglobin trees for both α and β sequences did not yield the monophyletic Archosauria (i.e., crocodilians + birds) and Lepidosauria (i.e., Sphenodon+ squamates) groups defined by the morphology tree. This discrepancy, as well as some other poorly resolved nodes, might be due to substantial heterogeneity in evolutionary rates among single hemoglobin lineages. Estimation of branch lengths based on uncorrected amino acid substitutions and on distances corrected for multiple substitutions (PAM distances) revealed that relative rates for squamate α^A and α^D chains and crocodilian β chains are at least twice as high as those of the rest of the chains considered. In contrast to these rate inequalities between reptilian orders, little variation was found within squamates, which allowed determination of absolute evolutionary rates for this subset of hemoglobins. Rate estimates for hemoglobins of lizards and snakes yielded 1.7 (α^A) and 3.3 (β) million years/PAM when calibrated with published divergence time vs. PAM distance correlates for several speciation events within snakes and for the squamate ↔ sphenodontid split. This suggests that hemoglobin chains of squamate reptiles evolved ∼3.5 (α^A) or ∼1.7 times (β) faster than their mammalian equivalents. These data also were used to obtain a first estimate of some intrasquamate divergence times. Received: 15 September 1997 / Accepted: 4 February 1998     

Rapid Diversification of Marine Picophytoplankton with Dissimilar Light-Harvesting Structures Inferred from Sequences of Prochlorococcus and Synechococcus (Cyanobacteria)

Cultured isolates of the unicellular planktonic cyanobacteria Prochlorococcus and marine Synechococcus belong to a single marine picophytoplankton clade. Within this clade, two deeply branching lineages of Prochlorococcus, two lineages of marine A Synechococcus and one lineage of marine B Synechococcus exhibit closely spaced divergence points with low bootstrap support. This pattern is consistent with a near-simultaneous diversification of marine lineages with divinyl chlorophyll b and phycobilisomes as photosynthetic antennae. Inferences from 16S ribosomal RNA sequences including data for 18 marine picophytoplankton clade members were congruent with results of psbB and petB and D sequence analyses focusing on five strains of Prochlorococcus and one strain of marine A Synechococcus. Third codon position and intergenic region nucleotide frequencies vary widely among members of the marine picophytoplankton group, suggesting that substitution biases differ among the lineages. Nonetheless, standard phylogenetic methods and newer algorithms insensitive to such biases did not recover different branching patterns within the group, and failed to cluster Prochlorococcus with chloroplasts or other chlorophyll b-containing prokaryotes. Prochlorococcus isolated from surface waters of stratified, oligotrophic ocean provinces predominate in a lineage exhibiting low G + C nucleotide frequencies at highly variable positions. Received: 18 January 1997 / Accepted: 18 May 1997     

Evolution of the nad3-rps12 Gene Cluster in Angiosperm Mitochondria: Comparison of Edited and Unedited Sequences

We have analyzed the nad3-rps12 locus for eight angiosperms in order to compare the utility of mitochondrial DNA and edited mRNA sequences in phylogenetic reconstruction. The two coding regions, containing from 25 to 35 editing sites in the various plants, have been concatenated in order to increase the significance of the analysis. Differing from the corresponding chloroplast sequences, unedited mitochondrial DNA sequences seem to evolve under a quasi-neutral substitution process which undifferentiates the nucleotide substitution rates for the three codon positions. By using complete gene sequences (all codon positions) we found that genomic sequences provide a classical angiosperm phylogenetic tree with a clear-cut grouping of monocotyledons and dicotyledons with Magnoliidae at the basal branch of the tree. Conversely, owing to their low nucleotide substitution rates, edited mRNA sequences were found not to be suitable for studying phylogenetic relationships among angiosperms. Received: 24 January 1996 / Accepted: 5 June 1996     

Analysis of Ribosomal RNA Genes Suggests That Trypanosomes Are Monophyletic

To further investigate the phylogeny of protozoa from the order Kinetoplastida we have sequenced the small subunit (SSU) and a portion of the large subunit (LSU) nuclear rRNA genes. The SSU and LSU sequences were determined from a lizard trypanosome, Trypanosoma scelopori and a bodonid, Rhynchobodo sp., and the LSU sequences were determined from an insect trypanosomatid, Crithidia oncopelti, and a bodonid, Dimastigella trypaniformis. Contrary to previous results, in which trypanosomes were found to be paraphyletic, with Trypanosoma brucei representing the earliest-diverging lineage, we have now found evidence for the monophyly of trypanosomes. Addition of new taxa which subdivide long branches (such as that of T. brucei) have helped to identify homoplasies responsible for the paraphyletic trees in previous studies. Although the monophyly of the trypanosome clade is supported in the bootstrap analyses for maximum likelihood at 97% and maximum parsimony at 92%, there is only a small difference in ln-likelihood value or tree length between the most optimal monophyletic tree and the best suboptimal paraphyletic tree. Within the trypanosomatid subtree, the clade of trypanosomes is a sister group to the monophyletic clade of the nontrypanosome genera. Different groups of trypanosomes group on the tree according to their mode of transmission. This suggests that the adaptation to invertebrate vectors plays a more important role in the trypanosome evolution than the adaptation to vertebrate hosts. Received: 5 July 1996 / Accepted: 26 September 1996     

Erratic Evolution of Glycerol-3-Phosphate Dehydrogenase in Drosophila,Chymomyza, and Ceratitis

We have studied the evolution of Gpdh in 18 fruitfly species by sequencing 1,077 nucleotides per species on average. The region sequenced includes four exons coding for 277 amino acids and three variable-length introns. Phylogenies derived by a variety of methods confirm that the nominal genus Zaprionus belongs within the genus Drosophila, whereas Scaptodrosophila and Chymomyza are outside. The rate of GPDH evolution is erratic. The rate of amino acid replacements in a lineage appears to be 1.0 × 10⁻¹⁰/site/year when Drosophila species are considered (diverged up to 55 million years ago), but becomes 2.3 × 10⁻¹⁰ when they are compared to Chymomyza species (divergence around 60 My ago), and 4.6 × 10⁻¹⁰ when species of those two genera are compared with the medfly Ceratitis capitata (divergence around 100 My ago). In order to account for these observations, the rate of amino acid replacement must have been 15 or more times greater in some lineages and at some times than in others. At the nucleotide level, however, Gpdh evolves in a fairly clockwise fashion. Received: 13 June 1996 / Accepted: 16 August 1996     

Molecular evidence for an independent origin of modern triserial planktonic foraminifera from benthic ancestors

Gallitellia vivans is the only Recent representative of the triserial planktonic foraminiferal family Guembelitriidae. The origin and evolution of this interesting albeit poorly known family are enigmatic. To elucidate the phylogenetic relationships between G. vivans and other planktonic foraminifera, we sequenced the small subunit ribosomal DNA (SSU rDNA) for comparison to our extensive database of planktonic and benthic species. Our analyses suggest that G. vivans represents a separate lineage of planktonic foraminifera, which branches close to the benthic rotaliids Stainforthia and Virgulinella. Both genera resemble Gallitellia in general morphological appearance, having elongate triserial tests at least in their early ontogenic stages. The divergence time of G. vivans is estimated at ca. 18 Ma (early Miocene), suggesting an origin independent from the Cretaceous and Paleogene triserial planktonic foraminifera. Our study thus indicates that modern triserial planktonic foraminifera are not related to the Cretaceous–Paleogene triserial species, and that the sporadic occurrences in the fossil record are not the result of poor preservation, but reflect multiple transitions from benthic to planktonic mode of life.