Differential CpG methylation at Nnat in the early establishment of beta cell heterogeneity

Aims/hypothesis Beta cells within the pancreatic islet represent a heterogenous population wherein individual sub-groups of cells make distinct contributions to the overall control of insulin secretion. These include a subpopulation of highly connected ‘hub’ cells, important for the propagation of intercellular Ca2+ waves. Functional subpopulations have also been demonstrated in human beta cells, with an altered subtype distribution apparent in type 2 diabetes. At present, the molecular mechanisms through which beta cell hierarchy is established are poorly understood. Changes at the level of the epigenome provide one such possibility, which we explore here by focusing on the imprinted gene Nnat (encoding neuronatin [NNAT]), which is required for normal insulin synthesis and secretion. Methods Single-cell RNA-seq datasets were examined using Seurat 4.0 and ClusterProfiler running under R. Transgenic mice expressing enhanced GFP under the control of the Nnat enhancer/promoter regions were generated for FACS of beta cells and downstream analysis of CpG methylation by bisulphite sequencing and RNA-seq, respectively. Animals deleted for the de novo methyltransferase DNA methyltransferase 3 alpha (DNMT3A) from the pancreatic progenitor stage were used to explore control of promoter methylation. Proteomics was performed using affinity purification mass spectrometry and Ca2+ dynamics explored by rapid confocal imaging of Cal-520 AM and Cal-590 AM. Insulin secretion was measured using homogeneous time-resolved fluorescence imaging. Results Nnat mRNA was differentially expressed in a discrete beta cell population in a developmental stage- and DNA methylation (DNMT3A)-dependent manner. Thus, pseudo-time analysis of embryonic datasets demonstrated the early establishment of Nnat-positive and -negative subpopulations during embryogenesis. NNAT expression is also restricted to a subset of beta cells across the human islet that is maintained throughout adult life. NNAT+ beta cells also displayed a discrete transcriptome at adult stages, representing a subpopulation specialised for insulin production, and were diminished in db/db mice. ‘Hub’ cells were less abundant in the NNAT+ population, consistent with epigenetic control of this functional specialisation. Conclusions/interpretation These findings demonstrate that differential DNA methylation at Nnat represents a novel means through which beta cell heterogeneity is established during development. We therefore hypothesise that changes in methylation at this locus may contribute to a loss of beta cell hierarchy and connectivity, potentially contributing to defective insulin secretion in some forms of diabetes. Data availability The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD048465. Graphical Abstract Supplementary Information The online version contains peer-reviewed but unedited supplementary material available at 10.1007/s00125-024-06123-6.

methylation at this locus may contribute to a loss of beta cell hierarchy and connectivity, potentially contributing to defective insulin secretion in some forms of diabetes.

Introduction
Insulin-producing beta cells are central to the modulation of glucose homeostasis, and their impaired function, loss of identity or lowered numbers result in type 2 diabetes [1].Previous studies have provided an understanding of the transcriptional machinery that orchestrates beta cell development from early pancreatic and endocrine precursors [2].

Abbreviations
To bolster these transcriptional programmes in vivo, chronic regulation via the epigenome appears to step in to maintain beta cell identity in the long term [3][4][5][6][7].
Imprinted genes are expressed from a single allele in a parent-of-origin-specific manner and their expression is also controlled via epigenetic modifications, notably DNA methylation.Imprinted genes often play key physiological roles, particularly in early (fetal and postnatal) growth and development, controlling a wide range of cellular processes.Thus, human imprinting disorders involving altered expression from specific imprinted loci are associated with severe childhood developmental and metabolic complications (reviewed in [31]).
Imprinted genes play key functional roles in pancreatic beta cells by modulating insulin secretory machinery or beta cell mass [32].Correspondingly, imprinted gene expression is dysregulated both in beta cells with diminished GSIS and in pancreatic islets from individuals with type 2 diabetes, and single nucleotide polymorphisms (SNPs) at imprinted loci are associated with type 2 diabetes risk [18].Monoallelic expression of imprinted genes is maintained transgenerationally by differential methylation between parental alleles at imprinted loci in the germline [33], with additional 'somatic' or 'secondary' differentially methylated regions (DMRs) also established post fertilisation [34].
Here, we explore the possibility that differential methylation of the Nnat gene contributes to the functional heterogeneity of embryonic and adult pancreatic beta cells.

Methods
For details, please refer to electronic supplementary material (ESM) Methods.
Pseudo-time analysis scRNA-seq datasets [39] were processed with Seurat v4.3 and clustered and processed for pseudo-time analysis, followed by integration with a dataset from adult mouse islets [41,42], according to the ESM Methods.

Study approval
All animal procedures were in accordance with the UK Animals (Scientific Procedures) Act 1986 and approved by the UK Home Office (Project licence (PPL) number: PP4712156).All animal studies have been approved by the Imperial College Animal Welfare and Ethical Review Body.
Animal models Transgenic mouse lines expressing Cre recombinase under the control of the rat insulin promoter (RIP) [43] and with a tdTomato reporter downstream of a stop codon flanked by loxP sites [44] have been described previously.Mice with Nnat-driven enhanced GFP (eGFP) expression were purchased from the Mutant Mouse Resource and Research Centre (MMRRC, USA) repository [Tg(Nnat-EGFP)EA106Gsat/Mmucd, stock no.010611-UCD].Mice with global deletion of Nnat [35] and conditional Dnmt3a null mice generated by crossing to Pdx1-Cre are described in the ESM Methods.db/db and control C57BLKS/J mice on the BKS background were obtained from The Jackson Laboratory, USA (JAX 000642 and JAX 000662, respectively).

Primary islet isolation and FACS
To isolate primary islets, the pancreas was inflated with Liberase TM (Roche, UK), digested at 37°C and purified using a Histopaque 1119/1083/1077 gradient (Sigma-Aldrich, UK), and islets were hand-picked.For FACS-based experiments, cells from purified islets were isolated using Accutase (Sigma-Aldrich) and sorted using a FACSAria III flow cytometer (BD Biosciences, UK), as detailed in the ESM Methods.Total insulin content was assessed using ultra-sensitive insulin homogeneous time-resolved fluorescence (HTRF) assay kits (Cisbio, France).
Intracellular calcium imaging Pancreatic islets from reporter mice expressing Nnat-eGFP were isolated as above, with Ca 2+ imaging of whole islets performed after loading the cytosol with 2 μmol/l Cal-590 AM (Stratech, UK).Images were captured on an Axiovert microscope (Zeiss, Germany) equipped with a ×10, 0.3-0.5 numerical aperture (NA) objective and a ImagEM camera (Hamamatsu, Japan) coupled to a Nipkow spinning-disk head (CSU-10, Yokogawa, Japan) and illuminated at 490 nm or 530 nm.
For experiments involving global Nnat null mice [35], islets were incubated with 4.5 μmol/l Cal-520 AM (Stratech), and imaging was performed on a Nikon (Japan) Eclipse Ti microscope equipped with a ×40/1.2NA oil objective and an ibidi heating system.Cal-520 AM was excited with a 491 nm laser line and emitted light filtered at 525/50 nm.Images were acquired with an ORCA-Flash 4.0 camera (Hamamatsu) and Metamorph software (Molecular Devices).Pearson-based connectivity and correlation analyses in an imaged islet were performed with Ca 2+ signals smoothed, binarised and analysed as described in the ESM Methods.

Histological techniques and immunofluorescence
Dissected tissues were fixed, cryoprotected and embedded in optimal cutting temperature (OCT) and stored at −80°C.Sections (10 μm) were cut using a CM1950 Cryostat (Leica, Germany) and immunostained using primary and secondary antibodies (listed in the ESM Methods), and imaged using a TCS SP5 confocal microscope (Leica).Immunostaining of paraffin-embedded pancreas from mice with beta cellselective deletion of Dnmt3a was performed according to references [6,22] as described in the ESM Methods.Human pancreases were processed, immunostained and imaged as described in the ESM Methods with the approval of the Ethics Committee of the University of Pisa, upon written consent of donors' next-of-kin.
Cell culture, RNA silencing and GSIS Human EndoC-βH1 and rat INS1E beta cells were cultured as previously described [35] and incubated with lentiviruses expressing an NNATtargeting short hairpin RNA (shRNA; Sigma-Aldrich) or Silencer Select siRNAs (Ambion, USA), respectively, as described in the ESM Methods.GSIS assays were performed using ultra-sensitive insulin HTRF assay kits (Cisbio) for insulin quantification.
Immunoprecipitation and mass spectrometry Cells were processed for immunoprecipitation and incubated with antibodies against NNAT (ab27266, Abcam, UK, 1:500) prior to mass spectrometry analysis [35] as described in the ESM Methods.
For RNA-seq, RNA from FACS-purified cells was quantified and assessed for integrity using a Bioanalyzer 2100 and an RNA 6000 Pico assay (Agilent Technologies).RNA was processed using a NEBNext Ultra II Directional RNA Library Prep Kit for Illumina paired with a Poly(A) mRNA magnetic isolation module and AMPure XP SPRIselect beads (Beckman Coulter, UK).Libraries were sequenced using a NextSeq 500 High Output sequencer (Illumina, USA), with 2×75 bp length at 50 million reads per sample.Further details can be found in the ESM Methods.
Bisulphite sequencing Genomic DNA was extracted from FACS-purified cells using Allprep DNA/RNA/protein mini kits from Qiagen.Bisulphite sequencing experiments were performed using an EZ DNA Methylation-Gold kit from Zymo Research, USA, an EpiTaq hot start kit (TaKaRa Bio, Japan) and a CloneJET PCR cloning kit (Life Technologies), with DNA purified using a Wizard SV 96 plasmid system (Promega, USA) before Sanger sequencing (Genewiz, now Azenta Life Sciences, UK), and then analysed using bisulphite sequencing DNA methylation analysis (BISMA) software.

Statistical analysis
Data are shown as mean±SEM in all figure panels.Data were assessed using GraphPad Prism 9.0, with <5% error probability considered significant (i.e.p<0.05).Further statistical information, such as n numbers and p values, are provided in the figure legends.

NNAT expression is heterogeneous across the islet and represents a highly differentiated beta cell subtype
We have previously shown that Nnat is crucial for insulin storage and GSIS in the mouse [35].NNAT expression is diminished in rodent models of obesity and diabetes including the Zucker diabetic fatty (ZDF) rat model (85% decrease, p=0.0023, n=5 [45]) and obese ob/ob mice (ESM Fig. 1a).Pancreatic sections of db/db mice displayed a loss of NNAT immunoreactivity across islets, and a significant reduction of NNAT + beta cells (ESM Fig. 1b,c).
Reanalysis of published scRNA-seq datasets from primary mouse beta cells at embryonic (embryonic day 17.5 [E17.5])[39,40] and adult [5,41,42] stages identified Nnat as amongst the most highly variable genes between individual cells (ESM Table 1).Pseudo-time analysis of embryonic (E12.5 to E17.5) mouse beta cells [39] clearly demarcated the differentiation of beta cells through embryonic cell states, with late embryonic (E17.5)beta cells overlapping with high levels of Nnat expression (Fig. 1a-c, ESM Figs 2a-f, 3a-k).Analysis of the proportions of different beta cell progenitors revealed that Nnat + cells become more prominent at later stages of development (Fig. 1d).Integrating cells from the embryonic beta cell trajectory with adult beta cells [41,42] allowed us to evaluate the progressive change in beta cell markers and Nnat expression from development to maturity (Fig. 1e,f, ESM Fig. 4a-i), and revealed large increases in Ins1 expression alongside Nnat downregulation in most beta cells.Nevertheless, a considerable number of adult beta cells (with high expression of Ins1/Ins2) also expressed Nnat (Fig. 1f).Thus, Nnat appears to mark latestage beta cell differentiation, peaking in expression around E17.5 and then gradually being downregulated across most beta cells in adulthood.
We confirmed some of these findings at the protein level using immunofluorescence in postnatal and adult mice (Fig. 3a,b).We noted that NNAT protein expression is highly dynamic in beta cells across the mouse islet throughout the postnatal stage, transitioning from expression throughout the majority (90.4±3.2%,n=15 mice) of beta cells at late embryogenesis to a subset (24.4±3.0%,n=5 mice) of beta cells by postnatal day 14 (P14) (Fig. 3a,b).NNAT + beta cells persisted throughout the islet into adulthood (Fig. 3a,b).NNAT expression was not apparent in other mouse islet cell types, including alpha and delta cells, at P14.However, NNAT expression was detectable in a small number of alpha and delta cells at E17.5 (Fig. 3c).
NNAT immunoreactivity is detected in a subset of beta cells in human islets and NNAT deficiency in human beta cells blunts GSIS While the role of NNAT has previously been examined in mouse beta cells [35], no data currently exist for human beta cells.We therefore next compared the functionality of human NNAT + and NNAT − beta cells.Transient silencing of NNAT in human EndoC-βH1 beta cells completely abrogated glucose-stimulated insulin release (Fig. 4a-c).As in mouse beta cells [35], NNAT interacted in cellulo with subunits of the SPC, SEC11A, SPCS1, SPCS2 and SPCS3 (ESM Fig. 8a), likely via SPCS1 (ESM Fig. 8b).

Heterogeneous NNAT expression in beta cells during postnatal development is associated with altered CpG methylation
To extend the findings above using an orthogonal approach, and to facilitate subsequent functional studies, we utilised a Bacterial Artificial Chromosome (BAC) reporter mouse line expressing eGFP under the control of the Nnat promoter and upstream enhancers (ESM Results, ESM Fig. 10a-e).Aligning with studies of the endogenous gene (Figs 1, 2, 3 and 4), adult reporter mice expressed eGFP in a subpopulation (~15% of total) of beta cells (Fig. 5a).
To isolate and purify beta cells based on NNAT levels, we crossed Nnat-eGFP reporter mice to animals expressing Cre recombinase under the RIP [43] and to transgenic mice expressing a tdTomato reporter downstream of a loxPflanked stop codon [44].Dispersion of primary adult islets into single cells, and subsequent FACS analysis, verified the presence of NNAT + and NNAT − beta cells, with the NNAT + fraction composing 15.3±1.8%(n=8 mice) of the total beta cell compartment (Fig. 5b,c, ESM Figs 11a,b, 12a-f).Collection of FACS-purified NNAT + and NNAT − beta cells and subsequent bisulphite sequencing analysis revealed that CpG methylation at the gametic DMR (known to control monoallelic Nnat expression; Introduction) was unchanged (ESM Fig. 11c-e).Nevertheless, CpG methylation at the Nnat promoter was significantly altered across this genomic region, with minimal CpG methylation observed in NNAT + beta cells (Fig. 5d-f).
Recent whole-genome methylation analysis in mouse germ cells has revealed that the Nnat promoter region is unmethylated in sperm, and this was confirmed by targeted bisulphite sequencing analysis (ESM Fig. 11f-h), whereas it is fully methylated in oocytes [47].Thus, the 'classical' imprinting of Nnat is unaltered between beta cell subtypes with near binary ('on/off') expression of NNAT.However, and overlaying this control, a second DMR at the Nnat promoter within the non-imprinted allele dictates beta cell subtype specificity of expression (Fig. 5g).

The de novo DNA methyltransferase DNMT3A establishes NNAT beta cell subtype specificity
To ask whether this apparent transition in Nnat expression in the early postnatal period is driven by de novo CpG methylation, we assessed endogenous NNAT beta cell immunoreactivity in postnatal mice conditionally deleted for the methyltransferase DNA methyltransferase 3 alpha (DNMT3A) at the pancreatic progenitor stage (using Pdx1-Cre) [6].NNAT staining in control mice at P6 demonstrated expression in a subpopulation of beta cells, whereas deletion of DNMT3A resulted in a loss of this heterogenous expression across the islet (Fig. 6a,b, ESM Fig. 13a).These findings demonstrate that de novo methylation is likely to drive NNAT restriction across the beta cell complement in the first few days of postnatal life.

NNAT + and NNAT − beta cells have distinct transcriptional signatures
In contrast to the embryonic populations described above, RNA-seq revealed that FACS-sorted adult (8 weeks) NNAT + and NNAT − beta cells displayed high transcriptional overlap and were similar to each other when compared with non-beta endocrine cells (Fig. 7a).Major beta cell identity markers such as Ins1, Ins2, Mafa, Slc2a2 (Glut2) and Nkx6.1 were not differentially expressed (ESM Table 5).
The NNAT + beta cell population shows impaired glucose-stimulated Ca 2+ dynamics and is de-enriched for highly connected 'hub' cells To determine whether NNAT may influence beta cell connectivity and membership of the 'hub' cell subgroup [14], we studied glucose-induced Ca 2+ dynamics in Nnat-deficient (Nnat +/−p ) islets [35].Ca 2+ increases in response to high (11 mmol/l) glucose were significantly higher in Nnat +/−p than control islets, and while the mean Pearson's coefficient of correlation and proportion of highly connected 'hub' cells were not significantly different between Nnat-deficient and control islets (ESM Results, ESM Fig. 14a-d), wild-type 'hub' cells were significantly more connected than Nnat-deficient 'hub' cells (ESM Fig. 14e,f).Thus, NNAT is a marker of less wellconnected cells.

Discussion
We describe here a novel subgroup of beta cells characterised by transient methylation of a second DMR (but not the gametic DMR) in the Nnat locus.This is consistent with previous findings that methylation at gametic DMRs is thought to be relatively stable [50], whereas 'secondary' imprinting regions have been shown to be more sensitive to nutrient-or physiology-based changes [6,51].We have previously described the importance of Nnat for beta cell insulin content and secretion [35] and therefore hypothesise that NNAT + and NNAT − cells display differences in secretory function.Correspondingly, transcriptomic analysis of purified adult NNAT + and NNAT − beta cells demonstrated that the NNAT + fraction is enriched for functional pathways including translation initiation, ETC/oxidative phosphorylation pathways and co-translational protein ER membrane targeting.
Extending these studies to humans, NNAT deficiency in human EndoC-βH1 beta cells severely blunted GSIS, likely via inhibition of signal peptidase-mediated processing [35].We show here that E17.5 Nnat + beta cells are enriched for expression of the SPC and translocon apparatus (Spcs1, Spcs2, Sec11a, Sec11c, Sec61) and that the SPC/NNAT interaction is via SPCS1, and likely the means through which NNAT influences cellular insulin content.Importantly, we demonstrate that heterogeneous expression of NNAT is also a feature of the adult human pancreatic islet.Finally, NNAT + cells were de-enriched for highly connected 'hubs' and, thus, are likelier to belong to the 'follower' population [14].Interestingly, Ca 2+ responses to 11 mmol/l glucose were significantly higher in Nnat-deficient vs wild-type mouse islets, in contrast to previous findings where differences were not observed in response to 16.7 mmol/l glucose [35].A similar strong tendency was also seen when comparing NNAT + vs NNAT − cells in the same islet in Nnat-eGFP reporter mice, and suggests that NNAT + cells are less responsive to metabolic stimulation by glucose, in line with their enrichment in a 'follower' subset of cells, and consistent with a role in insulin production rather than glucose detection [14,15].The molecular underpinnings of the weaker Ca 2+ responses are unclear, but do not appear to involve differences in the levels of transcripts encoding Slc2a2 (Glut2) or Gck.Moreover, the transcriptome of NNAT + cells does not show enrichment for genes enriched in 'hub' [15] or 'leader' cells [25].
Interestingly, bimodal expression of Nnat was already clearly evident at embryonic stages, with Nnat + beta cells enriched for markers of late-stage beta cell differentiation, as well as Ins1 and Ins2 mRNAs.These differences were, however, less marked in the adult islet, though CD24a [7] was more weakly expressed in the NNAT + population.De-enrichment of markers of other islet cell types (and the immaturity marker Cd81) was also a common feature of both embryonic and adult NNAT + beta cells.Interestingly, both embryonic and adult NNAT + beta cells were enriched for Npy, suggesting that these cells were not fully matured.NNAT + beta cells, however, had significantly higher insulin content, suggesting a possible, at least partial, overlap with the recently described 'CD63 hi ' population [23], and the 'β HI ' population described in reference [7].
Our work also provides evidence that the epigenome controls the fate of specific beta cell subpopulations.CpG methylation plays a crucial role in early beta cell developmental maturation, including the silencing of 'disallowed' genes such as Hk1, Mct1 (Slc16a1) and Ldha via DNMT3A [4,48].Here we show that islets transition from a state wherein the majority of beta cells are NNAT + in late embryogenesis to comprising a restricted subpopulation of NNAT + beta cells by P7.Whether this represents changes at the level of individual beta cells, or the turnover of the NNAT + population and replacement with a largely NNAT − population, was not determined, and would require investigation using other approaches, including fate mapping (lineage tracing).
Adult NNAT + beta cells are virtually unmethylated at the Nnat promoter whereas robust promoter methylation was apparent in NNAT − cells.These findings, and the fact that the Nnat promoter is differentially methylated between sperm and oocytes (see the Results), indicate that promoter methylation during this transition is likely to be acquired on the paternal allele (from which Nnat is selectively expressed via genomic imprinting).Our studies also show that specific deletion of the de novo methyltransferase DNMT3A at the pancreatic progenitor stage resulted in partial loss of NNAT-based beta cell heterogeneity.Thus, DNA methylation modulates restricted NNAT expression in specific beta cells during early maturation.
Might these findings be relevant for the pathogenesis of type 2 diabetes?Whether demethylation of the second DMR at the Nnat locus may occur in this disease is an interesting possibility which remains to be explored.In rodents, even mild hyperglycaemia deregulates Nnat expression alongside that of several other critical beta cell identity genes [52], and we describe here a significant reduction of NNAT + beta cells in db/db mice.Moreover, altered CpG methylation is a common observation in islets from individuals with diabetes at imprinted and non-imprinted loci (reviewed in [18,32]).
In conclusion, the present work demonstrates how an effector gene in pancreatic beta cells, NNAT/Nnat, is controlled at the level of the epigenome (DNA methylome), contributing to a functional hierarchy between cells.Chemical modification of the epigenome may in future provide an attractive therapeutic angle not only for beta cell replacement or regeneration [53] but also to modulate beta cell function and cell-cell connectivity in type 2 diabetes.

Fig. 1
Fig. 1 Nnat expression is a marker for late-stage beta cell differentiation during islet development.(a-c) UMAP projection of scRNA-seq data from embryonic mouse islets (E12.5-E17.5),with cells labelled by their cellular state (a), pseudo-time (scale bar represents abstract unit of time) (b) and corresponding Nnat expression (log-normalised expression level) (c).(d) Changes in the distribution of beta cell pre-

Fig. 2
Fig. 2 Transcriptomic analyses of islet cells at E17.5 reveal Nnat + cells to be a more differentiated beta cell subcluster.(a, b) Islet cells plotted in the UMAP space and coloured by Ins1 (a) and Nnat (b) expression, respectively.Individual cellular clusters are labelled.Beta cell clusters are highlighted with black circles.(c) Distribution of Nnat expression among beta cell clusters, along with the expression range used in defining Nnat + beta cells (highlighted by the red rectangle).(d) Volcano plot showing the log 2 fold change in gene expression between Nnat + and Nnat − beta cells against FDR-corrected p

Fig. 3 A
Fig. 3 A subpopulation of NNAT + beta cells develops during the early postnatal period in mice.(a) Representative confocal microscopy of pancreatic cryosections from wild-type mice on a C57BL/6J background of developmental stages from E17.5 through the postnatal period into adulthood.Sections were immunostained with antibodies against endogenous NNAT (green) and insulin (INS, red).Nuclei are visualised with DAPI and sections from P56 mice with constituent deletion of Nnat were used as an immunostaining control.Scale bar, 100 μm (inset, 10 μm).(b) Quantification of NNAT + beta cells from images shown in (a), expressed as NNAT/INS co-positive cells

Fig. 4
Fig. 4 NNAT is expressed in a subset of human beta cells and NNAT deficiency in human beta cells reduces insulin secretion.(a) GSIS analysis at low (3 mmol/l) and high (16.7 mmol/l) glucose in human EndoC-βH1 beta cells following 72 h of lentiviral-mediated shRNA knockdown of NNAT with data expressed as insulin secreted into culture media as a percentage of total cellular insulin content (n=8 independent cultures per group, two-way ANOVA with Sidak's multiple corrections test).(b, c) RT-PCR (n=8, unpaired Student's t test) and western blot (n=4, Mann-Whitney U test) analysis of NNAT expression in EndoC-βH1 beta cells after transient NNAT silencing as in (a).(c) Western blotting analysis of NNAT protein levels shown via a representative blot of two independent experiments.β-Tubulin was

Fig. 5
Fig. 5 Beta cell heterogeneity of NNAT expression is associated with changes in CpG methylation at the Nnat promoter.(a) Representative confocal microscopy of pancreatic cryosections from P56 mice with Nnat-driven eGFP expression from the paternal allele (Nnat WT/eGFPpat ) (n=7 mice).Sections were immunostained with antibodies against eGFP (green) and insulin (INS, grey).Nuclei are visualised with DAPI.Scale bar, 50 μm.(b) FACS separation of dispersed primary islet cells from reporter mice with insulindriven expression of tdTomato (to label beta cells) and Nnat-driven eGFP expression from the paternal (Nnat WT/eGFPpat ) or maternal (Nnat WT/eGFPmat ) allele or wild-type (Nnat WT/WT ) at this locus (representative FACS plot of the dispersed islet preparation from a single mouse shown) (n=8, 3 and 3 mice per genotype, respectively, Kruskal-Wallis test with Dunn's multiple comparisons).(c) Quan-

Fig. 6
Fig. 6 Postnatal restriction of NNAT in a subset of beta cells is at least partially driven by the de novo methyltransferase DNMT3A.(a) Representative confocal microscopy of pancreatic cryosections from mice with conditional deletion of DNMT3A under the control of the Pdx1 promoter (Pdx1-Cre + Dnmt3a fl/fl ) vs control (Pdx1-Cre − Dnmt3a fl/fl ) mice at P6.Sections were immunostained with antibodies

Fig. 7
Fig. 7 NNAT + adult beta cells are transcriptionally distinct and have significantly higher insulin content.(a) Correlation matrix of differentially expressed genes between NNAT + and NNAT − beta cells as assessed by RNA-seq analysis (n=4 FACS-purified populations from individual mouse islet preparations).Scale bar (0-255) gives an integer provided to the R function where the darker colour corresponds to higher sample correlation.(b) Heatmap of the top ten most differentially expressed genes reduced in NNAT + beta cells compared with NNAT − beta cells.(c, d) GSEA showing categories significantly enriched (c) and reduced (d) in NNAT + (vs NNAT − ) beta cells.(e, f) Representative confocal microscopy of pancreatic cryosections from P56 (8-week-old) wild-type mice on a C57BL/6J background

Fig. 8
Fig. 8 NNAT + beta cells show altered glucose-induced Ca 2+ dynamics and are de-enriched for highly connected 'hub' cells within individual islets.(a) Ca 2+ -bound Cal-590 AM fluorescence in response to high glucose (11 mmol/l) in NNAT − and NNAT + cells from primary islets from Nnat WT/eGFPpat reporter mice expressed as normalised intensity over time (F/F min ) (n=61 islets total from five mice per genotype; quantification of AUC on the right, p=0.063,Wilcoxon matched-pairs signed rank test).(b) Representative Cartesian map of beta cells with colour-coded lines connecting cells according to the strength of coactivation (colour-coded R values from 0 to 1, blue to red).Beta cells are represented by differently coloured nodes depending on their coactivity with the other beta cells, where black nodes indicate coactivity with ≥80% of the remaining beta cells, while grey, white and orange nodes represent coactivity with ≥60%, ≥40% and <40%, respectively.Nodes circled with a solid black line indicate