A novel disease mechanism leading to the expression of a disallowed gene in the pancreatic beta-cell identified by non-coding, regulatory mutations controlling HK1
===================================================================================================================================================================

* Matthew N. Wakeling
* Nick D. L. Owens
* Jessica R. Hopkinson
* Matthew B. Johnson
* Jayne A.L. Houghton
* Antonia Dastamani
* Christine S. Flaxman
* Rebecca C. Wyatt
* Thomas I. Hewat
* Jasmin J. Hopkins
* Thomas W. Laver
* Rachel Van Heugten
* Michael N. Weedon
* Elisa De Franco
* Kashyap A. Patel
* Sian Ellard
* Noel G. Morgan
* Edmund Cheesman
* Indraneel Banerjee
* Andrew T. Hattersley
* Mark J. Dunne
* International Congenital Hyperinsulinism Consortium
* Sarah J. Richardson
* Sarah E. Flanagan

## Abstract

Gene expression is tightly regulated with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function. This silencing is largely controlled by non-coding elements and their disruption might cause human disease. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding *de novo* mutations affecting a 42bp conserved region encompassed by a regulatory element in intron 2 of *Hexokinase 1* (*HK1*), a pancreatic beta-cell disallowed gene. We demonstrated that these mutations resulted in expression of *HK1* in the pancreatic beta-cells causing inappropriate insulin secretion and congenital hyperinsulinism. These mutations identify a regulatory region critical for cell-specific silencing. Importantly, this has revealed a new disease mechanism for non-coding mutations that cause inappropriate expression of a disallowed gene.

Genetic discovery in Mendelian disease has focused on identifying highly penetrant mutations affecting the function of genes expressed in clinically affected tissue(s). Whilst this approach has proven successful, the underlying aetiology of over 3000 presumed monogenic diseases remains undefined1-4, of which many have significant clinical and genetic heterogeneity5,6. Congenital Hyperinsulinism (CHI), is characterised by inappropriate insulin secretion during hypoglycaemia. It is a clinically and genetically heterogeneous disease where, despite extensive sequencing efforts, the underlying aetiology is not known in up to 50% of individuals7,8.

To identify new aetiologies, we performed whole-genome sequencing (WGS) on 135 individuals with biochemically confirmed, persistent CHI without an identified mutation in a known disease causing gene (**Supplemental Table 1**). We initially searched for genes containing novel coding *de novo* or biallelic variants in two or more individuals. When no such genes were found, we turned to the non-coding genome and searched for *de novo* copy number variants. This identified a single locus within intron 2 of Hexokinase 1 (*HK1*) containing two heterozygous ∼4.5Kb deletions with a 2,787bp overlap in two unrelated individuals (**Supplemental Figure 1**). Digital droplet PCR confirmed the deletions in both probands and an affected twin, in keeping with a germline mosaic mutation.

We next searched for *de novo* single nucleotide variants and indels within the minimal deleted region within our WGS discovery cohort and in a replication cohort of 27 individuals with CHI who had undergone pancreatectomy as part of routine clinical care (**Supplemental Table 1**). This identified 7 different *de novo* mutations in 12 probands. In two cases the mutation had been inherited by similarly affected offspring (**Figure 1**). All mutations were novel9,10 and within a 42bp conserved region constrained against variation in gnomAD.11

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/05/2021.12.03.21267240/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/12/05/2021.12.03.21267240/F1)

Figure 1. 
**a**) Schematic representation of the *HK1* gene (GRCh37/hg19 chr10:71,029,756-71,161,637). The full-length testis-specific isoform (ENST00000448642) and the shorter ubiquitously-expressed isoform (ENST00000359426) are depicted. The positions of the two ∼4.5Kb deletions within intron 2 of the shorter isoform and the 7 different heterozygous point mutations/indels within the minimal deleted region identified in 14 probands are shown. The number of probands with each mutation is provided. An asterisk (*) denotes a single proband with two *de novo* mutations (see **Supplementary Figure 7**). The mutation positions are given according to the GRCh37/hg19 genomic coordinates. b) Partial pedigrees showing inherited *HK1* mutations in 3 families. Filled symbols represent individuals with congenital hyperinsulinism. (‡) Haplotype analysis confirmed that the *HK1* mutation had arisen *de novo* in patient 3.2. M, *HK1* mutation; N, no mutation; NA, DNA not available. Pedigrees for the 11 probands with *de novo HK1* mutations are not shown.

Overall, the finding of 9 different *de novo* mutations in 14 probands, co-segregating with disease in three additional family members provides overwhelming evidence of disease-causality.

All 17 individuals with a *HK1* non-coding mutation had severe early-onset CHI (median age at diagnosis: Birth [IQR: 0-14 days]). These mutations were not associated with macrosomia at birth (median birthweight z score: 0.61 [IQR: 0.10-1.84]) suggesting that insulin secretion was not markedly increased *in utero*. In all cases the hyperinsulinism persisted, with the eldest individual still known to be affected being 18 years of age. These findings suggest that the mutations act to disrupt glucose-induced insulin secretion throughout post-natal life and into adulthood but have less impact during fetal development.

The *HK1* mutations caused a beta-cell specific defect with no common additional features between patients. In five individuals, medical management was ineffective, leading to pancreatic resection. Of these, the only case that had resolution of the CHI developed insulin-dependent diabetes following near-total pancreatectomy. Histopathological analysis of resected pancreatic tissue (n=2) demonstrated a discrete pathology, distinguishable from diffuse disease resulting from *ABCC8* mutations, the commonest cause of CHI (**Supplemental Figure 2**). An overview of the clinical features is provided in **Supplemental Table 2**.

*HK1* encodes a glycolytic enzyme that is silenced in the pancreas and liver but expressed in all other mature tissues (**Supplemental Figure 3a-c**) where it supports glucose metabolism, ensuring cell survival12. Within the pancreatic beta-cell, glucose is phosphorylated to glucose-6-phosphate by Glucokinase (GCK; Hexokinase 4) which, due to a low-binding affinity (Km ∼8mM) acts as the pancreatic glucose-sensor coupling insulin release to the prevailing glucose concentration13. Hexokinase 1 has a markedly higher affinity for glucose (Km <50μM)13. Silencing of *HK1* in favour of *GCK* in the beta-cell therefore ensures appropriate glucose-sensing, minimising insulin release at low glucose levels. Whilst bi-allelic loss-of-function *HK1* mutations have been reported to cause non-spherocytic haemolytic anaemia14, dominant or recessive coding mutations have not been described in individuals with defects in glucose homeostasis; in keeping with the absence of *HK1* protein in beta-cells.

In resected pancreatic tissue of individuals with non-coding mutations, we found that *HK1* was expressed and co-localised with insulin in the islets of affected tissue but not in unaffected controls (**Figure 2**). *HK1* did not co-localise with glucagon (secreted by pancreatic alpha-cells) suggesting that the impact of the mutations was beta-cell specific (**Supplemental Figure 4**). These results confirmed that the mutations cause *HK1* to be inappropriately expressed in the pancreatic beta-cell and explain why there is increased insulin secretion in these individuals during hypoglycaemia.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/05/2021.12.03.21267240/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2021/12/05/2021.12.03.21267240/F2)

Figure 2: *HK1* expression is present in the beta-cells of donors with *HK1*-CHI, but not controls.
**a**) Staining of HK1 (red); insulin (blue) and DAPI (white) in pancreatic tissue resected from one patient with a *HK1* mutation. Within this section HK1co-localises with insulin (purple staining) whereas, no *HK1* expression is observed in the non-hyperinsulinism (Non-HI) control donor. Scale bar – 100μm. **b**) HK1 expression (Median fluorescence intensity (MFI)) is significantly increased in the *HK1* HI donors (n=2 donors; 17015 beta-cells), when compared to Non-HI control donors (n=2 donors; 21408 beta-cells). Two-tailed Mann-Whitney (U=590830; \**\*|\*p<0.0001).

We identified an islet enhancer encompassing our critical region15 bound by a broad set of islet transcription factors, with NKX2-2 and FOXA2 most prominent (**Figure 3**). Single-nuclei ATAC-seq (snATAC-seq) data in human islets16 revealed a peak of open chromatin in beta-cells and showed that the relevant *HK1* promoter remains accessible (**Figure 3b**). Human islet Hi-C data17 confirmed that both our region and the relevant *HK1* promoter are contained within a well-insulated domain, without evidence of contact to distal loci (**Figure 3d**) further supporting direct regulation by the critical region on the promoter.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/05/2021.12.03.21267240/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/12/05/2021.12.03.21267240/F3)

Figure 3: 
**a**) Expression of *HK1* and GCK for comparison over the course of pancreatic differentiation from embryonic stem cells to maturing beta-cells, expression data over a beta-cell differentiation pseudotime18, Gaussian process regression median and 95% CI shown. DE – definitive endoderm, GT – gut tube, FG – foregut, PE – pancreatic endoderm, EP – endocrine precursors, IB – immature beta-cells. See **Supplemental Figure 3** for further *HK1* expression. **b**) Transcription factor binding in human islets15 and chromatin accessibility from snATAC-seq in the cluster of cells assigned to beta-cells16 over *HK1* locus, the region containing the mutations is marked as grey box. Chromatin accessibility in other endocrine cell types is provided in **Supplemental Figure 5. c**) Chromatin accessibility during pancreatic differentiation over critical region, stages: ES – embryonic stem cell, DE – definitive endoderm, GT – gut tube, Early/Late PP – pancreatic progenitor, describe ATAC-seq from in vitro differentiation. Beta – gives accessibility in snATAC-seq16 (as shown in B). **d)** HiC data in human islets17 reveals *HK1* is contained within a well-insulated domain, heatmap gives log scale contact frequencies, white triangle and circles mark chromatin loops called in same study. **e)** Transcription factor motif families disrupted by the mutations, shown are those motif families maximally disrupted by each mutation (red font). Notable secondary motif families given in **Supplemental Figure 7 and Supplemental Table 4**. Left shows scatter of normalised motif scores for the reference versus mutated score (motif score normalised by maximal motif score), black line gives 1-1 line of equal motif. Number gives final two digits of chromatin position, e.g. 45 – 71108645. Right gives sequence logo for the three maximally disrupted families. Members of each family are expressed in beta-cells **Supplemental Figure 6**.

During development, *HK1* is highly expressed in embryonic stem cells but is progressively downregulated during pancreatic cell-type differentiation becoming extinguished during beta-cell maturation18,19 (**Figure 3a, Supplemental Figure 3d**). GCK shows increased expression during differentiation as *HK1* is downregulated (**Figure 3a**). Despite being absent from beta-cells *HK1* is expressed in mature stellate- and duct-cells suggesting that it becomes de-repressed in other pancreatic cell-types (**Supplemental Figure 3e**).

Chromatin accessibility data shows our critical region remains closed over pancreatic-cell differentiation until pancreatic progenitor stages, when it becomes accessible and bound by a complement of factors in the pancreatic progenitor regulatory network20,21 (**Figure 3c, Supplemental Figure 5b-c**). This suggests the region is necessary for HK1 repression in late pancreatic development and mediating beta-cell specific control in adult cells but that it is not required for HK1 repression during early differentiation.

Analysis of histone modifications during pancreatic differentiation, in human islets20 and the EndoC-βH1 beta-cell line22, revealed that repression is actively maintained. The critical region is marked by a bivalent state encompassed by focal peaks of active enhancer marks H3K4me1 and H3K27ac and a broader polycomb H3K27me3 domain (**Supplemental Figure 5d**).

Analysis of transcription factor motif families revealed that mutations disrupt transcription factor binding sites for FOX, NKX2 and NFAT families (**Supplemental Table 4, Figure 3c**), which each have family members expressed in beta-cells (**Supplemental Figure 6**). This is supported by ChIP-seq data revealing that NKX2-2 and FOXA2 were prominently bound at this region in islets (**Figure 3b**). FOXA2 loss-of-function causes CHI23, and NKX2-2 is part of a large repressor complex that regulates beta-cell specification24 and could therefore have a role in maintaining the repressive epigenetic state at the HK1 locus. Multiple members of the NFAT family are expressed in beta-cells (**Supplemental Figure 6**), with NFATC2 inducing beta-cell proliferation in human islets25. Whilst NFATC2 is expressed at low levels in adult beta-cells it has greater expression in pancreatic progenitors (**Supplemental Figure 6**), supporting HK1 repression as mediated by loss of NFAT binding in late stages of beta-cell development and proliferation. Analysis of snATAC-data from human islets revealed endocrine cells present in hormone-high and hormone-low states16, with the former associated with increased promoter accessibility over secreted endocrine genes and the latter over cell-cycle genes. Our critical region is accessible in the hormone-high state of alpha, beta, and delta-cells, and only remains accessible in beta-cells in the hormone low state (**Supplemental Figure 5a**). Therefore, accessibility may be maintained at the critical region during cell proliferation, by factors including NFATC2, thereby ensuring HK1 is repressed in beta-cells.

Discovering non-coding regulatory mutations in *HK1* enabled us to identify the first regulatory element critical for the selective silencing of an otherwise ubiquitously expressed gene within a single cell-type. It is interesting to note that interruption of transcription factor binding is implicated in each of the 9 distinct mutations and the locus is decorated by both active and repressive epigenetic marks suggesting that continual transcription factor binding is necessary to maintain repression in beta-cells. Future work should prioritise regulatory regions surrounding disallowed genes with similar epigenetic marks and to explore whether different regulatory mechanisms exist to silence disallowed genes within the beta-cell.

Over 60 beta-cell disallowed genes have been described, although the mechanism(s) controlling cell-specific silencing in humans have not been fully determined26,27. Linkage analysis of *HK1* to CHI28 and increased HK1 expression29 have been reported in several patients, but disease causality and genetic aetiology were not established. Novel promoter variants in the beta-cell disallowed gene, *SLC16A1* have also been reported in two families with exercise–induced hyperinsulinism. However, these variants increased transcription across cell types leading to the hypothesis that differences in the post-transcriptional regulation of mRNA across tissues could explain the beta-cell-specific phenotype30.

The identification of non-coding *HK1* mutations in patients with CHI represents one of the first discoveries of regulatory non-coding mutations affecting a gene in which coding mutations do not cause the same phenotype31. These findings highlight a role for undiscovered regulatory mutations causing disease through inappropriate expression of a normally functioning protein in a specific cell-type.

Mutations affecting the regulation of *HK1* cause CHI through its aberrant expression in beta-cells. These findings are important for future efforts to discover non-coding regulatory mutations as they establish a critical role for disallowed genes in Mendelian disease.

## Supporting information

Supplemental Data [[supplements/267240_file03.pdf]](pending:yes)

Supplemental Table 4 [[supplements/267240_file04.xlsx]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors.

## Competing Financial Interests

The authors declare no competing financial interests.

## Author Contributions

S.E.F. designed the study. A.D, S.E., I.B., S.E.F. and the International Congenital Hyperinsulinism Consortium recruited patients to the study. Ma.N.W*, M.B.J., J.A.L.H., R.V-H., S.E. and S.E.F. performed the molecular genetic analysis and the interpretation of resulting data. Ma.N.W. and Mi.N.W. performed the bioinformatic analysis, T.I.H., J.H., I.B., A.T.H. and S.E.F. analysed the clinical data. E.C. and M.J.D. performed the pathological analysis of the pancreatic tissue, J.R.H, C.S.F., R.W., N.G.M. and S.J.R. designed and performed the immunohistochemistry studies, N.D.L.O. designed and performed the epigenomic data analysis, Ma.N.W., N.D.L.O., J.R.H, M.B.J., T.W.L., E.D-F., K.A.P., A.T.H., M.J.D., S.J.R and S.E.F. prepared the draft manuscript. All authors contributed to the discussion of the results and to manuscript preparation.

## ONLINE METHODS

### Subjects

Congenital hyperinsulinism (CHI) was defined as an inappropriately high level of plasma insulin at the time of hypoglycaemia associated with inappropriately supressed ketones and free fatty acids presenting within the first 12 months of life. The definition of hypoglycaemia was based on the recommendations of the Pediatric Endocrine Society (blood glucose <2.8mmol/L in the presence of detectable insulin)32. Clinical details of the subjects are provided in **Supplemental Table 1**. Subjects with CHI were recruited by their clinicians for molecular genetic analysis to The Exeter Genomics Laboratory. Disease-causing mutations in the known CHI genes had been excluded by targeted-next generation sequencing in all cases33.

### Ethical considerations

This study was approved by the North Wales Research Ethics Committee (517/WA/0327) and was conducted in accordance with the Declaration of Helsinki, with all subjects or their parents providing informed consent for genetic testing. All tissue samples were studied with full ethics approval (West of Scotland Research Ethics Committee, reference: 20/WS/0074 IRAS project ID: 283620 or nPOD)34.

### Whole-genome sequencing

Whole-genome sequencing of DNA extracted from peripheral blood leukocytes was completed on 135 probands with CHI (n=3 with Illumina HiSeq 2500, n=69 with Illumina HiSeq X10, n=63 with BGISeq-500). The mean read depth across the whole genome was 36.9 (stddev 4.9). An additional 191 family members were also sequenced. The sequence data was aligned using BWA MEM 0.7.15, and processed using a pipeline based on the GATK best practices (Picard version 2.7.1, GATK version 3.7). Variants were annotated using Alamut batch standalone version 1.11 (Rouen, France).

CNVs were called across the genome by SavvyCNV35 using a bin size of 2Kbp. We also developed a new tool (FindLargeInsertSizes, available at [https://github.com/rdemolgen/SavvySuite](https://github.com/rdemolgen/SavvySuite)) to detect small (≥1Kbp) deletions, insertions, inversions, and translocations using read pair information. We used this tool to screen all whole genome sequenced samples for structural variants within *HK1* intron 2.

When a mutation was identified, Sanger sequencing (for point mutations/indels) or digital droplet PCR (for deletions) was performed on samples from each available family member to confirm genome sequencing results.

### Digital droplet PCR (ddPCR)

was performed on leukocyte DNA from patients 1 and 2.1 and their unaffected parents to confirm the results of the whole-genome sequencing deletion analysis. ddPCR was also performed on the affected twin brother of patient 2.1 and the 27 individuals with CHI of unknown cause who had undergone pancreatectomy.

Relevant primers (**Supplemental Table 3**) were used in reactions for droplet generation, PCR and detection using the Bio-Rad QX200 ddPCR EvaGreen system (Hercules, CA, USA) to search for copy number changes. Reactions were performed as per the manufacturer’s recommendations with an annealing/extension temperature during PCR of 59°C. All data were analysed using Bio-Rad QuantaSoft software.

### Sanger sequencing of the *HK1* regulatory element

Sanger sequencing was performed on DNA extracted from peripheral blood leukocytes from 27 individuals with CHI of unknown cause who had undergone pancreatectomy. Briefly, a 397bp genomic region (GRCh37/hg19, Chr10:71,108,536-71,108,932) within intron 2 of *HK1* (NM_033497) was amplified by PCR (primer sequences are listed in **Supplemental Table 3**). PCR products were sequenced on an ABI3730 capillary machine (Applied Biosystems, Warrington, UK) and analysed using Mutation Surveyor v3.24 software (SoftGenetics, State College, PA, USA). When a mutation was identified, samples from family members were tested to investigate co-segregation and microsatellite analysis using the PowerPlex kit (Promega, Southampton, UK) was performed to confirm family relationships.

### Histopathology

Formalin-fixed paraffin-embedded pancreatic tissue was available for immediate analysis from two individuals with *HK1* mutations age-matched controls (n=4) and age-matched tissues from patients with *ABCC8* mutations (n=2). Immunohistochemistry was performed as described previously on 5-mm-thick sections of tissue36. For high-content quantification of tissue, sections were first digitized using the 3DHistech Pannoramic 250 Flash II slide scanner (3DHISTECH, Budapest, Hungary) and quantified using QuPath37.

### Immunofluorescence

After dewaxing and rehydration, pancreatic samples from two *HK1* patients and 4 age-matched control tissues (patients with *ABCC8* mutations (n=2); non-hyperinsulinism donors (n=2)) were subjected to heat–induced epitope retrieval (HIER) in 10mM citrate pH6 buffer for 20 minutes. After blocking, the sections were probed in a sequential manner with rabbit monoclonal anti-hexokinase 1 (Abcam ab150423; 1/100 overnight); mouse monoclonal anti-glucagon (Abcam ab10988; 1/2000 for 1h) and guinea-pig anti-insulin (Agilent; C#IR002; 1/5 for 1h). The relevant antigen-antibody complexes were detected using secondary antibodies conjugated with fluorescent dyes (Alexa Fluor™ anti-mouse 555, anti-rabbit 488, anti-guinea pig 647) (Invitrogen, Paisley, U.K). Cell nuclei were stained with DAPI. After mounting, the sections were imaged via a Leica DMi8 confocal microscope (Leica Microsystems UK, Milton Keynes, UK) and the distribution of HK1, insulin and glucagon examined in multiple islets. In addition, sections were scanned at 40x magnification using an Akoya Biosciences Vectra® Polaris™ Automated Quantitative Pathology Imaging System. The quantification of the HK1 expression in islets was determined using the Random Forest Classifier Module (Version 3.2.1851.354), DenseNet AI V2 and HighPlex FL v4.04 modules included in Indica Labs HALO Image analysis platform (Version 3.2.1851.354). The Random Forest Classifier was used to identify the islets, then the DenseNet AI V2 module was used to identify beta and alpha cells within the islets. This enabled exclusion of HK1 positive non-islet cells such as red blood cells, stellate cells and endothelial cells which are frequently found in proximity to the islet. The median fluorescence intensity (MFI) of HK1 expression in beta and alpha cells were then calculated using the HighPlex FL v4.04 module. This pipeline was applied to all corresponding tissue sections. GraphPad Prism V9.2.0 (332) was used to demonstrate the median and IQR of beta and alpha cells MFI in each of the donors assessed.

### Epigenomic analysis

#### Public Bulk ChIP and ATAC-seq

For ChIP-seq and ATAC-seq datasets reads were downloaded from accessions provided in (**Supplemental Table 4**), and then aligned with Bowtie238 to the GRC37/hg19 genome with Bowtie2 v2.3.5.1 with default parameters for single end reads and with additional options “- I 0 -X 1000 --no-discordant --no-mixed” for paired-end reads. Alignments were filtered for those with mapping quality > 30 and then reads with identical aligning coordinates were treated as duplicates and collapsed to a single alignment. All data is visualised as reads per million using [https://github.com/owensnick/GenomeFragments.jl](https://github.com/owensnick/GenomeFragments.jl). For human islet single-nuclei ATAC-seq data16 (GSE160472), reads from the authors’ combinatorial barcode approach were aligned as above and separated into cell types using author provided cluster labels: v[https://github.com/kjgaulton/pipelines/tree/master/islet\_snATAC_pipeline](https://github.com/kjgaulton/pipelines/tree/master/islet_snATAC_pipeline).

Human islet Hi-C data was obtained from experiment accession TSTSR043623 and file accession DFF064KIG (.hic file) and TSTFF938730 (bedpe file)17. The .hic file from experiment was downloaded and converted to .cool format using hic2cool [https://github.com/4dn-dcic/hic2cool](https://github.com/4dn-dcic/hic2cool) with default parameters. Contact matrices were obtained using cooler39 at 5kb resolution using KR balance and log-transformed for visualisation.

#### scRNA-seq Expression Datasets

scRNA-seq data collected over a time course of pancreatic differentiation18 projected onto a differentiation pseudotime was obtained from the reference’s supplementary table. We identity consistent temporal trends using Gaussian Process (GP) regression, following the approach we have previously applied40. Briefly, to stabilise variance we transform data log(αy+ β) where y is the gene expression for a given gene and α = 100, β = 1, then perform GP regression using Matern52 kernel and invert the transformation to report GP median and 95% confidence intervals. All GP regression was performed with GaussianProcesses.jl ([https://github.com/STOR-i/GaussianProcesses.jl](https://github.com/STOR-i/GaussianProcesses.jl); [https://arxiv.org/abs/1812.09064](https://arxiv.org/abs/1812.09064)).

For scRNA-seq data in human islets41 accession GSE101207 we gene counts per cell for the size healthy donors and normalised by depth per cell. In **Supplemental Figures 3, 6 and 7** we show boxplots of normalised counts for each gene over healthy donors.

Finally, we obtained data on gene expression during beta-cell maturation from the Broad Single Cell portal19. In **Supplemental Figures 3, 6 and 7** we show mean normalised counts for each gene over annotated cell types.

### Motif analysis

To assess motifs maximally disrupted by the mutations, we took all motifs in the non-redundant JASPAR database42 and found the maximal scoring match that spanned the position of each mutation for both reference and mutated sequences using [https://github.com/exeter-tfs/MotifScanner.jl](https://github.com/exeter-tfs/MotifScanner.jl). Motif scores are likelihood-ratio scores of motif PWM over a background of the A,C,G,T frequencies in the hg19 genome. We normalised all motif scores, by the maximal score for each PWM and to both primary and secondary candidate motifs that may be disrupted by the mutations. We considered two tiers; Tier 1 having normalised motif score ≥ 0.6, and Tier 2 having 0.45 ≤ normalised motif score < 0.45. For each tier we calculated a disruption score by subtracting the mutated sequence score from the reference sequence score and ranked these across all motifs. To report a single disrupted motif family, we grouped multiple motifs with overlapping alignments by family by removing any trailing numbers from the gene symbol (e.g. NKX2-2 → NKX2, HIC2 → HIC) and all FOX factors we grouped in the FOX family. Tier 1 disrupted motif families include NFAT, NKX2 and FOX (**Figure 3e**). One has two *de novo* mutations, the leftmost is shared with other patients and is a candidate for NFAT disruption, and the rightmost interrupts a HIC family motif (**Supplemental Figure 7a**). HIC family members include HIC1 and HIC2, with HIC2 the most greatly expressed in beta-cells (**Supplemental Figure 7b**). Interestingly, HIC2 is a transactivator of SIRT143 and the loss of SIRT1 impairs glucose sensing in beta-cells in mice44. Tier 2 disrupted motif families include TEAD and SMAD (**Supplemental Figure 7c**). The TEAD family motif shares TTCA consensus with the NFAT family and is an alternative candidate to NFAT, TEAD1 is expressed in beta-cells (**Supplemental Figure 7d**) and plays a critical role in pancreatic progenitors21, however it should be noted that TEAD1 does not bind the critical region in pancreatic progenitors when the region is bound by FOXA2 (**Figure 3b**). Multiple members of the SMAD family are expressed in beta-cells (**Supplemental Figure 7e**), SMAD factors are signal transducers of TGF-beta signalling and play an important role in beta-cell development, function, and proliferation45.

## Acknowledgements

SEF has a Sir Henry Dale Fellowship jointly funded by the Wellcome Trust and the Royal Society (105636/Z/14/Z). MaNW and MBJ are the recipients of an Independent Fellowship and TWL and NO a Lectureship from the Exeter Diabetes Centre of Excellence funded by Research England’s Expanding Excellence in England (E3) fund. EDF is the recipient of a Diabetes UK RD Lawrence fellowship (19/0005971). KAP has a Career Development fellowship funded by the Wellcome Trust (219606/Z/19/Z). IB and MJD were supported by the Northern Congenital Hyperinsulinism (NORCHI) charitable fund, by the Manchester Academic Health Sciences Centre and by The University of Manchester MRC Confidence in Concept (CiC) Award (MC_PC_18056). SJR and NGM are grateful to JDRF for funding (2-SRA-2018-474-S-B). Ivo Barić is member of the European Reference Network for Rare Hereditary Metabolic Disorders (MetabERN), Project ID No 739543.This research was performed with the support of the Network for Pancreatic Organ donors with Diabetes (nPOD; RRID: SCR_014641), a collaborative type 1 diabetes research project supported by JDRF (nPOD: 5-SRA-2018-557-Q-R) and The Leona M. & Harry B. Helmsley Charitable Trust (Grant#2018PG-T1D053, G-2108-04793). The content and views expressed are the responsibility of the authors and do not necessarily reflect the official view of nPOD. Organ Procurement Organizations (OPO) partnering with nPOD to provide research resources are listed at [http://www.jdrfnpod.org/for-partners/npod-partners/](http://www.jdrfnpod.org/for-partners/npod-partners/). This research has been conducted using the UK Biobank Resource. This work was carried out under UK Biobank project number 9055 and 9072. Targeted next-generation sequencing to exclude known disease-causing mutations was funded by Congenital Hyperinsulinism International (a 501(c)3 organisation) for 20 individuals within the discovery cohort.

## Footnotes

*   This research was funded in whole, or in part, by Wellcome [105636/Z/14/Z]. For the purpose of open access, the author has applied a CC BY public copyright licence to any Author accepted Manuscript version arising from this submission

*   Received December 3, 2021.
*   Revision received December 3, 2021.
*   Accepted December 5, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/)

## References

1.  1.Ferreira, C.R. The burden of rare diseases. Am J Med Genet A 179, 885–892 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.a.61124&link_type=DOI) 

2.  2.Online Mendelian Inheritance in Man, OMIM. McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University (Baltimore, MD), Accessed 18/11/2021. World Wide Web URL: [https://omim.org/](https://omim.org/).
    
    

3.  3.Turro, E. et al. Whole-genome sequencing of patients with rare diseases in a national health system. Nature 583, 96–102 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2434-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32581362&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

4.  4.Spielmann, M. & Mundlos, S. Looking beyond the genes: the role of non-coding variants in human disease. Hum Mol Genet 25, R157–R165 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddw205&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27354350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

5.  5.Shashi, V. et al. The utility of the traditional medical genetics diagnostic evaluation in the context of next-generation sequencing for undiagnosed genetic disorders. Genet Med 16, 176–82 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/gim.2013.99&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23928913&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

6.  6.Sawyer, S.L. et al. Utility of whole-exome sequencing for those near the end of the diagnostic odyssey: time to address gaps in care. Clin Genet 89, 275–84 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cge.12654&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26283276&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

7.  7.Kapoor, R.R. et al. Clinical and molecular characterisation of 300 patients with congenital hyperinsulinism. Eur J Endocrinol 168, 557–64 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZWplIjtzOjU6InJlc2lkIjtzOjk6IjE2OC80LzU1NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzA1LzIwMjEuMTIuMDMuMjEyNjcyNDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

8.  8.Snider, K.E. et al. Genotype and phenotype correlations in 417 children with congenital hyperinsulinism. J Clin Endocrinol Metab 98, E355–63 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2012-2169&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23275527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316270900022&link_type=ISI) 

9.  9.Halldorsson, B.V. et al. The sequences of 150,119 genomes in the UK biobank. bioRxiv, 2021.11.16.468246 (2021).
    
    

10. 10.Lek, M. et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 536, 285–91 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature19057&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27535533&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000381804900026&link_type=ISI) 

11. 11.Karczewski, K.J. et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

12. 12.Robey, R.B. & Hay, N. Mitochondrial hexokinases, novel mediators of the antiapoptotic effects of growth factors and Akt. Oncogene 25, 4683–96 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.onc.1209595&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16892082&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239687000006&link_type=ISI) 

13. 13.Becker, T.C., BeltrandelRio, H., Noel, R.J., Johnson, J.H. & Newgard, C.B. Overexpression of hexokinase I in isolated islets of Langerhans via recombinant adenovirus. Enhancement of glucose metabolism and insulin secretion at basal but not stimulatory glucose levels. J Biol Chem 269, 21234–8 (1994).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNjkvMzMvMjEyMzQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8wNS8yMDIxLjEyLjAzLjIxMjY3MjQwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

14. 14.Bianchi, M. & Magnani, M. Hexokinase mutations that produce nonspherocytic hemolytic anemia. Blood Cells Mol Dis 21, 2–8 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/bcmd.1995.0002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7655856&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RP69700001&link_type=ISI) 

15. 15.Pasquali, L. et al. Pancreatic islet enhancer clusters enriched in type 2 diabetes risk-associated variants. Nat Genet 46, 136–143 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2870&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24413736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

16. 16.Chiou, J. et al. Single-cell chromatin accessibility identifies pancreatic islet cell type-and state-specific regulatory programs of diabetes risk. Nat Genet 53, 455–466 (2021).
    
    

17. 17.Greenwald, W.W. et al. Pancreatic islet chromatin accessibility and conformation reveals distal enhancer networks of type 2 diabetes risk. Nat Commun 10, 2078 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09975-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31064983&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

18. 18.Weng, C. et al. Single-cell lineage analysis reveals extensive multimodal transcriptional control during directed beta-cell differentiation. Nat Metab 2, 1443–1458 (2020).
    
    

19. 19.Balboa, D. et al. Functional, metabolic and transcriptional maturation of stem cell derived beta cells. bioRxiv, 2021.03.31.437748 (2021).
    
    

20. 20.Geusz, R.J. et al. Pancreatic progenitor epigenome maps prioritize type 2 diabetes risk genes with roles in development. Elife 10(2021).
    
    

21. 21.Cebola, I. et al. TEAD and YAP regulate the enhancer network of human embryonic pancreatic progenitors. Nat Cell Biol 17, 615–626 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncb3160&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25915126&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

22. 22.Dhawan, S. et al. DNA methylation directs functional maturation of pancreatic beta cells. J Clin Invest 125, 2851–60 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI79956&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26098213&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

23. 23.Giri, D. et al. Novel FOXA2 mutation causes Hyperinsulinism, Hypopituitarism with Craniofacial and Endoderm-derived organ abnormalities. Hum Mol Genet 26, 4315–4326 (2017).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28973288&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

24. 24.Papizan, J.B. et al. Nkx2.2 repressor complex regulates islet beta-cell specification and prevents beta-to-alpha-cell reprogramming. Genes Dev 25, 2291–305 (2011).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6MTA6IjI1LzIxLzIyOTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8wNS8yMDIxLjEyLjAzLjIxMjY3MjQwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

25. 25.Keller, M.P. et al. The Transcription Factor Nfatc2 Regulates beta-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets. PLoS Genet 12, e1006466 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1006466&link_type=DOI) 

26. 26.Pullen, T.J. et al. Identification of genes selectively disallowed in the pancreatic islet. Islets 2, 89–95 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/isl.2.2.11025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21099300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000289687400005&link_type=ISI) 

27. 27.Lemaire, K., Thorrez, L. & Schuit, F. Disallowed and Allowed Gene Expression: Two Faces of Mature Islet Beta Cells. Annu Rev Nutr 36, 45–71 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-nutr-071715-050808&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27146011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

28. 28.Pinney, S.E. et al. Dominant form of congenital hyperinsulinism maps to HK1 region on 10q. Horm Res Paediatr 80, 18–27 (2013).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23859901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

29. 29.Henquin, J.C. et al. Congenital hyperinsulinism caused by hexokinase I expression or glucokinase-activating mutation in a subset of beta-cells. Diabetes 62, 1689–96 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjIvNS8xNjg5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMDUvMjAyMS4xMi4wMy4yMTI2NzI0MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

30. 30.Otonkoski, T. et al. Physical exercise-induced hypoglycemia caused by failed silencing of monocarboxylate transporter 1 in pancreatic beta cells. Am J Hum Genet 81, 467–74 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/520960&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701893&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249128200004&link_type=ISI) 

31. 31.Almeida, A.M. et al. Hypomorphic promoter mutation in PIGM causes inherited glycosylphosphatidylinositol deficiency. Nat Med 12, 846–51 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1410&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16767100&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238862800074&link_type=ISI) 

32. 32.Thornton, P.S. et al. Recommendations from the Pediatric Endocrine Society for Evaluation and Management of Persistent Hypoglycemia in Neonates, Infants, and Children. J Pediatr 167, 238–45 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpeds.2015.03.057&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25957977&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

33. 33.Ellard, S. et al. Improved genetic testing for monogenic diabetes using targeted next-generation sequencing. Diabetologia 56, 1958–63 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-013-2962-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23771172&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323059800014&link_type=ISI) 

34. 34.Campbell-Thompson, M. et al. Network for Pancreatic Organ Donors with Diabetes (nPOD): developing a tissue biobank for type 1 diabetes. Diabetes Metab Res Rev 28, 608–17 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/dmrr.2316&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22585677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

35. 35.Laver, T.W. et al. SavvyCNV: genome-wide CNV calling from off-target reads. bioRxiv, 617605 (2019).
    
    

36. 36.Han, B. et al. Enhanced Islet Cell Nucleomegaly Defines Diffuse Congenital Hyperinsulinism in Infancy but Not Other Forms of the Disease. Am J Clin Pathol 145, 757–68 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ajcp/aqw075&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27334808&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

37. 37.Bankhead, P. et al. QuPath: Open source software for digital pathology image analysis. Sci Rep 7, 16878 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-17204-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29203879&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

38. 38.Langmead, B. & Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat Methods 9, 357–9 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1923&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22388286&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302218500017&link_type=ISI) 

39. 39.Abdennur, N. & Mirny, L.A. Cooler: scalable storage for Hi-C data and other genomically labeled arrays. Bioinformatics 36, 311–316 (2020).
    
    

40. 40.Owens, N.D.L. et al. Measuring Absolute RNA Copy Numbers at High Temporal Resolution Reveals Transcriptome Kinetics in Development. Cell Rep 14, 632–647 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2015.12.050&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26774488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

41. 41.Fang, Z. et al. Single-Cell Heterogeneity Analysis and CRISPR Screen Identify Key beta-Cell-Specific Disease Genes. Cell Rep 26, 3132–3144 e7 (2019).
    
    

42. 42.Fornes, O. et al. JASPAR 2020: update of the open-access database of transcription factor binding profiles. Nucleic Acids Res 48, D87–D92 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkaa516&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31701148&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 

43. 43.Song, J.Y. et al. HIC2, a new transcription activator of SIRT1. FEBS Lett 593, 1763–1776 (2019).
    
    

44. 44.Luu, L. et al. The loss of Sirt1 in mouse pancreatic beta cells impairs insulin secretion by disrupting glucose sensing. Diabetologia 56, 2010–20 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-013-2946-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23783352&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F05%2F2021.12.03.21267240.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323059800020&link_type=ISI) 

45. 45.Lee, J.H., Lee, J.H. & Rane, S.G. TGF-beta Signaling in Pancreatic Islet beta Cell Development and Function. Endocrinology 162(2021).