The impact of CREBRF rs373863828 Pacific-variant on infant body composition

In Māori and Pacific adults, the CREBRF rs373863828 minor (A) allele is associated with increased body mass index (BMI) but reduced incidence of type-2 and gestational diabetes mellitus. In this prospective cohort study of Māori and Pacific infants, nested within a nutritional intervention trial for pregnant women with obesity and without pregestational diabetes, we investigated whether the rs373863828 A allele is associated with differences in growth and body composition from birth to 12–18 months’ corrected age. Infants with and without the variant allele were compared using generalised linear models adjusted for potential confounding by gestation length, sex, ethnicity and parity, and in a secondary analysis, additionally adjusted for gestational diabetes. Carriage of the rs373863828 A allele was not associated with altered growth and body composition from birth to 6 months. At 12–18 months, infants with the rs373863828 A allele had lower whole-body fat mass [FM 1.4 (0.7) vs. 1.7 (0.7) kg, aMD −0.4, 95% CI −0.7, 0.0, P = 0.05; FM index 2.2 (1.1) vs. 2.6 (1.0) kg/m2 aMD −0.6, 95% CI −1.2,0.0, P = 0.04]. However, this association was not significant after adjustment for gestational diabetes, suggesting that it may be mediated, at least in part, by the beneficial effect of CREBRF rs373863828 A allele on maternal glycemic status.

lean mass at 2 to 4 months 12 and a taller stature at 5 to 18 years of age 13 .Consistent with its protective role in diabetes [8][9][10] , the rs373863828 A allele was associated with a five-fold reduction in the likelihood of GDM in New Zealand Māori and Pacific pregnant women with obesity, without any apparent influence on gestational weight gain or maternal lipid concentrations 14 .
It is unclear if fetal carriage of the rs373863828 A allele alters neonatal and infant growth and body composition independent of its effects on maternal glucose homeostasis.The presence of such an association could have important implications for the interpretation of metabolic studies in pregnancy in Māori and Pacific populations.Therefore, among Māori and Pacific infants exposed to maternal obesity and a high prevalence of GDM, we investigated the relationship between CREBRF rs373863828 genotype and body size, composition and growth from birth to 12 to 18 months of age.We hypothesised that Māori and Pacific infants carrying the rs373863828 A allele, compared to those without the allele, would have increased growth in weight and length at 12 to 18 months but a leaner body composition.

Study population
This study was undertaken among infants born to women in the Healthy Mums and Babies (HUMBA) Trial (ACTRN12615000400561, first registration 29/4/2015) conducted in South Auckland, New Zealand, where more than half of the maternity population is of Māori or Pacific descent 15 .Ethical approval was obtained from the Southern Health and Disability ethics committee, New Zealand (14/STH/205) and the study was performed in accordance with the Declaration of Helsinki.The HUMBA Trial had a 2 × 2 factorial design and investigated in pregnant women with obesity (BMI ≥ 30 kg/m 2 ) whether excessive gestational weight gain and birthweight in their infants could be reduced by: (1) a multi-faceted dietary intervention provided by Community Health Workers compared with routine dietary advice; and/or (2) probiotics compared with placebo 15 .Women with preexisting diabetes [HbA 1c ≥ 50 mmol/mol (≥ 6.7%)] in early pregnancy were excluded.The dietary intervention, which consisted of four home-based education sessions, reinforced with behaviour change techniques, personalised pregnancy weight gain targets and motivational text messaging, had a modest effect on reducing gestational weight gain but did not affect infant birthweight.Probiotics did not alter either of the co-primary outcomes.
Infants were eligible for this study if their mother or father self-identified as being of Māori or Pacific (Polynesian) descent, and their parent/guardian provided written informed consent for genetic testing.Genomic deoxyribonucleic acid (DNA) was collected from buccal samples at follow-up at either of two time points, 5 to 6 months or 12 to 18 months, using the Oragene saliva kit (DNA Genotek Inc., Ottawa, Canada), and was extracted using the PureLink Genomic DNA Mini Kit (Invitrogen, USA).DNA samples were stored at −80 °C prior to genetic testing.
Anthropometric and body composition measurements were performed in infants at birth and 5 to 6 months and 12 to 18 months' corrected age, as previously described 15 .Briefly, weight was measured to 10 g accuracy by electronic scale; length to 1 mm by infantometer; circumferences for head, left arm (mid-acromio-radiale), chest and abdomen to 1 mm by lasso tape; and triceps and subscapular skinfold thickness to 0.2 mm by Holtain calipers.Whole-body fat and lean mass (FM, LM) were measured at birth and 5 to 6 months by air displacement plethysmography (PEA POD, Cosmed, Concord, USA) and at 12 to 18 months by hand-to-foot bioimpedance spectroscopy analysis (SFB7, Impedimed, Queensland, Australia).

Statistical analysis
Analysis was performed using SAS v9.4 (SAS Institute, Cary, NC, USA).Socioeconomic deprivation was defined as being in the lowest national quintile for New Zealand Deprivation Index, calculated using maternal address at study recruitment 16 .Customised birthweight centiles were calculated using the New Zealand version of GROW 17 .Sex-and corrected-age-specific z-scores for anthropometric measures at birth were calculated from the United Kingdom World Health Organization (UK-WHO) preterm population reference and subsequently from the WHO 2006 Child Growth Standard 18 .Small for gestational age (SGA) and LGA were defined as < 10th or > 90th birthweight centile, respectively, for customised and population references.Left arm muscle area was calculated from arm circumference and triceps skinfold thickness 19 .Whole-body FM and LM were converted to height indices (FMI, LMI kg/m 2 ) to account for the body size.Conditional growth in weight and length from birth and 5 months to 12 to 18 months' corrected age was assessed as conditional growth SS, calculated from the residuals of the regression of the current z-scores on previous z-scores 20 .
In the primary analysis (model 1), the carriage of the rs373863828 A allele was tested for association with body size and composition at birth, 5 to 6 months' and 12 to 18 months' corrected age, and growth in weight and length from birth to 12 to 18 months' using generalised linear models, adjusted for potential confounding by gestational age, sex, ethnicity and parity.Estimated exposure effects are presented as adjusted mean difference (aMD) or adjusted odds ratio (aOR).We estimated that a study of 70 infants with an allelic ratio of AAA/AG to GG of 1:2 would have 80% power to detect a 0.7 standard deviation (SD) difference in continuous measures.
In secondary analysis (model 2), regression models were additionally adjusted for potential mediation by maternal GDM status using the diagnostic thresholds of the International Association of Diabetes in Pregnancy Study Groups (IADPSG) 21 .For this analysis, only infants whose mothers completed an oral glucose tolerance test were included, to ensure accurate diagnosis of GDM.In post hoc exploratory analysis, the primary model was further adjusted for the HUMBA trial interventions.

Results
Of 140 HUMBA infants followed up at 12 to 18 months corrected age, 67 were of Māori or Pacific descent and had consent provided for genetic testing.Eighteen (27%) were carriers for the CREBRF rs373863828 A allele, with two (11%) being homozygous (Table1).Approximately half of infants with the rs373863828 A allele had a mother with the variant allele, consistent with Mendelian inheritance.Infants with and without the rs373863828 A allele were similar for maternal age, parity, body mass index and smoking; socioeconomic deprivation; sex; ethnicity; mode of birth and gestation length (Table 1).Although not statistically significant, mothers of infants carrying the rs373863828 A allele were three times less likely to have GDM (Table 1), consistent with our previous report 14 .www.nature.com/scientificreports/In primary analysis, infant carriage of the rs373863828 A allele was not significantly associated with body size at birth, including weight and length z-scores, BMI or being LGA, nor measures of body composition, including skinfold thickness, arm muscle area or whole-body FM and LM, and associated indices (Table 2).Similarly, at 5 to 6 months' corrected age, infant carriage of the rs373863828 A allele was not significantly associated with body size or composition (Table 3).In secondary post hoc analysis, results were not appreciably altered by adjusting the primary model for exposure to the HUMBA Trial interventions.
In secondary analysis, this association was not significant when additionally adjusted for maternal GDM (Table 4).In secondary post hoc analysis, results at 12 to 18 months were not appreciably altered by adjusting the primary model for exposure to the HUMBA Trial interventions.
Conditional growth in weight and length was not significantly different between infants with and without the rs373863828 A allele from birth to 5 to 6 months or 12 to 18 months, and from 5 to 6 months to 12 to 18 months (Table 5).However, conditional growth in length from 5 to 6 months was positive in infants with the rs373863828 A allele but negative in those without (Table 5).

Discussion
We found that HUMBA infants of Polynesian ancestry carrying the CREBRF rs373863828 A allele had similar body size and composition at birth and in early infancy compared to those without the variant allele.However, by 12 to 18 months' corrected age, infants carrying the rs373863828 A allele had lower whole-body FM by 0.4 kg or ~ 0.6 SD, with non-significant reductions in measures of abdominal but not subcutaneous fat, suggesting a possible beneficial effect on central adiposity.Despite the decrease in FM, BMI at 12 to 18 months was similar between infants with and without the rs373863828 A allele, due to a similar absolute but non-significant increase in whole-body LM.  www.nature.com/scientificreports/

Physiological role of the CREBRF gene
The CREBRF gene is widely expressed in all tissues and encodes for a transcription factor involved in basic cell functions, such as energy metabolism and cell differentiation 9 .The minor rs373863828 A allele, possibly positively selected for in Polynesian and Micronesian populations, is known to be associated with higher BMI but a leaner body composition 11,22 , lower fasting glucose concentrations, and decreased incidence of type 2 diabetes and GDM 23 .The rs373863828 A allele is associated with increased glucose-stimulated insulin release, without an effect on insulin sensitivity, suggesting that the variant may preserve β-cell function 24 .Postprandial insulin secretion is a biphasic process, with a rapid first-phase of release followed by a reduced, sustained second phase, which continues until normoglycaemia is restored 25 .In adults, carriage of the rs373863828 A allele seems to enhance the first phase of insulin release, resulting in better glycaemic control 24 .In addition, the rs373863828 A allele may increase the capacity for insulin-independent glucose clearance in skeletal muscle, either by increasing muscle mass, the most abundant component of lean mass 26 , or by upregulating the glycolysis pathway in muscle 27 .

Impact of the CREBRF rs373863828 A allele in neonates
The impact of the CREBRF rs373863828 A allele on infant growth has been less well studied.Arslanian et al.  assessed the relationship between the CREBRF rs373863828 A allele and neonatal body composition (mean age 6 days) measured by dual-energy X-ray 12 .Similar to our study, they did not find any difference between neonates with and without the rs373863828 A allele in body size, whole body fat and lean mass, and skinfold thickness, even after accounting for age at assessment, sex, breastfeeding status, and maternal BMI.Neonates with the rs373863828 A allele had significantly higher bone mass, although the difference was small (+ 3 g).We additionally adjusted body composition models for maternal GDM, but this did not influence results at birth.Together, these two studies suggest that fetal carriage of the CREBRF rs373863828 A allele is not associated with altered size at birth, regardless of any effect of maternal rs373863828 status on glucose homeostasis in pregnancy, and that studies investigating the effect of pregnancy interventions on neonatal body size and composition in Polynesian populations are unlikely to be confounded by the CREBRF rs373863828 A allele.

Impact of the CREBRF rs373863828 A allele in infancy
However, by 2 to 4 months of age, Arslanian et al. found that infants carrying the rs373863828 A allele had a modest increase in whole-body lean mass (+ 200 to 250 g) but no evidence of altered body fat or length 12 .We did not find any significant difference in whole-body LM at 5 to 6 months' corrected age, although the absolute difference between groups was similar (+ 300 g), suggestive of a type II error.The study by Arslanian et al. may have had greater precision due to a larger sample size and use of DEXA.Consistent with these findings, we found that by 12 to 18 months, carriage of the rs373863828 A allele was associated with a non-significant increase in whole-body LM of 400 g, and results were similar after adjusting for exposure to maternal GDM.In contrast to the data for LM, the association between the rs373863828 A allele and whole-body FM at 12 to 18 months was reduced and no longer significant after adjustment for maternal GDM.This may be explained by the fact that GDM has been associated with increased FM in infancy 28 but in this cohort maternal carriage of the CREBRF rs373863828 A allele was associated with reduced incidence of GDM 14 .Future studies should evaluate further whether any potential effect of the rs373863828 A allele on FM gain in infancy is mediated, at least in part, by maternal CREBRF rs373863828 and GDM status, and whether this apparent leaner body composition tracks through childhood and beyond.

Association of between the CREBRF rs373863828 A allele and linear growth
The timing of a possible effect of the CREBRF rs373863828 A allele on linear growth is also unclear.Carlson et al.  found that in Samoan children and teenagers from age 5 to 18 years, the CREBRF rs373863828 A allele was associated with a 0.4 increase in height z-score per allele, consistent with findings in Samoan adults 13 .Although we did not find significant differences in length in infancy, conditional growth in length from birth to 12 to 18 months was positive infants who carried the rs373863828 A allele but negative in those who did not, suggesting that the variant may be associated accelerated linear growth from late infancy.It is unclear if body segment proportions are altered, although accelerated growth before puberty usually affects appendicular more than axial growth.Further research is needed to understand the biological mechanism by which CREBRF rs373863828 A allele may affect skeletal growth, although it is possible that this may be partly mediated by reduced exposure to GDM as GDM has previously been associated with decreased linear growth in early childhood 28 .

Strengths and limitations
Strengths of our study include prospective data collection from early pregnancy, standardised assessment for GDM, detailed anthropometric and body composition assessment, use of conditional growth analysis and assessment of maternal genotype.The main limitation of our study is the modest sample size and the potential for type 2 error.Other limitations are that there were some missing maternal data on GDM for the secondary analysis, and that we could not explore the relationship between infant body composition and maternal BMI as all women in the HUMBA Trial were classified as having obesity.

Conclusion
In this cohort of Māori and Pacific infants exposed antenatally to maternal obesity, carriage of the CREBRF minor rs373863828 A allele did not alter body size and composition at birth and early infancy but was associated with lower FM at 12 to 18 months' corrected age, although this association may have been mediated by reduced exposure to GDM.

Table 1 .
Characteristics of infants and their mothers.Data are number (percent) or mean (standard deviation) a Homozygous (A/A) n = 2. b New Zealand Deprivation Index (2013) decile 8 to 1019.A CREBRF rs373863828 minor allele, BMI body mass index, GDM gestational diabetes mellitus, IADPSG International Association of Diabetes in Pregnancy Study Groups, OGTT oral glucose tolerance test.P value is for the comparison between genotype groups (A/G or A/A vs. GG) using Fisher exact test or Student t-test.To convert glucose concentration from mmol/L to mg/dL divide by 0.0555.Vol:.(1234567890)Scientific Reports | (2024) 14:8825 | https://doi.org/10.1038/s41598-024-59417-5

Table 2 .
Body size and composition at birth.Data are number (percent) or mean (standard deviation).Z-scores specific for gestation and sex.Model 1: adjusted for potential confounding by gestational age, sex, ethnicity and parity.Model 2: adjusted for potential confounding by gestational age, sex, ethnicity, parity and mediation by gestational diabetes (International Association of Diabetes in Pregnancy Study Groups diagnostic thresholds).FM fat mass, FMI fat mass index, LM lean mass, LMI lean mass index.

Table 3 .
Body size and composition at 5 to 6 months.Data are number (percent) or mean (standard deviation).Z-scores specific for gestation and sex.Model 1: adjusted for potential confounding by gestational age, sex, ethnicity and parity.Model 2: adjusted for potential confounding by gestational age, sex, ethnicity, parity and mediation by gestational diabetes (International Association of Diabetes in Pregnancy Study Groups diagnostic thresholds).FM fat mass, FMI fat mass index, LM lean mass, FMI fat mass index.

Table 5 .
Growth in infancy.Data are mean (standard deviation).SS, standard score.Model 1: adjusted for potential confounding by gestational age, sex, ethnicity and parity.Model 2: adjusted for potential confounding by gestational age, sex, ethnicity, parity and mediation by gestational diabetes (International Association of Diabetes in Pregnancy Study Groups diagnostic thresholds).