An unexpectedly large and loosely packed mitochondrial genome in the charophycean green alga Chlorokybus atmophyticus

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

Structure and Phylogeny of Chloroplast and Mitochondrial Genomes of a Chlorophycean Algae Pectinodesmus pectinatus (Scenedesmaceae, Sphaeropleales)

Pectinodesmus pectinatus is a green alga of commercial interest in sewage purification. Clarification of its organelle genomes is helpful for genetic manipulation, taxonomic revisions and evolutionary research. Here, de novo sequencing was used to determine chloroplast genome and mitochondrial genome of P. pectinatus strain F34. The chloroplast genome was composed of a large single-copy (LSC) region of 99,156 bp, a small single-copy (SSC) region of 70,665 bp, and a pair of inverted repeats (IRs) with a length of 13,494 bp each separated by LSC and SSC. The chloroplast genome contained 69 protein-coding genes, 25 transfer-RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes. The mitochondrial genome was 32,195 bp in length and consisted of 46 unique genes, including 16 protein-coding genes, 27 tRNA genes and 3 rRNA genes. The predominant mutations in organelle genomes were T/A to G/C transitions. Phylogenic analysis indicated P. pectinatus was a sister species to Tetradesmus obliquus and Hariotina sp. within the Pectinodesmus genus. In analysis with CGView Comparison Tool, P. pectinatus organelle genomes displayed the highest sequence similarity with that of T. obliquus. These findings advanced research on the taxonomy and phylogeny of Chlorophyceae algae and particularly revealed the role of P. pectinatus in microalgae evolution.

The complete mitochondrial genome of Isochrysis galbana harbors a unique repeat structure and a specific trans-spliced cox1 gene

The haptophyte Isochrysis galbana is considered as a promising source for food supplements due to its rich fucoxanthin and polyunsaturated fatty acids content. Here, the I. galbana mitochondrial genome (mitogenome) was sequenced using a combination of Illumina and PacBio sequencing platforms. This 39,258 bp circular mitogenome has a total of 46 genes, including 20 protein-coding genes, 24 tRNA genes and two rRNA genes. A large block of repeats (~12.7 kb) was segregated in one region of the mitogenome, accounting for almost one third of the total size. A trans-spliced gene cox1 was first identified in I. galbana mitogenome and was verified by RNA-seq and DNA-seq data. The massive expansion of tandem repeat size and cis- to trans-splicing shift could be explained by the high mitogenome rearrangement rates in haptophytes. Strict SNP calling based on deep transcriptome sequencing data suggested the lack of RNA editing in both organelles in this species, consistent with previous studies in other algal lineages. To gain insight into haptophyte mitogenome evolution, a comparative analysis of mitogenomes within haptophytes and among eight main algal lineages was performed. A core gene set of 15 energy and metabolism genes is present in haptophyte mitogenomes, consisting of 1 cob, 3 cox, 7 nad, 2 atp and 2 ribosomal genes. Gene content and order was poorly conserved in this lineage. Haptophyte mitogenomes have lost many functional genes found in many other eukaryotes including rps/rpl, sdh, tat, secY genes, which make it contain the smallest gene set among all algal taxa. All these implied the rapid-evolving and more recently evolved mitogenomes of haptophytes compared to other algal lineages. The phylogenetic tree constructed by cox1 genes of 204 algal mitogenomes yielded well-resolved internal relationships, providing new evidence for red-lineages that contained plastids of red algal secondary endosymbiotic origin. This newly assembled mitogenome will add to our knowledge of general trends in algal mitogenome evolution within haptophytes and among different algal taxa.

Group II intron and repeat-rich red algal mitochondrial genomes demonstrate the dynamic recent history of autocatalytic RNAs

Background

Group II introns are mobile genetic elements that can insert at specific target sequences, however, their origins are often challenging to reconstruct because of rapid sequence decay following invasion and spread into different sites. To advance understanding of group II intron spread, we studied the intron-rich mitochondrial genome (mitogenome) in the unicellular red alga, Porphyridium.

Results

Analysis of mitogenomes in three closely related species in this genus revealed they were 3–6-fold larger in size (56–132 kbp) than in other red algae, that have genomes of size 21–43 kbp. This discrepancy is explained by two factors, group II intron invasion and expansion of repeated sequences in large intergenic regions. Phylogenetic analysis demonstrates that many mitogenome group II intron families are specific to Porphyridium, whereas others are closely related to sequences in fungi and in the red alga-derived plastids of stramenopiles. Network analysis of intron-encoded proteins (IEPs) shows a clear link between plastid and mitochondrial IEPs in distantly related species, with both groups associated with prokaryotic sequences.

Conclusion

Our analysis of group II introns in Porphyridium mitogenomes demonstrates the dynamic nature of group II intron evolution, strongly supports the lateral movement of group II introns among diverse eukaryotes, and reveals their ability to proliferate, once integrated in mitochondrial DNA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01200-3.

Variation in protein gene and intron content among land plant mitogenomes

The functional content of the mitochondrial genome (mitogenome) is highly diverse across eukaryotes. Among land plants, our understanding of the variation in mitochondrial gene and intron content is improving from concerted efforts to densely sample mitogenomes from diverse land plants. Here I review the current state of knowledge regarding the diversity in content of protein genes and introns in the mitogenomes of all major land plant lineages. Mitochondrial protein gene content is largely conserved among mosses and liverworts, but it varies substantially among and within other land plant lineages due to convergent losses of genes encoding ribosomal proteins and, to a lesser extent, genes for proteins involved in cytochrome c maturation and oxidative phosphorylation. Mitochondrial intron content is fairly stable within each major land plant lineage, but highly variable among lineages, resulting from occasional gains and many convergent losses over time. Trans-splicing has evolved dozens of times in various vascular plant lineages, particularly those with relatively higher rates of mitogenomic rearrangement. Across eukaryotes, mitochondrial protein gene and intron content has been shaped massive convergent evolution.

Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution

The complete chloroplast and mitochondrial genomes of Charophyta have shed new light on land plant terrestrialization. Here, we report the organellar genomes of the Zygnema circumcarinatum strain UTEX 1559, and a comparative genomics investigation of 33 plastomes and 18 mitogenomes of Chlorophyta, Charophyta (including UTEX 1559 and its conspecific relative SAG 698-1a), and Embryophyta. Gene presence/absence was determined across these plastomes and mitogenomes. A comparison between the plastomes of UTEX 1559 (157 548 bp) and SAG 698-1a (165 372 bp) revealed very similar gene contents, but substantial genome rearrangements. Surprisingly, the two plastomes share only 85.69% nucleotide sequence identity. The UTEX 1559 mitogenome size is 215 954 bp, the largest among all sequenced Charophyta. Interestingly, this large mitogenome contains a 50 kb region without homology to any other organellar genomes, which is flanked by two 86 bp direct repeats and contains 15 ORFs. These ORFs have significant homology to proteins from bacteria and plants with functions such as primase, RNA polymerase, and DNA polymerase. We conclude that (i) the previously published SAG 698-1a plastome is probably from a different Zygnema species, and (ii) the 50 kb region in the UTEX 1559 mitogenome might be recently acquired as a mobile element.

Complete mitogenome of the streptophyte green alga Coleochaete scutata (Coleochaetophyceae)

We have sequenced the mitogenome of Coleochaete scutata strain SAG 110.80M. This mitogenome is the largest among the streptophyte green algae examined to date. At 242,024 bp, it is 4.3-fold larger than the mitogenome of Chaetosphaeridium globosum, the only other mitogenome available for the Coleochaetophyceae. This size difference is mainly explained by differences in the abundance of introns and in the length of intergenic regions containing vestiges of coding sequences thought to be of foreign origin. With 31 group I and 26 group II introns, the C. scutata mitogenome is the most intron-rich organelle genome known among streptophyte algae.

Dynamic evolution of mitochondrial genomes in Trebouxiophyceae, including the first completely assembled mtDNA from a lichen-symbiont microalga (Trebouxia sp. TR9)

Trebouxiophyceae (Chlorophyta) is a species-rich class of green algae with a remarkable morphological and ecological diversity. Currently, there are a few completely sequenced mitochondrial genomes (mtDNA) from diverse Trebouxiophyceae but none from lichen symbionts. Here, we report the mitochondrial genome sequence of Trebouxia sp. TR9 as the first complete mtDNA sequence available for a lichen-symbiont microalga. A comparative study of the mitochondrial genome of Trebouxia sp. TR9 with other chlorophytes showed important organizational changes, even between closely related taxa. The most remarkable change is the enlargement of the genome in certain Trebouxiophyceae, which is principally due to larger intergenic spacers and seems to be related to a high number of large tandem repeats. Another noticeable change is the presence of a relatively large number of group II introns interrupting a variety of tRNA genes in a single group of Trebouxiophyceae, which includes Trebouxiales and Prasiolales. In addition, a fairly well-resolved phylogeny of Trebouxiophyceae, along with other Chlorophyta lineages, was obtained based on a set of seven well-conserved mitochondrial genes.

The complete mitochondrial genome of the Caulerpa lentillifera (Ulvophyceae, Chlorophyta): Sequence, genome content, organization structure and phylogenetic consideration

The complete mitochondrial genome is greatly important for studies on genetic structure and phylogenetic relationship at various taxonomic levels. To obtain information about the evolutionary trends of mtDNA in the Ulvophyceae and also to gain insights into the phylogenetic relationships between ulvophytes and other chlorophytes, we determined the mtDNA sequence of Caulerpa lentillifera (sea grape) using de novo mitochondrial genome sequencing. The complete genomic DNA of C. lentillifera was circular and 209,034 bp in length, and it was the largest green-algal mitochondrial genome sequenced to date, with a low gene density of 65.2%, which is reminiscent of the "expanded" pattern of evolution exhibited by embryophyte mtDNAs. The C. lentillifera mtDNA consisted of a typical set of 17 protein-coding genes (PCGs), 20 transfer RNA (tRNA) genes, three ribosomal RNA (rRNA) genes, 42 putative open reading frames (ORFs) and 29 introns, which had homologs in green-algal mtDNAs displaying an "ancestral" or a "reduced-derived" pattern of evolution. The overall base composition of its mitochondrial genome was 24.19% for A, 24.94% for T, 25.80% for G, 25.07% for C and 50.87% for GC. The mitochondrial genome of C. lentillifera was characterized by numerous small intergenic regions and introns, which was clearly different from other green algae. With the exception of the NADH dehydrogenase subunit 6 (ND6), ND1, ATP and three tRNA genes (tRNA-His, tRNA-Thr and tRNA-Ala), all other mitochondrial genes were encoded on the heavy strand. All of the PCGs had ATG as their start codon and employed TAA, TGA or TAG as their termination codon. To gain insights into the evolutionary trends of mtDNA in the Ulvophyceae, we inferred the complete mtDNA sequence of C. lentillifera, an ulvophyte belonging to a distinct, early-diverging lineage. Taken together, our data offered useful information for the studies on phylogenetic hypotheses and phylogenetic relationships of C. lentillifera within the Chlorophyta.

New Taxa of Streptophyte Algae (Streptophyta) from Terrestrial Habitats Revealed Using an Integrative Approach

Two new genera (Streptosarcina and Streptofilum) and three new species (Streptosarcina arenaria, S. costaricana and Streptofilum capillatum) of streptophyte algae were detected in cultures isolated from terrestrial habitats of Europe and Central America and described using an integrative approach. Additionally, a strain isolated from soil in North America was identified as Hormidiella parvula and proposed as an epitype of this species. The molecular phylogeny based on 18S rRNA and rbcL genes, secondary structure of ITS-2, as well as the morphology of vegetative and reproductive stages, cell ultrastructure, ecology and distribution of the investigated strains were assessed. The new genus Streptosarcina forms a sister lineage to the genus Hormidiella (Klebsormidiophyceae). Streptosarcina is characterized by packet-like (sarcinoid) and filamentous thalli with true branching and a cell organization typical for Klebsormidiophyceae. Streptofilum forms a separate lineage within Streptophyta. This genus represents an easily disintegrating filamentous alga which exhibits a cell coverage of unique structure: layers of submicroscopic scales of piliform shape covering the plasmalemma and exfoliate inside the mucilage envelope surrounding cells. The implications of the discovery of the new taxa for understanding evolutionary tendencies in the Streptophyta, a group of great evolutionary interest, are discussed.

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C-169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.

Complete mitochondrial genomes of prasinophyte algae Pyramimonas parkeae and Cymbomonas tetramitiformis

Mitochondria are archetypal eukaryotic organelles that were acquired by endosymbiosis of an ancient species of alpha-proteobacteria by the last eukaryotic common ancestor. The genetic information contained within the mitochondrial genome has been an important source of information for resolving relationships among eukaryotic taxa. In this study, we utilized mitochondrial and chloroplast genomes to explore relationships among prasinophytes. Prasinophytes are represented by diverse early-diverging green algae whose physical structures and genomes have the potential to elucidate the traits of the last common ancestor of the Viridiplantae (or Chloroplastida). We constructed de novo mitochondrial genomes for two prasinophyte algal species, Pyramimonas parkeae and Cymbomonas tetramitiformis, representing the prasinophyte clade. Comparisons of genome structure and gene order between these species and to those of other prasinophytes revealed that the mitochondrial genomes of P. parkeae and C. tetramitiformis are more similar to each other than to other prasinophytes, consistent with other molecular inferences of the close relationship between these two species. Phylogenetic analyses using the inferred amino acid sequences of mitochondrial and chloroplast protein-coding genes resolved a clade consisting of P. parkeae and C. tetramitiformis; and this group (representing the prasinophyte clade I) branched with the clade II, consistent with previous studies based on the use of nuclear gene markers.

Mitochondrial Genome Structure of Photosynthetic Eukaryotes

Current ideas of plant mitochondrial genome organization are presented. Data on the size and structural organization of mtDNA, gene content, and peculiarities are summarized. Special emphasis is given to characteristic features of the mitochondrial genomes of land plants and photosynthetic algae that distinguish them from the mitochondrial genomes of other eukaryotes. The data published before the end of 2014 are reviewed.

Phylogenetic patterns of gene rearrangements in four mitochondrial genomes from the green algal family Hydrodictyaceae (Sphaeropleales, Chlorophyceae)

Background

The variability in gene organization and architecture of green algal mitochondrial genomes is only recently being studied on a finer taxonomic scale. Sequenced mt genomes from the chlorophycean orders Volvocales and Sphaeropleales exhibit considerable variation in size, content, and structure, even among closely related genera. However, sampling of mt genomes on a within-family scale is still poor and the sparsity of information precludes a thorough understanding of genome evolution in the green algae.

Methods

Genomic DNA of representative taxa were sequenced on an Illumina HiSeq2500 to produce 2x100 bp paired reads, and mitochondrial genomes were assembled and annotated using Geneious v.6.1.5. Phylogenetic analysis of 13 protein-coding mitochondrial genes spanning the Sphaeropleales was performed.

Results

This study presents one of the first within-family comparisons of mt genome diversity, and is the first to report complete mt genomes for the family Hydrodictyaceae (order Sphaeropleales). Four complete mt genomes representing three taxa and four phylogenetic groups, Stauridium tetras, Pseudopediastrum boryanum, and Pediastrum duplex, range in size from 37,723 to 53,560 bp. The size variability is primarily due to intergenic region expansion, and intron content is generally low compared with other mt genomes of Sphaeropleales.

Conclusions

Certain gene rearrangements appear to follow a phylogenetic pattern, and with a more thorough taxon sampling genome-level sequence may be useful in resolving systematic conundrums that plague this morphologically diverse family.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2056-5) contains supplementary material, which is available to authorized users.

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).

An intronic open reading frame was released from one of group II introns in the mitochondrial genome of the haptophyte Chrysochromulina sp. NIES-1333

Mitochondrial (mt) genome sequences, which often bear introns, have been sampled from phylogenetically diverse eukaryotes. Thus, we can anticipate novel insights into intron evolution from previously unstudied mt genomes. We here investigated the origins and evolution of three introns in the mt genome of the haptophyte Chrysochromulina sp. NIES-1333, which was sequenced completely in this study. All the three introns were characterized as group II, on the basis of predicted secondary structure, and the conserved sequence motifs at the 5′ and 3′ termini. Our comparative studies on diverse mt genomes prompt us to propose that the Chrysochromulina mt genome laterally acquired the introns from mt genomes in distantly related eukaryotes. Many group II introns harbor intronic open reading frames for the proteins (intron-encoded proteins or IEPs), which likely facilitate the splicing of their host introns. However, we propose that a “free-standing,” IEP-like protein, which is not encoded within any introns in the Chrysochromulina mt genome, is involved in the splicing of the first cox1 intron that lacks any open reading frames.

Tracing the evolution of streptophyte algae and their mitochondrial genome

Six monophyletic groups of charophycean green algae are recognized within the Streptophyta. Although incongruent with earlier studies based on genes from three cellular compartments, chloroplast and nuclear phylogenomic analyses have resolved identical relationships among these groups, placing the Zygnematales or the Zygnematales + Coleochaetales as sister to land plants. The present investigation aimed at determining whether this consensus view is supported by the mitochondrial genome and at gaining insight into mitochondrial DNA (mtDNA) evolution within and across streptophyte algal lineages and during the transition toward the first land plants. We present here the newly sequenced mtDNAs of representatives of the Klebsormidiales (Entransia fimbriata and Klebsormidium spec.) and Zygnematales (Closterium baillyanum and Roya obtusa) and compare them with their homologs in other charophycean lineages as well as in selected embryophyte and chlorophyte lineages. Our results indicate that important changes occurred at the levels of genome size, gene order, and intron content within the Zygnematales. Although the representatives of the Klebsormidiales display more similarity in genome size and intron content, gene order seems more fluid and gene losses more frequent than in other charophycean lineages. In contrast, the two members of the Charales display an extremely conservative pattern of mtDNA evolution. Collectively, our analyses of gene order and gene content and the phylogenies we inferred from 40 mtDNA-encoded proteins failed to resolve the relationships among the Zygnematales, Coleochaetales, and Charales; however, they are consistent with previous phylogenomic studies in favoring that the morphologically complex Charales are not sister to land plants.

From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes

BACKGROUND

Next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants (Viridiplantae). Although these data have helped resolve the phylogeny of numerous clades (e.g., green algae, angiosperms, and gymnosperms), their utility for inferring relationships across all green plants is uncertain. Viridiplantae originated 700-1500 million years ago and may comprise as many as 500,000 species. This clade represents a major source of photosynthetic carbon and contains an immense diversity of life forms, including some of the smallest and largest eukaryotes. Here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data.

RESULTS

We assembled protein-coding sequence data for 78 genes from 360 diverse green plant taxa with complete or nearly complete plastid genome sequences available from GenBank. Phylogenetic analyses of the plastid data recovered well-supported backbone relationships and strong support for relationships that were not observed in previous analyses of major subclades within Viridiplantae. However, there also is evidence of systematic error in some analyses. In several instances we obtained strongly supported but conflicting topologies from analyses of nucleotides versus amino acid characters, and the considerable variation in GC content among lineages and within single genomes affected the phylogenetic placement of several taxa.

CONCLUSIONS

Analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. This framework includes: i) the placement of Zygnematophyceace as sister to land plants (Embryophyta), ii) a clade of extant gymnosperms (Acrogymnospermae) with cycads + Ginkgo sister to remaining extant gymnosperms and with gnetophytes (Gnetophyta) sister to non-Pinaceae conifers (Gnecup trees), and iii) within the monilophyte clade (Monilophyta), Equisetales + Psilotales are sister to Marattiales + leptosporangiate ferns. Our analyses also highlight the challenges of using plastid genome sequences in deep-level phylogenomic analyses, and we provide suggestions for future analyses that will likely incorporate plastid genome sequence data for thousands of species. We particularly emphasize the importance of exploring the effects of different partitioning and character coding strategies.

The mitochondrial genome of the lycophyte Huperzia squarrosa: the most archaic form in vascular plants

Mitochondrial genomes have maintained some bacterial features despite their residence within eukaryotic cells for approximately two billion years. One of these features is the frequent presence of polycistronic operons. In land plants, however, it has been shown that all sequenced vascular plant chondromes lack large polycistronic operons while bryophyte chondromes have many of them. In this study, we provide the completely sequenced mitochondrial genome of a lycophyte, from Huperzia squarrosa, which is a member of the sister group to all other vascular plants. The genome, at a size of 413,530 base pairs, contains 66 genes and 32 group II introns. In addition, it has 69 pseudogene fragments for 24 of the 40 protein- and rRNA-coding genes. It represents the most archaic form of mitochondrial genomes of all vascular plants. In particular, it has one large conserved gene cluster containing up to 10 ribosomal protein genes, which likely represents a polycistronic operon but has been disrupted and greatly reduced in the chondromes of other vascular plants. It also has the least rearranged gene order in comparison to the chondromes of other vascular plants. The genome is ancestral in vascular plants in several other aspects: the gene content resembling those of charophytes and most bryophytes, all introns being cis-spliced, a low level of RNA editing, and lack of foreign DNA of chloroplast or nuclear origin.

The mitochondrial genomes of the early land plants Treubia lacunosa and Anomodon rugelii: dynamic and conservative evolution

Early land plant mitochondrial genomes captured important changes of mitochondrial genome evolution when plants colonized land. The chondromes of seed plants show several derived characteristics, e.g., large genome size variation, rapid intra-genomic rearrangement, abundant introns, and highly variable levels of RNA editing. On the other hand, the chondromes of charophytic algae are still largely ancestral in these aspects, resembling those of early eukaryotes. When the transition happened has been a long-standing question in studies of mitochondrial genome evolution. Here we report complete mitochondrial genome sequences from an early-diverging liverwort, Treubia lacunosa, and a late-evolving moss, Anomodon rugelii. The two genomes, 151,983 and 104,239 base pairs in size respectively, contain standard sets of protein coding genes for respiration and protein synthesis, as well as nearly full sets of rRNA and tRNA genes found in the chondromes of the liverworts Marchantia polymorpha and Pleurozia purpurea and the moss Physcomitrella patens. The gene orders of these two chondromes are identical to those of the other liverworts and moss. Their intron contents, with all cis-spliced group I or group II introns, are also similar to those in the previously sequenced liverwort and moss chondromes. These five chondromes plus the two from the hornworts Phaeoceros laevis and Megaceros aenigmaticus for the first time allowed comprehensive comparative analyses of structure and organization of mitochondrial genomes both within and across the three major lineages of bryophytes. These analyses led to the conclusion that the mitochondrial genome experienced dynamic evolution in genome size, gene content, intron acquisition, gene order, and RNA editing during the origins of land plants and their major clades. However, evolution of this organellar genome has remained rather conservative since the origin and initial radiation of early land plants, except within vascular plants.

Interorganellar DNA transfer in wheat: dynamics and phylogenetic origin

A homology search of wheat chloroplast (ct) and mitochondrial (mt) genomes identified 54 ctDNA segments that have homology with 66 mtDNA segments. The mtDNA segments were classified according to their origin: orthologs (prokaryotic origin), xenologs (interorganellar DNA transfer origin) and paralogs (intraorganellar DNA amplification origin). The 66 mtDNA sequences with homology to ctDNA segments included 14 paralogs, 18 orthologs and 34 xenologs. Analysis of the xenologs indicated that the DNA transfer occurred unidirectionally from the ct genome to the mt genome. The evolutionary timing of each interorganellar DNA transfer that generated a xenolog was estimated. This analysis showed that 2 xenologs originated early in green plant evolution, 4 in angiosperm evolution, 3 in monocotyledon evolution, 9 during cereal diversification and 8 in the evolution of wheat. Six other xenologs showed recurrent transfer from the ct to mt genomes in more than one taxon. The two remaining xenologs were uninformative on the evolutionary timing of their transfer. The wheat mt nad9 gene was found to be chimeric, consisting of the cereal nad9 gene and its 291 bp 5'-flanking region that included a 58 bp xenolog of the ct-ndhC origin.

Evolution and diversity of plant cell walls: from algae to flowering plants

All photosynthetic multicellular Eukaryotes, including land plants and algae, have cells that are surrounded by a dynamic, complex, carbohydrate-rich cell wall. The cell wall exerts considerable biological and biomechanical control over individual cells and organisms, thus playing a key role in their environmental interactions. This has resulted in compositional variation that is dependent on developmental stage, cell type, and season. Further variation is evident that has a phylogenetic basis. Plants and algae have a complex phylogenetic history, including acquisition of genes responsible for carbohydrate synthesis and modification through a series of primary (leading to red algae, green algae, and land plants) and secondary (generating brown algae, diatoms, and dinoflagellates) endosymbiotic events. Therefore, organisms that have the shared features of photosynthesis and possession of a cell wall do not form a monophyletic group. Yet they contain some common wall components that can be explained increasingly by genetic and biochemical evidence.

Lineage-specific group II intron gains and losses of the mitochondrial rps3 gene in gymnosperms

According to PCR assays and sequencing, we now report the shared presence of two rps3 introns, namely the rps3i74 and the rps3i249, in the mitochondria of all the classes representing the surviving lineages of gymnosperms, and unveil several lineages experiencing intron loss. Interestingly, the rps3 intron gains and losses within the four groups of gymnosperms let us sort out the Pinaceae and the non-Pinaceae into intron (+)- and intron (-)-lineages, respectively. Worthy of mention is also the finding that only Gnetum within the Gnetales harbours both the rps3 introns. This intron distribution pattern is consistent with the hypothesis that the two rps3 introns were likely present in the common ancestor of the seed plants and, then, independently lost in the non-Pinaceae during gymnosperm evolution. The derived secondary structural model of the novel group IIA intron improves our understanding of the significance and origin of the extraordinary length polymorphisms observed among rps3i249 orthologs. Despite the remarkable structural plasticity to adopt and reject introns, the rps3 mRNAs undergo accurate processing by splicing and extensive editing in gymnosperm mitochondria. This study provides additional insights into the evolutionarily high dynamics of mitochondrial introns which may come and go in closely related plant species. The turnover of the mitochondrial rps3 group II introns seen among lineages of seed plants further suggests that these introns might be an additional signature to discriminate between particularly cryptical taxonomic groups for which there is a need of a further evaluation of their evolutionary affiliation.

Uncovering the evolutionary origin of plant molecular processes: comparison of Coleochaete (Coleochaetales) and Spirogyra (Zygnematales) transcriptomes

The large and diverse land plant lineage is nested within a clade of fresh water green algae, the charophytes. Collection of genome-scale data for land plants and other organisms over the past decade has invigorated the field of evolutionary biology. One of the core questions in the field asks: how did a colonization event by a green algae over 450 mya lead to one of the most successful lineages on the tree of life? This question can best be answered using the comparative method, the first step of which is to gather genome-scale data across closely related lineages to land plants. Before sequencing an entire genome it is useful to first gather transcriptome data: it is less expensive, it targets the protein coding regions of the genome, and provides support for gene models for future genome sequencing. We built Expressed Sequence Tag (EST) libraries for two charophyte species, Coleochaete orbicularis (Coleochaetales) and Spirogyra pratensis (Zygnematales). We used both Sanger sequencing and next generation 454 sequencing to cover as much of the transcriptome as possible.

RESULTS

Our sequencing effort for Spirogyra pratensis yielded 9,984 5' Sanger reads plus 598,460 GS FLX Standard 454 sequences; Coleochaete orbicularis yielded 4,992 5' Sanger reads plus 673,811 GS FLX Titanium 454 sequences. After clustering S. pratensis yielded 12,000 unique transcripts, or unigenes, and C. orbicularis yielded 19,000. Both transcriptomes were very plant-like, i.e. most of the transcripts were more similar to streptophytes (land plants + charophyte green algae) than to other green algae in the sister group chlorophytes. BLAST results of several land plant genes hypothesized to be important in early land plant evolution resulted in high quality hits in both transcriptomes revealing putative orthologs ripe for follow-up studies.

CONCLUSIONS

Two main conclusions were drawn from this study. One illustrates the utility of next generation sequencing for transcriptome studies: larger scale data collection at a lower cost enabled us to cover a considerable portion of the transcriptome for both species. And, two, that the charophyte green algal transcriptoms are remarkably plant-like, which gives them the unique capacity to be major players for future evolutionary genomic studies addressing origin of land plant questions.

A deviant genetic code in the reduced mitochondrial genome of the picoplanktonic green alga Pycnococcus provasolii

Reduction in size of flagellated chlorophytes occurred multiple times during evolution, providing the opportunity to study the consequences of cell reduction on genome architecture. Recent investigations on the chloroplast genomes of the tiny prasinophyceans Ostreococcus tauri (Mamiellales), Micromonas sp. RCC299 (Mamiellales), and Pycnococcus provasolii (Pseudocourfieldiales) highlighted their extreme compaction and reduced gene repertoires. Genome compaction is also exemplified by the Ostreococcus and Micromonas mitochondrial DNAs (mtDNAs) although they have retained almost all of the about 65 genes presumably present in the mitochondria of ancestral prasinophyceans. In this study, the mitochondrial genome of Pycnococcus was sequenced and compared to those of previously examined chlorophytes. Our results document the first case where cellular reduction of a free-living alga was accompanied by marked reduction in gene content of both the mitochondrial and chloroplast genomes. At 24,321 bp, the intronless Pycnococcus mitochondrial genome falls within the lower size range displayed by green algal mtDNAs. The 36 conserved genes, specifying two rRNAs with conventional structures, 16 tRNAs and 18 proteins, are all encoded on the same DNA strand and represent 88% of the genome. Besides a pronounced codon bias, the protein-coding genes feature a variant genetic code characterized by the use of TGA (normally a stop codon) to code for tryptophan, and the unprecedented use of TTA and TTG (normally leucine codons) as stop codons. We conclude that substantial reduction of the mitochondrial genome occurred in at least three independent chlorophyte lineages and that this process entailed a number of convergent changes in these lineages.

MITOCHONDRIAL DNA IN THE OOGAMOCHLAMYS CLADE (CHLOROPHYCEAE): HIGH GC CONTENT AND UNIQUE GENOME ARCHITECTURE FOR GREEN ALGAE(1)

Most mitochondrial genomes in the green algal phylum Chlorophyta are AT-rich, circular-mapping DNA molecules. However, mitochondrial genomes from the Reinhardtii clade of the Chlorophyceae lineage are linear and sometimes fragmented into subgenomic forms. Moreover, Polytomella capuana, from the Reinhardtii clade, has an elevated GC content (57.2%). In the present study, we examined mitochondrial genome conformation and GC bias in the Oogamochlamys clade of the Chlorophyceae, which phylogenetic data suggest is closely related to the Reinhardtii clade. Total DNA from selected Oogamochlamys taxa, including four Lobochlamys culleus (H. Ettl) Pröschold, B. Marin, U. G. Schlöss. et Melkonian strains, Lobochlamys segnis (H. Ettl) Pröschold, B. Marin, U. G. Schlöss. et Melkonian, and Oogamochlamys gigantea (O. Dill) Pröschold, B. Marin, U. G. Schlöss. et Melkonian, was subjected to Southern blot analyses with cob and cox1 probes, and the results suggest that the mitochondrial genome of these taxa is represented by multiple-sized linear DNA fragments with overlapping homologies. On the basis of these data, we propose that linear mitochondrial DNA with a propensity to become fragmented arose in an ancestor common to the Reinhardtii and Oogamochlamys clades or even earlier in the evolutionary history of the Chlorophyceae. Analyses of partial cob and cox1 sequences from these Oogamochlamys taxa revealed an unusually high GC content (49.9%-65.1%) and provided evidence for the accumulation of cob and cox1 pseudogenes and truncated sequences in the mitochondrial genome of all L. culleus strains examined.

Ribosomal protein L10 is encoded in the mitochondrial genome of many land plants and green algae

Background

The mitochondrial genomes of plants generally encode 30-40 identified protein-coding genes and a large number of lineage-specific ORFs. The lack of wide conservation for most ORFs suggests they are unlikely to be functional. However, an ORF, termed orf-bryo1, was recently found to be conserved among bryophytes suggesting that it might indeed encode a functional mitochondrial protein.

Results

From a broad survey of land plants, we have found that the orf-bryo1 gene is also conserved in the mitochondria of vascular plants and charophycean green algae. This gene is actively transcribed and RNA edited in many flowering plants. Comparative sequence analysis and distribution of editing suggests that it encodes ribosomal protein L10 of the large subunit of the ribosome. In several lineages, such as crucifers and grasses, where the rpl10 gene has been lost from the mitochondrion, we suggest that a copy of the nucleus-encoded chloroplast-derived rpl10 gene may serve as a functional replacement.

Conclusion

Despite the fact that there are now over 20 mitochondrial genome sequences for land plants and green algae, this gene has remained unidentified and largely undetected until now because of the unlikely coincidence that most of the earlier sequences were from the few lineages that lack the intact gene. These results illustrate the power of comparative sequencing to identify novel genomic features.

The complete mitochondrial genome sequence of the hornwort Megaceros aenigmaticus shows a mixed mode of conservative yet dynamic evolution in early land plant mitochondrial genomes

Land plants possess some of the most unusual mitochondrial genomes among eukaryotes. However, in early land plants these genomes resemble those of green and red algae or early eukaryotes. The question of when during land plant evolution the dramatic change in mtDNAs occurred remains unanswered. Here we report the first completely sequenced mitochondrial genome of the hornwort, Megaceros aenigmaticus, a member of the sister group of vascular plants. It is a circular molecule of 184,908 base pairs, with 32 protein genes, 3 rRNA genes, 17 tRNA genes, and 30 group II introns. The genome contains many genes arranged in the same order as in those of a liverwort, a moss, several green and red algae, and Reclinomonas americana, an early-branching eukaryote with the most ancestral form of mtDNA. In particular, the gene order between mtDNAs of the hornwort and Physcomitrella patens (moss) differs by only 8 inversions and translocations. However, the hornwort mtDNA possesses 4 derived features relative to green alga mtDNAs--increased genome size, RNA editing, intron gains, and gene losses--which were all likely acquired during the origin and early evolution of land plants. Overall, this genome and those of other 2 bryophytes show that mitochondrial genomes in early land plants, unlike their seed plant counterparts, exhibit a mixed mode of conservative yet dynamic evolution.

Exploring plant biodiversity: the Physcomitrella genome and beyond

For decades, plant molecular biology has focused on only a few angiosperm species. Recently, the approximately 500mega base pairs (Mb) of the haploid Physcomitrella patens genome were sequenced and annotated. Mosses such as P. patens occupy a key evolutionary position halfway between green algae and flowering plants. This draft genome, in comparison to existing genome data from other plants, allows evolutionary insights into the conquest of land by plants and the molecular biodiversity that land plants exhibit. As a model organism, P. patens provides a well-developed molecular toolbox, including efficient gene targeting in combination with the morphologically simple moss tissues. We describe current as well as future tools for P. patens research and the prospects they offer for plant research in general.

Nucleotide diversity in the mitochondrial and nuclear compartments of Chlamydomonas reinhardtii: investigating the origins of genome architecture

Background

The magnitude of intronic and intergenic DNA can vary substantially both within and among evolutionary lineages; however, the forces responsible for this disparity in genome compactness are conjectural. One explanation, termed the mutational-burden hypothesis, posits that genome compactness is primarily driven by two nonadaptive processes: mutation and random genetic drift – the effects of which can be discerned by measuring the nucleotide diversity at silent sites (π_silent), defined as noncoding sites and the synonymous sites of protein-coding regions. The mutational-burden hypothesis holds that π_silent is negatively correlated to genome compactness. We used the model organism Chlamydomonas reinhardtii, which has a streamlined, coding-dense mitochondrial genome and an noncompact, intron-rich nuclear genome, to investigate the mutational-burden hypothesis. For measuring π_silent we sequenced the complete mitochondrial genome and portions of 7 nuclear genes from 7 geographical isolates of C. reinhardtii.

Results

We found significantly more nucleotide diversity in the nuclear compartment of C. reinhardtii than in the mitochondrial compartment: net values of π_silent for the nuclear and mitochondrial genomes were 32 × 10^-3 and 8.5 × 10^-3, respectively; and when insertions and deletions (indels) are factored in, these values become 49 × 10^-3 for the nuclear DNA and 11 × 10^-3 for the mitochondrial DNA (mtDNA). Furthermore, our investigations of C. reinhardtii revealed 4 previously undiscovered mitochondrial introns, one of which contains a fragment of the large-subunit (LSU) rRNA gene and another of which is found in a region of the LSU-rRNA gene not previously reported (for any taxon) to contain introns.

Conclusion

At first glance our results are in opposition to the mutational-burden hypothesis: π_silent was approximately 4 times greater in the nuclear compartment of C. reinhardtii relative to the mitochondrial compartment. However, when we consider the encumbrance of noncoding DNA in each of these C. reinhardtii compartments, we conclude that introns in the mtDNA impose a greater burden than those in the nuclear DNA and suggest that the same may be true for the intergenic regions. Overall, we cannot reject the mutational-burden hypothesis and feel that more data on nucleotide diversity from green algae and other protists are needed.

Order within a mosaic distribution of mitochondrial c-type cytochrome biogenesis systems?

Mitochondrial cytochromes c and c(1) are present in all eukaryotes that use oxygen as the terminal electron acceptor in the respiratory chain. Maturation of c-type cytochromes requires covalent attachment of the heme cofactor to the protein, and there are at least five distinct biogenesis systems that catalyze this post-translational modification in different organisms and organelles. In this study, we use biochemical data, comparative genomic and structural bioinformatics investigations to provide a holistic view of mitochondrial c-type cytochrome biogenesis and its evolution. There are three pathways for mitochondrial c-type cytochrome maturation, only one of which is present in prokaryotes. We analyze the evolutionary distribution of these biogenesis systems, which include the Ccm system (System I) and the enzyme heme lyase (System III). We conclude that heme lyase evolved once and, in many lineages, replaced the multicomponent Ccm system (present in the proto-mitochondrial endosymbiont), probably as a consequence of lateral gene transfer. We find no evidence of a System III precursor in prokaryotes, and argue that System III is incompatible with multi-heme cytochromes common to bacteria, but absent from eukaryotes. The evolution of the eukaryotic-specific protein heme lyase is strikingly unusual, given that this protein provides a function (thioether bond formation) that is also ubiquitous in prokaryotes. The absence of any known c-type cytochrome biogenesis system from the sequenced genomes of various trypanosome species indicates the presence of a third distinct mitochondrial pathway. Interestingly, this system attaches heme to mitochondrial cytochromes c that contain only one cysteine residue, rather than the usual two, within the heme-binding motif. The isolation of single-cysteine-containing mitochondrial cytochromes c from free-living kinetoplastids, Euglena and the marine flagellate Diplonema papillatum suggests that this unique form of heme attachment is restricted to, but conserved throughout, the protist phylum Euglenozoa.