Nicola L. Beer

Update on mutations in glucokinase (GCK), which cause maturity‐onset diabetes of the young, permanent neonatal diabetes, and hyperinsulinemic hypoglycemia

Glucokinase is a key regulatory enzyme in the pancreatic beta‐cell. It plays a crucial role in the regulation of insulin secretion and has been termed the glucose sensor in pancreatic beta‐cells. Given its central role in the regulation of insulin release it is understandable that mutations in the gene encoding glucokinase (GCK) can cause both hyper‐ and hypoglycemia. Heterozygous inactivating mutations in GCK cause maturity‐onset diabetes of the young (MODY) subtype glucokinase (GCK), characterized by mild fasting hyperglycemia, which is present at birth but often only detected later in life during screening for other purposes. Homozygous inactivating GCK mutations result in a more severe phenotype presenting at birth as permanent neonatal diabetes mellitus (PNDM). A growing number of heterozygous activating GCK mutations that cause hypoglycemia have also been reported. A total of 620 mutations in the GCK gene have been described in a total of 1,441 families. There are no common mutations, and the mutations are distributed throughout the gene. The majority of activating mutations cluster in a discrete region of the protein termed the allosteric activator site. The identification of a GCK mutation in patients with both hyper‐ and hypoglycemia has implications for the clinical course and clinical management of their disorder. Hum Mutat 30: 1–15, 2009.

Loss-of-function mutations in SLC30A8 protect against type 2 diabetes

Loss-of-function mutations protective against human disease provide in vivo validation of therapeutic targets, but none have yet been described for type 2 diabetes (T2D). Through sequencing or genotyping of ∼150,000 individuals across 5 ancestry groups, we identified 12 rare protein-truncating variants in SLC30A8, which encodes an islet zinc transporter (ZnT8) and harbors a common variant (p.Trp325Arg) associated with T2D risk and glucose and proinsulin levels. Collectively, carriers of protein-truncating variants had 65% reduced T2D risk (P = 1.7 × 10−6), and non-diabetic Icelandic carriers of a frameshift variant (p.Lys34Serfs*50) demonstrated reduced glucose levels (−0.17 s.d., P = 4.6 × 10−4). The two most common protein-truncating variants (p.Arg138* and p.Lys34Serfs*50) individually associate with T2D protection and encode unstable ZnT8 proteins. Previous functional study of SLC30A8 suggested that reduced zinc transport increases T2D risk, and phenotypic heterogeneity was observed in mouse Slc30a8 knockouts. In contrast, loss-of-function mutations in humans provide strong evidence that SLC30A8 haploinsufficiency protects against T2D, suggesting ZnT8 inhibition as a therapeutic strategy in T2D prevention.

The P446L variant in GCKR associated with fasting plasma glucose and triglyceride levels exerts its effect through increased glucokinase activity in liver

Genome-wide association studies have identified a number of signals for both Type 2 Diabetes and related quantitative traits. For the majority of loci, the transition from association signal to mutational mechanism has been difficult to establish. Glucokinase (GCK) regulates glucose storage and disposal in the liver where its activity is regulated by glucokinase regulatory protein (GKRP; gene name GCKR). Fructose-6 and fructose-1 phosphate (F6P and F1P) enhance or reduce GKRP-mediated inhibition, respectively. A common GCKR variant (P446L) is reproducibly associated with triglyceride and fasting plasma glucose levels in the general population. The aim of this study was to determine the mutational mechanism responsible for this genetic association. Recombinant human GCK and both human wild-type (WT) and P446L-GKRP proteins were generated. GCK kinetic activity was observed spectrophotometrically using an NADP+-coupled assay. WT and P446L-GKRP-mediated inhibition of GCK activity and subsequent regulation by phosphate esters were determined. Assays matched for GKRP activity demonstrated no difference in dose-dependent inhibition of GCK activity or F1P-mediated regulation. However, the response to physiologically relevant F6P levels was significantly attenuated with P446L-GKRP (n = 18; P ≤ 0.03). Experiments using equimolar concentrations of both regulatory proteins confirmed these findings (n = 9; P < 0.001). In conclusion, P446L-GKRP has reduced regulation by physiological concentrations of F6P, resulting indirectly in increased GCK activity. Altered GCK regulation in liver is predicted to enhance glycolytic flux, promoting hepatic glucose metabolism and elevating concentrations of malonyl-CoA, a substrate for de novo lipogenesis, providing a mutational mechanism for the reported association of this variant with raised triglycerides and lower glucose levels.

Assessing the phenotypic effects in the general population of rare variants in genes for a dominant Mendelian form of diabetes.

Genome sequencing can identify individuals in the general population who harbor rare coding variants in genes for Mendelian disorders and who may consequently have increased disease risk. Previous studies of rare variants in phenotypically extreme individuals display ascertainment bias and may demonstrate inflated effect-size estimates. We sequenced seven genes for maturity-onset diabetes of the young (MODY) in well-phenotyped population samples (n = 4,003). We filtered rare variants according to two prediction criteria for disease-causing mutations: reported previously in MODY or satisfying stringent de novo thresholds (rare, conserved and protein damaging). Approximately 1.5% and 0.5% of randomly selected individuals from the Framingham and Jackson Heart Studies, respectively, carry variants from these two classes. However, the vast majority of carriers remain euglycemic through middle age. Accurate estimates of variant effect sizes from population-based sequencing are needed to avoid falsely predicting a substantial fraction of individuals as being at risk for MODY or other Mendelian diseases.

Identification and Functional Characterization of G6PC2 Coding Variants Influencing Glycemic Traits Define an Effector Transcript at the G6PC2-ABCB11 Locus

Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P<5×10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights.

Correlation of rare coding variants in the gene encoding human glucokinase regulatory protein with phenotypic, cellular, and kinetic outcomes

Defining the genetic contribution of rare variants to common diseases is a major basic and clinical science challenge that could offer new insights into disease etiology and provide potential for directed gene- and pathway-based prevention and treatment. Common and rare nonsynonymous variants in the GCKR gene are associated with alterations in metabolic traits, most notably serum triglyceride levels. GCKR encodes glucokinase regulatory protein (GKRP), a predominantly nuclear protein that inhibits hepatic glucokinase (GCK) and plays a critical role in glucose homeostasis. The mode of action of rare GCKR variants remains unexplored. We identified 19 nonsynonymous GCKR variants among 800 individuals from the ClinSeq medical sequencing project. Excluding the previously described common missense variant p.Pro446Leu, all variants were rare in the cohort. Accordingly, we functionally characterized all variants to evaluate their potential phenotypic effects. Defects were observed for the majority of the rare variants after assessment of cellular localization, ability to interact with GCK, and kinetic activity of the encoded proteins. Comparing the individuals with functional rare variants to those without such variants showed associations with lipid phenotypes. Our findings suggest that, while nonsynonymous GCKR variants, excluding p.Pro446Leu, are rare in individuals of mixed European descent, the majority do affect protein function. In sum, this study utilizes computational, cell biological, and biochemical methods to present a model for interpreting the clinical significance of rare genetic variants in common disease.

Identification and functional characterisation of novel glucokinase mutations causing maturity-onset diabetes of the young in Slovakia.

Heterozygous glucokinase (GCK) mutations cause a subtype of maturity-onset diabetes of the young (GCK-MODY). Over 600 GCK mutations have been reported of which ∼65% are missense. In many cases co-segregation has not been established and despite the importance of functional studies in ascribing pathogenicity for missense variants these have only been performed for <10% of mutations. The aim of this study was to determine the minimum prevalence of GCK-MODY amongst diabetic subjects in Slovakia by sequencing GCK in 100 Slovakian probands with a phenotype consistent with GCK-MODY and to explore the pathogenicity of identified variants through family and functional studies. Twenty-two mutations were identified in 36 families (17 missense) of which 7 (I110N, V200A, N204D, G258R, F419S, c.580-2A>C, c.1113–1114delGC) were novel. Parental DNA was available for 22 probands (covering 14/22 mutations) and co-segregation established in all cases. Bioinformatic analysis predicted all missense mutations to be damaging. Nine (I110N, V200A, N204D, G223S, G258R, F419S, V244G, L315H, I436N) mutations were functionally evaluated. Basic kinetic analysis explained pathogenicity for 7 mutants which showed reduced glucokinase activity with relative activity indices (RAI) between 0.6 to <0.001 compared to wild-type GCK (1.0). For the remaining 2 mutants additional molecular mechanisms were investigated. Differences in glucokinase regulatory protein (GKRP) –mediated-inhibition of GCK were observed for both L315H & I436N when compared to wild type (IC50 14.6±0.1 mM & 20.3±1.6 mM vs.13.3±0.1 mM respectively [p<0.03]). Protein instability as assessed by thermal lability studies demonstrated that both L315H and I436N show marked thermal instability compared to wild-type GCK (RAI at 55°C 8.8±0.8% & 3.1±0.4% vs. 42.5±3.9% respectively [p<0.001]). The minimum prevalence of GCK-MODY amongst Slovakian patients with diabetes was 0.03%. In conclusion, we have identified 22 GCK mutations in 36 Slovakian probands and demonstrate that combining family, bioinformatic and functional studies can aid the interpretation of variants identified by molecular diagnostic screening.

Discovery of a Novel Site Regulating Glucokinase Activity following Characterization of a New Mutation Causing Hyperinsulinemic Hypoglycemia in Humans

Type 2 diabetes is a global problem, and current ineffective therapeutic strategies pave the way for novel treatments like small molecular activators targeting glucokinase (GCK). GCK activity is fundamental to beta cell and hepatocyte glucose metabolism, and heterozygous activating and inactivating GCK mutations cause hyperinsulinemic hypoglycemia (HH) and maturity onset diabetes of the young (MODY) respectively. Over 600 naturally occurring inactivating mutations have been reported, whereas only 13 activating mutations are documented to date. We report two novel GCK HH mutations (V389L and T103S) at residues where MODY mutations also occur (V389D and T103I). Using recombinant proteins with in vitro assays, we demonstrated that both HH mutants had a greater relative activity index than wild type (6.0 for V389L, 8.4 for T103S, and 1.0 for wild type). This was driven by an increased affinity for glucose (S0.5, 3.3 ± 0.1 and 3.5 ± 0.1 mm, respectively) versus wild type (7.5 ± 0.1 mm). Correspondingly, the V389D and T103I MODY mutants had markedly reduced relative activity indexes (<0.1). T103I had an altered affinity for glucose (S0.5, 24.9 ± 0.6 mm), whereas V389D also exhibited a reduced affinity for ATP and decreased catalysis rate (S0.5, 78.6 ± 4.5 mm; ATPKm, 1.5 ± 0.1 mm; Kcat, 10.3 ± 1.1s−1) compared with wild type (ATPKm, 0.4 ± <0.1; Kcat, 62.9 ± 1.2). Both Thr-103 mutants showed reduced inhibition by the endogenous hepatic inhibitor glucokinase regulatory protein. Molecular modeling demonstrated that Thr-103 maps to the allosteric activator site, whereas Val-389 is located remotely to this position and all other previously reported activating mutations, highlighting α-helix 11 as a novel region regulating GCK activity. Our data suggest that pharmacological manipulation of GCK activity at locations distal from the allosteric activator site is possible.

Insights into islet development and biology through characterization of a human iPSC-derived endocrine pancreas model

ABSTRACT Directed differentiation of stem cells offers a scalable solution to the need for human cell models recapitulating islet biology and T2D pathogenesis. We profiled mRNA expression at 6 stages of an induced pluripotent stem cell (iPSC) model of endocrine pancreas development from 2 donors, and characterized the distinct transcriptomic profiles associated with each stage. Established regulators of endodermal lineage commitment, such as SOX17 (log2 fold change [FC] compared to iPSCs = 14.2, p-value = 4.9 × 10−5) and the pancreatic agenesis gene GATA6 (log2 FC = 12.1, p-value = 8.6 × 10−5), showed transcriptional variation consistent with their known developmental roles. However, these analyses highlighted many other genes with stage-specific expression patterns, some of which may be novel drivers or markers of islet development. For example, the leptin receptor gene, LEPR, was most highly expressed in published data from in vivo-matured cells compared to our endocrine pancreas-like cells (log2 FC = 5.5, p-value = 2.0 × 10−12), suggesting a role for the leptin pathway in the maturation process. Endocrine pancreas-like cells showed significant stage-selective expression of adult islet genes, including INS, ABCC8, and GLP1R, and enrichment of relevant GO-terms (e.g. “insulin secretion”; odds ratio = 4.2, p-value = 1.9 × 10−3): however, principal component analysis indicated that in vitro-differentiated cells were more immature than adult islets. Integration of the stage-specific expression information with genetic data from T2D genome-wide association studies revealed that 46 of 82 T2D-associated loci harbor genes present in at least one developmental stage, facilitating refinement of potential effector transcripts. Together, these data show that expression profiling in an iPSC islet development model can further understanding of islet biology and T2D pathogenesis.

Insights Into the Pathogenicity of Rare Missense GCK Variants From the Identification and Functional Characterization of Compound Heterozygous and Double Mutations Inherited in Cis

OBJECTIVE To demonstrate the importance of using a combined genetic and functional approach to correctly interpret a genetic test for monogenic diabetes. RESEARCH DESIGN AND METHODS We identified three probands with a phenotype consistent with maturity-onset diabetes of the young (MODY) subtype GCK-MODY, in whom two potential pathogenic mutations were identified: [R43H/G68D], [E248 K/I225M], or [G261R/D217N]. Allele-specific PCR and cosegregation were used to determine phase. Single and double mutations were kinetically characterized. RESULTS The mutations occurred in cis (double mutants) in two probands and in trans in one proband. Functional studies of all double mutants revealed inactivating kinetics. The previously reported GCK-MODY mutations R43H and G68D were inherited from an affected father and unaffected mother, respectively. Both our functional and genetic studies support R43H as the cause of GCK-MODY and G68D as a neutral rare variant. CONCLUSIONS These data highlight the need for family/functional studies, even for previously reported pathogenic mutations.