METTL3-mediated m6A modification regulates cell cycle progression of dental pulp stem cells
Stem Cell Research & Therapy volume 12, Article number: 159 (2021)
Dental pulp stem cells (DPSCs) are a promising cell source in endodontic regeneration and tissue engineering with limited self-renewal and pluripotency capacity. N6-methyladenosine (m6A) is the most prevalent, reversible internal modification in RNAs associated with stem cell fate determination. In this study, we aim to explore the biological effect of m6A methylation in DPSCs.
m6A immunoprecipitation with deep sequencing (m6A RIP-seq) demonstrated the features of m6A modifications in DPSC transcriptome. Lentiviral vectors were constructed to knockdown or overexpress methyltransferase like 3 (METTL3). Cell morphology, viability, senescence, and apoptosis were analyzed by β-galactosidase, TUNEL staining, and flow cytometry. Bioinformatic analysis combing m6A RIP and shMETTL3 RNA-seq functionally enriched overlapped genes and screened target of METTL3. Cell cycle distributions were assayed by flow cytometry, and m6A RIP-qPCR was used to confirm METTL3-mediated m6A methylation.
Here, m6A peak distribution, binding area, and motif in DPSCs were first revealed by m6A RIP-seq. We also found a relatively high expression level of METTL3 in immature DPSCs with superior regenerative potential and METTL3 knockdown induced cell apoptosis and senescence. A conjoint analysis of m6A RIP and RNA sequencing showed METTL3 depletion associated with cell cycle, mitosis, and alteration of METTL3 resulted in cell cycle arrest. Furthermore, the protein interaction network of differentially expressed genes identified Polo-like kinase 1 (PLK1), a critical cycle modulator, as the target of METTL3-mediated m6A methylation in DPSCs.
These results revealed m6A methylated hallmarks in DPSCs and a regulatory role of METTL3 in cell cycle control. Our study shed light on therapeutic approaches in vital pulp therapy and served new insight into stem cell-based tissue engineering.
Adult stem cells residing in various tissues play an essential role in maintaining tissue homeostasis and ensuring regenerative needs. These cells may remain quiescent in a specialized microenvironment for a quite long time and can still differentiate into specific cell types once there are physical damages or activating signals . The regenerative and self-renewal capacity of stem cells highlight the therapeutic approaches not only as regenerative medicine in local injury repair but also in tissue engineering. Adult stem cells with diverse biological properties had been reported in lots of tissues including the bone marrow, blood, muscle, skin, liver, and also teeth. Among these, the dental pulp is rich in mesenchymal stem cells which can be easily obtained from permanent and primary teeth as cell resources. Dental pulp stem cells (DPSCs) derived from cranial neural crest cells are responsible for dentine-pulp complex regeneration and possess superior self-renewal and pluripotency capacity than other adult stem cells [2,3,4]. DPSCs can continuously grow up to thirty passages and still maintain pluripotent to differentiate into odontogenic, chondrogenic, adipogenic, and neurogenic cells .
As other adult stem cells, DPSCs are capable to repair mild injury while fail to cope with irreversible or severe damage due to its limited regenerative potential. The long lifespan and persistence of adult stem cells contributed to the accumulation of multiple damages which resulted in regenerative capacity reduction and insufficient injury repair . The reprogramming and self-renewal ability of adult stem cells peaked in immature states from young individuals, then declined along with age in older ones, and these results are consistent with clinical observation. In the dental pulp, researches had demonstrated a better success rate of vital-pulp preservation by indirect and direct pulp capping in young age especially immature teeth with incomplete root development [6, 7]. DPSCs of immature teeth exhibited superior self-renewal and regenerative potential; however, the stemness may gradually reduce over time in mature ones [8, 9]. Exploring the regulatory network and molecular mechanisms of DPSCs’ biological activities will greatly advance therapeutic approaches in vital pulp therapy and also stem cell-based tissue engineering.
N6-methyladenosine (m6A) is the most prevalent internal modification of messenger RNA and non-coding RNA in eukaryotes [10, 11]. Though m6A modification has been identified for decades, the contributions remained a mystery; until recently, a flurry of studies revealed its functional effects on gene expressions and diverse biological processes . Reversible m6A modification can dynamically regulate almost every step in RNA metabolisms including mRNA stability, translation, non-coding RNA processing, and alternative splicing according to signal cues, and participate in multiple biological processes as homeostasis maintenance, stem cell fate determination, and also embryonic development [12,13,14]. In mammals, m6A marks are catalyzed by methyltransferase complex consisted of enzymatic subunit methyltransferase like 3 (METTL3) together with METTL14 and Wilms tumor 1-associated protein (WTAP) and can be demethylases by fat mass and obesity-associated protein (FTO) and alkB homolog 5 (ALKBH5). As the main “writer,” METTL3-mediated m6A modification was reported to modulate stem cell pluripotency, self-renewal, sex determination, and meiosis initiation [15,16,17]. Until now, little was known about the potential effect of m6A methylation and the molecular mechanism in DPSCs.
In this study, we aim to figure out the potential involvement of m6A modification in dental pulp hemostasis and elucidate the underlining molecular mechanism. Our study first demonstrated the features and distribution of m6A hallmarks in DPSCs by m6A RIP-seq and uncover a regulatory role of m6A methyltransferase-METTL3 in cell senescence and apoptosis. Furthermore, we discovered METTL3 make an effect on cell cycle control via a m6A-dependent manner which provide new insights into the RNA epigenetic mechanism in DPSCs.
This entire study was approved by the human research committee of Shunde Hospital, Southern Medical University (the First People’s Hospital of Shunde), Foshan, China. The Shunde Hospital Southern Medical University Medical Ethics Committee reviewed all protocols. All subjects were informed and performed according to the guidelines after written informed consents were obtained.
Dental pulp stem cells culture and characterization
Healthy dental pulp tissues from three independent donors were gently removed from the permanent teeth, and DPSCs were cultured as described previously . Briefly, pulp tissues were isolated, then digested in collagenase type I (3 mg/ml) for 1 h to obtain single-cell suspensions. Cells were seeded into a 6-cm dish and cultured with Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin at 37 °C under 5% CO2 condition. The medium was changed every 2 days, and DPSCs at passages 2–4 were used for the following experiments.
Flow cytometry was performed to identify DPSC phenotypes by screening the surface markers against CD29, CD34, CD44, CD45, and CD90 (BD Biosciences). Meanwhile, DPSCs were induced to differentiate into odontoblasts or adipocytes by an odontogenic or adipogenic medium as described in previous studies [19, 20]. After 7 days of induction by odontogenic medium, cells were fixed and stained with alkaline phosphatase (ALP), while 14 days for alizarin red staining (ARS). For adipogenic induction, cells were subjected to oil red O staining after 3 weeks of culture. For chondrogenic differentiation, DPSCs were cultured in an induction medium (Invitrogen) for 21 days before subjected to alcian blue staining to detect the extracellular matrix of chondrocytes.
m6A-methylated RNA immunoprecipitation sequencing
The m6A-methylated RNA immunoprecipitation sequencing (m6A RIP-Seq) was used to observe the characterizations of m6A modification in DPSCs. Total RNA was isolated with TRIzol Reagent and enriched by Oligo(dT)-attached magnetic beads. The purified m6A RIPed RNA fragments of DPSCs were collected with the Millipore kit according to the protocol; purified mRNA was then fragmentized into about 100 nt for sequencing using fragmentation buffer. The Magna ChIP Protein A/G Magnetic Beads were incubated with m6A antibody in immunoprecipitation buffer. RNA fragmentation was enriched by magnetic beads and subjected to deep sequencing on an Illumina Novaseq™ 6000 platform at the LC-BIO Bio-Tech Ltd. as described previously .
Real-time polymerase chain reaction
Total RNA from DPSCs were isolated with an RNeasy Mini kit (Foregene Biotec) and then reverse-transcribed with RT reagent kit (Takara Biotechnology) as described previously to obtain complementary DNA (cDNA). Real-time quantitative PCR was performed with Advanced™ SYBR Green Super Mix (Bio-Rad) according to the standard protocol. Data were analyzed using the standard curve method and normalized to GAPDH mRNA levels. The primer sequences used in real-time PCR were summerized in supplementary table 1.
To investigate specific target genes of METTL3, the m6A enrichment of PLK1 in m6A-IP RNA and input RNA as described before was quantified with quantitative PCR analysis for m6A RIP-qPCR.
The gene expression array GSE52853 was obtained from the GEO database. GSE64392 which was submitted by Tamaoki et al. analyzed the gene expression profiles in DPSCs from mature and immature human permanent teeth . In this study, GEO2R was used to identify the expression level of m6A-related genes for further study.
Lentiviral vectors construction and infection
Two independent biological replicates of shRNA lentiviral vectors were used in this article. Number#1, vectors encoding homo METTL3 gene and green fluorescent protein (GFP) were pursued from Santa Cruz Company; number#2 vectors encoding METTL3 and negative control were constructed by Shanghai Genechem Company. The overexpression lentiviral vectors of homo METTL3 gene and GFP were also obtained from Genechem Company.
Western blot analysis
The protein of DPSCs was lysed with RIPA buffer and assayed by the BCA protein assay kit (Thermo Scientific). Ten to thirty μg of protein in each group was subjected to SDS polyacrylamide gel electrophoresis and then transferred to polyvinylidene fluoride membranes (Millipore). The transferred proteins were reacted with primary antibody overnight at 4 °C and then labeled with secondary antibody for 1 h at room temperature. Primary antibodies in this study include METTL3 (1:1000; Cell Signaling Technology), GAPDH (1:2000; Proteintech Group), phosphorylation-p53 and p53 (1:1000; Santa Cruz), and PLK1 (1:1000; Cell Signaling Technology). The protein complexes were visualized with Super Signal-enhanced (Thermo Scientific) chemiluminescence.
DPSCs transfected with shMETTL3 or shCTR lentivirus were fixed and washed with PBS before β-galactosidase (β-gal) staining to detect cell senescence . Cells were subjected to staining solution for 16–24 h at 37 °C according to the β-gal staining kit (Beyotime Biotechnology). Six randomly chosen images were obtained, and then β-gal-positive cells were counted.
Flow cytometric analysis
DPSCs transfected with shRNA or overexpression lentivirus were plated in a 6-well plate. Then, cells were subjected to APC-conjugated Annexin V/propidium iodide assay (KeyGEN Biotechnology) according to a previous study and analyzed using flow cytometry .
Apoptosis was also evaluated with the TUNEL Apoptosis Assay Kit (KeyGEN Biotechnology). Briefly, DPSCs transfected with shRNA lentivirus were fixed and stained with TUNEL reaction mixture as protocol. Six randomly fluorescence images were obtained, and percentage of red-labeled TUNEL-positive cell was analyzed by ImageJ.
RNA sequencing and profile analysis
For sequencing analysis, total RNA was extracted by TRIzol Reagent and purified with Oligo(dT)-attached magnetic beads. The purified RNA fragments were converted into double-stranded cDNA. Paired-end runs with 250–350 base pair (bp) read lengths were used for RNA sequencing (RNA-seq) by Novaseq™ 6000 following the manufacturer’s recommendations. FastQC and FASTX toolkit were used to control data quality and map to homo reference genomes by HISAT2. Differentially expressed genes with fold change ≥2.0 and p value < 0.05 were subjected to analysis.
The enrichment of Gene Ontology (GO) was performed with R package GOseq and the DAVID database. Then, we mapped the differentially expressed genes to the Search Tool for Retrieval of Interacting Genes (STRING) database to obtain the protein-protein interactions, and the networks were displayed by Cytoscape 3.5.1.
Cell cycle analysis
METTL3 knockdown and overexpression vectors were transduced in DPSCs, wihle PLK1 inhibitor BI2536 (1 nM, MedChemExpress)  or DMSO as negative control were used to inhibit PLK1 expression in shMETTL3-DPSCs. DPSCs were washed with PBS and then fixed in 70% cold ethanol overnight. At least 50,000 cells were subjected to propidium iodide staining with a Cell Cycle and Apoptosis Analysis Kit (Beyotime), and the DNA content was measured by fluorescence-activated cell sorting (FACS) instrument (BD Biosciences).
Experiments were carried out at least three times, and the data was presented as mean ± standard deviation. Statistical significance (p < 0.05) level was determined by analyzing comparisons of variance using Graphpad Prism 7.0.
Characterizations and m6A modification features of DPSCs
DPSCs appeared as spindle-like, elongated fibroblastic shape (Supplementary Figure S1). The surface markers of DPSCs were analyzed by flow cytometry; mesenchymal stem cell markers such as CD29, CD44, and CD90 were highly expressed with over 90% positive cells, while the expression of CD34 and CD45 (hematopoietic cell antigen) can be barely detected (Fig. 1a). After odontogenic medium induction for 3, 7, and 14 days, the early differentiation stage of DPSCs was evidenced by strong alkaline phosphatase (ALP) straining and upregulated osteo/odontogenesis-related markers such as RUNX2, ALP, and COL1 (Fig. 1b, c) while late differentiation stage identified by significant calcium mineralization in alizarin red staining (ARS) and mRNA expression of BSP, DSPP, and OCN (Fig. 1d, e). Meanwhile, accumulation of lipid vacuoles was detected by oil O staining in DPSCs after adipogenic induction, accompanying with increased expression of adipogenic-related markers as CREB, LPL (Fig. 1f, g). After chondrogenic differentiation, DPSCs were positive for alcian blue staining which indicated the glycosaminoglycan formation (Fig. 1h). These data suggested that DPSCs were capable of differentiating into odontoblasts and adipocytes which verified multiple differentiation potential and regenerative capacity did exist in DPSCs.
As the most prevalent RNA modification, m6A marks can dynamically regulate RNA metabolisms and participate in stem cell fate determination. To explore the binding area, motif, and distribution of m6A modifications in DPSCs, m6A immunoprecipitation with deep sequencing (m6A RIP-seq) was performed. We identified high similarity and overlapping of m6A peaks in two independent replicates and significantly enriched in GGm6ACA motif which is consistent with previous reports (Fig. 2a). A total of 13,191 m6A-modified transcripts were identified in m6A RIP-seq (p value < 0.05, fold enrichment> 1). The m6A peaks were abundant in 3′ untranslated regions (3′ UTR) 45.25%, exons 31.76%, 5′ UTR 22.99% (Fig. 2b, c), and transcription factor-binding sites were predominantly distributed within 100 kb related to transcription start sites (TSS) (Fig. 2d).
METTL3 knockdown induced cell apoptosis and senescence of DPSCs
The dental pulp in immature teeth with incomplete root development preserved better regenerative potential than mature teeth . We further investigated the gene expressions related to m6A modification including “writers” (METTL3, METTL14, WTAP), “erasers” (FTO, ALKBH5), and “readers” YTH domain-containing family protein (YTHDF1, YTHDF2, YTHDF3). The high expression levels of “writers” METTL3 and METTL14 in immature dental pulp tissues were observed compared with mature ones (Fig. 2e). In the global transcriptome profile of DPSCs in different developmental stages, the relative gene expressions of METTL3, WTAP, YTHDF1, and YTHDF2 were reduced in mature DPSCs compared to immature ones, while others had no significant differences (Fig. 2f).
As the only methyltransferase with catalyzing activity, METTL3 showed a remarkable high expression level in immature DPSCs. To investigate the biological effect of METTL3 inhibition in mature DPSCs, two independent lentiviral vectors were constructed to knockdown METTL3 expression. Both METTL3 shRNA-1 and shRNA-2 could inhibit its mRNA, protein expression in DPSCs (Fig. 3a, S2A), and so does lentiviral vector-overexpressing METTL3 (Fig. 3d, S2B). Representative images of cell morphology showed that shMETTL3-DPSCs exhibited elongate spindle shape and enlarged size compared with the control group (Fig. 3b). Also, β-galactosidase (β-gal) staining detected an increasing number of β-gal-positive cells after METTL3 inhibition which indicated DPSCs senescence (Fig. 3c, d). CCK8 assay showed no significant differences in cell viability in DPSCs (Fig. 3f) while flow cytometry demonstrated a significantly higher apoptosis rate after METTL3 knockdown (Fig. 3g). Additionally, phosphorylation of apoptotic-related p53 protein was upclined, and TUNEL assay showed increasing apoptosis by METTL3 knockdown (Fig. 3h, i). These data suggested that METTL3 inhibition in immature DPSCs induced cell senescence and apoptosis which could relate to the reduced self-renewal and regenerative potential.
Bioinformatic analysis of METTL3 knockdown in DPSCs
RNA-seq was carried out to screen the transcriptome profile of DPSCs transfected with two independent shMETTL3 vectors. A list of total 312 genes with GeneBank accession number from the National Center for Biotechnology Information (NCBI) database were identified as differentially expressed genes (DEGs) in DPSCs (p value < 0.05, FC > 2 or < − 2). Among these genes, 130 genes were upregulated while 182 ones were downregulated (Fig. 4a).
Bioinformatic analysis combing 13,191 m6A-methylated genes identified in m6A RIP-seq and 312 DEGs in shMETLL3 RNA-seq revealed that 216 genes were overlapped which might be regulated by METTL3-mediated m6A RNA methylation (Fig. 4b). Gene Ontology Consortium (GO) analysis was performed to functionally enrich the overlapped genes into specific biological contexts and summarize genes of related functions. For cell component (CC), genes were mainly enriched in the plasma membrane nucleus, while genes exhibited significant enrichments in protein binding for molecular function (MF) (Fig. 4c). Biological processes (BP) were essential to evaluate stem cell activities. Cell cycle, cell division, and positive regulation of transcription were the top 3 related processes which were related to the self-renewal and regenerative capacity of stem cells (Fig. 4c, d). In summary, the alternative gene expressions enriched in items as cell cycle, chromosome which indicated that cell cycle was the most relevant function resulted from METTL3 knockdown in DPSCs.
Alteration of METTL3 disrupted cell cycle via PLK1 m6A modulation
Flow cytometric analysis was performed to analyze the biological role of METTL3 in DPSCs cycle distribution. METTL3 knockdown by both shMETTL3 lentiviral vectors in DPSCs showed upregulated in the percentage of S phase while no significant difference was found in G2-M phase comparing to the negative control group which suggested METTL3 knockdown contributed to S phase cycle arrest in DPSCs (Fig. 5a, b). METLL3 overexpression led to different cycle distribution with a lightly upregulated in G2-M phase, and no variant in both G0-G1 and S phase (Fig. 5c, d).
To further explore the underlying mechanism of the impaired cell cycle regulated by METTL3, protein-protein interactions (PPIs) of DEGs were screened by the STRING database and network was then subjected to Cytoscape for further analysis (Fig. 6a). A critical cycle and mitosis modulator, Polo-like kinase 1 (PLK1), was identified as the central node of PPI network. qPCR and Western blot were used to confirm PLK1 expression after METTL3 alteration. mRNA and protein expression level of PLK1 were significantly upregulated after METTL3 knockdown while downregulated in METTL3 overexpression (Fig. 6b, c). Data visualization of m6A RIP-seq displayed the specific m6A peaks and distribution in PLK1 and m6A RIP-qPCR demonstrated m6A modification level fallen off a cliff after METTL3 inhibition in DPSCs (Fig. 6d, e). In addition, PLK1 inhibitor BI2536 partly rescued the S phase cycle arrest induced by METTL3 knockdown (Fig. 6f). These data indicated METTL3-meditated PLK1 m6A modification regulated the cell cycle progression of DPSCs.
As the most prevalent, reversible internal modification in nuclear RNAs, m6A modifications are involved in various biological processes including stem cell fate determination, embryonic development, and cell cycle control [25,26,27]. m6A methyltransferase complex and demethylases cooperated to influence different life stages of stem cells in post-translational levels . DPSCs served as a promising cell source in endodontic regeneration and tissue engineering; however, its regenerative capacity is limited and reduced in a long lifespan. Our data first identified the binding area, motif, and distribution of m6A peaks in DPSCs by m6A RIP-seq. The expression level of m6A methyltransferases-METTL3 was relatively high in immature DPSCs which are known to possess superior regenerative potential. METTL3 inhibition in immature DPSCs induced cell senescence and apoptosis. A conjoint analysis of RNA-seq and m6A RIP-seq showed METTL3 knockdown mainly associated cell cycle and chromosome stability that tightly related to stem cell self-renewal and fate decision. Mechanistically, the alteration of METTL3 led to cell cycle arrest in DPSCs by regulating m6A modification of PLK1. Our data indicated that METTL3 may support sufficient regenerative potential in DPSCs by cell cycle control.
Over 80% of transcripts were methylated during development which occur to effect the self-renewal and pluripotency of stem cells [28, 29]. m6A modification can influence stability, altered splicing, exporting, translation of mRNA in a fine, and cell-specific manner. All m6A marks share the redundant consensus RRm6ACH ([G/A/U][G>A]m6AC[U>A>C]) motif which is enriched in 3′ and 5′ UTR of RNA transcripts especially long internal exons . In our study, m6A RIP-seq revealed that in DPSCs, the consensus motif of m6A was GGACAG and mainly distributed in 3′ UTR, then 5′ UTR and the methylated level were higher in exons than introns which was in consistent with previous studies. The clustered m6A modifications in 3′ UTR around the stop codon influenced stability, localization, and translation of RNA while in 5′ UTR effected translation in a cap-independent manner [14, 31]. Methyltransferase complex is responsible to distinguish and target m6A sites, and METTL3 is the only catalytically active subunit which catalyzes the methylation of adenosines to form m6A modification . METTL3 can independently read, regulate m6A marks near 3′ UTR of mRNA and promote translation of this specific mRNA . Here, we identified METTL3 was the essential m6A-related gene differentially expressed in different stages of DPSCs which might relate to regenerative capacity alteration.
METTL3-mediated m6A methylation modified the transcripts in a cell-dependent manner which lead to divergent effects, somehow even conflicting . METTL3 knockdown and m6A depletion could limit self-renewal and promote cell differentiation in epiblast stem cells while improving pluripotency, block regeneration in embryonic stem cells [13, 26, 34]. A recent study suggested METTL3 knockout accelerated mesenchymal stem cell senescence by modulating m6A modification of MIS12 which supported our data . In our study, METTL3 inhibition in immature DPSCs promoted cell senescence and apoptosis which might further result in reduced regenerative potential in mature teeth. METTL3 depletion could make a contrary effect in apoptosis in different states. Some studies reported that knockdown METTL3 could rescue cell apoptosis induced by high glucose in human lens epithelial cells and also hypoxia/reoxygenation-treated cardiomyocytes [36, 37]. Meanwhile, METTL3 conditional knockout induced apoptosis of newborn granule cells which led to manifest cerebellar hypoplasia . Furthermore, METTL3 inhibition promoted cell apoptosis by modulating m6A modification in multiple cancers [17, 39]. We found that knockdown METTL3 induced DPSC apoptosis along with p53 pathway activation which could impair stem cell self-renewal and fate decision.
Functional stem cell senescence and apoptosis are cooperated with proper cell cycle progress which make essential effect in cellular homeostasis. Prime stem cells usually stayed quiescence as a cell resource; progenitor cell expansion needs a high cycling rate or shorter time, while some differentiated cells withdraw from the cell cycle could still reenter under certain signals [40, 41]. Dynamic activation of the cell cycle was responsible for bone marrow stem cell self-renewal and hematopoietic stem cell maturation . The transcriptome and epigenetic landscape alterations in cell cycle correspond with different status in development, maturation, and differentiation . m6A modifications were reported to participate in stem cell self-renew and pluripotency by cell cycle control [29, 43, 44]. In cortical neurogenesis, the reduction of m6A modification resulted in a longer cell cycle which related to insufficient differentiation and delayed generation . By bioinformatics analysis, we found inhibition of METTL3-meditated m6A mainly enriched in cycle-mediated factors and METTL3 alteration in DPSCs led to the cell cycle arrest in DPSCs. These findings were supported by recent studies which indicated a regulatory role of METTL3-meditated m6A methylation in cell cycle control. Methyltransferase complex consisted of METTL3, METTL14, and WTAP, inhibition of either could contribute to cell cycle arrest and disrupt the adipogenesis of adipocytes . As an oncogene, METTL3 depletion could result in cell cycle arrest in acute myeloid leukemia and bladder cancer in a m6A-dependent manner [46, 47]. METTL3 was also reported to be a tumor suppressor by regulating cell cycle and self-renewal in renal cell carcinoma and glioblastoma [48, 49].
Cell cycle progress is controlled by robust cell cycle-related genes and checkpoints which were essential for self-renewal and homeostasis of stem cells. As a critical mitotic cell cycle regulator, PLK1 make effect on centrosome maturation, cytokinesis execution, bipolar spindle formation, and chromosome segregation. The expression of PLK1 was tightly associated with cell cycle stages that it initiated at a relatively low level and accumulated during the S phase, then reach a peak in the G2/M phase [50, 51]. PLK1 was widely studied as a target in cancer therapy while residued relative low expression level in normal cells. PLK1 overexpression in mammary epithelial and murine cells contributed to mitosis defection and aberration in chromosome segregation, cytokinesis which supported unscheduled PLK1 expression disrupted mitosis, and defected genomic stability [52, 53]. Consistent with this study, we discovered the expression of PLK1 in DPSCs and elevated PLK1 expression following METTL3 inhibition resulted in the S phase arrested which suggested a regulatory role of PLK1 in cell cycle control of DPSCs. PLK1 was reported to mediate G2 And S-phase expressed 1 (GTSE1) phosphorylation and p53 inactivation which led to G2 checkpoint recovery, DNA damage checkpoints termination, and chromosome instability . In this study, our data showed that METTL3 regulated PLK1 expression in a m6A-dependent manner and participated in cell cycle control, which further effect DPSC senescence and apoptosis.
Taken together, we demonstrated that m6A modification and METTL3 alteration participated in different stages of DPSCs which would relate to stem cell fate determination. Impaired METTL3 expression in DPSCs led to increasing cell senescence and apoptosis by interfering with the mitotic cell cycle in a m6A-dependent manner. Our results first uncovered the m6A methylated hallmarks in DPSCs and revealed the potential link between METTL3-mediated m6A modification and cell cycle control. Further researches are still needed to confirm these findings in vivo and enrich the RNA epigenetic mechanisms in pluripotency and self-renewal of DPSCs which would advance our understanding and serve new strategies in endodontic regeneration and tissue engineering.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Dental pulp stem cells
- m6A RIP-seq:
m6A immunoprecipitation with deep sequencing
Methyltransferase like 3
Polo-like kinase 1
Methyltransferase like 14
Wilms tumor 1-associated protein
Fat mass and obesity-associated protein
alkB homolog 5
Alizarin red staining
- m6A RIP-Seq:
m6A-methylated RNA immunoprecipitation sequencing
Real-time polymerase chain reaction
Retrieval of Interacting Genes
YTH domain-containing family protein
National Center for Biotechnology Information
Transcription start sites
G2 And S-phase expressed 1
Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, et al. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318(5858):1917–20.
Sui B, Chen C, Kou X, Li B, Xuan K, Shi S, et al. Pulp stem cell-mediated functional pulp regeneration. J Dent Res. 2019;98(1):27–35.
Huang GT, Gronthos S, Shi S. Mesenchymal stem cells derived from dental tissues vs. those from other sources: their biology and role in regenerative medicine. J Dent Res. 2009;88(9):792–806.
Gronthos S, Brahim J, Li W, Fisher LW, Cherman N, Boyde A, et al. Stem cell properties of human dental pulp stem cells. J Dent Res. 2002;81(8):531–5.
Oh J, Lee YD, Wagers AJ. Stem cell aging: mechanisms, regulators and therapeutic opportunities. Nat Med. 2014;20(8):870–80.
Aguilar P, Linsuwanont P. Vital pulp therapy in vital permanent teeth with cariously exposed pulp: a systematic review. J Endod. 2011;37(5):581–7.
Lipski M, Nowicka A, Kot K, Postek-Stefanska L, Wysoczanska-Jankowicz I, Borkowski L, et al. Factors affecting the outcomes of direct pulp capping using Biodentine. Clin Oral Investig. 2018;22(5):2021–9.
Tamaoki N, Takahashi K, Aoki H, Iida K, Kawaguchi T, Hatakeyama D, et al. The homeobox gene DLX4 promotes generation of human induced pluripotent stem cells. Sci Rep. 2014;4:7283.
Macrin D, Alghadeer A, Zhao YT, Miklas JW, Hussein AM, Detraux D, et al. Metabolism as an early predictor of DPSCs aging. Sci Rep. 2019;9(1):2195.
Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18(1):31–42.
Narayan P, Rottman FM. An in vitro system for accurate methylation of internal adenosine residues in messenger RNA. Science. 1988;242(4882):1159–62.
Zhang C, Chen Y, Sun B, Wang L, Yang Y, Ma D, et al. m6A modulates haematopoietic stem and progenitor cell specification. Nature. 2017;549(7671):273–6.
Geula S, Moshitch-Moshkovitz S, Dominissini D, Mansour AA, Kol N, Salmon-Divon M, et al. Stem cells. m6A mRNA methylation facilitates resolution of naive pluripotency toward differentiation. Science. 2015;347(6225):1002–6.
Zhou J, Wan J, Shu XE, Mao Y, Liu XM, Yuan X, et al. N6-methyladenosine guides mRNA alternative translation during integrated stress response. Mol Cell. 2018;69(4):636–47.e7.
Xu K, Yang Y, Feng GH, Sun BF, Chen JQ, Li YF, et al. Mettl3-mediated m6A regulates spermatogonial differentiation and meiosis initiation. Cell Res. 2017;27(9):1100–14.
Wu Y, Xie L, Wang M, Xiong Q, Guo Y, Liang Y, et al. Mettl3-mediated m6A RNA methylation regulates the fate of bone marrow mesenchymal stem cells and osteoporosis. Nat Commun. 2018;9(1):4772.
Vu LP, Pickering BF, Cheng Y, Zaccara S, Nguyen D, Minuesa G, et al. The N6-methyladenosine (m6A)-forming enzyme METTL3 controls myeloid differentiation of normal hematopoietic and leukemia cells. Nat Med. 2017;23(11):1369–76.
Luo H, Wang C, Liu M, Yin B, AP, Huang D, et al. Inhibition of SOX9 promotes inflammatory and immune responses of dental pulp. J Endod. 2018;44(5):792–9.
Farahzadi R, Fathi E, Vietor I. Mesenchymal stem cells could be considered as a candidate for further studies in cell-based therapy of Alzheimer’s disease via targeting the signaling pathways. ACS Chem Neurosci. 2020;11(10):1424–35.
Liu W, Luo H, Wang R, Kang Y, Liao W, Sun Y, et al. Rapamycin-induced autophagy promotes the chondrogenic differentiation of synovium-derived mesenchymal stem cells in the temporomandibular joint in response to IL-1β. BioMed research international. 2020;2020:4035306.
Wang Y, Zheng Y, Guo D, Zhang X, Guo S, Hui T, et al. m6A methylation analysis of differentially expressed genes in skin tissues of coarse and fine type Liaoning cashmere goats. Front Genet. 2019;10:1318.
Farahzadi R, Fathi E, Mesbah-Namin SA, Zarghami N. Anti-aging protective effect of L-carnitine as clinical agent in regenerative medicine through increasing telomerase activity and change in the hTERT promoter CpG island methylation status of adipose tissue-derived mesenchymal stem cells. Tissue Cell. 2018;54:105–13.
Fathi E, Vietor I. Mesenchymal stem cells promote caspase expression in Molt-4 leukemia cells via GSK-3α/β and ERK1/2 signaling pathways as a therapeutic strategy. Curr Gene Ther. 2021;21(1):81-8.
Ueda A, Oikawa K, Fujita K, Ishikawa A, Sato E, Ishikawa T, et al. Therapeutic potential of PLK1 inhibition in triple-negative breast cancer. Lab Invest. 2019;99(9):1275–86.
Horiuchi K, Kawamura T, Iwanari H, Ohashi R, Naito M, Kodama T, et al. Identification of Wilms’ tumor 1-associating protein complex and its role in alternative splicing and the cell cycle. J Biol Chem. 2013;288(46):33292–302.
Batista PJ, Molinie B, Wang J, Qu K, Zhang J, Li L, et al. m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells. Cell Stem Cell. 2014;15(6):707–19.
Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, et al. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014;10(2):93–5.
Zhao BS, He C. Fate by RNA methylation: m6A steers stem cell pluripotency. Genome Biol. 2015;16:43.
Yoon KJ, Ringeling FR, Vissers C, Jacob F, Pokrass M, Jimenez-Cyrus D, et al. Temporal control of mammalian cortical neurogenesis by m6A methylation. Cell. 2017;171(4):877–89.e17.
Heck AM, Wilusz CJ. Small changes, big implications: the impact of m6A RNA methylation on gene expression in pluripotency and development. Biochim Biophys Acta Gene Regul Mech. 2019;1862(9):194402.
Meyer KD, Saletore Y, Zumbo P, Elemento O, Mason CE, Jaffrey SR. Comprehensive analysis of mRNA methylation reveals enrichment in 3′ UTRs and near stop codons. Cell. 2012;149(7):1635–46.
Wang X, Feng J, Xue Y, Guan Z, Zhang D, Liu Z, et al. Structural basis of N6-adenosine methylation by the METTL3-METTL14 complex. Nature. 2016;534(7608):575–8.
Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, et al. N6-methyladenosine modulates messenger RNA translation efficiency. Cell. 2015;161(6):1388–99.
Wang Y, Li Y, Toth JI, Petroski MD, Zhang Z, Zhao JC. N6-methyladenosine modification destabilizes developmental regulators in embryonic stem cells. Nat Cell Biol. 2014;16(2):191–8.
Wu Z, Shi Y, Lu M, Song M, Yu Z, Wang J, et al. METTL3 counteracts premature aging via m6A-dependent stabilization of MIS12 mRNA. Nucleic Acids Res. 2020;48(19):11083–96.
Yang J, Liu J, Zhao S, Tian F. N6-methyladenosine METTL3 modulates the proliferation and apoptosis of lens epithelial cells in diabetic cataract. Mol Ther Nucleic Acids. 2020;20:111–6.
Song H, Feng X, Zhang H, Luo Y, Huang J, Lin M, et al. METTL3 and ALKBH5 oppositely regulate m6A modification of TFEB mRNA, which dictates the fate of hypoxia/reoxygenation-treated cardiomyocytes. Autophagy. 2019;15(8):1419–37.
Wang CX, Cui GS, Liu X, Xu K, Wang M, Zhang XX, et al. METTL3-mediated m6A modification is required for cerebellar development. Plos Biol. 2018;16(6):e2004880.
Wang H, Xu B, Shi J. N6-methyladenosine METTL3 promotes the breast cancer progression via targeting Bcl-2. Gene. 2020;722:144076.
Orford KW, Scadden DT. Deconstructing stem cell self-renewal: genetic insights into cell-cycle regulation. Nat Rev Genet. 2008;9(2):115–28.
Boheler KR. Stem cell pluripotency: a cellular trait that depends on transcription factors, chromatin state and a checkpoint deficient cell cycle. J Cell Physiol. 2009;221(1):10–7.
Zape JP, Lizama CO, Cautivo KM, Zovein AC. Cell cycle dynamics and complement expression distinguishes mature haematopoietic subsets arising from hemogenic endothelium. Cell Cycle. 2017;16(19):1835–47.
Hirayama M, Wei FY, Chujo T, Oki S, Yakita M, Kobayashi D, et al. FTO demethylates cyclin D1 mRNA and controls cell-cycle progression. Cell Rep. 2020;31(1):107464.
Abakir A, Giles TC, Cristini A, Foster JM, Dai N, Starczak M, et al. N6-methyladenosine regulates the stability of RNA:DNA hybrids in human cells. Nat Genet. 2020;52(1):48–55.
Kobayashi M, Ohsugi M, Sasako T, Awazawa M, Umehara T, Iwane A, et al. The RNA Methyltransferase Complex of WTAP, METTL3, and METTL14 Regulates Mitotic Clonal Expansion in Adipogenesis. Molecular and cellular biology. 2018;38(16):e00116-18.
Barbieri I, Tzelepis K, Pandolfini L, Shi J, Millán-Zambrano G, Robson SC, et al. Promoter-bound METTL3 maintains myeloid leukaemia by m6A-dependent translation control. Nature. 2017;552(7683):126–31.
Han J, Wang JZ, Yang X, Yu H, Zhou R, Lu HC, et al. METTL3 promote tumor proliferation of bladder cancer by accelerating pri-miR221/222 maturation in m6A-dependent manner. Mol Cancer. 2019;18(1):110.
Li X, Tang J, Huang W, Wang F, Li P, Qin C, et al. The M6A methyltransferase METTL3: acting as a tumor suppressor in renal cell carcinoma. Oncotarget. 2017;8(56):96103–16.
Cui Q, Shi H, Ye P, Li L, Qu Q, Sun G, et al. m6A RNA methylation regulates the self-renewal and tumorigenesis of glioblastoma stem cells. Cell Rep. 2017;18(11):2622–34.
Bassermann F, Frescas D, Guardavaccaro D, Busino L, Peschiaroli A, Pagano M. The Cdc14B-Cdh1-Plk1 axis controls the G2 DNA-damage-response checkpoint. Cell. 2008;134(2):256–67.
Gheghiani L, Loew D, Lombard B, Mansfeld J, Gavet O. PLK1 activation in late G2 sets up commitment to mitosis. Cell Rep. 2017;19(10):2060–73.
Kong D, Farmer V, Shukla A, James J, Gruskin R, Kiriyama S, et al. Centriole maturation requires regulated Plk1 activity during two consecutive cell cycles. J Cell Biol. 2014;206(7):855–65.
de Cárcer G, Venkateswaran SV, Salgueiro L, El Bakkali A, Somogyi K, Rowald K, et al. Plk1 overexpression induces chromosomal instability and suppresses tumor development. Nat Commun. 2018;9(1):3012.
Cholewa BD, Liu X, Ahmad N. The role of polo-like kinase 1 in carcinogenesis: cause or consequence? Cancer Res. 2013;73(23):6848–55.
Thanks for all the support and contributions of participators and volunteers.
This work was supported by the National Natural Science Foundation of China (81900989, 81870786), China Postdoctoral Science Foundation (2020M672548, 2019M662993), and Medical Scientific Research Foundation of Guangdong Province, China (A2020227).
Ethics approval and consent to participate
This entire study was approved by the Human Research Committee of Shunde Hospital, Southern Medical University (the First People’s Hospital of Shunde), Foshan, China (scientific research ethical code 20200824). All subjects were informed and performed under the supervision of the Shunde Hospital, Southern Medical University Medical Ethics Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The sequence of primers used in PCR.
The images of DPSCs under microscope.
Quantitative analysis of protein expression bands after METTL3 inhibition and expression.
About this article
Cite this article
Luo, H., Liu, W., Zhang, Y. et al. METTL3-mediated m6A modification regulates cell cycle progression of dental pulp stem cells. Stem Cell Res Ther 12, 159 (2021). https://doi.org/10.1186/s13287-021-02223-x