Gly40Ser polymorphism of the glucagon receptor gene is associated with central adiposity in men

Determinants of plasma levels of proglucagon and the metabolic impact of glucagon receptor signalling: a UK Biobank study

AIMS/HYPOTHESES

Glucagon and glucagon-like peptide-1 (GLP-1) are derived from the same precursor; proglucagon, and dual agonists of their receptors are currently being explored for the treatment of obesity and metabolic dysfunction-associated steatotic liver disease (MASLD). Elevated levels of endogenous glucagon (hyperglucagonaemia) have been linked with hyperglycaemia in individuals with type 2 diabetes but are also observed in individuals with obesity and MASLD. GLP-1 levels have been reported to be largely unaffected or even reduced in similar conditions. We investigated potential determinants of plasma proglucagon and associations of glucagon receptor signalling with metabolic diseases based on data from the UK Biobank.

METHODS

We used exome sequencing data from the UK Biobank for ~410,000 white participants to identify glucagon receptor variants and grouped them based on their known or predicted signalling. Data on plasma levels of proglucagon estimated using Olink technology were available for a subset of the cohort (~40,000). We determined associations of glucagon receptor variants and proglucagon with BMI, type 2 diabetes and liver fat (quantified by liver MRI) and performed survival analyses to investigate if elevated proglucagon predicts type 2 diabetes development.

RESULTS

Obesity, MASLD and type 2 diabetes were associated with elevated plasma levels of proglucagon independently of each other. Baseline proglucagon levels were associated with the risk of type 2 diabetes development over a 14 year follow-up period (HR 1.13; 95% CI 1.09, 1.17; n=1562; p=1.3×10-12). This association was of the same magnitude across strata of BMI. Carriers of glucagon receptor variants with reduced cAMP signalling had elevated levels of proglucagon (β 0.847; 95% CI 0.04, 1.66; n=17; p=0.04), and carriers of variants with a predicted frameshift mutation had higher levels of liver fat compared with the wild-type reference group (β 0.504; 95% CI 0.03, 0.98; n=11; p=0.04).

CONCLUSIONS/INTERPRETATION

Our findings support the suggestion that glucagon receptor signalling is involved in MASLD, that plasma levels of proglucagon are linked to the risk of type 2 diabetes development, and that proglucagon levels are influenced by genetic variation in the glucagon receptor, obesity, type 2 diabetes and MASLD. Determining the molecular signalling pathways downstream of glucagon receptor activation may guide the development of biased GLP-1/glucagon co-agonist with improved metabolic benefits.

DATA AVAILABILITY

All coding is available through https://github.com/nicwin98/UK-Biobank-GCG.

Downregulation of the kidney glucagon receptor, essential for renal function and systemic homeostasis, contributes to chronic kidney disease

The glucagon receptor (GCGR) in the kidney is expressed in nephron tubules. In humans and animal models with chronic kidney disease, renal GCGR expression is reduced. However, the role of kidney GCGR in normal renal function and in disease development has not been addressed. Here, we examined its role by analyzing mice with constitutive or conditional kidney-specific loss of the Gcgr. Adult renal Gcgr knockout mice exhibit metabolic dysregulation and a functional impairment of the kidneys. These mice exhibit hyperaminoacidemia associated with reduced kidney glucose output, oxidative stress, enhanced inflammasome activity, and excess lipid accumulation in the kidney. Upon a lipid challenge, they display maladaptive responses with acute hypertriglyceridemia and chronic proinflammatory and profibrotic activation. In aged mice, kidney Gcgr ablation elicits widespread renal deposition of collagen and fibronectin, indicative of fibrosis. Taken together, our findings demonstrate an essential role of the renal GCGR in normal kidney metabolic and homeostatic functions. Importantly, mice deficient for kidney Gcgr recapitulate some of the key pathophysiological features of chronic kidney disease.

Naturally occurring mutations in G protein-coupled receptors associated with obesity and type 2 diabetes mellitus

G protein-coupled receptors (GPCRs) are the largest family of membrane receptors involved in the regulation of almost all known physiological processes. Dysfunctions of GPCR-mediated signaling have been shown to cause various diseases. The prevalence of obesity and type 2 diabetes mellitus (T2DM), two strongly associated disorders, is increasing worldwide, with tremendous economical and health burden. New safer and more efficacious drugs are required for successful weight reduction and T2DM treatment. Multiple GPCRs are involved in the regulation of energy and glucose homeostasis. Mutations in these GPCRs contribute to the development and progression of obesity and T2DM. Therefore, these receptors can be therapeutic targets for obesity and T2DM. Indeed some of these receptors, such as melanocortin-4 receptor and glucagon-like peptide 1 receptor, have provided important new drugs for treating obesity and T2DM. This review will focus on the naturally occurring mutations of several GPCRs associated with obesity and T2DM, especially incorporating recent large genomic data and insights from structure-function studies, providing leads for future investigations.

Molecular and in vivo phenotyping of missense variants of the human glucagon receptor

Naturally occurring missense variants of G protein–coupled receptors with loss of function have been linked to metabolic disease in case studies and in animal experiments. The glucagon receptor, one such G protein–coupled receptor, is involved in maintaining blood glucose and amino acid homeostasis; however, loss-of-function mutations of this receptor have not been systematically characterized. Here, we observed fewer glucagon receptor missense variants than expected, as well as lower allele diversity and fewer variants with trait associations as compared with other class B1 receptors. We performed molecular pharmacological phenotyping of 38 missense variants located in the receptor extracellular domain, at the glucagon interface, or with previously suggested clinical implications. These variants were characterized in terms of cAMP accumulation to assess glucagon-induced Gα_s coupling, and of recruitment of β-arrestin-1/2. Fifteen variants were impaired in at least one of these downstream functions, with six variants affected in both cAMP accumulation and β-arrestin-1/2 recruitment. For the eight variants with decreased Gα_s signaling (D63^ECDN, P86^ECDS, V96^ECDE, G125^ECDC, R225^3.30H, R308^5.40W, V368^6.59M, and R378^7.35C) binding experiments revealed preserved glucagon affinity, although with significantly reduced binding capacity. Finally, using the UK Biobank, we found that variants with wildtype-like Gα_s signaling did not associate with metabolic phenotypes, whereas carriers of cAMP accumulation-impairing variants displayed a tendency toward increased risk of obesity and increased body mass and blood pressure. These observations are in line with the essential role of the glucagon system in metabolism and support that Gα_s is the main signaling pathway effecting the physiological roles of the glucagon receptor.

Renal potassium handling in carriers of the Gly40Ser mutation of the glucagon receptor suggests a role for glucagon in potassium homeostasis

Plasma potassium concentration (PK ) is tightly regulated. Insulin is known to favor potassium entry into cells. But how potassium leaves the cells later on is not often considered. Previous studies in rats showed that glucagon infusion increased urinary potassium excretion dose-dependently and reversibly. This prompted us to investigate the possible influence of glucagon on potassium handling in humans. We took advantage of the Gly40Ser mutation of the glucagon receptor (GR) that results in a partial loss of function of the GR. In the Olivetti cohort (male workers), 25 subjects who carried this mutation were matched 1:4 to 100 noncarriers for age and weight. Estimated osmolarity of serum and 24-h urine (Sosm and Uosm, respectively) was calculated from the concentrations of the main solutes: [(Na+K)*2 + urea (+glucose for serum)]. Transtubular potassium gradient (TTKG), reflecting the intensity of K secretion in the distal nephron, was calculated as [(urine K/serum K)(Uosm /Sosm )]. There was no significant difference in serum K, or 24-h urine urea, Na and K excretion rates. But urine K concentration was significantly lower in carriers than in noncarriers. Means (interquartile range): 38 (34-43) versus 47 (43-51) mmol/L, P = 0.030. TTKG was also significantly lower in carriers: 4.2 (3.9-4.6) versus 5.0 (4.7-5.2), P = 0.015. This difference remained statistically significant after adjustments for serum insulin and 24-h Na and urea excretions. These results in humans suggest that glucagon stimulates K secretion in the distal nephron. Thus, in conjunction with insulin, glucagon may also participate in K homeostasis by promoting renal K excretion.

Differential impact of amino acid substitutions on critical residues of the human glucagon-like peptide-1 receptor involved in peptide activity and small-molecule allostery

The glucagon-like peptide-1 receptor (GLP-1R) is a class B G protein-coupled receptor that has a critical role in the regulation of glucose homeostasis, principally through the regulation of insulin secretion. The receptor system is highly complex, able to be activated by both endogenous [GLP-1(1-36)NH2, GLP-1(1-37), GLP-1(7-36)NH2, GLP-1(7-37), oxyntomodulin], and exogenous (exendin-4) peptides in addition to small-molecule allosteric agonists (compound 2 [6,7-dichloro-2-methylsulfonyl-3-tert-butylaminoquinoxaline], BETP [4-(3-benzyloxy)phenyl)-2-ethylsulfinyl-6-(trifluoromethyl)pyrimidine]). Furthermore, the GLP-1R is subject to single-nucleotide polymorphic variance, resulting in amino acid changes in the receptor protein. In this study, we investigated two polymorphic variants previously reported to impact peptide-mediated receptor activity (M149) and small-molecule allostery (C333). These residues were mutated to a series of alternate amino acids, and their functionality was monitored across physiologically significant signaling pathways, including cAMP, extracellular signal-regulated kinase 1 and 2 phosphorylation, and intracellular Ca(2+) mobilization, in addition to peptide binding and cell-surface expression. We observed that residue 149 is highly sensitive to mutation, with almost all peptide responses significantly attenuated at mutated receptors. However, most reductions in activity were able to be restored by the small-molecule allosteric agonist compound 2. Conversely, mutation of residue 333 had little impact on peptide-mediated receptor activation, but this activity could not be modulated by compound 2 to the same extent as that observed at the wild-type receptor. These results provide insight into the importance of residues 149 and 333 in peptide function and highlight the complexities of allosteric modulation within this receptor system.

Minireview: Signal bias, allosterism, and polymorphic variation at the GLP-1R: implications for drug discovery

The glucagon-like peptide-1 receptor (GLP-1R) controls the physiological responses to the incretin hormone glucagon-like peptide-1 and is a major therapeutic target for the treatment of type 2 diabetes, owing to the broad range of effects that are mediated upon its activation. These include the promotion of glucose-dependent insulin secretion, increased insulin biosynthesis, preservation of β-cell mass, improved peripheral insulin action, and promotion of weight loss. Regulation of GLP-1R function is complex, with multiple endogenous and exogenous peptides that interact with the receptor that result in the activation of numerous downstream signaling cascades. The current understanding of GLP-1R signaling and regulation is limited, with the desired spectrum of signaling required for the ideal therapeutic outcome still to be determined. In addition, there are several single-nucleotide polymorphisms (used in this review as defining a natural change of single nucleotide in the receptor sequence; clinically, this is viewed as a single-nucleotide polymorphism only if the frequency of the mutation occurs in 1% or more of the population) distributed within the coding sequence of the receptor protein that have the potential to produce differential responses for distinct ligands. In this review, we discuss the current understanding of GLP-1R function, in particular highlighting recent advances in the field on ligand-directed signal bias, allosteric modulation, and probe dependence and the implications of these behaviors for drug discovery and development.

The Human Obesity Gene Map: The 2003 Update

This is the tenth update of the human obesity gene map, incorporating published results up to the end of October 2003 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. Transgenic and knockout murine models relevant to obesity are also incorporated (N = 55). As of October 2003, 41 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. QTLs reported from animal models currently number 183. There are 208 human QTLs for obesity phenotypes from genome-wide scans and candidate regions in targeted studies. A total of 35 genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 272 studies reporting positive associations with 90 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, more than 430 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.

The Human Obesity Gene Map: The 2004 Update

This paper presents the eleventh update of the human obesity gene map, which incorporates published results up to the end of October 2004. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTLs) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2004, 173 human obesity cases due to single-gene mutations in 10 different genes have been reported, and 49 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 166 genes which, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 221. The number of human obesity QTLs derived from genome scans continues to grow, and we have now 204 QTLs for obesity-related phenotypes from 50 genome-wide scans. A total of 38 genomic regions harbor QTLs replicated among two to four studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably with 358 findings of positive associations with 113 candidate genes. Among them, 18 genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, >600 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful publications and genomic and other relevant sites can be found at http://obesitygene.pbrc.edu.

Genetic susceptibility to obstructive sleep apnea in the obese child

The etiology of obstructive sleep apnea (OSA) is multifactorial, consisting of a complex interplay between anatomic and neuromuscular factors and an underlying genetic predisposition toward this disease. Several of the factors that have been reported to play a role in the pathogenesis of OSA can serve as intermediate phenotypes in investigations targeting genetic susceptibility to OSA. A precise underpinning of the genetic basis of OSA has been thus far difficult because it is still unknown whether or not the recognized candidate genes for OSA are directly causal to the phenotype, or whether their effects on OSA are mediated through the intermediate phenotypes of OSA. Future studies utilizing phenotypically homogenous groups such as those with childhood OSA and technological advances such as haplotype analysis in a case control design are extremely promising. Developing predictive models that incorporate genetic and phenotypic markers will enable early diagnosis and, therefore, intervention, ultimately resulting in reduction of morbidity and of the public health concerns associated with OSA in obese children.

The human obesity gene map: the 2005 update

This paper presents the 12th update of the human obesity gene map, which incorporates published results up to the end of October 2005. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTL) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2005, 176 human obesity cases due to single-gene mutations in 11 different genes have been reported, 50 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 244 genes that, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 408. The number of human obesity QTLs derived from genome scans continues to grow, and we now have 253 QTLs for obesity-related phenotypes from 61 genome-wide scans. A total of 52 genomic regions harbor QTLs supported by two or more studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably, with 426 findings of positive associations with 127 candidate genes. A promising observation is that 22 genes are each supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. The electronic version of the map with links to useful publications and relevant sites can be found at http://obesitygene.pbrc.edu.

Comparison of Gene Expression in Intra-Abdominal and Subcutaneous Fat: A Study of Men with Morbid Obesity and Nonobese Men Using Microarray and Proteomics

Extent of intra-abdominal fat had significant linear relations with six metabolic coronary risk factors: systolic and diastolic blood pressure, fasting blood concentrations of glucose, high density lipoprotein (HDL) cholesterol, triglyceride, and plasminogen activator inhibitor-1. Tumor necrosis factor-alpha and adiponectin can be biological mediators from the intra-abdominal fat to the metabolic coronary risk factors. Complementarily, we describe a new study that will analyze the gene expression in intra-abdominal and subcutaneous fat on mRNA and protein level using high throughput methods. The study will elucidate further whether intra-abdominal obesity is the common denominator for the different components of the metabolic syndrome.

Genome-wide linkage analysis for severe obesity in french caucasians finds significant susceptibility locus on chromosome 19q

To ascertain whether distinct chromosomal loci existed that were linked to severe obesity, as well as to utilize the increased heritability of this excessive phenotype, we performed a genome-wide scan in severely obese French Caucasians. The 109 selected pedigrees, totaling 447 individuals, required both the proband and a sibling to be severely obese (BMI >or=35 kg/m(2)), and 84.8% of the nuclear families possessed >or=1 morbidly obese sibling (BMI >or=40). Severe and morbid obesity are still relatively rare in France, with rates of 2.5 and 0.6%, respectively. The initial genome scan consisted of 395 evenly spaced microsatellite markers. Six regions were found to have suggestive linkage on 4q, 6cen-q, 17q, and 19q for a BMI >or=35 phenotypic subset, and 5q and 10q for an inclusive BMI >or=27 group. The highest peak on chromosome 19q (logarithm of odds [LOD] = 3.59) was significant by genome scan simulation testing (P = 0.042). These regions then underwent second-stage mapping with an additional set of 42 markers. BMI >or=35 analysis defined regions on 17q23.3-25.1 and 19q13.33-13.43 with an maximum likelihood score LOD of 3.16 and 3.21, respectively. Subsequent pooled data analysis with an additional previous population of 66 BMI >or=35 sib-pairs led to a significant LOD score of 3.8 at the 19q locus (empirical P = 0.023). For more moderate obesity and overweight susceptibility loci, BMI >or=27 analysis confirmed suggestive linkage to chromosome regions 5q14.3-q21.3 (LOD = 2.68) and 10q24.32-26.2 (LOD = 2.47). Plausible positional candidate genes include NR1H2 and TULP2.

Association studies of genetic polymorphisms in central obesity: a critical review

During the past decade, mutations affecting liability to central obesity have been discovered at a phenomenal rate, and despite few consistently replicated findings, a number of intriguing results have emerged in the literature. Association studies have been proposed to identify the genetic determinants of complex traits such as central obesity. The advantages of the association method include its relative robustness to genetic heterogeneity and the ability to detect much smaller effect sizes than is detectable using feasible sample sizes in linkage studies. However, the current literature linking central obesity to genetic variants is teeming with reports of associations that either cannot be replicated or for which corroboration by linkage has been impossible to find. Explanations for this lack of reproducibility are well rehearsed, and typically include poor study design, incorrect assumptions about the underlying genetic architecture, and simple overinterpretation of data. These limitations create concern about the validity of association studies and cause problems in establishing robust criteria for undertaking association studies. In this article, the current status of the literature of association studies for genetic dissection of central obesity is critically reviewed.

The human obesity gene map: the 2002 update

This is the ninth update of the human obesity gene map, incorporating published results through October 2002 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and various animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. For the first time, transgenic and knockout murine models exhibiting obesity as a phenotype are incorporated (N = 38). As of October 2002, 33 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and the causal genes or strong candidates have been identified for 23 of these syndromes. QTLs reported from animal models currently number 168; there are 68 human QTLs for obesity phenotypes from genome-wide scans. Additionally, significant linkage peaks with candidate genes have been identified in targeted studies. Seven genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 222 studies reporting positive associations with 71 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. More than 300 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.