RNA editing in plant mitochondria: 20 years later

The assembly and comparative analysis of the first complete mitogenome of Lindera aggregata

Lindera aggregata, a member belongs to the genus Lindera of Lauraceae family. Its roots and leaves have been used as traditional Chinese medicine or functional food for thousands of years. However, its mitochondrial genome has not been explored. Our aim is to sequence and assemble the mitogenome of L. aggregata to elucidate the genetic mechanism and evolutionary pathway. The results had shown that the mitogenome was extremely complex and had a unique multi-branched conformation with total size of 912,473 bp. Comprehensive analysis of protein coding genes of 7 related species showed that there were 40 common genes in their mitogenome. Interestingly, positive selection had become an important factor in the evolution of ccmB, ccmFC, rps10, rps11 and rps7 genes. Furthermore, our data highlighted the repeated trend of homologous fragment migrations between chloroplast and mitochondrial organelles, and 38 homologous fragments were identified. Phylogenetic analysis identified a tree that was basically consistent with the phylogeny of Laurales species described in the APG IV system. To sum up, this study will be helpful to the study of population genetics and evolution of Lindera species.

Comprehensive Analysis of the Complete Mitochondrial Genome of Rehmannia chingii: An Autotrophic Species in the Orobanchaceae Family

Rehmannia chingii is an important medicinal plant with immense value in scientific research. However, its mitochondrial genome (mitogenome) has not yet been characterized. Herein, based on whole-genome Illumina short reads and PacBio HiFi reads, we obtained the complete mitogenome of R. chingii through a de novo assembly strategy. We carried out comparative genomic analyses and found that, in comparison with the plastid genome (plastome) showing a high degree of structural conservation, the R. chingii mitogenome structure is relatively complex, showing an intricate ring structure with 16 connections, owing to five repetitive sequences. The R. chingii mitogenome was 783,161 bp with a GC content of 44.8% and contained 77 genes, comprising 47 protein-coding genes (CDS), 27 tRNA genes, and 3 rRNA genes. We counted 579 RNA editing events in 47 CDS and 12,828 codons in all CDSs of the R. chingii mitogenome. Furthermore, 24 unique sequence transfer fragments were found between the mitogenome and plastome, comprising 8 mitogenome CDS genes and 16 plastome CDS genes, corresponding to 2.39% of the R. chingii mitogenome. Mitogenomes had shorter but more collinear regions, evidenced by a comparison of the organelles of non-parasitic R. chingii, hemiparasitic Pedicularis chinensis, and holoparasitic Aeginetia indica in the Orobanchaceae family. Moreover, from non-parasitic to holoparasitic species, the genome size in the mitogenomes of Orobanchaceae species did not decrease gradually. Instead, the smallest mitogenome was found in the hemiparasitic species P. chinensis, with a size of 225,612 bp. The findings fill the gap in the mitogenome research of the medicinal plant R. chingii, promote the progress of the organelle genome research of the Orobanchaceae family, and provide clues for molecular breeding.

Comparative analysis of mitochondrial genomes of Broussonetia spp. (Moraceae) reveals heterogeneity in structure, synteny, intercellular gene transfer, and RNA editing

The genus Broussonetia (Moraceae) is comprised of three non-hybrid recognized species that all produce high quality fiber essential in the development of papermaking and barkcloth-making technology. In addition, these species also have medicinal value in several countries. Despite their important economical, medicinal, and ecological values, the complete mitogenome of Broussonetia has not been reported and investigated, which would greatly facilitate molecular phylogenetics, species identification and understanding evolutionary processes. Here, we assembled the first-reported three complete Broussonetia (B. papyrifera, B. kaempferi, and B. monoica) mitochondrial genomes (mitogenome) based on a hybrid strategy using Illumina and Oxford Nanopore Technology sequencing data, and performed comprehensive comparisons in terms of their structure, gene content, synteny, intercellular gene transfer, phylogeny, and RNA editing. Our results showed their huge heterogeneities among the three species. Interestingly, the mitogenomes of B. monoica and B. papyrifera consisted of a single circular structure, whereas the B. kaempferi mitogenome was unique and consisted of a double circular structure. Gene content was consistent except for a few transfer RNA (tRNA) genes. The Broussonetia spp. mitogenomes had high sequence conservation but B. monoica with B. kaempferi contained more synteny blocks and were more related, a finding that was well-supported in organellar phylogeny. Fragments that had been transferred between mitogenomes were detected at plastome hotspots that had integrated under potential mediation of tRNA genes. In addition, RNA editing sites showed great differences in abundance, type, location and efficiency among species and tissues. The availability of these complete gap-free mitogenomes of Broussonetia spp. will provide a valuable genetic resource for evolutionary research and understanding the communications between the two organelle genomes.

Cloning and Organelle Expression of Bamboo Mitochondrial Complex I Subunits Nad1, Nad2, Nad4, and Nad5 in the Yeast Saccharomyces cerevisiae

Mitochondrial respiratory complex I catalyzes electron transfer from NADH to ubiquinone and pumps protons from the matrix into the intermembrane space. In particular, the complex I subunits Nad1, Nad2, Nad4, and Nad5, which are encoded by the nad1, nad2, nad4, and nad5 genes, reside at the mitochondrial inner membrane and possibly function as proton (H⁺) and ion translocators. To understand the individual functional roles of the Nad1, Nad2, Nad4, and Nad5 subunits in bamboo, each cDNA of these four genes was cloned into the pYES2 vector and expressed in the mitochondria of the yeast Saccharomyces cerevisiae. The mitochondrial targeting peptide mt gene (encoding MT) and the egfp marker gene (encoding enhanced green fluorescent protein, EGFP) were fused at the 5'-terminal and 3'-terminal ends, respectively. The constructed plasmids were then transformed into yeast. RNA transcripts and fusion protein expression were observed in the yeast transformants. Mitochondrial localizations of the MT-Nad1-EGFP, MT-Nad2-EGFP, MT-Nad4-EGFP, and MT-Nad5-EGFP fusion proteins were confirmed by fluorescence microscopy. The ectopically expressed bamboo subunits Nad1, Nad2, Nad4, and Nad5 may function in ion translocation, which was confirmed by growth phenotype assays with the addition of different concentrations of K⁺, Na⁺, or H⁺.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

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.

Cytonuclear Genetic Incompatibilities in Plant Speciation

Due to the endosymbiotic origin of organelles, a pattern of coevolution and coadaptation between organellar and nuclear genomes is required for proper cell function. In this review, we focus on the impact of cytonuclear interaction on the reproductive isolation of plant species. We give examples of cases where species exhibit barriers to reproduction which involve plastid-nuclear or mito-nuclear genetic incompatibilities, and describe the evolutionary processes at play. We also discuss potential mechanisms of hybrid fitness recovery such as paternal leakage. Finally, we point out the possible interplay between plant mating systems and cytonuclear coevolution, and its consequence on plant speciation.

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.

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.

Plant Mitochondrial Inner Membrane Protein Insertion

During the biogenesis of the mitochondrial inner membrane, most nuclear-encoded inner membrane proteins are laterally released into the membrane by the TIM23 and the TIM22 machinery during their import into mitochondria. A subset of nuclear-encoded mitochondrial inner membrane proteins and all the mitochondrial-encoded inner membrane proteins use the Oxa machinery-which is evolutionarily conserved from the endosymbiotic bacterial ancestor of mitochondria-for membrane insertion. Compared to the mitochondria from other eukaryotes, plant mitochondria have several unique features, such as a larger genome and a branched electron transport pathway, and are also involved in additional cellular functions such as photorespiration and stress perception. This review focuses on the unique aspects of plant mitochondrial inner membrane protein insertion machinery, which differs from that in yeast and humans, and includes a case study on the biogenesis of Cox2 in yeast, humans, two plant species, and an algal species to highlight lineage-specific similarities and differences. Interestingly, unlike mitochondria of other eukaryotes but similar to bacteria and chloroplasts, plant mitochondria appear to use the Tat machinery for membrane insertion of the Rieske Fe/S protein.

Evolutionary origin of RNA editing

The term "RNA editing" encompasses a wide variety of mechanistically and phylogenetically unrelated processes that change the nucleotide sequence of an RNA species relative to that of the encoding DNA. Two general classes of editing, substitution and insertion/deletion, have been described, with all major types of cellular RNA (messenger, ribosomal, and transfer) undergoing editing in different organisms. In cases where RNA editing is required for function (e.g., to generate a translatable open reading frame in a mRNA), editing is an obligatory step in the pathway of genetic information expression. How, when, and why individual RNA editing systems originated are intriguing biochemical and evolutionary questions. Here I review briefly what is known about the biochemistry, genetics, and phylogenetics of several very different RNA editing systems, emphasizing what we can deduce about their origin and evolution from the molecular machinery involved. An evolutionary model, centered on the concept of "constructive neutral evolution", is able to account in a general way for the origin of RNA editing systems. The model posits that the biochemical elements of an RNA editing system must be in place before there is an actual need for editing, and that RNA editing systems are inherently mutagenic because they allow potentially deleterious or lethal mutations to persist at the genome level, whereas they would otherwise be purged by purifying selection.

Characterization of the chloroplast genome sequence of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is an economically important crop, which is grown for oil production. To better understand the molecular basis of oil palm chloroplasts, we characterized the complete chloroplast (cp) genome sequence obtained from 454 pyrosequencing. The oil palm cp genome is 156,973 bp in length consisting of a large single-copy region of 85,192 bp flanked on each side by inverted repeats of 27,071 bp with a small single-copy region of 17,639 bp joining the repeats. The genome contains 112 unique genes: 79 protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. By aligning the cp genome sequence with oil palm cDNA sequences, we observed 18 non-silent and 10 silent RNA editing events among 19 cp protein-coding genes. Creation of an initiation codon by RNA editing in rpl2 has been reported in several monocots and was also found in the oil palm cp genome. Fifty common chloroplast protein-coding genes from 33 plant taxa were used to construct ML and MP phylogenetic trees. Their topologies are similar and strongly support for the position of E. guineensis as the sister of closely related species Phoenix dactylifera in Arecaceae (palm families) of monocot subtrees.

How a neutral evolutionary ratchet can build cellular complexity

Complex cellular machines and processes are commonly believed to be products of selection, and it is typically understood to be the job of evolutionary biologists to show how selective advantage can account for each step in their origin and subsequent growth in complexity. Here, we describe how complex machines might instead evolve in the absence of positive selection through a process of "presuppression," first termed constructive neutral evolution (CNE) more than a decade ago. If an autonomously functioning cellular component acquires mutations that make it dependent for function on another, pre-existing component or process, and if there are multiple ways in which such dependence may arise, then dependence inevitably will arise and reversal to independence is unlikely. Thus, CNE is a unidirectional evolutionary ratchet leading to complexity, if complexity is equated with the number of components or steps necessary to carry out a cellular process. CNE can explain "functions" that seem to make little sense in terms of cellular economy, like RNA editing or splicing, but it may also contribute to the complexity of machines with clear benefit to the cell, like the ribosome, and to organismal complexity overall. We suggest that CNE-based evolutionary scenarios are in these and other cases less forced than the selectionist or adaptationist narratives that are generally told.

The evolution of RNA editing and pentatricopeptide repeat genes

The pentatricopeptide repeat (PPR) is a degenerate 35-amino-acid structural motif identified from analysis of the sequenced genome of the model plant Arabidopsis thaliana. From the wealth of sequence information now available from plant genomes, the PPR protein family is now known to be one of the largest families in angiosperm species, as most genomes encode 400-600 members. As the number of PPR genes is generally only c. 10-20 in other eukaryotic organisms, including green algae, the family has obviously greatly expanded during land plant evolution. This provides a rare opportunity to study selection pressures driving a 50-fold expansion of a single gene family. PPR proteins are sequence-specific RNA-binding proteins involved in many aspects of RNA processing in organelles. In this review, we will summarize our current knowledge about the evolution of PPR genes, and will discuss the relevance of the dramatic expansion in the family to the functional diversification of plant organelles, focusing primarily on RNA editing.

Evolution of exon-intron structure and alternative splicing

Despite significant advances in high-throughput DNA sequencing, many important species remain understudied at the genome level. In this study we addressed a question of what can be predicted about the genome-wide characteristics of less studied species, based on the genomic data from completely sequenced species. Using NCBI databases we performed a comparative genome-wide analysis of such characteristics as alternative splicing, number of genes, gene products and exons in 36 completely sequenced model species. We created statistical regression models to fit these data and applied them to loblolly pine (Pinus taeda L.), an example of an important species whose genome has not been completely sequenced yet. Using these models, the genome-wide characteristics, such as total number of genes and exons, can be roughly predicted based on parameters estimated from available limited genomic data, e.g. exon length and exon/gene ratio.

Are substitution rates and RNA editing correlated?

Background

RNA editing is a post-transcriptional process that, in seed plants, involves a cytosine to uracil change in messenger RNA, causing the translated protein to differ from that predicted by the DNA sequence. RNA editing occurs extensively in plant mitochondria, but large differences in editing frequencies are found in some groups. The underlying processes responsible for the distribution of edited sites are largely unknown, but gene function, substitution rate, and gene conversion have been proposed to influence editing frequencies.

Results

We studied five mitochondrial genes in the monocot order Alismatales, all showing marked differences in editing frequencies among taxa. A general tendency to lose edited sites was observed in all taxa, but this tendency was particularly strong in two clades, with most of the edited sites lost in parallel in two different areas of the phylogeny. This pattern is observed in at least four of the five genes analyzed. Except in the groups that show an unusually low editing frequency, the rate of C-to-T changes in edited sites was not significantly higher that in non-edited 3^rd codon positions. This may indicate that selection is not actively removing edited sites in nine of the 12 families of the core Alismatales. In all genes but ccmB, a significant correlation was found between frequency of change in edited sites and synonymous substitution rate. In general, taxa with higher substitution rates tend to have fewer edited sites, as indicated by the phylogenetically independent correlation analyses. The elimination of edited sites in groups that lack or have reduced levels of editing could be a result of gene conversion involving a cDNA copy (retroprocessing). If so, this phenomenon could be relatively common in the Alismatales, and may have affected some groups recurrently. Indirect evidence of retroprocessing without a necessary correlation with substitution rate was found mostly in families Alismataceae and Hydrocharitaceae (e.g., groups that suffered a rapid elimination of all their edited sites, without a change in substitution rate).

Conclusions

The effects of substitution rate, selection, and/or gene conversion on the dynamics of edited sites in plant mitochondria remain poorly understood. Although we found an inverse correlation between substitution rate and editing frequency, this correlation is partially obscured by gene retroprocessing in lineages that have lost most of their edited sites. The presence of processed paralogs in plant mitochondria deserves further study, since most evidence of their occurrence is circumstantial.

RNA editing is absent in a single mitochondrial gene of Didymium iridis

An open reading frame (ORF) was found in the mitochondrial genome of the Pan2-16 strain of Didymium iridis that showed high similarity to the NADH dehydrogenase subunit 3 (nad3) gene in other organisms. So far all other typical mitochondrial genes identified in this organism require RNA editing to generate ORFs capable of directing protein synthesis. The D. iridis sequence was compared to the putative nad3 gene in the related myxomycete Physarum polycephalum, which would require editing. Based on this comparison, editing sites could be predicted for the P. polycelphalum gene that would result in the synthesis of a highly conserved ND3 protein between the two organisms. To determine the editing status of the nad3 gene in other D. iridis strains, PCR was used to amplify this region from eight other independent isolates of the A1 Central American interbreeding series. In each case a 378 base pair ORF was detected by PCR amplification and sequencing. Three patterns of sequence variation were observed; however all base substitutions were in the third codon position and silent with respect to the amino acids encoded. The distribution of the sequence variants was mapped geographically. The requirement for RNA editing in all other typical mitochondrial genes of D. iridis and P. polycephalum and the presence of RNA editing in the nad3 gene of P. polycephalum suggest that the D. iridis nad3 gene might have been edited at one time. We propose that the D. iridis nad3 gene may have lost the requirement for RNA editing by reverse transcription of an edited transcript that subsequently was inserted into the genome.

RNA editing in six mitochondrial ribosomal protein genes of Didymium iridis

Similarity searches with Didymium iridis mitochondrial genomic DNA identified six possible ribosomal protein-coding regions, however, each region contained stop codons that would need to be removed by RNA editing to produce functional transcripts. RT-PCR was used to amplify these regions from total RNA for cloning and sequencing. Six functional transcripts were verified for the following ribosomal protein genes: rpS12, rpS7, rpL2, rpS19, rpS3, and rpL16. The editing events observed, such as single C and U nucleotide insertions and a dinucleotide insertion, were consistent with previously observed editing patterns seen in D. iridis. Additionally, a new form of insertional editing, a single A insertion, was observed in a conserved region of the rpL16 gene. While the majority of codons created by editing specify hydrophobic amino acids, a greater proportion of the codons created in these hydrophilic ribosomal proteins called for positively charged amino acids in comparison to the previously characterized hydrophobic respiratory protein genes. This first report of edited soluble mitochondrial ribosomal proteins in myxomycetes expands upon the RNA editing patterns previously seen; there was: a greater proportion of created codons specifying positively charged amino acids, a shift in the codon position edited, and the insertion of single A nucleotides.