Comparative Analysis of the Complete Plastid Genome Sequence of the Red Alga Gracilaria tenuistipitata var. liui Provides Insights into the Evolution of Rhodoplasts and Their Relationship to Other Plastids

We sequenced to completion the circular plastid genome of the red alga Gracilaria tenuistipitata var. liui. This is the first plastid genome sequence from the subclass Florideophycidae (Rhodophyta). The genome is composed of 183,883 bp and contains 238 predicted genes, including a single copy of the ribosomal RNA operon. Comparisons with the plastid genome of Porphyra pupurea reveal strong conservation of gene content and order, but we found major genomic rearrangements and the presence of coding regions that are specific to Gracilaria. Phylogenetic analysis of a data set of 41 concatenated proteins from 23 plastid and two cyanobacterial genomes support red algal plastid monophyly and a specific evolutionary relationship between the Florideophycidae and the Bangiales. Gracilaria maintains a surprisingly ancient gene content in its plastid genome and, together with other Rhodophyta, contains the most complete repertoire of plastid genes known in photosynthetic eukaryotes.Supplementary material () is available for this article.[Reviewing Editor: Dr. W. Ford Doolittle]     

Parallel Loss of Plastid Introns and Their Maturase in the Genus Cuscuta

Plastid genome content and arrangement are highly conserved across most land plants and their closest relatives, streptophyte algae, with nearly all plastid introns having invaded the genome in their common ancestor at least 450 million years ago. One such intron, within the transfer RNA trnK-UUU, contains a large open reading frame that encodes a presumed intron maturase, matK. This gene is missing from the plastid genomes of two species in the parasitic plant genus Cuscuta but is found in all other published land plant and streptophyte algal plastid genomes, including that of the nonphotosynthetic angiosperm Epifagus virginiana and two other species of Cuscuta. By examining matK and plastid intron distribution in Cuscuta, we add support to the hypothesis that its normal role is in splicing seven of the eight group IIA introns in the genome. We also analyze matK nucleotide sequences from Cuscuta species and relatives that retain matK to test whether changes in selective pressure in the maturase are associated with intron deletion. Stepwise loss of most group IIA introns from the plastid genome results in substantial change in selective pressure within the hypothetical RNA-binding domain of matK in both Cuscuta and Epifagus, either through evolution from a generalist to a specialist intron splicer or due to loss of a particular intron responsible for most of the constraint on the binding region. The possibility of intron-specific specialization in the X-domain is implicated by evidence of positive selection on the lineage leading to C. nitida in association with the loss of six of seven introns putatively spliced by matK. Moreover, transfer RNA gene deletion facilitated by parasitism combined with an unusually high rate of intron loss from remaining functional plastid genes created a unique circumstance on the lineage leading to Cuscuta subgenus Grammica that allowed elimination of matK in the most species-rich lineage of Cuscuta.     

Sporadic Distribution of tRNACCUArg Introns among α-Purple Bacteria: Evidence for Horizontal Transmission and Transposition of a Group I Intron

A group I intron interrupts the tRNA_CCU^Arg gene of the α-purple bacterium Agrobacterium tumefaciens (B. Reinhold-Hurek and D. A. Shub, Nature [London] 357:173–176, 1992). In this study, we assess the distribution of the corresponding intron among 12 additional species of α-purple bacteria. Of 10 newly identified tRNA_CCU^Arg genes, we found only two that contained an intron homologous to that of the Agrobacterium tRNA_CCU^Arg intron. This restricted and scattered distribution of the tRNA_CCU^Arg intron among α-purple bacteria is consistent with a recent origin and horizontal transmission. Primary and secondary structural similarities between tRNA_UAA^Leu introns found in strains of the cyanobacterium Microcystis aeruginosa (K. Rudi and K. S. Jacobsen, FEMS Microbiol. Lett. 156:293–298, 1997) and α-purple tRNA_CCU^Arg introns suggest that these introns share a more recent common ancestor than either does with other known cyanobacterial tRNA_UAA^Leu introns.     

Evidence for Transitional Stages in the Evolution of Euglenid Group II Introns and Twintrons in the Monomorphina aenigmatica Plastid Genome

Background

Photosynthetic euglenids acquired their plastid by secondary endosymbiosis of a prasinophyte-like green alga. But unlike its prasinophyte counterparts, the plastid genome of the euglenid Euglena gracilis is riddled with introns that interrupt almost every protein-encoding gene. The atypical group II introns and twintrons (introns-within-introns) found in the E. gracilis plastid have been hypothesized to have been acquired late in the evolution of euglenids, implying that massive numbers of introns may be lacking in other taxa. This late emergence was recently corroborated by the plastid genome sequences of the two basal euglenids, Eutreptiella gymnastica and Eutreptia viridis, which were found to contain fewer introns.

Methodology/Principal Findings

To gain further insights into the proliferation of introns in euglenid plastids, we have characterized the complete plastid genome sequence of Monomorphina aenigmatica, a freshwater species occupying an intermediate phylogenetic position between early and late branching euglenids. The M. aenigmatica UTEX 1284 plastid genome (74,746 bp, 70.6% A+T, 87 genes) contains 53 intron insertion sites, of which 41 were found to be shared with other euglenids including 12 of the 15 twintron insertion sites reported in E. gracilis.

Conclusions

The pattern of insertion sites suggests an ongoing but uneven process of intron gain in the lineage, with perhaps a minimum of two bursts of rapid intron proliferation. We also identified several sites that represent intermediates in the process of twintron evolution, where the external intron is in place, but not the internal one, offering a glimpse into how these convoluted molecular contraptions originate.     

Evolution of the plastid ribosomal RNA operon in a nongreen parasitic plant: Accelerated sequence evolution,altered promoter structure,and tRNA pseudogenes

The nucleotide sequence of a 7.4 kb region containing the entire plastid ribosomal RNA operon of the nongreen parasitic plant Epifagus virginiana has been determined. Analysis of the sequence indicates that all four rRNA genes are intact and almost certainly functional. In contrast, the split genes for tRNA^Ile and tRNA^Ala present in the 16S-23S rRNA spacer region have become pseudogenes, and deletion upstream of the 16S rRNA gene has removed a tRNA^Val gene and most of the promoter region for the rRNA operon. The rate of nucleotide substitution in 16S and 23S rRNAs is several times higher in Epifagus than in tobacco, a related photosynthetic plant. Possible reasons for this, including relaxed translational constraints, are discussed.     

Piecing together the puzzle of parasitic plant plastome evolution

The importance of photosynthesis as a mode of energy production has put plastid genomes of plants under a constant purifying selection. This has shaped the characteristic features of plastid genomes across the entire spectrum of photosynthetic plants and has led to a highly uniform and conserved plastid genome with respect to structure, size, gene order, intron and editing site positions and coding capacity. Parasitic species that have dropped photosynthesis as the main energy provider share striking deviations from the plastid genome norm: multiple rearrangements within the circular chromosome, pseudogenization and gene deletions, promoter losses, intron losses as well as the extensive loss of mRNA editing competence have been reported. The collective loss of larger sets of functionally related genes like those for the plastid NADH–dehydrogenase complex and concomitant losses of RNA polymerase genes together with their target promoters point to “domino effects” where an initial loss might have triggered others. An example, which will be discussed in more detail, is the concomitant loss of the intron maturase gene matK and all introns that are supposedly subject to MatK-dependent splicing in two Cuscuta species.     

Organization of the Mitochondrial Genome of a Deep-Sea Fish, Gonostoma gracile (Teleostei: Stomiiformes): First Example of Transfer RNA Gene Rearrangements in Bony Fishes

We determined the complete nucleotide sequence of the mitochondrial genome (except for a portion of the putative control region) for a deep-sea fish, Gonostoma gracile. The entire mitochondrial genome was purified by gene amplification using long polymerase chain reaction (long PCR), and the products were subsequently used as templates for PCR with 30 sets of newly designed, fish-universal primers that amplify contiguous, overlapping segments of the entire genome. Direct sequencing of the PCR products showed that the genome contained the same 37 mitochondrial structural genes as found in other vertebrates (two ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes), with the order of all rRNA and protein-coding genes, and 19 tRNA genes being identical to that in typical vertebrates. The gene order of the three tRNAs (tRNA^Glu, tRNA^Thr, and tRNA^Pro) relative to cytochrome b, however, differed from that determined in other vertebrates. Two steps of tandem duplication of gene regions, each followed by deletions of genes, can be invoked as mechanisms generating such rearrangements of tRNAs. This is the first example of tRNA gene rearrangements in a bony fish mitochondrial genome. Received August 5, 1998; accepted February 19, 1999.     

The Mitochondrial Genome of the Prasinophyte Prasinoderma coloniale Reveals Two Trans-Spliced Group I Introns in the Large Subunit rRNA Gene

Organelle genes are often interrupted by group I and or group II introns. Splicing of these mobile genetic occurs at the RNA level via serial transesterification steps catalyzed by the introns''own tertiary structures and, sometimes, with the help of external factors. These catalytic ribozymes can be found in cis or trans configuration, and although trans-arrayed group II introns have been known for decades, trans-spliced group I introns have been reported only recently. In the course of sequencing the complete mitochondrial genome of the prasinophyte picoplanktonic green alga Prasinoderma coloniale CCMP 1220 (Prasinococcales, clade VI), we uncovered two additional cases of trans-spliced group I introns. Here, we describe these introns and compare the 54,546 bp-long mitochondrial genome of Prasinoderma with those of four other prasinophytes (clades II, III and V). This comparison underscores the highly variable mitochondrial genome architecture in these ancient chlorophyte lineages. Both Prasinoderma trans-spliced introns reside within the large subunit rRNA gene (rnl) at positions where cis-spliced relatives, often containing homing endonuclease genes, have been found in other organelles. In contrast, all previously reported trans-spliced group I introns occur in different mitochondrial genes (rns or coxI). Each Prasinoderma intron is fragmented into two pieces, forming at the RNA level a secondary structure that resembles those of its cis-spliced counterparts. As observed for other trans-spliced group I introns, the breakpoint of the first intron maps to the variable loop L8, whereas that of the second is uniquely located downstream of P9.1. The breakpoint In each Prasinoderma intron corresponds to the same region where the open reading frame (ORF) occurs when present in cis-spliced orthologs. This correlation between the intron breakpoint and the ORF location in cis-spliced orthologs also holds for other trans-spliced introns; we discuss the possible implications of this interesting observation for trans-splicing of group I introns.     

Organization of the Mitochondrial Genome of Antarctic Krill <Emphasis Type="Italic">Euphausia superba</Emphasis> (Crustacea: Malacostraca)

We determined the nearly complete DNA sequence of the mitochondrial genome of Antarctic krill Euphausia superba (Crustacea: Malacostraca), one of the most ecologically and commercially important zooplankters in Antarctic waters. All of the genome sequences were purified by gene amplification using long polymerase chain reaction (PCR), and the products were subsequently used as templates for either direct sequencing using a primer-walking strategy or nested PCR with crustacea-versatile primers. Although we were unable to determine a portion of the genome owing to technical difficulties, the sequenced position, 14,606 bp long, contained all of the 13 protein-coding genes, 19 of the 22 transfer RNA genes, and the large subunit as well as a portion of the small subunit ribosomal RNA genes. Gene rearrangement was observed for 3 transfer RNA genes (tRNA^Cys, tRNA^Tyr, and tRNA^Trp) and the 2 leucine tRNA genes.     

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA^His differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.     

Plant nuclear tRNAMet genes are ubiquitously interrupted by introns

We have isolated three independent clones for nuclear elongator tRNA^Met genes from an Arabidopsis DNA library using a tRNA^Met-specific probe generated by PCR. Each of the coding sequences for tRNA^Met in these clones is identical and is interrupted by an identical 11 bp long intervening sequence at the same position in the anticodon loop of the tRNA. Their sequences differ at two positions from the intron in a soybean counterpart. Southern analysis of Arabidopsis DNA demonstrates that a gene family coding for tRNA^Met is dispersed at at least eight loci in the genome. The unspliced precursor tRNA^Met intermediate was detected by RNA analysis using an oligonucleotide probe complementary to the putative intron sequence. In order to know whether introns commonly interrupt plant tRNA^Met genes, their coding sequences were PCR-amplified from the DNAs of eight phylogenetically separate plant species. All 53 sequences determined contain 10 to 13 bp long intervening sequences, always positioned one base downstream from the anticodon. They can all be potentially folded into the secondary structure characteristic for plant intron-containing precursor tRNAs. Surprisingly, GC residues are always present at the 5-distal end of each intron.     

The Plastid Genome of the Red Macroalga Grateloupia taiwanensis (Halymeniaceae)

The complete plastid genome sequence of the red macroalga Grateloupia taiwanensis S.-M.Lin & H.-Y.Liang (Halymeniaceae, Rhodophyta) is presented here. Comprising 191,270 bp, the circular DNA contains 233 protein-coding genes and 29 tRNA sequences. In addition, several genes previously unknown to red algal plastids are present in the genome of G. taiwanensis. The plastid genomes from G. taiwanensis and another florideophyte, Gracilaria tenuistipitata var. liui, are very similar in sequence and share significant synteny. In contrast, less synteny is shared between G. taiwanensis and the plastid genome representatives of Bangiophyceae and Cyanidiophyceae. Nevertheless, the gene content of all six red algal plastid genomes here studied is highly conserved, and a large core repertoire of plastid genes can be discerned in Rhodophyta.     

Ferritin gene organization: Differences between plants and animals suggest possible kingdom-specific selective constraints

Ferritin, a protein widespread in nature, concentrates iron ∼10¹¹–10¹²-fold above the solubility within a spherical shell of 24 subunits; it derives in plants and animals from a common ancestor (based on sequence) but displays a cytoplasmic location in animals compared to the plastid in contemporary plants. Ferritin gene regulation in plants and animals is altered by development, hormones, and excess iron; iron signals target DNA in plants but mRNA in animals. Evolution has thus conserved the two end points of ferritin gene expression, the physiological signals and the protein structure, while allowing some divergence of the genetic mechanisms. Comparison of ferritin gene organization in plants and animals, made possible by the cloning of a dicot (soybean) ferritin gene presented here and the recent cloning of two monocot (maize) ferritin genes, shows evolutionary divergence in ferritin gene organization between plants and animals but conservation among plants or among animals; divergence in the genetic mechanism for iron regulation is reflected by the absence in all three plant genes of the IRE, a highly conserved, noncoding sequence in vertebrate animal ferritin mRNA. In plant ferritin genes, the number of introns (n= 7) is higher than in animals (n= 3). Second, no intron positions are conserved when ferritin genes of plants and animals are compared, although all ferritin gene introns are in the coding region; within kingdoms, the intron positions in ferritin genes are conserved. Finally, secondary protein structure has no apparent relationship to intron/exon boundaries in plant ferritin genes, whereas in animal ferritin genes the correspondence is high. The structural differences in introns/exons among phylogenetically related ferritin coding sequences and the high conservation of the gene structure within plant or animal kingdoms suggest that kingdom-specific functional constraints may exist to maintain a particular intron/exon pattern within ferritin genes. In the case of plants, where ferritin gene intron placement is unrelated to triplet codons or protein structure, and where ferritin is targeted to the plastid, the selection pressure on gene organization may relate to RNA function and plastid/nuclear signaling. Received: 25 July 1995 / Accepted: 3 October 1995     

The Complete Chloroplast and Mitochondrial Genomes of the Green Macroalga Ulva sp. UNA00071828 (Ulvophyceae,Chlorophyta)

Sequencing mitochondrial and chloroplast genomes has become an integral part in understanding the genomic machinery and the phylogenetic histories of green algae. Previously, only three chloroplast genomes (Oltmannsiellopsis viridis, Pseudendoclonium akinetum, and Bryopsis hypnoides) and two mitochondrial genomes (O. viridis and P. akinetum) from the class Ulvophyceae have been published. Here, we present the first chloroplast and mitochondrial genomes from the ecologically and economically important marine, green algal genus Ulva. The chloroplast genome of Ulva sp. was 99,983 bp in a circular-mapping molecule that lacked inverted repeats, and thus far, was the smallest ulvophycean plastid genome. This cpDNA was a highly compact, AT-rich genome that contained a total of 102 identified genes (71 protein-coding genes, 28 tRNA genes, and three ribosomal RNA genes). Additionally, five introns were annotated in four genes: atpA (1), petB (1), psbB (2), and rrl (1). The circular-mapping mitochondrial genome of Ulva sp. was 73,493 bp and follows the expanded pattern also seen in other ulvophyceans and trebouxiophyceans. The Ulva sp. mtDNA contained 29 protein-coding genes, 25 tRNA genes, and two rRNA genes for a total of 56 identifiable genes. Ten introns were annotated in this mtDNA: cox1 (4), atp1 (1), nad3 (1), nad5 (1), and rrs (3). Double-cut-and-join (DCJ) values showed that organellar genomes across Chlorophyta are highly rearranged, in contrast to the highly conserved organellar genomes of the red algae (Rhodophyta). A phylogenomic investigation of 51 plastid protein-coding genes showed that Ulvophyceae is not monophyletic, and also placed Oltmannsiellopsis (Oltmannsiellopsidales) and Tetraselmis (Chlorodendrophyceae) closely to Ulva (Ulvales) and Pseudendoclonium (Ulothrichales).     

Complete Mitochondrial DNA Sequence of Conger myriaster (Teleostei: Anguilliformes): Novel Gene Order for Vertebrate Mitochondrial Genomes and the Phylogenetic Implications for Anguilliform Families

The complete nucleotide sequence of the mitochondrial genome was determined for a conger eel, Conger myriaster (Elopomorpha: Anguilliformes), using a PCR-based approach that employs a long PCR technique and many fish-versatile primers. Although the genome [18,705 base pairs (bp)] contained the same set of 37 mitochondrial genes [two ribosomal RNA (rRNA), 22 transfer RNA (tRNA), and 13 protein-coding genes] as found in other vertebrates, the gene order differed from that recorded for any other vertebrates. In typical vertebrates, the ND6, tRNA^Glu, and tRNA^Pro genes are located between the ND5 gene and the control region, whereas the former three genes, in C. myriaster, have been translocated to a position between the control region and the tRNA^Phe gene that are contiguously located at the 5′ end of the 12S rRNA gene in typical vertebrates. This gene order is similar to the recently reported gene order in four lineages of birds in that the latter lack the ND6, tRNA^Glu, and tRNA^Pro genes between the ND5 gene and the control region; however, the relative position of the tRNA^Pro to the ND6–tRNA^Glu genes in C. myriaster was different from that in the four birds, which presumably resulted from different patterns of tandem duplication of gene regions followed by gene deletions in two distantly related groups of organisms. Sequencing of the ND5–cyt b region in 11 other anguilliform species, representing 11 families, plus one outgroup species, revealed that the same gene order as C. myriaster was shared by another 4 families, belonging to the suborder Congroidei. Although the novel gene orders of four lineages of birds were indicated to have multiple independent origins, phylogenetic analyses using nucleotide sequences from the mitochondrial 12S rRNA and cyt b genes suggested that the novel gene orders of the five anguilliform families had originated in a single ancestral species. Received: 13 July 2000 / Accepted: 30 November 2000     

Complete mitochondrial DNA sequence and phylogenetic analysis of Zhikong scallop <Emphasis Type="Italic">Chlamys farreri</Emphasis> (Bivalvia: Pectinidae)

The complete mitochondrial genome of Zhikong scallop Chlamys farreri is 21,695 bp in length and contains 12 protein-coding genes (the atp8 gene is absent, as in most bivalves), 2 ribosomal RNA genes, and 22 transfer RNA genes. The heavy strand has an overall A+T content of 58.7%. GC and AT skews for the mt genome of C. farreri are 0.337 and ?0.184, respectively, indicating the nucleotide bias against C and A. The mitochondrial gene order of C. farreri differs drastically from the scallops Argopecten irradians, Mimachlamys nobilis and Placopecten magellanicus, which belong to the same family Pectinidae. 6623 bp non-coding nucleotides exist intergenically in the mitogenome of C. farreri, with a large continuous sequence (4763 bp) between tRNA ^Val and tRNA ^Asn . Two repeat families are found in the large continuous sequence, which seems to be a common feature of scallops. Phylogenetic analysis based on 12 concatenated amino acid sequences of protein-coding genes supports the monophyly of Pectinidae and paraphyletic Pteriomorphia with respect to Heteroconchia.     

Mitochondrial tRNA gene clusters in Aspergillus nidulans: Organization and nucleotide sequence

The organization of tRNA genes on the circular 32 kb mitochondrial genome of the ascomycete Aspergillus nidulans has been studied by gel transfer hybridization and by DNA sequencing. Most of the tRNA genes are tightly clustered within two regions (1 kb each) flanking the split gene for the large ribosomal subunit RNA. The upstream cluster contains nine genes, the downstream cluster eleven genes. The twenty tRNA genes are on the same strand as the two rRNA genes and are separated from each other by AT-rich spacer sequences, usually consisting of only a few nucleotides. Two tRNA genes (leul and ala) are joined end to end. The occurrence of two tRNA^Gty genes is the first exception to the observation that in mitochondria all four-codon families are read by a single tRNA. Both genes are adjacent and show extensive sequence homology, suggesting relatively recent origin by gene duplication. The product of glyl has a U in the wobble position as do all other tRNA gene products specific for four-codon families, whereas the gly2 product, which has a rare A in the same position, should read only the codon GGU. The products of metl and thr have an A and G in positions 18 and 55, respectively, like the mitochondrial tRNA^fMet and tRNA^Thr of Neurospora crassa. Other unusual features are the replacement of the invariant G-C pair at positions 53 and 61 by A-T in met2, glyl and gly2, the replacement of the invariant T at position 8 by A in phe and G in pro and the deletion of a nucleotide at position 9 in ser2.     

Lateral transfer of introns in the cryptophyte plastid genome

Cryptophytes are unicellular eukaryotic algae that acquired photosynthesis secondarily through the uptake and retention of a red-algal endosymbiont. The plastid genome of the cryptophyte Rhodomonas salina CCMP1319 was recently sequenced and found to contain a genetic element similar to a group II intron. Here, we explore the distribution, structure and function of group II introns in the plastid genomes of distantly and closely related cryptophytes. The predicted secondary structures of six introns contained in three different genes were examined and found to be generally similar to group II introns but unusually large in size (including the largest known noncoding intron). Phylogenetic analysis suggests that the cryptophyte group II introns were acquired via lateral gene transfer (LGT) from a euglenid-like species. Unexpectedly, the six introns occupy five distinct genomic locations, suggesting multiple LGT events or recent transposition (or both). Combined with structural considerations, RT–PCR experiments suggest that the transferred introns are degenerate ‘twintrons’ (i.e. nested group II/group III introns) in which the internal intron has lost its splicing capability, resulting in an amalgamation with the outer intron.