- Short report
- Open Access
Single cell transcriptomics suggest that human adipocyte progenitor cells constitute a homogeneous cell population
Stem Cell Research & Therapyvolume 8, Article number: 250 (2017)
Regulation of adipose tissue stem cells (ASCs) and adipogenesis impact the development of excess body fat-related metabolic complications. Animal studies have suggested the presence of distinct subtypes of ASCs with different differentiation properties. In addition, ASCs are becoming the biggest source of mesenchymal stem cells used in therapies, which requires deep characterization. Using unbiased single cell transcriptomics we aimed to characterize ASC populations in human subcutaneous white adipose tissue (scWAT). The transcriptomes of 574 single cells from the WAT total stroma vascular fraction (SVF) of four healthy women were analyzed by clustering and t-distributed stochastic neighbor embedding visualization. The identified cell populations were then mapped to cell types present in WAT using data from gene expression microarray profiling of flow cytometry-sorted SVF. Cells clustered into four distinct populations: three adipose tissue-resident macrophage subtypes and one large, homogeneous population of ASCs. While pseudotemporal ordering analysis indicated that the ASCs were in slightly different differentiation stages, the differences in gene expression were small and could not distinguish distinct ASC subtypes. Altogether, in healthy individuals, ASCs seem to constitute a single homogeneous cell population that cannot be subdivided by single cell transcriptomics, suggesting a common origin for human adipocytes in scWAT.
White adipose tissue (WAT) dysfunction is central to the pathologies associated with overweight/obesity such as insulin resistance, type 2 diabetes, dyslipidemia, and atherosclerosis. Obesity and insulin resistance are characterized by increased fat cell size, changes in lipid/glucose metabolism, as well as increased infiltration of leukocytes, primarily macrophages [1, 2].
Adipocytes develop from adipocyte stem cells (ASCs) in a process termed adipogenesis. Given that ~ 10% of the adipocyte pool is renewed annually in adult humans, altered adipogenesis might impact on adipose tissue function [3, 4]. This notion is supported by the observation that hypertrophic WAT (few, large fat cells), in comparison to a hyperplastic phenotype (many small fat cells), is closely linked to low adipocyte turnover and a pernicious metabolic profile [3,4,5]. The current view is that hypertrophic obesity develops when ASC differentiation is attenuated, leading to ectopic lipid deposition in peripheral tissues such as the liver, muscle, and vessels . Recent data in humans confirm that adipogenic markers in the entire ASC population correlate with fat cell size and poor metabolic measures . Altogether, dysregulation of ASCs and adipogenesis appear important in the development of metabolic complications to excess body fat.
ASCs are the most abundant cell type in the stroma vascular fraction (SVF) of WAT, and murine studies have suggested the presence of several ASC populations displaying different capacities to undergo adipogenesis [8, 9]. Some markers that enrich for ASCs with marked differentiation capacity have also been identified in humans, including CD34  and CD36 . However, it is still unknown whether different populations of human ASCs are present in vivo.
In addition, ASCs are becoming an important source of mesenchymal stem cells (MSCs) used in allogeneic cellular therapeutics . Therefore, characterization of heterogeneity of human ASCs is of great importance even for cellular therapies.
We aimed to map ASC populations in human subcutaneous WAT by single cell RNA sequencing of the total SVF of healthy individuals.
Materials and methods
Human subjects and adipose tissue dissociation
For single cell RNA sequencing, the SVF from subcutaneous (sc) WAT of four healthy individuals (Table 1) was isolated as described previously . For flow cytometry sorting, SVF from scWAT of six healthy individuals (Table 2) was isolated and prepared as described previously . The study was approved by the regional ethics board and all subjects provided their written informed consent.
Single cell capture and imaging
Amplification, tagmentation, and sequencing
RT and PCR mixes were added to the chip and samples were further processed using the C1 instrument script, which included lysis, reverse transcription, and amplification. cDNA quality was analyzed with an Agilent BioAnalyzer. All procedures including tagmentation and sequencing were as described previously  and in Additional file 1: Supplemental methods.
Single cell RNA-sequencing data from 574 cells were analyzed in a custom Python environment. The data analysis workflow was as described in detail previously . In brief, the following steps were performed: cell selection; clustering of all cells (first-level clustering); t-distributed stochastic neighbor embedding (t-SNE) visualization of all cells; identification of differential expressed genes in cell populations using negative binominal regression; clustering of ASCs (second-level clustering); rare cell detection; and pseudotemporal modeling. All procedures are described in detail in Additional file 1: Supplemental methods. The expression data were corrected for batch effects using ComBat  and normalized according to total molecule number before cubic spline fitting.
Flow cytometry sorting and RNA expression profiling by microarray
Flow cytometry sorting of human WAT SVF was performed as described previously . RNA was prepared from eight different cell WAT cell types (ASCs, total adipose tissue macrophages (ATMs), M1 ATMs, M2 ATMs, total T cells, CD4+ T cells, CD8+ T cells, and mature adipocytes). Ten nanograms of RNA was amplified using four cycles and loaded onto Clariom™D microarray chips. For details see Additional file 1: Supplemental methods. Microarray data have been published in GEO (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE100795; token qxuxgcoojfwdpgd).
To identify ASC subpopulations in human scWAT, we sequenced SVF-derived single cells; 574 cells passed quality control. Subsequently, most variable genes were selected (Additional file 2: Table S1). Cell and gene clustering, as well as heatmap analysis, could separate the cells into four groups, which were present in all individuals (Fig. 1a). t-SNE visualization also suggested four major cell populations (Fig. 1b). We identified the genes that best characterized these cell groups (Fig. 1c) and examined their expression in microarrays from FACS-sorted SVF of scWAT obtained from six different patients (Fig. 1d). This showed that the largest t-SNE population represented ASCs while the remaining three populations mapped to ATMs of M1, M2, and an intermediate subtype (Fig. 1d). Analysis of the single cell transcriptome for established markers specific for ASCs and macrophages confirmed the predicted populations (Fig. 1e).
Although the single cell transcriptomics revealed distinct ATM subtypes, the ASC cluster appeared homogeneous. To identify distinct ASC subtypes we selected the most variable genes within the ASC population (n = 381 cells, genes listed in Additional file 2: Table S2) and performed cell and gene clustering analysis summarized as a heatmap (Fig. 2a) and a t-SNE plot (Fig. 2b). However, no distinct clusters were found using these approaches. Furthermore, analysis designed to search for rare cell types  could also not reveal any subtypes (data not shown).
We considered that even though no specific ASC populations were detected, the cells might represent a gradient of differentiation stages. Therefore, we placed individual cells along a pseudotemporal trajectory in t-SNE space (Fig. 2b) and screened for pseudotime-dependent genes (n = 70 genes over significance cutoff point). Gene expression was visualized using a rolling-wave plot (Fig. 2c). Several collagen genes, CD55, and Thy1 were more highly expressed in the cells localizing to the beginning of the pseudotemporal ordering, while ribosomal genes and KLF4 were enriched toward the end. However, these genes showed only very minor differences in expression within ASCs (Fig. 2d). In conclusion, while we were able to divide ASCs into slightly different pseudotime-related states (Fig. 2e), we found no evidence of distinct ASC subtypes in human WAT.
To our knowledge, this is the first study reporting single cell transcriptomic data of total resident SVF cells from human WAT. Our main finding is that in healthy individuals ASCs seem to constitute a single and homogeneous population without evidence of any distinct subtypes, which is in agreement with earlier study using single cell PCR of presorted ASCs . In contrast, we were able to define several ATM populations which served as an indirect quality control for our data set.
By performing pseudotemporal ordering of the cells we observed a significant gradient of gene expression for 70 genes. Previously characterized markers for committed adipocyte progenitors, such as CD34, PDGFRA, CD29, and CD36 , were not among these. A few other markers known to be regulated by adipogenesis or cell commitment (CD55, SFRP4, SEMA3C) [19,20,21] appeared to be differentially expressed along the pseudotime axis. Unfortunately, expression differences were too small and variable to enable FACS sorting of the cells belonging to the early versus late phase of pseudotemporal order.
Markers for brown (MYF5, PAX5, MYOD1 ) or beige/bright (MYF11 ) adipocyte progenitors were absent in our data set. A few cells expressed low levels of CD24  while no cells expressed VSTM2A , suggesting that neither of these markers label ASCs in humans.
As already mentioned, several murine studies and few human studies have indicated that the ASC population is heterogeneous in the adipogenic capacity [8,9,10,11]; however, entire transcriptomes of the cells purified ex vivo have never been compared. We do not exclude that single markers might be connected with commitment to adipogenesis or osteogenesis. However, such commitment cannot be observed on the single cell transcriptome level, meaning that differences in the gene transcription level in the human ASC population are small. Clear detection of CD36 and CD34 expression as well as PDGFRA and PDGFRB expression suggests that we do not lack populations that were shown earlier to mark cells with different adipogenic capacity/commitment, but these surface markers do not correspond to specific transcriptomes of the cells.
Although our data suggest no major differences in the transcriptome of individual ASCs, we cannot exclude that differences and distinct populations may be present under specific conditions such as insulin resistance and/or detected using other approaches than those used herein (sequencing of SVF samples on the C1 Fluidigm system). In this study, we aimed to visualize heterogeneity of ASCs in healthy individuals and we cannot exclude that morbid obesity or diabetes might induce changes in the ASC population. Finally, because we only examined scWAT we cannot exclude that distinct progenitor populations may exist in visceral WAT or in brown/beige adipose tissue. Furthermore, we cannot exclude the possibility of variability occurring beyond the level of mRNA expression (e.g., epigenetics), or only being detected by gene-specific approaches. These questions were out of the scope of this report and will be addressed by future studies.
Taken together, our single cell transcriptomic approach suggests that ASCs in healthy individuals constitute a homogeneous cell population with only small variations in differentiation state.
Adipose tissue stem cell
Adipose tissue macrophage
Fluorescence-activated cell sorting
Mesenchymal stem cell
Stroma vascular fraction
t-Distributed stochastic neighbor embedding
White adipose tissue
Arner P. Human fat cell lipolysis: biochemistry, regulation and clinical role. Best Pract Res Clin Endocrinol Metab. 2005;19:471–82.
Guilherme A, Virbasius JV, Puri V, Czech MP. Adipocyte dysfunctions linking obesity to insulin resistance and type 2 diabetes. Nat Rev Mol Cell Biol. 2008;9:367–77.
Spalding KL, Arner E, Westermark PO, Bernard S, Buchholz BA, Bergmann O, et al. Dynamics of fat cell turnover in humans. Nature. 2008;453:783–7.
Arner E, Westermark PO, Spalding KL, Britton T, Ryden M, Frisen J, et al. Adipocyte turnover: relevance to human adipose tissue morphology. Diabetes. 2010;59:105–9.
Arner E, Ryden M, Arner P. Tumor necrosis factor alpha and regulation of adipose tissue. N Engl J Med. 2010;362:1151–3.
Gustafson B, Hammarstedt A, Hedjazifar S, Smith U. Restricted adipogenesis in hypertrophic obesity: the role of WISP2, WNT, and BMP4. Diabetes. 2013;62:2997–3004.
Acosta JR, Douagi I, Andersson DP, Backdahl J, Ryden M, Arner P, et al. Increased fat cell size: a major phenotype of subcutaneous white adipose tissue in non-obese individuals with type 2 diabetes. Diabetologia. 2015;59:560–70.
Rodeheffer MS, Birsoy K, Friedman JM. Identification of white adipocyte progenitor cells in vivo. Cell. 2008;135:240–9.
Secco B, Camire E, Briere MA, Caron A, Billong A, Gelinas Y, et al. Amplification of adipogenic commitment by VSTM2A. Cell Rep. 2017;18:93–106.
Sengenes C, Lolmede K, Zakaroff-Girard A, Busse R, Bouloumie A. Preadipocytes in the human subcutaneous adipose tissue display distinct features from the adult mesenchymal and hematopoietic stem cells. J Cell Physiol. 2005;205:114–22.
Gao H, Volat F, Sandhow L, Galitzky J, Nguyen T, Esteve D, et al. CD36 is a marker of human adipocyte progenitors with pronounced adipogenic and triglyceride accumulation potential. Stem Cells. 2017;35:1799-814.
Girdlestone J. Mesenchymal stromal cells with enhanced therapeutic properties. Immunotherapy. 2016;8:1405–16.
Pettersson AM, Stenson BM, Lorente-Cebrian S, Andersson DP, Mejhert N, Kratzel J, et al. LXR is a negative regulator of glucose uptake in human adipocytes. Diabetologia. 2013;56:2044–54.
Islam S, Zeisel A, Joost S, La Manno G, Zajac P, Kasper M, et al. Quantitative single-cell RNA-seq with unique molecular identifiers. Nat Methods. 2014;11:163–6.
Joost S, Zeisel A, Jacob T, Sun X, La Manno G, Lonnerberg P, et al. Single-cell transcriptomics reveals that differentiation and spatial signatures shape epidermal and hair follicle heterogeneity. Cell Syst. 2016;3:e229–37.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8:118–27.
Grun D, Lyubimova A, Kester L, Wiebrands K, Basak O, Sasaki N, et al. Single-cell messenger RNA sequencing reveals rare intestinal cell types. Nature. 2015;525:251–5.
Hardy WR, Moldovan NI, Moldovan L, Livak KJ, Datta K, Goswami C, et al. Transcriptional networks in single perivascular cells sorted from human adipose tissue reveal a hierarchy of mesenchymal stem cells. Stem Cells. 2017;35:1273–89.
Cawthorn WP, Scheller EL, MacDougald OA. Adipose tissue stem cells meet preadipocyte commitment: going back to the future. J Lipid Res. 2012;53:227–46.
Ehrlund A, Mejhert N, Lorente-Cebrian S, Astrom G, Dahlman I, Laurencikiene J, et al. Characterization of the Wnt inhibitors secreted frizzled-related proteins (SFRPs) in human adipose tissue. J Clin Endocrinol Metab. 2013;98:E503–8.
Mejhert N, Wilfling F, Esteve D, Galitzky J, Pellegrinelli V, Kolditz CI, et al. Semaphorin 3C is a novel adipokine linked to extracellular matrix composition. Diabetologia. 2013;56:1792–801.
Sanchez-Gurmaches J, Guertin DA. Adipocytes arise from multiple lineages that are heterogeneously and dynamically distributed. Nat Commun. 2014;5:4099.
Long JZ, Svensson KJ, Tsai L, Zeng X, Roh HC, Kong X, et al. A smooth muscle-like origin for beige adipocytes. Cell Metab. 2014;19:810–20.
The authors thank Sten Linnarsson and Amit Zeisel for advice when setting up the single cell sequencing experiments, Iyadh Douagi and the flow cytometry facility at HERM KI for help in FACS sorting, the core facility BEA for advice in microarray profiling, and Gaby Åström, Kerstin Wåhlen, and Elisabet Dungner for technical help.
This work was supported by the Swedish Research Council, Novo Nordisk Foundation including the Tripartite Immuno-metabolism Consortium (TrIC) Grant Number NNF15CC0018486 and the MSAM consortium NNF15SA0018346, the Swedish Diabetes Foundation, Stockholm County Council, and the Diabetes Research Program at Karolinska Institutet.
Availability of data and material
The microarray dataset was submitted to the GEO database and will be made publically available on the date of manuscript publication.
Single cell sequencing data are available from the corresponding author on reasonable request.
KK’s current address is Departments of Medicine and Genetics, Stanford University, Stanford, CA 94305, USA.
Ethics approval and consent to participate
The study was approved by the regional ethics board at Karolinska Institutet and all subjects provided their written informed consent.
Consent for publication
Written informed consent was obtained from the patients for recording/using their individual details (age, gender, weight, and height) for research purposes. Samples were coded at the collection point and cannot by any means be connected to personal identity. The consent form is held in the patients’ clinical notes and is available for review by the Editor-in-Chief.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.