Epitranscriptomic Analysis of m6A Methylome After Peripheral Nerve Injury

N6-methyladenosine (m6A) is one of the most plentiful internal RNA modifications, especially in eukaryotic messenger RNA (mRNA), which plays pivotal roles in the regulation of mRNA life cycle and nerve development. However, the mRNA m6A methylation pattern in peripheral nervous injury (PNI) has not been investigated. In this study, sciatic nerve samples were collected from 7 days after sciatic nerve injury (SNI) and control rats. Quantitative real-time PCR demonstrated that m6A-related methyltransferase/demethylase genes were remarkably upregulated in SNI group compared with control group. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) was performed to reveal the m6A methylation landscape. The results showed that 4,014 m6A peaks were significantly altered, including 2,144 upregulated and 1,870 downregulated m6A peaks, which were corresponded to 1,858 genes. Moreover, 919 differentially expressed genes were identified by the conjoint analysis of MeRIP-seq and RNA-seq. GO and KEGG pathway analyses were performed to determine the biological functions and signaling pathways of the m6A-modified genes. Notably, these genes were mainly related to the immune system process, cell activation, and nervous system development in GO analysis. KEGG pathway analysis revealed that these genes were involved in the cell cycle, B cell receptor signaling pathway, axon guidance pathway, and calcium signaling pathway. Furthermore, the m6A methylation and protein expression levels of autophagy-related gene (Atg7) were increased, together with the activation of autophagy. These findings shed some light on the epigenetic regulation of gene expression, which may provide a new opinion to promote functional recovery after PNI.

Introduction

Peripheral nerve repair following injury is a complicated process. The distal stump of injured nerve undergoes Wallerian degeneration due to the loss of nutrition and energy support (Waller, 1851; Stoll et al., 2002). Neuronal apoptosis/necrosis and axon disintegration may induce the dedifferentiation and reprogramming of Schwann cells (Arthur-Farraj et al., 2017; Clements et al., 2017). The injury-induced Schwann cells can proliferate to form Bungner’s bands that establish connections between broken stumps and guide the direction of axon regeneration (Waller, 1851). Schwann cells also inhibit the production of myelin protein and remove myelin debris, which contain axon growth inhibitory proteins (e.g., myelin-associated glycoprotein) (Shen et al., 1998). At the same time, Schwann cells secrete various trophic factors that accelerate axon regeneration coincided with the recruitment of immune cells (neutrophils, macrophages, and T lymphocytes) to the nerve injury sites (Gaudet et al., 2011; Richner et al., 2014). Hematogenous macrophages transport across the high permeability blood–nerve barrier to phagocytose myelin debris and enhance distal nerve regeneration (Chen et al., 2015). Apart from the above, Schwann cells activate autophagy to relieve local pressure and balance the nerve microenvironment during axonal regeneration (Pereira et al., 2012). On the 7th day after peripheral nerve injury (PNI), the number and proliferative activity of activated Schwann cells attain the highest level (Zhou et al., 2017). However, the recovery of motor and sensory functions is limited in some patients with severe PNI, and the molecular mechanisms underlying this phenomenon are still unclear.

The roles of epigenetic modifications in the peripheral nervous system are being actively studied, especially histone acetylation (Puttagunta et al., 2014), and DNA methylation (Loh et al., 2017). It has been proven that these epigenetic changes can promote the transcriptional activation of multiple regeneration-associated genes. However, the m6A modification in axon regeneration remains to be characterized. Similar to other epigenetic modifications, m6A methylation is also dynamic and reversible (Füllgrabe et al., 2014; Zhao et al., 2017; Weng et al., 2018). The formation of m6A is regulated by both methyltransferases (WTAP, METTL3, and METTL14; termed as “writers”) and demethylases (FTO and ALKBH5; termed as “erasers”) (Yeo and O’Rahilly, 2012; Zheng et al., 2013; Yue et al., 2015), and then m6A-recognizing proteins (YTHDF1/2, IGF2BP1, and HNRNPA2B1; termed as “readers”) bind to the specific site to mediate its function (Weng et al., 2018). Some studies have demonstrated that m6A methylation affects almost every stage in the messenger RNA (mRNA) cycle by altering RNA structure, regulating maturation, promoting translation, and accelerating mRNA decay (Zhao et al., 2017).

In this study, we first assessed the expression levels of methyltransferase and demethylase genes in sciatic nerve injury (SNI) and control groups. Then, we compared the m6A-tagged transcription profiles before and after SNI via methylated RNA immunoprecipitation sequencing (MeRIP-seq) and identified differentially methylated peaks. After that, the data were combined with RNA sequencing (RNA-seq) to further analyze the relationship between mRNA expression and m6A methylation. Finally, we found that the mRNA m6A methylation and protein expression levels of Atg7 were increased after SNI. ATG7 acts as the hub of microtubule-associated protein light chain 3 (LC3)-phosphatidylethanolamine and ATG12, thereby facilitating the expansion of autophagosomal membranes (Feng et al., 2015). This provides us with a new theoretical basis for PNI treatment.

Materials and Methods

Animal Model of Sciatic Nerve Injury

Female Wistar rats (100–120 g, 4 weeks old) were randomly divided into control and SNI groups at different time points. The animal model of SNI was constructed as described previously (Zhou et al., 2017). We collected sciatic nerve samples (about 1 cm) at the center of the injury as well as from the control rats, and then frozen at −80°C until further analysis. All animals were maintained under a standard condition (12:12-h light/dark cycle), with unlimited access to food and water. The experimental procedures involving animals were conducted in accordance with the Institutional Animal Care guidelines of the Radiation Medicine Chinese Academy of Medical Sciences and were approved by the Ethics Committee of the Radiation Medicine Chinese Academy of Medical Sciences.

RNA Isolation and qRT-PCR

The mRNA expression levels of m6A writers and erasers were detected by quantitative real-time PCR (qRT-PCR). Total RNA was isolated from the sciatic nerve tissues in control (uninjured) and SNI (7 days after injury) groups using a TRIzol Reagent (15596026, Invitrogen). The purity and quantity of RNA samples were assessed by NanoDrop ND-1000 instrument. cDNA was prepared by reverse transcription using a FastKing RT kit with gDNase (KR116, Tiangen, China). QRT-PCR was performed using the Ultra SYBR Mixture (CW0957, Cwbio, China) according to the manufacturer’s specifications. GAPDH was used as an internal control. Each experiment was repeated at least three times, and all data were presented as mean ± standard deviation (SD). The primer sequences are shown in Table 1.

TABLE 1

Table 1. Primers used for qRT-PCR analysis.

MeRIP-seq and RNA-seq

The sequencing experiments were conducted at Cloud-seq Biotech Inc. (Shanghai, China). First, we collected and quantified total RNA from the sciatic nerve tissues of control (uninjured) and SNI (7 days after injury) groups. Then, Arraystar Seq-StarTMpoly(A) mRNA Isolation Kit (Arraystar, MD, United States) was used to extract poly(A) mRNA from total RNA. The tubes with 1 μg/μL RNA and fragmentation buffer (metal-ion) were incubated at 70°C for 6 min in a preheated thermal cycler block. Subsequently, 30 μL of 3 M sodium acetate (NaOAc), 1 μL of glycoBlue, and 750 μL volumes of 100% ethanol were added, mixed and kept at −80°C overnight. After centrifugation and pellet resuspension, RNA fragmentation was conducted by running 0.5 μg of RNA through 1.5% (wt/v) agarose gel. The fragmented RNA was immunoprecipitated with magnetic beads containing m6A antibody (1 μL m6A antibody in 50 μL 1 × IP buffer) using the GenSeq^TM m6A RNA IP kit (GenSeq Inc., China). The input samples without immunoprecipitation were used as a control group. The following reagents were used for elution and purification in sequence: 1 × IP buffer, LB Buffer, HS Buffer, RLT Buffer, 100% and 75% ethanol.

MeRIP-seq libraries were constructed using the NEBNext^® Ultra II Directional RNA Library Prep Kit (New England Biolabs, Inc., United States). For RNA-seq, total RNA (1 μg) was used to remove the rRNAs via Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, United States). RNA-seq libraries were constructed using the TruSeq Stranded Total RNA Library Prep Kit (Illumina, San Diego, CA, United States). Then, both libraries were controlled and evaluated with the Bio Analyzer 2100 system (Agilent Technologies, United States). The 10 pM libraries were denatured as single-stranded DNA molecules and amplified in situ as clusters. MeRIP-seq and RNA-seq procedures were carried out on an Illumina NovaSeq 6000 instrument with 150-bp paired-end reads. The quality control metrics of the sequencing data are listed in Supplementary Table 1.

Sequencing Data Analysis

For MeRIP-seq data analysis, paired-end reads were conducted via Illumina NovaSeq 6000 sequencer and were quality checked with Q30. Cutadapt software (v1.9.3) was used to remove 3′ adaptor-trimming and low-quality reads (Martin, 2011). Clean reads of all libraries were mapped to the reference genome (Rat5) via Hisat2 software (v2.0.4) (Kim et al., 2015). Peaks with m6A modification were detected by MACS software (Zhang et al., 2008). Differentially methylated sites were identified by diffReps and visualized through IGV software (Shen et al., 2013). The peaks recognized by both software and overlapped with two or more exons were filtered out and annotated by homemade scripts. MEME software was used to perform motif-based sequence analysis (Bailey, 2011). GO and KEGG pathway analyses were conducted to identify the biological functions and signaling pathways of differentially m6A-modified genes.

For mRNA sequencing data analysis, the paired-end reads from Illumina NovaSeq 6000 sequencer were quality checked with Q30, and the low-quality reads were removed by Cutadapt software (v1.9.3) (Martin, 2011). Then, the high-quality clean reads from MeRIP-seq were aligned to the reference genome (Rat5) via Hisat2. Cuffdiff software (part of Cufflinks) was used to calculate the gene level FPKM (fragments per kilobase of exon per million fragments mapped reads) as the expression profiles of mRNA. Differentially expressed genes were identified based on the fold-change and p-value calculated from FPKM values (Trapnell et al., 2010). GO and KEGG pathway analyses were performed to assess the biological functions and signaling pathways of differentially expressed genes.

Gene-Specific m6A qPCR

Gene-specific m6A qPCR was carried out at Cloud-seq Biotech Inc. First, we obtained and quantified the total RNA from both control (uninjured) and SNI (7 days after injury) groups. Then, metal-ion was used to fragment the extracted RNA to approximately 100 nt in length and reserved one-tenth of fragmented RNA as INPUT control. Next, the fragmented RNA were subjected to immunoprecipitation with m6A antibody-coated magnetic beads. Methylated RNA bound to the magnetic beads was eluted and purified. Finally, the m6A modification of Atg7 gene was verified by qPCR analysis. The primer sequences are listed in Table 1.

Western Blot Analysis

The sciatic nerve samples obtained from control (uninjured) and SNI (2, 5, 7, and 10 days after injury) groups were ground using an automatic grinder, and then lysed with RIPA lysis buffer. BCA protein assay kit (P0012, Beyotime) was used to measure the concentration of total protein. An equal amount of protein sample was separated by 10% SDS-PAGE, and then transferred onto polyvinylidene difluoride membranes. After blocking with 5% non-fat milk for 1 h, the membranes were incubated with primary antibody overnight at 4°C. Subsequently, the membranes were incubated with the corresponding secondary antibody for 1 h. Enhanced chemiluminescent substrate (Bio-Rad) was used to visualize the protein bands. The antibodies used for Western blotting were as follows: ATG7 (1:1000, 8558S; Cell Signaling Technology), LC3 (1:1,000, L7543; Sigma-Aldrich), and GAPDH (1:2,000, ab181602; Abcam).

Transmission Electron Microscopy

The fresh samples of control and SNI (2 days after injury) groups were fixed with 2.5% glutaraldehyde overnight at 4°C and then with OsO4 at 4°C for 1 h. Different concentrations of ethanol solution (50, 70, 80, 90, 95, and 100%) were used for gradient dehydration. The samples were permeated with acetone and Epon embedding agent, sectioned using an ultramicrotome, and stained with lead citrate solution and uranyl acetate saturated in 50% ethanol. Autophagosomes were observed under a JEM-1010 electron microscope (Jeol, Tokyo, Japan) at × 3,000 magnification.

Statistical Analysis

All statistical analyses were performed using the GraphPad Prism software. The results were presented as mean ± SD. ImageJ software was used to calculate the gray values of Western blot images. Significant differences between two groups were assessed by Student’s t test, ANOVA, and post hoc analysis. Fisher’s exact test was used for all bioinformatic analyses.

Results

Ubiquitous Expression of m6A Writers and Erasers in Sciatic Nerve After Injury

To investigate the effects of epigenetic modulators on rat sciatic nerve (7 days after SNI), we determined the expression profiles of m6A writers (or known as methyltransferases; METTL3, METTL14, and WTAP) and erasers (or known as emethylases; FTO and ALKBH5). As shown in Figure 1, the mRNA expression levels of these methyltransferases and demethylases were remarkably higher in SNI group than in control group. Among them, the expression levels of FTO (26.1) and METTL14 (30.1) were more significantly upregulated compared with other key genes. The ubiquitous expression of these genes indicates that m6A plays a prominent role during SNI.

FIGURE 1

Figure 1. The expression levels of writers and erasers. Quantitative RT-PCR showing the upregulated levels of METTL3, METTL14, WTAP, FTO, and ALKBH5. Among them, the expression levels of METTL14 and FTO were the most significant. GAPDH was used as an internal control. SNI group (n = 3); control group (n = 3). “**” and “***” indicate p < 0.01 and p < 0.001, respectively, via Student’s t test.

Transcriptome Profile of m6A Methylation in SNI Samples

To analyze the level and distribution of m6A methylation in SNI and control groups, MeRIP-seq was performed to assess the epitranscriptome m6A-seq data (GEO accession number: GSE171866). In total, SNI samples had 2,144 and 1,870 significantly upregulated and downregulated m6A peaks, which corresponded to the transcripts of 1,083 and 775 genes, respectively (fold-change ≥ 2.0 and p < 0.00001; Figure 2A). Among them, the top 20 differentially methylated m6A peaks are shown in Table 2, and the remaining m6A peaks can be found in Supplementary Table 2.

FIGURE 2

Figure 2. The degree and pattern of m6A modification peaks in SNI and control groups. (A) Volcano plots indicating the distribution of differential m6A peaks after SNI. Fold-change ≥ 2.0, p < 0.00001, n = 3, via Fisher’s exact test. (B) Priority region and average distribution of m6A peaks. (C,D) The distribution of m6A peaks in SNI and control groups, respectively. (E) The top motifs enriched across differential m6A peaks. (F) Data visualization of the m6A-modified Gtse1 mRNA.

TABLE 2

Table 2. The top 20 altered methylated m6A mRNAs.

We further identified the preferential locations of m6A peaks in mRNA coding regions. The results demonstrated that the peaks were markedly correlated with the coding sequence (CDS), and there were significantly enriched peaks at the beginning of 3′ untranslated region (3′UTRs) and near the stop codon (Figures 2B–D). Motif enrichment analysis indicated that GGACU was the most enriched m6A sites of were in both SNI and control groups (Figure 2E), which provides further proof for the conservativeness of m6A methylation. Gtse1, a significantly hypermethylated peak on the CDS, was also detected (Figure 2F). These results strongly suggest that m6A methylation is dynamic and reversible.

Biological Functions and Signaling Pathways Associated With Differentially Methylated mRNA

To clarify the biological functions and signaling pathways of m6A modification during SNI, both GO enrichment and KEGG pathway analyses were carried out. In GO analysis, the biological functions of differentially methylated genes can be classified into three parts: biological processes (BP), cellular component (CC), and molecular function (MF). The top 10 most remarkably enriched BPs, CCs, and MFs of the genes with up/downregulated m6A peaks are listed in Figures 3A,B. Through the KEGG pathway analysis, the crucial pathways of differentially m6A-modified genes were identified. The results showed that the hypermethylated genes were involved in the cell cycle, B cell receptor signaling pathway, and Toll-like receptor (TLR) signaling pathway (Figure 3C). The hypomethylated genes were involved in the calcium signaling pathway, Rap1 signaling pathway, and axon guidance (Figure 3D). This indicates m6A methylation is a potential target to improve sciatic nerve repair after injury.

FIGURE 3

Figure 3. GO and KEGG pathway analyses of the differentially methylated genes. Significantly enriched GO terms of the hypermethylated (A) and hypomethylated genes (B). Significantly enriched KEGG pathways of the hypermethylated (C) and hypomethylated (D) genes.

Conjoint Analysis of MeRIP-seq and RNA-seq Data

To explore the association between gene expression and m6A methylation, RNA-seq data were used to determine the extent of altered mRNA expression (GEO accession number: GSE171866). The differentially expressed transcripts consisted of 2,055 upregulated transcripts and 785 downregulated transcripts (fold-change ≥ 2.0 and p < 0.05; Figure 4A). The top 20 differentially expressed genes are listed in Table 3, while the top 10 GO terms and KEGG pathways are presented in Supplementary Figure 2. Cluster analysis was performed to assess the expression patterns of differential methylated genes (Figure 4B). In the light of the conjoint analysis of MeRIP-seq and RNA-seq data, all differentially m6A peaks with differential mRNA levels were divided into four groups (Figure 4C). Notably, all 635 hypermethylated m6A peaks in mRNA transcripts were upregulated, while 283 of the 284 hypomethylated m6A peaks were downregulated in mRNA transcripts. These results suggest that there is a strong correlation between m6A methylation and gene expression.

FIGURE 4

Figure 4. Conjoint analysis of m6A methylation and mRNA expression after SNI. (A) Volcano plots and (B) Heatmap plots displaying the differentially expressed genes between SNI and control groups. (C) Four-quadrant graph showing the association between mRNA m6A methylation and transcriptome level.

TABLE 3

Table 3. The top 20 differently expressed mRNAs in SNI based on log2 (fold change).

The m6A Methylation and Protein Expression Level of Atg7 Are Increased in SNI Group

Autophagy plays a vital role in the repair of SNI. In this study, we found that the m6A abundance of Atg7 mRNA was obviously increased after SNI (Figure 5A). To further verify the expression of Atg7 mRNA methylation, gene-specific m6A qPCR was performed. Interestingly, the result of gene-specific m6A qPCR was consistent with that of m6A-seq (Figure 5B). The protein expression of Atg7 was also upregulated, indicating m6A methylation may serve as a potential mechanism to regulate autophagy during SNI (Figures 5C,E). In addition, the expression of the autophagy-related marker LC3 was increased (Figures 5D,E), and the autophagosomes were abundantly found in sciatic nerve after injury (Figure 5F). Altogether, our findings reveal a prominent association between the protein expression of Atg7 and its mRNA methylation.

FIGURE 5

Figure 5. The m6A methylation and protein expression of Atg7 gene are altered, and autophagy is activated following PNI. (A,B) Integrative genomics viewer and gene-specific m6A qPCR demonstrating the increased methylation of Atg7 gene after SNI. n = 3, via Student’s t test. (C) Upregulation of ATG7 protein level. n = 3, via Western blot analysis. (D) Expression level of the autophagy-related marker LC3 is increased. n = 3. (E) Quantification of Western blots for ATG7 and LC3II/LC3I. n = 3, via ANOVA and post hoc analysis. (F) Transmission electron microscopy showing autophagosome (yellow arrow) and lysosome (blue arrow) in the sciatic nerve tissues. “*,” “**,” and “***” indicate p < 0.05, p < 0.01, and p < 0.001, respectively.

Discussion

The roles of m6A modification in neurodevelopment have been verified in these aspects, including learning and memory (Li et al., 2019). Despite these efforts, the epigenetic changes in rat sciatic nerve remain poorly understood. In this study, we obtained a landscape of m6A methylation after SNI via MeRIP-Seq. There were 4,014 differentially m6A peaks, corresponding to 2,144 genes, and each gene was associated with at least one or more m6A peaks. The distribution of m6A peaks appeared in the CDS, which was also consistent with other diseases such as cancer (Zhang et al., 2020b), traumatic brain injury (Yu et al., 2020), and preeclampsia (Wang et al., 2020a). After PNI, these m6A peaks preferentially located at the start codons (19.17 vs. 17.49%), and there was a decreasing trend in the CDS (66.75 vs. 68.49%). The motif sequence GGACU displayed remarkably higher proportions than other motifs in both control and SNI groups, as well as in other species (e.g., human and mice). Besides, it is noteworthy that RRACH (where R denotes purine, A stands for m6A, and H represents a non-guanine base) is their consensus sequence (Dominissini et al., 2012; Zhang et al., 2020a).

In addition, we found that the mRNA expression levels of methylase and demethylase genes were upregulated in SNI group compared with control group. GAPDH is generally used as a housekeeping gene for the normalization of gene expression data. However, in our results, the Ct values of GAPDH in SNI and control groups were 16.9 and 19.3, respectively. This might be due to the fact that our samples were derived from rat sciatic nerve tissue. Compared with cells, tissue samples are more difficult to obtain the same Ct values. Some studies have pointed out that the mRNA expression of β-actin is increased after axon injury and may play a certain role in the acceleration of axonal outgrowth (Sántha et al., 2012; Donnelly et al., 2013).

The upregulated expression levels of FTO and METTL14 were more significant compared with other key genes. FTO and METTL14 act as eraser and writer, respectively, which regulate the dynamic levels of m6A RNA methylation. Some studies have suggested that FTO promotes the initiation of autophagy by reducing the mRNA methylation levels of Atg5, Atg7, and Ulk1 (Jin et al., 2018; Wang et al., 2020b). Autophagy is activated after PNI to promote nerve regeneration and motor recovery (Huang et al., 2016). However, the methylation level of Atg7 mRNA was increased in our study, and there were no significant differences in the mRNA methylation patterns of Atg5 and Ulk1 between SNI and control groups. In the next step, we will determine which methyltransferase regulates Atg7 mRNA and construct a conditional gene knockout model to verify its role in PNI. Perhaps that it is a new way to clarify the underlying mechanism of autophagy after PNI.

Fto depletion inhibits the length of dendrites and the number of intersections in primary neurons and dynamically regulates the proliferation and differentiation of adult neural stem cells (Cao et al., 2020). Another study reported that a conditional knockout of Fto in brain lipid inhibits adult neurogenesis by inducing cell apoptosis (Gao et al., 2020). FTO promotes neuropathic pain by elevating G9a expression in primary sensory neurons (Li et al., 2020). G9a is also involved in the regulation of autophagy via linking to the promoter of Atg gene (Talebian et al., 2020). On the other hand, FTO can be bound with a variety of reaction substrates (e.g., m6Am and m1A), suggesting the complexity and functionality of this gene (Wei et al., 2018).

METTL14 plays an important role in the formation of multicomponent methyltransferase complex as a reaction catalyst to activate and promote METTL3 activity (Wang et al., 2017). In the central nervous system, depletion of Mettl14 prevents the differentiation and maturation of oligodendrocytes, thus resulting in a decreased number of mature cells and hypomyelination (Xu et al., 2020). After PNI, the premyelinating phenotypes of Schwann cells are characterized by neural cell adhesion molecule, glial fibrillary acidic protein, and p75 neurotrophin receptor (Chen et al., 2007). Our results also showed that METTLE4 was upregulated, suggesting that it may have the potential to promote the development, maturation, and regeneration of Schwann cells. Besides, Mettl14 knockout in cortical neural progenitor cells can slow down the cell-cycle progression. Thus, it is speculated that METTL14 may also affect the development and regeneration of neural progenitor cells after PNI (Yoon et al., 2017).

KEGG pathway analysis revealed that the TLR signaling pathway was regulated by the genes with upregulated m6A modification sites during SNI. TLRs are dependent on pathogen-related molecular patterns to induce activities that stimulate immunity and host defense (Fitzgerald and Kagan, 2020). Recent findings have demonstrated that the expression of TLRs is upregulated in activated Schwann cells after nerve injury, which are involved in Wallerian degeneration and proinflammatory mediator generation, thereby recruiting more immune cells into the nerve injury sites (Thakur et al., 2017). TLR7 and TLR8 participate in neuropathic pain by activating NF-κB and Erk signaling pathways, respectively (Zhang et al., 2018; He et al., 2020). In Tlr2 knockout mice, thermal hyperalgesia and mechanical allodynia were attenuated compared with wild-type mice (Kim et al., 2007). According to our results, we believe that modulating m6A modifications of the TLR signaling pathway may be a potential therapeutic target for SNI in the near future.

Axon guidance pathway was regulated by the genes with downregulated m6A-modified sites during SNI. These genes included Ablim1, Cxcl12, Dpysl5, Ephb6, L1cam, Lrrc4, Ntn1, Plxnb1, Robo3, Sema4g, Sema5a, Sema6d, Slit2, and Slit3, all of which could be used as potential targets for regulating axon growth. CXCL12 is primarily responsible for Schwann cell migration and autophagy by regulating PI3K-AKT-mTOR signaling pathway (Gao et al., 2019); L1CAM promotes Schwann cell-neurite interactions and enhance neurite outgrowth (Romano et al., 2015); NTN1 could induce Schwann cell proliferation and migration to form the Büngner bands and extend toward the distal nerve stump (Dun and Parkinson, 2017).

Based on the conjoint analysis of MeRIP-seq and RNA-seq data, 635 hypermethylated m6A genes were upregulated, and 283 of 284 hypomethylated m6A genes were downregulated. There was a significant correlation between m6A methylation and mRNA expression level. Indeed, m6A not only affects mRNA folding and maturation to induce protein expression, but also regulates mRNA decay to curtail its lifetime (Zhao et al., 2017). Both methyltransferase and demethylase genes can determine the m6A status of a random mRNA; and in turn, the m6A-binding proteins can target and categorize the methylated mRNA into different functional groups. For example, YTHDF1 interacts with eIF3 to increase translation efficiency (Wang et al., 2015); while YTHDF2 binds to mRNA decay sites, thereby resulting in mRNA instability and degradation (Wang et al., 2014).

Among the 635 hypermethylated m6A genes, three genes (Atf3, Creb1, and Jun) were associated with functional recovery after PNI. In Atf3 mutant mice, the axon regeneration ability of DRG neurons is weakened (Gey et al., 2016). Overexpression of Creb enhances the regeneration of lesioned dorsal column axons (Gao et al., 2004). The loss of neuron Jun decreases the rate of axon regeneration after crushing, and blocks most severe axons from reconnecting to their targets (Ruff et al., 2012). Therefore, it can be inferred that m6A-modified genes play critical roles in the development of peripheral nervous system. However, detailed molecular mechanisms are still unknown, and this aspect deserves careful consideration in the near future.

Conclusion

This study reveals the m6A RNA methylation landscape via high-throughput sequencing, which sheds new light on the pivotal roles of epigenetic changes in SNI. Conjoint analysis of m6A-RIP-seq and RNA-seq indicated the effect of m6A methylation on regulating gene expression. The m6A methylation and protein expression level of Atg7 were increased, laying the foundation for understanding the epitranscriptomic mechanisms of autophagy. Further experiments should be carried out to verify the differential methylated genes in SNI patients.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI Gene Expression Omnibus, accession no: GSE171866.

Ethics Statement

The animal study was reviewed and approved by Ethics Committee of Radiation Medicine Chinese Academy of Medical Sciences.

Author Contributions

LZ and DH participated in the design of the research and the writing of the manuscript. PM, BM, and JQ conducted the experiments. GT and ZL analyzed the data. XZ revised the manuscript and were responsible for supervision of study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (81930070 and 81371957) and Tianjin key Research and Development Plan, Key Projects for Science and Technology Support(19YFZCSY00660).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.686000/full#supplementary-material

References

Arthur-Farraj, P. J., Morgan, C. C., Adamowicz, M., Gomez-Sanchez, J. A., Fazal, S. V., Beucher, A., et al. (2017). Changes in the coding and non-coding transcriptome and DNA methylome that define the schwann cell repair phenotype after nerve injury. Cell Rep. 20, 2719–2734. doi: 10.1016/j.celrep.2017.08.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Bailey, T. L. (2011). DREME: motif discovery in transcription factor ChIP-seq data. Bioinformatics 27, 1653–1659. doi: 10.1093/bioinformatics/btr261

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, Y., Zhuang, Y., Chen, J., Xu, W., Shou, Y., Huang, X., et al. (2020). Dynamic effects of Fto in regulating the proliferation and differentiation of adult neural stem cells of mice. Hum. Mol. Genet. 29, 727–735. doi: 10.1093/hmg/ddz274

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, P., Piao, X., and Bonaldo, P. (2015). Role of macrophages in Wallerian degeneration and axonal regeneration after peripheral nerve injury. Acta Neuropathol. 130, 605–618. doi: 10.1007/s00401-015-1482-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z. L., Yu, W. M., and Strickland, S. (2007). Peripheral regeneration. Annu. Rev. Neurosci. 30, 209–233. doi: 10.1146/annurev.neuro.30.051606.094337

PubMed Abstract | CrossRef Full Text | Google Scholar

Clements, M. P., Byrne, E., Camarillo Guerrero, L. F., Cattin, A. L., Zakka, L., Ashraf, A., et al. (2017). The wound microenvironment reprograms schwann cells to invasive mesenchymal-like cells to drive peripheral nerve regeneration. Neuron 96, 98–114.e117. doi: 10.1016/j.neuron.2017.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dominissini, D., Moshitch-Moshkovitz, S., Schwartz, S., Salmon-Divon, M., Ungar, L., Osenberg, S., et al. (2012). Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature 485, 201–206. doi: 10.1038/nature11112

PubMed Abstract | CrossRef Full Text | Google Scholar

Donnelly, C. J., Park, M., Spillane, M., Yoo, S., Pacheco, A., Gomes, C., et al. (2013). Axonally synthesized β-actin and GAP-43 proteins support distinct modes of axonal growth. J. Neurosci. 33, 3311–3322. doi: 10.1523/jneurosci.1722-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Dun, X. P., and Parkinson, D. B. (2017). Role of netrin-1 signaling in nerve regeneration. Int. J. Mol. Sci. 18:491. doi: 10.3390/ijms18030491

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, Y., Yao, Z., and Klionsky, D. J. (2015). How to control self-digestion: transcriptional, post-transcriptional, and post-translational regulation of autophagy. Trends Cell Biol. 25, 354–363. doi: 10.1016/j.tcb.2015.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzgerald, K. A., and Kagan, J. C. (2020). Toll-like receptors and the control of immunity. Cell 180, 1044–1066. doi: 10.1016/j.cell.2020.02.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Füllgrabe, J., Klionsky, D. J., and Joseph, B. (2014). The return of the nucleus: transcriptional and epigenetic control of autophagy. Nat. Rev. Mol. Cell. Biol. 15, 65–74. doi: 10.1038/nrm3716

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, D., Tang, T., Zhu, J., Tang, Y., Sun, H., and Li, S. (2019). CXCL12 has therapeutic value in facial nerve injury and promotes Schwann cells autophagy and migration via PI3K-AKT-mTOR signal pathway. Int. J. Biol. Macromol. 124, 460–468. doi: 10.1016/j.ijbiomac.2018.10.212

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, H., Cheng, X., Chen, J., Ji, C., Guo, H., Qu, W., et al. (2020). Fto-modulated lipid niche regulates adult neurogenesis through modulating adenosine metabolism. Hum. Mol. Genet. 29, 2775–2787. doi: 10.1093/hmg/ddaa171

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Y., Deng, K., Hou, J., Bryson, J. B., Barco, A., Nikulina, E., et al. (2004). Activated CREB is sufficient to overcome inhibitors in myelin and promote spinal axon regeneration in vivo. Neuron 44, 609–621. doi: 10.1016/j.neuron.2004.10.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaudet, A. D., Popovich, P. G., and Ramer, M. S. (2011). Wallerian degeneration: gaining perspective on inflammatory events after peripheral nerve injury. J. Neuroinflammation 8:110. doi: 10.1186/1742-2094-8-110

PubMed Abstract | CrossRef Full Text | Google Scholar

Gey, M., Wanner, R., Schilling, C., Pedro, M. T., Sinske, D., and Knöll, B. (2016). Atf3 mutant mice show reduced axon regeneration and impaired regeneration-associated gene induction after peripheral nerve injury. Open Biol. 6:160091. doi: 10.1098/rsob.160091

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., Han, G., Wu, S., Du, S., Zhang, Y., Liu, W., et al. (2020). Toll-like receptor 7 contributes to neuropathic pain by activating NF-κB in primary sensory neurons. Brain Behav. Immun. 87, 840–851. doi: 10.1016/j.bbi.2020.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, H. C., Chen, L., Zhang, H. X., Li, S. F., Liu, P., Zhao, T. Y., et al. (2016). Autophagy promotes peripheral nerve regeneration and motor recovery following sciatic nerve crush injury in rats. J. Mol. Neurosci. 58, 416–423. doi: 10.1007/s12031-015-0672-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, S., Zhang, X., Miao, Y., Liang, P., Zhu, K., She, Y., et al. (2018). m(6)a RNA modification controls autophagy through upregulating ULK1 protein abundance. Cell Res. 28, 955–957. doi: 10.1038/s41422-018-0069-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Kim, M. A., Cho, I. H., Kim, M. S., Lee, S., Jo, E. K., et al. (2007). A critical role of toll-like receptor 2 in nerve injury-induced spinal cord glial cell activation and pain hypersensitivity. J. Biol. Chem. 282, 14975–14983. doi: 10.1074/jbc.M607277200

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Yang, X., Qi, Z., Sang, Y., Liu, Y., Xu, B., et al. (2019). The role of mRNA m(6)A methylation in the nervous system. Cell Biosci. 9:66. doi: 10.1186/s13578-019-0330-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Guo, X., Sun, L., Xiao, J., Su, S., Du, S., et al. (2020). N(6)-methyladenosine demethylase FTO contributes to neuropathic pain by stabilizing G9a expression in primary sensory neurons. Adv. Sci. (Weinh) 7:1902402. doi: 10.1002/advs.201902402

PubMed Abstract | CrossRef Full Text | Google Scholar

Loh, Y., Koemeter-Cox, A., Finelli, M., Shen, L., Friedel, R., and Zou, H. (2017). Comprehensive mapping of 5-hydroxymethylcytosine epigenetic dynamics in axon regeneration. Epigenetics 12, 77–92. doi: 10.1080/15592294.2016.1264560

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. Martin 17:3. doi: 10.14806/ej.17.1.200

CrossRef Full Text | Google Scholar

Pereira, J. A., Lebrun-Julien, F., and Suter, U. (2012). Molecular mechanisms regulating myelination in the peripheral nervous system. Trends Neurosci. 35, 123–134. doi: 10.1016/j.tins.2011.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Puttagunta, R., Tedeschi, A., Sória, M., Hervera, A., Lindner, R., Rathore, K., et al. (2014). PCAF-dependent epigenetic changes promote axonal regeneration in the central nervous system. Nat. Commun. 5:3527. doi: 10.1038/ncomms4527

PubMed Abstract | CrossRef Full Text | Google Scholar

Richner, M., Ulrichsen, M., Elmegaard, S. L., Dieu, R., Pallesen, L. T., and Vaegter, C. B. (2014). Peripheral nerve injury modulates neurotrophin signaling in the peripheral and central nervous system. Mol. Neurobiol. 50, 945–970. doi: 10.1007/s12035-014-8706-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Romano, N. H., Madl, C. M., and Heilshorn, S. C. (2015). Matrix RGD ligand density and L1CAM-mediated schwann cell interactions synergistically enhance neurite outgrowth. Acta Biomater. 11, 48–57. doi: 10.1016/j.actbio.2014.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruff, C. A., Staak, N., Patodia, S., Kaswich, M., Rocha-Ferreira, E., Da Costa, C., et al. (2012). Neuronal c-Jun is required for successful axonal regeneration, but the effects of phosphorylation of its N-terminus are moderate. J. Neurochem. 121, 607–618. doi: 10.1111/j.1471-4159.2012.07706.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sántha, P., Pákáski, M., Fazekas, O., Szucs, S., Fodor, E. K., and Kálmán, J. Jr., et al. (2012). [Acute and chronic stress induced changes in gene transcriptions related to Alzheimer’s disease]. Ideggyogy Sz 65, 195–200.

Google Scholar

Shen, L., Shao, N. Y., Liu, X., Maze, I., Feng, J., and Nestler, E. J. (2013). diffReps: detecting differential chromatin modification sites from ChIP-seq data with biological replicates. PLoS One 8:e65598. doi: 10.1371/journal.pone.0065598

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y. J., DeBellard, M. E., Salzer, J. L., Roder, J., and Filbin, M. T. (1998). Myelin-associated glycoprotein in myelin and expressed by Schwann cells inhibits axonal regeneration and branching. Mol. Cell Neurosci. 12, 79–91. doi: 10.1006/mcne.1998.0700

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoll, G., Jander, S., and Myers, R. R. (2002). Degeneration and regeneration of the peripheral nervous system: from Augustus Waller’s observations to neuroinflammation. J. Peripher. Nerv. Syst. 7, 13–27. doi: 10.1046/j.1529-8027.2002.02002.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Talebian, S., Daghagh, H., Yousefi, B., Őzkul, Y., Ilkhani, K., Seif, F., et al. (2020). The role of epigenetics and non-coding RNAs in autophagy: a new perspective for thorough understanding. Mech. Ageing Dev. 190:111309. doi: 10.1016/j.mad.2020.111309

PubMed Abstract | CrossRef Full Text | Google Scholar

Thakur, K. K., Saini, J., Mahajan, K., Singh, D., Jayswal, D. P., Mishra, S., et al. (2017). Therapeutic implications of toll-like receptors in peripheral neuropathic pain. Pharmacol. Res. 115, 224–232. doi: 10.1016/j.phrs.2016.11.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

Waller, A. (1851). Experiments on the section of the glosso-pharyngeal and hypoglossal nerves of the frog, and observations of the alterations produced thereby in the structure of their primitive fibres. Edinb. Med. Surg. J. 76, 369–376.

Google Scholar

Wang, J., Gao, F., Zhao, X., Cai, Y., and Jin, H. (2020a). Integrated analysis of the transcriptome-wide m6A methylome in preeclampsia and healthy control placentas. PeerJ 8:e9880. doi: 10.7717/peerj.9880

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Feng, J., Xue, Y., Guan, Z., Zhang, D., Liu, Z., et al. (2017). Corrigendum: structural basis of N-adenosine methylation by the METTL3-METTL14 complex. Nature 542:260. doi: 10.1038/nature21073

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Lu, Z., Gomez, A., Hon, G. C., Yue, Y., Han, D., et al. (2014). N6-methyladenosine-dependent regulation of messenger RNA stability. Nature 505, 117–120. doi: 10.1038/nature12730

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Wu, R., Liu, Y., Zhao, Y., Bi, Z., Yao, Y., et al. (2020b). m(6)A mRNA methylation controls autophagy and adipogenesis by targeting Atg5 and Atg7. Autophagy 16, 1221–1235. doi: 10.1080/15548627.2019.1659617

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Zhao, B. S., Roundtree, I. A., Lu, Z., Han, D., Ma, H., et al. (2015). N(6)-methyladenosine modulates messenger RNA translation efficiency. Cell 161, 1388–1399. doi: 10.1016/j.cell.2015.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, J., Liu, F., Lu, Z., Fei, Q., Ai, Y., He, P. C., et al. (2018). Differential m(6)A, m(6)A(m), and m(1)a demethylation mediated by FTO in the cell nucleus and cytoplasm. Mol Cell 71, 973–985.e975. doi: 10.1016/j.molcel.2018.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, H., Huang, H., Wu, H., Qin, X., Zhao, B. S., Dong, L., et al. (2018). METTL14 inhibits hematopoietic stem/progenitor differentiation and promotes leukemogenesis via mRNA m(6)a modification. Cell Stem. Cell 22, 191–205.e199. doi: 10.1016/j.stem.2017.11.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H., Dzhashiashvili, Y., Shah, A., Kunjamma, R. B., Weng, Y. L., Elbaz, B., et al. (2020). m(6)a mRNA methylation is essential for oligodendrocyte maturation and CNS myelination. Neuron 105, 293–309.e295. doi: 10.1016/j.neuron.2019.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeo, G. S., and O’Rahilly, S. (2012). Uncovering the biology of FTO. Mol. Metab 1, 32–36. doi: 10.1016/j.molmet.2012.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, K. J., Ringeling, F. R., Vissers, C., Jacob, F., Pokrass, M., Jimenez-Cyrus, D., et al. (2017). Temporal control of mammalian cortical neurogenesis by m(6)a methylation. Cell 171, 877–889.e817. doi: 10.1016/j.cell.2017.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J., Zhang, Y., Ma, H., Zeng, R., Liu, R., Wang, P., et al. (2020). Epitranscriptomic profiling of N6-methyladenosine-related RNA methylation in rat cerebral cortex following traumatic brain injury. Mol. Brain 13:11. doi: 10.1186/s13041-020-0554-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, Y., Liu, J., and He, C. (2015). RNA N6-methyladenosine methylation in post-transcriptional gene expression regulation. Genes Dev. 29, 1343–1355. doi: 10.1101/gad.262766.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Shi, X., Huang, T., Zhao, X., Chen, W., Gu, N., et al. (2020a). Dynamic landscape and evolution of m6A methylation in human. Nucleic Acids Res. 48, 6251–6264. doi: 10.1093/nar/gkaa347

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Liu, T., Meyer, C. A., Eeckhoute, J., Johnson, D. S., Bernstein, B. E., et al. (2008). Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9:R137. doi: 10.1186/gb-2008-9-9-r137

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Wang, Q., Zhang, M., Zhang, W., Zhao, L., Yang, C., et al. (2020b). Comprehensive analysis of the transcriptome-wide m6A methylome in colorectal cancer by MeRIP sequencing. Epigenetics 16, 425–435. doi: 10.1080/15592294.2020.1805684

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z. J., Guo, J. S., Li, S. S., Wu, X. B., Cao, D. L., Jiang, B. C., et al. (2018). TLR8 and its endogenous ligand miR-21 contribute to neuropathic pain in murine DRG. J. Exp. Med. 215, 3019–3037. doi: 10.1084/jem.20180800

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, B. S., Roundtree, I. A., and He, C. (2017). Post-transcriptional gene regulation by mRNA modifications. Nat. Rev. Mol. Cell. Biol. 18, 31–42. doi: 10.1038/nrm.2016.132

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, G., Dahl, J. A., Niu, Y., Fedorcsak, P., Huang, C. M., Li, C. J., et al. (2013). ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol. Cell 49, 18–29. doi: 10.1016/j.molcel.2012.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X. H., Lin, W., Ren, Y. M., Liu, S., Fan, B. Y., Wei, Z. J., et al. (2017). Comparison of DNA methylation in schwann cells before and after peripheral nerve injury in rats. Biomed Res. Int. 2017:5393268. doi: 10.1155/2017/5393268

PubMed Abstract | CrossRef Full Text | Google Scholar