CURE-Chloroplast: a chloroplast C-to-U RNA editing predictor for seed plants

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

A-to-I RNA editing shows dramatic up-regulation in osteosarcoma and broadly regulates tumor-related genes by altering microRNA target regions

A-to-I RNA editing is a prevalent type of RNA modification in animals. The dysregulation of RNA editing has led to multiple human cancers. However, the role of RNA editing has never been studied in osteosarcoma, a complex bone cancer with unknown molecular basis. We retrieved the RNA-sequencing data from 24 primary osteosarcoma patients and 3 healthy controls. We systematically profiled the RNA editomes in these samples and quantitatively identified reliable differential editing sites (DES) between osteosarcoma and normal samples. RNA editing efficiency is dramatically increased in osteosarcoma, presumably due to the significant up-regulation of editing enzymes ADAR1 and ADAR2. Up-regulated DES in osteosarcoma are enriched in 3'UTRs. Strikingly, such 3'UTR sites are further enriched in microRNA binding regions of gene EMP2 and other oncogenes, abolishing the microRNA suppression on target genes. Accordingly, the expression of these tumor-promoting genes is elevated in osteosarcoma. There might be an RNA editing-dependent pathway leading to osteosarcoma. We expanded our knowledge on the potential roles of RNA editing in oncogenesis. Based on these molecular features, our work is valuable for future prognosis and diagnosis of osteosarcoma.

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.

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.

Comprehensive analysis of chloroplast genome of Albizia julibrissin Durazz. (Leguminosae sp.)

The Albizia julibrissin chloroplasts have a classical chloroplast genome structure, containing 93 coding genes and 34 non-coding genes. Our research provides basic data for plant phylogenetic evolutionary studies. There is limited genomic information available for the important Chinese herb Albizia julibrissin Durazz. In this study, we constructed the chloroplast (Cp) genome of A. julibrissin. The length of the assembled Cp genome was 175,922 bp consisting of four conserved regions: a 5145 bp small single-copy (SSC) region, a 91,323 bp large single-copy (LSC) region, and two identical length-inverted repeat (IR) regions (39,725 bp). This Cp genome included 34 non-coding RNAs and 93 unique genes, the former contains 30 transfer and 4 ribosomal RNA genes. Gene annotation indicated some of the coding genes (82) in the A. julibrissin Cp genome classified in the Leguminosae family, with some to other related families (11). The results show that low GC content (36.9%) and codon bias towards A- or T-terminal codons may affect the frequency of gene codon usage. The sequence analysis identified 30 forward, 18 palindrome, and 1 reverse repeat > 30 bp length, and 149 simple sequence repeats (SSR). Fifty-five RNA editing sites in the Cp of A. julibrissin were predicted, most of which are C-to-U conversions. Analysis of the reverse repeat expansion or contraction and divergence area between several species, including A. julibrissin, was performed. The phylogenetic tree revealed that A. julibrissin was most closely related to Albizia odoratissima and Albizia bracteata, followed by Samanea saman, forming an evolutionary branch with Mimosa pudica and Leucaena trichandra. The research results are helpful for breeding and genetic improvement of A. julibrissin, and also provide valuable information for understanding the evolution of this plant.

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis, identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa 'Evolvulus' could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01051-w.

Detection and Analysis of C-to-U RNA Editing in Rice Mitochondria-Encoded ORFs

Cytidine to uridine (C-to-U) RNA editing is an important type of substitutional RNA modification and is almost omnipresent in plant chloroplasts and mitochondria. In rice mitochondria, 491 C-to-U editing sites have been identified previously, and case studies have elucidated the function of several C-to-U editing sites in rice, but the functional consequence of most C-to-U alterations needs to be investigated further. Here, by means of Sanger sequencing and publicly available RNA-seq data, we identified a total of 569 C-to-U editing sites in rice mitochondria-encoded open reading frames (ORFs), 85.41% of these editing sites were observed on the first or the second base of a codon, resulting in the alteration of encoded amino acid. Moreover, we found some novel editing sites and several inaccurately annotated sites which may be functionally important, based on the highly conserved amino acids encoded by these edited codons. Finally, we annotated all 569 C-to-U RNA editing sites in their biological context. More precise information about C-to-U editing sites in rice mitochondria-encoded ORFs will facilitate our investigation on the function of C-to-U editing events in rice and also provide a valid benchmark from rice for the analysis of mitochondria C-to-U editing in other plant species.

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

Background

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

Results

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

Conclusions

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

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

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

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

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

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.

The Complete Chloroplast Genome Sequences of the Medicinal Plant Forsythia suspensa (Oleaceae)

Forsythia suspensa is an important medicinal plant and traditionally applied for the treatment of inflammation, pyrexia, gonorrhea, diabetes, and so on. However, there is limited sequence and genomic information available for F. suspensa. Here, we produced the complete chloroplast genomes of F. suspensa using Illumina sequencing technology. F. suspensa is the first sequenced member within the genus Forsythia (Oleaceae). The gene order and organization of the chloroplast genome of F. suspensa are similar to other Oleaceae chloroplast genomes. The F. suspensa chloroplast genome is 156,404 bp in length, exhibits a conserved quadripartite structure with a large single-copy (LSC; 87,159 bp) region, and a small single-copy (SSC; 17,811 bp) region interspersed between inverted repeat (IRa/b; 25,717 bp) regions. A total of 114 unique genes were annotated, including 80 protein-coding genes, 30 tRNA, and four rRNA. The low GC content (37.8%) and codon usage bias for A- or T-ending codons may largely affect gene codon usage. Sequence analysis identified a total of 26 forward repeats, 23 palindrome repeats with lengths >30 bp (identity > 90%), and 54 simple sequence repeats (SSRs) with an average rate of 0.35 SSRs/kb. We predicted 52 RNA editing sites in the chloroplast of F. suspensa, all for C-to-U transitions. IR expansion or contraction and the divergent regions were analyzed among several species including the reported F. suspensa in this study. Phylogenetic analysis based on whole-plastome revealed that F. suspensa, as a member of the Oleaceae family, diverged relatively early from Lamiales. This study will contribute to strengthening medicinal resource conservation, molecular phylogenetic, and genetic engineering research investigations of this species.

RNA editing and drug discovery for cancer therapy

RNA editing is vital to provide the RNA and protein complexity to regulate the gene expression. Correct RNA editing maintains the cell function and organism development. Imbalance of the RNA editing machinery may lead to diseases and cancers. Recently, RNA editing has been recognized as a target for drug discovery although few studies targeting RNA editing for disease and cancer therapy were reported in the field of natural products. Therefore, RNA editing may be a potential target for therapeutic natural products. In this review, we provide a literature overview of the biological functions of RNA editing on gene expression, diseases, cancers, and drugs. The bioinformatics resources of RNA editing were also summarized.

PREPACT 2.0: Predicting C-to-U and U-to-C RNA Editing in Organelle Genome Sequences with Multiple References and Curated RNA Editing Annotation

RNA editing is vast in some genetic systems, with up to thousands of targeted C-to-U and U-to-C substitutions in mitochondria and chloroplasts of certain plants. Efficient prognoses of RNA editing in organelle genomes will help to reveal overlooked cases of editing. We present PREPACT 2.0 (http://www.prepact.de) with numerous enhancements of our previously developed Plant RNA Editing Prediction & Analysis Computer Tool. Reference organelle transcriptomes for editing prediction have been extended and reorganized to include 19 curated mitochondrial and 13 chloroplast genomes, now allowing to distinguish RNA editing sites from “pre-edited” sites. Queries may be run against multiple references and a new “commons” function identifies and highlights orthologous candidate editing sites congruently predicted by multiple references. Enhancements to the BLASTX mode in PREPACT 2.0 allow querying of complete novel organelle genomes within a few minutes, identifying protein genes and candidate RNA editing sites simultaneously without prior user analyses.

Capturing the biofuel wellhead and powerhouse: the chloroplast and mitochondrial genomes of the leguminous feedstock tree Pongamia pinnata

Pongamia pinnata (syn. Millettia pinnata) is a novel, fast-growing arboreal legume that bears prolific quantities of oil-rich seeds suitable for the production of biodiesel and aviation biofuel. Here, we have used Illumina® 'Second Generation DNA Sequencing (2GS)' and a new short-read de novo assembler, SaSSY, to assemble and annotate the Pongamia chloroplast (152,968 bp; cpDNA) and mitochondrial (425,718 bp; mtDNA) genomes. We also show that SaSSY can be used to accurately assemble 2GS data, by re-assembling the Lotus japonicus cpDNA and in the process assemble its mtDNA (380,861 bp). The Pongamia cpDNA contains 77 unique protein-coding genes and is almost 60% gene-dense. It contains a 50 kb inversion common to other legumes, as well as a novel 6.5 kb inversion that is responsible for the non-disruptive, re-orientation of five protein-coding genes. Additionally, two copies of an inverted repeat firmly place the species outside the subclade of the Fabaceae lacking the inverted repeat. The Pongamia and L. japonicus mtDNA contain just 33 and 31 unique protein-coding genes, respectively, and like other angiosperm mtDNA, have expanded intergenic and multiple repeat regions. Through comparative analysis with Vigna radiata we measured the average synonymous and non-synonymous divergence of all three legume mitochondrial (1.59% and 2.40%, respectively) and chloroplast (8.37% and 8.99%, respectively) protein-coding genes. Finally, we explored the relatedness of Pongamia within the Fabaceae and showed the utility of the organellar genome sequences by mapping transcriptomic data to identify up- and down-regulated stress-responsive gene candidates and confirm in silico predicted RNA editing sites.

The comparative chloroplast genomic analysis of photosynthetic orchids and developing DNA markers to distinguish Phalaenopsis orchids

The chloroplast genome of Phalaenopsis equestris was determined and compared to those of Phalaenopsis aphrodite and Oncidium Gower Ramsey in Orchidaceae. The chloroplast genome of P. equestris is 148,959 bp, and a pair of inverted repeats (25,846 bp) separates the genome into large single-copy (85,967 bp) and small single-copy (11,300 bp) regions. The genome encodes 109 genes, including 4 rRNA, 30 tRNA and 75 protein-coding genes, but loses four ndh genes (ndhA, E, F and H) and seven other ndh genes are pseudogenes. The rate of inter-species variation between the two moth orchids was 0.74% (1107 sites) for single nucleotide substitution and 0.24% for insertions (161 sites; 1388 bp) and deletions (189 sites; 1393 bp). The IR regions have a lower rate of nucleotide substitution (3.5-5.8-fold) and indels (4.3-7.1-fold) than single-copy regions. The intergenic spacers are the most divergent, and based on the length variation of the three intergenic spacers, 11 native Phalaenopsis orchids could be successfully distinguished. The coding genes, IR junction and RNA editing sites are relatively more conserved between the two moth orchids than between those of Phalaenopsis and Oncidium spp.

Alignment free comparison: similarity distribution between the DNA primary sequences based on the shortest absent word

This work proposes an alignment free comparison model for the DNA primary sequences. In this paper, we treat the double strands of the DNA rather than single strand. We define the shortest absent word of the double strands between the DNA sequences and some properties are studied to speed up the algorithm for searching the shortest absent word. We present a novel model for comparison, in which the similarity distribution is introduced to describe the similarity between the sequences. A distance measure is deduced based on the Shannon entropy meanwhile is used in phylogenetic analysis. Some experiments show that our model performs well in the field of sequence analysis.

Chloroplast genome variation in upland and lowland switchgrass

Switchgrass (Panicum virgatum L.) exists at multiple ploidies and two phenotypically distinct ecotypes. To facilitate interploidal comparisons and to understand the extent of sequence variation within existing breeding pools, two complete switchgrass chloroplast genomes were sequenced from individuals representative of the upland and lowland ecotypes. The results demonstrated a very high degree of conservation in gene content and order with other sequenced plastid genomes. The lowland ecotype reference sequence (Kanlow Lin1) was 139,677 base pairs while the upland sequence (Summer Lin2) was 139,619 base pairs. Alignments between the lowland reference sequence and short-read sequence data from existing sequence datasets identified as either upland or lowland confirmed known polymorphisms and indicated the presence of other differences. Insertions and deletions principally occurred near stretches of homopolymer simple sequence repeats in intergenic regions while most Single Nucleotide Polymorphisms (SNPs) occurred in intergenic regions and introns within the single copy portions of the genome. The polymorphism rate between upland and lowland switchgrass ecotypes was found to be similar to rates reported between chloroplast genomes of indica and japonica subspecies of rice which were believed to have diverged 0.2–0.4 million years ago.

Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics

Here, we report a completely sequenced plastome using Illumina/Solexa sequencing-by-synthesis (SBS) technology. The plastome of Magnolia kwangsiensis Figlar & Noot. is 159 667 bp in length with a typical quadripartite structure: 88 030 bp large single-copy (LSC) and 18 669 bp small single-copy (SSC) regions, separated by two 26 484 bp inverted repeat (IR) regions. The overall predicted gene number is 129, among which 17 genes are duplicated in IR regions. The plastome of M. kwangsiensis is identical in its gene order to previously published plastomes of magnoliids. Furthermore, the C-to-U type RNA editing frequency of 114 seed plants is positively correlated with plastome GC content and plastome length, whereas plastome length is not correlated with GC content. A total of 16 potential putative barcoding or low taxonomic level phylogenetic study markers in Magnoliaceae were detected by comparing the coding and noncoding regions of the plastome of M. kwangsiensis with that of Liriodendron tulipifera L. At least eight markers might be applied not only to Magnoliaceae but also to other taxa. The 86 mononucleotide cpSSRs that distributed in single-copy noncoding regions are highly valuable to study population genetics and conservation genetics of this endangered rare species.

The ins and outs of editing and splicing of plastid RNAs: lessons from parasitic plants

In chloroplasts of higher plants, editing and splicing of transcripts is a prerequisite for the proper expression of the plastid genetic information and thereby for photosynthesis. Holoparasitic plants differ from photosynthetic plants in that they have abandoned a photoautotrophic life style, which has led to a reduction or loss of photosynthetic activity. The analysis of several parasitic plant plastid genomes revealed that coding capacities were reduced to different extent, encompassing genes that regulate plastid gene expression as well as photosynthesis genes. The reorganization of the plastid genome is also reflected in overall increases in point mutation rates that parallel the vanishing of RNA editing sites. Unprecedented in land plants is the parallel loss of the plastid gene coding for an intron maturase and all but one group IIa introns in two parasitic species. These observations highlight the plastome-wide effects that are associated with a relaxed evolutionary pressure in plants living a heterotrophic life style.

RESOPS: a database for analyzing the correspondence of RNA editing sites to protein three-dimensional structures

Transcripts from mitochondrial and chloroplast DNA of land plants often undergo cytidine to uridine conversion-type RNA editing events. RESOPS is a newly built database that specializes in displaying RNA editing sites of land plant organelles on protein three-dimensional (3D) structures to help elucidate the mechanisms of RNA editing for gene expression regulation. RESOPS contains the following information: unedited and edited cDNA sequences with notes for the target nucleotides of RNA editing, conceptual translation from the edited cDNA sequence in pseudo-UniProt format, a list of proteins under the influence of RNA editing, multiple amino acid sequence alignments of edited proteins, the location of amino acid residues coded by codons under the influence of RNA editing in protein 3D structures and the statistics of biased distributions of the edited residues with respect to protein structures. Most of the data processing procedures are automated; hence, it is easy to keep abreast of updated genome and protein 3D structural data. In the RESOPS database, we clarified that the locations of residues switched by RNA editing are significantly biased to protein structural cores. The integration of different types of data in the database also help advance the understanding of RNA editing mechanisms. RESOPS is accessible at http://cib.cf.ocha.ac.jp/RNAEDITING/.