Maturity Onset Diabetes of the Young (MODY) is a collection of monogenic forms of early-onset familial diabetes resulting from β-cell dysfunction. To date, mutations in 15 genes involved in β-cell development and function have been associated with MODY(1). Monogenic diabetes cases account for approximately 1%-5% of the total diabetic patient population(2, 3). However, as the phenotype of MODY overlaps with other forms of diabetes, many patients are misdiagnosed(4, 5). Data from monogenic diabetes registries also suggests that in the US, ∼70% of the registered cases (2680/3800) still do not have a known genetic cause, and similarly many patients in the UK (∼50%) with monogenic diabetes have MODY-causing mutations that remain to be determined(6, 7). Additionally, some MODY-associated mutations are only reported in single families, and thus require genetic and mechanistic confirmation (8, 9). We recently identified a mutation in KCNK16, the gene encoding TALK-1 channels (p. TALK-1 L114P), which co-segregates with MODY in a four-generation family(10). Importantly, this mutation causes a gain-of-function (GOF) in TALK-1 channel activity. As KCNK16 is the most islet-restricted and abundant β-cell K+ channel transcript, the TALK-1 L114P mutation is likely to perturb β-cell function and cause MODY(11, 12). This was further strengthened by a recent report of a Japanese family with the identical KCNK16 mutation (p. TALK-1 L114P) segregating with MODY(13). Thus, it is important to determine the mechanistic underpinnings of how this specific mutation in TALK-1 impacts islet function and results in glucose dyshomeostasis.

TALK-1 is a key regulator of β-cell membrane potential (Vm), glucose-stimulated Ca2+ influx and insulin secretion(14). A non-synonymous GOF polymorphism in KCNK16 (rs1535500; p. TALK-1 A277E) is also associated with an increased risk for type-2 diabetes (T2D)(15, 16). Moreover, we recently determined that the MODY-associated TALK-1 L114P mutation results in a significant GOF in TALK-1 K+ flux(10). Thus, there is strong genetic evidence that alterations in TALK-1 function (e.g., TALK-1 L114P and TALK-1 A277E) result in diabetic phenotypes. Interestingly, when heterologously expressed in β-cells, TALK-1 L114P almost completely inhibited glucose-stimulated Vm depolarization and Ca2+ influx in most β-cells. While these TALK-1 L114P-mediated drastic changes in electrical activity and Ca2+ would be predicted to fully inhibit glucose-stimulated insulin secretion (GSIS), they only resulted in partial blunting of β-cell GSIS. However, this is consistent with TALK-1 L114P MODY patients that required low-dose insulin therapy which also suggested that this mutation does not fully suppress GSIS. This differs from previous studies showing that inhibition of β-cell Ca2+ influx with either K+ channel pharmacological activation or GOF mutations (e.g., KATP R201H) results in complete inhibition of GSIS(17). Therefore, it is critical to determine if endogenous TALK-1 L114P expression also shows complete inhibition of β-cell Ca2+ influx and how this impacts GSIS. Another potential mechanism for the modest impact of TALK-1 L114P on GSIS could be due to its function in other islet cells. For example, TALK-1 is also expressed in δ-cells where its activity limits somatostatin secretion(18).

Because somatostatin exerts inhibitory tone on islet β- and α-cells, TALK-1 L114P mediated reductions in δ-cell somatostatin secretion would be predicted to increase both glucagon and insulin secretion. Due to the β-cell intrinsic role of TALK-1 L114P channels, the effect of somatostatin on GSIS might be limited, however it remains to be determined if the fasting hyperglycemia observed in KCNK16-MODY patients is due in part to elevated glucagon secretion.

The nature and severity of MODY phenotypes is dictated by how specific gene mutations affect β-cell function. KATP is the only ion channel besides TALK-1 to be linked to MODY; this is due to mutations in genes encoding the KATP channel complex (KCNJ11 and ABCC8) or mutations in genes affecting ATP synthesis (e.g., glucokinase)(19, 20). Interestingly, MODY-associated mutations in KCNJ11 and ABCC8 have also been found to cause other diabetic phenotypes including permanent or transient neonatal diabetes, and late-onset diabetes(21, 22). Although the two families with TALK-1 L114P primarily exhibit a MODY phenotype, it remains to be determined if KCNK16 mutations are associated with other diabetic phenotypes besides MODY. Indeed, it was originally predicted that GOF TALK-1 mutations could cause neonatal diabetes, however, the voltage dependence of activation of TALK-1 makes this less likely that GOF TALK-1 channels would be active at resting Vm (23). Thus, TALK-1 GOF would limit but not abolish GSIS during β-cell glucose-mediated Vm depolarization and cause MODY. As KATP channels show larger unitary conductance than TALK-1, and are active at resting Vm, severe GOF in KATP channels strongly suppress β-cell Ca2+ entry and GSIS (24, 25). The biophysical differences in KATP and TALK-1 especially in the context of nucleotide-gating predict differences in their corresponding MODY phenotypes. The phenotypes of KCNJ11- and ABCC8-MODY are well characterized; however, with only two known KCNK16-MODY families, it remains to be elucidated how endogenous TALK-1 L114P disrupts islet electrical activity and glucose homeostasis.

Here we developed a mouse model harboring the MODY-associated Kcnk16 variant p. TALK-1 L114P to investigate the impact of this mutation on glucose homeostasis and confirm its association with MODY. Interestingly, we observe increased neonatal lethality in heterozygous and homozygous Kcnk16 L114P mice in the C57BL/6J (B6) and the B6;CD-1(ICR) mixed genetic backgrounds, respectively, likely due to severe hyperglycemia and lack of GSIS. Whereas in adult mice, TALK-1-L114P blunts glucose-stimulated β-cell electrical activity, Ca2+ handling, and GSIS, thus significantly impairing glucose tolerance. These data strongly suggest that alterations in TALK-1 activity can disrupt islet hormone secretion and glucose homeostasis. Importantly, this study further confirms that the TALK-1-L114P mutation results in a MODY phenotype, but additionally predicts that GOF TALK-1 mutations cause transient neonatal diabetes.


Kcnk16 L114P neonates exhibit loss of glucose stimulated Ca2+ entry and insulin secretion leading to transient neonatal hyperglycemia and death.

To confirm the association of TALK-1 L114P in causing MODY, a mouse model harboring the Kcnk16 mutation was developed in the C57BL/6J background using CRISPR/spCas9 (B6 Kcnk16 L114P; Figure 1A and supplemental figures 1A and 1B). Surprisingly, heterozygous Kcnk16 L114P (L/P) mice exhibited neonatal lethality as indicated by extremely low number of Kcnk16 L114P (L/P) mice at weaning (Figure 1B). To increase the likelihood of survival of neonates, B6 Kcnk16 L114P (L/P) mice were crossed with the outbred CD-1(ICR) mouse strain resulting in a progeny in the mixed (50:50) background (B6;CD-1 Kcnk16 L114P; Figure 1A). Intriguingly, in the mixed background, neonatal lethality was observed in the homozygous B6;CD-1 Kcnk16 L114P (P/P) mice on ∼post-natal day 4 (P4), but not in the heterozygous B6;CD-1 Kcnk16 L114P (L/P) mice (Supplemental figure 1C). Lethality was likely independent of growth defects as body weight did not differ between genotypes (Supplemental figure 1D). To test if lethality occurs from GOF TALK-1 L114P mediated defect in neonatal islet function, we assessed glycemic control and insulin secretion on P4. Heterozygous and homozygous Kcnk16 L114P neonates showed severe hyperglycemia and a concurrent reduction in plasma insulin levels in a gene-dosage dependent manner compared to the control littermates (WT) on P4 (Figures 1C and 1D, and supplemental figure 1E). Additionally, TALK-1 L114P mutation blunted GSIS in islets from heterozygous neonates and decreased it further in islets from homozygous neonates (Figures 1E and 1F). These changes in GSIS and glucose homeostasis are likely due to TALK-1 L114P mediated alteration in islet function and not islet mass as pancreas weight was unchanged (Figure 1G and supplemental figure 1F). Furthermore, neonatal islet composition was unaltered except for a modest increase in the α-cell (glucagon+) area/islet in the Kcnk16 L114P(L/P) pancreas (Figures 1H-J). The loss of GSIS can be explained by a complete lack of glucose-stimulated Ca2+ entry in islets from Kcnk16 L114P (L/P and P/P) neonates (Figures 1K-N). However, these islets exhibited KCl-mediated β-cell Vm depolarization-induced Ca2+ entry; this indicates that TALK-1 L114P channels hyperpolarize β-cell Vm, thereby overriding the effect of KATP closure due to glucose metabolism and lead to insufficient GSIS and hyperglycemia. Importantly, insulin treatment (Glargine Lantus; 0.2U/Kg/day, once-daily subcutaneous injection) was able to extend the lifespan of Kcnk16 L114P (P/P) neonates suggesting that lethality results from hyperglycemia due to inadequate insulin secretion (Figure 1O).

Kcnk16 L114P neonates exhibit loss of glucose stimulated Ca2+ entry and insulin secretion leading to transient neonatal hyperglycemia and death.

A. Schematic of the Kcnk16 L114P mouse line generation in the C57BL/6J background and the C57BL/6J:CD-1(ICR) mixed background. B. χ2 analysis of the F1 progeny from C57BL/6J and heterozygous Kcnk16 L114P (L/P) crossings genotyped at weaning on post-natal day 21 (P21). C. Blood glucose measurements of male wildtype (WT; black; N=8), heterozygous Kcnk16 L114P (L/P; green; N=16), and homozygous Kcnk16 L114P (P/P; blue; N=9) mice on P4. D. Plasma insulin measurements performed on P4 in WT (N=4), Kcnk16 L114P (L/P; N=4), and Kcnk16 L114P (P/P; N=5) neonates. E. and F. In-vitro glucose-stimulated insulin secretion from P4 mouse islets stimulated with 2.8 mM glucose (G) or 20 mM G, respectively (WT; N=4, Kcnk16 L114P (L/P); N=8, and Kcnk16 L114P (P/P); N=6). G. Pancreas weight/ body weight measurements of P4 male mice (WT; N=6, Kcnk16 L114P (L/P); N=14, and Kcnk16 L114P (P/P); N=9). H. Representative immunostaining images of pancreas sections from P0 WT and Kcnk16 L114P (L/P) mice (N=3 mice/genotype), stained against insulin (INS; green), glucagon (GCG; red), somatostatin (SST; magenta), and Hoechst (blue); scale bar=100 μm. I. and J. Average islet area and area fraction of hormone staining per islet quantified using Fiji ImageJ software in P0 mouse pancreas sections. K. Representative glucose stimulated Ca2+ influx traces from P4 mouse islets sequentially stimulated with 2 mM glucose (G), 7 mM G, 20 mM G, and 20 mM G with 30 mM KCl (WT; N=6, Kcnk16 L114P (L/P); N=9, and Kcnk16 L114P (P/P); N=7). L-N. Average area under the curve (AUC) analysis of normalized Ca2+ during 7 mM G, 20 mM G, and 20 mM G+30 mM KCl stimulations. O. Survival curve for Kcnk16 L114P (P/P) mice treated with (N=6) or without (N=6) insulin (Lantus insulin glargine; 0.2 U/kg/day) starting at P0 up until death or weaning. Data are presented as mean±SEM. Data were analyzed using student’s t-test, one-way ANOVA, and two-way ANOVA as appropriate.

Adult Kcnk16 L114P mice exhibit fasting hyperglycemia and glucose intolerance

Because neonatal Kcnk16 L114P (L/P) mice showed transient hyperglycemia, we next assessed if TALK-1 L114P impairs glucose homeostasis in young-adulthood. Male Kcnk16 L114P (L/P) mice developed glucose intolerance as early as 11-weeks of age in the B6;CD-1 background and 8-weeks of age in the B6 background compared to their respective control littermates (data not shown). This defect in glucose homeostasis is maintained with age in the male Kcnk16 L114P (L/P) mice (Figures 2A-C and supplemental figures 2A-C). Furthermore, male Kcnk16 L114P (L/P) mice also exhibit fasting hyperglycemia indicating a likely non-β-cell-specific islet defect in this mouse model (Figure 2B). There were no body weight differences in males (Figure 2D). Interestingly, the female Kcnk16 L114P (L/P) mice develop moderate glucose intolerance compared to the controls (WT) in both B6 and B6;CD-1 genetic backgrounds (Figures 2E-G and supplemental figures 2D-F). Moreover, TALK-1 L114P mutation leads to greater body weight gain in females compared to the WT controls (Figure 2H); this increase is likely contributed by lean mass which might render a protective effect on glucose tolerance in Kcnk16 L114P (L/P) females compared to their male counterparts (Supplemental figures 3A-F). The T2D-associated SNP in KCNK16 (rs1535500) which causes GOF in TALK-1 activity, has shown a significant association with increased total cholesterol in a Han Chinese population (26). However, an increase in TALK-1 activity due to the Kcnk16 L114P mutation did not affect plasma- and liver-triglyceride and cholesterol levels in mice (Supplemental figures 3G-J).

Adult Kcnk16 L114P mice exhibit fasting hyperglycemia and glucose intolerance.

A. and B. Intraperitoneal glucose tolerance test (i.p. GTT) performed in 15-week-old and 18-week-old male mice following a 4-hour fast in response to 2mg/g glucose injection (WT; black; N=8-10 and Kcnk16 L114P (L/P); green; N=9-10). C. Average AUC of the 2-hr GTT excursion profiles from ages 7 weeks up to 18 weeks in male mice. D. Weekly body weight measurements of male WT; N=5 mice and Kcnk16 L114P (L/P); N=5 mice E. and F. I.P. GTT performed in 15-week-old and 18-week-old female mice following a 4-hour fast in response to 2mg/g glucose injection WT; N=9-11 and Kcnk16 L114P (L/P); N=10-11. G. Average AUC of the 2-hr GTT excursion profiles from ages 7 weeks up to 18 weeks in female mice. H. Body weight measurements of female WT (black; N=4) and Kcnk16 L114P (L/P; blue; N=5) mice. Data are presented as mean±SEM. Data were analyzed using student’s t-test.

Adult Kcnk16 L114P mice show disrupted islet hormone secretion and islet composition.

We then tested if the TALK-1-L114P-mediated impairment in glucose homeostasis arises from defective islet function. Both male and female Kcnk16 L114P (L/P) mice exhibited reduced plasma insulin levels compared to controls at 15 min and 30 min following an intraperitoneal glucose injection after a 4-hour fast (Figures 3A and 3B). Insulin sensitivity was not altered in the Kcnk16 L114P (L/P) mice indicating that the reduction in glucose tolerance is primarily due to an islet secretion defect (Figure 3C). This was confirmed by in vitro GSIS and glucose-inhibition of glucagon secretion (GIGS) measurements (Figures 3D and 3E). Islets from Kcnk16 L114P (L/P) mice exhibited reduced ability to secrete insulin in response to glucose and showed elevated glucagon secretion under low glucose (2 mM) and euglycemic (7 mM) conditions (Figures 3D and 3E). This indicates that the fasting hyperglycemia observed in male Kcnk16 L114P (L/P) mice (Figure 2B) is likely a result of the observed increase in glucagon secretion under fasting conditions. We next investigated if impaired islet hormone secretion results from changes in islet composition. Immunostaining analyses revealed an increase in glucagon-positive area/ islet and a concurrent modest reduction in insulin-positive area/ islet in Kcnk16 L114P (L/P) pancreata (Figures 3F-J). Together, these results indicate that TALK-1 L114P mutation leads to disruptions in both islet composition and hormone secretion giving rise to a MODY-like phenotype in adult mice.

Adult Kcnk16 L114P mice show disrupted islet hormone secretion and islet composition.

A. and B. Plasma insulin levels in male (A) and female (B) WT and Kcnk16 L114P(L/P) mice following a 4-hour fast at the indicated time points before and after a 2mg/g glucose injection. c. Glucose (%) change in response to an i.p. human insulin injection (0.75 UI/kg body weight) was measured as an indicator of insulin sensitivity in WT and Kcnk16 L114P male mice after a 4-hour fast. D. and E. In-vitro insulin secretion (N=5 mice/genotype) and glucagon secretion (N=3 mice/genotype) from male mice at the specified glucose concentrations F. Representative immunostaining images of pancreas sections from WT and Kcnk16 L114P(L/P) male mice (N=3/genotype) stained against insulin (INS; green), glucagon (GCG; red), somatostatin (SST; magenta), and Hoechst (blue); scale bar = 100 μm. G-J. Average islet area and area of hormone staining/ islet for β-cells (INS), α-cells (GCG), and δ-cells (SST). Data are presented as mean±SEM. Data were analyzed using student’s t-test or two-way ANOVA as appropriate.

Kcnk16 L114P blunts β-cell glucose-stimulated electrical excitability and increases whole-cell K2P currents.

We used patch-clamp electrophysiology to test if the reduction in GSIS results from decreased β-cell electrical excitability due to a GOF in TALK-1 activity in Kcnk16 L114P mouse islets. Perforated patch-clamp recordings revealed blunting of glucose-stimulated β-cell Vm depolarization and loss of action potential firing in islets from Kcnk16 L114P (L/P) mice (Figures 4A and 4B). The lack of islet Vm depolarization and action potentials was in part a result of increased whole-cell K2P currents at depolarized Vm (Figures 4C and 4D). However, TALK-1 L114P (L/P and P/P) channels did not show a massive GOF at hyperpolarized β-cell membrane potentials under low and high glucose where these K+ currents are typically active (Supplementary figures 4A and 4B). Notably, the heterologously expressed human TALK-1 L114P channels exhibited a drastic GOF (>300-fold)(10); this suggests the likelihood of unknown endogenous regulator(s) of TALK-1 in β-cells. It is also likely that recording conditions lead to poor isolation of the endogenous TALK-1 L114P currents. Thus, K2P currents were assessed by overexpressing mouse TALK-1 L114P channels in HEK293T cells under identical recording conditions as the human TALK-1 L114P study. Indeed, overexpression of the mouse mutant channel results in large K+ currents indicative of a similar GOF in mouse TALK-1 current as the human TALK-1 current due to the L114P mutation (Figures 4E and 4F).

Kcnk16 L114P blunts β-cell glucose-stimulated electrical excitability and increases whole-cell K2P currents.

A. Representative perforated patch-clamp Vm recordings in response to 2 mM G and 11 mM G in islets from WT and Kcnk16 L114P(L/P) mice (N=6-10 islets from 5 mice/genotype). B. Average resting islet Vm under 2 mM G and plateau islet Vm at 11 mM G. C. Representative whole-cell K2P current density (pA/pF) recorded using a voltage ramp (−120 mV to +60 mV) at 11 mM G in β-cells from WT and Kcnk16 L114P(L/P) mice. D. Average current density (pA/pF) at the specified membrane voltages during the voltage ramp recordings shown in panel C. (N=9-13 cells/genotype). E. Representative whole-cell K2P current (pA) recorded using a voltage ramp (−120 mV to +60 mV) in 11 mM G in HEK293T cells expressing either Kcnk16 WT or Kcnk16 L114P(L/P). F. Average current (pA) at the specified membrane voltages during the voltage ramp recordings shown in panel E. (N=9-11 cells/ group). Data are presented as mean±SEM. Data were analyzed using two-way ANOVA.

Kcnk16 L114P reduces glucose- and tolbutamide-stimulated Ca2+ entry and augments IP3-induced [Ca2+]ER release.

β-cell Ca2+ entry was monitored in response to increasing concentrations of glucose (G) in islets from WT and Kcnk16 L114P (L/P) mice to assess if TALK-1 L114P-mediated Vm hyperpolarization reduces glucose-stimulated Ca2+ influx (GSCI). TALK-1 L114P blunted islet GSCI in mice on both B6 and B6;CD-1 genetic backgrounds compared to WT controls (Figures 5A-F and supplemental figures 5A and 5B). However, islets from male Kcnk16 L114P (L/P) mice showed a larger reduction in normalized Ca2+ peak and normalized area under the curve (AUC) in response to glucose compared to their controls than islets from female Kcnk16 L114P (L/P) mice compared to their controls (Figures 5A-F). TALK-1 L114P channels did not affect the basal Ca2+ levels at 2 mM G; however, the activity of mutant channels resulted in increased glucose-stimulated phase 0 response compared to controls (Supplemental figures 6A and 6B). This suggests a reduction in endoplasmic reticulum (ER) Ca2+ ([Ca2+]ER) storage in Kcnk16 L114P islets. Furthermore, Kcnk16 L114P islets also showed a complete lack of glucose-stimulated Ca2+ oscillations monitored at 9 mM G (Supplemental figures 6C and 6D). We then assessed if the reduction in GSCI results from TALK-1 L114P overactivity on the β-cell plasma membrane by monitoring tolbutamide-stimulated Vm depolarization-induced Ca2+ entry. TALK-1 L114P mediated β-cell Vm hyperpolarization overrides the Vm depolarization caused by KATP closure by tolbutamide. Thus, Kcnk16 L114P islets do not exhibit Ca2+ entry in response to tolbutamide. Interestingly, islets from TALK-1 L114P mice showed equivalent KCl-stimulated Ca2+ entry compared to control islets; this suggests that the loss of GSCI in the TALK-1 L114P islets occurs due to enhanced K+ ion flux through GOF in TALK-1 activity (Figures 5G and 5H). Our previous studies showed that TALK-1 channels are expressed on the ER membrane where they function to provide a K+ countercurrent for Ca2+ release from the ER lumen(27). Thus, we tested if TALK-1 L114P regulates [Ca2+]ER homeostasis by stimulating Gq-signaling using acetylcholine which results in IP3 generation, subsequently causing IP3-induced [Ca2+]ER release. Indeed, IP3-induced [Ca2+]ER release was enhanced in islets from Kcnk16 L114P mice compared to controls suggesting that GOF TALK-1 channels are localized to the ER membrane where they facilitate increased Ca2+ release from the ER lumen (Figures 5I and 5J).

Kcnk16 L114P reduces glucose- and tolbutamide-stimulated islet Ca2+ entry and augments IP3-induced islet [Ca2+]ER release.

A. and B. Representative [Ca2+]c traces in islets from male (A) and female (B) WT and Kcnk16 L114P(L/P) mice in response to 2 mM G, 10 mM G, and 20 mM G. C-F. Average normalized Ca2+ peak (C. and E.) and total AUC (D. and F.) in response to the indicated glucose concentrations in islets from male and female WT and Kcnk16 L114P(L/P) mice (N=3-4/genotype). G. Representative [Ca2+]c traces in islets from WT and Kcnk16 L114P(L/P) male mice in response to 100 μM tolbutamide followed by 30 mM KCl stimulation. H. Average AUC of normalized [Ca2+]c during 100 μM tolbutamide or 100 μM tolbutamide with 30 mM KCl stimulation (N=3 mice/genotype). I. Representative [Ca2+]c traces in response to 100 μM acetylcholine in the absence of extracellular Ca2+ in islets from male WT and Kcnk16 L114P(L/P) mice. J. Average AUC of normalized [Ca2+]c following 100 μM acetylcholine stimulated [Ca2+]ER release (N=876 cells; WT, N=513 cells; Kcnk16 L114P (L/P)). Data are presented as mean±SEM. Data were analyzed using student’s t-test.

TALK-1 L114P islets exhibit altered expression of genes involved in β-cell identity, function, ion channel activity, hormone activity, inflammatory signaling, and extracellular matrix interaction pathways.

Despite the loss of glucose-stimulated islet electrical activity and Ca2+ entry, Kcnk16 L114P mice do not exhibit a drastic reduction in GSIS and overt diabetes. Therefore, compensatory mechanisms such as the amplification pathway for insulin secretion might be altered in Kcnk16 L114P islets (28). Moreover, prolonged hyperglycemia (although modest in this model), would also be predicted to result in long-term gene expression changes in Kcnk16 L114P islets that alter β-cell function. Thus, gene expression differences between WT and Kcnk16 L114P mouse islets were quantified with bulk-RNA sequencing and validated with qRT-PCR (Figures 6 A-D). We observed increased expression of many genes which regulate Ca2+-independent potentiation of insulin secretion such as Adcy5, Creb5, Adcyap1, and Adcyap1r1; genes that may promote Ca2+-independent secretion from Kcnk16 L114P islets(29, 30). Additionally, broad changes in ion channel activity genes (e.g., Cacna1g, Cacng8, Scn5a, Pkd1l1, Kcnk2, Gabrg1, and Fxyd3) may contribute to altered Ca2+ handling in Kcnk16 L114P islets. Reduced Ca2+ entry in Kcnk16 L114P islets results in decreased expression of genes previously shown to be elevated in chronically depolarized Abcc8−/− β-cells such as Serpina7, Asb11, Sall1, and Aldh1a3 (Figure 6A) (31). Finally, several stress-associated, fibrosis-related, and inflammatory signaling pathway genes are upregulated in Kcnk16 L114P islets likely due to prolonged hyperglycemic conditions. These include dedifferentiation markers Sox4, Sox6, Sox9, Hk2, Vim, and Cd36, extracellular matrix-interaction pathway genes Col1a1, Col1a2, Col3a1, Col14a1, Col6a1, and Dcn, and inflammatory signaling genes Cxcl1, Ccl2, Ccl11, Ccl22, Tgfb2, Il33, and Il6 (3237). The data show that Kcnk16 L114P islets exhibit gene expression differences in many pathways important for normal function. Taken together, this suggests that perturbed Kcnk16 L114P islet intrinsic (e.g., Ca2+ handling) and extrinsic (e.g., hyperglycemia) pathways result in direct as well as indirect disruption of islet function.

Kcnk16 L114P islets exhibit altered expression of genes involved in β-cell identity and function, ion channel activity, hormone activity, inflammatory signaling, and extracellular matrix interaction pathways.

A. Heatmap of a selected gene subsets showing differential gene expression in WT and Kcnk16 L114P islets. Normalized expression levels were scaled and centered by rows. B. Volcano plot displays genes differentially expressed between WT and Kcnk16 L114P samples. Differentially expressed genes are defined by FDR <0.05 and log2FC (≥1). C. Dotplot represents the top 10 most significantly (FDR <0.05) altered Gene Ontology (Molecular Function). GeneRatio represent (count of enriched genes) / (count of genes in the GO term). The color represents FDR adjusted p-values and the size of the dot represents the number of genes that are significant from the experimental dataset. D. qRT-PCR validation of the gene expression differences in WT and Kcnk16 L114P samples for the selected genes observed through bulk-sequencing.


GOF in TALK-1 activity is associated with diabetic phenotypes, suggesting a causal role for overactive TALK-1 in islet dysfunction and diabetes progression. Utilizing a GOF model of TALK-1 (L114P) this study uncovered that disrupted β-cell and α-cell function resulted in glucose intolerance in adolescent mice confirming the association of the TALK-1 L114P mutation with MODY-like diabetes. Importantly, our data also revealed that TALK-1 L114P mutation can cause transient neonatal diabetes. This finding suggests that transient neonatal diabetic patients with unknown genetic diagnosis should be screened for mutations in KCNK16. Moreover, our data provides further genetic evidence that TALK-1 is a potentially novel therapeutic target for diabetes treatment. Specifically, in mice heterozygous for Kcnk16 L114P mutation we observe neonatal hyperglycemia due to blunted GSIS, which can additionally result in neonatal death in mice homozygous for the mutant allele of Kcnk16. In young-adulthood, TALK-1 L114P causes glucose intolerance due to a reduction in GSIS mediated by enhanced β-cell Vm hyperpolarization and reduced glucose-stimulated Ca2+ entry. In addition to the β-cell intrinsic defect, the TALK-1 L114P mutation led to an increase in α-cell area fraction in islets and elevation of glucagon secretion under fasting conditions. Together these data highlight the crucial role of TALK-1 in β-cell function and glucose homeostasis and raise the possibility of TALK-1 inhibition as a druggable target for not only KCNK16-associated MODY but possibly for other forms of diabetes. Unlike the only other MODY-associated K+ channelopathy (KCNJ11 and ABCC8 i.e., KATP channel mutations), TALK-1 is not inhibited by sulfonylureas (14).

β-cell maturation and glucose-responsiveness of neonatal islets rapidly develops after birth and GSIS is required for efficient glucose uptake, which contributes to normal growth. Following birth, a shift to intermittent feeding and elevated plasma glucose requires an increase in β-cell insulin secretion for efficient nutrient absorption, as well as a suppression of insulin release during fasting to avoid hypoglycemia (38). Glucose-sensitivity in mouse β-cell Ca2+ handling and insulin secretion develops over the first 4 postnatal days due to changes in β-cell metabolism, KATP surface localization, and Ca2+-dependent secretory machinery (39). Thus, mice expressing ATP-insensitive KATP channels die shortly after birth due to hypoinsulinemia, severe hyperglycemia, and ketoacidosis (40). Kcnk16 L114P mice show a similar phenotype, with neonatal islets showing a complete loss of GSCI, a drastic reduction in GSIS, and severe hyperglycemia by P5. Interestingly, TALK-1 L114P-mediated susceptibility to neonatal lethality was dependent on genetic diversity in mouse strains. Heterozygous Kcnk16 L114P caused almost complete neonatal lethality in the B6 background but only hyperglycemia in the B6:CD-1 background, whereas homozygous Kcnk16 L114P led to neonatal lethality in the B6:CD-1 mice. Although the neonatal glycemic data from the KCNK16 L114P families does not exist, it may be that neonatal lethality is not observed in affected individuals. However, it is interesting to note that all currently identified individuals carrying the KCNK16 L114P mutation are females (N=8) (10, 13). The importance of hypoinsulinemia in causing neonatal diabetes and lethality in TALK-1 L114P mice was confirmed by insulin treatment which extended their lifespan. This resembles other monogenic forms of neonatal diabetes that require exogenous insulin treatment for survival (41). Indeed, it has been well established that insulin signaling is required for neonatal survival; for example, a similar neonatal lethality phenotype was observed in mice without insulin receptors (Insr-/-) where death results from hyperglycemia and diabetic ketoacidosis by P3 (42). Future studies are required to determine if TALK-1 GOF mutations result in human transient neonatal diabetes and are possibly followed by a MODY phenotype. Other monogenic diabetes mutations such as in genes encoding KATP channels result in transient neonatal diabetes and diabetes remits later in life (43).

In adolescence, the timeline of development of MODY-like diabetes in the Kcnk16 L114P (L/P) model is consistent with data from MODY patients. Similar to the timeline of disease progression in other MODY mouse models, Kcnk16 L114P mice developed glucose intolerance during adolescence (∼8 weeks in the B6 strain) (44, 45). The onset and severity of glucose intolerance in Kcnk16 L114P mice also recapitulates the data from probands from the Japanese and Australian KCNK16 L114P MODY families, who were diagnosed at 11 and 15 years of age and displayed an abnormal oral glucose tolerance test (blood glucose:19 and 19.6 mmol/L, 2 hours after 75 g oral glucose bolus, respectively). The disease severity is more prominent in the inbred B6 strain compared to the mixed B6:CD-1 strain likely owing to the outbred characteristics and genetic diversity of the CD-1 strain, which recapitulates the diverse nature of MODY manifestation observed in human patients. Intriguingly, the patients from the Japanese KCNK16 L114P MODY family showed a greater insulin requirement compared to the patients from Australian KCNK16 L114P MODY family (10, 13). Interestingly, male Kcnk16 L114P (L/P) mice exhibit more severe impairment in glucose homeostasis compared to female Kcnk16 L114P (L/P) mice. A variety of factors could be at play here which could explain the observed sex-differences, including the female hormone 17β-estradiol (E2) which is critical for protection against glucolipotoxicity and oxidative stress(46). Furthermore, for a given insulin action, comparing older men and women, women have higher insulin levels (47). Moreover, GSCI in islets from female Kcnk16 L114P mice is significantly greater than islets from male Kcnk16 L114P mice, which would be predicted to lead to greater insulin secretion and lesser impairment in glucose tolerance in female Kcnk16 L114P mice. Sexual dimorphism observed in GSCI could be mediated by differences in Ca2+ handling or differences in TALK-1 function. Whether sex-differences exist in humans with TALK-1 L114P mutation remains to be determined as all affected individuals in both KCNK16-MODY (p. TALK-1 L114P) families were females. Neonatal lethality was more penetrant in male than female Kcnk16 L114P (P/P) mice, thus, it will also be important to determine if KCNK16-MODY patients show neonatal phenotypes and display sexual dimorphism.

Similar to the proband from the human KCNK16 L114P MODY family who showed elevated fasting blood glucose (∼7mmol/L), Kcnk16 L114P mice also exhibit fasting hyperglycemia. While this was predicted to be due to decreased insulin secretion, Kcnk16 L114P islets showed equivalent insulin secretion under euglycemic conditions. However, these islets exhibited a significant elevation in fasting glucagon secretion which likely contributes to the fasting hyperglycemia. If TALK-1 L114P channels were expressed in α-cells, it would result in inhibition of α-cell Ca2+ entry and glucagon secretion which supports previously described lack of TALK-1 protein expression in α-cells (18). This suggests that hyperglucagonemia in TALK-1 L114P islets is likely due to loss of inhibitory paracrine signaling. As insulin secretion does not change under fasting and euglycemic conditions, hyperglucagonemia might primarily be mediated by reduced somatostatin secretion. This is supported by a significant decrease in Sst expression in Kcnk16 L114P islets (Figure 6), which is likely due to reduced δ-cell secretion. Also, our previous data in global TALK-1 KO mice show higher islet somatostatin secretion and lower glucagon secretion, thus δ-cell TALK-1 L114P would be predicted to limit Ca2+ influx and somatostatin secretion (18). Intriguingly, Kcnk16 L114P islets additionally exhibit increased α-cell area fraction and a concurrent modest reduction in β-cell area fraction compared to control islets. These changes in islet composition are consistent with both type-1 diabetes (T1D) and T2D data showing increased α-cell: β-cell ratio (48, 49), which show elevated islet glucagon secretion. α-cell hyperplasia may result from increased activity/secretion, which is supported by other mediators of α-cell secretion (e.g., amino acids) that cause hyperplasia as well (50). However, loss of somatostatin also elevates α-cell secretion without altering α-cell mass (51). Moreover, the increased somatostatin and insulin secretion only result in reduced α-cell secretion without altering α-cell number. Thus, the exact mechanism of how TALK-1 L114P mediates increased α-cell number is not fully understood. Taken together, the KcnK16 L114P mouse model shows disrupted glucagon and insulin secretion leading to fasting hyperglycemia and glucose intolerance, which provides confirmation that TALK-1 GOF mutations likely cause MODY.

One obstacle in determining how MODY-associated mutations result in β-cell dysfunction is the limited availability of primary islet tissue from MODY families. Due to this, the initial assessment of MODY-associated human TALK-1 L114P was performed in overexpression systems, which resulted in a drastic GOF. Similarly, overexpression of mouse Kcnk16 L114P also leads to a significant GOF in TALK-1 activity (7.96-fold at −30 mV and 6.13-fold at 0 mV compared to TALK-1 WT). However, surprisingly β-cells from Kcnk16 L114P (L/P and P/P) mice showed only a modest GOF in K2P currents. While the slight increase in β-cell K+ conductance from Kcnk16 L114P mice would be predicted to partially alter islet excitability, these islets exhibit a complete loss of glucose-stimulated Vm depolarization. This suggests that endogenous TALK-1 L114P polarizes plasma membrane potential, which is further supported by the robust KCl-induced Vm depolarization and Ca2+ entry observed in these islets. Because KCl shifts the reversal potential of K+ channels to a more depolarized Vm, the constant Vm hyperpolarization in TALK-1 L114P β-cells likely results from increased K+ conductance through these channels. Yet, K+ conductance in TALK-1 L114P β-cells was not significantly different at the hyperpolarized membrane potentials which these cells reside at under low and high glucose conditions. The recording conditions may lead to poor isolation of the endogenous TALK-1 L114P currents; however, overexpression of this mutant channel results in large K+ currents under identical recording conditions. The K+ conductance differences between heterologously versus endogenously expressed TALK-1 L114P channels points towards unidentified regulators of β-cell TALK-1 activity which could include endogenous ligands, protein interactions, and cellular localization of the channel. This is presumably not due to changes in TALK-1 protein levels because Kcnk16 mRNA expression was not altered in control and Kcnk16 (c. 337 T>C) islets. The subtle increase in β-cell K+ conductance correlates with a modest MODY-like phenotype and likely allows for incomplete suppression of β-cell function. This is also observed in KATP-MODY, where channel activity is only modestly increased (21). Importantly, our study suggests that drastic TALK-1 GOF mutations only lead to modest β-cell K+ conductance, which may explain why both families with TALK-1 MODY carry the same pore-domain mutation (L114P) in TALK-1. Additionally, it is likely that other less-severe GOF mutations in TALK-1 (e.g., A277E) result in a milder phenotype such as T2D. Thus, it will be important to establish endogenous modulators of TALK-1 channels (WT and L114P) and how they contribute to β-cell dysfunction.

Suppression of glucose-stimulated electrical activity and Ca2+ influx in Kcnk16 L114P islets would be predicted to cause a greater reduction of GSIS than that observed. However, a few β-cells from Kcnk16 L114P mice show modest Vm depolarization; this resulted in a slight increase in Ca2+ in response to glucose in a small subset of Kcnk16 L114P islets, which could result in some GSIS. Interestingly, Kcnk16 L114P islets also show reduced expression of Fxyd3, which encodes the auxiliary subunit of Na+/K+-ATPase and is a known negative regulator of GSIS in diabetic mice and humans (52). Other interesting gene expression differences that may increase GSIS in Kcnk16 L114P islets include elevated expression (2.13-fold) of Adcy5, a Ca2+-independent amplification pathway gene (29, 53). In humans, Adcy5 depletion impairs glucose-dependent elevation of cAMP and associated insulin secretion, thus elevation in Adcy5 expression would be predicted to increase GSIS via cAMP signaling(29). Expression levels of other cAMP-dependent pathway genes were also elevated in the Kcnk16 L114P islets including Creb5, Adcyap1, and Adcyap1r1(30). Moreover, Slit1 and Srgap3, part of the SLIT-ROBO signaling which enhances GSIS, were upregulated (54). SLIT-ROBO signaling regulates not only Ca2+ handling but also actin remodeling and thus, Ca2+-independent signaling pathways. These gene expression changes suggest an increase in GSCI-independent mechanism(s) of insulin secretion may compensate for TALK-1 L114P-mediated loss of β-cell electrical activity and Ca2+ entry.

Chronic hyperglycemia in diabetic patients and rodents results in islet dysfunction and destruction, thus hyperglycemia observed in the Kcnk16 L114P mice could exacerbate β-cell failure. Glucotoxicity results in numerous islet transcriptome changes that in-part contribute to dysfunction. Similar to islets from diabetic patients, Kcnk16 L114P islets exhibit an increase in Aldob (124.6-fold) and Nnat (1.7-fold) expression (55, 56). Additionally, Pdk4 expression, a marker for the shift from utilization of glucose to fatty acids as the primary fuel source, is higher in Kcnk16 L114P islets (2.78-fold). Elevated islet PDK4 expression is also observed in patients with T2D and in animals on a HFD (37). Although hyperglycemia would be predicted to negatively impact Kcnk16 L114P islet function, no overt changes in β-cell mass were observed in these mice. This differs from mice with islets expressing KATP GOF mutation, which show loss of β-cell mass (57). Kcnk16 L114P may not cause a significant β-cell destruction because these mice show sufficient GSIS to prevent overt diabetes.

In summary, we showed that the MODY-associated TALK-1 L114P mutation elevates α-cell glucagon secretion under fasting and euglycemic conditions and blunts glucose-stimulated β-cell electrical activity and Ca2+ entry leading to reduced insulin secretion. Together, elevated glucagon impairs fasting glycemia and reduced GSIS increases post-prandial glucose levels in adults. This phenotype was more prominent in the male Kcnk16 L114P mice compared to the female Kcnk16 L114P mice raising the question whether sex-differences translate in humans carrying this mutation. Surprisingly, the TALK-1 L114P mutation also resulted in severe transient neonatal diabetes which was lethal in the C57Bl/6J genetic background. Thus, these data hold potential clinical utility in that neonatal diabetes patients with unknown genetic linkage should be screened for mutations in KCNK16. Together, these data strengthen the rationale for designing TALK-1 inhibitors for use as a therapeutic modality to treat diabetes.


Chemicals and reagents

All research materials were purchased from Thermo-Fisher (Waltham, MA, USA) or Sigma-Aldrich (St. Louis, MO, USA) unless otherwise specified.

Mouse model generation and ethical approval

Neonatal mice used for the studies were postnatal day 0- to postnatal day 5-old, and adult mice used for the studies were 6- to 26-week-old, age- and gender-matched, bred in-house on a C57BL/6J or mixed C57BL/6J:CD-1(ICR) background. Animals were handled in compliance with guidelines approved by the Vanderbilt University Animal Care and Use Committee protocols (#M2200007-00). C57BL/6J.Kcnk16L114P mice (Kcnk16<em1Djaco>; MGI:7486559) were produced by the Vanderbilt Genome Editing Resource (Vanderbilt University, Nashville, TN. USA). Ribonucleoprotein complexes (RNPs) comprising chemically modified ctRNA (crRNA + tracrRNA) (50 ng/µl) and enhanced specificity SpCas9 protein (100 ng/µl), together with a 180-nucleotide single-stranded DNA (ssDNA) donor containing the Kcnk16 L114P mutation (50 ng/ul), were obtained from MilliporeSigma (Burlington, MA, USA). These components were diluted in 10 mM Tris, 0.1 mM EDTA, pH 7.6, sourced from Teknova (Half Moon Bay, California, USA), and administered via pronuclear injection into C57BL/6J embryos acquired from mice from Jackson Labs (Bar Harbor, Maine, USA). crRNA sequence: 5’CCCTGCAGGTTATGGAAACC. 180 nucleotide ssDNA sequence: 5’ CTAGAGCTGGTGGTTGGGGGTGGGAGCCAGTTCTGGGCTCTCTTTTCCCCGCATC TGCACACTCCCTTGCCCTGCAGGTTATGGGAATCCAGCCCCCAGCACGGAGGCAG GGCAGGTCTTCTGTGTCTTCTATGCTCTGATGGGGATCCCACTCAATGTGGTCTTC CTCAACCATCTGGG. Mosaic F0 animals were screened for the L114P point mutation by standard PCR followed by a restriction fragment length polymorphism assay for a de novo HinfI site incorporated with silent mutations into the ssDNA. Animals carrying the desired mutation were confirmed by Sanger sequencing. Founder Kcnk16 L114P (L/P) mouse was backcrossed onto the C57Bl/6J (B6) strain for 2 generations to obtain mice used for all studies performed on the B6 background. Additionally, B6 Kcnk16 L114P (L/P) male mice were crossed with the CD-1 (ICR) strain to obtain F1 mice on a hybrid B6;CD-1 background to reduce the incidence of neonatal lethality. Crossings of heterozygous F1 B6;CD-1 Kcnk16 L114P (L/P) mice were used for generating homozygous B6;CD-1 Kcnk16 L114P (P/P) mice. For all the studies, littermates expressing the wildtype Kcnk16 allele were used as controls (WT).


Mouse pancreata were fixed in 4% paraformaldehyde and embedded with paraffin. Rehydrated 5-μm sections were stained with primary antibodies against insulin (dilution 1:1000; Dako, Santa Clara, CA, USA), somatostatin (dilution 1:300, GeneTex., Irvine, CA, USA), and glucagon (dilution 1:100; Abcam, Cambridge, United Kingdom) followed by secondary antibodies (dilution 1:500; anti-guinea pig Alexa Fluor 488, dilution 1:500; anti-mouse, Alexa Fluor 647, and dilution 1:500; anti-rabbit Alexa Fluor 546) as previously described(14). Sections were imaged either with a Nikon Eclipse TE2000-U microscope or fluorescent ScanScope (Aperio).

Islet isolation

Islets from neonatal mice were isolated on post-natal day 4 using the protocol described by Huang and Gu (58). Briefly, pancreata were isolated and broken into ∼2mm pieces and digested in 200ul collagenase P (Roche, Basel, Switzerland). For digestion, the tube was left in a 37 °C incubator for up to 15 min and inverted two times every min. The lysate was spun at 500 x g for 10 s followed by three washes in RPMI, after which the islets were handpicked using a brightfield microscope.

Islets from adult mouse pancreata were isolated by collagenase P digestion and density gradient centrifugation as previously described(14). Following isolation, islets were either dispersed into clusters of cells or single cells with trituration in 0.005% trypsin or maintained as whole islets. Cells were cultured in RPMI 1640 supplemented with 15% FBS, 100IU/mL penicillin, 100mg/mL streptomycin, and 5.5 mM glucose (RPMI) in a humidified incubator at 37°C with an atmosphere of 95% air and 5% CO2.

Whole-cell K2P currents

TALK-1 L114P currents were monitored using the whole-cell patch clamp technique using an Axopatch 200B amplifier with pCLAMP10 software. Digidata 1440 was used to digitize currents that were low-pass-filtered at 1 kHz and sampled at 10 kHz. Cells were washed with the extracellular buffer (modified Krebs-Ringer-HEPES buffer (KRHB)) containing (mM) 119.0 NaCl, 2.0 CaCl2, 4.7 KCl, 25.0 HEPES, 1.2 MgSO4, 1.2 KH2PO4, and 11 mM glucose (pH 7.4 with NaOH). For isolation of K2P channel currents, KATP channels were blocked with 100 μM tolbutamide, voltage-gated K+ channels were blocked with 10 mM tetraethylammonium (TEA)(18, 27). Patch electrodes (3–5 MΩ) were backfilled with intracellular solution (IC) containing (mM) 140.0 KCl, 1.0 MgCl2, 10.0 EGTA, 10.0 HEPES, and 4.0 Mg-ATP (pH 7.25 with KOH). β-cell Vm was ramped from −120 mV to +60 mV from a holding potential of −80 mV to generate K2P currents. Currents were measured in single β-cells from WT or Kcnk16 L114P (L/P) mice, or in HEK293FT cells expressing either Kcnk16 WT or mouse Kcnk16 L114P channels. The whole-cell currents were analyzed using ClampFit (Molecular Devices) and Excel (Microsoft Corp., Redmond, WA, USA).

For K2P recordings in HEK293FT cells, the cells were grown to ∼80% confluency in Dulbecco’s Modified Eagle Media (DMEM) GlutaMax-I (Thermo Fisher Scientific) supplemented with 10% fetal bovine serum (FBS, Gibco), 100 IU·ml−1 penicillin (Gibco), and 100 mg·ml−1 streptomycin (Gibco) at 37°C, 5% CO2 in 100 mm tissue culture dishes. Cells were transfected with either pLV-CMV-mKcnk16:P2A:EGFP or pLV-CMV-mKcnk16 L114P:P2A:EGFP plasmids using Lipofectamine 3000 and P3000 (Thermo Fisher Scientific) in antibiotic-free Opti-MEM™ I Reduced Serum Medium as per the manufacturer’s protocol. K2P currents were only recorded from EGFP-positive cells.

β-cell Vm recordings

β-cell Vm was recorded by the perforated patch clamp technique using an Axopatch 200B amplifier with pCLAMP10 software. Whole islets were washed with KRHB with (mM)

119.0 NaCl, 2.0 CaCl2, 4.7 KCl, 25.0 HEPES, 1.2 MgSO4, 1.2 KH2PO4 (adjusted to pH

7.4 with NaOH) supplemented with 2 mM glucose and incubated in KRHB for 30 min at 37°C, 5% CO2. Patch electrodes (3-5 MΩ) were backfilled with IC containing (mM) 140.0 KCl, 1.0 MgCl2, and 5.0 HEPES (adjusted to pH 7.2 with KOH) supplemented with 20 μg/mL amphotericin B. Islets were perifused with KRHB supplemented with 2 mM glucose followed by KRHB with 10 mM glucose for monitoring Vm changes. β-cell Vm recordings were analyzed using ClampFit (Molecular Devices), Excel (Microsoft Corp., Redmond, WA, USA), and GraphPad Prism 8 (GraphPad Software Inc.).

Intracellular Ca2+ imaging

On the day of experiment, islets were incubated for 30 min in RPMI supplemented with 2 μM Fura-2, AM (Molecular Probes) and 2 mM glucose. Fura-2, AM fluorescence (Ratio 340Ex/380Ex-535Em; F340/F380) was measured every 5 s as an indicator of intracellular Ca2+ using a Nikon Eclipse Ti2 microscope equipped with a Photometrics Prime 95B 25mm sCMOS Camera(59). For [Ca2+]c measurements, β-cell glucose-stimulated Ca2+ influx was monitored in KRHB supplemented with the glucose concentrations specified in the figures. For IP3-induced [Ca2+]ER release measurements, islets were perifused in KRHB buffer containing 11 mM glucose, 100 μM diazoxide, without extracellular Ca2+. Fura-2, AM fluorescence was monitored as an indicator of IP3-mediated [Ca2+]ER release upon stimulation of muscarinic receptor signaling by 100μM acetylcholine. For all measurements, the cells were perifused at a flow rate of 2 mL/min. Ex; Excitation wavelength (nm), Em; Emission wavelength (nm).

Glucose homeostasis

Chow-diet fed male and female mice underwent glucose tolerance test (GTT) and insulin tolerance test (ITT) as previously described(14). Briefly, mice were fasted for 4 h and then intraperitoneally (i.p.) injected with either 2 mg dextrose/g body weight for GTTs, or 0.75 UI human recombinant insulin/kg body weight for ITTs (catalog no. 12585014, Gibco™). Tail glucose measurements were then taken at the indicated time points in the figures to measure glucose clearance.

Body composition, tissue and plasma triglyceride and cholesterol measurements

Measurement of lean tissue, fat and fluid in living mice was performed using Bruker’s minispec Body Composition Analyzer. Plasma samples and livers were collected from ad-lib fed mice for triglyceride and cholesterol measurements. Total cholesterol and triglyceride was measured using standard enzymatic assays by the Vanderbilt University Medical Center Lipid Core.

Plasma insulin, in vitro insulin, and glucagon secretion assays

Plasma insulin from the neonates was measured on post-natal day 4 using Ultrasensitive Mouse Insulin ELISA kit (Catalog no. 10-1249-01). In adult mice, plasma insulin measurements were conducted in mice fasted for ∼4 hours followed by intraperitoneal injection of 2mg dextrose/g body weight at 0-, 15-, and 30-min post injection. Tail blood samples were collected in Microvette® CB 300 K2 EDTA tubes (Catalog no. 16.444.100, Sarstedt) at the indicated time points and plasma insulin were measured using mouse ultrasensitive insulin ELISA kit (Mercodia Inc., Sweden). For in vitro insulin and glucagon secretion assays, islets were isolated from mice fed a standard chow diet and were incubated over-night in RPMI supplemented with 0.5mg/mL BSA. On the following day, islets were equilibrated in DMEM containing 0.5 mg/mL BSA, 0.5 mM CaCl2 and 10.0 mM HEPES (DMEM*) supplemented with 10% FBS and 5.5 mM glucose for 1 h at 37°C, 5% CO2. 20 islets/well were picked into 400 uL DMEM* without FBS at glucose concentrations specified in the figures in 24-well plate(s) and insulin or glucagon secretion was measured over 1 h at 37°C and stored at −20°C until analysis. Insulin secretion was measured using mouse insulin ELISA kit (Catalog no. 10-1247-01, Mercodia Inc, Sweden) and glucagon secretion was measured using the mouse glucagon ELISA kit (and 10-1281-01, Mercodia Inc, Sweden).

Bulk-RNA sequencing

RNA was isolated from islets from ∼15-week-old male mice (WT and Kcnk16 L114P (L/P)) using Maxwell® 16 LEV simplyRNA Purification Kits (Catalog no. AS1280, Promega, USA). RNA integrity was analyzed using an Agilent 2100 Bioanalyzer, and only those samples with an RNA integrity number of seven or above were used. The sequencing libraries were constructed using the NEBNext Ultra II Directional RNA Library Prep Kit for Illumina P/N: E7760S/L by the Vanderbilt Technologies for Advanced Genomics (VANTAGE) core. An Illumina NovaSeq 6000 instrument was used to produce paired-end, 150-nucleotide reads for each RNA sample. Paired-end RNA sequencing reads (150bp long) were trimmed and filtered for quality using Trimgalore v0.6.7(60). Trimmed reads were aligned and counted using Spliced Transcripts Alignment to a Reference (STAR) (61) v2.7.9a with the –quantMode GeneCounts parameter against the mm10 mouse genome and GENCODE comprehensive gene annotations (Release M23). ∼178 million uniquely mapped reads were acquired per sample. DESeq2 package v1.34.0 (62) was used to perform normalization and downstream differential expression. Features counted fewer than 5 times across at least 3 samples were removed. Freezing condition of the samples and sequencing batch was included as batch factor in DESeq2 design to increase the sensitivity for finding differences between Kcnk16 L114P v/s WT samples. Gene enrichment analysis implemented from Gene Ontology was applied using the clusterProfiler v4.2.2 package in R. Annotated gene sets GO was sourced from Genome wide annotation for Mouse (Carlson M (2019). Genome wide annotation for Mouse. R package version 3.8.2.). For GO, genes significantly up-or downregulated in different conditions were used as input. False discovery rate-adjusted P value < 0.05 and log2 fold change >1 was used to define differentially expressed genes.

Quantitative PCR

RNA was isolated from islets from 15-week-old, chow diet fed mice using Maxwell® 16 LEV simplyRNA Purification Kits (Catalog no. AS1280, Promega, USA). Reverse transcription was performed using a SuperScript™ IV First-Strand Synthesis System (Catalog no. 18091050, Invitrogen, Waltham, MA, USA). 20 ng cDNA was used for real-time qPCR reactions with KAPA SYBR FAST qPCR Kit (Catalog no. KK4618, Roche, Basel, Switzerland) using a CFX Opus Real-Time PCR System (Bio-Rad Laboratories). Primers used for qRT-PCR are listed in Supplemental table 1.

Statistical Analysis

Functional data were analyzed using Axon Clampfit (Molecular Devices), GraphPad Prism 8 (GraphPad Software Inc.), or Excel (Microsoft Corp., Redmond, WA, USA) and presented as mean ± standard error (SE) for the specified number of samples (N). Statistical significance was determined using two-tailed t-tests, one-way ANOVA, or two-way ANOVA as appropriate. P-value ≤0.05 was considered statistically significant.


The Kcnk16 L114P mouse model was developed by the Vanderbilt Genome Editing Resource Core (RRID:SCR018826) which is supported by the Diabetes Research and Training Center Grant (DK020593), and the Cancer Center Support Grant (CA68485), and the Vanderbilt Center for Stem Cell Biology. Paraffin-embedded pancreata were prepared and processed by The Vanderbilt Translational Pathology Shared Resource (2P30 CA068485-14). Immunostaining slide scanning was performed using the Islet and Pancreas Analysis (IPA) Core supported by the Vanderbilt Diabetes Research Center (DRTC; NIH grant DK20593). Bulk-RNA sequencing was performed by the Vanderbilt Technologies for Advanced Genomics (VANTAGE) core, part of the Vanderbilt University Medical Center and analysis of the bulk RNA-seq data was performed by Creative Data Solutions, part of the Vanderbilt Center for Stem Cell Biology. Body composition measurements were performed at the Vanderbilt Mouse Metabolic Phenotyping Center which is supported by the NIH grant 5U2CDK059637. Triglyceride and cholesterol measurements were performed by the Vanderbilt University Medical Center Lipid Core supported by the NIH DRTC grant DK020593. Research in the laboratory of DAJ was supported by NIH grants R01DK097392, R01DK129340, and R01DK115620. A few illustrations were created using

Author contributions

DAJ and AYN conceived the project. AYN, PKD, JK, SS, and JP conducted the experiments and/or analyzed/reviewed the data. LS and MAM developed the mouse model. AYN, SS, JP, LS, MAM and DAJ wrote and/or edited the manuscript.

Declaration of interests

All authors declare no competing interests.