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.

ß-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.

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.

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.

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.

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.

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.

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.

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.

Comparative genome-wide analysis of circular RNAs in Brassica napus L.: target-site versus non-target-site resistance to herbicide stress.

Circular RNAs (circRNAs), a class of non-coding RNA molecules, are recognized for their unique functions; however, their responses to herbicide stress in Brassica napus remain unclear. In this study, the role of circRNAs in response to herbicide treatment was investigated in two rapeseed cultivars: MH33, which confers non-target-site resistance (NTSR), and EM28, which exhibits target-site resistance (TSR). The genome-wide circRNA profiles of herbicide-stressed and non-stressed seedlings were analyzed. The findings indicate that NTSR seedlings exhibited a greater abundance of circRNAs, shorter lengths of circRNAs and their parent genes, and more diverse functions of parent genes compared with TSR seedlings. Compared to normal-growth plants, the herbicide-stressed group exhibited similar trends in the number of circRNAs, functions of parent genes, and differentially expressed circRNAs as observed in NTSR seedlings. In addition, a greater number of circRNAs that function as competing microRNA (miRNA) sponges were identified in the herbicide stress and NTSR groups compared to the normal-growth and TSR groups, respectively. The differentially expressed circRNAs were validated by qPCR. The differntially expressed circRNA-miRNA networks were predicted, and the mRNAs targeted by these miRNAs were annotated. Our results suggest that circRNAs play a crucial role in responding to herbicide stress, exhibiting distinct responses between NTSR and TSR in rapeseed. These findings offer valuable insights into the mechanisms underlying herbicide resistance in rapeseed.

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.

Insights into the Early Steps of the Symbiotic Interaction between Soybean (Glycine max) and Sinorhizobium fredii Symbiosis Using Transcriptome, Small RNA, and Degradome Sequencing.

Symbiotic nitrogen fixation carried out by the soybean-rhizobia symbiosis increases soybean yield and reduces the amount of nitrogen fertilizer that has been applied. MicroRNAs (miRNAs) are crucial in plant growth and development, prompting an investigation into their role in the symbiotic interaction of soybean with partner rhizobia. Through integrated small RNA, transcriptome, and degradome sequencing analysis, 1215 known miRNAs, 314 of them conserved, and 187 novel miRNAs were identified, with 44 differentially expressed miRNAs in soybean roots inoculated with Sinorhizobium fredii HH103 and a ttsI mutant. The study unveiled that the known miRNA gma-MIR398a-p5 was downregulated in the presence of the ttsI mutation, while the target gene of gma-MIR398a-p5, Glyma.06G007500, associated with nitrogen metabolism, was upregulated. The results of this study offer insights for breeding high-efficiency nitrogen-fixing soybean varieties, enhancing crop yield and quality.

Genome-wide characterization of drought-responsive long non-coding RNAs in sorghum (Sorghum bicolor).

Drought stress strongly affects crop yield. Although knowledge of long non-coding RNAs (lncRNAs) has been updated continuously and rapidly, information about lncRNAs in drought resistance regulation is extremely limited in sorghum. Here, lncRNA-sequencing was performed with seedlings of a sorghum cultivar (Jinza29) under three water control treatments to investigate the mechanism of lncRNAs responsible for drought resistance in sorghum. A total of 377 differentially expressed lncRNAs (DElncRNAs) were identified. We also predicted 4322 and 2827 transcripts as potential cis-target and trans-target genes for drought-responsive lncRNAs, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that those target genes exhibited marked enrichment into "oxidoreductase activity", "signal transducer activity", "DNA repair", "photosynthesis", "glutathione metabolism", and "phenylpropanoid biosynthesis" and other terms associated with abiotic stress resistance. Moreover, several lncRNAs were estimated to modulate the expression of other genes related to stress response and photosynthetic carbon metabolism. Additionally, we found 107 DElncRNAs that might be candidate target mimics for 56 miRNAs. LncRNAs play important roles in drought adaptation of sorghum through interacting with protein-encoding genes. The obtained results provided novel insights into the biological characteristics of lncRNAs and offered potential regulatory factors for genetically enhancing drought resistance in sorghum.

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.