RNA editing in plants and its evolution

Untranslated yet indispensable-UTRs act as key regulators in the environmental control of gene expression

To survive and thrive in a dynamic environment, plants must continuously monitor their surroundings and adjust their development and physiology accordingly. Changes in gene expression underlie these developmental and physiological adjustments, and are traditionally attributed to widespread transcriptional reprogramming. Growing evidence, however, suggests that post-transcriptional mechanisms also play a vital role in tailoring gene expression to a plant's environment. Untranslated regions (UTRs) act as regulatory hubs for post-transcriptional control, harbouring cis-elements that affect an mRNA's processing, localization, translation, and stability, and thereby tune the abundance of the encoded protein. Here, we review recent advances made in understanding the critical function UTRs exert in the post-transcriptional control of gene expression in the context of a plant's abiotic environment. We summarize the molecular mechanisms at play, present examples of UTR-controlled signalling cascades, and discuss the potential that resides within UTRs to render plants more resilient to a changing climate.

Learning from the Codon Table: Convergent Recoding Provides Novel Understanding on the Evolution of A-to-I RNA Editing

Adenosine-to-inosine (A-to-I) RNA editing recodes the genetic information. Apart from diversifying the proteome, another tempting advantage of RNA recoding is to correct deleterious DNA mutation and restore ancestral allele. Solid evidences for beneficial restorative editing are very rare in animals. By searching for "convergent recoding" under a phylogenetic context, we proposed this term for judging the potential restorative functions of particular editing site. For the well-known mammalian Gln>Arg (Q>R) recoding site, its ancestral state in vertebrate genomes was the pre-editing Gln, and all 470 available mammalian genomes strictly avoid other three equivalent ways to achieve Arg in protein. The absence of convergent recoding from His>Arg, or synonymous mutations on Gln codons, could be attributed to the strong maintenance on editing motif and structure, but the absence of direct A-to-G mutation is extremely unexpected. With similar ideas, we found cases of convergent recoding in Drosophila genus, reducing the possibility of their restorative function. In summary, we defined an interesting scenario of convergent recoding, the occurrence of which could be used as preliminary judgements for whether a recoding site has a sole restorative role. Our work provides novel insights to the natural selection and evolution of RNA editing.

The HPE1 RNA-binding protein modulates chloroplast RNA editing to promote photosynthesis under cold stress in Arabidopsis

Cold stress has severe negative consequences for plant growth and crop yield. Here, we report that an Arabidopsis thaliana mutant that lacks the HPE1 gene, which encodes an RNA-binding protein, maintains higher photosynthetic activity under cold stress, together with higher accumulation of thylakoid proteins. We showed that HPE1 interacts with MORF2 and MORF9 and thereby mediates RNA editing in chloroplasts. Loss of HPE1 function increased the editing efficiency at four RNA editing sites, rpoC-488, ndhB-149, ndhB-746 and matK-706, under cold stress and altered the expression of nuclear photosynthesis-related genes and cold-responsive genes. We propose that HPE1-mediated RNA editing acts as a trigger for retrograde signaling that affects photosynthesis under cold stress.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

Insights into structure, codon usage, repeats, and RNA editing of the complete mitochondrial genome of Perilla frutescens (Lamiaceae)

Perilla frutescens (L.) Britton, a member of the Lamiaceae family, stands out as a versatile plant highly valued for its unique aroma and medicinal properties. Additionally, P. frutescens seeds are rich in Îś-linolenic acid, holding substantial economic importance. While the nuclear and chloroplast genomes of P. frutescens have already been documented, the complete mitochondrial genome sequence remains unreported. To this end, the sequencing, annotation, and assembly of the entire Mitochondrial genome of P. frutescens were hereby conducted using a combination of Illumina and PacBio data. The assembled P. frutescens mitochondrial genome spanned 299,551 bp and exhibited a typical circular structure, involving a GC content of 45.23%. Within the genome, a total of 59 unique genes were identified, encompassing 37 protein-coding genes, 20 tRNA genes, and 2 rRNA genes. Additionally, 18 introns were observed in 8 protein-coding genes. Notably, the codons of the P. frutescens mitochondrial genome displayed a notable A/T bias. The analysis also revealed 293 dispersed repeat sequences, 77 simple sequence repeats (SSRs), and 6 tandem repeat sequences. Moreover, RNA editing sites preferentially produced leucine at amino acid editing sites. Furthermore, 70 sequence fragments (12,680 bp) having been transferred from the chloroplast to the mitochondrial genome were identified, accounting for 4.23% of the entire mitochondrial genome. Phylogenetic analysis indicated that among Lamiaceae plants, P. frutescens is most closely related to Salvia miltiorrhiza and Platostoma chinense. Meanwhile, inter-species Ka/Ks results suggested that Ka/Ks < 1 for 28 PCGs, indicating that these genes were evolving under purifying selection. Overall, this study enriches the mitochondrial genome data for P. frutescens and forges a theoretical foundation for future molecular breeding research.

MORF9-dependent specific plastid RNA editing inhibits root growth under sugar starvation in Arabidopsis

Multiple organellar RNA editing factor (MORF) complex was shown to be highly associated with C-to-U RNA editing of vascular plant editosome. However, mechanisms by which MORF9-dependent plastid RNA editing controls plant development and responses to environmental alteration remain obscure. In this study, we found that loss of MORF9 function impaired PSII efficiency, NDH activity, and carbohydrate production, rapidly promoted nuclear gene expression including sucrose transporter and sugar/energy responsive genes, and attenuated root growth under sugar starvation conditions. Sugar repletion increased MORF9 and MORF2 expression in wild-type seedlings and reduced RNA editing of matK-706, accD-794, ndhD-383 and ndhF-290 in the morf9 mutant. RNA editing efficiency of ndhD-383 and ndhF-290 sites was diminished in the gin2/morf9 double mutants, and that of matK-706, accD-794, ndhD-383 and ndhF-290 sites were significantly diminished in the snrk1/morf9 double mutants. In contrast, overexpressing HXK1 or SnRK1 promoted RNA editing rate of matK-706, accD-794, ndhD-383 and ndhF-290 in leaves of morf9 mutants, suggesting that HXK1 partially impacts MORF9 mediated ndhD-383 and ndhF-290 editing, while SnRK1 may only affect MORF9-mediated ndhF-290 site editing. Collectively, these findings suggest that sugar and/or its intermediary metabolites impair MORF9-dependent plastid RNA editing resulting in derangements of plant root development.

Is RNA editing truly absent in the complex thalloid liverworts (Marchantiopsida)? Evidence of extensive RNA editing from Cyathodium cavernarum

RNA editing is a crucial modification in plants' organellar transcripts that converts cytidine to uridine (C-to-U; and sometimes uridine to cytidine) in RNA molecules. This post-transcriptional process is controlled by the PLS-class protein with a DYW domain, which belongs to the pentatricopeptide repeat (PPR) protein family. RNA editing is widespread in land plants; however, complex thalloid liverworts (Marchantiopsida) are the only group reported to lack both RNA editing and DYW-PPR protein. The liverwort Cyathodium cavernarum (Marchantiopsida, Cyathodiaceae), typically found in cave habitats, was newly found to have 129 C-to-U RNA editing sites in its chloroplast and 172 sites in its mitochondria. The Cyathodium genus, specifically C. cavernarum, has a large number of PPR editing factor genes, including 251 DYW-type PPR proteins. These DYW-type PPR proteins may be responsible for C-to-U RNA editing in C. cavernarum. Cyathodium cavernarum possesses both PPR DYW proteins and RNA editing. Our analysis suggests that the remarkable RNA editing capability of C. cavernarum may have been acquired alongside the emergence of DYW-type PPR editing factors. These findings provide insight into the evolutionary pattern of RNA editing in land plants.

Comprehensive identification, characterization, and expression analysis of the MORF gene family in Brassica napus

BACKGROUND

RNA editing in chloroplast and mitochondrion transcripts of plants is an important type of post-transcriptional RNA modification in which members of the multiple organellar RNA editing factor gene family (MORF) play a crucial role. However, a systematic identification and characterization of MORF members in Brassica napus is still lacking.

RESULTS

In this study, a total of 43 MORF genes were identified from the genome of the Brassica napus cultivar "Zhongshuang 11". The Brassica napus MORF (BnMORF) family members were divided into three groups through phylogenetic analysis. BnMORF genes distributed on 14 chromosomes and expanded due to segmental duplication and whole genome duplication repetitions. The majority of BnMORF proteins were predicted to be localized to mitochondria and chloroplasts. The promoter cis-regulatory element analysis, spatial-temporal expression profiling, and co-expression network of BnMORF genes indicated the involvement of BnMORF genes in stress and phytohormone responses, as well as growth and development.

CONCLUSION

This study provides a comprehensive analysis of BnMORF genes and lays a foundation for further exploring their physiological functions in Brassica napus.

Exploring functional conservation in silico: a new machine learning approach to RNA-editing

Around 50 years ago, molecular biology opened the path to understand changes in forms, adaptations, complexity, or the basis of human diseases through myriads of reports on gene birth, gene duplication, gene expression regulation, and splicing regulation, among other relevant mechanisms behind gene function. Here, with the advent of big data and artificial intelligence (AI), we focus on an elusive and intriguing mechanism of gene function regulation, RNA editing, in which a single nucleotide from an RNA molecule is changed, with a remarkable impact in the increase of the complexity of the transcriptome and proteome. We present a new generation approach to assess the functional conservation of the RNA-editing targeting mechanism using two AI learning algorithms, random forest (RF) and bidirectional long short-term memory (biLSTM) neural networks with an attention layer. These algorithms, combined with RNA-editing data coming from databases and variant calling from same-individual RNA and DNA-seq experiments from different species, allowed us to predict RNA-editing events using both primary sequence and secondary structure. Then, we devised a method for assessing conservation or divergence in the molecular mechanisms of editing completely in silico: the cross-testing analysis. This novel method not only helps to understand the conservation of the editing mechanism through evolution but could set the basis for achieving a better understanding of the adenosine-targeting mechanism in other fields.

PPR596 Is Required for nad2 Intron Splicing and Complex I Biogenesis in Arabidopsis

Mitochondria are essential organelles that generate energy via oxidative phosphorylation. Plant mitochondrial genome encodes some of the respiratory complex subunits, and these transcripts require accurate processing, including C-to-U RNA editing and intron splicing. Pentatricopeptide repeats (PPR) proteins are involved in various organellar RNA processing events. PPR596, a P-type PPR protein, was previously identified to function in the C-to-U editing of mitochondrial rps3 transcripts in Arabidopsis. Here, we demonstrate that PPR596 functions in the cis-splicing of nad2 intron 3 in mitochondria. Loss of the PPR596 function affects the editing at rps3eU1344SS, impairs nad2 intron 3 splicing and reduces the mitochondrial complex I's assembly and activity, while inducing alternative oxidase (AOX) gene expression. This defect in nad2 intron splicing provides a plausible explanation for the slow growth of the ppr595 mutants. Although a few P-type PPR proteins are involved in RNA C-to-U editing, our results suggest that the primary function of PPR596 is intron splicing.

Detection and editing of the updated Arabidopsis plastid- and mitochondrial-encoded proteomes through PeptideAtlas

Arabidopsis (Arabidopsis thaliana) ecotype Col-0 has plastid and mitochondrial genomes encoding over 100 proteins. Public databases (e.g. Araport11) have redundancy and discrepancies in gene identifiers for these organelle-encoded proteins. RNA editing results in changes to specific amino acid residues or creation of start and stop codons for many of these proteins, but the impact of RNA editing at the protein level is largely unexplored due to the complexities of detection. Here, we assembled the nonredundant set of identifiers, their correct protein sequences, and 452 predicted nonsynonymous editing sites of which 56 are edited at lower frequency. We then determined accumulation of edited and/or unedited proteoforms by searching ∼259 million raw tandem MS spectra from ProteomeXchange, which is part of PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/). We identified all mitochondrial proteins and all except 3 plastid-encoded proteins (NdhG/Ndh6, PsbM, and Rps16), but no proteins predicted from the 4 ORFs were identified. We suggest that Rps16 and 3 of the ORFs are pseudogenes. Detection frequencies for each edit site and type of edit (e.g. S to L/F) were determined at the protein level, cross-referenced against the metadata (e.g. tissue), and evaluated for technical detection challenges. We detected 167 predicted edit sites at the proteome level. Minor frequency sites were edited at low frequency at the protein level except for cytochrome C biogenesis 382 at residue 124 (Ccb382-124). Major frequency sites (>50% editing of RNA) only accumulated in edited form (>98% to 100% edited) at the protein level, with the exception of Rpl5-22. We conclude that RNA editing for major editing sites is required for stable protein accumulation.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns

BACKGROUND

Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events.

RESULTS

In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species.

CONCLUSIONS

Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

YUCCA2 (YUC2)-Mediated 3-Indoleacetic Acid (IAA) Biosynthesis Regulates Chloroplast RNA Editing by Relieving the Auxin Response Factor 1 (ARF1)-Dependent Inhibition of Editing Factors in Arabidopsis thaliana

Although recent research progress on the abundant C-to-U RNA editing events in plant chloroplasts and mitochondria has uncovered many recognition factors and their molecular mechanisms, the intrinsic regulation of RNA editing within plants remains largely unknown. This study aimed to establish a regulatory relationship in Arabidopsis between the plant hormone auxin and chloroplast RNA editing. We first analyzed auxin response elements (AuxREs) present within promoters of chloroplast editing factors reported to date. We found that each has more than one AuxRE, suggesting a potential regulatory role of auxin in their expression. Further investigation unveiled that the depletion of auxin synthesis gene YUC2 reduces the expression of several editing factors. However, in yuc2 mutants, only the expression of CRR4, DYW1, ISE2, and ECD1 editing factors and the editing efficiency of their corresponding editing sites, ndhD-2 and rps14-149, were simultaneously suppressed. In addition, exogenous IAA and the overexpression of YUC2 enhanced the expression of these editing factors and the editing efficiency at the ndhD-2 and rps14-149 sites. These results suggested a direct effect of auxin upon the editing of the ndhD-2 and rps14-149 sites through the modulation of the expression of the editing factors. We further demonstrated that ARF1, a downstream transcription factor in the auxin-signaling pathway, could directly bind to and inactivate the promoters of CRR4, DYW1, and ISE2 in a dual-luciferase reporter system, thereby inhibiting their expression. Moreover, the overexpression of ARF1 in Arabidopsis significantly reduced the expression of the three editing factors and the editing efficiency at the ndhD-2 and rps14-149 sites. These data suggest that YUC2-mediated auxin biosynthesis governs the RNA-editing process through the ARF1-dependent signal transduction pathway.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

The emerging role of epitranscriptome in shaping stress responses in plants

KEY MESSAGE

RNA modifications and editing changes constitute 'epitranscriptome' and are crucial in regulating the development and stress response in plants. Exploration of the epitranscriptome and associated machinery would facilitate the engineering of stress tolerance in crops. RNA editing and modifications post-transcriptionally decorate almost all classes of cellular RNAs, including tRNAs, rRNAs, snRNAs, lncRNAs and mRNAs, with more than 170 known modifications, among which m6A, Ψ, m5C, 8-OHG and C-to-U editing are the most abundant. Together, these modifications constitute the "epitranscriptome", and contribute to changes in several RNA attributes, thus providing an additional structural and functional diversification to the "cellular messages" and adding another layer of gene regulation in organisms, including plants. Numerous evidences suggest that RNA modifications have a widespread impact on plant development as well as in regulating the response of plants to abiotic and biotic stresses. High-throughput sequencing studies demonstrate that the landscapes of m6A, m5C, Am, Cm, C-to-U, U-to-G, and A-to-I editing are remarkably dynamic during stress conditions in plants. GO analysis of transcripts enriched in Ψ, m6A and m5C modifications have identified bonafide components of stress regulatory pathways. Furthermore, significant alterations in the expression pattern of genes encoding writers, readers, and erasers of certain modifications have been documented when plants are grown in challenging environments. Notably, manipulating the expression levels of a few components of RNA editing machinery markedly influenced the stress tolerance in plants. We provide updated information on the current understanding on the contribution of RNA modifications in shaping the stress responses in plants. Unraveling of the epitranscriptome has opened new avenues for designing crops with enhanced productivity and stress resilience in view of global climate change.

Comparative Analysis of Tylosema esculentum Mitochondrial DNA Revealed Two Distinct Genome Structures

Tylosema esculentum, commonly known as the marama bean, is an underutilized legume with nutritious seeds, holding potential to enhance food security in southern Africa due to its resilience to prolonged drought and heat. To promote the selection of this agronomically valuable germplasm, this study assembled and compared the mitogenomes of 84 marama individuals, identifying variations in genome structure, single-nucleotide polymorphisms (SNPs), insertions/deletions (indels), heteroplasmy, and horizontal transfer. Two distinct germplasms were identified, and a novel mitogenome structure consisting of three circular molecules and one long linear chromosome was discovered. The structural variation led to an increased copy number of specific genes, nad5, nad9, rrnS, rrn5, trnC, and trnfM. The two mitogenomes also exhibited differences at 230 loci, with only one notable nonsynonymous substitution in the matR gene. Heteroplasmy was concentrated at certain loci on chromosome LS1 (OK638188). Moreover, the marama mitogenome contained an over 9 kb insertion of cpDNA, originating from chloroplast genomes, but had accumulated mutations and lost gene functionality. The evolutionary and comparative genomics analysis indicated that mitogenome divergence in marama might not be solely constrained by geographical factors. Additionally, marama, as a member from the Cercidoideae subfamily, tends to possess a more complete set of mitochondrial genes than Faboideae legumes.

Chloroplast gene expression: Recent advances and perspectives

Chloroplasts evolved from an ancient cyanobacterial endosymbiont more than 1.5 billion years ago. During subsequent coevolution with the nuclear genome, the chloroplast genome has remained independent, albeit strongly reduced, with its own transcriptional machinery and distinct features, such as chloroplast-specific innovations in gene expression and complicated post-transcriptional processing. Light activates the expression of chloroplast genes via mechanisms that optimize photosynthesis, minimize photodamage, and prioritize energy investments. Over the past few years, studies have moved from describing phases of chloroplast gene expression to exploring the underlying mechanisms. In this review, we focus on recent advances and emerging principles that govern chloroplast gene expression in land plants. We discuss engineering of pentatricopeptide repeat proteins and its biotechnological effects on chloroplast RNA research; new techniques for characterizing the molecular mechanisms of chloroplast gene expression; and important aspects of chloroplast gene expression for improving crop yield and stress tolerance. We also discuss biological and mechanistic questions that remain to be answered in the future.

The DYW domain of RARE1 plays an indispensable role in regulating accD-C794 RNA editing in Arabidopsis thaliana

The Arabidopsis pentatricopeptide repeat (PPR) proteins, required for accD RNA editing 1 (RARE1) and early chloroplast biogenesis 2 (AtECB2), each contain a DYW domain deemed essential for cytosine deamination at the accD-C794 RNA editing site in chloroplasts. Complementation assays using the rare1 mutant investigate the correlation between these PPRs and their respective DYW domain functions in RNA editing of accD-C794. The results demonstrate that the coding sequence of AtECB2 cannot replace that of RARE1. Moreover, rare1 mutants complemented with DYW-deleted RARE1 failed to recover the RNA editing of accD-C794 even in the presence of the highly similar DYW domain of the AtECB2 protein. These findings indicate that RARE1 and AtECB2 possess divergent roles in RNA editing, with specificity for accD-C794 directly attributable to DYW domain within RARE1. Structural modeling data suggest this functioning pertains to a local α-helical motif that residues slightly N-terminal to the consensus glutamate and CXXCH motif in the DYW domain for cytidine deamination during C-to-U editing by RARE1 that is absent within AtECB2.

Benchmarking RNA Editing Detection Tools

RNA, like DNA and proteins, can undergo modifications. To date, over 170 RNA modifications have been identified, leading to the emergence of a new research area known as epitranscriptomics. RNA editing is the most frequent RNA modification in mammalian transcriptomes, and two types have been identified: (1) the most frequent, adenosine to inosine (A-to-I); and (2) the less frequent, cysteine to uracil (C-to-U) RNA editing. Unlike other epitranscriptomic marks, RNA editing can be readily detected from RNA sequencing (RNA-seq) data without any chemical conversions of RNA before sequencing library preparation. Furthermore, analyzing RNA editing patterns from transcriptomic data provides an additional layer of information about the epitranscriptome. As the significance of epitranscriptomics, particularly RNA editing, gains recognition in various fields of biology and medicine, there is a growing interest in detecting RNA editing sites (RES) by analyzing RNA-seq data. To cope with this increased interest, several bioinformatic tools are available. However, each tool has its advantages and disadvantages, which makes the choice of the most appropriate tool for bench scientists and clinicians difficult. Here, we have benchmarked bioinformatic tools to detect RES from RNA-seq data. We provide a comprehensive view of each tool and its performance using previously published RNA-seq data to suggest recommendations on the most appropriate for utilization in future studies.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

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

Deciphering the enigma of RNA editing in the ATP1_alpha subunit of ATP synthase in Triticum aestivum

There is evidence that RNA editing is related to plant cellular stress as well as electron transport organelles, such as mitochondria. The mitochondrial atp1 gene encodes the alpha-subunit of Atp synthase. Control as well as two periods of drought stress treatments were analyzed in the cDNAs generated from the mitochondrial atp1 gene of two cultivars of Triticum aestivum [Giza 168 (G168) and Gemmiza 10 (GM10)]. Following RNA-seq data assembly, atp1 cDNAs from the control (acc. no. OQ129415), 2-hour (acc. no. OQ129416), and 12-hour (acc. no. OQ129417) time points of the T. aestivum cultivar G168 were obtained. Control (acc. no. OQ129419), 2-hour (acc. no. OQ129420), and 12-hour (acc. no. OQ129421) samples all included reconstructed atp1 transcripts from Gemmiza 10. Atp1 transcripts were assembled using the wheat atp1 gene (acc. no. NC_036024). RNA-seq raw data was utilized to identify 11 RNA editing sites in atp1 in the tolerant cultivar Giza168 and 6 in the sensitive cultivar Gemmiza10. The significant difference in RNA editing observed between control and drought stress conditions in sites led to synonymous amino acids. This led to no change in tertiary structure between tolerant and sensitive cultivars. But the change was focused between produced protein and its correspondence sequence on DNA.

Progress, challenge and prospect of plant plastome annotation

The plastome (plastid genome) represents an indispensable molecular data source for studying phylogeny and evolution in plants. Although the plastome size is much smaller than that of nuclear genome, and multiple plastome annotation tools have been specifically developed, accurate annotation of plastomes is still a challenging task. Different plastome annotation tools apply different principles and workflows, and annotation errors frequently occur in published plastomes and those issued in GenBank. It is therefore timely to compare available annotation tools and establish standards for plastome annotation. In this review, we review the basic characteristics of plastomes, trends in the publication of new plastomes, the annotation principles and application of major plastome annotation tools, and common errors in plastome annotation. We propose possible methods to judge pseudogenes and RNA-editing genes, jointly consider sequence similarity, customed algorithms, conserved domain or protein structure. We also propose the necessity of establishing a database of reference plastomes with standardized annotations, and put forward a set of quantitative standards for evaluating plastome annotation quality for the scientific community. In addition, we discuss how to generate standardized GenBank annotation flatfiles for submission and downstream analysis. Finally, we prospect future technologies for plastome annotation integrating plastome annotation approaches with diverse evidences and algorithms of nuclear genome annotation tools. This review will help researchers more efficiently use available tools to achieve high-quality plastome annotation, and promote the process of standardized annotation of the plastome.

The ATP Synthase γ Subunit ATPC1 Regulates RNA Editing in Chloroplasts

RNA editing is the process of modifying RNA molecules by inserting, deleting, or substituting nucleotides. In flowering plants, RNA editing occurs predominantly in RNAs encoded by the organellar genomes of mitochondria and chloroplasts, and the main type of editing involves the substitution of cytidine with uridine at specific sites. Abnormal RNA editing in plants can affect gene expression, organelle function, plant growth, and reproduction. In this study, we report that ATPC1, the gamma subunit of ATP synthase in Arabidopsis chloroplasts, has an unexpected role in the regulation of editing at multiple sites of plastid RNAs. The loss of function of ATPC1 severely arrests chloroplast development, causing a pale-green phenotype and early seedling lethality. Disruption of ATPC1 increases the editing of matK-640, rps12-i-58, atpH-3'UTR-13210, and ycf2-as-91535 sites while decreasing the editing of rpl23-89, rpoA-200, rpoC1-488, and ndhD-2 sites. We further show that ATPC1 participates in RNA editing by interacting with known multiple-site chloroplast RNA editing factors, including MORFs, ORRM1, and OZ1. The transcriptome in the atpc1 mutant is profoundly affected, with a pattern of defective expression of chloroplast development-related genes. These results reveal that the ATP synthase γ subunit ATPC1 is involved in multiple-site RNA editing in Arabidopsis chloroplasts.

Transcriptomics and Metabolomics Analysis Provides Insight into Leaf Color and Photosynthesis Variation of the Yellow-Green Leaf Mutant of Hami Melon (Cucumis melo L.)

Leaf color mutants are ideal materials for studying the regulatory mechanism of chloroplast development and photosynthesis. We isolated a cucumis melo spontaneous mutant (MT), which showed yellow-green leaf phenotype in the whole growing period and could be inherited stably. We compared its leaves with the wild type (WT) in terms of cytology, physiology, transcriptome and metabolism. The results showed that the thylakoid grana lamellae of MT were loosely arranged and fewer in number than WT. Physiological experiments also showed that MT had less chlorophyll content and more accumulation of reactive oxygen species (ROS) than WT. Furthermore, the activity of several key enzymes in C₄ photosynthetic carbon assimilation pathway was more enhanced in MT than WT. Transcriptomic and metabolomic analyses showed that differential expression genes and differentially accumulated metabolites in MT were mainly co-enriched in the pathways related to photosystem-antenna proteins, central carbon metabolism, glutathione metabolism, phenylpropanoid biosynthesis and flavonoid metabolism. We also analyzed several key proteins in photosynthesis and chloroplast transport by Western blot. In summary, the results may provide a new insight into the understanding of how plants respond to the impaired photosynthesis by regulating chloroplast development and photosynthetic carbon assimilation pathways.

Young Leaf White Stripe encodes a P-type PPR protein required for chloroplast development

Pentatricopeptide repeat (PPR) proteins function in post-transcriptional regulation of organellar gene expression. Although several PPR proteins are known to function in chloroplast development in rice (Oryza sativa), the detailed molecular functions of many PPR proteins remain unclear. Here, we characterized a rice young leaf white stripe (ylws) mutant, which has defective chloroplast development during early seedling growth. Map-based cloning revealed that YLWS encodes a novel P-type chloroplast-targeted PPR protein with 11 PPR motifs. Further expression analyses showed that many nuclear- and plastid-encoded genes in the ylws mutant were significantly changed at the RNA and protein levels. The ylws mutant was impaired in chloroplast ribosome biogenesis and chloroplast development under low-temperature conditions. The ylws mutation causes defects in the splicing of atpF, ndhA, rpl2, and rps12, and editing of ndhA, ndhB, and rps14 transcripts. YLWS directly binds to specific sites in the atpF, ndhA, and rpl2 pre-mRNAs. Our results suggest that YLWS participates in chloroplast RNA group II intron splicing and plays an important role in chloroplast development during early leaf development.

Low copy numbers for mitochondrial DNA moderates the strength of nuclear-cytoplasmic incompatibility in plants

Plant cells contain only small amounts of mitochondrial DNA (mtDNA), with the genomic information shared among multiple mitochondria. The biological relevance and molecular mechanism underlying this hallmark of plant cells has been unclear. Here, we report that Arabidopsis thaliana plants exhibited significantly reduced growth and mitochondrial dysfunction when the mtDNA copy number was increased to the degree that each mitochondrion possessed DNA. The amounts of mitochondrion-encoded transcripts increased several fold in the presence of elevated mtDNA levels. However, the efficiency of RNA editing decreased with this excess of mitochondrion-encoded transcripts, resulting in impaired assembly of mitochondrial complexes containing mtDNA-encoded subunits, such as respiratory complexes I and IV. These observations indicate the occurrence of nuclear-mitochondrial incompatibility in the cells with increased amounts of mtDNA and provide an initial answer to the fundamental question of why plant cells have much lower mtDNA levels than animal cells. We propose that keeping mtDNA levels low moderates nuclear-mitochondrial incompatibility and that this may be a crucial factor driving plant cells to restrict the copy numbers of mtDNA.

On the origin and evolution of RNA editing in metazoans

Extensive adenosine-to-inosine (A-to-I) editing of nuclear-transcribed mRNAs is the hallmark of metazoan transcriptional regulation. Here, by profiling the RNA editomes of 22 species that cover major groups of Holozoa, we provide substantial evidence supporting A-to-I mRNA editing as a regulatory innovation originating in the last common ancestor of extant metazoans. This ancient biochemistry process is preserved in most extant metazoan phyla and primarily targets endogenous double-stranded RNA (dsRNA) formed by evolutionarily young repeats. We also find intermolecular pairing of sense-antisense transcripts as an important mechanism for forming dsRNA substrates for A-to-I editing in some but not all lineages. Likewise, recoding editing is rarely shared across lineages but preferentially targets genes involved in neural and cytoskeleton systems in bilaterians. We conclude that metazoan A-to-I editing might first emerge as a safeguard mechanism against repeat-derived dsRNA and was later co-opted into diverse biological processes due to its mutagenic nature.

Polymeric structure of the Cannabis sativa L. mitochondrial genome identified with an assembly graph model

Cannabis sativa L. belongs to the family Cannabaceae in Rosales. It has been widely used as medicines, building materials, and textiles. Elucidating its genome is critical for molecular breeding and synthetic biology study. Many studies have shown that the mitochondrial genomes (mitogenomes) and even chloroplast genomes (plastomes) had complex polymeric structures. Using the Nanopore sequencing platform, we sequenced, assembled, and analyzed its mitogenome and plastome. The resulting unitig graph suggested that the mitogenome had a complex polymeric structure. However, a gap-free, circular sequence was further assembled from the unitig graph. In contrast, a circular sequence representing the plastome was obtained. The mitogenome major conformation was 415,837 bp long, and the plastome was 153,927 bp long. To test if the repeat sequences promote recombination, which corresponds to the branch points in the structure, we tested the sequences around repeats by long-read mapping. Among 208 pairs of predicted repeats, the mapping results supported the presence of cross-over around 25 pairs of repeats. Subsequent PCR amplification confirmed the presence of cross-over around 15 of the 25 repeats. By comparing the mitogenome and plastome sequences, we identified 19 mitochondria plastid DNAs, including seven complete genes (trnW-CCA, trnP-UGG, psbJ, trnN-GUU, trnD-GUC, trnH-GUG, trnM-CAU) and nine gene fragments. Furthermore, the selective pressure analysis results showed that five genes (atp1, ccmB, ccmC, cox1, nad7) had 19 positively selected sites. Lastly, we predicted 28 RNA editing sites. A total of 8 RNA editing sites located in the coding regions were successfully validated by PCR amplification and Sanger sequencing, of which four were synonymous, and four were nonsynonymous. In particular, the RNA editing events appeared to be tissue-specific in C. sativa mitogenome. In summary, we have confirmed the major confirmation of C. sativa mitogenome and characterized its structural features in detail. These results provide critical information for future variety breeding and resource development for C. sativa.

SlRIP1b is a global organellar RNA editing factor, required for normal fruit development in tomato plants

RNA editing in plant organelles involves numerous C-U conversions, which often restore evolutionarily conserved codons and may generate new translation initiation and termination codons. These RNA maturation events rely on a subset of nuclear-encoded protein cofactors. Here, we provide evidence of the role of SlRIP1b on RNA editing of mitochondrial transcripts in tomato (Solanum lycopersicum) plants. SlRIP1b is a RIP/MORF protein that was originally identified as an interacting partner of the organellar editing factor SlORRM4. Mutants of SlRIP1b, obtained by CRISPR/Cas9 strategy, exhibited abnormal carpel development and grew into fruit with more locules. RNA-sequencing revealed that SlRIP1b affects the C-U editing of numerous mitochondrial pre-RNA transcripts and in particular altered RNA editing of various cytochrome c maturation (CCM)-related genes. The slrip1b mutants display increased H₂ O₂ and aberrant mitochondrial morphologies, which are associated with defects in cytochrome c biosynthesis and assembly of respiratory complex III. Taken together, our results indicate that SlRIP1b is a global editing factor that plays a key role in CCM and oxidative phosphorylation system biogenesis during fruit development in tomato plants. These data provide important insights into the molecular roles of organellar RNA editing factors during fruit development.

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.

Lactobacillus for ribosome peptide editing cancer

This study reviews newly discovered insect peptide point mutations as new possible cancer research targets. To interpret newly discovered peptide point mutations in insects as new possible cancer research targets, we focused on the numerous peptide changes found in the 'CSP' family on the sex pheromone gland of the female silkworm moth Bombyx mori. We predict that the Bombyx peptide modifications will have a significant effect on cancer CUP (cancers of unknown primary) therapy and that bacterial peptide editing techniques, specifically Lactobacillus combined to CRISPR, will be used to regulate ribosomes and treat cancer in humans.

Isolation and Characterization of SPOTTED LEAF42 Encoding a Porphobilinogen Deaminase in Rice

The formation and development of chloroplasts play a vital role in the breeding of high-yield rice (Oryza sativa L.). Porphobilinogen deaminases (PBGDs) act in the early stage of chlorophyll and heme biosynthesis. However, the role of PBGDs in chloroplast development and chlorophyll production remains elusive in rice. Here, we identified the spotted leaf 42 (spl42) mutant, which exhibited a reddish-brown spotted leaf phenotype. The mutant showed a significantly lower chlorophyll content, abnormal thylakoid morphology, and elevated activities of reactive oxygen species (ROS)-scavenging enzymes. Consistently, multiple genes related to chloroplast development and chlorophyll biosynthesis were significantly down-regulated, whereas many genes involved in leaf senescence, ROS production, and defense responses were upregulated in the spl42 mutant. Map-based cloning revealed that SPL42 encodes a PBGD. A C-to-T base substitution occurred in spl42, resulting in an amino acid change and significantly reduced PBGD enzyme activity. SPL42 targets to the chloroplast and interacts with the multiple organelle RNA editing factors (MORFs) OsMORF8-1 and OsMORF8-2 to affect RNA editing. The identification and characterization of spl42 helps in elucidating the molecular mechanisms associated with chlorophyll synthesis and RNA editing in rice.

Characterization and phylogenetic analysis of the complete mitochondrial genome sequence of Photinia serratifolia

Plant mitochondrial genomes (mitogenomes) are a valuable source of genetic information for a better understanding of phylogenetic relationships. However, no mitogenome of any species in the genus of Photinia has been reported. In this study, using NGS sequencing, we reported the mitogenome assembly and annotation of Photinia serratifolia, which is 473,579 bp in length, contains 38 protein-coding genes, 23 tRNAs, and 6 rRNAs, with 61 genes have no introns. The rps2 and rps11 genes are missing in the P. serratifolia mitogenome. Although there are more editing sites (488) in the P. serratifolia mitogenome than in most angiosperms, fewer editing types were found in the P. serratifolia mitogenome, showing a clear bias in RNA-editing. Phylogenetic analysis based on the mitogenomes of P. serratifolia and 8 other taxa of the Rosaceae family reflected the exact evolutionary and taxonomic status of P. serratifolia. However, Ka/Ks analysis revealed that 72.69% of the protein-coding genes in the P. serratifolia mitogenome had undergone negative selections, reflecting the importance of those genes in the P. serratifolia mitogenome. Collectively, these results will provide valuable information for the evolution of P. serratifolia and provide insight into the evolutionary relationships within Photinia and the Rosaceae family.

Occurrence of RNA post-transcriptional modifications in plant viruses and viroids and their correlation with structural and functional features

Post-transcriptional modifications of RNA bases are widespread across all the tree of life and have been linked to RNA maturation, stability, and molecular interactions. RNA modifications have been extensively described in endogenous eukaryotic mRNAs, however, little is known about the presence of RNA modifications in plant viral and subviral RNAs. Here, we used a computational approach to infer RNA modifications in plant-pathogenic viruses and viroids using high-throughput annotation of modified ribonucleotides (HAMR), a software that predicts modified ribonucleotides using high-throughput RNA sequencing data. We analyzed datasets from representative members of different plant viruses and viroids and compared them to plant-endogenous mRNAs. Our approach was able to predict potential RNA chemical modifications (RCMs) in all analyzed pathogens. We found that both DNA and RNA viruses presented a wide range of RCM proportions while viroids had lowest values. Furthermore, we found that for viruses with segmented genomes, some genomic RNAs had a higher proportion of RCM. Interestingly, nuclear-replicating viroids showed most of the predicted modifications located in the pathogenesis region, pointing towards a possible functional role of RCMs in their infectious cycle. Thus, our results strongly suggest that plant viral and subviral RNAs might contain a variety of previously unreported RNA modifications, thus opening a new perspective in the multifaceted process of plant-pathogen interactions.

Genome-wide identification and expression analysis of peach multiple organellar RNA editing factors reveals the roles of RNA editing in plant immunity

BACKGROUND

Multiple organellar RNA editing factor (MORF) genes play key roles in chloroplast developmental processes by mediating RNA editing of Cytosine-to-Uracil conversion. However, the function of MORF genes in peach (Prunus persica), a perennial horticultural crop species of Rosaceae, is still not well known, particularly the resistance to biotic and abiotic stresses that threaten peach yield seriously.

RESULTS

In this study, to reveal the regulatory roles of RNA editing in plant immunity, we implemented genome-wide analysis of peach MORF (PpMORF) genes in response to biotic and abiotic stresses. The chromosomal and subcellular location analysis showed that the identified seven PpMORF genes distributed on three peach chromosomes were mainly localized in the mitochondria and chloroplast. All the PpMORF genes were classified into six groups and one pair of PpMORF genes was tandemly duplicated. Based on the meta-analysis of two types of public RNA-seq data under different treatments (biotic and abiotic stresses), we observed down-regulated expression of PpMORF genes and reduced chloroplast RNA editing, especially the different response of PpMORF2 and PpMORF9 to pathogens infection between resistant and susceptible peach varieties, indicating the roles of MORF genes in stress response by modulating the RNA editing extent in plant immunity. Three upstream transcription factors (MYB3R-1, ZAT10, HSFB3) were identified under both stresses, they may regulate resistance adaption by modulating the PpMORF gene expression.

CONCLUSION

These results provided the foundation for further analyses of the functions of MORF genes, in particular the roles of RNA editing in plant immunity. In addition, our findings will be conducive to clarifying the resistance mechanisms in peaches and open up avenues for breeding new cultivars with high resistance.

RNA Editing in Chloroplast: Advancements and Opportunities

Many eukaryotic and prokaryotic organisms employ RNA editing (insertion, deletion, or conversion) as a post-transcriptional modification mechanism. RNA editing events are common in these organelles of plants and have gained particular attention due to their role in the development and growth of plants, as well as their ability to cope with abiotic stress. Owing to rapid developments in sequencing technologies and data analysis methods, such editing sites are being accurately predicted, and many factors that influence RNA editing are being discovered. The mechanism and role of the pentatricopeptide repeat protein family of proteins in RNA editing are being uncovered with the growing realization of accessory proteins that might help these proteins. This review will discuss the role and type of RNA editing events in plants with an emphasis on chloroplast RNA editing, involved factors, gaps in knowledge, and future outlooks.

Exploring the RNA Editing Events and Their Potential Regulatory Roles in Tea Plant (Camellia sinensis L.)

RNA editing is a post-transcriptional modification process that alters the RNA sequence relative to the genomic blueprint. In plant organelles (namely, mitochondria and chloroplasts), the most common type is C-to-U, and the absence of C-to-U RNA editing results in abnormal plant development, such as etiolation and albino leaves, aborted embryonic development and retarded seedling growth. Here, through PREP, RES-Scanner, PCR and RT-PCR analyses, 38 and 139 RNA editing sites were identified from the chloroplast and mitochondrial genomes of Camellia sinensis, respectively. Analysis of the base preference around the RNA editing sites showed that in the -1 position of the edited C had more frequent occurrences of T whereas rare occurrences of G. Three conserved motifs were identified at 25 bases upstream of the RNA editing site. Structural analyses indicated that the RNA secondary structure of 32 genes, protein secondary structure of 37 genes and the three-dimensional structure of 5 proteins were altered due to RNA editing. The editing level analysis of matK and ndhD in six tea cultivars indicated that matK-701 might be involved in the color change of tea leaves. Furthermore, 218 PLS-CsPPR proteins were predicted to interact with the identified RNA editing sites. In conclusion, this study provides comprehensive insight into RNA editing events, which will facilitate further study of the RNA editing phenomenon of the tea plant.

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

A chloroplast-localized pentatricopeptide repeat protein involved in RNA editing and splicing and its effects on chloroplast development in rice

BACKGROUND

The chloroplast is the organelle responsible for photosynthesis in higher plants. The generation of functional chloroplasts depends on the precise coordination of gene expression in the nucleus and chloroplasts and is essential for the development of plants. However, little is known about nuclear-plastid regulatory mechanisms at the early stage of chloroplast generation in rice.

RESULTS

In this study, we identified a rice (Oryza sativa) mutant that exhibited albino and seedling-lethal phenotypes and named it ssa1(seedling stage albino1). Transmission electron microscopy (TEM) analysis indicated that the chloroplasts of ssa1 did not have organized thylakoid lamellae and that the chloroplast structure was destroyed. Genetic analysis revealed that the albino phenotypes of ssa1 were controlled by a pair of recessive nuclear genes. Map-based cloning experiments found that SSA1 encoded a pentapeptide repeat (PPR) protein that was allelic to OSOTP51,which was previously reported to participate in Photosystem I (PSI) assembly. The albino phenotype was reversed to the wild type (WT) phenotype when the normal SSA1 sequence was expressed in ssa1 under the drive of the actin promoter. Knockout experiments further created mutants ssa1-2/1-9, which had a phenotype similar to that of ssa1. SSA1 consisted of 7 pentatricopeptide repeat domains and two C-terminal LAGLIDADG tandem sequence motifs and was located in the chloroplast. GUS staining and qRT-PCR analysis showed that SSA1 was mainly expressed in young leaves and stems. In the ssa1 mutants, plastid genes transcribed by plastid-encoded RNA polymerase decreased, while those transcribed by nuclear-encoded RNA polymerase increased at the mRNA level. Loss-of-function SSA1 destroys RNA editing of ndhB-737 and intron splicing of atpF and ycf3-2 in the plastid genome. Yeast two-hybrid and BiFC assays revealed that SSA1 physically interacted with two new RNA editing partners, OsMORF8 and OsTRXz, which have potential functions in RNA editing and chloroplast biogenesis.

CONCLUSIONS

Rice SSA1 encodes a pentatricopeptide repeat protein, which is targeted to the chloroplast. SSA1 regulates early chloroplast development and plays a critical role in RNA editing and intron splicing in rice. These data will facilitate efforts to further elucidate the molecular mechanism of chloroplast biogenesis.

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.

Evolution of mitochondrial RNA editing in extant gymnosperms

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

An Overview of Pentatricopeptide Repeat (PPR) Proteins in the Moss Physcomitrium patens and Their Role in Organellar Gene Expression

Pentatricopeptide repeat (PPR) proteins are one type of helical repeat protein that are widespread in eukaryotes. In particular, there are several hundred PPR members in flowering plants. The majority of PPR proteins are localized in the plastids and mitochondria, where they play a crucial role in various aspects of RNA metabolism at the post-transcriptional and translational steps during gene expression. Among the early land plants, the moss Physcomitrium (formerly Physcomitrella) patens has at least 107 PPR protein-encoding genes, but most of their functions remain unclear. To elucidate the functions of PPR proteins, a reverse-genetics approach has been applied to P. patens. To date, the molecular functions of 22 PPR proteins were identified as essential factors required for either mRNA processing and stabilization, RNA splicing, or RNA editing. This review examines the P. patens PPR gene family and their current functional characterization. Similarities and a diversity of functions of PPR proteins between P. patens and flowering plants and their roles in the post-transcriptional regulation of organellar gene expression are discussed.

Fujian cytoplasmic male sterility and the fertility restorer gene OsRf19 provide a promising breeding system for hybrid rice

Although hybrid rice has been widely utilized for nearly half a century, tremendously improving rice productivity worldwide, the breeding of hybrids has been difficult because of genetic complications in male sterility and fertility-restoring systems currently available in rice. Here, we characterized Fujian Abortive cytoplasmic male sterility (CMS-FA) rice, which has shown stable male sterility controlled by the mitochondrial gene FA182; a single nuclear gene, OsRf19, completely restores fertility. This single-gene inheritance has greatly eased the breeding process. By converting CMS-WA hybrids with the CMS-FA system, we developed six hybrids that showed equivalent or better performance relative to their CMS-WA counterparts. CMS-FA/OsRf19 provides a promising system for future hybrid rice breeding.

Cytoplasmic male sterility (CMS) determined by mitochondrial genes and restorer of fertility (Rf) controlled by nuclear-encoded genes provide the breeding systems of many hybrid crops for the utilization of heterosis. Although several CMS/Rf systems have been widely exploited in rice, hybrid breeding using these systems has encountered difficulties due to either fertility instability or complications of two-locus inheritance or both. In this work, we characterized a type of CMS, Fujian Abortive cytoplasmic male sterility (CMS-FA), with stable sporophytic male sterility and a nuclear restorer gene that completely restores hybrid fertility. CMS is caused by the chimeric open reading frame FA182 that specifically occurs in the mitochondrial genome of CMS-FA rice. The restorer gene OsRf19 encodes a pentatricopeptide repeat (PPR) protein targeted to mitochondria, where it mediates the cleavage of FA182 transcripts, thus restoring male fertility. Comparative sequence analysis revealed that OsRf19 originated through a recent duplication in wild rice relatives, sharing a common ancestor with OsRf1a/OsRf5, a fertility restorer gene for Boro II and Hong-Lian CMS. We developed six restorer lines by introgressing OsRf19 into parental lines of elite CMS-WA hybrids; hybrids produced from these lines showed equivalent or better agronomic performance relative to their counterparts based on the CMS-WA system. These results demonstrate that CMS-FA/OsRf19 provides a highly promising system for future hybrid rice breeding.

Comparative chloroplast genome and transcriptome analysis on the ancient genus Isoetes from China

Isoetes is a famous living fossil that plays a significant role in the evolutionary studies of the plant kingdom. To explore the adaptive evolution of the ancient genus Isoetes from China, we focused on Isoetes yunguiensis (Q.F. Wang and W.C. Taylor), I. shangrilaensis (X. Li, Y.Q. Huang, X.K. Dai & X. Liu), I. taiwanensis (DeVol), I. sinensis (T.C. Palmer), I. hypsophila_GHC (Handel-Mazzetti), and I. hypsophila_HZS in this study. We sequenced, assembled, and annotated six individuals' chloroplast genomes and transcriptomes, and performed a series of analyses to investigate their chloroplast genome structures, RNA editing events, and adaptive evolution. The six chloroplast genomes of Isoetes exhibited a typical quadripartite structure with conserved genome sequence and structure. Comparative analyses of Isoetes species demonstrated that the gene organization, genome size, and GC contents of the chloroplast genome are highly conserved across the genus. Besides, our positive selection analyses suggested that one positively selected gene was statistically supported in Isoetes chloroplast genomes using the likelihood ratio test (LRT) based on branch-site models. Moreover, we detected positive selection signals using transcriptome data, suggesting that nuclear-encoded genes involved in the adaption of Isoetes species to the extreme environment of the Qinghai-Tibetan Plateau (QTP). In addition, we identified 291-579 RNA editing sites in the chloroplast genomes of six Isoetes based on transcriptome data, well above the average of angiosperms. RNA editing in protein-coding transcripts results from amino acid changes to increase their hydrophobicity and conservation in Isoetes, which may help proteins form functional three-dimensional structure. Overall, the results of this study provide comprehensive transcriptome and chloroplast genome resources and contribute to a better understanding of adaptive evolutionary and molecular biology in Isoetes.