ß-Aminobutyric acid promotes stress tolerance, physiological adjustments, as well as broad epigenetic changes at DNA and RNA nucleobases in field elms (Ulmus minor).

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.

AmiR-P3: An AI-based microRNA prediction pipeline in plants.

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. CONCLUSION: 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.

The interaction between miR165/166 and miR160 regulates Arabidopsis thaliana seed size, weight, and number in a ROS-dependent manner.

MAIN CONCLUSION: 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.

Identification and Characterization of miRNAs and lncRNAs Associated with Salinity Stress in Rice Panicles.

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.

Comparative analyses suggest a link between mRNA splicing, stability, and RNA covalent modifications in flowering plants.

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.

RiceMetaSys: Drought-miR, a one-stop solution for drought responsive miRNAs-mRNA module in rice.

MicroRNAs are key players involved in stress responses in plants and reports are available on the role of miRNAs in drought stress response in rice. This work reports the development of a database, RiceMetaSys: Drought-miR, based on the meta-analysis of publicly available sRNA datasets. From 28 drought stress-specific sRNA datasets, we identified 216 drought-responsive miRNAs (DRMs). The major features of the database include genotype-, tissue- and miRNA ID-specific search options and comparison of genotypes to identify common miRNAs. Co-localization of the DRMs with the known quantitative trait loci (QTLs), i.e., meta-QTL regions governing drought tolerance in rice pertaining to different drought adaptive traits, narrowed down this to 37 promising DRMs. To identify the high confidence target genes of DRMs under drought stress, degradome datasets and web resource on drought-responsive genes (RiceMetaSys: DRG) were used. Out of the 216 unique DRMs, only 193 had targets with high stringent parameters. Out of the 1081 target genes identified by Degradome datasets, 730 showed differential expression under drought stress in at least one accession. To retrieve complete information on the target genes, the database has been linked with RiceMetaSys: DRG. Further, we updated the RiceMetaSys: DRGv1 developed earlier with the addition of DRGs identified from RNA-seq datasets from five rice genotypes. We also identified 759 putative novel miRNAs and their target genes employing stringent criteria. Novel miRNA search has all the search options of known miRNAs and additionally, it gives information on their in silico validation features. Simple sequence repeat markers for both the miRNAs and their target genes have also been designed and made available in the database. Network analysis of the target genes identified 60 hub genes which primarily act through abscisic acid pathway and jasmonic acid pathway. Co-localization of the hub genes with the meta-QTL regions governing drought tolerance narrowed down this to 16 most promising DRGs. Database URL: http://14.139.229.201/RiceMetaSys_miRNA Updated database of RiceMetaSys URL: http://14.139.229.201/RiceMetaSysA/Drought/.

m6A modification plays an integral role in mRNA stability and translation during pattern-triggered immunity.

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.

Plant Transcriptome Analysis with HISAT-StringTie-Ballgown and TopHat-Cufflinks Pipelines.

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.

Long-read RNA-Seq for the discovery of long noncoding and antisense RNAs in plant organelles.

Plant organelle transcription has been studied for decades. As techniques advanced, so did the fields of mitochondrial and plastid transcriptomics. The current view is that organelle genomes are pervasively transcribed, irrespective of their size, content, structure, and taxonomic origin. However, little is known about the nature of organelle noncoding transcriptomes, including pervasively transcribed noncoding RNAs (ncRNAs). Next-generation sequencing data have uncovered small ncRNAs in the organelles of plants and other organisms, but long ncRNAs remain poorly understood. Here, we argue that publicly available third-generation long-read RNA sequencing data from plants can provide a fine-tuned picture of long ncRNAs within organelles. Indeed, given their bloated architectures, plant mitochondrial genomes are well suited for studying pervasive transcription of ncRNAs. Ultimately, we hope to showcase this new avenue of plant research while also underlining the limitations of the proposed approach.

Regulatory Peptide Encoded by the Primary Transcript of miR396a Influences Gene Expression and Root Development in Solanum lycopersicum.

MicroRNAs (miRNAs) are the processing products of primary miRNAs (pri-miRNAs) that regulate the expression of target genes. Recent studies have demonstrated that some pri-miRNAs can encode small peptides (miPEPs) that perform significant biological functions. The function of miPEPs in tomatoes, an important model horticultural crop, remains to be investigated. Here, we characterized the primary sequence of tomato miR396a using 5' RACE and confirmed the presence of miPEP396a in tomato by verifying the translational activity of the start codon. It primarily resides in the nucleus to exert its function and additionally regulates the expression of pri-miR396a, miR396a, and its target genes. Transcriptomic and metabolomic analyses showed that in vitro synthesis of miPEP396a significantly increased the expression of genes related to phenylpropanoid biosynthesis and hormones in tomato. Meanwhile, our in vitro application of miPEP396a in tomato significantly inhibited the elongation of tomato primary roots. In conclusion, our results indicate that miPEP396a regulates root growth in tomato by specifically promoting miR396a expression, provide insight into the function of miPEPs in tomato and potential applications.

MicroRNA166: Old Players and New Insights into Crop Agronomic Traits Improvement.

MicroRNA (miRNA), a type of non-coding RNA, is crucial for controlling gene expression. Among the various miRNA families, miR166 stands out as a highly conserved group found in both model and crop plants. It plays a key role in regulating a wide range of developmental and environmental responses. In this review, we explore the diverse sequences of MIR166s in major crops and discuss the important regulatory functions of miR166 in plant growth and stress responses. Additionally, we summarize how miR166 interacts with other miRNAs and highlight the potential for enhancing agronomic traits by manipulating the expression of miR166 and its targeted HD-ZIP III genes.

Advances in CircRNAs in the Past Decade: Review of CircRNAs Biogenesis, Regulatory Mechanisms, and Functions in Plants.

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.

The Biosynthesis Process of Small RNA and Its Pivotal Roles in Plant Development.

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.

Elucidating the regulatory role of long non-coding RNAs in drought stress response during seed germination in leaf mustard.

Leaf mustard (Brassica juncea L. Czern & Coss), an important vegetable crop, experiences pronounced adversity due to seasonal drought stress, particularly at the seed germination stage. Although there is partial comprehension of drought-responsive genes, the role of long non-coding RNAs (lncRNAs) in adjusting mustard's drought stress response is largely unexplored. In this study, we showed that the drought-tolerant cultivar 'Weiliang' manifested a markedly lower base water potential (-1.073 MPa vs -0.437 MPa) and higher germination percentage (41.2% vs 0%) than the drought-susceptible cultivar 'Shuidong' under drought conditions. High throughput RNA sequencing techniques revealed a significant repertoire of lncRNAs from both cultivars during germination under drought stress, resulting in the identification of 2,087 differentially expressed lncRNAs (DELs) and their correspondingly linked 12,433 target genes. It was noted that 84 genes targeted by DEL exhibited enrichment in the photosynthesis pathway. Gene network construction showed that MSTRG.150397, a regulatory lncRNA, was inferred to potentially modulate key photosynthetic genes (Psb27, PetC, PetH, and PsbW), whilst MSTRG.107159 was indicated as an inhibitory regulator of six drought-responsive PIP genes. Further, weighted gene co-expression network analysis (WGCNA) corroborated the involvement of light intensity and stress response genes targeted by the identified DELs. The precision and regulatory impact of lncRNA were verified through qPCR. This study extends our knowledge of the regulatory mechanisms governing drought stress responses in mustard, which will help strategies to augment drought tolerance in this crop.

LncRNA cis- and trans-regulation provides new insight into drought stress responses in wild barley.

Drought is one of the most common abiotic stresses that affect barley productivity. Long noncoding RNA (lncRNA) has been reported to be widely involved in abiotic stress, however, its function in the drought stress response in wild barley remains unclear. In this study, RNA sequencing was performed to identify differentially expressed lncRNAs (DElncRNA) among two wild barley and two cultivated barley genotypes. Then, the cis-regulatory networks were according to the chromosome position and the expression level correction. The GO annotation indicates that these cis-target genes are mainly involved in "ion transport transporter activity" and "metal ion transport transporter activity". Through weighted gene co-expression network analysis (WGCNA), 10 drought-related modules were identified to contract trans-regulatory networks. The KEGG annotation demonstrated that these trans-target genes were enriched for photosynthetic physiology, brassinosteroid biosynthesis, and flavonoid metabolism. In addition, we constructed the lncRNA-mediated ceRNA regulatory network by predicting the microRNA response elements (MREs). Furthermore, the expressions of lncRNAs were verified by RT-qPCR. Functional verification of a candidate lncRNA, MSTRG.32128, demonstrated its positive role in drought response and root growth and development regulation. Hormone content analysis provided insights into the regulatory mechanisms of MSTRG.32128 in root development, revealing its involvement in auxin and ethylene signal transduction pathways. These findings advance our understanding of lncRNA-mediated regulatory mechanisms in barley under drought stress. Our results will provide new insights into the functions of lncRNAs in barley responding to drought stress.

A Novel Method for In Situ RT-PCR Based on Capsules from Centrifuge Tubes, Ideal for Transcripts Detection in Plant Tissues.

In situ RT-PCR presents advantages over other expression analysis methods due to its rapid processing and low-cost equipment. However, this technique is not without its challenges. A protocol based on a capsule made from centrifuge tubes that offers advantages over slides is presented. This capsule protects histological sections from drying out, and its easy assembly reduces time pauses between incubations. In addition, the container size where the sample is deposited allows the addition and withdrawal of the different solutions. The capsule does not need previous sealing after each incubation, and, above all, it is a low-cost and accessible material. A guideline for tissue sectioning using a cryostat that offers advantages over other sectioning methods is also described.

High-throughput single-molecule long-read RNA sequencing analysis of tissue-specific genes and isoforms in lettuce (Lactuca sativa L.).

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.

lncRNAs and epigenetics regulate plant's resilience against biotic stresses.

With the advent of transcriptomic techniques involving single-stranded RNA sequencing and chromatin isolation by RNA purification-based sequencing, transcriptomic studies of coding and non-coding RNAs have been executed efficiently. These studies acknowledged the role of non-coding RNAs in modulating gene expression. Long non-coding RNAs (lncRNAs) are a kind of non-coding RNAs having lengths of >200 nucleotides, playing numerous roles in plant developmental processes such as photomorphogenesis, epigenetic changes, reproductive tissue development, and in regulating biotic and abiotic stresses. Epigenetic changes further control gene expression by changing their state to "ON-OFF" and also regulate stress memory and its transgenerational inheritance. With well-established regulatory mechanisms, they act as guides, scaffolds, signals, and decoys to modulate gene expression. They act as a major operator of post-transcriptional modifications such as histone and epigenetic modifications, and DNA methylations. The review elaborates on the roles of lncRNAs in plant immunity and also discusses how epigenetic markers alter gene expression in response to pest/pathogen attack and influences chromatin-associated stress memory as well as transgenerational inheritance of epigenetic imprints in plants. The review further summarizes some research studies on how histone modifications and DNA methylations resist pathogenic and pest attacks by activating defense-related genes.