Pathogenic monoallelic variants in GLIS3 increase type 2 diabetes risk and identify a subgroup of patients sensitive to sulfonylureas.

AIMS/HYPOTHESIS: GLIS3 encodes a transcription factor involved in pancreatic beta cell development and function. Rare pathogenic, bi-allelic mutations in GLIS3 cause syndromic neonatal diabetes whereas frequent SNPs at this locus associate with common type 2 diabetes risk. Because rare, functional variants located in other susceptibility genes for type 2 diabetes have already been shown to strongly increase individual risk for common type 2 diabetes, we aimed to investigate the contribution of rare pathogenic GLIS3 variants to type 2 diabetes. METHODS: GLIS3 was sequenced in 5471 individuals from the Rare Variants Involved in Diabetes and Obesity (RaDiO) study. Variant pathogenicity was assessed following the criteria established by the American College of Medical Genetics and Genomics (ACMG). To address the pathogenic strong criterion number 3 (PS3), we conducted functional investigations of these variants using luciferase assays, focusing on capacity of GLIS family zinc finger 3 (GLIS3) to bind to and activate the INS promoter. The association between rare pathogenic or likely pathogenic (P/LP) variants and type 2 diabetes risk (and other metabolic traits) was then evaluated. A meta-analysis combining association results from RaDiO, the 52K study (43,125 individuals) and the TOPMed study (44,083 individuals) was finally performed. RESULTS: Through targeted resequencing of GLIS3, we identified 105 rare variants that were carried by 395 participants from RaDiO. Among them, 49 variants decreased the activation of the INS promoter. Following ACMG criteria, 18 rare variants were classified as P/LP, showing an enrichment in the last two exons compared with the remaining exons (p<5×10-6; OR>3.5). The burden of these P/LP variants was strongly higher in individuals with type 2 diabetes (p=3.0×10-3; OR 3.9 [95% CI 1.4, 12]), whereas adiposity, age at type 2 diabetes diagnosis and cholesterol levels were similar between variant carriers and non-carriers with type 2 diabetes. Interestingly, all carriers with type 2 diabetes were sensitive to oral sulfonylureas. A total of 7 P/LP variants were identified in both 52K and TOPMed studies. The meta-analysis of association studies obtained from RaDiO, 52K and TOPMed showed an enrichment of P/LP GLIS3 variants in individuals with type 2 diabetes (p=5.6×10-5; OR 2.1 [95% CI 1.4, 2.9]). CONCLUSIONS/INTERPRETATION: Rare P/LP GLIS3 variants do contribute to type 2 diabetes risk. The variants located in the distal part of the protein could have a direct effect on its functional activity by impacting its transactivation domain, by homology with the mouse GLIS3 protein. Furthermore, rare P/LP GLIS3 variants seem to have a direct clinical effect on beta cell function, which could be improved by increasing insulin secretion via the use of sulfonylureas.

Control of human pancreatic beta cell kinome by glucagon-like peptide-1 receptor biased agonism.

AIM: To determine the kinase activity profiles of human pancreatic beta cells downstream of glucagon-like peptide-1 receptor (GLP-1R) balanced versus biased agonist stimulations. MATERIALS AND METHODS: This study analysed the kinomic profiles of human EndoC-ßh1 cells following vehicle and GLP-1R stimulation with the pharmacological agonist exendin-4, as well as exendin-4-based biased derivatives exendin-phe1 and exendin-asp3 for acute (10-minute) versus sustained (120-minute) responses, using PamChip protein tyrosine kinase and serine/threonine kinase assays. The raw data were filtered and normalized using BioNavigator. The kinase analyses were conducted with R, mainly including kinase-substrate mapping and Kyoto Encyclopedia of Genes and Genomes pathway analysis. RESULTS: The present analysis reveals that kinomic responses are distinct for acute versus sustained GLP-1R agonist exposure, with individual responses associated with agonists presenting specific bias profiles. According to pathway analysis, several kinases, including JNKs, PKCs, INSR and LKB1, are important GLP-1R signalling mediators, constituting potential targets for further research on biased GLP-1R downstream signalling. CONCLUSION: The results from this study suggest that differentially biased exendin-phe1 and exendin-asp3 can modulate distinct kinase interaction networks. Further understanding of these mechanisms will have important implications for the selection of appropriate anti-type 2 diabetes therapies with optimized downstream kinomic profiles.

The Map3k12 (Dlk)/JNK3 signaling pathway is required for pancreatic beta-cell proliferation during postnatal development.

Unveiling the key pathways underlying postnatal beta-cell proliferation can be instrumental to decipher the mechanisms of beta-cell mass plasticity to increased physiological demand of insulin during weight gain and pregnancy. Using transcriptome and global Serine Threonine Kinase activity (STK) analyses of islets from newborn (10 days old) and adult rats, we found that highly proliferative neonatal rat islet cells display a substantially elevated activity of the mitogen activated protein 3 kinase 12, also called dual leucine zipper-bearing kinase (Dlk). As a key upstream component of the c-Jun amino terminal kinase (Jnk) pathway, Dlk overexpression was associated with increased Jnk3 activity and was mainly localized in the beta-cell cytoplasm. We provide the evidence that Dlk associates with and activates Jnk3, and that this cascade stimulates the expression of Ccnd1 and Ccnd2, two essential cyclins controlling postnatal beta-cell replication. Silencing of Dlk or of Jnk3 in neonatal islet cells dramatically hampered primary beta-cell replication and the expression of the two cyclins. Moreover, the expression of Dlk, Jnk3, Ccnd1 and Ccnd2 was induced in high replicative islet beta cells from ob/ob mice during weight gain, and from pregnant female rats. In human islets from non-diabetic obese individuals, DLK expression was also cytoplasmic and the rise of the mRNA level was associated with an increase of JNK3, CCND1 and CCND2 mRNA levels, when compared to islets from lean and obese patients with diabetes. In conclusion, we find that activation of Jnk3 signalling by Dlk could be a key mechanism for adapting islet beta-cell mass during postnatal development and weight gain.

Pathogenic, Total Loss-of-Function DYRK1B Variants Cause Monogenic Obesity Associated With Type 2 Diabetes.

OBJECTIVE: Rare variants in DYRK1B have been described in some patients with central obesity, type 2 diabetes, and early-onset coronary disease. Owing to the limited number of conducted studies, the broader impact of DYRK1B variants on a larger scale has yet to be investigated. RESEARCH DESIGN AND METHODS: DYRK1B was sequenced in 9,353 participants from a case-control study for obesity and type 2 diabetes. Each DYRK1B variant was functionally assessed in vitro. Variant pathogenicity was determined using criteria from the American College of Medical Genetics and Genomics (ACMG). The effect of pathogenic or likely pathogenic (P/LP) variants on metabolic traits was assessed using adjusted mixed-effects score tests. RESULTS: Sixty-five rare, heterozygous DYRK1B variants were identified and were not associated with obesity or type 2 diabetes. Following functional analyses, 20 P/LP variants were pinpointed, including 6 variants that exhibited a fully inhibitory effect (P/LP-null) on DYRK1B activity. P/LP and P/LP-null DYRK1B variants were associated with increased BMI and obesity risk; however, the impact was notably more pronounced for the P/LP-null variants (effect of 8.0 ± 3.2 and odds ratio of 7.9 [95% CI 1.2-155]). Furthermore, P/LP-null variants were associated with higher fasting glucose and type 2 diabetes risk (effect of 2.9 ± 1.0 and odds ratio of 4.8 [95% CI 0.85-37]), while P/LP variants had no effect on glucose homeostasis. CONCLUSIONS: P/LP, total loss-of-function DYRK1B variants cause monogenic obesity associated with type 2 diabetes. This study underscores the significance of conducting functional assessments in order to accurately ascertain the tangible effects of P/LP DYRK1B variants.

Exploring the role of purinergic receptor P2RY1 in type 2 diabetes risk and pathophysiology: Insights from human functional genomics.

OBJECTIVE: Human functional genomics has proven powerful in discovering drug targets for common metabolic disorders. Through this approach, we investigated the involvement of the purinergic receptor P2RY1 in type 2 diabetes (T2D). METHODS: P2RY1 was sequenced in 9,266 participants including 4,177 patients with T2D. In vitro analyses were then performed to assess the functional effect of each variant. Expression quantitative trait loci (eQTL) analysis was performed in pancreatic islets from 103 pancreatectomized individuals. The effect of P2RY1 on glucose-stimulated insulin secretion was finally assessed in human pancreatic beta cells (EndoCßH5), and RNA sequencing was performed on these cells. RESULTS: Sequencing P2YR1 in 9,266 participants revealed 22 rare variants, seven of which were loss-of-function according to our in vitro analyses. Carriers, except one, exhibited impaired glucose control. Our eQTL analysis of human islets identified P2RY1 variants, in a beta-cell enhancer, linked to increased P2RY1 expression and reduced T2D risk, contrasting with variants located in a silent region associated with decreased P2RY1 expression and increased T2D risk. Additionally, a P2RY1-specific agonist increased insulin secretion upon glucose stimulation, while the antagonist led to decreased insulin secretion. RNA-seq highlighted TXNIP as one of the main transcriptomic markers of insulin secretion triggered by P2RY1 agonist. CONCLUSION: Our findings suggest that P2RY1 inherited or acquired dysfunction increases T2D risk and that P2RY1 activation stimulates insulin secretion. Selective P2RY1 agonists, impermeable to the blood-brain barrier, could serve as potential insulin secretagogues.

Knocking Down CDKN2A in 3D hiPSC-Derived Brown Adipose Progenitors Potentiates Differentiation, Oxidative Metabolism and Browning Process.

Human induced pluripotent stem cells (hiPSCs) have the potential to be differentiated into any cell type, making them a relevant tool for therapeutic purposes such as cell-based therapies. In particular, they show great promise for obesity treatment as they represent an unlimited source of brown/beige adipose progenitors (hiPSC-BAPs). However, the low brown/beige adipocyte differentiation potential in 2D cultures represents a strong limitation for clinical use. In adipose tissue, besides its cell cycle regulator functions, the cyclin-dependent kinase inhibitor 2A (CDKN2A) locus modulates the commitment of stem cells to the brown-like type fate, mature adipocyte energy metabolism and the browning of adipose tissue. Here, using a new method of hiPSC-BAPs 3D culture, via the formation of an organoid-like structure, we silenced CDKN2A expression during hiPSC-BAP adipogenic differentiation and observed that knocking down CDKN2A potentiates adipogenesis, oxidative metabolism and the browning process, resulting in brown-like adipocytes by promoting UCP1 expression and beiging markers. Our results suggest that modulating CDKN2A levels could be relevant for hiPSC-BAPs cell-based therapies.

The transcription factor E2F1 controls the GLP-1 receptor pathway in pancreatic ß cells.

The glucagon-like peptide 1 (Glp-1) has emerged as a hormone with broad pharmacological potential in type 2 diabetes (T2D) treatment, notably by improving ß cell functions. The cell-cycle regulator and transcription factor E2f1 is involved in glucose homeostasis by modulating ß cell mass and function. Here, we report that ß cell-specific genetic ablation of E2f1 (E2f1ß-/-) impairs glucose homeostasis associated with decreased expression of the Glp-1 receptor (Glp1r) in E2f1ß-/- pancreatic islets. Pharmacological inhibition of E2F1 transcriptional activity in nondiabetic human islets decreases GLP1R levels and blunts the incretin effect of GLP1R agonist exendin-4 (ex-4) on insulin secretion. Overexpressing E2f1 in pancreatic ß cells increases Glp1r expression associated with enhanced insulin secretion mediated by ex-4. Interestingly, ex-4 induces retinoblastoma protein (pRb) phosphorylation and E2f1 transcriptional activity. Our findings reveal critical roles for E2f1 in ß cell function and suggest molecular crosstalk between the E2F1/pRb and GLP1R signaling pathways.

Histone deacetylase 9 promoter hypomethylation associated with adipocyte dysfunction is a statin-related metabolic effect.

BACKGROUND: Adipogenesis, the process whereby preadipocytes differentiate into mature adipocytes, is crucial for maintaining metabolic homeostasis. Cholesterol-lowering statins increase type 2 diabetes (T2D) risk possibly by affecting adipogenesis and insulin resistance but the (epi)genetic mechanisms involved are unknown. Here, we characterised the effects of statin treatment on adipocyte differentiation using in vitro human preadipocyte cell model to identify putative effective genes. RESULTS: Statin treatment during adipocyte differentiation caused a reduction in key genes involved in adipogenesis, such as ADIPOQ, GLUT4 and ABCG1. Using Illumina's Infinium '850K' Methylation EPIC array, we found a significant hypomethylation of cg14566882, located in the promoter of the histone deacetylase 9 (HDAC9) gene, in response to two types of statins (atorvastatin and mevastatin), which correlates with an increased HDAC9 mRNA expression. We confirmed that HDAC9 is a transcriptional repressor of the cholesterol efflux ABCG1 gene expression, which is epigenetically modified in obesity and prediabetic states. Thus, we assessed the putative impact of ABCG1 knockdown in mimicking the effect of statin in adipogenesis. ABCG1 KD reduced the expression of key genes involved in adipocyte differentiation and decreased insulin signalling and glucose uptake. In human blood cells from two cohorts, ABCG1 expression was impaired in response to statins, confirming that ABCG1 is targeted in vivo by these drugs. CONCLUSIONS: We identified an epigenetic link between adipogenesis and adipose tissue insulin resistance in the context of T2D risk associated with statin use, which has important implications as HDAC9 and ABCG1 are considered potential therapeutic targets for obesity and metabolic diseases.

Loss-of-function mutations in MRAP2 are pathogenic in hyperphagic obesity with hyperglycemia and hypertension.

The G-protein-coupled receptor accessory protein MRAP2 is implicated in energy control in rodents, notably via the melanocortin-4 receptor1. Although some MRAP2 mutations have been described in people with obesity1-3, their functional consequences on adiposity remain elusive. Using large-scale sequencing of MRAP2 in 9,418 people, we identified 23 rare heterozygous variants associated with increased obesity risk in both adults and children. Functional assessment of each variant shows that loss-of-function MRAP2 variants are pathogenic for monogenic hyperphagic obesity, hyperglycemia and hypertension. This contrasts with other monogenic forms of obesity characterized by excessive hunger, including melanocortin-4 receptor deficiency, that present with low blood pressure and normal glucose tolerance4. The pleiotropic metabolic effect of loss-of-function mutations in MRAP2 might be due to the failure of different MRAP2-regulated G-protein-coupled receptors in various tissues including pancreatic islets.

Loss-of-function mutations in ADCY3 cause monogenic severe obesity.

Study of monogenic forms of obesity has demonstrated the pivotal role of the central leptin-melanocortin pathway in controlling energy balance, appetite and body weight 1 . The majority of loss-of-function mutations (mostly recessive or co-dominant) have been identified in genes that are directly involved in leptin-melanocortin signaling. These genes, however, only explain obesity in <5% of cases, predominantly from outbred populations 2 . We previously showed that, in a consanguineous population in Pakistan, recessive mutations in known obesity-related genes explain ~30% of cases with severe obesity3-5. These data suggested that new monogenic forms of obesity could also be identified in this population. Here we identify and functionally characterize homozygous mutations in the ADCY3 gene encoding adenylate cyclase 3 in children with severe obesity from consanguineous Pakistani families, as well as compound heterozygous mutations in a severely obese child of European-American descent. These findings highlight ADCY3 as an important mediator of energy homeostasis and an attractive pharmacological target in the treatment of obesity.

Increased Hepatic PDGF-AA Signaling Mediates Liver Insulin Resistance in Obesity-Associated Type 2 Diabetes.

In type 2 diabetes (T2D), hepatic insulin resistance is strongly associated with nonalcoholic fatty liver disease (NAFLD). In this study, we hypothesized that the DNA methylome of livers from patients with T2D compared with livers of individuals with normal plasma glucose levels can unveil some mechanism of hepatic insulin resistance that could link to NAFLD. Using DNA methylome and transcriptome analyses of livers from obese individuals, we found that hypomethylation at a CpG site in PDGFA (encoding platelet-derived growth factor α) and PDGFA overexpression are both associated with increased T2D risk, hyperinsulinemia, increased insulin resistance, and increased steatohepatitis risk. Genetic risk score studies and human cell modeling pointed to a causative effect of high insulin levels on PDGFA CpG site hypomethylation, PDGFA overexpression, and increased PDGF-AA secretion from the liver. We found that PDGF-AA secretion further stimulates its own expression through protein kinase C activity and contributes to insulin resistance through decreased expression of insulin receptor substrate 1 and of insulin receptor. Importantly, hepatocyte insulin sensitivity can be restored by PDGF-AA-blocking antibodies, PDGF receptor inhibitors, and by metformin, opening therapeutic avenues. Therefore, in the liver of obese patients with T2D, the increased PDGF-AA signaling contributes to insulin resistance, opening new therapeutic avenues against T2D and possibly NAFLD.

Expression and functional assessment of candidate type 2 diabetes susceptibility genes identify four new genes contributing to human insulin secretion.

OBJECTIVES: Genome-wide association studies (GWAS) have identified >100 loci independently contributing to type 2 diabetes (T2D) risk. However, translational implications for precision medicine and for the development of novel treatments have been disappointing, due to poor knowledge of how these loci impact T2D pathophysiology. Here, we aimed to measure the expression of genes located nearby T2D associated signals and to assess their effect on insulin secretion from pancreatic beta cells. METHODS: The expression of 104 candidate T2D susceptibility genes was measured in a human multi-tissue panel, through PCR-free expression assay. The effects of the knockdown of beta-cell enriched genes were next investigated on insulin secretion from the human EndoC-ßH1 beta-cell line. Finally, we performed RNA-sequencing (RNA-seq) so as to assess the pathways affected by the knockdown of the new genes impacting insulin secretion from EndoC-ßH1, and we analyzed the expression of the new genes in mouse models with altered pancreatic beta-cell function. RESULTS: We found that the candidate T2D susceptibility genes' expression is significantly enriched in pancreatic beta cells obtained by laser capture microdissection or sorted by flow cytometry and in EndoC-ßH1 cells, but not in insulin sensitive tissues. Furthermore, the knockdown of seven T2D-susceptibility genes (CDKN2A, GCK, HNF4A, KCNK16, SLC30A8, TBC1D4, and TCF19) with already known expression and/or function in beta cells changed insulin secretion, supporting our functional approach. We showed first evidence for a role in insulin secretion of four candidate T2D-susceptibility genes (PRC1, SRR, ZFAND3, and ZFAND6) with no previous knowledge of presence and function in beta cells. RNA-seq in EndoC-ßH1 cells with decreased expression of PRC1, SRR, ZFAND6, or ZFAND3 identified specific gene networks related to T2D pathophysiology. Finally, a positive correlation between the expression of Ins2 and the expression of Prc1, Srr, Zfand6, and Zfand3 was found in mouse pancreatic islets with altered beta-cell function. CONCLUSIONS: This study showed the ability of post-GWAS functional studies to identify new genes and pathways involved in human pancreatic beta-cell function and in T2D pathophysiology.

Impaired histone deacetylases 5 and 6 expression mimics the effects of obesity and hypoxia on adipocyte function.

OBJECTIVE: The goal of the study was to investigate the role of histone deacetylases (HDACs) in adipocyte function associated with obesity and hypoxia. METHODS: Total proteins and RNA were prepared from human visceral adipose tissues (VAT) of human obese and normal weight subjects and from white adipose tissue (WAT) of C57Bl6-Rj mice fed a normal or high fat diet (HFD) for 16 weeks. HDAC activity was measured by colorimetric assay whereas the gene and protein expression were monitored by real-time PCR and by western blotting, respectively. RNA interference (RNAi) was used to silence the expression of genes in 3T3-L1 adipocytes. RESULTS: Total HDAC activity was decreased in VAT and WAT from obese individuals and from mice fed a HFD, respectively. The HDAC activity reduction was associated with decreased HDAC5/Hdac5 and HDAC6/Hdac6 expression in human and mice adipocyte fraction. Similarly, hypoxia hampered total Hdac activity and reduced the expression of Hdac5 and Hdac6 in 3T3-L1 adipocytes. The decrease of both Hdac5 and Hdac6 by hypoxia was associated with altered expression of adipokines and of the inducible cAMP early repressor (Icer), a key repressor that is defective in human and mice obesity. Silencing of Icer in adipocytes reproduced the changes in adipokine levels under hypoxia and obesity, suggesting a causative effect. Finally, modeling the defect of the two Hdacs in adipocytes by RNAi or selective inhibitors mimicked the effects of hypoxia on the expression of Icer, leading to impairment of insulin-induced glucose uptake. CONCLUSION: Hdac5 and Hdac6 expression are required for the adequate expression of Icer and adipocyte function. Altered adipose expression of the two Hdacs in obesity by hypoxia may contribute to the development of metabolic abnormalities.

JNK3 is required for the cytoprotective effect of exendin 4.

Preservation of beta cell against apoptosis is one of the therapeutic benefits of the glucagon-like peptide-1 (GLP1) antidiabetic mimetics for preserving the functional beta cell mass exposed to diabetogenic condition including proinflammatory cytokines. The mitogen activated protein kinase 10 also called c-jun amino-terminal kinase 3 (JNK3) plays a protective role in insulin-secreting cells against death caused by cytokines. In this study, we investigated whether the JNK3 expression is associated with the protective effect elicited by the GLP1 mimetic exendin 4. We found an increase in the abundance of JNK3 in isolated human islets and INS-1E cells cultured with exendin 4. Induction of JNK3 by exendin 4 was associated with an increased survival of INS-1E cells. Silencing of JNK3 prevented the cytoprotective effect of exendin 4 against apoptosis elicited by culture condition and cytokines. These results emphasize the requirement of JNK3 in the antiapoptotic effects of exendin 4.