Identifying RNA Modifications by Direct RNA Sequencing Reveals Complexity of Epitranscriptomic Dynamics in Rice

Simultaneous nanopore profiling of mRNA m6A and pseudouridine reveals translation coordination

N6-methyladenosine (m6A) and pseudouridine (Ψ) are the two most abundant modifications in mammalian messenger RNA, but the coordination of their biological functions remains poorly understood. We develop a machine learning-based nanopore direct RNA sequencing method (NanoSPA) that simultaneously analyzes m6A and Ψ in the human transcriptome. Applying NanoSPA to polysome profiling, we reveal opposing transcriptomic co-occurrence of m6A and Ψ and synergistic, hierarchical effects of m6A and Ψ on the polysome.

Nanopore sequencing: flourishing in its teenage years

Over the past decade, nanopore sequencing has experienced significant advancements and changes, transitioning from an initially emerging technology to a significant instrument in the field of genomic sequencing. However, as advancements in next-generation sequencing technology persist, nanopore sequencing also improves. This paper reviews the developments, applications, and outlook on nanopore sequencing technology. Currently, nanopore sequencing supports both DNA and RNA sequencing, making it widely applicable in areas such as telomere-to-telomere (T2T) genome assembly, direct RNA sequencing (DRS), and metagenomics. The openness and versatility of nanopore sequencing have established it as a preferred option for an increasing number of research teams, signaling a transformative influence on life science research. As the nanopore sequencing technology advances, it provides a faster, more cost-effective approach with extended read lengths, demonstrating the significant potential for complex genome assembly, pathogen detection, environmental monitoring, and human disease research, offering a fresh perspective in sequencing technologies.

Nanopore direct RNA sequencing reveals N6-methyladenosine and polyadenylation landscapes on long non-coding RNAs in Arabidopsis thaliana

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in various biological processes, including stage development in plants. N6-methyladenosine (m6A) modification and polyadenylation are noteworthy regulatory processes that impact transcript functions by modulating their abundance. However, the specific landscapes of m6A modification and polyadenylation on lncRNAs remain largely unexplored. The advent of nanopore direct RNA sequencing (DRS) provides unprecedented opportunities for directly detecting m6A modifications and estimating polyadenine (poly[A]) tail lengths on individual RNA molecules.

RESULTS

Here we utilized nanopore DRS to identify lncRNAs and map the transcriptome-wide m6A modification and polyadenylation landscapes in the model plant Arabidopsis thaliana. Leveraging the Low-abundance Aware Full-length Isoform clusTEr (LAFITE) assembly pipeline, we identified 1149 novel lncRNAs in seventeen nanopore DRS datasets from the wild-type Columbia-0. Through the precise detection of 2381 m6A modification sites on lncRNAs at single-base resolution, we observed that lncRNAs exhibited lower methylation levels compared to protein-coding RNAs, and m6A modification facilitated lncRNA abundance. Additionally, we estimated the poly(A) tail lengths of individual lncRNAs and found that poly(A) tails contributed to lncRNA stability, while their effect was not length-dependent. Furthermore, by comparing lncRNA abundance between 2-week seedlings and 5-week floral buds, we revealed the dynamic expression patterns of lncRNAs during the transition from the vegetative stage to the reproductive stage. These observations provided insights into their potential roles in specific tissues or stages in Arabidopsis, including regulating stage development. Moreover, by integrating information on m6A modification, we unveiled a positive correlation between methylation variances and differential expressions of lncRNAs during stage development.

CONCLUSIONS

These findings highlighted the significance of epigenetic modification and post-transcriptional processing in shaping lncRNA expression and their functions during Arabidopsis stage development, contributing to the growing field of lncRNA research in plants.

CLINICAL TRIAL NUMBER

Not applicable.

Transcriptome-wide N6-methyladenosine profiling reveals growth-defense trade-offs in the response of rice to brown planthopper (Nilaparvata lugens) infestation

BACKGROUND

N6-Methyladenosine (m6A) is a common messenger RNA (mRNA) modification that affects various physiological processes in stress responses. However, the role of m6A modifications in plants responses to herbivore stress remains unclear.

RESULTS

Here, we found that an infestation of brown planthopper (Nilaparvata lugens) female adults enhanced the resistance of rice to N. lugens. The m6A methylome analysis of N. lugens-infested and uninfested rice samples was performed to explore the interaction between rice and N. lugens. The m6A methylation mainly occurred in genes that were actively expressed in rice following N. lugens infestation, while an analysis of the whole-genomic mRNA distribution of m6A showed that N. lugens infestation caused an overall decrease in the number of m6A methylation sites across the chromosomes. The m6A methylation of genes involved in the m6A modification machinery and several defense-related phytohormones (jasmonic acid and salicylic acid) pathways was increased in N. lugens-infested rice compared to that in uninfested rice. In contrast, m6A modification levels of growth-related phytohormone (auxin and gibberellin) biosynthesis-related genes were significantly attenuated during N. lugens infestation, accompanied by the down-regulated expression of these transcripts, indicating that rice growth was restricted during N. lugens attack to rapidly optimize resource allocation for plant defense. Integrative analysis of the differential patterns of m6A methylation and the corresponding transcripts showed a positive correlation between m6A methylation and transcriptional regulation.

CONCLUSION

The m6A modification is an important strategy for regulating the expression of genes involved in rice defense and growth during rice-N. lugens interactions. These findings provide new ideas for formulating strategies to control herbivorous pests. © 2024 Society of Chemical Industry.

OsMBF1a Facilitates Seed Germination by Regulating Biosynthesis of Gibberellic Acid and Abscisic Acid in Rice

Seed germination is a pivotal stage in the plant life cycle, orchestrated by a myriad of internal and external factors, notably plant hormones. The underlying molecular mechanisms governing rice seed germination remain largely unelucidated. Herein, we uncover OsMBF1a as a crucial regulatory factor that employs a dual strategy to promote seed germination: positively activating genes involved in gibberellin (GA) biosynthesis pathways, while negatively regulating key genes responsible for abscisic acid (ABA) synthesis. Furthermore, OsMBF1a modulates the endogenous levels of ABA and GA in rice seeds, reinforcing its central role in the germination process. The expression of ZmMBF1a and ZmMBF1b, the homologous genes in maize, in rice seeds similarly affects germination, indicating the conserved functionality of MBF1 family genes in regulating seed germination. This study provides novel insights into the molecular mechanisms underlying rice seed germination and underscores the significance of MBF1 family genes in plant growth and development.

Nanopore Direct RNA Sequencing Reveals the Short-Term Salt Stress Response in Maize Roots

Transcriptome analysis, relying on the cutting-edge sequencing of cDNA libraries, has become increasingly prevalent within functional genome studies. However, the dependence on cDNA in most RNA sequencing technologies restricts their ability to detect RNA base modifications. To address this limitation, the latest Oxford Nanopore Direct RNA Sequencing (ONT DRS) technology was employed to investigate the transcriptome of maize seedling roots under salt stress. This approach aimed to unveil both the RNA transcriptional profiles and alterations in base modifications. The analysis of the differential expression revealed a total of 1398 genes and 2223 transcripts that exhibited significant variation within the maize root system following brief exposure to salt stress. Enrichment analyses, such as the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway assessments, highlighted the predominant involvement of these differentially expressed genes (DEGs) in regulating ion homeostasis, nitrogen metabolism, amino acid metabolism, and the phytohormone signaling pathways. The protein-protein interaction (PPI) analysis showed the participation of various proteins related to glycolytic metabolism, nitrogen metabolism, amino acid metabolism, abscisic acid signaling, and the jasmonate signaling pathways. It was through this intricate molecular network that these proteins collaborated to safeguard root cells against salt-induced damage. Moreover, under salt stress conditions, the occurrence of variable shear events (AS) in RNA modifications diminished, the average length of poly(A) tails underwent a slight decrease, and the number of genes at the majority of the variable polyadenylation (APA) sites decreased. Additionally, the levels of N5-methylcytosine (m5C) and N6-methyladenosine (m6A) showed a reduction. These results provide insights into the mechanisms of early salt tolerance in maize.