Abstract
Aims/hypothesis
The A (minor) allele of CREBRF rs373863828 has been associated with increased BMI and reduced risk of type 2 diabetes in the Samoan populations of Samoa and American Samoa. Our aim was to test rs373863828 for associations with BMI and the odds of type 2 diabetes, gout and chronic kidney disease (CKD) in Māori and Pacific (Polynesian) people living in Aotearoa/New Zealand.
Methods
Linear and logistic regression models were used to analyse the association of the A allele of CREBRF rs373863828 with BMI, log-transformed BMI, waist circumference, type 2 diabetes, gout and CKD in 2286 adults. The primary analyses were adjusted for age, sex, the first four genome-wide principal components and (where appropriate) BMI, waist circumference and type 2 diabetes. The primary analysis was conducted in ancestrally defined groups and association effects were combined using meta-analysis.
Results
For the A allele of rs373863828, the effect size was 0.038 (95% CI 0.022, 0.055, p = 4.8 × 10−6) for log-transformed BMI, with OR 0.59 (95% CI 0.47, 0.73, p = 1.9 × 10−6) for type 2 diabetes. There was no evidence for an association of genotype with variance in BMI (p = 0.13), and nor was there evidence for associations with serum urate (β = 0.012 mmol/l, pcorrected = 0.10), gout (OR 1.00, p = 0.98) or CKD (OR 0.91, p = 0.59).
Conclusions/interpretation
Our results in New Zealand Polynesian adults replicate, with very similar effect sizes, the association of the A allele of rs373863828 with higher BMI but lower odds of type 2 diabetes among Samoan adults living in Samoa and American Samoa.
Similar content being viewed by others
Introduction
A missense variant (rs373863828, p.Arg457Gln) in CREBRF (encoding cAMP-responsive element binding protein 3 [CREB3] regulatory factor) has been associated with BMI in the Samoan and American Samoan populations [1]. The minor allele (c.1370A, p.457Gln) was associated with a 1.36 kg/m2 greater BMI (~4 kg in body weight at 1.7 m tall) [1, 2]. In a small sample from the Kingdom of Tonga (n = 171) the A (minor) allele was associated with an increase of 3.1 kg/m2 in BMI [3]. The high A allele frequency of rs373863828 among Samoans in Samoa and American Samoa (minor allele frequency [MAF] 0.26) and in Tonga (MAF 0.15) [1, 3] is in contrast to an exceedingly rare frequency in other populations in the Genome Aggregation Database (http://gnomad.broadinstitute.org, accessed 1 October 2017; MAF 5.3 × 10−5 in East Asian individuals, 3.3 × 10−5 in South Asian individuals, 4.0 × 10−5 in European individuals and absent in 7651 African individuals) [4], and supports the hypothesis that rs373863828 is an important risk factor for obesity that is unique to Samoan and Tongan populations and possibly other Polynesian populations.
Unlike FTO-IRX3 and many other obesity risk variants in other populations, the BMI-increasing A allele of rs373863828 has been associated with lower odds of type 2 diabetes (OR 0.59–0.74 after adjustment for BMI) [1]. This contrasts with the established observational positive association between type 2 diabetes and increased BMI. ‘Favourable adiposity’ genetic variants have recently been identified, which have been associated with a higher BMI but lower risk of type 2 diabetes, hypertension and coronary artery disease in European populations [5,6,7]. For example, in a study that used data from the UK Biobank, a genetic score of 11 such alleles was associated with a higher BMI (+0.12 kg/m2) and higher body fat percentage (+0.30%) in the top 50% of people carrying the most favourable adiposity alleles, compared to the bottom 50%. However, for a given BMI, individuals carrying these alleles had reduced odds of type 2 diabetes (OR 0.84), hypertension (OR 0.94) and heart disease (OR 0.92) [8]. Such genetic variants are hypothesised to promote energy storage by expanding favourable subcutaneous adipose tissue depots, resulting in less accumulation of fat in non-adipose tissues and a consequent reduction in the risk of type 2 diabetes and other metabolic consequences of insulin resistance [8, 9]. In contrast, aside from increasing BMI and decreasing both fasting glucose and the risk of type 2 diabetes, the CREBRF rs373863828 A allele has not been associated with other features of insulin resistance such as hypertension, altered lipid levels or a higher HOMA-IR [1].
In Aotearoa/New Zealand, obesity and type 2 diabetes are both very prevalent in Māori and Pacific people [9] and are associated with other very prevalent metabolic-based conditions, in particular gout and chronic kidney disease (CKD) [10,11,12]. Moreover, increased BMI is a cause of increased urate levels [13], and diabetic kidney disease has been ascribed as the cause of renal failure in 69% of Polynesian people [14]. CREBRF rs373863828 has not previously been tested for an association with gout or CKD [1]. Polynesian populations include those from West Polynesia (Samoa, Tonga, Niue and Tokelau) and East Polynesia (Aotearoa/New Zealand Māori and Cook Island Māori). Given different pathogenic allele frequencies and different linkage disequilibrium structures between East and West Polynesian populations [15, 16], investigating CREBRF rs373863828 in other Polynesian population groups may provide further insights into its association with increasing BMI, yet apparent reduction in the risk of type 2 diabetes.
In this study, we tested the associations of CREBRF rs373863828 with BMI, waist circumference, type 2 diabetes, gout and CKD in people of Polynesian ancestry living in Aotearoa/New Zealand. In addition, based on the association of the FTO genotype with BMI variance in Europeans [17], we tested the association of rs373863828 with variance in log-transformed BMI to detect possible underlying genetic and/or environmental influences on phenotypic variability.
Methods
Study population
Individuals aged ≥16 years, primarily from the Auckland, Waikato and Christchurch regions of Aotearoa/New Zealand, were recruited to the Genetics of Gout, Diabetes and Kidney Disease in Aotearoa New Zealand Study [18]. A separate Māori sample set from the rohe (area) of the Ngāti Porou iwi (tribe) of the Tairāwhiti (east coast of the North Island) region was also included in the Aotearoa/New Zealand Māori analysis. This sample set was recruited in collaboration with the Ngāti Porou Hauora (Health Service) Charitable Trust. A Pukapuka Island sample set was recruited in collaboration with the Pukapuka Community of New Zealand Inc. in Mangere, South Auckland.
Information obtained at recruitment included age, sex, height, weight and waist circumference, as measured by trained assessors. BMI was calculated by dividing an individual’s weight by the square of their height in metres. Participants were also asked about the ancestry of each of their grandparents. Type 2 diabetes was ascertained by physician diagnosis, participant report or the use of glucose-lowering therapy. Blood biochemical measurements were performed at the Southern Community Laboratories (www.sclabs.co.nz). eGFR was derived from participants’ serum creatinine, age and sex using the Chronic Kidney Disease Epidemiology Collaboration equation [19]. Stage 4 and 5 CKD were defined by an eGFR of <30 ml min−1 [1.73 m]−2. Obesity was defined as a BMI of >32 kg/m2 [20]. Ethical approval for this study was given by the NZ Multi-Region Ethics Committee (MEC/05/10/130; MEC/10/09/092; MEC/11/04/036) and the Northern Y Region Health Research Ethics Committee (Ngāti Porou Hauora Charitable Trust study; NTY07/07/074). All participants provided written informed consent for the collection of samples and subsequent analysis.
Participants were separated into sample sets based on the self-reported Pacific nation of ancestry of their grandparents. Those participants who also reported non-Polynesian ancestry were grouped according to their Polynesian ancestry. This resulted in seven sample sets: Aotearoa/New Zealand Māori (n = 1296, including 270 people from the Ngāti Porou Hauora Charitable Trust study), Cook Island Māori (n = 205), Samoan (n = 387), Tongan (n = 181), Niuean (n = 47), Pukapukan (n = 75) and an ‘other’ Polynesian group (n = 271), which included individuals of Tahitian (n = 3), Tokelauan (n = 6) and Tuvaluan (n = 5) ancestry, along with individuals who self-reported grandparental ancestry from more than one Pacific nation (n = 257). Pukapuka is part of the Cook Islands, but is geographically located within West Polynesia. These analysis groups were further refined based on clustering of genome-wide principal component vectors one to four (details of calculation below), resulting in the exclusion of 176 people (ESM Fig. 1) for one of three reasons: (1) a clear mismatch between self-reported and principal component-defined ancestry (n = 16; potentially indicating an error in DNA or phenotype sample ID alignment during processing); (2) clustering outside of any of the principal component clusters (n = 26); and (3) individuals self-reported as Aotearoa/New Zealand Māori and Cook Island Māori who clustered within the European and between the European and Eastern Polynesian clusters (n = 135). The final groups used in all analyses totalled 2286 and were: Aotearoa/New Zealand Māori (n = 1154), Cook Island Māori (n = 197), Samoan (n = 378), Tongan (n = 175), Niuean (n = 47), Pukapukan (n = 70) and ‘other’ Polynesian (n = 265). Baseline characteristics for the final groupings are shown in Table 1.
Whole-genome microarray analysis
The Illumina Infinium CoreExome v24 BeadChip platform (Illumina, San Diego, CA, USA) was used to genotype participants for approximately 500,000 variants across the whole genome. Genotyping was performed at the University of Queensland (Centre for Clinical Genomics) for the Genetics of Gout, Diabetes and Kidney Disease in Aotearoa/New Zealand Study cohorts and at AgResearch (Invermay, Dunedin, New Zealand) for the Ngāti Porou Hauora Charitable Trust cohort. BeadChip genotyping batches were auto-clustered using GenomeStudio v2011.1 software (Illumina, San Diego, CA, USA). The Illumina GenomeStudio best practice guidelines and other quality-control protocols were applied [21, 22]. The genotyping batches were then merged and relevant quality-control steps were repeated for the full dataset.
Determination of principal components
Whole-genome principal component analysis vectors were calculated using a subset of 2858 ancestry-informative markers (as identified by Illumina) extracted from the CoreExome whole-genome genotypes. The SmartPCA (EIGENSOFT v6.0.1; www.hsph.harvard.edu/alkes-price/software/) [23] program was used, with an output of ten eigenvectors, no outlier removal and no population size limit. Individuals of non-Polynesian ancestry were included, and the first four vectors were plotted against each other to view the clustering of ancestral groupings (Asian, European, Eastern Polynesian and Western Polynesian) (ESM Fig. 1). The first four vectors, which explained 97.1% of the proportion of the variability within the sample sets, were included as covariates in the linear and logistic regression analyses to account for population stratification and cryptic relatedness.
CREBRF rs373863828 genotyping
rs373863828 was directly genotyped because this variant was not present on the CoreExome platform and we were unable to impute the region owing to the unavailability of Māori and Pacific reference haplotypes. A custom-designed TaqMan probe-set (Applied Biosystems, Foster City, CA, USA) was created for rs373863828 using a custom Python script (snp_design; DOI: https://doi.org/10.5281/zenodo.56250) to annotate the human genome build 37 reference sequence (ftp://ftp.ensembl.org/pub/grch37; accessed 1 August 2016) with rs373863828 and any surrounding SNPs (obtained from the NCBI dbSNP build 147 common SNP list; ftp://ftp.ncbi.nlm.nih.gov/snp). Forward primer: CAAGAGAGGATGCTGAGACCAT; reverse primer: ACCATGATGTAAGCCATTTTTCTGATACA; probe 1 (VIC): TGAGTGGAACCGAGATAC probe 2 (FAM): AGTGGAACCAAGATAC. Genotyping was performed using the LightCycler 480 Real-Time Polymerase Chain Reaction System (Roche Applied Science, Indianapolis, IN, USA) in 384 well plates. There was a 99% successful genotyping call rate. Re-genotyping of 25% of the sample set demonstrated 100% concordance.
Association testing
Analyses were performed using the R v3.3.2 statistical software (within RStudio v0.99.902; www.rstudio.com). A multivariable linear regression model was used to test for associations between the rs373863828 minor allele (c.1370A, p.457Gln) and the continuous variables log-transformed BMI (a Box-Cox normality plot for the BMI data yielded λ = −0.39), untransformed BMI, waist circumference and serum urate, with the β coefficient representing the estimated effect of each copy of the rs373863828 A allele. For binary outcomes (obesity, type 2 diabetes, gout and CKD), a multivariable binomial logistic regression model was used in a similar manner, with the allelic OR representing the estimated effect of each copy of the A allele. All analyses were adjusted for sex, age, the first four whole-genome principal component vectors and BMI, type 2 diabetes or waist circumference, where appropriate. Analyses were also performed in groups stratified by type 2 diabetes status. Each Polynesian population sample set was analysed separately, and the effects were combined using an inverse-variance-weighted fixed-effect meta-analysis. Heterogeneity among sample sets was assessed during the meta-analysis using Cochran’s heterogeneity (Q) statistic, with a random-effects analysis used when there was evidence of heterogeneity (p < 0.05). For the BMI, waist circumference and type 2 diabetes association analyses, p < 0.05 was set as the significance value, given the prior probability of detecting associations [1]. For the other outcomes (n = 4; serum urate, gout, CKD and variance in BMI), p < 0.0125 was set as the significance value, accounting for the fact that none of these outcomes had previously been tested for associations with rs373863828 [1].
The association analyses for the major outcomes (BMI, log e -transformed BMI and type 2 diabetes) were repeated using two different approaches to account for relatedness within the analysis group. First, a kinship coefficient matrix was calculated in PLINK v1.9 [24] using 40,156 independent autosomal markers in the pooled dataset. Participant pairs with a kinship coefficient (\( \widehat{\pi} \)) ≥0.125 (equivalent to first cousins or closer) were identified and one individual was excluded (at random) to create a maximal set of unrelated individuals (n = 1665) for analysis. Multivariate linear (or logistic) regression analyses testing the CREBRF rs373863828 association were performed in this group, adjusting for sex, age, the first four principal component vectors, membership of the seven previously described ancestral groups and BMI, type 2 diabetes or waist circumference, where appropriate. Second, for BMI, a linear mixed-model regression analysis, which included the kinship coefficient matrix calculated above to adjust for relatedness, was performed in the entire analysis group (not excluding related individuals) using GenABEL (www.rdocumentation.org/packages/GenABEL/versions/1.8-0) [25, 26] and the same adjusting variables as above.
Models of inheritance were investigated by formulating the genotype predictor in the linear and logistic regression models (adjusting for age, sex and principal components 1–4) in different ways: additive model (0, 1, 2); dominant model (0, 1, 1); or recessive model (0, 0, 1). For each model the coded genotypes were treated as a continuous variable, producing a single effect estimate (β or OR). A model-selection tool (Akaike’s information criterion [AIC] [27]) was used to select the most likely model. Smaller AIC values indicate better models, but where the difference was <2 the simplest model was chosen.
Power
Based on previously reported estimates [1] and α = 0.05, the power to detect an effect size of 1.36 kg/m2 per A allele was >80% in the combined Māori and Pacific Island sample set for a MAF of ≥0.15 (ESM Fig. 2) The power to detect a moderate protective effect for type 2 diabetes (OR 0.59) of the minor allele was >90% for a MAF of ≥0.10 (ESM Fig. 2). Power calculations for gout and CKD as outcomes with α = 0.0125 (to account for serum urate, gout, CKD and variance in BMI, not previously tested [1]) showed that the power was adequate only to detect effect sizes of OR ≥1.75 for gout and CKD and ≥0.032 mmol/l for serum urate (ESM Fig. 2).
Testing for an association of rs373863828 with variance in log-transformed BMI
The association of a genetic variant with variance in phenotype can detect a locus interacting in a non-additive way without prior knowledge of the interacting factor (environmental, intrinsic, genetic). Testing for an association of rs373863828 with variance in log-transformed BMI was performed as previously described [28]. The variable used as a measure of variance was produced from residuals obtained from cohort-specific analyses regressing age, age2 and age by sex. An independent ranked inverse normal transformation of absolute residuals generated Z scores, with squared Z scores (Z2) being the variance variable. To account for the influence of rs373863828 mean effect on the variance, the mean log-transformed BMI (per genotype) was subtracted from the log-transformed BMI of each participant and the Z scores were recalculated. Linear models associating CREBRF genotype with both the unadjusted and mean-effect adjusted variance Z scores were used. This analysis was also performed on the Samoan rs373863828 genotype data [1], with age, age2, age by sex and polity as the adjusting variables in the Z score calculation steps. The Aotearoa/New Zealand sample sets and the two Samoan sample sets [1] were combined using an inverse-variance weighted fixed-effect meta-analysis [28].
Results
Prevalence of rs373863828
The allele and genotype frequencies of rs373863828 in each Polynesian sample set are presented in Table 1 and Fig. 1. There was no evidence of deviation from Hardy–Weinberg equilibrium in the various sample sets (Table 1; p for all datasets was ≥0.049). The relative frequencies of the minor (c.1370A) allele differed among the Polynesian groups, with the Samoan (MAF 0.236) and Pukapukan (MAF 0.243) groups exhibiting the highest frequencies and the Niuean (MAF 0.096) group exhibiting the lowest frequency (Table 1).
Association analysis of rs373863828 with adiposity measures
A fixed-effect meta-analysis of the Polynesian samples showed a significant association of rs373863828 with log-transformed BMI (β = 0.038, p = 4.8 × 10−6), with no evidence of heterogeneity among the sample sets (p = 0.19) (Table 2, Fig. 2; ESM Table 1). Association analysis with untransformed BMI revealed similar results (Table 2, Fig. 2; ESM Table 2). Analysis of the pooled dataset removing related individuals yielded very similar effect sizes (β = 0.037 [95% CI 0.017, 0.057], p = 2.3 × 10−4 for log-transformed BMI and β = 1.38 kg/m2 [95% CI 0.66, 2.09], p = 1.6 × 10−4 for untransformed BMI). Analysis adjusting for relatedness using GenABEL also yielded very similar effect sizes (β = 0.038 [95% CI 0.021, 0.054], p = 1.2 × 10−5 for log-transformed BMI and β = 1.34 kg/m2 [95% CI 0.73, 1.94], p = 1.6 × 10−5 for untransformed BMI). The proportion of variance explained by rs373863828 in the pooled dataset was 1.1%. One copy of the minor allele was sufficient to confer the effect (Table 2; β = 1.80 for the heterozygote group and β = 1.49 for the A allele homozygote group compared with major allele homozygotes in untransformed BMI analysis). This was supported by AIC analysis, where the dominant model had a difference of 3.3 less than the additive model. In the combined group there was an association with higher odds of obesity (Table 2; OR 1.33, p = 8 × 10−4 for BMI >32 kg/m2 and OR 1.54, p = 1 × 10−5 for BMI >40 kg/m2). There was a significant association with untransformed BMI in both the type 2 diabetes-positive and -negative groups (β = 1.62 kg/m2, p = 2.4 × 10−3 and β = 1.21 kg/m2, p = 5.0 × 10−4, respectively, with no evidence of heterogeneity in effect size [pHet = 0.48]), indicating no evidence of a type 2 diabetes-specific effect on BMI. There was no indication of sex-specific effects (Table 2, Fig. 3) and a sex-by-rs373863828 interaction analysis of the pooled sample set showed no evidence of an interaction with BMI (p = 0.67).
A fixed-effect meta-analysis showed evidence of an association between rs373863828 and increased waist circumference in the full Polynesian sample set (β = 2.98 cm per minor allele, p = 1.3 × 10−5), with no evidence of heterogeneity among sample sets (p = 0.089) (Fig. 2; ESM Table 3). Adjustment of the waist circumference analysis for BMI abrogated the association with waist circumference (Table 2; β decreased from 2.98 to 0.66, p = 0.092).
Association analysis of rs373863828 with type 2 diabetes
A fixed-effect meta-analysis of the various sample sets revealed a significant association with reduced odds of type 2 diabetes (OR 0.65, p = 3.4 × 10−5) that was strengthened after adjustment for BMI (OR 0.59, p = 1.9 × 10−6), with no evidence of heterogeneity among sample sets (p = 0.089) (Table 2, Fig. 4; ESM Table 4). Association analysis of the pooled dataset with related individuals removed yielded very similar effect sizes (OR 0.56 [95% CI 0.43, 0.72], p = 7.9 × 10−6). Similar to the observation for BMI, one copy of the minor allele appeared to be sufficient to confer the effect (Table 2: OR 0.55 for the heterozygote group compared with the major allele homozygote group). The AIC analysis showed the dominant model to be similar to the additive model (0.21 less). Therefore, the mode of inheritance for type 2 diabetes was concluded to be more consistent with an additive model.
We combined all available data. A fixed-effect meta-analysis that included the data presented here and data from the previously published Samoan [1] and Tongan [3] studies yielded β = 1.43 kg/m2 (95% CI 1.17, 1.68), p = 3.8 × 10−28 for BMI and OR 0.62 (95% CI 0.55, 0.70), p = 1.7 × 10−14 for type 2 diabetes, with no evidence of heterogeneity.
Association analysis of rs373863828 with serum urate, gout and CKD
We tested for associations with serum urate (Table 2; ESM Fig. 3). Unadjusted, there was evidence of an association of the A allele with higher levels (β = 0.015 mmol/l, p = 0.005, pcorrected = 0.020); however, this effect was mitigated after adjusting for BMI (β = 0.012 mmol/l, p = 0.026, Pcorrected = 0.10). Similarly there was no evidence of an association with gout (OR 1.00, p = 0.98) (Table 2; ESM Table 5, ESM Fig. 3). There was no statistically significant (p < 0.0125) evidence of an association with CKD either before (OR 0.72, p = 0.030, pcorrected = 0.12) or after adjustment for type 2 diabetes and BMI (OR 0.91, p = 0.59) (Table 2; ESM Table 6, ESM Fig. 3). Finally, analysis in type 2 diabetes-positive and -negative groups yielded evidence of an association only for serum urate in the type 2 diabetes-negative group (Table 2). However there was no evidence of a difference in effect size between the type 2 diabetes-positive and -negative groups for serum urate (p > 0.13).
Association analysis of rs373863828 with variance in phenotype
The rs373863828 variant was not significantly associated with log-transformed BMI variance at the CREBRF locus for any of the Aotearoa/New Zealand sample sets or for the two previously published Samoan cohorts [1]. Fixed-effect meta-analysis of the Aotearoa/New Zealand (n = 2282), 1990s Samoan (n = 1020) and discovery Samoan (n = 1876) cohorts showed no evidence of an association of rs373863828 with variance in log-transformed BMI in either the unadjusted (β = 0.065, p = 0.05) or adjusted models (β = 0.050, p = 0.13) (ESM Table 7, ESM Fig. 4).
Discussion
In this study, we confirmed the association of CREBRF rs373863828 with a higher BMI (1.38 kg/m2 per minor allele) and larger waist circumference (2.98 cm), but lower risk of type 2 diabetes (OR 0.59) in adults of Polynesian (Māori and Pacific) ancestry living in Aotearoa/New Zealand. These results are very similar to those from Samoans living in Samoa and American Samoa, where each copy of the minor allele was associated with a 1.36 kg/m2 increase in BMI and an OR for type 2 diabetes of 0.59 [1]. The per cent variation in BMI explained was also similar (1.9% in the Samoan discovery sample and 1.1% in the Samoan replication sample, vs 1.1% in our Aotearoa/New Zealand sample), suggesting that the main effect of the rs373863828 A allele is relatively impervious to environment. Consistent also with the Samoan data [1], adjustment for BMI strengthened the association with type 2 diabetes in the Aotearoa/New Zealand sample set (OR decreased from 0.65 to 0.59 [Table 2]; Samoan discovery OR decreased from 0.64 to 0.59 [1]; Samoan replication OR decreased from 0.83 to 0.74 [1]). However, adjustment for BMI diminished the association with waist circumference (Table 2), indicating that the waist circumference association was driven by overall body mass distribution rather than central adiposity.
Precisely how these genetic epidemiological findings relate to the actual CREBRF-mediated molecular pathogenesis of obesity and type 2 diabetes is unclear in the absence of detailed knowledge of the molecular pathways involving CREBRF and in the absence of genetic association data with detailed body-composition measures as outcomes. Most population genetic variants associated with generalised obesity are also associated with insulin resistance, hypertension, dyslipidaemia and type 2 diabetes, compatible with the degree of adiposity. However, some genetic variants are associated with a higher BMI and percentage body fat, but lower odds of type 2 diabetes along with lower insulin resistance, hypertension, circulating triacylglycerols and LDL-cholesterol levels. These variants are known as ‘favourable adiposity’ variants [8] due to higher subcutaneous-to-visceral adipose tissue that suggests preferential fat storage away from visceral organs. The association of the A allele of rs373863828 with a higher BMI but reduced odds of type 2 diabetes is not entirely compatible with ‘favourable adiposity’ due to the lack of an association with hypertension or lipids [1]. Consistent with the reduced odds of type 2 diabetes, there was a weak association of the A allele with increased insulin sensitivity by HOMA-IR in the Samoan and American Samoan populations [1]; HOMA-IR data were not available for the current study. Cellular bioenergetics models have shown that the rs373863828 A allele promotes lipid and triacylglycerol storage in adipocytes at a reduced energy cost, suggesting that the metabolic activity of CREBRF in fat is important [1]. Detailed clinical studies are required to clarify whether visceral and subcutaneous body fat storage depots are altered among carriers of this variant.
It is notable that, from what is currently understood about the physiological role of CREBRF, it has no obvious role in regulating appetite, which is seen at FTO-IRX3 and other loci regulating BMI in Europeans. However, CREBRF is widely expressed (www.gtexportal.org), including throughout the brain, and CREBRF is known to regulate the CREB3/Luman protein, which is localised to the endoplasmic reticulum and plays an important role in axonal regeneration [29]. Interestingly, the CREB3/Luman protein was identified through its association with herpes simplex virus-related host cell factor 1, which has led to the hypothesis that CREB3/Luman may play a role in emergence of the virus from latency [30]. It will be important to explore the relationship between CREBRF expression in the hypothalamic nuclei and the role of CREB3/Luman in intra-axonal translation and retrograde trafficking to promote neuronal survival in response to viral stimuli, and how this might relate to BMI and type 2 diabetes.
A recent study reported an association of the CREBRF rs373863828 A allele with higher BMI in a pooled multi-ethnic sample of 4572 New Zealand children of Māori, Pacific, European and Asian ethnicity [31]. Given that the association analysis was pooled and population-specific association analyses were not reported, we are unable to directly compare our data with that study. The A allele was reported at a prevalence of 0.015 in European children and 0.011 in Asian children [31], in contrast to the very low prevalence reported in the Genome Aggregation Database (MAF 5.3 × 10−5 in East Asian individuals, 3.3 × 10−5 in South Asian individuals and 4.0 × 10−5 in European individuals).
Our study has confirmed that the presence of each additional CREBRF rs373863828 A allele is associated with increased BMI yet reduced odds of type 2 diabetes in adults of Māori and Pacific (Polynesian) ancestry living in Aotearoa/New Zealand. While the prevalence of both obesity and type 2 diabetes is increased among New Zealand Māori and those of Pacific ethnicity [9] compared with New Zealand Europeans, our study confirms that Pacific-specific genetic variation underpins some of the interindividual heterogeneity observed in the discordant manifestations of obesity and type 2 diabetes. This study supports the need to conduct comprehensive gene–phenotype studies in populations currently under-represented in genomic studies, among which different genetically segregating pathways linking obesity and type 2 diabetes clearly exist. Such studies are not only important for these populations per se, but are also important in illuminating the molecular biology of the pathogenesis of metabolic disease in the wider human population and have the potential to lead to novel clinical interventions.
Data availability
The dataset generated during the current study is not publicly available owing to consent restrictions, but can be requested from the corresponding author under an appropriate arrangement.
Abbreviations
- AIC:
-
Akaike’s information criterion
- CKD:
-
Chronic kidney disease
- CREB3:
-
cAMP-responsive element binding protein 3
- CREBRF:
-
cAMP-responsive element binding protein 3 regulatory factor
- MAF:
-
Minor allele frequency
References
Minster RL, Hawley NL, Su CT et al (2016) A thrifty variant in CREBRF strongly influences body mass index in Samoans. Nat Genet 48:1049–1054
Loos RJ (2016) CREBRF variant increases obesity risk and protects against diabetes in Samoans. Nat Genet 48:976–978
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:847–849
Lek M, Karczewski KJ, Minikel EV et al (2016) Analysis of protein-coding genetic variation in 60,706 humans. Nature 536:285–291
Kilpelainen TO, Zillikens MC, Stancakova A et al (2011) Genetic variation near IRS1 associates with reduced adiposity and an impaired metabolic profile. Nat Genet 43:753–760
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:4369–4377
Scott RA, Fall T, Pasko D et al (2014) Common genetic variants highlight the role of insulin resistance and body fat distribution in type 2 diabetes, independent of obesity. Diabetes 63:4378–4387
Yaghootkar H, Lotta LA, Tyrrell J et al (2016) Genetic evidence for a link between favorable adiposity and lower risk of type 2 diabetes, hypertension, and heart disease. Diabetes 65:2448–2460
Virtue S, Vidal-Puig A (2010) Adipose tissue expandability, lipotoxicity and the metabolic syndrome—an allostatic perspective. Biochim Biophys Acta 1801:338–349
Winnard D, Wright C, Jackson G et al (2012) Gout, diabetes and cardiovascular disease in the Aotearoa New Zealand adult population: co-prevalence and implications for clinical practice. N Z Med J 126:53–64
Winnard D, Wright C, Taylor WJ et al (2012) National prevalence of gout derived from administrative health data in Aotearoa New Zealand. Rheumatology (Oxford) 51:901–909
Hollis-Moffatt JE, Phipps-Green AJ, Chapman B et al (2012) The renal urate transporter SLC17A1 locus: confirmation of association with gout. Arthritis Res Ther 14:R92
Robinson PC, Choi HK, Do R, Merriman TR (2016) Insight into rheumatological cause and effect through the use of Mendelian randomization. Nat Rev Rheumatol 12:486–496
Australia and New Zealand Dialysis & Transplant Registry (2015) Available from www.anzdata.org.au/v1/report_2015.html. Accessed 26 Sept 2017
Friedlaender JS, Friedlaender FR, Reed FA et al (2008) The genetic structure of Pacific islanders. PLoS Genet 4:e19
Terrell JE (2010) Social network analysis of the genetic structure of Pacific islanders. Ann Hum Genet 74:211–232
Yang J, Loos RJ, Powell JE et al (2012) FTO genotype is associated with phenotypic variability of body mass index. Nature 490:267–272
Phipps-Green AJ, Merriman ME, Topless R et al (2016) Twenty-eight loci that influence serum urate levels: analysis of association with gout. Ann Rheum Dis 75:124–130
Levey AS, Stevens LA, Schmid CH et al (2009) A new equation to estimate glomerular filtration rate. Ann Intern Med 150:604–612
Swinburn BA, Ley SJ, Carmichael HE, Plank LD (1999) Body size and composition in Polynesians. Int J Obes Relat Metab Disord 23:1178–1183
Guo Y, He J, Zhao S et al (2014) Illumina human exome genotyping array clustering and quality control. Nat Protoc 9:2643–2662
Illumina (2014) Infinium genotyping data analysis: a guide for analyzing infinium genotyping data using the GenomeStudio genotyping module. Available from https://www.illumina.com/Documents/products/technotes/technote_infinium_genotyping_data_analysis.pdf
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38:904–909
Purcell S, Neale B, Todd-Brown K et al (2007) PLINK: a toolset for whole-genome association and population-based linkage analysis. Am J Hum Genet 81:559–575
Aulchenko YS, Ripke S, Isaacs A, van Dujin CM (2007) GenABEL: an R library for genome-wide association analysis. Bioinformatics 23:1294–1296
Chen WM, Abecasis GR (2007) Family-based association tests for genome-wide association scans. Am J Hum Genet 81:913–926
Akaike H (1974) A new look at the statistical model identification. IEEE Trans Autom Control 19:716–723
Topless RK, Flynn TJ, Cadzow M et al (2015) Association of SLC2A9 genotype with phenotypic variability of serum urate in pre-menopausal women. Front Genet 6:313
Ying Z, Misra V, Verge VM (2014) Sensing nerve injury at the axonal ER: activated Luman/CREB3 serves as a novel axonally synthesized retrograde regeneration signal. Proc Natl Acad Sci U S A 111:16142–16147
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:7999–8010
Berry S, Walker C, Ly K et al (2017) Widespread prevalence of a CREBRF variant amongst Māori and Pacific children is associated with weight and height in early childhood. Int J Obes. https://doi.org/10.1038/ijo.2017.230
Kampstra P (2008) Beanplot: a boxplot alternative for visual comparison of distributors. J Stat Softw 28:1–9
Acknowledgements
The authors sincerely thank the participants for generously donating their time and information to this study. The authors would like to thank J. Drake (Department of Rheumatology, Canterbury District Health Board, Christchurch, New Zealand), J. de Kwant, R. Laurence, C. Franklin and M. House (all Department of Medicine, University of Auckland, Auckland, New Zealand), N. Aupouri, R. Akuhata and C. Ford (all of Ngāti Porou Hauora Charitable Trust, Te Puia Springs, New Zealand) and G. Sexton (Counties Manukau District Health Board, Auckland, New Zealand) for recruitment.
Funding
The Health Research Council of New Zealand (grant no. 08/075, 10/548, 11/1075, 14/527) and the Maurice Wilkins Centre funded the New Zealand component of this study, and the National Institutes of Health funded the Samoa and American Samoa components (grant no. R01-HL093093 and R01-HL133040).
Author information
Authors and Affiliations
Contributions
MK, TJM, PRS, RM and TRM contributed to the design of the study. OD, LM, JdZ, LKS, ND, JHH, NR, TN, MSR, RD, STM and SV contributed to data collection, and RKT, LY, JMDT, WWHE, DEW, RLM, PW, DG and ANS contributed to data analysis and interpretation. MK, RM and TRM drafted the manuscript and all of the other authors reviewed it. The manuscript was approved by all authors. TRM is the guarantor of this work.
Corresponding author
Ethics declarations
The authors declare that there is no duality of interest associated with this manuscript.
Electronic supplementary material
ESM
(PDF 872 kb)
Rights and permissions
About this article
Cite this article
Krishnan, M., Major, T.J., Topless, R.K. et al. Discordant association of the CREBRF rs373863828 A allele with increased BMI and protection from type 2 diabetes in Māori and Pacific (Polynesian) people living in Aotearoa/New Zealand. Diabetologia 61, 1603–1613 (2018). https://doi.org/10.1007/s00125-018-4623-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00125-018-4623-1