Transcriptome-Wide Analysis of RNA m6A Methylation and Gene Expression Changes Among Two Arabidopsis Ecotypes and Their Reciprocal Hybrids

Analyzing the relationship of RNA and DNA methylation with gene expression

BACKGROUND

DNA 5-methylcytosine (5mC) and RNA N6-methyladenosine (m6A) methylation are prevalent modifications in eukaryotes, both playing crucial roles in gene regulation. Recent studies have explored their crosstalk and impact on transcription. However, the intricate relationships among 5mC, m6A, and gene expression remain incompletely elucidated.

RESULTS

We collect data on 5mC, m6A, and gene expression from samples from three tissues from each of four pregnant cattle and sheep. We construct a comprehensive genome-wide self-interaction (same gene) and across-interaction (across genes) network of 5mC and m6A within gene-bodies or promoters and gene expression in both species. Qualitative analysis identifies uniquely expressed genes with specific m6A methylation in each tissue from both species. A quantitative comparison of gene expression ratio between methylated and unmethylated genes for m6A within gene body and promoter, and 5mC within gene body and promoter confirms the positive effect of RNA methylation on gene expression. Importantly, the influence of RNA methylation on gene expression is stronger than that of DNA methylation. The predominant self- and across-interactions are between RNA methylation within gene bodies and gene expression, as well as between RNA methylation within promoters and gene expression in both species.

CONCLUSIONS

RNA methylation has a stronger effect on gene expression than does DNA methylation within gene bodies and promoters. DNA and RNA methylation in gene-bodies has a greater impact on gene expression than those in promoters. These findings deepen comprehension of the dynamics and complex relationships among the epigenome, epitranscriptome, and transcriptome, offering fresh insights for advancing epigenetics research.

PEA-m6A: an ensemble learning framework for accurately predicting N6-methyladenosine modifications in plants

N 6-methyladenosine (m6A), which is the mostly prevalent modification in eukaryotic mRNAs, is involved in gene expression regulation and many RNA metabolism processes. Accurate prediction of m6A modification is important for understanding its molecular mechanisms in different biological contexts. However, most existing models have limited range of application and are species-centric. Here we present PEA-m6A, a unified, modularized and parameterized framework that can streamline m6A-Seq data analysis for predicting m6A-modified regions in plant genomes. The PEA-m6A framework builds ensemble learning-based m6A prediction models with statistic-based and deep learning-driven features, achieving superior performance with an improvement of 6.7% to 23.3% in the area under precision-recall curve compared with state-of-the-art regional-scale m6A predictor WeakRM in 12 plant species. Especially, PEA-m6A is capable of leveraging knowledge from pretrained models via transfer learning, representing an innovation in that it can improve prediction accuracy of m6A modifications under small-sample training tasks. PEA-m6A also has a strong capability for generalization, making it suitable for application in within- and cross-species m6A prediction. Overall, this study presents a promising m6A prediction tool, PEA-m6A, with outstanding performance in terms of its accuracy, flexibility, transferability, and generalization ability. PEA-m6A has been packaged using Galaxy and Docker technologies for ease of use and is publicly available at https://github.com/cma2015/PEA-m6A.

Principles, functions, and biological implications of m6A in plants

Over the past decade, N 6-methyladenosine (m6A) has emerged as a prevalent and dynamically regulated modification across the transcriptome; it has been reversibly installed, removed, and interpreted by specific binding proteins, and has played crucial roles in molecular and biological processes. Within this scope, we consolidate recent advancements of m6A research in plants regarding gene expression regulation, diverse physiologic and pathogenic processes, as well as crop trial implications, to guide discussions on challenges associated with and leveraging epitranscriptome editing for crop improvement.

Current technical advancements in plant epitranscriptomic studies

The growth and development of plants are the result of the interplay between the internal developmental programming and plant-environment interactions. Gene expression regulations in plants are made up of multi-level networks. In the past few years, many studies were carried out on co- and post-transcriptional RNA modifications, which, together with the RNA community, are collectively known as the "epitranscriptome." The epitranscriptomic machineries were identified and their functional impacts characterized in a broad range of physiological processes in diverse plant species. There is mounting evidence to suggest that the epitranscriptome provides an additional layer in the gene regulatory network for plant development and stress responses. In the present review, we summarized the epitranscriptomic modifications found so far in plants, including chemical modifications, RNA editing, and transcript isoforms. The various approaches to RNA modification detection were described, with special emphasis on the recent development and application potential of third-generation sequencing. The roles of epitranscriptomic changes in gene regulation during plant-environment interactions were discussed in case studies. This review aims to highlight the importance of epitranscriptomics in the study of gene regulatory networks in plants and to encourage multi-omics investigations using the recent technical advancements.

Epigenetic control of plant senescence and cell death and its application in crop improvement

Plant senescence is the last stage of plant development and a type of programmed cell death, occurring at a predictable time and cell. It involves the functional conversion from nutrient assimilation to nutrient remobilization, which substantially impacts plant architecture and plant biomass, crop quality, and horticultural ornamental traits. In past two decades, DNA damage was believed to be a main reason for cell senescence. Increasing evidence suggests that the alteration of epigenetic information is a contributing factor to cell senescence in organisms. In this review, we summarize the current research progresses of epigenetic and epitranscriptional mechanism involved in cell senescence of plant, at the regulatory level of DNA methylation, histone methylation and acetylation, chromatin remodeling, non-coding RNAs and RNA methylation. Furthermore, we discuss their molecular genetic manipulation and potential application in agriculture for crop improvement. Finally we point out the prospects of future research topics.

Examining parent-of-origin effects on transcription and RNA methylation in mediating aggressive behavior in honey bees (Apis mellifera)

Conflict between genes inherited from the mother (matrigenes) and the father (patrigenes) is predicted to arise during social interactions among offspring if these genes are not evenly distributed among offspring genotypes. This intragenomic conflict drives parent-specific transcription patterns in offspring resulting from parent-specific epigenetic modifications. Previous tests of the kinship theory of intragenomic conflict in honey bees (Apis mellifera) provided evidence in support of theoretical predictions for variation in worker reproduction, which is associated with extreme variation in morphology and behavior. However, more subtle behaviors - such as aggression - have not been extensively studied. Additionally, the canonical epigenetic mark (DNA methylation) associated with parent-specific transcription in plant and mammalian model species does not appear to play the same role as in honey bees, and thus the molecular mechanisms underlying intragenomic conflict in this species is an open area of investigation. Here, we examined the role of intragenomic conflict in shaping aggression in honey bee workers through a reciprocal cross design and Oxford Nanopore direct RNA sequencing. We attempted to probe the underlying regulatory basis of this conflict through analyses of parent-specific RNA m6A and alternative splicing patterns. We report evidence that intragenomic conflict occurs in the context of honey bee aggression, with increased paternal and maternal allele-biased transcription in aggressive compared to non-aggressive bees, and higher paternal allele-biased transcription overall. However, we found no evidence to suggest that RNA m6A or alternative splicing mediate intragenomic conflict in this species.

m6A mRNA modification promotes chilling tolerance and modulates gene translation efficiency in Arabidopsis

N 6-methyladenosine (m6A), the most prevalent mRNA modification in eukaryotes, is an emerging player of gene regulation at transcriptional and translational levels. Here, we explored the role of m6A modification in response to low temperature in Arabidopsis (Arabidopsis thaliana). Knocking down mRNA adenosine methylase A (MTA), a key component of the modification complex, by RNA interference (RNAi) led to drastically reduced growth at low temperature, indicating a critical role of m6A modification in the chilling response. Cold treatment reduced the overall m6A modification level of mRNAs especially at the 3' untranslated region. Joint analysis of the m6A methylome, transcriptome and translatome of the wild type (WT) and the MTA RNAi line revealed that m6A-containing mRNAs generally had higher abundance and translation efficiency than non-m6A-containing mRNAs under normal and low temperatures. In addition, reduction of m6A modification by MTA RNAi only moderately altered the gene expression response to low temperature but led to dysregulation of translation efficiencies of one third of the genes of the genome in response to cold. We tested the function of the m6A-modified cold-responsive gene ACYL-COA:DIACYLGLYCEROL ACYLTRANSFERASE 1 (DGAT1) whose translation efficiency but not transcript level was reduced in the chilling-susceptible MTA RNAi plant. The dgat1 loss-of-function mutant exhibited reduced growth under cold stress. These results reveal a critical role of m6A modification in regulating growth under low temperature and suggest an involvement of translational control in chilling responses in Arabidopsis.

Detection, regulation, and functions of RNA N6-methyladenosine modification in plants

N⁶-Methyladenosine (m⁶A) is the most abundant internal chemical modification in eukaryotic mRNA and plays important roles in gene expression regulation, including transcriptional and post-transcriptional regulation. m⁶A is a reversible modification that is installed, removed, and recognized by methyltransferases (writers), demethylases (erasers), and m⁶A-binding proteins (readers), respectively. Recently, the breadth of research on m⁶A in plants has expanded, and the vital roles of m⁶A in plant development, biotic and abiotic stress responses, and crop trait improvement have been investigated. In this review, we discuss recent developments in research on m⁶A and highlight the detection methods, distribution, regulatory proteins, and molecular and biological functions of m⁶A in plants. We also offer some perspectives on future investigations, providing direction for subsequent research on m⁶A in plants.

N6 -Methyladenosine modification of mRNA contributes to the transition from 2D to 3D growth in the moss Physcomitrium patens

Plants colonized the land approximately 470 million years ago, coinciding with the development of apical cells that divide in three planes. The molecular mechanisms that underly the development of the 3D growth pattern are poorly understood, mainly because 3D growth in seed plants starts during embryo development. In contrast, the transition from 2D to 3D growth in the moss Physcomitrium patens has been widely studied, and it involves a large turnover of the transcriptome to allow the establishment of stage-specific transcripts that facilitate this developmental transition. N⁶ -Methyladenosine (m⁶ A) is the most abundant, dynamic and conserved internal nucleotide modification present on eukaryotic mRNA and serves as a layer of post-transcriptional regulation directly affecting several cellular processes and developmental pathways in many organisms. In Arabidopsis, m⁶ A has been reported to be essential for organ growth and determination, embryo development and responses to environmental signals. In this study, we identified the main genes of the m⁶ A methyltransferase complex (MTC), MTA, MTB and FIP37, in P. patens and demonstrate that their inactivation leads to the loss of m⁶ A in mRNA, a delay in the formation of gametophore buds and defects in spore development. Genome-wide analysis revealed several transcripts affected in the Ppmta background. We demonstrate that the PpAPB1-PpAPB4 transcripts, encoding central factors orchestrating the transition from 2D to 3D growth in P. patens, are modified by m⁶ A, whereas in the Ppmta mutant the lack of the m⁶ A marker is associated with a corresponding decrease in transcript accumulation. Overall, we suggest that m⁶ A is essential to enable the proper accumulation of these and other bud-specific transcripts directing the turnover of stage-specific transcriptomes, and thus promoting the transition from protonema to gametophore buds in P. patens.

Ribosome profiling reveals the translational landscape and allele-specific translational efficiency in rice

Translational regulation is a critical step in the process of gene expression and governs the synthesis of proteins from mRNAs. Many studies have revealed translational regulation in plants in response to various environmental stimuli. However, there have been no studies documenting the comprehensive landscape of translational regulation and allele-specific translational efficiency in multiple plant tissues, especially those of rice, a main staple crop that feeds nearly half of the world's population. Here we used RNA sequencing and ribosome profiling data to analyze the transcriptome and translatome of an elite hybrid rice, Shanyou 63 (SY63), and its parental varieties Zhenshan 97 and Minghui 63. The results revealed that gene expression patterns varied more among tissues than among varieties at the transcriptional and translational levels. We identified 3392 upstream open reading frames (uORFs), and the uORF-containing genes were enriched in transcription factors. Only 668 of 13 492 long non-coding RNAs could be translated into peptides. Finally, we discovered numerous genes with allele-specific translational efficiency in SY63 and demonstrated that some cis-regulatory elements may contribute to allelic divergence in translational efficiency. Overall, these findings may improve our understanding of translational regulation in rice and provide information for molecular breeding research.

Recent advances in the plant epitranscriptome

Chemical modifications of RNAs, known as the epitranscriptome, are emerging as widespread regulatory mechanisms underlying gene regulation. The field of epitranscriptomics advances recently due to improved transcriptome-wide sequencing strategies for mapping RNA modifications and intensive characterization of writers, erasers, and readers that deposit, remove, and recognize RNA modifications, respectively. Herein, we review recent advances in characterizing plant epitranscriptome and its regulatory mechanisms in post-transcriptional gene regulation and diverse physiological processes, with main emphasis on N⁶-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C). We also discuss the potential and challenges for utilization of epitranscriptome editing in crop improvement.

EGFR and ERK activation resists flavonoid quercetin-induced anticancer activities in human cervical cancer cells in vitro

In the present study, due to the complex and numerous targets of Sarcandrae Herb (also known as Zhong Jie Feng), network pharmacology was performed to analyze its therapeutic effect on 2 cervical cancer cell lines, which could assist with the development of novel therapies. The results suggested that the natural flavonoid quercetin (Que), the effective antitumor ingredient in SH, which is widely present in a variety of plants, may depend on the target, EGFR. Previous studies have shown that EGFR serves a crucial role in the occurrence and development of cervical cancer, but its downstream molecules and regulatory mechanisms remain unknown. The anti-cervical cancer cell properties of Que, which are present in ubiquitous plants, were examined in vitro to identify the association between Que and its underlying pathway using MTT assays, flow cytometry, western blot analysis and Transwell assays. It was found that Que reduced cervical cancer cell viability, promoted G₂/M phase cell cycle arrest and cell apoptosis, as well as inhibited cell migration and invasion. The Tyr1068 phosphorylation site of EGFR and the corresponding ERK target were also examined and the 2 kinases were markedly activated by Que. Furthermore, the EGFR inhibitor, afatinib and the ERK inhibitor, U0126 blocked the increase of EGFR and ERK phosphorylation, and resulted in a notable enhancement of apoptosis and cell cycle arrest. Therefore, to the best of our knowledge, the current results provided the first evidence that EGFR and ERK activation induced by Que could resist Que-induced anticancer activities. On this basis, the present study determined the role of EGFR and the underlying signaling pathways involved in the anti-cervical cancer malignant behavior induced by Que and identified the negative regulatory association.