Identification and functional characterization of G6PC2 coding variants influencing glycemic traits define an effector transcript at the G6PC2-ABCB11 locus

Mendelian randomization study on the causal relationships among fasting blood glucose, plasma proteins, and squamous cell lung cancer

BACKGROUND

Squamous cell lung cancer (SQCLC) represents the second most common subtype of lung cancer (LC) with characteristics of treatment resistance. Plasma proteins often influence levels of fasting blood glucose (FBG), consequently impacting LC. However, the precise role of FBG in this association remains unclear.

OBJECTIVE

To investigate the causal relationships of FBG with LC and its subtypes, plasma proteins, and SQCLC, as well as the mediating role of FBG.

METHODS

Mendelian randomization (MR) analysis was employed to assess the causal associations of FBG with LC and its subtypes, plasma proteins and SQCLC, and plasma proteins and FBG, using the two-step MR approach with the primary method being Inverse Variance Weighted (IVW). Protein-Protein Interaction (PPI) network was utilized to identify hub genes of plasma proteins causally linked to SQCLC.

RESULTS

FBG was a risk factor for SQCLC (OR: 1.376, 95% CI 1.017-1.862, P = 0.038) but had no significant causal associations with LC and other subtypes (P > 0.05). Furthermore, 54 plasma proteins had significant causal associations with SQCLC (P < 0.05). EEF2 K (OR: 1.111, 95% CI 1.015-1.216, P = 0.023) and SSR1 (OR: 0.546, 95% CI 0.487-0.613, P < 0.001) were identified as a risk and protective factor for FBG, respectively. Mediation analysis indicated a significant negative mediating effect of FBG in the causal relationship between SSR1 and SQCLC (B = - 0.193, 95% CI - 0.312-0.074, P = 0.001), with a mediation proportion of 44.4%.

CONCLUSION

Our study revealed FBG as a risk factor for SQCLC and demonstrated the mediating role of FBG in the causal association between SSR1 and SQCLC.

Efficient storage and regression computation for population-scale genome sequencing studies

MOTIVATION

The growing availability of large-scale population biobanks has the potential to significantly advance our understanding of human health and disease. However, the massive computational and storage demands of whole genome sequencing (WGS) data pose serious challenges, particularly for underfunded institutions or researchers in developing countries. This disparity in resources can limit equitable access to cutting-edge genetic research.

RESULTS

We present novel algorithms and regression methods that dramatically reduce both computation time and storage requirements for WGS studies, with particular attention to rare variant representation. By integrating these approaches into PLINK 2.0, we demonstrate substantial gains in efficiency without compromising analytical accuracy. In an exome-wide association analysis of 19.4 million variants for the body mass index phenotype in 125 077 individuals (AllofUs project data), we reduced runtime from 695.35 min (11.5 h) on a single machine to 1.57 min with 30 GB of memory and 50 threads (or 8.67 min with 4 threads). Additionally, the framework supports multi-phenotype analyses, further enhancing its flexibility.

AVAILABILITY AND IMPLEMENTATION

Our optimized methods are fully integrated into PLINK 2.0 and can be accessed at: https://www.cog-genomics.org/plink/2.0/.

The role of glucose-6-phosphatase activity in glucose homeostasis and its potential for diabetes therapy

Glucose-6-phosphatase catalytic subunit (G6PC)1 and G6PC2 are crucial for glucose metabolism, regulating processes like glycolysis, gluconeogenesis, and glycogenolysis. Despite their structural and functional similarities, G6PC1 and G6PC2 exhibit distinct tissue-specific expression patterns, G6P hydrolysis kinetics, and physiological functions. This review provides a comprehensive overview of their enzymology and distinct roles in glucose homeostasis. We examine how inactivating mutations in G6PC1 lead to glycogen storage disease, and how elevated G6PC1 and G6PC2 expression can affect the incidence of diabetic complications, risk for type 2 diabetes mellitus (T2DM) and various cancers. We also discuss the potential of inhibiting G6PC1 and G6PC2 to protect against complications from elevated blood glucose levels, and highlight drug development efforts targeting G6PC1 and G6PC2, and the therapeutic potential of inhibitors for disease prevention.

ZHX3 interacts with CEBPB to repress hepatic gluconeogenic gene expression and uric acid secretion

ZHX3, which encodes for a transcriptional repressor, is associated with fasting blood glucose (FBG) levels and increased type 2 diabetes (T2D) risk but its role in cell types involved in glucose metabolism is not well understood. Here, we show that the deletion of ZHX3 in the human pancreatic β-cell line EndoC-βH1 did not impair glucose-stimulated insulin secretion (GSIS) nor perturb its transcriptome. On the other hand, we found that ZHX3 represses the expression of gluconeogenic genes PCK1 and G6PC1 in the human hepatoma line HepG2. Transcriptomic analysis of ZHX3-deficient HepG2 cells revealed that the uric acid transporter gene SLC17A1 was up-regulated, which consequentially led to increased uric acid secretion. High levels of uric acid could then impair GSIS in EndoC-βH1 cells. Subsequently, in-depth co-immunoprecipitation followed by mass spectrometry analysis of ZHX3 in HepG2 cells identified transcription factor CEBPB as its binding partner, required to repress the transcription of PCK1, G6PC1, and partially SLC17A1 in HepG2 cells. Overall, our study uncovered the role of ZHX3 in regulating glucose metabolism in hepatocytes, thereby influencing FBG levels and their association with T2D risk.

Efficient storage and regression computation for population-scale genome sequencing studies

In the era of big data in human genetics, large-scale biobanks aggregating genetic data from diverse populations have emerged as important for advancing our understanding of human health and disease. However, the computational and storage demands of whole genome sequencing (WGS) studies pose significant challenges, especially for researchers from underfunded institutions or developing countries, creating a disparity in research capabilities. We introduce new approaches that significantly enhance computational efficiency and reduce data storage requirements for WGS studies. By developing algorithms for compressed storage of genetic data, focusing particularly on optimizing the representation of rare variants, and designing regression methods tailored for the scale and complexity of WGS data, we significantly lower computational and storage costs. We integrate our approach into PLINK 2.0. The implementation demonstrates considerable reductions in storage space and computational time without compromising analytical accuracy, as evidenced by the application to the AllofUs project data. We optimized the runtime of an exome-wide association analysis involving 19.4 million variants and the body mass index phenotype of 125,077 individuals, reducing it from 695.35 minutes (approximately 11.5 hours) on a single machine to just 1.57 minutes using 30 GB of memory and 50 threads (or 8.67 minutes with 4 threads). Additionally, we extended this approach to support multi-phenotype analyses. We anticipate that our approach will enable researchers across the globe to unlock the potential of population biobanks, accelerating the pace of discoveries that can improve our understanding of human health and disease.

Understanding the Genetic Landscape of Gestational Diabetes: Insights into the Causes and Consequences of Elevated Glucose Levels in Pregnancy

Background/Objectives: During pregnancy, physiological changes in maternal circulating glucose levels and its metabolism are essential to meet maternal and fetal energy demands. Major changes in glucose metabolism occur throughout pregnancy and consist of higher insulin resistance and a compensatory increase in insulin secretion to maintain glucose homeostasis. For some women, this change is insufficient to maintain normoglycemia, leading to gestational diabetes mellitus (GDM), a condition characterized by maternal glucose intolerance and hyperglycaemia first diagnosed during the second or third trimester of pregnancy. GDM is diagnosed in approximately 14.0% of pregnancies globally, and it is often associated with short- and long-term adverse health outcomes in both mothers and offspring. Although recent studies have highlighted the role of genetic determinants in the development of GDM, research in this area is still lacking, hindering the development of prevention and treatment strategies. Methods: In this paper, we review recent advances in the understanding of genetic determinants of GDM and glycaemic traits during pregnancy. Results/Conclusions: Our review highlights the need for further collaborative efforts as well as larger and more diverse genotyped pregnancy cohorts to deepen our understanding of the genetic aetiology of GDM, address research gaps, and further improve diagnostic and treatment strategies.

Causal link between metabolic related factors and osteoarthritis: a Mendelian randomization investigation

Introduction

Metabolic syndrome (MetS) is significantly associated with osteoarthritis (OA), especially in MetS patients with blood glucose abnormalities, such as elevated fasting blood glucose (FG), which may increase OA risk. Dietary modifications, especially the intake of polyunsaturated fatty acids (PUFAs), are regarded as a potential means of preventing MetS and its complications. However, regarding the effects of FG, Omega-3s, and Omega-6s on OA, the research conclusions are conflicting, which is attributed to the complexity of the pathogenesis of OA. Therefore, it is imperative to thoroughly evaluate multiple factors to fully understand their role in OA, which needs further exploration and clarification.

Methods

Two-sample univariable Mendelian randomization (UVMR) and multivariable Mendelian randomization (MVMR) were employed to examine the causal effect of metabolic related factors on hip OA (HOA) or knee OA (KOA). The exposure and outcome datasets were obtained from Open GWAS IEU. All cases were independent European ancestry data. Three MR methods were performed to estimate the causal effect: inverse-variance weighting (IVW), weighted median method (WMM), and MR-Egger regression. Additionally, the intercept analysis in MR-Egger regression is used to estimate pleiotropy, and the IVW method and MR-Egger regression are used to test the heterogeneity.

Results

The UVMR analysis revealed a causal relationship between FG and HOA. By MVMR analysis, the study discovered a significant link between FG (OR = 0.79, 95%CI: 0.64∼0.99, p = 0.036) and KOA after accounting for body mass index (BMI), age, and sex hormone-binding globulin (SHBG). However, no causal effects of FG on HOA were seen. Omega-3s and Omega-6s did not have a causal influence on HOA or KOA. No significant evidence of pleiotropy was identified.

Discussion

The MR investigation showed a protective effect of FG on KOA development but no causal relationship between FG and HOA. No causal effect of Omega-3s and Omega-6s on HOA and KOA was observed. Shared genetic overlaps might also exist between BMI and age, SHBG and PUFAs for OA development. This finding offers a novel insight into the treatment and prevention of KOA from glucose metabolism perspective. The FG cutoff value should be explored in the future, and consideration should be given to demonstrating the study in populations other than Europeans.

Biochemical and metabolic characterization of a G6PC2 inhibitor

Three glucose-6-phosphatase catalytic subunits, that hydrolyze glucose-6-phosphate (G6P) to glucose and inorganic phosphate, have been identified, designated G6PC1-3, but only G6PC1 and G6PC2 have been implicated in the regulation of fasting blood glucose (FBG). Elevated FBG has been associated with multiple adverse clinical outcomes, including increased risk for type 2 diabetes and various cancers. Therefore, G6PC1 and G6PC2 inhibitors that lower FBG may be of prophylactic value for the prevention of multiple conditions. The studies described here characterize a G6PC2 inhibitor, designated VU0945627, previously identified as Compound 3. We show that VU0945627 preferentially inhibits human G6PC2 versus human G6PC1 but activates human G6PC3. VU0945627 is a mixed G6PC2 inhibitor, increasing the Km but reducing the Vmax for G6P hydrolysis. PyRx virtual docking to an AlphaFold2-derived G6PC2 structural model suggests VU0945627 binds two sites in human G6PC2. Mutation of residues in these sites reduces the inhibitory effect of VU0945627. VU0945627 does not inhibit mouse G6PC2 despite its 84% sequence identity with human G6PC2. Mutagenesis studies suggest this lack of inhibition of mouse G6PC2 is due, in part, to a change in residue 318 from histidine in human G6PC2 to proline in mouse G6PC2. Surprisingly, VU0945627 still inhibited glucose cycling in the mouse islet-derived βTC-3 cell line. Studies using intact mouse liver microsomes and PyRx docking suggest that this observation can be explained by an ability of VU0945627 to also inhibit the G6P transporter SLC37A4. These data will inform future computational modeling studies designed to identify G6PC isoform-specific inhibitors.

MultiSuSiE improves multi-ancestry fine-mapping in All of Us whole-genome sequencing data

Leveraging data from multiple ancestries can greatly improve fine-mapping power due to differences in linkage disequilibrium and allele frequencies. We propose MultiSuSiE, an extension of the sum of single effects model (SuSiE) to multiple ancestries that allows causal effect sizes to vary across ancestries based on a multivariate normal prior informed by empirical data. We evaluated MultiSuSiE via simulations and analyses of 14 quantitative traits leveraging whole-genome sequencing data in 47k African-ancestry and 94k European-ancestry individuals from All of Us. In simulations, MultiSuSiE applied to Afr47k+Eur47k was well-calibrated and attained higher power than SuSiE applied to Eur94k; interestingly, higher causal variant PIPs in Afr47k compared to Eur47k were entirely explained by differences in the extent of LD quantified by LD 4th moments. Compared to very recently proposed multi-ancestry fine-mapping methods, MultiSuSiE attained higher power and/or much lower computational costs, making the analysis of large-scale All of Us data feasible. In real trait analyses, MultiSuSiE applied to Afr47k+Eur94k identified 579 fine-mapped variants with PIP > 0.5, and MultiSuSiE applied to Afr47k+Eur47k identified 44% more fine-mapped variants with PIP > 0.5 than SuSiE applied to Eur94k. We validated MultiSuSiE results for real traits via functional enrichment of fine-mapped variants. We highlight several examples where MultiSuSiE implicates well-studied or biologically plausible fine-mapped variants that were not implicated by other methods.

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.

Rare Heterozygous Loss-of-Function Variants in the Human GLP-1 Receptor Are Not Associated With Cardiometabolic Phenotypes

CONTEXT

Lost glucagon-like peptide 1 receptor (GLP-1R) function affects human physiology.

OBJECTIVE

This work aimed to identify coding nonsynonymous GLP1R variants in Danish individuals to link their in vitro phenotypes and clinical phenotypic associations.

METHODS

We sequenced GLP1R in 8642 Danish individuals with type 2 diabetes or normal glucose tolerance and examined the ability of nonsynonymous variants to bind GLP-1 and to signal in transfected cells via cyclic adenosine monophosphate (cAMP) formation and β-arrestin recruitment. We performed a cross-sectional study between the burden of loss-of-signaling (LoS) variants and cardiometabolic phenotypes in 2930 patients with type 2 diabetes and 5712 participants in a population-based cohort. Furthermore, we studied the association between cardiometabolic phenotypes and the burden of the LoS variants and 60 partly overlapping predicted loss-of-function (pLoF) GLP1R variants found in 330 566 unrelated White exome-sequenced participants in the UK Biobank cohort.

RESULTS

We identified 36 nonsynonymous variants in GLP1R, of which 10 had a statistically significant loss in GLP-1-induced cAMP signaling compared to wild-type. However, no association was observed between the LoS variants and type 2 diabetes, although LoS variant carriers had a minor increased fasting plasma glucose level. Moreover, pLoF variants from the UK Biobank also did not reveal substantial cardiometabolic associations, despite a small effect on glycated hemoglobin A1c.

CONCLUSION

Since no homozygous LoS nor pLoF variants were identified and heterozygous carriers had similar cardiometabolic phenotype as noncarriers, we conclude that GLP-1R may be of particular importance in human physiology, due to a potential evolutionary intolerance of harmful homozygous GLP1R variants.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Omics biomarkers and an approach for their practical implementation to delineate health status for personalized nutrition strategies

Personalized nutrition (PN) has gained much attention as a tool for empowerment of consumers to promote changes in dietary behavior, optimizing health status and preventing diet related diseases. Generalized implementation of PN faces different obstacles, one of the most relevant being metabolic characterization of the individual. Although omics technologies allow for assessment the dynamics of metabolism with unprecedented detail, its translatability as affordable and simple PN protocols is still difficult due to the complexity of metabolic regulation and to different technical and economical constrains. In this work, we propose a conceptual framework that considers the dysregulation of a few overarching processes, namely Carbohydrate metabolism, lipid metabolism, inflammation, oxidative stress and microbiota-derived metabolites, as the basis of the onset of several non-communicable diseases. These processes can be assessed and characterized by specific sets of proteomic, metabolomic and genetic markers that minimize operational constrains and maximize the information obtained at the individual level. Current machine learning and data analysis methodologies allow the development of algorithms to integrate omics and genetic markers. Reduction of dimensionality of variables facilitates the implementation of omics and genetic information in digital tools. This framework is exemplified by presenting the EU-Funded project PREVENTOMICS as a use case.

Loss of RREB1 in pancreatic beta cells reduces cellular insulin content and affects endocrine cell gene expression

AIMS/HYPOTHESIS

Genome-wide studies have uncovered multiple independent signals at the RREB1 locus associated with altered type 2 diabetes risk and related glycaemic traits. However, little is known about the function of the zinc finger transcription factor Ras-responsive element binding protein 1 (RREB1) in glucose homeostasis or how changes in its expression and/or function influence diabetes risk.

METHODS

A zebrafish model lacking rreb1a and rreb1b was used to study the effect of RREB1 loss in vivo. Using transcriptomic and cellular phenotyping of a human beta cell model (EndoC-βH1) and human induced pluripotent stem cell (hiPSC)-derived beta-like cells, we investigated how loss of RREB1 expression and activity affects pancreatic endocrine cell development and function. Ex vivo measurements of human islet function were performed in donor islets from carriers of RREB1 type 2 diabetes risk alleles.

RESULTS

CRISPR/Cas9-mediated loss of rreb1a and rreb1b function in zebrafish supports an in vivo role for the transcription factor in beta cell mass, beta cell insulin expression and glucose levels. Loss of RREB1 also reduced insulin gene expression and cellular insulin content in EndoC-βH1 cells and impaired insulin secretion under prolonged stimulation. Transcriptomic analysis of RREB1 knockdown and knockout EndoC-βH1 cells supports RREB1 as a novel regulator of genes involved in insulin secretion. In vitro differentiation of RREB1KO/KO hiPSCs revealed dysregulation of pro-endocrine cell genes, including RFX family members, suggesting that RREB1 also regulates genes involved in endocrine cell development. Human donor islets from carriers of type 2 diabetes risk alleles in RREB1 have altered glucose-stimulated insulin secretion ex vivo, consistent with a role for RREB1 in regulating islet cell function.

CONCLUSIONS/INTERPRETATION

Together, our results indicate that RREB1 regulates beta cell function by transcriptionally regulating the expression of genes involved in beta cell development and function.

Bidirectional Mendelian randomization study of insulin-related traits and risk of ovarian cancer

BACKGROUND

It is well known that the occurrence and development of ovarian cancer are closely related to the patient's weight and various endocrine factors in the body.

AIM

Mendelian randomization (MR) was used to analyze the bidirectional relationship between insulin related characteristics and ovarian cancer.

METHODS

The data on insulin related characteristics are from up to 5567 diabetes free patients from 10 studies, mainly including fasting insulin level, insulin secretion rate, peak insulin response, etc. For ovarian cancer, UK Biobank data just updated in 2021 was selected, of which the relevant gene data was from 199741 Europeans. Mendelian randomization method was selected, with inverse variance weighting (IVW) as the main estimation, while MR Pleiotropy, MR Egger, weighted median and other methods were used to detect the heterogeneity of data and whether there was multi validity affecting conclusions.

RESULTS

Among all insulin related indicators (fasting insulin level, insulin secretion rate, peak insulin response), the insulin secretion rate was selected to have a causal relationship with the occurrence of ovarian cancer (IVW, P < 0.05), that is, the risk of ovarian cancer increased with the decrease of insulin secretion rate. At the same time, we tested the heterogeneity and polymorphism of this indicator, and the results were non-existent, which ensured the accuracy of the analysis results. Reverse causal analysis showed that there was no causal effect between the two (P>0.05).

CONCLUSION

The impairment of the insulin secretion rate has a causal effect on the risk of ovarian cancer, which was confirmed by Mendel randomization. This suggests that the human glucose metabolism cycle represented by insulin secretion plays an important role in the pathogenesis of ovarian cancer, which provides a new idea for preventing the release of ovarian cancer.

Identification of Drought-Tolerance Genes in the Germination Stage of Soybean

Drought stress influences the vigor of plant seeds and inhibits seed germination, making it one of the primary environmental factors adversely affecting food security. The seed germination stage is critical to ensuring the growth and productivity of soybeans in soils prone to drought conditions. We here examined the genetic diversity and drought-tolerance phenotypes of 410 accessions of a germplasm diversity panel for soybean and conducted quantitative genetics analyses to identify loci associated with drought tolerance of seed germination. We uncovered significant differences among the diverse genotypes for four growth indices and five drought-tolerance indices, which revealed abundant variation among genotypes, upon drought stress, and for genotype × treatment effects. We also used 158,327 SNP markers and performed GWAS for the drought-related traits. Our data met the conditions (PCA + K) for using a mixed linear model in TASSEL, and we thus identified 26 SNPs associated with drought tolerance indices for germination stage distributed across 10 chromosomes. Nine SNP sites, including, for example, Gm20_34956219 and Gm20_36902659, were associated with two or more phenotypic indices, and there were nine SNP markers located in or adjacent to (within 500 kb) previously reported drought tolerance QTLs. These SNPs led to our identification of 41 candidate genes related to drought tolerance in the germination stage. The results of our study contribute to a deeper understanding of the genetic mechanisms underlying drought tolerance in soybeans at the germination stage, thereby providing a molecular basis for identifying useful soybean germplasm for breeding new drought-tolerant varieties.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

The contribution of common and rare genetic variants to variation in metabolic traits in 288,137 East Asians

Metabolic traits are heritable phenotypes widely-used in assessing the risk of various diseases. We conduct a genome-wide association analysis (GWAS) of nine metabolic traits (including glycemic, lipid, liver enzyme levels) in 125,872 Korean subjects genotyped with the Korea Biobank Array. Following meta-analysis with GWAS from Biobank Japan identify 144 novel signals (MAF ≥ 1%), of which 57.0% are replicated in UK Biobank. Additionally, we discover 66 rare (MAF < 1%) variants, 94.4% of them co-incident to common loci, adding to allelic series. Although rare variants have limited contribution to overall trait variance, these lead, in carriers, substantial loss of predictive accuracy from polygenic predictions of disease risk from common variant alone. We capture groups with up to 16-fold variation in type 2 diabetes (T2D) prevalence by integration of genetic risk scores of fasting plasma glucose and T2D and the I349F rare protective variant. This study highlights the need to consider the joint contribution of both common and rare variants on inherited risk of metabolic traits and related diseases.

Insulin Response to Oral Glucose and Cardiometabolic Disease: A Mendelian Randomization Study to Assess Potential Causality

Mendelian randomization (MR) suggests that postprandial hyperinsulinemia (unadjusted for plasma glucose) increases BMI, but its impact on cardiometabolic disease, a leading cause for mortality and morbidity in people with obesity, is not established. Fat distribution i.e., increased centripetal and/or reduced femoro-gluteal adiposity, is causally associated with and better predicts cardiometabolic disease than BMI. We therefore undertook bidirectional MR to assess the effect of corrected insulin response (CIR) (insulin 30 min after a glucose challenge adjusted for plasma glucose) on BMI, waist-to-hip ratio (WHR), leg fat, type 2 diabetes (T2D), triglyceride (TG), HDL, liver fat, hypertension (HTN), and coronary artery disease (CAD) in people of European descent. Inverse variance-weighted MR suggests a potential causal association between increased CIR and increased BMI (b = 0.048 ± 0.02, P = 0.03), increased leg fat (b = 0.029 ± 0.012, P = 0.01), reduced T2D (b = -0.73 ± 0.15, P = 6 × 10-7, odds ratio [OR] 0.48 [95% CI 0.36-0.64]), reduced TG (b = -0.07 ± 0.02, P = 0.003), and increased HDL (b = 0.04 ± 0.01, P = 0.006) with some evidence of horizontal pleiotropy. CIR had neutral effects on WHR (b = 0.009 ± 0.02, P = 0.69), liver fat (b = -0.08 ± 0.04, P = 0.06), HTN (b = -0.001 ± 0.004, P = 0.7, OR 1.00 [95% CI 0.99-1.01]), and CAD (b = -0.002 ± 0.002, P = 0.48, OR 0.99 [95% CI 0.81-1.21]). T2D decreased CIR (b -0.22 ± 0.04, P = 1.3 × 10-7), with no evidence that BMI, TG, HDL, liver fat, HTN, and CAD modulate CIR. In conclusion, we did not find evidence that increased CIR increases cardiometabolic disease. It might increase BMI with favorable fat distribution, reduce T2D, and improve lipids.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

Mutations in G6PC2 gene with increased risk for development of type 2 diabetes: Understanding via computational approach

An increase in the fast blood glucose (FBG) levels has been linked to an increased risk of developing a chronic condition, type 2 diabetes (T2D). The mutation in the G6PC2 gene was identified to have a lead role in the modulation of FBG levels. The abnormal regulation of this enzyme influences glucose-stimulated insulin secretion (GSIS), which controls the insulin levels corresponding to the system's glucose level. This study focuses on the mutations at the G6PC2 gene, which cause the variation from normal expression levels and increase the risk of T2D. We examined the non-synonymous single nucleotide polymorphisms (nsSNPs) present in the G6PC2 and subjected them to pathogenicity, stability, residue conservation, and membrane simulation. The individual representation of surrounding amino acids in the mutant (I63T) model showed the loss of hydrophobic interactions compared to the native G6PC2. In addition, the trajectory results from the membrane simulation exhibited reduced stability, and the least compactness was identified for the I63T mutant model. Our study shed light on the structural and conformational changes at the transmembrane region due to the I63T mutation in G6PC2. Additionally, the Gibbs free energy landscape analysis against the two principal components showed structural differences and decreased the conformational stability of the I63T mutant model compared to the native. Like those presented in this study, dynamical simulations may indeed be crucial to comprehending the structural insights of G6PC2 mutations in cardiovascular-associated mortality and T2D.

Empowering consumers to PREVENT diet-related diseases through OMICS sciences (PREVENTOMICS): protocol for a parallel double-blinded randomised intervention trial to investigate biomarker-based nutrition plans for weight loss

INTRODUCTION

Personalised nutrition holds immense potential over conventional one-size-fits-all approaches for preventing and treating diet-related diseases, such as obesity. The current study aims to examine whether a personalised nutritional plan produces more favourable health outcomes than a standard approach based on general dietary recommendations in subjects with overweight or obesity and elevated waist circumference.

METHODS AND ANALYSIS

This project is a 10-week parallel, double-blinded randomised intervention trial. We plan to include 100 adults aged 18-65 years interested in losing weight, with body mass index ≥27 but<40 kg/m² and elevated waist circumference (males >94 cm; females >80 cm). Participants will be categorised into one of five predefined 'clusters' based on their individual metabolic biomarker profile and genetic background, and will be randomised in a 1:1 ratio to one of two groups: (1) personalised plan group that will receive cluster-specific meals every day for 6 days a week, in conjunction with a personalised behavioural change programme via electronic push notifications; or (2) control group that will receive meals following the general dietary recommendations in conjunction with generic health behaviour prompts. The primary outcome is the difference between groups (personalised vs control) in the change in fat mass from baseline. Secondary outcomes include changes in weight and body composition, fasting blood glucose and insulin, lipid profile, adipokines, inflammatory biomarkers, and blood pressure. Other outcomes involve measures of physical activity and sleep patterns, health-related quality of life, dietary intake, eating behaviour, and biomarkers of food intake. The effect of the intervention on the primary outcome will be analysed by means of linear mixed models.

ETHICS AND DISSEMINATION

The protocol has been approved by the Ethics Committee of the Capital Region, Copenhagen, Denmark. Study findings will be disseminated through peer-reviewed publications, conference presentations and media outlets.

TRIAL REGISTRATION NUMBER

NCT04590989.

A robust method for collider bias correction in conditional genome-wide association studies

Estimated genetic associations with prognosis, or conditional on a phenotype (e.g. disease incidence), may be affected by collider bias, whereby conditioning on the phenotype induces associations between causes of the phenotype and prognosis. We propose a method, 'Slope-Hunter', that uses model-based clustering to identify and utilise the class of variants only affecting the phenotype to estimate the adjustment factor, assuming this class explains more variation in the phenotype than any other variant classes. Simulation studies show that our approach eliminates the bias and outperforms alternatives even in the presence of genetic correlation. In a study of fasting blood insulin levels (FI) conditional on body mass index, we eliminate paradoxical associations of the underweight loci: COBLLI; PPARG with increased FI, and reveal an association for the locus rs1421085 (FTO). In an analysis of a case-only study for breast cancer mortality, a single region remains associated with more pronounced results.

Nonsynonymous single-nucleotide polymorphisms in the G6PC2 gene affect protein expression, enzyme activity, and fasting blood glucose

G6PC2 encodes a glucose-6-phosphatase (G6Pase) catalytic subunit that modulates the sensitivity of insulin secretion to glucose and thereby regulates fasting blood glucose (FBG). A common single-nucleotide polymorphism (SNP) in G6PC2, rs560887 is an important determinant of human FBG variability. This SNP has a subtle effect on G6PC2 RNA splicing, which raises the question as to whether nonsynonymous SNPs with a major impact on G6PC2 stability or enzyme activity might have a broader disease/metabolic impact. Previous attempts to characterize such SNPs were limited by the very low inherent G6Pase activity and expression of G6PC2 protein in islet-derived cell lines. In this study, we describe the use of a plasmid vector that confers high G6PC2 protein expression in islet cells, allowing for a functional analysis of 22 nonsynonymous G6PC2 SNPs, 19 of which alter amino acids that are conserved in mouse G6PC2 and the human and mouse variants of the related G6PC1 isoform. We show that 16 of these SNPs markedly impair G6PC2 protein expression (>50% decrease). These SNPs have variable effects on the stability of human and mouse G6PC1, despite the high sequence homology between these isoforms. Four of the remaining six SNPs impaired G6PC2 enzyme activity. Electronic health record-derived phenotype analyses showed an association between high-impact SNPs and FBG, but not other diseases/metabolites. While homozygous G6pc2 deletion in mice increases the risk of hypoglycemia, these human data reveal no evidence that the beneficial use of partial G6PC2 inhibitors to lower FBG would be associated with unintended negative consequences.

Dissecting the IL-6 pathway in cardiometabolic disease: a Mendelian randomization study on both IL6 and IL6R

BACKGROUND

Chronic inflammation is a risk factor for cardiovascular disease. IL-6 signaling perturbation through IL-6 or IL-6R blockade may have potential benefit on cardiovascular risk. It is unknown whether targeting either IL-6 or IL-6 receptor may result in similar effects on CVD and adverse events. We compared the anticipated effects of targeting IL-6 and IL-6 receptor on cardiometabolic risk and potential side effects.

METHODS

We constructed four instruments: two main instruments with genetic variants in the IL6 and IL6R loci weighted for their association with CRP, and two after firstly filtering variants for their association with IL-6 or IL-6R expression. Analyses were performed for coronary artery disease (CAD), ischemic stroke, atrial fibrillation (AF), heart failure, type 2 diabetes (T2D), rheumatoid arthritis (RA), infection endpoints, and quantitative hematological, metabolic, and anthropometric parameters.

RESULTS

A 1 mg/L lower CRP by the IL6 instrument was associated with lower CAD (OR 0.86, 95% CI 0.77;0.96), AF, and T2D risk. A 1mg/L lower CRP by the IL6R instrument was associated with lower CAD (OR 0.90, 95% CI 0.86;0.95), any stroke and ischemic stroke, AF, RA risk and higher pneumonia risk. The eQTL filtered results were in concordance with the main results, but with wider confidence intervals.

CONCLUSIONS

IL-6 signalling perturbation by either IL6 or IL6R genetic instruments is associated with a similar risk reduction for multiple cardiometabolic diseases, suggesting that both IL-6 and IL-6R are potential therapeutic targets to lower CVD. Moreover, IL-6 rather than IL-6R inhibition might have a more favorable pneumonia risk.

Haplotypes of the genes (GCK and G6PC2) underlying the glucose/glucose-6-phosphate cycle are associated with pancreatic beta cell glucose sensitivity in patients with newly diagnosed type 2 diabetes from the VNDS study (VNDS 11)

BACKGROUND

Elevated fasting plasma glucose has been associated with increased risk for development of type 2 diabetes (T2D). The balance between glucokinase (GCK) and glucose-6-phosphate catalytic subunit 2 (G6PC2) activity are involved in glucose homeostasis through glycolytic flux, and subsequent insulin secretion.

AIM

In this study, we evaluated the association between the genetic variability of G6PC2 and GCK genes and T2D-related quantitative traits.

METHODS

In 794 drug-naïve, GADA-negative, newly diagnosed T2D patients (VNDS; NTC01526720) we performed: genotyping of 6 independent tag-SNPs within GCK gene and 5 tag-SNPs within G6PC2 gene; euglycaemic insulin clamp to assess insulin sensitivity; OGTT to estimate beta-cell function (derivative and proportional control; DC, PC) by mathematical modeling. Genetic association analysis has been conducted using Plink software.

RESULTS

Two SNPs within GCK gene (rs882019 and rs1303722) were associated to DC in opposite way (both p < 0.004). Two G6PC2 variants (rs13387347 and rs560887) were associated to both parameters of insulin secretion (DC and PC) and to fasting C-peptide levels (all p < 0.038). Moreover, subjects carrying the A allele of rs560887 showed higher values of 2h-plasma glucose (2hPG) (p = 0.033). Haplotype analysis revealed that GCK (AACAAA) haplotype was associated to decreased fasting C-peptide levels, whereas, the most frequent haplotype of G6PC2 (GGAAG) was associated with higher fasting C-peptide levels (p = 0.001), higher PC (β = 6.87, p = 0.022) and the lower 2hPG (p = 0.012).

CONCLUSION

Our findings confirmed the role of GCK and G6PC2 in regulating the pulsatility in insulin secretion thereby influencing insulin-signaling and leading to a gradual modulation in glucose levels in Italian patients with newly diagnosed T2D.

Genome-wide association study identifying novel variant for fasting insulin and allelic heterogeneity in known glycemic loci in Chilean adolescents: The Santiago Longitudinal Study

BACKGROUND

The genetic underpinnings of glycemic traits have been understudied in adolescent and Hispanic/Latino (H/L) populations in comparison to adults and populations of European ancestry.

OBJECTIVE

To identify genetic factors underlying glycemic traits in an adolescent H/L population.

METHODS

We conducted a genome-wide association study (GWAS) of fasting glucose (FG) and fasting insulin (FI) in H/L adolescents from the Santiago Longitudinal Study.

RESULTS

We identified one novel variant positioned in the CSMD1 gene on chromosome 8 (rs77465890, effect allele frequency = 0.10) that was associated with FI (β = -0.299, SE = 0.054, p = 2.72×10-8 ) and was only slightly attenuated after adjusting for body mass index z-scores (β = -0.252, SE = 0.047, p = 1.03×10-7 ). We demonstrated directionally consistent, but not statistically significant results in African and Hispanic adults of the Population Architecture Using Genomics and Epidemiology Consortium. We also identified secondary signals for two FG loci after conditioning on known variants, which demonstrate allelic heterogeneity in well-known glucose loci.

CONCLUSION

Our results exemplify the importance of including populations with diverse ancestral origin and adolescent participants in GWAS of glycemic traits to uncover novel risk loci and expand our understanding of disease aetiology.

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.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Progress in Defining the Genetic Contribution to Type 2 Diabetes in Individuals of East Asian Ancestry

PURPOSE OF REVIEW

Prevalence of type 2 diabetes (T2D) and progression of complications differ between worldwide populations. While obesity is a major contributing risk factor, variations in physiological manifestations, e.g., developing T2D at lower body mass index in some populations, suggest other contributing factors. Early T2D genetic associations were mostly discovered in European ancestry populations. This review describes the progression of genetic discoveries associated with T2D in individuals of East Asian ancestry in the last 10 years and highlights the shared genetic susceptibility between the population groups and additional insights into genetic contributions to T2D.

RECENT FINDINGS

Through increased sample size and power, new genetic associations with T2D were discovered in East Asian ancestry populations, often with higher allele frequencies than European ancestry populations. As we continue to generate maps of T2D-associated variants across diverse populations, there will be a critical need to expand and diversify other omics resources to enable integration for clinical translation.

Human Pluripotent Stem Cells to Model Islet Defects in Diabetes

Diabetes mellitus is characterized by elevated levels of blood glucose and is ultimately caused by insufficient insulin production from pancreatic beta cells. Different research models have been utilized to unravel the molecular mechanisms leading to the onset of diabetes. The generation of pancreatic endocrine cells from human pluripotent stem cells constitutes an approach to study genetic defects leading to impaired beta cell development and function. Here, we review the recent progress in generating and characterizing functional stem cell-derived beta cells. We summarize the diabetes disease modeling possibilities that stem cells offer and the challenges that lie ahead to further improve these models.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

Nutrient consumption-dependent association of a glucagon-like peptide-1 receptor gene polymorphism with insulin secretion

Since type 2 diabetes (DM) is a life-style related disease, life-style should be considered when association between genetic factors and DM are examined. However, most studies did not examine genetic associations in consideration with lifestyle. Glucagon-like peptide-1 (GLP-1) receptor (GLP1R) mediates the insulinotropic action of GLP-1 in β-cells. We here examined the association while taking into consideration of interactions between the gene polymorphism and various nutrient factors. Participants from the population-based Iwaki study of Japanese subjects held in 2014–2017 with information on nutritional intake evaluated by self-administered dietary history questionnaire, and GLP1R genotype (rs3765467: A/G), were included (n = 1,560). Although not significant, insulin secretion indices assessed by homeostasis model assessment of β-cell function (HOMA-β) in subjects with the GG genotype tended to be lower than in those with the AA+AG genotypes in most groups stratified into tertiles based on daily nutrient consumptions (high, middle, and low). Stratification also showed that the GG genotype was a significant risk for decreased insulin secretion (HOMA-β ≤ 30) even after adjustment for multiple factors (age, body mass index, alcohol consumption), but only in the highest tertiles of energy, protein and carbohydrate consumption in men [odds ratios (95% confidence interval) 3.95 (1.03–15.1), 15.83 (1.58–158.9), and 4.23 (1.10–11.2), respectively]. A polymorphism of the GLP1R gene was associated with decreased insulin secretion in a nutrient consumption-dependent manner in Japanese men, indicating an interaction between GLP1R and nutritional factors in the pathophysiology of DM.

Convex combination sequence kernel association test for rare-variant studies

We propose a novel variant set test for rare-variant association studies, which leverages multiple single-nucleotide variant (SNV) annotations. Our approach optimizes a convex combination of different sequence kernel association test (SKAT) statistics, where each statistic is constructed from a different annotation and combination weights are optimized through a multiple kernel learning algorithm. The combination test statistic is evaluated empirically through data splitting. In simulations, we find our method preserves type I error at α = 2.5 × 1 0 - 6 and has greater power than SKAT(-O) when SNV weights are not misspecified and sample sizes are large ( N ≥ 5 , 000 ). We utilize our method in the Framingham Heart Study (FHS) to identify SNV sets associated with fasting glucose. While we are unable to detect any genome-wide significant associations between fasting glucose and 4-kb windows of rare variants ( p < 1 0 - 7 ) in 6,419 FHS participants, our method identifies suggestive associations between fasting glucose and rare variants near ROCK2 ( p = 2.1 × 1 0 - 5 ) and within CPLX1 ( p = 5.3 × 1 0 - 5 ). These two genes were previously reported to be involved in obesity-mediated insulin resistance and glucose-induced insulin secretion by pancreatic beta-cells, respectively. These findings will need to be replicated in other cohorts and validated by functional genomic studies.

Genetic variants associated with glycemic response to treatment with dipeptidylpeptidase 4 inhibitors

Aim: We examined associations of eight SNPs in/near seven candidate genes with glycemic response to 6 month treatment with DPP4 inhibitors. Patients & methods: 206 patients with type 2 diabetes (116 men and 90 women) were treated with sitagliptin or vildagliptin (both 100 mg/day) in combination with metformin or metformin/sulphonylurea over 6 months, and the reduction in glycated hemoglobin (HbA1c) was measured. Results: Rs6923761 in GLP1R was significantly associated with a reduction in HbA1c (adjusted p = 0.006). Homozygotes for the minor A allele had smaller reduction in HbA1c by 0.4% (4 mmol/mol) than the G allele carriers (p = 0.016). Conclusion: The missense variant rs6923761 in the GLP1R gene was associated with a smaller glycemic response to 6 month gliptin therapy in diabetic patients of central European origin.

From Genetic Association to Molecular Mechanisms for Islet-cell Dysfunction in Type 2 Diabetes

Genome-wide association studies (GWAS) have identified over 400 signals robustly associated with risk for type 2 diabetes (T2D). At the vast majority of these loci, the lead single nucleotide polymorphisms (SNPs) reside in noncoding regions of the genome, which hampers biological inference and translation of genetic discoveries into disease mechanisms. The study of these T2D risk variants in normoglycemic individuals has revealed that a significant proportion are exerting their disease risk through islet-cell dysfunction. The central role of the islet is also demonstrated by numerous studies, which have shown an enrichment of these signals in islet-specific epigenomic annotations. In recent years the emergence of authentic human beta-cell lines, and advances in genome-editing technologies coupled with improved protocols differentiating human pluripotent stem cells into beta-like cells has opened up new opportunities for T2D disease modeling. Here we review the current understanding on the genetic basis of T2D focusing on approaches, which have facilitated the identification of causal variants and their effector transcripts in human islets. We will present examples of functional studies based on animal and conventional cellular systems and highlight the potential of novel stem cell-based T2D disease models.

Transcription Factor RREB1: from Target Genes towards Biological Functions

The Ras-responsive element binding protein 1(RREB1) is a member of zinc finger transcription factors, which is widely involved in biological processes including cell proliferation, transcriptional regulation and DNA damage repair. New findings reveal RREB1 functions as both transcriptional repressors and transcriptional activators for transcriptional regulation of target genes. The activation of RREB1 is regulated by MAPK pathway. We have summarized the target genes of RREB1 and discussed RREB1 roles in the cancer development. In addition, increasing evidences suggest that RREB1 is a potential risk gene for type 2 diabetes and obesity. We also review the current clinical application of RREB1 as a biomarker for melanoma detection. In conclusion, RREB1 is a promising diagnostic biomarker or new drug target for cancers and metabolic diseases.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

The Contribution of Low-Frequency and Rare Coding Variation to Susceptibility to Type 2 Diabetes

PURPOSE OF REVIEW

Soon after the first genome-wide association study (GWAS) for type 2 diabetes (T2D) was published, it was hypothesized that rare and low-frequency variants might explain a substantial proportion of disease risk. Rare coding variants in particular were emphasized given their large expected role in disease. This review summarizes the extent to which recent T2D genetic studies provide evidence for or against this hypothesis.

RECENT FINDINGS

Following a comprehensive study of T2D genetic architecture using three sequencing and genotyping technologies, four even larger studies have provided a yet higher resolution view of the role of rare and low-frequency coding variation in T2D susceptibility. Empirical evidence strongly suggests that common regulatory variants are the dominant contributor to T2D heritability. However, rare coding variants may nonetheless be pervasive across T2D-relevant genes. A strategy using common variants to map disease genes, and rare coding variants to link molecular gene perturbations to cellular and phenotypic effects, may be an effective means to investigate T2D pathogenesis and potential new therapies.

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.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Pharmacogenetics of type 2 diabetes mellitus, the route toward tailored medicine

Type 2 diabetes mellitus (T2DM) is a chronic disease that has reached the levels of a global epidemic. In order to achieve optimal glucose control, it is often necessary to rely on combination therapy of multiple drugs or insulin because uncontrolled glucose levels result in T2DM progression and enhanced risk of complications and mortality. Several antihyperglycemic agents have been developed over time, and T2DM pharmacotherapy should be prescribed based on suitability for the individual patient's characteristics. Pharmacogenetics is the branch of genetics that investigates how our genome influences individual responses to drugs, therapeutic outcomes, and incidence of adverse effects. In this review, we evaluated the pharmacogenetic evidences currently available in the literature, and we identified the top informative genetic variants associated with response to the most common anti-diabetic drugs: metformin, DPP-4 inhibitors/GLP1R agonists, thiazolidinediones, and sulfonylureas/meglitinides. Overall, we found 40 polymorphisms for each drug class in a total of 71 loci, and we examined the possibility of encouraging genetic screening of these variants/loci in order to critically implement decision-making about the therapeutic approach through precision medicine strategies. It is possible then to anticipate that when the clinical practice will take advantage of the genetic information of the diabetic patients, this will provide a useful resource for the prevention of T2DM progression, enabling the identification of the precise drug that is most likely to be effective and safe for each patient and the reduction of the economic impact on a global scale.