Identification and expression profiling of microRNAs in leaf tissues of Foeniculum vulgare Mill. under salinity stress.

Foeniculum vulgare Mill. commonly known as fennel, is a globally recognized aromatic medicinal plant and culinary herb with widespread popularity due to its antimicrobial, antioxidant, carminative, and diuretic properties, among others. Although the phenotypic effects of salinity stress have been previously explored in fennel, the molecular mechanisms underlying responses to elevated salinity in this plant remain elusive. MicroRNAs (miRNAs) are tiny, endogenous, and extensively conserved non-coding RNAs (ncRNAs) typically ranging from 20 to 24 nucleotides (nt) in length that play a major role in a myriad of biological functions. In fact, a number of miRNAs have been extensively associated with responses to abiotic stress in plants. Consequently, employing computational methodologies and rigorous filtering criteria, 40 putative miRNAs belonging to 25 different families were characterized from fennel in this study. Subsequently, employing the psRNATarget tool, a total of 67 different candidate target transcripts for the characterized fennel miRNAs were predicted. Additionally, the expression patterns of six selected fennel miRNAs (i.e. fvu-miR156a, fvu-miR162a-3p, fvu-miR166a-3p, fvu-miR167a-5p, fvu-miR171a-3p, and fvu-miR408-3p) were analyzed under salinity stress conditions via qPCR. This article holds notable significance as it identifies not only 40 putative miRNAs in fennel, a non-model plant, but also pioneers the analysis of their expression under salinity stress conditions.

Characterisation of the Arabidopsis thaliana telomerase TERT-TR complex.

Most eukaryotic organisms employ a telomerase complex for the maintenance of chromosome ends. The core of this complex is composed of telomerase reverse transcriptase (TERT) and telomerase RNA (TR) subunits. The TERT reverse transcriptase (RT) domain synthesises telomeric DNA using the TR template sequence. The other TERT domains contribute to this process in different ways. In particular, the TERT RNA-binding domain (TRBD) interacts with specific TR motif(s). Using a yeast 3-hybrid system, we show the critical role of Arabidopsis thaliana (At) TRBD and embryophyta-conserved KRxR motif in the unstructured linker preceding the TRBD domain for binding to the recently identified AtTR subunit. We also show the essential role of the predicted P4 stem and pseudoknot AtTR structures and provide evidence for the binding of AtTRBD to pseudoknot and KRxR motif stabilising interaction with the P4 stem structure. Our results thus provide the first insight into the core part of the plant telomerase complex.

Straightforward and affordable agroinfiltration with RUBY accelerates RNA silencing research.

Transient expression and induction of RNA silencing by agroinfiltration is a fundamental method in plant RNA biology. Here, we introduce a new reporter assay using RUBY, which encodes three key enzymes of the betalain biosynthesis pathway, as a polycistronic mRNA. The red pigmentation conferred by betalains allows visual confirmation of gene expression or silencing levels without tissue disruption, and the silencing levels can be quantitatively measured by absorbance in as little as a few minutes. Infiltration of RUBY in combination with p19, a well-known RNA silencing suppressor, induced a fivefold higher accumulation of betalains at 7 days post infiltration compared to infiltration of RUBY alone. We demonstrated that co-infiltration of RUBY with two RNA silencing inducers, targeting either CYP76AD1 or glycosyltransferase within the RUBY construct, effectively reduces RUBY mRNA and betalain levels, indicating successful RNA silencing. Therefore, compared to conventional reporter assays for RNA silencing, the RUBY-based assay provides a simple and rapid method for quantitative analysis without the need for specialized equipment, making it useful for a wide range of RNA silencing studies.

Noncoding RNAs in regulation of plant secondary metabolism.

Plant secondary metabolites (PSMs) are a large class of structurally diverse molecules, mainly consisting of terpenoids, phenolic compounds, and nitrogen-containing compounds, which play active roles in plant development and stress responses. The biosynthetic processes of PSMs are governed by a sophisticated regulatory network at multiple levels. Noncoding RNAs (ncRNAs) such as microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs) may serve as post-transcriptional regulators for plant secondary metabolism through acting on genes encoding either transcription factors or participating enzymes in relevant metabolic pathways. High-throughput sequencing technologies have facilitated the large-scale identifications of ncRNAs potentially involved in plant secondary metabolism in model plant species as well as certain species with enriched production of specific types of PSMs. Moreover, a series of miRNA-target modules have been functionally characterized to be responsible for regulating PSM biosynthesis and accumulation in plants under abiotic or biotic stresses. In this review, we will provide an overview of current findings on the ncRNA-mediated regulation of plant secondary metabolism with special attention to its participation in plant stress responses, and discuss possible issues to be addressed in future fundamental research and breeding practice.

Mitochondrial RNA Helicases: Key Players in the Regulation of Plant Organellar RNA Splicing and Gene Expression.

Mitochondrial genomes of land plants are large and exhibit a complex mode of gene organization and expression, particularly at the post-transcriptional level. The primary organellar transcripts in plants undergo extensive maturation steps, including endo- and/or exo-nucleolytic cleavage, RNA-base modifications (mostly C-to-U deaminations) and both 'cis'- and 'trans'-splicing events. These essential processing steps rely on the activities of a large set of nuclear-encoded factors. RNA helicases serve as key players in RNA metabolism, participating in the regulation of transcription, mRNA processing and translation. They unwind RNA secondary structures and facilitate the formation of ribonucleoprotein complexes crucial for various stages of gene expression. Furthermore, RNA helicases are involved in RNA metabolism by modulating pre-mRNA maturation, transport and degradation processes. These enzymes are, therefore, pivotal in RNA quality-control mechanisms, ensuring the fidelity and efficiency of RNA processing and turnover in plant mitochondria. This review summarizes the significant roles played by helicases in regulating the highly dynamic processes of mitochondrial transcription, RNA processing and translation in plants. We further discuss recent advancements in understanding how dysregulation of mitochondrial RNA helicases affects the splicing of organellar genes, leading to respiratory dysfunctions, and consequently, altered growth, development and physiology of land plants.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function.

Reproductive phasiRNAs (phased, small interfering RNAs) are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt (nucleotides) phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at the premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 (Dicer-like 5) for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely i) not triggered by microRNAs, ii) not loaded by AGO18 proteins, and iii) not capable of mediating PHAS precursor cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

A conserved Pol II elongator SPT6L mediates Pol V transcription to regulate RNA-directed DNA methylation in Arabidopsis.

In plants, the plant-specific RNA polymerase V (Pol V) transcripts non-coding RNAs and provides a docking platform for the association of accessory proteins in the RNA-directed DNA methylation (RdDM) pathway. Various components have been uncovered that are involved in the process of DNA methylation, but it is still not clear how the transcription of Pol V is regulated. Here, we report that the conserved RNA polymerase II (Pol II) elongator, SPT6L, binds to thousands of intergenic regions in a Pol II-independent manner. The intergenic enrichment of SPT6L, interestingly, co-occupies with the largest subunit of Pol V (NRPE1) and mutation of SPT6L leads to the reduction of DNA methylation but not Pol V enrichment. Furthermore, the association of SPT6L at Pol V loci is dependent on the Pol V associated factor, SPT5L, rather than the presence of Pol V, and the interaction between SPT6L and NRPE1 is compromised in spt5l. Finally, Pol V RIP-seq reveals that SPT6L is required to maintain the amount and length of Pol V transcripts. Our findings thus uncover the critical role of a Pol II conserved elongator in Pol V mediated DNA methylation and transcription, and shed light on the mutual regulation between Pol V and II in plants.

Exploring the dynamic adaptive responses of Epimedium pubescens to phosphorus deficiency by Integrated transcriptome and miRNA analysis.

Phosphorus, a crucial macronutrient essential for plant growth and development. Due to widespread phosphorus deficiency in soils, phosphorus deficiency stress has become one of the major abiotic stresses that plants encounter. Despite the evolution of adaptive mechanisms in plants to address phosphorus deficiency, the specific strategies employed by species such as Epimedium pubescens remain elusive. Therefore, this study observed the changes in the growth, physiological reponses, and active components accumulation in E. pubescensunder phosphorus deficiency treatment, and integrated transcriptome and miRNA analysis, so as to offer comprehensive insights into the adaptive mechanisms employed by E. pubescens in response to phosphorus deficiency across various stages of phosphorus treatment. Remarkably, our findings indicate that phosphorus deficiency induces root growth stimulation in E. pubescens, while concurrently inhibiting the growth of leaves, which are of medicinal value. Surprisingly, this stressful condition results in an augmented accumulation of active components in the leaves. During the early stages (30 days), leaves respond by upregulating genes associated with carbon metabolism, flavonoid biosynthesis, and hormone signaling. This adaptive response facilitates energy production, ROS scavenging, and morphological adjustments to cope with short-term phosphorus deficiency and sustain its growth. As time progresses (90 days), the expression of genes related to phosphorus cycling and recycling in leaves is upregulated, and transcriptional and post-transcriptional regulation (miRNA regulation and protein modification) is enhanced. Simultaneously, plant growth is further suppressed, and it gradually begins to discard and decompose leaves to resist the challenges of long-term phosphorus deficiency stress and sustain survival. In conclusion, our study deeply and comprehensively reveals adaptive strategies utilized by E. pubescens in response to phosphorus deficiency, demonstrating its resilience and thriving potential under stressful conditions. Furthermore, it provides valuable information on potential target genes for the cultivation of E. pubescens genotypes tolerant to low phosphorus.

Maize miRNAs and their putative target genes involved in chilling stress response in 5-day old seedlings.

BACKGROUND: In the context of early sowing of maize as a promising adaptation strategy that could significantly reduce the negative effects of climate change, an in-depth understanding of mechanisms underlying plant response to low-temperature stress is demanded. Although microRNAs (miRNAs) have been recognized as key regulators of plant stress response, research on their role in chilling tolerance of maize during early seedling stages is scarce. Therefore, it is of great significance to explore chilling-responsive miRNAs, reveal their expression patterns and associated target genes, as well as to examine the possible functions of the conserved and novel miRNAs. In this study, the role of miRNAs was examined in 5d-old maize seedlings of one tolerant and one sensitive inbred line exposed to chilling (10/8 °C) stress for 6 h and 24 h, by applying high throughput sequencing. RESULTS: A total of 145 annotated known miRNAs belonging to 30 families and 876 potentially novel miRNAs were identified. Differential expression (DE) analysis between control and stress conditions identified 98 common miRNAs for both genotypes at one time point and eight miRNAs at both time points. Target prediction and enrichment analysis showed that the DE zma-miR396, zma-miR156, zma-miR319, and zma-miR159 miRNAs modulate growth and development. Furthermore, it was found that several other DE miRNAs were involved in abiotic stress response: antioxidative mechanisms (zma-miR398), signal transduction (zma-miR156, zma-miR167, zma-miR169) and regulation of water content (zma-miR164, zma-miR394, zma-miR396). The results underline the zma-miRNAs involvement in the modulation of their target genes expression as an important aspect of the plant's survival strategy and acclimation to chilling stress conditions. CONCLUSIONS: To our understanding, this is the first study on miRNAs in 5-d old seedlings' response to chilling stress, providing data on the role of known and novel miRNAs post-transcriptional regulation of expressed genes and contributing a possible platform for further network and functional analysis.

Improving plant miRNA-target prediction with self-supervised k-mer embedding and spectral graph convolutional neural network.

Deciphering the targets of microRNAs (miRNAs) in plants is crucial for comprehending their function and the variation in phenotype that they cause. As the highly cell-specific nature of miRNA regulation, recent computational approaches usually utilize expression data to identify the most physiologically relevant targets. Although these methods are effective, they typically require a large sample size and high-depth sequencing to detect potential miRNA-target pairs, thereby limiting their applicability in improving plant breeding. In this study, we propose a novel miRNA-target prediction framework named kmerPMTF (k-mer-based prediction framework for plant miRNA-target). Our framework effectively extracts the latent semantic embeddings of sequences by utilizing k-mer splitting and a deep self-supervised neural network. We construct multiple similarity networks based on k-mer embeddings and employ graph convolutional networks to derive deep representations of miRNAs and targets and calculate the probabilities of potential associations. We evaluated the performance of kmerPMTF on four typical plant datasets: Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, and Prunus persica. The results demonstrate its ability to achieve AUPRC values of 84.9%, 91.0%, 80.1%, and 82.1% in 5-fold cross-validation, respectively. Compared with several state-of-the-art existing methods, our framework achieves better performance on threshold-independent evaluation metrics. Overall, our study provides an efficient and simplified methodology for identifying plant miRNA-target associations, which will contribute to a deeper comprehension of miRNA regulatory mechanisms in plants.

Tracking the messengers: Emerging advances in mRNA-based plant communication.

Messenger RNAs (mRNAs) are the templates for protein translation but can also act as non-cell-autonomous signaling molecules. Plants input endogenous and exogenous cues to mobile mRNAs and output them to local or systemic target cells and organs to support specific plant responses. Mobile mRNAs form ribonucleoprotein (RNP) complexes with proteins during transport. Components of these RNP complexes could interact with plasmodesmata (PDs), a major mediator of mRNA transport, to ensure mRNA mobility and transport selectivity. Based on advances in the last two to three years, this review summarizes mRNA transport mechanisms in local and systemic signaling from the perspective of RNP complex formation and PD transport. We also discuss the physiological roles of endogenous mRNA transport and the recently revealed roles of non-cell-autonomous mRNAs in inter-organism communication.

Dissecting the molecular puzzle of the editosome core in Arabidopsis organelles.

Over the last decade, the composition of the C-to-U RNA editing complex in embryophyte organelles has turned out to be much more complex than first expected. While PPR proteins were initially thought to act alone, significant evidences have clearly depicted a sophisticated mechanism with numerous protein-protein interaction involving PPR and non-PPR proteins. Moreover, the identification of specific functional partnership between PPRs also suggests that, in addition to the highly specific PPRs directly involved in the RNA target recognition, non-RNA-specific ones are required. Although some of them, such as DYW1 and DYW2, were shown to be the catalytic domains of the editing complex, the molecular function of others, such as NUWA, remains elusive. It was suggested that they might stabilize the complex by acting as a scaffold. We here performed functional complementation of the crr28-2 mutant with truncated CRR28 proteins mimicking PPR without the catalytic domain and show that they exhibit a specific dependency to one of the catalytic proteins DYW1 or DYW2. Moreover, we also characterized the role of the PPR NUWA in the editing reaction and show that it likely acts as a scaffolding factor. NUWA is no longer required for efficient editing of the CLB19 editing sites once this RNA specific PPR is fused to the DYW catalytic domain of its partner DYW2. Altogether, our results strongly support a flexible, evolutive and resilient editing complex in which RNA binding activity, editing activity and stabilization/scaffolding function can be provided by one or more PPRs.

Molecular basis of A. thaliana KEOPS complex in biosynthesizing tRNA t6A.

In archaea and eukaryotes, the evolutionarily conserved KEOPS is composed of four core subunits-Kae1, Bud32, Cgi121 and Pcc1, and a fifth Gon7/Pcc2 that is found in fungi and metazoa. KEOPS cooperates with Sua5/YRDC to catalyze the biosynthesis of tRNA N6-threonylcarbamoyladenosine (t6A), an essential modification needed for fitness of cellular organisms. Biochemical and structural characterizations of KEOPSs from archaea, yeast and humans have determined a t6A-catalytic role for Kae1 and auxiliary roles for other subunits. However, the precise molecular workings of KEOPSs still remain poorly understood. Here, we investigated the biochemical functions of A. thaliana KEOPS and determined a cryo-EM structure of A. thaliana KEOPS dimer. We show that A. thaliana KEOPS is composed of KAE1, BUD32, CGI121 and PCC1, which adopts a conserved overall arrangement. PCC1 dimerization leads to a KEOPS dimer that is needed for an active t6A-catalytic KEOPS-tRNA assembly. BUD32 participates in direct binding of tRNA to KEOPS and modulates the t6A-catalytic activity of KEOPS via its C-terminal tail and ATP to ADP hydrolysis. CGI121 promotes the binding of tRNA to KEOPS and potentiates the t6A-catalytic activity of KEOPS. These data and findings provide insights into mechanistic understanding of KEOPS machineries.

ARGONAUTE1-binding Tudor domain proteins function in small interfering RNA production for RNA-directed DNA methylation.

Transposable elements (TEs) contribute to plant evolution, development, and adaptation to environmental changes, but the regulatory mechanisms are largely unknown. RNA-directed DNA methylation (RdDM) is 1 TE regulatory mechanism in plants. Here, we identified that novel ARGONAUTE 1 (AGO1)-binding Tudor domain proteins Precocious dissociation of sisters C/E (PDS5C/E) are involved in 24-nt siRNA production to establish RdDM on TEs in Arabidopsis thaliana. PDS5 family proteins are subunits of the eukaryote-conserved cohesin complex. However, the double mutant lacking angiosperm-specific subfamily PDS5C and PDS5E (pds5c/e) exhibited different developmental phenotypes and transcriptome compared with those of the double mutant lacking eukaryote-conserved subfamily PDS5A and PDS5B (pds5a/b), suggesting that the angiosperm-specific PDS5C/E subfamily has a unique function in angiosperm plants. Proteome and imaging analyses revealed that PDS5C/E interact with AGO1. The pds5c/e double mutant had defects in 24-nt siRNA accumulation and CHH DNA methylation on TEs. In addition, some lncRNAs that accumulated in the pds5c/e mutant were targeted by AGO1-loading 21-nt miRNAs and 21-nt siRNAs. These results indicate that PDS5C/E and AGO1 participate in 24-nt siRNA production for RdDM in the cytoplasm. These findings indicate that angiosperm plants evolved a new regulator, the PDS5C/E subfamily, to control the increase in TEs during angiosperm evolution.

Pleiotropic Effects of miR5504 Underlying Plant Height, Grain Yield and Quality in Rice.

MicroRNAs (miRNAs) are known to play critical roles in regulating rice agronomic traits through mRNA cleavage or translational repression. Our previous study indicated that miR5504 regulates plant height by affecting cell proliferation and expansion. Here, the two independent homozygous mir5504 mutants (CR1 and CR2) and overexpression lines (OE1 and OE2) were further used to investigate the functions of miR5504. The panicle length, 1000-grain weight and grain yield per plant of miR5504-OE lines were identical to those of Nipponbare (NIP), but the 1000-grain weight of mir5504 mutants was reduced by about 10% and 9%, respectively. Meanwhile, the grain width and thickness of mir5504 mutants decreased significantly by approximately 10% and 11%, respectively. Moreover, the cytological results revealed a significant decrease in cell number along grain width direction and cell width in spikelet in mir5504, compared with those in NIP. In addition, several major storage substances of the rice seeds were measured. Compared to NIP, the amylose content of the mir5504 seeds was noticeably decreased, leading to an increase of nearly 10 mm in gel consistency (GC) in mir5504 lines. Further investigation confirmed that LOC_Os08g16914 was the genuine target of miR5504: LOC_Os08g16914 over-expression plants phenocopied the mir5504 mutants. This study provides insights into the role of miR5504 in rice seed development.

A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants.

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

Plant non-coding RNAs: The new frontier for the regulation of plant development and adaptation to stress.

Most plant transcriptomes constitute functional non-coding RNAs (ncRNAs) that lack the ability to encode proteins. In recent years, more research has demonstrated that ncRNAs play important regulatory roles in almost all plant biological processes by modulating gene expression. Thus, it is important to study the biogenesis and function of ncRNAs, particularly in plant growth and development and stress tolerance. In this review, we systematically explore the process of formation and regulatory mechanisms of ncRNAs, particularly those of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, we provide a comprehensive overview of the recent advancements in ncRNAs research, including their regulation of plant growth and development (seed germination, root growth, leaf morphogenesis, floral development, and fruit and seed development) and responses to abiotic and biotic stress (drought, heat, cold, salinity, pathogens and insects). We also discuss research challenges and provide recommendations to advance the understanding of the roles of ncRNAs in agronomic applications.

A long noncoding RNA functions in pumpkin fruit development through S-adenosyl-L-methionine synthetase.

Long noncoding RNAs (lncRNAs) play important roles in various biological processes. However, the regulatory roles of lncRNAs underlying fruit development have not been extensively studied. The pumpkin (Cucurbita spp.) is a preferred model for understanding the molecular mechanisms regulating fruit development because of its variable shape and size and large inferior ovary. Here, we performed strand-specific transcriptome sequencing on pumpkin (Cucurbita maxima "Rimu") fruits at 6 developmental stages and identified 5,425 reliably expressed lncRNAs. Among the 332 lncRNAs that were differentially expressed during fruit development, the lncRNA MSTRG.44863.1 was identified as a negative regulator of pumpkin fruit development. MSTRG.44863.1 showed a relatively high expression level and an obvious period-specific expression pattern. Transient overexpression and silencing of MSTRG.44863.1 significantly increased and decreased the content of 1-aminocyclopropane carboxylic acid (a precursor of ethylene) and ethylene production, respectively. RNA pull-down and microscale thermophoresis assays further revealed that MSTRG.44863.1 can interact with S-adenosyl-L-methionine synthetase (SAMS), an enzyme in the ethylene synthesis pathway. Considering that ethylene negatively regulates fruit development, these results indicate that MSTRG.44863.1 plays an important role in the regulation of pumpkin fruit development, possibly through interacting with SAMS and affecting ethylene synthesis. Overall, our findings provide a rich resource for further study of fruit-related lncRNAs while offering insights into the regulation of fruit development in 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.