Lysosomal protein surface expression discriminates fat- from bone-forming human mesenchymal precursor cells

  1. Jiajia Xu
  2. Yiyun Wang
  3. Ching-Yun Hsu
  4. Stefano Negri
  5. Robert J Tower
  6. Yongxing Gao
  7. Ye Tian
  8. Takashi Sono
  9. Carolyn A Meyers
  10. Winters R Hardy
  11. Leslie Chang
  12. Shuaishuai Hu
  13. Nusrat Kahn
  14. Kristen Broderick
  15. Bruno Péault  Is a corresponding author
  16. Aaron W James  Is a corresponding author
  1. Departments of Pathology, Johns Hopkins University, United States
  2. Departments of Orthopaedics, Johns Hopkins University, United States
  3. Department of Oral and Maxillofacial Surgery, School of Stomatology, China Medical University, China
  4. UCLA and Orthopaedic Hospital Department of Orthopaedic Surgery and the Orthopaedic Hospital Research Center, United States
  5. Departments of Plastic Surgery, Johns Hopkins University, United States
  6. Center For Cardiovascular Science and Center for Regenerative Medicine, University of Edinburgh, United Kingdom

Abstract

Tissue resident mesenchymal stem/stromal cells (MSCs) occupy perivascular spaces. Profiling human adipose perivascular mesenchyme with antibody arrays identified 16 novel surface antigens, including endolysosomal protein CD107a. Surface CD107a expression segregates MSCs into functionally distinct subsets. In culture, CD107alow cells demonstrate high colony formation, osteoprogenitor cell frequency, and osteogenic potential. Conversely, CD107ahigh cells include almost exclusively adipocyte progenitor cells. Accordingly, human CD107alow cells drove dramatic bone formation after intramuscular transplantation in mice, and induced spine fusion in rats, whereas CD107ahigh cells did not. CD107a protein trafficking to the cell surface is associated with exocytosis during early adipogenic differentiation. RNA sequencing also suggested that CD107alow cells are precursors of CD107ahigh cells. These results document the molecular and functional diversity of perivascular regenerative cells, and show that relocation to cell surface of a lysosomal protein marks the transition from osteo- to adipogenic potential in native human MSCs, a population of substantial therapeutic interest.

Introduction

Within mammalian white adipose tissue (WAT), a perivascular population of mesenchymal progenitor cells is well documented with respect to multipotency and tissue renewal capabilities (Corselli et al., 2012; Kramann et al., 2016). The ability of human WAT-resident perivascular cells to differentiate into bone-forming osteoblasts and incite or participate in bone repair is also well known (Askarinam et al., 2013; Chung et al., 2014; James et al., 2012a; James et al., 2012b; James et al., 2012c; Lee et al., 2015; Meyers et al., 2018b; Tawonsawatruk et al., 2016) (see (James et al., 2017; James and Péault, 2019) for reviews). The bulk of WAT-resident perivascular cells with mesenchymal progenitor cell attributes reside in the tunica adventitia – the outer collagen-rich sheath of blood vessels (Corselli et al., 2012; James et al., 2012a; West et al., 2016). Microvascular pericytes, although less frequent in absolute numbers, also demonstrate progenitor cell attributes (Chen et al., 2013; Crisan et al., 2009; Crisan et al., 2008). With several recent studies from our group in human (Ding et al., 2019; Hardy et al., 2017) and mouse WAT (Wang et al., 2020), it is clear that perivascular cells, including those found within the tunica adventitia (adventitial cells or adventicytes), demonstrate more phenotypic and functional diversity than previously understood.

CD107a (lysosome-associated membrane protein-1, LAMP1) is a member of a family of structurally related type one membrane proteins predominantly expressed in lysosomes and other intracellular vesicles (Carlsson et al., 1988; de Saint-Vis et al., 1998; Defays et al., 2011; Ramprasad et al., 1996). CD107a is far less frequently expressed on the cell surface, which is the result of both trafficking of nascent protein to the plasma membrane as well as the fusion of late endosomes and lysosomes to the cell membrane (Akasaki et al., 1993; Dell'Angelica et al., 2000). In inflammatory cells, surface CD107a reflects the state of activation (Janvier and Bonifacino, 2005) and has been implicated in cell adhesion (Kannan et al., 1996; Min et al., 2013). In separate reports, CD107a has been described in intracellular vesicles in both osteoblasts and adipocytes (Bandeira et al., 2018; Solberg et al., 2015), yet beyond this, essentially nothing is known regarding CD107a in mesenchymal stem cell fate or differentiation decisions.

Here, antibody array screening of FACS-defined stromal vascular fraction (SVF) perivascular cells identified several novel cell surface antigens, including CD107a, enriched within subpopulations of human adventicytes and pericytes. Flow cytometry and immunohistochemical analyses confirmed that cells with membranous surface CD107a expression reside in a perivascular microanatomical niche within WAT. CD107ahigh cells represent an adipocyte precursor cell, while CD107alow cells represent progenitors with increased osteoblast potential. CD107a trafficking to the cell surface was observed to occur during early adipocyte differentiation – results confirmed by single-cell RNA sequencing datasets from mouse and human adipose tissues. Upon transplantation into immunocompromised rodents, CD107alow cells robustly induce bone formation, both in an intramuscular ectopic ossicle assay in mice and a lumbar spine fusion rat model. These results suggest that cell surface CD107a divides osteoblast from adipocyte perivascular precursors within human tissues.

Results

Identification of CD107a as a novel cell surface antigen expressed among adipose tissue (AT)-resident perivascular stem cells

To identify new markers that may define subsets of perivascular cells, a cell surface antigen screen (Lyoplate) was performed on previously defined perivascular cell fractions (Crisan et al., 2008; James et al., 2012c; Xu et al., 2019), including CD34+CD146- adventitial cells and CD146+CD34- pericytes after exclusion of non-viable, endothelial, and hematopoietic cells (PI+ CD31+ or CD45+ fractions) (Table 1). Several markers were confirmed to be highly expressed among both adventicytes and pericytes, including, for example, the progenitor cell and MSC marker CD90 (Thy-1) and the perivascular cell antigen CD140b (PDGFRβ). Novel markers to divide perivascular progenitors ranged broadly, including the endolysosomal protein CD107a (32% and 82% expression among adventicytes and pericytes, respectively). Another endolysosomal protein, CD107b, was also present on each perivascular cell fraction (13% and 46% expressing adventicytes and pericytes, respectively). Other markers noted to be expressed differentially in subsets of adventicytes and pericytes included CD98, CD140a (PDGFRα), CD142, CD165, CD200, and CD271 (NGFR) (Table 1).

Table 1
Surface antigens expressed within human adventitial cells versus pericytes.

Results derived from Lyoplate analysis of CD34+CD146-CD45-CD31- adventitial cells or CD146+CD34-CD45-CD31- pericytes.

CD markerProtein nameFrequency in adventitial cells
(CD34+CD146-CD45-CD31-)
Frequency in pericytes
(CD146+CD34-CD45-CD31-)
CD90Thy-197%70%
CD91Low-density lipoprotein-related receptor97%61%
CD95Fas receptor (TNFRSF6)42%22%
CD98Large neutral amino acid transporter (LAT1)17%65%
CD105Endoglin47%14%
CD107aLysosomal-associated membrane protein 1 (LAMP1)32%82%
CD107bLysosomal-associated membrane protein 2 (LAMP2)13%46%
CD130Interleukin six beta transmembrane protein39%61%
CD140aPlatelet-derived growth factor receptor alpha (PDGFRA)82%13%
CD140bPlatelet-derived growth factor receptor beta (PDGFRB)97%34%
CD142Tissue factor, PTF, Factor III, or thromboplastin47%75%
CD147Basigin (BSG)99%99%
CD151Raph blood group71%100%
CD164Sialomucin core protein 24 or endolyn91%97%
CD165AD277%21%
CD271Nerve growth factor receptor (NGFR)64%10%

Next, previously derived transcriptomics data on human perivascular cells were analyzed to confirm LAMP1 gene expression, encoding CD107a. Using WAT-derived pericytes (n = 3 samples, GEO dataset: GSE125545) or adventitial perivascular stem cells (n = 3 samples, GEO dataset: GSE130086) (Xu et al., 2019), high expression of the LAMP1 gene was confirmed (mean FPKM values of 9.576 and 9.619, respectively).

Spatial localization of CD107a was next assessed by immunostaining of subcutaneous WAT (Figure 1A–E, n = 3 samples). CD107a immunoreactivity was found most frequently within the outermost layers of larger arteries (Figure 1B) and veins (Figure 1C). Within arteries, the outer tunica adventitia showed a high frequency of CD107a+CD34+ cells (Figure 1B1), while the inner adventitia showed predominantly CD107a-CD34+ cells (Figure 1B2). The smooth muscle media largely did not show CD107a immunoreactivity (Figure 1B3), which was confirmed by dual immunohistochemistry for CD107a and αSMA (Figure 1—figure supplement 1A). Co-expression with the recently described adventitial marker Gli1 was assessed (Kramann et al., 2016), which showed little overlap with CD107a immunoreactive adventitial cells (Figure 1—figure supplement 1B). Smaller caliber arteries (Figure 1D) and veins showed a high frequency of dual expressing CD107a+CD34+ cells within the adventitia. Capillaries within WAT showed some CD107a immunoreactive pericytes, which co-expressed CD146 but not CD31, and were present in an abluminal location (Figure 1E, appearing yellow). CD107a immunoreactivity within the perivascular mesenchymal niche was confirmed across other AT depots, including pericardial, perigonadal, perirenal, and omental human fat (Figure 1—figure supplement 2, n = 3 samples per depot). Small and medium caliber vessels showed perivascular immunoreactivity across all adipose depots.

Figure 1 with 3 supplements see all
Perivascular CD107a expression typifies a subset of perivascular cells within human subcutaneous white adipose tissue (WAT).

Immunofluorescent staining of CD107a (green) and CD34 (red) in human adipose tissue. (A) Tile scan. (B) Larger artery in cross-section, including (B1) outer tunica adventitia, (B2) inner tunica adventitia, and (B3) tunica media and intima. (C) Larger vein in cross-section, including (C1) high magnification of vessel wall. (D) Smaller caliber artery in cross-section, including (D1) high magnification of vessel wall. (E) Capillary in longitudinal cross-section, and (E1) high magnification. (F) Representative FlowJo plot to demonstrate partitioning CD107alowCD31-CD45- and CD107ahighCD31-CD45- fractions from human stromal vascular fraction (SVF). Frequency of CD107alow/high cells across samples is shown in Supplementary file 1 (N = 8 samples). (G) Confirmatory immunofluorescent staining of FAC-sorted CD107alow and CD107ahigh mesenchymal cells (CD31-CD45-PI- cells). (H,I) Representative flow cytometry analysis of freshly isolated CD107alow and CD107ahigh mesenchymal cells, including CD107a, CD34, and CD146. Frequency of expression is shown in relation to isotype control (blue vs. red lines). Frequency of CD34+ and CD146+ cells across all samples is shown in Supplementary file 2, 3 (N = 4 samples). Scale bars: 500 μm (A), 50 μm (B,C,G) and 10 μm (B1–B3,C1,D,D1,E,E1).

Next, flow cytometry demonstrated a spectrum of CD107a membranous staining across the viable, non-endothelial/noninflammatory cells of human WAT (Figure 1F). The PI-CD31-CD45- component of SVF was divided by FACS into CD107alow and CD107ahigh cell populations for further analysis (Figure 1F,G). Mean frequency of CD107alow mesenchymal cells represented 33.75% of viable SVF, while mean frequency of CD107ahigh mesenchymal cells represented 5.20% (Supplementary file 1). Flow cytometry analysis was next performed within CD107ahighCD31-CD45- and CD107alowCD31-CD45- mesenchymal populations (Figure 1H,I). High expression of CD107a was confirmed by flow cytometry among freshly sorted CD107ahigh cell preparations (Figure 1I). Concordant with histological observations, CD34+ and CD146+ cells were identified in both CD107alow and CD107ahigh cell fractions (Figure 1H,I, mean frequencies reported in Supplementary file 2, 3). These data confirmed that surface CD107a expression is present in both perivascular native MSC niches within WAT, and that surface CD107a could be used to prospectively purify mesenchymal cell subpopulations with disparate staining intensities.

Next, canonical markers of culture-expanded human MSCs were examined by flow cytometry in freshly sorted CD107alow and CD107ahigh cells, including CD44, CD73, CD90, and CD105 (Supplementary file 4, Figure 1—figure supplement 3). With the exception of CD105, all markers showed overall similar expression patterns across CD107alow and CD107ahigh cell populations (n = 3 samples per group). CD105 expression was disproportionately present within CD107ahigh mesenchymal cells (mean frequency 0.44% and 8.12% among CD107alow and CD107ahigh cells, respectively).

CD107alow AT-derived stromal cells represent osteoblast precursor cells

CD107alow and CD107ahigh cells were again derived from the CD31-CD45- fraction of adipose tissue samples, and in vitro properties examined (Figure 2). Morphology of adherent CD107alow and CD107ahigh cells was broadly similar, with a fibroblastic shape (Figure 2—figure supplement 1). CD107alow cells demonstrated a higher proliferative rate in comparison to CD107ahigh cells (Figure 2A). Among freshly isolated cells, the vast majority of colony forming units-fibroblast (CFU-F) was identified within the CD107alow cell fraction (Figure 2B,C). Among equivalent cells at passage 4, an enrichment in CFU-F was still observed in CD107alow cells (Figure 2D). CFUs-osteoblast (CFU-OB) likewise showed a similar enrichment among CD107alow cells. CFU-OB assays performed in growth medium showed alkaline phosphatase (ALP)+ colonies among CD107alow cells only (Figure 2E,F). The same experiment performed in osteogenic differentiation medium showed an enrichment in CFU-OB among CD107alow cells (Figure 2G,H). Among passaged cells in sub-confluent monolayer, osteogenic differentiation was next examined (Figure 2I–O). ALP staining and quantification demonstrated an enrichment among CD107alow cells (Figure 2I,J). Bone nodule formation was likewise increased in CD107alow cells as compared to CD107ahigh counterparts (Figure 2K,L). Osteogenic gene expression across timepoints of differentiation likewise showed an enrichment for RUNX2 (Runt related transcription factor 2), ALPL, and osteopontin (SPP1) (Figure 2M–O). Thus, the CD107alow mesenchymal component of human WAT contains a precursor cell population with high osteoblastogenic potential.

Figure 2 with 1 supplement see all
Stem/osteoprogenitor cell identity of CD107alow mesenchymal cells.

CD107alowCD31-CD45- and CD107ahighCD31-CD45- cells derived from the same sample of human subcutaneous WAT were exposed to the indicated growth or osteogenic conditions. (A) Cell proliferation among CD107alow and CD107ahigh mesenchymal cells, by MTS assays at 48 hr. (B–D) Fibroblastic colony formation frequency (CFU-F) among human CD107alow and CD107ahigh mesenchymal cells, shown by (B) representative images among freshly isolated cells, (C) CFU-F quantification among freshly isolated cells, (D) CFU-F quantification among passage 4 cells. Whole well images are shown. (E–H) Osteoblastic colony formation frequency (CFU-OB) detected in human CD107alow and CD107ahigh cells. Experiments performed in growth medium (GM) (E,F) or osteogenic differentiation medium (ODM) (G,H). Whole well images are shown. (I,J) Alkaline phosphatase (ALP) staining and photometric quantification at d 10 of osteogenic differentiation among human CD107alow and CD107ahigh cells. Representative whole well and high magnification images are shown. (K,L) Alizarin red (AR) staining and photometric quantification at d 10 of osteogenic differentiation among human CD107alow and CD107ahigh cells. Representative whole well and high magnification images are shown. (M–O) Osteogenic gene expression among human CD107alow and CD107ahigh cells by qRT-PCR at d 3, 7, and 10 of differentiation, including (M) Runt related transcription factor 2 (RUNX2), (N) ALPL, and (O) Osteopontin (SPP1). Osteogenic differentiation examined in N = 3 human cell preparations, and at least experimental triplicate. Dots in scatterplots represent values from individual wells, while mean and one SD are indicated by crosshairs and whiskers. In column graphs, mean values and one SD are shown. *p<0.05; **p<0.01; ***p<0.01 in relation to corresponding CD107alow cell population. Statistical analysis was performed using a two-tailed Student t-test (A–L) or two-way ANOVA followed by Sidak’s multiple comparisons test (M–O). Experiments performed in at least biologic triplicate. Scale bars: 200 μm.

CD107ahigh AT-derived stromal cells represent adipocyte precursor cells

Converse experiments to assay adipogenesis were next performed among cell subsets with differential CD107a expression (Figure 3). Among freshly isolated cells, essentially all CFU-adipocyte (CFU-AD) were found within the CD107ahigh cell population (Figure 3A,B). Next, sub-confluent CD107alow and CD107ahigh cells were propagated under adipogenic differentiation conditions (Figure 3C–G). Oil red O staining was significantly more abundant among the CD107ahigh cells (Figure 3C,D). Adipogenic gene expression was next assessed along adipogenic differentiation (Figure 3E–G). All marker gene transcripts showed significant enrichment among CD107ahigh cells in comparison to CD107alow counterparts, including peroxisome proliferator-activated receptor gamma (PPARG), lipoprotein lipase (LPL), and fatty acid-binding protein 4 (FABP4).

Figure 3 with 4 supplements see all
Adipoprogenitor cell identity of CD107ahigh human mesenchymal cells and correlation to exocytosis during early adipogenic differentiation.

(A–G) CD107alowCD31-CD45- and CD107ahighCD31-CD45- cells derived from human subcutaneous WAT were exposed to adipogenic differentiation conditions. (A,B) Adipocyte colony formation frequency (CFU-AD) detected in human CD107alow and CD107ahigh cells. Whole well images shown. (C,D) Oil red O (ORO) staining and photometric quantification at d7 of adipogenic differentiation among human CD107alow and CD107ahigh cells. Representative whole well and high magnification images shown. (E–G) Adipogenic gene expression among human CD107alow and CD107ahigh cells by qRT-PCR at d 3, 7, and 10 of differentiation, including (E) Peroxisome proliferator-activated receptor-γ (PPARG), (F) Lipoprotein lipase (LPL), and (G) Fatty acid-binding protein 4 (FABP4). (H) Immunocytochemical staining of membranous CD107a in the presence of growth medium (GM) or adipogenic differentiation medium (ADM) after 3 d using human, culture-defined adipose-derived stem cells (ASCs). CD107a immunoreactivity appears red, nuclear counterstain appears blue. (I) Photographic quantification of membranous CD107a immunofluorescence under GM or ADM conditions. (J–L) Induction of membranous CD107a expression after adipogenic differentiation across cell types, including (J) culture-defined human ASCs, (K) FACS-purified human perivascular stem cells (PSC), and (L) culture-defined human BMSCs, assessed by flow cytometry after 3 d under GM or ADM conditions. (M,N) Trafficking of CD107a to the cell surface during adipogenesis was inhibited after treatment with Vacuolin-1 (Vac-1, 1 μM), assessed by CD107a immunostaining (M) and flow cytometry (N) after 3 d under GM and ADM conditions. The cell membrane was labeled using Wheat Germ Agglutinin Conjugates (red), while overlap with CD107a immunostaining appears yellow, and DAPI nuclear counterstain appears blue. (O) Dimensional reduction and unsupervised clustering of human stromal vascular fraction (SVF) adipogenic lineage from subcutaneous WAT revealed three cell groups. (P) Trajectory analyses of human SVF adipogenic lineage, colored based on their unsupervised clustering identity. DPP4 (early) and GGT5 (late) expression were used to identify trajectory origin. (Q) Pseudotemporal cell ordering along differentiation trajectories. Pseudotime is depicted from red to purple. (R) Expression heatmap across pseudotime of genes associated with exocytosis. (S) Combined, normalized expression of exocytosis genes shows enrichment (>1) in early progenitors (green shaded area), while more differentiated cells show reduced average expression. Dots in scatterplots represent values from individual wells, while mean and one SD are indicated by crosshairs and whiskers. In column graphs, mean values and one SD are shown. *p<0.05; **p<0.01; ***p<0.001. Statistical analysis was performed using a two-tailed Student t-test (B,D,I–L) or two-way ANOVA followed by Sidak’s multiple comparisons test (E–G,N). Experiments performed in at least biologic triplicate. Black scale bar: 100 μm; white scale bar: 20 μm.

Finally, to confirm that CD107alow and CD107ahigh mesenchymal cell fractions both represented multipotent precursor cells, parallel chondrogenic differentiation assays were performed in three dimensional micromass culture (Figure 3—figure supplement 1). Here, both cell populations demonstrated a progressive increase in cartilage associated gene expression after 7 d in chondrogenic differentiation conditions. No significant differences in cartilage gene expression were observed between CD107alow and CD107ahigh cell fractions (Figure 3—figure supplement 1A). Alcian blue stained sections of micromass cultures at 21 d of differentiation likewise showed a similar appearance between CD107alow and CD107ahigh cells (Figure 3—figure supplement 1B). Thus, CD107alow and CD107ahigh subpopulations of human SVF both house multipotent mesenchymal cells, but with considerably divergent osteoblastic and adipocytic differentiation potentials.

CD107a traffics to the cell surface during early adipocyte differentiation

Cell surface expression of CD107a results predominantly from trafficking of endolysosomal CD107a+ vesicles to the cell surface. To investigate, unpurified ASCs were exposed to growth conditions or adipogenic differentiation conditions and cell surface CD107a was assessed by immunocytochemistry or flow cytometry (Figure 3H–J). After 3 d exposure to adipogenic conditions, a 12.03 fold increase in immunostaining intensity and a 253.5% increase in the number of CD107ahigh cells were noted. Flow cytometry across several other human cell types confirmed this finding, including FACS-purified perivascular stem cells (PSCs) (Figure 3K) and culture-defined human bone marrow mesenchymal stem cells (BMSCs) (Figure 3L) (188.1–455.4% increase in CD107ahigh cell frequency). Parallel experiments were performed under osteogenic differentiation conditions, which found no significant increase in CD107a staining intensity (Figure 3—figure supplement 2). The increase in membranous CD107a during early adipogenesis was reversed by vacuolin-1 (Vac-1), an inhibitor of Ca2+-dependent fusion of lysosomes with the cell membrane. Treatment with Vac-1 significantly decreased the frequency of CD107ahigh ASCs, and prevented an increase of cell surface CD107a by adipogenic conditions (Figure 3M,N).

In order to confirm that exocytosis is a common feature of early adipogenic differentiation, existing single-cell RNA sequencing datasets of human and mouse AT-derived cells were re-assessed (Merrick et al., 2019). Among human AT-derived cells, re-clustering and cell trajectory analysis identified early progenitor cells (expressing DPP4 and CD55), late progenitor cells (expressing GGT5 and F3), as well as an intermediate cell type (middle progenitor cells) (Figure 3O–Q). KEGG terms for exocytosis showed enrichment within early and middle progenitor cells, as visualized by heatmaps across pseudotime (Figure 3R) and normalized expression of overall exocytosis gene activation (Figure 3S). Similar results linking activation of exocytosis to early adipogenic differentiation were obtained from mouse subcutaneous WAT (Figure 3—figure supplement 3; Merrick et al., 2019). Here, after re-clustering and cell trajectory analysis (Figure 3—figure supplement 3A–C), normalized expression of exocytosis gene activation identified similar trends as in human cells (Figure 3—figure supplement 3D,E). Normalized exocytosis gene activation scores were again enriched within early Dpp4-expressing progenitor cells in comparison to more mature Dlk1-expressing pre-adipocytes.

Knockdown (KD) experiments in human ASCs and CD107alow/high cells did not identify a significant functional role for CD107a in osteo/adipogenesis (Figure 3—figure supplement 4). SiRNA-mediated KD of LAMP1 (encoding CD107a) showed no appreciable effect on the osteogenic differentiation of human ASCs (Figure 3—figure supplement 4A–F). During adipogenesis, LAMP1 KD ASCs and CD107alow adventicytes showed a modest increase in lipid droplet accumulation (Figure 3—figure supplement 4G,K,L,P,Q) and likewise a modest increase in expression of adipogenic marker genes (Figure 3—figure supplement 4H–J,M–O). The high adipogenic differentiation potential of CD107ahigh cells was not significantly altered with LAMP1 KD (Figure 3—figure supplement 4P,Q). Thus, high expression of membranous CD107a, rather than having vital function in cellular differentiation, correlates with exocytosis during early adipogenic differentiation.

Transcriptomic analysis suggests a progenitor cell phenotype for CD107alow cells

Differences in differentiation potential were next investigated using transcriptomic analysis of CD107alow and CD107ahigh stromal cells. RNA sequencing comparative analysis was performed across these two stromal cells. Clear separation between gene expression profiles was observed when comparing CD107alow with CD107ahigh stromal cells, as assessed by principal component analysis (Figure 4A). Further confirming our FACS purification, endothelial and inflammatory marker genes were rarely or not expressed among CD107alow and CD107ahigh stromal cells (Figure 4—figure supplement 1). Progenitor cell markers were expressed among both CD107alow and CD107ahigh stromal cells, with some subtle differences noted (Figure 4B). Transcripts of MYC, LEPR (Leptin receptor), MCAM (CD146), and PDGFRA (Platelet-derived growth factor receptor α) while expressed across all samples were enriched among CD107alow cells. Likewise, NES (Nestin), THY1 (CD90), PDGFRB (Platelet-derived growth factor receptor β) and TBX18 (T-Box transcription factor 18) were expressed across all samples, but more highly among CD107ahigh cells. Other typical MSC markers were more evenly distributed across cell preparations, including CD44, and NT5E (CD73). Consistent with in vitro differentiation potential, genes associated with adipogenic differentiation were highly expressed among CD107ahigh stromal cells, such as FABP4 (Fatty acid-binding protein 4), LPL (Lipoprotein lipase), PPARGC1A (PPARG coactivator 1 α), and CEBPA (CCAAT enhancer binding protein α). In addition, negative regulators of adipogenesis were increased among CD107alow stromal cells, such as KLF2 (Krüppel-like factor 2), KLF3, SIRT1 (Sirtuin 1), and DDIT3 (DNA Damage Inducible Transcript 3) (Figure 4C; Banerjee et al., 2003; Pereira et al., 2004; Sue et al., 2008; Zhou et al., 2015).

Figure 4 with 2 supplements see all
Bulk RNA sequencing among uncultured CD107alow and CD107ahigh mesenchymal cells and relationship to putative adipose cell hierarchy.

(A–F) Total RNA sequencing comparison of CD107alowCD31-CD45- and CD107ahighCD31-CD45- mesenchymal cells from a single human subcutaneous WAT sample. (A) Principal component analysis among CD107alowCD31-CD45- and CD107ahighCD31-CD45- cells. (B) Heat map demonstrating mRNA expression levels of stemness-related markers and perivascular cell markers among CD107alowCD31-CD45- and CD107ahighCD31-CD45- mesenchymal cells. (C) Expression of adipogenic gene markers among CD107alowCD31-CD45- and CD107ahighCD31-CD45- mesenchymal cells, shown in heat map. (D) Ingenuity pathway analysis (IPA) identified representative pathways that were upregulated (Z-score >0; red color) or downregulated (Z-score <0; blue color) in CD107ahighCD31-CD45- compared with CD107alowCD31-CD45- mesenchymal cells. (E,F) Comparison of CD107ahigh/low bulk sequencing data to human SVF single-cell sequencing data (see again Figure 3O–Q). (E) Pathways enriched in both CD107alow bulk RNA-seq and early pseudotime genes derived from scRNA-seq. (F) Pathways enriched in both CD107ahigh bulk RNA-seq and late pseudotime genes derived from scRNA-seq.

QIAGEN Ingenuity Pathway Analysis (IPA) showed that the activated pathways in CD107ahigh stromal cells are associated with the positive regulation of adipogenesis, including, for example, white adipose tissue browning pathway and Sirtuin signaling (Z scores 1.342 and 1.387; CD107ahigh compared with CD107alow stromal cells; Figure 4D; Kurylowicz, 2019). Conversely, upregulated signaling pathways in CD107alow stromal cells included Wnt/β-catenin signaling as well as pathways associated with cellular respiration and metabolism, including Oxidative Phosphorylation and Glutathione metabolism (Z scores −1; CD107ahigh compared with CD107alow stromal cells; Figure 4D and Figure 4—figure supplement 2). In order to further evaluate differences in CD107a cell fractions, pathway analyses were next cross-referenced with prior AT-derived single-cell RNA sequencing data (see again Figure 3O). Highly enriched GO terms among CD107alow stromal cells were likewise found to be enriched among DPP4+ cell fractions during ‘early’ pseudotime (Figure 4E). This included terms associated with Wnt signaling as well as energy metabolism. Conversely, GO terms enriched within CD107ahigh stromal cells were likewise enriched among GGT5+ cell fractions within ‘late’ pseudotime (Figure 4F). This included terms associated with the regulation of cellular proliferation, cell adhesion, as well as remodeling of extracellular matrix.

CD107alow rather than CD107ahigh cells induce ectopic bone formation

We next sought to extend our findings to xenotransplantation studies. If a CD107alow mesenchymal cell population over-represents stem/osteoblast precursor cells, we hypothesized that CD107alow cells would preferentially form ectopic bone within an intramuscular transplantation model (Figure 5; James et al., 2012c; Meyers et al., 2018b). First, CD107alow and CD107ahigh cell subsets were derived from the same patient sample and mixed mechanically with a demineralized bone matrix putty carrier (DBX Putty, MTF Biologics) before intramuscular implantation in NOD-SCID mice. The carrier without cells was used as an acellular control. Details of cell implant composition and animal allocation are summarized in Supplementary file 5. Intramuscular implants were imaged by micro-computed tomography (μCT) at 8 weeks, demonstrating an accumulation of bone tissue among CD107alow laden implants in relation to either CD107ahigh implants or acellular control (Figure 5A). Quantitative μCT analysis demonstrated a significant increase in bone volume (BV, 97.7% increase), fractional bone volume (BV/TV, 73.4% increase), and bone surfaces (BS, 91.4% increase) among CD107alow as compared to CD107ahigh implants (Figure 5B–D). Albeit to a lesser degree, CD107ahigh cells did exhibit bone-forming potential in comparison to acellular control (292.1–473.3% increase in μCT quantitative metrics, Figure 5B–D). Histologic analysis revealed conspicuous areas of woven bone among CD107alow laden implants, which were not commonly seen among CD107ahigh implants (Figure 5E). Bone histomorphometric analysis confirmed these observations, demonstrating significantly increased osteoblast number (N.Ob, 237.5% increase), increased osteoblast number per bone surface (N.Ob/BS, 232.5% increase), and osteocyte number (N.Ot, 460.3% increase) (Figure 5F–H). ALP staining and semi-quantitative analysis confirmed an overall increase in serial sections of CD107alow treated implants (Figure 5I,J, 14.3% increase among CD107alow implant sites). Enrichment in the terminal osteogenic differentiation marker osteocalcin (OCN) was also confirmed among CD107alow implants, shown by immunostaining and semi-quantitative analysis (Figure 5K,L, 345.3% increase among CD107alow implant sites). Detection of human nuclear antigen (HNA) among implant sites confirmed the persistence of human cells across both groups, which were overall similar in frequency (Figure 5—figure supplement 1).

Figure 5 with 1 supplement see all
CD107alow mesenchymal cells promote ectopic bone formation in vivo.

(A) FACS-purified CD107alowCD31-CD45- and CD107ahighCD31-CD45- mesenchymal cells from the same human subcutaneous WAT sample were implanted intramuscularly in equal numbers in the hindlimbs of NOD-SCID mice. A demineralized bone matrix (DBM) carrier was used, and an acellular control used as a further comparison. Bone formation was assayed after eight wks. Further details on implant composition and animal allocation are found in Supplementary file 5. (B) Representative micro-computed tomography (μCT) reconstruction images of the implant site among control (DBM only), CD107alow, and CD107ahigh cell grafts. Mineralized bone appears gray. (C–E) μCT based quantification of ectopic bone formation, including (C) Bone volume (BV), (D) fractional Bone volume (BV/TV), and (E) bone surface (BS). (F) Representative histologic appearance by routine H and E of the implant sites among control (DBM only), CD107alow, and CD107ahigh cell grafts. (G–I) Bone histomorphometric measurements among each treatment group, including (G) osteoblast number (N.Ob), (H) osteoblast number per bone surface (N.Ob/BS), and (I) osteocyte number (N.Ot). (J,K) Representative alkaline phosphatase (ALP) staining appearing blue (J), and photographic quantification (K). (L,M) Representative Osteocalcin (OCN) immunohistochemical staining (L), and photographic quantification (M). OCN immunostaining appears red, while DAPI nuclear counterstain appears blue. Dots in scatterplots represent values from individual implants, while mean and one SD are indicated by crosshairs and whiskers. M, muscle. *p<0.05; **p<0.01; ***p<0.001. Statistical analysis was performed using a one-way ANOVA followed by Tukey’s post hoc test. N = 8 implants per group. Black and white scale bars: 50 μm.

CD107alow rather than CD107ahigh mesenchymal cells induce spine fusion

Having observed that CD107alow cell preparations demonstrate enhanced ectopic bone formation, we next challenged these cells to a posterolateral lumbar spine fusion model within athymic rats (Figure 6; Chung et al., 2014; Lee et al., 2015). CD107alow and CD107ahigh cell subsets from patient-identical samples were implanted bilaterally in an L4-L5 spine fusion model (Figure 6—figure supplement 1). Details of cell implant composition and animal allocation are summarized in Supplementary file 6. A qualitative increase in radiodensity was observed among CD107alow treated animals within the spine fusion bed over the post-operative period by high-resolution roentgenography (Figure 6—figure supplement 2). Progressive increase in density of the implant sites was confirmed by dual-energy X-ray absorptiometry (DXA)-based quantification, with a gradual and significant increase in CD107alow treated spinal implants in comparison to either CD107ahigh treated cells or control (Figure 6A). Fusion rate was next assessed by a validated manual palpation scoring (Figure 6B; Grauer et al., 2001). Consistent with prior studies (Chung et al., 2014), after 8 weeks acellular control-treated animals showed 14.3% fusion (1/7 animals). Analyses performed after 8 weeks demonstrated 62.5% spine fusion among CD107alow treated animals (6/8 animals). In comparison, CD107ahigh treated animals showed 37.5% fusion (3/8 animals). μCT imaging and reconstructions demonstrated lack of bone bridging within the spinal fusion segments of control-treated and CD107ahigh treated implant sites (Figure 6C). In comparison, more robust evidence of bone bridging was observed among CD107alow spine fusion segments (Figure 6C). Quantitative μCT analysis demonstrated a significant increase in bone volume (BV), fractional bone volume (BV/TV), and bone surfaces (BS) among CD107alow implant sites in comparison to acellular control (Figure 6D–F, 58.6–80.7% increase across μCT metrics). In contrast, CD107ahigh spine fusion segments demonstrate no statistically significant change in μCT assessments in comparison to acellular control (Figure 6D–F, 15.3–25.2% change in comparison to acellular control). These findings were confirmed using histologic and histomorphometric assessments of the spinal fusion segment across treatment groups (Figure 6G–J). Histologic analysis revealed conspicuous areas of woven bone among CD107alow implants, which were not commonly seen among CD107ahigh implants (Figure 6G). Bone histomorphometric analysis confirmed these observations, demonstrating significantly increased osteoblast number (172.2% increase among CD107alow implant sites in comparison to CD107ahigh implant sites), increased osteoblast number per bone surface (183.5% increase), and osteocyte number (357.1% increase) (Figure 6H–J). ALP enzymatic staining and OCN immunohistochemical staining confirmed the above findings (Figure 6K–N). In summary, CD107alow but not CD107ahigh mesenchymal cell demonstrate improvements in bone-forming potential across two orthopaedic models.

Figure 6 with 2 supplements see all
CD107alow mesenchymal cells induce spine fusion in vivo.

FACS-purified CD107alowCD31-CD45- and CD107ahighCD31-CD45- mesenchymal cells from the same human subcutaneous WAT sample were implanted in equal numbers in a posterolateral spinal fusion model in athymic rats. A demineralized bone matrix (DBM) carrier was used, and an acellular control used as a further comparison. Animals were analyzed at up to eight wks post-operatively. (A) Bone mineral density (BMD) assessed by DXA (dual-energy X-ray absorptiometry) within the lumbar implantation site, at 0, 4, and 8 wks. (B) Spine fusion rate, assessed by manual palpation after eight wks. *: CD107alow compared with acellular control. (C) Representative micro-computed tomography (μCT) reconstruction images of the spine fusion site among CD107alow and CD107ahigh treated samples, in comparison to acellular control. Images are shown from the dorsal aspect. (C1) Corresponding high magnification μCT reconstruction of the fusion site. (C2) Corresponding coronal μCT cross-sectional image. (D–F) μCT-based quantification of bone formation within the spine fusion site, including (D) Bone volume (BV), (E) fractional Bone volume (BV/TV), and (F) bone surface (BS). (G) Representative histologic appearance by routine H and E of the implant sites among control (DBM only), CD107alow, and CD107ahigh cell grafts within the spine fusion site. (H–J) Bone histomorphometric measurements among each treatment group, including (H) osteoblast number (N.Ob), (I) osteoblast number per bone surface (N.Ob/BS), and (J) osteocyte number (N.Ot). (K,L) Representative alkaline phosphatase (ALP) staining appearing blue (K), and photographic quantification within the spine fusion site (L). (M) Representative Osteocalcin (OCN) immunohistochemical staining (M), and photographic quantification within the spine fusion site (N). OCN immunostaining appears red, while DAPI nuclear counterstain appears blue. Dots in scatterplots represent values from individual animal measurements, while mean and one SD are indicated by crosshairs and whiskers. #p<0.05 and ###p<0.001 in relation to corresponding 0 wk timepoint; *p<0.05; **p<0.01; ***p<0.001. Statistical analysis was performed using a two-way ANOVA (A) or one-way ANOVA followed by Tukey’s post hoc test (D–N). N = 6–8 animals per group. Black and white scale bars: 50 μm.

Discussion

Mesenchymal progenitor cells are broadly distributed in post-natal organs, where they are concentrated principally in perivascular areas. Microvascular pericytes were first recognized to include such progenitors since they grow into mesenchymal stem cells in culture (Crisan et al., 2008). The outermost layer enwrapping arteries and veins, or tunica adventitia, that used to be considered as a mere fibroblast-populated collagen sheath anchoring vessels within tissues, is also home to presumptive MSCs (Corselli et al., 2012; Kramann et al., 2016). Here, we have used an antibody array targeting all human CD surface markers to identify several novel antigens expressed by human adipose tissue-resident perivascular cells. We found, among other surface antigens, that CD107a, aka LAMP-1, is expressed at the surface of subsets of adventitial cells and pericytes, which was confirmed in terms of gene expression on RNA sequencing libraries, and corroborated by immunohistochemistry, surface CD107a being co-expressed with canonical markers of pericytes and adventitial cells. Further, flow cytometry analysis of the total stromal vascular fraction extracted from adipose tissue showed a continuum of non-endothelial, non-hematopoietic CD107a expressing cells that could be gated back to pericytes and adventitial cells. Altogether, these results confirmed unequivocally that surface CD107a/LAMP-1, used generally as a marker of NK cell activity (Bryceson et al., 2005), is present on subsets of human MSC-related perivascular cells, and established the conditions for the purification of these cells by flow cytometry.

Albeit not described before at the surface of human pericytes and other perivascular cells, CD107a has been earlier detected on human bone marrow and dental conventionally cultured MSCs, where it binds the enamel matrix protein amelogenin and in turn induces cell proliferation (Huang et al., 2010). Besides, CD146+CD107a+ human bone marrow MSCs have been recently described as endowed with the highest immunomodulatory and secretory, hence therapeutic, potential in experimental joint inflammation (Bowles et al., 2019). Besides suggested association of surface CD107a with progenitor cell proliferation and immunomodulation related tissue repair (Bowles et al., 2019; Huang et al., 2010), lysosomal CD107a has also been linked to neural stem cell potential (Yagi et al., 2010). Here, we show that surface expression of CD107a divides adipocyte- from osteoblast precursors within human perivasculature. Adipocyte progenitor distribution was virtually confined to the CD107ahigh cell subset; in agreement, adipogenic potential in culture was restricted to this cell compartment. However, knockdown of CD107a in ASCs slightly promoted adipogenic differentiation, suggesting that CD107a can be used for identification of functionally relevant subsets, which is likely not explained by intrinsic function of CD107a protein. CD107a is thought to be responsible for maintaining the structural integrity of the lysosomal compartment (Eskelinen, 2006). Lysosomes provide the degradative enzymes for autophagy and are involved in autophagy regulation, primarily through its relationship with the master kinase complex, mTORC1 (Mrschtik and Ryan, 2015; Yim and Mizushima, 2020). Therefore, CD107a may also play an essential role in autophagy, which itself has been shown crucial for adipocyte differentiation (Guo et al., 2013). Whether CD107ahigh cells enhance adipogenesis by promoting autophagy is an interesting and worthy follow-up topic. Conversely, osteogenic ability in vitro was almost totally restricted to purified CD107alow cells, which also exhibited the highest CFU-F potential. We confirmed the exclusive bone-forming potential of CD107alow cells in vivo, in situations of ectopic intramuscular ossification and lumbar spine fusion. Thus, CD107alow cells likely represent a more progenitor cell, while CD107ahigh cells are a more mature, differentiated cell population. Similar results are found in the hematopoietic system, where the most primitive hematopoietic stem cells are Thy-1low, whereas Thy-1high cells belong to a well-defined blood cell lineage (Spangrude et al., 1988).

The natural function of these ubiquitous perivascular mesenchymal progenitor cells remains a dominant, yet unanswered question. While subsets of mouse pericytes are known, in contexts of experimental disease or injury, to give rise in situ to white adipocytes, follicular dendritic cells, myoblasts, and myofibroblasts (Dellavalle et al., 2011; Dulauroy et al., 2012; Krautler et al., 2012; Murray et al., 2017; Tang et al., 2008), the osteogenic and chondrogenic potentials present in these microvascular cells are unlikely to be ever used in the turnover and repair of soft tissues. On the other hand, a contribution to blood vessel pathologic remodeling of perivascular presumptive MSCs has been recognized, as osteoblast and smooth muscle cell forerunners for calcific arteriosclerosis and atherosclerosis, respectively (Kramann et al., 2016). How such pathogenic side effects exerted by perivascular progenitors are counterbalanced by beneficial contributions to tissue homeostasis is not known, but might be related to differential recruitment of uni- or multipotent progenitors, as well as activation in these cells of mechanisms independent of cell differentiation such as growth factor production or cell contact-mediated immunomodulation (Pittenger et al., 2019). A precise phenotypic and functional characterization of perivascular progenitor cells is a requisite for the understanding of the role of these cells in developmental, regenerative, and pathophysiologic processes. In particular, unwrapping the intrinsic heterogeneity of these cells is a prioritized, ongoing process.

Although abundant evidence of the osteogenic potential of perivascular cells exists in humans and mice (Askarinam et al., 2013; Chung et al., 2014; Corselli et al., 2012; Crisan et al., 2008; James et al., 2017; James and Péault, 2019; James et al., 2012a; James et al., 2012a; James et al., 2012b; James et al., 2012c; Kramann et al., 2016; Lee et al., 2015; Meyers et al., 2018b; Meyers et al., 2018b; Tawonsawatruk et al., 2016; Wang et al., 2020; Xu et al., 2019), there is recent proof that this competence is restricted to discrete cell subsets. For instance, CD10 expression marks a subset of human adipose tissue adventitial cells with higher bone-forming potential (Ding et al., 2019), and the highest calcification potential is attributed to mouse adventitial cells co-expressing CD34 and PDGFRα (Wang et al., 2020). On the other hand, human perivascular cells expressing the ROR2 Wnt receptor exhibit stronger chondrogenic ability than ROR2 negative counterparts (Dickinson et al., 2017). Therefore, the concept is emerging of a developmental micro-heterogeneity of perivascular cells, the surrounding niches being suggested to host a whole hierarchy of mesenchymal progenitor cells already documented to include adipocyte-, chondrocyte-, and osteoblast progenitors. Although the data herein were obtained on abdominal subcutaneous adipose tissue, our preliminary results indicate that the same cell partition can be achieved in other fat depots. Moreover, we have observed the same restriction of adipogenic and osteogenic potentials, respectively, to CD107ahigh and CD107alow perivascular cells sorted from the human placenta, which suggests the broad and possibly ubiquitous distribution of these functionally distinct subpopulations. Strikingly, these tissues are not natural sites of ossification, which raises the recurring question of the physiologic significance of these unrelated developmental potentials. It is conceivable that mesodermal cell turnover and regeneration in adult tissues be mediated exclusively by multipotent, MSC like cells, irrelevant potentials in a given tissue being always repressed. It is more difficult to justify that progenitors committed to a given cell compartment be maintained in a tissue devoid of this very cell lineage, such as osteogenic dedicated progenitors in adipose tissue, at the expense of adipogenic cells. Since the hypothesis that these unrelated progenitors can be mobilized through blood circulation to drive regeneration in other organs is not supported for the moment, we can only speculate that such atypical differentiation potentials are irreversibly associated with other tissue repair mechanisms of broader applicability, which remain to be identified. Further characterization of the role of CD107a at the cell surface may contribute to clarifying this issue.

Our in vitro studies, coupled with single-cell transcriptomics, suggest that CD107a, an endolysosome transmembrane protein, traffics to the cell surface during early adipogenesis, suggesting a specific function in this cell lineage. Future studies will tell whether CD107a/LAMP-1, a unique novel marker of the perivascular mesenchymal stem cell hierarchy, will also shed light on the tissue regeneration mechanisms initiated in this niche.

Materials and methods

Immunohistochemistry and microscopy

Request a detailed protocol

All human tissues were obtained under Johns Hopkins University institutional IRB approval with a waiver of informed consent. For histology, human subcutaneous fat tissue was obtained from three anonymized female donors from the abdominal or thigh area. Human fat tissue was embedded in optimal cutting temperature compound (OCT) (Sakura, Torrance, CA), and cryo-sectioned at 30 μm thickness. For immunofluorescent staining, all sections were blocked with 5% goat serum in PBS for 1 hr at room temperature (RT). The following primary antibodies were used: anti-αSMA (RRID:AB_1951138, 1:100), anti-CD107a (RRID:AB_1727417/RRID:AB_10719137/RRID:AB_470708/RRID:AB_449893, 1:100), anti-CD31 (RRID:AB_448167/RRID:AB_726362, 1:100), anti-CD34 (RRID:AB_1640331, 1:100), anti-CD146 (RRID:AB_2143375, 1:100), or anti-Gli1 (RRID:AB_880198, 1:100; see antibody details in Supplementary file 7) for overnight incubation at 4o C. Next, anti-rabbit Alexa Fluor 647-conjugated (RRID:AB_2722623, 1:200), anti-mouse Alexa Fluor 488-conjugated (RRID:AB_2688012, 1:200), anti-goat Alexa Fluor 647-conjugated (RRID:AB_2687955, 1:200), or anti-rat Alexa Fluor 647-conjugated secondary antibodies (RRID:AB_2864291, 1:200, Abcam, San Francisco, CA) were used (incubation 2 hr at RT). DAPI mounting medium was used (RRID:AB_2336788, Vector laboratories, Burlingame, CA), and visualized using a Zeiss 800 confocal microscope (Zeiss, Thornwood, NY). For colorimetric immunohistochemistry staining, sections were blocked with 2.5% horse serum for 20 min at RT. Anti-CD107a primary antibody (RRID:AB_470708, Abcam, 1:100) was added and incubated overnight at 4° C. Next, incubation with alkaline phosphatase (AP) polymer anti-mouse IgG reagent was performed for 30 min (RRID:AB_2336535, MP-5402, Vector laboratories), followed by AP substrate solution (RRID:AB_2336847, SK-5100, Vector laboratories), followed by hematoxylin counterstain and microscopic imaging using a Leica DM6 B microscope (Leica Microsystems Inc, Wetzlar, Germany).

Adipose-derived stromal cells (ASCs) isolation and FACS isolation of human AT cell populations

Request a detailed protocol

For cell isolation, human lipoaspirate was obtained from healthy adult donors and was stored no more than 72 hr at 4°C before processing with some modifications from prior protocols (Wang et al., 2019; Xu et al., 2019). Patient gender and donor area are shown in Supplementary file 1. Equal volume of phosphate-buffered saline (PBS) was used to wash lipoaspirates. The washed lipoaspirate was digested with 1 mg/ml type II collagenase (Washington Biochemical; Lakewood, NJ) in Dulbecco's Modified Eagle's Medium (DMEM) containing 0.5% bovine serum albumin (Sigma-Aldrich, St. Louis, MO) at 37°C for 60 min under agitation, followed by centrifugation to remove adipocytes. The cell pellet was resuspended and incubated in red cell lysis buffer (155 mM NH4Cl, 10 mM KHCO3, and 0.1 mM EDTA) at RT for 10 min. After centrifugation, the stromal vascular fraction (SVF) was resuspended in PBS and filtered at 40 μm. In select studies, SVF was culture- expanded as adipose-derived stromal cells (ASCs) for further evaluation. Human ASCs were cultured and used for experiments at passage 2–6. Fluorescence activated cell sorting (FACS) was next performed using a Beckman MoFlo (Beckman, Indianapolis, IN), with analysis performed using the FlowJo software (RRID:SCR_008520, BD Biosciences, San Jose, CA). Cells were incubated with anti-CD45-allophycocyanin-cyanin 7 (RRID:AB_396891, 1:30; BD Pharmingen, San Diego, CA), anti-CD31-allophycocyanin-cyanin 7 (RRID:AB_2738350, 1:100, BD Pharmingen), anti-CD34 (RRID:AB_11154586, 1:60, BD Pharmingen), anti-CD146 (RRID:AB_324069, 1:100, Bio-Rad, Hercules, CA), and/or anti-CD107a-allophycocyanin (RRID:AB_1727417, 1:20; BD Pharmingen), for 20 min on ice. Propidium iodide (PI) staining solution (BD Pharmingen) was used to gate out non-viable cells. See Supplementary file 7 for a list of antibodies used. Gating was performed to isolate either CD146+ pericytes (CD146+CD34-CD31-CD45-), CD34+ adventicytes (CD34+CD146-CD31-CD45-), CD107alow cells (CD107alowCD31-CD45-) or CD107ahigh cells (CD107ahighCD31-CD45-). Cell surface markers were analyzed using either Lyoplate (BD Biosciences) or flow cytometry. For flow cytometry, cells were incubated with the following antibodies for 20 min on ice: anti-CD34 PE-CF594 (RRID:AB_11154586), anti-CD146 FITC (RRID:AB_324069), anti-CD44 Alexa Fluor 700 (RRID:AB_10645788), anti-CD73 PE (RRID:AB_2033967), anti-CD90 FITC (RRID:AB_395969), anti-CD105 PE-CF594 (RRID:AB_11154054), and anti-CD107a APC (RRID:AB_1727417, Supplementary file 7). For select studies using culture-expanded cells, flow cytometry was performed after trypsinization and cell re-suspension in HBSS (Life Technologies, Gaithersburg, MD) with 0.5% bovine serum albumin (Sigma-Aldrich). All cells were cultured at 37°C in a humidified atmosphere containing 95% air and 5% CO2. Unless otherwise stated, cells were cultured in Endothelial Cell Growth Medium-2 (EGM-2; Lonza, Gaithersburg, MD).

Human bone marrow mesenchymal stem cell isolation

Request a detailed protocol

Bone marrow mesenchymal stem cells (BMSCs) from anonymized human femur and tibia were isolated using previously reported methods (Xu et al., 2019). Marrow cells were flushed with PBS and passed through a 70 μm cell strainer (BD Biosciences) to obtain single cells, which were seeded into T75 flasks. Non-adherent cells were removed after 5 d and medium was changed every 3 d. BMSCs were cultured in growth medium consisting of DMEM, 15% fetal bovine serum (FBS; Gibco, Grand Island, NY), 1% penicillin/streptomycin (Gibco).

Identification of novel human perivascular cell markers using lyoplate

Request a detailed protocol

The BD Lyoplate Human Cell Surface Marker Screening Panel contains 242 purified and lyophilized monoclonal antibodies to cell surface markers, along with AlexaFluor 647 conjugated goat anti-mouse Ig and goat anti-rat Ig secondary antibodies, distributed in three 96-well plates, as well as mouse and rat isotype controls for assessing isotype-specific background. The Lyoplate array was used according to manufacturer’s instructions. Aspirated human subcutaneous fat was digested with collagenase, washed by centrifugation and the SVF recovered as described above. After washing and red cell lysis, SVF cells were stained with the following reagents: DAPI 1/100 (1 μg/ml final), FITC-CD146 1/40, PE-CD45 1/20, PE-Cy7-CD34 1/33. Using a multi-channel pipette, 100 μl aliquots of antibody stained SVF (500,000 to 1 million cells) were distributed in the wells, and 50 μl of reconstituted Lyoplate antibody solution were added to each well according to the template. Plates were incubated on ice in the dark for 30 min, then cells were washed twice by adding 100 μl of staining solution to each well and spinning at 300xg for 5 min. 100 μl of 4% paraformaldehyde (PFA) were added to each well and incubated at RT for 30 min. Labeled cells were washed again and either stored at 4°C or analyzed directly on a LSR II flow cytometer (BD Biosciences).

Transcriptomics

Request a detailed protocol

In select experiments, global gene expression analysis of CD107alowCD31-CD45- and CD107ahighCD31-CD45- cells from adipose tissue was performed. The RNA content of CD107alow and CD107ahigh cells was detected by total RNA sequencing. Briefly, total RNA was extracted from CD107alow and CD107ahigh cells with Trizol (Life Technologies). After purification and reverse transcription, cDNA samples were sent to the JHMI Transcriptomics and Deep Sequencing Core and quantified by deep sequencing with the Illumina NextSeq 500 platform (Illumina, San Diego, CA). Data were analyzed using software packages including Partek Genomics Suite (RRID:SCR_011860), Spotfire DecisionSite with Functional Genomics (RRID:SCR_008858), and QIAGEN Ingenuity Pathway Analysis (RRID:SCR_008653).

Single-cell RNA sequencing (scRNA-seq)

Request a detailed protocol

ScRNA-seq data were obtained from the Gene Expression Omnibus (GEO) repository, accession number GSE128889 (GSM3717979, GSM3717977). Initial quality control removed cells expressing >200 and<6000 genes and a mitochondrial content >5%. Data normalization, dimensional reduction and clustering were conducted in Seurat (RRID:SCR_016341) as previously described in the original publication with the exception of altered clustering resolutions. Trajectory plots were generated in Monocle (RRID:SCR_018685) as previously described. For exocytosis pathway activation, a gene list of exocytosis activating/promoting genes was generated from the previously annotated KEGG pathway. Gene lists were filtered for genes that met the Monocle cutoff criteria as an expressed gene (expressed in a minimum of 10 cells). Expression levels were normalized first for individual cell UMI counts, then to the average gene expression across the whole sample population. This normalized expression across pseudotime was displayed following average with nearest neighbors. Values above one indicate enriched expression of genes associated with exocytosis while values below one indicate expression below population averages. Pathway analyses were conducted on CD107alow and CD107ahigh cells from bulk RNA-seq experiments, as well as on early- and late-expressing, pseudotemporally-regulated genes derived from scRNA-seq data.

Osteogenic differentiation assay and ALP/Alizarin red staining

Request a detailed protocol

For osteogenic differentiation, cells were seeded at the density of 2.5 × 105/ml in 12- or 24-well plates. Upon confluency, medium was changed to osteogenic differentiation medium, composed of DMEM, 10% FBS (Gibco), 1% penicillin/streptomycin, with 50 μM ascorbic acid, 10 mM β-glycerophosphate, and 100 nM dexamethasone (Sigma-Aldrich) (Xu et al., 2019). Medium was changed every 3 d. Leukocyte Alkaline Phosphatase Kit (Sigma-Aldrich) was used for alkaline phosphatase staining at 3 or 7 d of differentiation. Image J (RRID:SCR_003070) was used to detect the integrated density of ALP staining. Alizarin red S (Sigma-Aldrich) was used to stain cultures at 7 or 10 d of differentiation to detect mineralization. Next, calcium precipitate was dissolved with 0.1N sodium hydroxide and quantified by absorbance at 548 nm. Experiments were performed in N = 3 biological replicates, and in experimental triplicates in each case.

Adipogenic differentiation assay and oil red O staining

Request a detailed protocol

For adipogenic differentiation, cells were seeded at the density of 2.5 × 105/ml in 12- or 24-well plates. Upon subconfluency, adipogenic differentiation medium (DMEM, 1% penicillin/streptomycin, 10% FBS with 1 μM dexamethasone, 200 μM indomethacin, 10 μg/ml insulin, and 500 μM 3-isobutyl-1-methylxanthine (Sigma-Aldrich)) was used. Medium was changed every three d. Cells were fixed with 4% PFA and stained with Oil Red O solution at 5 to 7 d of differentiation (Meyers et al., 2018a). Experiments were performed in N = 3 biological replicates, and in experimental triplicates in each case.

Chondrogenic differentiation assay and Alcian blue staining

Request a detailed protocol

For chondrogenic differentiation, cells were seeded at high-density micromass environment (1 × 107/ml, 10 μl/drop) in 12-well plates and cultured in 37°C. After 4 hr, chondrogenic differentiation medium (DMEM, 1% penicillin/streptomycin, 10% FBS with 10 ng/ml transforming growth factor-β3 (R and D Systems, Minneapolis, MN), 100x ITS+ Premix (Corning Incorporated, Corning, NY), 50 μg/ml ascorbic acid, 40 μg/ml proline, 100 μg/ml pyruvate, and 100 nM dexamethasone (Sigma-Aldrich)) was added. Medium was changed every three d. Cells were fixed with 4% PFA and embedded in OCT for cryosectioning at 18 μm thickness. Slides were stained with Alcian Blue and Fast Red.

Proliferation assay

Request a detailed protocol

Proliferation assays were performed in 96 well plates (2 × 103 cells/well) and measured at 72 hr using the CellTiter96 AQueous One Solution Cell Proliferation Assay kit (MTS, G358A; Promega, Madison, WI). Briefly, 20 μl of MTS solution was added to each well and incubated for 1 hr at 37°C. The absorbance was assayed at 490 nm using Epoch microspectrophotometer (Bio-Tek, Winooski, VT).

CFU assay

Request a detailed protocol

For all CFU assays, 1,000 cells / well were seeded in 6-well plates. For CFU-F analysis, cells were cultivated for 14 d in growth medium, fixed with 100% methanol and stained with 0.5% crystal violet. For CFU-OB assays, cells were cultivated for 7 d in growth medium, followed by 3 d culture in osteogenic differentiation medium followed by alkaline phosphatase staining. For CFU-AD assays, cells were cultured for 8 d in growth medium, followed by 8 d in adipogenic differentiation medium, followed by Oil red O staining. For quantification, the total number of positive colonies was calculated per well. All CFU assays were performed with N = 6 wells per group.

Immunocytochemistry

Request a detailed protocol

Cells were seeded at the density of 1.5 × 105/ml in EZ SLIDE (Merck Millipore, Billerica, MA). After confluence, medium was changed to osteogenic differentiation medium or adipogenic differentiation medium for 3 d. To visualize the membranous expression of CD107a, cells were directly stained with anti-CD107a-allophycocyanin (RRID:AB_1727417, 1:20; BD Pharmingen) for 20 min at 4℃. Then cells were washed with PBS and fixed with 4% PFA, followed by DAPI mounting medium (RRID:AB_2336788, Vector laboratories).

Vacuolin-1 treatment

Request a detailed protocol

To inhibit exocytosis, cells were pre-treated with vacuolin-1 (1 μM; Sigma-Aldrich) for 24 hr before adipogenic differentiation medium was added. Cells were incubated in the presence or absence of vacuolin-1 (1 μM) in adipogenic differentiation medium for 3 d. Surface CD107a expression was detected by flow cytometry or confocal microscopy (Zeiss 800). For immunofluorescence staining, the cell membrane was labeled using Wheat Germ Agglutinin Conjugates (Thermo Fisher Scientific, Waltham, MA).

siRNA knockdown

Request a detailed protocol

In select experiments, siRNA-mediated knockdown of LAMP1 was performed among primary human ASCs prior to osteogenic or adipogenic differentiation. LAMP1 siRNA (Cat# s8080 and s8082) and negative control siRNA (Cat# 4390843) were obtained from Thermo Fisher Scientific. TransIT-LT1 Transfection Reagent (Mirus Bio, Madison, WI) was used as described by the manufacturer. The medium was changed after 4 hr. Validation by qRT-PCR was performed.

Real-time polymerase chain reaction

Request a detailed protocol

Gene expression was analyzed with quantitative real-time polymerase chain reaction (qRT-PCR) (Xu et al., 2019). Total RNA was isolated from cells using TRIzol (Life Technologies). Next, RNA was reverse transcribed into cDNA by iScript cDNA Synthesis Kit (Bio-Rad) following manufacturer's instructions. SYBR Green PCR Master Mix (Life Technologies) was used for RT-PCR. See Supplementary file 8 for primer information. N = 3 wells per group were used, with all studies performed in biologic triplicates.

Western blot

Request a detailed protocol

Proteins were extracted from cultured cells following lysis in ice cold RIPA buffer (Thermo Scientific) with protease inhibitor cocktail (Cell Signaling Technology, Danvers, MA, USA). Proteins were separated by SDS–polyacrylamide gel electrophoresis and transferred onto a nitrocellulose membrane. The blotted nitrocellulose membranes were blocked with 5% bovine serum albumin for 1 hr and then probed with primary antibodies at 4°C overnight. Finally, membranes were incubated with a horseradish-peroxidase (HRP)-conjugated secondary antibody and detected by ChemiDoc XRS+ System (Bio-rad). Quantification of protein bands was performed using Image J software (RRID:SCR_003070).

Intramuscular implantation

Request a detailed protocol

Animals were housed and experiments were performed in accordance with institutional guidelines at Johns Hopkins University under ACUC approval. A DBX putty (courtesy of Musculoskeletal Transplant Foundation, Edison, NJ) was used for ectopic bone formation in mice. Briefly, CD107alow or CD107ahigh cells derived from the same human WAT sample at passage five were prepared at a density of 3 million total cells in 40 μl PBS and mechanically mixed with 45 mg DBX putty. DBX alone was used as an acellular control. The cell preparation was then implanted intramuscularly into the thigh muscle pouch of 8-week-old male NOD-SCID mice (RRID:IMSR_JAX:001303, The Jackson Laboratory, Bar Harbor, ME) as previously described with some modifications (James et al., 2012c). Mice were anesthetized by isoflurane inhalation and premedicated with buprenorphine. Incisions in the hindlimbs were made, and pockets were cut in the biceps femoris muscles by blunt dissection, parallel to the muscle fiber long axis. Dissection methods and the surgical manipulation of tissues were kept as constant as possible across animals. The muscle and skin were each closed with 4–0 Vicryl*Plus sutures (Ethicon Endo-Surgery, Blue Ash, OH). See Supplementary file 5 for an outline of animals per experimental group. Surgical implantations and subsequent analyses were performed blinded to treatment group.

Lumbar spine fusion

Request a detailed protocol

Posterolateral lumbar spinal fusion was performed on 23-week-old male athymic rats (RRID:RGD_2312499) as previously described (Chung et al., 2014). Posterior midline incisions were made over the caudal portion of the lumbar spine, and two separate fascial incisions were made 4 mm bilaterally from the midline. L4 and L5 lumbar spines were exposed by blunt muscle splitting technique and decorticated using a low speed burr and micro-drill (Roboz Surgical Instrument Co., Gaithersburg, MD). Next, DBX (300 μl per side) mixed with CD107alow or CD107ahigh cells from the same AT sample (1.5 million cells, P3-5 passage) or DBX alone were implanted between the transverse processes bilaterally into the paraspinal muscle bed. Finally, the fasciae and skin were each closed using continuous suture (4–0 Vicryl*Plus, Ethicon Endo-Surgery). In vivo imaging was performed by a combination X-ray/DXA (Faxitron Bioptics, Tucson, AZ) at 0, 4, and 8 weeks after surgery. Rats were sacrificed 8 weeks after surgery, and the spines were harvested for further analysis. See Supplementary file 6 for animals per experimental group. Surgical procedures and subsequent analyses were performed blinded to treatment group.

Post mortem analyses

Request a detailed protocol

Samples were fixed in 4% PFA for 24–48 hr and evaluated using a high-resolution micro-computed tomography (μCT) imaging system (SkyScan 1275; Bruker MicroCT N.V, Kontich, Belgium). For intramuscular implants, scans were obtained at an image resolution of 15 μm with a 1 mm of aluminum filter (X-ray voltage of 65 kVP, anode current of 153 uA, exposure time of 218 ms). For spine fusion samples, scans were obtained at an image resolution of 22 μm with a 1 mm of aluminum filter (X-ray voltage of 55 kVP, anode current of 181 uA, exposure time of 218 ms). NRecon software (SkyScan, Bruker) was used to reconstruct images from the 2D X-ray projections. For the 3D morphometric analyses of images, CTVox and CTAn were used (SkyScan, Bruker). For muscle pouch implantation and spine fusion analysis, volumes of interest were shaped to encompass all the implant and exclude native bone, with a threshold value of 65.

After radiographic imaging, samples were transferred to 14% EDTA for decalcification for 28–60 d. Samples were then embedded in OCT for cryosectioning at 18 μm thickness. H and E and ALP staining were performed on serial sections (Leukocyte Alkaline Phosphatase Kit, Sigma-Aldrich). For immunofluorescent staining, sections were washed with PBS three times and blocked with 5% goat serum in PBS for 1 hr at 25° C. Antigen retrieval was performed by trypsin enzymatic antigen retrieval solution (ab970; Abcam) for 10 min at 25° C. Primary anti-osteocalcin antibody (RRID:AB_10675660, 1:100, Abcam) was added to each section and incubated at 4° C overnight. Next, an Alexa Fluor 647 goat anti-rabbit IgG (H+L) polyclonal (RRID:AB_2722623, 1:200, Abcam) was used as the secondary antibody. Sections were counterstained with DAPI mounting medium (RRID:AB_2336788, Vector laboratories).

Statistical analysis

Request a detailed protocol

Quantitative data are expressed as mean ±one SD. Statistical analyses were performed using the SPSS16.0 software (RRID:SCR_002865) or GraphPad Prism (RRID:SCR_002798, Version 7.0). Our in vitro studies comparing CD107alow to CD107ahigh cells resulted in effect size of 3.75. Based on α = 0.05 and power = 0.8, statistical significance should be observed with N = 6 animals per group assuming a one-way ANOVA with a 0.05 significance level. Student’s t-test was used for two-group comparisons, and one-way ANOVA test was used for comparisons of three or more groups, followed by Tukey’s post hoc test. Two-way ANOVA test was used for comparisons of two or three groups with different time points. Differences were considered significant with *p<0.05, **p<0.01, and ***p<0.001.

Appendix 1

Appendix 1—key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional
information
Strain, strain background (mouse)NOD-SCIDJackson LaboratoryStrain #001303, RRID:IMSR_JAX:001303Male, 8 week old
Strain, strain background (rat)RNU Nude RatCharles River Laboratories IncStrain #316, RRID:RGD_2312499Male, 23 week old
Transfected construct (human)Negative Control siRNAThermoFisher ScientificCat#4390843
Transfected construct (human)LAMP1 siRNAThermoFisher ScientificCat#4392420
Assay ID s8082
Transfected construct (human)LAMP1 siRNA 2#ThermoFisher ScientificCat#4392420
Assay ID s8080
Antibodyanti-Human CD31-APC-Cy7 (Mouse monoclonal)BD PharmingenCat# 563653, RRID:AB_2738350FACS/Flow cytometry (1:100)
Antibodyanti-Human CD31 (Mouse monoclonal)AbcamCat# ab24590, RRID:AB_448167Immunofluorescent staining (1:100)
Antibodyanti-Human CD31 (Rabbit polyclonal)AbcamCat# ab28364, RRID:AB_726362Immunofluorescent staining (1:100)
Antibodyanti-Human CD34-PE-CF594 (Mouse monoclonal)BD PharmingenCat# 562383, RRID:AB_11154586Flow cytometry (1:60)
Antibodyanti-Human CD34 (Rabbit monoclonal)AbcamCat# ab81289, RRID:AB_1640331Immunofluorescent staining (1:100)
Antibodyanti-Human CD44-AF700 (Mouse monoclonal)BD PharmingenCat# 561289, RRID:AB_10645788Flow cytometry (1:20)
Antibodyanti-Human CD45-APC-Cy7 (Mouse monoclonal)BD PharmingenCat# 557833, RRID:AB_396891FACS/Flow cytometry (1:30)
Antibodyanti-Human CD73-PE (Mouse monoclonal)BD PharmingenCat# 561014, RRID:AB_2033967Flow cytometry (1:5)
Antibodyanti-Human CD90-FITC (Mouse monoclonal)BD PharmingenCat# 555595, RRID:AB_395969Flow cytometry (1:20)
Antibodyanti-Human CD105-PE-CF594 (Mouse monoclonal)BD PharmingenCat# 562380, RRID:AB_11154054Flow cytometry (1:20)
Antibodyanti-Human CD107a-APC (Mouse monoclonal)BD PharmingenCat# 560664, RRID:AB_1727417FACS/Flow cytometry (1:20)/Immunocytochemistry (1:100)
Antibodyanti-Human CD107a (Mouse monoclonal)R and D SystemsCat# MAB4800, RRID:AB_10719137Immunofluorescent staining (1:100)
Antibodyanti-Human CD107a (Mouse monoclonal)AbcamCat# ab25630, RRID:AB_470708Immunohistochemistry (1:100)/Western blot (1:1000)
Antibodyanti-Human CD107a (Rat monoclonal)AbcamCat# ab25245, RRID:AB_449893Immunofluorescent staining (1:100)
Antibodyanti-Human CD146-FITC (Mouse monoclonal)Bio-RadCat# MCA2141F, RRID:AB_324069Flow cytometry (1:100)
Antibodyanti-Human CD146 (Rabbit monoclonal)AbcamCat# ab75769, RRID:AB_2143375Immunofluorescent staining (1:100)
Antibodyanti-human GAPDH (Rabbit monoclonal)Cell Signaling TechnologyCat# 5174, RRID:AB_10622025Western blot (1:1000)
Antibodyanti-human Gli1 (Rabbit polyoclonal)AbcamCat# ab49314, RRID:AB_880198Immunofluorescent staining (1:100)
Antibodyanti-human Nuclei (Mouse monoclonal)Sigma-AldrichCat# MAB1281, RRID:AB_94090Immunofluorescent staining (1:500)
Antibodyanti-human Osteocalcin (Rabbit polyclonal)AbcamCat# ab93876, RRID:AB_10675660Immunofluorescent staining (1:100)
Antibodyanti-human αSMA (Rabbit polyclonal)AbcamCat# ab21027, RRID:AB_1951138Immunofluorescent staining (1:100)
AntibodyAnti-rabbit IgG, HRP-linked (Goat polyclonal)Cell Signaling TechnologyCat# 7074, RRID:AB_2099233Western blot (1:5000)
AntibodyAnti-mouse IgG, HRP-linked (Horse polyclonal)Cell Signaling TechnologyCat# 7076, RRID:AB_330924Western blot (1:5000)
Antibodyanti-mouse IgG H and L-AF488 (Goat polyclonal)AbcamCat# ab150117, RRID:AB_2688012Immunofluorescent staining (1:200)
Antibodyanti-rabbit IgG H and L-AF488 (Goat polyclonal)AbcamCat# ab150077, RRID:AB_2630356Immunofluorescent staining (1:200)
Antibodyanti-rabbit IgG H+L-DyLight 594 (Goat)Vector LaboratoriesCat# DI-1594, RRID:AB_2336413Immunofluorescent staining (1:200)
Antibodyanti-goat IgG H and L-AF647 (Donkey polyclonal)AbcamCat# ab150135, RRID:AB_2687955Immunofluorescent staining (1:200)
Antibodyanti-mouse IgG H and L-AF647 (Goat polyclonal)AbcamCat# ab150119, RRID:AB_2811129Immunofluorescent staining (1:200)
Antibodyanti-rabbit IgG H and L-AF647 (Goat polyclonal)AbcamCat# ab150079, RRID:AB_2722623Immunofluorescent staining (1:200)
Antibodyanti-rat IgG H and L-AF647 (Goat polyclonal)AbcamCat# ab150167, RRID:AB_2864291Immunofluorescent staining (1:200)
AntibodyImmPRESS-AP Anti-Mouse Reagent antibody (Horse)Vector LaboratoriesCat# MP-5402, RRID:AB_2336535Immunohistochemistry (200 μl)
Sequence-based reagentACAN_FThis paperPCR primersAGGCTGGGGAGAGAACTGAAAAG
Sequenced-based reagentACAN_RThis paperPCR primersGCTCACAATGGGGTATCTGACAG
Sequenced-based reagentACTB_FThis paperPCR primersCTGGAACGGTGAAGGTGACA
Sequenced-based reagentACTB_RThis paperPCR primersAAGGGACTTCCTGTAACAATGCA
Sequenced-based reagentALPL_FThis paperPCR primersACCACCACGAGAGTGAACCA
Sequenced-based reagentALPL_RThis paperPCR primersCGTTGTCTGAGTACCAGTCCC
Sequenced-based reagentCOL2A1_FThis paperPCR primersCCGCGGTGAGCCATGATTCG
Sequenced-based reagentCOL2A1_RThis paperPCR primersCAGGCCCAGGAGGTCCTTTGGG
Sequenced-based reagentCOMP_FPMID:23382851PCR primersCAACTGTCCCCAGAAGAGCAA
Sequenced-based reagentCOMP_RPMID:23382851PCR primersTGGTAGCCAAAGATGAAGCCC
Sequenced-based reagentFABP4_FThis paperPCR primersACGAGAGGATGATAAACTGGTGG
Sequenced-based reagentFABP4_RThis paperPCR primersGCGAACTTCAGTCCAGGTCAAC
Sequenced-based reagentGAPDH_FPMID:31482845PCR primersCTGGGCTACACTGAGCACC
Sequenced-based reagentGAPDH_RPMID:31482845PCR primersAAGTGGTCGTTGAGGGCAATG
Sequenced-based reagentLAMP1_FThis paperPCR primersGTCTTCTTCGTGCCGGCGT
Sequenced-based reagentLAMP1_RThis paperPCR primersGCAGGTCAAAGGTCATGTTCTT
Sequenced-based reagentLPL_FThis paperPCR primersTTGCAGAGAGAGGACTCGGA
Sequenced-based reagentLPL_RThis paperPCR primersGGAGTTGCACCTGTATGCCT
Sequenced-based reagentSPP1_FThis paperPCR primersCCTCCTAGGCATCACCTGTG
Sequenced-based reagentSPP1_RThis paperPCR primersCCACACTATCACCTCGGCC
Sequenced-based reagentRUNX2_FPMID:31482845PCR primersTGGTTACTGTCATGGCGGGTA
Sequenced-based reagentRUNX2_RPMID:31482845PCR primersTCTCAGATCGTTGAACCTTGCTA
Sequenced-based reagentPPARG_FThis paperPCR primersGACAGGAAAGACAACAGACAAATC
Sequenced-based reagentPPARG_RThis paperPCR primersGGGGTGATGTGTTTGAACTTG
Sequenced-based reagentRUNX2_FThis paperPCR primersTGGTTACTGTCATGGCGGGTA
Sequenced-based reagentRUNX2_RThis paperPCR primersTCTCAGATCGTTGAACCTTGCTA
Sequenced-based reagentSOX9_FThis paperPCR primersGAGGAAGTCGGTGAAGAACG
Sequenced-based reagentSOX9_RThis paperPCR primersATCGAAGGTCTCGATGTTGG
Commercial assay or kitLeukocyte Alkaline Phosphatase KitSigma-AldrichCat#85L2-1KT
Chemical compound, drugVacuolin-1Sigma-AldrichCat#673000
Chemical compound, drugTransIT-LT1 Transfection ReagentMirus BioCat#MIR2300Transfection
Software, algorithmFlowJoFlowJoRRID:SCR_008520
Software, algorithmImageJNIHRRID:SCR_003070
Software, algorithmPrismGraphPadRRID:SCR_002798
Software, algorithmSeuratSeuratRRID:SCR_016341
Software, algorithmMonocle 3MonocleRRID:SCR_018685
Software, algorithmPartek Genomics SuitePartekRRID:SCR_011860
Software, algorithmSpotfireSpotfireRRID:SCR_008858
Software, algorithmIngenuity Pathway AnalysisQIAGENRRID:SCR_008653
Software, algorithmSPSSSPSSRRID:SCR_002865
OtherDAPI stainVector LaboratoriesCat# H-1500, RRID:AB_2336788
OtherVector Red Substrate KitVector LaboratoriesCat# SK-5100, RRID:AB_2336847

Data availability

Expression data that support the findings of this study have been deposited in Gene Expression Omnibus (GEO) with the accession codes GSE148519 and GSE128889 (GSM3717979, GSM3717977).

The following data sets were generated
    1. James AW
    2. Xu J
    (2020) NCBI Gene Expression Omnibus
    ID GSE148519. Expression data of CD107aLow and CD107aHigh cells isolated from human adipose tissue.
The following previously published data sets were used
    1. Seale P
    2. Merrick D
    3. Sakers A
    (2019) NCBI Gene Expression Omnibus
    ID GSE128889. Identification of a mesenchymal progenitor cell hierarchy in adipose tissue.

References

    1. Carlsson SR
    2. Roth J
    3. Piller F
    4. Fukuda M
    (1988)
    Isolation and characterization of human lysosomal membrane glycoproteins, h-lamp-1 and h-lamp-2. major sialoglycoproteins carrying polylactosaminoglycan
    The Journal of Biological Chemistry 263:18911–18919.

Decision letter

  1. Clifford J Rosen
    Senior Editor; Maine Medical Center Research Institute, United States
  2. Subburaman Mohan
    Reviewing Editor; Loma Linda University, United States
  3. Yuji Mishina
    Reviewer; University of Michigan School of Dentistry, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Acceptance summary:

The reviewers and the editors are in general agreement that this is a well-designed and executed study that evaluated the usefulness of CD170a as a surrogate marker to distinguish between osteoblast precursors and adipocytes within the human perivascular tissue. The authors have clearly demonstrated the osteogenic potential of CD107alow expressing cells, both in vitro and in vivo. Furthermore, they have carefully and thoroughly addressed the reviewer suggestions and made appropriate amendments in the revised paper. The precise phenotypic and functional characterization of CD107a perivascular progenitor cells will provide a foundation for future understanding of the role of these cells in developmental, regeneration and pathological process.

Decision letter after peer review:

Thank you for submitting your article "Lysosomal protein surface expression discriminates fat- from bone-forming human mesenchymal precursor cells" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Clifford Rosen as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Yuji Mishina (Reviewer #2).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

As the editors have judged that your manuscript is of interest, but as described below that additional experiments are required before it is published, we would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is "in revision at eLife". Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)

Summary:

In this study James et al. have performed extensive analysis of CD107a expression in ASCs and determined that CD107alow expressing cells exhibit osteogenic potential both in vitro and in vivo. While all three reviewers are in agreement that the manuscript is well written and the establishment that CD107a serves as a surrogate marker to distinguish between adipogenic and osteogenic human ASCs would provide an important contribution to the field, they also have raised a number of questions (see below) that relate to the role of CD107a in cell fate. The reviewers and the editors realize the difficulty in performing the suggested additional in vivo bone regeneration experiments in a timely manner due to current Covid19 situation. However, the authors could use the in vitro model system to complement in vivo data to provide additional data to support the claims made in the paper. For example, the CD107a knockdown experiments using purified populations would provide more convincing data for direct involvement of CD107a in cell fate specification (or just an innocent bystander). Please consider clarifying some aspects noted below that were raised in the reviewer comments without additional experiments.

Please let us know if you have any questions.

Reviewer #1:

In this study, Peault B et al. claim that the cell surface expression of lysosomal associated protein-1 (LAMP1 or CD107a) serves as a surrogate marker to distinguish between osteogenic and adipogenic human adipose tissue derived stem cells (ASCs).

Using BD Lyoplate assay, authors identified CD107a as a cell surface protein differentially expressed between pericytes (CD146+CD34-) and adventicytes (CD146-CD34+). Authors recently demonstrated that these cells work together in bone formation (Wang et al., 2019). In this paper, authors show that CD107alow ASC are more osteogenic relative to CD107high cells, whereas the latter is more adipogenic than the former in vitro. in vitro experiments further demonstrate that CD107 cell surface expression is induced during adipogenesis. Gene silencing of CD107, however, does not interfere with adipogenic potential of ASCs. Consistent with in vitro findings, CD107ahigh ASCs demonstrate highly osteogenic potential in ectopic bone formation and spinal fusion mouse models.

The potential significance of CD107a as a surrogate marker to enrich adipogenic ASCs is shown in vitro and the functional impact of CD107alow ASCs in bone formation in vivo was well documented in this paper. However, the biological significance of lysosomal trafficking in ASCs is unclear.

1) Increased CD107a expression on the cell surface may simply reflect increased exocytosis process in differentiating preadipocytes. A subset of ASCs might be resistant to adipogenesis, which results in reduced CD107a expression on cell surface. From this study alone we do not know whether the intrinsic ability of lysosomal trafficking confers adipogenic versus osteogenic capacity of ASCs. Addressing the functional significance of lysosomal trafficking in adipogenesis may help us better understand the cellular mechanism of human ASC differentiation.

2) In vivo experiment of ectopic bone formation or bone fusion did not compare CD107alow cells versus unsorted ASCs. If the focus of this paper is to use CD107a to enrich osteogenic ASCs for bone regeneration, comparing between enriched and unsorted cell populations will be needed.

3) Adipogenic potential of CD107ahigh ASCs in vivo is not shown; as such, the in vivo significance of CD107ahigh cells is unclear.

4) As discussed, developmental and pathological significance of CD107a in ASCs is unclear unless authors demonstrate the cell surface expression of this lysosomal protein in development and certain disease process, such as obesity.

5) Increased expression of CD107a in pericytes (CD146+) relative to adventicytes (CD34+) is intriguing. As the roles of pericytes versus fibroblasts in adipogenesis during development and in obesity are actively investigated at least in mice (1, 2), comparing the adipogenic potential of human pericytes relative to adventicytes may help us better understand human adipose tissue biology.

References:

1) Hepler C, Shan B, Zhang Q, Henry GH, Shao M, Vishvanath L, Ghaben AL, Mobley AB, Strand D, Hon GC, Gupta RK. Identification of functionally distinct fibro-inflammatory and adipogenic stromal subpopulations in visceral adipose tissue of adult mice. eLife. 2018;7:e39636. doi: 10.7554/eLife.39636.

2) Cattaneo P, Mukherjee D, Spinozzi S, Zhang L, Larcher V, Stallcup WB, Kataoka H, Chen J, Dimmeler S, Evans SM, Guimarães-Camboa N. Parallel Lineage-Tracing Studies Establish Fibroblasts as the Prevailing in vivo Adipocyte Progenitor. Cell Reports. 2020;30(2):571-82.e2. doi: https://doi.org/10.1016/j.celrep.2019.12.046.

Reviewer #2:

The authors identified endolysosomal protein CD107a as a novel cell surface antigen expressed among adipose tissue (AT)-resident perivascular stem cells. CD107alow and CD107ahigh AT-derived stromal cells represent osteoblast and adipocyte precursor cells, respectively. They also confirmed human CD107alow cells drive dramatic bone formation using two mouse models. They demonstrated that CD107a correlates with exocytosis during early adipogenic differentiation, rather than having vital function in cellular differentiation. They further systematically compared the transcriptome of both cell types and identified CD107alow cells are precursors of CD107ahigh cells using bioinformatics analysis.

The authors did extensive and systemic analysis to compare CD107alow and CD107ahigh cell populations. The paper is well written and the identification of this new stem cell marker is an important contribution to the field.

This reviewer does not have any major issues.

Reviewer #3:

This is an interesting paper on the identification of CD107a as a novel marker of MSC to determine a subpopulation of the progenitors with specific differentiation and biological function in human stem cell biology. However, I have several comments that will require further attention:

1) In the Introduction, the authors mentioned the basic function of CD107a (Lamp1) marker as a lysosomal membrane protein. However, they did not mention its important role in autophagy, which should be included as it can make relevant explanation for their observation as the process of autophagy is crucial for regulation of adipogenesis (DOI: 10.1128/MCB.00193-13). The authors should add this information as it can be further discussed in relation to their findings.

2) Regarding the expression profile of CD107a in human adipose tissue, did the authors check the proportion of the cells and cell types with CD107a expression to identify the heterogeneity of the cells expressing this surface marker in AT and also if it correlates with number of lysosomes or their state of the activation?

3) Further, the authors mentioned that CD107a distinguish between AD and OB progenitors among MSC derived from AT. Did the authors check the expression profile of CD107a during AD and OB differentiation if its expression profile changes during these processes? They should add these data and comment on it in relation to CD107a function in MSC. Also, if there are known other factors e.g. inflammatory/metabolic molecules (e.g. TNFa, TGFb, or lipids) that could stimulate expression of CD107a?

4) Regarding the source of human AT samples, the authors should add the information on the age and BMI or possible medication used in the subjects as these are the major factors that can affect the quality and cellular composition of the tissue.

5) For the cell sorting, how stable was the expression of CD107a after sorting in the culture, did the CD107ahigh and low cells maintain their difference in CD107a expression through the passaging? It is very common that MSC after sorting change their expression profile through the sub-culturing. The authors should add these data and comment on them.

6) Did the authors measure proliferation rate of CD107ahigh and low expressed cells? Did they measure ALP activity or inflammatory properties (e.g. expression of cytokines, growth factors) in these cells? These are functional assays that would add more information about their cellular properties and determination.

7) Regarding the CFU-AD/OB data, how did the authors evaluate this parameter, i.e. they calculated only positive colonies or the ratio of positive colonies to total number of colonies? This should be specified to know if there was change in determination or switch in differentiation potential of the cells.

8) Regarding the different source of MSC on CD107a expression, the authors mentioned some data on CD107a expression in different fat depots or placenta (Discussion, fourth paragraph). They should add this information, also if they could add the data on BMSC and their expression of CD107a as these cells are more prone to osteogenesis. It will be interesting to compare different sources of MSC and the CD107a expression if it can be used for the specification of particular cell subpopulation.

9) Regarding the RNA-seq data between CD107alow and high (Figure 4C), is the expression of AD genes significant between groups? It does not seem to be much different, e.g. CEBPb, PPARg, Sirt1 indeed seem to be more expressed in CD107alow. The authors should more carefully interpret these data and correct their postulation.

10) Regarding the major outcome of this study, it is not clear what the authors would like to point out with their findings. Does it mean that MSC with lower CD107a expression have higher multipotent characteristics and better regeneration properties? They should elaborate more on these aspects of their findings and discuss it more with other literature. Also, if the testing regenerative properties of CD107alow and high cells should be attributed, the authors should use better model for tissue regeneration for example calvarial defect regeneration model or monocortical bone defect model. As these models represent more appropriate system with injury stimuli for activation of the regenerative properties of the cells. The authors should work more on this part of the results.

https://doi.org/10.7554/eLife.58990.sa1

Author response

Summary:

In this study James et al. have performed extensive analysis of CD107a expression in ASCs and determined that CD107alow expressing cells exhibit osteogenic potential both in vitro and in vivo. While all three reviewers are in agreement that the manuscript is well written and the establishment that CD107a serves as a surrogate marker to distinguish between adipogenic and osteogenic human ASCs would provide an important contribution to the field, they also have raised a number of questions (see below) that relate to the role of CD107a in cell fate. The reviewers and the editors realize the difficulty in performing the suggested additional in vivo bone regeneration experiments in a timely manner due to current Covid19 situation. However, the authors could use the in vitro model system to complement in vivo data to provide additional data to support the claims made in the paper. For example, the CD107a knockdown experiments using purified populations would provide more convincing data for direct involvement of CD107a in cell fate specification (or just an innocent bystander). Please consider clarifying some aspects noted below that were raised in the reviewer comments without additional experiments.

We thank the editor for this suggestion. To further pursue this, we examined CD107a knockdown in FACS isolated CD107alow and CD107ahigh cells. Our results showed that knockdown of CD107a in CD107alow cells slightly promoted adipogenic differentiation, which was in similarity to unsorted cells. In CD107ahigh cells, knockdown had no effect (revised Figure 3—figure supplement 4P, Q). We speculate that as CD107ahigh cells have high intrinsic adipogenic potential, the effects of knockdown were not apparent in this group. In aggregate, these data suggest that CD107a is a marker (or innocent bystander), rather than directly involved in cell fate specification.

Reviewer #1:

In this study, Peault B et al. claim that the cell surface expression of lysosomal associated protein-1 (LAMP1 or CD107a) serves as a surrogate marker to distinguish between osteogenic and adipogenic human adipose tissue derived stem cells (ASCs).

Using BD Lyoplate assay, authors identified CD107a as a cell surface protein differentially expressed between pericytes (CD146+CD34-) and adventicytes (CD146-CD34+). Authors recently demonstrated that these cells work together in bone formation (1). In this paper, authors show that CD107alow ASC are more osteogenic relative to CD107high cells, whereas the latter is more adipogenic than the former in vitro. in vitro experiments further demonstrate that CD107 cell surface expression is induced during adipogenesis. Gene silencing of CD107, however, does not interfere with adipogenic potential of ASCs. Consistent with in vitro findings, CD107ahigh ASCs demonstrate highly osteogenic potential in ectopic bone formation and spinal fusion mouse models.

The potential significance of CD107a as a surrogate marker to enrich adipogenic ASCs is shown in vitro and the functional impact of CD107alow ASCs in bone formation in vivo was well documented in this paper. However, the biological significance of lysosomal trafficking in ASCs is unclear.

1) Increased CD107a expression on the cell surface may simply reflect increased exocytosis process in differentiating preadipocytes. A subset of ASCs might be resistant to adipogenesis, which results in reduced CD107a expression on cell surface. From this study alone we do not know whether the intrinsic ability of lysosomal trafficking confers adipogenic versus osteogenic capacity of ASCs. Addressing the functional significance of lysosomal trafficking in adipogenesis may help us better understand the cellular mechanism of human ASC differentiation.

We thank the reviewer for this suggestion. As you point out, it is likely that increased CD107a expression reflects increased exocytosis process in differentiating preadipocytes. In our past observations, CD107a knockdown led to a modest but significant increase in adipogenesis. To further expand on these results, adipogenic differentiation of ASCs was performed with Vac-1 treatment, the inhibitor of exocytosis which we previously observed prevented CD107a translocation to the cell membrane. Consistent with siRNA knockdown results, Vac-1 treatment also led to a modest increase in adipogenic differentiation (see Author response image 1). These results, along with new siRNA experiments using CD107alow and high cells (Figure 3—figure supplement 4P, Q), suggest that partitioning allows for identification of functionally relevant subsets which is likely not explained by intrinsic function of CD107a protein. Changes to the Discussion to make this point more clear have been added (second paragraph).

Author response image 1
Inhibition of exocytosis with Vacuolin-1 moderately promotes adipogenic differentiation of human ASCs.

(A) Oil red O (ORO) staining (left) and photometric quantification (right) at d7 of adipogenic differentiation with or without Vacuolin-1 treatment (Vac-1, 1 μM). Black scale bar: 50 μm. (B-D) Adipogenic gene expression with or without Vac-1 (1 μM) treatment, including (B) Peroxisome proliferator-activated receptor-γ (PPARG), (C) Lipoprotein lipase (LPL), and (D) Fatty acid binding protein 4 (FABP4). Gene expression assessed at 7 d of adipogenic differentiation. Data repeated in experimental triplicate. *P<0.05; **P<0.01. Statistical analysis was performed using a two-tailed Student t-test.

2) In vivo experiment of ectopic bone formation or bone fusion did not compare CD107alow cells versus unsorted ASCs. If the focus of this paper is to use CD107a to enrich osteogenic ASCs for bone regeneration, comparing between enriched and unsorted cell populations will be needed.

We thank the reviewer for the excellent suggestion. The main focus of the paper was to better understand the cellular heterogeneity of human progenitor cells within the perivascular niche. As pointed out by the editor, timely additional in vivo experiments are not currently feasible due to COVID19 related shutdowns in our area. However, we agree that the future use of CD107a cell subsets for tissue engineering efforts in the bone or fat field would require clear evidence of improvement in comparison to unsorted ASCs.

3) Adipogenic potential of CD107ahigh ASCs in vivo is not shown; as such, the in vivo significance of CD107ahigh cells is unclear.

Thank you for the excellent point. From our data and across 8 patient samples, we did find that CD107ahigh stromal cells represent adipocyte precursor cells in vitro. Just as perivascular cells in PDGFRA1 or PDGFRB2 reporter animals have been observed to replenish adipocytes among either homeostasis or high fat diet conditions, it is quite possible that CD107a-expressing perivascular cells perform a similar function in vivo. Our studies are confined to human tissues, and future studies examining the dynamic cell fate of CD107ahigh cells would need to be performed in the mouse in the context of lineage tracing. Of note, we did return to our muscle pouch and spine fusion assays and performed Oil red O staining to look for frank adipocytes within these tissues. No definite adipocytes were observed across any samples, suggesting that these particular carriers are not permissive to adipocyte formation. Future tissue engineering studies with CD107ahigh cells for adipose tissue formation are indeed needed.

4) As discussed, developmental and pathological significance of CD107a in ASCs is unclear unless authors demonstrate the cell surface expression of this lysosomal protein in development and certain disease process, such as obesity.

Thank you for the good point. Mice deficient in either CD107a or CD107b are viable and fertile. However, mice deficient in both CD107a and CD107b have an embryonic lethal phenotype. In mouse embryonic fibroblasts, mutual disruption of both CD107s is associated with an increased accumulation of autophagic vacuoles and unesterified cholesterol3. In our preliminary data, CD107b is also expressed in the human adipose stromal cells. However, using FACS we found that only a small fraction of cells (0.39% in total ASCs) co-expressed surface CD107a and CD107b, making functional redundancy of these proteins in ASCs less likely. Any link between CD107a/CD107b and either obesity or a skeletal phenotype in mouse or humans is not known. In pilot studies, we did examine CD107a expression in a tissue microarray of healthy and diseased human coronary arteries. Although still preliminary data, adventitial and plaque associated CD107a immunostaining was present, particularly in diseased vessels (Author response image 2). It remains to be seen the extent to which CD107a-expressing adventitial cells are involved in pathologic vascular remodeling.

Author response image 2
CD107a immunohistochemical staining in human atherosclerotic vessel.

CD107a positive staining appears red. Red scale bar: 200 μm; black scale bar: 50 μm.

5) Increased expression of CD107a in pericytes (CD146+) relative to adventicytes (CD34+) is intriguing. As the roles of pericytes versus fibroblasts in adipogenesis during development and in obesity are actively investigated at least in mice (1, 2), comparing the adipogenic potential of human pericytes relative to adventicytes may help us better understand human adipose tissue biology.

Thank you for the interesting point. In our experience, adventitial cells rather than pericytes have a higher intrinsic adipogenic potential, which most likely correlates with their theoretically higher ‘stem-like’ identity4. In terms of frequency, and across 8 patient samples, membranous CD107a was present in CD146+ pericytes, CD34+ adventicytes and also a fraction of as of yet undefined mesenchymal cells (C146-CD34-CD31-CD45- cells). Interestingly, the percentage of pericytes and adventicytes that made up the CD107ahigh cell population varied to some degree across samples (see Supplementary files 2 and 3). In the future, membranous CD107a among mesenchymal cells or a subfraction of perivascular mesenchymal cells could be used prospectively as a biomarker. Nevertheless, our study was a first step in analyzing a completely unknown marker in MSC biology. The further functional segregation of CD107a-expressing cell subsets (pericytes and adventitial cells) is a fascinating area of research, but lies considerably outside the scope of the present work.

References:

1) Hepler C, Shan B, Zhang Q, Henry GH, Shao M, Vishvanath L, Ghaben AL, Mobley AB, Strand D, Hon GC, Gupta RK. Identification of functionally distinct fibro-inflammatory and adipogenic stromal subpopulations in visceral adipose tissue of adult mice. eLife. 2018;7:e39636. doi: 10.7554/eLife.39636.

2) Cattaneo P, Mukherjee D, Spinozzi S, Zhang L, Larcher V, Stallcup WB, Kataoka H, Chen J, Dimmeler S, Evans SM, Guimarães-Camboa N. Parallel Lineage-Tracing Studies Establish Fibroblasts as the Prevailing in vivo Adipocyte Progenitor. Cell Reports. 2020;30(2):571-82.e2. doi: https://doi.org/10.1016/j.celrep.2019.12.046.

Reviewer #3:

This is an interesting paper on the identification of CD107a as a novel marker of MSC to determine a subpopulation of the progenitors with specific differentiation and biological function in human stem cell biology. However, I have several comments that will require further attention:

1) In the Introduction, the authors mentioned the basic function of CD107a (Lamp1) marker as a lysosomal membrane protein. However, they did not mention its important role in autophagy, which should be included as it can make relevant explanation for their observation as the process of autophagy is crucial for regulation of adipogenesis (DOI: 10.1128/MCB.00193-13). The authors should add this information as it can be further discussed in relation to their findings.

Thank you for the excellent point, which we have added to the Discussion (second paragraph). CD107a is responsible for maintaining lysosomal integrity, pH and catabolism3,5. Lysosomes provide the degradative enzymes for autophagy and participate in autophagy regulation6, and CD107a may have an essential role in the process of autophagy. Also as you point out, the activation of autophagy is crucial for adipogenic differentiation7. Whether CD107ahigh cells enhance adipogenesis by promoting autophagy is an interesting and worthy follow-up topic. This additional discussion has been added to the resubmission (see the aforementioned paragraph).

2) Regarding the expression profile of CD107a in human adipose tissue, did the authors check the proportion of the cells and cell types with CD107a expression to identify the heterogeneity of the cells expressing this surface marker in AT and also if it correlates with number of lysosomes or their state of the activation?

Thank you for the excellent suggestion. The proportion of cell types present within CD107ahigh and low cell preparations is reported in Supplementary files 2 and 3. Briefly, both CD107ahigh and low cell populations showed relatively equivalent numbers of CD146+ pericytes (6.94% in CD107alow cells and 4.47% in CD107ahigh cells), while CD107alow cells contained a slightly greater population of CD34+ adventitial cells (70.88% in CD107alow cells and 52.45% in CD107ahigh cells). According to your suggestion, we quantified the number of lysosomes among the CD107alow and CD107ahigh cells using Lysosomal Staining Kit (Abcam) and found no difference in lysosome staining intensity (see Author response image 3).

Author response image 3
No change in lysosome immunohistochemical staining among human CD107alow and CD107ahigh cells.

Lysosome positive staining appears green. White scale bar: 20 μm.

3) Further, the authors mentioned that CD107a distinguish between AD and OB progenitors among MSC derived from AT. Did the authors check the expression profile of CD107a during AD and OB differentiation if its expression profile changes during these processes? They should add these data and comment on it in relation to CD107a function in MSC. Also, if there are known other factors e.g. inflammatory/metabolic molecules (e.g. TNFa, TGFb, or lipids) that could stimulate expression of CD107a?

Thank you for the excellent suggestions. We observed that CD107a was increased prominently during early adipogenic differentiation (Figure 3H-J). Briefly, after 3 d exposure to adipogenic conditions, a 12.03-fold increase in immunostaining intensity and a 253.5% increase in the number of CD107ahigh cells were noted. In contrast, no significant change in membranous CD107a expression profile was found during osteogenic differentiation (Figure 3—figure supplement 2). As to your good point regarding inflammatory mediators of CD107a expression, this is well reported in other cell types. For example, IL-2 stimulates CD107a expression on NK and cytotoxic T cells8. In addition, CD107a production is increased with the stimulation of IL-15, IFN-α, and IFN-β in NK cells9. Linking perivascular inflammation to mesenchymal CD107a expression in obesity or other pathophysiology would be an interesting subject for the future.

4) Regarding the source of human AT samples, the authors should add the information on the age and BMI or possible medication used in the subjects as these are the major factors that can affect the quality and cellular composition of the tissue.

Thank you for the suggestions. Details of age and BMI have been added in Supplementary file 1. Unfortunately, data on medication usage was not available for these de-identified samples.

5) For the cell sorting, how stable was the expression of CD107a after sorting in the culture, did the CD107ahigh and low cells maintain their difference in CD107a expression through the passaging? It is very common that MSC after sorting change their expression profile through the sub-culturing. The authors should add these data and comment on them.

Thank you for the good point. The frequency of membranous CD107a in CD107ahigh cells did not change after three passages. However, the frequency of membranous CD107a among passaged CD107alow cells increased slightly over three passages (1.41%).

6) Did the authors measure proliferation rate of CD107ahigh and low expressed cells? Did they measure ALP activity or inflammatory properties (e.g. expression of cytokines, growth factors) in these cells? These are functional assays that would add more information about their cellular properties and determination.

Thank you for the suggestion. Proliferation assays have been added, in which we observed that CD107alow cells have a higher proliferative index in comparison to CD107ahigh cells (Figure 2A, subsection “CD107alow AT-derived stromal cells represent osteoblast precursor cells!). As suggested, we observed that ALP activity was quite different among CD107alow and high cells (Figure 2I, J), in which ALP staining and quantification demonstrated an enrichment among CD107alow cells. According to your suggestions, we re-analyzed existing RNA sequencing data and found that CD107alow and high cells expressed few inflammatory cytokines, such as IL-1B, IL-6, and IL-33. In addition, both cell populations produced common growth factors such as FGF2, BMP2, PDGFA, PDGFB, VEGFA and IGF1, and at comparable levels (see Author response image 4).

Author response image 4
Expression of inflammatory cytokines and common growth and differentiation factors among human CD107alowCD31-CD45- and CD107ahighCD31-CD45- mesenchymal cells.

Data shown in volcano plot. X-axis represents Log2 fold change for each gene. Y-axis represents -Log10 p value. Red dots indicate >2SD increase among CD107ahigh mesenchymal cells. Blue dots indicate >2SD increase among CD107alow mesenchymal cells. Inflammatory cytokines are shown in purple box, while growth factors are shown in green box.

7) Regarding the CFU-AD/OB data, how did the authors evaluate this parameter, i.e. they calculated only positive colonies or the ratio of positive colonies to total number of colonies? This should be specified to know if there was change in determination or switch in differentiation potential of the cells.

Thank you for the good point. For CFU-OB/AD quantification, we calculated the total number of positive colonies in each well. This has been added in the Materials and methods section (subsection “Colony-forming unit (CFU) assay”).

8) Regarding the different source of MSC on CD107a expression, the authors mentioned some data on CD107a expression in different fat depots or placenta (Discussion, fourth paragraph). They should add this information, also if they could add the data on BMSC and their expression of CD107a as these cells are more prone to osteogenesis. It will be interesting to compare different sources of MSC and the CD107a expression if it can be used for the specification of particular cell subpopulation.

Thank you for the suggestion. CD107a expression within different fat depots is shown in Figure 1—figure supplement 2, in which CD107a immunoreactivity within the perivascular mesenchymal niche was noted, including pericardial, perigonadal, perirenal and omental human fat. In addition, we have begun to investigate disparate tissues such as placenta and fetal skeletal muscle. An example of the perivascular staining patterns are shown in Author response image 5. However, this distinct project is being performed by another graduate student, and is reserved for a distinct thesis project. Finally, we did observe that a portion of BMSC expressed membranous CD107a (2.54%, Figure 3L). Interestingly, BMSC showed a similar increase in CD107a expression during adipogenesis as did ASC (Figure 3L). We hope in future work to more thoroughly understand how CD107a identifies perivascular cell subsets across and between tissues.

Author response image 5
CD107a immunohistochemical staining in human fetal skeletal muscle.

CD107a positive staining appears red, while CD31 positive staining appears black.

9) Regarding the RNA-seq data between CD107alow and high (Figure 4C), is the expression of AD genes significant between groups? It does not seem to be much different, e.g. CEBPb, PPARg, Sirt1 indeed seem to be more expressed in CD107alow. The authors should more carefully interpret these data and correct their postulation.

Thank you for the suggestion. LRP5 (3.76-fold increase; p value 0.0028), FABP4 (1.40-fold increase; p value 0.0160), NCOR2 (4.69-fold increase; p value 0.0060), and ANGPT2 (2.26-fold increase; p value 0.0220) were significantly increased in CD107ahigh stromal cells. CEBPA (1.08-fold increase; p value 0.8488), CEBPD (1.48-fold increase; p value 0.0873), and LPL (1.11-fold increase; p value 0.3596) trended towards an upregulation among in CD107ahigh cells, but this did not reach statistical significance. Some significantly upregulated genes in CD107alow stromal cells were negative regulators of adipogenesis, such as KLF2, KLF3, SIRT1, and DDIT310-13 (subsection “Transcriptomic analysis suggests a progenitor cell phenotype for CD107alow cells”, first paragraph).

10) Regarding the major outcome of this study, it is not clear what the authors would like to point out with their findings. Does it mean that MSC with lower CD107a expression have higher multipotent characteristics and better regeneration properties? They should elaborate more on these aspects of their findings and discuss it more with other literature. Also, if the testing regenerative properties of CD107alow and high cells should be attributed, the authors should use better model for tissue regeneration for example calvarial defect regeneration model or monocortical bone defect model. As these models represent more appropriate system with injury stimuli for activation of the regenerative properties of the cells. The authors should work more on this part of the results.

Thank you for the excellent suggestions and the opportunity to summarize. Our main findings were that CD107alow cells demonstrate high colony formation, osteoprogenitor cell frequency, and osteogenic potential in vitro and in vivo. Conversely, CD107ahigh cells include almost exclusively adipocyte progenitor cells and exhibit high fat cell forming potential. Overall this suggests a more ‘stem-like’ identity for CD107alow cells, and that CD107ahigh cell populations have a high pre-adipocyte content which is functionally associated with exocytosis during early adipogenic differentiation.

According to your earlier suggestion, we have reviewed in more detail the literature regarding CD107a regulation of autophagy (Discussion, second paragraph), as well as inflammatory inducers of CD107a expression (see response to point 3). As you point out, there are many potential bone injury models that could be chosen, including monocortical bone defects or calvarial bone defects. In our experience with the mentioned models, these models rely on both paracrine stimulation of endogenous cells and direct ossification of the implanted cells. In contrast, the spine fusion model relies more on direct ossification of the implanted cells, which was one of the primary factors which influenced our choice of models.

References:

1) Wang Y, Xu J, Meyers CA, et al. PDGFRalpha marks distinct perivascular populations with different osteogenic potential within adipose tissue. Stem Cells 38, 276-290, doi:10.1002/stem.3108 (2020).

2) Gao Z, Daquinag, AC, Su F, Snyder B and Kolonin MG. PDGFRalpha/PDGFRbeta signaling balance modulates progenitor cell differentiation into white and beige adipocytes. Development 145, doi:10.1242/dev.155861 (2018).

3) Eeva-Liisa E. Roles of LAMP-1 and LAMP-2 in lysosome biogenesis and autophagy. Mol Aspects Med 27(5-6), 495-502, doi: 10.1016/j.mam.2006.08.005 (2006).

4) Hardy WR, Moldovan N, Moldovan L, et al. Transcriptional Networks in Single Perivascular Cells Sorted from Human Adipose Tissue Reveal a Hierarchy of Mesenchymal Stem Cells. Stem Cells 35(5), 1273-1289, doi: 10.1002/stem.2599 (2017).

5) Andrejewski N, Punnonen EL, Guhde G, et al. Normal lysosomal morphology and function in LAMP-1-deficient mice. The Journal of Biological Chemistry 274(18), 12692–12701, doi: 10.1074/jbc.274.18.12692 (1999).

6) Yim WW, Mizushima N. Lysosome biology in autophagy. Cell Discovery 6, doi: 10.1038/s41421-020-0141-7 (2020).

7). Guo L, Huang J, Liu Y, et al. Transactivation of Atg4b by C/EBPβ promotes autophagy to facilitate adipogenesis. Mol Cell Biol 33(16), 3180-3190, doi: 10.1128/MCB.00193-13 (2013).

8) Aktas E, Kucuksezer UC, Bilgic S, et al. Relationship between CD107a expression and cytotoxic activity. Cellular Immunology 254(2), 149-154, doi: 10.1016/j.cellimm.2008.08.007 (2009).

9) Xie Z, Zheng J, Wang Y, et al. Deficient IL-2 produced by activated CD56+ T cells contributes to impaired NK cell-mediated ADCC function in chronic HIV-1 infection. Front Immunol 10, 1647, doi: 10.3389/fimmu.2019.01647 (2019).

10) Banerjee S, Feinberg M, Watanabe M, et al. The Krüppel-like factor KLF2 inhibits peroxisome proliferator-activated receptor-γ expression and adipogenesis. J Biol Chem 278(4), 2581-2584, doi: 10.1074/jbc.M210859200 (2003).

11) Sue N, Jack B, Eaton S, et al. Targeted disruption of the basic Krüppel-like factor gene (Klf3) reveals a role in adipogenesis. Mol Cell Biol 28(12), 3967-3978, doi: 10.1128/MCB.01942-07 (2008).

12) Zhou Y, Zhou Z, Zhang W, et al. SIRT1 inhibits adipogenesis and promotes myogenic differentiation in C3H10T1/2 pluripotent cells by regulating Wnt signaling. Cell Biosci 5, 61, doi: 10.1186/s13578-015-0055-5 (2015).

13) Pereira R, Delany A, Canalis E. CCAAT/enhancer binding protein homologous protein (DDIT3) induces osteoblastic cell differentiation. Endocrinology 145(4), 1952-1960, doi: 10.1210/en.2003-0868 (2004).

https://doi.org/10.7554/eLife.58990.sa2

Article and author information

Author details

  1. Jiajia Xu

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Conceptualization, Data curation, Validation, Investigation, Methodology, Writing - original draft, Writing - review and editing
    Contributed equally with
    Yiyun Wang
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6084-2029
  2. Yiyun Wang

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Validation, Investigation, Writing - original draft
    Contributed equally with
    Jiajia Xu
    Competing interests
    No competing interests declared
  3. Ching-Yun Hsu

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Software, Investigation, Methodology
    Competing interests
    No competing interests declared
  4. Stefano Negri

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Data curation, Validation, Investigation, Methodology
    Competing interests
    No competing interests declared
  5. Robert J Tower

    1. Departments of Pathology, Johns Hopkins University, Baltimore, United States
    2. Departments of Orthopaedics, Johns Hopkins University, Baltimore, United States
    Contribution
    Software, Validation, Investigation
    Competing interests
    No competing interests declared
  6. Yongxing Gao

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Validation, Investigation, Writing - original draft
    Competing interests
    No competing interests declared
  7. Ye Tian

    1. Departments of Pathology, Johns Hopkins University, Baltimore, United States
    2. Department of Oral and Maxillofacial Surgery, School of Stomatology, China Medical University, Shenyang, China
    Contribution
    Validation, Investigation, Methodology
    Competing interests
    No competing interests declared
  8. Takashi Sono

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Investigation, Methodology
    Competing interests
    No competing interests declared
  9. Carolyn A Meyers

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Data curation, Validation, Investigation
    Competing interests
    No competing interests declared
  10. Winters R Hardy

    1. Departments of Pathology, Johns Hopkins University, Baltimore, United States
    2. UCLA and Orthopaedic Hospital Department of Orthopaedic Surgery and the Orthopaedic Hospital Research Center, Los Angeles, United States
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
  11. Leslie Chang

    Departments of Pathology, Johns Hopkins University, Baltimore, United States
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
  12. Shuaishuai Hu

    UCLA and Orthopaedic Hospital Department of Orthopaedic Surgery and the Orthopaedic Hospital Research Center, Los Angeles, United States
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
  13. Nusrat Kahn

    UCLA and Orthopaedic Hospital Department of Orthopaedic Surgery and the Orthopaedic Hospital Research Center, Los Angeles, United States
    Contribution
    Validation, Investigation
    Competing interests
    No competing interests declared
  14. Kristen Broderick

    Departments of Plastic Surgery, Johns Hopkins University, Baltimore, United States
    Contribution
    Resources
    Competing interests
    No competing interests declared
  15. Bruno Péault

    1. UCLA and Orthopaedic Hospital Department of Orthopaedic Surgery and the Orthopaedic Hospital Research Center, Los Angeles, United States
    2. Center For Cardiovascular Science and Center for Regenerative Medicine, University of Edinburgh, Edinburgh, United Kingdom
    Contribution
    Conceptualization, Supervision, Investigation, Visualization, Methodology, Project administration, Writing - review and editing
    For correspondence
    BPeault@mednet.ucla.edu
    Competing interests
    BP is the inventor of perivascular stem cell related patents held by the UC Regents (Patent number: 9730965) and is on the editorial board of Stem Cells.
  16. Aaron W James

    1. Departments of Pathology, Johns Hopkins University, Baltimore, United States
    2. UCLA and Orthopaedic Hospital Department of Orthopaedic Surgery and the Orthopaedic Hospital Research Center, Los Angeles, United States
    Contribution
    Conceptualization, Supervision, Funding acquisition, Validation, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    awjames@jhmi.edu
    Competing interests
    AWJ is a paid consultant for Novadip. This arrangement has been reviewed and approved by the JHU in accordance with its conflict of interest polices. AWJ receives funding for unrelated research from MTF Biologics and Novadip, and is on the editorial board of American Journal of Pathology and Bone Research. AWJ is the inventor of methods to purify CD107a progenitor cells held by the Johns Hopkins University.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2002-622X

Funding

National Institute of Arthritis and Musculoskeletal and Skin Diseases (R01 AR070773)

  • Aaron W James

National Institute of Arthritis and Musculoskeletal and Skin Diseases (K08 AR068316)

  • Aaron W James

National Institute of Dental and Craniofacial Research (R21 DE027922)

  • Aaron W James

Department of Defense (W81XWH-18-1-0121)

  • Aaron W James

American Cancer Society (Research Scholar Grant)

  • Aaron W James

American Cancer Society (RSG-18-027-01-CSM)

  • Aaron W James

Maryland Stem Cell Research Fund

  • Aaron W James

MTF Biologics

  • Aaron W James

British Israel Research Academic Exchange Partnership

  • Bruno Péault

British Heart Foundation

  • Bruno Péault

Department of Defense (W81XWH-18-1-0121)

  • Aaron W James

Department of Defense (W81XWH-18-1-0336)

  • Aaron W James

Department of Defense (W81XWH-18-10613)

  • Aaron W James

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank the JHU microscopy core facility, JHMI deep sequencing and microarray core facility, and Hao Zhang within the JHU Bloomberg Flow Cytometry and Immunology Core. The slides of human adipose tissue depots used in Figure 1—figure supplement 2 were donated by Drs. Trivia Frazier and Jeffrey Gimble at Obatala Sciences, Inc (New Orleans LA). AWJ was supported by the NIH/NIAMS (R01 AR070773, K08 AR068316), NIH/NIDCR (R21 DE027922), Department of Defense (W81XWH-18-1-0121, W81XWH-18-1-0336, W81XWH-18–10613), American Cancer Society (Research Scholar Grant, RSG-18-027-01-CSM), the Maryland Stem Cell Research Foundation, and MTF Biologics. In addition, MTF Biologics donated reagents for the study. This work was also supported by grants from BIRAX (the Britain Israel Research and Academic Exchange Partnership) and the British Heart Foundation (BHF) to BP. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Health, Department of Defense, or U.S. Army.

Ethics

Human subjects: All human tissues were obtained under Johns Hopkins University institutional IRB approval with a waiver of informed consent (Approval No. IRB00119905 and IRB00137530).

Animal experimentation: All animal experiments were performed according to the approved protocol of the Animal Care and Use Committee (ACUC) at Johns Hopkins University (Approval No. MO19M266 and RA19M268).

Senior Editor

  1. Clifford J Rosen, Maine Medical Center Research Institute, United States

Reviewing Editor

  1. Subburaman Mohan, Loma Linda University, United States

Reviewer

  1. Yuji Mishina, University of Michigan School of Dentistry, United States

Version history

  1. Received: May 15, 2020
  2. Accepted: September 25, 2020
  3. Version of Record published: October 12, 2020 (version 1)

Copyright

© 2020, Xu et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,622
    Page views
  • 163
    Downloads
  • 10
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Jiajia Xu
  2. Yiyun Wang
  3. Ching-Yun Hsu
  4. Stefano Negri
  5. Robert J Tower
  6. Yongxing Gao
  7. Ye Tian
  8. Takashi Sono
  9. Carolyn A Meyers
  10. Winters R Hardy
  11. Leslie Chang
  12. Shuaishuai Hu
  13. Nusrat Kahn
  14. Kristen Broderick
  15. Bruno Péault
  16. Aaron W James
(2020)
Lysosomal protein surface expression discriminates fat- from bone-forming human mesenchymal precursor cells
eLife 9:e58990.
https://doi.org/10.7554/eLife.58990

Further reading

    1. Medicine
    Luyang Cao, Lixiang Ma ... Jingsong Xu
    Research Article

    Billions of apoptotic cells are removed daily in a human adult by professional phagocytes (e.g. macrophages) and neighboring nonprofessional phagocytes (e.g. stromal cells). Despite being a type of professional phagocyte, neutrophils are thought to be excluded from apoptotic sites to avoid tissue inflammation. Here, we report a fundamental and unexpected role of neutrophils as the predominant phagocyte responsible for the clearance of apoptotic hepatic cells in the steady state. In contrast to the engulfment of dead cells by macrophages, neutrophils burrowed directly into apoptotic hepatocytes, a process we term perforocytosis, and ingested the effete cells from the inside. The depletion of neutrophils caused defective removal of apoptotic bodies, induced tissue injury in the mouse liver, and led to the generation of autoantibodies. Human autoimmune liver disease showed similar defects in the neutrophil-mediated clearance of apoptotic hepatic cells. Hence, neutrophils possess a specialized immunologically silent mechanism for the clearance of apoptotic hepatocytes through perforocytosis, and defects in this key housekeeping function of neutrophils contribute to the genesis of autoimmune liver disease.

    1. Medicine
    Hong Zheng, Qianjin Li ... Cheng-Kui Qu
    Research Article

    While mitochondria in different tissues have distinct preferences for energy sources, they are flexible in utilizing competing substrates for metabolism according to physiological and nutritional circumstances. However, the regulatory mechanisms and significance of metabolic flexibility are not completely understood. Here, we report that the deletion of Ptpmt1, a mitochondria-based phosphatase, critically alters mitochondrial fuel selection – the utilization of pyruvate, a key mitochondrial substrate derived from glucose (the major simple carbohydrate), is inhibited, whereas the fatty acid utilization is enhanced. Ptpmt1 knockout does not impact the development of the skeletal muscle or heart. However, the metabolic inflexibility ultimately leads to muscular atrophy, heart failure, and sudden death. Mechanistic analyses reveal that the prolonged substrate shift from carbohydrates to lipids causes oxidative stress and mitochondrial destruction, which in turn results in marked accumulation of lipids and profound damage in the knockout muscle cells and cardiomyocytes. Interestingly, Ptpmt1 deletion from the liver or adipose tissue does not generate any local or systemic defects. These findings suggest that Ptpmt1 plays an important role in maintaining mitochondrial flexibility and that their balanced utilization of carbohydrates and lipids is essential for both the skeletal muscle and the heart despite the two tissues having different preferred energy sources.