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

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.

Visual observation of polystyrene nano-plastics in grape seedlings of Thompson Seedless and assessing their effects via transcriptomics and metabolomics

Micro/nano-plastics (MNPs) are emerging non-point source pollutants that have garnered increasing attention owing to their threat to ecosystems. Studies on the effects of MNPs on horticultural crops are scarce. Specifically, whether MNPs can be absorbed and transported by grapevines have not been reported. To fill this gap, we added polystyrene nanoplastics (PS-NPs, 100 nm) to a hydroponic environment and observed their distribution in grape seedlings of Thompson Seedless (TS, Vitis vinifera L.). After 15 d of exposure, plastic nanospheres were detected on the cell walls of the roots, stems, and leaves using confocal microscopy and scanning electron microscopy. This indicated that PS-NPs can also be absorbed by the root system through the epidermis-cortex interface in grapevines and transported upward along the xylem conduit. Furthermore, we analyzed the molecular response mechanisms of TS grapes to the PS-NPs. Through the measurement of relevant indicators and combined omics analysis, we found that plant hormone signal transduction, flavonoid and flavonol biosynthesis, phenylpropanoid biosynthesis, and MAPK signaling pathway biosynthesis played crucial roles in its response to PS-NPs. The results not only revealed the potential risk of MNPs being absorbed by grapevines and eventually entering the food chain but also provided valuable scientific evidence and data for the assessment of plant health and ecological risk.

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.

Transcriptome-wide m6A methylation profile reveals its potential role underlying drought response in wheat (Triticum aestivum L.)

MAIN CONCLUSION

This study revealed the transcriptome-wide m6A methylation profile under drought stress and found that TaETC9 might regulate drought tolerance through mediating RNA methylation in wheat. Drought is one of the most destructive environmental constraints limiting crop growth and development. N6-methyladenosine (m6A) is a prevalent and important post-transcriptional modification in various eukaryotic RNA molecules, playing the crucial role in regulating drought response in plants. However, the significance of m6A in wheat (Triticum aestivum L.), particularly its involvment in drought response, remains underexplored. In this study, we investigated the transcriptome-wide m6A profile under drought stress using parallel m6A immunoprecipitation sequencing (MeRIP-seq). Totally, 4221 m6A peaks in 3733 m6A-modified genes were obtained, of which 373 methylated peaks exhibited differential expression between the control (CK) and drought-stressed treatments. These m6A loci were significantly enriched in proximity to stop codons and within the 3'-untranslated region. Integration of MeRIP-seq and RNA-seq revealed a positive correlation between m6A methylation and mRNA abundance and the genes displaying both differential methylation and expression were obtained. Finally, qRT-PCR analyses were further performed and the results found that the m6A-binding protein (TaETC9) showed significant up-regulation, while the m6A demethylase (TaALKBH10B) was significantly down-regulated under drought stress, contributing to increased m6A levels. Furthermore, the loss-of-function mutant of TaECT9 displayed significantly higher drought sensitivity compared to the wild type, highlighting its role in regulating drought tolerance. This study reported the first wheat m6A profile associated with drought stress, laying the groundwork for unraveling the potential role of RNA methylation in drought responses and enhancing stress tolerance in wheat through epigenetic approaches.

N6-Methyladenosine (m6A) Sequencing Reveals Heterodera glycines-Induced Dynamic Methylation Promoting Soybean Defense

Unraveling the intricacies of soybean cyst nematode (Heterodera glycines) race 4 resistance and susceptibility in soybean breeding lines-11-452 (highly resistant) and Dongsheng1 (DS1, highly susceptible)-was the focal point of this study. Employing cutting-edge N6-methyladenosine (m6A) and RNA sequencing techniques, we delved into the impact of m6A modification on gene expression and plant defense responses. Through the evaluation of nematode development in both resistant and susceptible roots, a pivotal time point (3 days postinoculation) for m6A methylation sequencing was identified. Our sequencing data exhibited robust statistics, successful soybean genome mapping, and prevalent m6A peak distributions, primarily in the 3' untranslated region and stop codon regions. Analysis of differential methylation peaks and differentially expressed genes revealed distinctive patterns between resistant and susceptible genotypes. In the highly resistant line (11-452), key resistance and defense-associated genes displayed increased expression coupled with inhibited methylation, encompassing crucial players such as R genes, receptor kinases, and transcription factors. Conversely, the highly susceptible DS1 line exhibited heightened expression correlated with decreased methylation in genes linked to susceptibility pathways, including Mildew Locus O-like proteins and regulatory elements affecting defense mechanisms. Genome-wide assessments, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, and differential methylation peak/differentially expressed gene overlap emphasized the intricate interplay of m6A modifications, alternative splicing, microRNA, and gene regulation in plant defense. Protein-protein interaction networks illuminated defense-pivotal genes, delineating divergent mechanisms in resistant and susceptible responses. This study sheds light on the dynamic correlation between methylation, splicing, and gene expression, providing profound insights into plant responses to nematode infection.

Cadmium exacerbates liver injury by remodeling ceramide metabolism: Multiomics and laboratory evidence

Cadmium (Cd) is a toxic heavy metal that primarily targets the liver. Cd exposure disrupts specific lipid metabolic pathways; however, the underlying mechanisms remain unclear. This study aimed to investigate the lipidomic characteristics of rat livers after Cd exposure as well as the potential mechanisms of Cd-induced liver injury. Our analysis of established Cd-exposed rat and cell models showed that Cd exposure resulted in liver lipid deposition and hepatocyte damage. Lipidomic detection, transcriptome sequencing, and experimental analyses revealed that Cd mainly affects the sphingolipid metabolic pathway and that the changes in ceramide metabolism are the most significant. In vitro experiments revealed that the inhibition of ceramide synthetase activity or activation of ceramide decomposing enzymes ameliorated the proapoptotic and pro-oxidative stress effects of Cd, thereby alleviating liver injury. In contrast, the exogenous addition of ceramide aggravated liver injury. In summary, Cd increased ceramide levels by remodeling ceramide synthesis and catabolism, thereby promoting hepatocyte apoptosis and oxidative stress and ultimately aggravating liver injury. Reducing ceramide levels can serve as a potential protective strategy to mitigate the liver toxicity of Cd. This study provides new evidence for understanding Cd-induced liver injury at the lipidomic level and insights into the health risks and toxicological mechanisms associated with Cd.

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.

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.

Genome-wide identification of the N6-methyladenosine regulatory genes reveals NtFIP37B increases drought resistance of tobacco (Nicotiana tabacum L.)

BACKGROUND

N6-methyladenosine (m6A) is one of the common internal RNA modifications found in eukaryotes. The m6A modification can regulate various biological processes in organisms through the modulation of alternative splicing, alternative polyadenylation, folding, translation, localization, transport, and decay of multiple types of RNA, without altering the nucleotide sequence. The three components involved in m6A modification, namely writer, eraser, and reader, mediate the abundance of RNA m6A modification through complex collaborative actions. Currently, research on m6A regulatory genes in plants is still in its infancy.

RESULTS

In this study, we identified 52 candidate m6A regulatory genes in common tobacco (Nicotiana tabacum L.). Gene structure, conserved domains, and motif analysis showed structural and functional diversity among different subgroups of tobacco m6A regulatory genes. The amplification of m6A regulatory genes were mainly driven by polyploidization and dispersed duplication, and duplicated genes evolved through purified selection. Based on the potential regulatory network and expression pattern analysis of m6A regulatory genes, a significant number of m6A regulatory genes might play important roles in growth, development, and stress response processes. Furthermore, we have confirmed the critical role of NtFIP37B, an m6A writer gene in tobacco, in enhancing drought resistance.

CONCLUSIONS

This study provides useful information for better understanding the evolution of m6A regulatory genes and the role of m6A modification in tobacco stress response, and lays the foundation for further elucidating the function of m6A regulatory genes in tobacco.

Role of RNA Modifications, Especially m6A, in Aflatoxin Biosynthesis of Aspergillus flavus

RNA modifications play key roles in eukaryotes, but the functions in Aspergillus flavus are still unknown. Temperature has been reported previously to be a critical environmental factor that regulates the aflatoxin production of A. flavus, but much remains to be learned about the molecular networks. Here, we demonstrated that 12 kinds of RNA modifications in A. flavus were significantly changed under 29 °C compared to 37 °C incubation; among them, m6A was further verified by a colorimetric method. Then, the transcriptome-wide m6A methylome and m6A-altered genes were comprehensively illuminated through methylated RNA immunoprecipitation sequencing and RNA sequencing, from which 22 differentially methylated and expressed transcripts under 29 °C were screened out. It is especially notable that AFCA_009549, an aflatoxin biosynthetic pathway gene (aflQ), and the m6A methylation of its 332nd adenine in the mRNA significantly affect aflatoxin biosynthesis in A. flavus both on media and crop kernels. The content of sterigmatocystin in both ΔaflQ and aflQA332C strains was significantly higher than that in the WT strain. Together, these findings reveal that RNA modifications are associated with secondary metabolite biosynthesis of A. flavus.

Comprehensive analysis of aflatoxin B1 biosynthesis in Aspergillus flavus via transcriptome-wide m6A methylome response to cycloleucine

Aspergillus flavus and its toxic aflatoxins secondary metabolites contaminate food and grains, posing a severe threat to human health and leading to liver cancer. Here, we demonstrated that cycloleucine blocked aflatoxin B1 synthesis by inhibiting N6-methyladenosine (m6A) methylation modification of messenger RNA (mRNA). m6A Methylation Immunoprecipitation Sequencing (m6A MeRIP-Seq)-based comprehensive transcriptome-wide m6A profiling identified 102 differentially expressed genes that underwent m6A modification, of which 22 hypermethylated genes were downregulated and 49 hypomethylated genes were upregulated, suggesting a negative correlation between m6A methylation and gene expression. Notably, cycloleucine inhibited aflatoxin B1 production via multiple targets. The m6A sites of several key genes involved in the aflatoxin B1 biosynthesis pathway were significantly enriched in the coding sequence and around the stop codon, resulting in their downregulation. Furthermore, m6A methylation on genes related to the aflatoxin B1 biosynthesis pathway led to reduced mRNA stability. Cycloleucine inhibition of aflatoxin B1 production highlights its potential as an agent for removing mycotoxins in environmental pollution. ENVIRONMENTAL IMPLICATION: Aflatoxins, highly carcinogenic secondary metabolites produced by Aspergillus flavus, frequently contaminate crops such as peanut, corn, wheat and sesame leading to irreversible loss in the quality and yield of agricultural products and posing serious threats to food safety. Aflatoxins has also been linked to developmental delays and liver cancer in humans. In our study, 'monitoring aflatoxin concentrations and its bioaccumulation in organisms' has been conducted. The results demonstrated that aflatoxin production in A. flavus was completely blocked after cycloleucine treatment. Additionally, we demonstrated that inhibition of aflatoxin was linked to N6-methyladenosine methylation of multiple genes in aflatoxin biosynthesis pathway.

The role of epigenetic and epitranscriptomic modifications in plants exposed to non-essential metals

Contamination of the soil with non-essential metals and metalloids is a serious problem in many regions of the world. These non-essential metals and metalloids are toxic to all organisms impacting crop yields and human health. Crop plants exposed to high concentrations of these metals leads to perturbed mineral homeostasis, decreased photosynthesis efficiency, inhibited cell division, oxidative stress, genotoxic effects and subsequently hampered growth. Plants can activate epigenetic and epitranscriptomic mechanisms to maintain cellular and organism homeostasis. Epigenetic modifications include changes in the patterns of cytosine and adenine DNA base modifications, changes in cellular non-coding RNAs, and remodeling histone variants and covalent histone tail modifications. Some of these epigenetic changes have been shown to be long-lasting and may therefore contribute to stress memory and modulated stress tolerance in the progeny. In the emerging field of epitranscriptomics, defined as chemical, covalent modifications of ribonucleotides in cellular transcripts, epitranscriptomic modifications are postulated as more rapid modulators of gene expression. Although significant progress has been made in understanding the plant's epigenetic changes in response to biotic and abiotic stresses, a comprehensive review of the plant's epigenetic responses to metals is lacking. While the role of epitranscriptomics during plant developmental processes and stress responses are emerging, epitranscriptomic modifications in response to metals has not been reviewed. This article describes the impact of non-essential metals and metalloids (Cd, Pb, Hg, Al and As) on global and site-specific DNA methylation, histone tail modifications and epitranscriptomic modifications in plants.

Role of Epigenetic Factors in Response to Stress and Establishment of Somatic Memory of Stress Exposure in Plants

All species are well adapted to their environment. Stress causes a magnitude of biochemical and molecular responses in plants, leading to physiological or pathological changes. The response to various stresses is genetically predetermined, but is also controlled on the epigenetic level. Most plants are adapted to their environments through generations of exposure to all elements. Many plant species have the capacity to acclimate or adapt to certain stresses using the mechanism of priming. In most cases, priming is a somatic response allowing plants to deal with the same or similar stress more efficiently, with fewer resources diverted from growth and development. Priming likely relies on multiple mechanisms, but the differential expression of non-coding RNAs, changes in DNA methylation, histone modifications, and nucleosome repositioning play a crucial role. Specifically, we emphasize the role of BRM/CHR17, BRU1, FGT1, HFSA2, and H2A.Z proteins as positive regulators, and CAF-1, MOM1, DDM1, and SGS3 as potential negative regulators of somatic stress memory. In this review, we will discuss the role of epigenetic factors in response to stress, priming, and the somatic memory of stress exposures.

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.

Genome-Wide Identification, Expression and Evolution Analysis of m6A Writers, Readers and Erasers in Aegilops_tauschii

N6-methyladenosine modifications (m6A) is one of the most abundant and prevalent post-transcriptional RNA modifications in plants, playing the crucial role in plant growth and development and stress adaptation. However, the m6A regulatory machinery in Aegilops_tauschii, the D genome progenitor of common wheat, is not well understood at present. Here, we systematically identified the m6A-related genes in Aegilops with a genome-wide search approach. In total, 25 putative m6A genes composed of 5 writers, 13 readers and 7 erasers were obtained. A phylogenetic analysis clearly grouped them into three subfamilies with the same subfamily showing similar gene structures and conserved domains. These m6A genes were found to contain a large number of cis-acting elements associating with plant hormones, regulation of growth and development as well as stress response, suggesting their widespread regulation function. Furthermore, the expression profiling of them was investigated using RNA-seq data to obtain stress-responsive candidates, of which 5 were further validated with a qPCR analysis. Finally, the genetic variation of m6A-related genes was investigated between Aegilops and D subgenome of wheat based on re-sequencing data, and an obvious genetic bottleneck occurred on them during the wheat domestication process. The promising haplotype association with domestication and agronomic traits was also detected. This study provided some insights on the genomic organization and evolutionary features of m6A-related genes in Aegilops, which will facilitate the further functional study and also contribute to broaden the genetic basis for genetic improvement in wheat and other crops.

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.

Transcriptome-wide m6A methylation in natural yellow leaf of Catalpa fargesii

N6-methyladenosine (m6A) is the most abundant internal modification in eukaryotic messenger RNA, and involved in various biological processes in plants. However, the distribution features and functions of mRNA m6A methylation have been poorly explored in woody perennial plants. In this study, a new natural variety with yellow-green leaves, named Maiyuanjinqiu, was screened from the seedlings of Catalpa fargesii. Based on the preliminary experiment, the m6A methylation levels in the leaves of Maiyuanjinqiu were significantly higher than those in C. fargesii. Furthermore, a parallel analysis of m6A-seq and RNA-seq was carried out in different leaf color sectors. The result showed that m6A modification were mostly identified around the 3'-untranslated regions (3'-UTR), which was slightly negatively correlated with the mRNA abundance. KEGG and GO analyses showed that m6A methylation genes were associated with photosynthesis, pigments biosynthesis and metabolism, oxidation-reduction and response to stress, etc. The overall increase of m6A methylation levels in yellow-green leaves might be associated with the decreased the expression of RNA demethylase gene CfALKBH5. The silencing of CfALKBH5 caused a chlorotic phenotype and increased m6A methylation level, which further confirmed our hypothesis. Our results suggested that mRNA m6A methylation could be considered as a vital epigenomic mark and contribute to the natural variations in plants.

Investigation of the role of TmMYB16/123 and their targets (TmMTP1/11) in the tolerance of Taxus media to cadmium

The toxicity and stress caused by heavy metal contamination has become an important constraint to the growth and flourishing of trees. In particular, species belonging to the genus Taxus, which are the only natural source for the anti-tumor medicine paclitaxel, are known to be highly sensitive to environmental changes. To investigate the response of Taxus spp. to heavy metal stress, we analyzed the transcriptomic profiles of Taxus media trees exposed to cadmium (Cd2+). In total, six putative genes from the metal tolerance protein (MTP) family were identified in T. media, including two Cd2+ stress inducible TMP genes (TmMTP1, TmMTP11 and Taxus media). Secondary structure analyses predicted that TmMTP1 and TmMTP11, which are members of the Zn-CDF and Mn-CDF subfamily proteins, respectively, contained six and four classic transmembrane domains, respectively. The introduction of TmMTP1/11 into the ∆ycf1 yeast cadmium-sensitive mutant strain showed that TmMTP1/11 might regulate the accumulation of Cd2+ to yeast cells. To screen the upstream regulators, partial promoter sequences of the TmMTP1/11 genes were isolated using the chromosome walking method. Several myeloblastosis (MYB) recognition elements were identified in the promoters of these genes. Furthermore, two Cd2+-induced R2R3-MYB TFs, TmMYB16 and TmMYB123, were identified. Both in vitro and in vivo assays confirmed that TmMTB16/123 play a role in Cd2+ tolerance by activating and repressing the expression of TmMTP1/11 genes. The present study elucidated new regulatory mechanisms underlying the response to Cd stress and can contribute to the breeding of Taxus species with high environmental adaptability.

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.

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.

Transcriptome Analysis Reveals Differentially Expressed Genes Involved in Cadmium and Arsenic Accumulation in Tea Plant (Camellia sinensis)

Tea (Camellia sinensis) is the second most consumed drink in the world. Rapid industrialization has caused various impacts on nature and increased pollution by heavy metals. However, the molecular mechanisms of cadmium (Cd) and arsenic (As) tolerance and accumulation in tea plants are poorly understood. The present study focused on the effects of heavy metals Cd and As on tea plants. Transcriptomic regulation of tea roots after Cd and As exposure was analyzed to explore the candidate genes involved in Cd and As tolerance and accumulation. In total, 2087, 1029, 1707, and 366 differentially expressed genes (DEGs) were obtained in Cd1 (with Cd treatment for 10 days) vs. CK (without Cd treatment), Cd2 (with Cd treatment for 15 days) vs. CK, As1 (with As treatment for 10 days) vs. CK (without Cd treatment), and As2 (with As treatment for 15 days) vs. CK, respectively. Analysis of DEGs showed that a total of 45 DEGs with the same expression patterns were identified in four pairwise comparison groups. One ERF transcription factor (CSS0000647) and six structural genes (CSS0033791, CSS0050491, CSS0001107, CSS0019367, CSS0006162, and CSS0035212) were only increased at 15 d of Cd and As treatments. Using weighted gene co-expression network analysis (WGCNA) revealed that the transcription factor (CSS0000647) was positively correlated with five structural genes (CSS0001107, CSS0019367, CSS0006162, CSS0033791, and CSS0035212). Moreover, one gene (CSS0004428) was significantly upregulated in both Cd and As treatments, suggesting that these genes might play important roles in enhancing the tolerance to Cd and As stresses. These results provide candidate genes to enhance multi-metal tolerance through the genetic engineering technology.

Transcriptional Regulatory Network of Plant Cadmium Stress Response

Cadmium (Cd) is a non-essential heavy metal with high toxicity to plants. Plants have acquired specialized mechanisms to sense, transport, and detoxify Cd. Recent studies have identified many transporters involved in Cd uptake, transport, and detoxification. However, the complex transcriptional regulatory networks involved in Cd response remain to be elucidated. Here, we provide an overview of current knowledge regarding transcriptional regulatory networks and post-translational regulation of the transcription factors involved in Cd response. An increasing number of reports indicate that epigenetic regulation and long non-coding and small RNAs are important in Cd-induced transcriptional responses. Several kinases play important roles in Cd signaling that activate transcriptional cascades. We also discuss the perspectives to reduce grain Cd content and improve crop tolerance to Cd stress, which provides a theoretical reference for food safety and the future research of plant varieties with low Cd accumulation.

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.

N6-methyladenosine (m6A) RNA modification as a metabolic switch between plant cell survival and death in leaf senescence

Crop losses caused by climate change and various (a)biotic stressors negatively affect agriculture and crop production. Therefore, it is vital to develop a proper understanding of the complex response(s) to (a)biotic stresses and delineate them for each crop plant as a means to enable translational research. In plants, the improvement of crop quality by m⁶A editing is believed to be a promising strategy. As a reaction to environmental changes, m⁶A modification showed a high degree of sensitivity and complexity. We investigated differences in gene medleys between dark-induced leaf senescence (DILS) and developmental leaf senescence in barley, including inter alia RNA modifications active in DILS. The identified upregulated genes in DILS include RNA methyltransferases of different RNA types, embracing enzymes modifying mRNA, tRNA, and rRNA. We have defined a decisive moment in the DILS model which determines the point of no return, but the mechanism of its control is yet to be uncovered. This indicates the possibility of an unknown additional switch between cell survival and cell death. Discoveries of m⁶A RNA modification changes in certain RNA species in different stages of leaf senescence may uncover the role of such modifications in metabolic reprogramming. Nonetheless, there is no such data about the process of leaf senescence in plants. In this scope, the prospect of finding connections between the process of senescence and m⁶A modification of RNA in plants seems to be compelling.

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.

Development of human retinal organoid models for bisphenol toxicity assessment

Bisphenols, including Bisphenol A (BPA), Tetrabromobisphenol A (TBBPA), and Tetrabromobisphenol S (TBBPS), have been widely applied in the production of polycarbonate plastics and epoxy resins and have been detected in the environment worldwide. The frequent detection of bisphenols in maternal and fetal samples has raised concerns about their toxic effects on human embryonic development, especially on the development of the central nervous system. However, the effect of bisphenols on human retinal development is still unknown. In this study, to evaluate the toxicity of bisphenols on early retinal development, human embryonic stem cells were induced to differentiate into retinal organoids that responded to BPA, TBBPA, and TBBPS, at human exposure relevant concentrations. The global gene expression of retinal organoids was analyzed by RNA sequencing (RNA-seq). A set of retinal development-related biological processes, including neuron differentiation, phototransduction, axon guidance, and retina layer formation, were identified in retinal organoids corresponding to different developmental stages. The RNA-seq data also showed that BPA, TBBPA, and TBBPS influenced retinal development by interfering with the Cytokine-cytokine receptor interaction pathway. HSPA6, HIF1A-AS3, CDC20B, IL19, OAS1, HSPA7, and RN7SK were dysregulated by these chemicals. Additionally, BPA, TBBPA, and TBBPS exhibited different toxic effects on neural retina development, with TBBPA appearing to exert more toxicity than BPA and TBBPS. Furthermore, three bisphenols exhibited different effects at different stages of neural retina development. The sensitivity of retinal development to bisphenols depends on their developmental stage. This study provides new insights into the deep dissection of retinotoxicity after prenatal bisphenol exposure.

Transcriptome-wide N6-methyladenosine (m6A) methylation in soybean under Meloidogyne incognita infection

N6-methyladenosine (m6A) is a reversible epigenetic modification of mRNA and other RNAs that plays a significant role in regulating gene expression and biological processes. However, m6A abundance, dynamics, and transcriptional regulatory mechanisms remain unexplored in the context of soybean resistance to Meloidogyne incognita. In this study, we performed a comparative analysis of transcriptome-wide m6A and metabolome profiles of soybean root tissues with and without M. incognita infection. Global m6A hypermethylation was widely induced in response to M. incognita infection and was enriched around the 3' end of coding sequences and in 3' UTR regions. There were 2069 significantly modified m6A sites, 594 differentially expressed genes, and 103 differentially accumulated metabolites between infected and uninfected roots, including coumestrol, psoralidin, and 2-hydroxyethylphosphonate. Among 101 m6A-modified DEGs, 34 genes were hypomethylated and upregulated, and 39 genes were hypermethylated and downregulated, indicating a highly negative correlation between m6A methylation and gene transcript abundance. A number of these m6A-modified DEGs, including WRKY70, ERF60, POD47 and LRR receptor-like serine/threonine-protein kinases, were involved in plant defense responses. Our study provides new insights into the critical role of m6A modification in early soybean responses to M. incognita.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00077-2.

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.

Profiling of Transcriptome-Wide N6-Methyladenosine (m6A) Modifications and Identifying m6A Associated Regulation in Sperm Tail Formation in Anopheles sinensis

Recent discoveries of reversible N6-methyladenosine (m6A) methylation on messenger RNA (mRNA) and mapping of m6A methylomes in many species have revealed potential regulatory functions of this RNA modification by m6A players—writers, readers, and erasers. Here, we first profile transcriptome-wide m6A in female and male Anopheles sinensis and reveal that m6A is also a highly conserved modification of mRNA in mosquitoes. Distinct from mammals and yeast but similar to Arabidopsis thaliana, m6A in An. sinensis is enriched not only around the stop codon and within 3′-untranslated regions but also around the start codon and 5′-UTR. Gene ontology analysis indicates the unique distribution pattern of m6A in An. sinensis is associated with mosquito sex-specific pathways such as tRNA wobble uridine modification and phospholipid-binding in females, and peptidoglycan catabolic process, exosome and signal recognition particle, endoplasmic reticulum targeting, and RNA helicase activity in males. The positive correlation between m6A deposition and mRNA abundance indicates that m6A can play a role in regulating gene expression in mosquitoes. Furthermore, many spermatogenesis-associated genes, especially those related to mature sperm flagellum formation, are positively modulated by m6A methylation. A transcriptional regulatory network of m6A in An. sinensis is first profiled in the present study, especially in spermatogenesis, which may provide a new clue for the control of this disease-transmitting vector.