Variations in the G6PC2/ABCB11 genomic region are associated with fasting glucose levels

Multiple promoter and enhancer differences likely contribute to augmented G6PC2 expression in human versus mouse pancreatic islet alpha cells

G6PC2 encodes a glucose-6-phosphatase catalytic subunit that opposes the action of glucokinase in pancreatic islets, thereby modulating the sensitivity of insulin and glucagon secretion to glucose. In mice, G6pc2 is expressed at ~20-fold higher levels in β-cells than in α-cells, whereas in humans G6PC2 is expressed at only ~5-fold higher levels in β-cells. We therefore hypothesize that G6PC2 likely influences glucagon secretion to a greater degree in humans. With a view to generating a humanized mouse that recapitulates augmented G6PC2 expression levels in α-cells, we sought to identify the genomic regions that confer differential mouse G6pc2 expression in α-cells versus β-cells as well as the evolutionary changes that have altered this ratio in humans. Studies in islet-derived cell lines suggest that the elevated G6pc2 expression in mouse β-cells versus α-cells is mainly due to a difference in the relative activity of the proximal G6pc2 promoter in these cell types. Similarly, the smaller difference in G6PC2 expression between α-cells and β-cells in humans is potentially explained by a change in relative proximal G6PC2 promoter activity. However, we show that both glucocorticoid levels and multiple differences in the relative activity of eight transcriptional enhancers between mice and humans likely contribute to differential G6PC2 expression. Finally, we show that a mouse-specific non-coding RNA, Gm13613, whose expression is controlled by G6pc2 enhancer I, does not regulate G6pc2 expression, indicating that altered expression of Gm13613 in a humanized mouse that contains both the human promoter and enhancers should not affect G6PC2 function.

Identification of structural motifs critical for human G6PC2 function informed by sequence analysis and an AlphaFold2-predicted model

G6PC2 encodes a glucose-6-phosphatase (G6Pase) catalytic subunit, primarily expressed in pancreatic islet β cells, which modulates the sensitivity of insulin secretion to glucose and thereby regulates fasting blood glucose (FBG). Mutational analyses were conducted to validate an AlphaFold2 (AF2)-predicted structure of human G6PC2 in conjunction with a novel method to solubilize and purify human G6PC2 from a heterologous expression system. These analyses show that residues forming a predicted intramolecular disulfide bond are essential for G6PC2 expression and that residues forming part of a type 2 phosphatidic acid phosphatase (PAP2) motif are critical for enzyme activity. Additional mutagenesis shows that residues forming a predicted substrate cavity modulate enzyme activity and substrate specificity and residues forming a putative cholesterol recognition amino acid consensus (CRAC) motif influence protein expression or enzyme activity. This CRAC motif begins at residue 219, the site of a common G6PC2 non-synonymous single-nucleotide polymorphism (SNP), rs492594 (Val219Leu), though the functional impact of this SNP is disputed. In microsomal membrane preparations, the L219 variant has greater activity than the V219 variant, but this difference disappears when G6PC2 is purified in detergent micelles. We hypothesize that this was due to a differential association of the two variants with cholesterol. This concept was supported by the observation that the addition of cholesteryl hemi-succinate to the purified enzymes decreased the Vmax of the V219 and L219 variants ∼8-fold and ∼3 fold, respectively. We anticipate that these observations should support the rational development of G6PC2 inhibitors designed to lower FBG.

An Enhancer Within Abcb11 Regulates G6pc2 in C57BL/6 Mouse Pancreatic Islets

G6PC2 is predominantly expressed in pancreatic islet β-cells where it encodes a glucose-6-phosphatase catalytic subunit that modulates the sensitivity of insulin secretion to glucose by opposing the action of glucokinase, thereby regulating fasting blood glucose (FBG). Prior studies have shown that the G6pc2 promoter alone is unable to confer sustained islet-specific gene expression in mice, suggesting the existence of distal enhancers that regulate G6pc2 expression. Using information from both mice and humans and knowledge that single nucleotide polymorphisms (SNPs) both within and near G6PC2 are associated with variations in FBG in humans, we identified several putative enhancers 3' of G6pc2. One region, herein referred to as enhancer I, resides in the 25th intron of Abcb11 and binds multiple islet-enriched transcription factors. CRISPR-mediated deletion of enhancer I in C57BL/6 mice had selective effects on the expression of genes near the G6pc2 locus. In isolated islets, G6pc2 and Spc25 expression were reduced ∼50%, and Gm13613 expression was abolished, whereas Cers6 and nostrin expression were unaffected. This partial reduction in G6pc2 expression enhanced islet insulin secretion at basal glucose concentrations but did not affect FBG or glucose tolerance in vivo, consistent with the absence of a phenotype in G6pc2 heterozygous C57BL/6 mice.

ARTICLE HIGHLIGHTS

Large-scale exome array summary statistics resources for glycemic traits to aid effector gene prioritization

Background

Genome-wide association studies for glycemic traits have identified hundreds of loci associated with these biomarkers of glucose homeostasis. Despite this success, the challenge remains to link variant associations to genes, and underlying biological pathways.

Methods

To identify coding variant associations which may pinpoint effector genes at both novel and previously established genome-wide association loci, we performed meta-analyses of exome-array studies for four glycemic traits: glycated hemoglobin (HbA1c, up to 144,060 participants), fasting glucose (FG, up to 129,665 participants), fasting insulin (FI, up to 104,140) and 2hr glucose post-oral glucose challenge (2hGlu, up to 57,878). In addition, we performed network and pathway analyses.

Results

Single-variant and gene-based association analyses identified coding variant associations at more than 60 genes, which when combined with other datasets may be useful to nominate effector genes. Network and pathway analyses identified pathways related to insulin secretion, zinc transport and fatty acid metabolism. HbA1c associations were strongly enriched in pathways related to blood cell biology.

Conclusions

Our results provided novel glycemic trait associations and highlighted pathways implicated in glycemic regulation. Exome-array summary statistic results are being made available to the scientific community to enable further discoveries.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic hESC-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined β-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-β cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.

Type 2 Diabetes-Related Variants Influence the Risk of Developing Prostate Cancer: A Population-Based Case-Control Study and Meta-Analysis

In this study, we have evaluated whether 57 genome-wide association studies (GWAS)-identified common variants for type 2 diabetes (T2D) influence the risk of developing prostate cancer (PCa) in a population of 304 Caucasian PCa patients and 686 controls. The association of selected single nucleotide polymorphisms (SNPs) with the risk of PCa was validated through meta-analysis of our data with those from the UKBiobank and FinnGen cohorts, but also previously published genetic studies. We also evaluated whether T2D SNPs associated with PCa risk could influence host immune responses by analysing their correlation with absolute numbers of 91 blood-derived cell populations and circulating levels of 103 immunological proteins and 7 steroid hormones. We also investigated the correlation of the most interesting SNPs with cytokine levels after in vitro stimulation of whole blood, peripheral mononuclear cells (PBMCs), and monocyte-derived macrophages with LPS, PHA, Pam3Cys, and Staphylococcus Aureus. The meta-analysis of our data with those from six large cohorts confirmed that each copy of the FTO_rs9939609A, HNF1B_rs7501939T, HNF1B_rs757210T, HNF1B_rs4430796G, and JAZF1_rs10486567A alleles significantly decreased risk of developing PCa (p = 3.70 × 10^-5, p = 9.39 × 10^-54, p = 5.04 × 10^-54, p = 1.19 × 10^-71, and p = 1.66 × 10^-18, respectively). Although it was not statistically significant after correction for multiple testing, we also found that the NOTCH2_rs10923931T and RBMS1_rs7593730 SNPs associated with the risk of developing PCa (p = 8.49 × 10^-4 and 0.004). Interestingly, we found that the protective effect attributed to the HFN1B locus could be mediated by the SULT1A1 protein (p = 0.00030), an arylsulfotransferase that catalyzes the sulfate conjugation of many hormones, neurotransmitters, drugs, and xenobiotic compounds. In addition to these results, eQTL analysis revealed that the HNF1B_rs7501939, HNF1B_rs757210, HNF1B_rs4430796, NOTCH2_rs10923931, and RBMS1_rs7593730 SNPs influence the risk of PCa through the modulation of mRNA levels of their respective genes in whole blood and/or liver. These results confirm that functional TD2-related variants influence the risk of developing PCa, but also highlight the need of additional experiments to validate our functional results in a tumoral tissue context.

Multi-omics analysis identifies CpGs near G6PC2 mediating the effects of genetic variants on fasting glucose

AIMS/HYPOTHESIS

An elevated fasting glucose level in non-diabetic individuals is a key predictor of type 2 diabetes. Genome-wide association studies (GWAS) have identified hundreds of SNPs for fasting glucose but most of their functional roles in influencing the trait are unclear. This study aimed to identify the mediation effects of DNA methylation between SNPs identified as significant from GWAS and fasting glucose using Mendelian randomisation (MR) analyses.

METHODS

We first performed GWAS analyses for three cohorts (Taiwan Biobank with 18,122 individuals, the Healthy Aging Longitudinal Study in Taiwan with 1989 individuals and the Stanford Asia-Pacific Program for Hypertension and Insulin Resistance with 416 individuals) with individuals of Han Chinese ancestry in Taiwan, followed by a meta-analysis for combining the three GWAS analysis results to identify significant and independent SNPs for fasting glucose. We determined whether these SNPs were methylation quantitative trait loci (meQTLs) by testing their associations with DNA methylation levels at nearby CpG sites using a subsample of 1775 individuals from the Taiwan Biobank. The MR analysis was performed to identify DNA methylation with causal effects on fasting glucose using meQTLs as instrumental variables based on the 1775 individuals. We also used a two-sample MR strategy to perform replication analysis for CpG sites with significant MR effects based on literature data.

RESULTS

Our meta-analysis identified 18 significant (p < 5 × 10^-8) and independent SNPs for fasting glucose. Interestingly, all 18 SNPs were meQTLs. The MR analysis identified seven CpGs near the G6PC2 gene that mediated the effects of a significant SNP (rs2232326) in the gene on fasting glucose. The MR effects for two CpGs were replicated using summary data based on the European population, using an exonic SNP rs2232328 in G6PC2 as the instrument.

CONCLUSIONS/INTERPRETATION

Our analysis results suggest that rs2232326 and rs2232328 in G6PC2 may affect DNA methylation at CpGs near the gene and that the methylation may have downstream effects on fasting glucose. Therefore, SNPs in G6PC2 and CpGs near G6PC2 may reside along the pathway that influences fasting glucose levels. This is the first study to report CpGs near G6PC2, an important gene for regulating insulin secretion, mediating the effects of GWAS-significant SNPs on fasting glucose.

Understanding the genetic architecture of the metabolically unhealthy normal weight and metabolically healthy obese phenotypes in a Korean population

Understanding the mechanisms underlying the metabolically unhealthy normal weight (MUHNW) and metabolically healthy obese (MHO) phenotypes is important for developing strategies to prevent cardiometabolic diseases. Here, we conducted genome-wide association studies (GWASs) to identify the MUHNW and MHO genetic indices. The study dataset comprised genome-wide single-nucleotide polymorphism genotypes and epidemiological data from 49,915 subjects categorised into four phenotypes-metabolically healthy normal weight (MHNW), MUHNW, MHO, and metabolically unhealthy obese (MUHO). We conducted two GWASs using logistic regression analyses and adjustments for confounding variables (model 1: MHNW versus MUHNW and model 2: MHO versus MUHO). GCKR, ABCB11, CDKAL1, LPL, CDKN2B, NT5C2, APOA5, CETP, and APOC1 were associated with metabolically unhealthy phenotypes among normal weight individuals (model 1). LPL, APOA5, and CETP were associated with metabolically unhealthy phenotypes among obese individuals (model 2). The genes common to both models are related to lipid metabolism (LPL, APOA5, and CETP), and those associated with model 1 are related to insulin or glucose metabolism (GCKR, CDKAL1, and CDKN2B). This study reveals the genetic architecture of the MUHNW and MHO phenotypes in a Korean population-based cohort. These findings could help identify individuals at a high metabolic risk in normal weight and obese populations and provide potential novel targets for the management of metabolically unhealthy phenotypes.

G6PC2 confers protection against hypoglycemia upon ketogenic diet feeding and prolonged fasting

Objective

G6PC2 is predominantly expressed in pancreatic islet beta cells. G6PC2 hydrolyzes glucose-6-phosphate to glucose and inorganic phosphate, thereby creating a futile substrate cycle that opposes the action of glucokinase. This substrate cycle determines the sensitivity of glucose-stimulated insulin secretion to glucose and hence regulates fasting blood glucose (FBG) but not fasting plasma insulin (FPI) levels. Our objective was to explore the physiological benefit this cycle confers.

Methods

We investigated the response of wild type (WT) and G6pc2 knockout (KO) mice to changes in nutrition.

Results

Pancreatic G6pc2 expression was little changed by ketogenic diet feeding but was inhibited by 24 hr fasting and strongly induced by high fat feeding. When challenged with either a ketogenic diet or 24 hr fasting, blood glucose fell to 70 mg/dl or less in G6pc2 KO but not WT mice, suggesting that G6PC2 may have evolved, in part, to prevent hypoglycemia. Prolonged ketogenic diet feeding reduced the effect of G6pc2 deletion on FBG. The hyperglycemia associated with high fat feeding was partially blunted in G6pc2 KO mice, suggesting that under these conditions the presence of G6PC2 is detrimental. As expected, FPI changed but did not differ between WT and KO mice in response to fasting, ketogenic and high fat feeding.

Conclusions

Since elevated FBG levels are associated with increased risk for cardiovascular-associated mortality (CAM), these studies suggest that, while G6PC2 inhibitors would be useful for lowering FBG and the risk of CAM, partial inhibition will be important to avoid the risk of hypoglycemia.

•

G6pc2 deletion lowers fasting blood glucose (FBG) in chow and high fat fed mice.

•

Elevated FBG increases the risk of cardiovascular-associated mortality (CAM).

•

G6pc2 deletion results in hypoglycemia in mice on a ketogenic diet.

•

G6pc2 deletion results in hypoglycemia in mice following prolonged fasting.

•

G6PC2 inhibitors may prevent CAM but increase risk of hypoglycemia.

Pancreatic islet beta cell-specific deletion of G6pc2 reduces fasting blood glucose

The G6PC1, G6PC2 and G6PC3 genes encode distinct glucose-6-phosphatase catalytic subunit (G6PC) isoforms. In mice, germline deletion of G6pc2 lowers fasting blood glucose (FBG) without affecting fasting plasma insulin (FPI) while, in isolated islets, glucose-6-phosphatase activity and glucose cycling are abolished and glucose-stimulated insulin secretion (GSIS) is enhanced at submaximal but not high glucose. These observations are all consistent with a model in which G6PC2 regulates the sensitivity of GSIS to glucose by opposing the action of glucokinase. G6PC2 is highly expressed in human and mouse islet beta cells however, various studies have shown trace G6PC2 expression in multiple tissues raising the possibility that G6PC2 also affects FBG through non-islet cell actions. Using real-time PCR we show here that expression of G6pc1 and/or G6pc3 are much greater than G6pc2 in peripheral tissues, whereas G6pc2 expression is much higher than G6pc3 in both pancreas and islets with G6pc1 expression not detected. In adult mice, beta cell-specific deletion of G6pc2 was sufficient to reduce FBG without changing FPI. In addition, electronic health record-derived phenotype analyses showed no association between G6PC2 expression and phenotypes clearly unrelated to islet function in humans. Finally, we show that germline G6pc2 deletion enhances glycolysis in mouse islets and that glucose cycling can also be detected in human islets. These observations are all consistent with a mechanism by which G6PC2 action in islets is sufficient to regulate the sensitivity of GSIS to glucose and hence influence FBG without affecting FPI.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Obesity genetics and cardiometabolic health: Potential for risk prediction

The increasing burden of obesity worldwide and its effect on cardiovascular disease (CVD) risk is an opportunity for evaluation of preventive approaches. Both obesity and CVD have a genetic background and polymorphisms within genes which enhance expression of variant proteins that influence CVD in obesity. Genome-based prediction may therefore be a feasible strategy, but the identification of genetically driven risk factors for CVD manifesting as clinically recognized phenotypes is a major challenge. Clusters of such risk factors include hyperglycaemia, hypertension, ectopic liver fat, and inflammation. All involve multiple genetic pathways having complex interactions with variable environmental influences. The factors that make significant contributions to CVD risk include altered carbohydrate homeostasis, ectopic deposition of fat in muscle and liver, and inflammation, with contributions from the gut microbiome. A futuristic model depends on harnessing the predictive power of plausible genetic variants, phenotype reversibility, and effective therapeutic choices based on genotype-phenotype interactions. Inverting disease phenotypes into ideal cardiovascular health metrics could improve genetic and epigenetic assessment, and form the basis of a future model for risk detection and early intervention.

The Genetic Basis of Metabolic Disease

Recent developments in genetics and genomics are providing a detailed and systematic characterization of the genetic underpinnings of common metabolic diseases and traits, highlighting the inherent complexity within systems for homeostatic control and the many ways in which that control can fail. The genetic architecture underlying these common metabolic phenotypes is complex, with each trait influenced by hundreds of loci spanning a range of allele frequencies and effect sizes. Here, we review the growing appreciation of this complexity and how this has fostered the implementation of genome-scale approaches that deliver robust mechanistic inference and unveil new strategies for translational exploitation.

An in vivo screen for neuronal genes involved in obesity identifies Diacylglycerol kinase as a regulator of insulin secretion

Objective

Obesity is a complex disorder involving many genetic and environmental factors that are required to maintain energy homeostasis. While studies in human populations have led to significant progress in the generation of an obesity gene map and broadened our understanding of the genetic basis of common obesity, there is still a large portion of heritability and etiology that remains unknown. Here, we have used the genetically tractable fruit fly, Drosophila melanogaster, to identify genes/pathways that function in the nervous system to regulate energy balance.

Methods

We performed an in vivo RNAi screen in Drosophila neurons and assayed for obese or lean phenotypes by measuring changes in levels of stored fats (in the form of triacylglycerides or TAG). Three rounds of screening were performed to verify the reproducibility and specificity of the adiposity phenotypes. Genes that produced >25% increase in TAG (206 in total) underwent a second round of screening to verify their effect on TAG levels by retesting the same RNAi line to validate the phenotype. All remaining hits were screened a third time by testing the TAG levels of additional RNAi lines against the genes of interest to rule out any off-target effects.

Results

We identified 24 genes including 20 genes that have not been previously associated with energy homeostasis. One identified hit, Diacylglycerol kinase (Dgk), has mammalian homologues that have been implicated in genome-wide association studies for metabolic defects. Downregulation of neuronal Dgk levels increases TAG and carbohydrate levels and these phenotypes can be recapitulated by reducing Dgk levels specifically within the insulin-producing cells that secrete Drosophila insulin-like peptides (dILPs). Conversely, overexpression of kinase-dead Dgk, but not wild-type, decreased circulating dILP2 and dILP5 levels resulting in lower insulin signalling activity. Despite having higher circulating dILP levels, Dgk RNAi flies have decreased pathway activity suggesting that they are insulin-resistant.

Conclusion

Altogether, we have identified several genes that act within the CNS to regulate energy homeostasis. One of these, Dgk, acts within the insulin-producing cells to regulate the secretion of dILPs and energy homeostasis in Drosophila.

•

RNAi screen in neurons identifies 24 regulators of energy homeostasis.

•

One of the hits, Dgk, affects lipid and carbohydrate homeostasis.

•

Dgk acts within the IPCs to regulate dILP secretion and insulin signalling activity.

Identification and functional analysis of glycemic trait loci in the China Health and Nutrition Survey

To identify genetic contributions to type 2 diabetes (T2D) and related glycemic traits (fasting glucose, fasting insulin, and HbA1c), we conducted genome-wide association analyses (GWAS) in up to 7,178 Chinese subjects from nine provinces in the China Health and Nutrition Survey (CHNS). We examined patterns of population structure within CHNS and found that allele frequencies differed across provinces, consistent with genetic drift and population substructure. We further validated 32 previously described T2D- and glycemic trait-loci, including G6PC2 and SIX3-SIX2 associated with fasting glucose. At G6PC2, we replicated a known fasting glucose-associated variant (rs34177044) and identified a second signal (rs2232326), a low-frequency (4%), probably damaging missense variant (S324P). A variant within the lead fasting glucose-associated signal at SIX3-SIX2 co-localized with pancreatic islet expression quantitative trait loci (eQTL) for SIX3, SIX2, and three noncoding transcripts. To identify variants functionally responsible for the fasting glucose association at SIX3-SIX2, we tested five candidate variants for allelic differences in regulatory function. The rs12712928-C allele, associated with higher fasting glucose and lower transcript expression level, showed lower transcriptional activity in reporter assays and increased binding to GABP compared to the rs12712928-G, suggesting that rs12712928-C contributes to elevated fasting glucose levels by disrupting an islet enhancer, resulting in reduced gene expression. Taken together, these analyses identified multiple loci associated with glycemic traits across China, and suggest a regulatory mechanism at the SIX3-SIX2 fasting glucose GWAS locus.

A common variant within the HNF1B gene is associated with overall survival of multiple myeloma patients: results from the IMMEnSE consortium and meta-analysis

Diabetogenic single nucleotide polymorphisms (SNPs) have recently been associated with multiple myeloma (MM) risk but their impact on overall survival (OS) of MM patients has not been analysed yet. In order to investigate the impact of 58 GWAS-identified variants for type 2 diabetes (T2D) on OS of patients with MM, we analysed genotyping data of 936 MM patients collected by the International Multiple Myeloma rESEarch (IMMENSE) consortium and an independent set of 700 MM patients recruited by the University Clinic of Heidelberg. A meta-analysis of the cox regression results of the two sets showed that rs7501939 located in the HNF1B gene negatively impacted OS (HRRec= 1.44, 95% CI = 1.18-1.76, P = 0.0001). The meta-analysis also showed a noteworthy gender-specific association of the SLC30A8rs13266634 SNP with OS. The presence of each additional copy of the minor allele at rs13266634 was associated with poor OS in men whereas no association was seen in women (HRMen-Add = 1.32, 95% CI 1.13-1.54, P = 0.0003). In conclusion, these data suggest that the HNF1Brs7501939 SNP confers poor OS in patients with MM and that a SNP in SLC30A8 affect OS in men.

Genetic aetiology of glycaemic traits: approaches and insights

Glycaemic traits such as fasting and post-challenge glucose and insulin measures, as well as glycated haemoglobin (HbA1c), are used to diagnose and monitor diabetes. These traits are risk factors for cardiovascular disease even below the diabetic threshold, and their study can additionally yield insights into the pathophysiology of type 2 diabetes. To date, a diverse set of genetic approaches have led to the discovery of over 97 loci influencing glycaemic traits. In this review, we will focus on recent advances in the genetic aetiology of glycaemic traits, and the resulting biological insights. We will provide a brief overview of results ranging from common, to low- and rare-frequency variant-trait association studies, studies leveraging the diversity across populations, and studies harnessing the power of genetic and genomic approaches to gain insights into the biological underpinnings of these traits.

Meta-analyses of the association of G6PC2 allele variants with elevated fasting glucose and type 2 diabetes

Objective

To collectively evaluate the association of glucose-6-phosphatase catalytic unit 2 (G6PC2) allele variants with elevated fasting glucose (FG) and type 2 diabetes (T2D).

Design

Meta-analysis

Data sources

PubMed, Web of Knowledge and Embase databases.

Study selection

Full text articles of studies that identified an association of G6PC2 with T2D and elevated FG.

Patient involvement

There was no T2D patient involvement in the analyses on the association of FG with G6PC2, there were T2D patients and non-diabetes patient involvement in the analyses on the association of T2D with G6PC2.

Statistical analysis

Random-effects meta-analyses were used to calculate the pool effect sizes. I² metric and H² tests were used to calculate the heterogeneity. Begg's funnel plot and Egger’s linear regression test were done to assess publication bias.

Results

Of the 423 studies identified, 21 were eligible and included. Data on three loci (rs560887, rs16856187 and rs573225) were available. The G allele at rs560887 in three ethnicities, the C allele at rs16856187 and the A allele at rs573225 all had a positive association with elevated FG. Per increment of G allele at rs560887 and A allele at rs573225 resulted in a FG 0.070 mmol/l and 0.075 mmol/l higher (ß (95% CI) = 0.070 (0.060, 0.079), p = 4.635e-50 and 0.075 (0.065, 0.085), p = 5.856e-48, respectively). With regard to the relationship of rs16856187 and FG, an increase of 0.152 (95% CI: 0.034–0.270; p = 0.011) and 0.317 (95% CI: 0.193–0.442, p = 6.046e-07) was found in the standardized mean difference (SMD) of FG for the AC and CC genotypes, respectively, when compared with the AA reference genotype. However, the G-allele of rs560887 in Caucasians under the additive model and the C-allele of rs16856187 under the allele and dominant models were associated with a decreased risk of T2D (OR (95% CI) = 0.964 (0.947, 0.981), p = 0.570e-4; OR (95% CI) = 0.892 (0.832, 0.956), p = 0.001; and OR (95% CI) = 0.923(0.892, 0.955), p = 5.301e-6, respectively).

Conclusions

Our meta-analyses demonstrate that all three allele variants of G6PC2 (rs560887, rs16856187 and rs573225) are associated with elevated FG, with two variants (rs560887 in the Caucasians subgroup and rs16856187 under the allele and dominant model) being associated with T2D as well. Further studies utilizing larger sample sizes and different ethnic populations are needed to extend and confirm these findings.

Human genetics as a model for target validation: finding new therapies for diabetes

Type 2 diabetes is a global epidemic with major effects on healthcare expenditure and quality of life. Currently available treatments are inadequate for the prevention of comorbidities, yet progress towards new therapies remains slow. A major barrier is the insufficiency of traditional preclinical models for predicting drug efficacy and safety. Human genetics offers a complementary model to assess causal mechanisms for target validation. Genetic perturbations are 'experiments of nature' that provide a uniquely relevant window into the long-term effects of modulating specific targets. Here, we show that genetic discoveries over the past decades have accurately predicted (now known) therapeutic mechanisms for type 2 diabetes. These findings highlight the potential for use of human genetic variation for prospective target validation, and establish a framework for future applications. Studies into rare, monogenic forms of diabetes have also provided proof-of-principle for precision medicine, and the applicability of this paradigm to complex disease is discussed. Finally, we highlight some of the limitations that are relevant to the use of genome-wide association studies (GWAS) in the search for new therapies for diabetes. A key outstanding challenge is the translation of GWAS signals into disease biology and we outline possible solutions for tackling this experimental bottleneck.

Effects of G6pc2 deletion on body weight and cholesterol in mice

Genome-wide association study (GWAS) data have linked the G6PC2 gene to variations in fasting blood glucose (FBG). G6PC2 encodes an islet-specific glucose-6-phosphatase catalytic subunit that forms a substrate cycle with the beta cell glucose sensor glucokinase. This cycle modulates the glucose sensitivity of insulin secretion and hence FBG. GWAS data have not linked G6PC2 to variations in body weight but we previously reported that female C57BL/6J G6pc2-knockout (KO) mice were lighter than wild-type littermates on both a chow and high-fat diet. The purpose of this study was to compare the effects of G6pc2 deletion on FBG and body weight in both chow-fed and high-fat-fed mice on two other genetic backgrounds. FBG was reduced in G6pc2 KO mice largely independent of gender, genetic background or diet. In contrast, the effect of G6pc2 deletion on body weight was markedly influenced by these variables. Deletion of G6pc2 conferred a marked protection against diet-induced obesity in male mixed genetic background mice, whereas in 129SvEv mice deletion of G6pc2 had no effect on body weight. G6pc2 deletion also reduced plasma cholesterol levels in a manner dependent on gender, genetic background and diet. An association between G6PC2 and plasma cholesterol was also observed in humans through electronic health record-derived phenotype analyses. These observations suggest that the action of G6PC2 on FBG is largely independent of the influences of environment, modifier genes or epigenetic events, whereas the action of G6PC2 on body weight and cholesterol are influenced by unknown variables.

Susceptibility to type 2 diabetes may be modulated by haplotypes in G6PC2, a target of positive selection

BACKGROUND

The endoplasmic reticulum enzyme glucose-6-phosphatase catalyzes the common terminal reaction in the gluconeogenic/glycogenolytic pathways and plays a central role in glucose homeostasis. In most mammals, different G6PC subunits are encoded by three paralogous genes (G6PC, G6PC2, and G6PC3). Mutations in G6PC and G6PC3 are responsible for human mendelian diseases, whereas variants in G6PC2 are associated with fasting glucose (FG) levels.

RESULTS

We analyzed the evolutionary history of G6Pase genes. Results indicated that the three paralogs originated during early vertebrate evolution and that negative selection was the major force shaping diversity at these genes in mammals. Nonetheless, site-wise estimation of evolutionary rates at corresponding sites revealed weak correlations, suggesting that mammalian G6Pases have evolved different structural features over time. We also detected pervasive positive selection at mammalian G6PC2. Most selected residues localize in the C-terminal protein region, where several human variants associated with FG levels also map. This region was re-sequenced in ~560 subjects from Saudi Arabia, 185 of whom suffering from type 2 diabetes (T2D). The frequency of rare missense and nonsense variants was not significantly different in T2D and controls. Association analysis with two common missense variants (V219L and S342C) revealed a weak but significant association for both SNPs when analyses were conditioned on rs560887, previously identified in a GWAS for FG. Two haplotypes were significantly associated with T2D with an opposite effect direction.

CONCLUSIONS

We detected pervasive positive selection at mammalian G6PC2 genes and we suggest that distinct haplotypes at the G6PC2 locus modulate susceptibility to T2D.

The Metabolic Syndrome in Men study: a resource for studies of metabolic and cardiovascular diseases

The Metabolic Syndrome in Men (METSIM) study is a population-based study including 10,197 Finnish men examined in 2005-2010. The aim of the study is to investigate nongenetic and genetic factors associated with the risk of T2D and CVD, and with cardiovascular risk factors. The protocol includes a detailed phenotyping of the participants, an oral glucose tolerance test, fasting laboratory measurements including proton NMR measurements, mass spectometry metabolomics, adipose tissue biopsies from 1,400 participants, and a stool sample. In our ongoing follow-up study, we have, to date, reexamined 6,496 participants. Extensive genotyping and exome sequencing have been performed for essentially all METSIM participants, and >2,000 METSIM participants have been whole-genome sequenced. We have identified several nongenetic markers associated with the development of diabetes and cardiovascular events, and participated in several genetic association studies to identify gene variants associated with diabetes, hyperglycemia, and cardiovascular risk factors. The generation of a phenotype and genotype resource in the METSIM study allows us to proceed toward a "systems genetics" approach, which includes statistical methods to quantitate and integrate intermediate phenotypes, such as transcript, protein, or metabolite levels, to provide a global view of the molecular architecture of complex traits.

Single-cell transcriptomes identify human islet cell signatures and reveal cell-type-specific expression changes in type 2 diabetes

Blood glucose levels are tightly controlled by the coordinated action of at least four cell types constituting pancreatic islets. Changes in the proportion and/or function of these cells are associated with genetic and molecular pathophysiology of monogenic, type 1, and type 2 (T2D) diabetes. Cellular heterogeneity impedes precise understanding of the molecular components of each islet cell type that govern islet (dys)function, particularly the less abundant delta and gamma/pancreatic polypeptide (PP) cells. Here, we report single-cell transcriptomes for 638 cells from nondiabetic (ND) and T2D human islet samples. Analyses of ND single-cell transcriptomes identified distinct alpha, beta, delta, and PP/gamma cell-type signatures. Genes linked to rare and common forms of islet dysfunction and diabetes were expressed in the delta and PP/gamma cell types. Moreover, this study revealed that delta cells specifically express receptors that receive and coordinate systemic cues from the leptin, ghrelin, and dopamine signaling pathways implicating them as integrators of central and peripheral metabolic signals into the pancreatic islet. Finally, single-cell transcriptome profiling revealed genes differentially regulated between T2D and ND alpha, beta, and delta cells that were undetectable in paired whole islet analyses. This study thus identifies fundamental cell-type–specific features of pancreatic islet (dys)function and provides a critical resource for comprehensive understanding of islet biology and diabetes pathogenesis.

G6PC2 Modulates the Effects of Dexamethasone on Fasting Blood Glucose and Glucose Tolerance

The glucose-6-phosphatase catalytic subunit 2 (G6PC2) gene encodes an islet-specific glucose-6-phosphatase catalytic subunit. G6PC2 forms a substrate cycle with glucokinase that determines the glucose sensitivity of insulin secretion. Consequently, deletion of G6pc2 lowers fasting blood glucose (FBG) without affecting fasting plasma insulin. Although chronic elevation of FBG is detrimental to health, glucocorticoids induce G6PC2 expression, suggesting that G6PC2 evolved to transiently modulate FBG under conditions of glucocorticoid-related stress. We show, using competition and mutagenesis experiments, that the synthetic glucocorticoid dexamethasone (Dex) induces G6PC2 promoter activity through a mechanism involving displacement of the islet-enriched transcription factor MafA by the glucocorticoid receptor. The induction of G6PC2 promoter activity by Dex is modulated by a single nucleotide polymorphism, previously linked to altered FBG in humans, that affects FOXA2 binding. A 5-day repeated injection paradigm was used to examine the chronic effect of Dex on FBG and glucose tolerance in wild-type (WT) and G6pc2 knockout mice. Acute Dex treatment only induces G6pc2 expression in 129SvEv but not C57BL/6J mice, but this chronic treatment induced G6pc2 expression in both. In 6-hour fasted C57BL/6J WT mice, Dex treatment lowered FBG and improved glucose tolerance, with G6pc2 deletion exacerbating the decrease in FBG and enhancing the improvement in glucose tolerance. In contrast, in 24-hour fasted C57BL/6J WT mice, Dex treatment raised FBG but still improved glucose tolerance, with G6pc2 deletion limiting the increase in FBG and enhancing the improvement in glucose tolerance. These observations demonstrate that G6pc2 modulates the complex effects of Dex on both FBG and glucose tolerance.

Functional Analysis of Mouse G6pc1 Mutations Using a Novel In Situ Assay for Glucose-6-Phosphatase Activity and the Effect of Mutations in Conserved Human G6PC1/G6PC2 Amino Acids on G6PC2 Protein Expression

Elevated fasting blood glucose (FBG) has been associated with increased risk for development of type 2 diabetes. Single nucleotide polymorphisms (SNPs) in G6PC2 are the most important common determinants of variations in FBG in humans. Studies using G6pc2 knockout mice suggest that G6pc2 regulates the glucose sensitivity of insulin secretion. G6PC2 and the related G6PC1 and G6PC3 genes encode glucose-6-phosphatase catalytic subunits. This study describes a functional analysis of 22 non-synonymous G6PC2 SNPs, that alter amino acids that are conserved in human G6PC1, mouse G6pc1 and mouse G6pc2, with the goal of identifying variants that potentially affect G6PC2 activity/expression. Published data suggest strong conservation of catalytically important amino acids between all four proteins and the related G6PC3 isoform. Because human G6PC2 has very low glucose-6-phosphatase activity we used an indirect approach, examining the effect of these SNPs on mouse G6pc1 activity. Using a novel in situ functional assay for glucose-6-phosphatase activity we demonstrate that the amino acid changes associated with the human G6PC2 rs144254880 (Arg79Gln), rs149663725 (Gly114Arg) and rs2232326 (Ser324Pro) SNPs reduce mouse G6pc1 enzyme activity without affecting protein expression. The Arg79Gln variant alters an amino acid mutation of which, in G6PC1, has previously been shown to cause glycogen storage disease type 1a. We also demonstrate that the rs368382511 (Gly8Glu), rs138726309 (His177Tyr), rs2232323 (Tyr207Ser) rs374055555 (Arg293Trp), rs2232326 (Ser324Pro), rs137857125 (Pro313Leu) and rs2232327 (Pro340Leu) SNPs confer decreased G6PC2 protein expression. In summary, these studies identify multiple G6PC2 variants that have the potential to be associated with altered FBG in humans.

Type 2 diabetes-related variants influence the risk of developing multiple myeloma: results from the IMMEnSE consortium

Type 2 diabetes (T2D) has been suggested to be a risk factor for multiple myeloma (MM), but the relationship between the two traits is still not well understood. The aims of this study were to evaluate whether 58 genome-wide-association-studies (GWAS)-identified common variants for T2D influence the risk of developing MM and to determine whether predictive models built with these variants might help to predict the disease risk. We conducted a case-control study including 1420 MM patients and 1858 controls ascertained through the International Multiple Myeloma (IMMEnSE) consortium. Subjects carrying the KCNQ1rs2237892T allele or the CDKN2A-2Brs2383208G/G, IGF1rs35767T/T and MADDrs7944584T/T genotypes had a significantly increased risk of MM (odds ratio (OR)=1.32-2.13) whereas those carrying the KCNJ11rs5215C, KCNJ11rs5219T and THADArs7578597C alleles or the FTOrs8050136A/A and LTArs1041981C/C genotypes showed a significantly decreased risk of developing the disease (OR=0.76-0.85). Interestingly, a prediction model including those T2D-related variants associated with the risk of MM showed a significantly improved discriminatory ability to predict the disease when compared to a model without genetic information (area under the curve (AUC)=0.645 vs AUC=0.629; P=4.05×10(-) (06)). A gender-stratified analysis also revealed a significant gender effect modification for ADAM30rs2641348 and NOTCH2rs10923931 variants (Pinteraction=0.001 and 0.0004, respectively). Men carrying the ADAM30rs2641348C and NOTCH2rs10923931T alleles had a significantly decreased risk of MM whereas an opposite but not significant effect was observed in women (ORM=0.71 and ORM=0.66 vs ORW=1.22 and ORW=1.15, respectively). These results suggest that TD2-related variants may influence the risk of developing MM and their genotyping might help to improve MM risk prediction models.

Novel stable isotope analyses demonstrate significant rates of glucose cycling in mouse pancreatic islets

A polymorphism located in the G6PC2 gene, which encodes an islet-specific glucose-6-phosphatase catalytic subunit, is the most important common determinant of variations in fasting blood glucose (FBG) levels in humans. Studies of G6pc2 knockout (KO) mice suggest that G6pc2 represents a negative regulator of basal glucose-stimulated insulin secretion (GSIS) that acts by hydrolyzing glucose-6-phosphate (G6P), thereby reducing glycolytic flux. However, this conclusion conflicts with the very low estimates for the rate of glucose cycling in pancreatic islets, as assessed using radioisotopes. We have reassessed the rate of glucose cycling in pancreatic islets using a novel stable isotope method. The data show much higher levels of glucose cycling than previously reported. In 5 mmol/L glucose, islets from C57BL/6J chow-fed mice cycled ∼16% of net glucose uptake. The cycling rate was further increased at 11 mmol/L glucose. Similar cycling rates were observed using islets from high fat-fed mice. Importantly, glucose cycling was abolished in G6pc2 KO mouse islets, confirming that G6pc2 opposes the action of the glucose sensor glucokinase by hydrolyzing G6P. The demonstration of high rates of glucose cycling in pancreatic islets explains why G6pc2 deletion enhances GSIS and why variants in G6PC2 affect FBG in humans.

Genetic modifiers of response to glucose-insulin-potassium (GIK) infusion in acute coronary syndromes and associations with clinical outcomes in the IMMEDIATE trial

Modifiers of response to glucose, insulin and potassium (GIK) infusion may affect clinical outcomes in acute coronary syndromes (ACS). In an Immediate Myocardial Metabolic Enhancement During Initial Assessment And Treatment In Emergency Care (IMMEDIATE) trial's sub-study (n = 318), we explored effects of 132,634 genetic variants on plasma glucose and potassium response to 12-h GIK infusion. Associations between metabolite-associated variants and infarct size (n = 84) were assessed. The 'G' allele of rs12641551, near ACSL1, as well as the 'A' allele of XPO4 rs2585897 were associated with a differential glucose response (P for 2 degrees of freedom test, P2df ⩽ 4.75 × 10(-7)) and infarct size with GIK (P2df < 0.05). Variants within or near TAS1R3, LCA5, DNAH5, PTPRG, MAGI1, PTCSC3, STRADA, AKAP12, ARFGEF2, ADCYAP1, SETX, NDRG4 and ABCB11 modified glucose response, and near CSF1/AHCYL1 potassium response (P2df ⩽ 4.26 × 10(-7)), but not outcomes. Gene variants may modify glucose and potassium response to GIK infusion, contributing to cardiovascular outcomes in ACS.

Low-frequency and rare exome chip variants associate with fasting glucose and type 2 diabetes susceptibility

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol l(-1), P=3.4 × 10(-12)), T2D risk (OR[95%CI]=0.86[0.76-0.96], P=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose(-1), P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l(-1), P=4.3 × 10(-4)). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8 × 10(-6)) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l(-1), P=1.3 × 10(-8)). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

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.

Understanding how FI and FG levels are regulated is important because their derangement is a feature of T2D. Despite recent success from GWAS in identifying regions of the genome influencing glycemic traits, collectively these loci explain only a small proportion of trait variance. Unlocking the biological mechanisms driving these associations has been challenging because the vast majority of variants map to non-coding sequence, and the genes through which they exert their impact are largely unknown. In the current study, we sought to increase our understanding of the physiological pathways influencing both traits using exome-array genotyping in up to 33,231 non-diabetic individuals to identify coding variants and consequently genes associated with either FG or FI levels. We identified novel association signals for both traits including the receptor for GLP-1 agonists which are a widely used therapy for T2D. Furthermore, we identified coding variants at several GWAS loci which point to the genes underlying these association signals. Importantly, we found that multiple coding variants in G6PC2 result in a loss of protein function and lower fasting glucose levels.

Genome-wide association study identifies two novel Loci with sex-specific effects for type 2 diabetes mellitus and glycemic traits in a korean population

BACKGROUND

Until recently, genome-wide association study (GWAS)-based findings have provided a substantial genetic contribution to type 2 diabetes mellitus (T2DM) or related glycemic traits. However, identification of allelic heterogeneity and population-specific genetic variants under consideration of potential confounding factors will be very valuable for clinical applicability. To identify novel susceptibility loci for T2DM and glycemic traits, we performed a two-stage genetic association study in a Korean population.

METHODS

We performed a logistic analysis for T2DM, and the first discovery GWAS was analyzed for 1,042 cases and 2,943 controls recruited from a population-based cohort (KARE, n=8,842). The second stage, de novo replication analysis, was performed in 1,216 cases and 1,352 controls selected from an independent population-based cohort (Health 2, n=8,500). A multiple linear regression analysis for glycemic traits was further performed in a total of 14,232 nondiabetic individuals consisting of 7,696 GWAS and 6,536 replication study participants. A meta-analysis was performed on the combined results using effect size and standard errors estimated for stage 1 and 2, respectively.

RESULTS

A combined meta-analysis for T2DM identified two new (rs11065756 and rs2074356) loci reaching genome-wide significance in CCDC63 and C12orf51 on the 12q24 region. In addition, these variants were significantly associated with fasting plasma glucose and homeostasis model assessment of β-cell function. Interestingly, two independent single nucleotide polymorphisms were associated with sex-specific stratification in this study.

CONCLUSION

Our study showed a strong association between T2DM and glycemic traits. We further observed that two novel loci with multiple diverse effects were highly specific to males. Taken together, these findings may provide additional insights into the clinical assessment or subclassification of disease risk in a Korean population.

Insights into the genetic susceptibility to type 2 diabetes from genome-wide association studies of glycaemic traits

Over the past 8 years, the genetics of complex traits have benefited from an unprecedented advancement in the identification of common variant loci for diseases such as type 2 diabetes (T2D). The ability to undertake genome-wide association studies in large population-based samples for quantitative glycaemic traits has permitted us to explore the hypothesis that models arising from studies in non-diabetic individuals may reflect mechanisms involved in the pathogenesis of diabetes. Amongst 88 T2D risk and 72 glycaemic trait loci, only 29 are shared and show disproportionate magnitudes of phenotypic effects. Important mechanistic insights have been gained regarding the physiological role of T2D loci in disease predisposition through the elucidation of their contribution to glycaemic trait variability. Further investigation is warranted to define causal variants within these loci, including functional characterisation of associated variants, to dissect their role in disease mechanisms and to enable clinical translation.

Moving on from GWAS: functional studies on the G6PC2 gene implicated in the regulation of fasting blood glucose

Genome-wide association studies (GWAS) have shown that single-nucleotide polymorphisms (SNPs) in G6PC2 are the most important common determinants of variations in fasting blood glucose (FBG) levels. Molecular studies examining the functional impact of these SNPs on G6PC2 gene transcription and splicing suggest that they affect FBG by directly modulating G6PC2 expression. This conclusion is supported by studies on G6pc2 knockout (KO) mice showing that G6pc2 represents a negative regulator of basal glucose-stimulated insulin secretion that acts by hydrolyzing glucose-6-phosphate, thereby reducing glycolytic flux and opposing the action of glucokinase. Suppression of G6PC2 activity might, therefore, represent a novel therapy for lowering FBG and the risk of cardiovascular-associated mortality. GWAS and G6pc2 KO mouse studies also suggest that G6PC2 affects other aspects of beta cell function. The evolutionary benefit conferred by G6PC2 remains unclear, but it is unlikely to be related to its ability to modulate FBG.

Identification of HKDC1 and BACE2 as genes influencing glycemic traits during pregnancy through genome-wide association studies

Maternal metabolism during pregnancy impacts the developing fetus, affecting offspring birth weight and adiposity. This has important implications for metabolic health later in life (e.g., offspring of mothers with pre-existing or gestational diabetes mellitus have an increased risk of metabolic disorders in childhood). To identify genetic loci associated with measures of maternal metabolism obtained during an oral glucose tolerance test at ∼28 weeks' gestation, we performed a genome-wide association study of 4,437 pregnant mothers of European (n = 1,367), Thai (n = 1,178), Afro-Caribbean (n = 1,075), and Hispanic (n = 817) ancestry, along with replication of top signals in three additional European ancestry cohorts. In addition to identifying associations with genes previously implicated with measures of glucose metabolism in nonpregnant populations, we identified two novel genome-wide significant associations: 2-h plasma glucose and HKDC1, and fasting C-peptide and BACE2. These results suggest that the genetic architecture underlying glucose metabolism may differ, in part, in pregnancy.

A genome-wide association study of chronic otitis media with effusion and recurrent otitis media identifies a novel susceptibility locus on chromosome 2

Chronic otitis media with effusion (COME) and recurrent otitis media (ROM) have been shown to be heritable, but candidate gene and linkage studies to date have been equivocal. Our aim was to identify genetic susceptibility factors using a genome-wide association study (GWAS). We genotyped 602 subjects from 143 families with 373 COME/ROM subjects using the Illumina Human CNV370-Duo DNA Bead Chip (324,748 SNPs). We carried out the GWAS scan and imputed SNPs at the regions with the most significant associations. Replication genotyping in an independent family-based sample was conducted for 53 SNPs: the 41 most significant SNPs with P < 10(-4) and 12 imputed SNPs with P < 10(-4) on chromosome 15 (near the strongest signal). We replicated the association of rs10497394 (GWAS discovery P = 1.30 × 10(-5)) on chromosome 2 in the independent otitis media population (P = 4.7 × 10(-5); meta-analysis P = 1.52 × 10(-8)). Three additional SNPs had replication P values < 0.10. Two were on chromosome 15q26.1 including rs1110060, the strongest association with COME/ROM in the primary GWAS (P = 3.4 ×10(-7)) in KIF7 intron 7 (P = 0.072), and rs10775247, a non-synonymous SNP in TICRR exon 2 (P = 0.075). The third SNP rs386057 was on chromosome 5 in TPPP intron 1 (P = 0.045). We have performed the first GWAS of COME/ROM and have identified a SNP rs10497394 on chromosome 2 is significantly associated with COME/ROM susceptibility. This SNP is within a 537 kb intergenic region, bordered by CDCA7 and SP3. The genomic and functional significance of this newly identified locus in COME/ROM pathogenesis requires additional investigation.

Multiple functional polymorphisms in the G6PC2 gene contribute to the association with higher fasting plasma glucose levels

AIMS/HYPOTHESIS

We previously identified the G6PC2 locus as a strong determinant of fasting plasma glucose (FPG) and showed that a common G6PC2 intronic single nucleotide polymorphism (SNP) (rs560887) and two common G6PC2 promoter SNPs (rs573225 and rs13431652) are highly associated with FPG. However, these promoter SNPs have complex effects on G6PC2 fusion gene expression, and our data suggested that only rs13431652 is a potentially causative SNP. Here we examine the effect of rs560887 on G6PC2 pre-mRNA splicing and the contribution of an additional common G6PC2 promoter SNP, rs2232316, to the association signal.

METHODS

Minigene analyses were used to characterise the effect of rs560887 on G6PC2 pre-mRNA splicing. Fusion gene and gel retardation analyses characterised the effect of rs2232316 on G6PC2 promoter activity and transcription factor binding. The genetic association of rs2232316 with FPG variation was assessed using regression adjusted for age, sex and BMI in 4,220 Europeans with normal FPG.

RESULTS

The rs560887-G allele was shown to enhance G6PC2 pre-mRNA splicing, whereas the rs2232316-A allele enhanced G6PC2 transcription by promoting Foxa2 binding. Genetic analyses provide evidence for association of the rs2232316-A allele with increased FPG (β = 0.04 mmol/l; p = 4.3 × 10(-3)) as part of the same signal as rs560887, rs573225 and rs13431652.

CONCLUSIONS/INTERPRETATION

As with rs13431652, the in situ functional data with rs560887 and rs2232316 are in accord with the putative function of G6PC2 in pancreatic islets, and suggest that all three are potentially causative SNPs that contribute to the association between G6PC2 and FPG.

G6PC2: a negative regulator of basal glucose-stimulated insulin secretion

Elevated fasting blood glucose (FBG) is associated with increased risk for the development of type 2 diabetes and cardiovascular-associated mortality. Genome-wide association studies (GWAS) have linked polymorphisms in G6PC2 with variations in FBG and body fat, although not insulin sensitivity or glucose tolerance. G6PC2 encodes an islet-specific, endoplasmic reticulum-resident glucose-6-phosphatase catalytic subunit. A combination of in situ perfused pancreas, in vitro isolated islet, and in vivo analyses were used to explore the function of G6pc2 in mice. G6pc2 deletion had little effect on insulin sensitivity and glucose tolerance, whereas body fat was reduced in female G6pc2 knockout (KO) mice on both a chow and high-fat diet, observations that are all consistent with human GWAS data. G6pc2 deletion resulted in a leftward shift in the dose-response curve for glucose-stimulated insulin secretion (GSIS). As a consequence, under fasting conditions in which plasma insulin levels were identical, blood glucose levels were reduced in G6pc2 KO mice, again consistent with human GWAS data. Glucose-6-phosphatase activity was reduced, whereas basal cytoplasmic calcium levels were elevated in islets isolated from G6pc2 KO mice. These data suggest that G6pc2 represents a novel, negative regulator of basal GSIS that acts by hydrolyzing glucose-6-phosphate, thereby reducing glycolytic flux.

Combined genome-wide linkage and association analyses of fasting glucose level in healthy twins and families of Korea

This study was undertaken to identify genetic polymorphisms that are associated with the risk of an elevated fasting glucose (FG) level using genome-wide analyses. We explored a quantitative trait locus (QTL) for FG level in a genome-wide study from a Korean twin-family cohort (the Healthy Twin Study) using a combined linkage and family-based association analysis approach. We investigated 1,754 individuals, which included 432 families and 219 pairs of monozygotic twins. Regions of chromosomes 2q23.3-2q31.1, 15q26.1-15q26.3, 16p12.1, and 20p13-20p12.2, were found to show evidence of linkage with FG level, and several markers in these regions were found to be significantly associated with FG level using family-based or general association tests. In particular, a single-nucleotide polymorphism (rs6138953) on the PTPRA gene in the 20p13 region (combined P = 1.8 × 10(-6)) was found to be associated with FG level, and the PRKCB1 gene (in 16p12.1) to be possibly associated with FG level. In conclusion, multiple regions of chromosomes 2q23.3-2q31.1, 15q26.1-15q26.3, 16p12.1, and 20p13-20p12.2 are associated with FG level in our Korean twin-family cohort. The combined approach of genome-wide linkage and family-based association analysis is useful to identify novel or known genetic regions concerning FG level in a family cohort study.

A Drosophila functional evaluation of candidates from human genome-wide association studies of type 2 diabetes and related metabolic traits identifies tissue-specific roles for dHHEX

Background

Genome-wide association studies (GWAS) identify regions of the genome that are associated with particular traits, but do not typically identify specific causative genetic elements. For example, while a large number of single nucleotide polymorphisms associated with type 2 diabetes (T2D) and related traits have been identified by human GWAS, only a few genes have functional evidence to support or to rule out a role in cellular metabolism or dietary interactions. Here, we use a recently developed Drosophila model in which high-sucrose feeding induces phenotypes similar to T2D to assess orthologs of human GWAS-identified candidate genes for risk of T2D and related traits.

Results

Disrupting orthologs of certain T2D candidate genes (HHEX, THADA, PPARG, KCNJ11) led to sucrose-dependent toxicity. Tissue-specific knockdown of the HHEX ortholog dHHEX (CG7056) directed metabolic defects and enhanced lethality; for example, fat-body-specific loss of dHHEX led to increased hemolymph glucose and reduced insulin sensitivity.

Conclusion

Candidate genes identified in human genetic studies of metabolic traits can be prioritized and functionally characterized using a simple Drosophila approach. To our knowledge, this is the first large-scale effort to study the functional interaction between GWAS-identified candidate genes and an environmental risk factor such as diet in a model organism system.

Genetic determinants of cardiometabolic risk: a proposed model for phenotype association and interaction

This review provides a translational and unifying summary of metabolic syndrome genetics and highlights evidence that genetic studies are starting to unravel and untangle origins of the complex and challenging cluster of disease phenotypes. The associated genes effectively express in the brain, liver, kidney, arterial endothelium, adipocytes, myocytes, and β cells. Progression of syndrome traits has been associated with ectopic lipid accumulation in the arterial wall, visceral adipocytes, myocytes, and liver. Thus, it follows that the genetics of dyslipidemia, obesity, and nonalcoholic fatty liver disease are central in triggering progression of the syndrome to overt expression of disease traits and have become a key focus of interest for early detection and for designing prevention and treatments. To support the "birds' eye view" approach, we provide a road-map depicting commonality and interrelationships between the traits and their genetic and environmental determinants based on known risk factors, metabolic pathways, pharmacologic targets, treatment responses, gene networks, pleiotropy, and association with circadian rhythm. Although only a small portion of the known heritability is accounted for and there is insufficient support for clinical application of gene-based prediction models, there is direction and encouraging progress in a rapidly moving field that is beginning to show clinical relevance.

The hydrodynamic tail vein assay as a tool for the study of liver promoters and enhancers

The hydrodynamic tail vein injection is a technique that is used to deliver nucleic acids into live mice. Delivery through this method results in the in vivo transfection of foreign DNA primarily in the liver. Here, we describe the use of this technique to test for regulatory activity of liver promoters and enhancers, using a dual luciferase reporter system as the readable/measureable output and how this application can be used for pharmacogenomic studies.

Effects of genetic variants previously associated with fasting glucose and insulin in the Diabetes Prevention Program

Common genetic variants have been recently associated with fasting glucose and insulin levels in white populations. Whether these associations replicate in pre-diabetes is not known. We extended these findings to the Diabetes Prevention Program, a clinical trial in which participants at high risk for diabetes were randomized to placebo, lifestyle modification or metformin for diabetes prevention. We genotyped previously reported polymorphisms (or their proxies) in/near G6PC2, MTNR1B, GCK, DGKB, GCKR, ADCY5, MADD, CRY2, ADRA2A, FADS1, PROX1, SLC2A2, GLIS3, C2CD4B, IGF1, and IRS1 in 3,548 Diabetes Prevention Program participants. We analyzed variants for association with baseline glycemic traits, incident diabetes and their interaction with response to metformin or lifestyle intervention. We replicated associations with fasting glucose at MTNR1B (P<0.001), G6PC2 (P = 0.002) and GCKR (P = 0.001). We noted impaired β-cell function in carriers of glucose-raising alleles at MTNR1B (P<0.001), and an increase in the insulinogenic index for the glucose-raising allele at G6PC2 (P<0.001). The association of MTNR1B with fasting glucose and impaired β-cell function persisted at 1 year despite adjustment for the baseline trait, indicating a sustained deleterious effect at this locus. We also replicated the association of MADD with fasting proinsulin levels (P<0.001). We detected no significant impact of these variants on diabetes incidence or interaction with preventive interventions. The association of several polymorphisms with quantitative glycemic traits is replicated in a cohort of high-risk persons. These variants do not have a detectable impact on diabetes incidence or response to metformin or lifestyle modification in the Diabetes Prevention Program.

Abcb11 deficiency induces cholestasis coupled to impaired β-fatty acid oxidation in mice

The bile salt export pump (BSEP) is an ATP-binding cassette transporter that serves as the primary system for removing bile salts from the liver. In humans, deficiency of BSEP, which is encoded by the ABCB11 gene, causes severe progressive cholestatic liver disease from early infancy. In previous studies of Abcb11 deficiency in mice generated on a mixed genetic background, the animals did not recapitulate the human disease. We reasoned that ABCB11 deficiency may cause unique changes in hepatic metabolism that are predictive of liver injury. To test this possibility, we first determined that Abcb11 knock-out (KO) C57BL/6J mice recapitulate human deficiency. Before the onset of cholestasis, Abcb11 KO mice have altered hepatic lipid metabolism coupled with reduced expression of genes important in mitochondrial fatty acid oxidation. This was associated with increased serum free-fatty acids, reduced total white adipose, and marked impairment of long-chain fatty acid β-oxidation. Importantly, metabolomic analysis confirmed that Abcb11 KO mice have impaired mitochondrial fatty acid β-oxidation with the elevated fatty acid metabolites phenylpropionylglycine and phenylacetylglycine. These metabolic changes precede cholestasis but may be of relevance to cholestatic disease progression because altered fatty acid metabolism can enhance reactive oxygen species that might exacerbate cholestatic liver damage.

Fasting glucose GWAS candidate region analysis across ethnic groups in the Multiethnic Study of Atherosclerosis (MESA)

Genetic variants associated with fasting glucose in European ancestry populations are increasingly well understood. However, the nature of the associations between these single nucleotide polymorphisms (SNPs) and fasting glucose in other racial and ethnic groups is unclear. We sought to examine regions previously identified to be associated with fasting glucose in Caucasian genome-wide association studies (GWAS) across multiple ethnicities in the Multiethnic Study of Atherosclerosis (MESA). Nondiabetic MESA participants with fasting glucose measured at the baseline exam and with GWAS genotyping were included; 2,349 Caucasians, 664 individuals of Chinese descent, 1,366 African Americans, and 1,171 Hispanics. Genotype data were generated from the Affymetrix 6.0 array and imputation in IMPUTE. Fasting glucose was regressed on SNP dosage data in each ethnic group adjusting for age, gender, MESA study center, and ethnic-specific principal components. SNPs from the three gene regions with the strongest associations to fasting glucose in previous Caucasian GWAS (MTNR1B / GCK / G6PC2) were examined in depth. There was limited power to replicate associations in other ethnic groups due to smaller allele frequencies and limited sample size; SNP associations may also have differed across ethnic groups due to differing linkage disequilibrium patterns with causal variants. rs10830963 in MTNR1B and rs4607517 in GCK demonstrated consistent magnitude and direction of association with fasting glucose across ethnic groups, although the associations were often not nominally significant. In conclusion, certain SNPs in MTNR1B and GCK demonstrate consistent effects across four racial and ethnic groups, narrowing the putative region for these causal variants.

Type 2 Diabetes Genetics: Beyond GWAS

The global epidemic of type 2 diabetes mellitus (T2D) is one of the most challenging problems of the 21(st) century leading cause of and the fifth death worldwide. Substantial evidence suggests that T2D is a multifactorial disease with a strong genetic component. Recent genome-wide association studies (GWAS) have successfully identified and replicated nearly 75 susceptibility loci associated with T2D and related metabolic traits, mostly in Europeans, and some in African, and South Asian populations. The GWAS serve as a starting point for future genetic and functional studies since the mechanisms of action by which these associated loci influence disease is still unclear and it is difficult to predict potential implication of these findings in clinical settings. Despite extensive replication, no study has unequivocally demonstrated their clinical role in the disease management beyond progression to T2D from impaired glucose tolerance. However, these studies are revealing new molecular pathways underlying diabetes etiology, gene-environment interactions, epigenetic modifications, and gene function. This review highlights evolving progress made in the rapidly moving field of T2D genetics that is starting to unravel the pathophysiology of a complex phenotype and has potential to show clinical relevance in the near future.

Genetic determinants of variability in glycated hemoglobin (HbA(1c)) in humans: review of recent progress and prospects for use in diabetes care

Glycated hemoglobin A(1c) (HbA(1c)) indicates the percentage of total hemoglobin that is bound by glucose, produced from the nonenzymatic chemical modification by glucose of hemoglobin molecules carried in erythrocytes. HbA(1c) represents a surrogate marker of average blood glucose concentration over the previous 8 to 12 weeks, or the average lifespan of the erythrocyte, and thus represents a more stable indicator of glycemic status compared with fasting glucose. HbA(1c) levels are genetically determined, with heritability of 47% to 59%. Over the past few years, inroads into understanding genetic predisposition by glycemic and nonglycemic factors have been achieved through genomewide analyses. Here I review current research aimed at discovering genetic determinants of HbA(1c) levels, discussing insights into biologic factors influencing variability in the general and diabetic population, and across different ethnicities. Furthermore, I discuss briefly the relevance of findings for diabetes monitoring and diagnosis.

Deletion of the G6pc2 gene encoding the islet-specific glucose-6-phosphatase catalytic subunit-related protein does not affect the progression or incidence of type 1 diabetes in NOD/ShiLtJ mice

OBJECTIVE

Islet-specific glucose-6-phosphatase catalytic subunit-related protein (IGRP), now known as G6PC2, is a major target of autoreactive T cells implicated in the pathogenesis of type 1 diabetes in both mice and humans. This study aimed to determine whether suppression of G6p2 gene expression might therefore prevent or delay disease progression.

RESEARCH DESIGN AND METHODS

G6pc2(-/-) mice were generated on the NOD/ShiLtJ genetic background, and glycemia was monitored weekly up to 35 weeks of age to determine the onset and incidence of diabetes. The antigen specificity of CD8(+) T cells infiltrating islets from NOD/ShiLtJ G6pc2(+/+) and G6pc2(-/-) mice at 12 weeks was determined in parallel.

RESULTS

The absence of G6pc2 did not affect the time of onset, incidence, or sex bias of type 1 diabetes in NOD/ShiLtJ mice. Insulitis was prominent in both groups, but whereas NOD/ShiLtJ G6pc2(+/+) islets contained CD8(+) T cells reactive to the G6pc2 NRP peptide, G6pc2 NRP-reactive T cells were absent in NOD/ShiLtJ G6pc2(-/-) islets.

CONCLUSIONS

These results demonstrate that G6pc2 is an important driver for the selection and expansion of islet-reactive CD8(+) T cells infiltrating NOD/ShiLtJ islets. However, autoreactivity to G6pc2 is not essential for the emergence of autoimmune diabetes. The results remain consistent with previous studies indicating that insulin may be the primary autoimmune target, at least in NOD/ShiLtJ mice.

Association of genetic Loci with glucose levels in childhood and adolescence: a meta-analysis of over 6,000 children

OBJECTIVE

To investigate whether associations of common genetic variants recently identified for fasting glucose or insulin levels in nondiabetic adults are detectable in healthy children and adolescents.

RESEARCH DESIGN AND METHODS

A total of 16 single nucleotide polymorphisms (SNPs) associated with fasting glucose were genotyped in six studies of children and adolescents of European origin, including over 6,000 boys and girls aged 9-16 years. We performed meta-analyses to test associations of individual SNPs and a weighted risk score of the 16 loci with fasting glucose.

RESULTS

Nine loci were associated with glucose levels in healthy children and adolescents, with four of these associations reported in previous studies and five reported here for the first time (GLIS3, PROX1, SLC2A2, ADCY5, and CRY2). Effect sizes were similar to those in adults, suggesting age-independent effects of these fasting glucose loci. Children and adolescents carrying glucose-raising alleles of G6PC2, MTNR1B, GCK, and GLIS3 also showed reduced β-cell function, as indicated by homeostasis model assessment of β-cell function. Analysis using a weighted risk score showed an increase [β (95% CI)] in fasting glucose level of 0.026 mmol/L (0.021-0.031) for each unit increase in the score.

CONCLUSIONS

Novel fasting glucose loci identified in genome-wide association studies of adults are associated with altered fasting glucose levels in healthy children and adolescents with effect sizes comparable to adults. In nondiabetic adults, fasting glucose changes little over time, and our results suggest that age-independent effects of fasting glucose loci contribute to long-term interindividual differences in glucose levels from childhood onwards.

Associations of common genetic variants with age-related changes in fasting and postload glucose: evidence from 18 years of follow-up of the Whitehall II cohort

OBJECTIVE

In the general, nondiabetic population, fasting glucose increases only slightly over time, whereas 2-h postload glucose shows a much steeper age-related rise. The reasons underlying these different age trajectories are unknown. We investigated whether common genetic variants associated with fasting and 2-h glucose contribute to age-related changes of these traits.

RESEARCH DESIGN AND METHODS

We studied 5,196 nondiabetic participants of the Whitehall II cohort (aged 40-78 years) attending up to four 5-yearly oral glucose tolerance tests. A genetic score was calculated separately for fasting and 2-h glucose, including 16 and 5 single nucleotide polymorphisms, respectively. Longitudinal modeling with age centered at 55 years was used to study the effects of each genotype and genetic score on fasting and 2-h glucose and their interactions with age, adjusting for sex and time-varying BMI.

RESULTS

The fasting glucose genetic score was significantly associated with fasting glucose with a 0.029 mmol/L (95% CI 0.023-0.034) difference (P = 2.76 × 10(-21)) per genetic score point, an association that remained constant over time (age interaction P = 0.17). Two-hour glucose levels differed by 0.076 mmol/L (0.047-0.105) per genetic score point (P = 3.1 × 10(-7)); notably, this effect became stronger with increasing age by 0.006 mmol/L (0.003-0.009) per genetic score point per year (age interaction P = 3.0 × 10(-5)), resulting in diverging age trajectories by genetic score.

CONCLUSIONS

Common genetic variants contribute to the age-related rise of 2-h glucose levels, whereas associations of variants for fasting glucose are constant over time, in line with stable age trajectories of fasting glucose.