Association study of CREBRF missense variant (rs373863828:G > A; p.Arg457Gln) with levels of serum lipid profile in the Pacific populations

Associations between fasting glucose rate-of-change and the missense variant, rs373863828, in an adult Samoan cohort

BACKGROUND

The A allele of rs373863828 in CREB3 regulatory factor is associated with high Body Mass Index, but lower odds of type 2 diabetes. These associations have been replicated elsewhere, but to date all studies have been cross-sectional. Our aims were (1) to describe the development of type 2 diabetes and change in fasting glucose between 2010 and 2018 among a longitudinal cohort of adult Samoans without type 2 diabetes or who were not using diabetes medications at baseline, and (2) to examine associations between fasting glucose rate-of-change (mmol/L per year) and the A allele of rs373863828.

METHODS

We describe and test differences in fasting glucose, the development of type 2 diabetes, body mass index, age, smoking status, physical activity, urbanicity of residence, and household asset scores between 2010 and 2018 among a cohort of n = 401 adult Samoans, selected to have a ~2:2:1 ratio of GG:AG: AA rs373863828 genotypes. Multivariate linear regression was used to test whether fasting glucose rate-of-change was associated with rs373863828 genotype, and other baseline variables.

RESULTS

By 2018, fasting glucose and BMI significantly increased among all genotype groups, and a substantial portion of the sample developed type 2 diabetes mellitus. The A allele was associated with a lower fasting glucose rate-of-change (β = -0.05 mmol/L/year per allele, p = 0.058 among women; β = -0.004 mmol/L/year per allele, p = 0.863 among men), after accounting for baseline variables. Mean fasting glucose and mean BMI increased over an eight-year period and a substantial number of individuals developed type 2 diabetes by 2018. However, fasting glucose rate-of-change, and type 2 diabetes development was lower among females with AG and AA genotypes.

CONCLUSIONS

Further research is needed to understand the effect of the A allele on fasting glucose and type 2 diabetes development. Based on our observations that other risk factors increased over time, we advocate for the continued promotion for diabetes prevention and treatment programming, and the reduction of modifiable risk factors, in this setting.

Multivariate analysis of a missense variant in CREBRF reveals associations with measures of adiposity in people of Polynesian ancestries

The minor allele of rs373863828, a missense variant in CREB3 Regulatory Factor, is associated with several cardiometabolic phenotypes in Polynesian peoples. To better understand the variant, we tested the association of rs373863828 with a panel of correlated phenotypes (body mass index [BMI], weight, height, HDL cholesterol, triglycerides, and total cholesterol) using multivariate Bayesian association and network analyses in a Samoa cohort (n = 1632), Aotearoa New Zealand cohort (n = 1419), and combined cohort (n = 2976). An expanded set of phenotypes (adding estimated fat and fat-free mass, abdominal circumference, hip circumference, and abdominal-hip ratio) was tested in the Samoa cohort (n = 1496). In the Samoa cohort, we observed significant associations (log10 Bayes Factor [BF] ≥ 5.0) between rs373863828 and the overall phenotype panel (8.81), weight (8.30), and BMI (6.42). In the Aotearoa New Zealand cohort, we observed suggestive associations (1.5 < log10 BF < 5) between rs373863828 and the overall phenotype panel (4.60), weight (3.27), and BMI (1.80). In the combined cohort, we observed concordant signals with larger log10 BFs. In the Samoa-specific expanded phenotype analyses, we also observed significant associations between rs373863828 and fat mass (5.65), abdominal circumference (5.34), and hip circumference (5.09). Bayesian networks provided evidence for a direct association of rs373863828 with weight and indirect associations with height and BMI.

The protective effect of rs373863828 on type 2 diabetes does not operate through a body composition pathway in adult Samoans

OBJECTIVE

The aim of this study was to understand whether the paradoxical association of missense variant rs373863828 in CREB3 regulatory factor (CREBRF) with higher BMI but lower odds of diabetes is explained by either metabolically favorable body fat distribution or greater fat-free mass.

METHODS

This study explored the association of the minor allele with dual-energy x-ray absorptiometry-derived body composition in n = 421 Samoans and used path analysis to examine the mediating role of fat and fat-free mass on the relationship between rs373863828 and fasting glucose.

RESULTS

Among females, the rs373863828 minor A allele was associated with greater BMI. There was no association of genotype with percent body fat, visceral adiposity, or fat distribution in either sex. In both females and males, lean mass was greater with each A allele: 2.16 kg/copy (p = 0.0001) and 1.73 kg/copy (p = 0.02), respectively. Path analysis showed a direct negative effect of rs373863828 genotype on fasting glucose (p = 0.004) consistent with previous findings, but also an indirect positive effect on fasting glucose operating through fat-free mass (p = 0.027).

CONCLUSIONS

The protective effect of rs373863828 in CREBRF, common among Pacific Islanders, on type 2 diabetes does not operate through body composition. Rather, the variant's effects on body size/composition and fasting glucose likely operate via different, tissue-specific mechanisms.

Generation of a genetically modified pig model with CREBRFR457Q variant

Obesity is among the strongest risk factors for type 2 diabetes (T2D). The CREBRF missense allele rs373863828 (p. Arg457Gln, p. R457Q) is associated with increased body mass index but reduced risk of T2D in people of Pacific ancestry. To investigate the functional consequences of the CREBRF variant, we introduced the corresponding human mutation R457Q into the porcine genome. The CREBRF^R457Q pigs displayed dramatically increased fat deposition, which was mainly distributed in subcutaneous adipose tissue other than visceral adipose tissue. The CREBRF^R457Q variant promoted preadipocyte differentiation. The increased differentiation capacity of precursor adipocytes conferred pigs the unique histological phenotype that adipocytes had a smaller size but a greater number in subcutaneous adipose tissue (SAT) of CREBRF^R457Q variant pigs. In addition, in SAT of CREBRF^R457Q pigs, the contents of the peroxidative metabolites 4-hydroxy-nonenal and malondialdehyde were significantly decreased, while the activity of antioxidant enzymes, such as glutathione peroxidase, superoxide dismutase, and catalase, was increased, which was in accordance with the declined level of the reactive oxygen species (ROS) in CREBRF^R457Q pigs. Together, these data supported a causal role of the CREBRF^R457Q variant in the pathogenesis of obesity, partly via adipocyte hyperplasia, and further suggested that reduced oxidative stress in adipose tissue may mediate the relative metabolic protection afforded by this variant despite the related obesity.

Population-specific reference panels are crucial for genetic analyses: an example of the CREBRF locus in Native Hawaiians

Statistical imputation applied to genome-wide array data is the most cost-effective approach to complete the catalog of genetic variation in a study population. However, imputed genotypes in underrepresented populations incur greater inaccuracies due to ascertainment bias and a lack of representation among reference individuals, further contributing to the obstacles to study these populations. Here we examined the consequences due to the lack of representation by genotyping in a large number of self-reported Native Hawaiians (N = 3693) a functionally important, Polynesian-specific variant in the CREBRF gene, rs373863828. We found the derived allele was significantly associated with several adiposity traits with large effects (e.g. ~ 1.28 kg/m2 per allele in body mass index as the most significant; P = 7.5 × 10-5), consistent with the original findings in Samoans. Due to the current absence of Polynesian representation in publicly accessible reference sequences, rs373863828 or its proxies could not be tested through imputation using these existing resources. Moreover, the association signals at the entire CREBRF locus could not be captured by alternative approaches, such as admixture mapping. In contrast, highly accurate imputation can be achieved even if a small number (<200) of internally constructed Polynesian reference individuals were available; this would increase sample size and improve the statistical evidence of associations. Taken together, our results suggest the alarming possibility that lack of representation in reference panels could inhibit discovery of functionally important loci such as CREBRF. Yet, they could be easily detected and prioritized with improved representation of diverse populations in sequencing studies.

A missense variant in CREBRF, rs373863828, is associated with fat-free mass, not fat mass in Samoan infants

BACKGROUND/OBJECTIVES

In Samoa, where 80% of the adult population is living with obesity, understanding the determinants of adiposity and growth during infancy may inform prevention efforts. We examined the association of a missense variant, rs373863828, in the CREBRF gene with body composition in Samoan infants. Adults with one or more copies of the rs373863828 minor allele (A) have higher odds of obesity, based on body-mass index (BMI), but paradoxically decreased odds of diabetes compared to those without the allele. Our study may offer novel insight into the natural history and pathogenesis of this unexpected relationship.

SUBJECTS/METHODS

In a prospective study, we measured body composition in early infancy, and at 2- and 4-months of age using anthropometry and dual-energy x-ray absorptiometry (DXA). We genotyped subjects at the CREBRF rs373863828 locus and compared infants with (AA/AG) and without (GG) the variant. In longitudinal analyses, we calculated the absolute change in each outcome from the early infant to the 4-month assessment, adjusting for baseline and other covariates.

RESULTS

In cross-sectional analyses, there was no significant difference in infant BMI or fat mass by genotype. After adjusting for covariates, infants with the variant had 4.0 ± 1.8 g more bone mass (p = 0.026) and 210.9 ± 79.6 g more lean mass (p = 0.009) at 4-months and accumulated 176.9 ± 73.0 g more lean mass between the early infant and 4-month assessment (p = 0.017).

CONCLUSIONS

The CREBRF rs373863828 minor allele (A) was not associated with increased BMI or adiposity in Samoan infants, but instead with increased lean and bone mass. Our findings suggest that lean (i.e., muscle) and bone mass accretion should be explored as pathways to explain the "protective" effect of the CREBRF variant against diabetes.

Exploring the Paradoxical Relationship of a Creb 3 Regulatory Factor Missense Variant With Body Mass Index and Diabetes Among Samoans: Protocol for the Soifua Manuia (Good Health) Observational Cohort Study

BACKGROUND

The prevalence of obesity and diabetes in Samoa, like many other Pacific Island nations, has reached epidemic proportions. Although the etiology of these conditions can be largely attributed to the rapidly changing economic and nutritional environment, a recently identified genetic variant, rs373863828 (CREB 3 regulatory factor, CREBRF: c.1370G>A p.[R457Q]) is associated with increased odds of obesity, but paradoxically, decreased odds of diabetes.

OBJECTIVE

The overarching goal of the Soifua Manuia (Good Health) study was to precisely characterize the association of the CREBRF variant with metabolic (body composition and glucose homeostasis) and behavioral traits (dietary intake, physical activity, sleep, and weight control behaviors) that influence energy homeostasis in 500 adults.

METHODS

A cohort of adult Samoans who participated in a genome-wide association study of adiposity in Samoa in 2010 was followed up, based on the presence or absence of the CREBRF variant, between August 2017 and March 2019. Over a period of 7-10 days, each participant completed the main study protocol, which consisted of anthropometric measurements (weight, height, circumferences, and skinfolds), body composition assessment (bioelectrical impedance and dual-energy x-ray absorptiometry), point-of-care glycated hemoglobin measurement, a fasting blood draw and oral glucose tolerance test, urine collection, blood pressure measurement, hand grip strength measurement, objective physical activity and sleep apnea monitoring, and questionnaire measures (eg, health interview, cigarette and alcohol use, food frequency questionnaire, socioeconomic position, stress, social support, food and water insecurity, sleep, body image, and dietary preferences). In January 2019, a subsample of the study participants (n=118) completed a buttock fat biopsy procedure to collect subcutaneous adipose tissue samples.

RESULTS

Enrollment of 519 participants was completed in March 2019. Data analyses are ongoing, with results expected in 2020 and 2021.

CONCLUSIONS

While the genetic variant rs373863828, in CREBRF, has the largest known effect size of any identified common obesity gene, very little is currently understood about the mechanisms by which it confers increased odds of obesity but paradoxically lowered odds of type 2 diabetes. The results of this study will provide insights into how the gene functions on a whole-body level, which could provide novel targets to prevent or treat obesity, diabetes, and associated metabolic disorders. This study represents the human arm of a comprehensive and integrated approach involving humans as well as preclinical models that will provide novel insights into metabolic disease.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

RR1-10.2196/17329.

The Genetic Basis of Metabolic Disease

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