Distribution of introns in the mitochondrial gene nad1 in land plants: phylogenetic and molecular evolutionary implications

Comprehensive analysis of the Lycopodium japonicum mitogenome reveals abundant tRNA genes and cis-spliced introns in Lycopodiaceae species

Lycophytes and ferns represent one of the earliest-diverging lineages of vascular plants, with the Lycopodiaceae family constituting the basal clade among lycophytes. In this research, we successfully assembled and annotated the complete Lycopodium japonicum Thunb. (L. japonicum) mitochondrial genome (mitogenome) utilizing PacBio HiFi sequencing data, resulting in a single circular molecule with a size of 454,458 bp. 64 unique genes were annotated altogether, including 34 protein-coding genes, 27 tRNAs and 3 rRNAs. It also contains 32 group II introns, all of which undergo cis-splicing. We identified 195 simple sequence repeats, 1,948 dispersed repeats, and 92 tandem repeats in the L. japonicum mitogenome. Collinear analysis indicated that the mitogenomes of Lycopodiaceae are remarkably conserved compared to those of other vascular plants. We totally identified 326 RNA editing sites in 31 unique protein-coding genes with 299 sites converting cytosine to uracil and 27 sites the reverse. Notably, the L. japonicum mitogenome has small amounts foreign DNA from plastid or nuclear origin, accounting for only 2.81% of the mitogenome. The maximum likelihood phylogenetic analysis based on 23 diverse land plant mitogenomes and plastid genomes supports the basal position of lycophytes within vascular plants and they form a sister clade to all other vascular lineages, which is consistent with the PPG I classification system. As the first reported mitogenome of Lycopodioideae subfamily, this study enriches our understanding of Lycopodium mitogenomes, and sets the stage for future research on mitochondrial diversity and evolution within the lycophytes and ferns.

The first mitogenome of Lauraceae (Cinnamomum chekiangense)

•The first reported mitochondrial genome (Cinnamomum chekiangense) of the Lauraceae family.•The mitogenome of C. chekiangense retains almost all of the ancestral protein-coding genes and has the highest RNA editing number in angiosperms.•Both of the plastid and mitochondrial phylogenetic trees support the magnoliids as a sister group of monocots and eudicots.

Assembly and comparative analysis of the complete mitochondrial genome of Trigonella foenum-graecum L

BACKGROUND

Trigonella foenum-graecum L. is a Leguminosae plant, and the stems, leaves, and seeds of this plant are rich in chemical components that are of high research value. The chloroplast (cp) genome of T. foenum-graecum has been reported, but the mitochondrial (mt) genome remains unexplored.

RESULTS

In this study, we used second- and third-generation sequencing methods, which have the dual advantage of combining high accuracy and longer read length. The results showed that the mt genome of T. foenum-graecum was 345,604 bp in length and 45.28% in GC content. There were 59 genes, including: 33 protein-coding genes (PCGs), 21 tRNA genes, 4 rRNA genes and 1 pseudo gene. Among them, 11 genes contained introns. The mt genome codons of T. foenum-graecum had a significant A/T preference. A total of 202 dispersed repetitive sequences, 96 simple repetitive sequences (SSRs) and 19 tandem repetitive sequences were detected. Nucleotide diversity (Pi) analysis counted the variation in each gene, with atp6 being the most notable. Both synteny and phylogenetic analyses showed close genetic relationship among Trifolium pratense, Trifolium meduseum, Trifolium grandiflorum, Trifolium aureum, Medicago truncatula and T. foenum-graecum. Notably, in the phylogenetic tree, Medicago truncatula demonstrated the highest level of genetic relatedness to T. foenum-graecum, with a strong support value of 100%. The interspecies non-synonymous substitutions (Ka)/synonymous substitutions (Ks) results showed that 23 PCGs had Ka/Ks < 1, indicating that these genes would continue to evolve under purifying selection pressure. In addition, setting the similarity at 70%, 23 homologous sequences were found in the mt genome of T. foenum-graecum.

CONCLUSIONS

This study explores the mt genome sequence information of T. foenum-graecum and complements our knowledge of the phylogenetic diversity of Leguminosae plants.

A mitochondrial pentatricopeptide repeat protein enhances cold tolerance by modulating mitochondrial superoxide in rice

Cold stress affects rice growth and productivity. Defects in the plastid-localized pseudouridine synthase OsPUS1 affect chloroplast ribosome biogenesis, leading to low-temperature albino seedlings and accumulation of reactive oxygen species (ROS). Here, we report an ospus1-1 suppressor, sop10. SOP10 encodes a mitochondria-localized pentatricopeptide repeat protein. Mutations in SOP10 impair intron splicing of the nad4 and nad5 transcripts and decrease RNA editing efficiency of the nad2, nad6, and rps4 transcripts, resulting in deficiencies in mitochondrial complex I, thus decrease ROS generation and rescuing the albino phenotype. Overexpression of different compartment-localized superoxide dismutases (SOD) genes in ospus1-1 reverses the ROS over-accumulation and albino phenotypes to various degrees, with Mn-SOD reversing the best. Mutation of SOP10 in indica rice varieties enhances cold tolerance with lower ROS levels. We find that the mitochondrial superoxide plays a key role in rice cold responses, and identify a mitochondrial superoxide modulating factor, informing efforts to improve rice cold tolerance.

Categorizing 161 plant (streptophyte) mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

Group II introns are common in the two endosymbiotic organelle genomes of the plant lineage. Chloroplasts harbor 22 positionally conserved group II introns whereas their occurrence in land plant (embryophyte) mitogenomes is highly variable and specific for the seven major clades: liverworts, mosses, hornworts, lycophytes, ferns, gymnosperms and flowering plants. Each plant group features "signature selections" of ca. 20-30 paralogues from a superset of altogether 105 group II introns meantime identified in embryophyte mtDNAs, suggesting massive intron gains and losses along the backbone of plant phylogeny. We report on systematically categorizing plant mitochondrial group II introns into "families", comprising evidently related paralogues at different insertion sites, which may even be more similar than their respective orthologues in phylogenetically distant taxa. Including streptophyte (charophyte) algae extends our sampling to 161 and we sort 104 streptophyte mitochondrial group II introns into 25 core families of related paralogues evidently arising from retrotransposition events. Adding to discoveries of only recently created intron paralogues, hypermobile introns and twintrons, our survey led to further discoveries including previously overlooked "fossil" introns in spacer regions or e.g., in the rps8 pseudogene of lycophytes. Initially excluding intron-borne maturase sequences for family categorization, we added an independent analysis of maturase phylogenies and find a surprising incongruence between intron mobility and the presence of intron-borne maturases. Intriguingly, however, we find that several examples of nuclear splicing factors meantime characterized simultaneously facilitate splicing of independent paralogues now placed into the same intron families. Altogether this suggests that plant group II intron mobility, in contrast to their bacterial counterparts, is not intimately linked to intron-encoded maturases.

Mitochondrial Genome Sequence of Salvia officinalis (Lamiales: Lamiaceae) Suggests Diverse Genome Structures in Cogeneric Species and Finds the Stop Gain of Genes through RNA Editing Events

Our previous study was the first to confirm that the predominant conformation of mitochondrial genome (mitogenome) sequence of Salvia species contains two circular chromosomes. To further understand the organization, variation, and evolution of Salvia mitogenomes, we characterized the mitogenome of Salvia officinalis. The mitogenome of S. officinalis was sequenced using Illumina short reads and Nanopore long reads and assembled using a hybrid assembly strategy. We found that the predominant conformation of the S. officinalis mitogenome also had two circular chromosomes that were 268,341 bp (MC1) and 39,827 bp (MC2) in length. The S. officinalis mitogenome encoded an angiosperm-typical set of 24 core genes, 9 variable genes, 3 rRNA genes, and 16 tRNA genes. We found many rearrangements of the Salvia mitogenome through inter- and intra-specific comparisons. A phylogenetic analysis of the coding sequences (CDs) of 26 common protein-coding genes (PCGs) of 11 Lamiales species and 2 outgroup taxa strongly indicated that the S. officinalis was a sister taxon to S. miltiorrhiza, consistent with the results obtained using concatenated CDs of common plastid genes. The mapping of RNA-seq data to the CDs of PCGs led to the identification of 451 C-to-U RNA editing sites from 31 PCGs of the S. officinalis mitogenome. Using PCR amplification and Sanger sequencing methods, we successfully validated 113 of the 126 RNA editing sites from 11 PCGs. The results of this study suggest that the predominant conformation of the S. officinalis mitogenome are two circular chromosomes, and the stop gain of rpl5 was found through RNA editing events of the Salvia mitogenome.

Assembly of the complete mitochondrial genome of an endemic plant, Scutellaria tsinyunensis, revealed the existence of two conformations generated by a repeat-mediated recombination

We assembled the complete mitochondrial genome of Scutellaria tsinyunensis in this study. Repeat-mediated recombination resulted in the formation of two conformations of the mitochondrial genome in S. tsinyunensis. Scutellaria tsinyunensis belongs to the family Lamiaceae, distributed only in the Jinyun Mountain, Chongqing, China. As a valuable endemic and small population species, it is regarded as a natural resource potentially with significant economic and ecological importance. In this study, we assembled a complete and gap-free mitochondrial genome of S. tsinyunensis. This genome had a length of 354,073 bp and the base composition of the genome was A (27.44%), T (27.30%), C (22.58%), and G (22.68%). This genome encodes 59 genes, including 32 protein-coding genes, 24 tRNA genes, and 3 rRNA genes. The Sanger sequencing and Oxford Nanopore sequencing confirmed a pair of direct repeats had mediated genome recombination, resulting in the formation of two conformations. The gene conversation between plastome and mitochondrial genome was also observed in S. tsinyunensis by detecting gene migration, including six tRNA genes (namely, trnW-CCA, trnI-CAU, trnH-UUU, trnD-GUC, trnN-GUU, and trnM-CAU), five protein-coding gene fragments, and the fragments from 2 rRNA genes. Moreover, the dN/dS analysis revealed the atp9 gene had undergone strong negative selection, and four genes (atp4, mttB, ccmFc, and ccmB) probably had undergone positive selection during evolution in Lamiales. This work reported the first mitochondrial genome of S. tsinyunensis, which could be used as a reference genome for the important medicinal plants of the genus Scutellaria, and also provide much-desired information for molecular breeding.

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.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

Towards a plant model for enigmatic U-to-C RNA editing: the organelle genomes, transcriptomes, editomes and candidate RNA editing factors in the hornwort Anthoceros agrestis

Hornworts are crucial to understand the phylogeny of early land plants. The emergence of 'reverse' U-to-C RNA editing accompanying the widespread C-to-U RNA editing in plant chloroplasts and mitochondria may be a molecular synapomorphy of a hornwort-tracheophyte clade. C-to-U RNA editing is well understood after identification of many editing factors in models like Arabidopsis thaliana and Physcomitrella patens, but there is no plant model yet to investigate U-to-C RNA editing. The hornwort Anthoceros agrestis is now emerging as such a model system. We report on the assembly and analyses of the A. agrestis chloroplast and mitochondrial genomes, their transcriptomes and editomes, and a large nuclear gene family encoding pentatricopeptide repeat (PPR) proteins likely acting as RNA editing factors. Both organelles in A. agrestis feature high amounts of RNA editing, with altogether > 1100 sites of C-to-U and 1300 sites of U-to-C editing. The nuclear genome reveals > 1400 genes for PPR proteins with variable carboxyterminal DYW domains. We observe significant variants of the 'classic' DYW domain, in the meantime confirmed as the cytidine deaminase for C-to-U editing, and discuss the first attractive candidates for reverse editing factors given their excellent matches to U-to-C editing targets according to the PPR-RNA binding code.

Comparative Analysis of Four Calypogeia Species Revealed Unexpected Change in Evolutionarily-Stable Liverwort Mitogenomes

Liverwort mitogenomes are considered to be evolutionarily stable. A comparative analysis of four Calypogeia species revealed differences compared to previously sequenced liverwort mitogenomes. Such differences involve unexpected structural changes in the two genes, cox1 and atp1, which have lost three and two introns, respectively. The group I introns in the cox1 gene are proposed to have been lost by two-step localized retroprocessing, whereas one-step retroprocessing could be responsible for the disappearance of the group II introns in the atp1 gene. These cases represent the first identified losses of introns in mitogenomes of leafy liverworts (Jungermanniopsida) contrasting the stability of mitochondrial gene order with certain changes in the gene content and intron set in liverworts.

Group II introns in wheat mitochondria have degenerate structural features and varied splicing pathways

Mitochondrial introns in flowering plant genes are virtually all classified as members of the group II ribozyme family although certain structural features have degenerated to varying degrees over evolutionary time. We are interested in the impact that unconventional intron architecture might have on splicing biochemistry in vivo and we have focused in particular on intronic domains V and VI, which for self-splicing introns provide a key component of the catalytic core and the bulged branchpoint adenosine, respectively. Notably, the two transesterification steps in classical group II splicing are the same as for nuclear spliceosomal introns and release the intron as a lariat. Using RT-PCR and circularized RT-PCR, we had previously demonstrated that several wheat mitochondrial introns which lack a branchpoint adenosine have atypical splicing pathways, and we have now extended this analysis to the full set of wheat introns, namely six trans-splicing and sixteen cis-splicing ones. A number of introns are excised using non-lariat pathways and interestingly, we find that several introns which do have a conventional domain VI also use pathways that appear to exploit other internal or external nucleophiles, with the lariat form being relatively minor. Somewhat surprisingly, several introns with weakly-structured domain V/VI helices still exhibit classical lariat splicing, suggesting that accessory factors aid in restoring a splicing-competent conformation. Our observations illustrate that the loss of conventional group II features during evolution is correlated with altered splicing biochemistry in an intron-distinctive manner.

Complete mitochondrial genomes from the ferns Ophioglossum californicum and Psilotum nudum are highly repetitive with the largest organellar introns

Currently, complete mitochondrial genomes (mitogenomes) are available from all major land plant lineages except ferns. Sequencing of fern mitogenomes could shed light on the major evolutionary transitions that established mitogenomic diversity among extant lineages. In this study, we generated complete mitogenomes from the adder's tongue fern (Ophioglossum californicum) and the whisk fern (Psilotum nudum). The Psilotum mitogenome (628 kb) contains a rich complement of genes and introns, some of which are the largest of any green plant organellar genome. In the Ophioglossum mitogenome (372 kb), gene and intron content is slightly reduced, including the loss of all four mitochondrial ccm genes. Transcripts of nuclear Ccm genes also were not detected, suggesting loss of the entire mitochondrial cytochrome c maturation pathway from Ophioglossum. Both fern mitogenomes are highly repetitive, yet they show extremely low levels of active recombination. Transcriptomic sequencing uncovered ˜1000 sites of C-to-U RNA editing in both species, plus a small number (< 60) of U-to-C edit sites. Overall, the first mitochondrial genomes of ferns show a mix of features shared with lycophytes and/or seed plants and several novel genomic features, enabling a robust reconstruction of the mitogenome in the common ancestor of vascular plants.

Loss of a Trans-Splicing nad1 Intron from Geraniaceae and Transfer of the Maturase Gene matR to the Nucleus in Pelargonium

The mitochondrial nad1 gene of seed plants has a complex structure, including four introns in cis or trans configurations and a maturase gene (matR) hosted within the final intron. In the geranium family (Geraniaceae), however, sequencing of representative species revealed that three of the four introns, including one in a trans configuration and another that hosts matR, were lost from the nad1 gene in their common ancestor. Despite the loss of the host intron, matR has been retained as a freestanding gene in most genera of the family, indicating that this maturase has additional functions beyond the splicing of its host intron. In the common ancestor of Pelargonium, matR was transferred to the nuclear genome, where it was split into two unlinked genes that encode either its reverse transcriptase or maturase domain. Both nuclear genes are transcribed and contain predicted mitochondrial targeting signals, suggesting that they express functional proteins that are imported into mitochondria. The nuclear localization and split domain structure of matR in the Pelargonium nuclear genome offers a unique opportunity to assess the function of these two domains using transgenic approaches.

Analysis of the Complete Mitochondrial Genome Sequence of the Diploid Cotton Gossypium raimondii by Comparative Genomics Approaches

Cotton is one of the most important economic crops and the primary source of natural fiber and is an important protein source for animal feed. The complete nuclear and chloroplast (cp) genome sequences of G. raimondii are already available but not mitochondria. Here, we assembled the complete mitochondrial (mt) DNA sequence of G. raimondii into a circular genome of length of 676,078 bp and performed comparative analyses with other higher plants. The genome contains 39 protein-coding genes, 6 rRNA genes, and 25 tRNA genes. We also identified four larger repeats (63.9 kb, 10.6 kb, 9.1 kb, and 2.5 kb) in this mt genome, which may be active in intramolecular recombination in the evolution of cotton. Strikingly, nearly all of the G. raimondii mt genome has been transferred to nucleus on Chr1, and the transfer event must be very recent. Phylogenetic analysis reveals that G. raimondii, as a member of Malvaceae, is much closer to another cotton (G. barbadense) than other rosids, and the clade formed by two Gossypium species is sister to Brassicales. The G. raimondii mt genome may provide a crucial foundation for evolutionary analysis, molecular biology, and cytoplasmic male sterility in cotton and other higher plants.

Monilophyte mitochondrial rps1 genes carry a unique group II intron that likely originated from an ancient paralog in rpl2

Intron patterns in plant mitochondrial genomes differ significantly between the major land plant clades. We here report on a new, clade-specific group II intron in the rps1 gene of monilophytes (ferns). This intron, rps1i25g2, is strikingly similar to rpl2i846g2 previously identified in the mitochondrial rpl2 gene of seed plants, ferns, and the lycophyte Phlegmariurus squarrosus Although mitochondrial ribosomal protein genes are frequently subject to endosymbiotic gene transfer among plants, we could retrieve the mitochondrial rps1 gene in a taxonomically wide sampling of 44 monilophyte taxa including basal lineages such as the Ophioglossales, Psilotales, and Marattiales with the only exception being the Equisetales (horsetails). Introns rps1i25g2 and rpl2i846g2 were likewise consistently present with only two exceptions: Intron rps1i25g2 is lost in the genus Ophioglossum and intron rpl2i846g2 is lost in Equisetum bogotense Both intron sequences are moderately affected by RNA editing. The unprecedented primary and secondary structure similarity of rps1i25g2 and rpl2i846g2 suggests an ancient retrotransposition event copying rpl2i846g2 into rps1, for which we suggest a model. Our phylogenetic analysis adding the new rps1 locus to a previous data set is fully congruent with recent insights on monilophyte phylogeny and further supports a sister relationship of Gleicheniales and Hymenophyllales.

Convergent Evolution of Fern-Specific Mitochondrial Group II Intron atp1i361g2 and Its Ancient Source Paralogue rps3i249g2 and Independent Losses of Intron and RNA Editing among Pteridaceae

Mitochondrial intron patterns are highly divergent between the major land plant clades. An intron in the atp1 gene, atp1i361g2, is an example for a group II intron specific to monilophytes (ferns). Here, we report that atp1i361g2 is lost independently at least 4 times in the fern family Pteridaceae. Such plant organelle intron losses have previously been found to be accompanied by loss of RNA editing sites in the flanking exon regions as a consequence of genomic recombination of mature cDNA. Instead, we now observe that RNA editing events in both directions of pyrimidine exchange (C-to-U and U-to-C) are retained in atp1 exons after loss of the intron in Pteris argyraea/biaurita and in Actiniopteris and Onychium We find that atp1i361g2 has significant similarity with intron rps3i249g2 present in lycophytes and gymnosperms, which we now also find highly conserved in ferns. We conclude that atp1i361g2 may have originated from the more ancestral rps3i249g2 paralogue by a reverse splicing copy event early in the evolution of monilophytes. Secondary structure elements of the two introns, most characteristically their domains III, show strikingly convergent evolution in the monilophytes. Moreover, the intron paralogue rps3i249g2 reveals relaxed evolution in taxa where the atp1i361g2 paralogue is lost. Our findings may reflect convergent evolution of the two related mitochondrial introns exerted by co-evolution with an intron-binding protein simultaneously acting on the two paralogues.

Organellar phylogenomics of an emerging model system: Sphagnum (peatmoss)

BACKGROUND AND AIMS

Sphagnum-dominated peatlands contain approx. 30 % of the terrestrial carbon pool in the form of partially decomposed plant material (peat), and, as a consequence, Sphagnum is currently a focus of studies on biogeochemistry and control of global climate. Sphagnum species differ in ecologically important traits that scale up to impact ecosystem function, and sequencing of the genome from selected Sphagnum species is currently underway. As an emerging model system, these resources for Sphagnum will facilitate linking nucleotide variation to plant functional traits, and through those traits to ecosystem processes. A solid phylogenetic framework for Sphagnum is crucial to comparative analyses of species-specific traits, but relationships among major clades within Sphagnum have been recalcitrant to resolution because the genus underwent a rapid radiation. Herein a well-supported hypothesis for phylogenetic relationships among major clades within Sphagnum based on organellar genome sequences (plastid, mitochondrial) is provided.

METHODS

We obtained nucleotide sequences (273 753 nucleotides in total) from the two organellar genomes from 38 species (including three outgroups). Phylogenetic analyses were conducted using a variety of methods applied to nucleotide and amino acid sequences. The Sphagnum phylogeny was rooted with sequences from the related Sphagnopsida genera, Eosphagnum and Flatbergium

KEY RESULTS

Phylogenetic analyses of the data converge on the following subgeneric relationships: (Rigida (((Subsecunda) (Cuspidata)) ((Sphagnum) (Acutifolia))). All relationships were strongly supported. Species in the two major clades (i.e. Subsecunda + Cuspidata and Sphagnum + Acutifolia), which include >90 % of all Sphagnum species, differ in ecological niches and these differences correlate with other functional traits that impact biogeochemical cycling. Mitochondrial intron presence/absence are variable among species and genera of the Sphagnopsida. Two new nomenclatural combinations are made, in the genera Eosphagnum and Flatbergium

CONCLUSIONS

Newly resolved relationships now permit phylogenetic analyses of morphological, biochemical and ecological traits among Sphagnum species. The results clarify long-standing disagreements about subgeneric relationships and intrageneric classification.

Multiple splicing pathways of group II trans-splicing introns in wheat mitochondria

Trans-splicing of discontinuous introns in plant mitochondria requires the assembly of independently-transcribed precursor RNAs into splicing-competent structures, and they are expected to be excised as Y-branched molecules ("broken lariats") because these introns belong to the group II ribozyme family. We now demonstrate that this is just one of several trans-splicing pathways for wheat mitochondrial nad1 intron 4 and nad5 intron 2; they also use a hydrolytic pathway and the liberated 5'-half-intron linear molecules are unexpectedly abundant in the RNA population. We also observe a third productive splicing pathway for nad5 intron 2 that yields full-length excised introns in which the termini are joined in vivo and possess non-encoded nucleotides. In the case of trans-splicing nad1 intron 1, which has a weakly-structured and poorly-conserved core sequence, excision appears to be solely through a hydrolytic pathway. When wheat embryos are germinated in the cold rather than at room temperature, an increased complexity in trans-splicing products is seen for nad1 intron 4, suggesting that there can be environmental effects on the RNA folding of bipartite introns. Our observations provide insights into intron evolution and the complexity of RNA processing events in plant mitochondria.

Dynamic evolution of Geranium mitochondrial genomes through multiple horizontal and intracellular gene transfers

The exchange of genetic material between cellular organelles through intracellular gene transfer (IGT) or between species by horizontal gene transfer (HGT) has played an important role in plant mitochondrial genome evolution. The mitochondrial genomes of Geraniaceae display a number of unusual phenomena including highly accelerated rates of synonymous substitutions, extensive gene loss and reduction in RNA editing. Mitochondrial DNA sequences assembled for 17 species of Geranium revealed substantial reduction in gene and intron content relative to the ancestor of the Geranium lineage. Comparative analyses of nuclear transcriptome data suggest that a number of these sequences have been functionally relocated to the nucleus via IGT. Evidence for rampant HGT was detected in several Geranium species containing foreign organellar DNA from diverse eudicots, including many transfers from parasitic plants. One lineage has experienced multiple, independent HGT episodes, many of which occurred within the past 5.5 Myr. Both duplicative and recapture HGT were documented in Geranium lineages. The mitochondrial genome of Geranium brycei contains at least four independent HGT tracts that are absent in its nearest relative. Furthermore, G. brycei mitochondria carry two copies of the cox1 gene that differ in intron content, providing insight into contrasting hypotheses on cox1 intron evolution.

Chloroplast RNA editing going extreme: more than 3400 events of C-to-U editing in the chloroplast transcriptome of the lycophyte Selaginella uncinata

RNA editing in chloroplasts and mitochondria of land plants differs significantly in abundance. For example, some 200-500 sites of cytidine-to-uridine RNA editing exist in flowering plant mitochondria as opposed to only 30-50 such C-to-U editing events in their chloroplasts. In contrast, we predicted significantly more chloroplast RNA editing for the protein-coding genes in the available complete plastome sequences of two species of the spike moss genus Selaginella (Lycopodiophyta). To evaluate these predictions we investigated the Selaginella uncinata chloroplast transcriptome. Our exhaustive cDNA studies identified the extraordinary number of 3415 RNA-editing events, exclusively of the C-to-U type, in the 74 mRNAs encoding intact reading frames in the S. uncinata chloroplast. We find the overwhelming majority (61%) of the 428 silent editing events leaving codon meanings unaltered directly neighboring other editing events, possibly suggesting a sterically more flexible RNA-editing deaminase activity in Selaginella. No evidence of RNA editing was found for tRNAs or rRNAs but we identified a total of 74 editing sites in cDNA sequences of four group II introns (petBi6g2, petDi8g2, ycf3i124g2, and ycf3i354g2) retained in partially matured transcripts, which strongly contribute to improved base-pairing in the intron secondary structures as a likely prerequisite for their splicing.

Complete sequences of organelle genomes from the medicinal plant Rhazya stricta (Apocynaceae) and contrasting patterns of mitochondrial genome evolution across asterids

Background

Rhazya stricta is native to arid regions in South Asia and the Middle East and is used extensively in folk medicine to treat a wide range of diseases. In addition to generating genomic resources for this medicinally important plant, analyses of the complete plastid and mitochondrial genomes and a nuclear transcriptome from Rhazya provide insights into inter-compartmental transfers between genomes and the patterns of evolution among eight asterid mitochondrial genomes.

Results

The 154,841 bp plastid genome is highly conserved with gene content and order identical to the ancestral organization of angiosperms. The 548,608 bp mitochondrial genome exhibits a number of phenomena including the presence of recombinogenic repeats that generate a multipartite organization, transferred DNA from the plastid and nuclear genomes, and bidirectional DNA transfers between the mitochondrion and the nucleus. The mitochondrial genes sdh3 and rps14 have been transferred to the nucleus and have acquired targeting presequences. In the case of rps14, two copies are present in the nucleus; only one has a mitochondrial targeting presequence and may be functional. Phylogenetic analyses of both nuclear and mitochondrial copies of rps14 across angiosperms suggests Rhazya has experienced a single transfer of this gene to the nucleus, followed by a duplication event. Furthermore, the phylogenetic distribution of gene losses and the high level of sequence divergence in targeting presequences suggest multiple, independent transfers of both sdh3 and rps14 across asterids. Comparative analyses of mitochondrial genomes of eight sequenced asterids indicates a complicated evolutionary history in this large angiosperm clade with considerable diversity in genome organization and size, repeat, gene and intron content, and amount of foreign DNA from the plastid and nuclear genomes.

Conclusions

Organelle genomes of Rhazya stricta provide valuable information for improving the understanding of mitochondrial genome evolution among angiosperms. The genomic data have enabled a rigorous examination of the gene transfer events. Rhazya is unique among the eight sequenced asterids in the types of events that have shaped the evolution of its mitochondrial genome. Furthermore, the organelle genomes of R. stricta provide valuable genomic resources for utilizing this important medicinal plant in biotechnology applications.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-405) contains supplementary material, which is available to authorized users.

nMAT4, a maturase factor required for nad1 pre-mRNA processing and maturation, is essential for holocomplex I biogenesis in Arabidopsis mitochondria

Group II introns are large catalytic RNAs that are found in bacteria and organellar genomes of lower eukaryotes, but are particularly prevalent within mitochondria in plants, where they are present in many critical genes. The excision of plant mitochondrial introns is essential for respiratory functions, and is facilitated in vivo by various protein cofactors. Typical group II introns are classified as mobile genetic elements, consisting of the self-splicing ribozyme and its own intron-encoded maturase protein. A hallmark of maturases is that they are intron-specific, acting as cofactors that bind their intron-containing pre-RNAs to facilitate splicing. However, the degeneracy of the mitochondrial introns in plants and the absence of cognate intron-encoded maturase open reading frames suggest that their splicing in vivo is assisted by 'trans'-acting protein factors. Interestingly, angiosperms harbor several nuclear-encoded maturase-related (nMat) genes that contain N-terminal mitochondrial localization signals. Recently, we established the roles of two of these paralogs in Arabidopsis, nMAT1 and nMAT2, in the splicing of mitochondrial introns. Here we show that nMAT4 (At1g74350) is required for RNA processing and maturation of nad1 introns 1, 3 and 4 in Arabidopsis mitochondria. Seed germination, seedling establishment and development are strongly affected in homozygous nmat4 mutants, which also show modified respiration phenotypes that are tightly associated with complex I defects.

Impact of low temperature on splicing of atypical group II introns in wheat mitochondria

To investigate the impact of cold on group II intron splicing, we compared the physical forms of excised mitochondrial introns from wheat embryos germinated at room temperature and 4°C. For introns which deviate from the conventional branchpoint structure, we observed predominantly heterogeneous circularized introns in the cold rather than linear polyadenylated forms arising from a hydrolytic pathway as seen at room temperature. In addition, intron-containing precursors are elevated relative to mature mRNAs upon cold treatment. Our findings indicate that low temperature growth not only reduces splicing efficiency, but also shifts the splicing biochemistry of atypical group II introns to novel, yet productive, pathways.

Complete plastid genomes from Ophioglossum californicum, Psilotum nudum, and Equisetum hyemale reveal an ancestral land plant genome structure and resolve the position of Equisetales among monilophytes

BACKGROUND

Plastid genome structure and content is remarkably conserved in land plants. This widespread conservation has facilitated taxon-rich phylogenetic analyses that have resolved organismal relationships among many land plant groups. However, the relationships among major fern lineages, especially the placement of Equisetales, remain enigmatic.

RESULTS

In order to understand the evolution of plastid genomes and to establish phylogenetic relationships among ferns, we sequenced the plastid genomes from three early diverging species: Equisetum hyemale (Equisetales), Ophioglossum californicum (Ophioglossales), and Psilotum nudum (Psilotales). A comparison of fern plastid genomes showed that some lineages have retained inverted repeat (IR) boundaries originating from the common ancestor of land plants, while other lineages have experienced multiple IR changes including expansions and inversions. Genome content has remained stable throughout ferns, except for a few lineage-specific losses of genes and introns. Notably, the losses of the rps16 gene and the rps12i346 intron are shared among Psilotales, Ophioglossales, and Equisetales, while the gain of a mitochondrial atp1 intron is shared between Marattiales and Polypodiopsida. These genomic structural changes support the placement of Equisetales as sister to Ophioglossales + Psilotales and Marattiales as sister to Polypodiopsida. This result is augmented by some molecular phylogenetic analyses that recover the same relationships, whereas others suggest a relationship between Equisetales and Polypodiopsida.

CONCLUSIONS

Although molecular analyses were inconsistent with respect to the position of Marattiales and Equisetales, several genomic structural changes have for the first time provided a clear placement of these lineages within the ferns. These results further demonstrate the power of using rare genomic structural changes in cases where molecular data fail to provide strong phylogenetic resolution.

Loss of two introns from the Magnolia tripetala mitochondrial cox2 gene implicates horizontal gene transfer and gene conversion as a novel mechanism of intron loss

Intron loss is often thought to occur through retroprocessing, which is the reverse transcription and genomic integration of a spliced transcript. In plant mitochondria, several unambiguous examples of retroprocessing are supported by the parallel loss of an intron and numerous adjacent RNA edit sites, but in most cases, the evidence for intron loss via retroprocessing is weak or lacking entirely. To evaluate mechanisms of intron loss, we designed a polymerase chain reaction (PCR)-based assay to detect recent intron losses from the mitochondrial cox2 gene within genus Magnolia, which was previously suggested to have variability in cox2 intron content. Our assay showed that all 22 examined species have a cox2 gene with two introns. However, one species, Magnolia tripetala, contains an additional cox2 gene that lacks both introns. Quantitative PCR showed that both M. tripetala cox2 genes are present in the mitochondrial genome. Although the intronless gene has lost several ancestral RNA edit sites, their distribution is inconsistent with retroprocessing models. Instead, phylogenetic and gene conversion analyses indicate that the intronless gene was horizontally acquired from a eudicot and then underwent gene conversion with the native intron-containing gene. The models are presented to summarize the roles of horizontal gene transfer and gene conversion as a novel mechanism of intron loss.

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.

Genes and processed paralogs co-exist in plant mitochondria

RNA-mediated gene duplication has been proposed to create processed paralogs in the plant mitochondrial genome. A processed paralog may retain signatures left by the maturation process of its RNA precursor, such as intron removal and no need of RNA editing. Whereas it is well documented that an RNA intermediary is involved in the transfer of mitochondrial genes to the nucleus, no direct evidence exists for insertion of processed paralogs in the mitochondria (i.e., processed and un-processed genes have never been found simultaneously in the mitochondrial genome). In this study, we sequenced a region of the mitochondrial gene nad1, and identified a number of taxa were two different copies of the region co-occur in the mitochondria. The two nad1 paralogs differed in their (a) presence or absence of a group II intron, and (b) number of edited sites. Thus, this work provides the first evidence of co-existence of processed paralogs and their precursors within the plant mitochondrial genome. In addition, mapping the presence/absence of the paralogs provides indirect evidence of RNA-mediated gene duplication as an essential process shaping the mitochondrial genome in plants.

Group I intron-mediated trans-splicing in mitochondria of Gigaspora rosea and a robust phylogenetic affiliation of arbuscular mycorrhizal fungi with Mortierellales

Gigaspora rosea is a member of the arbuscular mycorrhizal fungi (AMF; Glomeromycota) and a distant relative of Glomus species that are beneficial to plant growth. To allow for a better understanding of Glomeromycota, we have sequenced the mitochondrial DNA of G. rosea. A comparison with Glomus mitochondrial genomes reveals that Glomeromycota undergo insertion and loss of mitochondrial plasmid-related sequences and exhibit considerable variation in introns. The gene order between the two species is almost completely reshuffled. Furthermore, Gigaspora has fragmented cox1 and rns genes, and an unorthodox initiator tRNA that is tailored to decoding frequent UUG initiation codons. For the fragmented cox1 gene, we provide evidence that its RNA is joined via group I-mediated trans-splicing, whereas rns RNA remains in pieces. According to our model, the two cox1 precursor RNA pieces are brought together by flanking cox1 exon sequences that form a group I intron structure, potentially in conjunction with the nad5 intron 3 sequence. Finally, we present analyses that address the controversial phylogenetic association of Glomeromycota within fungi. According to our results, Glomeromycota are not a separate group of paraphyletic zygomycetes but branch together with Mortierellales, potentially also Harpellales.

Exclusive conservation of mitochondrial group II intron nad4i548 among liverworts and its use for phylogenetic studies in this ancient plant clade

Liverworts occupy a pivotal position in land plant (embryophyte) phylogeny as the presumed earliest-branching major clade, sister to all other land plants, including the mosses, hornworts, lycophytes, monilophytes and seed plants. Molecular support for this earliest dichotomy in land plant phylogeny comes from strikingly different occurrences of introns in mitochondrial genes distinguishing liverworts from all other embryophytes. Exceptionally, however, the nad5 gene--the mitochondrial locus hitherto used most widely to elucidate early land plant phylogeny--carries a group I type intron that is shared between liverworts and mosses. We here explored whether a group II intron, the other major type of organellar intron, would similarly be conserved in position across the entire diversity of extant liverworts and could be of use for phylogenetic analyses in this supposedly most ancient embryophyte clade. To this end, we investigated the nad4 gene as a candidate locus possibly featuring different introns in liverworts as opposed to the non-liverwort embryophyte (NLE) lineage. We indeed found group II intron nad4i548 universally conserved in a wide phylogenetic sampling of 55 liverwort taxa, confirming clade specificity and surprising evolutionary stability of plant mitochondrial introns. As expected, intron nad4i548g2 carries phylogenetic information in its variable sequences, which confirms and extends previous cladistic insights on liverwort evolution. We integrate the new nad4 data with those of the previously established mitochondrial nad5 and the chloroplast rbcL and rps4 genes and present a phylogeny based on the fused datasets. Notably, the phylogenetic analyses suggest a reconsideration of previous phylogenetic and taxonomic assignments for the genera Calycularia and Mylia and resolve a sister group relationship of Ptilidiales and Porellales.

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.

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.

Purification, characterization, and molecular cloning of a novel antifungal lectin from the roots of Ophioglossum pedunculosum

A novel mannan-specific lectin was isolated from the roots of a traditional Chinese herbal medicine, Ophioglossum pedunculosum through ion-exchange chromatography and gel filtration. With a molecular mass of 19,835.7 Da demonstrated by MALDI-TOF analysis, this novel agglutinin was designated as O. pedunculosum agglutinin (OPA), specifically agglutinating human O erythrocytes and rabbit erythrocytes. The hemagglutination could be strongly inhibited by mannan and thyroglobulin, the activity of which was stable in pH range of 4.0-8.0 and at temperatures below 50 °C. Chemical modification studies indicated that tryptophan and arginine residues were essential for its hemagglutinating activity. Meanwhile, it showed antifungal activities toward Sclerotium rolfsii and Fusarium graminearum. In addition, to amplify cDNA of OPA by 3'/5'-rapid amplification of cDNA ends (RACE), the N-terminal 30 amino acids sequence of OPA was determined, and degenerate primers were designed. The obtained full-length cDNA of OPA contained 885 bp with an open-reading frame of 600 bp encoding a precursor protein of 199 amino acids, while the mature protein had 170 amino acids.

Extreme RNA editing in coding islands and abundant microsatellites in repeat sequences of Selaginella moellendorffii mitochondria: the root of frequent plant mtDNA recombination in early tracheophytes

Using an independent fosmid cloning approach and comprehensive transcriptome analysis to complement data from the Selaginella moellendorffii genome project, we determined the complete mitochondrial genome structure of this spikemoss. Numerous recombination events mediated mainly via long sequence repeats extending up to 7kbp result in a complex mtDNA network structure. Peculiar features associated with the repeat sequences are more than 80 different microsatellite sites (predominantly trinucleotide motifs). The S. moellendorffii mtDNA encodes a plant-typical core set of a twin-arginine translocase (tatC), 17 respiratory chain subunits, and 2 rRNAs but lacks atp4 and any tRNA genes. As a further novelty among plant chondromes, the nad4L gene is encoded within an intron of the nad1 gene. A total of 37 introns occupying the 20 mitochondrial genes (four of which are disrupted into trans-splicing arrangements including two novel instances of trans-splicing introns) make the S. moellendorffii chondrome the intron-richest and gene-poorest plant mtDNA known. Our parallel transcriptome analyses demonstrates functional splicing of all 37 introns and reveals a new record amount of plant organelle RNA editing with a total of 2,139 sites in mRNAs and 13 sites in the two rRNAs, all of which are exclusively of the C-to-U type.

A multi-locus molecular phylogeny of the Lepidoziaceae: laying the foundations for a stable classification

The Lepidoziaceae, with over 700 species in 30 genera, is one of the largest leafy liverwort families. Despite receiving considerable attention, the composition of subfamilies and genera remains unsatisfactorily resolved. In this study, 10 loci (one nuclear 26S, two mitochondrial nad1 and rps3, and seven chloroplast atpB, psbA, psbT-psbH, rbcL, rps4, trnG and trnL-trnF) are used to estimate the phylogeny of 93 species of Lepidoziaceae. These molecular data provide strong evidence against the monophyly of three subfamilies; Lepidozioideae, Lembidioideae and Zoopsidoideae, and seven of the 20 sampled genera; Lepidozia, Telaranea, Kurzia, Zoopsis, Lembidium, Paracromastigum and Chloranthelia. Several robust clades are recognised that might provide the basis for a revised subfamily circumscription including a narrower circumscription of the Lepidozioideae and a more inclusive Lembidioideae. Neogrollea notabilis is returned to the Lepidoziaceae and Megalembidium insulanum is placed in the Lembidioideae.

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.

Introducing intron locus cox1i624 for phylogenetic analyses in Bryophytes: on the issue of Takakia as sister genus to all other extant mosses

Liverworts are well supported as the sister group to all other land plants (embryophytes) by molecular data. Observations strongly supporting this earliest dichotomy in embryophyte evolution are the strikingly different introns occurring in the mitochondrial DNAs of liverworts versus non-liverwort embryophytes (NLE), including the mosses. A final conclusion on the most basal lineages of mosses, for which genera such as Sphagnum and Takakia are the most likely candidates, is lacking. We have now investigated cox1i624, a mitochondrial group I intron conserved between the moss Physcomitrella patens and the liverwort Marchantia polymorpha. Focusing on a sampling of liverwort and moss genera, which had previously been identified as early branching taxa in their respective clades, we find that group I intron cox1i624 is universally conserved in all 33 mosses and 11 liverworts investigated. The group I intron core secondary structure is well conserved between the two ancient land plant clades. However, whereas dramatic size reductions are seen in the moss phylogeny, exactly the opposite is observed for liverworts. The cox1i624g1 locus was used for phylogenetic tree reconstruction also in combination with data sets of nad5i753g1 as well as chloroplast loci rbcL and rps4. The phylogenetic analyses revealed (i) very good support for the Treubiopsida as sister clade to all other liverworts, (ii) a sister group relationship of the nematodontous Tetraphidopsida and Polytrichopsida and (iii) two rivalling hypotheses about the basal-most moss genus with mitochondrial loci suggesting an isolated Takakia as sister to all other mosses and chloroplast loci indicating a Takakia-Sphagnum clade.

Extensive mis-splicing of a bi-partite plant mitochondrial group II intron

Expression of the seed plant mitochondrial nad5 gene involves two trans-splicing events that remove fragmented group II introns and join the small, central exon c to exons b and d. We show that in both monocot and eudicot plants, extensive mis-splicing of the bi-partite intron 2 takes place, resulting in the formation of aberrantly spliced products in which exon c is joined to various sites within exon b. These mis-spliced products accumulate to levels comparable to or greater than that of the correctly spliced mRNA. We suggest that mis-splicing may result from folding constraints imposed on intron 2 by base-pairing between exon a and a portion of the bi-partite intron 3 downstream of exon c. Consistent with this hypothesis, we find that mis-splicing does not occur in Oenothera mitochondria, where intron 3 is further fragmented such that the predicted base-pairing region is not covalently linked to exon c. Our findings suggest that intron fragmentation may lead to mis-splicing, which may be corrected by further intron fragmentation.

The complete mitochondrial genome sequence of the liverwort Pleurozia purpurea reveals extremely conservative mitochondrial genome evolution in liverworts

Plant mitochondrial genomes have been known to be highly unusual in their large sizes, frequent intra-genomic rearrangement, and generally conservative sequence evolution. Recent studies show that in early land plants the mitochondrial genomes exhibit a mixed mode of conservative yet dynamic evolution. Here, we report the completely sequenced mitochondrial genome from the liverwort Pleurozia purpurea. The circular genome has a size of 168,526 base pairs, containing 43 protein-coding genes, 3 rRNA genes, 25 tRNA genes, and 31 group I or II introns. It differs from the Marchantia polymorpha mitochondrial genome, the only other liverwort chondriome that has been sequenced, in lacking two genes (trnRucg and trnTggu) and one intron (rrn18i1065gII). The two genomes have identical gene orders and highly similar sequences in exons, introns, and intergenic spacers. Finally, a comparative analysis of duplicated trnRucu and other trnR genes from the two liverworts and several other organisms identified the recent lateral origin of trnRucg in Marchantia mtDNA through modification of a duplicated trnRucu. This study shows that the mitochondrial genomes evolve extremely slowly in liverworts, the earliest-diverging lineage of extant land plants, in stark contrast to what is known of highly dynamic evolution of mitochondrial genomes in seed plants.

A trans-splicing group I intron and tRNA-hyperediting in the mitochondrial genome of the lycophyte Isoetes engelmannii

Plant mitochondrial genomes show much more evolutionary plasticity than those of animals. We analysed the first mitochondrial DNA (mtDNA) of a lycophyte, the quillwort Isoetes engelmannii, which is separated from seed plants by more than 350 million years of evolution. The Isoetes mtDNA is particularly rich in recombination events, and chloroplast as well as nuclear DNA inserts document the incorporation of foreign sequences already in this most ancestral vascular plant lineage. On the other hand, particularly small group II introns and short intergenic regions reveal a tendency of evolution towards a compact mitochondrial genome. RNA editing reaches extreme levels exceeding 100 pyrimidine exchanges in individual mRNAs and, hitherto unobserved in such frequency, also in tRNAs with 18 C-to-U conversions in the tRNA for proline. In total, some 1500 sites of RNA editing can be expected for the Isoetes mitochondrial transcriptome. As a unique molecular novelty, the Isoetes cox1 gene requires trans-splicing via a discontinuous group I intron demonstrating disrupted, but functional, RNAs for yet another class of natural ribozymes.