REDO: RNA Editing Detection in Plant Organelles Based on Variant Calling Results

Complete mitochondrial genome of Agropyron cristatum reveals gene transfer and RNA editing events

BACKGROUND

As an important forage in arid and semi-arid regions, Agropyron cristatum provides livestock with exceptionally high nutritional value. Additionally, A. cristatum exhibits outstanding genetic characteristics to endure drought and disease. Therefore, rich genetic diversity serves as a cornerstone for the improvement of major food crops. The purposes of this study were to systematically describe mitogenome of A.cristatum and preliminarily analyze its internal variations.

RESULT

The A. cristatum mitogenome was a single-ring molecular structure of 381,065 bp that comprised 52 genes, including 35 protein-coding, 3 rRNA and 14 tRNA genes. Among these, two pseudoprotein-coding genes and multiple copies of tRNA genes were observed. A total of 320 repetitive sequences was found to cover more than 10% of the mitogenome (105 simple sequences, 185 dispersed and 30 tandem repeats), which led to a large number of fragment rearrangements in the mitogenome of A. cristatum. Leucine was the most frequent amino acid (n = 1087,10.8%) in the protein-coding genes of A. cristatum mitogenome, and the highest usage codon was ATG (initiation codon). The number of A/T changes at the third base of the codon was much higher than that of G/C. Among 23 PCGs, the range of Pi values is from 0.0021 to 0.0539, with an average of 0.013. Additionally, 81 RNA editing sites were predicted, which were considerably fewer than those reported in other plant mitogenomes. Most of the RNA editing site base positions were concentrated at the first and second codon bases, which were C to T transitions. Moreover, we identified 95 sequence fragments (total length of 34, 343 bp) that were transferred from the chloroplast to mitochondria genes, introns, and intergenic regions. The stability of the tRNA genes was maintained during this process. Selection pressure analysis of 23 protein-coding genes shared by 15 Poaceae plants, showed that most genes were subjected to purifying selection during evolution, whereas rps4, cob, mttB, and ccmB underwent positive selection in different plants. Finally, a phylogenetic tree was constructed based on 22 plant mitogenomes, which showed that Agropyron plants have a high degree of independent heritability in Triticeae.

CONCLUSION

The findings of this study provide new data for a better understanding of A. cristatum genes, and demonstrate that mitogenomes are suitable for the study of plant classifications, such as those of Agropyron. Moreover, it provides a reference for further exploration of the phylogenetic relationships within Agropyron species, and establishes a theoretical basis for the subsequent development and utilization of A. cristatum plant germplasm resources.

Complete mitochondrial genome assembly of Zizania latifolia and comparative genome analysis

Zizania latifolia (Griseb.) Turcz. ex Stapf has been cultivated as a popular aquatic vegetable in China due to its important nutritional, medicinal, ecological, and economic values. The complete mitochondrial genome (mitogenome) of Z. latifolia has not been previously studied and reported, which has hindered its molecular systematics and understanding of evolutionary processes. Here, we assembled the complete mitogenome of Z. latifolia and performed a comprehensive analysis including genome organization, repetitive sequences, RNA editing event, intercellular gene transfer, phylogenetic analysis, and comparative mitogenome analysis. The mitogenome of Z. latifolia was estimated to have a circular molecule of 392,219 bp and 58 genes consisting of three rRNA genes, 20 tRNA genes, and 35 protein-coding genes (PCGs). There were 46 and 20 simple sequence repeats (SSRs) with different motifs identified from the mitogenome and chloroplast genome of Z. latifolia, respectively. Furthermore, 49 homologous fragments were observed to transfer from the chloroplast genome to the mitogenome of Z. latifolia, accounting for 47,500 bp, presenting 12.1% of the whole mitogenome. In addition, there were 11 gene-containing homologous regions between the mitogenome and chloroplast genome of Z. latifolia. Also, approximately 85% of fragments from the mitogenome were duplicated in the Z. latifolia nuclear genome. Selection pressure analysis revealed that most of the mitochondrial genes were highly conserved except for ccmFc, ccmFn, matR, rps1, and rps3. A total of 93 RNA editing sites were found in the PCGs of the mitogenome. Z. latifolia and Oryza minuta are the most closely related, as shown by collinear analysis and the phylogenetic analysis. We found that repeat sequences and foreign sequences in the mitogenomes of Oryzoideae plants were associated with genome rearrangements. In general, the availability of the Z. latifolia mitogenome will contribute valuable information to our understanding of the molecular and genomic aspects of Zizania.

Comparative Chloroplast Genomes Analysis Provided Adaptive Evolution Insights in Medicago ruthenica

A perennial leguminous forage, Medicago ruthenica has outstanding tolerance to abiotic stresses. The genome of Medicago ruthenica is large and has a complex genetic background, making it challenging to accurately determine genetic information. However, the chloroplast genome is widely used for researching issues related to evolution, genetic diversity, and other studies. To better understand its chloroplast characteristics and adaptive evolution, chloroplast genomes of 61 Medicago ruthenica were assembled (including 16 cultivated Medicago ruthenica germplasm and 45 wild Medicago ruthenica germplasm). These were used to construct the pan-chloroplast genome of Medicago ruthenica, and the chloroplast genomes of cultivated and wild Medicago ruthenica were compared and analyzed. Phylogenetic and haplotype analyses revealed two main clades of 61 Medicago ruthenica germplasm chloroplast genomes, distributed in eastern and western regions. Meanwhile, based on chloroplast variation information, 61 Medicago ruthenica germplasm can be divided into three genetic groups. Unlike the phylogenetic tree constructed from the chloroplast genome, a new intermediate group has been identified, mainly consisting of samples from the eastern region of Inner Mongolia, Shanxi Province, and Hebei Province. Transcriptomic analysis showed that 29 genes were upregulated and three genes were downregulated. The analysis of these genes mainly focuses on enhancing plant resilience and adapting adversity by stabilizing the photosystem structure and promoting protein synthesis. Additionally, in the analysis of adaptive evolution, the accD, clpP and ycf1 genes showed higher average Ka/Ks ratios and exhibited significant nucleotide diversity, indicating that these genes are strongly positively selected. The editing efficiency of the ycf1 and clpP genes significantly increases under abiotic stress, which may positively contribute to plant adaptation to the environment. In conclusion, the construction and comparative analysis of the complete chloroplast genomes of 61 Medicago ruthenica germplasm from different regions not only revealed new insights into the genetic variation and phylogenetic relationships of Medicago ruthenica germplasm, but also highlighted the importance of chloroplast transcriptome analysis in elucidating the model of chloroplast responses to abiotic stress. These provide valuable information for further research on the adaptive evolution of Medicago ruthenica.

Transcriptome analysis of Citrus sinensis reveals potential responsive events triggered by Candidatus Liberibacter asiaticus

Citrus Huanglongbing (HLB), caused by Candidatus Liberibacter asiaticus (CLas), is a devastating immune-mediated disorder that has a detrimental effect on the citrus industry, with the distinguishing feature being an eruption of reactive oxygen species (ROS). This study explored the alterations in antioxidant enzyme activity, transcriptome, and RNA editing events of organelles in C. sinensis during CLas infection. Results indicated that there were fluctuations in the performance of antioxidant enzymes, such as ascorbate peroxidase (APX), catalase (CAT), glutathione reductase (GR), peroxidase (POD), and superoxide dismutase (SOD), in plants affected by HLB. Transcriptome analysis revealed 3604 genes with altered expression patterns between CLas-infected and healthy samples, including those associated with photosynthesis, biotic interactions, and phytohormones. Samples infected with CLas showed a decrease in the expression of most genes associated with photosynthesis and gibberellin metabolism. It was discovered that RNA editing frequency and the expression level of various genes in the chloroplast and mitochondrion genomes were affected by CLas infection. Our findings provide insights into the inhibition of photosynthesis, gibberellin metabolism, and antioxidant enzymes during CLas infection in C. sinensis.

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

Complete mitochondrial genome of the endangered Prunus pedunculata (Prunoideae, Rosaceae) in China: characterization and phylogenetic analysis

Introduction

Prunus pedunculata (Prunoideae: Rosaceae), a relic shrub with strong resistance and multiple application values, is endangered in China. Extensive research had been devoted to gene expression, molecular markers, plastid genome analysis, and genetic background investigations of P. pedunculata. However, the mitochondrial genome of this species has not been systematically described, owing to the complexity of the plant mitogenome.

Methods

In the present research, the complete mitochondrial genome of P. pedunculata was assembled, annotated, and characterized. The genomic features, gene content and repetitive sequences were analyzed. The genomic variation and phylogenetic analysis have been extensively enumerated.

Results and discussion

The P. pedunculata mitogenome is a circular molecule with a total length of 405,855 bp and a GC content of 45.63%, which are the smallest size and highest GC content among the known Prunus mitochondrial genomes. The mitogenome of P. pedunculata encodes 62 genes, including 34 unique protein-coding genes (PCGs, excluding three possible pseudogenes), three ribosomal RNA genes, and 19 transfer RNA genes. The mitogenome is rich in repetitive sequences, counting 112 simple sequence repeats, 15 tandem repeats, and 50 interspersed repetitive sequences, with a total repeat length of 11,793 bp, accounting for 2.91% of the complete genome. Leucine (Leu) was a predominant amino acid in PCGs, with a frequency of 10.67%, whereas cysteine (Cys) and tryptophan (Trp) were the least adopted. The most frequently used codon was UUU (Phe), with a relative synonymous codon usage (RSCU) value of 1.12. Selective pressure was calculated based on 20 shared PCGs in the mitogenomes of the 32 species, most of which were subjected to purifying selection (Ka/Ks < 1), whereas ccmC and ccmFn underwent positive selection. A total of 262 potential RNA editing sites in 26 PCGs were identified. Furthermore, 56 chloroplast-derived fragments were ascertained in the mitogenome, ranging from 30 to 858 bp, and were mainly located across IGS (intergenic spacer) regions or rRNA genes. These findings verify the occurrence of intracellular gene transfer events from the chloroplast to the mitochondria. Furthermore, the phylogenetic relationship of P. pedunculata was supported by the mitogenome data of 30 other taxa of the Rosaceae family. Understanding the mitochondrial genome characteristics of P. pedunculata is of great importance to promote comprehension of its genetic background and this study provides a basis for the genetic breeding of Prunus.

A Genome-Wide Analysis of the Pentatricopeptide Repeat Protein Gene Family in Two Kiwifruit Species with an Emphasis on the Role of RNA Editing in Pathogen Stress

Kiwifruit is a perennial fruit tree with high nutritional and economic value; however, various pathogen stresses have resulted in reductions in its yield and quality. Pentatricopeptide repeat proteins (PPRs), characterized by tandem repetitions of 35 amino acid motifs, play roles in RNA editing, mRNA stability, and splicing. They may also regulate plant development and growth. Nevertheless, the roles of PPRs in plant development and disease resistance remain unclear. In this study, we focused on the roles of PPRs in the fruit development and pathogen stress of kiwifruit and conducted a series of analyses of the PPR gene family in two representative kiwifruit species (Actinidia chinensis (Ach) and Actinidia eriantha (Ace)) with markedly different degrees of disease resistance. A total of 497 and 499 PPRs were identified in Ach and Ace, respectively. All the kiwifruit PPRs could be phylogenetically divided into four subfamilies. There were about 40.68% PPRs predicted to be localized to mitochondria or chloroplasts. A synteny analysis showed that the expansion of the kiwifruit PPRs mainly originated from segmental duplication. Based on RNA-seq data from the fruit over 12 periods of development and maturity, a weighted correlation network analysis suggested that two PPRs, Actinidia20495.t1 and Actinidia15159.t1, may be involved in fruit development and maturation. In addition, we observed different responses with respect to the expression of PPRs and RNA editing between resistant and susceptible kiwifruits following infection with pathogenic bacteria, indicating the regulatory role of PPRs in the stress response via the modulation of RNA editing. The differentially expressed upstream transcription factors of the PPRs were further identified; they may regulate resistance adaption by modulating the expression of the PPRs. Collectively, these results suggest that PPRs play roles in the development and disease resistance of kiwifruit and provide candidate genes for further clarifying the resistance mechanisms in kiwifruits.

Comparative analysis of the chloroplast genomes of Rosa species and RNA editing analysis

BACKGROUND

The genus Rosa (Rosaceae) contains approximately 200 species, most of which have high ecological and economic values. Chloroplast genome sequences are important for studying species differentiation, phylogeny, and RNA editing.

RESULTS

In this study, the chloroplast genomes of three Rosa species, Rosa hybrida, Rosa acicularis, and Rosa rubiginosa, were assembled and compared with other reported Rosa chloroplast genomes. To investigate the RNA editing sites in R. hybrida (commercial rose cultivar), we mapped RNA-sequencing data to the chloroplast genome and analyzed their post-transcriptional features. Rosa chloroplast genomes presented a quadripartite structure and had highly conserved gene order and gene content. We identified four mutation hotspots (ycf3-trnS, trnT-trnL, psbE-petL, and ycf1) as candidate molecular markers for differentiation in the Rosa species. Additionally, 22 chloroplast genomic fragments with a total length of 6,192 bp and > 90% sequence similarity with their counterparts were identified in the mitochondrial genome, representing 3.96% of the chloroplast genome. Phylogenetic analysis including all sections and all subgenera revealed that the earliest divergence in the chloroplast phylogeny roughly distinguished species of sections Pimpinellifoliae and Rosa and subgenera Hulthemia. Moreover, DNA- and RNA-sequencing data revealed 19 RNA editing sites, including three synonymous and 16 nonsynonymous, in the chloroplast genome of R. hybrida that were distributed among 13 genes.

CONCLUSIONS

The genome structure and gene content of Rosa chloroplast genomes are similar across various species. Phylogenetic analysis based on the Rosa chloroplast genomes has high resolution. Additionally, a total of 19 RNA editing sites were validated by RNA-Seq mapping in R. hybrida. The results provide valuable information for RNA editing and evolutionary studies of Rosa and a basis for further studies on genomic breeding of Rosa species.

Mitogenomic Research of Silverleaf Sunflower (Helianthus argophyllus) and Its Interspecific Hybrids

Interspecific hybridization is widespread for sunflowers, both in wild populations and commercial breeding. One of the most common species that can efficiently cross with Helianthus annuus is the Silverleaf sunflower-Helianthus argophyllus. The current study carried out structural and functional organization analyses of mitochondrial DNA in H. argophyllus and the interspecific hybrid, H. annuus (VIR114A line) × H. argophyllus. The complete mitogenome of H. argophyllus counts 300,843 bp, has a similar organization to the mitogenome of cultivated sunflower, and holds SNPs typical for wild sunflowers. RNA editing analysis predicted 484 sites in H. argophyllus mitochondrial CDS. The mitochondrial genome of the H. annuus × H. argophyllus hybrid is identical to the maternal line (VIR114A). We expected that significant rearrangements in the mitochondrial DNA of the hybrid would occur, due to the frequent recombination. However, the hybrid mitogenome lacks rearrangements, presumably due to the preservation of nuclear-cytoplasmic interaction paths.

Chloroplast Pan-Genomes and Comparative Transcriptomics Reveal Genetic Variation and Temperature Adaptation in the Cucumber

Although whole genome sequencing, genetic variation mapping, and pan-genome studies have been done on a large group of cucumber nuclear genomes, organelle genome information is largely unclear. As an important component of the organelle genome, the chloroplast genome is highly conserved, which makes it a useful tool for studying plant phylogeny, crop domestication, and species adaptation. Here, we have constructed the first cucumber chloroplast pan-genome based on 121 cucumber germplasms, and investigated the genetic variations of the cucumber chloroplast genome through comparative genomic, phylogenetic, haplotype, and population genetic structure analysis. Meanwhile, we explored the changes in expression of cucumber chloroplast genes under high- and low-temperature stimulation via transcriptome analysis. As a result, a total of 50 complete chloroplast genomes were successfully assembled from 121 cucumber resequencing data, ranging in size from 156,616-157,641 bp. The 50 cucumber chloroplast genomes have typical quadripartite structures, consisting of a large single copy (LSC, 86,339-86,883 bp), a small single copy (SSC, 18,069-18,363 bp), and two inverted repeats (IRs, 25,166-25,797 bp). Comparative genomic, haplotype, and population genetic structure results showed that there is more genetic variation in Indian ecotype cucumbers compared to other cucumber cultivars, which means that many genetic resources remain to be explored in Indian ecotype cucumbers. Phylogenetic analysis showed that the 50 cucumber germplasms could be classified into 3 types: East Asian, Eurasian + Indian, and Xishuangbanna + Indian. The transcriptomic analysis showed that matK were significantly up-regulated under high- and low-temperature stresses, further demonstrating that cucumber chloroplasts respond to temperature adversity by regulating lipid metabolism and ribosome metabolism. Further, accD has higher editing efficiency under high-temperature stress, which may contribute to the heat tolerance. These studies provide useful insight into genetic variation in the chloroplast genome, and established the foundation for exploring the mechanisms of temperature-stimulated chloroplast adaptation.

A Bibliometric Study for Plant RNA Editing Research: Trends and Future Challenges

RNA editing is a post-transcriptional process that introduces changes in RNA sequences encoded by nuclear, mitochondrial, or plastid genomes. To understand the research progress of plant RNA editing, we comprehensively analyze the articles on plant RNA editing from 2001 to 2022 through bibliometric methods. Nucleic Acids Research, Plant Journal and Plant cell are the journals that deserve attention with their high production, total local citation scores (TLCS), and h-indexes. The USA, China, and Germany are the top three countries with highly productive publications. Ulm University, Cornell University, and Chinese Acad Sci are excellent cooperative institutions with a high level of influence in the field, and KNOOP V and TAKENAKA M are good partnership. Plant RNA editing researches concentrate on the subject categories of Biochemistry & Molecular Biology, Plant Sciences, Genetics & Heredity, etc. Plant mitochondria, genome editing and messenger-RNA may be the research hotspots in the future. The main plant RNA editing research tools are JACUSA, SPRINT, and REDO, and the main databases are REDIdb, PED, and dbRES. At present, the research streams are (1) RNA editing sites; (2) Pentapeptide repeat protein (PPR) involved in RNA editing; (3) RNA editing factors. Overall, this article summarizes the research overview of plant RNA editing until 2022 and provides theoretical implications for its possible future directions.

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.

Genome-Wide Analysis of Multiple Organellar RNA Editing Factor (MORF) Family in Kiwifruit (Actinidia chinensis) Reveals Its Roles in Chloroplast RNA Editing and Pathogens Stress

Kiwifruit (Actinidia chinensis) is well known for its high vitamin C content and good taste. Various diseases, especially bacterial canker, are a serious threat to the yield of kiwifruit. Multiple organellar RNA editing factor (MORF) genes are pivotal factors in the RNA editosome that mediates Cytosine-to-Uracil RNA editing, and they are also indispensable for the regulation of chloroplast development, plant growth, and response to stresses. Although the kiwifruit genome has been released, little is known about MORF genes in kiwifruit at the genome-wide level, especially those involved in the response to pathogens stress. In this study, we identified ten MORF genes in the kiwifruit genome. The genomic structures and chromosomal locations analysis indicated that all the MORF genes consisted of three conserved motifs, and they were distributed widely across the seven linkage groups and one contig of the kiwifruit genome. Based on the structural features of MORF proteins and the topology of the phylogenetic tree, the kiwifruit MORF gene family members were classified into six groups (Groups A-F). A synteny analysis indicated that two pairs of MORF genes were tandemly duplicated and five pairs of MORF genes were segmentally duplicated. Moreover, based on analysis of RNA-seq data from five tissues of kiwifruit, we found that both expressions of MORF genes and chloroplast RNA editing exhibited tissue-specific patterns. MORF2 and MORF9 were highly expressed in leaf and shoot, and may be responsible for chloroplast RNA editing, especially the ndhB genes. We also observed different MORF expression and chloroplast RNA editing profiles between resistant and susceptible kiwifruits after pathogen infection, indicating the roles of MORF genes in stress response by modulating the editing extend of mRNA. These results provide a solid foundation for further analyses of the functions and molecular evolution of MORF genes, in particular, for clarifying the resistance mechanisms in kiwifruits and breeding new cultivars with high resistance.

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.

Assembly and Analysis of the Complete Mitochondrial Genome of Capsella bursa-pastoris

Shepherd's purse (Capsella bursa-pastoris) is a cosmopolitan annual weed and a promising model plant for studying allopolyploidization in the evolution of angiosperms. Though plant mitochondrial genomes are a valuable source of genetic information, they are hard to assemble. At present, only the complete mitogenome of C. rubella is available out of all species of the genus Capsella. In this work, we have assembled the complete mitogenome of C. bursa-pastoris using high-precision PacBio SMRT third-generation sequencing technology. It is 287,799 bp long and contains 32 protein-coding genes, 3 rRNAs, 25 tRNAs corresponding to 15 amino acids, and 8 open reading frames (ORFs) supported by RNAseq data. Though many repeat regions have been found, none of them is longer than 1 kbp, and the most frequent structural variant originated from these repeats is present in only 4% of the mitogenome copies. The mitochondrial DNA sequence of C. bursa-pastoris differs from C. rubella, but not from C. orientalis, by two long inversions, suggesting that C. orientalis could be its maternal progenitor species. In total, 377 C to U RNA editing sites have been detected. All genes except cox1 and atp8 contain RNA editing sites, and most of them lead to non-synonymous changes of amino acids. Most of the identified RNA editing sites are identical to corresponding RNA editing sites in A. thaliana.

The Amount of RNA Editing Sites in Liverwort Organellar Genes Is Correlated with GC Content and Nuclear PPR Protein Diversity

RNA editing occurs in the organellar mRNAs of all land plants but the marchantioid liverworts, making liverworts a perfect group for studying the evolution of RNA editing. Here, we profiled the RNA editing of 42 exemplars spanning the ordinal phylogenetic diversity of liverworts, and screened for the nuclear-encoded pentatricopeptide repeat (PPR) proteins in the transcriptome assemblies of these taxa. We identified 7,428 RNA editing sites in 128 organellar genes from 31 non-marchantioid liverwort species, and characterized 25,059 PPR protein sequences. The abundance of organellar RNA editing sites varies greatly among liverwort lineages, genes, and codon positions, and shows strong positive correlations with the GC content of protein-coding genes, and the diversity of the PLS class of nuclear PPR proteins.

Dynamic response of RNA editing to temperature in grape by RNA deep sequencing

RNA editing is a post-transcriptional process of modifying genetic information on RNA molecules, which provides cells an additional level of gene expression regulation. Unlike mammals, in land plants, RNA editing converts C-to-U residues in organelles. However, its potential roles in response to different stressors (heat, salt, and so on) remains unclear. Grape is one of the most popular and economically important fruits in the world, and its production, like other crops, must deal with abiotic and biotic stresses, which cause reductions in yield and fruit quality. In our study, we tested the influence of the environmental factor temperature on RNA editing process in the whole mRNA from grape organelle. In total, we identified 122 and 627 RNA editing sites in chloroplast and mitochondria respectively with the average editing efficiency nearly ~ 60%. The analyses revealed that number of non-synonymous editing sites were higher than that of synonymous editing sites, and the amino acid substitution type tends to be hydrophobic. Additionally, the overall editing level decreased with the temperature rises, especially for several gene transcripts in chloroplast and mitochondria (matK, ndhB, etc.). We also found that the expression level of most PPR genes decreased with the temperature rises, which may contribute to the decline of RNA editing efficiency at high temperature. Our findings suggested that the RNA editing events were very sensitive to heat stress; the changes of amino acid in RNA editing genes may contribute to the stress adaption for grape.

RNA editing in plants: A comprehensive survey of bioinformatics tools and databases

RNA editing is a widespread epitranscriptomic mechanism by which primary RNAs are specifically modified through insertions/deletions or nucleotide substitutions. In plants, RNA editing occurs in organelles (plastids and mitochondria), involves the cytosine to uridine modification (rarely uridine to cytosine) within protein-coding and non-protein-coding regions of RNAs and affects organelle biogenesis, adaptation to environmental changes and signal transduction. High-throughput sequencing technologies have dramatically improved the detection of RNA editing sites at genomic scale. Consequently, different bioinformatics resources have been released to discovery and/or collect novel events. Here, we review and describe the state-of-the-art bioinformatics tools devoted to the characterization of RNA editing in plant organelles with the aim to improve our knowledge about this fascinating but yet under investigated process.

High Level of Conservation of Mitochondrial RNA Editing Sites Among Four Populus Species

RNA editing occurs in the endosymbiont organelles of higher plants as C-to-U conversions of defined nucleotides. The availability of large quantities of RNA sequencing data makes it possible to identify RNA editing sites and to quantify their editing extent. We have investigated RNA editing in 34 protein-coding mitochondrial transcripts of four Populus species, a genus noteworthy for its remarkably small number of RNA editing sites compared to other angiosperms. 27 of these transcripts were subject to RNA editing in at least one species. In total, 355 RNA editing sites were identified with high confidence, their editing extents ranging from 10 to 100%. The most heavily edited transcripts were ccmB with the highest density of RNA editing sites (53.7 sites / kb) and ccmFn with the highest number of sites (39 sites). Most of the editing events are at position 1 or 2 of the codons, usually altering the encoded amino acid, and are highly conserved among the species, also with regard to their editing extent. However, one SNP was found in the newly sequenced and annotated mitochondrial genome of P. alba resulting in the loss of an RNA editing site compared to P. tremula and P. davidiana This SNP causes a C-to-T transition and an amino acid exchange from Ser to Phe, highlighting the widely discussed role of RNA editing in compensating mutations.