Brown adipose tissue (BAT) dissipates chemical energy and generates heat to protect animals from cold and obesity. Rodents possess two types of UCP-1 positive brown adipocytes arising from distinct developmental lineages: "classical" brown adipocytes develop during the prenatal stage whereas "beige" or "brite" cells that reside in white adipose tissue (WAT) develop during the postnatal stage in response to chronic cold or PPARγ agonists. Beige cells' inducible characteristics make them a promising therapeutic target for obesity treatment, however, the relevance of this cell type in humans remains unknown. In the present study, we determined the gene signatures that were unique to classical brown adipocytes and to beige cells induced by a specific PPARγ agonist rosiglitazone in mice. Subsequently we applied the transcriptional data to humans and examined the molecular signatures of human BAT isolated from multiple adipose depots. To our surprise, nearly all the human BAT abundantly expressed beige cell-selective genes, but the expression of classical brown fat-selective genes were nearly undetectable. Interestingly, expression of known brown fat-selective genes such as PRDM16 was strongly correlated with that of the newly identified beige cell-selective genes, but not with that of classical brown fat-selective genes. Furthermore, histological analyses showed that a new beige cell marker, CITED1, was selectively expressed in the UCP1-positive beige cells as well as in human BAT. These data indicate that human BAT may be primary composed of beige/brite cells.
Brown adipose tissue (BAT) dissipates chemical energy and generates heat to protect animals from cold and obesity. Rodents possess two types of UCP-1 positive brown adipocytes arising from distinct developmental lineages: "classical" brown adipocytes develop during the prenatal stage whereas "beige" or "brite" cells that reside in white adipose tissue (WAT) develop during the postnatal stage in response to chronic cold or PPARγ agonists. Beige cells' inducible characteristics make them a promising therapeutic target for obesity treatment, however, the relevance of this cell type in humans remains unknown. In the present study, we determined the gene signatures that were unique to classical brown adipocytes and to beige cells induced by a specific PPARγ agonist rosiglitazone in mice. Subsequently we applied the transcriptional data to humans and examined the molecular signatures of human BAT isolated from multiple adipose depots. To our surprise, nearly all the human BAT abundantly expressed beige cell-selective genes, but the expression of classical brown fat-selective genes were nearly undetectable. Interestingly, expression of known brown fat-selective genes such as PRDM16 was strongly correlated with that of the newly identified beige cell-selective genes, but not with that of classical brown fat-selective genes. Furthermore, histological analyses showed that a new beige cell marker, CITED1, was selectively expressed in the UCP1-positive beige cells as well as in human BAT. These data indicate that human BAT may be primary composed of beige/brite cells.
Brown adipose tissue (BAT) is specialized to generate heat by dissipating chemical energy as a defense against cold and obesity. Thermogenesis in the BAT is mediated through a brown fat-specific mitochondrial protein, uncoupling protein 1 (UCP1), which plays an important role in the control of energy homeostasis. Indeed, loss of UCP1 causes cold-intolerance as well as obesity at thermoneutality in mice [1], [2]. Importantly, recent studies using non-invasive imaging technologies such as 18fluoro-labeled 2-deoxy-glucose positron emission tomography (18FDG-PET) scanning and MRI scanning clearly demonstrated that adult humans have significant amounts of active brown fat deposits [3], [4], [5], [6], [7], [8], [9]. A histological study also identified brown fat precursors in the BAT from adult humans [10]. Furthermore, the glucose uptake capacity in the BAT correlates inversely with adiposity, indicating that variation in the amount and/or thermogenic activity of BAT may contribute to the propensity for weight gain in humans [11], [12].Current evidence suggests that there are at least two types of brown adipocytes that arise from distinct developmental lineages in rodents. “Classical” brown adipocytes that reside in the interscapular BAT depot and perirenal regions develop during the prenatal stage from Myf-5 positive myoblast precursors [13]. Global gene expression analyses indicate that the brown fat precursors in the interscapular BAT show a gene signature that resembles that of skeletal muscle cells [14]. In addition, a proteome study has shown that mitochondria from the interscapular BAT are most similar to their counterparts in skeletal muscle at transcriptional and protein levels [15]. These Myf-5 derived precursors differentiate into brown adipocytes through the action of two transcriptional factors, PRDM16 and C/EBP-β [13], [16], although it remains unknown whether the Myf-5 positive cells clonally give rise to brown adipocytes and myocytes.Another type of brown adipocytes is sporadically found as a copious cluster in the white adipose of adult animals that have been chronically exposed to cold or to PPARγ ligands (so called beige or brite cells). These “inducible” brown adipocytes possess the morphological and biochemical characteristics of classical brown adipocytes, including the expression of UCP1 and multilocular lipid droplet [17], however, they arise from a non-Myf-5 lineage [13]. A recent study has shown that PDGFRa-positive bi-potent progenitors in the abdominal WAT give rise to brown adipocytes in response to beta-adrenergic stimulation in vivo
[18]. Cannon and Nedergaard’s group has also shown that these brite cells express several well-known brown fat-selective genes such as Ucp1, but do not express myocyte-enriched genes, indicating that these two types of brown adipocytes have distinct expression profiles [19], [20]. Of note, Kozak’s group has implied that signaling pathways to control these beige cells in the WAT are different from classical brown adipocytes in the interscapular BAT depots [21], [22], [23]. Importantly, increases in number of beige cells in WAT are associated well with a protection against diet induced obesity and metabolic diseases [24], [25], [26]. Hence, an emerging interest in the field of brown fat biology is to characterize the molecular and functional signatures of the beige cells in vivo. In particular, the relevance of the beige cells in humans must be now addressed.In this study, we asked if humans possess beige cells. To this end, we first employed microarray analysis to define the unique gene sets that are distinctively expressed in the Myf-5 derived classical brown adipocytes and in the non-Myf-5 beige/brite cells in mice. Subsequently we applied the transcriptional data to humanadipose tissues and determined the molecular signatures of human BAT. Surprisingly, nearly all the human BAT isolated from multiple locations abundantly expressed genes unique to beige cells, correlating very well with several known brown fat-selective genes such as PRDM16 and PGC1α. On the contrary, the classical brown adipocyte-selective genes were nearly undetected in human BAT. Our transcriptional and histological analyses indicate that human BAT possesses a molecular signature of beige cells.
Materials and Methods
Animals
All animal experiments were performed according to procedures approved by UCSF’s Institutional Animal Care and Use Committee for animal care and handling. C57BL/6 mice (Jackson Laboratories) were used for the animal experiments. To induce the browning of white fat in mice, CL316243 at a dose of 1 mg/kg or saline control was injected intraperitoneally for 8 days into 6–8 week-old mice. Inguinal adipose tissue was isolated, fixed in 4% paraformaldehyde, and embedded in paraffin for histological analysis.
Immunohistochemistry
Paraffin-embedded sections were cut to 7 µm sections. They were subjected to citrate-based antigen retrieval using Antigen Retrieval Citra Solution (Biogenex) and incubated for 1 hour at room temperature with rabbit polyclonal antibodies against CITED1 at a dilution of 1∶500 (Abcam) or UCP1 at a dilution of 1∶1000 (Abcam) or IgG as negative control. Secondary detection was performed for 30 minutes at room temperature using anti-rabbitAlexa 488 and Alexa 594 (Invitrogen). Images were captured using a Zeiss Axioskop 2 Widefield microscope.
Human BAT Samples and MRI Imaging
The protocol for the human study was approved by the Institutional Review Board of Children’s Hospital Los Angeles (CHLA), where the imaging and postmortem examinations were performed. Written informed consent was obtained from all parent/legal guardian(s) in the presence of a witness. MRI scanning was performed in the patients to identify the BAT as previously described [9], [27]. Briefly, 13 post-mortem examinations of children ranging from 3-day to 18-years-old were scanned using MRI implemented on a 3 T whole-body human platform (Achieva, R2.6.3, Philips Healthcare, Cleveland, OH). The mDIXON (Philips Healthcare, Cleveland, Ohio) chemical-shift water-fat multi-echo pulse sequence provided by the manufacturer was utilized in this work using an 8- or 16-channel coil. A consistent acquisition voxel size of 1 mm × 1 mm × 2 mm was maintained across all subjects. For the purpose of this study, the fat fraction (FF) of adipose tissue in the chest and the abdomen was calculated using multi-echo chemical-shift water-fat MRI algorithms as previously described [9], [27]. Image analysis and quantification of FF maps were performed with SliceOmatic (Tomovision, Inc., Magog, Canada) segmentation software. Based on the MRI FF, biopsies of adipose tissue were obtained within 24 hours after death by a pediatric pathologist at sites highly suggestive of BAT. Routine histological and immunohistochemical assays and molecular analysis of tissue samples were performed at the Department of Pathology at CHLA and at UCSF.
Cell Culture
Primary stromal vascular (SV) cells were isolated from inguinal WAT or from interscapular BAT of the same C57BL/6 mice using Collagenase D (1.5 u/ml) and Dispase II (2.4 u/ml). To clearly define the molecular signatures of PPARγ ligand-induced beige cells, we used C57BL/6 mice that possess less beige cells in the absence of environmental stimuli. SV cells were plated in collagen coated culture dishes in DMEM/F12 medium (D-glucose 17.51 mM). The confluent cells were subsequently differentiated in DMEM/F12 medium containing 10% FBS, 0.5 mM isobutylmethylxanthine, 125 nM indomethacin, 1 mM dexamethasone, 850 nM insulin and 1 nM T3. To induce beige cell differentiation, a specific PPARγ agonist rosiglitazone at 1 µM were added in the medium [28]. Two days after induction, cells were switched to the maintenance medium containing 10% FBS, 850 nM insulin and 1 nM T3 for additional 4–5 days until cells were fully differentiated to mature adipocytes. R2F primary skin fibroblasts isolated from human newborn foreskin were cultured following the methods described elsewhere [29].
RNA Isolation and Gene Expression Analysis
Total RNA was isolated from tissues using the Trizol reagent (Invitrogen) or RiboZol reagent (AMRESCO) following the manufacturer’s protocol. Quality of RNA from all the samples was checked by spectrophotometer. Samples that did not clear the quality control (O.D.260/280<1.8) were excluded from the analyses. Reverse transcription reactions were performed using IScript cDNA synthesis kit (Bio-Rad). The sequences of primers used in this study can be found in Table S1. Quantitative real-time PCR (qRT-PCR) was performed with SYBR green fluorescent dye using an ABI ViiA™7 PCR machine. Relative mRNA expression was determined by the ΔΔ-Ct method using TATA-binding protein (TBP) as an endogenous control to normalize samples. For a microarray analysis, Affymetrix GeneChip Mouse Genome 430 2.0 array was used according to established methods [30]. The array data were analyzed using the DNA-Chip Analyzer (dChip) software. The statistical significance of differences in gene expression was assessed by unpaired t-test (P<0.05). Microarray data has been deposited in Gene Expression Omnibus (GEO): GSE35011. Patient and sample information for the tissues used in the expression analyses were provided in Table S2.
Heatmap Visualization
The measured mRNA expression were clustered according to the gene group and visualized as a heatmap representation using Multi Experiment Viewer [31]. Expression levels of all genes were visualized in the same green-black-red scheme in order to evaluate relative expression levels in all samples.
Statistical Analyses
Hierarchical clustering was performed with average linkage method using Pearson correlation as distance metric. Integration of correlation coefficient among gene group was done using average of weighted Fisher’s z-score. We transformed each coefficient into z-score and computed average and then inverse transformed. Test of population correlation coefficient was conducted by Z-test. To detect outliers, we used box-and-whisker plot. Samples beyond the whiskers of the plot were considered outlier. All these statistical analysis were performed using the JMP software version 9 (SAS Institute) or R version 2.11.1.
Results
Distinct Molecular Signatures between Myf-5 Derived Brown Adipocytes and Non-Myf-5 Derived Beige Cells
Classical brown adipocytes in interscapular BAT (Myf-5 derived) and “inducible” beige cells in WAT (non-Myf-5 derived) have distinct developmental origins, although both cell types have morphological and biochemical characteristics of brown fat such as the expression of UCP1 [13]. This raises an important question as to how similar the two types of brown adipocytes are at molecular and functional levels. To this end, we first employed microarray analysis to systematically determine the transcriptional signatures unique to each cell type. To avoid contamination from non-adipocytes, the stromal vascular (SV) fraction was isolated from interscapular BAT or inguinal WAT and fully differentiated in vitro. A subset of the SV cells from inguinal WAT were differentiated in the presence of rosiglitazone at 1 µM to induce the development of beige cells [28]. At least 95% of the cells differentiated into mature adipocytes under this culture condition.Subsequently, mRNA expression of genes was compared between classical brown adipocytes, inguinal white adipocytes and rosiglitazone-treated adipocytes (beige cells). As shown in Fig.1A, we identified a total of 496 genes (>2 fold, P<0.05) that were differentially expressed among the three groups. These genes could be clustered into four groups as follows: Group A included 32 genes whose expression was enriched in inguinal white adipocytes; Group B contained 141 genes whose expression was enriched in rosiglitazone-treated adipocytes (beige cells); Group C included 29 genes enriched both in classical brown adipocytes and in beige cells; Group D contained 294 genes enriched in classical brown adipocytes. For instance, as verified by qRT-PCR, well-known brown fat-selective genes, such as Ucp1, Cidea, Cox8b and Pgc1α were expressed at similar levels both in classical brown adipocytes and in rosiglitazone-treated adipocytes (Fig.1B, upper left). Such “common” brown fat-selective genes for the classical brown adipocytes and rosiglitazone-treated adipocytes (beige cells) represent only a small portion of the total (9%, 29/323 genes in classical brown adipocytes or 17%, 29/170 genes in beige cells). In contrast, group B were genes whose expression was enriched only in the rosiglitazone-induced beige cells (listed in Table S3). This includes fibroblast growth factor 21 (Fgf21), carbonic anhydrase 4 (Car4), and Cbp/p300-interacting transactivator with Glu/Asp-rich carboxy-terminal domain 1 (Cited1). Expression of homeo box gene 9 (Hox9a) was enriched in both white adipocytes and beige cells (Fig.1B, upper middle). Group D includes genes that are highly enriched only in classical brown adipocytes (listed in Table S4). This group contained zinc finger protein of the cerebellum 1 (Zic1), LIM homeobox protein 8 (Lhx8), and epithelial stromal interaction 1 (Epstl1) (Fig.2B, upper right). This is consistent with the previous studies by Cannon and Nedergaard’s group [19], [20]. Expression of angiotensinogen (Ang), resistin, WDNM1-like protein (Wdnm1), and Serpina3a, as found in Group A and reported in previous studies [32], [33], was enriched in white adipocytes (Fig.1B, lower left). On the other hand, expression of general adipogenic genes shared in differentiated white adipocytes, beige cells and brown adipocytes, such as Pparg and Adiponectin (Adipoq) were expressed at similar levels, suggesting that adipocyte differentiation was similar in three groups. These data clearly suggest that the two types of brown adipocytes have distinct characteristics at the molecular levels.
Figure 1
Distinct molecular signatures in Myf-5 derived classical brown adipocytes and non-Myf-5 derived beige cells.
(A) Microarray analysis was performed in fully differentiated brown adipocytes from interscapular BAT (classical brown adipocytes), white adipocytes isolated from inguinal WAT and rosiglitazone-treated adipocytes (beige cells). The color scale shows the fold changes in mRNA expression levels of the genes in blue (down-regulation)-white-red (up-regulation) scheme. (B) mRNA levels for the indicated genes were analyzed by qRT-PCR. Relative expression levels in the brown adipocyte were set as 1.0. * P<0.05, **P<0.01. Data are expressed as means ± SD.
Figure 2
Isolation of human BAT from multiple BAT depots.
(A) Axial (top) and coronal (bottom) MRI scanning images illustrating bilateral supraclavicular BAT depots from chemical-shift water-fat MRI in a 4-months-old female. BAT depots are denoted by arrows and subcutaneous WAT depots are denoted by arrowheads. The first column shows the water images. The second column shows the fat images. The third column shows the fat fraction map of water content, which translates to a lower fat fraction in BAT than WAT. The coronal images were obtained anteriorly at the level of the humeral epiphysis. (B) The isolated BAT depots were fixed and stained by H&E. The sample was obtained from the left supraclavicular area of the same 4-month old female shown in the MR images in Fig. 2A. Note that nearly all the isolated BAT samples consisted of brown adipocytes containing multilocular lipid droplets.
Distinct molecular signatures in Myf-5 derived classical brown adipocytes and non-Myf-5 derived beige cells.
(A) Microarray analysis was performed in fully differentiated brown adipocytes from interscapular BAT (classical brown adipocytes), white adipocytes isolated from inguinal WAT and rosiglitazone-treated adipocytes (beige cells). The color scale shows the fold changes in mRNA expression levels of the genes in blue (down-regulation)-white-red (up-regulation) scheme. (B) mRNA levels for the indicated genes were analyzed by qRT-PCR. Relative expression levels in the brown adipocyte were set as 1.0. * P<0.05, **P<0.01. Data are expressed as means ± SD.
Isolation of Human BAT from Multiple Depots
Recent advance in MRI scanning technology allowed us to identify the BAT in multiple adipose depots in humans [9]. To isolate the BAT from multiple locations in given individuals, we first scanned 13 patients by MRI. Fig2A illustrates axial (top panels) and coronal (bottom panels) images of MRI results in a 4-months-old female. As shown in the first column (water image), BAT is hyperintense, or brighter image than WAT, since BAT contains greater intracellular water content than WAT. On the contrary, in the second column (fat image), BAT is hypointense or darker, than WAT, due to its lower intracellular fat content. Collectively, as shown in the third column, the difference in fat and water content translates to a lower fat fraction in BAT than WAT with a fraction of about 30–60%. Of note, nearly all of the subcutaneous compartments were red and had a fat fraction between 80–100%, which was indicative of white adipose tissue. In contrast, the adipose tissues in the supraclavicular and intra-abdominal areas showed fat fraction characteristics of brown fat. The BAT depots were located in the subcutaneous supraclavicular areas, posterior mediastinum, retroperitoneal, intraabdominal and mesenteric depots, with thigh tissue samples for 4 of the 13 patients. In contrast to rodents in which BAT is preferentially located in the interscapular area, this tissue in humans is primarily recognized in the cervical-supraclavicular area, where it is present across a wide range of ages [3], [4], [34], [35]. On occasion, newborns and very young infants present BAT in the interscapular area, but this depot is replaced by white fat shortly after.
Isolation of human BAT from multiple BAT depots.
(A) Axial (top) and coronal (bottom) MRI scanning images illustrating bilateral supraclavicular BAT depots from chemical-shift water-fat MRI in a 4-months-old female. BAT depots are denoted by arrows and subcutaneous WAT depots are denoted by arrowheads. The first column shows the water images. The second column shows the fat images. The third column shows the fat fraction map of water content, which translates to a lower fat fraction in BAT than WAT. The coronal images were obtained anteriorly at the level of the humeral epiphysis. (B) The isolated BAT depots were fixed and stained by H&E. The sample was obtained from the left supraclavicular area of the same 4-month old female shown in the MR images in Fig. 2A. Note that nearly all the isolated BAT samples consisted of brown adipocytes containing multilocular lipid droplets.Subsequently, the identified BAT depots were biopsied for histological analysis. As represented in Fig. 2B and also previously reported (Gilsanz et al. 2012), H&E staining analyses showed that most of these tissues in supraclavicular areas, posterior mediastinum, retroperitoneal, intraabdominal and mesenteric depots consisted of bona-fide brown adipocytes containing multilocular lipid droplets.
Human BAT Expresses a Large Array of Beige Cell-selective Genes
To examine the molecular signatures of human BAT, total RNA was isolated from the human BAT samples isolated from multiple locations, including subcutaneous supraclavicular areas, posterior mediastinum, retroperitoneal, intraabdominal, mesenteric depots and thigh tissues. As controls, total RNA was also isolated from white adipose tissue and from smooth muscle and included in the analysis. Subsequently, we measured the expression of classical brown adipocyte-selective or beige cell-selective genes that were identified in Fig.1. A recent study by Spiegelman’s group has reported that three genes, CD137, TMEM26, and TBX1 were enriched in the beige cells but not in the classical brown adipocytes nor in the white adipocytes [36]. These genes were therefore included in the expression analyses.Fig.3A shows the expression profiles of the common brown fat-selective genes (Group A), classical brown fat-selective genes (Group B), and beige cell-selective genes (Group C). Expression of PPARγ, a master regulator of adipocyte development, was also measured to validate that these samples were indeed adipose. To our surprise, expression of classical brown fat markers was nearly undetectable, whereas beige cell markers were abundantly expressed in early all human BAT. Notably, beige-selective transcription factors (HoxC8, HoxC9 and CITED1) were highly expressed in human BAT. The trend of expression profile was consistent after normalizing the each value with that in non-adipocyte cells such as primary skin fibroblasts (Fig. S1).
Figure 3
Transcriptional analyses of human BAT.
(A) Expression profiles of the common brown fat-selective genes (Group A), classical brown fat-selective genes (Group B), and beige cell-selective genes (Group C) in human BAT from multiple adipose depots. The color scale shows the mRNA levels of the genes in green (low or no-expression)-black-red (high-expression) scheme. All genes are shown in the identical scale. WAT: white adipose tissue. M: smooth muscle. Each alphabet indicates the tissues from the same patient. (B) Correlation matrix of beige-selective and classical brown-selective gene expression. The color scale shows the Pearson correlation between each gene mRNA expression in blue (negative correlation, −1.0)-gray (no correlation, 0)-red (positive correlation, 1.0) scheme. (C) Strength of correlation with representative group A gene (PRDM16 and PGC1a) and classical-brown-selective genes or beige cell-selective genes. **P<0.01 relative to representative group B gene. (D) Correlation between the mRNA expression of PRDM16 and PGCa1 (left), CITED1 (middle), and ZIC1 (right). These genes represent each gene class. Dotted lines denote density ellipse (95% confidence interval of plot). (E) Box-and-whisker plot (upper) and histogram (bottom) to graphically summarize the distribution of each gene expression levels. The lines extending from end of the box (quartile) are whiskers, which edge is the outermost data point(s) that fall within the distance defined by quartile ±1.5*(IQR). Values beyond the whiskers, denoted by a red point indicate outlier samples.
Transcriptional analyses of human BAT.
(A) Expression profiles of the common brown fat-selective genes (Group A), classical brown fat-selective genes (Group B), and beige cell-selective genes (Group C) in human BAT from multiple adipose depots. The color scale shows the mRNA levels of the genes in green (low or no-expression)-black-red (high-expression) scheme. All genes are shown in the identical scale. WAT: white adipose tissue. M: smooth muscle. Each alphabet indicates the tissues from the same patient. (B) Correlation matrix of beige-selective and classical brown-selective gene expression. The color scale shows the Pearson correlation between each gene mRNA expression in blue (negative correlation, −1.0)-gray (no correlation, 0)-red (positive correlation, 1.0) scheme. (C) Strength of correlation with representative group A gene (PRDM16 and PGC1a) and classical-brown-selective genes or beige cell-selective genes. **P<0.01 relative to representative group B gene. (D) Correlation between the mRNA expression of PRDM16 and PGCa1 (left), CITED1 (middle), and ZIC1 (right). These genes represent each gene class. Dotted lines denote density ellipse (95% confidence interval of plot). (E) Box-and-whisker plot (upper) and histogram (bottom) to graphically summarize the distribution of each gene expression levels. The lines extending from end of the box (quartile) are whiskers, which edge is the outermost data point(s) that fall within the distance defined by quartile ±1.5*(IQR). Values beyond the whiskers, denoted by a red point indicate outlier samples.Next we systematically asked, based on the gene expression profile, how similar human BAT was to the classical brown adipocytes or to the beige cells. First, we used Pearson correlation to evaluate the similarities of expression pattern among classical brown- and beige-selective genes. As shown in Fig. 3B, all the beige selective genes except TBX1 showed high (r>0.87) cross-correlations, justifying the analysis based on a set of gene group designation. We then integrated correlation coefficient between established BAT markers such as PRDM16 and each classical brown- or beige-selective gene into a “gene-set correlation coefficient” according to the groupings; this allowed us to analyze the expression data as “gene-sets” rather than single-gene-based analysis that is sensitive to outlier samples inherent in human postmortem samples. As shown in Fig. 3C (left), the correlation coefficient between PRDM16 and beige cell-selective genes was significantly higher, as compared with those between PRDM16 and classical brown-selective genes (P<0.01, Z-test). Similarly, this difference was highly significant in PGC1a, a well-known common brown adipocyte-enriched gene (P<0.01, Fig. 3C, right). Furthermore, the expression profiles of Group B and C was clearly different at each individual gene level; As shown in Fig.3D, mRNA expression of PRDM16 correlated very well with that of PGC1a (left), and also with that of beige cell-selective genes including CITED1 (middle), but not with that of classical brown fat genes such as ZIC1 (right). mRNA expression of UCP1 was not included in this analysis because outlier plot analysis detected too many outlier samples in UCP1 mRNA levels, partly due to insufficient number of samples. Therefore this did not allow us to run an appropriate correlation analysis using UCP1 mRNA levels. On the other hand, most samples fell within the range of the whisker in PRDM16 and PGC1a mRNA levels. (Fig. 3E). Taken together, these data indicate that molecular signature of human BAT resembles beige cells.
CITED1 is a New Marker for Beige Cells in Mice and in Humans
Although the H&E staining showed that the human BAT samples mostly consisted of brown adipocytes containing multilocular lipid droplets (Fig.2B), qRT-PCR based expression analyses could not exclude the possibility that expression of the beige-selective genes might come from non-brown adipocytes in a given tissue. To address this, we performed immunohistochemistry of a new beige marker, CITED1 in mice and in humans. UCP1 immunohistochemistry was also performed in serial sections to identify the brown adipocytes. CITED1 was originally identified in a murinemelanoma cell line [37] and in humanthyroid carcinoma [38], [39]. It has been reported that CITED1 acts as a transcriptional co-activator for estrogen receptor [40]. As shown in Fig.4A, CITED1 was selectively expressed in the UCP1-positive beige cells in the inguinal WAT in mice chronically treated with a β3 agonist CL316243 (upper panels). In contrast, there were almost no UCP1-positive or CITED1-positive cells in the inguinal WAT in mice treated with saline control (Fig.4A lower panels). Consistent with this, CITED1 was well expressed in the UCP1-positive brown adipocytes in human BAT (Fig. 4B, upper panels). Negative controls for the immunofluorescent stainings of UCP1 and CITED1 (Cy3 and FITC, respectively) confirmed the specificity of the staining (Fig. 4B, lower panels). These data indicate that CITED1 serves as a marker for beige cells both in mice as well as in humans.
Figure 4
Expression of a new beige cell marker CITED1 in mice and in humans.
(A) Immunohistochemistry of UCP1 (left) and CITED1 (right) in serial sections of WAT in mice that were treated with CL316243 at a dose of 1 mg/kg (upper panels) or saline control (lower panels) for 8 days. Scale bar, 50 µm. (B) Immunohistochemistry of UCP1 (left) and CITED1 (right) in serial sections of human BAT. Lower panels show IgG negative controls of immunohistochemistry for UCP1 (Cy3) and for CITED1 (FITC), respectively in human BAT. Scale bar, 50 µm.
Expression of a new beige cell marker CITED1 in mice and in humans.
(A) Immunohistochemistry of UCP1 (left) and CITED1 (right) in serial sections of WAT in mice that were treated with CL316243 at a dose of 1 mg/kg (upper panels) or saline control (lower panels) for 8 days. Scale bar, 50 µm. (B) Immunohistochemistry of UCP1 (left) and CITED1 (right) in serial sections of human BAT. Lower panels show IgG negative controls of immunohistochemistry for UCP1 (Cy3) and for CITED1 (FITC), respectively in human BAT. Scale bar, 50 µm.
Discussion
Two types of brown adipocytes with distinct developmental origins has been previously described in mice: Myf5-derived pre-exiting brown adipocytes in the interscapular BAT and non-Myf5-derived beige cells that emerge in the subcutaneous WAT in response to chronic cold exposure or to PPARγ agonists. The present study characterized the unique molecular signatures for the two types of brown adipocytes, which allowed us to ask whether humans possess beige cells. To our surprise, a number of beige cell-selective genes were abundantly expressed in most, if not all, human BAT isolated from multiple adipose depots, correlating very well with the expression of known common brown fat-selective genes such as PRDM16 and PGC1α, but not with that of classical brown fat-selective genes. Although the current study is limited to molecular analyses, the present data indicate that human BAT may resemble beige cells. This view is in consistent with the recent report by Spiegelman’s group showing that three beige cell markers CD137, TMEM26, and TBX1 are enriched in the humanUCP1-positive brown adipocytes [36]. To be noted, histological analyses of human BAT showed that UCP1-positive brown adipocytes formed a copious cluster and resided together with white adipocytes. These clusters are a morphological characteristic of beige cells in the mouse WAT [3], [4], [5], [6], [41]. Global analyses of human brown adipocytes will further clarify its cellular identity and developmental origin in detail. The number of patient examined in this study was not sufficient to critically examine the age-related changes in the expression of brown fat-selective genes, however, we will investigate this point in our future studies.We found strong correlations among beige cell-selective gene expressions in human BAT, indicating that there is a dominant regulatory pathway for beige cell development. To be noted, expressions of PRDM16 correlated very well with those of most beige cell genes. PRDM16 is a 140 kDa zinc-finger protein containing PR-(PRD1-BF-1-RIZ1 homologous) domain that was originally identified at a chromosomal breakpoint in humanmyeloid leukemia cells [42], [43]. Indeed, it has been shown that PRDM16 expression in primary and immortalized white adipocytes is sufficient to activate a large array of brown-fat selective genes via interacting with several transcriptional factors such as PPARγ and PGC1α [26], [32], [44]. In addition, PRDM16 expression together with C/EBPβ is able to convert human foreskin-derived fibroblasts into UCP1-positive brown adipocytes [16]. Our recent study has also shown that PRDM16 is required for the development of beige cells (so called “browning” of white fat) in response to PPARγ agonists in mice [28]. These studies indicate that PRDM16 is one of the dominant factors for the formation of human BAT, however, the requirement of PRDM16 in humans remains to be tested.The development of beige cells can occur in adult mice and is highly inducible in response to a variety of environmental cues and endogenous hormones, therefore, this cell type possesses a tremendous potential as a therapeutic target. The present study indicates the relevance of this cell type in humans, further supporting this idea. One of the emerging questions would be to understand the developmental processes and molecular control of the beige cells. Wu et al., have recently cloned beige-cell lines from mouse inguinal WAT and showed that beige-precursors expressed unique sets of genes that could distinguish beige cells from white or classical brown adipocytes [36]. Meanwhile, Cinti’s group has shown that the browning of white fat in response to cold is mainly due to the transdifferentiation of mature white adipocytes into brown adipocytes in mice [45]. It needs to be addressed how proliferation of beige cell precursors and white-to-brown adipocyte transdifferentiation contributes to the development of human BAT.Previous attempts to activate thermogenesis in existing BAT by beta3-adrenergic receptor agonists have been unsuccessful in humans [46]. Because human BAT possesses molecular characteristics of beige cells, a future therapeutic strategy for humanobesity would be to induce the browning of white fat. For example, inhibition of TGF-β signaling by neutralizing antibodies sufficiently increases whole body energy expenditure and counteracts diet-induced obesity, at least in part, by browning white fat in mice [47], [48]. Alternatively, using endogenous hormones such as FGF21 [49], natriuretic peptides [50], prostaglandin [51], [52], or Irisin [53] to promote browning of white fat would be a plausible approach.Expression profile after normalizing the each gene mRNA level with that of primary skin fibroblasts.(TIF)Click here for additional data file.Primer sequences used in quantitative real-time PCR.(PDF)Click here for additional data file.Patient and sample information for the tissues used in the expression analyses.(PDF)Click here for additional data file.Genes enriched in rosiglitazone-inducible beige/brite cells.(PDF)Click here for additional data file.Genes enriched in classical brown adipocytes(PDF)Click here for additional data file.
Authors: G Barbatelli; I Murano; L Madsen; Q Hao; M Jimenez; K Kristiansen; J P Giacobino; R De Matteis; S Cinti Journal: Am J Physiol Endocrinol Metab Date: 2010-03-30 Impact factor: 4.310
Authors: Wouter D van Marken Lichtenbelt; Joost W Vanhommerig; Nanda M Smulders; Jamie M A F L Drossaerts; Gerrit J Kemerink; Nicole D Bouvy; Patrick Schrauwen; G J Jaap Teule Journal: N Engl J Med Date: 2009-04-09 Impact factor: 91.245
Authors: Aaron M Cypess; Sanaz Lehman; Gethin Williams; Ilan Tal; Dean Rodman; Allison B Goldfine; Frank C Kuo; Edwin L Palmer; Yu-Hua Tseng; Alessandro Doria; Gerald M Kolodny; C Ronald Kahn Journal: N Engl J Med Date: 2009-04-09 Impact factor: 91.245
Authors: Shingo Kajimura; Patrick Seale; Takuya Tomaru; Hediye Erdjument-Bromage; Marcus P Cooper; Jorge L Ruas; Sherry Chin; Paul Tempst; Mitchell A Lazar; Bruce M Spiegelman Journal: Genes Dev Date: 2008-05-15 Impact factor: 11.361
Authors: Pontus Boström; Jun Wu; Mark P Jedrychowski; Anisha Korde; Li Ye; James C Lo; Kyle A Rasbach; Elisabeth Almer Boström; Jang Hyun Choi; Jonathan Z Long; Shingo Kajimura; Maria Cristina Zingaretti; Birgitte F Vind; Hua Tu; Saverio Cinti; Kurt Højlund; Steven P Gygi; Bruce M Spiegelman Journal: Nature Date: 2012-01-11 Impact factor: 49.962
Authors: Andrew L Carey; Renata Pajtak; Melissa F Formosa; Bruce Van Every; David A Bertovic; Mitchell J Anderson; Nina Eikelis; Gavin W Lambert; Victor Kalff; Stephen J Duffy; Martin H Cherk; Bronwyn A Kingwell Journal: Diabetologia Date: 2015-03-01 Impact factor: 10.122
Authors: Paul Cohen; Julia D Levy; Yingying Zhang; Andrea Frontini; Dmitriy P Kolodin; Katrin J Svensson; James C Lo; Xing Zeng; Li Ye; Melin J Khandekar; Jun Wu; Subhadra C Gunawardana; Alexander S Banks; João Paulo G Camporez; Michael J Jurczak; Shingo Kajimura; David W Piston; Diane Mathis; Saverio Cinti; Gerald I Shulman; Patrick Seale; Bruce M Spiegelman Journal: Cell Date: 2014-01-16 Impact factor: 41.582
Authors: Eun Hee Koh; Yong Chen; David A Bader; Mark P Hamilton; Bin He; Brian York; Shingo Kajimura; Sean E McGuire; Sean M Hartig Journal: J Biol Chem Date: 2016-10-10 Impact factor: 5.157
Authors: Laurent Vergnes; Graeme R Davies; Jason Y Lin; Michael W Yeh; Masha J Livhits; Avital Harari; Michael E Symonds; Harold S Sacks; Karen Reue Journal: J Clin Endocrinol Metab Date: 2016-08-30 Impact factor: 5.958