CREBRF rs373863828 minimally regulates gene expression in subcutaneous adipose tissue of Samoan adults

A population-specific missense variant, rs373863828 (G > A, p.Arg457Gln), has been identified in Pacific Islanders and is associated with increased body mass index and lower odds of type 2 diabetes. However, the mechanism underlying these associations is unknown. The variant is located in exon 5 of CREBRF, a gene which is rapidly upregulated in response to nutrient deprivation. This study aimed to assess the effect of rs373863828 on gene expression in adipose tissue from Samoan adults from the Soifua Manuia (''Good Health'') study. Needle subcutaneous adipose (upper buttock) biopsies were obtained in Apia, Samoa, from fasting Samoans without diabetes. RNA isolated from adipose tissue (n=91, 65.9% female) was subjected to bulk paired-end mRNA-sequencing (2x150 base pair). Differential gene expression among the three rs373863828 genotypes (GG, AG, AA) was assessed using a likelihood-ratio test implemented in DESeq2, and overrepresentation analysis for significant differentially expressed genes (DEGs) (padj < 0.1) was performed with clusterProfiler and oRG.DB. Bulk expression was deconvoluted to determine the cell type composition of the adipose tissue. Twenty-two DEGs were associated with rs373863828 genotype including CBX7 (padj=0.026 ) and SELENOM (padj= 0.004), both of which have roles in adipogenesis and metabolism. A single expression profile comprising 13 of the 22 DEGs corresponding to CREBRF genotype was identified, and overrepresentation analysis supported this finding, with terms such as ''keratinocyte differentiation'' and ''skin development'' enriched among significant DEGs. However, PERMANOVA of adipose cell type proportions showed no difference (p > 0.05) among the three genotypes. rs373863828 genotype was associated with expression of 22 genes in whole subcutaneous adipose tissue of this sample of fasting Samoan adults without diabetes. Future work should assess cell- and tissue-specific rs373863828 genotype effects using a combination of single-cell approaches and human cell models with perturbations associated with nutritional stress.