Skip to main content
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care
  • Subscribe
  • Log in
  • My Cart
  • Follow ada on Twitter
  • RSS
  • Visit ada on Facebook
Diabetes

Advanced Search

Main menu

  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care

User menu

  • Subscribe
  • Log in
  • My Cart

Search

  • Advanced search
Diabetes
  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
Islet Studies

m6A mRNA Methylation Controls Functional Maturation in Neonatal Murine β-Cells

  1. Yanqiu Wang1,
  2. Jiajun Sun1,
  3. Zhen Lin2,
  4. Weizhen Zhang3,
  5. Shu Wang1,
  6. Weiqing Wang1⇑,
  7. Qidi Wang1,4⇑ and
  8. Guang Ning1⇑
  1. 1Shanghai National Clinical Research Center for Endocrine and Metabolic Diseases, Key Laboratory for Endocrine and Metabolic Diseases of the National Health Commission of the PR China, Shanghai Institute of Endocrine and Metabolic Diseases, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
  2. 2State Key Laboratory of Molecular Biology, Shanghai Key Laboratory of Molecular Andrology, CAS Center for Excellence in Molecular Cell Science, Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences, University of Chinese Academy of Sciences, Shanghai, China
  3. 3Department of Physiology and Pathophysiology, School of Basic Science, Peking University Health Science Center, Beijing, China
  4. 4Sino-French Research Center for Life Sciences and Genomics, Ruijin Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China
  1. Corresponding author: Guang Ning, gning{at}sibs.ac.cn, Qidi Wang, wqd11094{at}rjh.com.cn, or Weiqing Wang, wqingw{at}hotmail.com
  1. Y.W. and J.S. contributed equally to this work.

Diabetes 2020 Aug; 69(8): 1708-1722. https://doi.org/10.2337/db19-0906
PreviousNext
  • Article
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF
Loading

Abstract

The N6-methyladenosine (m6A) RNA modification is essential during embryonic development of various organs. However, its role in embryonic and early postnatal islet development remains unknown. Mice in which RNA methyltransferase-like 3/14 (Mettl3/14) were deleted in Ngn3+ endocrine progenitors (Mettl3/14nKO) developed hyperglycemia and hypoinsulinemia at 2 weeks after birth. We found that Mettl3/14 specifically regulated both functional maturation and mass expansion of neonatal β-cells before weaning. Transcriptome and m6A methylome analyses provided m6A-dependent mechanisms in regulating cell identity, insulin secretion, and proliferation in neonatal β-cells. Importantly, we found that Mettl3/14 were dispensable for β-cell differentiation but directly regulated essential transcription factor MafA expression at least partially via modulating its mRNA stability. Failure to maintain this modification impacted the ability to fulfill β-cell functional maturity. In both diabetic db/db mice and patients with type 2 diabetes (T2D), decreased Mettl3/14 expression in β-cells was observed, suggesting its possible role in T2D. Our study unraveled the essential role of Mettl3/14 in neonatal β-cell development and functional maturation, both of which determined functional β-cell mass and glycemic control in adulthood.

Introduction

Mature β-cells have the unique ability to secrete insulin in response to extracellular glucose concentrations (1). Defects in either β-cell number and/or function can lead to loss of functional β-cell mass and, eventually, to diabetes (2). Establishment of adequate functional β-cell mass is essential for glycemic control in adulthood, which is mainly accomplished during neonatal period before weaning by β-cell differentiation, self-replication, and functional maturation (3,4). Neonatal β-cells are immature and are unable to secrete insulin appropriately in response to a glucose challenge (5). During the early postnatal period (P0-P14), neonatal β-cells undergo dramatic changes to acquire glucose-responsive insulin secretion ability and gain functional maturity (3,6,7). Previous studies have identified orchestrated mechanisms including molecular pathways (8), transcriptional signals (9–12), and epigenetic regulators (5,13) in shaping the transcriptional networks that reinforce the functional identity of mature β-cells. Epigenetic regulators, such as DNA methylation and histone modifications, have profound effects on β-cell development and functional maturation (5,13). We and others have reported that DNMT3A directs functional maturation in murine β-cells (5,8). Meanwhile, genome-wide analysis of histone marks during β-cell development further underscored the importance of histone modifications in β-cell maturation (13). N6-methyladenosine (m6A) modification of mRNA is emerging as an important regulator of gene expression (14,15). However, its role in neonatal β-cells is still unknown.

Being the most prevalent internal modification in mRNA, the m6A mRNA modification functionally modulates mRNA splicing, export, localization, translation, and stability and has crucial roles in various normal and pathological processes (16,17). m6A methylation is dynamic and is accomplished by the orchestrated action of methyltransferase complex methyltransferase-like 3 (Mettl3), methyltransferase-like 14 (Mettl14), and Wilms tumor 1-associated protein (WTAP), while it is removed by demethylases obesity-associated protein (FTO) and AlkB homolog 5 (ALKBH5) (18). In particular, reversible m6A mRNA methylation plays critical roles in modulating dynamic transcriptome switching in response to diverse signals during embryonic development of various organs and organisms (17). As the key writers of m6A, Mettl3 and/or Mettl14 have been identified to regulate Drosophila neuronal function and sex determination (19), mouse stem cell pluripotency (20), naive T cell differentiation (21), cortical neurogenesis (22), spermatogenesis (23,24), and hematopoietic stem/progenitor cell specification (25). Mettl3 and or Mettl14 were also reported to participate in the pathology and progression of multiple human cancer diseases, i.e., myeloid leukemia (26), liver cancer (27), lung cancer (28), and breast cancer (29). Very recently, m6A modification was shown to regulate mature β-cell insulin secretion and survival (30,31). It is currently unknown whether and how m6A RNA modifications regulate embryonic and early postnatal islet development.

To investigate the role of Mettl3/14-mediated m6A modifications in islet development, we generated mice in which Mettl3 and/or Mettl14 were deleted in Ngn3+ endocrine progenitors. Loss of Mettl3/14 in endocrine progenitors (Mettl3/14nKO) caused hyperglycemia at 2 weeks of age, with reduced cell number and impaired functional maturation of neonatal β-cells. Transcriptome and m6A methylome analyses of primary islets from P14 wild-type (WT) and Mettl3/14nKO mice identified m6A-dependent mechanisms in regulating neonatal β-cell identity and insulin secretion. Loss of Mettl3/14 silenced MafA expression, possibly via modulating its mRNA stability and, thus prevented neonatal β-cell to fulfill functional maturation. Importantly, as early as P0, reduced MafA protein abundance and impaired glucose-stimulated insulin secretion (GSIS) were already present in mutant islets, which were prior to β-cell loss and changes in other transcription factors, i.e., Pdx1 and Nkx6.1. The findings that Mettl3 and Mettl14 expression were significantly decreased in pancreatic β-cells of db/db mice and patients with type 2 diabetes (T2D) further suggested the possible role of m6A in T2D in both rodents and humans.

Research Design and Methods

Mice

Mettl3- and Mettl14-floxed mice (Mettl3flox/flox and Mettl14flox/flox) were generated and provided by Prof. Ming-Han Tong (University of Chinese Academy of Sciences) (24). Ngn3-cre mice were a kind gift from Prof. Wei-Zhen Zhang. The Ngn3-cre;Mettl3flox/flox (Mettl3nKO) or Ngn3-cre;Mettl14flox/flox (Mettl14nKO) mice were generated by crossing Mettl3flox/flox or Mettl14flox/flox mice with Ngn3-Cre mice. Mettl3- and Mettl14- double knockout (KO) mice (Mettl3/14nKO) were generated by crossing Mettl3flox/flox Mettl14flox/flox mice with Ngn3-cre;Mettl3flox/+ Mettl14flox/+ mice. Ngn3-Cre, Mettl3flox/flox, Mettl14flox/flox, or Mettl3flox/flox Mettl14flox/flox mice were used as their littermate controls. All mice were housed in the animal facility on a 12-h/12-h light/dark cycle. Normal chow and water were available ad libitum. Blood glucose concentrations were measured by glucometers, and plasma insulin concentrations were determined by using an ELISA kit (Mouse Ultrasensitive Insulin ELISA kit; Alpco). Intraperitoneal glucose tolerance tests were performed on 8-week-old mice after overnight fasting as previously described (8). Insulin tolerance tests were performed on 4-week-old Mettl3/14nKO mice and WT mice after 6 h of fasting. Male mice were used in all the experiments in the current study, unless otherwise stated. All animal experiments were approved by the Animal Care Committee of Shanghai Jiao Tong University.

Human Subjects

Paraffin sections of pancreas far from the margin of the pancreatectomy were collected from our previous research (32). In brief, all patients with partial pancreatectomy performed in Ruijin Hospital between 2013 and 2017 were enrolled. Those who had been reported as having a malignant tumor were excluded, and then five patients with T2D and five age- and BMI-matched subjects who did not have diabetes (nondiabetic [ND]) were finally included in this study. Detailed information and clinical characteristics for each patient were listed in Supplementary Table 1. This study was approved by the Institutional Review Board of the Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine and was in accordance with the principles of the Declaration of Helsinki II.

Cell Culture and Lentivirus Infection

The mouse insulinoma cells (MIN6 cell line) purchased from CAMS Cell Culture Center (Beijing, China) were grown in DMEM medium (Gibco) containing 25.0 mmol/L glucose, 15% FBS, 100 IU/mL penicillin, 100 μg/mL streptomycin, 10.2 mmol/L L-glutamine, and 2.5 mmol/L β-mercaptoethanol at 37°C in a humidified 5% CO2 atmosphere. For knocking down Mettl3 and Mettl14, shRNA lentiviruses that targeted Mettl3, Mettl14, or a control lentivirus were constructed, packaged, purified, and titrated at GeneChem Co. Ltd. MIN6 cells were infected with purified lentivirus at 50 multiplicity of infection for 48 h. After infection, cells were harvested for RNA extraction and further analysis.

Immunostaining Analysis

Pancreas sections were dissected, fixed, and processed as described before (33). For the immunostaining analysis, entire pancreatic tissues were continuously sectioned at 5-μm thickness. Immunochemistry procedures were performed on continuous sections (selected every 150-μm apart, 10–12 sections per animal) to obtain representative β-cell mass information of the whole pancreas. Immunochemistry staining of insulin for β-cell mass analysis was performed using a diaminobenzidine peroxidase substrate kit (Vector Laboratories, Burlingame, CA) counterstaining with eosin. Digital images of whole pancreas were captured by MZ 100 microscope (Nikon Corp., Tokyo, Japan). Total pancreatic and insulin-positive areas of each section were measured using Meta-Morph version 7.1 (Molecular Devices, Sunnyvale, CA).

Immunoblot Analysis

P0 mouse islets and MIN6 cells were lysed, quantified, and blotted as described before. Primary antibodies are listed as following: rabbit anti-MAFA (1:1,000; Bethyl); rabbit anti-METTL3 (1:1,000; Proteintech); rabbit anti-METTL14 (1:1,000; Proteintech); rabbit anti-Pdx1 (1:1,000; CST); mouse anti-Nkx6.1 (1:1,000; DSHB); and mouse anti-GAPDH (1:10,000; Proteintech) was used as an internal control to normalize band intensity.

Extraction of RNA and Quantitative Real-Time PCR Analysis

Total cell RNA was extracted using the TRIzol reagent (Invitrogen). Reverse transcription and quantitative real-time PCR were performed as previously described (8). PCRs were performed in duplicate. The expression levels were normalized to individual β-actin. Primers used in this study were listed in Supplementary Table 2.

m6A Quantification

Islets were isolated from 8-week-old Mettl3nKO, Mettl14nKO, Mettl3/14nKO, and WT mice, and then, total RNA was extracted using a RNeasy Micro kit (Qiagen) following the manufacturer’s protocol. The change of global m6A levels in mRNA was measured using EpiQuik m6A RNA Methylation Quantification Kit (Colorimetric) (Epigentek) following the manufacturer’s protocol; 200-ng poly-A-purified RNA was used for each sample analysis.

mRNA Stability Assay

The mRNA half-life measurements were performed according to previous publications (34). In brief, MIN6 cells with or without Mettl3/14 knockdown were treated with 5-μg/mL actinomycin D (Sigma) for 1 or 3 h at the end of culture and then collected for RNA extraction and real-time PCR analysis. Since actinomycin D treatment results in transcription stalling, the change of mRNA concentration at a given time (dC/dt) is proportional to the constant of mRNA decay (Kdecay) and mRNA concentration (C).

m6A MeRIP Sequencing, RNA Sequencing, and Data Analysis

Islet total RNA was prepared as previously described. The m6A mRNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were performed by Cloudseq Biotech Inc. (Shanghai, China) according to published procedures. Briefly, m6A RNA immunoprecipitation was performed with the GenSeq m6A-MeRIP Kit (GenSeq Inc., Shanghai, China) following the manufacturer’s instructions. Both the input sample without immunoprecipitation and the m6A IP samples were used for library generation with NEBNext Ultra II Directional RNA Library Prep Kit (New England Biolabs, Inc.). The library quality was evaluated with BioAnalyzer 2100 system (Agilent Technologies, Inc.). Library sequencing was performed on an illumina Hiseq instrument with 150 base pair (bp) paired-end reads. Paired-end reads were harvested from the Illumina HiSEq 4000 sequencer, and they were quality controlled by Q30 after 3′ adaptor trimming and removal of low-quality reads by cutadapt software (v1.9.3). First, clean reads of all libraries were aligned to the reference genome (MM10) by Hisat2 software (v2.0.4). Methylated sites on RNAs (peaks) were identified by MACS software. Differentially methylated sites were identified by diffReps. These peaks identified by both software overlapping with exons of mRNA were figured out and chosen by homemade scripts. For RNA-seq, total RNA was used for removing the rRNAs with NEBNext rRNA Depletion Kit (New England Biolabs, Inc.) following the manufacturer’s instructions. RNA libraries were constructed, and libraries were controlled for quality and quantified. Library sequencing was performed on an illumina Hiseq instrument with 150 bp paired end reads. Then, guided by the Ensembl gtf gene annotation file, cuffdiff software (part of cufflinks) was used to get the gene-level FPKM (fragments per kilobase of exon model per million mapped reads) as the expression profiles of mRNA. Fold change and P value were calculated based on FPKM.

Statistics Analyses

The exact sample size for each experiment was indicated in the figure legends. All statistics comparing two groups used two-sided Student t tests. ANOVA was used for multiple groups. Statistical analyses were performed with GraphPad Prism 7. P < 0.05 was considered as statistically significant.

Data and Resource Availability

The P0 and P14 RNA-seq and m6A MeRIP-seq data are deposited at the Gene Expression Omnibus (GEO) database under the accession numbers: GSE149193, GSE132323, and GSE132319.

Results

Loss of m6A in Ngn3+ Endocrine Progenitors Leads to Hypoinsulinemia and Severe Diabetes

To examine whether Mettl3/14-mediated m6A mRNA methylation participates in β-cell development, we first checked the expression patterns of Mettl3 and Mettl14 in mouse pancreas at E17.5, P0, P4, P8, P14, and P56 by immunostaining (Fig. 1A). We found relatively weak Mettl3 and Mettl14 expressions in β-cells at E17.5 and P0, whereas the immunofluorescence intensities of both the two methyltransferases increased dramatically during β-cell maturation window P4-P8, reached maximal levels at P14, and maintained at high levels throughout adulthood (Fig. 1B and C). These data suggested that Mettl3- and Mettl14-mediated m6A mRNA methylation might participate in postnatal β-cell functional maturation.

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Loss of m6A in Ngn3+ endocrine progenitors leads to hypoinsulinemia and severe diabetes. A: Representative pancreatic sections from E17.5-P14 WT mice were immunostained for Mettl3 or Mettl14 (red) with insulin (green). B and C: Relative immune fluorescence intensity of Mettl3 (B) and Mettl14 (C) inside islets were determined at indicated times (n = 3). D: Model illustration of the generation of Mettl3nKO, Mettl14nKO, and Mettl3/14nKO mice. E: Representative pancreatic sections from P14 Mettl3/14nKO and WT mice were immunostained for Mettl3 or Mettl14 (red) with Ins (white) and Gcg (green). F: Relative m6A amounts relative to adenosine (A) in mRNA extracted from the Mettl3nKO, Mettl14nKO, Mettl3/14nKO, and WT islets were quantified (n = 3). G: Random blood glucose levels of Mettl3nKO, Mettl14nKO, Mettl3/14nKO, and WT mice were monitored weekly (n = 4–9). H and I: Intraperitoneal glucose tolerance tests were performed on 8-week-old WT, Mettl3nKO (H), and Mettl14nKO (I) mice after overnight fasting (n = 5). J: Plasma insulin levels of P14 Mettl3/14nKO and WT mice were determined (n = 3–4). K: Body weights of Mettl3nKO, Mettl14nKO, Mettl3/14nKO, and WT mice were monitored weekly (n = 4–9). Data are presented as mean ± SEM of independent experiment indicated as above. *P < 0.05, **P < 0.01, ***P < 0.001 by Student t test. Nuclei were counterstained with DAPI (blue). Scale bars, 20 μm.

We then generated mice in which Mettl3 (Ngn3-cre;Mettl3flox/flox, named as Mettl3nKO), Mettl14 (Ngn3-cre;Mettl14flox/flox, named as Mettl14nKO) or both Mettl3/Mettl14 (Ngn3-cre;Mettl3flox/flox/Mettl14flox/flox, named as Mettl3/14nKO) were specifically ablated in Ngn3+ endocrine progenitors (Fig. 1D). The Ngn3-Cre, Mettl3flox/flox, Mettl14flox/flox, and Mettl3flox/flox/Mettl14flox/flox mice were all healthy and phenotypically normal, and thus, they were jointly used as WT littermate controls in the following experiments. Successful KO of Mettl3 and Mettl14 were confirmed by immunostaining: both METTL3 and METTL14 proteins were selectively absent in endocrine cells of P14 Mettl3/14nKO mice (Fig. 1E). Moreover, m6A/A% levels in RNA extracted from Mettl3nKO and Mettl14nKO islets were significantly reduced compared with WT controls, and this decrease was more pronounced in Mettl3/14nKO islets (Fig. 1F).

At 2 weeks after birth, Mettl3nKO and Mettl14nKO had comparable random blood glucose levels as WT controls (Fig. 1G). However, Mettl3nKO and Mettl14nKO began to show significant increases in random blood glucose levels at 4 weeks (14.34 ± 3.26 vs. 9.95 ± 1.18 mmol/L in WT, P < 0.01) and 3 weeks (13.22 ± 1.28 vs. 9.77 ± 1.70 mmol/L in WT, P < 0.01) of age, respectively; and their blood glucose rose gradually with age (Fig. 1G). At 8 weeks of age, both Mettl3nKO and Mettl14nKO showed a dramatic increase in glycemia after intraperitoneal glucose injection (Fig. 1H and I). Notably, a more severe phenotype was observed in double KO mice: at 2 weeks of age, Mettl3/14nKO mice already showed a slight increase in random blood glucose levels (9.00 ± 1.45 vs. 7.88 ± 1.22 mmol/L in WT, P = 0.07), and they developed severe hyperglycemia as early as 3 weeks of age (random blood glucose levels >20 mmol/L) and had persistently severe hyperglycemia (>30 mmol/L) thereafter (Fig. 1G). Moreover, plasma insulin levels in 2-week-old Mettl3/14nKO mice was reduced by 70%, suggesting β-cell failure predominated in these diabetic mutants (Fig. 1J). None of the mutants showed differences in body weight, compared with their age-matched controls (Fig. 1K). Our data also showed that insulin sensitivity was not impaired in 4-week-old male Mettl3/14nKO mice compared with WT (Supplementary Fig. 1A). Female Mettl3/14nKO mice were also diabetic at 2 weeks of age and manifested hyperglycemia thereafter (Supplementary Fig. 1B).

Mettl3/14 Preferentially Regulates Postnatal β-Cell Mass Establishment

We performed immunostaining against four endocrine hormones including insulin (β), glucagon (α), somatostatin (δ), and pancreatic polypeptide (PP) on P0 pancreas from WT and Mettl3/14nKO. We found similar cord-like islet structures in newborn Mettl3/14nKO and WT mice (Fig. 2A). Then, we calculated the number of four different islet endocrine cells per section and found the absolute number of β, α, δ, and PP cells remained constant at birth between the two groups (Fig. 2B). Interestingly, after 14 days, Mettl3/14nKO islets exhibited significantly reduced β-cell number per section (369.6 ± 45.1 vs. 819.4 ± 217.2 in WT, P < 0.05) (Fig. 2C and D). On the contrary, no differences were detected in the number of α, δ and PP cells per section between the two groups (Fig. 2C and D). The above data indicated a preferential role of Mettl3/14 in β-cell development during early postnatal period (P0-P14). In parallel, β-cell mass was 30% lower in Mettl3/14nKO mice than in WT, while α-cell mass remained unchanged (Fig. 2E and F). The reduction in β-cell mass was attributed to changes in both proliferation and apoptosis: 40% reduction in the proportion of Ki67+insulin+ cells (Fig. 2G and Supplementary Fig. 2A) and a twofold increase in the percentage of TUNEL+insulin+ cells (Fig. 2H and Supplementary Fig. 2B) were found in 2-week-old Mettl3/14nKO. β-Cell size was unchanged in mutant islets (Fig. 2I and Supplementary Fig. 2C). Taken together, Mettl3/Mettl14 specifically regulated β-cell number via controlling neonatal β-cell proliferation and survival and, thus, determined postnatal β-cell mass formation during P0-P14.

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2

Loss of m6A preferentially regulates postnatal β-cell mass establishment. A: Representative pancreatic sections immunostained for Gcg (red), Sst (red), PP (pancreatic polypeptide [Ppy], red) with Ins (green) in Mettl3/14nKO and WT pancreas at P0. B: The absolute cell number of α, β, δ, and PP cells per pancreatic section of P0 Mettl3/14nKO and WT mice were determined (n = 3). C: Representative pancreatic sections immunostained for Gcg (red), Sst (red), PP (Ppy, red) with Ins (green) in Mettl3/14nKO and WT pancreas at P14. D: The absolute cell number of α, β, δ, and PP cells per pancreatic section at P14 Mettl3/14nKO and WT were determined (n = 3). E–I: β-Cell mass (E), α-cell mass (F), the percentage of Ki67+/insulin+ β-cells (G), the percentage of TUNEL+/insulin+ β-cells (H), and β-cell size (I) in P14 WT and Mettl3/14nKO mice were determined (n = 3). Data are presented as mean ± SEM of independent experiment indicated as above, *P < 0.05 by Student t test. Nuclei were counterstained with DAPI (blue). Scale bars, 20 μm.

Loss of Cell Identity and Function in Neonatal Mettl3/14nKO β-Cells

Neonatal β-cells are immature and undergo a dramatic change in the early postnatal period to become functional mature and acquire glucose responsiveness (3). We then checked expression levels of the critical transcription factors required for β-cell identity and function in WT and Mettl3/14nKO at P14. We found that Pdx1, MafA, and Nkx6.1 were absent in a vast majority of Mettl3/14-deficient β-cells at P14 (Fig. 3A–C). We then calculated the percentages of Pdx1+insulin+, MafA+insulin+, and Nkx6.1+insulin+ β-cells in P14 Mettl3/14nKO and WT pancreas and found that they decreased about 63.3%, 49.0%, and 78.8%, respectively, in P14 mutants (Fig. 3D–F). In parallel, the expression of Urocortin 3 (Ucn3) (3), a molecular marker of mature β-cells, was present in a vast majority of insulin+ β-cells in WT, but it was remarkably decreased in P14 mutant β-cells (Fig. 3G). We also detected reduced expression of GLUT2 (35) in P14 Mettl3/14nKO (Fig. 3H). Then, we isolated primary islets from P14 WT and Mettl3/14nKO to evaluate their function in vitro. P14 Mettl3/14nKO islets displayed impaired insulin release at both basal (2.8 mmol/L) and high glucose (16.7 mmol/L) levels (Fig. 3I). We further studied β-cell ultrastructure in P14 Mettl3/14nKO and WT pancreas using transmission electron microscopy (Fig. 3J). On the basis of the recognized criterion (36), we quantified the secretory vesicle subsets including mature (with dense core granules, red arrow), immature (with light core granules, green arrow), and empty vesicles (with no granules, blue arrow) (Fig. 3J). The proportion of mature vesicles was significantly decreased and the ratios of immature and empty vesicles were significantly increased in Mettl3/14-deficient β-cells (Fig. 3K). These results indicated that Mettl3/14 deficiency led to the loss of β-cell identity, impaired glucose-induced insulin release, and defections in secretory vesicle maturation in neonatal pancreatic β-cells. Moreover, we detected substantial Aldh1a3+ cells in 8-week-old Mettl3/14nKO islets (Fig. 3L). Meanwhile, we also performed immunostaining against chromogranin A (CGA) and four endocrine cell types (insulin [ins], glucagon [Gcg], somatostatin [Sst], and PP) in 8-week-old Mettl3/14nKO and WT pancreas (Fig. 3M). We observed a dramatic increase in the percentage of dedifferentiated cells (islet hormone-negative and CGA-positive cells: 30.8 ± 9.0% vs. 3.1 ± 1.5%, Mettl3/14nKO vs. WT, P < 0.001, yellow arrows in Fig. 3M) in 8-week-old Mettl3/14nKO pancreas sections. These results indicated that β-cells underwent dedifferentiation in 8-week-old Mettl3/14nKO mice.

Figure 3
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3

Loss of β-cell identity and function in P14 Mettl3/14nko mice. A–C: Representative pancreatic sections from P14 WT and Mettl3/14nKO were coimmunostained for Ins (green) and β-cell identity genes: Pdx1(red) (A); MafA (red) (B) and Nkx6.1 (red) (C). D–F: The percentages of Pdx1+, MafA+, and Nkx6.1+ β-cells in P14 WT and Mettl3/14nKO were determined (n = 3). G: Representative pancreatic sections from P14 WT and Mettl3/14nKO were coimmunostained for Ins (green) and β-cell maturation marker Ucn3 (red). H: Representative pancreatic sections from P14 WT and Mettl3/14nKO mice were coimmunostained for Ins (green) and β-cell GLUT2 (red). I: Isolated P14 islets were incubated at 2.8 mmol/L or 16.7 mmol/L glucose for 1 h. Secreted insulin levels were measured and were normalized to total insulin in the islets (n = 3). J and K: Representative transmission electron microscopy of islets from P14 WT and Mettl3/14nKO mice. Red and green arrows point to typical vesicles containing mature and immature granules. Blue arrows point to typical empty vesicles (J). Analysis of the percentage of mature, immature, and empty vesicles from P14 WT and Mettl3/14nKO mice (K) (n = 3). L: Representative immunofluorescence stainings against Aldh1a3 (red) and insulin (green) in 8-week-old WT and Mettl3/14nKO pancreas sections are shown. M: Pancreatic sections of 8-week-old WT and Mettl3/14nKO immunostained with insulin (green), endocrine cocktail (Gcg/Sst/PP, red), and CGA (white) are shown. Yellow arrows indicated the dedifferentiated cells (endocrine hormone–negative/CGA-positive). N: The ratio of dedifferentiated (Dediff) cells was determined (n = 3). Data are presented as mean ± SEM of independent experiment indicated as above. **P < 0.01, ***P < 0.001 by Student t test. Nuclei were counterstained with DAPI (blue). Scale bars, 20 μm or 500 nm.

Transcriptome Profile of P14 Mettl3/14nKO Islets

To investigate the underlying mechanisms, we performed RNA-seq analysis on P14 WT and Mettl3/14nKO islets. We detected significant differences in the expression levels of 544 genes (fold change >1.5, P < 0.05), among which 294 were upregulated and 250 were downregulated (Fig. 4A). We further analyzed alterations of known pathways using the ingenuity pathway analysis and listed the top downregulated and upregulated signaling pathways in P14 Mettl3/14nKO islets (Fig. 4B). Pathways related to insulin secretion, ion transportation, and response to glucose were among the top downregulated cellular biological processes, while the Wnt-signaling pathway, negative regulation of intracellular signal transduction, and negative regulation of cell proliferation pathway were significantly upregulated in mutant islets (Fig. 4B and C). Interestingly, several crucial β-cell transcription factors Pdx1, Nkx6.1, and MafA were downregulated at P14, among which MafA was the most significantly repressed gene after loss of Mettl3/14 (Fig. 4C). Moreover, a cluster of genes involved in ion transportation and insulin secretion were reduced in mutant islets (i.e., Fkbp1b, Ero1lb, Sytl4, Maob, Nnat, and Trpm5) (Fig. 4C). Loss of Mettl3/14 also affected expressions of genes important for β-cell proliferation, including downregulation of several positive proliferation regulators (i.e., Nkx6.1, Rtkn2, Nasp, and Plk3) and upregulation of some proliferation inhibitors (i.e., Sfrp5, lims2, and Klf11) (Fig. 4C). The dramatic reduction in Pdx1, Nkx6.1, MafA, and Trpm5 mRNA levels were still observed in 8-week-old Mettl3/14nKO islets (Supplementary Fig. 3A). The above results reinforced that methyltransferase-like Mettl3/14 indeed regulated the expression of those critical β-cell genes.

Figure 4
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4

Transcriptome and m6A methylome profile of P14 Mettl3/14nKO islets. A: Scatter plots of relative expression levels of 544 differentially expressed genes identified from RNA-seq of P14 WT and Mettl3/14nKO islets. Upregulated genes are marked in red, and downregulated genes are marked in blue (fold change >1.5, P < 0.05). B: GO analysis of differentially expressed genes as identified by RNA-seq of P14 WT and Mettl3/14nKO associated with β-cell function. TOR, target of rapamycin. C: The heatmap shows relative expression levels of genes critical for β-cell identity, insulin secretion, ion transportation, and proliferation inhibition between P14 WT and Mettl3/14nKO. D: Sequence motif identified within m6A peaks in P14 islets by HOMER database. E: Pie chart depicting the fraction of m6A peaks in five transcript segments in P14 WT islets. TSS, transcriptional start site. F: GO analysis of enriched signaling pathways of m6A target genes in P14 islets. G: Venn diagram showed the overlap between genes with m6A modifications and genes that were differentially expressed in RNA-seq of P14 Mettl3/14nkO. H: Top five downregulated and upregulated genes in the intersection of m6A-modified genes and differentially expressed genes in RNA-seq are presented in the table, respectively. I: Integrative Genomics Viewer tracks showed RNA-seq reads distribution in MafA mRNA of P14 WT and Mettl3/14nKO (upper panel) and showed MeRIP-seq reads distribution in MafA mRNA (lower panel). J: MIN6 cells were transfected with Lv-ShMettl3, Lv-ShMettl14, Lv-ShMettl3/14, or control virus for 48 h; then, MafA protein abundance was assayed by immunoblot. K: The MafA mRNA remaining after actinomycin D (ActD) treatment with or without ShMettl3/14 in MIN6 cells is shown. (n = 3) Data are presented as mean ± SEM of independent experiments indicated as above. *P < 0.05. **P < 0.01 by one-way ANOVA. CDS, coding sequences.

To identify the direct effects of Mettl3/14, we treated MIN6 cells with Mettl3 and Mettl14 shRNA lentiviruses. Successful knockdown of Mettl3 and Mettl14 significantly reduced mRNA expression of β-cell transcription factors Pdx1, Nkx6.1, and MafA (Supplementary Fig. 4A), increased expression of proliferation inhibitors (Sfrp5, Tlr2, Frzb, Lims2, and Klf11) (Supplementary Fig. 4B), diminished expression of genes related to ion transportation and insulin secretion (Ero1lb, Nnat, Pcsk9, Trpm5, Fkbp1b, Arrb1, Syt3, Flna, Sytl4, and Kcng3) (Supplementary Fig. 4C). Taken together, both in vivo and in vitro, we showed that loss of Mettl3/14 impaired genes important for β-cell identity, insulin secretion, and proliferation.

Identification of m6A Targets in Neonatal Mettl3/14nKO Islets

To obtain insights into the transcriptome-wide m6A distribution in neonatal β-cells, we performed MeRIP-seq analysis using P14 mouse islets. Compared with the unbound fractions, 8,987 methylation sites (corresponding to 4,873 transcripts; fold change >2) were significantly enriched in m6A antibody-bound fractions and, thus, were identified as high-confidence m6A targets. m6A peaks in islet transcripts were significantly enriched in GGACU motif (Fig. 4D). m6A was distributed throughout islet mRNA transcripts and were abundant predominantly in the coding sequences (Fig. 4E). Gene Ontology (GO) analysis showed m6A-containing mRNAs were enriched in biological processes related to RNA Polymerase II Transcription, cell projection organization, RNA Metabolism, protein ubiquitination, mRNA metabolic process, and cell morphogenesis involved in differentiation (Fig. 4F).

To identify the direct m6A-modified targets, we correlated MeRIP-seq results with RNA-seq data identified from P14 islets. The correlation analysis eventually identified 170 differentially expressed genes with m6A modifications (Fig. 4G). Several m6A-modified genes which were important for β-cell function, such as MafA (11), Maob (37), Fkbp1b (38), and Kcng3 (39) were identified as direct m6A targets, and their expressions were significantly downregulated in neonatal Mettl3/14nKO islets (Fig. 4H and Supplementary Fig. 5A–C). On the contrary, Sfrp5 (40), a well-established proliferation inhibitor in β-cells, was upregulated in mutants (Supplementary Fig. 5D).

Mettl3/14-Mediated m6A Modification Regulates MafA Expression via Modulating mRNA Stability

We listed the top five downregulated and upregulated m6A-modified genes that were preferentially changed in neonatal Mettl3/14nKO mice (Fig. 4H), among which MafA was identified as the most significantly reduced gene (Fig. 4H). Integrative Genomics Viewer tracks showed MafA transcripts had abundantly enriched m6A peaks spread its mRNA, which indeed had the highest methylation score among all enriched transcripts (Fig. 4H and I). On the contrary, the other two preferentially changed transcription factors, Pdx1 and Nkx6.1, showed no m6A modifications in murine β-cells as identified in our MeRIP-seq.

To determine the direct effect of Mettl3/14 on MafA protein expression, we then treated MIN6 cells with ShMettl3/14 virus for 48 h and found reduced MafA protein abundance after loss of Mettl3 and/or Mettl14 in MIN6 cells (Fig. 4J). These results supported that Mettl3/14 directly regulated MafA expression in β-cells.

It is known that m6A modifications on mRNA transcripts might affect mRNA stability and translation (17). We then treated control or Mettl3/14-deficient MIN6 cells with transcription inhibitor actinomycin D and checked the role of m6A on MafA mRNA decay. We found that knockdown of Mettl3/14 significantly decreased the mRNA stability of MafA (Fig. 4K), but had no effect on Pdx1, Nkx6.1 mRNA stability (Supplementary Fig. 6A and B). These findings indicated that Mettl3/14 directly regulates MafA expression in neonatal β-cells, at least partly, via modulating its mRNA stability.

Transcriptome Profile and Functional Analysis on P0 Mettl3/14nKO Islets

To exclude the possibility that Mettl3/Mettl14 inactivation compromised β-cell differentiation that contributed to later β-cell functional immaturity, comprehensive gene expression analyses on P0 islets were performed. At birth, mutant mice showed comparable blood glucose levels (5.37 ± 0.98 vs. 4.23 ± 0.92 mmol/L in WT, P = 0.17) and islet composition (Fig. 2A and B) compared with WT. RNA-seq on isolated P0 islets identified 945 differentially expressed genes (fold change >1.5, P < 0.05), among which 191 were upregulated and 754 were downregulated (Fig. 5A). Pathway analysis identified that pathways related to the response to external stimulus, developmental process, metabolic process, and regulation of secretion were downregulated, while pathways like response to stress, defense response, and the apoptosis process were significantly upregulated (Fig. 5B). We further compared these genes with the MeRIP-seq database and identified 326 differentially expressed genes with m6A modifications (Fig. 5C and D). Among them, genes related to ion transportation and insulin secretion (i.e., Fxyd5, Slc3a2, Vamp8, and Syt9) were significantly downregulated (Fig. 5D).

Figure 5
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5

Mettl3/14 regulates MafA expression to drive β-cell functional maturation. A: Scatter plots of relative expression levels of 945 differentially expressed genes identified from RNA-seq of P0 WT and Mettl3/14nKO islets (fold change >1.5, P < 0.05). B: GO analysis of differentially expressed genes as identified by RNA-seq of P0 WT and Mettl3/14nKO. C: Venn diagram shows the overlap between genes with m6A modifications and genes that were differentially expressed in RNA-seq of P0 Mettl3/14nkO. D: Heatmap shows relative expression levels of these 326 overlapped genes between P0 WT and Mettl3/14nKO. E: Representative pancreatic sections from P0 WT and Mettl3/14nKO mice were coimmunostained for Ins (green) and β-cell identity genes: MafA/Pdx1/Nkx6.1 (red), Ins (green), and GLUT2 (red). F and G: P0 WT and Mettl3/14nKO islets were isolated, and then Mettl3, Mettl14, and MafA (F) and Pdx1 and Nkx6.1 (G) protein abundances were assayed by immunoblot. H: Band intensities of MafA, Pdx1, and Nkx6.1 proteins were normalized based on the corresponding GAPDH intensity (n = 3). I: P0 WT and Mettl3/14nKO islets were isolated and incubated with 2.8/16.7 mmol/L glucose for 1 h. Secreted insulin was normalized to total insulin in the islets (n = 4–5). Data are presented as mean ± SEM, **P < 0.01 by Student t test. Nuclei were counterstained with DAPI (blue). Scale bars, 20 μm.

Importantly, although similar expression levels of MafA, Nkx6.1, and Pdx1 mRNA were observed in RNA-seq, by in situ immunostaining, we still detected significantly reduced MafA protein abundance in insulin-positive cells of P0 Mettl3/14nKO islets, whereas no reductions in protein levels were detected in Pdx1 and Nkx6.1 at this early time point (Fig. 5E). We also found slightly diminished GLUT2 expression in P0 Mettl3/14nKO β-cells (Fig. 5E). However, Ucn3 was rarely expressed in new-born β-cells (Supplementary Fig. 7A). Moreover, we performed Western blot on isolated islets from P0 WT and mutant mice. Again, we found reduced MafA protein abundance in P0 Mettl3/14-deficient islets (Fig. 6F and H), whereas the expression levels of Pdx1 and Nkx6.1 were comparable to that of WT islets (Fig. 6G and H). The above data indicated that m6A-induced MafA change was prior to that of Pdx1 and Nkx6.1 at P0. Indeed, we found impaired glucose responsive insulin secretion ability in P0 Mettl3/14nKO islets. The mutant islet displayed comparable basal insulin secretion at 2.8 mmol/L glucose, but they had a trend toward reduced GSIS response under high glucose (16.7 mmol/L) stimulation (P = 0.059) (Fig. 5I).

Figure 6
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6

Reduced Mettl3 and Mettl14 expression in db/db mice and patients with T2D. A: Representative pancreatic sections from db/db and lean control mice were coimmunostained for Ins (white), Gcg (green), Mettl3 (red), or Mettl14 (red). B and C: The percentages of Mettl3 or Mettl14-positive/insulin-positive cells in db/db and lean control mice were calculated (n = 3). D: Representative images of pancreatic sections from patients with T2D and those who were ND stained for Ins (white), Gcg (green), Mettl3 (red), or Mettl14 (red) are shown. BMI values are given in kg/m2. E and F: The ratio of Mettl3- or Mettl14-positive/insulin-positive cells in patients with T2D and in patients who were ND were determined (n = 5). Data are presented as mean ± SEM. **P < 0.01 by Student t test. Nuclei were counterstained with DAPI (blue). Scale bars, 20 μm.

Reduced Mettl3 and Mettl14 Expression in db/db Mice and Patients With T2D

To evaluate a possible role of Mettl3/14-mediated m6A RNA methylation in the pathology process of T2D, we checked the expression of these two methyltransferases in islets of lean control and db/db mice at 8 weeks of age (Fig. 6A). The random blood glucose level of db/db mice was significantly higher than lean controls (23.2 ± 3.5 vs. 8.2 ± 0.6 mmol/L, db/db vs. WT, P < 0.05). We found Mettl3 and Mettl14 were abundantly expressed in both endocrine and exocrine cells in pancreas of lean control mice; however, by counting the percentage of Mettl3+ or Mettl14+ β-cells in diabetic db/db mice, both of their expressions in β-cells were significantly reduced by 31.3% and 27.1%, respectively, compared with lean control mice (Fig. 6A–C).

We also evaluated the expression of METTL3 or METTL14 in β-cells in pancreas sections of five patients with T2D and five age- and BMI-matched ND patients who had performed partial pancreatectomy (Supplementary Table 1). According to American Diabetes Association guidelines, all the five patients with T2D were under adequate glucose control with antidiabetic therapy (Supplementary Table 1). We found METTL3 and METTL14 were present in islet endocrine cells in the ND group, including both α-cells and β-cells. In contrast, islets from patients with T2D showed dramatic reductions in METTL3 and METTL14 expression (Fig. 6D). We further calculated the percentage of METTL3+ or METTL14+ β-cells in the two groups: in diabetic islets, the ratio of METTL3+/insulin+ cells (50.1 ± 16.2% vs. 95.6 ± 0.6%, T2D group vs. ND group, P < 0.01) (Fig. 6E) and METTL14+/insulin+ cells (55.3 ± 9.9% vs. 94.2 ± 3.0%, T2D group vs. ND group, P < 0.01) (Fig. 6F) declined substantially. These observations indicated that loss of Mettl3/14 might be involved in pathogenesis of T2D in both rodents and human.

Discussion

Recent studies have reported the role of posttranscriptional mRNA modifications on organ development, cellular function, and human diseases, thus highlighting its involvement in both physiological and pathological status (19,23–25,41). Dynamic m6A modification was known to affect development and fate determination of various cell types (i.e., hematopoietic stem/progenitor cells or spermatogenesis) (23–25). Germ cell–specific deletion of the m6A RNA methyltransferase Mettl3 or Mettl14 disrupted spermiogenesis through inhibiting diploid spermatogonia stem cells proliferation/differentiation (24). Neural stem cells lacking Mettl14 displayed markedly decreased proliferation and premature differentiation (42). During zebrafish embryogenesis, Mettl3-depletion led to continuous activation of the Notch signaling pathway, ultimately resulting in the blockage of endothelial-to-hematopoietic transition to specify the earliest hematopoietic stem/progenitor cells (25). Very recently, Men et al. (31) reported that Mettl14 deficiency in β-cells resulted in defects of β-cell survival and insulin secretion (30). However, whether m6A modification was implicated in neonatal islet development and the detailed regulation mechanisms remained unknown.

In the current study, by using transgenic mice in which Mettl3/14 were specifically ablated in Ngn3+ endocrine progenitors, we observed that the composition of the four endocrine cell types (α, β, δ, and PP cells) remained unchanged in WT and Mettl3/14nKO at birth. Strikingly, a dramatic drop in β-cell number and mass during the weaning period (P0-P14) led to the conclusion that m6A is essential for neonatal β-cell mass establishment. Interestingly, this regulation seems to be specific to β-cells, since Mettl3/14 deficiency had no effect on α, δ, and PP cell formation and expansion. Since Ngn3+ cells would also give rise to enteroendocrine cells, the interaction between enteroendocrine cells and endocrine cells should also be considered.

The postnatal weaning period is pivotal for immature β-cells to acquire their glucose-responsive insulin-producing mature phenotype (3). Previous studies have provided evidence on transcriptional factors (i.e., Nkx6.1, MafA, and Pdx1) (9–11) and cellular signals (i.e., mechanistic target of rapamycin [mTOR], connective tissue growth factor [CTGF], and estrogen-related receptor γ [ERRγ]) (8,43,44) in modulating β-cell functional maturation and/or postnatal expansion. In recent years, epigenetic modifiers were confirmed to be involved in the functional maturation of β-cells in the postnatal period to promote adequate metabolic reprogramming (5,8,13). For the first time, we demonstrated that Mettl3 and Mettl14 were essential for functional maturation of neonatal murine β-cells. First, the loss of Mettl3/14 compromised expression of critical β-cell transcription factors Pdx1, MafA, and Nkx6.1, which were reported to modulate β-cell maturation and function. Second, Mettl3/14-deficient β-cells exhibited reduced expression of mature β-cell marker Ucn3 and GLUT2. Third, impaired glucose-induced insulin release was detected in Mettl3/14nKO β-cells, with defects in genes responsible for glucose sensing (loss of GLUT2) (35), insulin biosynthesis (loss of Ero1lb) (45), and granule docking (loss of Sytl4) (46). Fourth, the secretory vesicles in Mettl3/14nKO β-cells were mostly immature. Interestingly, the pattern of Mettl3 and Mettl14 expression during the physiological maturation window in murine β-cells further supported our observation. The relatively low levels of Mettl3/14 at E17.5 and P0 might explain that they are dispensable for cell lineage formation during differentiation. The strong induction of Mettl3/14 during early postnatal period P0-P14 and high expression level throughout adulthood reinforced the importance of Mettl3/14 in driving β-cell functional maturation and maintaining β-cell function.

Combined analysis of P14 RNA-seq and MeRIP-seq identified that MafA was the direct target of Mettl3/14 for identity maintenance and functional maturation of neonatal β-cells. In Mettl3/14nKO pancreas, remarkably diminished MafA expression was found as early as P0, before changes of other transcriptional factors (Pdx1 and Nkx6.1) and islet morphology. Moreover, our MeRIP-seq analysis of neonatal islets demonstrated abundant m6A modifications in MafA mRNA, which were mostly located in the exon region, whereas Pdx1 and Nkx6.1 mRNA had rather weak or no m6A modifications. Mechanical studies proved that knockdown of Mettl3/14 reduced MafA mRNA stability and shortened its mRNA half-life, while had no effect on that of Pdx1 and Nkx6.1. These observations suggest m6A/MafA serves as a direct trigger of functional maturation process in neonatal β-cells.

To exclude the possibility that Mettl3/Mettl14 inactivation compromised β-cell differentiation that contributed to later β-cell functional immaturity, comprehensive gene expression analyses in P0 β-cells were performed. In our P0 RNA-seq, we did not detect significant differences in MafA mRNA, which might be caused by their relatively low expression level at this early time point. However, we confirmed MafA downregulation at the protein level in P0 mutant islets by both immunofluorescence staining and Western blot. The downregulation in MafA expression was accompanied with reduced GLUT2 expression and impaired GSIS in P0 mutant islets, when the islet composition and blood glucose level were unaffected. Interestingly, m6A-induced MafA alteration at P0 preceded and was followed by the changes in Nkx6.1 and Pdx1. It is well accepted that Pdx1 and Nkx6.1 are essential for β-cell differentiation (47–50). On the contrary, MafA is not involved in endocrine specification during embryonic development but is critical for β-cell functional maturation (51,52). After birth, MafA expression remains low, and it increases gradually with age (53). The role of MafA on regulating replication/survival and function of β-cells is gradually increased after birth, but it is not conspicuous at birth (54). Our findings provide evidence that Mettl3/14 regulates MafA expression and, thus, controls β-cell functional maturation after birth.

Indeed, our Mettl3/14nKO mice were to some extent phenotypically similar to MafAKO mice (11): both of them showed normal islet morphology and comparable β-cell numbers at birth but a progressive loss of β-cell identity and function with age. The link between m6A/MafA modulation in neonatal β-cells indicates that this finely tuned regulation network is crucial for proper β-cell maturation during a critical time window. Interestingly, like other diabetic mice models, we observed a dramatic increase in the percentage of dedifferentiated β-cells in 8-week-old Mettl3/14nKO islets. Considering the high percentage of dedifferentiation (30%) and the relatively short hyperglycemic period (8 weeks old), the possible role of Mettl3/14 on β-cell dedifferentiation, independent of hyperglycemia, needs to be further clarified. Moreover, it would be more ideal to perform RNA-seq analyses on purified β-cells from earlier developmental times instead of P14 islets to fully clarify the identical regulatory mechanisms in Mettl3/14 KO β-cells.

Using diabetic mouse models, we and others have shown that loss of Mettl3/14 in β-cells is one of the critical mechanisms in β-cell failure during the pathology of diabetes (30,31). The current study provided observational evidence on remarkably reduced METTL3 and METTL14 expression in islets from patients with T2D compared with patients who were ND, revealing the significance of Mettl3/14-mediated m6A signaling in the occurrence and development of human diabetes. This is consistent with a very recent finding by De Jesus et al. (55), who reported that m6A mRNA methylation was involved in human β-cell biology. A deeper understanding of the role of m6A mRNA methylation in β-cell biology at different stages of development or disease would help to find therapeutic strategies for maintaining or regenerating functional β-cell mass. Taken together, our data showed, for the first time, that as the most abundant mRNA modification in mammalian cells, Mettl3/14-mediated m6A modification played critical roles in governing neonatal β-cell functional maturation and identity maintenance. Loss of Mettl3/14 in neonatal islets led to overt neonatal diabetes due to incapability of establishing adequate functional β-cell mass after birth. Our results highlight the important role of Mettl3/14-mediated m6A modification in driving functional maturation in neonatal β-cells, which might be a potential target for diabetes therapy.

Article Information

Acknowledgments. The authors thank Prof. Ming-Han Tong (University of Chinese Academy of Sciences, Shanghai, China) for providing the Mettl3flox/flox and Mettl14flox/flox mice and for his valuable discussion and comments for the study. The authors thank Prof. Wei-Zhen Zhang (Department of Physiology and Pathophysiology, School of Basic Science, Peking University Health Science Center, Beijing, China) for providing the Ngn3-cre mice.

Funding. This work was supported by National Natural Sciences Foundation of China grants 81670700 and 81870527.

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. Y.W. and J.S. performed all the experiments and analyzed the data. Y.W., J.S., and Q.W. wrote the manuscript. Z.L., W.Z., and S.W. contributed to the data discussion. W.W., Q.W., and G.N. designed the project, supervised research, and coordinated the execution of the experimental plan. G.N. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Footnotes

  • This article contains supplementary material online at https://doi.org/10.2337/figshare.12272627.

  • Received September 10, 2019.
  • Accepted May 8, 2020.
  • © 2020 by the American Diabetes Association
https://www.diabetesjournals.org/content/license

Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at https://www.diabetesjournals.org/content/license.

References

  1. ↵
    1. Prentki M,
    2. Nolan CJ
    . Islet beta cell failure in type 2 diabetes. J Clin Invest 2006;116:1802–1812
    OpenUrlCrossRefPubMedWeb of Science
  2. ↵
    1. Meier JJ,
    2. Bonadonna RC
    . Role of reduced β-cell mass versus impaired β-cell function in the pathogenesis of type 2 diabetes. Diabetes Care 2013;36(Suppl. 2):S113–S119
    OpenUrlFREE Full Text
  3. ↵
    1. Blum B,
    2. Hrvatin S,
    3. Schuetz C,
    4. Bonal C,
    5. Rezania A,
    6. Melton DA
    . Functional beta-cell maturation is marked by an increased glucose threshold and by expression of urocortin 3. Nat Biotechnol 2012;30:261–264
    OpenUrlCrossRefPubMed
  4. ↵
    1. Wang Y,
    2. Sun J,
    3. Ni Q, et al
    . Dual effect of Raptor on neonatal β-cell proliferation and identity maintenance. Diabetes 2019;68:1950–1964
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Dhawan S,
    2. Tschen SI,
    3. Zeng C, et al
    . DNA methylation directs functional maturation of pancreatic β cells. J Clin Invest 2015;125:2851–2860
    OpenUrlCrossRefPubMed
  6. ↵
    1. Rorsman P,
    2. Arkhammar P,
    3. Bokvist K, et al
    . Failure of glucose to elicit a normal secretory response in fetal pancreatic beta cells results from glucose insensitivity of the ATP-regulated K+ channels. Proc Natl Acad Sci U S A 1989;86:4505–4509
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Otonkoski T,
    2. Andersson S,
    3. Knip M,
    4. Simell O
    . Maturation of insulin response to glucose during human fetal and neonatal development. Studies with perifusion of pancreatic isletlike cell clusters. Diabetes 1988;37:286–291
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Ni Q,
    2. Gu Y,
    3. Xie Y, et al
    . Raptor regulates functional maturation of murine beta cells. Nat Commun 2017;8:15755
    OpenUrlCrossRefPubMed
  9. ↵
    1. Gao T,
    2. McKenna B,
    3. Li C, et al
    . Pdx1 maintains β cell identity and function by repressing an α cell program. Cell Metab 2014;19:259–271
    OpenUrlCrossRefPubMed
    1. Taylor BL,
    2. Liu FF,
    3. Sander M
    . Nkx6.1 is essential for maintaining the functional state of pancreatic beta cells. Cell Rep 2013;4:1262–1275
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    1. Nishimura W,
    2. Takahashi S,
    3. Yasuda K
    . MafA is critical for maintenance of the mature beta cell phenotype in mice. Diabetologia 2015;58:566–574
    OpenUrlCrossRefPubMed
  11. ↵
    1. Gu C,
    2. Stein GH,
    3. Pan N, et al
    . Pancreatic beta cells require NeuroD to achieve and maintain functional maturity. Cell Metab 2010;11:298–310
    OpenUrlCrossRefPubMedWeb of Science
  12. ↵
    1. Lu TT,
    2. Heyne S,
    3. Dror E, et al
    . The polycomb-dependent epigenome controls β cell dysfunction, dedifferentiation, and diabetes. Cell Metab 2018;27:1294–1308.e7
    OpenUrlCrossRef
  13. ↵
    1. Holoch D,
    2. Moazed D
    . RNA-mediated epigenetic regulation of gene expression. Nat Rev Genet 2015;16:71–84
    OpenUrlCrossRefPubMed
  14. ↵
    1. Meyer KD,
    2. Jaffrey SR
    . The dynamic epitranscriptome: N6-methyladenosine and gene expression control. Nat Rev Mol Cell Biol 2014;15:313–326
    OpenUrlCrossRefPubMed
  15. ↵
    1. Fu Y,
    2. Dominissini D,
    3. Rechavi G,
    4. He C
    . Gene expression regulation mediated through reversible m⁶A RNA methylation. Nat Rev Genet 2014;15:293–306
    OpenUrlCrossRefPubMed
  16. ↵
    1. Roundtree IA,
    2. Evans ME,
    3. Pan T,
    4. He C
    . Dynamic RNA modifications in gene expression regulation. Cell 2017;169:1187–1200
    OpenUrlCrossRefPubMed
  17. ↵
    1. Lee M,
    2. Kim B,
    3. Kim VN
    . Emerging roles of RNA modification: m(6)A and U-tail. Cell 2014;158:980–987
    OpenUrlCrossRefPubMed
  18. ↵
    1. Lence T,
    2. Akhtar J,
    3. Bayer M, et al
    . m6A modulates neuronal functions and sex determination in Drosophila. Nature 2016;540:242–247
    OpenUrlCrossRefPubMed
  19. ↵
    1. Geula S,
    2. Moshitch-Moshkovitz S,
    3. Dominissini D, et al
    . Stem cells. m6A mRNA methylation facilitates resolution of naïve pluripotency toward differentiation. Science 2015;347:1002–1006
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Li HB,
    2. Tong J,
    3. Zhu S, et al
    . m6A mRNA methylation controls T cell homeostasis by targeting the IL-7/STAT5/SOCS pathways. Nature 2017;548:338–342
    OpenUrlCrossRefPubMed
  21. ↵
    1. Yoon KJ,
    2. Ringeling FR,
    3. Vissers C, et al
    . Temporal control of mammalian cortical neurogenesis by m6A methylation. Cell 2017;171:877–889.e17
    OpenUrlCrossRefPubMed
  22. ↵
    1. Xu K,
    2. Yang Y,
    3. Feng GH, et al
    . Mettl3-mediated m6A regulates spermatogonial differentiation and meiosis initiation. Cell Res 2017;27:1100–1114
    OpenUrlCrossRefPubMed
  23. ↵
    1. Lin Z,
    2. Hsu PJ,
    3. Xing X, et al
    . Mettl3-/Mettl14-mediated mRNA N6-methyladenosine modulates murine spermatogenesis. Cell Res 2017;27:1216–1230
    OpenUrlCrossRefPubMed
  24. ↵
    1. Zhang C,
    2. Chen Y,
    3. Sun B, et al
    . m6A modulates haematopoietic stem and progenitor cell specification. Nature 2017;549:273–276
    OpenUrlCrossRefPubMed
  25. ↵
    1. Barbieri I,
    2. Tzelepis K,
    3. Pandolfini L, et al
    . Promoter-bound METTL3 maintains myeloid leukaemia by m6A-dependent translation control. Nature 2017;552:126–131
    OpenUrlCrossRefPubMed
  26. ↵
    1. Chen M,
    2. Wei L,
    3. Law CT, et al
    . RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology 2018;67:2254–2270
    OpenUrlCrossRefPubMed
  27. ↵
    1. Lin S,
    2. Choe J,
    3. Du P,
    4. Triboulet R,
    5. Gregory RI
    . The m(6)A methyltransferase METTL3 promotes translation in human cancer cells. Mol Cell 2016;62:335–345
    OpenUrlCrossRefPubMed
  28. ↵
    1. Cai X,
    2. Wang X,
    3. Cao C, et al
    . HBXIP-elevated methyltransferase METTL3 promotes the progression of breast cancer via inhibiting tumor suppressor let-7g. Cancer Lett 2018;415:11–19
    OpenUrl
  29. ↵
    1. Liu J,
    2. Luo G,
    3. Sun J, et al
    . METTL14 is essential for β-cell survival and insulin secretion. Biochim Biophys Acta Mol Basis Dis 2019;1865:2138–2148
    OpenUrl
  30. ↵
    1. Men L,
    2. Sun J,
    3. Luo G,
    4. Ren D
    . Acute deletion of METTL14 in β-cells of adult mice results in glucose intolerance. Endocrinology 2019;160:2388–2394
    OpenUrl
  31. ↵
    1. Sun J,
    2. Ni Q,
    3. Xie J, et al
    . β-cell dedifferentiation in patients with T2D with adequate glucose control and nondiabetic chronic pancreatitis. J Clin Endocrinol Metab 2019;104:83–94
    OpenUrl
  32. ↵
    1. Li W,
    2. Zhang H,
    3. Nie A, et al
    . mTORC1 pathway mediates beta cell compensatory proliferation in 60 % partial-pancreatectomy mice. Endocrine 2016;53:117–128
    OpenUrl
  33. ↵
    1. Chen CY,
    2. Ezzeddine N,
    3. Shyu AB
    . Messenger RNA half-life measurements in mammalian cells. Methods Enzymol 2008;448:335–357
    OpenUrlCrossRefPubMedWeb of Science
  34. ↵
    1. Rorsman P,
    2. Renström E
    . Insulin granule dynamics in pancreatic beta cells. Diabetologia 2003;46:1029–1045
    OpenUrlCrossRefPubMedWeb of Science
  35. ↵
    1. Pictet RL,
    2. Clark WR,
    3. Williams RH,
    4. Rutter WJ
    . An ultrastructural analysis of the developing embryonic pancreas. Dev Biol 1972;29:436–467
    OpenUrlCrossRefPubMedWeb of Science
  36. ↵
    1. Ganic E,
    2. Johansson JK,
    3. Bennet H,
    4. Fex M,
    5. Artner I
    . Islet-specific monoamine oxidase A and B expression depends on MafA transcriptional activity and is compromised in type 2 diabetes. Biochem Biophys Res Commun 2015;468:629–635
    OpenUrl
  37. ↵
    1. Tonne JM,
    2. Sakuma T,
    3. Deeds MC, et al
    . Global gene expression profiling of pancreatic islets in mice during streptozotocin-induced β-cell damage and pancreatic Glp-1 gene therapy. Dis Model Mech 2013;6:1236–1245
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Ottschytsch N,
    2. Raes A,
    3. Van Hoorick D,
    4. Snyders DJ
    . Obligatory heterotetramerization of three previously uncharacterized Kv channel alpha-subunits identified in the human genome. Proc Natl Acad Sci U S A 2002;99:7986–7991
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Guan B,
    2. Li W,
    3. Li F, et al
    . Sfrp5 mediates glucose-induced proliferation in rat pancreatic β-cells. J Endocrinol 2016;229:73–83
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Wang CX,
    2. Cui GS,
    3. Liu X, et al
    . METTL3-mediated m6A modification is required for cerebellar development. PLoS Biol 2018;16:e2004880
    OpenUrlCrossRef
  41. ↵
    1. Wang Y,
    2. Li Y,
    3. Yue M, et al
    . N6-methyladenosine RNA modification regulates embryonic neural stem cell self-renewal through histone modifications. Nat Neurosci 2018;21:195–206
    OpenUrlCrossRefPubMed
  42. ↵
    1. Yoshihara E,
    2. Wei Z,
    3. Lin CS, et al
    . ERRγ is required for the metabolic maturation of therapeutically functional glucose-responsive β cells. Cell Metab 2016;23:622–634
    OpenUrlCrossRefPubMed
  43. ↵
    1. Riley KG,
    2. Pasek RC,
    3. Maulis MF, et al
    . Connective tissue growth factor modulates adult β-cell maturity and proliferation to promote β-cell regeneration in mice. Diabetes 2015;64:1284–1298
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Zito E,
    2. Chin KT,
    3. Blais J,
    4. Harding HP,
    5. Ron D
    . ERO1-beta, a pancreas-specific disulfide oxidase, promotes insulin biogenesis and glucose homeostasis. J Cell Biol 2010;188:821–832
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Mizuno K,
    2. Fujita T,
    3. Gomi H,
    4. Izumi T
    . Granuphilin exclusively mediates functional granule docking to the plasma membrane. Sci Rep 2016;6:23909
    OpenUrlCrossRefPubMed
  46. ↵
    1. Offield MF,
    2. Jetton TL,
    3. Labosky PA, et al
    . PDX-1 is required for pancreatic outgrowth and differentiation of the rostral duodenum. Development 1996;122:983–995
    OpenUrlAbstract
    1. Ahlgren U,
    2. Jonsson J,
    3. Jonsson L,
    4. Simu K,
    5. Edlund H
    . beta-cell-specific inactivation of the mouse Ipf1/Pdx1 gene results in loss of the beta-cell phenotype and maturity onset diabetes. Genes Dev 1998;12:1763–1768
    OpenUrlAbstract/FREE Full Text
    1. Jonsson J,
    2. Carlsson L,
    3. Edlund T,
    4. Edlund H
    . Insulin-promoter-factor 1 is required for pancreas development in mice. Nature 1994;371:606–609
    OpenUrlCrossRefPubMedWeb of Science
  47. ↵
    1. Sander M,
    2. Sussel L,
    3. Conners J, et al
    . Homeobox gene Nkx6.1 lies downstream of Nkx2.2 in the major pathway of beta-cell formation in the pancreas. Development 2000;127:5533–5540
    OpenUrlAbstract
  48. ↵
    1. Nishimura W,
    2. Bonner-Weir S,
    3. Sharma A
    . Expression of MafA in pancreatic progenitors is detrimental for pancreatic development. Dev Biol 2009;333:108–120
    OpenUrlCrossRefPubMed
  49. ↵
    1. El Khattabi I,
    2. Sharma A
    . Proper activation of MafA is required for optimal differentiation and maturation of pancreatic β-cells. Best Pract Res Clin Endocrinol Metab 2015;29:821–831
    OpenUrl
  50. ↵
    1. Aguayo-Mazzucato C,
    2. Koh A,
    3. El Khattabi I, et al
    . Mafa expression enhances glucose-responsive insulin secretion in neonatal rat beta cells. Diabetologia 2011;54:583–593
    OpenUrlCrossRefPubMedWeb of Science
  51. ↵
    1. Nishimura W,
    2. Kondo T,
    3. Salameh T, et al
    . A switch from MafB to MafA expression accompanies differentiation to pancreatic beta-cells. Dev Biol 2006;293:526–539
    OpenUrlCrossRefPubMedWeb of Science
  52. ↵
    1. De Jesus DF,
    2. Zhang Z,
    3. Kahraman S, et al
    . m6A mRNA methylation regulates human β-cell biology in physiological states and in type 2 diabetes. Nat Metab 2019;1:765–774
    OpenUrl
PreviousNext
Back to top
Diabetes: 69 (8)

In this Issue

August 2020, 69(8)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by Author
  • Masthead (PDF)
Sign up to receive current issue alerts
View Selected Citations (0)
Print
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word about Diabetes.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
m6A mRNA Methylation Controls Functional Maturation in Neonatal Murine β-Cells
(Your Name) has forwarded a page to you from Diabetes
(Your Name) thought you would like to see this page from the Diabetes web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
m6A mRNA Methylation Controls Functional Maturation in Neonatal Murine β-Cells
Yanqiu Wang, Jiajun Sun, Zhen Lin, Weizhen Zhang, Shu Wang, Weiqing Wang, Qidi Wang, Guang Ning
Diabetes Aug 2020, 69 (8) 1708-1722; DOI: 10.2337/db19-0906

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Add to Selected Citations
Share

m6A mRNA Methylation Controls Functional Maturation in Neonatal Murine β-Cells
Yanqiu Wang, Jiajun Sun, Zhen Lin, Weizhen Zhang, Shu Wang, Weiqing Wang, Qidi Wang, Guang Ning
Diabetes Aug 2020, 69 (8) 1708-1722; DOI: 10.2337/db19-0906
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Research Design and Methods
    • Results
    • Discussion
    • Article Information
    • Footnotes
    • References
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Pancreatic β-Cell–Specific Deletion of VPS41 Causes Diabetes Due to Defects in Insulin Secretion
  • Retinol-Binding Protein 4 Activates STRA6, Provoking Pancreatic β-Cell Dysfunction in Type 2 Diabetes
  • Loss of Furin in β-Cells Induces an mTORC1-ATF4 Anabolic Pathway That Leads to β-Cell Dysfunction
Show more Islet Studies

Similar Articles

Subjects

  • Islet Biology-Beta Cell-Development and Postnatal Growth

Navigate

  • Current Issue
  • Online Ahead of Print
  • Scientific Sessions Abstracts
  • Collections
  • Archives
  • Submit
  • Subscribe
  • Email Alerts
  • RSS Feeds

More Information

  • About the Journal
  • Instructions for Authors
  • Journal Policies
  • Reprints and Permissions
  • Advertising
  • Privacy Policy: ADA Journals
  • Copyright Notice/Public Access Policy
  • Contact Us

Other ADA Resources

  • Diabetes Care
  • Clinical Diabetes
  • Diabetes Spectrum
  • Scientific Sessions Abstracts
  • Standards of Medical Care in Diabetes
  • BMJ Open - Diabetes Research & Care
  • Professional Books
  • Diabetes Forecast

 

  • DiabetesJournals.org
  • Diabetes Core Update
  • ADA's DiabetesPro
  • ADA Member Directory
  • Diabetes.org

© 2021 by the American Diabetes Association. Diabetes Print ISSN: 0012-1797, Online ISSN: 1939-327X.