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     

Extensive Intraindividual Variation in Plastid rDNA Sequences from the Holoparasite <Emphasis Type="Italic">Cynomorium coccineum</Emphasis> (Cynomoriaceae)

Ribosomal genes are considered to have a high degree of sequence conservation between species and also at higher taxonomic levels. In this paper we document a case where a single individual of Cynomorium coccineum (Cynomoriaceae), a nonphotosynthetic holoparasitic plant, contains highly divergent plastid ribosomal genes. PCR amplification a nearly complete ribosomal DNA cistron was performed using genomic DNA, the products cloned, and the 23S rDNA genes were sequenced from 19 colonies. Of these, five distinct types were identified. Fifteen of the sequences were nearly identical (11 or fewer differences) and these were designated Type I. The remaining types (II–V) were each represented by a single clone and differed from Type I by 93 to 255 changes. Compared with green vascular plants, we found that there are more substitutional differences in the 23S rDNA sequences within a single individual of Cynomorium than among all sequenced photosynthetic vascular plants. Several trends of molecular evolution observed in 16S rDNA from other holoparasitic angiosperms and heterotrophic green algae have been also observed in Cynomorium 23S rDNA. Higher-order structures were constructed for representatives of the five clone types, and in all cases these possessed complete complements of the major structural elements present in functional plastid 23S rRNAs. These data indicate that such molecules may be subject to purifying selection, thus providing indirect evidence that they have retained some degree of functionality. This intraindividual polymorphism is probably a case of plastid heteroplasmy but translocation of ribosomal cistrons to the nucleus or mitochondria has not been tested and therefore cannot be ruled out.This article contains online supplementary material.     

The effect of relaxed functional constraints on the photosynthetic gene rbcL in photosynthetic and nonphotosynthetic parasitic plants

The photosynthetic gene rbcL has been lost or dramatically altered in some lineages of nonphotosynthetic parasitic plants, but the dynamics of these events following loss of photosynthesis and whether rbcL has sustained functionally significant changes in photosynthetic parasitic plants are unknown. To assess the changes to rbcL associated with the loss of functional constraints for photosynthesis, nucleotide sequences from nonparasitic and parasitic plants of Scrophulariales were used for phylogeny reconstruction and character analysis. Plants in this group display a broad range of parasitic abilities, from photosynthetic ("hemiparasites") to nonphotosynthetic ("holoparasites"). With the exception of Conopholis (Orobanchaceae), the rbcL locus is present in all parasitic plants of Scrophulariales examined. Several holoparasitic genera included in this study, including Boschniakia, Epifagus, Orobanche, and Hyobanche, have rbcL pseudogenes. However, the holoparasites Alectra orobanchoides, Harveya capensis, Harveya purpurea, Lathraea clandestina, Orobanche corymbosa, O. fasciculata, and Striga gesnerioides have intact open reading frames (ORFs) for the rbcL gene. Phylogenetic hypotheses based on rbcL are largely in agreement with those based on sequences of the nonphotosynthetic genes rps2 and matK and show a single origin of parasitism, and loss of photosynthesis and pseudogene formation have been independently derived several times in Scrophulariales. The mutations in rbcL in nonparasitic and hemiparasitic plants would result in largely conservative amino acid substitutions, supporting the hypothesis that functional proteins can experience only a limited range of changes, even in minimally photosynthetic plants. In contrast, ORFs in some holoparasites had many previously unobserved missense substitutions at functionally important amino acid residues, suggesting that rbcL genes in these plants have evolved under relaxed or altered functional constraints.      

Rate accelerations in nuclear 18S rDNA of mycoheterotrophic and parasitic angiosperms

Rate variation in genes from all three genomes has been observed frequently in plant lineages with a parasitic and mycoheterotrophic mode of life. While the loss of photosynthetic ability leads to a relaxation of evolutionary constraints in genes involved in the photosynthetic apparatus, it remains to be determined how prevalent increased substitution rates are in nuclear DNA of non-photosynthetic angiosperms. In this study we infer rates of molecular evolution of 18S rDNA of all parasitic and mycoheterotorphic plant families (except Lauraceae and Polygalaceae) using relative rate tests. In several holoparasitic and mycoheterotrophic plant lineages extremely high substitution rates are observed compared to other photosynthetic angiosperms. The position and frequency of these substitutions have been identified to understand the mutation dynamics of 18S rRNA in achlorophyllous plants. Despite the presence of significantly elevated substitution rates, very few mutations occur in major functional and structural regions of the small ribosomal molecule, providing evidence that the efficiency of the translational apparatus in non-photosynthetic plants has not been affected.     

Structural analyses of plastid-derived 16S rRNAs in holoparasitic angiosperms

Higher-order structures have been constructed for plastid-encoded small-subunit (SSU, 16S), rRNAs from representatives of seven nonphotosynthetic holoparasitic angiosperm families: Apodanthaceae, Cynomoriaceae, Cytinaceae, Balanophoraceae, Hydnoraceae, Mitrastemonaceae, and Rafflesiaceae. Whereas most pairwise comparisons among angiosperms differ by 2–3% in substitutions, the 16S rRNAs of the holoparasites show an increasingly greater number of mutations: Cynomorium (7.3%), Cytinus (8.0%), Bdallophyton (12.7%), Mitrastema (14.9%), Hydnora (19.4%), Pilostyles (30.4%) and Corynaea (35.9%). Despite this high level of sequence variation, SSU structures constructed for all species except Pilostyles possess the typical complement of 50 helices (that contain numerous compensatory mutations) thereby providing indirect evidence supporting their functionality. Pilostyles, likely with the most unusual plastid 16S rRNA yet documented, lacks four major helices and contains lengthy insertions for four others. Sequences of products generated via RT-PCR show that these structural modifications are present on a mature (transcribed) rRNA. The trend toward increasing numbers of base substitutions in the holoparasites is accompanied by a marked increase in AA+U content of the rRNA. This A/T drift phenomenon of rDNA is especially apparent in Corynaea whose SSU rDNA sequence is 72% A+T. A comparison of Cytinus to tobacco showed that substitution rates appear to be dependent upon the composition of neighboring bases. Transversions represented 26% of the mutations when flanking bases were G or C whereas transversions increased to 36% when the flanking bases were A to T. The underlying molecular mechanism associated with these high substitution rates is presently unknown, however, relaxation of selection pressure on ribosome function resulting in altered DNA replication and/or repair systems may be involved.     

Origin and Evolution of Paralogous rRNA Gene Clusters Within the Flatworm Family Dugesiidae (Platyhelminthes, Tricladida)

Analysis of the 18S rDNA sequences of five species of the family Dugesiidae (phylum Platyhelminthes, suborder Tricladida, infraorder Paludicola) and eight species belonging to families Dendrocoelidae and Planaridae and to the infraorder Maricola showed that members of the family Dugesiidae have two types of 18S rDNA genes, while the rest of the species have only one. The duplication event also affected the ITS-1, 5.8S, ITS-2 region and probably the 28S gene. The mean divergence value between the type I and the type II sequences is 9% and type II 18S rDNA genes are evolving 2.3 times more rapidly than type I. The evolutionary rates of type I and type II genes were calibrated from biogeographical data, and an approximate date for the duplication event of 80–120 million years ago was calculated. The type II gene was shown, by RT-PCR, to be transcribed in adult individuals of Schmidtea polychroa, though at very low levels. This result, together with the fact that most of the functionally important positions for small-subunit rRNA in prokaryotes have been conserved, indicates that the type II gene is probably functional. Received: 24 March 1998 / Accepted: 17 March 1999     

Analysis of the Primary Sequence and Secondary Structure of the Unusually Long SSU rRNA of the Soil Bug, Armadillidium vulgare

The complete nucleotide sequence of the SSU rRNA gene from the soil bug, Armadillidium vulgare (Crustacea, Isopoda), was determined. It is 3214 bp long, with a GC content of 56.3%. It is not only the longest SSU rRNA gene among Crustacea but also longer than any other SSU rRNA gene except that of the strepsipteran insect, Xenos vesparum (3316 bp). The unusually long sequence of this species is explained by the long sequences of variable regions V4 and V7, which make up more than half of the total length. RT-PCR analysis of these two regions showed that the long sequences also exist in the mature rRNA and sequence simplicity analysis revealed the presence of slippage motifs in these two regions. The putative secondary structure of the rRNA is typical for eukaryotes except for the length and shape variations of the V2, V4, V7, and V9 regions. Each of the V2, V4, and V7 regions was elongated, while the V9 region was shortened. In V2, two bulges, located between helix 8 and helix 9 and between helix 9 and helix 10, were elongated. In V4, stem E23-3 was dramatically expanded, with several small branched stems. In V7, stem 43 was branched and expanded. Comparisons with the unusually long SSU rRNAs of other organisms imply that the increase in total length of SSU rRNA is due mainly to expansion in the V4 and V7 regions. Received: 2 March 1999 / Accepted: 22 July 1999     

Equally Parsimonious Pathways Through an RNA Sequence Space Are Not Equally Likely

An experimental system for determining the potential ability of sequences resembling 5S ribosomal RNA (rRNA) to perform as functional 5S rRNAs in vivo in the Escherichia coli cellular environment was devised previously. Presumably, the only 5S rRNA sequences that would have been fixed by ancestral populations are ones that were functionally valid, and hence the actual historical paths taken through RNA sequence space during 5S rRNA evolution would have most likely utilized valid sequences. Herein, we examine the potential validity of all sequence intermediates along alternative equally parsimonious trajectories through RNA sequence space which connect two pairs of sequences that had previously been shown to behave as valid 5S rRNAs in E. coli. The first trajectory requires a total of four changes. The 14 sequence intermediates provide 24 apparently equally parsimonious paths by which the transition could occur. The second trajectory involves three changes, six intermediate sequences, and six potentially equally parsimonious paths. In total, only eight of the 20 sequence intermediates were found to be clearly invalid. As a consequence of the position of these invalid intermediates in the sequence space, seven of the 30 possible paths consisted of exclusively valid sequences. In several cases, the apparent validity/invalidity of the intermediate sequences could not be anticipated on the basis of current knowledge of the 5S rRNA structure. This suggests that the interdependencies in RNA sequence space may be more complex than currently appreciated. If ancestral sequences predicted by parsimony are to be regarded as actual historical sequences, then the present results would suggest that they should also satisfy a validity requirement and that, in at least limited cases, this conjecture can be tested experimentally. Received: 27 August 1996 / Accepted: 14 April 1997     

Do nonasterid holoparasitic flowering plants have plastid genomes?

Past work involving the plastid genome (plastome) of holoparasitic plants has been confined to Scrophulariaceae (or Orobanchaceae) which have truncated plastomes owing to loss of photosynthetic and other genes. Nonasterid holoparasites from Balanophoraceae (Corynaea), Hydnoraceae (Hydnora) and Cytinaceae (Cytinus) were tested for the presence of plastid genes and a plastome. Using PCR, plastid 16S rDNA was successfully amplified and sequenced from the above three holoparasites. The sequence of Cytinus showed 121 single base substitutions relative to Nicotiana (8% of the molecule) whereas higher sequence divergence was observed in Hydnora and Corynaea (287 and 513 changes, respectively). Secondary structural models for these 16S rRNAs show that most changes are compensatory, thus suggesting they are functional. Probes constructed for 16S rDNA and for four plastid-encoded ribosomal protein genes (rps2, rps4, rps7 and rpl16) were used in Southern blots of digested genomic DNA from the three holoparasites. Positive hybridizations were obtained using each of the five probes only for Cytinus. For SmaI digests, all plastid gene probes hybridized to a common fragment ca. 20 kb in length in this species. Taken together, these data provide preliminary evidence suggestive of the retention of highly diverged and truncated plastid genome in Cytinus. The greater sequence divergence for 16S rDNA and the negative hybridization results for Hydnora and Corynaea suggests two possibilities: the loss of typically conserved elements of their plastomes or the complete absence of a plastome.     

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     

High rates of nucleotide substitution in nuclear small-subunit (18S) rDNA from holoparasitic flowering plants

Relative rate tests, using Gnetum as a reference taxon, were conducted on nuclear 18S rRNA sequences from 10 angiosperms including autotrophic nonparasites (Arabidopsis, Asarum, Glycine, Malpighia, and Zea), a chlorophyllous hemiparasite (Arceuthobium—Viscaceae), and achlorophyllous holoparasites (Balanophora—Balanophoraceae, Prosopanche—Hydnoraceae, and Rafflesia and Rhizanthes—Rafflesiaceae). Compared with Glycine, the mean number of substitutions per site (K) for five autotrophic angiosperms is 0.036 whereas for the holoparasites K = 0.126, i.e., 3.5 times higher. Comparisons of autotrophic species with short and long generation times showed no differences in K; hence, divergent rRNA sequences in the holoparasites are likely attributable to other mechanisms. These might include genetic bottlenecks, effective population size, and/or molecular drive. High substitution rates appear to be associated only with those parasitic angiosperms that have developed a highly modified haustorial system and extreme nutritional dependence upon the host. At present, high substitution rates in these parasites confound attempts to determine their phylogenetic position relative to other angiosperms. Correspondence to: D.L. Nickrent     

A Group I Intron in the 18S Ribosomal DNA from the Parasitic Fungus Isaria japonica

The nucleotide sequence of the 18S rDNA coding gene in the ascomycetes parasitic fungus Isaria japonica contains a group I intron with a length of 379 nucleotides. The identification of the DNA sequence as a group I intron is based on its sequence homology to other fungal group I introns. Its group I intron contained the highly conserved sequence elements P, Q, R, and S found in other group I introns. Surprisingly, the intron sequence of I. japonica is more similar to that of Ustilago maydis than to the one found in Sclerotinia sclerotiorum. This is in contrast to the sequence identity found on the neighboring rDNA. This is an interesting finding and suggests a horizontal transfer of group I intron sequences. Received: 19 September 1997 / Accepted: 10 September 1998     

Phylogenetic Position of the Phacotaceae Within the Chlamydophyceae as Revealed by Analysis of 18S rDNA and rbcL Sequences

Four genera of the Phacotaceae (Phacotus, Pteromonas, Wislouchiella, Dysmorphococcus), a family of loricated green algal flagellates within the Volvocales, were investigated by means of transmission electron microscopy and analysis of the nuclear encoded small-subunit ribosomal RNA (18S rRNA) genes and the plastid-encoded rbcL genes. Additionally, the 18S rDNA of Haematococcus pluvialis and the rbcL sequences of Chlorogonium elongatum, C. euchlorum, Dunaliella parva, Chloromonas serbinowii, Chlamydomonas radiata, and C. tetragama were determined. Analysis of ultrastructural data justified the separation of the Phacotaceae into two groups. Phacotus, Pteromonas, and Wislouchiella generally shared the following characters: egg-shaped protoplasts, a single pyrenoid with planar thylakoid double-lamellae, three-layered lorica, flagellar channels as part of the central lorica layer, mitochondria located in the central cytoplasm, lorica development that occurs in mucilaginous zoosporangia that are to be lysed, and no acid-resistant cell walls. Dysmorphococcus was clearly different in each of the characters mentioned. Direct comparison of sequences of Phacotus lenticularis, Pteromonas sp., Pteromonas protracta, and Wislouchiella planctonica revealed DNA sequence homologies of ≥98.0% within the 18S gene and 93.9% within the rbcL gene. D. globosus was quite different from these species, with a maximum of 92.9% homology in the 18S rRNA and ≤86.6% in the rbcL gene. It showed major similarities to the 18S rDNA of Dunaliella salina, with 95.3%, and to the rbcL sequence of Chlamydomonas tetragama, with 90.3% sequence homology. Additionally, the Phacotaceae sensu stricto exclusively shared 10 (rbcL: 4) characters which were present neither in other Chlamydomonadales nor in Dysmorphococcus globosus. Different phylogenetic analysis methods confirmed the hypothesis that the Phacotaceae are polyphyletic. The Phacotaceae sensu stricto form a stable cluster with affinities to the Dunaliellaes and possibly Haematococcus pluvialis. Dysmorphococcus globosus represented an independent lineage that is possibly related to Chlamydomonas moewusii and C. tetragama. Received: 9 June 1997 / Accepted: 17 October 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     

Phylogeny of Ultra-Rapidly Evolving Dinoflagellate Chloroplast Genes: A Possible Common Origin for Sporozoan and Dinoflagellate Plastids

Complete chloroplast 23S rRNA and psbA genes from five peridinin-containing dinoflagellates (Heterocapsa pygmaea, Heterocapsa niei, Heterocapsa rotun-data, Amphidinium carterae, and Protoceratium reticulatum) were amplified by PCR and sequenced; partial sequences were obtained from Thoracosphaera heimii and Scrippsiella trochoidea. Comparison with chloroplast 23S rRNA and psbA genes of other organisms shows that dinoflagellate chloroplast genes are the most divergent and rapidly evolving of all. Quartet puzzling, maximum likelihood, maximum parsimony, neighbor joining, and LogDet trees were constructed. Intersite rate variation and invariant sites were allowed for with quartet puzzling and neighbor joining. All psbA and 23S rRNA trees showed peridinin-containing dinoflagellate chloroplasts as monophyletic. In psbA trees they are related to those of chromists and red algae. In 23S rRNA trees, dinoflagellates are always the sisters of Sporozoa (apicomplexans); maximum likelihood analysis of Heterocapsa triquetra 16S rRNA also groups the dinoflagellate and sporozoan sequences, but the other methods were inconsistent. Thus, dinoflagellate chloroplasts may actually be related to sporozoan plastids, but the possibility of reproducible long-branch artifacts cannot be strongly ruled out. The results for all three genes fit the idea that dinoflagellate chloroplasts originated from red algae by a secondary endosymbiosis, possibly the same one as for chromists and Sporozoa. The marked disagreement between 16S rRNA trees using different phylogenetic algorithms indicates that this is a rather poor molecule for elucidating overall chloroplast phylogeny. We discuss possible reasons why both plastid and mitochondrial genomes of alveolates (Dinozoa, Sporozoa and Ciliophora) have ultra-rapid substitution rates and a proneness to unique genomic rearrangements. Received: 27 December 1999 / Accepted: 24 March 2000     

Horizontal transfer of genes coding for the photosynthetic reaction centers of purple bacteria

Phylogenetic trees were drawn and analyzed based on the nucleotide sequences of the 1.5-kb gene fragment coding for the L and M subunits of the photochemical reaction center of various purple photosynthetic bacteria. These trees are mostly consistent with phylogenetic trees based on 16S rRNA and soluble cytochrome c, but differ in some significant details. This inconsistency implies horizontal transfer of the genes that code for the photosynthetic apparatus in purple bacteria. Possibilities of similar transfers of photosynthesis genes during the evolution of photosynthesis are discussed especially for the establishment of oxygenic photosynthesis. Received: 8 July 1996 / Accepted: 12 March 1997     

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     

Anoxygenic phototrophy across the phylogenetic spectrum: current understanding and future perspectives

The phylogenetic heterogeneity of anoxygenic phototrophic bacteria has been revealed by 16S rRNA sequence analysis, the results of which have led to extensive taxonomic rearrangements within previously defined taxa of phototrophs and stimulated interest in this group of organisms. Anoxygenic photosynthetic bacteria can be found within 4 of the 12 phylogenetic lineages, and in some cases are highly related to non-photosynthetic members of these groups. The largest number of phototrophs are found in the class Proteobacteria. Comparative phylogenetic analysis using 23S rDNA sequences generally supports the topology obtained from 16S rDNA sequences. The photosynthetic reaction centers are conserved in all photosynthetic bacteria, and are of two types. One is shared by the Proteobacteria and Chloroflexus aurantiacus and is similar to Photosystem II of cyanobacteria, while heliobacteria and Chlorobium and relatives possess a reaction center similar to the cyanobacterial Photosystem I. These similarities are supported by sequence analysis of core reaction center peptides, but contradict phylogenies reconstructed from rRNA sequence analysis. Genome analysis by means of physical mapping has been performed for only three species of anoxygenic phototrophs. Some conservation of operon structure and gene sequence has been found within the Proteobacteria, but does not extend to other phototrophs. Received: 29 December 1995 / Accepted: 19 July 1996