The process of RNA editing in plant mitochondria

Emerging roles of the C-to-U RNA editing in plant stress responses

RNA editing is an important post-transcriptional event in all living cells. Within chloroplasts and mitochondria of higher plants, RNA editing involves the deamination of specific cytosine (C) residues in precursor RNAs to uracil (U). An increasing number of recent studies detail specificity of C-to-U RNA editing as an essential prerequisite for several plant stress-related responses. In this review, we summarize the current understanding of responses and functions of C-to-U RNA editing in plants under various stress conditions to provide theoretical reference for future research.

Analysis of the Rhodomyrtus tomentosa mitochondrial genome: Insights into repeat-mediated recombination and intra-cellular DNA transfer

Plant mitochondrial genomes participate in encoding proteins crucial to the major producers of ATP in the cell and replication and heredity of their own DNA. The sequences and structure of the plant mitochondrial genomes profoundly impact these fundamental processes, and studies of plant mitochondrial genomes are needed. We reported the complete sequences of the Rhodomyrtus tomentosa mitochondrial genome here, totaling 400,482 bp. Nanopore ONT reads and PCR amplification provided evidence for recombination mediated by the eight repeat pairs for the R. tomentosa mitochondrial genome. Thirty-eight genes were identified in the R. tomentosa mitochondrial genome. Comparative analyses of the mitochondrial genome and plastome and PCR amplification suggest that five fragments of mitochondrial plastid DNA were unfunctional sequences resulting from intracellular gene transfer. Phylogenetic analysis based on each and all of the 27 mitochondrial protein-coding genes of nine Myrtales species revealed that R. tomentosa always clustered with other species of Myrtaceae. This study uncovered the enormous complexity of the R. tomentosa mitochondrial genome, the active repeat-mediated recombinations, the presence of mitochondrial plastid DNAs, and the topological incongruence of Myrtales among the single-gene trees.

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

The first mitochondrial genome of Calophyllum soulattri Burm.f

Calophyllum soulattri Burm.f. is traditionally used to treat skin infections and reduce rheumatic pain, yet genetic and genomic studies are still limited. Here, we present the first complete mitochondrial genome of C. soulattri. It is 378,262 bp long with 43.97% GC content, containing 55 genes (30 protein-coding, 5 rRNA, and 20 tRNA). Repeat analysis of the mitochondrial genome revealed 194 SSRs, mostly mononucleotides, and 266 pairs of dispersed repeats ( ≥ 30 bp) that were predominantly palindromic. There were 23 homologous fragments found between the mitochondrial and plastome genomes. We also predicted 345 C-to-U RNA editing sites from 30 protein-coding genes (PCGs) of the C. soulatrii mitochondrial genome. These RNA editing events created the start codon of nad1 and the stop codon of ccmFc. Most PCGs of the C. soulattri mitochondrial genome underwent negative selection, but atp4 and ccmB experienced positive selection. Phylogenetic analyses showed C. soulattri is a sister taxon of Garcinia mangostana. This study has shed light on C. soulattri's evolution and Malpighiales' phylogeny. As the first complete mitochondrial genome in Calophyllaceae, it can be used as a reference genome for other medicinal plant species within the family for future genetic studies.

Assembly and comparative analysis of the complete mitochondrial genome of Pinellia ternata

Pinellia ternata is an important natural medicinal herb in China. However, it is susceptible to withering when exposed to high temperatures during growth, which limits its tuber production. Mitochondria usually function in stress response. The P . ternata mitochondrial (mt) genome has yet to be explored. Therefore, we integrated PacBio and Illumina sequencing reads to assemble and annotate the mt genome of P . ternata . The circular mt genome of P . ternata is 876 608bp in length and contains 38 protein-coding genes (PCGs), 20 tRNA genes and three rRNA genes. Codon usage, sequence repeats, RNA editing and gene migration from chloroplast (cp) to mt were also examined. Phylogenetic analysis based on the mt genomes of P . ternata and 36 other taxa revealed the taxonomic and evolutionary status of P . ternata . Furthermore, we investigated the mt genome size and GC content by comparing P . ternata with the other 35 species. An evaluation of non-synonymous substitutions and synonymous substitutions indicated that most PCGs in the mt genome underwent negative selection. Our results provide comprehensive information on the P . ternata mt genome, which may facilitate future research on the high-temperature response of P . ternata and provide new molecular insights on the Araceae family.

Withania somnifera mitochondrial atp4 gene editing alters the ATP synthase b subunit, independent of salt stress

Numerous studies have shown that stress in plant cells and organelles with transport electron chains is related to RNA editing. The ATP synthase complex present in mitochondria plays a crucial role in cellular respiration and consists of several subunits. Among them is the b subunit, which is encoded by the mitochondrial atp4 gene. Computing-based analysis of the effects of RNA editing of the Withania somnifera atp4 gene in mitochondria leading to alterations in the b subunit of ATP synthase. Using the CLC Genomic Workbench 3, RNA editing analysis between the control and salt stress conditions was not significantly different. Depending on RNA editing, the tertiary structure model revealed a change in the states of the b subunit, reflecting differences in the central stalk and F1-catalytic domain. The study found that polar edits in the N-terminus of the b subunit allow for efficient H + ion selectivity and introduce a new coiled-coil alpha-helical structure that may help stabilize the complex. The most noteworthy finding of this study was the strong impact of these editing events on the tertiary structure of the b subunit, which has the potential to affect the ATPase activity and indicate that the editing in this subunit aimed to restore the original active protein and not as a response to salt stress.

Effect of salinity on ccmfn gene RNA editing of mitochondria in wild barley and uncommon types of RNA editing

The primary function of mitochondria is cellular respiration and energy production. Cytochrome C complex is an essential complex that transports electrons in the respiratory chain between complex III and complex IV. One of this complex's main subunits is CcmFN, which is believed to be crucial for holocytochrome assembly. In wild-type plant Hordeum vulgare subsp. spontaneum, four ccmfn cDNAs are subjected to high salt stress (500 mM salinity), 0 h (or control) (GenBank accession no. ON764850), after 2 h (GenBank accession no. ON7648515), after 12 h (GenBank accession no. ON764852), and after 24 h (GenBank accession no. ON764853) and mtDNA of ccmfn gene (GenBank accession no. ON764854). Using raw data from RNA-seq, 47 sites with nucleotide and amino acid modifications were detected. There were ten different RNA editing types, with most of them are C to U. Unusual editing types in plants have also been found, such as A to C, C to A, A to G, A to U, T to A, T to C, C to G, G to C, and T to G. High levels of editing were observed in control as well as treatments of salinity stress. Amino acid changes were found in 43 sites; nearly all showed hydrophilic to hydrophilic alterations. Only C749 showed regulation under salinity stress.

C-to-U RNA Editing: A Site Directed RNA Editing Tool for Restoration of Genetic Code

The restoration of genetic code by editing mutated genes is a potential method for the treatment of genetic diseases/disorders. Genetic disorders are caused by the point mutations of thymine (T) to cytidine (C) or guanosine (G) to adenine (A), for which gene editing (editing of mutated genes) is a promising therapeutic technique. In C-to-Uridine (U) RNA editing, it converts the base C-to-U in RNA molecules and leads to nonsynonymous changes when occurring in coding regions; however, for G-to-A mutations, A-to-I editing occurs. Editing of C-to-U is not as physiologically common as that of A-to-I editing. Although hundreds to thousands of coding sites have been found to be C-to-U edited or editable in humans, the biological significance of this phenomenon remains elusive. In this review, we have tried to provide detailed information on physiological and artificial approaches for C-to-U RNA editing.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Drought Response in Wheat

RNA/DNA difference (RDD) is a post-transcriptional RNA modification to enrich genetic information, widely involved in regulating diverse biological processes in eukaryotes. RDDs in the wheat nuclear genome, especially those associated with drought response or tolerance, were not well studied up to now. In this study, we investigated the RDDs related to drought response based on the RNA-seq data of drought-stressed and control samples in wheat. In total, 21,782 unique RDDs were identified, of which 265 were found to be drought-induced, representing the first drought-responsive RDD landscape in the wheat nuclear genome. The drought-responsive RDDs were located in 69 genes, of which 35 were differentially expressed under drought stress. Furthermore, the effects of RNA/DNA differences were investigated, showing that they could result in changes of RNA secondary structure, miRNA-target binding as well as protein conserved domains in the RDD-containing genes. In particular, the A to C mutation in TraesCS2A02G053100 (orthology to OsRLCK) led to the loss of tae-miR9657b-5p targeting, indicating that RNA/DNA difference might mediate miRNA to regulate the drought-response process. This study reported the first drought-responsive RDDs in the wheat nuclear genome. It sheds light on the roles of RDD in drought tolerance, and may also contribute to wheat genetic improvement based on epi-transcriptome methods.

OTP970 Is Required for RNA Editing of Chloroplast ndhB Transcripts in Arabidopsis thaliana

RNA editing is essential for compensating for defects or mutations in haploid organelle genomes and is regulated by numerous trans-factors. Pentatricopeptide repeat (PPR) proteins are the prime factors that are involved in RNA editing; however, many have not yet been identified. Here, we screened the plastid-targeted PLS-DYW subfamily of PPR proteins belonging to Arabidopsis thaliana and identified ORGANELLE TRANSCRIPT PROCESSING 970 (OTP970) as a key player in RNA editing in plastids. A loss-of-function otp970 mutant was impaired in RNA editing of ndhB transcripts at site 149 (ndhB-C149). RNA-immunoprecipitation analysis indicated that OTP970 was associated with the ndhB-C149 site. The complementation of the otp970 mutant with OTP970 lacking the DYW domain (OTP970^∆DYW) failed to restore the RNA editing of ndhB-C149. ndhB gene encodes the B subunit of the NADH dehydrogenase-like (NDH) complex; however, neither NDH activity and stability nor NDH-PSI supercomplex formation were affected in otp970 mutant compared to the wild type, indicating that alteration in amino acid sequence is not necessary for NdhB function. Together, these results suggest that OTP970 is involved in the RNA editing of ndhB-C149 and that the DYW domain is essential for its function.

The first report of RNA U to C or G editing in the mitochondrial NADH dehydrogenase subunit 5 (Nad5) transcript of wild barley

BACKGROUND

Nad dehydrogenase complex in mtDNA has a significant role in cellular respiration. One of the largest subunits in the complex is subunit 5 (Nad5).

METHODS AND RESULTS

Four cDNAs of the Hordeum vulgare subsp. spontaneum nad5 gene have been characterized and subjected to four phases of 0.5 M salinity, at 0 h (control, accession no. MT235236), after 2 h (acc. no. MT235237), after 12 h (acc. no. MT235238) and after 24 h (acc. no. MT235239). Utilizing raw data from RNA-seq, ten RNA editing sites were reported. Seven sites have common editing from C to U in positions (C1490, C1859, C1895, C1900, C1901, C1916, C1918). A rare editing event U to C was detected in two positions (U1650 and U1652) and a novel editing event U to G was for the first time in positions nad5-U231. The highest editing level was shown in 2 and 12 h after salinity exposure. After 24 h, these edits were disrupted, possibly due to the launch of the programed cell death mechanism. However, the RNA editing in positions U1650, U1652 and U231 was fixed at all exposure times.

CONCLUSIONS

Although study clarified the role of salinity stress in nad5 RNA editing sites, the main achievements are first report of U to G RNA editing in plants at position U231 and first report of U to C editing in the nad5 gene at U1650 and U1652.

Differential RNA Editing and Intron Splicing in Soybean Mitochondria during Nodulation

Nitrogen fixation in soybean consumes a tremendous amount of energy, leading to substantial differences in energy metabolism and mitochondrial activities between nodules and uninoculated roots. While C-to-U RNA editing and intron splicing of mitochondrial transcripts are common in plant species, their roles in relation to nodule functions are still elusive. In this study, we performed RNA-seq to compare transcript profiles and RNA editing of mitochondrial genes in soybean nodules and roots. A total of 631 RNA editing sites were identified on mitochondrial transcripts, with 12% or 74 sites differentially edited among the transcripts isolated from nodules, stripped roots, and uninoculated roots. Eight out of these 74 differentially edited sites are located on the matR transcript, of which the degrees of RNA editing were the highest in the nodule sample. The degree of mitochondrial intron splicing was also examined. The splicing efficiencies of several introns in nodules and stripped roots were higher than in uninoculated roots. These include nad1 introns 2/3/4, nad4 intron 3, nad5 introns 2/3, cox2 intron 1, and ccmFc intron 1. A greater splicing efficiency of nad4 intron 1, a higher NAD4 protein abundance, and a reduction in supercomplex I + III2 were also observed in nodules, although the causal relationship between these observations requires further investigation.

Light/heat effects on RNA editing in chloroplast NADH-plastoquinone oxidoreductase subunit 2 (ndhB) gene of Calotropis (Calotropis procera)

Background

RNA editing is common in terrestrial plants, especially in mitochondria and chloroplast. In the photosynthesis process, NAD dehydrogenase plays a very important role. Subunit 2 of NADH-dehydrogenase is one of the major subunits in NAD dehydrogenase complex. Using desert plant Calotropis (Calotropis procera), this study focuses on the RNA editing activity of ndhB based on light time.

Results

NdhB (NADH-dehydrogenase subunit 2) gene accession no. MK144329 was isolated from Calotropis procera genomic data (PRJNA292713). Additionally, using RNA-seq data, the cDNA of the ndhB gene of C. procera was isolated at three daylight periods, i.e., dawn (accession no. MK165161), at midday (accession no. MK165160), and pre-dusk (accession no. MK165159). Seven RNA editing sites have been found in several different positions (nucleotide no. C467, C586, C611, C737, C746, C830, and C1481) within the ndhB coding region. The rate of these alterations was deferentially edited across the three daylight periods. RNA editing rate of ndhB gene was highest at dawn, (87.5, 79.6, 78.5, 76, 68.6, 39.3, and 96.9%, respectively), less in midday (74.8, 54.1, 62.6, 47.4, 45.5, 47.4, and 93.4%, respectively), and less at pre-dusk (67, 52.6, 56.9, 40.1, 40.7, 33.2, and 90%, respectively), also all these sites were validated by qRT-PCR.

Conclusion

The differential editing of chloroplast ndhB gene across light periods may be led to a somehow relations between the RNA editing and control of photosynthesis.

Systematic analysis reveals cis and trans determinants affecting C-to-U RNA editing in Arabidopsis thaliana

Background

C-to-U RNA editing is prevalent in the mitochondrial and chloroplast genes in plants. The biological functions of a fraction of C-to-U editing sites are continuously discovered by case studies. However, at genome-wide level, the cis and trans determinants affecting the occurrence or editing levels of these C-to-U events are relatively less studied. What is known is that the PPR (pentatricopeptide repeat) proteins are the main trans-regulatory elements responsible for the C-to-U conversion, but other determinants especially the cis-regulatory elements remain largely uninvestigated.

Results

By analyzing the transcriptome and translatome data in Arabidopsis thaliana roots and shoots, combined with RNA-seq data from hybrids of Arabidopsis thaliana and Arabidopsis lyrata, we perform genome-wide investigation on the cis elements and trans-regulatory elements that potentially affect C-to-U editing events. An upstream guanosine or double-stranded RNA (dsRNA) regions are unfavorable for editing events. Meanwhile, many genes including the transcription factors may indirectly play regulatory roles in trans.

Conclusions

The 5-prime thymidine facilitates editing and dsRNA structures prevent editing in cis. Many transcription factors affect editing in trans. Although the detailed molecular mechanisms underlying the cis and trans regulation remain to be experimentally verified, our findings provide novel aspects in studying the botanical C-to-U RNA editing events.

Topologies of N6 -adenosine methylation (m6 A) in land plant mitochondria and their putative effects on organellar gene expression

Mitochondria serve as major sites of ATP production and play key roles in many other metabolic processes that are critical to the cell. As relicts of an ancient bacterial endosymbiont, mitochondria contain their own hereditary material (i.e. mtDNA, or mitogenome) and a machinery for protein biosynthesis. The expression of the mtDNA in plants is complex, particularly at the post-transcriptional level. Following transcription, the polycistronic pre-RNAs undergo extensive modifications, including trimming, splicing and editing, before being translated by organellar ribosomes. Our study focuses on N⁶ -methylation of adenosine ribonucleotides (m⁶ A-RNA) in plant mitochondria. m⁶ A is a prevalent modification in nuclear-encoded mRNAs. The biological significance of this dynamic modification is under investigation, but it is widely accepted that m⁶ A mediates structural switches that affect RNA stability and/or activity. Using m⁶ A-pulldown/RNA-seq (m⁶ A-RIP-seq) assays of Arabidopsis and cauliflower mitochondria, we provide information on the m⁶ A-RNA landscapes in Arabidopsis thaliana and Brassica oleracea mitochondria. The results show that m⁶ A targets different types of mitochondrial transcripts, including known genes, mtORFs, as well as non-coding (transcribed intergenic) RNA species. While ncRNAs undergo multiple m⁶ A modifications, N⁶ -methylation of adenosine residues with mRNAs seem preferably positioned near start codons and may modulate their translatability.

Reduced C-to-U RNA editing rates might play a regulatory role in stress response of Arabidopsis

C-to-U RNA editing is prevalent in the mitochondrial and chloroplast genes in plants. The C-to-U editing rates are constantly very high. During genome evolution, those edited cytidines are likely to be replaced with thymidines at the DNA level. C-to-U editing events are suggested to be designed for reversing the unfavorable T-to-C DNA mutations. Despite the existing theory showing the importance of editing mechanisms, few studies have investigated the genome-wide adaptive signals of the C-to-U editome or the potential function of C-to-U editing events in the stress response. By analyzing the transcriptome and translatome data of normal and heat-shocked Arabidopsis thaliana and the RNA-seq from cold-stressed plants, combined with genome-wide comparison of mitochondrial/chloroplast genes and nuclear genes from multiple aspects, we present the conservational and translational features of each gene and depict the dynamic mitochondrial/chloroplast C-to-U RNA editome. We found that the tAI (tRNA adaptation index) and basic translation levels are lower for mitochondrial/chloroplast genes than for nuclear genes. Interestingly, although we found adaptive signals for the global C-to-U RNA editome in mitochondrial/chloroplast genes, the C-to-U (T) alteration would usually cause a reduction in the codon tAI value. Moreover, the C-to-U editing rates are significantly reduced under heat or cold stress when compared to the normal condition. This reduction is irrelevant to the temperature-sensitive RNA structures. Several cases have illustrated that under heat stress, the reduced C-to-U editing rates alleviate ribosome stalling and consequently facilitate the local translation. Our study reveals that in Arabidopsis thaliana the mitochondrial/chloroplast C-to-U RNA editing rates are reduced under heat or cold stress. This reduction is associated with the alleviation of decreased tAI/translation rate of edited codons. The regulation of C-to-U editing rates could be the tradeoff between quantity and quality. We profile the dynamic change of C-to-U RNA editome under heat stress and propose a potential role of editing sites in the heat response. Our work should be appealing to the plant physiologists as well as the RNA editing community.

Evolutionary Model of Plastidial RNA Editing in Angiosperms Presumed from Genome-Wide Analysis of Amborella trichopoda

Amborella trichopoda is placed close to the base of the angiosperm lineage (basal angiosperm). By genome-wide RNA sequencing, we identified 184C-to-U RNA editing sites in the plastid genome of Amborella. This number is much higher than that observed in other angiosperms including maize (44 sites), rice (39 sites) and grape (115 sites). Despite the high frequency of RNA editing, the biased distribution of RNA editing sites in the genome, target codon preference and nucleotide preference adjacent to the edited cytidine are similar to that in other angiosperms, suggesting a common editing machinery. Consistent with this idea, the Amborella nuclear genome encodes 2-3 times more of the E- and DYW-subclass members of pentatricopeptide repeat proteins responsible for RNA editing site recognition in plant organelles. Among 165 editing sites in plastid protein coding sequences in Amborella, 100 sites were conserved at least in one out of 38 species selected to represent key branching points of the angiosperm phylogenetic tree. We assume these 100 sites represent at least a subset of the sites in the plastid editotype of ancestral angiosperms. We then mapped the loss and gain of editing sites on the phylogenetic tree of angiosperms. Our results support the idea that the evolution of angiosperms has led to the loss of RNA editing sites in plastids.

Complete Sequence, Multichromosomal Architecture and Transcriptome Analysis of the Solanum tuberosum Mitochondrial Genome

Mitochondrial genomes (mitogenomes) in higher plants can induce cytoplasmic male sterility and be somehow involved in nuclear-cytoplasmic interactions affecting plant growth and agronomic performance. They are larger and more complex than in other eukaryotes, due to their recombinogenic nature. For most plants, the mitochondrial DNA (mtDNA) can be represented as a single circular chromosome, the so-called master molecule, which includes repeated sequences that recombine frequently, generating sub-genomic molecules in various proportions. Based on the relevance of the potato crop worldwide, herewith we report the complete mtDNA sequence of two S. tuberosum cultivars, namely Cicero and Désirée, and a comprehensive study of its expression, based on high-coverage RNA sequencing data. We found that the potato mitogenome has a multi-partite architecture, divided in at least three independent molecules that according to our data should behave as autonomous chromosomes. Inter-cultivar variability was null, while comparative analyses with other species of the Solanaceae family allowed the investigation of the evolutionary history of their mitogenomes. The RNA-seq data revealed peculiarities in transcriptional and post-transcriptional processing of mRNAs. These included co-transcription of genes with open reading frames that are probably expressed, methylation of an rRNA at a position that should impact translation efficiency and extensive RNA editing, with a high proportion of partial editing implying frequent mis-targeting by the editing machinery.

The chloroplast and mitochondrial C-to-U RNA editing in Arabidopsis thaliana shows signals of adaptation

C-to-U RNA editing is the conversion from cytidine to uridine at RNA level. In plants, the genes undergo C-to-U RNA modification are mainly chloroplast and mitochondrial genes. Case studies have identified the roles of C-to-U editing in various biological processes, but the functional consequence of the majority of C-to-U editing events is still undiscovered. We retrieved the deep sequenced transcriptome data in roots and shoots of Arabidopsis thaliana and profiled their C-to-U RNA editomes and gene expression patterns. We investigated the editing level and conservation pattern of these C-to-U editing sites. The levels of nonsynonymous C-to-U editing events are higher than levels of synonymous events. The fraction of nonsynonymous editing sites is higher than neutral expectation. Highly edited cytidines are more conserved at DNA level, and the gene expression levels are correlated with C-to-U editing levels. Our results demonstrate that the global C-to-U editome is shaped by natural selection and that many nonsynonymous C-to-U editing events are adaptive. The editing mechanism might be positively selected and maintained and could have profound effects on the modified RNAs.

Potential of Transcript Editing Across Mitogenomes of Early Land Plants Shows Novel and Familiar Trends

RNA editing alters the identity of nucleotides in an RNA sequence so that the mature transcript differs from the template defined in the genome. This process has been observed in chloroplasts and mitochondria of both seed and early land plants. However, the frequency of RNA editing in plant mitochondria ranges from zero to thousands of editing sites. To date, analyses of RNA editing in mitochondria of early land plants have been conducted on a small number of genes or mitochondrial genomes of a single species. This study provides an overview of the mitogenomic RNA editing potential of the main lineages of these two groups of early land plants by predicting the RNA editing sites of 33 mitochondrial genes of 37 species of liverworts and mosses. For the purpose of the research, we newly assembled seven mitochondrial genomes of liverworts. The total number of liverwort genera with known complete mitogenome sequences has doubled and, as a result, the available complete mitogenome sequences now span almost all orders of liverworts. The RNA editing site predictions revealed that C-to-U RNA editing in liverworts and mosses is group-specific. This is especially evident in the case of liverwort lineages. The average level of C-to-U RNA editing appears to be over three times higher in liverworts than in mosses, while the C-to-U editing frequency of the majority of genes seems to be consistent for each gene across bryophytes.

FEDRO: a software tool for the automatic discovery of candidate ORFs in plants with c →u RNA editing

BACKGROUND

RNA editing is an important mechanism for gene expression in plants organelles. It alters the direct transfer of genetic information from DNA to proteins, due to the introduction of differences between RNAs and the corresponding coding DNA sequences. Software tools successful for the search of genes in other organisms not always are able to correctly perform this task in plants organellar genomes. Moreover, the available software tools predicting RNA editing events utilise algorithms that do not account for events which may generate a novel start codon.

RESULTS

We present FEDRO, a Java software tool implementing a novel strategy to generate candidate Open Reading Frames (ORFs) resulting from Cytidine to Uridine (c→u) editing substitutions which occur in the mitochondrial genome (mtDNA) of a given input plant. The goal is to predict putative proteins of plants mitochondria that have not been yet annotated. In order to validate the generated ORFs, a screening is performed by checking for sequence similarity or presence in active transcripts of the same or similar organisms. We illustrate the functionalities of our framework on a model organism.

CONCLUSIONS

The proposed tool may be used also on other organisms and genomes. FEDRO is publicly available at http://math.unipa.it/rombo/FEDRO .

EMP18 functions in mitochondrial atp6 and cox2 transcript editing and is essential to seed development in maize

RNA editing plays an important role in organellar gene expression in plants, and pentatricopeptide repeat (PPR) proteins are involved in this function. Because of its large family size, many PPR proteins are not known for their function and roles in plant growth and development. Through genetic and molecular analyses of the empty pericarp18 (emp18) mutant in maize (Zea mays), we cloned the Emp18 gene, revealed its molecular function, and defined its role in the mitochondrial complex assembly and seed development. Emp18 encodes a mitochondrial-localized DYW-PPR protein. Null mutation of Emp18 arrests embryo and endosperm development at an early stage in maize, resulting in embryo lethality. Mutants are deficient in the cytidine (C)-to-uridine (U) editing at atp6-635 and cox2-449, which converts a Leu to Pro in ATP6 and a Met to Thr in Cox2. The atp6 gene encodes the subunit a of F₁ F_o -ATPase. The Leu to Pro alteration disrupts an α-helix of subunit a, resulting in a dramatic reduction in assembly and activity of F₁ F_o -ATPase holoenzyme and an accumulation of free F₁ -subcomplex. These results demonstrate that EMP18 functions in the C-to-U editing of atp6 and cox2, and is essential to mitochondrial biogenesis and seed development in maize.

Molecular and Functional Diversity of RNA Editing in Plant Mitochondria

RNA editing is a fundamental biochemical process relating to the modification of nucleotides in messenger RNAs of functional genes in cells. RNA editing leads to re-establishment of conserved amino acid residues for functional proteins in nuclei, chloroplasts, and mitochondria. Identification of RNA editing factors that contributes to target site recognition increases our understanding of RNA editing mechanisms. Significant progress has been made in recent years in RNA editing studies for both animal and plant cells. RNA editing in nuclei and mitochondria of animal cells and in chloroplast of plant cells has been extensively documented and reviewed. RNA editing has been also extensively documented on plant mitochondria. However, functional diversity of RNA editing factors in plant mitochondria is not overviewed. Here, we review the biological significance of RNA editing, recent progress on the molecular mechanisms of RNA editing process, and function diversity of editing factors in plant mitochondrial research. We will focus on: (1) pentatricopeptide repeat proteins in Arabidopsis and in crop plants; (2) the progress of RNA editing process in plant mitochondria; (3) RNA editing-related RNA splicing; (4) RNA editing associated flower development; (5) RNA editing modulated male sterile; (6) RNA editing-regulated cell signaling; and (7) RNA editing involving abiotic stress. Advances described in this review will be valuable in expanding our understanding in RNA editing. The diverse functions of RNA editing in plant mitochondria will shed light on the investigation of molecular mechanisms that underlies plant development and abiotic stress tolerance.

The Pentatricopeptide Repeat Protein MEF31 is Required for Editing at Site 581 of the Mitochondrial tatC Transcript and Indirectly Influences Editing at Site 586 of the Same Transcript

Pentatricopeptide repeat (PPR) proteins constitute the largest family of proteins in angiosperms, and most members are predicted to play roles in the maturation of organellar RNAs. Here we describe the novel mitochondrial editing factor 31 (MEF31), an E-PPR protein involved in editing at two close sites in the same transcript encoding subunit C of the twin-arginine translocation (tat) pathway. MEF31 is essential for editing at site tatC-581 and application of the recently proposed amino acid code for RNA recognition by PPR proteins supports the view that MEF31 directly targets this site by recognizing its cis sequence. In contrast, editing at site tatC-586 five nucleotides downstream is only partially affected in plants lacking MEF31, being restored to wild-type levels in complemented plants. Application of the amino acid code and analysis of individual RNA molecules for editing at sites 581 and 586 suggest that MEF31 does not directly target site tatC-586, and only indirectly influences editing at this site. It is likely that editing at site tatC-581 improves recognition of the site tatC-586 cis sequence by a second unknown PPR protein.

Rapid evolutionary divergence of diploid and allotetraploid Gossypium mitochondrial genomes

Background

Cotton (Gossypium spp.) is commonly grouped into eight diploid genomic groups and an allotetraploid genomic group, AD. The mitochondrial genomes supply new information to understand both the evolution process and the mechanism of cytoplasmic male sterility. Based on previously released mitochondrial genomes of G. hirsutum (AD₁), G. barbadense (AD₂), G. raimondii (D₅) and G. arboreum (A₂), together with data of six other mitochondrial genomes, to elucidate the evolution and diversity of mitochondrial genomes within Gossypium.

Results

Six Gossypium mitochondrial genomes, including three diploid species from D and three allotetraploid species from AD genome groups (G. thurberi D₁, G. davidsonii D_3-d and G. trilobum D₈; G. tomentosum AD₃, G. mustelinum AD₄ and G. darwinii AD₅), were assembled as the single circular molecules of lengths about 644 kb in diploid species and 677 kb in allotetraploid species, respectively. The genomic structures of mitochondrial in D group species were identical but differed from the mitogenome of G. arboreum (A₂), as well as from the mitogenomes of five species of the AD group. There mainly existed four or six large repeats in the mitogenomes of the A + AD or D group species, respectively. These variations in repeat sequences caused the major inversions and translocations within the mitochondrial genome. The mitochondrial genome complexity in Gossypium presented eight unique segments in D group species, three specific fragments in A + AD group species and a large segment (more than 11 kb) in diploid species. These insertions or deletions were most probably generated from crossovers between repetitive or homologous regions. Unlike the highly variable genome structure, evolutionary distance of mitochondrial genes was 1/6th the frequency of that in chloroplast genes of Gossypium. RNA editing events were conserved in cotton mitochondrial genes. We confirmed two near full length of the integration of the mitochondrial genome into chromosome 1 of G. raimondii and chromosome A03 of G. hirsutum, respectively, with insertion time less than 1.03 MYA.

Conclusion

Ten Gossypium mitochondrial sequences highlight the insights to the evolution of cotton mitogenomes.

Electronic supplementary material

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

Transcriptomic analysis of salt stress responsive genes in Rhazya stricta

Rhazya stricta is an evergreen shrub that is widely distributed across Western and South Asia, and like many other members of the Apocynaceae produces monoterpene indole alkaloids that have anti-cancer properties. This species is adapted to very harsh desert conditions making it an excellent system for studying tolerance to high temperatures and salinity. RNA-Seq analysis was performed on R. stricta exposed to severe salt stress (500 mM NaCl) across four time intervals (0, 2, 12 and 24 h) to examine mechanisms of salt tolerance. A large number of transcripts including genes encoding tetrapyrroles and pentatricopeptide repeat (PPR) proteins were regulated only after 12 h of stress of seedlings grown in controlled greenhouse conditions. Mechanisms of salt tolerance in R. stricta may involve the upregulation of genes encoding chaperone protein Dnaj6, UDP-glucosyl transferase 85a2, protein transparent testa 12 and respiratory burst oxidase homolog protein b. Many of the highly-expressed genes act on protecting protein folding during salt stress and the production of flavonoids, key secondary metabolites in stress tolerance. Other regulated genes encode enzymes in the porphyrin and chlorophyll metabolic pathway with important roles during plant growth, photosynthesis, hormone signaling and abiotic responses. Heme biosynthesis in R. stricta leaves might add to the level of salt stress tolerance by maintaining appropriate levels of photosynthesis and normal plant growth as well as by the participation in reactive oxygen species (ROS) production under stress. We speculate that the high expression levels of PPR genes may be dependent on expression levels of their targeted editing genes. Although the results of PPR gene family indicated regulation of a large number of transcripts under salt stress, PPR actions were independent of the salt stress because their RNA editing patterns were unchanged.

Salt stress affects mRNA editing in soybean chloroplasts

Soybean, a crop known by its economic and nutritional importance, has been the subject of several studies that assess the impact and the effective plant responses to abiotic stresses. Salt stress is one of the main environmental stresses and negatively impacts crop growth and yield. In this work, the RNA editing process in the chloroplast of soybean plants was evaluated in response to a salt stress. Bioinformatics approach using sRNA and mRNA libraries were employed to detect specific sites showing differences in editing efficiency. RT-qPCR was used to measure editing efficiency at selected sites. We observed that transcripts of NDHA, NDHB, RPS14 and RPS16 genes presented differences in coverage and editing rates between control and salt-treated libraries. RT-qPCR assays demonstrated an increase in editing efficiency of selected genes. The salt stress enhanced the RNA editing process in transcripts, indicating responses to components of the electron transfer chain, photosystem and translation complexes. These increases can be a response to keep the homeostasis of chloroplast protein functions in response to salt stress.

Identification of Symmetrical RNA Editing Events in the Mitochondria of Salvia miltiorrhiza by Strand-specific RNA Sequencing

Salvia miltiorrhiza is one of the most widely-used medicinal plants. Here, we systematically analyzed the RNA editing events in its mitochondria. We developed a pipeline using REDItools to predict RNA editing events from stand-specific RNA-Seq data. The predictions were validated using reverse transcription, RT-PCR amplification and Sanger sequencing experiments. Putative sequences motifs were characterized. Comparative analyses were carried out between S. miltiorrhiza, Arabidopsis thaliana and Oryza sativa. We discovered 1123 editing sites, including 225 “C to U” sites in the protein-coding regions. Fourteen of sixteen (87.5%) sites were validated. Three putative DNA motifs were identified around the predicted sites. The nucleotides on both strands at 115 of the 225 sites had undergone RNA editing, which we called symmetrical RNA editing (SRE). Four of six these SRE sites (66.7%) were experimentally confirmed. Re-examination of strand-specific RNA-Seq data from A. thaliana and O. sativa identified 327 and 369 SRE sites respectively. 78, 20 and 13 SRE sites were found to be conserved among A. thaliana, O. sativa and S. miltiorrhiza respectively. This study provides a comprehensive picture of RNA editing events in the mitochondrial genome of S. miltiorrhiza. We identified SREs for the first time, which may represent a universal phenomenon.

Synonymous codon usage bias in plant mitochondrial genes is associated with intron number and mirrors species evolution

Synonymous codon usage bias (SCUB) is a common event that a non-uniform usage of codons often occurs in nearly all organisms. We previously found that SCUB is correlated with both intron number and exon position in the plant nuclear genome but not in the plastid genome; SCUB in both nuclear and plastid genome can mirror the evolutionary specialization. However, how about the rules in the mitochondrial genome has not been addressed. Here, we present an analysis of SCUB in the mitochondrial genome, based on 24 plant species ranging from algae to land plants. The frequencies of NNA and NNT (A- and T-ending codons) are higher than those of NNG and NNC, with the strongest preference in bryophytes and the weakest in land plants, suggesting an association between SCUB and plant evolution. The preference for NNA and NNT is more evident in genes harboring a greater number of introns in land plants, but the bias to NNA and NNT exhibits even among exons. The pattern of SCUB in the mitochondrial genome differs in some respects to that present in both the nuclear and plastid genomes.

RNA sociology: group behavioral motifs of RNA consortia

RNA sociology investigates the behavioral motifs of RNA consortia from the social science perspective. Besides the self-folding of RNAs into single stem loop structures, group building of such stem loops results in a variety of essential agents that are highly active in regulatory processes in cellular and non-cellular life. RNA stem loop self-folding and group building do not depend solely on sequence syntax; more important are their contextual (functional) needs. Also, evolutionary processes seem to occur through RNA stem loop consortia that may act as a complement. This means the whole entity functions only if all participating parts are coordinated, although the complementary building parts originally evolved for different functions. If complementary groups, such as rRNAs and tRNAs, are placed together in selective pressure contexts, new evolutionary features may emerge. Evolution initiated by competent agents in natural genome editing clearly contrasts with statistical error replication narratives.

The pentatricopeptide repeat protein MEF26 participates in RNA editing in mitochondrial cox3 and nad4 transcripts

In angiosperms most members of the large nuclear-encoded family of pentatricopeptide repeat (PPR) proteins are predicted to play relevant roles in the maturation of organellar RNAs. Here we report the novel Mitochondrial Editing Factor 26, a DYW-PPR protein involved in RNA editing at two sites. While at one site, cox3-311, editing is abolished in the absence of MEF26, the other site, nad4-166, is still partially edited. These sites share similar cis-elements and application of the recently proposed amino acid code for RNA recognition by PPR proteins ranks them at first and second positions of the most probable targets.

RNA editing in plants and its evolution

RNA editing alters the identity of nucleotides in RNA molecules such that the information for a protein in the mRNA differs from the prediction of the genomic DNA. In chloroplasts and mitochondria of flowering plants, RNA editing changes C nucleotides to U nucleotides; in ferns and mosses, it also changes U to C. The approximately 500 editing sites in mitochondria and 40 editing sites in plastids of flowering plants are individually addressed by specific proteins, genes for which are amplified in plant species with organellar RNA editing. These proteins contain repeat elements that bind to cognate RNA sequence motifs just 5' to the edited nucleotide. In flowering plants, the site-specific proteins interact selectively with individual members of a different, smaller family of proteins. These latter proteins may be connectors between the site-specific proteins and the as yet unknown deaminating enzymatic activity.

RNA editing in plant mitochondria—connecting RNA target sequences and acting proteins

RNA editing changes several hundred cytidines to uridines in the mRNAs of mitochondria in flowering plants. The target cytidines are identified by a subtype of PPR proteins characterized by tandem modules which each binds with a specific upstream nucleotide. Recent progress in correlating repeat structures with nucleotide identities allows to predict and identify target sites in mitochondrial RNAs. Additional proteins have been found to play a role in RNA editing; their precise function still needs to be elucidated. The enzymatic activity performing the C to U reaction may reside in the C-terminal DYW extensions of the PPR proteins; however, this still needs to be proven. Here we update recent progress in understanding RNA editing in flowering plant mitochondria.

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.

The life of plant mitochondrial complex I

The mitochondrial NADH dehydrogenase complex (complex I) of the respiratory chain has several remarkable features in plants: (i) particularly many of its subunits are encoded by the mitochondrial genome, (ii) its mitochondrial transcripts undergo extensive maturation processes (e.g. RNA editing, trans-splicing), (iii) its assembly follows unique routes, (iv) it includes an additional functional domain which contains carbonic anhydrases and (v) it is, indirectly, involved in photosynthesis. Comprising about 50 distinct protein subunits, complex I of plants is very large. However, an even larger number of proteins are required to synthesize these subunits and assemble the enzyme complex. This review aims to follow the complete "life cycle" of plant complex I from various molecular perspectives. We provide arguments that complex I represents an ideal model system for studying the interplay of respiration and photosynthesis, the cooperation of mitochondria and the nucleus during organelle biogenesis and the evolution of the mitochondrial oxidative phosphorylation system.

Group II intron splicing factors in plant mitochondria

Group II introns are large catalytic RNAs (ribozymes) which are found in bacteria and organellar genomes of several lower eukaryotes, but are particularly prevalent within the mitochondrial genomes (mtDNA) in plants, where they reside in numerous critical genes. Their excision is therefore essential for mitochondria biogenesis and 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 intron-encoded maturase protein. A hallmark of maturases is that they are intron specific, acting as cofactors which bind their own cognate containing pre-mRNAs to facilitate splicing. However, the plant organellar introns have diverged considerably from their bacterial ancestors, such as they lack many regions which are necessary for splicing and also lost their evolutionary related maturase ORFs. In fact, only a single maturase has been retained in the mtDNA of various angiosperms: the matR gene encoded in the fourth intron of the NADH-dehydrogenase subunit 1 (nad1 intron 4). Their degeneracy and the absence of cognate ORFs suggest that the splicing of plant mitochondria introns is assisted by trans-acting cofactors. Interestingly, in addition to MatR, the nuclear genomes of angiosperms also harbor four genes (nMat 1-4), which are closely related to maturases and 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. In addition to the nMATs, genetic screens led to the identification of other genes encoding various factors, which are required for the splicing and processing of mitochondrial introns in plants. In this review we will summarize recent data on the splicing and processing of mitochondrial introns and their implication in plant development and physiology, with a focus on maturases and their accessory splicing cofactors.

Deletions in cox2 mRNA result in loss of splicing and RNA editing and gain of novel RNA editing sites

As previously demonstrated, the maize cox2 RNA is fully edited in cauliflower mitochondria. Use of constructs with a deleted cox2 intron, however, led to a loss of RNA editing at almost all editing sites, with only a few sites still partially edited. Likewise, one deletion in exon 1 and three in exon 2 abolish RNA editing at all cox2 sites analyzed. Furthermore, intron splicing is abolished using these deletions. Mutation of a cytosine residue, which is normally edited and localized directly adjacent to the intron, to thymidine did not result in restoration of splicing, indicating that the loss of splicing was not due to loss of RNA editing. One deletion in exon 2 did not lead to loss of splicing. Instead, most editing sites were found to be edited, only three were not edited. Unexpectedly, we observed additional RNA editing events at new sites. Thus it appears that deletions in the cox2 RNA sequence can have a strong effect on RNA processing, leading to loss of splicing, loss of editing at all sites, or even to a gain of new editing sites. As these effects are not limited to the vicinity of the respective deletions, but appear to be widespread or even affect all editing sites, they may not be explained by the loss of PPR binding sites. Instead, it appears that several parts of the cox2 transcript are required for proper RNA processing. This indicates the roles of the RNA sequence and structural elements in the recognition of the editing sites.

MASCP gator: an overview of the Arabidopsis proteomic aggregation portal

A key challenge in the area of bioinformatics in the coming decades is the ability to manage the wealth of information that is being generated from the variety of high throughput methodologies currently being undertaken in laboratories across the world. While these approaches have made available large volumes of data to the research community, less attention has been given to the problem of how to intuitively present the data to enable greater biological insights. Recently, an attempt was made to tackle this problem in the area of Arabidopsis proteomics. The model plant has been the target of countless proteomics surveys producing an exhaustive array of data and online repositories. The MASCP Gator is an aggregation portal for proteomic data currently being produced by the community and unites a large collection of specialized resources to a single portal (http://gator.masc-proteomics.org/). Here we describe the latest additions, upgrades and features to this resource further expanding its role into protein modifications and genome sequence variations.

The complete sequence of the mitochondrial genome of Butomus umbellatus--a member of an early branching lineage of monocotyledons

In order to study the evolution of mitochondrial genomes in the early branching lineages of the monocotyledons, i.e., the Acorales and Alismatales, we are sequencing complete genomes from a suite of key taxa. As a starting point the present paper describes the mitochondrial genome of Butomus umbellatus (Butomaceae) based on next-generation sequencing data. The genome was assembled into a circular molecule, 450,826 bp in length. Coding sequences cover only 8.2% of the genome and include 28 protein coding genes, four rRNA genes, and 12 tRNA genes. Some of the tRNA genes and a 16S rRNA gene are transferred from the plastid genome. However, the total amount of recognized plastid sequences in the mitochondrial genome is only 1.5% and the amount of DNA transferred from the nucleus is also low. RNA editing is abundant and a total of 557 edited sites are predicted in the protein coding genes. Compared to the 40 angiosperm mitochondrial genomes sequenced to date, the GC content of the Butomus genome is uniquely high (49.1%). The overall similarity between the mitochondrial genomes of Butomus and Spirodela (Araceae), the closest relative yet sequenced, is low (less than 20%), and the two genomes differ in size by a factor 2. Gene order is also largely unconserved. However, based on its phylogenetic position within the core alismatids Butomus will serve as a good reference point for subsequent studies in the early branching lineages of the monocotyledons.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

The genome and transcriptome of perennial ryegrass mitochondria

Background

Perennial ryegrass (Lolium perenne L.) is one of the most important forage and turf grass species of temperate regions worldwide. Its mitochondrial genome is inherited maternally and contains genes that can influence traits of agricultural importance. Moreover, the DNA sequence of mitochondrial genomes has been established and compared for a large number of species in order to characterize evolutionary relationships. Therefore, it is crucial to understand the organization of the mitochondrial genome and how it varies between and within species. Here, we report the first de novo assembly and annotation of the complete mitochondrial genome from perennial ryegrass.

Results

Intact mitochondria from perennial ryegrass leaves were isolated and used for mtDNA extraction. The mitochondrial genome was sequenced to a 167-fold coverage using the Roche 454 GS-FLX Titanium platform, and assembled into a circular master molecule of 678,580 bp. A total of 34 proteins, 14 tRNAs and 3 rRNAs are encoded by the mitochondrial genome, giving a total gene space of 48,723 bp (7.2%). Moreover, we identified 149 open reading frames larger than 300 bp and covering 67,410 bp (9.93%), 250 SSRs, 29 tandem repeats, 5 pairs of large repeats, and 96 pairs of short inverted repeats. The genes encoding subunits of the respiratory complexes – nad1 to nad9, cob, cox1 to cox3 and atp1 to atp9 – all showed high expression levels both in absolute numbers and after normalization.

Conclusions

The circular master molecule of the mitochondrial genome from perennial ryegrass presented here constitutes an important tool for future attempts to compare mitochondrial genomes within and between grass species. Our results also demonstrate that mitochondria of perennial ryegrass contain genes crucial for energy production that are well conserved in the mitochondrial genome of monocotyledonous species. The expression analysis gave us first insights into the transcriptome of these mitochondrial genes in perennial ryegrass.

Evolutionary aspects of a unique internal mitochondrial targeting signal in nuclear-migrated rps19 of sugar beet (Beta vulgaris L.)

The endosymbiotic theory postulates that many genes migrated from endosymbionts to the nuclear genomes of their hosts. Some migrated genes lack presequences directing proteins to mitochondria, and their mitochondrial targeting signals appear to be inscribed in the core coding regions as internal targeting signals (ITSs). ITSs may have evolved after sequence transfer to nuclei or ITSs may have pre-existed before sequence transfer. Here, we report the molecular cloning of a sugar beet gene for ribosomal protein S19 (Rps19; the first letter is capitalized when the gene is a nuclear gene). We show that sugar beet Rps19 (BvRps19) is an ITS-type gene. Based on amino-acid sequence comparison, dicotyledonous rps19s (the first letter is lower-cased when the gene is a mitochondrial gene), such as tobacco rps19 (Ntrps19), resemble an ancestral form of BvRps19. We investigated whether differences in amino-acid sequences between BvRps19 and Ntrps19 were involved in ITS evolution. Analyses of the intracellular localization of chimaeric GFP-fusion proteins that were transiently expressed in Welsh onion cells showed that Ntrps19-gfp was not localized in mitochondria. When several BvRps19-type amino acid substitutions, none of which was seen in any other angiosperm rps19, were introduced into Ntrps19-gfp, the modified Ntrps19-gfp became localized in mitochondria, supporting the notion that an ITS in BvRps19 evolved following sequence transfer to nuclei. Not all of these substitutions were seen in other ITS-type Rps19s, suggesting that the ITSs of Rps19 are diverse.

Nucleo-cytoplasmic interactions affect RNA editing of cox2, atp6 and atp9 in alloplasmic male-sterile rice (Oryza sativa L.) lines

RNA editing plays an important role in the regulation of mitochondrial gene expression in flowering plants. In this study, we examined RNA editing of the mitochondrial genes cox2, atp6 and atp9 in five isonuclear alloplasmic male-sterile lines (IAMSLs) of rice to investigate whether different cytoplasmic types affect RNA editing. Although many editing sites were conserved among the three genes, we found that the editing efficiency of certain sites was significantly different between different IAMSLs or between IAMSLs and their corresponding cytoplasmic donor CMS lines. Furthermore, several editing sites were found to be either present or absent in certain IAMSLs and their corresponding CMS lines. These results indicate that nuclear loci, as well as unknown editing factors within the mitochondria of different cytoplasmic types, may be involved in RNA editing, and they suggest that RNA editing in plant mitochondria is affected by nucleo-cytoplasmic interactions.

Darned in 2013: inclusion of model organisms and linking with Wikipedia

DARNED (DAtabase of RNa EDiting, available at http://darned.ucc.ie) is a centralized repository of reference genome coordinates corresponding to RNA nucleotides having altered templated identities in the process of RNA editing. The data in DARNED are derived from published datasets of RNA editing events. RNA editing instances have been identified with various methods, such as bioinformatics screenings, deep sequencing and/or biochemical techniques. Here we report our current progress in the development and expansion of the DARNED. In addition to novel database features the DARNED update describes inclusion of Drosophila melanogaster and Mus musculus RNA editing events and the launch of a community-based annotation in the RNA WikiProject.

The mitochondrial genome of an aquatic plant, Spirodela polyrhiza

BACKGROUND

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

METHODS

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

CONCLUSIONS

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.

Two related RNA-editing proteins target the same sites in mitochondria of Arabidopsis thaliana

The facilitators for specific cytosine-to-uridine RNA-editing events in plant mitochondria and plastids are pentatricopeptide repeat (PPR)-containing proteins with specific additional C-terminal domains. Here we report the related PPR proteins mitochondrial editing factor 8 (MEF8) and MEF8S with only five such repeats each to be both involved in RNA editing at the same two sites in mitochondria of Arabidopsis thaliana. Mutants of MEF8 show diminished editing in leaves but not in pollen, whereas mutants of the related protein MEF8S show reduced RNA editing in pollen but not in leaves. Overexpressed MEF8 or MEF8S both increase editing at the two target sites in a mef8 mutant. Double mutants of MEF8 and MEF8S are not viable although both identified target sites are in mRNAs for nonessential proteins. This suggests that MEF8 and MEF8S may have other essential functions beyond these two editing sites in complex I mRNAs.