The CREBRF diabetes-protective rs373863828-A allele is associated with enhanced early insulin release in men of Māori and Pacific ancestry

Hannah J. Burden, Shannon Adams, Braydon Kulatea, Morag Wright-McNaughton, Danielle Sword, Jennifer J. Ormsbee, Conor Watene-O’Sullivan, Tony R. Merriman, Jennifer L. Knopp, J. Geoffrey Chase, Jeremy D. Krebs, Rosemary M. Hall, Lindsay D. Plank, Rinki Murphy, Peter R. Shepherd, Troy L. Merry
doi: https://doi.org/10.1101/2021.02.27.21252567

Abstract

Aim The minor A allele of rs373863828 (CREBRF p.Arg457Gln) is associated with increased body mass index (BMI), but reduced risk of type 2 and gestational diabetes in Polynesian (Pacific peoples and Aotearoa New Zealand Māori) populations. This study investigates the effect of the A allele on insulin release and sensitivity in overweight/obese men without diabetes.

Methods A mixed meal tolerance test was completed by 172 men (56 with the A allele) of Māori or Pacific ancestry, and 44 (24 with the A allele) had a frequently sampled intravenous glucose tolerance test and hyperinsulinemic-euglycemic clamp. Mixed linear models with covariates age, ancestry and BMI were used to analyse the association between the A allele of rs373863828 and markers of insulin release and blood glucose regulation.

Results The A allele of rs373863828 is associated with a greater increase in plasma insulin 30 min following a meal challenge without affecting the elevation in plasma glucose or incretins GLP-1 or GIP. Consistent with this point, following an intravenous infusion of a glucose bolus, participants with an A allele had higher early (p<0.05 at 2 and 4 min) plasma insulin and C-peptide concentrations for a similar elevation in blood glucose as those homozygous for the major (G) allele. Despite increased plasma insulin, rs373863828 genotype was not associated with a significant difference (p>0.05) in insulin sensitivity index or glucose disposal during hyperinsulinemic-euglycemic clamp.

Conclusion/interpretation rs373863828-A allele associates with increased glucose-stimulated insulin release without affecting insulin sensitivity, suggesting that CREBRF p.Arg457Gln may increase maximal ability for β-cells to release insulin to reduce the risk of type 2 diabetes.

Introduction

Type 2 diabetes (T2D) is a heterogeneous polygenic disease resulting from a complex interaction with environmental and genetic factors [1]. Population-specific genetic variation impacts T2D risk [2-4]. The A allele of rs373863828, a Arg457Gln missense variant in the CREBRF gene, has been associated with a 30-40% reduction in risk of T2D in people of Polynesian (Aotearoa NZ Māori and Pacific peoples) [4-6] and Oceanic (Marianas and Micronesian) [7] ancestries, as well as reduced risk of gestational diabetes in Māori and Pacific women with obesity [8]. CREBRF rs373863828 has a minor allele frequency of 5–27% in people of Polynesian and Oceanic ancestry but is rare (<0.01%) in non-Pacific populations [4, 5, 9, 10]. The mechanism by which this genetic variant protects from diabetes is not known.

At a population level body mass index (BMI), a widely used marker for adiposity and obesity, positively correlates with risk of T2D [11] and higher BMI associates with poorer metabolic outcomes in Māori and Pacific peoples [12-14]. Paradoxically, the rs373863828-A allele not only associates with reduced risk of T2D, but also an increase in BMI of 1.36-1.45 kg/m2 (p ≤ 1 × 10−5) per allele [4, 5], which is the largest effect currently known of a common single variant on BMI [15]. T2D is characterised by defects in insulin secretion, insulin action, or both, resulting in chronic hyperglycemia [16]. Excess adipose tissue can reduce insulin sensitivity and impair the ability of insulin to act on peripheral tissues to dispose of glucose and supress hepatic glucose production following a meal [17].

The Gln (A) allele of CREBRF p.Arg457Gln promotes fat cell development and adiposity [4, 9], suggesting that the protective effect of the A allele on T2D may relate to preferential lipid storage in subcutaneous fat depots over metabolically more harmful visceral or ectopic fat deposition [15]. While genetically-driven metabolically-favourable fat distribution phenotypes have been reported [18, 19], there is little evidence of CREBRF Arg457Gln being a factorable adiposity allele. Indeed, the A allele associates with greater (rather than reduced) abdominal/waist circumference [4, 5], a surrogate measure for centrally located and thus metabolically unfavourable fat mass.

An alternative hypothesis for the inverse risk relationship between T2D and BMI is that higher relative muscle mass may partially account for the increased BMI, increasing the ability to dispose of glucose. Consistent with this, the A allele of rs373863828 associates with greater lean mass in infants [20]. Both higher muscle mass and metabolically-favourable adipose tissue distribution could improve insulin sensitivity to protect from diabetes. However, whether the CREBRF Arg457Gln affects insulin sensitivity is not clear, since in people without diabetes the A allele associates with similar [7] or increased fasted insulin resistance [4] as assessed by HOMA-IR. Interestingly, however, Hanson et al. [7] reported 20% higher HOMA-β per A allele in Oceanic peoples without diabetes. HOMA-β is a marker of pancreatic β-cell function, and a decline in functional β-cell mass is a defining characteristic of T2D pathogenesis [16]. Therefore, it is possible that the CREBRF Arg457Gln variant promotes insulin sensitivity or β-cell function to maintain glucose homeostasis to slow or prevent the development of T2D.

Here we investigated the effect of the A allele of rs373863828 on insulin sensitivity and insulin secretion in men of Māori and Pacific ancestry. We report the A allele is not associated with changes in insulin sensitivity assessed by hyperinsulinemic-euglycemic clamp. However, in response to a mixed meal challenge and intravenous glucose bolus, those with an A allele had higher plasma insulin concentration. These results suggest the CREBRF Arg457Gln variant may enhance β-cell function to protect from diabetes.

Materials and Methods

Participants

Healthy men of Te Moana Nui a Kiwa [Pacific Ocean; Aotearoa New Zealand Māori and/or Pacific peoples (Polynesian)] were recruited to take part in the study in either Auckland or Wellington, Aotearoa New Zealand. Informed written consent was provided, and this study was approved by the Health and Disability Ethics Committee of New Zealand (17/STH/79). Inclusion criteria were being male, aged 18-45 years, having two or more grandparents who self-report as being of Māori and/or Pacific descent [21], and being overweight or obese (BMI between 25-45kg/m2). Exclusion criteria included being diagnosed with a chronic illness, including cardiovascular, diabetes, and or related metabolic disease.

Mixed Meal Tolerance Test (MMTT)

Participants were required to report to the laboratory between 0800 and 0900 hrs following a 10 hour overnight fast. They were asked to refrain from exercise, smoking and caffeine for the 24 hours prior to their visit. A cannula was placed in an antecubital vein for blood sampling, and a baseline blood sample collected prior to participants consuming a standardised mixed meal consisting of a 98 g English muffin, 13 g of butter, 23 g raspberry jam and a chocolate nutritional drink (Complan, London UK; 55 g powder in 300 ml of water) with a macronutrient break down of 71 g carbohydrate, 23 g fat and 19 g protein. Subsequent blood samples were collected at t= 30, 60, 90, 120 and 150 minutes post meal.

Intravenous glucose tolerance test and hyperinsulinemic-euglycaemic clamp

Forty-four participants were recalled based on genotype (20 GG for rs373863828, and 24 AG or AA) at least two weeks following the mixed meal tolerance test to undertake a frequently sampled intravenous glucose tolerance test (FS-IVGTT) and hyperinsulinemic-euglycemic clamp [22]. All visits commenced at 0800 hrs and participants were asked to fast for 10 hrs overnight and refrained from exercise and caffeine for the previous 24 hours. Body composition was determined by dual energy x-ray absorptiometry (DXA, model iDXA, GE Healthcare, Madison, WI) and a cannula was then placed into the antecubital vein for infusion of 20% dextrose and insulin while another was placed in a wrist vein of the opposite arm for blood sampling. The blood sampling hand and wrist were placed in a heat box set at 55°C for 1 min prior to samples being taken to achieve venous blood arterialisation. Following baseline sampling a 20% dextrose bolus (0.3 g/kg, to a maximum of 30 g) was administered through the antecubital cannula within two minutes and blood samples were taken at t= 0, 2, 4, 6, 8, 10, 30 and 60 minutes. Following the 60 min blood sample, an insulin and 20% dextrose infusion were started. The priming insulin infusion dose consisted of 3 mins at 90, 75 and 60 mU/m2 of body surface area/min, after which insulin infusion rate was set to 45 mU/m2 of body surface area/min for the remainder of the clamp period. Body surface area was calculated using the Mosteller formula [23]. Blood samples were taken every 5 mins for the period of the clamp to test blood glucose (HemoCue Glucose 201+ system, Angelholm, Sweden) and the 20% dextrose infusion rate was adjusted to maintain blood glucose at 5 mmol/L. Blood samples were taken and plasma recovered at t=160, 170 and 180 mins (steady state phase) for analysis.

Blood collection and analysis

Blood samples were collected into either EDTA or p800 (for incretins) vacutainers. EDTA tubes were centrifuged at 2000 rpm, 4°C for 10 min and p800 at 1300 rpm at 4°C for 20 mins, and the plasma was recovered and stored at -20°C until analysis. Plasma glucose, insulin, C-peptide, non-esterified fatty acids (NEFA) and triglycerides (TAG), high (HDL) and low density lipoprotein (LDL), alanine aminotransferase (ALT), C-reactive protein (CRP) and cholesterol (CHOL) were measured using Roche COBAS e 411 and c 311 analysers and commercially available COBAS elecsys assays (Roche, Mannheim, Germany). Glucagon-like peptide-1 (GLP-1) and gastric inhibitory polypeptide (GIP) were analysed using a Millipore Human Metabolic Hormone Magnetic Bead Panel (Merck, Darmstadt, Germany). A whole blood sample was taken at baseline for analysis of HbA1c (Afinion™ HbA1c, Abbott, US) and extraction of DNA. DNA was isolated using a commercially available kit (DNA extraction minikit, Qiagen, Hilden Germany). rs373863828 genotype was determined using a custom-designed Taqman probe-set (Applied Biosystems, Foster City, CA, USA) as described previously [5, 8].

Insulin sensitivity and secretion

Homeostatic model assessment of insulin resistance 2 (HOMA2-IR) was calculated via the Diabetes Trials Unit (Headington, Oxford, UK) online calculator described by Levy et al. [24]. Hyperinsulinemic-euglycemic clamp glucose disposal rate (GDR; glucose infusion rate (mL/hr)/ 0.3/weight (kg)) and insulin sensitivity index (ISI = GDR/((steady state insulin-fasting insulin) x steady state blood glucose) were calculated as described by Tripathy et al. [22]. Steady state blood glucose was defined as the average blood glucose at 160-180 min of the hyperinsulinemic-euglycemic clamp. Van Cauter et al. [25] kinetic parameters were used to calculate an insulin secretion curve for each subject using their measured C-peptide from the FS-IVGTT. The total area under the insulin secretion curve (AUC) from t=0 minutes to t=30 minutes was calculated as a measure of total units of secretion, the 1st phase insulin secretion response from t=0 to t=4 minutes, and the 2nd phase from t=5 to t=30 minutes. Disposition index was calculated as ISI from hyperinsulinemic-euglycemic clamp x average of plasma insulin at t= 2 min and 4 min of FS-IVGTT [26].

Statistical analysis

All statistical analysis was undertaken using SPSS version 26 (IBM Corp, Chicago, IL, USA). Multivariable linear regression models were used to test for associations between rs373863828 A allele and all parameters presented in Tables 1 and 2 with age, proportion of self-reported grandparents of Polynesian ancestry (ANC), and BMI as covariates as indicated. Statistical significance was accepted as p<0.05. Repeated measures mixed linear models were used to test associations between the A allele and all variables measured over time course experiments (MMTT and FS-IVGTT), with age, ANC and BMI as covariates. Where a statistically significant (p<0.05) main effect was seen, Sidak post-hoc analysis was undertaken, and statistical outliers were identified and removed using the Grubbs’ test.

View this table:
Table 1.

Characteristics and insulin secretion and insulin sensitivity markers in men of Polynesian ancestry with (AX) or without (GG) a A allele for rs373863828 who underwent the mixed meal tolerance test. Association significance was assessed via multivariable linear regression with age, ancestry and body mass index as covariates.

View this table:
Table 2.

Characteristics and markers of inulin release and insulin sensitivity in men of Polynesian ancestry with (AX) or without (GG) a A allele for rs373863828 and who underwent the intravenous glucose tolerance test and hyperinsulinemic-euglycemic clamp. Association significance was assessed via multivariable linear regression with age, ancestry and body mass index as covariates.

Results

One hundred and seventy-two participants took part in the mixed meal tolerance test (MMTT), of which 116 had a GG (CREBRF Arg457) genotype and 56 had one or two A alleles (GA or AA; CREBRF Gln457) for rs373863828. Due to limited numbers of participants who were homozygous for the A allele (n=6), these participants were grouped together with those who were heterozygous (AG) and this grouping is indicated by the genotype AX. By design, GG and AX genotype groups were similar in age and BMI (Table 1). In the fasted state, those participants with the GG or AX genotypes had similar plasma glucose, insulin, HbA1c, triglyceride and cholesterol concentrations. Similarly, the presence of the A allele did not associate with significantly different circulating concentrations of the inflammatory marker CRP, the liver health marker ALT or HOMA2-IR or HOMA2-β (Table 1).

rs373863828-A allele is associated with increased postprandial plasma insulin

There was a significant increase in blood glucose at 30 minutes following consumption of a standardised mixed meal, which returned to fasted levels by 60 minutes and was not affected by rs373863828 genotype (Fig. 1a). Despite the similar blood glucose response to a meal, the AX genotype group had higher plasma insulin at 30 minutes following the meal (p=0.006) (Fig. 1b). Meal responses in plasma C-peptide, non-esterified fatty acids and triglycerides (Fig. 1c-e) did not differ by rs373863828 genotype. Incretin hormones GLP-1 and GIP can potentiate β-cell insulin concentrations in response to a rise in blood glucose [27]. Given the elevated circulating insulin concentrations associated with the A allele during the mixed meal tolerance test, we investigated whether it might be due to increased concentrations of GLP-1 or GIP. While plasma GLP-1 and GIP increased (p<0.001) in response to the mixed meal challenge, it was not affected by rs373863828 genotype (Fig. 2 a & b).

Figure 1.
Figure 1. The A allele of rs373863828 is associated with higher postprandial plasma insulin.

Plasma glucose (a), insulin (b), c-peptide (c), non-esterified fatty acid (NEFA) (d) and triglyceride (e) concentration before (0 min) and following a standardised meal in men of Polynesian ancestry with (AX, n=56) or without (GG, n=113) a A allele for rs373863828. Data is presented as mean ± standard error of mean (SEM). Geno, Time and Inter indicate p-values for genotype, time and interaction (genotype x time) effect were assessed by mixed liner model with age, ancestry and body mass index as covariates. ** p<0.01 vs. GG at the indicated timepoint as determined by Sidak post-hoc analysis.

Figure 2.
Figure 2. rs373863828 genotype does not affect postprandial plasma incretins.

Plasma glucagon-like peptide-1 (GLP-1) (a) and gastric inhibitor peptide (GIP) (b) concentration before (0 min) and following a standardised meal in men of Polynesian ancestry with (AX, n=53) or without (GG, n=80) an A allele for rs373863828. Data is presented as mean ± standard error of mean (SEM). Geno, Time and Inter indicate p-values for genotype, time and interaction (genotype x time) effect as assessed by repeated measures mixed liner model with age, ancestry and body mass index as covariates.

rs373863828-A allele associates with increased insulin secretion in response to IV-glucose bolus

The subgroup who participated in the intravenous glucose tolerance test and hyperinsulinemic-euglycemic clamp included 20 participants with the GG and 24 with the AX genotype, who were similar in age, BMI, relative fat and lean mass (Table 2). The FS-IVGTT allows assessment of glucose-stimulated insulin secretion independent of oral absorption effects. Both genotypes had a similar peak in blood glucose concentrations of approximately 15 mmol/L during the first 2 minutes of the FS-IVGTT, with the AX genotype group having lower blood glucose at 30 min post IV-glucose bolus compared to the GG genotype group (Fig. 3a). The A allele was also associated with higher glucose-induced increase in plasma insulin and C-peptide (interaction effect p<0.05) during the FS-IVGTT, which was most evident in the 6 min following glucose bolus (Fig. 3b-c). As such, we modelled 1st (0-4 min) and 2nd (5-30 min) phase insulin secretion AUC using the Van Cauter et al. [25] method. Both genotype groups had similar total and 2nd phase insulin secretion (p>0.50). The A allele was not significantly associated with 1st phase insulin secretion (p=0.19) or disposition index (which assess the relationship between acute insulin response and insulin sensitivity [26, 28]; p=0.13), despite both variables tending to be approximately 1.4-fold higher in participants of the AX genotype (Table 2).

Figure 3.
Figure 3. The A allele of rs373863828 is associated with higher insulin secretion following a intravenous bolus of glucose.

Blood glucose (a), and plasma insulin (b) and c-peptide before (0 min) and following an intravenous glucose bolus (0.3mg/kg) in men of Polynesian ancestry with (AX, n=24) or without (GG, n=20) a A allele for rs373863828. Data is presented as mean ± standard error of mean (SEM). Geno, Time and Inter indicate p-values for genotype, time and interaction (genotype x time) effect as assessed by repeated measures mixed liner model with age, ancestry and body mass index as covariates. *p<0.05, **p <0.01 vs. GG at the indicated timepoint as determined by Sidak post-hoc analysis.

rs373863828-A allele does not associate with insulin sensitivity

Since elevated insulin concentrations for the same blood glucose concentrations in the mixed meal tolerance test can indicate insulin resistance, the hyperinsulinemic-euglycemic clamp was used to assess insulin sensitivity [22]. During the steady-state portion of the hyperinsulinemic-euglycemic clamp the two genotype groups had similar blood glucose, plasma insulin and glucose infusion rate (Table 2). Furthermore, rs373863828 genotype did not affect whole body glucose disposal rate (GDR) or insulin sensitivity index during the clamp (Table 2).

Discussion

The CREBRF Arg457Gln allele (rs373863828-A) is associated with reduced risk of T2D despite increased BMI [4, 5]. To investigate how the minor (A) allele of rs373863828 may protect from T2D, we assessed whole body glucose homeostasis of men of Māori and Pacific ancestry using a mixed meal challenge, intravenous glucose tolerance test and hyperinsulinemic-euglycemic clamp. We found insulin sensitivity to be similar between genotype groups, however, the A allele was associated with indices of increased early insulin release in response to an IV-glucose bolus and meal challenge. This result is consistent with a report that the A allele associates with increased HOMA-β in non-diabetic Oceanic peoples [7] and suggests that CREBRF rs373863828-A may improve β-cell function without altering insulin sensitivity.

Glucose-stimulated insulin secretion from pancreatic β-cells is biphasic, with a rapid 1st phase insulin secretion from a readily releasable pool of insulin followed by a reduced, but sustained 2nd phase insulin release which continues until normoglycemia is returned [29]. Blunting of 1st phase insulin response promotes postprandial hyperglycaemia and is a determinant of developing T2D [29]. While the participants in this study were normoglycemic, as assessed by HbA1c and fasting blood glucose, they were also overweight/obese (average BMI>32 mg/m2) [30] and had low level insulin resistance with an average HOMA2-IR index of ∼1.8 [31]. As such, it is possible this participant group had blunted early glucose-stimulated insulin secretion, and this blunting was attenuated in those with the A allele for rs373863828. This outcome would be consistent with the A allele protecting from T2D since preserving early phase insulin secretion is thought to be a promising T2D prevention strategy [29, 32], and restoration of early insulin release in response to weight loss is a predictor of T2D remission [33].

Plasma insulin concentrations were increased in the AX genotype group during the mixed meal tolerance test independent of changes in plasma GLP-1 and GIP, suggesting that the A allele may affect insulin release independent of gut-derived insulinotropic hormones. This outcome was supported by intravenously-infused glucose (bypassing absorption induced-gut hormone release) resulting in a greater elevation in plasma C-peptide and insulin in those with an A allele at rs373863828. Consistent with secretion, rather than clearance, elevating plasma insulin the AX group also had a more pronounced elevation in plasma C-peptide. C-peptide is produced when proinsulin is cleaved to insulin and is not as readily cleared from the circulation as insulin, making it a more consistent marker of secretion [34]. It is not entirely clear why insulin but not C-peptide was elevated at 30 min following a meal in the AX group. This may be the result of not capturing insulin secretion during the first 10 min following the meal, which where the genotype effect on insulin and C-peptide was most apparent during IV-GTT. Alternatively it may potentially reflect lower hepatic insulin clearance contributing to the elevated circulation insulin in the AX group following a meal [35]. However, insulin levels were similar between genotypes during the hyperinsulinemic-euglycemic clamp suggesting that any effect of the rs373863828-A allele on liver insulin clearance rate would most likely be subtle and/or potentially relate to a mechanism associated dependent on oral, as opposed to intravenous, glucose load [36].

Modelling of the insulin secretion profile did not detect significant differences in genotype, despite 1st phase insulin secretion appearing to be 1.4-fold greater in participants with the A allele. This may be the results of a ceiling effect since we would not expect overt β-cell in this normoglycemia cohort, and therefore it is important to next assess whether rs373863828 affects 1st phase insulin secretion under conditions where there is normally a substantial suppression in early insulin release such as in pre-diabetes or early stages of T2D [29]. This may also explain why a significant genotype difference in disposition index (insulin sensitivity multiplied by the amount of insulin secreted) was not evident. The hyperbolic relationship between insulin sensitivity and insulin secretion defines the disposition index, where changes in insulin sensitivity are compensated for by inverse changes in insulin secretion and a low disposition index associates with increased risk of T2D [37-39]. However, our cohort did not fit this hyperbolic model largely owning to a limited range of insulin sensitivity and insulin secretion profiles making it difficult determine whether rs373863828 genotype influences disposition index. It is also possible that the increased insulin levels associated with the A allele is a compensatory response to reduced peripheral tissue insulin sensitivity [16]. However, it seems unlikely given the hyperinsulinemic-euglycemic clamp showed the rs373863828 genotype does not affect insulin-induced glucose disposal and HOMA2-IR was similar between genotypes. Furthermore, the increased insulin response in the AX genotype group associated with a greater decline in blood glucose during IV-GTT, providing further evidence that the A allele of rs373863828 is not associating with insulin resistance.

The glucose bolus given in the IV-glucose tolerance test resulted in a >3-fold increase in blood glucose above fasting, which would likely result in near maximal glucose-stimulated insulin secretion. This may suggest that the A allele is associated with improved maximal ability for β-cells to release insulin, rather than responding to a greater glucose stimulus for release. Increased β-cell insulin secretion could be the result of improved secretory capacity or elevated β-cell mass [40]. While the effect of rs373863828 genotype on either β-cell development, mass or secretory capacity has not been investigated, β-cell replication during infancy is a major indicator of β-cell mass in adults [41].

Recently, the A allele of rs373863828 has been shown to associate with increased growth of lean and bone mass in the first 4-months of life in Samoan infants [20]. Whether β-cell mass follows this same trajectory has yet to be determined, but it is interesting to speculate that the A allele promotes greater β-cell mass by adulthood, and thus a greater capacity to produce insulin to control blood glucose concentration when faced with environmental factors that would otherwise promote β-cell loss and T2D. Alternatively, the A allele may promote β-cell resilience to metabolic stress-induced impairment in the internal mechanisms governing insulin secretion, dedifferentiation or death, all of which have been proposed to play a causal role in the transition from prediabetes to overt-T2D [42, 43].

CREBRF (also known as luman regulatory factor) encodes CREB3 regulatory factor protein that represses the activity of CREB3, a transmembrane transcription factor involved in coordinating the unfolded protein response (UPR) during endoplasmic reticulum (ER) stress [44-46]. While it is not known if the CREBRF p.Gln457 allele affects ER stress in β-cells, overexpression of CREBRF p.Gln457 does attenuate starvation-induced death of 3T3-L1 adipocytes [4] and elevated ER stress is one of the mechanisms through which starvation can initiate cell apoptosis [47]. Because pancreatic β-cells are inherently susceptible to ER stress-induced apoptosis and dysfunction [48], lowering ER-stress response to metabolic overload could prolong the duration in which β-cells can optimally release insulin and delay the onset of T2D.

There are a number of limitations when interpreting findings in the current study. Firstly, only men were included. Since the CREBRF p.Arg457Gln appears to increase height to a greater extent in New Zealand men of Polynesian ancestry than women [6], it is important to next assess whether this variant influences insulin secretion in females. Secondly, to avoid age and long-term disease as confounding factors, we recruited participants <45 years old and free of long-term disease. While this strategy is important to understanding whether CREBRF p.Arg457Gln alters glucose homeostasis prior to disease development, it is possible that the impact of the A allele changes with age or with development of age-related metabolic diseases and related medications. Indeed, when body composition was assessed by DXA scan, the A allele of rs373863828 associated with increased fat mass in old age [9], but, in contrast, also increased lean mass in infants [20]. As alluded to earlier it may be under conditions where there is increased pressure on the β-cells to produce insulin (such as in prediabetes) the effect of CREBRF p.Arg457Gln on β-cells function becomes more apparent.

Collectively, our findings provide the first evidence that the Oceanic-specific A allele of rs373863828 is associated with increased insulin release in response to a glucose challenge in overweight/obese young males of Māori and Pacific ancestry, and this appears to be independent of insulin sensitivity. We speculate that enhanced β-cell function or mass could underpin the T2D protective effect of the A allele of rs373863828.

Data Availability

Considered upon request

Competing interest

The authors declare no competing interests or disclosures.

Acknowledgements

He mihi nui tēnei ki ngā kai tuku taonga (we would like to thank the study participants for donating their time and tissue samples). We would like to thank Prof. Vicky Cameron and Sara Raudsepp for analysing plasma incretins, and Ngati Porou Hauora for their assistance in initiating this study. This study was funded by the Auckland Medical Research Foundation (AMRF), the Health Research Council (HRC) New Zealand and the Maurice Wilkins Center (MWC). TLM is supported by a Rutherford Discovery Fellowship.

Reference

  1. [1].↵
    Ali O (2013) Genetics of type 2 diabetes. World J Diabetes 4(4): 114–123. doi:10.4239/wjd.v4.i4.114
    OpenUrlCrossRefPubMed
  2. [2].↵
    Moltke I, Grarup N, Jorgensen ME, et al. (2014) A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes. Nature 512(7513): 190–193. doi:10.1038/nature13425
    OpenUrlCrossRefPubMed
  3. [3].
    Corona E, Chen R, Sikora M, et al. (2013) Analysis of the genetic basis of disease in the context of worldwide human relationships and migration. PLoS Genet 9(5): e1003447. doi:10.1371/journal.pgen.1003447
    OpenUrlCrossRefPubMed
  4. [4].↵
    Minster RL, Hawley NL, Su CT, et al. (2016) A thrifty variant in CREBRF strongly influences body mass index in Samoans. Nat Genet 48(9): 1049–1054. doi:10.1038/ng.3620
    OpenUrlCrossRefPubMed
  5. [5].↵
    Krishnan M, Major TJ, Topless RK, et al. (2018) Discordant association of the CREBRF rs373863828 A allele with increased BMI and protection from type 2 diabetes in Maori and Pacific (Polynesian) people living in Aotearoa/New Zealand. Diabetologia 61(7): 1603–1613. doi:10.1007/s00125-018-4623-1
    OpenUrlCrossRefPubMed
  6. [6].↵
    Metcalfe LK, Krishnan M, Turner N, et al. (2020) The Maori and Pacific specific CREBRF variant and adult height. Int J Obes (Lond) 44(3): 748–752. doi:10.1038/s41366-019-0437-6
    OpenUrlCrossRef
  7. [7].↵
    Hanson RL, Safabakhsh S, Curtis JM, et al. (2019) Association of CREBRF variants with obesity and diabetes in Pacific Islanders from Guam and Saipan. Diabetologia 62(9): 1647–1652. doi:10.1007/s00125-019-4932-z
    OpenUrlCrossRef
  8. [8].↵
    Krishnan M, Murphy R, Okesene-Gafa KAM, et al. (2020) The Pacific-specific CREBRF rs373863828 allele protects against gestational diabetes mellitus in Maori and Pacific women with obesity. Diabetologia 63(10): 2169–2176. doi:10.1007/s00125-020-05202-8
    OpenUrlCrossRef
  9. [9].↵
    Lin M, Caberto C, Wan P, et al. (2020) Population-specific reference panels are crucial for genetic analyses: an example of the CREBRF locus in Native Hawaiians. Hum Mol Genet 29(13): 2275–2284. doi:10.1093/hmg/ddaa083
    OpenUrlCrossRef
  10. [10].↵
    Naka I, Furusawa T, Kimura R, et al. (2017) A missense variant, rs373863828-A (p.Arg457Gln), of CREBRF and body mass index in Oceanic populations. J Hum Genet 62(9): 847–849. doi:10.1038/jhg.2017.44
    OpenUrlCrossRef
  11. [11].↵
    Colditz GA, Willett WC, Rotnitzky A, Manson JE (1995) Weight gain as a risk factor for clinical diabetes mellitus in women. Ann Intern Med 122(7): 481–486. doi:10.7326/0003-4819-122-7-199504010-00001
    OpenUrlCrossRefPubMedWeb of Science
  12. [12].↵
    Taylor RW, Brooking L, Williams SM, et al. (2010) Body mass index and waist circumference cutoffs to define obesity in indigenous New Zealanders. Am J Clin Nutr 92(2): 390–397. doi:10.3945/ajcn.2010.29317
    OpenUrlAbstract/FREE Full Text
  13. [13].
    McAuley KA, Williams SM, Mann JI, Goulding A, Murphy E (2002) Increased risk of type 2 diabetes despite same degree of adiposity in different racial groups. Diabetes Care 25(12): 2360–2361. doi:10.2337/diacare.25.12.2360
    OpenUrlFREE Full Text
  14. [14].↵
    Linhart C, Naseri T, Lin S, et al. (2018) Continued increases in blood pressure over two decades in Samoa (1991-2013); around one-third of the increase explained by rising obesity levels. BMC Public Health 18(1): 1122. doi:10.1186/s12889-018-6016-2
    OpenUrlCrossRef
  15. [15].↵
    Loos RJ (2016) CREBRF variant increases obesity risk and protects against diabetes in Samoans. Nat Genet 48(9): 976–978. doi:10.1038/ng.3653
    OpenUrlCrossRef
  16. [16].↵
    DeFronzo RA (2004) Pathogenesis of type 2 diabetes mellitus. Med Clin North Am 88(4): 787-835, ix. doi:10.1016/j.mcna.2004.04.013
    OpenUrlCrossRefPubMedWeb of Science
  17. [17].↵
    Frayn KN (2001) Adipose tissue and the insulin resistance syndrome. Proc Nutr Soc 60(3): 375–380. doi:10.1079/pns200195
    OpenUrlCrossRefPubMedWeb of Science
  18. [18].↵
    Yaghootkar H, Scott RA, White CC, et al. (2014) Genetic evidence for a normal-weight “metabolically obese” phenotype linking insulin resistance, hypertension, coronary artery disease, and type 2 diabetes. Diabetes 63(12): 4369–4377. doi:10.2337/db14-0318
    OpenUrlAbstract/FREE Full Text
  19. [19].↵
    Huang LO, Rauch A, Mazzaferro E, et al. (2021) Genome-wide discovery of genetic loci that uncouple excess adiposity from its comorbidities. Nat Metab 3(2): 228–243. doi:10.1038/s42255-021-00346-2
    OpenUrlCrossRef
  20. [20].↵
    Arslanian KJ, Fidow UT, Atanoa T, et al. (2021) A missense variant in CREBRF, rs373863828, is associated with fat-free mass, not fat mass in Samoan infants. Int J Obes (Lond) 45(1): 45–55. doi:10.1038/s41366-020-00659-4
    OpenUrlCrossRef
  21. [21].↵
    Tanner C, Boocock J, Stahl EA, et al. (2017) Population-Specific Resequencing Associates the ATP-Binding Cassette Subfamily C Member 4 Gene With Gout in New Zealand Maori and Pacific Men. Arthritis Rheumatol 69(7): 1461–1469. doi:10.1002/art.40110
    OpenUrlCrossRef
  22. [22].↵
    Tripathy D, Wessman Y, Gullstrom M, Tuomi T, Groop L (2003) Importance of obtaining independent measures of insulin secretion and insulin sensitivity during the same test: results with the Botnia clamp. Diabetes Care 26(5): 1395–1401. doi:10.2337/diacare.26.5.1395
    OpenUrlAbstract/FREE Full Text
  23. [23].↵
    Mosteller RD (1987) Simplified calculation of body-surface area. N Engl J Med 317(17): 1098. doi:10.1056/NEJM198710223171717
    OpenUrlCrossRefPubMedWeb of Science
  24. [24].↵
    Levy JC, Matthews DR, Hermans MP (1998) Correct homeostasis model assessment (HOMA) evaluation uses the computer program. Diabetes Care 21(12): 2191–2192. doi:10.2337/diacare.21.12.2191
    OpenUrlFREE Full Text
  25. [25].↵
    Van Cauter E, Mestrez F, Sturis J, Polonsky KS (1992) Estimation of insulin secretion rates from C-peptide levels. Comparison of individual and standard kinetic parameters for C-peptide clearance. Diabetes 41(3): 368–377. doi:10.2337/diab.41.3.368
    OpenUrlAbstract/FREE Full Text
  26. [26].↵
    Lorenzo C, Wagenknecht LE, Rewers MJ, et al. (2010) Disposition index, glucose effectiveness, and conversion to type 2 diabetes: the Insulin Resistance Atherosclerosis Study (IRAS). Diabetes Care 33(9): 2098–2103. doi:10.2337/dc10-0165
    OpenUrlAbstract/FREE Full Text
  27. [27].↵
    Drucker DJ (2006) The biology of incretin hormones. Cell Metab 3(3): 153–165. doi:10.1016/j.cmet.2006.01.004
    OpenUrlCrossRefPubMedWeb of Science
  28. [28].↵
    Bergman RN, Phillips LS, Cobelli C (1981) Physiologic evaluation of factors controlling glucose tolerance in man: measurement of insulin sensitivity and beta-cell glucose sensitivity from the response to intravenous glucose. J Clin Invest 68(6): 1456–1467. doi:10.1172/jci110398
    OpenUrlCrossRefPubMedWeb of Science
  29. [29].↵
    Curry DL, Bennett LL, Grodsky GM (1968) Dynamics of insulin secretion by the perfused rat pancreas. Endocrinology 83(3): 572–584. doi:10.1210/endo-83-3-572
    OpenUrlCrossRefPubMedWeb of Science
  30. [30].↵
    Swinburn BA, Ley SJ, Carmichael HE, Plank LD (1999) Body size and composition in Polynesians. Int J Obes Relat Metab Disord 23(11): 1178–1183. doi:10.1038/sj.ijo.0801053
    OpenUrlCrossRefPubMedWeb of Science
  31. [31].↵
    Geloneze B, Vasques AC, Stabe CF, et al. (2009) HOMA1-IR and HOMA2-IR indexes in identifying insulin resistance and metabolic syndrome: Brazilian Metabolic Syndrome Study (BRAMS). Arq Bras Endocrinol Metabol 53(2): 281–287. doi:10.1590/s0004-27302009000200020
    OpenUrlCrossRefPubMed
  32. [32].↵
    Bruttomesso D, Pianta A, Mari A, et al. (1999) Restoration of early rise in plasma insulin levels improves the glucose tolerance of type 2 diabetic patients. Diabetes 48(1): 99–105. doi:10.2337/diabetes.48.1.99
    OpenUrlAbstract
  33. [33].↵
    Taylor R, Al-Mrabeh A, Zhyzhneuskaya S, et al. (2018) Remission of Human Type 2 Diabetes Requires Decrease in Liver and Pancreas Fat Content but Is Dependent upon Capacity for beta Cell Recovery. Cell Metab 28(4): 547–556 e543. doi:10.1016/j.cmet.2018.07.003
    OpenUrlCrossRefPubMed
  34. [34].↵
    Polonsky KS, Rubenstein AH (1984) C-peptide as a measure of the secretion and hepatic extraction of insulin. Pitfalls and limitations. Diabetes 33(5): 486–494. doi:10.2337/diab.33.5.486
    OpenUrlFREE Full Text
  35. [35].↵
    Bojsen-Moller KN, Lundsgaard AM, Madsbad S, Kiens B, Holst JJ (2018) Hepatic Insulin Clearance in Regulation of Systemic Insulin Concentrations-Role of Carbohydrate and Energy Availability. Diabetes 67(11): 2129–2136. doi:10.2337/db18-0539
    OpenUrlAbstract/FREE Full Text
  36. [36].↵
    Nauck MA, Homberger E, Siegel EG, et al. (1986) Incretin effects of increasing glucose loads in man calculated from venous insulin and C-peptide responses. J Clin Endocrinol Metab 63(2): 492–498. doi:10.1210/jcem-63-2-492
    OpenUrlCrossRefPubMedWeb of Science
  37. [37].↵
    Kahn SE, Prigeon RL, McCulloch DK, et al. (1993) Quantification of the relationship between insulin sensitivity and beta-cell function in human subjects. Evidence for a hyperbolic function. Diabetes 42(11): 1663–1672. doi:10.2337/diab.42.11.1663
    OpenUrlAbstract/FREE Full Text
  38. [38].
    Bergman RN, Finegood DT, Kahn SE (2002) The evolution of beta-cell dysfunction and insulin resistance in type 2 diabetes. Eur J Clin Invest 32 Suppl 3: 35–45. doi:10.1046/j.1365-2362.32.s3.5.x
    OpenUrlCrossRefPubMedWeb of Science
  39. [39].↵
    Bergman RN, Ader M, Huecking K, Van Citters G (2002) Accurate assessment of beta-cell function: the hyperbolic correction. Diabetes 51 Suppl 1: S212–220. doi:10.2337/diabetes.51.2007.s212
    OpenUrlAbstract/FREE Full Text
  40. [40].↵
    Tiedge M (2014) Inside the pancreas: progress and challenges of human beta cell mass quantification. Diabetologia 57(5): 856–859. doi:10.1007/s00125-014-3206-z
    OpenUrlCrossRef
  41. [41].↵
    Meier JJ, Butler AE, Saisho Y, et al. (2008) Beta-cell replication is the primary mechanism subserving the postnatal expansion of beta-cell mass in humans. Diabetes 57(6): 1584–1594. doi:10.2337/db07-1369
    OpenUrlAbstract/FREE Full Text
  42. [42].↵
    Talchai C, Xuan S, Lin HV, Sussel L, Accili D (2012) Pancreatic beta cell dedifferentiation as a mechanism of diabetic beta cell failure. Cell 150(6): 1223–1234. doi:10.1016/j.cell.2012.07.029
    OpenUrlCrossRefPubMedWeb of Science
  43. [43].↵
    Porte D, Jr.., Kahn SE (2001) beta-cell dysfunction and failure in type 2 diabetes: potential mechanisms. Diabetes 50 Suppl 1: S160–163. doi:10.2337/diabetes.50.2007.s160
    OpenUrlFREE Full Text
  44. [44].↵
    DenBoer LM, Hardy-Smith PW, Hogan MR, Cockram GP, Audas TE, Lu R (2005) Luman is capable of binding and activating transcription from the unfolded protein response element. Biochem Biophys Res Commun 331(1): 113–119. doi:10.1016/j.bbrc.2005.03.141
    OpenUrlCrossRefPubMedWeb of Science
  45. [45].
    Liang G, Audas TE, Li Y, et al. (2006) Luman/CREB3 induces transcription of the endoplasmic reticulum (ER) stress response protein Herp through an ER stress response element. Mol Cell Biol 26(21): 7999–8010. doi:10.1128/MCB.01046-06
    OpenUrlAbstract/FREE Full Text
  46. [46].↵
    Audas TE, Li Y, Liang G, Lu R (2008) A novel protein, Luman/CREB3 recruitment factor, inhibits Luman activation of the unfolded protein response. Mol Cell Biol 28(12): 3952–3966. doi:10.1128/MCB.01439-07
    OpenUrlAbstract/FREE Full Text
  47. [47].↵
    Sano R, Reed JC (2013) ER stress-induced cell death mechanisms. Biochim Biophys Acta 1833(12): 3460–3470. doi:10.1016/j.bbamcr.2013.06.028
    OpenUrlCrossRefPubMedWeb of Science
  48. [48].↵
    Fonseca SG, Gromada J, Urano F (2011) Endoplasmic reticulum stress and pancreatic beta-cell death. Trends Endocrinol Metab 22(7): 266–274. doi:10.1016/j.tem.2011.02.008
    OpenUrlCrossRefPubMed
Back to top
Posted March 01, 2021.
Download PDF
Data/Code
Email
The CREBRF diabetes-protective rs373863828-A allele is associated with enhanced early insulin release in men of Māori and Pacific ancestry
Hannah J. Burden, Shannon Adams, Braydon Kulatea, Morag Wright-McNaughton, Danielle Sword, Jennifer J. Ormsbee, Conor Watene-O’Sullivan, Tony R. Merriman, Jennifer L. Knopp, J. Geoffrey Chase, Jeremy D. Krebs, Rosemary M. Hall, Lindsay D. Plank, Rinki Murphy, Peter R. Shepherd, Troy L. Merry
medRxiv 2021.02.27.21252567; doi: https://doi.org/10.1101/2021.02.27.21252567
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area