Coordination of m6A mRNA methylation and gene transcriptome in rice response to cadmium stress

Exogenous MgH2-derived hydrogen alleviates cadmium toxicity through m6A RNA methylation in rice

Cadmium (Cd) contamination poses a substantial threat to crop yields and human health. While magnesium hydride (MgH2) has been reported as a hydrogen (H2) donor that promotes plant growth under heavy metal contamination, its role in rice remains elusive. Herein, seedlings of Oryza sativa L. Japonica variety Zhonghua 11 (ZH11) were selected and exposed to 20 µL of 1-mol/L cadmium chloride (CdCl2) solution via hydroponics to simulate Cd stress. Meanwhile, 0.1 mg of MgH2 was used to slow-release H2 to the experimental group to explore its potential effects on rice over a 2-week period. The results indicated that Cd exposure severely inhibited the growth and development of ZH11 rice seedlings. However, the exogenous slow-release of H2 from MgH2 effectively mitigated this inhibitory effect by restoring the balance of reactive oxygen species (ROS), maintaining endogenous H2 homeostasis, and supporting the photosynthetic system. High-performance liquid chromatography analysis revealed that exogenous H2 reduces m6A RNA methylation levels in mRNA under Cd stress. Consequently, MeRIP-seq was conducted to investigate the effect of Cd exposure in rice in the presence and absence of H2. The m6A modifications were enriched at the start codon, stop codon, and 3' UTR. By integrating RNA-seq data, 118 transcripts were identified as differentially methylated and expressed genes under Cd stress. These gene annotations were associated with ROS, biological stress, and hormonal responses. Notably, 297 differentially methylated and expressed genes were identified under Cd stress in the presence of H2, linked to heavy metals, protein kinases, and calcium signaling regulation. Cd strongly activates the MAPK pathway in response to stress. Exogenous H2 reduces Cd accumulation as well as enhances plant tolerance and homeostasis by lowering m6A levels, thereby decreasing the mRNA stability of these genes. Our findings indicate that MgH2, by supplying H2, regulates gene expression through m6A RNA methylation and confers Cd tolerance in rice. This study provides potential candidate genes for studying the remediation of heavy metal pollution in plants.

An mRNA methylase and demethylase regulate sorghum salt tolerance by mediating N6-methyladenosine modification

N 6-methyladenosine (m6A) modification is a crucial and widespread molecular mechanism governing plant development and stress tolerance. The specific impact of m6A regulation on plants with inherently high salt tolerance remains unclear. Existing research primarily focuses on the overexpression or knockout of individual writer or eraser components to alter m6A levels. However, a comprehensive study simultaneously altering overall m6A modification levels within the same experiment is lacking. Such an investigation is essential to determine whether opposing changes in m6A modification levels exert entirely different effects on plant salt tolerance. In this study, we identified the major writer member mRNA adenosine methylase A (SbMTA) in sorghum (Sorghum bicolor) as critical for sorghum survival. The sbmta mutant exhibits a phenotype characterized by reduced overall m6A, developmental arrest, and, ultimately, lethality. Overexpression of SbMTA increased m6A levels and salt tolerance, while overexpression of the m6A eraser alkylated DNA repair protein AlkB homolog 10B (SbALKBH10B) in sorghum showed the opposite phenotype. Comparative analyses between sorghum with different m6A levels reveal that SbMTA- and SbALKBH10B-mediated m6A alterations significantly impact the stability and expression levels of genes related to the abscisic acid signaling pathway and growth under salt stress. In summary, this study unveils the intricate relationship between m6A modifications and salt tolerance in sorghum, providing valuable insights into how m6A modification levels on specific transcripts influence responses to salt stress.

Multi-kernel feature extraction with dynamic fusion and downsampled residual feature embedding for predicting rice RNA N6-methyladenine sites

RNA N$^{6}$-methyladenosine (m$^{6}$A) is a critical epigenetic modification closely related to rice growth, development, and stress response. m$^{6}$A accurate identification, directly related to precision rice breeding and improvement, is fundamental to revealing phenotype regulatory and molecular mechanisms. Faced on rice m$^{6}$A variable-length sequence, to input into the model, the maximum length padding and label encoding usually adapt to obtain the max-length padded sequence for prediction. Although this can retain complete sequence information, resulting in sparse information and invalid padding, reducing feature extraction accuracy. Simultaneously, existing rice-specific m$^{6}$A prediction methods are still at an early stage. To address these issues, we develop a new end-to-end deep learning framework, MFDm$^{6}$ARice, for predicting rice m$^{6}$A sites. In particular, to alleviate sparseness, we construct a multi-kernel feature fusion module to mine essential information in max-length padded sequences by multi-kernel feature extraction function and effectively transfer information through global-local dynamic fusion function. Concurrently, considering the complexity and computational efficiency of high-dimensional features caused by invalid padding, we design a downsampling residual feature embedding module to optimize feature space compression and achieve accurate feature expression and efficient computational performance. Experiments show that MFDm$^{6}$ARice outperforms comparison methods in cross-validation, same- and cross-species independent test sets, demonstrating good robustness and generalization. The application on maize m$^{6}$A indicates the MFDm$^{6}$ARice's scalability. Further investigations have shown that combining different kernel features, focusing on global channel-local spatial, and employing reasonable downsampling and residual connections can improve feature representation and extraction, ensure effective information transfer, and significantly enhance model performance.

Comprehensive analysis of epigenetic modifications in alfalfa under cadmium stress

Epigenetics plays an important role in plant growth and development and in environmental adaptation. Alfalfa, an important forage crop, is rich in nutrients. However, little is known about the molecular regulatory mechanisms underlying the response of alfalfa to cadmium (Cd) stress. Here, we performed DNA methylation (5mC), RNA methylation (m6A) and transcriptomic sequencing analyses of alfalfa roots under Cd stress. Whole-genome methylation sequencing and transcriptomic sequencing revealed that Cd stress reduced DNA methylation levels. Moreover, a reduced 5mC methylation level was associated with decreased expression of several DNA methyltransferase genes. Compared with those under normal (CK) conditions, the m6A modification levels under Cd stress were greater and were positively correlated with gene expression in alfalfa roots. We also found a negative correlation between the 5mC level and the m6A level, especially in CG and CHG contexts. In yeast, the overexpression of MsNARMP5 (natural resistance-associated macrophage protein) and MsPCR2 (plant cadmium resistance 2), which are modified by 5mC or m6A, significantly increased Cd stress tolerance. These results provide candidate genes for future studies on the mechanism of Cd stress tolerance in alfalfa roots and valuable information for studying heavy metal stress in alfalfa breeding.

ZnO nanoparticles enhances cadmium tolerance by modulating N6-methyladenosine (m6A) level of stress-responsive genes NRT1 and GM35E in vegetable soybean

Cadmium (Cd) is a hazardous heavy metal pollutant that poses significant risks to agricultural production and human health. Nanoparticles (NPs) can alleviate the effects of cadmium on crops by regulating the expression of stress-responsive genes, however, the mechanism of regulation is unknown. N6-methyladenosine (m6A) is a prevalent RNA modification, which determines the expression level of RNA. In this study, we performed m6A methylome analysis and gene editing to investigate the regulation of stress-responsive genes by m6A under Cd stress with ZnO nanoparticles (ZnO NPs). Firstly, we identified 16 differentially expressed genes (DEGs) with differential m6A-modification by m6A methylome seq, which included 6 stress-responsive genes. ZnO NPs treatment reduced the m6A level and stabilized the mRNAs of these stress-responsive genes. Then, we utilized the in vivo m6A modification tool, Plant m6A Editors (PMEs), specifically reduced the m6A level of NRT1 and GM35E in transgenic plants. The expression of NRT1 and GM35E was increased, and plants carrying PMEs-NRT1 and PMEs-GM35E showed stronger Cd tolerance than control. Therefore, we can conclude that NPs enhance Cd tolerance in plants by regulating the m6A methylation level of stress-responsive genes such as NRT1 and GM35E, which ultimately affects the expression of these genes. This study proposed a post-transcriptional regulatory network in vegetable soybean roots responding to Cd with ZnO NPs, potentially offering new insights into gene manipulation for controlling low Cd accumulation in crops.

Transcriptome-wide methylated RNA immunoprecipitation sequencing profiling reveals m6A modification involved in response to heat stress in Apostichopus japonicus

BACKGROUND

Global warming-induced environmental stresses have diverse effects on gene expression and regulation in the life processes of various aquatic organisms. N6 adenylate methylation (m6A) modifications are known to influence mRNA transcription, localization, translation, stability, splicing, and nuclear export, which are pivotal in mediating stress responses. Apostichopus japonicus is a significant species in aquaculture and a representative of benthic organisms in ecosystems, thus there is a growing need for research on its heat stress mechanism.

RESULTS

In this study, m6A-modified whole transcriptome profiles of the respiratory tree tissues of A. japonicus in the control (T18) and high-temperature stress (T32) groups were obtained using MeRIP-seq technology. The results showed that 7,211 common m6A peaks, and 9,459 genes containing common m6A were identified in three replicates T18 and T32 groups. The m6A peaks were found to be highly enriched in the 3' untranslated region, and the common sequence of the m6A peak was also enriched, which was shown as RRACH (R = G or A; H = A, C, or U). A total of 1,200 peaks were identified as significantly differentially enriched in the T32 group compared with the T18 group. Among them, 245 peaks were upregulated and 955 were downregulated, which indicated that high temperature stress significantly altered the methylation pattern of m6A, and there were more demethylation sites in the T32 group. Conjoint analysis of the m6A methylation modification and the transcript expression level (the MeRIP-seq and RNA-seq data) showed co-differentiated 395 genes were identified, which were subsequently divided into four groups with a predominant pattern that more genes with decreased m6A modification and up-regulated expression, including HSP70IV, EIF2AK1, etc. GO enrichment and KEGG analyses of differential m6A peak related genes and co-differentiated genes showed the genes were significantly associated with transcription process and pathways such as protein processing in the endoplasmic reticulum, Wnt signaling pathway, and mTOR signaling pathway, etc. CONCLUSION: The comparisons of m6A modification patterns and conjoint analyses of m6A modification and gene expression profiles suggest that m6A modification was involved in the regulation of heat stress-responsive genes and important functional pathways in A. japonicus in response to high-temperature stress. The study will contribute to elucidate the regulatory mechanism of m6A modification in the response of A. japonicus to environmental stress, as well as the conservation and utilization of sea cucumber resources in the context of environmental changes.

Epigenetic control of plant abiotic stress responses

On top of genetic information, organisms have evolved complex and sophisticated epigenetic regulation to adjust gene expression in response to developmental and environmental signals. Key epigenetic mechanisms include DNA methylation, histone modifications and variants, chromatin remodeling, and chemical modifications of RNAs. Epigenetic control of environmental responses is particularly important for plants, which are sessile and unable to move away from adverse environments. Besides enabling plants to rapidly respond to environmental stresses, some stress-induced epigenetic changes can be maintained, providing plants with a pre-adapted state to recurring stresses. Understanding these epigenetic mechanisms offers valuable insights for developing crop varieties with enhanced stress tolerance. Here, we focus on abiotic stresses and summarize recent progress in characterizing stress-induced epigenetic changes and their regulatory mechanisms and roles in plant abiotic stress resistance.

Transcriptome-wide m6A methylation profile reveals tissue specific regulatory networks in switchgrass (Panicum virgatum L.) under cadmium stress

The heavy metal cadmium (Cd), known for its high toxicity, poses a grave threat to human health through the food chain. N6-methyladenosine (m6A), the most abundant internal modification, regulates plant adaptation to various adversities, yet the panorama of m6A modifications in switchgrass under cadmium stress remains elusive. This study examines the physiological responses of switchgrass roots and shoots exposed to 50 μM CdCl2, alongside an overview of transcriptome-wide m6A methylation patterns. After cadmium treatment, methylation modifications are primarily enriched near stop codons and the 3'UTR region, with a negative correlation between m6A modification and gene expression levels. In shoots, approximately 58 % of DEGs with m6A modifications show upregulation in expression and decrease in m6A peaks, including zinc transporter 4-like (ZIP4). In roots, about 43 % of DEGs with m6A modifications exhibit downregulation in expression and increase in m6A peaks, such as the ABC transporter family member (ABCG25). We further validate the m6A enrichment, gene expression and mRNA stability of ZIP4 in response to Cd treatment. The results suggest that the negative correlation of m6A enrichment and gene expression is due to altered mRNA stability. Our study establishes an m6A regulatory network governing cadmium transport in switchgrass roots and shoots, offering new avenues for candidate gene manipulation in phytoremediation applications of heavy metal pollution.

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.

Cadmium exposure elicited dynamic RNA m6A modification and epi-transcriptomic regulation in the Pacific whiteleg shrimp Litopenaeus vannamei

N6-methyladenosine (m6A) methylation is the most prevalent post-transcriptional RNA modification in eukaryotic organisms, but its roles in the regulation of physiological resistance of marine crustaceans to heavy metal pollutants are poorly understood. In this study, the transcriptome-wide m6A RNA methylation profiles and dynamic m6A changes induced by acute Cd2+ exposure in the the pacific whiteleg shrimp Litopenaeus vannamei were comprehensively analyzed. Cd2+ toxicity caused a significant reduction in global RNA m6A methylation level, with major m6A regulators including the m6A methyltransferase METTL3 and the m6A binding protein YTHDF2 showing declined expression. Totally, 11,467 m6A methylation peaks from 6415 genes and 17,291 peaks within 7855 genes were identified from the Cd2+ exposure group and the control group, respectively. These m6A peaks were predominantly enriched in the 3' untranslated region (UTR) and around the start codon region of the transcripts. 7132 differentially expressed genes (DEGs) and 7382 differentially m6A-methylated genes (DMGs) were identified. 3186 genes showed significant changes in both gene expression and m6A methylation levels upon cadmium exposure, and they were related to a variety of biological processes and gene pathways. Notably, an array of genes associated with antioxidation homeostasis, transmembrane transporter activity and intracellular detoxification processes were significantly enriched, demonstrating that m6A modification may mediate the physiological responses of shrimp to cadmium toxicity via regulating ROS balance, Cd2+ transport and toxicity mitigation. The study would contribute to a deeper understanding of the evolutionary and functional significance of m6A methylation to the physiological resilience of decapod crustaceans to heavy metal toxicants.

Detection, distribution, and functions of RNA N 6-methyladenosine (m6A) in plant development and environmental signal responses

The epitranscriptomic mark N 6-methyladenosine (m6A) is the most common type of messenger RNA (mRNA) post-transcriptional modification in eukaryotes. With the discovery of the demethylase FTO (FAT MASS AND OBESITY-ASSOCIATED PROTEIN) in Homo Sapiens, this modification has been proven to be dynamically reversible. With technological advances, research on m6A modification in plants also rapidly developed. m6A modification is widely distributed in plants, which is usually enriched near the stop codons and 3'-UTRs, and has conserved modification sequences. The related proteins of m6A modification mainly consist of three components: methyltransferases (writers), demethylases (erasers), and reading proteins (readers). m6A modification mainly regulates the growth and development of plants by modulating the RNA metabolic processes and playing an important role in their responses to environmental signals. In this review, we briefly outline the development of m6A modification detection techniques; comparatively analyze the distribution characteristics of m6A in plants; summarize the methyltransferases, demethylases, and binding proteins related to m6A; elaborate on how m6A modification functions in plant growth, development, and response to environmental signals; and provide a summary and outlook on the research of m6A in plants.

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.

Transcriptome-wide m6A methylation profiling identifies GmAMT1;1 as a promoter of lead and cadmium tolerance in soybean nodules

Lead (Pb) and cadmium (Cd) are common heavy metal pollutants that are often found in the soil in soybean agricultural production, adversely impacting symbiotic nitrogen fixation in soybean nodules. In this study, the exposure of soybean nodules to Pb and Cd stress was found to reduce nitrogenase activity. Shifts in the RNA methylation profiles of nodules were subsequently examined by profiling the differential expression of genes responsible for regulating m6A modifications and conducting transcriptome-wide analyses of m6A methylation profiles under Pb and Cd stress condition. Differentially methylated genes (DMGs) that were differentially expressed were closely related to reactive oxygen species activity and integral membrane components. Overall, 19 differentially expressed DMGs were ultimately determined to be responsive to both Pb and Cd stress, including Glyma.20G082450, which encodes GmAMT1;1 and was confirmed to be a positive regulator of nodules tolerance to Pb and Cd. Together, these results are the first published data corresponding to transcriptome-wide m6A methylation patterns in soybean nodules exposed to Cd and Pb stress, and provide novel molecular insight into the regulation of Pb and Cd stress responses in nodules, highlighting promising candidate genes related to heavy metal tolerance, that may also be amenable to application in agricultural production. ENVIRONMENTAL IMPLICATIONS: Lead (Pb) and cadmium (Cd) are prevalent heavy metal pollutants in soil, and pose a major threat to crop production, food security and human health. Here, MeRIP-seq approach was employed to analyze the regulatory network activated in soybean nodules under Pb and Cd stress, ultimately leading to the identification of 19 shared differentially expressed DMGs. When overexpressed, GmATM1;1 was found to enhance the Pb and Cd tolerance of soybean nodules. These results provide a theoretical basis for studies on tolerance to heavy metals in symbiotic nitrogen fixation, and provide an approach to enhancing Pb and Cd tolerance in soybean production.

Exploring N6-methyladenosine (m6A) modification in tree species: opportunities and challenges

N 6-methyladenosine (m6A) in eukaryotes is the most common and widespread internal modification in mRNA. The modification regulates mRNA stability, translation efficiency, and splicing, thereby fine-tuning gene regulation. In plants, m6A is dynamic and critical for various growth stages, embryonic development, morphogenesis, flowering, stress response, crop yield, and biomass. Although recent high-throughput sequencing approaches have enabled the rapid identification of m6A modification sites, the site-specific mechanism of this modification remains unclear in trees. In this review, we discuss the functional significance of m6A in trees under different stress conditions and discuss recent advancements in the quantification of m6A. Quantitative and functional insights into the dynamic aspect of m6A modification could assist researchers in engineering tree crops for better productivity and resistance to various stress conditions.

Molecular-Assisted Breeding of Cadmium Pollution-Safe Cultivars

Cadmium (Cd) contamination in edible agricultural products, especially in crops intended for consumption, has raised worldwide concerns regarding food safety. Breeding of Cd pollution-safe cultivars (Cd-PSCs) is an effective solution to preventing the entry of Cd into the food chain from contaminated agricultural soil. Molecular-assisted breeding methods, based on molecular mechanisms for cultivar-dependent Cd accumulation and bioinformatic tools, have been developed to accelerate and facilitate the breeding of Cd-PSCs. This review summarizes the recent progress in the research of the low Cd accumulation traits of Cd-PSCs in different crops. Furthermore, the application of molecular-assisted breeding methods, including transgenic approaches, genome editing, marker-assisted selection, whole genome-wide association analysis, and transcriptome, has been highlighted to outline the breeding of Cd-PSCs by identifying critical genes and molecular biomarkers. This review provides a comprehensive overview of the development of Cd-PSCs and the potential future for breeding Cd-PSC using modern molecular technologies.

Coordination of m6A mRNA Methylation and Gene Transcriptome in Sugarcane Response to Drought Stress

The N6-methyladenosine (m6A) methylation of mRNA is involved in biological processes essential for plant growth. To explore the m6A modification of sugarcane and reveal its regulatory function, methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to construct the m6A map of sugarcane. In this study, m6A sites of sugarcane transcriptome were significantly enriched around the stop codon and within 3'-untranslated regions (3'UTR). Gene ontology (GO) analysis showed that the m6A modification genes are associated with metabolic biosynthesis. In addition, the m6A modification of drought-resistant transcript mRNA increased significantly under drought (DR) treatment, resulting in enhanced mRNA stability, which is involved in regulating sugarcane drought resistance. GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results showed that differentially methylated peak (DMP) modification of differentially expressed genes (DEGs) in DR were particularly associated with abscisic acid (ABA) biosynthesis. The upregulated genes were significantly enriched in the ABA metabolism, ethylene response, fatty acid metabolism, and negative regulation of the abscisic acid activation signaling pathway. These findings provide a basis and resource for sugarcane RNA epigenetic studies and further increase our knowledge of the functions of m6A modifications in RNA under abiotic stress.

The emerging role of epitranscriptome in shaping stress responses in plants

KEY MESSAGE

RNA modifications and editing changes constitute 'epitranscriptome' and are crucial in regulating the development and stress response in plants. Exploration of the epitranscriptome and associated machinery would facilitate the engineering of stress tolerance in crops. RNA editing and modifications post-transcriptionally decorate almost all classes of cellular RNAs, including tRNAs, rRNAs, snRNAs, lncRNAs and mRNAs, with more than 170 known modifications, among which m6A, Ψ, m5C, 8-OHG and C-to-U editing are the most abundant. Together, these modifications constitute the "epitranscriptome", and contribute to changes in several RNA attributes, thus providing an additional structural and functional diversification to the "cellular messages" and adding another layer of gene regulation in organisms, including plants. Numerous evidences suggest that RNA modifications have a widespread impact on plant development as well as in regulating the response of plants to abiotic and biotic stresses. High-throughput sequencing studies demonstrate that the landscapes of m6A, m5C, Am, Cm, C-to-U, U-to-G, and A-to-I editing are remarkably dynamic during stress conditions in plants. GO analysis of transcripts enriched in Ψ, m6A and m5C modifications have identified bonafide components of stress regulatory pathways. Furthermore, significant alterations in the expression pattern of genes encoding writers, readers, and erasers of certain modifications have been documented when plants are grown in challenging environments. Notably, manipulating the expression levels of a few components of RNA editing machinery markedly influenced the stress tolerance in plants. We provide updated information on the current understanding on the contribution of RNA modifications in shaping the stress responses in plants. Unraveling of the epitranscriptome has opened new avenues for designing crops with enhanced productivity and stress resilience in view of global climate change.

METTL3-mediated m6A mRNA modification was involved in cadmium-induced liver injury

Cadmium is an environmental pollutant that has extensive deleterious effects. However, the mechanisms underlying the hepatotoxicity induced by long-term exposure to cadmium remained undefined. In the present study, we explored the role of m6A methylation in the development of cadmium-induced liver disease. We showed a dynamic change of RNA methylation in liver tissue from mice administrated with cadmium chloride (CdCl2) for 3, 6 and 9 months, respectively. Particularly, the METTL3 expression was declined in a time-dependent manner, associated with the degree of liver injury, indicating the involvement of METTL3 in hepatotoxicity induced by CdCl2. Moreover, we established a mouse model with liver-specific over-expression of Mettl3 and administrated these mice with CdCl2 for 6 months. Notably, METTL3 highly expressed in hepatocytes attenuated CdCl2-induced steatosis and liver fibrosis in mice. In vitro assay also showed METTL3 overexpression ameliorated the CdCl2-induced cytotoxicity and activation of primary hepatic stellate cells. Furthermore, transcriptome analysis identified 268 differentially expressed genes both in mice liver tissue treated with CdCl2 for 3 months and 9 months. Among them, 115 genes were predicted to be regulated by METTL3 determined by m6A2Target database. Further analysis revealed the perturbation of metabolic pathway, glycerophospholipid metabolism, ErbB signaling pathway, Hippo signaling pathway, and choline metabolism in cancer, and circadian rhythm, led to hepatotoxicity induced by CdCl2. Collectively, our findings reveal new insight into the crucial role of epigenetic modifications in hepatic diseases caused by long-term exposure to cadmium.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

An analysis of differentially expressed and differentially m6A-modified transcripts in soybean roots treated with lead

Lead is one of the most common toxic heavy metal pollutants in nature, and exposure to lead can cause serious toxicity to many organisms. In this study, we collected root growth data from soybean plants exposed to lead for seven days and confirmed that lead significantly inhibited root growth. We performed a transcriptome-wide m6A methylation analysis to study the response of soybean RNA methylation groups to lead. The m6A modified regions were enriched near the 3'UTR region and stop codon, and m6A methylation was positively correlated with transcript abundance. In the presence of lead, the transcriptome range of m6A RNA methylation peaks increased, and we identified 1144 m6A modification peaks and 1094 differentially expressed genes. The integration of m6A methylation and transcriptomic results enabled us to identify 16 candidate genes whose transcripts were differentially methylated and differentially expressed under lead stress. Annotation results suggest that these genes may promote abiotic stress tolerance by impacting lead uptake, transport, and accumulation through ROS pathways, enzymes, transporters, and hormones. These results provide candidate genes for future studies of lead stress tolerance mechanisms in soybean roots and provide genetic resources for studying plant heavy metal stress in soybean breeding.

N6-methyladenosine RNA modification regulates cotton drought response in a Ca2+ and ABA-dependent manner

N⁶ -methyladenosine (m⁶ A) is the most prevalent internal modification present in mRNAs, and is considered to participate in a range of developmental and biological processes. Drought response is highly regulated at the genomic, transcriptional and post-transcriptional levels. However, the biological function and regulatory mechanism of m⁶ A modification in the drought stress response is still poorly understood. We generated a transcriptome-wide m⁶ A map using drought-resistant and drought-sensitive varieties of cotton under different water deficient conditions to uncover patterns of m⁶ A methylation in cotton response to drought stress. The results reveal that m⁶ A represents a common modification and exhibit dramatic changes in distribution during drought stress. More 5'UTR m⁶ A was deposited in the drought-resistant variety and was associated with a positive effect on drought resistance by regulating mRNA abundance. Interestingly, we observed that increased m⁶ A abundance was associated with increased mRNA abundance under drought, contributing to drought resistance, and vice versa. The demethylase GhALKBH10B was found to decrease m⁶ A levels, facilitating the mRNA decay of ABA signal-related genes (GhZEP, GhNCED4 and GhPP2CA) and Ca²⁺ signal-related genes (GhECA1, GhCNGC4, GhANN1 and GhCML13), and mutation of GhALKBH10B enhanced drought resistance at seedling stage in cotton. Virus-induced gene silencing (VIGS) of two Ca²⁺ -related genes, GhECA1 and GhCNGC4, reduced drought resistance with the decreased m⁶ A enrichment on silenced genes in cotton. Collectively, we reveal a novel mechanism of post-transcriptional modification involved in affecting drought response in cotton, by mediating m⁶ A methylation on targeted transcripts in the ABA and Ca²⁺ signalling transduction pathways.

The covalent nucleotide modifications within plant mRNAs: What we know, how we find them, and what should be done in the future

Although covalent nucleotide modifications were first identified on the bases of transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs), a number of these epitranscriptome marks have also been found to occur on the bases of messenger RNAs (mRNAs). These covalent mRNA features have been demonstrated to have various and significant effects on the processing (e.g. splicing, polyadenylation, etc.) and functionality (e.g. translation, transport, etc.) of these protein-encoding molecules. Here, we focus our attention on the current understanding of the collection of covalent nucleotide modifications known to occur on mRNAs in plants, how they are detected and studied, and the most outstanding future questions of each of these important epitranscriptomic regulatory signals.

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.

Changes in the m6A RNA methylome accompany the promotion of soybean root growth by rhizobia under cadmium stress

Cadmium (Cd) is the most widely distributed heavy metal pollutant in soil and has significant negative effects on crop yields and human health. Rhizobia can enhance soybean growth in the presence of heavy metals, and the legume-rhizobia symbiosis has been used to promote heavy-metal phytoremediation, but much remains to be learned about the molecular networks that underlie these effects. Here, we demonstrated that soybean root growth was strongly suppressed after seven days of Cd exposure but that the presence of rhizobia largely eliminated this effect, even prior to nodule development. Moreover, rhizobia did not appear to promote root growth by limiting plant Cd uptake: seedlings with and without rhizobia had similar root Cd concentrations. Previous studies have demonstrated a role for m6A RNA methylation in the response of rice and barley to Cd stress. We therefore performed transcriptome-wide m6A methylation profiling to investigate changes in the soybean RNA methylome in response to Cd with and without rhizobia. Here, we provide some of the first data on transcriptome-wide m6a RNA methylation patterns in soybean; m6A modifications were concentrated at the 3' UTR of transcripts and showed a positive relationship with transcript abundance. Transcriptome-wide m6A RNA methylation peaks increased in the presence of Cd, and the integration of m6A methylome and transcriptome results enabled us to identify 154 genes whose transcripts were both differentially methylated and differentially expressed in response to Cd stress. Annotation results suggested that these genes were associated with Ca²⁺ homeostasis, ROS pathways, polyamine metabolism, MAPK signaling, hormones, and biotic stress responses. There were 176 differentially methylated and expressed transcripts under Cd stress in the presence of rhizobia. In contrast to the Cd-only gene set, they were also enriched in genes related to auxin, jasmonic acid, and brassinosteroids, as well as abiotic stress tolerance. They contained fewer genes related to Ca²⁺ homeostasis and also included candidates with known functions in the legume-rhizobia symbiosis. These findings offer new insights into how rhizobia promote soybean root growth under Cd stress; they provide candidate genes for research on plant heavy metal responses and for the use of legumes in phytoremediation.

Transcriptome-wide profiling of mRNA N6-methyladenosine modification in rice panicles and flag leaves

N⁶-methyladenosine (m⁶A) modification is essential for plant growth and development. Exploring m⁶A methylation patterns in rice tissues is fundamental to understanding the regulatory effects of this modification. Here, we profiled the transcriptome-wide m⁶A landscapes of rice panicles at the booting stage (PB) and flowering stage (PF), and of flag leaves at the flowering stage (LF). The global m⁶A level differed significantly among the three tissues and was closely associated with the expression of writer and eraser genes. The methylated gene ratio was higher in the flag leaves than in the panicles. Compared with commonly methylated genes, tissue-specific methylated genes showed lower levels of both m⁶A modification and expression, and a preference for m⁶A deposition in the coding sequence region. The m⁶A profiles of the two organs had more distinct differences than the profiles of the same organ at different stages. A negative correlation between m⁶A levels and gene expression was observed in PF vs. PB but not in PF vs. LF, indicting the complicated regulatory effect of m⁶A on gene expression. The distinct expression patterns of m⁶A reader genes in different tissues indicate that readers may affect gene stability through binding. Overall, our findings demonstrated that m⁶A modification influences tissue function by regulating gene expression. Our findings provide valuable insights on the regulation and biological functions of m⁶A modifications in rice.

N6-methyladenosine mRNA methylation is important for the light response in soybean

N6-methyladenosine (m⁶A) modification of messenger RNA (mRNA) is the most prevalent and abundant modification in eukaryotic mRNA and posttranscriptionally modulates the transcriptome at almost all stages of mRNA metabolism. In plants, m⁶A is crucial for embryonic-phase growth, flowering time control, microspore generation and fruit maturation. However, the role of m⁶A in plant responses to light, the most important environmental stimulus, remains unexplored. Here, we profile the m⁶A transcriptome of Williams 82, a soybean cultivar, and reveal that m⁶A is highly conserved and plays an important role in the response to light stimuli in soybean. Similar to the case in Arabidopsis, m⁶A in soybean is enriched not only around the stop codon and within the 3'UTR but also around the start codon. Moreover, genes with methylation occurring in the 3'UTR have higher expression levels and are more prone to alternative splicing. The core genes in the light signaling pathway, GmSPA1a, GmPRR5e and GmBIC2b, undergo changes in methylation modification and transcription levels in response to light. KEGG pathway analysis revealed that differentially expressed genes with differential m⁶A peaks were involved in the "photosynthesis" and "circadian rhythm" pathways. Our results highlight the important role played by epitranscriptomic mRNA methylation in the light response in soybean and provide a solid basis for determining the functional role of light on RNA m⁶A modification in this plant.

Epitranscriptomic mRNA modifications governing plant stress responses: underlying mechanism and potential application

Plants inevitably encounter environmental adversities, including abiotic and biotic stresses, which significantly impede plant growth and reduce crop yield. Thus, fine-tuning the fate and function of stress-responsive RNAs is indispensable for plant survival under such adverse conditions. Recently, post-transcriptional RNA modifications have been studied as a potent route to regulate plant gene expression under stress. Among over 160 mRNA modifications identified to date, N⁶ -methyladenosine (m⁶ A) in mRNAs is notable because of its multifaceted roles in plant development and stress response. Recent transcriptome-wide mapping has revealed the distribution and patterns of m⁶ A in diverse stress-responsive mRNAs in plants, building a foundation for elucidating the molecular link between m⁶ A and stress response. Moreover, the identification and characterization of m⁶ A writers, readers and erasers in Arabidopsis and other model crops have offered insights into the biological roles of m⁶ A in plant abiotic stress responses. Here, we review the recent progress of research on mRNA modifications, particularly m⁶ A, and their dynamics, distribution, regulation and biological functions in plant stress responses. Further, we posit potential strategies for breeding stress-tolerant crops by engineering mRNA modifications and propose the future direction of research on RNA modifications to gain a much deeper understanding of plant stress biology.

Mechanism of macroalgae Gracilaria bailiniae responding to cadmium and lanthanum

Macroalgae can accumulate a wide array of metals, leading to their appliance as biomonitors of aquatic environments. With the rapid development of industrial and agricultural-based activities, Cd pollution in aquatic environments is considered an increasingly severe problem worldwide. Although La could alleviate the Cd stress in higher terrestrial plants, the response mechanisms of macroalgae to Cd and La are unknown. Along these lines, in this work, Cd significantly affected the growth, internal cellular structure, photosynthesis, pigment content, antioxidant enzyme activity, and lipid peroxidation level of G. bailiniae. However, the presence of La alleviated these adverse effects from Cd. Furthermore, the response mechanism of G. bailiniae to Cd was attributed to the self-antioxidant ability enhancement, membrane defense, and programmed-cellular regulation. However, the presence of La mediated the biosynthesis of both flavonoids and lipids, which inhibited the Cd accumulation, modulated algal stress signalling networks, renewed the impaired chlorophyll molecule, maintained the activity of the crucial enzyme, enhanced antioxidant ability, and maintained the stabilization of redox homeostasis, alleviating the adverse impact from Cd and improve the growth of G. bailiniae. The experimental results successfully demonstrate a new detoxicant to alleviate Cd stress, promoting a more comprehensive array of macroalgal applications.

Transcriptome-wide analyses of RNA m6A methylation in hexaploid wheat reveal its roles in mRNA translation regulation

N6-methyladenosine (m6A) is the most abundant RNA modification in eukaryotic messenger RNAs. m6A was discovered in wheat about 40 years ago; however, its potential roles in wheat remain unknown. In this study, we profiled m6As in spikelets transcriptome at the flowering stage of hexaploid wheat and found that m6As are evenly distributed across the A, B, and D subgenomes but their extents and locations vary across homeologous genes. m6As are enriched in homeologous genes with close expression levels and the m6A methylated genes are more conserved. The extent of m6A methylation is negatively correlated with mRNA expression levels and its presence on mRNAs has profound impacts on mRNA translation in a location-dependent manner. Specifically, m6As within coding sequences and 3'UTRs repress the translation of mRNAs while the m6As within 5'UTRs and start codons could promote it. The m6A-containing mRNAs are significantly enriched in processes and pathways of "translation" and "RNA transport," suggesting the potential role of m6As in regulating the translation of genes involved in translation regulation. Our data also show a stronger translation inhibition by small RNAs (miRNA and phasiRNA) than by m6A methylation, and no synergistical effect between the two was observed. We propose a secondary amplification machinery of translation regulation triggered by the changes in m6A methylation status. Taken together, our results suggest translation regulation as a key role played by m6As in hexaploid wheat.

m6 A-mediated regulation of crop development and stress responses

Dynamic chemical modifications in eukaryotic messenger RNAs (mRNAs) constitute an essential layer of gene regulation, among which N⁶ -methyladenosine (m⁶ A) was unveiled to be the most abundant. m⁶ A functionally modulates important biological processes in various mammals and plants through the regulation of mRNA metabolism, mainly mRNA degradation and translation efficiency. Physiological functions of m⁶ A methylation are diversified and affected by intricate sequence contexts and m⁶ A machineries. A number of studies have dissected the functional roles and the underlying mechanisms of m⁶ A modifications in regulating plant development and stress responses. Recently, it was demonstrated that the human FTO-mediated plant m⁶ A removal caused dramatic yield increases in rice and potato, indicating that modulation of m⁶ A methylation could be an efficient strategy for crop improvement. In this review, we summarize the current progress concerning the m⁶ A-mediated regulation of crop development and stress responses, and provide an outlook on the potential application of m⁶ A epitranscriptome in the future improvement of crops.

Exploring epitranscriptomics for crop improvement and environmental stress tolerance

Climate change and stressful environmental conditions severely hamper crop growth, development and yield. Plants respond to environmental perturbations, through their plasticity provided by key-genes, governed at post-/transcriptional levels. Gene-regulation in plants is a multilevel process controlled by diverse cellular entities that includes transcription factors (TF), epigenetic regulators and non-coding RNAs beside others. There are successful studies confirming the role of epigenetic modifications (DNA-methylation/histone-modifications) in gene expression. Recent years have witnessed emergence of a highly specialized field the "Epitranscriptomics". Epitranscriptomics deals with investigating post-transcriptional RNA chemical-modifications present across the life forms that change structural, functional and biological characters of RNA. However, deeper insights on of epitranscriptomic modifications, with >140 types known so far, are to be understood fully. Researchers have identified epitranscriptome marks (writers, erasers and readers) and mapped the site-specific RNA modifications (m6A, m⁵C, 3' uridylation, etc.) responsible for fine-tuning gene expression in plants. Simultaneous advancement in sequencing platforms, upgraded bioinformatic tools and pipelines along with conventional labelled techniques have further given a statistical picture of these epitranscriptomic modifications leading to their potential applicability in crop improvement and developing climate-smart crops. We present herein the insights on epitranscriptomic machinery in plants and how epitranscriptome and epitranscriptomic modifications underlying plant growth, development and environmental stress responses/adaptations. Third-generation sequencing technology, advanced bioinformatics tools and databases being used in plant epitranscriptomics are also discussed. Emphasis is given on potential exploration of epitranscriptome engineering for crop-improvement and developing environmental stress tolerant plants covering current status, challenges and future directions.

Epitranscriptomics: An Additional Regulatory Layer in Plants' Development and Stress Response

Epitranscriptomics has added a new layer of regulatory machinery to eukaryotes, and the advancement of sequencing technology has revealed more than 170 post-transcriptional modifications in various types of RNAs, including messenger RNA (mRNA), transfer RNA (tRNA), ribosomal RNA (rRNA), and long non-coding RNA (lncRNA). Among these, N6-methyladenosine (m6A) and N5-methylcytidine (m5C) are the most prevalent internal mRNA modifications. These regulate various aspects of RNA metabolism, mainly mRNA degradation and translation. Recent advances have shown that regulation of RNA fate mediated by these epitranscriptomic marks has pervasive effects on a plant's development and responses to various biotic and abiotic stresses. Recently, it was demonstrated that the removal of human-FTO-mediated m6A from transcripts in transgenic rice and potatoes caused a dramatic increase in their yield, and that the m6A reader protein mediates stress responses in wheat and apple, indicating that regulation of m6A levels could be an efficient strategy for crop improvement. However, changing the overall m6A levels might have unpredictable effects; therefore, the identification of precise m6A levels at a single-base resolution is essential. In this review, we emphasize the roles of epitranscriptomic modifications in modulating molecular, physiological, and stress responses in plants, and provide an outlook on epitranscriptome engineering as a promising tool to ensure food security by editing specific m6A and m5C sites through robust genome-editing technology.

Transcriptome-wide m6A methylation profile reveals regulatory networks in roots of barley under cadmium stress

Cadmium (Cd) pollutants restrict crop yield and food security in long-term agricultural activities. Crops have evolved adaptive strategies under Cd condition, however, the transcriptional regulatory mechanism of Cd-tolerant genes remains to be largely illustrated. In this study, barley roots were exposed to 5 µM CdCl₂ for physiological response and transcriptome-wide m6A methylation profile. Cd stress inhibited root growth after 7 d Cd treatment, which is mainly associated with inhibited absorption of Mn. After Cd treatment, 8151 significantly modified m6A sites and 3920 differentially expressed genes were identified. Transcriptome-wide m6A hypermethylation was widely induced by Cd stress and enriched near the stop codon and 3' UTR regions. Among 435 m6A modified DEGs, 319 hypermethylated genes were up-regulated and 84 hypomethylated genes were down-regulated, respectively, indicating a positive correlation of m6A methylation and expression. But well-known Cd transporter genes (HvNramp5, HvIRT1, HvHMA3, etc.) were not modified by m6A methylation, except for ABC transporters. We further found key Cd-responding regulatory genes were positively modulated with m6A methylation, including MAPK, WRKY and MYB members. This study proposed a transcriptional regulatory network of Cd stress response in barley roots, which may provide new insight into gene manipulation of controlling low Cd accumulation for crops.