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     

Mechanisms of Functional and Physical Genome Reduction in Photosynthetic and Nonphotosynthetic Parasitic Plants of the Broomrape Family

Nonphotosynthetic plants possess strongly reconfigured plastomes attributable to convergent losses of photosynthesis and housekeeping genes, making them excellent systems for studying genome evolution under relaxed selective pressures. We report the complete plastomes of 10 photosynthetic and nonphotosynthetic parasites plus their nonparasitic sister from the broomrape family (Orobanchaceae). By reconstructing the history of gene losses and genome reconfigurations, we find that the establishment of obligate parasitism triggers the relaxation of selective constraints. Partly because of independent losses of one inverted repeat region, Orobanchaceae plastomes vary 3.5-fold in size, with 45 kb in American squawroot (Conopholis americana) representing the smallest plastome reported from land plants. Of the 42 to 74 retained unique genes, only 16 protein genes, 15 tRNAs, and four rRNAs are commonly found. Several holoparasites retain ATP synthase genes with intact open reading frames, suggesting a prolonged function in these plants. The loss of photosynthesis alters the chromosomal architecture in that recombinogenic factors accumulate, fostering large-scale chromosomal rearrangements as functional reduction proceeds. The retention of DNA fragments is strongly influenced by both their proximity to genes under selection and the co-occurrence with those in operons, indicating complex constraints beyond gene function that determine the evolutionary survival time of plastid regions in nonphotosynthetic plants.     

The complete chloroplast genome of the chlorarachniophyte Bigelowiella natans: evidence for independent origins of chlorarachniophyte and euglenid secondary endosymbionts

Chlorarachniophytes are amoeboflagellate cercozoans that acquired a plastid by secondary endosymbiosis. Chlorarachniophytes are the last major group of algae for which there is no completely sequenced plastid genome. Here we describe the 69.2-kbp chloroplast genome of the model chlorarachniophyte Bigelowiella natans. The genome is highly reduced in size compared with plastids of other photosynthetic algae and is closer in size to genomes of several nonphotosynthetic plastids. Unlike nonphotosynthetic plastids, however, the B. natans chloroplast genome has not sustained a massive loss of genes, and it retains nearly all of the functional photosynthesis-related genes represented in the genomes of other green algae. Instead, the genome is highly compacted and gene dense. The genes are organized with a strong strand bias, and several unusual rearrangements and inversions also characterize the genome; notably, an inversion in the small-subunit rRNA gene, a translocation of 3 genes in the major ribosomal protein operon, and the fragmentation of the cluster encoding the large photosystem proteins PsaA and PsaB. The chloroplast endosymbiont is known to be a green alga, but its evolutionary origin and relationship to other primary and secondary green plastids has been much debated. A recent hypothesis proposes that the endosymbionts of chlorarachniophytes and euglenids share a common origin (the Cabozoa hypothesis). We inferred phylogenies using individual and concatenated gene sequences for all genes in the genome. Concatenated gene phylogenies show a relationship between the B. natans plastid and the ulvophyte-trebouxiophyte-chlorophyte clade of green algae to the exclusion of Euglena. The B. natans plastid is thus not closely related to that of Euglena, which suggests that plastids originated independently in these 2 groups and the Cabozoa hypothesis is false.     

Functional gene losses occur with minimal size reduction in the plastid genome of the parasitic liverwort Aneura mirabilis

Aneura mirabilis is a parasitic liverwort that exploits an existing mycorrhizal association between a basidiomycete and a host tree. This unusual liverwort is the only known parasitic seedless land plant with a completely nonphotosynthetic life history. The complete plastid genome of A. mirabilis was sequenced to examine the effect of its nonphotosynthetic life history on plastid genome content. Using a partial genomic fosmid library approach, the genome was sequenced and shown to be 108,007 bp with a structure typical of green plant plastids. Comparisons were made with the plastid genome of Marchantia polymorpha, the only other liverwort plastid sequence available. All ndh genes are either absent or pseudogenes. Five of 15 psb genes are pseudogenes, as are 2 of 6 psa genes and 2 of 6 pet genes. Pseudogenes of cysA, cysT, ccsA, and ycf3 were also detected. The remaining complement of genes present in M. polymorpha is present in the plastid of A. mirabilis with intact open reading frames. All pseudogenes and gene losses co-occur with losses detected in the plastid of the parasitic angiosperm Epifagus virginiana, though the latter has functional gene losses not found in A. mirabilis. The plastid genome sequence of A. mirabilis represents only the second liverwort, and first mycoheterotroph, to have its plastid genome sequenced. We observed a pattern of genome evolution congruent with functional gene losses in parasitic angiosperms but suggest that its plastid genome represents a genome in the early stages of decay following the relaxation of selection pressures.     

Nucleus-encoded,plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH) indicates a single origin for chromalveolate plastids

Plastids (the photosynthetic organelles of plants and algae) originated through endosymbiosis between a cyanobacterium and a eukaryote and subsequently spread to other eukaryotes by secondary endosymbioses between two eukaryotes. Mounting evidence favors a single origin for plastids of apicomplexans, cryptophytes, dinoflagellates, haptophytes, and heterokonts (together with their nonphotosynthetic relatives, termed chromalveolates), but so far, no single molecular marker has been described that supports this common origin. One piece of evidence comes from plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which originated by a gene duplication of the cytosolic form. However, no plastid GAPDH has been characterized from haptophytes, leaving an important piece of the puzzle missing. We have sequenced genes encoding cytosolic, mitochondrion-targeted, and plastid-targeted GAPDH proteins from a number of haptophytes and heterokonts and found haptophyte homologs that branch within a strongly supported clade of chromalveolate plastid-targeted genes, being more closely related to an apicomplexan homolog than was expected. The evolution of plastid-targeted GAPDH supports red algal ancestry of apicomplexan plastids and raises a number of questions about the importance of plastid loss and the possibility of cryptic plastids in nonphotosynthetic lineages such as ciliates.     

EVOLUTIONARY CONSEQUENCES OF THE LOSS OF PHOTOSYNTHESIS IN CHLAMYDOMONADACEAE: PHYLOGENETIC ANALYSIS OF Rrn18 (18S rDNA) IN 13 POLYTOMA STRAINS (CHLOROPHYTA)1

Complete sequences of the Rrn 18 genes were obtained from 13 strains of the nonphotosynthetic algal genus Polytoma. Phylogenetic analyses showed that these strains formed two clades. One clade shows only modest sequence diversity but is represented by strains collected at widely dispersed sites in Europe and America. The other clade consists of a single isolate from the Canary Islands. Both clades lie well within the extended clade that includes all species of Chlamydomonas for which sequence data are available. The two Polytoma clades are separated from each other by several green species, suggesting that the extant nonphotosynthetic Chlamydomonadaceae arose from photosynthetic ancestors at least twice. These results suggest that nonphotosynthetic mutants are capable of establishing lineages that can spread widely but have a higher probability of extinction than their photosynthetic congeners.     

Functional EF-Tu with large C-terminal extensions in an E coli strain with a precise deletion of both chromosomal tuf genes

An Escherichia coli strain was constructed in which both chromosomal genes encoding elongation factor (EF)-Tu (tufA and tufB) have been inactivated with precise coding sequence replacements. A tufA gene in an expression vector is supplied as the sole EF-Tu source. By using plasmid replacement, based on plasmid incompatibility, mutant EF-Tu variants with a large C'-terminal extension up to 270 amino acids were studied and proved to be functional in a strain lacking the chromosomal tufA and tufB genes.     

Characterization of cytoplasmic and nuclear genomes in the colorless alga Polytoma. I. Ultrastructural analysis of organelles

Electron microscope studies have been made on the fine structure of the colorless biflagellate, Polytoma obtusum, with main emphasis on the structural organization of the mitochondria and the leucoplast. Both organelles have been demonstrated to contain DNA aggregates as well as ribosomal particles within their matrix material. Reconstructions from serial showed that (a) the mitochondria were highly convoluted and irregular in shape and size, and (b) the leucoplast was a single cup-shaped entity, with large starch grains, localized at the posterior end, and multiple sites of DNA aggregates. The starch-containing compartments appeared to be interconnected by narrow tubular or sheetlike bridges. Cytoplasmic invaginations into the plastid, often containing mitochondria, were of frequent occurrence, and membranes of mitochondria and the leucoplast appeared to be closely apposed. Membranes elements, both sheetlike and vesicular, were also present in the matrix. The Polytoma leucoplast was, in certain respects, morphologically similar to the plastids of various photosynthetic mutants of Chlamydomonas, most of which show Menedelian segregation. It is suggested that Polytoma arose from a Chlamydomonas-like ancestor, possibly through combined mutational processess of both chloroplast and nuclear genomes. Since Polytoma leucoplasts contain both DNA and ribosomal particles, it is probable that these organelles still possess semiautonomy and limited ability for protein synthesis.     

A subset of conserved tRNA genes in plastid DNA of nongreen plants.

The plastid genome of the nonphotosynthetic parasitic plant Epifagus virginiana contains only 17 of the 30 tRNA genes normally found in angiosperm plastid DNA. Although this is insufficient for translation, the genome is functional, so import of cytosolic tRNAs into plastids has been suggested. This raises the question of whether the tRNA genes that remain in E. virginiana plastid DNA are active or have just fortuitously escaped deletion. We report the sequences of 20 plastid tRNA loci from Orobanche minor, which shares a nonphotosynthetic ancestor with E. virginiana. The two species have 9 intact tRNA genes in common, the others being defunct in one or both species. The intron-containing trnLUAA gene is absent from E. virginiana, but it is intact, transcribed, and spliced in O. minor. The shared intact genes are better conserved than intergenic sequences, which indicates that these genes are being maintained by natural selection and, therefore, must be functional. For the most part, the tRNA species conserved in nonphotosynthetic plastids are also those that have never been found to be imported in plant mitochondria, which suggests that the same rules may govern tRNA import in the two organelles. A small photosynthesis gene, psbI, is still intact in O. minor, and computer simulations show that some small nonessential genes have an appreciable chance of escaping deletion.     

Rapid evolution of the plastid translational apparatus in a nonphotosynthetic plant: Loss or accelerated sequence evolution of tRNA and ribosomal protein genes

Summary The vestigial plastid genome of Epifagus virginiana (beechdrops), a nonphotosynthetic parasitic flowering plant, is functional but lacks six ribosomal protein and 13 tRNA genes found in the chloroplast DNAs of photosynthetic flowering plants. Import of nuclear gene products is hypothesized to compensate for many of these losses. Codon usage and amino acid usage patterns in Epifagus plastic genes have not been affected by the tRNA gene losses, though a small shift in the base composition of the whole genome (toward A + T -richness) is apparent. The ribosomal protein and tRNA genes that remain have had a high rate of molecular evolution, perhaps due to relaxation of constraints on the translational apparatus. Despite the compactness and extensive gene loss, one translational gene (infA, encoding initiation factor 1) that is a pseudogene in tobacco has been maintained intact in Epifagus.Offprint requests to: J.D. Palmer     

Distribution and Evolution of Pseudogenes, Gene Losses, and a Gene Rearrangement in the Plastid Genome of the Nonphotosynthetic Liverwort, Aneura mirabilis (Metzgeriales, Jungermanniopsida)

The plastid genome sequence of the parasitic liverwort Aneura mirabilis revealed the loss of five chlororespiration (ndh) genes. Additionally, six ndh genes, subunits of photosystem I, photosystem II, and the cytochrome b6f complex were inferred to be pseudogenes. Pseudogenes of cysA, cyst, ccsA, and ycf3, an inversion of psbE and petL, were also detected. The designation of pseudogenes was made using comparisons with the distantly related liverwort Marchantia polymorpha. We sampled several populations of A. mirabilis and its photosynthetic sister groups to correlate functional gene losses with the evolution of a achlorophylly. The gene losses, pseudogenes, or the psbE-petL inversion were never detected in a photosynthetic Aneura but were detected in every population of A. mirabilis. One population of A. mirabilis revealed a unique deletion of 541 bp in the psbE-petL region; another is characterized by a unique deletion of 471 bp in the trnV(UAC)-ndhC region. The ratio of synonymous-to-nonsynonymous substitution rates (omega) was estimated for eight pseudogenes and six ORFs to detect relaxed purifying selection. A significant increase in omega for the nonphotosynthetic liverwort was detected in six pseudogenes. Relaxation purifying selection, determined by a significant increase in omega, was detected for three intact ORFs: psbA, psbM, and rbcL.     

Comparison of the complete sequence of the str operon in Salmonella typhimurium and Escherichia coli.

The nucleotide (nt) sequences of the str operon in Escherichia coli K-12 and Salmonella typhimurium LT2 were completed and compared at the nt and amino acid (aa) level. The order of conservation at the nt and aa level is rpsL greater than tufA greater than rpsG greater than f usA. A striking difference is that the rpsG-encoded ribosomal protein, S7, in E. coli K-12 is 23 aa longer than in S. typhimurium. The very low (0.18) codon adaptation index of this part of the E. coli K-12-encoding gene and the unusual stop codon (UGA) suggest that this is a relatively recent extension. A trend towards a higher G+C content in fusA (gene encoding elongation factor (EF)-G) and tufA (gene encoding EF-Tu) in S. typhimurium is noted. In fusA, nt substitutions at all three positions in a codon occur at a much higher frequency than expected from the number of nt substitutions in the gene, assuming they are random and independent events. An analysis of substitutions in this and other genes suggests that the triple substitutions in fusA, and some other genes, are the result of the sequential accumulation of individual mutations, probably driven by selection pressure for particular codons or aa.     

Characterization of cytoplasmic and nuclear genomes in the colorless alga Polytoma. III. Ribosomal RNA cistrons of the nucleus and leucoplast

The colorless alga Polytoma obtusum has been found to possess leucoplasts, and two kinds of ribosomes with sedimentation values of 73S and 79S. The ribosomal RNA (rRNA) of the 73S but not the 79S ribosomes was shown to hybridize with the leucoplast DNA (rho - 1.682 g/ml). Nuclear DNA of Polytoma (rho = 1.711) showed specific hybridization with rRNA from the 79S ribosomes. Saturation hybridization indicated that only one copy of the rRNA cistrons was present per leucoplast genome, with an average buoyant density of rho = 1.700. On the other hand, about 750 copies of the cytoplasmic rRNA cistrons were present per nuclear genome with a density of rho = 1.709. Heterologous hybridization studies with Chlamydomonas reinhardtii rRNAs showed an estimated 80% homology between the two cytoplasmic rRNAs, but only a 50% homology between chloroplast and leucoplast rRNAs of the two species. We conclude that the leucoplasts of Polytoma derive from chloroplasts of a Chlamydomonas-like ancestor, but that the leucoplast rRNA cistrons have diverged in evolution more extensively than the cistrons for cytoplasmic rRNA.     

UTILITY OF THE MITOCHONDRIAL nad4L GENE FOR ALGAL AND PROTISTAN PHYLOGENETIC ANALYSIS1

Mitochondrial gene sequences, in contrast to their plastid counterparts, have been rarely exploited for use in determining algal phylogenetic relationships. Their utilization would prove advantageous particularly for the assessment of evolutionary relationships among protistan lineages that include photo synthetic species and their colorless heterotrophic dependents that contain only a reduced plastid genome as well as ancestral nonphotosynthetic relatives. Recent studies that have examined mitochondrial protein (coxl, cox2, cob, and nad5) rather than ribosomal RNA genes confirm that these algal sequences are sufficiently conserved to be used routinely in algal systematics. In this investigation, we sequenced the mitochondrial nad4L gene, encoding a small (ca. 99-amino acid) polypeptide subunit of the NADH dehydrogenase complex, from two chromophytic algae and an oomycete to determine its suitability for use in determining algal and protistan evolutionary relationships. Our analysis demonstrated that nad4L could clearly resolve these three organisms as a distinct clade, the stramenopiles, to the exclusion of terrestrial plants/chlorophytes and eufungi. However, the position of deeply rooted species within the stramenopiles proved tenuous and would best be examined in conjunction with the analysis of longer mitochondrial protein gene sequences.     

Studies on the metabolism of the protozoa. The molecular structure of a starch-type polysaccharide from Polytoma uvella

1. Polytoma uvella, when grown in an acetate-containing medium, synthesizes a starch-type polysaccharide. 2. The starch differs from normal plant starches in having a lower iodine affinity, and in being relatively insoluble; this latter property makes fractionation difficult. 3. The starch contains about 16% of amylose, and on fractionation yields a branched amylopectin component that is similar in structure to a typical plant starch amylopectin.     

Codon Usage in Plastid Genes Is Correlated with Context, Position Within the Gene, and Amino Acid Content

Highly expressed plastid genes display codon adaptation, which is defined as a bias toward a set of codons which are complementary to abundant tRNAs. This type of adaptation is similar to what is observed in highly expressed Escherichia coli genes and is probably the result of selection to increase translation efficiency. In the current work, the codon adaptation of plastid genes is studied with regard to three specific features that have been observed in E. coli and which may influence translation efficiency. These features are (1) a relatively low codon adaptation at the 5′ end of highly expressed genes, (2) an influence of neighboring codons on codon usage at a particular site (codon context), and (3) a correlation between the level of codon adaptation of a gene and its amino acid content. All three features are found in plastid genes. First, highly expressed plastid genes have a noticeable decrease in codon adaptation over the first 10–20 codons. Second, for the twofold degenerate NNY codon groups, highly expressed genes have an overall bias toward the NNC codon, but this is not observed when the 3′ neighboring base is a G. At these sites highly expressed genes are biased toward NNT instead of NNC. Third, plastid genes that have higher codon adaptations also tend to have an increased usage of amino acids with a high G + C content at the first two codon positions and GNN codons in particular. The correlation between codon adaptation and amino acid content exists separately for both cytosolic and membrane proteins and is not related to any obvious functional property. It is suggested that at certain sites selection discriminates between nonsynonymous codons based on translational, not functional, differences, with the result that the amino acid sequence of highly expressed proteins is partially influenced by selection for increased translation efficiency. Received: 21 July 1999 / Accepted: 5 November 1999     

Nonessential plastid-encoded ribosomal proteins in tobacco: a developmental role for plastid translation and implications for reductive genome evolution

Plastid genomes of higher plants contain a conserved set of ribosomal protein genes. Although plastid translational activity is essential for cell survival in tobacco (Nicotiana tabacum), individual plastid ribosomal proteins can be nonessential. Candidates for nonessential plastid ribosomal proteins are ribosomal proteins identified as nonessential in bacteria and those whose genes were lost from the highly reduced plastid genomes of nonphotosynthetic plastid-bearing lineages (parasitic plants, apicomplexan protozoa). Here we report the reverse genetic analysis of seven plastid-encoded ribosomal proteins that meet these criteria. We have introduced knockout alleles for the corresponding genes into the tobacco plastid genome. Five of the targeted genes (ribosomal protein of the large subunit22 [rpl22], rpl23, rpl32, ribosomal protein of the small subunit3 [rps3], and rps16) were shown to be essential even under heterotrophic conditions, despite their loss in at least some parasitic plastid-bearing lineages. This suggests that nonphotosynthetic plastids show elevated rates of gene transfer to the nuclear genome. Knockout of two ribosomal protein genes, rps15 and rpl36, yielded homoplasmic transplastomic mutants, thus indicating nonessentiality. Whereas Δrps15 plants showed only a mild phenotype, Δrpl36 plants were severely impaired in photosynthesis and growth and, moreover, displayed greatly altered leaf morphology. This finding provides strong genetic evidence that chloroplast translational activity influences leaf development, presumably via a retrograde signaling pathway.