Frontiers in plant RNA research in ICAR2023: from lab to innovative agriculture

The recent growth in global warming, soil contamination, and climate instability have widely disturbed ecosystems, and will have a significant negative impact on the growth of plants that produce grains, fruits and woody biomass. To conquer this difficult situation, we need to understand the molecular bias of plant environmental responses and promote development of new technologies for sustainable maintenance of crop production. Accumulated molecular biological data have highlighted the importance of RNA-based mechanisms for plant stress responses. Here, we report the most advanced plant RNA research presented in the 33rd International Conference on Arabidopsis Research (ICAR2023), held as a hybrid event on June 5-9, 2023 in Chiba, Japan, and focused on "Arabidopsis for Sustainable Development Goals". Six workshops/concurrent sessions in ICAR2023 targeted plant RNA biology, and many RNA-related topics could be found in other sessions. In this meeting report, we focus on the workshops/concurrent sessions targeting RNA biology, to share what is happening now at the forefront of plant RNA research.

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.

Profiling of RNA editing events in plant organellar transcriptomes with high-throughput sequencing

RNA editing is a crucial post-transcriptional modification process in plant organellar RNA metabolism. rRNA removal-based total RNA-seq is one of the most common methods to study this event. However, the lack of commercial kits to remove rRNAs limits the usage of this method, especially for non-model plant species. DSN-seq is a transcriptome sequencing method utilizing duplex-specific nuclease (DSN) to degrade highly abundant cDNA species especially those from rRNAs while keeping the robustness of transcript levels of the majority of other mRNAs, and has not been applied to study RNA editing in plants before. In this study, we evaluated the capability of DSN-seq to reduce rRNA content and profile organellar RNA editing events in plants, as well we used commercial Ribo-off-seq and standard mRNA-seq as comparisons. Our results demonstrated that DSN-seq efficiently reduced rRNA content and enriched organellar transcriptomes in rice. With high sensitivity to RNA editing events, DSN-seq and Ribo-off-seq provided a more complete and accurate RNA editing profile of rice, which was further validated by Sanger sequencing. Furthermore, DSN-seq also demonstrated efficient organellar transcriptome enrichment and high sensitivity for profiling RNA editing events in Arabidopsis thaliana. Our study highlights the capability of rRNA removal-based total RNA-seq for profiling RNA editing events in plant organellar transcriptomes and also suggests DSN-seq as a widely accessible RNA editing profiling method for various plant species.

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

MicroRNA482/2118 is lineage-specifically involved in gibberellin signalling via the regulation of GID1 expression by targeting noncoding PHAS genes and subsequently instigated phasiRNAs

MicroRNA482/2118 (miR482/2118) is a 22-nt miRNA superfamily, with conserved functions in disease resistance and plant development. It usually instigates the production of phased small interfering RNAs (phasiRNAs) from its targets to expand or reinforce its silencing effect. Using a new high-quality reference genome sequence and comprehensive small RNA profiling, we characterized a newly evolved regulatory pathway of miR482/2118 in litchi. In this pathway, miR482/2118 cleaved a novel noncoding trans-acting gene (LcTASL1) and triggered phasiRNAs to regulate the expression of gibberellin (GA) receptor gene GIBBERELLIN INSENSITIVE DWARF1 (GID1) in trans; another trans-acting gene LcTASL2, targeted by LcTASL1-derived phasiRNAs, produced phasiRNAs as well to target LcGID1 to reinforce the silencing effect of LcTASL1. We found this miR482/2118-TASL-GID1 pathway was likely involved in fruit development, especially the seed development in litchi. In vivo construction of the miR482a-TASL-GID1 pathway in Arabidopsis could lead to defects in flower and silique development, analogous to the phenotype of gid1 mutants. Finally, we found that a GA-responsive transcription factor, LcGAMYB33, could regulate LcMIR482/2118 as a feedback mechanism of the sRNA-silencing pathway. Our results deciphered a lineage-specifically evolved regulatory module of miR482/2118, demonstrating the high dynamics of miR482/2118 function in plants.

Mating system is associated with seed phenotypes upon loss of RNA-directed DNA methylation in Brassicaceae

In plants, de novo DNA methylation is guided by 24-nt short interfering (si)RNAs in a process called RNA-directed DNA methylation (RdDM). Primarily targeted at transposons, RdDM causes transcriptional silencing and can indirectly influence expression of neighboring genes. During reproduction, a small number of siRNA loci are dramatically upregulated in the maternally derived seed coat, suggesting that RdDM might have a special function during reproduction. However, the developmental consequence of RdDM has been difficult to dissect because disruption of RdDM does not result in overt phenotypes in Arabidopsis (Arabidopsis thaliana), where the pathway has been most thoroughly studied. In contrast, Brassica rapa mutants lacking RdDM have a severe seed production defect, which is determined by the maternal sporophytic genotype. To explore the factors that underlie the different phenotypes of these species, we produced RdDM mutations in 3 additional members of the Brassicaceae family: Camelina sativa, Capsella rubella, and Capsella grandiflora. Among these 3 species, only mutations in the obligate outcrosser, C. grandiflora, displayed a seed production defect similar to Brassica rapa mutants, suggesting that mating system is a key determinant for reproductive phenotypes in RdDM mutants.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

KEY MESSAGE

Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

RIP-Seq analysis of non-PPR chloroplast editing factors reveals broad RNA interactions and enrichment of less efficiently translated RNAs by OZ1 and ORRM1 complexes

C-to-U RNA editing in angiosperm chloroplasts requires a large suite of proteins bound together in the editosome. The editosome is comprised of PPR proteins, RIP/MORFs, OZ proteins, and ORRM proteins that physically interact in high molecular weight complexes. The specific functions of non-PPR editing factors in the editosome are unclear, however, specific subsets of editing sites are affected by absence of non-PPR editing factors. Unlike the PPR components of editosomes that have predictable nucleotide specificities, domains present in non-PPR editing factors make RNA associations difficult to predict. In this study, chloroplast extracts were isolated from juvenile maize seedlings. RNAs were immunoprecipitated using polyclonal antibodies targeting non-PPR editing factors RIP9, OZ1, and ORRM1. RNA libraries from duplicate experiments were compared. RIP9 was associated with most of the non-ribosomal RNA content of chloroplasts, consistent with a general binding function to PPR L-motifs and tethering of large ribonucleoprotein complexes. The breadth of RNA associations was greater than predicted and include mRNAs without predicted editing sites, tRNA sequences, and introns. OZ1 and ORRM1 were associated with a highly similar pool of RNAs that have a bias toward lower translational efficiency values in mature chloroplasts. Lower translational efficiency was also associated with the pool of edited RNAs compared to RNAs without editing sites. The unexpected breadth of interactions by non-PPR editing factors suggests the editosome is large, diverse, and associated with RNAs with lower relative translational efficiency in mature chloroplasts.

Systematic identification and characterization of microRNAs with target genes involved in high night temperature stress at the filling stage of rice

High night temperature stress is one of the main environmental factors affecting rice yield and quality. More and more evidence shows that microRNA (miRNA) plays an important role in various abiotic stresses. However, the molecular network of miRNA regulation on rice tolerance to high night temperatures remains unclear. Here, small RNA, transcriptome and degradome sequencing were integrated to identify differentially expressed miRNAs, genes, and key miRNA-target gene pairs in rice heat-sensitive and heat-tolerant lines at the filling stage suffering from high night temperature stress. It was discovered that there were notable differences in the relative expression of 102 miRNAs between the two rice lines under stress. Meanwhile, 5263 and 5405 mRNAs were differentially expressed in the heat-sensitive line and heat-tolerant line, and functional enrichment analysis revealed that these genes were involved in heat-related processes and pathways. The miRNAs-mRNAs target relationship was further verified by degradome sequencing. Eventually, 49 miRNAs-222 mRNAs target pairs with reverse expression patterns showed significant relative expression changes between the heat-tolerant and the heat-sensitive line, being suggested to be responsible for the heat tolerance difference of these two rice lines. Functional analysis of these 222 mRNA transcripts showed that high night temperature-responsive miRNAs targeted these mRNAs involved in many heat-related biological processes, such as transcription regulation, chloroplast regulation, mitochondrion regulation, protein folding, hormone regulation and redox process. This study identified possible miRNA-mRNA regulation relationships in response to high night temperature stress in rice and potentially contributed to heat resistance breeding of rice in the future.

Molecular mechanisms of the RNA polymerases in plant RNA-directed DNA methylation

In plants, two atypical DNA-dependent RNA polymerases, RNA polymerase IV (Pol IV) and Pol V, and an RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) together produce noncoding RNAs (ncRNAs) to guide the plant-specific RNA-directed DNA methylation (RdDM). Although both Pol IV and Pol V have evolved from the canonical Pol II, they have adapted to different roles in RdDM. The mechanisms of their adaptation are key to understanding plant DNA methylation and the divergent evolution of polymerases. In this review, we summarize insights that have emerged from recent structural studies of Pol IV, Pol V, and RDR2 and discuss their structural features critical for efficient ncRNA production in RdDM.

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.

Harnessing the potential of non-coding RNA: An insight into its mechanism and interaction in plant biotic stress

Plants have developed diverse physical and chemical defence mechanisms to ensure their continued growth and well-being in challenging environments. Plants also have evolved intricate molecular mechanisms to regulate their responses to biotic stress. Non-coding RNA (ncRNA) plays a crucial role in this process that affects the expression or suppression of target transcripts. While there have been numerous reviews on the role of molecules in plant biotic stress, few of them specifically focus on how plant ncRNAs enhance resistance through various mechanisms against different pathogens. In this context, we explored the role of ncRNA in exhibiting responses to biotic stress endogenously as well as cross-kingdom regulation of transcript expression. Furthermore, we address the interplay between ncRNAs, which can act as suppressors, precursors, or regulators of other ncRNAs. We also delve into the regulation of ncRNAs in response to attacks from different organisms, such as bacteria, viruses, fungi, nematodes, oomycetes, and insects. Interestingly, we observed that diverse microorganisms interact with distinct ncRNAs. This intricacy leads us to conclude that each ncRNA serves a specific function in response to individual biotic stimuli. This deeper understanding of the molecular mechanisms involving ncRNAs in response to biotic stresses enhances our knowledge and provides valuable insights for future research in the field of ncRNA, ultimately leading to improvements in plant traits.

Nonsense-mediated mRNA decay pathway in plants under stress: general gene regulatory mechanism and advances

MAIN CONCLUSION

Nonsense-mediated mRNA decay in eukaryotes is vital to cellular homeostasis. Further knowledge of its putative role in plant RNA metabolism under stress is pivotal to developing fitness-optimizing strategies. Nonsense-mediated mRNA decay (NMD), part of the mRNA surveillance pathway, is an evolutionarily conserved form of gene regulation in all living organisms. Degradation of mRNA-bearing premature termination codons and regulation of physiological RNA levels highlight NMD's role in shaping the cellular transcriptome. Initially regarded as purely a tool for cellular RNA quality control, NMD is now considered to mediate various aspects of plant developmental processes and responses to environmental changes. Here we offer a basic understanding of NMD in eukaryotes by explaining the concept of premature termination codon recognition and NMD complex formation. We also provide a detailed overview of the NMD mechanism and its role in gene regulation. The potential role of effectors, including ABCE1, in ribosome recycling during the translation process is also explained. Recent reports of alternatively spliced variants of corresponding genes targeted by NMD in Arabidopsis thaliana are provided in tabular format. Detailed figures are also provided to clarify the NMD concept in plants. In particular, accumulating evidence shows that NMD can serve as a novel alternative strategy for genetic manipulation and can help design RNA-based therapies to combat stress in plants. A key point of emphasis is its function as a gene regulatory mechanism as well as its dynamic regulation by environmental and developmental factors. Overall, a detailed molecular understanding of the NMD mechanism can lead to further diverse applications, such as improving cellular homeostasis in living organisms.

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

PRMD: an integrated database for plant RNA modifications

The scope and function of RNA modifications in model plant systems have been extensively studied, resulting in the identification of an increasing number of novel RNA modifications in recent years. Researchers have gradually revealed that RNA modifications, especially N6-methyladenosine (m6A), which is one of the most abundant and commonly studied RNA modifications in plants, have important roles in physiological and pathological processes. These modifications alter the structure of RNA, which affects its molecular complementarity and binding to specific proteins, thereby resulting in various of physiological effects. The increasing interest in plant RNA modifications has necessitated research into RNA modifications and associated datasets. However, there is a lack of a convenient and integrated database with comprehensive annotations and intuitive visualization of plant RNA modifications. Here, we developed the Plant RNA Modification Database (PRMD; http://bioinformatics.sc.cn/PRMD and http://rnainformatics.org.cn/PRMD) to facilitate RNA modification research. This database contains information regarding 20 plant species and provides an intuitive interface for displaying information. Moreover, PRMD offers multiple tools, including RMlevelDiff, RMplantVar, RNAmodNet and Blast (for functional analyses), and mRNAbrowse, RNAlollipop, JBrowse and Integrative Genomics Viewer (for displaying data). Furthermore, PRMD is freely available, making it useful for the rapid development and promotion of research on plant RNA modifications.

Fluorescent In Situ Detection of Small RNAs in Plants Using sRNA-FISH

Plant small RNAs are 21-24 nucleotide, noncoding RNAs that function as regulators in plant growth and development. Colorimetric detection of plant small RNAs was made possible with the introduction of locked nucleic acid probes. However, fluorescent detection of plant small RNAs has been challenging due to the high autofluorescence from plant tissue. Here we report a fluorescent in situ detection method for plant small RNAs. This method can be applied to most plant samples and tissue types and also can be adapted for single-molecule detection of small RNAs with super-resolution microscopy.

RNA Isolation Method in Marginal Crops with High Agronomic Potential

Ribonucleic Acid (RNA) isolation is a basic technique in the field of molecular biology. The purpose of RNA isolation is to acquire pure and complete RNA that can be used to evaluate gene expression. Many methods can be used to perform RNA isolation, all of them based on the chemical properties of nucleic acids. However, some of them do not achieve high RNA yields and purity levels when used in a number of marginally studied crops of agronomic importance, such as grain and vegetable amaranth plants. In the method described here, the use of guanidinium thiocyanate and two additional precipitation steps with different reagents designed to obtain high yields and RNA purity levels from diverse plant species employed for plant functional genomics studies is described.

A universal protocol for high-quality DNA and RNA isolation from diverse plant species

Next-generation sequencing demands high-quality nucleic acid, yet isolating DNA and RNA is often challenging, particularly from plant tissues. Despite advances in developing various kits and reagents, these products are tailored to isolation of nucleic acid from model plant tissues. Here we introduce a universal lysis buffer to separate nucleic acid from various plant species, including recalcitrant plants, to facilitate molecular analyses, such as quantitative PCR (qPCR), transcriptomics, and whole-genome sequencing (WGS). The protocol is a modification of the original CTAB methods, which leads to nucleic acid isolation from many plant species, including monocots and eudicots. The lysis buffer consists of hexadecyltrimethylammonium bromide (CTAB), sodium chloride (NaCl), Tris base, ethylenediaminetetraacetic acid (EDTA) and β-mercaptoethanol (βME). The modified CTAB method enables the isolation of nucleic acid from small amounts of plant tissues (e.g., 15-100 mg) in a timely manner, which is well-suited for a large number of samples and also when adequate sample collection is a limiting factor. The protocol isolates not only DNA from various plant species but also RNA. This makes it highly effective for molecular analyses compared to previously described CTAB methods optimised for DNA isolation. The appropriate concentration of the components enables high-quality DNA and RNA isolation from plant tissues simultaneously. Additionally, this protocol is compatible with commercially available columns. For DNA and RNA to be qualified for next-generation sequencing platforms, the protocol is supplemented with columns to purify either DNA or RNA from the same tissue to meet high standards for sequencing analyses. This protocol provides an ideal approach to overcome potential obstacles in isolating high-quality DNA or RNA from a wide range of plant species for downstream molecular analysis.

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

Identification of microRNAs involved in ectomycorrhizal formation in Populus tomentosa

The majority of woody plants are able to form ectomycorrhizal (ECM) symbioses with fungi. During symbiotic development, plants undergo a complex re-programming process involving a series of physiological and morphological changes. MicroRNAs (miRNAs) are important components of the regulatory network underlying symbiotic development. To elucidate the mechanisms of miRNAs and miRNA-mediated mRNA cleavage during symbiotic development, we conducted high-throughput sequencing of small RNAs and degradome tags from roots of Populus tomentosa inoculated with Cenococcum geophilum. This process led to the annotation of 51 differentially expressed miRNAs between non-mycorrhizal and mycorrhizal roots of P. tomentosa, including 13 novel miRNAs. Increased or decreased accumulation of several novel and conserved miRNAs in ECM roots, including miR162, miR164, miR319, miR396, miR397, miR398, novel-miR44 and novel-miR47, suggests essential roles for these miRNAs in ECM formation. The degradome analysis identified root transcripts as miRNA-mediated mRNA cleavage targets, which was confirmed using real-time quantitative PCR. Several of the identified miRNAs and corresponding targets are involved in arbuscular mycorrhizal symbioses. In summary, increased or decreased accumulation of specific miRNAs and miRNA-mediated cleavage of symbiosis-related genes indicate that miRNAs play important roles in the regulatory network underlying symbiotic development.

PhasiHunter: a robust phased siRNA regulatory cascade mining tool based on multiple reference sequences

SUMMARY

In recent years, phased small interfering RNA has been found to play crucial roles in many biological processes in plants. However, efficiently predicting phasiRNA regulatory cascades with computational methods is still challenging. Here, we introduce PhasiHunter, a phasiRNA regulatory network prediction tool that has several distinctive features compared to existing tools: (i) PhasiHunter employs two major phasiRNA prediction algorithms, namely phase score and hypergeometric distribution-based methods, to ensure the integrity and accuracy of prediction; (ii) PhasiHunter can identify phasiRNAs and their regulatory networks based on multiple reference sequences and the predicted results can be automatically integrated; (iii) PhasiHunter can efficiently identify the phasiRNAs generated through alternative splicing events; and (iv) the excellent data structure and parallel computing architecture allow PhasiHunter to predict phasiRNAs and their regulatory pathways with high efficiency.

AVAILABILITY AND IMPLEMENTATION

PhasiHunter is an open-source tool that is available at https://github.com/HuangLab-CBI/PhasiHunter.

Plant immunity suppressor SKRP encodes a novel RNA-binding protein that targets exon 3' end of unspliced RNA

The regulatory roles of RNA splicing in plant immunity are emerging but still largely obscure. We reported previously that Phytophthora pathogen effector Avr3c targets a soybean protein SKRP (serine/lysine/arginine-rich protein) to impair soybean basal immunity by regulating host pre-mRNA alternative splicing, while the biochemical nature of SKRP remains unknown. Here, by using Arabidopsis as a model, we studied the mechanism of SKRP in regulating pre-mRNA splicing and plant immunity. AtSKRP confers impaired plant immunity against Phytophthora capsici and associates with spliceosome component PRP8 and splicing factor SR45, which positively and negatively regulate plant immunity, respectively. Enhanced crosslinking and immunoprecipitation followed by high-throughput sequencing (eCLIP-seq) showed AtSKRP is a novel RNA-binding protein that targets exon 3' end of unspliced RNA. Such position-specific binding of SKRP is associated with its activity in suppressing intron retention, including at positive immune regulatory genes UBP25 and RAR1. In addition, we found AtSKRP self-interact and forms oligomer, and these properties are associated with its function in plant immunity. Overall, our findings reveal that the immune repressor SKRP is a spliceosome-associated protein that targets exon 3' end to regulate pre-mRNA splicing in Arabidopsis.

Small RNA sequencing provides insights into molecular mechanism of flower development in Rhododendron pulchrum Sweet

Rhododendron pulchrum sweet, a member of the Ericaceae family possessing valuable horticultural properties, is widely distributed in the temperate regions. Though serving as bioindicator of metal pollution, the molecular mechanism regulating flowering in R. pulchrum is very limited. Illumina sequencing was performed to identify critical miRNAs in the synthesis of flavonoids at different developmental stages. Totally, 722 miRNAs belonging to 104 families were screened, and 84 novel mature miRNA sequences were predicted. The miR166, miR156, and miR167-1 families were dominant. In particular, 126 miRNAs were significantly differentially expressed among four different flowering stages. Totally, 593 genes were differentially regulated by miRNAs during the flower development process, which were mostly involved in "metabolic pathways", "plant hormone signal transduction", and "mitosis and regulation of biosynthetic processes". In pigment biosynthesis and signal transduction processes, gra-miR750 significantly regulated the expression of flavonoid 3',5'-hydroxylase; aof-miR171a, aof-miR171b, aof-miR171c, cas-miR171a-3p, and cas-miR171c-3p could regulate the expression of DELLA protein; aof-miR390, aof-miR396b, ath-miR3932b-5p, cas-miR171a-3p, aof-miR171a, and aof-miR171b regulated BAK1 expression. This research showed great potentials for genetic improvement of flower color traits for R. pulchrum and other Rhododendron species.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

An mTRAN-mRNA interaction mediates mitochondrial translation initiation in plants

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that "mitochondrial translation factors" mTRAN1 and mTRAN2 are land plant-specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)-like RNA binding proteins of the mitoribosomal "small" subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5' regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

GmGAMYB-BINDING PROTEIN 1 promotes small auxin-up RNA gene transcription to modulate soybean maturity and height

Flowering time, maturity, and plant height are crucial agronomic traits controlled by photoperiod that affect soybean (Glycine max [L.] Merr.) yield and regional adaptability. It is important to cultivate soybean cultivars of earlier maturity that adapt to high latitudes. GAMYB-binding protein 1 (GmGBP1), a member of the SNW/SKIP family of transcriptional coregulators in soybean, is induced by short days and interacts with transcription factor GAMYB (GmGAMYB) during photoperiod control of flowering time and maturity. In the present study, GmGBP1:GmGBP1 soybean showed the phenotypes of earlier maturity and higher plant height. Chromatin immunoprecipitation sequencing (ChIP-seq) assays of GmGBP1-binding sites and RNA sequencing (RNA-seq) of differentially expressed transcripts in GmGBP1:GmGBP1 further identified potential targets of GmGBP1, including small auxin-up RNA (GmSAUR). GmSAUR:GmSAUR soybean also showed earlier maturity and higher plant height. GmGBP1 interacted with GmGAMYB, bound to the promoter of GmSAUR and promoted the expression of FLOWER LOCUS T homologs 2a (GmFT2a) and FLOWERING LOCUS D LIKE 19 (GmFDL19). Flowering repressors such as GmFT4 were negatively regulated, resulting in earlier flowering and maturity. Furthermore, the interaction of GmGBP1 with GmGAMYB increased the gibberellin (GA) signal to promote height and hypocotyl elongation by activating GmSAUR and GmSAUR bound to the promoter of the GA-positive activating regulator gibberellic acid-stimulated Arabidopsis 32 (GmGASA32). These results suggested a photoperiod regulatory pathway in which the interaction of GmGBP1 with GmGAMYB directly activated GmSAUR to promote earlier maturity and plant height in soybean.

HEN1 SUPPRESSOR1 stabilizes polymerase IV RNAs via uridylation in Arabidopsis

3′-end uridylation of P4RNA directed by HEN1 SUPPRESSOR1 is critical for its stabilization.

SERRATE: a key factor in coordinated RNA processing in plants

The SERRATE (SE) protein is involved in the processing of RNA polymerase II (RNAPII) transcripts. It is associated with different complexes engaged in different aspects of plant RNA metabolism, including assemblies involved in transcription, splicing, polyadenylation, miRNA biogenesis, and RNA degradation. SE stability and interactome properties can be influenced by phosphorylation. SE exhibits an intriguing liquid-liquid phase separation property that may be important in the assembly of different RNA-processing bodies. Therefore, we propose that SE seems to participate in the coordination of different RNA-processing steps and can direct the fate of transcripts, targeting them for processing or degradation when they cannot be properly processed or are synthesized in excess.

Integrated analysis of smRNAome, transcriptome, and degradome data to decipher microRNAs regulating costunolide biosynthesis in Saussurea lappa

Saussurea lappa (S. lappa) has been known to synthesize medicinally important, costunolide. Due to its immense therapeutic importance, understanding of regulatory mechanism associated with its biosynthesis is crucial. The identification of genes and transcription factors (TFs) in S. lappa, created a clear picture of costunolide biosynthesis pathways. Further to understand the regulation of costunolide biosynthesis by miRNAs, an integrated study of transcriptome, miRNAs, and degradome was performed. Identified candidate miRNAs and associated feed-forward loops (FFLs) illustrates their regulatory role in secondary metabolite biosynthesis. Small RNA and degradome sequencing were performed for leaf and root tissues to determine miRNAs-targets pairs. A total of 711 and 525 such targets were obtained for novel and known miRNAs respectively. This data was used to generate costunolide-specific miRNA-TF-gene interactome to perform systematic analyses through graph theoretical approach. Interestingly, miR171c.1 and sla-miR121 were identified as key regulators to connect and co-regulate both mevalonate and sesquiterpenoid pathways to bio-synthesize costunolide. Tissue-specific FFLs were identified to be involved in costunolide biosynthesis which further suggests the evolutionary co-relation of root-specific networks in synthesis of secondary metabolites in addition to leaf-specific networks. This integrative approach allowed us to determine candidate miRNAs and associated tissue-specific motifs involved in the diversification of secondary metabolites. MiRNAs identified in present study can provide alternatives for bioengineering tool to enhance the synthesis of costunolide and other secondary metabolites in S. lappa.

Structural basis for RNA slicing by a plant Argonaute

Argonaute (AGO) proteins use small RNAs to recognize transcripts targeted for silencing in plants and animals. Many AGOs cleave target RNAs using an endoribonuclease activity termed 'slicing'. Slicing by DNA-guided prokaryotic AGOs has been studied in detail, but structural insights into RNA-guided slicing by eukaryotic AGOs are lacking. Here we present cryogenic electron microscopy structures of the Arabidopsis thaliana Argonaute10 (AtAgo10)-guide RNA complex with and without a target RNA representing a slicing substrate. The AtAgo10-guide-target complex adopts slicing-competent and slicing-incompetent conformations that are unlike known prokaryotic AGO structures. AtAgo10 slicing activity is licensed by docking target (t) nucleotides t9-t13 into a surface channel containing the AGO endoribonuclease active site. A β-hairpin in the L1 domain secures the t9-t13 segment and coordinates t9-t13 docking with extended guide-target pairing. Results show that prokaryotic and eukaryotic AGOs use distinct mechanisms for achieving target slicing and provide insights into small interfering RNA potency.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Structural insight into the activation of an Arabidopsis organellar C-to-U RNA editing enzyme by active site complementation

RNA-binding pentatricopeptide repeat (PPR) proteins catalyze hundreds of cytidine to uridine RNA editing events in plant organelles; these editing events are essential for proper gene expression. More than half of the PPR-type RNA editing factors, however, lack the DYW cytidine deaminase domain. Genetic analyses have suggested that their cytidine deaminase activity arises by association with a family of DYW1-like proteins that contain an N-terminally truncated DYW domain, but their molecular mechanism has been unclear. Here, we report the crystal structure of the Arabidopsis thaliana DYW1 deaminase domain at 1.8 Å resolution. DYW1 has a cytidine deaminase fold lacking the PG box. The internal insertion within the deaminase fold shows an α-helical fold instead of the β-finger reported for the gating domain of the A. thaliana ORGANELLE TRANSCRIPT PROCESSING 86. The substrate-binding pocket is incompletely formed and appears to be complemented in the complex by the E2 domain and the PG box of the interacting PPR protein. In vivo RNA editing assays corroborate the activation model for DYW1 deaminase. Our study demonstrates the common activation mechanism of the DYW1-like proteins by molecular complementation of the DYW domain and reconstitution of the substrate-binding pocket.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

Beyond transcription: compelling open questions in plant RNA biology

The study of RNAs has become one of the most influential research fields in contemporary biology and biomedicine. In the last few years, new sequencing technologies have produced an explosion of new and exciting discoveries in the field but have also given rise to many open questions. Defining these questions, together with old, long-standing gaps in our knowledge, is the spirit of this article. The breadth of topics within RNA biology research is vast, and every aspect of the biology of these molecules contains countless exciting open questions. Here, we asked 12 groups to discuss their most compelling question among some plant RNA biology topics. The following vignettes cover RNA alternative splicing; RNA dynamics; RNA translation; RNA structures; R-loops; epitranscriptomics; long non-coding RNAs; small RNA production and their functions in crops; small RNAs during gametogenesis and in cross-kingdom RNA interference; and RNA-directed DNA methylation. In each section, we will present the current state-of-the-art in plant RNA biology research before asking the questions that will surely motivate future discoveries in the field. We hope this article will spark a debate about the future perspective on RNA biology and provoke novel reflections in the reader.

A chaperonin containing T-complex polypeptide-1 facilitates the formation of the PbWoxT1-PbPTB3 ribonucleoprotein complex for long-distance RNA trafficking in Pyrus betulaefolia

Numerous plant endogenous mRNAs move via phloem and thus affect the growth and development of long-distant organs. mRNAs are transported with RNA-binding proteins forming a ribonucleoprotein complex. However, it remains elusive how such RNP complex assembles and facilitates mRNA trafficking. Protease digestion and RNA immunoprecipitation were used to investigate the RNP assembly function of the complete Chaperonin Containing T-complex Polypeptide-1. In situ hybridization, hairy root transformation, microprojectile bombardment, and grafting experiments demonstrate the role of CCT complex in the transport of a PbWoxT1-PbPTB3 RNP complex in Pyrus betulaefolia. PbCCT5 silenced caused defective movement of GFP-PbPTB3 and GFP-PbWoxT1 from hairy roots to new leaves via the phloem. PbCCT5 is shown to interact with PbPTB3. PbCCT complex enhanced PbPTB3 stabilization and permitted assembly of PbWoxT1 and PbPTB3 into an RNP complex. Furthermore, silencing of individual CCT subunits inhibited the intercellular movement of GFP-PbPTB3 and long-distance trafficking of PbWoxT1 and PbPTB3 in grafted plants. Taken together, the CCT complex assembles PbPTB3 and PbWoxT1 into an RNP complex in the phloem in order to facilitate the long-distance trafficking of PbWoxT1 in P. betulaefolia. This study therefore provides important insights into the mechanism of RNP complex formation and transport.

RETINOBLASTOMA-RELATED interactions with key factors of the RNA-directed DNA methylation (RdDM) pathway and its influence on root development

MAIN CONCLUSION

Our study presents evidence for a novel mechanism for RBR function in transcriptional gene silencing by interacting with key players of the RdDM pathway in Arabidopsis and several plant clades. Transposable elements and other repetitive elements are silenced by the RNA-directed DNA methylation pathway (RdDM). In RdDM, POLIV-derived transcripts are converted into double-stranded RNA (dsRNA) by the activity of RDR2 and subsequently processed into 24 nucleotide short interfering RNAs (24-nt siRNAs) by DCL3. 24-nt siRNAs serve as guides to direct AGO4-siRNA complexes to chromatin-bound POLV-derived transcripts generated from the template/target DNA. The interaction between POLV, AGO4, DMS3, DRD1, RDM1 and DRM2 promotes DRM2-mediated de novo DNA methylation. The Arabidopsis Retinoblastoma protein homolog (RBR) is a master regulator of the cell cycle, stem cell maintenance, and development. We in silico predicted and explored experimentally the protein-protein interactions (PPIs) between RBR and members of the RdDM pathway. We found that the largest subunits of POLIV and POLV (NRPD1 and NRPE1), the shared second largest subunit of POLIV and POLV (NRPD/E2), RDR1, RDR2, DCL3, DRM2, and SUVR2 contain canonical and non-canonical RBR binding motifs and several of them are conserved since algae and bryophytes. We validated experimentally PPIs between Arabidopsis RBR and several of the RdDM pathway proteins. Moreover, seedlings from loss-of-function mutants in RdDM and RBR show similar phenotypes in the root apical meristem. We show that RdDM and SUVR2 targets are up-regulated in the 35S:AmiGO-RBR background.

C-to-U and U-to-C: RNA editing in plant organelles and beyond

The genomes in the two energy-converting organelles of plant cells, chloroplasts and mitochondria, contain numerous 'errors' that are corrected at the level of RNA transcript copies. The genes encoded in the two endosymbiotic organelles would not function properly if their transcripts were not altered by site-specific cytidine-to-uridine (C-to-U) exchanges and by additional reverse U-to-C exchanges in hornworts, lycophytes, and ferns. These peculiar processes of plant RNA editing, re-establishing genetic information that could alternatively be present at the organelle genome level, has spurred much research over >30 years. Lately new studies have revealed numerous interesting insights, notably on the biochemical machinery identifying specific pyrimidine nucleobases for conversion from C to U and vice versa. Here, I will summarize prominent research findings that lately have contributed to our better understanding of these phenomena introducing an added layer of information processing in plant cells. Some of this recent progress is based on the successful functional expression of plant RNA editing factors in bacteria and mammalian cells. These research approaches have recapitulated natural processes of horizontal gene transfer through which some protist lineages seem to have acquired plant RNA editing factors and adapted them functionally for their own purposes.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

miRador: a fast and precise tool for the prediction of plant miRNAs

Plant microRNAs (miRNAs) are short, noncoding RNA molecules that restrict gene expression via posttranscriptional regulation and function in several essential pathways, including development, growth, and stress responses. Accurately identifying miRNAs in populations of small RNA sequencing libraries is a computationally intensive process that has resulted in the misidentification of inaccurately annotated miRNA sequences. In recent years, criteria for miRNA annotation have been refined with the aim to reduce these misannotations. Here, we describe miRador, a miRNA identification tool that utilizes the most up-to-date, community-established criteria for accurate identification of miRNAs in plants. We combined target prediction and Parallel Analysis of RNA Ends (PARE) data to assess the precision of the miRNAs identified by miRador. We compared miRador to other commonly used miRNA prediction tools and found that miRador is at least as precise as other prediction tools while being substantially faster than other tools. miRador should be broadly useful for the plant community to identify and annotate miRNAs in plant genomes.

Tweaking the Small Non-Coding RNAs to Improve Desirable Traits in Plant

Plant transcriptome contains an enormous amount of non-coding RNAs (ncRNAs) that do not code for proteins but take part in regulating gene expression. Since their discovery in the early 1990s, much research has been conducted to elucidate their function in the gene regulatory network and their involvement in plants' response to biotic/abiotic stresses. Typically, 20-30 nucleotide-long small ncRNAs are a potential target for plant molecular breeders because of their agricultural importance. This review summarizes the current understanding of three major classes of small ncRNAs: short-interfering RNAs (siRNAs), microRNA (miRNA), and transacting siRNAs (tasiRNAs). Furthermore, their biogenesis, mode of action, and how they have been utilized to improve crop productivity and disease resistance are discussed here.

Novel insights on genes and pathways involved in Pinus elliottii response to resinosis

Pinus elliottii, an important coniferous timber species, has recently become one of the most popular sources of resin in China. Resinosis is a common disease that may negatively affect pine tree growth and production. In this study, we used single-molecule real-time sequencing and Illumina RNA sequencing to generate an accurate transcriptome for P. elliottii. The transcriptome included 90,026 transcripts, 5160 long non-coding RNAs and 7710 transcription factors. We then analyzed RNA-sequencing, small RNA-sequencing and degradome data to identify genes, miRNAs and key miRNA-target pairs involved in response to resinosis in P. elliottii. We identified 1305 genes and 1151 miRNAs exhibiting significant differential expression in response to resinosis. According to the degradome sequencing analysis, 318 differentially expressed transcripts were targets of 14 differentially expressed miRNAs. Our study has provided resources for further functional characterization of genes and miRNAs involved in resinosis in P. elliottii, which should aid the future disease-resistance breeding of this species.

Enzymatic reactions of AGO4 in RNA-directed DNA methylation: siRNA duplex loading, passenger strand elimination, target RNA slicing, and sliced target retention

RNA-directed DNA methylation in plants is guided by 24-nt siRNAs generated in parallel with 23-nt RNAs of unknown function. We show that 23-nt RNAs function as passenger strands during 24-nt siRNA incorporation into AGO4. The 23-nt RNAs are then sliced into 11- and 12-nt fragments, with 12-nt fragments remaining associated with AGO4. Slicing recapitulated with recombinant AGO4 and synthetic RNAs reveals that siRNAs of 21-24 nt, with any 5'-terminal nucleotide, can guide slicing, with sliced RNAs then retained by AGO4. In vivo, RdDM target locus RNAs that copurify with AGO4 also display a sequence signature of slicing. Comparing plants expressing slicing-competent versus slicing-defective AGO4 shows that slicing elevates cytosine methylation levels at virtually all RdDM loci. We propose that siRNA passenger strand elimination and AGO4 tethering to sliced target RNAs are distinct modes by which AGO4 slicing enhances RNA-directed DNA methylation.

Chloroplast C-to-U RNA editing in vascular plants is adaptive due to its restorative effect: testing the restorative hypothesis

The adaptiveness of nonsynonymous RNA editing (recoding) could be conferred by the flexibility of the temporal-spatially controllable proteomic diversity, or by its restorative effect which fixes unfavorable genomic mutations at the RNA level. These two complementary hypotheses, namely, the diversifying hypothesis and the restorative hypothesis, have distinct predictions on the landscape of RNA editing sites. We collected the chloroplast C-to-U RNA editomes of 21 vascular plants (11 angiosperms, four gymnosperms, and six ferns) from a previous study, aiming to testify whether the plant editomes typically conform to the restorative hypothesis. All predictions made by the restorative hypothesis are verified: (i) nonsynonymous editing sites are more frequent and have higher editing levels than synonymous sites; (ii) nonsynonymous editing levels are extremely high and show weak tissue-specificity in plants; (iii) on the inferred genomic sites with recent T-to-C mutations, nonsynonymous sites but not synonymous sites are compensated by C-to-U RNA editing. In conclusion, nonsynonymous C-to-U RNA editing in plants is adaptive due to its restorative effects. The recoding levels are high and are constantly required across the whole plant so that the recoding events could perfectly mimic DNA mutations. The evolutionary significance of plant RNA editing is systematically demonstrated at the genome-wide level.

Changes in the m6A RNA methylome accompany the promotion of soybean root growth by rhizobia under cadmium stress

Cadmium (Cd) is the most widely distributed heavy metal pollutant in soil and has significant negative effects on crop yields and human health. Rhizobia can enhance soybean growth in the presence of heavy metals, and the legume-rhizobia symbiosis has been used to promote heavy-metal phytoremediation, but much remains to be learned about the molecular networks that underlie these effects. Here, we demonstrated that soybean root growth was strongly suppressed after seven days of Cd exposure but that the presence of rhizobia largely eliminated this effect, even prior to nodule development. Moreover, rhizobia did not appear to promote root growth by limiting plant Cd uptake: seedlings with and without rhizobia had similar root Cd concentrations. Previous studies have demonstrated a role for m6A RNA methylation in the response of rice and barley to Cd stress. We therefore performed transcriptome-wide m6A methylation profiling to investigate changes in the soybean RNA methylome in response to Cd with and without rhizobia. Here, we provide some of the first data on transcriptome-wide m6a RNA methylation patterns in soybean; m6A modifications were concentrated at the 3' UTR of transcripts and showed a positive relationship with transcript abundance. Transcriptome-wide m6A RNA methylation peaks increased in the presence of Cd, and the integration of m6A methylome and transcriptome results enabled us to identify 154 genes whose transcripts were both differentially methylated and differentially expressed in response to Cd stress. Annotation results suggested that these genes were associated with Ca²⁺ homeostasis, ROS pathways, polyamine metabolism, MAPK signaling, hormones, and biotic stress responses. There were 176 differentially methylated and expressed transcripts under Cd stress in the presence of rhizobia. In contrast to the Cd-only gene set, they were also enriched in genes related to auxin, jasmonic acid, and brassinosteroids, as well as abiotic stress tolerance. They contained fewer genes related to Ca²⁺ homeostasis and also included candidates with known functions in the legume-rhizobia symbiosis. These findings offer new insights into how rhizobia promote soybean root growth under Cd stress; they provide candidate genes for research on plant heavy metal responses and for the use of legumes in phytoremediation.

Coexpression network and trans-activation analyses of maize reproductive phasiRNA loci

The anther-enriched phased, small interfering RNAs (phasiRNAs) play vital roles in sustaining male fertility in grass species. Their long non-coding precursors are synthesized by RNA polymerase II and are likely regulated by transcription factors (TFs). A few putative transcriptional regulators of the 21- or 24-nucleotide phasiRNA loci (referred to as 21- or 24-PHAS loci) have been identified in maize (Zea mays), but whether any of the individual TFs or TF combinations suffice to activate any PHAS locus is unclear. Here, we identified the temporal gene coexpression networks (modules) associated with maize anther development, including two modules highly enriched for the 21- or 24-PHAS loci. Comparisons of these coexpression modules and gene sets dysregulated in several reported male sterile TF mutants provided insights into TF timing with regard to phasiRNA biogenesis, including antagonistic roles for OUTER CELL LAYER4 and MALE STERILE23. Trans-activation assays in maize protoplasts of individual TFs using bulk-protoplast RNA-sequencing showed that two of the TFs coexpressed with 21-PHAS loci could activate several 21-nucleotide phasiRNA pathway genes but not transcription of 21-PHAS loci. Screens for combinatorial activities of these TFs and, separately, the recently reported putative transcriptional regulators of 24-PHAS loci using single-cell (protoplast) RNA-sequencing, did not detect reproducible activation of either 21-PHAS or 24-PHAS loci. Collectively, our results suggest that the endogenous transcriptional machineries and/or chromatin states in the anthers are necessary to activate reproductive PHAS loci.

miR319 and its target TCP4 involved in plant architecture regulation in Brassica napus

Plant architecture is a collection of genetically controlled crop productivity and adaptation. MicroRNAs (miRNAs) have been proved to function in various biological processes, but little is known about how miRNA regulates plant architecture in rapeseed (Brassica napus L.). In this study, four small RNA libraries and two degradome libraries from shoot apex of normal and rod-like plants were sequenced. A total of 639 miRNA precursors and 16 differentially expressed miRNAs were identified in this study. In addition, 322 targets were identified through degradome sequencing. Among them, 14 targets were further validated via RNA ligase-mediated 5' rapid amplification of cDNA ends. Transgenic approach showed that increased TCP4 activity in Arabidopsis resulted in premature onset of maturation and reduced plant size along with early flowering and shortened flowering time. miR319-OE lines in Brassica napus exhibited serrated leaves and abnormal development of shoot apical meristem (SAM), which led to the deformed growth of stem and reduced plant height. In conclusion, our study lays the foundation for elucidating miRNA regulate plant architecture and provides new insight into the miR319/TCP4 module regulates plant architecture in rapeseed.