Next Article in Journal
Effect of Selected Silyl Groups on the Anticancer Activity of 3,4-Dibromo-5-Hydroxy-Furan-2(5H)-One Derivatives
Next Article in Special Issue
Histamine H3 Receptor Ligands—KSK-59 and KSK-73—Reduce Body Weight Gain in a Rat Model of Excessive Eating
Previous Article in Journal
A Comprehensive Review of HLA and Severe Cutaneous Adverse Drug Reactions: Implication for Clinical Pharmacogenomics and Precision Medicine
Previous Article in Special Issue
Withasomniferol D, a New Anti-Adipogenic Withanolide from the Roots of Ashwagandha (Withania somnifera)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pharmaceuticals
Volume 14
Issue 11
10.3390/ph14111078
pharmaceuticals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

BMP7 Increases UCP1-Dependent and Independent Thermogenesis with a Unique Gene Expression Program in Human Neck Area Derived Adipocytes

by
Abhirup Shaw
1,2,
Beáta B. Tóth
1,
Rini Arianti
1,2,
István Csomós
3,
Szilárd Póliska
4,
Attila Vámos
1,2,
Zsolt Bacso
3,5,
Ferenc Győry
6,
László Fésüs
1,*,† and
Endre Kristóf
1,*,†
1
Laboratory of Cell Biochemistry, Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, H-4032 Debrecen, Hungary
2
Doctoral School of Molecular Cell and Immune Biology, University of Debrecen, H-4032 Debrecen, Hungary
3
Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, H-4032 Debrecen, Hungary
4
Genomic Medicine and Bioinformatics Core Facility, Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, H-4032 Debrecen, Hungary
5
Faculty of Pharmacy, University of Debrecen, H-4032 Debrecen, Hungary
6
Department of Surgery, Faculty of Medicine, University of Debrecen, H-4032 Debrecen, Hungary
*
Authors to whom correspondence should be addressed.
†
These authors have contributed equally to this work and share last authorship.
Pharmaceuticals 2021, 14(11), 1078; https://doi.org/10.3390/ph14111078
Submission received: 29 September 2021 / Revised: 19 October 2021 / Accepted: 21 October 2021 / Published: 25 October 2021
(This article belongs to the Special Issue Searching for New Therapeutic Targets with Anti-obesity Potential)
Browse Figures
Review Reports Versions Notes

Abstract

:
White adipocytes contribute to energy storage, accumulating lipid droplets, whereas brown and beige adipocytes mainly function in dissipating energy as heat primarily via the action of uncoupling protein 1 (UCP1). Bone morphogenic protein 7 (BMP7) was shown to drive brown adipocyte differentiation in murine interscapular adipose tissue. Here, we performed global RNA-sequencing and functional assays on adipocytes obtained from subcutaneous (SC) and deep-neck (DN) depots of human neck and differentiated with or without BMP7. We found that BMP7 did not influence differentiation but upregulated browning markers, including UCP1 mRNA and protein in SC and DN derived adipocytes. BMP7 also enhanced mitochondrial DNA content, levels of oxidative phosphorylation complex subunits, along with PGC1α and p-CREB upregulation, and fragmentation of mitochondria. Furthermore, both UCP1-dependent proton leak and UCP1-independent, creatine-driven substrate cycle coupled thermogenesis were augmented upon BMP7 addition. The gene expression analysis also shed light on the possible role of genes unrelated to thermogenesis thus far, including ACAN, CRYAB, and ID1, which were among the highest upregulated ones by BMP7 treatment in both types of adipocytes. Together, our study shows that BMP7 strongly upregulates thermogenesis in human neck area derived adipocytes, along with genes, which might have a supporting role in energy expenditure.

1. Introduction

Two functionally distinct types of fat are present in mammals: white and brown adipose tissue. Healthy human adults possess thermogenic brown adipose tissue (BAT) that can dissipate energy as heat under sub-thermal conditions [1,2]. These depots are primarily located in six anatomical regions: supraclavicular, axillary, mediastinal, cervical, paravertebral, and abdominal [3,4,5]. Supraclavicular, paravertebral, and deep-neck (DN) regions contain the highest amounts of thermogenic adipose tissue. Together, these depots account for almost 5% of the basal metabolic rate, thereby playing an important role in combating weight gain and type 2 diabetes [6,7]. The thermogenic adipocytes in rodents are either classical brown or beige cells distinguished by their origin and distribution [8,9,10,11]. Non-shivering thermogenesis of interscapular BAT and inguinal white adipose tissue (WAT) depots of rodents is majorly mediated by uncoupling protein 1 (UCP1), which is mitochondria resident and functions via generating proton leak in the inner mitochondrial membrane thereby uncoupling the mitochondrial respiratory chain from ATP synthesis [8,12,13,14]. A UCP1-independent mechanism for thermogenesis has been recently discovered in murine beige fat, which is mediated via the futile creatine phosphate cycle [15,16,17]. In this process, creatine kinase first phosphorylates creatine using ATP. Creatine phosphate is immediately dephosphorylated, most probably by tissue-nonspecific alkaline phosphatase (TNAP), generating heat [18]. Novel strategies for specific stimulation of beneficial fat burning by the metabolically active BAT in humans can aid weight reduction and decrease insulin resistance in individuals with obesity [6,7].
Bone morphogenic proteins (BMPs) act via heterotetrametric complexes of transmembrane receptors, which are type I or type II serine/threonine kinase receptors [19,20,21]. Type I BMP receptors (BMPR) are divided into BMPR1A, BMPR1B, Activin A Receptor (ACVR)-like (L) 1, and ACVR1. BMP7 has been known to interact with BMPR1A, BMPR1B, and ACVR1. It also interacts with three distinct type II receptors: BMPR2, ACVR2A, and ACVR2B [22,23,24]. BMP7 acts as an auto/paracrine mediator that drives classical brown adipocyte differentiation in mice; it promotes differentiation of brown preadipocytes even in the absence of an adipogenic induction cocktail [25]. BMP7 increases the expression of UCP1, peroxisome proliferator-activated receptor-gamma coactivator (PGC)1α, and PRDM16 via the activation of p38 mitogen activated protein (MAP) kinase pathway [25]. BMP7 has also been shown to upregulate UCP1 expression in selected clones of immortalized white and brown preadipocytes derived from human neck [26]. Further studies are necessary to better understand the molecular mechanisms by which BMP7 can exert its thermogenic effects on distinct types of human adipocytes.
Recently, we screened and compared global gene expression patterns by RNA-sequencing of human primary differentiated adipocytes derived from progenitors of DN and subcutaneous neck (SC) origins [27,28]. Consistent with previous studies [26,29,30,31,32], DN derived adipocytes displayed higher browning features than the SC derived ones. The purpose of the present study was to learn whether BMP7 can enhance thermogenic differentiation of human neck area derived adipocytes. Our results indicate that BMP7 caused a significant upregulation of UCP1 protein expression in both SC and DN adipocytes that was further confirmed by immunostaining. BMP7 treatment enhanced mitochondrial biogenesis, fragmentation, and UCP1-dependent oxygen consumption rate (OCR). We also detected an increase in creatine driven substrate cycle coupled OCR and induction of genes related to this futile cycle, which generates heat without UCP1 [15,18,33]. The global gene expression data revealed pathways thus far unrelated to thermogenesis but upregulated as a result of BMP7 administration.

2. Results

2.1. SC and DN Derived Preadipocytes Are Differentiated to Adipocytes to the Same Extent in the Absence or Presence of BMP7

Adipose tissue biopsies of SC and DN origins were obtained in pairs from nine independent donors, followed by isolation and cultivation of preadipocytes according to previously described protocols [27]. White adipogenic differentiation medium was applied to the preadipocytes for 14 days, with or without BMP7. A heatmap prepared from data of global RNA-sequencing for illustrating the expression of general adipocyte marker genes (e.g., LEP, FABP4, ADIPOQ, PPARG) showed similarly elevated expression in BMP7 treated and untreated differentiating adipocytes (Figure 1A). Laser-scanning cytometry based quantification of adipogenic differentiation rate showed that approximately 60% of the cells were differentiated upon administration of the white adipogenic differentiation medium; the presence of BMP7 did not influence the differentiation rate (Figure 1B) [34,35]. Gene expression of BMPR subunits, such as BMPR1A, BMPR1B, ACVR1, BMPR2, ACVR2A, and ACVR2B were abundant at the preadipocyte stage in both SC and DN derived adipocytes, their expression was not affected by BMP7 significantly (Figure 1C–H).

2.2. BMP7 Upregulated Browning Marker Genes Including UCP1 in SC and DN Derived Differentiated Adipocytes

A recent in-depth analysis of white, brown, and beige adipocyte transcriptomes of murine and human origins has been utilized to develop an algorithm named, BATLAS that can effectively calculate brown/beige adipocyte content in cell and biopsy samples [36]. BATLAS analysis of gene expression data from the nine independent donor samples revealed that brown/beige content of DN samples was significantly higher than of SC ones, which was further elevated upon BMP7 treatment (Figure 2A). The presence of BMP7 also increased the browning capacity of SC derived adipocytes. Next, texture sum variance was quantified using laser-scanning cytometry [34,35], which revealed that BMP7 treatment slightly decreased the size of lipid droplets in both SC and DN derived adipocytes, that suggests increased capacity for lipolysis and energy expenditure (Figure 2B). Gene expression of UCP1 tended to elevate upon BMP7 treatment in both types of adipocytes as determined by RNA-sequencing (Figure 2C). This upregulation was further confirmed by RT-qPCR analysis and was found to be statistically significant in the case of DN adipocytes (Figure 2D). Consistent with previous results, DN derived adipocytes expressed more UCP1 protein than SC derived ones, which was significantly upregulated in response to BMP7 in adipocytes of both origins (Figure 2E). Additionally, immunostaining clearly revealed an increase of the UCP1 intensity upon BMP7 treatment in adipocytes differentiated from both types of progenitors (Figure 2F; see secondary antibody control in Figure S1A, Supplementary Materials). Laser-scanning cytometry assisted quantification of UCP1 immunostaining intensity showed a significant increase in DN adipocytes as a result of BMP7 (Figure 2G) [34,35]. Together these data indicate that BMP7 upregulates browning marker genes, including UCP1, in primary human neck area derived adipocytes.

2.3. BMP7 Facilitated Mitochondrial Biogenesis Leading to Increased Protein Expression of Mitochondrial OXPHOS Complex Subunits

Next, we addressed the question how BMP7 treatment influences mitochondrial biogenesis and content. The protein expression of PGC1α, one of the principal regulators of mitochondrial biogenesis [37,38], was found to be significantly upregulated upon BMP7 treatment in both SC and DN derived adipocytes (Figure 3A). PGC1α can also be induced by the transcription factor CREB, which can bind to a functional CRE in the PGC-1 promoter [37]. A significant induction of CREB phosphorylation in DN adipocytes further suggested an increase in mitochondrial biogenesis upon BMP7 treatment via upregulation of CREB pathway (Figure 3B). The amount of mitochondrial respiratory chain complex subunits was also elevated upon BMP7 treatment. Complex I, II, and III subunits showed significant upregulation upon BMP7 treatment in both types of adipocytes, while Complex IV was significantly upregulated in case of DN adipocytes. Complex V showed an increasing trend, although statistically not significant, in both SC and DN derived adipocytes (Figure 3C). Collectively, these data clearly demonstrate that BMP7 upregulates mitochondrial biogenesis in SC and DN derived adipocytes.

2.4. BMP7 Treatment Increased Mitochondrial Fragmentation and Elevated Cellular Respiration with Enhanced Proton Leak

UCP1-enriched mitochondria are mostly fragmented [39]. Immunostaining of translocase of outer mitochondrial membrane 20 (TOM20) was performed to visualize mitochondrial morphology. The content of fragmented mitochondria was significantly higher in differentiated DN derived adipocytes compared to SC derived ones, which corresponds to their higher UCP1 content and thermogenic capacity. BMP7 significantly increased the fragmented mitochondria content in both SC and DN derived adipocytes (Figure 4A,B; see Figure S1B, Supplementary Materials, for secondary antibody control). Total mitochondrial DNA content was also increased significantly upon BMP7 treatment in both SC and DN derived adipocytes (Figure 4C). These data suggest that mitochondria of adipocytes differentiated in the presence of BMP7 possess an enhanced thermogenic potential.
Next, we intended to investigate the functional consequences of the increased mitochondrial content and fragmentation. Monitoring of OCR revealed that both basal and 3′,5′-cyclic adenosine monophosphate (cAMP) stimulated OCR was significantly elevated upon BMP7 treatment in SC and DN derived adipocytes. Importantly, proton leak related OCR, revealed by oligomycin addition, was significantly elevated, which indicates increased thermogenesis mediated by UCP1 (Figure 5A). Basal extracellular acidification rate (ECAR) showed an increasing trend upon BMP7 treatment, while cAMP stimulated ECAR was significantly elevated in DN derived adipocytes differentiated in the presence of BMP7 (Figure 5B).

2.5. Creatine Driven Substrate Cycle Related Thermogenesis Was Facilitated upon BMP7 Treatment

Recent studies in mice have shown that a futile creatine phosphate cycle plays an important role in thermogenic metabolism of beige adipocytes [15]. Our research group demonstrated the contribution of this cycle to thermogenesis in ex vivo models of human subcutaneous and neck area adipocytes [40,41,42]. We could observe significant increase of creatine driven substrate cycle related OCR, revealed by applying the creatine transport inhibitor, β-guanidinopropionic acid (β-GPA) [11], in both SC and DN derived adipocytes differentiated in the presence of BMP7 (Figure 6A). RNA expression of mitochondrial creatine kinase 2 (CKMT2), one of the possible kinases acting in the futile cycle, was found slightly elevated in the generated RNA-sequencing dataset, and this was confirmed by RT-qPCR analysis in SC and DN adipocyte samples upon BMP7 treatment. Importantly, the CKMT2 protein level was also significantly higher in these adipocytes (Figure 6B). A recent publication has shown that cytosolic creatine kinase B (CKB) is targeted to mitochondria and is indispensable for thermogenesis by the creatine phosphate cycle [33]. Evaluation of RNA-sequencing data revealed an increasing trend for CKB expression upon BMP7 treatment. BMP7 caused significant upregulation of CKB protein expression in DN derived adipocytes (Figure 6C). The phosphatase TNAP has been identified to hydrolyze phosphocreatine to creatine, thereby driving the futile creatine phosphate cycle in mitochondria of thermogenic adipocytes [18]. TNAP was found to be expressed at a high extent in both SC and DN derived adipocytes, however, its gene expression remained unchanged upon BMP7 treatment (Figure 6D). Together, these data suggest that BMP7 upregulates thermogenesis not only via UCP1 but also by the futile creatine phosphate cycle in human neck area adipocytes.

2.6. BMP7 Upregulates Genes, Including Aggrecan (ACAN), Crystallin Alpha B (CRYAB), and Inhibitor of Differentiation 1 (ID1), Thus Far Not Linked to Thermogenesis

Next, we aimed to explore whether BMP7 treatment could upregulate pathways in neck area derived adipocytes, which have not been related to thermogenesis thus far. Further analysis of RNA-sequencing data revealed that 121 and 60 genes were upregulated in SC and DN derived adipocytes in response to BMP7 treatment, while 190 (such as COL6A6, MMP27, MGAT4C, CCL11, and LINC01028) and 87 (such as KRT1, MYH8, CCL11, MATN4, and FAM180B) genes were downregulated, respectively. Genes, up- or downregulated by BMP7, were visualized by Volcano Plot (Figure 7A) and clustered by a heatmap (Figure S2, Supplementary Materials). The lists of all the genes upregulated (Table S1, Supplementary Materials) and downregulated (Table S2, Supplementary Materials) are provided in descending order of their fold change. 38 genes were commonly upregulated between SC and DN derived adipocytes, while 45 genes were commonly downregulated, respectively (Figure 7B). Gephi illustration of Panther pathway analysis illustrates that signaling by BMP involving SMAD group of transcription factors (SMAD6, SMAD7, SMAD9) was commonly upregulated in both SC and DN adipocytes (Table 1, Figure S2, Supplementary Materials), which is consistent with previously published results [43]. Integrin cell surface interactions, involving ITGA9, COMP, or ITGA8 and GPCR ligand binding, involving ADRA2A, ADRA2C, FZD1, FZD5, or ACKR1 were among the significantly elevated pathways only in SC derived adipocytes upon BMP7 treatment (Table 1, Figure S3, Supplementary Materials). Extracellular matrix proteins, involving ACAN were also found to be upregulated only in SC derived adipocytes (Table 1).
According to the RNA-sequencing analysis, ACAN, CRYAB, and ID1 were strongly upregulated upon BMP7 treatment in both types of adipocytes (Figure 7A). Gene expression of ACAN was confirmed to be significantly higher in SC adipocytes upon BMP7 treatment by RT-qPCR analysis, protein expression of ACAN also followed a similar pattern (Figure 7C). Gene expression of CRYAB, a member of small heat shock protein (HSP) 20 family, was found to be elevated upon BMP7 treatment by RNA-sequencing and the upregulation was further confirmed by RT-qPCR in both SC and DN adipocytes. CRYAB protein expression was similarly elevated (Figure 7D). ID1 gene expression was also significantly higher in SC and DN derived adipocytes differentiated in the presence of BMP7 as compared to untreated ones quantified by both RNA-sequencing and RT-qPCR analysis. ID1 protein expression was significantly elevated in SC derived adipocytes, while in DN adipocytes it showed a trend of upregulation (Figure 7E).

3. Discussion

BMP7 functions as an autocrine/paracrine mediator that promotes classical brown and beige adipocyte differentiation in mice [25,44]. Pre-treatment with BMP7 of immortalized human neck derived white and brown preadipocyte clones, followed by the administration of an adipogenic induction medium, also significantly increased UCP1 gene expression [26]. Our research group observed that BMP7 can significantly upregulate gene expression of UCP1 and the classical brown adipocyte marker, ZIC1 in differentiated Simpson–Golabi–Behmel syndrome (SGBS) adipocytes [42]. BMP7 has also shown to improve insulin signal transduction in cultured human hepatocytes [45]. BMP7 exerts its effect via upregulation of p38 MAP kinase and PGC1α, which leads to upregulation of mitochondrial biogenesis and UCP1 expression [25]. BMP7 can also act via phosphorylation of SMAD 1/5/8 [43].
In our study, the receptors involved in BMP signaling [19,20,22,23] were found to be abundantly expressed in the preadipocyte stage irrespective of the anatomical origin of the progenitors isolated from human neck area (Figure 1). BMP7 significantly upregulated UCP1-dependent thermogenesis when applied on top of a white adipogenic differentiation medium in human primary SC and DN derived adipocytes (Figure 2). BMP7 treatment led to increased protein expression of mitochondrial oxidative phosphorylation complex subunits (Figure 3C) and total mitochondrial DNA content (Figure 4C) in parallel with the upregulation of PGC1α and p-CREB pathways (Figure 3A,B), which is consistent with previously published studies [25]. Our study first showed that BMP7 elevated fragmented mitochondria content in both SC and DN derived adipocytes (Figure 4A,B), which indicates an increased thermogenic potential [39]. This was further proved by a significant increase in basal, stimulated, and proton leak OCR in both types of adipocytes differentiated in the presence of BMP7 (Figure 5).
Kazak et al. discovered that a creatine-mediated futile cycle is responsible for enhanced mitochondrial respiration in beige fat [15]. A recent publication from the same group showed that an adipose tissue specific knockout of glycine amidinotransferase, the rate limiting enzyme for creatine biosynthesis, made mice prone to diet-induced obesity [46]. This clearly illustrated the importance of creatine metabolism in energy expenditure in vivo. Another study in mice showed that creatine driven thermogenesis plays an important role in both UCP1 positive and negative beige adipocytes [47]. Mitochondrial creatine kinases, CKMT1 A/B, CKMT2, and most recently CKB have been postulated to play an important role in the UCP1-independent thermogenesis mediated by the futile cycle [15,33]. Previously, our research group has shown that BMP7 treatment in SGBS adipocytes significantly elevated creatine cycle related OCR [42]. We also demonstrated that clozapine-induced human browning adipocytes of abdominal origin increased their energy expenditure via the futile creatine cycle [41]. Our previous study showed that mitochondrial creatine kinases, CKMT1 A/B were expressed at a greater extent in DN compared to SC adipocytes; creatine cycle related OCR followed a similar pattern as well [27]. We report for the first time that BMP7 treatment significantly increases creatine phosphate cycle related thermogenesis, which was evident from a strong increase in creatine cycle related OCR (Figure 6A). We also observed an increased gene and protein expression of mitochondrial CKMT2 (Figure 6B) and CKB (Figure 6C), which further supports this conclusion.
It is well characterized that long-term administration of rosiglitazone induces browning of human adipocytes ex vivo [27,34,48] and increases beiging of WAT in mice [49]. In SC adipocytes, only seven genes (AF131216.3, TMEM132C, SLC7A10, STOX1, RIMS4, PLPR4, and CMYA5) were commonly upregulated and ten (TMEM204, TMEM158, EPHB6, PLA2G5, RSPO1, DPT, TMC2, GPM6B, TRPV6, and HPSE) were downregulated, respectively, by both BMP7 and rosiglitazone [27]. Out of the commonly upregulated ones, SLC7A10 encodes the alanine-serine-cysteine transporter-1 (ASC-1), which also showed higher expression in DN derived adipocytes as compared to SC derived ones [40]. Expression of SLC7A10 in human WAT inversely correlates with visceral obesity and insulin resistance [50]. ASC-1 promoted mitochondrial respiration of adipocytes, especially when they were stimulated for thermogenesis by a cell permeable cAMP analogue [40,50]. In DN adipocytes, only RSPO3 was commonly upregulated and nine genes (WISP1, CLSTN2, LTBP2, DPT, IFITM1, INMT, TMEM158, TNC, and RSPO1) were downregulated, respectively, by both treatments [27]. Our presented results have also revealed certain genes, such as ACAN, CRYAB, or ID1, unrelated to human adipocyte browning thus far, which show a strong correlation with increased thermogenesis upon BMP7 treatment suggesting that they might have a positive effect on thermogenesis. Further research is necessary to unravel their mechanism of action and molecular targets with respect to thermogenesis.
Aggrecan protein, encoded by ACAN, is a chondroitin sulphated proteoglycan, which functions as a critical structural component of cartilage [51]. Aggrecan is also found in the brain almost exclusively in the perineuronal net and is postulated to play an important role in its formation and function [52,53,54,55] and a similar task might be also needed in the innervation of browning adipose tissue. Gene and protein expression of ACAN was significantly upregulated in BMP7 treated SC derived adipocytes (Figure 7C).
CRYAB is one of the major structural proteins of eye lens and is found to be expressed in several other tissues. It belongs to the small HSP family and behaves as a chaperone, thereby protecting against apoptosis and oxidative stress [56,57,58,59,60,61]. In the lens, CRYAB is a substrate for transglutaminase 2 (TGM2)-mediated cross-linking [62,63]. Our group found that gonadal WAT of TGM2−/− mice exhibited a lower expression of beige marker genes, such as UCP1, TBX1, and TNFRFS9; therefore, these mice were cold intolerant [64]. CRYAB has also been identified as a novel adipokine and its protein expression is strongly induced during adipogenesis, reaching a 10-times higher level in mature human white adipocytes than in preadipocytes [65,66]. Gene expression of CRYAB was higher in visceral WAT of children with obesity as compared to normal-weight [67]. CRYAB protein expression was found to be higher in subcutaneous WAT of old as compared to young individuals with obesity [68]. CRYAB mRNA had high expression in thermogenically more potent DN as compared to SC adipocytes and was less expressed in adipocytes that carried the FTO obesity-risk genotype [27]. In addition, gene and protein expression of CRYAB increased significantly upon BMP7 treatment in SC and DN derived adipocytes (Figure 7D) in parallel to increased thermogenic potential pointing to a thus far not recognized significance of CRYAB in this process.
The ID proteins (ID 1–4) are a subfamily of helix–loop–helix (HLH) transcription factors lacking a basic DNA binding domain. ID proteins act predominantly via dimerization with other transcriptional regulators, mostly by basic-HLH (bHLH) factors, which fail to bind to DNA and thereby function as dominant negative regulators of bHLH proteins [69]. A recent study that screened protein expression patterns of different metabolically active tissues showed that ID1 is highly expressed in both WAT and BAT of adult mice, with BAT showing the highest expression [70]. However, the same study found that an adipose tissue specific overexpression of ID1 resulted in high fat diet-induced obesity in mice. Depletion of ID1 increased UCP1 gene expression in mouse WAT upon cold exposure [71]. Another study indicated that ID1 deficiency improves glucose tolerance and lowers serum insulin levels in mouse WAT [72,73]. Our study here shows that BMP7 significantly upregulated ID1 gene and protein expression in SC and DN derived adipocytes (Figure 7E). The increased ID1 expression correlated with increased UCP1 expression, mitochondrial biogenesis, and content. Functional experiments are required in the future to explore the direct effect of ID1 on human adipocyte thermogenesis.

4. Materials and Methods

4.1. Materials

All chemicals were obtained from Sigma Aldrich (Munich, Germany) unless otherwise stated.

4.2. Isolation, Maintenance, Differentiation, and Treatment of Human Neck Derived Human Adipose-Derived Stromal Cells (hASCs)

hASCs were isolated from the stromal-vascular fractions of SC and DN tissues of patients (35–75 years) undergoing a planned elective surgery. Biopsies were obtained in pairs from SC and DN regions to avoid inter-individual variations. Patients with known diabetes, abnormal thyroid hormone levels or malignant tumors were not included in this study. Written informed consent was obtained from all participants prior surgery [27,28,74,75].
hASCs were isolated, maintained, and checked for mycoplasma contamination as previously described [27,74]. Cells were grown to 90% confluency and differentiated following a reported white adipogenic differentiation protocol, with or without BMP7 (R&D Systems, Minneapolis, MN, USA, 354-BP) at 50 ng/mL [27,34]. Media was replaced at an interval of 4 days, followed by collection of cells on day 14.

4.3. RNA Isolation, RT-qPCR, and RNA-Sequencing

Cells were collected in Trizol reagent (Thermo Fisher Scientific, Waltham, MA, USA), followed by manual RNA isolation by chloroform extraction and isopropanol precipitation. Global transcriptome data were obtained via high throughput mRNA-sequencing [27]. Pathways were obtained using Panther Reactome Pathways (http://www.pantherdb.org/, access date 20 September 2021). Heatmaps were constructed on the Morpheus web tool (https://software.broadinstitute.org/morpheus, access date 20 September 2021) with Pearson correlation of rows and complete linkage based on calculated z-score of DESeq normalized data after log2 transformation [27]. Gephi images were constructed as previously described [27]. Browning content was calculated using the BATLAS Webtool (http://green-l-12.ethz.ch:3838//BATLAS, access date 20 September 2021) [36].
RNA quality was evaluated by Nanodrop (Thermo Fisher Scientific), followed by generation of cDNA by TaqMan reverse transcription reagent kit (Thermo Fisher Scientific) and qPCR analysis [40]. Normalized gene expression was quantified using the probes listed in Table S3 (Supplementary Materials).

4.4. Western Blot, Immunoblotting and Antibodies

Samples were separated by SDS-PAGE, transferred to PVDF membrane, and blocked by 5% skimmed milk solution as previously described [40,42,76]. The following primary antibodies were used overnight in 1% skimmed milk solution: anti-UCP1 (1:750, R&D Systems, MAB6158), anti-pCREB (1:1000, Merck-Millipore, Burlington, MA, USA, 05-667), anti-CREB (1:1000, Abcam, Cambridge, UK, ab31387), anti-PGC1α (1:1000, Santa Cruz Biotechnology, Dallas, TX, USA, H-300), anti-CKMT2 (1:1000, Novus Biologicals, Centennial, CO, USA, NBP2-13841), anti-CKB (1:1000, Novus Biologicals, NBP1-84460), anti-CRYAB (1:1000, Novus Biologicals, NB100-2519), anti-β actin (1:5000, Novus Biologicals, A2066), anti-OXPHOS (1:1000, Abcam, ab110411), anti-Aggrecan (1:1000, Novus Biologicals, NB100-74350), and anti-ID1 (1:1000, Novus Biologicals, JM92-13). HRP-conjugated goat anti-rabbit (1:10,000, Advansta, San Jose, CA, USA, R-05072-500) or anti-mouse (1:5000, Advansta, R-05071-500) IgG were used as secondary antibodies, respectively. Immunoreactive proteins were visualized by Immobilion Western chemiluminescence substrate (Merck-Millipore). Densitometry was carried out by FIJI.

4.5. Immunostaining and Image Analysis

SC and DN derived hASCs were plated and differentiated in 8-well Ibidi μ-chambers (Ibidi GmbH, Gräfelfing, Germany). Cells were washed with PBS prior to fixation by 4% paraformaldehyde. Cells were permeabilized with 0.1% saponin, followed by blocking with 5% skimmed milk as per previously described protocols [76]. The cells were incubated overnight with the following primary antibodies: anti-TOM20 (1:75, WH0009804M1) or anti-UCP1 (1:200, U6382) at room temperature. Incubation with the following secondary antibodies was kept for 3 h at room temperature: Alexa 647 goat anti-mouse IgG (1:1000, Thermo Fischer Scientific, A21236) or Alexa 488 goat anti-rabbit IgG (1:1000, Thermo Fischer Scientific, A11034). Nuclei were labeled with Propidium Iodide (1.5 μg/mL, 1 h). Images were obtained by Olympus FluoView 1000 confocal microscope. Fragmented mitochondria were defined by the size (pixel2) 0–100 AU and counted using FIJI(ImageJ) [76]. Adipogenic differentiation rate, texture sum variance, and UCP1 immunostaining intensity were quantified using laser-scanning cytometry as previously described [34,35].

4.6. Determination of Cellular OCR and ECAR

OCR and ECAR were measured by XF96 oximeter (Seahorse Biosciences, North Billerica, MA, USA). Cells were seeded and differentiated on XF96 assay plates with or without BMP7. During measurement, baseline respiration and acidification were measured for 30 min followed by stimulated OCR and ECAR. A single bolus dose of dibutyryl-cAMP (at 500 μM final concentration) was added to mimic adrenergic stimulation leading to the stimulated OCR and ECAR, which were measured at 30-min intervals. 3 h post treatment, oligomycin at 2 μM final concentration was added to block ATP synthase activity to measure proton leak OCR. For measuring creatine cycle related OCR, β-GPA was added at 2 mM concentration after cAMP treatment. Antimycin A at 10 μM final concentration was added at the end for baseline correction. The OCR was normalized to protein content [40].

4.7. Statistics and Image Preparation

Results are expressed as mean ± SD for the number of independent repetitions indicated. One-way ANOVA with Tukey’s post hoc test were used for multiple comparisons of groups. Two-tailed paired t test was used to compare between two individual groups. Graphpad Prism 9 was used to visualize the graphs and evaluation of statistics.

5. Conclusions

As a conclusion, to the best of our knowledge, we showed for the first time that BMP7 can significantly facilitate thermogenesis in human primary SC and DN derived ex vivo differentiating adipocytes. The increased thermogenesis was mediated primarily via UCP1 dependent proton leak and, to some extent, by the futile creatine phosphate substrate cycle. We demonstrated that both fragmented mitochondria content and mitochondrial OXPHOS complex subunits were elevated upon BMP7 treatment. The global RNA-sequencing analysis gave clue to novel genes, such as ID1, ACAN, and CRYAB, which might play a supportive role in thermogenesis. Further functional experiments are necessary to validate and confirm their possible action with this regard.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ph14111078/s1, Figure S1: Secondary antibody control images illustrating the specificity of the antibodies used for UCP1 (A) and TOM20 (B) immunostaining, Figure S2: Heatmap illustrating the up- and downregulated genes upon BMP7 treatment in Subcutaneous (SC) and Deep-neck (DN) derived differentiated adipocytes, Table S1: List of genes upregulated by BMP7 treatment in SC and DN derived differentiated adipocytes arranged in descending order of their log2fold change, Table S2: List of genes downregulated by BMP7 treatment in SC and DN derived differentiated adipocytes arranged in descending order of their log2 fold change, Figure S3: Gephi diagrams illustrating pathways and gene interaction networks upregulated by BMP7 treatment in SC and DN derived differentiated adipocytes, Table S3: Table listing gene expression assays used in the study.

Author Contributions

Conceptualization, A.S., B.B.T., L.F. and E.K.; methodology, A.S., I.C., S.P., Z.B. and E.K.; software, B.B.T., R.A. and Z.B.; validation, A.S., B.B.T. and E.K.; formal analysis, A.S., B.B.T., R.A., S.P., A.V., Z.B. and E.K.; investigation, A.S., R.A., S.P., A.V. and E.K.; resources, F.G.; data curation, A.S., B.B.T. and Z.B.; writing—original draft preparation, A.S. and E.K.; writing—review and editing, A.S., B.B.T., L.F. and E.K.; visualization, A.S., B.B.T. and R.A.; supervision, L.F. and E.K.; project administration, A.S. and E.K.; funding acquisition, L.F. and E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union and the European Regional Development Fund (GINOP-2.3.2-15-2016-00006) and the National Research, Development and Innovation Office (NKFIH-FK131424 and K129139) of Hungary. EK was supported by the János Bolyai Fellowship of the Hungarian Academy of Sciences.

Institutional Review Board Statement

The study was conducted according to guidelines of the Declaration of Helsinki, and approved by the Medical Research Council of Hungary (20571-2/2017/EKU; 28.04.2017). Experiments were performed strictly in accordance with the approved ethical regulations and guidelines.

Informed Consent Statement

Written informed consent was obtained from all participants prior the collection of tissue samples.

Data Availability Statement

RNA-sequencing data were deposited to Sequence Read Archive (SRA) database [https://www.ncbi.nlm.nih.gov/sra] (access date 30 September 2021) under accession number PRJNA607438. Other data that support the findings of this study are available from the corresponding authors ([email protected], [email protected]) upon reasonable request.

Acknowledgments

We thank Jennifer Nagy for technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cypess, A.M.; Lehman, S.; Williams, G.; Tal, I.; Rodman, D.; Goldfine, A.B.; Kuo, F.C.; Palmer, E.L.; Tseng, Y.H.; Doria, A.; et al. Identification and importance of brown adipose tissue in adult humans. N. Engl. J. Med. 2009, 360, 1509–1517. [Google Scholar] [CrossRef] [Green Version]
  2. Leitner, B.P.; Huang, S.; Brychta, R.J.; Duckworth, C.J.; Baskin, A.S.; McGehee, S.; Tal, I.; Dieckmann, W.; Gupta, G.; Kolodny, G.M.; et al. Mapping of human brown adipose tissue in lean and obese young men. Proc. Natl. Acad. Sci. USA 2017, 114, 8649–8654. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Van Marken Lichtenbelt, W.D.; Vanhommerig, J.W.; Smulders, N.M.; Drossaerts, J.M.; Kemerink, G.J.; Bouvy, N.D.; Schrauwen, P.; Teule, G.J. Cold-activated brown adipose tissue in healthy men. N. Engl. J. Med. 2009, 360, 1500–1508. [Google Scholar] [CrossRef] [Green Version]
  4. Virtanen, K.A.; Lidell, M.E.; Orava, J.; Heglind, M.; Westergren, R.; Niemi, T.; Taittonen, M.; Laine, J.; Savisto, N.J.; Enerbäck, S.; et al. Functional brown adipose tissue in healthy adults. N. Engl. J. Med. 2009, 360, 1518–1525. [Google Scholar] [CrossRef] [PubMed]
  5. Saito, M.; Okamatsu-Ogura, Y.; Matsushita, M.; Watanabe, K.; Yoneshiro, T.; Nio-Kobayashi, J.; Iwanaga, T.; Miyagawa, M.; Kameya, T.; Nakada, K.; et al. High incidence of metabolically active brown adipose tissue in healthy adult humans: Effects of cold exposure and adiposity. Diabetes 2009, 58, 1526–1531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Van Marken Lichtenbelt, W.D.; Schrauwen, P. Implications of nonshivering thermogenesis for energy balance regulation in humans. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2011, 301, R285–R296. [Google Scholar] [CrossRef] [Green Version]
  7. Jimenez-Munoz, C.M.; López, M.; Albericio, F.; Makowski, K. Targeting Energy Expenditure-Drugs for Obesity Treatment. Pharmaceuticals 2021, 14, 435. [Google Scholar] [CrossRef] [PubMed]
  8. Kajimura, S.; Spiegelman, B.M.; Seale, P. Brown and Beige Fat: Physiological Roles beyond Heat Generation. Cell Metab. 2015, 22, 546–559. [Google Scholar] [CrossRef] [Green Version]
  9. Rosen, E.D.; Spiegelman, B.M. What we talk about when we talk about fat. Cell 2014, 156, 20–44. [Google Scholar] [CrossRef] [Green Version]
  10. Cohen, P.; Kajimura, S. The cellular and functional complexity of thermogenic fat. Nat. Rev. Mol. Cell Biol. 2021, 22, 393–409. [Google Scholar] [CrossRef]
  11. Hussain, M.F.; Roesler, A.; Kazak, L. Regulation of adipocyte thermogenesis: Mechanisms controlling obesity. FEBS J. 2020, 287, 3370–3385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Chouchani, E.T.; Kazak, L.; Spiegelman, B.M. New Advances in Adaptive Thermogenesis: UCP1 and Beyond. Cell Metab. 2019, 29, 27–37. [Google Scholar] [CrossRef] [PubMed]
  13. Fedorenko, A.; Lishko, P.V.; Kirichok, Y. Mechanism of fatty-acid-dependent UCP1 uncoupling in brown fat mitochondria. Cell 2012, 151, 400–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ikeda, K.; Yamada, T. UCP1 Dependent and Independent Thermogenesis in Brown and Beige Adipocytes. Front. Endocrinol. 2020, 11, 498. [Google Scholar] [CrossRef] [PubMed]
  15. Kazak, L.; Chouchani, E.T.; Jedrychowski, M.P.; Erickson, B.K.; Shinoda, K.; Cohen, P.; Vetrivelan, R.; Lu, G.Z.; Laznik-Bogoslavski, D.; Hasenfuss, S.C.; et al. A creatine-driven substrate cycle enhances energy expenditure and thermogenesis in beige fat. Cell 2015, 163, 643–655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Roesler, A.; Kazak, L. UCP1-independent thermogenesis. Biochem. J. 2020, 477, 709–725. [Google Scholar] [CrossRef] [PubMed]
  17. Kazak, L.; Cohen, P. Creatine metabolism: Energy homeostasis, immunity and cancer biology. Nat. Rev. Endocrinol. 2020, 16, 421–436. [Google Scholar] [CrossRef] [PubMed]
  18. Sun, Y.; Rahbani, J.F.; Jedrychowski, M.P.; Riley, C.L.; Vidoni, S.; Bogoslavski, D.; Hu, B.; Dumesic, P.A.; Zeng, X.; Wang, A.B.; et al. Mitochondrial TNAP controls thermogenesis by hydrolysis of phosphocreatine. Nature 2021, 593, 580–585. [Google Scholar] [CrossRef] [PubMed]
  19. Gámez, B.; Rodriguez-Carballo, E.; Ventura, F. BMP signaling in telencephalic neural cell specification and maturation. Front. Cell. Neurosci. 2013, 7, 87. [Google Scholar] [CrossRef] [Green Version]
  20. Mueller, T.D.; Nickel, J. Promiscuity and specificity in BMP receptor activation. FEBS Lett. 2012, 586, 1846–1859. [Google Scholar] [CrossRef] [PubMed]
  21. Katagiri, T.; Watabe, T. Bone Morphogenetic Proteins. Cold Spring Harb. Perspect. Biol. 2016, 8, a021899. [Google Scholar] [CrossRef] [Green Version]
  22. Hinck, A.P. Structural studies of the TGF-βs and their receptors—Insights into evolution of the TGF-β superfamily. FEBS Lett. 2012, 586, 1860–1870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Liu, F.; Ventura, F.; Doody, J.; Massagué, J. Human type II receptor for bone morphogenic proteins (BMPs): Extension of the two-kinase receptor model to the BMPs. Mol. Cell. Biol. 1995, 15, 3479–3486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gomez-Puerto, M.C.; Iyengar, P.V.; García de Vinuesa, A.; Ten Dijke, P.; Sanchez-Duffhues, G. Bone morphogenetic protein receptor signal transduction in human disease. J. Pathol. 2019, 247, 9–20. [Google Scholar] [CrossRef]
  25. Tseng, Y.H.; Kokkotou, E.; Schulz, T.J.; Huang, T.L.; Winnay, J.N.; Taniguchi, C.M.; Tran, T.T.; Suzuki, R.; Espinoza, D.O.; Yamamoto, Y.; et al. New role of bone morphogenetic protein 7 in brown adipogenesis and energy expenditure. Nature 2008, 454, 1000–1004. [Google Scholar] [CrossRef] [PubMed]
  26. Xue, R.; Lynes, M.D.; Dreyfuss, J.M.; Shamsi, F.; Schulz, T.J.; Zhang, H.; Huang, T.L.; Townsend, K.L.; Li, Y.; Takahashi, H.; et al. Clonal analyses and gene profiling identify genetic biomarkers of the thermogenic potential of human brown and white preadipocytes. Nat. Med. 2015, 21, 760–768. [Google Scholar] [CrossRef] [PubMed]
  27. Tóth, B.B.; Arianti, R.; Shaw, A.; Vámos, A.; Veréb, Z.; Póliska, S.; Győry, F.; Bacso, Z.; Fésüs, L.; Kristóf, E. FTO Intronic SNP Strongly Influences Human Neck Adipocyte Browning Determined by Tissue and PPARγ Specific Regulation: A Transcriptome Analysis. Cells 2020, 9, 987. [Google Scholar] [CrossRef] [Green Version]
  28. Shaw, A.; Tóth, B.B.; Király, R.; Arianti, R.; Csomós, I.; Póliska, S.; Vámos, A.; Korponay-Szabó, I.R.; Bacso, Z.; Győry, F.; et al. Irisin Stimulates the Release of CXCL1 From Differentiating Human Subcutaneous and Deep-Neck Derived Adipocytes via Upregulation of NFκB Pathway. Front. Cell Dev. Biol. 2021, 9, 737872. [Google Scholar] [CrossRef]
  29. Cypess, A.M.; White, A.P.; Vernochet, C.; Schulz, T.J.; Xue, R.; Sass, C.A.; Huang, T.L.; Roberts-Toler, C.; Weiner, L.S.; Sze, C.; et al. Anatomical localization, gene expression profiling and functional characterization of adult human neck brown fat. Nat. Med. 2013, 19, 635–639. [Google Scholar] [CrossRef]
  30. Wu, J.; Boström, P.; Sparks, L.M.; Ye, L.; Choi, J.H.; Giang, A.H.; Khandekar, M.; Virtanen, K.A.; Nuutila, P.; Schaart, G.; et al. Beige adipocytes are a distinct type of thermogenic fat cell in mouse and human. Cell 2012, 150, 366–376. [Google Scholar] [CrossRef] [Green Version]
  31. Jespersen, N.Z.; Larsen, T.J.; Peijs, L.; Daugaard, S.; Homøe, P.; Loft, A.; de Jong, J.; Mathur, N.; Cannon, B.; Nedergaard, J.; 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. [Google Scholar] [CrossRef] [Green Version]
  32. Shinoda, K.; Luijten, I.H.; Hasegawa, Y.; Hong, H.; Sonne, S.B.; Kim, M.; Xue, R.; Chondronikola, M.; Cypess, A.M.; Tseng, Y.H.; et al. Genetic and functional characterization of clonally derived adult human brown adipocytes. Nat. Med. 2015, 21, 389–394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Rahbani, J.F.; Roesler, A.; Hussain, M.F.; Samborska, B.; Dykstra, C.B.; Tsai, L.; Jedrychowski, M.P.; Vergnes, L.; Reue, K.; Spiegelman, B.M.; et al. Creatine kinase B controls futile creatine cycling in thermogenic fat. Nature 2021, 590, 480–485. [Google Scholar] [CrossRef] [PubMed]
  34. Kristóf, E.; Doan-Xuan, Q.M.; Bai, P.; Bacso, Z.; Fésüs, L. Laser-scanning cytometry can quantify human adipocyte browning and proves effectiveness of irisin. Sci. Rep. 2015, 5, 12540. [Google Scholar] [CrossRef] [PubMed]
  35. Doan-Xuan, Q.M.; Sarvari, A.K.; Fischer-Posovszky, P.; Wabitsch, M.; Balajthy, Z.; Fesus, L.; Bacso, Z. High content analysis of differentiation and cell death in human adipocytes. Cytometry A 2013, 83, 933–943. [Google Scholar] [CrossRef]
  36. Perdikari, A.; Leparc, G.G.; Balaz, M.; Pires, N.D.; Lidell, M.E.; Sun, W.; Fernandez-Albert, F.; Müller, S.; Akchiche, N.; Dong, H.; et al. BATLAS: Deconvoluting Brown Adipose Tissue. Cell Rep. 2018, 25, 784–797.e4. [Google Scholar] [CrossRef] [Green Version]
  37. Kelly, D.P.; Scarpulla, R.C. Transcriptional regulatory circuits controlling mitochondrial biogenesis and function. Genes Dev. 2004, 18, 357–368. [Google Scholar] [CrossRef] [Green Version]
  38. Puigserver, P.; Adelmant, G.; Wu, Z.; Fan, M.; Xu, J.; O’Malley, B.; Spiegelman, B.M. Activation of PPARγ coactivator-1 through transcription factor docking. Science 1999, 286, 1368–1371. [Google Scholar] [CrossRef] [PubMed]
  39. Pisani, D.F.; Barquissau, V.; Chambard, J.C.; Beuzelin, D.; Ghandour, R.A.; Giroud, M.; Mairal, A.; Pagnotta, S.; Cinti, S.; Langin, D.; et al. Mitochondrial fission is associated with UCP1 activity in human brite/beige adipocytes. Mol. Metab. 2018, 7, 35–44. [Google Scholar] [CrossRef] [PubMed]
  40. Arianti, R.; Vinnai, B.; Tóth, B.B.; Shaw, A.; Csősz, É.; Vámos, A.; Győry, F.; Fischer-Posovszky, P.; Wabitsch, M.; Kristóf, E.; et al. ASC-1 transporter-dependent amino acid uptake is required for the efficient thermogenic response of human adipocytes to adrenergic stimulation. FEBS Lett. 2021, 595, 2085–2098. [Google Scholar] [CrossRef]
  41. Kristóf, E.; Doan-Xuan, Q.M.; Sárvári, A.K.; Klusóczki, Á.; Fischer-Posovszky, P.; Wabitsch, M.; Bacso, Z.; Bai, P.; Balajthy, Z.; Fésüs, L. Clozapine modifies the differentiation program of human adipocytes inducing browning. Transl. Psychiatry 2016, 6, e963. [Google Scholar] [CrossRef] [PubMed]
  42. Klusóczki, Á.; Veréb, Z.; Vámos, A.; Fischer-Posovszky, P.; Wabitsch, M.; Bacso, Z.; Fésüs, L.; Kristóf, E. Differentiating SGBS adipocytes respond to PPARγ stimulation, irisin and BMP7 by functional browning and beige characteristics. Sci. Rep. 2019, 9, 5823. [Google Scholar] [CrossRef]
  43. Li, S.N.; Wu, J.F. TGF-β/SMAD signaling regulation of mesenchymal stem cells in adipocyte commitment. Stem Cell Res. Ther. 2020, 11, 41. [Google Scholar] [CrossRef]
  44. Schulz, T.J.; Huang, P.; Huang, T.L.; Xue, R.; McDougall, L.E.; Townsend, K.L.; Cypess, A.M.; Mishina, Y.; Gussoni, E.; Tseng, Y.H. Brown-fat paucity due to impaired BMP signalling induces compensatory browning of white fat. Nature 2013, 495, 379–383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Ma, H.; Yuan, J.; Ma, J.; Ding, J.; Lin, W.; Wang, X.; Zhang, M.; Sun, Y.; Wu, R.; Liu, C.; et al. BMP7 improves insulin signal transduction in the liver via inhibition of mitogen-activated protein kinases. J. Endocrinol. 2019, 243, 97–110. [Google Scholar] [CrossRef] [PubMed]
  46. Kazak, L.; Chouchani, E.T.; Lu, G.Z.; Jedrychowski, M.P.; Bare, C.J.; Mina, A.I.; Kumari, M.; Zhang, S.; Vuckovic, I.; Laznik-Bogoslavski, D.; et al. Genetic Depletion of Adipocyte Creatine Metabolism Inhibits Diet-Induced Thermogenesis and Drives Obesity. Cell Metab. 2017, 26, 660–671.e3. [Google Scholar] [CrossRef] [Green Version]
  47. Bertholet, A.M.; Kazak, L.; Chouchani, E.T.; Bogaczyńska, M.G.; Paranjpe, I.; Wainwright, G.L.; Bétourné, A.; Kajimura, S.; Spiegelman, B.M.; Kirichok, Y. Mitochondrial Patch Clamp of Beige Adipocytes Reveals UCP1-Positive and UCP1-Negative Cells Both Exhibiting Futile Creatine Cycling. Cell Metab. 2017, 25, 811–822.e4. [Google Scholar] [CrossRef] [Green Version]
  48. Elabd, C.; Chiellini, C.; Carmona, M.; Galitzky, J.; Cochet, O.; Petersen, R.; Pénicaud, L.; Kristiansen, K.; Bouloumié, A.; Casteilla, L.; et al. Human multipotent adipose-derived stem cells differentiate into functional brown adipocytes. Stem Cells 2009, 27, 2753–2760. [Google Scholar] [CrossRef]
  49. Petrovic, N.; Walden, T.B.; Shabalina, I.G.; Timmons, J.A.; Cannon, B.; Nedergaard, J. Chronic peroxisome proliferator-activated receptor gamma (PPARγ) 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, 7153–7164. [Google Scholar] [CrossRef] [Green Version]
  50. Jersin, R.; Tallapragada, D.S.P.; Madsen, A.; Skartveit, L.; Fjære, E.; McCann, A.; Dyer, L.; Willems, A.; Bjune, J.I.; Bjune, M.S.; et al. Role of the Neutral Amino Acid Transporter SLC7A10 in Adipocyte Lipid Storage, Obesity, and Insulin Resistance. Diabetes 2021, 70, 680–695. [Google Scholar] [CrossRef]
  51. Morawski, M.; Brückner, G.; Arendt, T.; Matthews, R.T. Aggrecan: Beyond cartilage and into the brain. Int. J. Biochem. Cell Biol. 2012, 44, 690–693. [Google Scholar] [CrossRef] [PubMed]
  52. Brückner, G.; Brauer, K.; Härtig, W.; Wolff, J.R.; Rickmann, M.J.; Derouiche, A.; Delpech, B.; Girard, N.; Oertel, W.H.; Reichenbach, A. Perineuronal nets provide a polyanionic, glia-associated form of microenvironment around certain neurons in many parts of the rat brain. Glia 1993, 8, 183–200. [Google Scholar] [CrossRef] [PubMed]
  53. Matthews, R.T.; Kelly, G.M.; Zerillo, C.A.; Gray, G.; Tiemeyer, M.; Hockfield, S. Aggrecan glycoforms contribute to the molecular heterogeneity of perineuronal nets. J. Neurosci. 2002, 22, 7536–7547. [Google Scholar] [CrossRef] [PubMed]
  54. Morawski, M.; Brückner, G.; Jäger, C.; Seeger, G.; Matthews, R.T.; Arendt, T. Involvement of perineuronal and perisynaptic extracellular matrix in Alzheimer’s disease neuropathology. Brain Pathol. 2012, 22, 547–561. [Google Scholar] [CrossRef] [Green Version]
  55. Giamanco, K.A.; Morawski, M.; Matthews, R.T. Perineuronal net formation and structure in aggrecan knockout mice. Neuroscience 2010, 170, 1314–1327. [Google Scholar] [CrossRef] [PubMed]
  56. Iwaki, T.; Kume-Iwaki, A.; Liem, R.K.; Goldman, J.E. Alpha B-crystallin is expressed in non-lenticular tissues and accumulates in Alexander’s disease brain. Cell 1989, 57, 71–78. [Google Scholar] [CrossRef]
  57. Kato, K.; Shinohara, H.; Kurobe, N.; Goto, S.; Inaguma, Y.; Ohshima, K. Immunoreactive alpha A crystallin in rat non-lenticular tissues detected with a sensitive immunoassay method. Biochim. Biophys. Acta 1991, 1080, 173–180. [Google Scholar] [CrossRef]
  58. Dubin, R.A.; Wawrousek, E.F.; Piatigorsky, J. Expression of the murine alpha B-crystallin gene is not restricted to the lens. Mol. Cell. Biol. 1989, 9, 1083–1091. [Google Scholar] [CrossRef] [Green Version]
  59. Horwitz, J. Alpha-crystallin can function as a molecular chaperone. Proc. Natl. Acad. Sci. USA 1992, 89, 10449–10453. [Google Scholar] [CrossRef] [Green Version]
  60. Rao, P.V.; Horwitz, J.; Zigler, J.S. Alpha-crystallin, a molecular chaperone, forms a stable complex with carbonic anhydrase upon heat denaturation. Biochem. Biophys. Res. Commun. 1993, 190, 786–793. [Google Scholar] [CrossRef]
  61. Raman, B.; Rao, C.M. Chaperone-like activity and quaternary structure of alpha-crystallin. J. Biol. Chem. 1994, 269, 27264–27268. [Google Scholar] [CrossRef]
  62. Shridas, P.; Sharma, Y.; Balasubramanian, D. Transglutaminase-mediated cross-linking of alpha-crystallin: Structural and functional consequences. FEBS Lett. 2001, 499, 245–250. [Google Scholar] [CrossRef] [Green Version]
  63. Groenen, P.J.; Bloemendal, H.; de Jong, W.W. The carboxy-terminal lysine of alpha B-crystallin is an amine-donor substrate for tissue transglutaminase. Eur. J. Biochem. 1992, 205, 671–674. [Google Scholar] [CrossRef] [PubMed]
  64. Mádi, A.; Cuaranta-Monroy, I.; Lénárt, K.; Pap, A.; Mezei, Z.A.; Kristóf, E.; Oláh, A.; Vámosi, G.; Bacsó, Z.; Bai, P.; et al. Browning deficiency and low mobilization of fatty acids in gonadal white adipose tissue leads to decreased cold-tolerance of transglutaminase 2 knock-out mice. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2017, 1862, 1575–1586. [Google Scholar] [CrossRef] [Green Version]
  65. Lehr, S.; Hartwig, S.; Lamers, D.; Famulla, S.; Müller, S.; Hanisch, F.G.; Cuvelier, C.; Ruige, J.; Eckardt, K.; Ouwens, D.M.; et al. Identification and validation of novel adipokines released from primary human adipocytes. Mol. Cell. Proteom. 2012, 11, M111.010504. [Google Scholar] [CrossRef] [Green Version]
  66. Hartwig, S.; De Filippo, E.; Göddeke, S.; Knebel, B.; Kotzka, J.; Al-Hasani, H.; Roden, M.; Lehr, S.; Sell, H. Exosomal proteins constitute an essential part of the human adipose tissue secretome. Biochim. Biophys. Acta Proteins Proteom. 2019, 1867, 140172. [Google Scholar] [CrossRef]
  67. Aguilera, C.M.; Gomez-Llorente, C.; Tofe, I.; Gil-Campos, M.; Cañete, R.; Gil, Á. Genome-wide expression in visceral adipose tissue from obese prepubertal children. Int. J. Mol. Sci. 2015, 16, 7723–7737. [Google Scholar] [CrossRef] [PubMed]
  68. Alfadda, A.A.; Benabdelkamel, H.; Masood, A.; Moustafa, A.; Sallam, R.; Bassas, A.; Duncan, M. Proteomic analysis of mature adipocytes from obese patients in relation to aging. Exp. Gerontol. 2013, 48, 1196–1203. [Google Scholar] [CrossRef]
  69. Benezra, R.; Davis, R.L.; Lockshon, D.; Turner, D.L.; Weintraub, H. The protein Id: A negative regulator of helix-loop-helix DNA binding proteins. Cell 1990, 61, 49–59. [Google Scholar] [CrossRef]
  70. Patil, M.; Sharma, B.K.; Satyanarayana, A. Id transcriptional regulators in adipogenesis and adipose tissue metabolism. Front. Biosci. Landmark Ed. 2014, 19, 1386–1397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Patil, M.; Sharma, B.K.; Elattar, S.; Chang, J.; Kapil, S.; Yuan, J.; Satyanarayana, A. Id1 Promotes Obesity by Suppressing Brown Adipose Thermogenesis and White Adipose Browning. Diabetes 2017, 66, 1611–1625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Zhao, Y.; Ling, F.; Griffin, T.M.; He, T.; Towner, R.; Ruan, H.; Sun, X.H. Up-regulation of the Sirtuin 1 (Sirt1) and peroxisome proliferator-activated receptor γ coactivator-1α (PGC-1α) genes in white adipose tissue of Id1 protein-deficient mice: Implications in the protection against diet and age-induced glucose intolerance. J. Biol. Chem. 2014, 289, 29112–29122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Satyanarayana, A.; Klarmann, K.D.; Gavrilova, O.; Keller, J.R. Ablation of the transcriptional regulator Id1 enhances energy expenditure, increases insulin sensitivity, and protects against age and diet induced insulin resistance, and hepatosteatosis. FASEB J. 2012, 26, 309–323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Kristóf, E.; Klusóczki, Á.; Veress, R.; Shaw, A.; Combi, Z.S.; Varga, K.; Győry, F.; Balajthy, Z.; Bai, P.; Bacso, Z.; et al. Interleukin-6 released from differentiating human beige adipocytes improves browning. Exp. Cell Res. 2019, 377, 47–55. [Google Scholar] [CrossRef]
  75. Sárvári, A.K.; Doan-Xuan, Q.M.; Bacsó, Z.; Csomós, I.; Balajthy, Z.; Fésüs, L. Interaction of differentiated human adipocytes with macrophages leads to trogocytosis and selective IL-6 secretion. Cell Death Dis. 2015, 6, e1613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Szatmári-Tóth, M.; Shaw, A.; Csomós, I.; Mocsár, G.; Fischer-Posovszky, P.; Wabitsch, M.; Balajthy, Z.; Lányi, C.; Győry, F.; Kristóf, E.; et al. Thermogenic Activation Downregulates High Mitophagy Rate in Human Masked and Mature Beige Adipocytes. Int. J. Mol. Sci. 2020, 21, 6640. [Google Scholar] [CrossRef]
Pharmaceuticals 14 01078 g001 550
Figure 1. Preadipocytes derived from subcutaneous (SC) and deep-neck (DN) depots of human neck widely expressed Bone morphogenic protein (BMP) 7 receptors and differentiated equally in the presence or absence of BMP7. SC and DN derived preadipocytes (Pre) were differentiated for two weeks following the white adipogenic differentiation protocol, BMP7 (50 ng/mL) was added where indicated. (A) Heatmap showing the expression pattern of general adipogenic differentiation marker genes and (B) quantification of differentiation rate assisted by laser-scanning cytometry. Quantification of gene expression of (C) BMPR1A, (D) BMPR1B, (E) ACVR1, (F) BMPR2, (G) ACVR2A, and (H) ACVR2B by RNA-sequencing. Data expressed as mean ± SD. n = 9.
Figure 1. Preadipocytes derived from subcutaneous (SC) and deep-neck (DN) depots of human neck widely expressed Bone morphogenic protein (BMP) 7 receptors and differentiated equally in the presence or absence of BMP7. SC and DN derived preadipocytes (Pre) were differentiated for two weeks following the white adipogenic differentiation protocol, BMP7 (50 ng/mL) was added where indicated. (A) Heatmap showing the expression pattern of general adipogenic differentiation marker genes and (B) quantification of differentiation rate assisted by laser-scanning cytometry. Quantification of gene expression of (C) BMPR1A, (D) BMPR1B, (E) ACVR1, (F) BMPR2, (G) ACVR2A, and (H) ACVR2B by RNA-sequencing. Data expressed as mean ± SD. n = 9.
Pharmaceuticals 14 01078 g001
Pharmaceuticals 14 01078 g002 550
Figure 2. BMP7 upregulated browning markers in both SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1. (A) Browning content as estimated by BATLAS Webtool (n = 9); (B) texture sum variance as quantified by laser-scanning cytometry (n = 5); quantification of UCP1 gene expression by (C) RNA-sequencing (n = 9) and (D) RT-qPCR normalized to GAPDH (n = 5). (E) UCP1 protein expression normalized to β-actin (n = 5); (F) representative UCP1 immunostaining images visualized by confocal microscopy. Nucleus was stained using propidium iodide (PI), BF represents bright field image. Scale bars represent 10 μm. (G) Quantification of UCP1 immunostaining intensity assisted by laser-scanning cytometry (n = 4000 cells of 4 donors). Data expressed as mean ± SD. * p < 0.05, ** p < 0.01, *** p < 0.001. Statistics: Paired t-test (A,B,G), GLM (C), and one-way ANOVA with Tukey’s post-test (D,E).
Figure 2. BMP7 upregulated browning markers in both SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1. (A) Browning content as estimated by BATLAS Webtool (n = 9); (B) texture sum variance as quantified by laser-scanning cytometry (n = 5); quantification of UCP1 gene expression by (C) RNA-sequencing (n = 9) and (D) RT-qPCR normalized to GAPDH (n = 5). (E) UCP1 protein expression normalized to β-actin (n = 5); (F) representative UCP1 immunostaining images visualized by confocal microscopy. Nucleus was stained using propidium iodide (PI), BF represents bright field image. Scale bars represent 10 μm. (G) Quantification of UCP1 immunostaining intensity assisted by laser-scanning cytometry (n = 4000 cells of 4 donors). Data expressed as mean ± SD. * p < 0.05, ** p < 0.01, *** p < 0.001. Statistics: Paired t-test (A,B,G), GLM (C), and one-way ANOVA with Tukey’s post-test (D,E).
Pharmaceuticals 14 01078 g002
Pharmaceuticals 14 01078 g003 550
Figure 3. BMP7 treatment increased mitochondrial biogenesis and protein expression of mitochondrial OXPHOS complex subunits in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1 and Figure 2. (A) Protein expression of PGC1α normalized to β-actin, (B) Ratio of p-CREB to CREB protein expression, and (C) protein expression of mitochondrial OXPHOS complexes (I, II, III, IV, and V) normalized to β-actin. Data expressed as mean ± SD. n = 5. * p < 0.05, ** p < 0.01. Statistics: one-way ANOVA with Tukey’s post-test.
Figure 3. BMP7 treatment increased mitochondrial biogenesis and protein expression of mitochondrial OXPHOS complex subunits in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1 and Figure 2. (A) Protein expression of PGC1α normalized to β-actin, (B) Ratio of p-CREB to CREB protein expression, and (C) protein expression of mitochondrial OXPHOS complexes (I, II, III, IV, and V) normalized to β-actin. Data expressed as mean ± SD. n = 5. * p < 0.05, ** p < 0.01. Statistics: one-way ANOVA with Tukey’s post-test.
Pharmaceuticals 14 01078 g003
Pharmaceuticals 14 01078 g004 550
Figure 4. Fragmented and total mitochondrial content were increased upon BMP7 treatment in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2 and Figure 3. (A) Representative confocal microscopy images of TOM20 immunostaining. PI was used to stain the nucleus. Scale bars represent 10 μm. (B) Quantification of fragmented mitochondrial content normalized to per nuclei (n = 55 cells of 3 donors); (C) total mitochondrial DNA content quantified by qPCR (n = 5). Data expressed as mean ± SD. * p < 0.05, ** p < 0.01, **** p < 0.0001. Statistics: one-way ANOVA with Tukey’s post-test.
Figure 4. Fragmented and total mitochondrial content were increased upon BMP7 treatment in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2 and Figure 3. (A) Representative confocal microscopy images of TOM20 immunostaining. PI was used to stain the nucleus. Scale bars represent 10 μm. (B) Quantification of fragmented mitochondrial content normalized to per nuclei (n = 55 cells of 3 donors); (C) total mitochondrial DNA content quantified by qPCR (n = 5). Data expressed as mean ± SD. * p < 0.05, ** p < 0.01, **** p < 0.0001. Statistics: one-way ANOVA with Tukey’s post-test.
Pharmaceuticals 14 01078 g004
Pharmaceuticals 14 01078 g005 550
Figure 5. BMP7 treatment resulted in increased cellular respiration, thermogenic response, and extracellular acidification in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2, Figure 3 and Figure 4. (A) Representative time lapse curve for measured oxygen consumption rate (OCR) after indicated treatments, followed by quantification of basal, stimulated, and proton leak OCR; (B) representative time lapse curve for measured extracellular acidification rate (ECAR) after the indicated treatment, followed by quantification of basal and stimulated ECAR. Data expressed as mean ± SD. n = 4. * p < 0.05, ** p < 0.01. Statistics: one-way ANOVA with Tukey’s post-test.
Figure 5. BMP7 treatment resulted in increased cellular respiration, thermogenic response, and extracellular acidification in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2, Figure 3 and Figure 4. (A) Representative time lapse curve for measured oxygen consumption rate (OCR) after indicated treatments, followed by quantification of basal, stimulated, and proton leak OCR; (B) representative time lapse curve for measured extracellular acidification rate (ECAR) after the indicated treatment, followed by quantification of basal and stimulated ECAR. Data expressed as mean ± SD. n = 4. * p < 0.05, ** p < 0.01. Statistics: one-way ANOVA with Tukey’s post-test.
Pharmaceuticals 14 01078 g005
Pharmaceuticals 14 01078 g006 550
Figure 6. BMP7 treatment upregulated creatine futile cycle-mediated thermogenesis in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. (A) Representative time lapse OCR curve after indicated treatments, followed by quantification of creatine cycle related (β-GPA inhibited) OCR (n = 5). (B) Quantification of gene expression of CKMT2 by RNA-sequencing (n = 9) and RT-qPCR (n = 5), followed by protein expression (n = 5). (C) Quantification of gene expression of CKB by RNA-sequencing (n = 9), followed by protein expression (n = 4). (D) Quantification of gene expression of TNAP by RNA-sequencing (n = 9). Data expressed as mean ± SD. RT-qPCR was normalized to GAPDH, protein expression was normalized to β-actin. * p < 0.05, ** p < 0.01. Statistics: one-way ANOVA with Tukey’s post-test.
Figure 6. BMP7 treatment upregulated creatine futile cycle-mediated thermogenesis in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. (A) Representative time lapse OCR curve after indicated treatments, followed by quantification of creatine cycle related (β-GPA inhibited) OCR (n = 5). (B) Quantification of gene expression of CKMT2 by RNA-sequencing (n = 9) and RT-qPCR (n = 5), followed by protein expression (n = 5). (C) Quantification of gene expression of CKB by RNA-sequencing (n = 9), followed by protein expression (n = 4). (D) Quantification of gene expression of TNAP by RNA-sequencing (n = 9). Data expressed as mean ± SD. RT-qPCR was normalized to GAPDH, protein expression was normalized to β-actin. * p < 0.05, ** p < 0.01. Statistics: one-way ANOVA with Tukey’s post-test.
Pharmaceuticals 14 01078 g006
Pharmaceuticals 14 01078 g007 550
Figure 7. Genes upregulated upon BMP7 treatment in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. (A) Volcano plot illustrating each of the upregulated genes in SC (red) and DN (blue) depots upon BMP7 treatment. (B) Venn-diagram showing the genes commonly up- or downregulated by BMP7 in SC and DN derived adipocytes. Gene expression of (C) ACAN, (D) CRYAB, and (E) ID1 as determined by RNA-sequencing (n = 9) and RT-qPCR (n = 5), followed by their respective protein expression (n = 5). Data presented as mean ± SD. RT-qPCR was normalized to GAPDH, protein expression was normalized to β-actin. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Statistics: GLM (RNA-sequencing), one-way ANOVA with Tukey’s test (RT-qPCR and densitometry of western blots).
Figure 7. Genes upregulated upon BMP7 treatment in SC and DN derived differentiated adipocytes. SC and DN preadipocytes were differentiated and treated as in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. (A) Volcano plot illustrating each of the upregulated genes in SC (red) and DN (blue) depots upon BMP7 treatment. (B) Venn-diagram showing the genes commonly up- or downregulated by BMP7 in SC and DN derived adipocytes. Gene expression of (C) ACAN, (D) CRYAB, and (E) ID1 as determined by RNA-sequencing (n = 9) and RT-qPCR (n = 5), followed by their respective protein expression (n = 5). Data presented as mean ± SD. RT-qPCR was normalized to GAPDH, protein expression was normalized to β-actin. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Statistics: GLM (RNA-sequencing), one-way ANOVA with Tukey’s test (RT-qPCR and densitometry of western blots).
Pharmaceuticals 14 01078 g007
Table
Table 1. Pathways of genes significantly upregulated upon BMP7 treatment in SC and DN derived differentiated adipocytes. Genes commonly upregulated in SC and DN derived adipocytes are marked red. False discovery rate (FDR) values are indicated.
Table 1. Pathways of genes significantly upregulated upon BMP7 treatment in SC and DN derived differentiated adipocytes. Genes commonly upregulated in SC and DN derived adipocytes are marked red. False discovery rate (FDR) values are indicated.
SC BMP7 Upregulated
Reactome PathwaysGene NameFDR
Signaling by BMPNOG, SMAD9, SMAD6, SMAD72.49 × 10−2
ECM proteoglycansITGA9, AL645608, ITGA8, ACAN, COMP, COL9A3, COL9A21.87 × 10−3
Integrin cell surface interactionsITGA9, AL645608, ITGA8, COMP, COL9A3, COL9A22.3 × 10−2
GPCR ligand bindingPLPPR4, NMUR1, ADRA2C, FZD5, ACKR1, ADRA2A, PTH1R, NPY1R, CHRM4, FZD1, EDNRA, PTGDR, CNR12.93 × 10−2
DN BMP7 Upregulated
Reactome PathwaysGene NameFDR
Signaling by BMPNOG, SMAD9, SMAD6, SMAD76.05 × 10−3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shaw, A.; Tóth, B.B.; Arianti, R.; Csomós, I.; Póliska, S.; Vámos, A.; Bacso, Z.; Győry, F.; Fésüs, L.; Kristóf, E. BMP7 Increases UCP1-Dependent and Independent Thermogenesis with a Unique Gene Expression Program in Human Neck Area Derived Adipocytes. Pharmaceuticals 2021, 14, 1078. https://doi.org/10.3390/ph14111078

AMA Style

Shaw A, Tóth BB, Arianti R, Csomós I, Póliska S, Vámos A, Bacso Z, Győry F, Fésüs L, Kristóf E. BMP7 Increases UCP1-Dependent and Independent Thermogenesis with a Unique Gene Expression Program in Human Neck Area Derived Adipocytes. Pharmaceuticals. 2021; 14(11):1078. https://doi.org/10.3390/ph14111078

Chicago/Turabian Style

Shaw, Abhirup, Beáta B. Tóth, Rini Arianti, István Csomós, Szilárd Póliska, Attila Vámos, Zsolt Bacso, Ferenc Győry, László Fésüs, and Endre Kristóf. 2021. "BMP7 Increases UCP1-Dependent and Independent Thermogenesis with a Unique Gene Expression Program in Human Neck Area Derived Adipocytes" Pharmaceuticals 14, no. 11: 1078. https://doi.org/10.3390/ph14111078

APA Style

Shaw, A., Tóth, B. B., Arianti, R., Csomós, I., Póliska, S., Vámos, A., Bacso, Z., Győry, F., Fésüs, L., & Kristóf, E. (2021). BMP7 Increases UCP1-Dependent and Independent Thermogenesis with a Unique Gene Expression Program in Human Neck Area Derived Adipocytes. Pharmaceuticals, 14(11), 1078. https://doi.org/10.3390/ph14111078

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop