Comprehensive transcriptomic and proteomic
characterization of human mesenchymal stem cells
reveals source specific cellular markers

Anja M Billing, Hisham Ben Hamidane, Shaima S Dib, Richard J Cotton, Aditya M Bhagwat, Pankaj Kumar, Shahina Hayat, Noha A Yousri, Neha Goswami, Karsten Suhre, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Anja M Billing, Hisham Ben Hamidane, Shaima S Dib, Richard J Cotton, Aditya M Bhagwat, et al. Comprehensive transcriptomic and proteomic characterization of human mesenchymal stem cells reveals source specific cellular markers. Scientific Reports, 2016, 6, pp.21507. 10.1038/srep21507. hal-03125470

HAL Id: hal-03125470

https://hal.umontpellier.fr/hal-03125470

Submitted on 29 Jan 2021

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Comprehensive transcriptomic and proteomic characterization of human mesenchymal stem cells reveals source specific cellular markers

Received: 10 December 2015

Accepted: 26 January 2016

Published: 09 February 2016

Anja M. Billing, Hisham Ben Hamidane, Shaima S. Dib, Richard J. Cotton, Aditya M. Bhagwat, Pankaj Kumar, Shahina Hayat, Noha A. Yousri, Neha Goswami, Karsten Suhre, Arash Rafii & Johannes Graumann

Mesenchymal stem cells (MSC) are multipotent cells with great potential in therapy, reflected by more than 500 MSC-based clinical trials registered with the NIH. MSC are derived from multiple tissues but require invasive harvesting and imply donor-to-donor variability. Embryonic stem cell-derived MSC (ESC-MSC) may provide an alternative, but how similar they are to *ex vivo* MSC is unknown. Here we performed an in depth characterization of human ESC-MSC, comparing them to human bone marrow-derived MSC (BM-MSC) as well as human embryonic stem cells (hESC) by transcriptomics (RNA-seq) and quantitative proteomics (nanoLC-MS/MS using SILAC). Data integration highlighted and validated a central role of vesicle-mediated transport and exosomes in MSC biology and also demonstrated, through enrichment analysis, their versatility and broad application potential. Particular emphasis was placed on comparing profiles between ESC-MSC and BM-MSC and assessing their equivalency. Data presented here shows that differences between ESC-MSC and BM-MSC are similar in magnitude to those reported for MSC of different origin and the former may thus represent an alternative source for therapeutic applications. Finally, we report an unprecedented coverage of MSC CD markers, as well as membrane associated proteins which may benefit immunofluorescence-based applications and contribute to a refined molecular description of MSC.

Mesenchymal stem cells (MSC) are self-renewing multipotent cells which hold great potential in regenerative medicine and tissue engineering. Their versatility, homing potential to injured tissue, and the limited risk of tumorigenesis¹ render them an interesting tool in cell-based personalized therapy. Leveraging differentiation ability, strong immunomodulatory capacity and paracrine support, they have been used in the experimental treatment of a plethora of diseases, evident in the more than 500 active MSC-based clinical trials worldwide currently registered with the NIH. Studies have shown that treatment with MSC is beneficial to patients with Graft versus Host Disease²,³, autoimmune diseases⁴⁻⁶, cardiovascular⁷ and liver disease⁸. In MSC-based therapies, the actual engraftment of MSC is considered to be very low and not persistent⁹,¹⁰. Tissue regeneration is initiated and supported by MSC at the site of tissue damage mainly by immunomodulation, through both cell-free (secretion) or cell-based mechanisms. Importantly, allogenic MSC do not induce a strong host response due to their low expression of MHC class II molecules and lack of co-stimulatory receptors, rendering them transplantable with comparative ease¹¹.

MSC were first isolated from bone-marrow¹² but actually reside perivascular in a wide range of organs and tissues such as adipose tissue, skin, muscle, tendon, lungs, dental pulp, umbilical cord, placenta, with each population displaying individual differentiation potential and phenotype¹³⁻¹⁵. In general, MSC can be differentiated *in vitro* into cells of the mesodermal lineage such as osteocytes, adipocytes and chondrocytes¹⁶,¹⁷. A smaller number of studies have also shown differentiation into myocytes, cardiomyocytes and even transdifferentiation into

Research Division, Weill Cornell Medicine-Qatar, Doha, State of Qatar. Correspondence and requests for materials should be addressed to J.G. (email: jog2030@qatar-med.cornell.edu)
cells of the endoderm or ectoderm<sup>17,18</sup>. Adult tissue derived MSCs can be expanded *in vitro* but lose proliferation, differentiation and immunomodulation potential gradually over time<sup>18,19</sup>. Recently described embryonic stem cell-derived MSC (ESC-MSC), however, are highly proliferative and an interesting source for clinical use. ESC-MSC carry typical MSC markers and have been shown to have similar *in vitro* and *in vivo* capabilities as BM-MSC, including homing, immunomodulation, regeneration and differentiation<sup>20-26</sup>. One *in vivo* study even showed the superiority of ESC-MSC in comparison to BM-MSC in the treatment of an experimental autoimmune encephalomyelitis (EAE) model of Multiple Sclerosis<sup>27</sup>.

To date, no surface markers have been exclusively associated with MSC. Based on a proposal by the International Society for Cell Therapy, cells may be classified as MSC if they adhere to plastic, carry a minimal subset of characteristic surface markers (CD73, CD90, CD105), are negative for hematopoietic lineage markers (CD14, CD19, CD45) and present the potential to differentiate into bone, fat and cartilage<sup>28</sup>.

For the development of targeted MSC-based therapy, an unambiguous definition of MSC is highly desirable and comprehensive molecular characterizations, as delivered by “Omics” methodologies, may be one approach on the path towards that goal<sup>29</sup>.

Given their donor-derived nature and limited expansion capacity in cell culture, the production of adult MSC at the scale required for efficient therapeutic use is challenging. MSC derived from embryonic stem cells (ESC-derived MSC; ESC-MSC), in comparison, are an alternative source with additional potential benefits with respect to biosafety<sup>30</sup>. A future goal may be the use of induced pluripotent stem cells to produce personalized MSC at large scale prior to treatment.

Establishing molecular equivalency (or the lack thereof) between differently sourced MSC is a central question for future clinical use of ESC-MSC. We present a comprehensive transcriptomic and proteomic characterization of ESC-MSC and their comparison to bone marrow-derived MSC as well as to embryonic stem cells (ESC). The ESC-MSC used were derived through a protocol potentially adaptable to clinical use, which enriches a uniform MSC population from spontaneously differentiated ESC<sup>31</sup>. The rationale for selecting this method was based on its simplicity, robustness and the adaptability to clinical applications. To our knowledge, this is the first large-scale proteomic characterization of ESC-MSC. Two global profiling techniques were combined for a comprehensive analysis: RNA deep sequencing as well as quantitative high resolution nano-LC-MS/MS based on stable isotope labeling with amino acid in cell culture (SILAC)<sup>32</sup>. This extensive comparison delivers unprecedented coverage of cluster of differentiation (CD) markers, which may contribute to a more stringent definition of MSC as a cell type.

## Methods

### Cell culture and SILAC labeling of embryonic stem cells (ESC) and mesenchymal stem cells (ESC-MSC, BM-MSC).
ESC-derived MSC were generated according to a published protocol<sup>31</sup>. Briefly, ESC seeded in 6-well plates were cultured under feeder-free condition on growth factor reduced Matrigel (BD Invitrogen) with daily changes of mTeSR1 medium (Stem Cell Technologies) supplemented with 1% penicillin/streptomycin. Three days after passage with accutase, medium was changed to DMEM/low glucose supplemented with 10% FBS, GlutaMax and 1% penicillin/streptomycin. Cells were left undisturbed for 6 days with one medium change after 3 days. Differentiating ESC were passaged 1:3 using trypsin and seeded on plastic. After another 6 days of culture with one medium change after three days, cells were passaged again with trypsin and seeded on plastic. From this passage onwards, cells are fully differentiated into ESC-MSC. ESC-derived MSC have been phenotyped and show classic MSC behavior. They are plastic-adherent, show trilineage differentiation and express MSC-typical surface markers (CD73, CD90, CD105)<sup>31</sup>. ESC-MSC used for transcriptomic and proteomic profiling were generated in three independent experiments. Bone marrow-derived mesenchymal stem cells (Lonza, StemCell Technologies) were purchased with their mesenchymal characteristics verified. Details of the BM-MSC donors: 40y/m (StemCell, MSC-001F, lot#BM2893), 39y/m (Lonza, PT2505, lot#1F3422), 27y/m (Lonza, PT2505, lot#318006), 20y/m (Lonza, PT2505, lot#8F3520). Since the origin ESC line is male, only male BM-MSC donors were selected. For proteomic and transcriptomic profiling, BM-MSC of three (40y, 39y, 27y) and four different donors (40y, 39y, 27y, 20y) were included, respectively. Representative pictures for ESC, ESC-MSC and BM-MSC are shown (Figure S1A). All cells were cultured at 37°C in an atmosphere supplemented to 5% CO₂.

Permission to use the human embryonic stem cell line (hESC) ES04 (WiCell institute) was obtained from the Cornell/Rockefeller/Sloan Kettering tri-institutional ESC research oversight committee. Funding was secured from nonfederal, US-external funding sources. ESCs were maintained under feeder-free condition on growth factor reduced Matrigel (BD Invitrogen) in mTeSR1 medium (Stem cell technologies) supplemented with 1% (v/v) penicillin/streptomycin. Medium was changed daily. ESC-MSC and BM-MSC were maintained in DMEM low glucose supplemented with 10% FBS, GlutaMax and 1% penicillin/streptomycin. ESC-MSC were sub-cultured every 3-4 days, BM-MSC every 7 days. Cells were harvested at 80% confluence at the following passages: ESC at passages 89-91, BM-MSC at passages 5-6 and ESC-MSC at passages 7-9.

SILAC labeling was performed as described<sup>33</sup>. Arginine/lysine-free medium (custom made-mTeSR1 for ESC, DMEM low glucose with dialyzed FBS for MSC) was supplemented either with L-lysine – U-¹³C₆¹⁵N₂ (Lys8) and L-arginine – U-¹³C₆¹⁵N₄ (Arg10) for “heavy” labeling, or with L-lysine – U-¹³C₄¹⁵N₀ (Lys4) and L-arginine – U-¹³C₆¹⁵N₀ (Arg6) for “medium” labeling, or with standard amino acids for “light” labeling at a concentration matching the original medium. Labeled amino acids were purchased from Cambridge Isotope Laboratories. Cells were labeled for >5 doublings and incorporation was verified prior to sample mixing. Labeled protein samples were combined as shown in Figure S4A.

### Immunoblotting.
Proteins (5 µg) were separated by SDS-PAGE (Bis-Tris NuPAGE gradient gels 4-12%, Lifescience technologies) and transferred onto PVDF membranes (HybondLF, GE Healthcare). The membranes were blocked in TBS/0.05% Tween20 containing 3% BSA and were incubated overnight with primary antibodies (Cell Signaling: αLSD1 #2139, αMEK1/2 #4694, αhistone H2A.Z #2718) followed by horse radish peroxidase
conjugated secondary antibodies (SantaCruz Technologies: α-rabbit #sc-2317, α-mouse #sc-2314, Promega: α-mouse #W4028). Membranes were incubated with the ECL™ Prime Western Blotting Detection Reagent (GE Healthcare) and chemiluminescence was detected with a CCD camera (SnapGene v.7.04).

**Next generation RNA sequencing.** Total RNA was extracted with TRIZOL and further cleaned on RNAeasy columns (Qiagen). RNA quality was evaluated by the Agilent Bioanalyzer 2100 system. Only preparations with RIN (RNA integrity number) values above 8 were considered.

Following RNA isolation, 100 ng of total RNA was converted to cDNA using the Ovation RNA-seq System V2 (Nugen Technologies, San Carlos, CA). 2µg of the amplified cDNA was sheared to 150–200 bp size distribution by Adaptive Focused Acoustics using a Covaris E220 instrument (Covaris, Woburn, MA) under the following conditions: 50µl total volume, 20% duty cycle, 175 intensity, 200 cycles per burst, 165 sec in frequency sweeping mode. The remainder of the library preparation followed manufacturer's protocol as described in Encore Rapid IL Multiplex System (Nugen Technologies). Briefly, the sheared cDNA was end-repaired to generate blunt ends then ligated to Illumina compatible adaptors with indexing tags, followed by 1× AMPure XP beads purification. The final NGS libraries were quantified using Agilent Bioanalyzer DNA Chip 1000 then pooled; 11 libraries per pool. Paired-end 100 bp deep sequencing was carried out on HiSeq 2500 (Illumina).

Using a customized analysis pipeline (Figure S1B) all samples were run together on a virtual machine. The analysis pipeline implemented as a bash shell script is composed of TopHat³⁴, Picard, (http://picard.sourceforge.net/index.shtml), Samtools³⁵ and Cuffdiff³⁴. TopHat aligns RNA-seq reads to the reference genome (GRCh37). Picard removes PCR duplicates and Samtools creates indices for bam files. We ran the Cuffdiff function from the Cufflinks package v2.1.1.³⁴ to find differentially expressed genes and transcripts in the samples using the default settings. Results from Cuffdiff were plotted in R using the Bioconductor CummeRbund package³⁴.

Information for mapped, unmapped and duplicate reads was extracted using Samtools. The length of the transcript in base pairs was calculated by summing up depth per genome position. The covered fraction of transcript was computed by dividing the total number of bases in mapped reads with the length of the transcript. Depth per genome position was calculated using GenomeCov function from BedTools³⁶. RNA quality control was performed with RNA-SeQC³⁷. Read metrics are summarized in Table S1B.

**Sample preparation for mass spectrometry.** Subcellular fractionation was performed using the Nuclear Extract Kit (ActiveMotif) resulting in cytosolic (CYT), nuclear (NUC) and chromatin-bound (CH) proteins. All extraction buffers were supplemented with protease inhibitors (Complete EDTA-free, Roche) and phosphatase inhibitors (PhosStop, Roche). Protein was methanol/chloroform precipitated³⁸ and resuspended in urea/thiourea buffer (6 M/2 M, 30 mM HEPES, pH 8). Protein concentration was determined by Bradford. Unlabeled and SILAC (M, H) samples were mixed at a ratio 1:1:1 using 60 µg of protein per sample. After reduction (5 mM DTT, 30 min) and alkylation (10 mM iodoacetamide, 20 min, in the dark), proteins were digested for 3 h by Lys-C (Wako) at a protein:enzyme ratio of 100:1. Samples were diluted with 10 mM ammonium bicarbonate to 2 M urea prior to overnight trypsin digestion (Sequencing grade, Promega) at a protein:enzyme ratio of 100:1.

**Isoelectric focusing.** Peptides (180 µg) were separated by in-solution isoelectric focusing (OFFGEL fractionator, Agilent) into 12 fractions over a pH range of 3–10. Fractionation was performed according to manufacturer's protocol with adaptation. Glycerol in the running buffer was reduced to 0.3% (original: 6%) and the ampholytes to 0.1% (original: 1%). Peptides were focused for 20 kVh. Peptides were harvested and wells incubated with 50 µl 50:49:1 methanol:H₂O:TFA for 15 min, pooling the wash with the corresponding fractions. Fractionated peptides were dried down with a speedvac concentrator and kept at –80 °C until further use. Peptides were solubilized by 1% acetonitrile/0.05% TFA and desalted on C₁₈ STAGETips³⁹. Briefly, STAGETips were packed with 3 C₁₈ disks and activated with methanol. After 2 washes of 2% acetonitrile/0.1% TFA, peptides were loaded. STAGETips were washed with 0.1 % acetic acid and 2% acetonitrile/0.1% TFA before elution with 60% acetonitrile/0.1% TFA. Eluates were dried down completely in a speedvac concentrator and reconstituted in 10 µl of mobile phase A (0.5% acetic acid) prior LC-MS/MS analysis.

**Mass spectrometry.** *Nano LC-MS/MS analysis.* Of the resolubilized peptide solution, 6–10 µl were subsequently analyzed by liquid chromatography (LC) using an EASY nLC-II coupled to a Q Exactive mass spectrometer (Thermo Scientific, Bremen, Germany)³³. Emitter columns were packed in-house with ReproSil-Pur 120 C18-AQ 3 µm diameter beads (Dr. Maisch GmbH, Germany). Liquid chromatography analysis was performed with a segmented gradient from 10 to 60% mobile phase B over 100 min at a flow rate of 250 nl/min. Mobile phase composition was as follows: mobile phase A, H₂O with 0.5% acetic acid; mobile phase B, H₂O: acetonitrile (20:80 volume ratio with 0.5% acetic acid). A data dependent mass spectrometric method was employed, where the 10 most intense ions, excluding unassigned charge states and singly charged ions, detected in the preceding full scan are isolated (3% isolation width) and fragmented using higher energy collisional dissociation (HCD) (normalized collision energy 25). Precursor scans (MS¹ level) were acquired at a resolution of 70,000 (m/z 300) and an AGC target value of 3,000,000 charges (maximum ion injection time 20 ms). Fragmentation spectra were acquired at a resolution of 17,500 (m/z 300) and an AGC target value of 100,000 charges (maximum ion injection time 120 ms). All scan events were recorded in profile mode. A dynamic exclusion list of 25 s was employed and the exclude isotopes functionality was activated.

*Data analysis.* MS data was analyzed by MaxQuant suite of algorithms version 1.4.1.2⁴⁰ using a Homo Sapiens database downloaded from UniprotKB on the 27th of November 2013 and comprising 88,473 protein isoforms entries. The following default search parameters were employed: first search mass accuracy tolerance 20 ppm, main search mass accuracy tolerance 4.5 ppm, FTMS MS/MS tolerance 20 ppm, minimum peptide length of 7
amino acids, peptide spectrum match FDR and protein FDR both set to 0.01 as calculated by the revert database approach⁴¹. Protein quantification was based on a minimum of two ratio counts, originating from unique or razor peptides only. Additionally, unless explicitly stated otherwise, other parameters of the data analysis were not changed from their MaxQuant 1.4.1.2 default value.

Multiplicity was set to 3 with Arg6/Lys4 and Arg10/Lys8 for medium and heavy labels respectively as required by the triple SILAC approach and the search was conducted using Trypsin/P enzyme specificity, allowing for maximum two missed cleavages, N-terminal acetylation and methionine oxidation as variable modifications and cysteine carbamidomethylation as fixed modification. Both re-quantify and match between runs (match time: 1 min, alignment time window: 20 min) functionalities were enabled. MaxQuant analysis was performed both on the combined (36 MS runs per sample) and separate organelle fractions (12 MS runs per sample). Despite reducing the ampholyte concentration during in-solution isoelectric focusing (0.1% instead of the recommended 1%) and additional washing on C₁₈ STAGETips, ampholytes eluted in the beginning of the chromatogram, interfering with timely MaxQuant processing. Raw files were thus preprocessed with the RecalOffline application (ThermoFisher Scientific) to exclude the first 35 min of acquisition.

**Bioinformatics and statistical analysis.** Proteomics data sets were analyzed with the empirical Bayes moderated T test implemented by the limma bioconductor package in the R environment. P values were corrected for multiple hypothesis testing using the Benjamini-Hochberg method (FDR ≤ 0.05). Differential expression was calculated on normalized log₂ ratios (nano LC-MS/MS). Differential expression for RNA-sequencing data was performed with the Cufflinks package³⁴ using the Cuffdiff function. PCA analysis was performed with the FactoMineR package v.1.29 within the R environment on FPKM values for RNA and on non-normalized intensities derived from mass spectrometry for PROT. Differentially expressed genes and proteins were processed using the in-house developed Comics R package for enrichment analysis.

**Data availability.** The mass spectrometry proteomics data, including all raw files and data analysis tables produced by MaxQuant analysis, have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository⁴² with the dataset identifier PXD001856. RNA sequencing data can be accessed via the Sequence Reads Archive by NCBI (BioProject ID: PRJNA277616).

# Results and Discussion

**Deep transcriptome analysis by RNA sequencing.** As a first step in the comprehensive characterization of mesenchymal stem cells (ESC and bone marrow derived) we comparatively analyzed the transcriptome by deep sequencing. Currently, a single study is published characterizing ESC-MSC with deep RNA sequencing. The cells used, however, were derived by co-culture with OP9 cells⁴³ and are thus not compatible with clinical usage.

The RNA-seq pipeline is depicted in Figure S1B. Differential expression analysis by CuffDiff³⁴ led to 7262, 7073 and 2521 modulated genes (FDR ≤ 0.05) for BM-MSC vs ESC (BM-E), ESC-MSC vs ESC (EM-E), and BM-MSC vs ESC-MSC (BM-EM), respectively. Among those are 3578, 3394 and 1001 modulated isoforms (FDR ≤ 0.05) corresponding to 2452, 2285 and 868 genes for BM-E, EM-E and BM-EM (Fig. 1; Figure S2A,B; Table S1). As a control, a set of pluripotency and MSC markers were analyzed and as expected, separation of the cell types (ESC and MSC) was observed (Figure S2C). A Pearson correlation matrix clearly demonstrates the differences between ESC and MSC, but also shows that the two MSC subpopulations (ESC-MSC and BM-MSC) segregate. ESC-MSC are closer to ESC (0.53–0.56) than BM-MSC to ESC (0.44–0.54). The two younger BM-MSC donors (27y, 20y–BM-MSC 3 and 4 respectively) are closer to ESC (0.5–0.54) than the other two (40y, 39y–BM-MSC 1 and 2 respectively) (0.44–0.49). This may be interpreted as ESC-MSC representing a particularly young subclass of MSC. Correlations between both MSC populations range from 0.86–0.89, whereas intra-group correlations range from 0.95–0.96 (EM) and 0.93–0.95 (BM) (Fig. 1c). ESC to MSC differences may be explained in PCA analysis by the first two components covering 94% of variance. Component 3 separates the two MSC populations, explaining less than 5% of variance. As in the correlation matrix, BM-MSC donors show a larger intra-group variation than their ESC-MSC counterparts (Fig. 1d).

The observed differences between the two MSC populations are similar in scope to the variability recently observed between differentially sourced MSC⁴⁴,⁴⁵. They may thus be reflective of ESC-MSC representing yet another sourcing-equivalent MSC subtype or indicate that ESC-MSC are closer to adult MSC from a different source, sharing expression differences and differentiation potential. This hypothesis is further strengthened in a direct side-by-side comparison with a published data set covering fetal (placenta-derived) and adult (bone marrow-derived) MSC⁴⁵ using our bioinformatics pipeline. A similar distance between the two MSC subtypes was observed (Figure S3).

When comparing up-regulated genes in MSC versus ESC (BM > E, EM > E), only minor changes can be observed between the two MSC subpopulations. Top enriched biological processes (GOBP) in MSC are vesicle-mediated transport, regulation of signal transduction and ECM organization. Enriched terms of the cellular compartment category (GOCC) include cytoplasm, vesicle, exosomes, organelle membrane, focal adhesion and Golgi apparatus. Tissue regeneration initiated by MSC is known to happen in a paracrine manner and the observed strong enrichment of exosomes as well as vesicle-mediated transport further confirms the central role of this process in regeneration support. Notably, developmental processes with a wide range of different tissues are enriched in MSC (Figure S6A, left, Table S3). Those include blood vessels, neurons, skeleton, muscle, heart, connective tissue, cartilage, kidney, gland, epithelium, urogenital tract, eye and others, reflecting the wide range of MSC application in clinical and laboratory research.

MSC from different sources, although expressing the classical surface markers (CD73, CD90, CD105) and exhibiting similar trilineage differentiation potential, show distinct molecular signatures, as previously reported⁴⁴⁻⁴⁶. Most significant changes between both MSC populations were observed for membrane-bound
<figure><img src="image_4.png" /><figcaption>Figure 1. Transcriptome analysis (RNA) of ESC-MSC (EM), BM-MSC (BM) and ESC (E) by RNA-seq.</figcaption></figure>

(a) Volcano plots for BM-MSC vs ESC (left panel), ESC-MSC vs ESC (middle panel) and BM-MSC vs ESC-MSC (right panel). Numbers are given for the total of differentially expressed genes as well as for the up-regulated per cell type. (b) Venn diagram of SDEs per comparisons (BM-MSC vs ESC, ESC-MSC vs ESC and BM-MSC vs ESC-MSC). (c) Pearson correlation matrix for all replicates. (d) PCA analysis based on FPKM values; principal components 1 to 3 were plotted against each other. (e) GO-term enrichment analysis for biological process (GOBP) and cellular component (GOCC). Bar charts represent the most significant top 20 terms of each category for each cell type sorted by mean −log₁₀ p values.

glycosylated or secreted proteins, indicating differential signaling or homing properties. Functionally, those proteins from both MSC populations are enriched for ECM organization, cell differentiation, and developmental
processes (Fig. 1e, Figure S6A). Uniquely enriched among up-regulated genes of BM-MSC are genes associated with mitotic cell cycle and regulation of proliferation, which is counterintuitive since BM-MSC are slower in their proliferation rate than ESC-MSC. However, upon applying a fold change cutoff of 2 log₂, all cell cycle related functions disappear, indicating the relative low amplitude of the change, possibly originating from overall low signal to noise of the associated probes. Transcripts present exclusively in BM-MSC (n = 71) are enriched in extracellular vesicle proteins. Among the RNA transcripts present only in ESC-MSC (n = 19) are four transcription factors or regulators involved in developmental processes (HOXD1, NKX2-5, LHX2, FGF12), concordant with our previous study where we described the increased differentiation ability of these MSC⁴⁷. The majority of ESC-MSC specific transcripts (n = 12) is unknown, assigned to neither genes nor noncoding RNAs.

**LC-MS/MS-based proteomics combined with SILAC.** For the quantitative proteomic analysis, experiments were designed to use triple SILAC labeling with ESC, ESC-MSC and BM-MSC labeled as “light”, “medium” and “heavy”, respectively. Three biological replicates were generated including label swapping. Biological triplicates for BM-MSC originated from three donors, whereas ESC-MSC triplicates were derived from independent differentiation experiments. For increased proteome coverage, subcellular fractionation into cytosol (CYT), nucleus (NUC), and chromatin (CH) was performed⁴⁸ (Figure S4A,B,C). Three markers have been selected to verify the enrichment of cytosolic, nuclear and chromatin-bound proteins by western blot (Figure S4D). To validate the fractionation, the dual specificity mitogen-activated protein kinase kinase 2 (gene: MAP2K2) was selected as a representative cytosolic protein, histone H2AZ (gene: H2AFV) as exclusively chromatin-bound, and lysine-specific histone demethylase 1 (gene: KDM1A) as present in both nuclear and chromatin-bound fractions. All three proteins were measured by mass spectrometry and their intensity profiles (Figure S4E) validated the western blot data and confirmed the specificity of the organelle fractionation technique employed. In the combined analysis (CYT, NUC and CH together), 8408 proteins were identified with 7129, 7023 and 7136 quantified proteins per replicate (FDR ≤ 1% on both peptide and protein level) (Figure S4F, Table S2). As an internal control, the data set was screened for proteins expected in ESC and MSC, respectively. Pluripotency-associated proteins identified include SOX2, POU5F1, SALL4, UTF1, DNMT3B, DPPA4, LIN28A, and GRB7. As classical MSC phenotyping surface markers⁴⁹ we identified ALCAM (CD116), ITGB1 (CD29), CD44, ENG (CD105), NT5E (CD73), and THY1 (CD90) (Figure S4G). As expected, expression patterns observed followed the cell type origin. For instance, bone marrow stromal cell antigen 1 (BST1, CD157), found to be only enhanced in BM-MSC, is a MSC marker involved in self-renewal, migration and osteogenic differentiation⁵⁰.

Differential expression of the combined analysis was assessed by the empirical Bayes moderate T test (limma, FDR ≤ 0.05) with 3733, 3401 and 34 proteins differentially expressed between BM-E, EM-E and BM-EM, respectively (Fig. 2a,b). Unlike RNA sequencing, only a low number of proteins are differentially expressed between BM-MSC and ESC-MSC. When reporting correlation coefficients for all samples based on their non-normalized log₁₀ raw intensities across all experiments (Fig. 2c) ESC and MSC clearly separate. Correlation coefficients between ESC and MSC range from 0.86–0.89 for ESC-MSC and 0.83–0.87 for BM-MSC, while the two MSC subtypes correlate in the range 0.94–0.96. The close relationship between ESC-MSC and BM-MSC is also validated by PCA analysis, where approximately 95% of the variance is explained by the first two components separating ESC and MSC (Fig. 2d).

Comparing our proteomic dataset on ESC and MSC to the two most comprehensive human proteome maps published to date, which contain together 19608 proteins and cover a wide range of tissues from adult and fetal origin⁵¹,⁵², we uniquely identified 54 proteins and 5409 peptides, which map to 3024 proteins (Figure S5A). Our study thus provides new additional data to the human proteome. Given their likely stem cell specificity and the origin of the data in the maps, the 54 proteins are unsurprisingly largely uncharacterized and thus lack functional data.

In comparison to the most comprehensive studies available on human ESC⁵³,⁵⁴ (Figure S5B), our data set contains 4500 proteins commonly expressed within the three studies which may be seen as an ESC core proteome. The data set adds 1661 unique identifications, expanding the ESC proteome and contributing to ESC biology. Finally, comparing our data set to the most comprehensive MSC proteomics study available to date⁵⁵, our relative contribution to the reported MSC proteome is even larger adding 3331 unique identifications (Figure S5C).

Among all proteins identified, 62 were classified as uncharacterized proteins including KIAAs and ORFs, from which 20 were significantly differentially expressed, hinting at a role in stem cell biology. A large proportion of transcription-related proteins in our dataset are zinc finger proteins (ZNF) (207) which is believed to be the largest family of transcription factors in mammals. Most of these ZNF remain uncharacterized⁵⁶, but some have been shown to be implicated in developmental processes⁵⁷,⁵⁸, with recent examples involving chondrogenesis⁵⁹ (ZNF449) and adipogenesis⁶⁰ (ZNF395). Therefore, the differentially expressed ZNF proteins observed here in ESC and MSC strongly suggests they also play an important role in stem cell biology. The highest up-regulated protein in ESC-MSC is ZNF322 (177 fold), associated with MAPK signaling and heart development⁶¹. ZFHX4, up-regulated in both ESC-MSC and BM-MSC is involved in neuronal differentiation⁶². FHL2, up-regulated in both MSC populations as well, controls osteogenic differentiation⁶³. Based on bioinformatics analysis we observed an enrichment of zinc fingers of the LIM-type in both MSC types. Functionally, these are associated with cytoskeleton organization and focal adhesion, but also tissue-specific differentiation⁶⁴.

The human proteome map published by the Pandey group⁵¹ includes tissues and cell types of fetal and adult origin. Among the 30 tissues analyzed, 2350 proteins were found to be expressed ubiquitously and can be considered house keeping proteins or core proteome. Filtering those out of protein lists should thus facilitate the selection of signature proteins that define cell type functionality. GO term enrichment analysis was thus performed on the 200 most intense proteins of ESC-MSC and BM-MSC after removing the core proteome (Fig. 2e). The top biological processes enriched in the 200 most intense MSC proteins include extracellular matrix organization,
<figure><img src="image_5.png" /><figcaption>**Figure 2. Proteome analysis (PROT) of ESC-MSC (EM), BM-MSC (BM) and ESC (E) by nano LC-MS/MS.**
(a) Volcano plots for BM-MSC vs ESC, ESC-MSC vs ESC and BM-MSC vs ESC-MSC. Numbers are given for the total of differentially expressed proteins as well as for the up-regulated per cell type. (b) Venn diagram of SDEs per comparisons (BM-MSC vs ESC, ESC-MSC vs ESC and BM-MSC vs ESC-MSC). (c) Pearson correlation matrix for all replicates. (d) PCA analysis based on non-normalized log₁₀ intensities was performed; principal components 1 to 3 were plotted against each other. (e) GO-term enrichment analysis for biological process (GOBP) and cellular component (GOCC). Bar charts represent the most significant top 20 terms of each category for each cell type sorted by mean −log₁₀ p values.</figcaption></figure>
adhesion and developmental processes (vasculature and endoderm). As seen with RNA-seq, vesicle-mediated transport, focal adhesion and vesicles/exosomes are among the most enriched terms in MSC when comparing to ESC (Fig. 2e).

**Data integration of the two large-scale techniques.** A large overlap of identifications can be obtained (6959) when comparing RNA-seq and LC-MS/MS data (Fig. 3a), extending to those measured as differentially expressed. 60% of proteomic (LC-MS/MS) changes were validated by RNA-seq (Fig. 3b). Both techniques give very similar results in terms of global trends of expression as emphasized by the Pearson correlation coefficient of 0.66 in agreement with the literature⁶⁵.

**Fundamental differences between ESC and MSC.** ESC and MSC differ extensively as demonstrated by PCA analysis for both techniques (Figs 1d, 2d). RNA-seq and proteomics validated genes/proteins differentially expressed between ESC and MSC (Fig. 3c) were subjected to enrichment analysis (Fig. 3d). Upon performing the analysis on cellular component associated terms (GOCC): cytoplasm, endomembrane system, vesicle, extracellular exosomes, focal adhesion, and Golgi apparatus were associated to proteins/genes up-regulation in MSC. Consistently, the biological processes (GOBP) found enriched in MSC by both techniques included vesicle-mediated transport, cellular component organization, extracellular matrix organization and intracellular transport. When refining the GOBP analysis by filtering for developmental terms, the top 20 enriched processes in MSC included: axon, neuron projection, neuron, vasculature, cardiovascular system, skeletal and muscle development. In total, 45 developmental terms (p val ≤ 0.05) were commonly enriched for both MSC subtypes (ESC-MSC and BM-MSC) by RNA and PROT, including for instance, and in addition to the above mentioned, heart, cartilage, mesenchyme, bone, and endoderm development. If considering developmental terms found exclusively by RNA-seq (p val ≤ 0.05), then kidney, gland, liver, eye, and urogenital system development may be added to the list.

Comparatively, upon conducting the analysis on cellular component associated terms (GOCC) for proteins/ genes up-regulated in ESC, nucleus, nuclear lumen, nucleoplasm, nucleolus, and chromosome showed significant enrichment. The corresponding top biological processes included chromosome organization, RNA processing, RNA splicing and cell cycle. These findings reflect the highly proliferative nature of ESC (cell division time of 15–16 h as compared to 2–3 d in MSC⁶⁶) and high relative contribution of the nucleus to total cell volume in this cell type.

**How do ESC-MSC compare to BM-MSC.** Central to answering the question on usability of ESC-MSC in clinical applications is their quantitative comparison to the currently employed standard: BM-MSC. RNA-sequencing reported a total of 2521 characteristic features differentially expressed between the two MSC populations, with 1501 and 1020 respectively up- and down-regulated when comparing BM-MSC to ESC-MSC (Fig. 1a). However, LC-MS/MS based analysis reported minimal differences with only 34 significantly differentially expressed proteins, 24 up- and 10 down-regulated, when comparing BM-MSC to ESC-MSC (Fig. 2a). This is presumably due to a lower analytical depth of PROT compared to RNA as well as an increased stringency of the MS-based statistical analysis, but also reflects changes at the RNA level that did not systematically translate to the protein level (Fig. 3b).

Applying significance B statistical analysis⁴⁰ instead of limma allowed for a more extensive comparison between PROT and RNA-seq results. A total of 1739 significantly differentially expressed proteins, 887 overexpressed in BM-MSC and 852 overexpressed in ESC-MSC, were reported using this approach (Fig. 4a). SDE proteins in at least two replicates (Fig. 4a intersections, Fig. 4b) and in any replicate (Fig. 4a union, Fig. 4c) were compared to RNA-seq data. Based on all intersecting significance B SDE proteins (240 for BM-MSC, 146 for ESC-MSC): 137 and 70 candidates (for BM-MSC and ESC-MSC respectively) were validated by RNA-seq. When considering the union of all significance B SDE proteins instead (887 for BM-MSC, 852 for ESC-MSC) those numbers increase to 321 and 187 validated proteins for BM-MSC and ESC-MSC, respectively.

GO term enrichment analyzes were performed for validated proteins (RNA/PROT overlap) of significance B intersection (Fig. 4d, left panels) and union (Fig. 4d, middle panels) as well as for all significance B SDE proteins measured by LC-MS/MS (Fig. 4d, right panels). Based on RNA validated intersecting significance B proteins, ECM organization and extracellular exosome pathways were found to be enriched in MSC. Upon decreasing stringency by including RNA validated proteins from the significance B union, similar functions were enriched at low p values, specifically developmental processes and exosomes. Furthermore, when considering all significance B SDE proteins identified in PROT only, a remarkably high enrichment of extracellular vesicle exosome (ESC-MSC: 212, p val 9.92E-16; BM-MSC: 307, p val 5.80E-52), ECM organization (ESC-MSC: 53, p val 1.13E-11; BM-MSC: 63, p val 1.71E-16) and developmental process (ESC-MSC: 320, p val 4.31E-06; BM-MSC: 348, p val 9.55E-09) was observed. These results consistently point towards exosome production playing a central role in both MSC subpopulations, and overall close resemblance of both MSC subtypes. A finer grained comparison between ESC-MSC and BM-MSC highlights biological processes expression levels differences relating to their source: as expected cell cycle and nuclear division related processes are enriched in ESC-MSC, consistent with their higher proliferation rate, whereas BM-MSC show enhanced extracellular organization related processes, consistent with previously reported findings.

Filtering the biological process enrichment for developmental terms, as done for the ESC to MSC comparison, hints at other specificities distinguishing both MSC subtypes (Fig. 4e). Interestingly, ESC-MSC show higher enrichment for neuron and axon development whereas BM-MSC are more enriched for vasculature development. Additionally and as expected, ESC-MSC are also exclusively enriched in stem cell and embryo development terms. Overall, developmental terms related functions appear very similar in both MSC populations.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Combined analysis (RNA, PROT) of MSC and ESC. (a) Comparison of identified features based on gene name. (b) Comparison between significant differentially expressed proteins (SDEs) (FDR ≤ 0.05) found by RNA and PROT for all six cases: EM > E, BM > E, BM > EM, EM < E, BM < E, BM < EM; comparisons are based on common genes/proteins identified by RNA and PROT. (c) Venn diagrams display SDEs of genes/proteins up-regulated in MSC (left) and ESC (right) and validated by both techniques (RNA and PROT). (d) GO-term enrichment analysis on RNA and PROT validated genes/proteins up-regulated in MSC (top) and ESC (bottom). Bar charts represent the most significant top 20 terms of each category per cell type and sorted by mean −log₁₀ p values. In addition to GOCC (left panels) and GOBP (middle panels), an enrichment analysis for GOBP development terms is presented (right panels).</figcaption></figure>

**Surface markers: CDs and other plasma membrane proteins.** Previously, several studies were conducted to catalogue MSC surface markers⁶⁷,⁶⁸, the most comprehensive to date including 116 CD markers found
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Comparison of BM and EM based on significance B. (a) Venn diagrams display PROT SDEs for BM > EM (top) and BM < EM (bottom) by significance B (p value ≤ 0.05) per replicate. Comparison to RNA results based on (b) intersecting proteins or (c) the union. (d) Enrichment analysis on BM vs EM based on significance B analysis of RNA validated candidates for proteins measured by PROT in at least 2 replicates (SigB_i/RNA, left) or 1 replicate (SigB_u/RNA, middle) as well as all significance B SDE proteins measured by PROT (SigB_u, right). Bar charts represent the most significant top 20 terms of GOBP (top) and GOCC (bottom) per cell type and sorted by mean −log₁₀ p values. (e) GOBP development terms specific enrichment analysis for GOBP terms enriched for significance B SDE proteins measured by PROT.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** New potential MSC markers among CD molecules. 37 CD markers were found to be up-regulated in MSC when compared to ESC by RNA, PROT or both (Table S3). Their normalized log₁₀ protein expression levels are shown across a wide range of tissues (top panel) and cell lines (bottom panel) according to ProteomicsDB. CD molecules up-regulated solely in ESC-MSC are highlighted in purple, whereas those up-regulated only in BM-MSC are highlighted in red. CD molecules potentially best suited as MSC markers are comprised in the horizontal rectangle frame.</figcaption></figure>

on bone marrow derived MSC⁶⁷. Using a combination of RNA-seq and quantitative proteomics, we identified 137 surface markers of the CD family (Table S4). As expected, common MSC markers (CD29, CD36, CD44, CD73, CD90, CD105, CD106, CD166) were highly expressed, whereas hematopoietic markers (CD19, CD20, CD34, CD45) were absent. High expression levels of THY-1 (CD90), VCAM1 (CD106), and LNGFR (CD271) on MSC have been associated with self-renewal and robust lineage differentiation⁶⁹ and were proposed as a quality control for clinical use of MSC. Among the SDE proteins validated by both techniques, 9 CD markers were found up-regulated in ESC (Table S4), and 28 CD markers were found up-regulated in both MSC subtypes; 3 and 6 specifically for ESC-MSC (LIFR, PTGFRN, PVR) and BM-MSC (VCAM1, CD151, LAMP1, ANPEP, PRNP, LAMP2), respectively (Fig. 5, Table S4). Among the 28 CD markers enriched in MSC were several tetraspanins known to be exosome markers such as CD9P, CD63, CD81, and CD151⁷⁰. Protein expression profiles of aforementioned MSC CD markers are shown for a broad range of tissues and cell lines accessible through ProteomicsDB⁵² (Fig 5). The potentially best suited MSC markers are TNFSF1A, LIFR, JAG1 and PVRL3, when considering tissues, or TNFRSF10D, TNFRSF1A, CD248 and LIFR, when considering cell lines. In general, most discriminative are TNFRSF1A (CD120a) and TNFRSF10D (CD264), and their use, or more generally our findings, may prove helpful in refining the definition of MSC as given by⁶⁷,⁶⁸.

To broaden our search for potential MSC surface markers, our dataset was further filtered for proteins integral to plasma membrane; 293 candidates identified by nano LC-MS/MS in both MSC subtypes were thus obtained.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** New potential MSC markers among plasma membrane proteins. Proteins identified in MSC (by proteomics) annotated as "integral to plasma membrane" (IPM) were screened for overall tissue expression based on ProteomicsDB. Heatmaps of the 60 most intense IPM proteins for EM and BM are shown, sorted by tissue type and normalized log₁₀ protein expression and separated in 3 panels as follows: IPM proteins 1 to 20 (upper panel), 21 to 40 (middle panel) and 41 to 60 (lower panel). Potential MSC markers are comprised in the horizontal rectangle frames.</figcaption></figure>
Tissue expression profiles for the 60 most intense candidates among those membrane proteins were then gener-
ated using ProteomicsDB⁵² (Fig. 6). Out of those 60 most intense plasma membrane proteins, 20 displayed tissue
specificity, defined as detectable expression levels in a maximum of 4 out of 56 representative tissues. On the basis
of this definition, those membrane proteins could be considered additional potential MSC markers and include:
ACVR1B, BACE1, CD8B, EPHA3, GABRB1, GPR64, IL1R1, LGR5, NRG1, P01762, SLC2A6, SLC7A6, SLC11A2,
SLC16A4, SLC28A1, SLC29A2, SLC38A5, TNFRSF1A, TRO, and TYRO3. Among the non-CD surface markers
with high intensity and the most specific expression (one out of 56 tissues) were ACVR1B, GABRB1, GPR64,
NGR1, P01762 and SLC28A1.

Conclusions

Using a combination of RNA sequencing and quantitative MS-based proteomics, we report here the most com-
prehensive proteome to date for both embryonic and mesenchymal stem cells, quantifying transcripts deriving
from more than 15000 genes and quantitatively profiling more than 8400 proteins.

Both techniques demonstrate high similarity of expression profiles for the two MSC subtypes and their
marked difference to the ESC proteome. RNA sequencing allowed to identify approximately 2500 significantly
differentially expressed (SDE) genes when comparing BM-MSC and ESC-MSC, while MS-based proteomics
reported 40 to 200 SDE proteins between the MSC subtypes (depending on the statistical analysis stringency). It
is worth noting that while the sensitivity advantage of nucleic acid based techniques may lead to overestimation
of differences through emphasis of biological noise, true biological differences may also be masked in mass spec-
trometry by dynamic range limitations. The higher number of SDE features measured by RNA sequencing may
also relate to translational and post-translational regulation not propagating into protein abundance differences.
Nonetheless, a 60% correlation coefficient between the two techniques was achieved, validating literature reports
and confirming the quality of the presented dataset.

Based on PCA and correlation analysis, the presented data suggests equivalency between bone marrow- and
embryonic stem cell-derived MSC, with differences similar to previous reports on donor-derived MSC from
different tissues. This interpretation is further supported by the observation that expression profiles of BM-MSC
from younger donors correlate better to those of ESC derived MSC. The reported differences between ESC and
MSC reflect the expected cell behavior: cell cycle centered for the former versus extracellular processes driven for
the latter. Proteins related to vesicle-mediated transport and exosomes were found to be central in MSC biology,
supporting the paracrine functionality reported by others. The striking enrichment in both MSC subtypes for
developmental terms points to the versatility of MSC, as reflected by their wide range of usage in trials. Apparent
from enrichment analysis, ESC-MSC may support more neurogenesis related processes, whereas BM-MSC more
vasculature development, suggesting the potential use of differentially sourced MSC to address more specific
medical conditions. Overall and keeping in mind that the results may be specific to the isolation and culturing
conditions used, we believe these findings support future evaluation of the clinical use of ESC-derived MSC as a
viable substitute for their bone marrow derived counterparts.

Finally, we reported a total of 137 surface markers of the CD family expressed in ESC and MSC identified jointly in RNA-seq and LC-MS/MS. Among those, 28 were found up-regulated in both MSC subtypes including exosome markers (CD9P, CD63, CD81, and CD151), further confirming the central importance of extracellular processes in MSC biology. Comparing protein expression profiles for the reported MSC CD markers across a variety of cell lines and tissues accessible through ProteomicsDB highlighted potential novel MSC markers: TNFSF1A, TNFRSF10D, LIFR, JAG1, PVRL3 and CD248. Among those, LIFR displays ESC-MSC specificity which may prove of use for quality control towards clinical application. Additionally, we have considered membrane associated proteins as potential cellular markers and have identified 20 tissue specific candidates from the 60 most intense plasma membrane proteins, among which ACVR1B, GABRB1, GPR64, NGR1, P01762 and SLC28A1 were the most promising. We believe MSC definition and, in particular, immunofluorescence-based applications may greatly benefit from the presented list of specific CD markers and membrane associated proteins.

References

1. Hass, R., Kasper, C., Böhm, S. & Jacobs, R. Different populations and sources of human mesenchymal stem cells (MSC): A comparison of adult and neonatal tissue-derived MSC. *Cell Commun. Signal*. doi: 10.1186/1478-811X-9-12 (2011).

2. Le Blanc, K. et al. Mesenchymal stem cells for treatment of steroid-resistant, severe, acute graft-versus-host disease: a phase II study. *The Lancet* **371**, 1579–1586 (2008).

3. Ringdén, O. et al. Mesenchymal stem cells for treatment of therapy-resistant graft-versus-host disease. *Transplantation* **81**, 1390–1397 (2006).

4. Sun, L. et al. Umbilical cord mesenchymal stem cell transplantation in severe and refractory systemic lupus erythematosus. *Arthritis Rheum*. **62**, 2467–2475 (2010).

5. Tyndall, A. & van Laar, J. M. Stem cells in the treatment of inflammatory arthritis. *Best Pract. Res. Clin. Rheumatol*. **24**, 565–574 (2010).

6. Wang, D. et al. Umbilical cord mesenchymal stem cell transplantation in active and refractory systemic lupus erythematosus: a multicenter clinical study. *Arthritis Res. Ther*. **16**, R79 (2014).

7. Hare, J. M. et al. Comparison of allogeneic vs autologous bone marrow-derived mesenchymal stem cells delivered by transendocardial injection in patients with ischemic cardiomyopathy: the POSEIDON randomized trial. *JAMA* **308**, 2369–2379 (2012).

8. Zhang, Z. et al. Human umbilical cord mesenchymal stem cells improve liver function and ascites in decompensated liver cirrhosis patients. *J. Gastroenterol. Hepatol*. **27**, 112–120 (2012).

9. von Bahr, L. et al. Analysis of Tissues Following Mesenchymal Stromal Cell Therapy in Humans Indicates Limited Long-Term Engraftment and No Ectopic Tissue Formation. *STEM CELLS* **30**, 1575–1578 (2012).

10. Prockop, D. J., Kota, D. J., Bazhanov, N. & Reger, R. L. Evolving paradigms for repair of tissues by adult stem/progenitor cells (MSCs). *J. Cell. Mol. Med*. **14**, 2190–2199 (2010).

11. Uccelli, A. & Prockop, D. J. Why should mesenchymal stem cells (MSCs) cure autoimmune diseases? *Curr. Opin. Immunol*. **22**, 768–774 (2010).
12. Friedenstein, A. J., Piatetzky-Shapiro, I. I. & Petrakova, K. V. Osteogenesis in transplants of bone marrow cells. *J. Embryol. Exp. Morphol.* **16**, 381–390 (1966).

13. Kern, S., Eichler, H., Stoeve, J., Klüter, H. & Bieback, K. Comparative Analysis of Mesenchymal Stem Cells from Bone Marrow, Umbilical Cord Blood, or Adipose Tissue. *STEM CELLS* **24**, 1294–1301 (2006).

14. Pevsner-Fischer, M., Levin, S. & Zipori, D. The origins of mesenchymal stromal cell heterogeneity. *Stem Cell Rev.* **7**, 560–568 (2011).

15. Davies, O. G., Cooper, P. R., Shelton, R. M., Smith, A. J. & Scheven, B. A. A comparison of the *in vitro* mineralisation and dentinogenic potential of mesenchymal stem cells derived from adipose tissue, bone marrow and dental pulp. *J. Bone Miner. Metab.* **33**, 371–382 (2014).

16. Friedenstein, A. J., Petrakova, K. V., Kurolesova, A. I. & Frolova, G. P. Heterotopic of bone marrow. Analysis of precursor cells for osteogenic and hematopoietic tissues. *Transplantation* **6**, 230–247 (1968).

17. Pittenger, M. F. et al. Multilineage potential of adult human mesenchymal stem cells. *Science* **284**, 143–147 (1999).

18. Bianco, P., Robey, P. G. & Simmons, P. J. Mesenchymal stem cells: revisiting history, concepts, and assays. *Cell Stem Cell* **2**, 313–319 (2008).

19. Wagner, W. et al. Aging and Replicative Senescence Have Related Effects on Human Stem and Progenitor Cells. *PLoS ONE* doi: 10.1371/journal.pone.0005846 (2009).

20. Boyd, N. L., Robbins, K. R., Dhara, S. K., West, F. D. & Stice, S. L. Human Embryonic Stem Cell-Derived Mesoderm-like Epithelium Transitions to Mesenchymal Progenitor Cells. *Tissue Eng. Part A* **15**, 1897–1907 (2009).

21. Hematti, P. Human embryonic stem cell-derived mesenchymal stromal cells. *Transfusion (Paris)* **51**, 138S–144S (2011).

22. de Peppo, G. M. et al. Human Embryonic Mesodermal Progenitors Highly Resemble Human Mesenchymal Stem Cells and Display High Potential for Tissue Engineering Applications. *Tissue Eng. Part A* **16**, 2161–2182 (2010).

23. Vodyanik, M. A. et al. A Mesoderm-Derived Precursor for Mesenchymal Stem and Endothelial Cells. *Cell Stem Cell* **7**, 718–729 (2010).

24. Kimbrel, E. A. et al. Mesenchymal Stem Cell Population Derived from Human Pluripotent Stem Cells Displays Potent Immunomodulatory and Therapeutic Properties. *Stem Cells Dev.* **23**, 1611–1624 (2014).

25. Sánchez, L. et al. Enrichment of Human ESC-Derived Multipotent Mesenchymal Stem Cells with Immunosuppressive and Anti-Inflammatory Properties Capable to Protect Against Experimental Inflammatory Bowel Disease. *STEM CELLS* **29**, 251–262 (2011).

26. Hajizadeh-Saffar, E. et al. Inducible VEGF Expression by Human Embryonic Stem Cell-Derived Mesenchymal Stromal Cells Reduces the Minimal Islet Mass Required to Reverse Diabetes. *Sci. Rep.* doi: 10.1038/srep09322 (2015).

27. Wang, X. et al. Human ESC-Derived MSCs Outperform Bone Marrow MSCs in the Treatment of an EAE Model of Multiple Sclerosis. *Stem Cell Rep.* **3**, 115–130 (2014).

28. Dominici, M. et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. *Cytotherapy* **8**, 315–317 (2006).

29. Raynaud, C. M. & Rafii, A. The Necessity of a Systematic Approach for the Use of MSCs in the Clinical Setting. *Stem Cells Int.* doi: 10.1155/2013/892340 (2013).

30. Trivedi, P. & Hematti, P. Derivation and immunological characterization of mesenchymal stromal cells from human embryonic stem cells. *Exp. Hematol.* **36**, 350–359 (2008).

31. Raynaud, C. M. et al. Human Embryonic Stem Cell Derived Mesenchymal Progenitors Express Cardiac Markers but Do Not Form Contractile Cardiomyocytes. *PLoS ONE* **8**, e54524 (2013).

32. Ong, S.-E. et al. Stable Isotope Labeling by Amino Acids in Cell Culture, SILAC, as a Simple and Accurate Approach to Expression Proteomics. *Mol. Cell. Proteomics* **1**, 376–386 (2002).

33. Liberski, A. R. et al. Adaptation of a Commonly Used, Chemically Defined Medium for Human Embryonic Stem Cells to Stable Isotope Labeling with Amino Acids in Cell Culture. *J. Proteome Res.* **12**, 3233–3245 (2013).

34. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. *Nat. Protoc.* **7**, 562–578 (2012).

35. Li, H. et al. The Sequence Alignment/Map format and SAMtools. *Bioinforma. Oxf. Engl.* **25**, 2078–2079 (2009).

36. Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. *Bioinformatics* **26**, 841–842 (2010).

37. DeLuca, D. S. et al. RNA-SeQC: RNA-seq metrics for quality control and process optimization. *Bioinformatics* **28**, 1530–1532 (2012).

38. Wessel, D. & Flügge, U. I. A method for the quantitative recovery of protein in dilute solution in the presence of detergents and lipids. *Anal. Biochem.* **138**, 141–143 (1984).

39. Rappsilber, J., Ishihama, Y. & Mann, M. Stop and go extraction tips for matrix-assisted laser desorption/ionization, nanoelectrospray, and LC/MS sample pretreatment in proteomics. *Anal. Chem.* **75**, 663–670 (2003).

40. Cox, J. & Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. *Nat. Biotechnol.* **26**, 1367–1372 (2008).

41. Elias, J. E. & Gygi, S. P. Target-decoy search strategy for mass spectrometry-based proteomics. *Methods Mol. Biol. Clifton NJ* **604**, 55–71 (2010).

42. Vizcaíno, J. A. et al. The Proteomics Identifications (PRIDE) database and associated tools: status in 2013. *Nucleic Acids Res.* **41**, D1063–D1069 (2013).

43. Xie, W. et al. Epigenomic analysis of multilineage differentiation of human embryonic stem cells. *Cell* **153**, 1134–1148 (2013).

44. Miranda, H. C. et al. A quantitative proteomic and transcriptomic comparison of human mesenchymal stem cells from bone marrow and umbilical cord vein. *PROTEOMICS* **12**, 2607–2617 (2012).

45. Roson-Burgo, B., Sanchez-Guijo, F., Del Cañizo, C. & De Las Rivas, J. Transcriptomic portrait of human Mesenchymal Stromal/Stem cells isolated from bone marrow and placenta. *BMC Genomics* **15**, 910 (2014).

46. Tsai, M.-S. et al. Functional Network Analysis of the Transcriptomes of Mesenchymal Stem Cells Derived from Amniotic Fluid, Amniotic Membrane, Cord Blood, and Bone Marrow. *STEM CELLS* **25**, 2511–2523 (2007).

47. Raynaud, C. M. et al. Endothelial cells provide a niche for placental hematopoietic stem/progenitor cell expansion through broad transcriptomic modification. *Stem Cell Res.* **11**, 1074–1090 (2013).

48. Graumann, J. et al. Stable Isotope Labeling by Amino Acids in Cell Culture (SILAC) and Proteome Quantitation of Mouse Embryonic Stem Cells to a Depth of 5,111 Proteins. *Mol. Cell. Proteomics* **7**, 672–683 (2008).

49. Calloni, R., Cordero, E. A. A., Henriques, J. A. P. & Bonatto, D. Reviewing and Updating the Major Molecular Markers for Stem Cells. *Stem Cells Dev.* **22**, 1455–1476 (2013).

50. Aomatsu, E. et al. Novel SCRG1/BST1 axis regulates self-renewal, migration, and osteogenic differentiation potential in mesenchymal stem cells. *Sci. Rep.* **4**, 1–9 (2014).

51. Kim, M.-S. et al. A draft map of the human proteome. *Nature* **509**, 575–581 (2014).

52. Wilhelm, M. et al. Mass-spectrometry-based draft of the human proteome. *Nature* **509**, 582–587 (2014).

53. Phanstiel, D. H. et al. Proteomic and phosphoproteomic comparison of human ES and iPS cells. *Nat. Methods* **8**, 821–827 (2011).

54. Rigbolt, K. T. G. et al. System-Wide Temporal Characterization of the Proteome and Phosphoproteome of Human Embryonic Stem Cell Differentiation. *Sci. Signal.* doi: 10.1126/scisignal.2001570 (2011).

55. Holley, R. J. et al. Comparative Quantification of the Surfaceome of Human Multipotent Mesenchymal Progenitor Cells. *Stem Cell Rep.* **4**, 473–488 (2015).
56. Vaquerizas, J. M., Kummerfeld, S. K., Teichmann, S. A. & Luscombe, N. M. A census of human transcription factors: function, expression and evolution. *Nat. Rev. Genet.* **10**, 252–263 (2009).

57. Aruga, J. The role of Zic genes in neural development. *Mol. Cell. Neurosci.* **26**, 205–221 (2004).

58. Ganss, B. & Jheon, A. Zinc Finger Transcription Factors in Skeletal Development. *Crit. Rev. Oral Biol. Med.* **15**, 282–297 (2004).

59. Okada, K. et al. Identification of SCAN Domain Zinc-Finger Gene ZNF449 as a Novel Factor of Chondrogenesis. *PLoS ONE* doi: 10.1371/journal.pone.0115169 (2014).

60. Hasegawa, R. et al. Identification of ZNF395 as a novel modulator of adipogenesis. *Exp. Cell Res.* **319**, 68–76 (2013).

61. Li, Y. et al. ZNF322, a novel human C2H2 Kruppel-like zinc-finger protein, regulates transcriptional activation in MAPK signaling pathways. *Biochem. Biophys. Res. Commun.* **325**, 1383–1392 (2004).

62. Hemmi, K. et al. A Homeodomain-Zinc Finger Protein, ZFHX4, Is Expressed in Neuronal Differentiation Manner and Suppressed in Muscle Differentiation Manner. *Biol. Pharm. Bull.* **29**, 1830–1835 (2006).

63. Brun, J. et al. The LIM-only protein FHL2 controls mesenchymal cell osteogenic differentiation and bone formation through Wnt5a and Wnt10b. *Bone* **53**, 6–12 (2013).

64. Kadrmas, J. L. & Beckerle, M. C. The LIM domain: from the cytoskeleton to the nucleus. *Nat. Rev. Mol. Cell Biol.* **5**, 920–931 (2004).

65. Benevento, M. et al. Proteome adaptation in cell reprogramming proceeds via distinct transcriptional networks. *Nat. Commun.* doi: 10.1038/ncomms6613 (2014).

66. Griffiths, S., Baraniak, P. R., Copland, I. B., Nerem, R. M. & McDevitt, T. C. Human platelet lysate stimulates high-passage and senescent human multipotent mesenchymal stromal cell growth and rejuvenation in *vitro*. *Cytotherapy* **15**, 1469–1483 (2013).

67. Mindaye, S. T., Ra, M., Lo Surdo, J. L., Bauer, S. R. & Alterman, M. A. Global proteomic signature of undifferentiated human bone marrow stromal cells: Evidence for donor-to-donor proteome heterogeneity. *Stem Cell Res.* **11**, 793–805 (2013).

68. Niehage, C. et al. The Cell Surface Proteome of Human Mesenchymal Stromal Cells. *PLoS ONE* doi: 10.1371/journal.pone.0020399 (2011).

69. Mabuchi, Y. et al. LNGFR+THY-1 + VCAM-1hi + Cells Reveal Functionally Distinct Subpopulations in Mesenchymal Stem Cells. *Stem Cell Rep.* **1**, 152–165 (2013).

70. Andreu, Z. & Yáñez-Mó, M. Tetraspanins in extracellular vesicle formation and function. *Immunother. Vaccines* (doi: 10.3389/fimmu.2014.00442 2014).

## Acknowledgements

J.G. and the Proteomics Core at WCMC-Q are supported by “Biomedical Research Program” funds at Weill Cornell Medical College in Qatar, a program funded by Qatar Foundation. This work was supported by a grant from Qatar National Research Fund’s National Priority Research Program (4–1267–1–194). The statements made herein are solely the responsibility of the authors.

## Author Contributions

J.G., A.R. and K.S. directed the study. A.M.B. designed the experiments. A.M.B. and S.S.D. carried out the experiments. A.M.B. and A.B. analyzed proteomics data and performed bioinformatics analysis. H.B.H. and N.G. measured samples by mass spectrometry. S.H. and P.K. performed R.N.A. sequencing data processing. A.M.B., S.H., P.K. and N.Y. performed statistics. A.M.B., R.J.C. and A.B. created figures. A.M.B., H.B.H. and J.G. drafted the manuscript.

## Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing financial interests: The authors declare no competing financial interests.

**How to cite this article:** Billing, A. M. *et al.* Comprehensive transcriptomic and proteomic characterization of human mesenchymal stem cells reveals source specific cellular markers. *Sci. Rep.* **6**, 21507; doi: 10.1038/srep21507 (2016).

<figure><img src="image_10.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/