- Open Access
METTL1-mediated m7G methylation maintains pluripotency in human stem cells and limits mesoderm differentiation and vascular development
Stem Cell Research & Therapy volume 11, Article number: 306 (2020)
7-Methylguanosine (m7G) is one of the most conserved modifications in nucleosides within tRNAs and rRNAs. It plays essential roles in the regulation of mRNA export, splicing, and translation. Recent studies highlighted the importance of METTL1-mediated m7G tRNA methylome in the self-renewal of mouse embryonic stem cells (mESCs) through its ability to regulate mRNA translation. However, the exact mechanisms by which METTL1 regulates pluripotency and differentiation in human induced pluripotent stem cells (hiPSCs) remain unknown. In this study, we evaluated the functions and underlying molecular mechanisms of METTL1 in regulating hiPSC self-renewal and differentiation in vivo and in vitro.
By establishing METTL1 knockdown (KD) hiPSCs, gene expression profiling was performed by RNA sequencing followed by pathway analyses. Anti-m7G northwestern assay was used to identify m7G modifications in tRNAs and mRNAs. Polysome profiling was used to assess the translation efficiency of the major pluripotent transcription factors. Moreover, the in vitro and in vivo differentiation capacities of METTL1-KD hiPSCs were assessed in embryoid body (EB) formation and teratoma formation assays.
METTL1 silencing resulted in alterations in the global m7G profile in hiPSCs and reduced the translational efficiency of stem cell marker genes. METTL1-KD hiPSCs exhibited reduced pluripotency with slower cell cycling. Moreover, METTL1 silencing accelerates hiPSC differentiation into EBs and promotes the expression of mesoderm-related genes. Similarly, METTL1 knockdown enhances teratoma formation and mesoderm differentiation in vivo by promoting cell proliferation and angiogenesis in nude mice.
Our findings provided novel insight into the critical role of METTL1-mediated m7G modification in the regulation of hiPSC pluripotency and differentiation, as well as its potential roles in vascular development and the treatment of vascular diseases.
Somatic cells can be transformed into induced pluripotent stem cells (iPSCs) functionally equivalent to embryonic stem cells (ESCs) by simultaneously introducing four transcription factors (TFs) Oct4, Sox2, Klf4, and c-Myc [1,2,3]. The critical role of these factors in maintaining pluripotency is highlighted by the fact that cell differentiation is enhanced when these pluripotent factors are suppressed [4, 5]. Recent studies have led to the identification of additional TFs and microRNAs (miRNAs) that can also affect the pluripotency and self-renewal capacity of iPSCs [6, 7].
Now, several studies have indicated that epigenetic modifications regulate stem cell fate decision and development [8,9,10]. One of the hotspots in epigenetic regulation, N6-methyladenosine (m6A) mRNA methylation, has been observed during embryogenesis and is believed to be involved in the regulation of stem cell development . However, the role of post-transcriptional modifications of RNA, including 7-methylguanosine (m7G), in the regulation of hiPSCs fate remains unclear. METTL1 is an m7G RNA methyltransferase, and mutations in the human m7G methyltransferase complex METTL1/WDR4 can lead to primordial dwarfism and brain malformation . Our previous study found that METTL1-mediated m7G tRNA modification regulates mRNA translation and the self-renewal, pluripotency, and neural lineage differentiation in mouse embryonic stem cells (mESCs) . However, the differentiation of hiPSC is accompanied by the translational regulation of multiple genes [14, 15]. Thus, we hypothesize that METTL1 may also play an important role in hiPSC pluripotency and differentiation.
In this study, we investigated whether METTL1-mediated m7G methylation regulates pluripotency of hiPSCs. We find that METTL1 knockdown (KD) significantly decreased cell viability and changed the cell cycle of hiPSCs. Moreover, METTL1 KD significantly accelerates hiPSC differentiation toward the mesoderm fate and suppressed neuroectoderm differentiation in vitro. Most importantly, compared with control hiPSCs, METTL1-KD hiPSCs significantly enhance teratoma formation in vivo by promoting cell proliferation and angiogenesis in nude mice. Therefore, our study revealed that METTL1-mediated m7G modification regulates hiPSC pluripotency and differentiation, providing novel insights into vascular development and the treatment of vascular diseases.
Materials and methods
Chemical reagents and antibodies
HiPSCs and hESCs were maintained in matrigel-coated plates in mTeSR1 medium and under standard cell culture conditions. Cells were passaged every 4 to 5 days. 293T cells were maintained in DMEM supplemented with 10% FBS. All cells were maintained at 37 °C in a humidified 5% CO2 incubator. The detailed list of cell lines and mediums used for cell cultures was provided in the Supplementary Table S1 and Table S3.
Knockdown of METTL1 in hiPSCs with shRNA lentivirus
The lentivirus expressing METTL1-targeting shRNAs has been previously described . Briefly, 293T cells were transfected with the appropriate lentiviral vector using the Lipofectamine 3000 reagent (Invitrogen) according to the manufacturer’s instructions. The pLKO.1-shMETTL1-1 and pLKO.1-shMETTL1-5 vectors were used for METTL1 knockdown. The target and packaging plasmids TAT, HEPM2, RAII, and VSVG were used at a 4:1:1:1:1 ratio. Cells transduced with pLKO.1-shGFP viruses were used for control. Cell transduction was performed using polybrene. The transduced cells were selected with 1 μg/ml puromycin for 5 days.
HiPSC self-renewal and proliferation analysis
To determine the proliferation capacity of hiPSCs, cells (2 × 104) were seeded in 6-well plates, and the cell numbers were quantified on days 3, 5, and 7. For the colony formation assay, 500 cells were seeded in 6-well plates, followed by colony size and number quantification on day 7. Alkaline phosphatase staining was performed using an alkaline phosphatase detection kit (Merck Millipore, SCR004) following the manufacturer’s instructions. Briefly, 4–5 days after culturing cells on coverslips, cells were fixed with 4% paraformaldehyde for 1 min at room temperature. After washing with PBS, the cells were incubated with working buffer 15–45 min at room temperature and protected from light, followed by microscopy.
RT-PCR and quantitative PCR analysis of gene expression
Total RNA was isolated using TRIzol, and cDNA synthesis was performed using ~ 600 ng of RNA and the all-in-one RT MasterMix kit (TransGen Biotech). For the qRT-PCR, the Step One Plus PCR system and SYBR Green Master Mix were used. All the relative gene expression was calculated using the function of comparative cycle threshold (Ct) and analyzed by 2−ΔΔCt method. The primer sequences used in this study are listed in the Supplementary Table S4.
Cells were lysed in RIPA buffer containing protease inhibitors, and ~ 40 μg of protein was resolved by 10% SDS-PAGE. Proteins were then transferred onto PVDF membranes and blocked with 5% BSA in TBS-T. Membranes were probed overnight at 4 °C with primary antibodies against METTL1, FLAG, OCT4, SOX2, Nanog, CyclinD1, CyclinD2, CDK2, and ACTB. After washing, the membranes were incubated with the appropriate secondary antibody, and protein bands were visualized using the ECL system (Thermo Scientific).
Determination of cell cycle analysis by flow cytometry
Cells were harvested using Accutase and fixed in ice-cold 70% ethanol overnight. Cells were stained with PI for 30 min, followed by flow cytometry analysis (Beckman Coulter, MoFlo Astrios EQ). Flow cytometry data were analyzed using Summit v. 5.3.
Determination of m7G modification by Northwestern blot
Northwestern blot was performed to assess m7G modifications among METTL1-KD and control hiPSCs. RNA (~ 2 μg) was incubated in 2× TBE loading buffer at 95 °C for 5 min and then loaded onto 10% polyacrylamide-TBE urea gels. After electrophoresis, RNAs were transferred onto positively charged nylon membranes and UV-crosslinked. Membranes were probed with anti-m7G antibodies, and signal was detected using the ECL system.
Determination of translation efficiency by sucrose-gradient centrifugation and polysome fractionation
Sucrose-gradient centrifugation and polysome fractionation were used to assess the impact of METTL1 silencing on the translation efficiency of stem cell markers. Cells (2 × 107) were lysed in polysome cell extraction buffer on ice (the components are listed in the Supplementary Table S5), followed by centrifugation at 13,000×g for 15 min at 4 °C. Samples were loaded onto a 10–50% sucrose gradient and centrifuged at 36,000×g and 4 °C for 2.5 h. The polysome fraction was collected for qRT-PCR analysis. The relative translation efficiency was calculated by normalizing the mRNA levels in polysome fraction to the mRNA of gene expression.
Induction of embryoid bodies (EBs)
To estimate the role of METTL1 in hiPSC differentiation in vitro, EB formation assay was performed as previously described [16, 17]. Briefly, hiPSCs were digested with Accutase, and 4 × 103 cells/well were seeded in 96-well spheroid microplates (Corning, Cat #4515) in KO-DMEM medium supplemented with 10% KSR, 1% NEAA, 1 mM L-Glutamine, 50 μM 2-ME, and 10 μM Y-27632. Cells were incubated for 6 days, and the medium was refreshed daily.
Induction of teratoma formation
To evaluate the role of METTL1 in hiPSC differentiation in vivo, animal experiments were performed in accordance with the guidelines provided by the First Affiliated Hospital of Sun Yat-sen University. Four-week-old BALB/C nude male mice were purchased from GemPharmatech Co. Ltd. (Nanjing, China). Mice were subcutaneously injected in the flank with 200 μl of PBS containing hiPSCs (6 × 106 cells). At 6 weeks post-injection, teratomas were extracted from the nude mice for histological analyses.
Determination of the expression of CD31 and SM22α by immunofluorescent staining
Tissue slides were fixed using 4% paraformaldehyde for 10 min and incubated with the rabbit anti-CD31 and mouse anti-SM22α antibodies overnight at 4 °C. Subsequently, sections were incubated with Alexa Fluor 488-conjugated donkey anti-rabbit IgG (1:400) or Alexa Fluor 594-conjugated donkey anti-mouse IgG (1:400) antibody for 1 h. Nuclei were counterstained with DAPI for 5 min. Coverslips were sealed and visualized under a confocal microscope (Zeiss).
Determination of three layers and generation of teratomas by immunohistochemistry
Teratoma tissues were fixed in 4% paraformaldehyde at 4 °C overnight, and then embedded in paraffin and sectioned (7-μm thickness). Hematoxylin and eosin staining was performed on the sections. Tissue sections were also subjected to immunohistochemistry for Pan-CK, NeuN, vimentin, CD 31, and Ki67.
Identification of METTL1-regulated genes by RNA sequencing (RNA-seq) analysis
A total of 2 μg RNA per sample was used for RNA sequencing (at least two replicates per sample). A complementary DNA library was prepared, and sequencing was performed by Beijing Annoroad Gene Technology Co. Ltd. After adaptor trimming and low-quality sequence filtering, the reads were mapped to the human reference genome version. Bowtie2 v2.2.3 was used for building the genome index, and sequences were then aligned to the reference genome using HISAT2 v2.1.0. Read counts were determined by HTSeq v0.6.0, and FPKM (fragments per kilobase million mapped reads) was then calculated to estimate the expression level of genes in each sample. The raw counts were then converted to reads per kilobase per million mapped reads (RPKM) using edegR. Genes with q ≤ 0.05 and |log2_ratio|≥ 1 were identified as differentially expressed genes (DEGs).
Gene ontology (GO) and KEGG enrichment analysis
Functions of METTL1-regulated DEGs were analyzed by GO and KEGG in the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/summary.jsp)  . GO enrichment analysis can predict the functional roles of DEGs associated with METTL1 knockdown in hiPSCs. GO enrichment in DEGs was assessed using the hypergeometric test, which uses the genes in the whole genome as background. GO terms with q < 0.05 were considered to be significantly enriched. KEGG analysis can define the pathways related to the DEGs associated with METTL1 knockdown. It was performed using the hypergeometric test, in which P values were adjusted for multiple comparisons, providing q values. KEGG terms with q < 0.05 were considered to be significantly enriched.
Quantification and statistical analysis
All statistical analyses were performed using the GraphPad Prism 7.0 software (GraphPad Software, La Jolla, CA, USA). All results are representative of at least three independent experiments unless stated otherwise. The normality of the data was assessed using the Shapiro-Wilk test, while equality of group variance was assessed using the Brown-Forsythe test; these tests were performed using SigmaPlot 14 (Systat Software, San Jose, CA). Data are presented as the mean ± SEM of the biological replicates. Comparisons between two groups were performed by unpaired, two-tailed t test, while more than two groups were compared using one-way ANOVA followed by Bonferroni’s post hoc test or by two-way ANOVA. P values were two-tailed, and values < 0.05 were considered as statistically significant. *P < 0.05; **P < 0.01; ***P < 0.001.
Gene expression profiling in METTL1 knockdown hiPSCs
To study the function of the tRNA m7G methyltransferase METTL1 in hiPSCs, METTL1 expression was silenced in hiPSCs using lentivirus expressing METTL1-targeting shRNAs (METTL1-KD). Pooled clones were selected, and knockdown efficiency of METTL1 was confirmed by qRT-PCR and western blot (Fig. S1a).
RNA-seq analysis in hiPSCs led to the identification of 6426 DEGs, 2301 (6%) of which were upregulated, whilst 4125 (11%) were downregulated when METTL1 was silenced (Fig. 1a, b). Furthermore, KEGG pathway analysis showed that numerous pathways associated with stem cell regulation were significantly downregulated, including p53 signaling pathway, cell cycle, and the signaling pathways regulating stem cell pluripotency in METTL1-KD hiPSCs (Fig. 1e). Consistent with the pathway analysis, gene ontology (GO) enrichment analysis of the genes downregulated in METTL1-KD hiPSCs showed enrichment in regulation of transcription (DNA-templated) and cell division (Fig. 1f). Otherwise, in the upregulated pathways enriched by KEGG, we found that it is mainly related to immune responses, including systemic lupus erythematosus, graft-versus-host disease, and antigen processing presentation (Fig. S1c). qRT-PCR was performed to assess the expression of 5 representative genes differentially expressed following METTL1 silencing, which confirmed the reproducibility of our RNA-seq findings (Fig. 1c, d).
Additionally, a heatmap of downregulated genes revealed enrichment in neuroectoderm markers that are involved in neuronal cell development and differentiation (Fig. S1b). Consistent with the RNA-seq data, qRT-PCR verification showed that the expression of neuroectoderm-specific genes was downregulated in hiPSCs upon METTL1 silencing, whereas mesoderm- and endoderm-specific genes were expressed at higher levels (Fig. 1g). Of note, the downregulation of ISL1 is consistent with its expression pattern in endoderm and mesoderm/cardiogenic lineages as well as in neural stem cells [19, 20]. Taken together, these data suggested a role of METTL1 in maintaining hiPSC self-renewal and pluripotency.
METTL1 promotes proliferation required for hiPSC self-renewal by regulating cell cycle
We then assessed the relevance of METTL1-mediated m7G tRNA modifications in hiPSC proliferation and self- renewal using two different single-cell METTL1-KD hiPSC clones (Fig. 2a, b). We found that METTL1-KD cells proliferated slower compared with control hiPSCs (Fig. 2c). Moreover, cell cycle analysis revealed that the percentage of cells in the G2 phase was significantly higher in METTL1-KD hiPSCs compared with control cells (Fig. 2d, e). We also found that the mRNA and protein levels of critical cell cycle regulators were significantly altered in METTL1-KD hiPSCs (Fig. S2a, b).
Morphological analysis revealed that silencing of METTL1 resulted in flat colonies with increased surface area, as well as reduced the number of colonies (Fig. 2f, g). Notably, the cells found in the periphery of METTL1-KD hiPSC colonies appeared to be partially differentiated (blue arrows in Fig. 2f). A similar phenomenon was also observed in METTL1-KD hESCs (Fig. S2c), suggesting that METTL1 knockdown leads to a loss of self-renewal and promotes stem cell differentiation.
Since high levels of alkaline phosphatase (AP) can be used as a marker of undifferentiated state in stem cells , we sought to perform AP staining assays. We found a significant decrease in the AP-positive population in METTL1-KD hiPSCs (Fig. 2h, i). These data confirmed the critical role of METTL1 in maintaining stem cell self-renewal and preventing spontaneous differentiation.
METTL1-mediated m7G methylome is required for hiPSC pluripotency by regulating the translation of stem cell markers
Since RNA-seq analysis indicated that signaling pathways regulating stem cell pluripotency were downregulated in METTL1-KD hiPSCs, the expression of pluripotency marker genes following METTL1 knockdown was investigated. We found that METTL1 downregulation resulted in only minor changes in the mRNA levels of stem cell markers Oct4, Nanog, and Sox2. However, protein levels of these markers were drastically decreased in METTL1-KD cells (Fig. 3a, b). Similar results were observed in two different hESC lines (Fig. S3a, b), suggesting that METTL1-mediated m7G RNA modifications might regulate the translation of stem cell markers.
We next sought to assess the impact of METTL1 silencing on the translation efficiency of stem cell markers. To this end, control and METTL1-KD hiPSCs were subjected to polysome fractionation via sucrose gradient centrifugation (Fig. 3c) and the mRNA levels of stem cell markers in each of the fractions were measured by qRT-PCR. We found that the polysome-bound mRNA levels of these markers were significantly decreased in METTL1-KD hiPSCs (Fig. 3d), indicating that the translation efficiency for pluripotency genes, including OCT4, Nanog, and SOX2, was remarkably decreased when METTL1 was silenced.
We next performed an anti-m7G northwestern assay to assess whether m7G modifications were altered in METTL1-KD hiPSCs. We observed that the m7G signal was lost in the RNAs of METTL1-KD cells (Fig. 3e, f). We also assessed whether the exogenous expression of METTL1 in METTL1-KD cells could rescue the defects in RNA methylation (Fig. 3g). We found that transfer of METTL1-WT, but not a catalytic inactive mutant, into METTL1-KD hiPSCs rescued m7G modifications on RNAs (Fig. 3h, i). Moreover, METTL1-WT also rescued the expression of OCT4, Nanog, and SOX2 in METTL1-KD cells (Fig. 3j). These findings indicated that METTL1 deficiency impacted hiPSC pluripotency by impairing the translation of three key pluripotency genes.
METTL1 silencing accelerates hiPSC differentiation into embryoid bodies and promotes the expression of angiogenesis-related genes
HiPSC cells can spontaneously aggregate and form embryoid bodies (EBs) in vitro which comprised all three embryonic germ layers. Since our results suggested a potential role of METTL1 in hiPSC differentiation, we conducted an EB formation assay. To ensure the homogeneity of EBs, we used ultra-low attachment 96-well spheroid microplates to induce EB formation (Fig. 4a). HiPSCs formed EBs 6 days after seeding, and the EBs formed by METTL1-KD hiPSCs were significantly larger than those formed by control hiPSCs (Fig. 4b). We also assessed the expression of different lineage-specific markers and found that METTL1 knockdown resulted in increased expression of mesoderm- and endoderm-lineage markers, while the expression of ectoderm markers was reduced (Fig. 4c).
We also performed RNA-seq analysis in day 6 METTL1-KD-derived and control EBs. GO enrichment analysis of DEGs revealed enrichment of neuroectoderm markers, as well as genes involved in neuron differentiation and nervous system development among the downregulated in METTL1-KD EBs genes (Fig. 4d). Interestingly, a significant enrichment of genes involved in the regulation of angiogenesis and skeletal system development, which arises from mesoderm, was observed in the upregulated genes (Fig. 4e). Consistent with the RNA-seq findings, qRT-PCR confirmed the differential expression of mRNA expression of angiogenesis- and vasculogenesis-related genes (Fig. 4f). Moreover, we compared the change in gene expression levels for key pluripotent markers (Fig. S4a). We found that there were significant changes in multiple pluripotent factors regulating the stem cells, including SOX2/OCT4 and NANOG (Fig. S4b). According to literature reports, the upregulation of OCT4 and the downregulation of Nanog and SOX2 would force the stem cells to differentiate toward mesoderm [22,23,24,25,26], which coincide with our results. Gene set enrichment analysis (GSEA) indicated that mesoderm development and cardiovascular system development gene signatures were enriched in genes upregulated in METTL1-KD EBs (Fig. 4g). Furthermore, KEGG pathway analysis showed that numerous pathways associated with angiogenesis and vasculogenesis were significantly differentially regulated (Fig. 4h).
Overall, our data reveal that METTL1 silencing promotes embryonic differentiation in vitro and induces the expression of mesodermal markers, polarizing stem cells differentiation toward mesoderm lineages and promoting angiogenesis and vasculogenesis.
METTL1 knockdown enhances teratoma formation in vivo by promoting cell proliferation and angiogenesis in nude mice
We then sought to assess whether teratoma formation in vivo was affected by METTL1 silencing. To this end, we injected control and METTL1-KD hiPSCs into the posterior limbs of BALB/C nude mice (Fig. S5a). Teratoma formation was observed at 6 weeks post-injection; we observed that METTL1-KD-derived teratomas had a significantly increased volume and weight compared to control teratomas (Fig. 5a, b; Fig. S5b).
We also investigated if the enhanced teratoma formation was due to increased differentiation of hiPSCs, as observed in the EB formation assay. H&E staining revealed that both the METTL1-KD-derived and control teratomas were fully differentiated into the three germ layers (Fig. S5c). These results indicated that hiPSCs maintained their ability to differentiate into cell types derived from all three primary germ layers despite METTL1 silencing. However, immunohistochemical analyses revealed a decreased number of cells expressing the ectodermal marker pan-cytokeratin METTL1-KD-derived teratomas, while the number of cells expressing the mesodermal marker vimentin was increased (Fig. 5c, d). Ki67 expression was also increased in METTL1-KD-derived teratomas, suggesting a higher number of proliferating cells. Moreover, a lower number of cells expressing the neural differentiation marker NeuN were detected in METTL1-KD-derived teratomas. In contrast, METTL1-KD-derived teratomas contained increased levels of the vascular endothelial cell marker CD31. qRT-PCR confirmed the increased levels of CD31, vWF, VEGFR2, and SM22α expression in METTL1-KD hiPSC-derived teratomas compared to control hiPSC-derived teratomas (Fig. 5e, f).
METTL1 knockdown promotes angiogenesis in vivo
To further confirm the increased neovascularization observed in METTL1-KD-derived teratomas, the expression of the surface markers of neovascular endothelial cells and vascular smooth muscle cells (VSMCs) CD31 and SM22α was quantified in sections from METTL1-KD-derived and control teratomas; the number of cavities per high power field (hpf) and the relative tube area were also calculated (Fig. 6a). Consistent with our previous findings, we found that METTL1 silencing enhanced the formatting of new blood vessels by endothelial cells and VSMCs. Furthermore, the numbers of CD31-positive ECs and SMA-positive VSMCs were significantly increased in METTL1-KD-derived EBs compared to the control EBs (Fig. 6b). These results indicated that METTL1 suppresses angiogenesis in vivo, as the formation of a vascular network was profoundly enhanced in METTL1-KD-derived EBs.
The use of iPSCs can overcome the ethical issues of ESCs and is considered as a promising alternative for drug screening, disease modeling, and development of clinical treatments [27,28,29]. The molecular mechanisms underlying self-renewal and maintenance of pluripotency in iPSCs have been extensively studied. The stem cell genome undergoes reversible epigenetic modifications that regulate gene expression and dictate cell development and differentiation [8, 11, 30, 31]. In our previous study, we reported that METTL1-deficient mouse ESCs display defective self-renewal and neural differentiation . However, the exact molecular mechanisms by which METTL1 regulates pluripotency and differentiation of stem cells remained unknown.
RNA-seq analysis in METTL1-KD hiPSCs revealed an apparent downregulation of genes involved in maintaining stem cell self-renewal and pluripotency, including BMP4, FGF2, and WNT2. Simultaneously, various pathways associated with stem cell regulation were significantly downregulated [32, 33]. P53 signaling pathway has been considered to play an active role in promoting cell apoptosis and also plays a critical role in inhibiting the cell cycle of ESCs . The upregulated p53 signaling pathway will lead to ESC differentiation [35, 36]. In our study, we found that although p53 signaling pathway is significantly inhibited in METTL1-KD hiPSCs, the self-renewal and stability of METTL1-KD hiPSCs were impaired. This may be due to the extensive influence of METTL1 on other signaling pathways that regulate hiPSCs. For example, the METTL1 regulated FGF/TGF-β/Wnt signalings play critical roles in maintaining the stability of stem cells [37, 38].
Changes in cell cycle regulation are also closely associated with stem cell fate [39, 40]. The cell cycle in stem cells is characterized by a faster G1 phase and extended S phase [41, 42], and hiPSCs maintain pluripotency through rapid self-renewal and inhibition of cell differentiation. Our findings suggested that METTL1 silencing resulted in cell cycle arrest in hiPSCs, with a potential impact on their self-renewal, which was in accordance with the previous studies that have demonstrated METTL1’s involvement in cell cycle regulation and pluripotency in mESCs .
Recent studies have shown that METTL1-mediated m7G methylation is not restricted to tRNA but also occurs on mRNAs and miRNAs, increasing translation efficiency and miRNA processing, respectively [43, 44]. Our polysome profiling revealed that suppression of METTL1-mediated m7G tRNA modifications in hiPSCs resulted in a marked decrease in translation efficiency of the stem cell transcription factors OCT4, SOX2, and Nanog. We showed that METTL1 knockdown reduced m7G tRNA modification and decreased the mRNA translation of stem cell transcription factors. However, whether m7G modification on RNAs other than tRNAs contributes to this process requires further investigation.
Next, we verify the role of METTL1 in differentiation in hiPSCs; we performed EB formation assay in vitro. The results showed that METTL1-KD hiPSCs tended to differentiate into mesoderm including neovascularization while inhibiting neuroectoderm differentiation. The latest research of our team also showed that METTL1 can promote the ability and function of hiPSCs to differentiate into endothelial progenitor cells (EPC) in vitro . In addition, hiPSCs with pluripotency can differentiate into teratoma with three germ layers in mice [46, 47]. We found that teratomas derived from METTL1-KD hiPSCs exhibited lower levels of ectoderm markers and increased levels of mesoderm markers, consistent with our EB formation results and further supporting the role of METTL1 in the regulation of germ layer differentiation. Consistent with increased polarization toward mesoderm fate after METTL1 silencing, we observed increased vascular CD31+ blood vessels in METTL1-KD teratomas. The increased vascularization might explain the fact that the downregulation of METTL1 in hiPSCs enhanced in vivo tumor formation and progression. Angiogenesis is essential for embryonic development, as well as for tumor initiation and progression. Increased angiogenesis is often associated with enhanced cancer cell proliferation and tumor growth , which might explain why METTL1-KD teratomas exhibited increased cell proliferation and enhanced growth.
Taken together, the results of our study highlighted the role of METTL1 in the processes of pluripotency maintenance, self-renewal, and differentiation of hiPSCs. Our findings indicated, for the first time, the physiological function of METTL1 in hiPSCs and that METTL1-mediated m7G methylation is essential for translation of pluripotency transcription factors. Our findings also suggested that METTL1 suppresses angiogenesis by inhibiting mesoderm differentiation (Fig. 7). More importantly, since hiPSCs offer a powerful tool to model human ontogenetic processes in vivo and in vitro [46, 49], the findings of this study have strong implications in human embryo development research. However, the molecular mechanisms underlying the regulation of METTL1 during these processes require further investigation. Understanding how m7G on mRNAs and tRNAs are dynamically regulated in stem cell differentiation may provide improved strategies for the treatment of neurological and vascular diseases.
This study aimed at assessing the biological relevance and exact mechanisms by which METTL1 regulates pluripotency and differentiation in human induced pluripotent stem cells (hiPSCs). In summary, the results of this study make a significant contribution to the literature for our findings demonstrated that METTL1 is essential in maintaining pluripotency in human stem cells. We found METTL1 silencing resulted in alterations in the global m7G profile in hiPSCs, as well as reduced the translational efficiency of stem cell marker genes. METTL1-KD hiPSCs exhibited reduced pluripotency with slower cell cycling. METTL1 silencing also skewed hiPSC differentiation toward the mesoderm fate and suppressed neuroectoderm differentiation both in vivo and in vitro.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Fetal bovine serum
Differentially expressed genes
Embryonic stem cells
False discovery rate
Reads per kilobase of transcript per million mapped reads
Gene set enrichment analysis
Hematoxylin and eosin staining
Kyoto Encyclopedia of Genes and Genomes
Non-essential amino acid
Methylated RNA immunoprecipitation sequence
Quantitative real-time PCR
Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131(5):861–72.
Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, Nie J, Jonsdottir GA, Ruotti V, Stewart R, et al. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318(5858):1917–20.
Wu SM, Hochedlinger K. Harnessing the potential of induced pluripotent stem cells for regenerative medicine. Nat Cell Biol. 2011;13(5):497–505.
Stadtfeld M, Hochedlinger K. Induced pluripotency: history, mechanisms, and applications. Genes Dev. 2010;24(20):2239–63.
Thomson M, Liu SJ, Zou LN, Smith Z, Meissner A, Ramanathan S. Pluripotency factors in embryonic stem cells regulate differentiation into germ layers. Cell. 2011;145(6):875–89.
Kuppusamy KT, Sperber H, Ruohola-Baker H. MicroRNA regulation and role in stem cell maintenance, cardiac differentiation and hypertrophy. Curr Mol Med. 2013;13(5):757–64.
Warren L, Lin C. mRNA-based genetic reprogramming. Mol Ther. 2019;27(4):729–34.
Jambhekar A, Dhall A, Shi Y. Roles and regulation of histone methylation in animal development. Nat Rev Mol Cell Biol. 2019;20(10):625–41.
Srinageshwar B, Maiti P, Dunbar GL, et al. Role of epigenetics in stem cell proliferation and differentiation: implications for treating neurodegenerative diseases. Int J Mol Sci. 2016;17(2):199-214.
Ren R, Ocampo A, Liu GH, Izpisua Belmonte JC. Regulation of stem cell aging by metabolism and epigenetics. Cell Metab. 2017;26(3):460–74.
Wu Y, Zhou C, Yuan Q. Role of DNA and RNA N6-adenine methylation in regulating stem cell fate. Curr Stem Cell Res Ther. 2018;13(1):31–8.
Shaheen R, Abdel-Salam GM, Guy MP, Alomar R, Abdel-Hamid MS, Afifi HH, Ismail SI, Emam BA, Phizicky EM, Alkuraya FS. Mutation in WDR4 impairs tRNA m(7) G46 methylation and causes a distinct form of microcephalic primordial dwarfism. Genome Biol. 2015;16:210.
Lin S, Liu Q, Lelyveld VS, Choe J, Szostak JW, Gregory RI. Mettl1/Wdr4-mediated m(7) G tRNA methylome is required for normal mRNA translation and embryonic stem cell self-renewal and differentiation. Mol Cell. 2018;71(2):244–55 e245.
Ingolia NT, Lareau LF, Weissman JS. Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes. Cell. 2011;147(4):789–802.
Sampath P, Pritchard DK, Pabon L, Reinecke H, Schwartz SM, Morris DR, Murry CE. A hierarchical network controls protein translation during murine embryonic stem cell self-renewal and differentiation. Cell Stem Cell. 2008;2(5):448–60.
Subramanian A, Guo B, Marsden MD, Galic Z, Kitchen S, Kacena A, Brown HJ, Cheng G, Zack JA. Macrophage differentiation from embryoid bodies derived from human embryonic stem cells. J Stem Cells. 2009;4(1):29–45.
Itskovitz-Eldor J, Schuldiner M, Karsenti D, Eden A, Yanuka O, Amit M, Soreq H, Benvenisty N. Differentiation of human embryonic stem cells into embryoid bodies compromising the three embryonic germ layers. Mol Med. 2000;6(2):88–95.
da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Qu Q, Li D, Louis KR, Li X, Yang H, Sun Q, Crandall SR, Tsang S, Zhou J, Cox CL, et al. High-efficiency motor neuron differentiation from human pluripotent stem cells and the function of Islet-1. Nat Commun. 2014;5:3449.
Laugwitz KL, Moretti A, Lam J, Gruber P, Chen Y, Woodard S, Lin LZ, Cai CL, Lu MM, Reth M, et al. Postnatal isl1+ cardioblasts enter fully differentiated cardiomyocyte lineages. Nature. 2005;433(7026):647–53.
Lutolf MP, Gilbert PM, Blau HM. Designing materials to direct stem-cell fate. Nature. 2009;462(7272):433–41.
Ding J, Xu H, Faiola F, Ma'ayan A, Wang J. Oct4 links multiple epigenetic pathways to the pluripotency network. Cell Res. 2012;22(1):155–67.
Lengner CJ, Welstead GG, Jaenisch R. The pluripotency regulator Oct4: a role in somatic stem cells? Cell Cycle. 2008;7(6):725–8.
Loh YH, Wu Q, Chew JL, Vega VB, Zhang W, Chen X, Bourque G, George J, Leong B, Liu J, et al. The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells. Nat Genet. 2006;38(4):431–40.
Masui S, Nakatake Y, Toyooka Y, Shimosato D, Yagi R, Takahashi K, Okochi H, Okuda A, Matoba R, Sharov AA, et al. Pluripotency governed by Sox2 via regulation of Oct3/4 expression in mouse embryonic stem cells. Nat Cell Biol. 2007;9(6):625–35.
Simandi Z, Horvath A, Wright LC, Cuaranta-Monroy I, De Luca I, Karolyi K, Sauer S, Deleuze JF, Gudas LJ, Cowley SM, et al. OCT4 acts as an integrator of pluripotency and signal-induced differentiation. Mol Cell. 2016;63(4):647–61.
Salaris F, Rosa A. Construction of 3D in vitro models by bioprinting human pluripotent stem cells: challenges and opportunities. Brain Res. 2019;1723:146393.
Lee CT, Bendriem RM, Wu WW, Shen RF. 3D brain organoids derived from pluripotent stem cells: promising experimental models for brain development and neurodegenerative disorders. J Biomed Sci. 2017;24(1):59.
Flaherty EK, Brennand KJ. Using hiPSCs to model neuropsychiatric copy number variations (CNVs) has potential to reveal underlying disease mechanisms. Brain Res. 2017;1655:283–93.
Ankam S, Rovini A, Baheti S, Hrstka R, Wu Y, Schmidt K, Wang H, Madigan N, Koenig LS, Stelzig K, et al. DNA methylation patterns in human iPSC-derived sensory neuronal differentiation. Epigenetics. 2019;14(9):927–37.
Dovey OM, Foster CT, Cowley SM. Histone deacetylase 1 (HDAC1), but not HDAC2, controls embryonic stem cell differentiation. Proc Natl Acad Sci U S A. 2010;107(18):8242–7.
Wu M, Chen G, Li YP. TGF-beta and BMP signaling in osteoblast, skeletal development, and bone formation, homeostasis and disease. Bone Res. 2016;4:16009.
Mohammadnia A, Yaqubi M, Pourasgari F, Neely E, Fallahi H, Massumi M. Signaling and gene regulatory networks governing definitive endoderm derivation from pluripotent stem cells. J Cell Physiol. 2016;231(9):1994–2006.
Hong H, Takahashi K, Ichisaka T, Aoi T, Kanagawa O, Nakagawa M, Okita K, Yamanaka S. Suppression of induced pluripotent stem cell generation by the p53-p21 pathway. Nature. 2009;460(7259):1132–5.
Lee DF, Su J, Ang YS, Carvajal-Vergara X, Mulero-Navarro S, Pereira CF, Gingold J, Wang HL, Zhao R, Sevilla A, et al. Regulation of embryonic and induced pluripotency by aurora kinase-p53 signaling. Cell Stem Cell. 2012;11(2):179–94.
Han MK, Song EK, Guo Y, Ou X, Mantel C, Broxmeyer HE. SIRT1 regulates apoptosis and Nanog expression in mouse embryonic stem cells by controlling p53 subcellular localization. Cell Stem Cell. 2008;2(3):241–51.
Van Camp JK, Beckers S, Zegers D, Van Hul W. Wnt signaling and the control of human stem cell fate. Stem Cell Rev Rep. 2014;10(2):207–29.
Coutu DL, Galipeau J. Roles of FGF signaling in stem cell self-renewal, senescence and aging. Aging (Albany NY). 2011;3(10):920–33.
Calder A, Roth-Albin I, Bhatia S, Pilquil C, Lee JH, Bhatia M, Levadoux-Martin M, McNicol J, Russell J, Collins T, et al. Lengthened G1 phase indicates differentiation status in human embryonic stem cells. Stem Cells Dev. 2013;22(2):279–95.
Sela Y, Molotski N, Golan S, Itskovitz-Eldor J, Soen Y. Human embryonic stem cells exhibit increased propensity to differentiate during the G1 phase prior to phosphorylation of retinoblastoma protein. Stem Cells. 2012;30(6):1097–108.
Pauklin S, Vallier L. The cell-cycle state of stem cells determines cell fate propensity. Cell. 2013;155(1):135–47.
Dalton S. Linking the cell cycle to cell fate decisions. Trends Cell Biol. 2015;25(10):592–600.
Zhou M, Wang H, Zeng X, Yin P, Zhu J, Chen W, Li X, Wang L, Wang L, Liu Y, et al. Mortality, morbidity, and risk factors in China and its provinces, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2019;394(10204):1145–58.
Pandolfini L, Barbieri I, Bannister AJ, Hendrick A, Andrews B, Webster N, Murat P, Mach P, Brandi R, Robson SC, et al. METTL1 promotes let-7 MicroRNA processing via m7G methylation. Mol Cell. 2019;74(6):1278–90 e1279.
Deng Y, Zhou Z, Lin S, Yu B. METTL1 limits differentiation and functioning of EPCs derived from human-induced pluripotent stem cells through a MAPK/ERK pathway. Biochem Biophys Res Commun. 2020;527(3):791-8.
Nelakanti RV, Kooreman NG, Wu JC: Teratoma formation: a tool for monitoring pluripotency in stem cell research. Curr Protoc Stem Cell Biol 2015, 32:4A 8 1–17.
Kooreman NG, Wu JC. Tumorigenicity of pluripotent stem cells: biological insights from molecular imaging. J R Soc Interface. 2010;7(Suppl 6):S753–63.
Carmeliet P, Jain RK. Molecular mechanisms and clinical applications of angiogenesis. Nature. 2011;473(7347):298–307.
Bulic-Jakus F, Katusic Bojanac A, Juric-Lekic G, Vlahovic M, Sincic N. Teratoma: from spontaneous tumors to the pluripotency/malignancy assay. Wiley Interdiscip Rev Dev Biol. 2016;5(2):186–209.
We thank the members of our laboratory for providing technical advice and encouragement.
This work was supported by the National Key Research and Development Program of China (2016YFC1300600), the National Natural Science Foundation of China (U1601219 and 8197042012), Guangdong Natural Science Foundation Outstanding Youth Scholar Project (2019B151502011), and Guangzhou People’s Livelihood Science and Technology Project (201903010006).
Ethics approval and consent to participate
All animal procedures were performed in accordance with the protocols approved by the Institutional Animal Care and Use Committee of Sun Yat-Sen University.
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.
Identification of METTL1-modulated genes by RNA-seq. a qRT-PCR and WB analyses of METTL1 expression levels in KD, WT, and control cells. b Heatmap showing differentially expressed genes involved in neurectoderm development. Each lane corresponds to an independent biological sample. Upregulated and downregulated genes are shown in red and blue, respectively. c KEGG pathway analysis of genes upregulated in METTL1-KD hiPSCs.
METTL1 regulates cell cycle progression in hiPSCs and morphology of ESCs. a qRT-PCR analysis assessing the mRNA levels of cell cycle markers in METTL1-KD and control hiPSCs. b Western blot assessing the protein levels of cell cycle regulators. Images are representative of three independent experiments. c Colony-formation assay using METTL1-KD and control hESCs. Representative images and colony quantification data are shown. Data are presented as mean ± SEM; n = 3; *, P < 0.05; **, P < 0.01; ***, P < 0.001.
METTL1-mediated m7G methylation is required for OCT4, Nanog, and SOX2 expression in ESCs. (a- b) qRT-PCR and WB assessing the mRNA and protein levels of Oct4, Nanog, and Sox2 in METTL1-KD and control H1ESCs (a) and H9ESCs (b) cells. Data are presented as mean ± SEM; n = 3; ***, P < 0.001.
Expression of pluripotency genes on METTL1-KD and control hiPSCs derived EBs. (a) Heatmap showing DEGs involved in stem cells pluripotency pathways. Each lane corresponds to an independent biological sample. Scale bar: log2 FPKM. (b) qRT-PCR analysis assessing the mRNA levels of OCT4, Nanog, and SOX2 in control and METTL1-KD derived EBs. Data are presented as mean ± SEM; n = 3; **, P < 0.01; ***, P < 0.001.
METTL1 knockdown promotes teratoma differentiation in nude mice. a Schematic diagram summarizing the procedure of subcutaneous injection of matrigel solution containing METTL1-KD hiPSCs or control hiPSCs into 5-week-old BALB/C nude mice. b Tumors were isolated 6 weeks after injection. c Haematoxylin and eosin (HE) staining of the teratoma tissues. The development of three germ layers (ectoderm, mesoderm, and endoderm) was evident in the teratoma tissue.
List of culture mediums and chemical reagents used in this study. Table S2. List of primary and secondary antibodies used. Table S3. Cell lines used in this study. Table S4. Sequences of the primes used in this study. Table S5. Components of polysome cell extraction buffer.
About this article
Cite this article
Deng, Y., Zhou, Z., Ji, W. et al. METTL1-mediated m7G methylation maintains pluripotency in human stem cells and limits mesoderm differentiation and vascular development. Stem Cell Res Ther 11, 306 (2020). https://doi.org/10.1186/s13287-020-01814-4
- N7-methylguanosine (m7G)
- Human induced pluripotent stem cells (hiPSCs)