An adaptive variant of TRIB2, rs1057001, is associated with higher expression levels of thermogenic genes in human subcutaneous and visceral adipose tissues
Journal of Physiological Anthropology volume 36, Article number: 16 (2017)
An obesity-related single-nucleotide polymorphism (SNP) of the Tribbles pseudokinase 2 gene (TRIB2) was shown to have underwent adaptive evolution in the last glacial period, suggesting a selective advantage of this SNP in human populations in cold environments. In order to verify this hypothesis, the effect of the TRIB2 SNP on the expression of genes involved in adaptive thermogenesis was tested using messenger RNAs prepared from adipose tissues of Japanese adults.
Complementary DNA was prepared from subcutaneous adipose tissues (SAT) and visceral adipose tissues (VAT) obtained from 48 Japanese adults. Transcript levels of 15 selected genes, including five genes that are upregulated in development of thermogenic adipocytes, were measured by using real-time polymerase chain reaction. Differences in transcript levels between the TRIB2 SNP genotype groups (AA genotype versus AT + TT genotype) were assessed using t test.
Of the five thermogenic genes, DIO2, CIDEA, PPARGC1A, and PRDM16 showed significantly higher transcript levels in SAT of individuals with the AA genotype relative to those with the AT + TT genotype (P < 0.05). However, only 2 out of the 10 non-thermogenic genes exhibited differences in transcript levels according to genotype. Additionally, in silico prediction indicated that this SNP likely affects the expression of nearby genes including TRIB2.
The higher expression levels of thermogenic genes in individuals homozygous for the adaptive variant of TRIB2 SNP suggest a greater propensity for induction of thermogenesis in adipose tissues in cold environments.
The primary cause of obesity is an imbalance between caloric intake and energy expenditure. Susceptibility to obesity in present human populations is thought to be shaped by past genetic adaptation to famine [1, 2]. Additionally, lower ambient temperatures may exert selective pressure on genetic variations that influence susceptibility to obesity, as cold-induced thermogenesis substantially increases energy expenditure . It may be hypothesized that genotypes linked to higher thermogenic capacity and leaner phenotypes were adaptive in archaic human populations during glacial periods . A single-nucleotide polymorphism (SNP) in the uncoupling protein 1 gene, an essential gene for thermogenesis in brown adipocytes, has been shown to support the role of adaptation to cold climates in shaping the susceptibility to obesity [5, 6].
A previous study reported that a functional SNP in the 3′ untranslated region of the tribbles pseudokinase 2 gene (TRIB2) strongly influences visceral fat accumulation in Japanese adults. The obesity-resistance-associated allele of TRIB2 underwent positive natural selection in East Asians during the last glacial maximum, suggesting that this TRIB2 variant links past adaptation to present resistance to obesity . TRIB2 suppresses the differentiation of adipocytes by promoting the degradation of the CCAAT/enhancer binding protein beta (CEBPB), a transcription factor that acts during early stages of brown adipocyte development [8, 9]. Human brown adipocytes are thought to be present only during early life in humans; however, recent studies have indicated the existence of active brown adipocytes in adult humans . Moreover, it has been shown that energy-storing white adipocytes (or/and their progenitor cells) may be transformed into thermogenic beige/brite adipocytes via a series of stimuli; these thermogenic adipocytes resemble brown adipocytes but possess distinctive developmental lineages . The adaptive TRIB2 allele, which is associated with the reduction of visceral fat, may contribute to resistance to cold environments by increasing the activity of these thermogenic adipocytes. In the present study, we examined the effects of the previously reported TRIB2 SNP on the expression of genes involved in thermogenesis in the adipose tissues of Japanese adults.
Subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) samples were obtained from 48 Japanese adults. These participants were enrolled from among patients admitted to the Jichi Medical University Hospital for gastrointestinal surgery or gynecological surgery. All participants provided written informed consent. SAT was collected from the abdominal skin incision and VAT from the extirpated organ, omentum, or mesenterium during surgical procedures. Complementary DNA (cDNA) and genomic DNA were prepared from adipose tissue specimens. Details of participants and samples were described in the previous paper . The TRIB2 SNP (dbSNP: rs1057001) was genotyped using TaqMan Genotyping Assays and ABI PRISM 7900HT . Transcription levels of selected genes were measured using TaqMan Gene Expression Assays or RT-PCR with SYBR Green on an ABI PRISM 7900HT. The list of target genes is shown in Table 1. Transcript levels of the target genes were normalized to that of the actin beta gene (measured by the TaqMan method) using the delta Ct method. Statistical tests were performed using SPSS version 23. Epigenetic modification patterns near rs1057001 and tightly linked variants (r 2 > 0.8 in the 1000 Genome Project East Asians) were retrieved using Haploreg4.1 . The design of the present study was approved by the ethical committee of Jichi Medical University.
Results and discussion
The genotypes of rs1057001 in participants were AA = 26, AT = 20, and TT = 2. The proportion of the three genotypes was in Hardy-Weinberg equilibrium status (P = 0.421, chi square test, degrees of freedom = 1). Individuals with TT genotype were infrequent and thus combined with individuals with AT genotype for further statistical analyses. Genotype groups did not show differences in terms of the proportion of males and females (P = 0.066, Fisher’s exact test), age (P = 0.253, t test), and body mass index (P = 0.960, t test).
The transcript levels of tested genes are shown in Table 2. We tested CIDEA, DIO2, PPARGC1A, PRDM16, and UCP1, which were previously linked to brown adipose tissue activity as well as induction of beige/brite adipocytes in mice and humans [14–18]. The UCP1 transcript was undetectable in the majority of the adipose tissue samples. Transcripts of the other four genes were present at higher levels in individuals with the AA genotype than in those with the AT + TT genotype (P < 0.05, t test). The level of PPARGC1A expression in SAT was significant, at 5% after the Bonferroni correction (numbers of successfully measured gene = 14). The genotypic differences in the transcript levels of thermogenic genes were more notable in SAT, in which induction of beige/brite adipocytes was observed in mice and humans . Furthermore, CIDEA and PRDM16 showed higher expression levels in SAT than in VAT (P < 0.05, after the Bonferroni correction was applied), supporting the greater prepotency of SAT for thermogenesis. We additionally tested the effect of TRIB2 genotypes on the expression of genes for de novo lipogenesis in white adipocytes (ACC, DGAT1, FASN, SCD1, SREBP1) and genes commonly expressed in white and brown adipocytes (ADIPOQ, ADRB3, FABP4, GLUT4, and LEP). Of the 10 genes, only ADIPOQ and GLUT4 in SAT showed significant differences in expression between the genotype groups (P < 0.05, t test).
The previous study showed that the A allele of rs1057001 is linked to lower transcription levels of TRIB2 in adipose tissues . We additionally investigated epigenetic modifications of rs1057001 and short nucleotide sequence variants with tight linkage disequilibrium status with rs1057001 (Fig. 1). Several of these variants were found in genome sequences with signatures of gene regulatory elements, including histone enhancer marks, DNase hypersensitive sites, and alteration of protein-binding motifs. These variants are considered to coordinately alter expression levels of TRIB2 in various tissues. The transcription level of a microRNA-encoding gene (MIR3125) as well as of TRIB2 showed significant correlations with genotypes of rs1057001 or the tightly linked variants. Although the physiological function of MIR3125 is still unknown, this gene is thought to contribute to the genotype-related effects of the TRIB2 variant of the thermogenic genes.
Higher expression levels of genes involved in the differentiation and function of brown and beige/brite adipocytes may indicate a greater propensity for thermogenesis in individuals with the AA genotype in response to cold stimuli. This observation may support the hypothesis that the A allele contributed to cold adaptation in ancestors of East Asians in the last glacial period . The suppression of Trib2 has been shown to reduce the degradation of CEBPB, a key transcription factor involved in brown adipocyte differentiation in murine models [8, 9]. This suggests that the AA genotype, which is linked with suppressed expression of TRIB2, is associated with more potent differentiation of thermogenic adipocytes. Transcription levels of UCP1, an important component of the thermogenic machinery of brown and beige/brite adipocytes, were very low in the present adipose tissue samples. UCP1 expression is known to occur at very low levels in non-stimulated white adipose tissues in mice and humans [18, 19]. The present adipose tissue samples, which were obtained from patients who were not subject to cold simulation, did not exhibit high levels of UCP1 expression.
The present study has several potential limitations, which should be addressed in future work; (1) the observed gene expression pattern may be confounded by non-adipocyte cell populations, mainly stromal vascular cells, in the adipose tissues. (2) It is unclear whether the abdominal SAT mirrors the white adipose tissues in the supraclavicular region, which represents the main site of expression of thermogenic genes in human beige/brite adipocytes . (3) The molecular mechanisms underlying the modulation of thermogenic gene expression by TRIB2 remain unknown.
In conclusion, the present study provides evidence for the role of TRIB2 in energy expenditure in adipose tissues. These findings are believed to explain the previously proposed positive natural selection for the TRIB2 cold-adaptive variant during the last glacial period.
Adrenoceptor beta 3
- CEBPB CCAAT:
Enhancer binding proteinbeta
Cell death-inducing DFFA-like effector A
Diacylglycerol O-acyltransferase 1
Iodothyronine deiodinase 2
Fatty acid-binding protein 4
Fatty acid synthase
Facilitated glucose transporter, member 4
Polymerase chain reaction
PR/SET domain 16
PPARG coactivator 1 alpha
Subcutaneous adipose tissues
Stearoyl-CoA desaturase 1
Sterol regulatory element binding protein1
Tribbles pseudokinase 2
Uncoupling protein 1
Visceral adipose tissues
Neel JV. Diabetes mellitus: a “thrifty” genotype rendered detrimental by “progress”? Am J Hum Genet. 1962;14:353–62.
Prentice AM, Hennig BJ, Fulford AJ. Evolutionary origins of the obesity epidemic: natural selection of thrifty genes or genetic drift following predation release? Int J Obesity. 2008;32:1607–10.
Nishimura T, Motoi M, Egashira Y, Choi D, Aoyagi K, Watanuki K. Seasonal variation of non-shivering thermogenesis (NST) during mild cold exposure. J Physiol Anthropol. 2015;34:11.
Sellayah D, Cagampang FR, Cox RD. On the evolutionary origins of obesity: a new hypothesis. Endocrinology. 2014;155:1573–88.
Hancock AM, Clark VJ, Qian Y, Di Rienzo A. Population genetic analysis of the uncoupling proteins supports a role for UCP3 in human cold resistance. Mol Biol Evol. 2011;28:601–14.
Yoneshiro T, Ogawa T, Okamoto N, Matsushita M, Aita S, Kameya T, et al. Impact of UCP1 and beta3AR gene polymorphisms on age-related changes in brown adipose tissue and adiposity in humans. Int J Obesity. 2013;37:993–8.
Nakayama K, Ogawa A, Miyashita H, Tabara Y, Igase M, Kohara K, et al. Positive natural selection of TRIB2, a novel gene that influences visceral fat accumulation, in East Asia. Hum Genet. 2013;132:201–17.
Naiki T, Saijou E, Miyaoka Y, Sekine K, Miyajima A. TRB2, a mouse Tribbles ortholog, suppresses adipocyte differentiation by inhibiting AKT and C/EBPbeta. J Biol Chem. 2007;282:24075–82.
Kajimura S, Seale P, Kubota K, Lunsford E, Frangioni JV, Gygi SP, et al. Initiation of myoblast to brown fat switch by a PRDM16-C/EBP-beta transcriptional complex. Nature. 2009;460:1154–8.
Yoneshiro T, Saito M. Activation and recruitment of brown adipose tissue as anti-obesity regimens in humans. Ann Med. 2015;47:133–41.
Wang W, Seale P. Control of brown and beige fat development. Nat Rev Mol Cell Biol. 2016;17:691–702.
Munkhtulga L, Nagashima S, Nakayama K, Utsumi N, Yanagisawa Y, Gotoh T, et al. Regulatory SNP in the RBP4 gene modified the expression in adipocytes and associated with BMI. Obesity. 2010;18:1006–14.
Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012;40:D930–4.
Bostrom P, Wu J, Jedrychowski MP, Korde A, Ye L, Lo JC, et al. A PGC1-alpha-dependent myokine that drives brown-fat-like development of white fat and thermogenesis. Nature. 2012;481:463–8.
Fisher FM, Kleiner S, Douris N, Fox EC, Mepani RJ, Verdeguer F, et al. FGF21 regulates PGC-1alpha and browning of white adipose tissues in adaptive thermogenesis. Genes Dev. 2012;26:271–81.
Petrovic N, Walden TB, Shabalina IG, Timmons JA, Cannon B, Nedergaard J. Chronic peroxisome proliferator-activated receptor gamma (PPARgamma) activation of epididymally derived white adipocyte cultures reveals a population of thermogenically competent, UCP1-containing adipocytes molecularly distinct from classic brown adipocytes. J Biol Chem. 2010;285(10):7153–64.
Seale P, Kajimura S, Yang W, Chin S, Rohas LM, Uldry M, et al. Transcriptional control of brown fat determination by PRDM16. Cell Metab. 2007;6:38–54.
Seale P, Conroe HM, Estall J, Kajimura S, Frontini A, Ishibashi J, et al. Prdm16 determines the thermogenic program of subcutaneous white adipose tissue in mice. J Clin Invest. 2011;121:96–105.
Jespersen NZ, Larsen TJ, Peijs L, Daugaard S, Homoe P, Loft A, et al. A classical brown adipose tissue mRNA signature partly overlaps with brite in the supraclavicular region of adult humans. Cell Metab. 2013;17:798–805.
We are grateful to the participants of the present study.
This work was partly supported by the MEXT KAKENHI Grant Numbers 26291096 and 26251050.
Availability of data and materials
The datasets during and/or analyzed during the current study are available from the corresponding author on reasonable request.
KN conceived this study. SI provided the DNA and cDNA samples. KN performed the experiments. KN wrote the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The design of this study was approved by the Ethical Committee of the Jichi Medical University. All participants provided written informed consent.
About this article
Cite this article
Nakayama, K., Iwamoto, S. An adaptive variant of TRIB2, rs1057001, is associated with higher expression levels of thermogenic genes in human subcutaneous and visceral adipose tissues. J Physiol Anthropol 36, 16 (2017). https://doi.org/10.1186/s40101-017-0132-z