BACKGROUND: β-Aminobutyric acid (BABA) has been successfully used to prime stress resistance in numerous plant species; however, its effectiveness in forest trees has been poorly explored thus far. This study aimed to investigate the influence of BABA on morphological, physiological, and epigenetic parameters in field elms under various growth conditions. Epigenetic changes were assessed in both DNA and RNA through the use of reversed-phase ultra-performance liquid chromatography (UPLC) coupled with sensitive mass spectrometry.
RESULTS: The presented results confirm the influence of BABA on the development, physiology, and stress tolerance in field elms. However, the most important findings are related to the broad epigenetic changes promoted by this amino acid, which involve both DNA and RNA. Our findings confirm, for the first time, that BABA influences not only well-known epigenetic markers in plants, such as 5-methylcytosine, but also several other non-canonical nucleobases, such as 5-hydroxymethyluracil, 5-formylcytosine, 5-hydroxymethylcytosine, N6-methyladenine, uracil (in DNA) and thymine (in RNA). The significant effect on the levels of N6-methyladenine, the main bacterial epigenetic marker, is particularly noteworthy. In this case, the question arises as to whether this effect is due to epigenetic changes in the microbiome, the plant genome, or both.
CONCLUSIONS: The plant phenotype is the result of complex interactions between the plant\'s DNA, the microbiome, and the environment. We propose that different types of epigenetic changes in the plant and microbiome may play important roles in the largely unknown memory process that enables plants to adapt faster to changing environmental conditions.

CONCLUSIONS: Our data link the miR165/166- and miR160-mediated regulatory modules to ROS and seed formation. Trade-offs of seed size, weight, and number probably require control of the expression of miR165/166 by miR160, modulation of ROS metabolism by miR165/166, and miR160 abundance by ROS-induced oxidative modifications The cycle of plant life and its yield productivity depends fundamentally on the establishment of the trade-offs of seed size, weight, and number. For annual plants, seed number should simply be a positive function of vegetative biomass and a negative function of seed size and/or weight. However, extensive natural variation within species is observed for these traits, for which an optimal solution is environmentally dependent. Understanding the miRNA-mediated post-transcriptional regulation of gene expression determining seed phenotype and number is crucial from both an evolutionary and applied perspective. Although extensive research has concentrated on the individual roles of miRNAs in plant life, fewer studies have centred on their functional interactions, hence this study aimed to examine whether the module of miR165/miR166 and/or miR160 interactions is involved in forming Arabidopsis thaliana seeds, and/or has an impact on their features. Considering that reactive oxygen species (ROS) are among key players in seed-related processes, it was also intriguing to verify if the mechanism of action of these miRNAs is associated with the ROS pathway. The plant material used in this study consisted of flower buds, green siliques, and freshly harvested seeds, of wild type (WT), and STTM165/166 and STTM160 × 165/166 mutants of A. thaliana plants which are powerful tools for functional analysis of miRNAs in plants. The novel results obtained during physiological phenotyping together with two-tailed qRT-PCR analysis of mature miR165, miR166, miR160, and spectrofluorimetric measurement of apoplastic hydrogen peroxide (H2O2) for the first time revealed that interaction between miR165/miR166 and miR160 may regulate seed size, weight and number in ROS-dependent manner.

BACKGROUND: In recent years, covalent modifications on RNA nucleotides have emerged as pivotal moieties influencing the structure, function, and regulatory processes of RNA Polymerase II transcripts such as mRNAs and lncRNAs. However, our understanding of their biological roles and whether these roles are conserved across eukaryotes remains limited.
RESULTS: In this study, we leveraged standard polyadenylation-enriched RNA-sequencing data to identify and characterize RNA modifications that introduce base-pairing errors into cDNA reads. Our investigation incorporated data from three Poaceae (Zea mays, Sorghum bicolor, and Setaria italica), as well as publicly available data from a range of stress and genetic contexts in Sorghum and Arabidopsis thaliana. We uncovered a strong enrichment of RNA covalent modifications (RCMs) deposited on a conserved core set of nuclear mRNAs involved in photosynthesis and translation across these species. However, the cohort of modified transcripts changed based on environmental context and developmental program, a pattern that was also conserved across flowering plants. We determined that RCMs can partly explain accession-level differences in drought tolerance in Sorghum, with stress-associated genes receiving a higher level of RCMs in a drought tolerant accession. To address function, we determined that RCMs are significantly enriched near exon junctions within coding regions, suggesting an association with splicing. Intriguingly, we found that these base-pair disrupting RCMs are associated with stable mRNAs, are highly correlated with protein abundance, and thus likely associated with facilitating translation.
CONCLUSIONS: Our data point to a conserved role for RCMs in mRNA stability and translation across the flowering plant lineage.

Salinity is a common abiotic stress that limits crop productivity. Although there is a wealth of evidence suggesting that miRNA and lncRNA play important roles in the response to salinity in rice seedlings and reproductive stages, the mechanism by which competing endogenous RNAs (ceRNAs) influence salt tolerance and yield in rice has been rarely reported. In this study, we conducted full whole-transcriptome sequencing of rice panicles during the reproductive period to clarify the role of ceRNAs in the salt stress response and yield. A total of 214 lncRNAs, 79 miRNAs, and 584 mRNAs were identified as differentially expressed RNAs under salt stress. Functional analysis indicates that they play important roles in GO terms such as response to stress, biosynthesis processes, abiotic stimuli, endogenous stimulus, and response to stimulus, as well as in KEGG pathways such as secondary metabolite biosynthesis, carotenoid biosynthesis, metabolic pathways, and phenylpropanoid biosynthesis. A ceRNA network comprising 95 lncRNA-miRNA-mRNA triplets was constructed. Two lncRNAs, MSTRG.51634.2 and MSTRG.48576.1, were predicted to bind to osa-miR172d-5p to regulate the expression of OsMYB2 and OsMADS63, which have been reported to affect salt tolerance and yield, respectively. Three lncRNAs, MSTRG.30876.1, MSTRG.44567.1, and MSTRG.49308.1, may bind to osa-miR5487 to further regulate the expression of a stress protein (LOC_Os07g48460) and an aquaporin protein (LOC_Os02g51110) to regulate the salt stress response. This study is helpful for understanding the underlying molecular mechanisms of ceRNA that drive the response of rice to salt stress and provide new genetic resources for salt-resistant rice breeding.

Plants employ distinct mechanisms to respond to environmental changes. Modification of mRNA by N 6-methyladenosine (m6A), known to affect the fate of mRNA, may be one such mechanism to reprogram mRNA processing and translatability upon stress. However, it is difficult to distinguish a direct role from a pleiotropic effect for this modification due to its prevalence in RNA. Through characterization of the transient knockdown-mutants of m6A writer components and mutants of specific m6A readers, we demonstrate the essential role that m6A plays in basal resistance and pattern-triggered immunity (PTI). A global m6A profiling of mock and PTI-induced Arabidopsis plants as well as formaldehyde fixation and cross-linking immunoprecipitation-sequencing of the m6A reader, EVOLUTIONARILY CONSERVED C-TERMINAL REGION2 (ECT2) showed that while dynamic changes in m6A modification and binding by ECT2 were detected upon PTI induction, most of the m6A sites and their association with ECT2 remained static. Interestingly, RNA degradation assay identified a dual role of m6A in stabilizing the overall transcriptome while facilitating rapid turnover of immune-induced mRNAs during PTI. Moreover, polysome profiling showed that m6A enhances immune-associated translation by binding to the ECT2/3/4 readers. We propose that m6A plays a positive role in plant immunity by destabilizing defense mRNAs while enhancing their translation efficiency to create a transient surge in the production of defense proteins.

BACKGROUND: MicroRNAs (miRNAs) are small noncoding RNAs that play important post-transcriptional regulatory roles in animals and plants. Despite the importance of plant miRNAs, the inherent complexity of miRNA biogenesis in plants hampers the application of standard miRNA prediction tools, which are often optimized for animal sequences. Therefore, computational approaches to predict putative miRNAs (merely) from genomic sequences, regardless of their expression levels or tissue specificity, are of great interest.
RESULTS: Here, we present AmiR-P3, a novel ab initio plant miRNA prediction pipeline that leverages the strengths of various utilities for its key computational steps. Users can readily adjust the prediction criteria based on the state-of-the-art biological knowledge of plant miRNA properties. The pipeline starts with finding the potential homologs of the known plant miRNAs in the input sequence(s) and ensures that they do not overlap with protein-coding regions. Then, by computing the secondary structure of the presumed RNA sequence based on the minimum free energy, a deep learning classification model is employed to predict potential pre-miRNA structures. Finally, a set of criteria is used to select the most likely miRNAs from the set of predicted miRNAs. We show that our method yields acceptable predictions in a variety of plant species.
CONCLUSIONS: AmiR-P3 does not (necessarily) require sequencing reads and/or assembled reference genomes, enabling it to identify conserved and novel putative miRNAs from any genomic or transcriptomic sequence. Therefore, AmiR-P3 is suitable for miRNA prediction even in less-studied plants, as it does not require any prior knowledge of the miRNA repertoire of the organism. AmiR-P3 is provided as a docker container, which is a portable and self-contained software package that can be readily installed and run on any platform and is freely available for non-commercial use from: https://hub.docker.com/r/micrornaproject/amir-p3.

Lettuce is one of the most widely cultivated and consumed dicotyledonous vegetables globally. Despite the availability of its reference genome sequence, lettuce gene annotation remains incomplete, impeding comprehensive research and the broad application of genomic resources. Long-read RNA isoform sequencing (Iso-Seq) offers substantial advantages for analyzing RNA alternative splicing and aiding gene annotation, yet it faces throughput limitations. We present the HIT-ISOseq method tailored for bulk sample analysis, significantly enhancing RNA sequencing throughput on the PacBio platform by concatenating cDNA. Here we show, HIT-ISOseq generates 3-4 cDNA molecules per CCS read in lettuce, yielding 15.7 million long reads per PacBio Sequel II SMRT Cell 8 M. We validate its effectiveness in analyzing six lettuce tissue samples, including roots, stems, and leaves, revealing tissue-specific gene expression patterns and RNA isoforms. Leveraging diverse tissue long-read RNA sequencing, we refine the transcript annotation of the lettuce reference genome, expanding its GO and KEGG annotation repertoire. Collectively, this study serves as a foundational reference for genome annotation and the analysis of multi-sample isoform expression, utilizing high-throughput long-read transcriptome sequencing.

The development of next-generation sequencing technology has led to a burst of data in a single assay. Management of a large dataset requires high demands on bioinformatic skills and computing resources. Here we present two popular pipelines for RNA-seq data analysis, using open-source software tools HISAT-StringTie-Ballgown and TopHat-Cufflinks. To meet the need of plant scientist, we describe in detail how to perform such comprehensive analysis beginning with raw RNA-seq reads and available reference genome. It allows biologists to align short reads to a reference genome, measure the transcript abundance, and analyze gene differential expression under two or more conditions. We also discuss other RNA-seq tools that are comparable or alternative to this protocol.

In the realm of plant biology, small RNAs (sRNAs) are imperative in the orchestration of gene expression, playing pivotal roles across a spectrum of developmental sequences and responses to environmental stressors. The biosynthetic cascade of sRNAs is characterized by an elaborate network of enzymatic pathways that meticulously process double-stranded RNA (dsRNA) precursors into sRNA molecules, typically 20 to 30 nucleotides in length. These sRNAs, chiefly microRNAs (miRNAs) and small interfering RNAs (siRNAs), are integral in guiding the RNA-induced silencing complex (RISC) to selectively target messenger RNAs (mRNAs) for post-transcriptional modulation. This regulation is achieved either through the targeted cleavage or the suppression of translational efficiency of the mRNAs. In plant development, sRNAs are integral to the modulation of key pathways that govern growth patterns, organ differentiation, and developmental timing. The biogenesis of sRNA itself is a fine-tuned process, beginning with transcription and proceeding through a series of processing steps involving Dicer-like enzymes and RNA-binding proteins. Recent advances in the field have illuminated the complex processes underlying the generation and function of small RNAs (sRNAs), including the identification of new sRNA categories and the clarification of their involvement in the intercommunication among diverse regulatory pathways. This review endeavors to evaluate the contemporary comprehension of sRNA biosynthesis and to underscore the pivotal role these molecules play in directing the intricate performance of plant developmental processes.

Circular RNA (circRNA) is a type of non-coding RNA with multiple biological functions. Whole circRNA genomes in plants have been identified, and circRNAs have been demonstrated to be widely present and highly expressed in various plant tissues and organs. CircRNAs are highly stable and conserved in plants, and exhibit tissue specificity and developmental stage specificity. CircRNAs often interact with other biomolecules, such as miRNAs and proteins, thereby regulating gene expression, interfering with gene function, and affecting plant growth and development or response to environmental stress. CircRNAs are less studied in plants than in animals, and their regulatory mechanisms of biogenesis and molecular functions are not fully understood. A variety of circRNAs in plants are involved in regulating growth and development and responding to environmental stress. This review focuses on the biogenesis and regulatory mechanisms of circRNAs, as well as their biological functions during growth, development, and stress responses in plants, including a discussion of plant circRNA research prospects. Understanding the generation and regulatory mechanisms of circRNAs is a challenging but important topic in the field of circRNAs in plants, as it can provide insights into plant life activities and their response mechanisms to biotic or abiotic stresses as well as new strategies for plant molecular breeding and pest control.