Skip to main content
Download PDF

Norepinephrine triggers an immediate-early regulatory network response in primary human white adipocytes

BMC Genomics volume 19, Article number: 794 (2018) Cite this article

Abstract

Background

Norepinephrine (NE) signaling has a key role in white adipose tissue (WAT) functions, including lipolysis, free fatty acid liberation and, under certain conditions, conversion of white into brite (brown-in-white) adipocytes. However, acute effects of NE stimulation have not been described at the transcriptional network level.

Results

We used RNA-seq to uncover a broad transcriptional response. The inference of protein-protein and protein-DNA interaction networks allowed us to identify a set of immediate-early genes (IEGs) with high betweenness, validating our approach and suggesting a hierarchical control of transcriptional regulation. In addition, we identified a transcriptional regulatory network with IEGs as master regulators, including HSF1 and NFIL3 as novel NE-induced IEG candidates. Moreover, a functional enrichment analysis and gene clustering into functional modules suggest a crosstalk between metabolic, signaling, and immune responses.

Conclusions

Altogether, our network biology approach explores for the first time the immediate-early systems level response of human adipocytes to acute sympathetic activation, thereby providing a first network basis of early cell fate programs and crosstalks between metabolic and transcriptional networks required for proper WAT function.

Background

White adipose tissue (WAT) is a multifunctional organ that governs energy storage, endocrine functions and signaling for maintaining energy homeostasis in the body [1]. WAT also has the capacity to expand in response to caloric intake, hormones, and in the aging process [2]. As a counterpart to WAT, brown adipose tissue (BAT), which is located predominantly in the interscapular area in human adults, maintains thermoregulation of the body during acute or prolonged cold exposure. It is equipped with a high density of mitochondria and multilocular lipid droplets [3]. The thermogenic function of BAT relies mainly on UCP1, a mitochondrial protein that uncouples oxidative phosphorylation from ATP synthesis, leading to energy dissipation [4]. Recently, the remodeling of white adipocytes to brown-like fat cells (brite) has been reported, with characteristics such as UCP1 expression and thermogenesis that resemble BAT [5]. The brite adipocyte upraise has been observed upon exposure to cold or in response to different stimuli, such as chronic exposition to rosiglitazone (Rosi), a PPARγ agonist [6], Celastrol, a plant-derived triterpene [7], and norephineprhine (NE) [8].

NE is an integral part of the sympathetic nervous system and mediates essential physiologic responses, including increased heart rate and blood pressure, mobilization of energy stores and control of core body temperature [9]. NE exerts its effects by binding to adrenergic receptors α and β, linked to GS proteins [10], which in turn are linked to an adenylate cyclase. NE binding thus causes a rise in the intracellular concentration of cyclic AMP (cAMP). Downstream effectors of cAMP include cAMP dependent protein kinase (PKA), which mediates most of the so far known intracellular events following NE binding [11]. In WAT, it is well known that adrenergic stimuli trigger lipolysis and mobilization of free fatty acids [12]. In BAT, NE-activated PKA can also phosphorylate p38 MAPK, which activates the transcription of UCP1 by phosphorylating the PPARγ coregulator 1α (PGC1α) codified by the gene PPARGC1A, and the transcription factor ATF [13,14,15].

It is reasonable to conclude that, in response to cold, signaling by NE can induce well-differentiated cellular programs in different adipose tissues. In BAT, NE triggers a thermogenic program dependent on a signaling cascade that leads to elevated UCP1 expression and activation of UCP1 by lipolysis [16]. In WAT two programs are carried out: lipolysis leads to fat mobilization, immediately executed and dependent almost exclusively on a phosphorylation cascade of existing proteins [17]. In addition, conversion of WAT to the brown-like phenotype implies a major transcriptional shift, caused by a profound remodeling of the superenhancers responsible for the maintenance of adipogenesis [18], the down-regulation of pro-adipogenic transcription factors (TFs) such as PPARγ, and the initiation of the transcriptional program that triggers thermogenesis [19]. However, while long-term NE-triggered phenotypical and physiological effects are known, immediate-early transcriptional responses to NE are still obscure.

Cell-extrinsic signals can activate a specific set of genes, called immediate-early genes (IEGs), which are transcribed within minutes after stimulation, are expressed in waves without the need for de novo protein synthesis [20, 21], and are transcription factors (TFs) that can control target gene expression. IEGs have emerged to respond to a variety of extrinsic stimuli in multiple cell types, indicating a common response mechanism probably comprising several hundred genes [22]. Nevertheless, little is known about IEGs that are triggered by NE in adipocytes.

The coordination of regulatory mechanisms becomes critical for an accurate gene expression pattern in biological processes, and this transcriptional regulation is structured into a hierarchical organization with regulators at different levels exhibiting unique characteristics [23]. Standard transcriptome analyses typically uncover expression changes for hundreds or thousands of genes. However, without additional system level approaches, they are not able to reveal the complex mechanisms behind changes in gene expression. Simplistic approaches have failed so far to characterize the complex cellular response defined by the collective contribution of regulatory and signaling pathways [24, 25]. As no single TF is sufficient to drive a complex regulatory process alone, changes in cellular states are determined by complex networks, involving both positive and negative regulatory interactions with a substantial number of TFs [26]. In the context of adipocytes, the early complex cellular response to NE and its underlying transcriptional regulatory network still awaits to be studied.

In order to get insight into this complexity, it is necessary to study molecules in a network context, including protein-protein interactions (PPI) and metabolic, signaling and transcriptional regulatory networks (TRN) [27]. The sum of all these interactions, the cellular network, allows to elucidate and visualize complex interactions and their information flow, where molecules are represented as nodes and their interaction as edges [28]. The analysis of structural network components and their behavior using topology offers a quantifiable measurement that provides insight into biological functions of cellular networks. The analysis of betweenness centrality, for example, allows the identification of highly interconnected nodes known as “bottlenecks”, which may point toward their essentiality [29, 30]. Therefore, the complex relationships between components can only be elucidated from a network perspective [31]. Moreover, integration of functional clustering and analysis of network topology paves the way to disclose new associations among genes cooperating in undisclosed, not yet annotated biological processes [32, 33].

Here, using these tools and concepts, we investigated the early transcriptional response of human white adipocytes to 3 h of NE stimulation. First, we found a broad transcriptional response with more than 2,100 differentially expressed genes assigned to a wide range of pathways. Second, we inferred the transcriptional regulatory network of these NE-responsive genes, which indicated that the known and novel immediate-early genes might serve as master and local regulators. Third, we performed a functional enrichment analysis that suggested the participation of several genes as a “bridge” between functional clusters. These transcriptional and post-transcriptional regulators and their target genes provide a novel insight into the circuitry and functional principles of the acute response to NE in human white adipocytes.

Results

NE stimulation triggers an acute and broad transcriptional response in primary human white adipocytes

In order to investigate the early response of primary human adipocytes to NE stimulation, we obtained human primary adipose-derived stem cells (hpASCs) from four female individuals. These hpASCs were exposed to an adipogenic cocktail, permitting the development of terminally differentiated white adipocytes within 9 days. These mature adipocytes were then treated with NE (i.e. stimulated) or vehicle (i.e. unstimulated), harvested after 3 h, and used to perform RNA-seq analysis (Fig. 1a).

Fig. 1
figure1

NE stimulation triggers an acute and broad transcriptional response in human adipocytes. a Workflow diagram for sample preparation. Adipocyte precursor cells were isolated from donors (n = 4), differentiated and treated with NE or vehicle for 3 h, harvested and subsequently analyzed by RNA-seq (The female human silhouette was modified from a clipart with a public domain license creative commons CC0). b Principal component analysis (PCA) and unsupervised hierarchical clustering (c) showed clustering between stimulated (orange) and unstimulated (blue) adipocytes in two distinct groups, suggesting that NE stimulation causes deep changes in the transcriptome. Columns are biological replicates and rows represent differentially expressed genes. d Reactome pathway enrichment for up- and e down-regulated DEGs after acute NE stimulation

Full size image

We performed a principal component analysis (PCA) of the RNA-seq data, which revealed that the expression signatures from the stimulated adipocytes separate from the unstimulated ones in each biological replicate (Fig. 1b). Differential gene expression analysis of stimulated compared to unstimulated adipocytes showed 2,160 differentially expressed genes (DEGs) that were significantly up-regulated (1,239) or down-regulated (921) (Fig. 1c). Among the up-regulated genes we found CREM [34] (1.5 log2FC) and PPARGC1A [35] (2.5 log2FC), both known to be up-regulated in response to NE, thereby validating our experimental approach.

In order to ascertain which pathways were acutely stimulated in white adipocytes upon NE treatment, we performed a pathway enrichment analysis using the Reactome database. Up-regulated genes were enriched for circadian cycle and gene-related signaling pathways (Fig. 1d), primarily nerve growth factor (NGF), platelet-derived growth factor (PDGF), VEGF, and WNT/Beta-catenin. On the other hand, enriched pathways for the down-regulated genes comprise p53 downstream targets, chromatin organization, apoptosis regulators, GAB1 signalosome super pathway (which includes PI3K/AKT/LAT2/LAB pathways), signaling by TGFβ receptor complex, and targets of SMAD heterodimers (Fig. 1e). The large number of alterations in gene expression suggests an acute transcriptional response, while the finding of many different modulated pathways postulates a broad transcriptional response to NE in primary human white adipocytes.

Acute NE responsive network has immediate-early genes as nodes with high betweenness

In order to gain insight into the molecular network triggered acutely by NE stimulation, we first constructed the experimentally verified interactions within up- and down-regulated DEGs by inferring their PPI and Protein-DNA networks. Second, we explored the connection between network structure and global network topology. The betweenness of all the nodes in the network was calculated using the Cytoscape [36] plug-in cytoHubba [37]. The nodes with the highest betweenness values in the network of up-regulated genes were JUN, FOS, NCOR2, FOXO1, CEBPB, TLE1, CREM, ATF3, NR4A1, EPHB2, HIPK2 and RHOB (Fig. 2a, Additional file 1: Table S1). JUN and FOS are described as IEGs [38], known to be transcribed in other models, during adipogenesis [39] and upon adrenergic stimulation [40]. However, the participation of these TFs in acute NE response in mature adipocytes was not yet known until now. In the network of down-regulated genes, the top nodes – ranked by betweenness – were SMAD3, FYN, RUNX2 and SOS1 (Fig. 2a, Additional file 1: Table S1). SMAD3 [41] and RUNX2 [42] are known to participate in adipocyte differentiation from precursors but not in early NE stimulation.

Fig. 2
figure2

Acute NE responsive network has immediate-early genes as nodes with high betweenness. a The interaction network reconstructed from DEGs shows nodes with high betweenness. Edges represent experimentally verified protein-protein or protein-DNA interactions. Up-regulated DEGs are in red and down-regulated DEGs in green. Size node is related to betweenness value. b Known immediate-early genes were found to be up-regulated after NE stimulation. Expression values are shown in log2 counts per million (CPM). Expression differences between conditions were analyzed with a Kruskall-Wallis Test, P < 0.0001, Dunn’s Multiple Comparison Test P < 0.05; ET Test, P < 0.05 (*** p < 0.0001, ** p < 0.001)

Full size image

Since JUN and other IEGs were up-regulated genes with higher betweenness within the network, we measured expression levels of JUN and other known IEGs such as JUND, ATF3, FOS and FOSL2 in detail. These genes were up-regulated in adipocytes stimulated with NE, demonstrating that NE stimulation activates IEG transcription (Fig. 2b).

All together, these results shows that several genes with high betweenness are IEGs. These findings suggest a transcriptional regulatory network organized by highly connected genes, identified as IEGs in the very early response to NE stimulation.

Transcriptional regulatory network analysis reveals immediate-early genes responsive to NE stimulation

Since we observed altered transcription of known IEGs, we were interested in getting insight into the emergent properties of the acute NE-treated adipocyte regulation by looking into the TRN. Therefore, we classified the TFs according to hierarchy as global (“master”) regulators (MRs) or mid-level (“local”) regulators (LRs). MRs are genes at the top of the gene regulatory hierarchy that regulate multiple downstream genes either directly or through a cascade of gene expression changes, and have the ability to re-direct the fate of cells [43]. LRs are defined as those TFs that are regulated by MRs, and are mostly dedicated to regulate a specific set of downstream targets [44]. In order to identify the regulatory hierarchy in our TRN, a TF motif activity prediction was performed on the networks of our significantly up- and down-regulated genes using iRegulon [45]. TFs with high, normalized enrichment scores (NES ≥ 3) were identified in each network, corresponding to an estimated false discovery rate of less than 0.01, which we used as threshold for our dataset. In the transcriptional regulatory network of up-regulated TFs we found a total of 147 overrepresented TFs, 60 of them with differentially expressed targets, comprising 10 MRs and 50 LRs (Additional file 1: Table S2). In the regulatory network of down-regulated TFs, we found a total of 143 TFs, 81 of them with differentially expressed targets, comprising four MRs and 77 LRs (Additional file 1: Table S3).

Due to the fact that IEGs do not require de novo protein synthesis for their expression, that they acutely respond to NE stimulation, and that they are TFs [22], we assume that our MRs may be IEGs. Several MRs, which we identified in our dataset as up-regulated, are widely acknowledged as IEGs in general, like FOS, JUN, JUND, ATF3, CREM, CEBPA, and CEBPB (previous section and Fig. 3a), thus serving as validation of our approach. Several of these MRs are known to be involved in WAT conversion to the brite phenotype, such as CREM, CEBPA, and CEBPB [46, 47]. Moreover, in our TRN we also identified TFs which are not yet recognized as IEGs, but were acutely responsive to NE, such as the heat shock transcription factor 1 (HSF1) and the nuclear factor, interleukin 3-regulated (NFIL3). Interestingly, the highest ranked MR (NES = 5.5) was HSF1, a TF involved in the transcriptional activation of the heat shock response (HSR) [48] and known to regulate energy expenditure through activation of a PGC1α-dependent metabolic program [7]. Another highly ranked MR is NFIL3, a transcriptional regulator known for its role in circadian rhythm [49].

Fig. 3
figure3

Regulatory network analysis reveals a set of immediate-early genes that are acutely responsive to NE stimulation. Integrated view of relationships between master (1) and local (2) regulators, their expression (3) and their predicted targets (4) mapped to the chromosomal context (5). Regulatory edges link the master regulators to their targets, which were colored as indicated in the figure. Up-regulated (a) and down-regulated (b) regulatory networks, based on acute responses to NE stimulation, revealed master regulators with different targets, as shown in histograms

Full size image

Moreover, the TRN of down-regulated genes identified four TFs as MRs, with ZBTB14 as the highest ranked (NES = 4.6), followed by FOXP1, SP1 and FOXO3 (Fig. 3b). Several pro-adipogenic factors, like PPARG, HOXC8, HOXB3, HOXA9, RARG and KLF5, are LRs downregulated upon NE treatment, suggesting a shutdown of the adipogenic program.

Thus, our results indicate that several TFs, including MRs and LRs, form a TRN that participates in the regulatory program of white adipocytes, with a shutdown of pro-adipogenic TFs and the initiation of pleiotropic genes capable of altering cell fate, including a number of factors known to be involved in the remodeling of white adipocytes. Moreover, our results pinpoint HSF1 and NFIL3 as promising novel IEGs acutely triggered by NE in primary human white adipocytes.

Functional enrichment analysis elucidates a complex response to NE in metabolic and signaling pathways

Our results revealed that NE stimulation up-regulated several IEGs that we also identified as MRs and LRs with pleiotropic function. Thus, we speculated whether they participate in the crosstalk between pathways and performed a functional gene enrichment analysis of the up-regulated DEGs using the R/Bioconductor package FGNet [50], reconstructing the functional gene network from up-regulated genes. FGNet which also organized the network into metagroups that condense multiple KEGG pathways, Reactome pathways, and GO terms (Fig. 4). The genes clustered in four metagroups. These metagroups contain genes classified according to their functions and seven clusters of common genes that are part of two or more metagroups. In each cluster, there are TFs, adaptor proteins, kinases as well DNA binding proteins (Table 1).

Fig. 4
figure4

Functional enrichment analysis suggests a complex changes in metabolic and signaling pathways upon acute NE stimulation. Functional gene network analysis of up-regulated genes showing metagroup enrichment in distinct cellular functions as represented by colored nodes and clusters. White nodes represent genes that connect two or more functional modules. Metagroup compositions are depicted in Table 1

Full size image
Table 1 Functional module composition. Each metagroup represents cellular functions, each represented by a specific color and composed of several gene clusters. A total of four metagroups and seven clusters were identified
Full size table

Because numerous genes tightly interconnect the four metagroups, our results suggest a crosstalk between pathways after NE stimulation. In conclusion, our approach provides hypotheses regarding the connections between functional metagroups. These connections can be used as the basis for investigating the mechanisms that link these functionally related genes. We expect that as more data become available, these networks will be further refined and expanded to provide a more in-depth insight into the regulatory network that drives biological function.

Discussion

Despite the key role of NE in the sympathetic nervous system, acute transcriptional processes unleashed by NE in human models have been barely studied. In the present work, we investigated the acute response to NE stimulation in the transcriptional landscape of human primary white adipocytes with the aim to identify immediate-early genes and their downstream network, as well as their relationship to biological pathways that may ultimately explain phenotypical changes that have been described as response to NE in previous studies.

Primary adipose-derived stem cells from four donors have been studied. In order to diminish this limitation in available donor samples, hpASC from donors were randomized and treated as repetitions, to dilute the batch effect and the differences between individuals. Our approach allowed us to identify significant early responders to NE, with some of them being known responders that serve as positive controls, and with some of them identified in this study as novel responders to NE. However, because adipocytes from female individuals only were studied, conclusions on sex-dependent responses to NE are beyond the scope of this study. Despite these limitations, the results of our study suggest important changes on the cellular transcriptional network, and we identified known and novel IEGs and pathways to be involved in early NE response.

Another issue is the contamination of adipocyte populations with macrophages when adipocytes are directly isolated from tissue. In order to avoid this problem, we isolated primary adipose-derived stem cells (hpASCs) that were subsequently differentiated into white adipocytes to ensure a majority population of adipocytes. Moreover, after sequencing, we searched for the expression of several macrophage markers, but their expression levels were magnitudes lower than those of adipocyte markers (Additional file 2: Figure S1) so that we can assume that the contribution of macrophages to the expression signature in our dataset is negligible in comparison to adipocytes.

Over the last decades, several studies focused on lipid and glucose mobilization triggered by NE in WAT through the ADRβ/cAMP/PKA/PI3K axis [51], and it was assumed that these functions were performed mostly without the need for “de novo” transcription. However, our results show that several TFs are acutely transcribed upon NE treatment, together with more than 2,000 DEGs, indicating an acute shift in gene regulation and cell fate after NE stimulation.

Transcriptional regulation in adipocytes has been intensely studied in recent years on brite adipocytes [52, 53]. Several lines of evidence indicate in mouse [54, 55] and human [56] that chronic adrenergic stimulation or agonistic activity triggers conversion of white adipocytes to a brown-like phenotype with thermogenic activity due to high expression and functionality of UCP1 [57, 58]. Most adipocyte studies have analyzed adipocyte responses to stimulation by NE or agonistic activity [59] after 24 to 72 h in mouse or cell systems, and after weeks in humans [55, 56], but not after 3 h as we did. Therefore, it was surprising that we found some of the considered brite adipocyte marker genes, such as CITED1, HOXA9, DIO2 or PPARGC1A [60] already being up-regulated at such an early stage of NE stimulation. In line with that, transcription of several white adipocyte marker genes, such as PPARG, RUNX1T1 or HOXC8 [61, 62] were repressed, as well as genes involved in osteogenesis, such as SMAD3, 5 and 9 [63, 64] (Additional file 3: Table S4). These results suggest, provocatively, that NE-driven conversion starts already at an earlier stage than expected.

At pathway level, we found an up-regulation of growth factors such as NGF, PDGF, VEGF, WNT and circadian clock factors. NGF, PDGF and VEGF are well-known growth factors that share a common signaling cascade, the Akt/mTOR/MAPK axis [65]. Intriguingly, there are lines of evidence which indicate that these factors are critical for cell fate and tissue dynamics [66,67,68].

On the other side of the spectrum, NE stimulation led to a down-regulation of TP53, chromatin organization, and programmed cell death pathways. In principle, these findings are not surprising, because it is necessary to shut down the cell padlocks in order to initiate a change in the cellular program [69, 70] and to avoid the induction of cell death [71].

The fact that not a large number of DEGs in our study could be assigned to gene and pathway annotations, coincides with an increasing body of evidence suggesting that canonical pathways are incomplete and inaccurate models to study the complex interplay of signal transduction, transcriptional, post-transcriptional, metabolic and other regulatory events that determine cell behavior [72]. A possible and widely used solution to analyze information flow in systems biology is network construction, since networks act as a paradigm for data integration and analysis, providing a systems-level understanding of the mechanisms underlying cell biology and disease [73, 74]. The interaction network of DEGs upon NE stimulation pinpointed several highly connected nodes based on betweenness. We could identify genes with several functions: those codifying for signaling transduction proteins (MAPK14, SOS1, ABL1), structural proteins (ACTB), cell cycle (CDK1, BCL2), metabolic signaling proteins (PIK3CA, AKT, mTOR), post-translational modifiers (UBE2D1, YWAHZ), and specific TFs (NCOR2, RELA, SMAD3, FOXO1, JUN, FOS, YAP1). Nodes with high betweenness have been demonstrated to play a central role in biological networks [75,76,77], and are frequently recognized as pleiotropic or essential genes, and their differential expression can predict information flow [78].

Regulatory network analysis based on TF motif activity has been used before ranking and dissecting the influence of TFs in biological networks [79]. Just as in interaction networks, bottlenecks in regulatory networks are highly related to their role as essential proteins with the addition of having implicit flow, since their edges are directed. Our study identified several TFs that could play a role as master and local regulators in white adipocytes upon NE stimulation, many of them being well-known IEGs (for instance, FOS, JUN, JUND, ATF3 or CEBPA/B). FOS and JUN are responsive to adrenergic stimulation in rat cardiomyocytes [40] suggesting a conserved regulatory response to NE between rat and human, in cardiomyocytes and adipocytes. Nevertheless, the bulk of knowledge about IEGs and their regulation is based on studies of neurons and the nervous system [80,81,82]. In those models, the main functions described are proliferation, differentiation, survival [83], and stem cells fate [84].

An unexpected result of our TRN analysis was the identification of HSF1 and NFIL3 as novel IEGs and MRs in white adipocytes. HSF1 is a TF that plays a central role in the transcriptional activation of the heat shock response (HSR), leading to the expression of heat shock proteins. Notably, HSF1 has previously been identified as an important factor in white to brite conversion after seven days of stimulation with celastrol [7], and its ablation impairs energy metabolism [85, 86]. Moreover, it has been described as a central regulator of cellular bioenergetics and protein homeostasis in liver [87]. Thus, our work pinpoints HSF1 as a novel IEG candidate that could regulate adipocyte energetics upon NE stimulation at very early stages. NFIL3 is a transcriptional regulator, mainly known for being involved in circadian regulation [88] and shutting down apoptotic signals promoting survival and regeneration [89]. While it is a known responder to adrenergic stress in fibroblasts, it has also been demonstrated to be induced by activation of the α1-adrenergic receptor and does not involve the cAMP pathway [90, 91]. Taken together, our results suggest that the acute response to NE stimulation is coordinated by several adrenergic receptors and pathways.

We performed a functional enrichment network analysis in order to identify specific functional modules and a crosstalk between identified pathways. The identification of genes that work as “bridges” or linkers between the four identified metagroups reveal not only TFs but also other molecules, such as kinases, receptors and even structural genes that can be critical regulators of cell processes. Although it is true that we do not have a quantitative way to support whether these metagroups have an impact to drive the phenotype, biological evidence shows that several genes participating in crosstalk make sense in cellular reprogramming and cell fate, as BCR on B cells [92], CCND1 on epidermal cells [93] or mTOR on pancreatic cells [94]. High redundancy in the signals across several functions, like immune response or signaling events was observed, where diverse clusters of transcriptional and post-transcriptional regulators seem to control a robust network. Our results not only allow the inference of transcriptional regulation, but also of post-transcriptional effectors, such as the YWHAZ, a gene encoding for the 14–3-3ζ protein, a signal adapter capable of functioning as a kinase and allowing the flow of information in signaling cascades on effectors like PI3K [95], mTOR or AKT [96].

Conclusions

We identified a broad transcriptional response of primary human white adipocytes to acute NE stimulation. With our reconstruction of functional, interaction, and transcriptional networks, we identified novel NE-responsive bottleneck genes with high betweenness, novel NE-responsive immediate-early gene candidates, and a complex response to NE in metabolic and signaling pathways. Our study provides a basis for hypothesis-driven studies on how the activation of IEGs orchestrates the shift in the transcriptional program after acute NE stimulation, and the amplitude of its effects in the cellular network of adipocytes. In future studies, the use of a cohort can help to infer possible sex differences and stratification to NE response between populations. Studies with human samples offer particular challenges; nevertheless, functional genomic assays could lead to an impact on a variety of practical applications such as the identification of novel candidate genes capable of modulating cell fate in response to extrinsic signals, the discovery and exploration of non-canonical pathways, and a deeper insight into crosstalks between metabolic and transcriptional networks

Methods

Isolation and culture of hpASCs

Human primary adipose-derived stem cells (hpASCs) were isolated from subcutaneous lipoaspirates from healthy female donors (n = 4) as previously described [97]. Cells at P0 were thawed, cultivated in EGM-2 Medium (Lonza), and used after 1–3 passages. For adipocyte differentiation, cells were seeded in 6-well plates (85000 cells / well) in EGM-2 Medium. Cells reached 100% confluence after 3 days, when the medium was replaced with fresh EGM-2 Medium. 2 days later (= day 0), adipocyte differentiation was induced by changing the medium to adipocyte differentiation (AD) medium (DMEM/Ham’s F12 (50:50), 5 mM HEPES, 2 mM L-glutamine, 100 μg/ml normocin, 860 nM insulin, 10 μg/ml apo-transferrin, 100 nM rosiglitazone, 0.2 nM triiodothyronin) supplemented with 100 μM 3-isobutyl-1-methylxanthine (IBMX), and 1 μM dexamethasone (Dex). Medium was replaced at days 2 and 5 with AD medium. At day 7, another medium change was performed using AD medium without insulin. Adipogenesis was monitored and confirmed based on the appearance of lipid droplets by light microscopy.

Acute norepinephrine (NE) stimulation

Based on literature, we chose the most widely used NE concentration and the time point which demonstrated an acute reaction to NE [98]. At day 9 of differentiation, norepinephrine (NE) stimulation was performed by changing the medium to AD medium without insulin, but supplemented with 25 μg/mL ascorbate and either 1 μM norepinephrine (NE; dissolved in 10 mM HCl) or vehicle (VE, 10 mM HCl).

RNA extraction and sequencing

Cells were harvested with TRIzol reagent (Invitrogen) 3 h after stimulation with NE, and RNA isolation was performed according to the manufacturer’s protocol. Quality of RNA was assayed by an RNA Nano chip using the BioAnalyzer 2100 (Agilent); all samples had RIN values ≥8.5. Four μg total RNA per sample were used for the TruSeq Stranded mRNA LT Sample Prep Kit (Illumina) to generate cDNA libraries according to the manufacturer’s protocol. Single read sequencing was carried out using Illumina/Solexa HiSeq 2000. High-throughput sequencing was conducted by the Biomedical Sequencing Facility (BSF) at CeMM in Vienna.

RNA-seq alignment and differential expression analysis

Raw RNA sequencing reads were aligned against the human hg38 genome using STAR aligner with default parameters [99]. The mapped reads were assigned to genes using featureCount from the bioconductor package Rsubread [100]. All annotated genes were quantified across each condition, corresponding to the University of California at Santa Cruz (UCSC) GRCh38.84 annotation. Normalization and differential expression analysis were performed using the R/Bioconductor package DESeq2 [101]. The effect size / differential expression threshold was situated in an absolute log2 fold change of at least 0.19 (1.17 FC) and an adjusted p-value (padj) < 0.01.

Principal component analysis (PCA) and hierarchical clustering

Principal component analysis (PCA) was performed using the procedure implemented in the R function prcomp in the normalized RNA-seq counts. Unsupervised hierarchical clustering of the RNA-seq data was performed with a centered Pearson correlation coefficient algorithm and a complete linkage method using the R function Heatmap.

Pathway analysis

Pathway enrichment analysis was performed using the R/Bioconductor package ReactomePA, version 1.22.0 [102]. P-values were adjusted for multiple comparisons using the Benjamini-Hochberg procedure.

Network reconstruction

Network reconstruction was performed in Cytoscape with the BisoGenet plugin [36, 103], using DE genes as bait nodes. Edges from experimentally validated interactions were added using the following parameters: Organism > Homo sapiens, gene identifiers only; Data Settings > protein-protein interactions; DIP, BIOGRID, HPRD, INTACT, MINT and BIND databases and all experimental methods; Protein-DNA interaction from experimentally validated interactions; BIND and ENCODE databases, all experimental methods; Method > Criteria to build the network, connecting input nodes with the “By adding edges” option, and as Output > Genes.

Betweenness analysis

The Cytoscape plugin Cyto-Hubba [37] was used to analyze the network topologies and to calculate betweenness centrality. DEGs receiving the highest scores in betweenness were selected as the highest-ranked genes in the network and visualized in their network context using atlas force layout.

Transcription factor and regulatory network analysis

We performed a gene based motif enrichment analysis of our network by using the Cytoscape iRegulon plugin [45]. The following parameters were used: Species> Homo sapiens; Search space> gene-based; Motif collection> 10 k; Track collection> 1120 ChIP-Seq tracks; Putative regulatory region> 20 kb centered around TSS; Motif ranking database> 20 kb centered around TSS (10 species); Track ranking database> 20 kb centered around TSS (ChIP-Seq derived); Region-based parameters>default; Recovery parameters: Enrichment score threshold> 3.0; ROC threshold for AUC calculation > 0.03; Rank threshold > 5000; Transcription factor prediction: Minimum identity between orthologous genes > 0.0 and Maximum FDR on motif similarity > 0.001. Circos plots of predicted transcriptional networks were created using Circa software (OMGGenomics, 2017).

Functional network construction

We constructed a network derived from a Functional Enrichment Analysis (FEA) using DAVID Functional Annotation Clustering (DAVID-FAC), clustering genes in groups of highly related terms [104]. Visualization of the functional gene network and metagroup inference were performed by the R/Bioconductor package FGNet [50].

Statistical analysis

Individual changes in immediate-early gene expression after NE stimulation were determined using RNA-seq read counts normalized by counts per million (CPM). Comparisons among groups were performed using a Kruskall-Wallis test followed by a post-hoc Dunn’s multiple comparison test. All the statistical analyses and graphics were performed using R, version 3.3.3 (2017-03-06).

Abbreviations

cAMP:

Cyclic adenosine monophosphate

CPM:

Counts per million

DEG:

Differentially expressed genes

hpASCs:

Human primary adipose-derived stem cells

IEGs:

Immediate-early genes

log2FC:

Logarithm 2-fold change

LRs:

Local regulators

MRs:

Master regulators

NE:

Norepinephrine

NES:

Normalized enrichment score

NGF:

Nerve growth factor

PCA:

Principal component analysis

PDGF:

Platelet-derived growth factor

PKA:

Protein kinase

PPI:

Protein-protein interaction

TFs:

Transcription factors

TRN:

Transcriptional regulatory network

WAT:

White adipose tissue

References

  1. 1.

    Trayhurn P, Beattie JH. Physiological role of adipose tissue: white adipose tissue as an endocrine and secretory organ. Proc Nutr Soc. 2001;60(3):329–39.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Mueller E. Browning and Graying: novel transcriptional regulators of Brown and Beige fat tissues and aging. Front Endocrinol (Lausanne). 2016;7:19.

    Google Scholar 

  3. 3.

    Nedergaard J, Bengtsson T, Cannon B. Unexpected evidence for active brown adipose tissue in adult humans. Am J Physiol Endocrinol Metab. 2007;293(2):E444–52.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Cannon B, Nedergaard J. Brown adipose tissue: function and physiological significance. Physiol Rev. 2004;84(1):277–359.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Pisani DF, Djedaini M, Beranger GE, Elabd C, Scheideler M, Ailhaud G, Amri EZ. Differentiation of human adipose-derived stem cells into "Brite" (Brown-in-white) adipocytes. Front Endocrinol (Lausanne). 2011;2:87.

    CAS  Google Scholar 

  6. 6.

    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.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Ma X, Xu L, Alberobello AT, Gavrilova O, Bagattin A, Skarulis M, Liu J, Finkel T, Mueller E. Celastrol protects against obesity and metabolic dysfunction through activation of a HSF1-PGC1alpha transcriptional Axis. Cell Metab. 2015;22(4):695–708.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Rosenwald M, Perdikari A, Rulicke T, Wolfrum C. Bi-directional interconversion of brite and white adipocytes. Nat Cell Biol. 2013;15(6):659–67.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Tank AW, Lee Wong D. Peripheral and central effects of circulating catecholamines. Compr Physiol. 2015;5(1):1–15.

    PubMed  Google Scholar 

  10. 10.

    Chen-Izu Y, Xiao RP, Izu LT, Cheng H, Kuschel M, Spurgeon H, Lakatta EG. G(i)-dependent localization of beta(2)-adrenergic receptor signaling to L-type ca(2+) channels. Biophys J. 2000;79(5):2547–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Fredriksson JM, Lindquist JM, Bronnikov GE, Nedergaard J. Norepinephrine induces vascular endothelial growth factor gene expression in brown adipocytes through a beta -adrenoreceptor/cAMP/protein kinase a pathway involving Src but independently of Erk1/2. J Biol Chem. 2000;275(18):13802–11.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Kuo A, Lee MY, Yang K, Gross RW, Sessa WC. Caveolin-1 regulates lipid droplet metabolism in endothelial cells via autocrine prostacyclin stimulated cAMP-mediated lipolysis. J Biol Chem. 2017;293(3):973-83.

    PubMed  Article  Google Scholar 

  13. 13.

    Robidoux J, Cao W, Quan H, Daniel KW, Moukdar F, Bai X, Floering LM, Collins S. Selective activation of mitogen-activated protein (MAP) kinase kinase 3 and p38alpha MAP kinase is essential for cyclic AMP-dependent UCP1 expression in adipocytes. Mol Cell Biol. 2005;25(13):5466–79.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Dong M, Lin J, Lim W, Jin W, Lee HJ. Role of brown adipose tissue in metabolic syndrome, aging, and cancer cachexia. Front Med. 2017;12(2):130-8.

    PubMed  Article  Google Scholar 

  15. 15.

    Loh RKC, Kingwell BA, Carey AL. Human brown adipose tissue as a target for obesity management; beyond cold-induced thermogenesis. Obes Rev. 2017;18(11):1227–42.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Au-Yong IT, Thorn N, Ganatra R, Perkins AC, Symonds ME. Brown adipose tissue and seasonal variation in humans. Diabetes. 2009;58(11):2583–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Collins S, Cao W, Robidoux J. Learning new tricks from old dogs: beta-adrenergic receptors teach new lessons on firing up adipose tissue metabolism. Mol Endocrinol. 2004;18(9):2123–31.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Loft A, Forss I, Siersbaek MS, Schmidt SF, Larsen AS, Madsen JG, Pisani DF, Nielsen R, Aagaard MM, Mathison A, et al. Browning of human adipocytes requires KLF11 and reprogramming of PPARgamma superenhancers. Genes Dev. 2015;29(1):7–22.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  19. 19.

    Loft A, Forss I, Mandrup S. Genome-wide insights into the development and function of thermogenic adipocytes. Trends Endocrinol Metab. 2017;28(2):104–20.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Herschman HR. Primary response genes induced by growth factors and tumor promoters. Annu Rev Biochem. 1991;60:281–319.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Bahrami S, Drablos F. Gene regulation in the immediate-early response process. Adv Biol Regul. 2016;62:37–49.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Fowler T, Sen R, Roy AL. Regulation of primary response genes. Mol Cell. 2011;44(3):348–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Jiang M, Gao Z, Wang J, Nurminsky DI. Evidence for a hierarchical transcriptional circuit in Drosophila male germline involving testis-specific TAF and two gene-specific transcription factors, mod and Acj6. FEBS Lett. 2017;12(2):130-8.

  24. 24.

    Natarajan M, Lin KM, Hsueh RC, Sternweis PC, Ranganathan R. A global analysis of cross-talk in a mammalian cellular signalling network. Nat Cell Biol. 2006;8(6):571–80.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Ledezma-Tejeida D, Ishida C, Collado-Vides J. Genome-wide mapping of transcriptional regulation and metabolism describes information-processing units in Escherichia coli. Front Microbiol. 2017;8:1466.

    PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Amit I, Garber M, Chevrier N, Leite AP, Donner Y, Eisenhaure T, Guttman M, Grenier JK, Li W, Zuk O, et al. Unbiased reconstruction of a mammalian transcriptional network mediating pathogen responses. Science. 2009;326(5950):257–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Barabasi AL, Oltvai ZN. Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004;5(2):101–13.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Aittokallio T, Schwikowski B. Graph-based methods for analysing networks in cell biology. Brief Bioinform. 2006;7(3):243–55.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Joy MP, Brock A, Ingber DE, Huang S. High-betweenness proteins in the yeast protein interaction network. J Biomed Biotechnol. 2005;2005(2):96–103.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. 30.

    Higareda-Almaraz JC, Enriquez-Gasca Mdel R, Hernandez-Ortiz M, Resendis-Antonio O, Encarnacion-Guevara S. Proteomic patterns of cervical cancer cell lines, a network perspective. BMC Syst Biol. 2011;5:96.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Rhee DY, Cho DY, Zhai B, Slattery M, Ma L, Mintseris J, Wong CY, White KP, Celniker SE, Przytycka TM, et al. Transcription factor networks in Drosophila melanogaster. Cell Rep. 2014;8(6):2031–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Katsogiannou M, Andrieu C, Baylot V, Baudot A, Dusetti NJ, Gayet O, Finetti P, Garrido C, Birnbaum D, Bertucci F, et al. The functional landscape of Hsp27 reveals new cellular processes such as DNA repair and alternative splicing and proposes novel anticancer targets. Mol Cell Proteomics. 2014;13(12):3585–601.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Bartolini D, Galli F. The functional interactome of GSTP: a regulatory biomolecular network at the interface with the Nrf2 adaption response to oxidative stress. J Chromatogr B Analyt Technol Biomed Life Sci. 2016;1019:29–44.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Stehle JH, Foulkes NS, Molina CA, Simonneaux V, Pevet P, Sassone-Corsi P. Adrenergic signals direct rhythmic expression of transcriptional repressor CREM in the pineal gland. Nature. 1993;365(6444):314–20.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Wan Y. Bone marrow mesenchymal stem cells: fat on and blast off by FGF21. Int J Biochem Cell Biol. 2013;45(3):546–9.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11.

    PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Melko M, Nguyen LS, Shaw M, Jolly L, Bardoni B, Gecz J. Loss of FMR2 further emphasizes the link between deregulation of immediate early response genes FOS and JUN and intellectual disability. Hum Mol Genet. 2013;22(15):2984–91.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Sarjeant K, Stephens JM. Adipogenesis. Cold Spring Harb Perspect Biol. 2012;4(9):a008417.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Iwaki K, Sukhatme VP, Shubeita HE, Chien KR. Alpha- and beta-adrenergic stimulation induces distinct patterns of immediate early gene expression in neonatal rat myocardial cells. Fos/Jun expression is associated with sarcomere assembly; Egr-1 induction is primarily an alpha 1-mediated response. J Biol Chem. 1990;265(23):13809–17.

    CAS  PubMed  Google Scholar 

  41. 41.

    Guo W, Flanagan J, Jasuja R, Kirkland J, Jiang L, Bhasin S. The effects of myostatin on adipogenic differentiation of human bone marrow-derived mesenchymal stem cells are mediated through cross-communication between Smad3 and Wnt/beta-catenin signaling pathways. J Biol Chem. 2008;283(14):9136–45.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Zhang L, Su P, Xu C, Chen C, Liang A, Du K, Peng Y, Huang D. Melatonin inhibits adipogenesis and enhances osteogenesis of human mesenchymal stem cells by suppressing PPARgamma expression and enhancing Runx2 expression. J Pineal Res. 2010;49(4):364–72.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Chan SS, Kyba M. What is a master regulator? J Stem Cell Res Ther. 2013;3:114.

  44. 44.

    Martinez-Antonio A, Collado-Vides J. Identifying global regulators in transcriptional regulatory networks in bacteria. Curr Opin Microbiol. 2003;6(5):482–9.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Janky R, Verfaillie A, Imrichova H, Van de Sande B, Standaert L, Christiaens V, Hulselmans G, Herten K, Naval Sanchez M, Potier D, et al. iRegulon: from a gene list to a gene regulatory network using large motif and track collections. PLoS Comput Biol. 2014;10(7):e1003731.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  46. 46.

    Kim SW, Choi JH, Mukherjee R, Hwang KC, Yun JW. Proteomic identification of fat-browning markers in cultured white adipocytes treated with curcumin. Mol Cell Biochem. 2016;415(1–2):51–66.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    Pardo R, Blasco N, Vila M, Beiroa D, Nogueiras R, Canas X, Simo R, Sanchis D, Villena JA. EndoG knockout mice show increased Brown adipocyte recruitment in white adipose tissue and improved glucose homeostasis. Endocrinology. 2016;157(10):3873–87.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Green M, Schuetz TJ, Sullivan EK, Kingston RE. A heat shock-responsive domain of human HSF1 that regulates transcription activation domain function. Mol Cell Biol. 1995;15(6):3354–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Ueda HR, Hayashi S, Chen W, Sano M, Machida M, Shigeyoshi Y, Iino M, Hashimoto S. System-level identification of transcriptional circuits underlying mammalian circadian clocks. Nat Genet. 2005;37(2):187–92.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Aibar S, Fontanillo C, Droste C, De Las Rivas J. Functional gene networks: R/bioc package to generate and analyse gene networks derived from functional enrichment and clustering. Bioinformatics. 2015;31(10):1686–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Chernogubova E, Cannon B, Bengtsson T. Norepinephrine increases glucose transport in brown adipocytes via beta3-adrenoceptors through a cAMP, PKA, and PI3-kinase-dependent pathway stimulating conventional and novel PKCs. Endocrinology. 2004;145(1):269–80.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Cinti S. Transdifferentiation properties of adipocytes in the adipose organ. Am J Physiol Endocrinol Metab. 2009;297(5):E977–86.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Barquissau V, Beuzelin D, Pisani DF, Beranger GE, Mairal A, Montagner A, Roussel B, Tavernier G, Marques MA, Moro C, et al. White-to-brite conversion in human adipocytes promotes metabolic reprogramming towards fatty acid anabolic and catabolic pathways. Mol Metab. 2016;5(5):352–65.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    Okamatsu-Ogura Y, Fukano K, Tsubota A, Nio-Kobayashi J, Nakamura K, Morimatsu M, Sakaue H, Saito M, Kimura K. Cell-cycle arrest in mature adipocytes impairs BAT development but not WAT browning, and reduces adaptive thermogenesis in mice. Sci Rep. 2017;7(1):6648.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  55. 55.

    Lim S, Honek J, Xue Y, Seki T, Cao Z, Andersson P, Yang X, Hosaka K, Cao Y. Cold-induced activation of brown adipose tissue and adipose angiogenesis in mice. Nat Protoc. 2012;7(3):606–15.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Sidossis LS, Porter C, Saraf MK, Borsheim E, Radhakrishnan RS, Chao T, Ali A, Chondronikola M, Mlcak R, Finnerty CC, et al. Browning of subcutaneous white adipose tissue in humans after severe adrenergic stress. Cell Metab. 2015;22(2):219–27.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Feldmann HM, Golozoubova V, Cannon B, Nedergaard J. UCP1 ablation induces obesity and abolishes diet-induced thermogenesis in mice exempt from thermal stress by living at thermoneutrality. Cell Metab. 2009;9(2):203–9.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Wu J, Bostrom P, Sparks LM, Ye L, Choi JH, Giang AH, Khandekar M, Virtanen KA, Nuutila P, Schaart G, et al. Beige adipocytes are a distinct type of thermogenic fat cell in mouse and human. Cell. 2012;150(2):366–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Merlin J, Sato M, Chia LY, Fahey R, Pakzad M, Nowell CJ, Summers RJ, Bengtsson T, Evans BA, Hutchinson DS. Rosiglitazone and a beta3-adrenoceptor agonist are both required for functional Browning of white adipocytes in culture. Front Endocrinol (Lausanne). 2018;9:249.

    Article  Google Scholar 

  60. 60.

    Shinoda K, Luijten IH, Hasegawa Y, Hong H, Sonne SB, Kim M, Xue R, Chondronikola M, Cypess AM, Tseng YH, et al. Genetic and functional characterization of clonally derived adult human brown adipocytes. Nat Med. 2015;21(4):389–94.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Nakagami H. The mechanism of white and brown adipocyte differentiation. Diabetes Metab J. 2013;37(2):85–90.

    PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Baboota RK, Singh DP, Sarma SM, Kaur J, Sandhir R, Boparai RK, Kondepudi KK, Bishnoi M. Capsaicin induces "brite" phenotype in differentiating 3T3-L1 preadipocytes. PLoS One. 2014;9(7):e103093.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  63. 63.

    Tsukamoto S, Mizuta T, Fujimoto M, Ohte S, Osawa K, Miyamoto A, Yoneyama K, Murata E, Machiya A, Jimi E, et al. Smad9 is a new type of transcriptional regulator in bone morphogenetic protein signaling. Sci Rep. 2014;4:7596.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Vanhatupa S, Ojansivu M, Autio R, Juntunen M, Miettinen S. Bone morphogenetic Protein-2 induces donor-dependent osteogenic and Adipogenic differentiation in human adipose stem cells. Stem Cells Transl Med. 2015;4(12):1391–402.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Malemud CJ. The PI3K/Akt/PTEN/mTOR pathway: a fruitful target for inducing cell death in rheumatoid arthritis? Future Med Chem. 2015;7(9):1137–47.

    CAS  PubMed  Article  Google Scholar 

  66. 66.

    Noskovicova N, Petrek M, Eickelberg O, Heinzelmann K. Platelet-derived growth factor signaling in the lung. From lung development and disease to clinical studies. Am J Respir Cell Mol Biol. 2015;52(3):263–84.

    CAS  PubMed  Article  Google Scholar 

  67. 67.

    Lattanzi W, Geloso MC. Editorial: crosstalk between the osteogenic and neurogenic stem cell niches: how far are they from each other? Front Cell Neurosci. 2015;9:504.

    PubMed  Article  Google Scholar 

  68. 68.

    Sharma AD, Wiederin J, Uz M, Ciborowski P, Mallapragada SK, Gendelman HE, Sakaguchi DS. Proteomic analysis of mesenchymal to Schwann cell transdifferentiation. J Proteome. 2017;165:93–101.

    CAS  Article  Google Scholar 

  69. 69.

    Feng Z. p53 regulation of the IGF-1/AKT/mTOR pathways and the endosomal compartment. Cold Spring Harb Perspect Biol. 2010;2(2):a001057.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  70. 70.

    Feng Z, Zhang H, Levine AJ, Jin S. The coordinate regulation of the p53 and mTOR pathways in cells. Proc Natl Acad Sci U S A. 2005;102(23):8204–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Gonzales KA, Liang H, Lim YS, Chan YS, Yeo JC, Tan CP, Gao B, Le B, Tan ZY, Low KY, et al. Deterministic restriction on pluripotent state dissolution by cell-cycle pathways. Cell. 2015;162(3):564–79.

    CAS  PubMed  Article  Google Scholar 

  72. 72.

    Califano A, Butte AJ, Friend S, Ideker T, Schadt E. Leveraging models of cell regulation and GWAS data in integrative network-based association studies. Nat Genet. 2012;44(8):841–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

    Furlong LI. Human diseases through the lens of network biology. Trends Genet. 2013;29(3):150–9.

    CAS  PubMed  Article  Google Scholar 

  74. 74.

    Carter H, Hofree M, Ideker T. Genotype to phenotype via network analysis. Curr Opin Genet Dev. 2013;23(6):611–21.

    CAS  PubMed  Article  Google Scholar 

  75. 75.

    Higareda-Almaraz JC, Valtierra-Gutierrez IA, Hernandez-Ortiz M, Contreras S, Hernandez E, Encarnacion-Guevara S. Analysis and prediction of pathways in HeLa cells by integrating biological levels of organization with systems-biology approaches. PLoS One. 2013;8(6):e65433.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    Xiao Q, Wang J, Peng X, Wu FX, Pan Y. Identifying essential proteins from active PPI networks constructed with dynamic gene expression. BMC Genomics. 2015;16(Suppl 3):S1.

    PubMed  PubMed Central  Article  Google Scholar 

  77. 77.

    Mistry D, Wise RP, Dickerson JA. DiffSLC: a graph centrality method to detect essential proteins of a protein-protein interaction network. PLoS One. 2017;12(11):e0187091.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  78. 78.

    Yu H, Kim PM, Sprecher E, Trifonov V, Gerstein M. The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput Biol. 2007;3(4):e59.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  79. 79.

    Jung H, Han S, Kim S. The construction of regulatory network for insulin-mediated genes by integrating methods based on transcription factor binding motifs and gene expression variations. Genomics Inform. 2015;13(3):76–80.

    PubMed  PubMed Central  Article  Google Scholar 

  80. 80.

    Heroux NA, Osborne BF, Miller LA, Kawan M, Buban KN, Rosen JB, Stanton ME. Differential expression of the immediate early genes c-Fos, arc, Egr-1, and Npas4 during long-term memory formation in the context preexposure facilitation effect (CPFE). Neurobiol Learn Mem. 2017;147:128-38.

    CAS  PubMed  Article  Google Scholar 

  81. 81.

    Ugajin A, Uchiyama H, Miyata T, Sasaki T, Yajima S, Ono M. Identification and initial characterization of novel neural immediate early genes possibly differentially contributing to foraging-related learning and memory processes in the honeybee. Insect Mol Biol. 2017;27(2):154-65.

    PubMed  Article  CAS  Google Scholar 

  82. 82.

    Muniz JA, Prieto JP, Gonzalez B, Sosa MH, Cadet JL, Scorza C, Urbano FJ, Bisagno V. Cocaine and caffeine effects on the conditioned place preference test: concomitant changes on early genes within the mouse prefrontal cortex and nucleus Accumbens. Front Behav Neurosci. 2017;11:200.

    PubMed  PubMed Central  Article  Google Scholar 

  83. 83.

    O'Donnell A, Odrowaz Z, Sharrocks AD. Immediate-early gene activation by the MAPK pathways: what do and don't we know? Biochem Soc Trans. 2012;40(1):58–66.

    CAS  PubMed  Article  Google Scholar 

  84. 84.

    Obier N, Cauchy P, Assi SA, Gilmour J, Lie ALM, Lichtinger M, Hoogenkamp M, Noailles L, Cockerill PN, Lacaud G, et al. Cooperative binding of AP-1 and TEAD4 modulates the balance between vascular smooth muscle and hemogenic cell fate. Development. 2016;143(23):4324–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Ma X, Lee P, Chisholm DJ, James DE. Control of adipocyte differentiation in different fat depots; implications for pathophysiology or therapy. Front Endocrinol (Lausanne). 2015;6:1.

    Google Scholar 

  86. 86.

    Jin X, Qiao A, Moskophidis D, Mivechi NF. WITHDRAWN: abrogation of heat shock factor 1 (Hsf1) phosphorylation deregulates its activity and lowers activation threshold, leading to obesity in mice. J Biol Chem. July 2017..

  87. 87.

    Qiao A, Jin X, Pang J, Moskophidis D, Mivechi NF. The transcriptional regulator of the chaperone response HSF1 controls hepatic bioenergetics and protein homeostasis. J Cell Biol. 2017;216(3):723–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. 88.

    Ohno T, Onishi Y, Ishida N. A novel E4BP4 element drives circadian expression of mPeriod2. Nucleic Acids Res. 2007;35(2):648–55.

    CAS  PubMed  Article  Google Scholar 

  89. 89.

    Junghans D, Chauvet S, Buhler E, Dudley K, Sykes T, Henderson CE. The CES-2-related transcription factor E4BP4 is an intrinsic regulator of motoneuron growth and survival. Development. 2004;131(18):4425–34.

    CAS  PubMed  Article  Google Scholar 

  90. 90.

    Hirai T, Tanaka K, Togari A. alpha1-adrenergic receptor signaling in osteoblasts regulates clock genes and bone morphogenetic protein 4 expression through up-regulation of the transcriptional factor nuclear factor IL-3 (Nfil3)/E4 promoter-binding protein 4 (E4BP4). J Biol Chem. 2014;289(24):17174–83.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  91. 91.

    Morris DP, Lei B, Longo LD, Bomsztyk K, Schwinn DA, Michelotti GA. Temporal dissection of rate limiting transcriptional events using pol II ChIP and RNA analysis of adrenergic stress gene activation. PLoS One. 2015;10(8):e0134442.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  92. 92.

    Rothstein TL, Guo B. Receptor crosstalk: reprogramming B cell receptor signalling to an alternate pathway results in expression and secretion of the autoimmunity-associated cytokine, osteopontin. J Intern Med. 2009;265(6):632–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    Zhao A, Yang L, Ma K, Sun M, Li L, Huang J, Li Y, Zhang C, Li H, Fu X. Overexpression of cyclin D1 induces the reprogramming of differentiated epidermal cells into stem cell-like cells. Cell Cycle. 2016;15(5):644–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  94. 94.

    Solloway MJ, Madjidi A, Gu C, Eastham-Anderson J, Clarke HJ, Kljavin N, Zavala-Solorio J, Kates L, Friedman B, Brauer M, et al. Glucagon couples hepatic amino acid catabolism to mTOR-dependent regulation of alpha-cell mass. Cell Rep. 2015;12(3):495–510.

    CAS  PubMed  Article  Google Scholar 

  95. 95.

    Neal CL, Xu J, Li P, Mori S, Yang J, Neal NN, Zhou X, Wyszomierski SL, Yu D. Overexpression of 14-3-3zeta in cancer cells activates PI3K via binding the p85 regulatory subunit. Oncogene. 2012;31(7):897–906.

    CAS  PubMed  Article  Google Scholar 

  96. 96.

    Mohammad DK, Nore BF, Hussain A, Gustafsson MO, Mohamed AJ, Smith CI. Dual phosphorylation of Btk by Akt/protein kinase b provides docking for 14-3-3zeta, regulates shuttling, and attenuates both tonic and induced signaling in B cells. Mol Cell Biol. 2013;33(16):3214–26.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  97. 97.

    Karbiener M, Pisani DF, Frontini A, Oberreiter LM, Lang E, Vegiopoulos A, Mossenbock K, Bernhardt GA, Mayr T, Hildner F, et al. MicroRNA-26 family is required for human adipogenesis and drives characteristics of brown adipocytes. Stem Cells. 2014;32(6):1578–90.

    CAS  PubMed  Article  Google Scholar 

  98. 98.

    Herron D, Rehnmark S, Nechad M, Loncar D, Cannon B, Nedergaard J. Norepinephrine-induced synthesis of the uncoupling protein thermogenin (UCP) and its mitochondrial targeting in brown adipocytes differentiated in culture. FEBS Lett. 1990;268(1):296–300.

    CAS  PubMed  Article  Google Scholar 

  99. 99.

    Dobin A, Gingeras TR. Mapping RNA-seq Reads with STAR. Curr Protoc Bioinformatics. 2015;51:11 14 11–19.

    PubMed  PubMed Central  Google Scholar 

  100. 100.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.

    CAS  PubMed  Article  Google Scholar 

  101. 101.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  102. 102.

    Yu G, He QY. ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualization. Mol BioSyst. 2016;12(2):477–9.

    CAS  PubMed  Article  Google Scholar 

  103. 103.

    Martin A, Ochagavia ME, Rabasa LC, Miranda J, Fernandez-de-Cossio J, Bringas R. BisoGenet: a new tool for gene network building, visualization and analysis. BMC Bioinformatics. 2010;11:91.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  104. 104.

    Huang da W, Sherman BT, Zheng X, Yang J, Imamichi T, Stephens R, Lempicki RA. Extracting biological meaning from large gene lists with DAVID. Curr Protoc Bioinformatics. 2009, Chapter 13:Unit;13:11.

    PubMed  Google Scholar 

Download references

Acknowledgements

We thank Florian Hildner and Heinz Redl, Austrian Cluster for Tissue Regeneration, Linz, Austria, for providing primary preadipocytes. We are grateful to Ilse A. Valtierra-Gutierrez, Rocío Enríquez-Gasca, Rocío Dominguez-Vidaña, for proofreading, observations and suggestions; to Alejandra Medina-Rivera and Daniela Ledezma-Tejeida for proofreading; and to Anne Loft, Anastasia Georgiadi, and Frits Mattijssen for helpful discussions.

Funding

This work was funded by the German Centre for Diabetes Research (DZD) and the Austrian Science Fund (FWF, P25729-B19). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

Processed RNA-seq data are included as Additional file 3: Table S4. RNA-seq raw data have been deposited in the ArrayExpress database at EMBL-EBI (http://www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-6345.

Author information

Affiliations

  1. Institute for Diabetes and Cancer (IDC), Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany

    Juan Carlos Higareda-Almaraz, Maude Giroud, Stephan Herzig & Marcel Scheideler

  2. Joint Heidelberg-IDC Translational Diabetes Program, Heidelberg University Hospital, Heidelberg, Germany

    Juan Carlos Higareda-Almaraz, Maude Giroud, Stephan Herzig & Marcel Scheideler

  3. Molecular Metabolic Control, Medical Faculty, Technical University, Munich, Germany

    Juan Carlos Higareda-Almaraz, Maude Giroud, Stephan Herzig & Marcel Scheideler

  4. German Center for Diabetes Research (DZD), Neuherberg, Germany

    Juan Carlos Higareda-Almaraz, Maude Giroud, Stephan Herzig & Marcel Scheideler

  5. Department of Phoniatrics, ENT University Hospital, Medical University of Graz, Graz, Austria

    Michael Karbiener

  6. CeMM, Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria

    Florian M. Pauler

  7. Present Address: Institute of Science and Technology (IST) Austria, Klosterneuburg, Austria

    Florian M. Pauler & Teresa Gerhalter

  8. NMR laboratory, Institute of Myology, Hopital Universitaire Pitie Salpetriere, Paris, France

    Juan Carlos Higareda-Almaraz & Marcel Scheideler

Authors
  1. Juan Carlos Higareda-Almaraz
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Michael Karbiener
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Maude Giroud
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Florian M. Pauler
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Teresa Gerhalter
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Stephan Herzig
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Marcel Scheideler
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JCHA, MK, SH and MS made the conception and design of the study. MG, FMP and TG made substantial contributions to the acquisition of data, while JCHA made a major contribution to the analysis and interpretation of data. All authors have sufficiently participated in this work, were involved in drafting the manuscript or revising it critically for important intellectual content, gave final approval of the version to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Corresponding authors

Correspondence to Juan Carlos Higareda-Almaraz or Marcel Scheideler.

Ethics declarations

Ethics approval and consent to participate

Donor samples and all experiments were performed in accordance with protocols approved by the Austrian Cluster for Tissue Regeneration, Linz, Austria. Samples were collected from each patient after obtaining their written informed consent. This study was performed according to the Declaration of Helsinki and approved by the local ethics commission (Ethikkommission des Landes Oberoesterreich). The analysis was conducted anonymously, without public disclosure of patient data.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Genes with highest betweenness. Table S2. Local and Master regulators up-regulated. Table S3. Local and Master regulators down-regulated (XLSX 23 kb).

Additional file 2:

Figure S1. Adipocyte and macrophage marker gene expression in white adipocytes in the unstimulated and 3 h NE-stimulated state. Y-axis indicates the transcript abundance in Transcript Per Million (TPM) (PDF 41 kb).

Additional file 3:

Table S4. Counts per Gene per Million Reads Mapped. (XLSX 2751 kb).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Higareda-Almaraz, J.C., Karbiener, M., Giroud, M. et al. Norepinephrine triggers an immediate-early regulatory network response in primary human white adipocytes. BMC Genomics 19, 794 (2018). https://doi.org/10.1186/s12864-018-5173-0

Download citation

Keywords

  • Norepinephrine stimulation
  • White adipocyte
  • Early cell fate
  • Network biology
  • Immediate-early gene
  • Transcriptional regulatory network

BMC Genomics

ISSN: 1471-2164

Contact us