Quantitative Proteomic Analysis of Human Embryonic Stem Cell Differentiation by 8-Plex iTRAQ Labelling

Abstract

Analysis of gene expression to define molecular mechanisms and pathways involved in human embryonic stem cells (hESCs) proliferation and differentiations has allowed for further deciphering of the self-renewal and pluripotency characteristics of hESC. Proteins associated with hESCs were discovered through isobaric tags for relative and absolute quantification (iTRAQ). Undifferentiated hESCs and hESCs in different stages of spontaneous differentiation by embryoid body (EB) formation were analyzed. Using the iTRAQ approach, we identified 156 differentially expressed proteins involved in cell proliferation, apoptosis, transcription, translation, mRNA processing, and protein synthesis. Proteins involved in nucleic acid binding, protein synthesis, and integrin signaling were downregulated during differentiation, whereas cytoskeleton proteins were upregulated. The present findings added insight to our understanding of the mechanisms involved in hESC proliferation and differentiation.

About the Authors:

Mahdieh Jadaliha

Contributed equally to this work with: Mahdieh Jadaliha, Hyoung-Joo Lee

Affiliations Department of Molecular Systems Biology, Cell Science Research Center, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran, Iran, Department of Biotechnology, College of Science, University of Tehran, Tehran, Iran

Hyoung-Joo Lee

Contributed equally to this work with: Mahdieh Jadaliha, Hyoung-Joo Lee

Affiliation: Department of Biochemistry, Yonsei Proteome Research Center and Biomedical Proteome Research Center, Yonsei University, Sudaemoon-Ku, Seoul, Korea

Mohammad Pakzad

Affiliation: Department of Stem Cells and Developmental Biology, Cell Science Research Center, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran, Iran

Ali Fathi

Affiliation: Department of Molecular Systems Biology, Cell Science Research Center, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran, Iran

Seul-Ki Jeong

Affiliation: Department of Biochemistry, Yonsei Proteome Research Center and Biomedical Proteome Research Center, Yonsei University, Sudaemoon-Ku, Seoul, Korea

Sang-Yun Cho

Affiliation: Department of Biochemistry, Yonsei Proteome Research Center and Biomedical Proteome Research Center, Yonsei University, Sudaemoon-Ku, Seoul, Korea

Hossein Baharvand

Affiliations Department of Stem Cells and Developmental Biology, Cell Science Research Center, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran, Iran, Department of Developmental Biology, University of Science and Culture, ACECR, Tehran, Iran

Young-Ki Paik

* E-mail: [email protected] (GSH); [email protected] (Y-KP)

Affiliation: Department of Biochemistry, Yonsei Proteome Research Center and Biomedical Proteome Research Center, Yonsei University, Sudaemoon-Ku, Seoul, Korea

Ghasem Hosseini Salekdeh

* E-mail: [email protected] (GSH); [email protected] (Y-KP)

Affiliations Department of Molecular Systems Biology, Cell Science Research Center, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran, Iran, Department of Systems Biology, Agricultural Biotechnology Research Institute of Iran, Karaj, Iran

Introduction

Human embryonic stem cells (hESCs) are pluripotent cells that have the potential to form any cell type and can be propagated in an undifferentiated state in vitro. Their exceptional properties mean they have tremendous potential for developmental biology, drug screening, functional genomics, and regenerative medicine. Developing reliable and reproducible protocols to differentiate hESCs into specific cell types and their transplantation into humans will require a detailed understanding of the molecular mechanisms that maintain the undifferentiated and pluripotent nature of hESCs.

Factors involved in hESCs self-renewal and pluripotency have been described [1], [2], [3], [4], [5], [6], and the overexpression of some of these factors in somatic cells has reprogrammed them into induced pluripotent stem cells (iPSCs), similar to hESCs. However, the genes and mechanisms that maintain the undifferentiated and pluripotent nature of hESCs are still largely unknown.

The proliferation and differentiation of embryonic stem cells (ESCs) are highly coordinated events that involve myriads of genes. A comprehensive profile of ESC gene expression can provide novel insights into the biology of these cells. Transcriptome analysis (i.e., microarray analysis) has shown tremendous potential for the analysis of ESC function and differentiation [7], [8], [9], but gene expression at the transcript level may not correlate well with its expression at the protein level due to alternative splicing, mRNA degradation, and posttranslational modifications such as phosphorylation and protein degradation. These concerns suggest that proteome analysis of ESCs can provide invaluable insights into pathways activated during ESC proliferation and differentiation.

Proteomic tools are valuable in studying ESC differentiation and elucidating the underlying molecular mechanisms [10], [11], [12]. Recent advances in state-of-the-art mass spectrometry (MS) techniques have demonstrated that MS-based quantitative proteomics approaches can significantly contribute to identifying proteins involved in ESC proliferation and differentiation.

Stable isotope labeling with amino acids in cell culture (SILAC) is also a powerful approach for quantitative proteomics [13] that has been used to comprehensively analyze self-renewing versus differentiating cells of two distinct hESC lines [14].

Another popular in vitro labeling method is the isobaric tag for relative and absolute quantitation (iTRAQ) reagent. A 4-plex iTRAQ has been used to analyze hESCs, mouse ES, and EC cells during differentiation by a 4-plex iTRAQ [15], [16]. Neural development from hESCs has been studied using an 8-plex iTRAQ reagent. A study of the progression of neural development from hESCs generated a catalog of approximately 1200 proteins and their relative quantitative expression patterns, which included several that changed expression levels during differentiation [17].

This study applied the 8-plex iTRAQ system to analyze hESCs during the differentiation of embryoid bodies (EB), as this system has the capability to compare several time points during a single experiment. The expression profiles for 1032 proteins during EB differentiation were analyzed and we identified 156 proteins that exhibited statistically significant changes in expression levels during differentiation.

Materials and Methods

Cell Culture

The hESC line, Royan H5, with normal karyotypes (46 XX) at passages 40–50 was used in this experiment. Briefly, the cells were first cultured on mouse embryonic fibroblasts (MEF), inactivated by mitomycin C (Sigma; M0503) [18]. The hESCs were then passaged and maintained under feeder-free conditions for 25–30 passages as described previously [19]. Then, ideal colonies were mechanically dissected into small pieces and replated on matrigel-coated dishes containing hESC medium (2 mM L-glutamine (Gibco; 25030-024), 0.1 mM β–mercaptoethanol (Sigma; M7522), DMEM/F12 medium (Gibco;21331-020) supplemented with 20% knock-out serum replacement (KSR, Gibco; 10828-028), 1% nonessential amino acid (Gibco; 11140-035), 100 units/ml penicillin and 100 µg/ml streptomycin (Gibco; 15070-063), 100 ng/mL basic-fibroblast growth factor (bFGF, Sigma; F0291)). The cells were grown in 5% CO2 at 95% humidity and the hESC medium was changed every day. The cells were further passaged as small clumps (100–500 cells) every 6–7 days after enzymatic treatment (2 mg/ml of dispase) and mechanical dissociation using a cell scraper by gently pipetting. To promote differentiation, hESCs were first cultured in suspension in ESC medium without KSR and containing fetal bovine serum (FBS) (ES-qualified; Gibco 16141-079), where they developed into multicellular aggregates called embryoid bodies (EBs). The EBs were cultured in suspension for 12 days and then plated onto gelatin-coated dishes for 8 days in the same medium to form a pool of spontaneously differentiated cells. Samples from undifferentiated hESCs and EBs at days 6 (EB6), 12 (EB12), and 20 (EB20) were collected for proteomics analysis.

To evaluate the percentage of undifferentiated hESCs, we analyzed the expression of key hESC markers including Nanog, Oct-4, SSEA-4, Tra-1-60, Tra-1-81, and SSEA-4 and Tra-1-60 or Tra-1-81 using two-color flow cytometry as previously described [11]. The analysis was performed by BD-FACS Caliber Flow Cytometer (Becton Dickinson) using following primary antibodies: anti-SSEA-4 (1∶50, ChemiconMAB4304) hOct-4 (1∶50, R&D MAB1759), Tra-1-60 (1∶20, Chemicon MAB4360), Tra-1-81 (1∶20, Chemicon MAB4381) and Nanog (1∶100, R&D MAB1994). Data from three independent replicate were analyzed by WinMDI software (version 2.8). Karyotype analysis and alkaline phosphatase staining was performed as described [11], [20].

Protein Preparation and iTRAQ Isobaric Labeling

Samples of hESCs with at least 6×106 cells in each of the three replicates from hESCs and differentiated derivatives at 6, 12, and 20 days after the initiation of differentiation were homogenized by sonication in 400 µl lysis buffer that consisted of 8 M urea, 4% w/v CHAPS, and one protease inhibitor tablet per 50 ml (Complete Protease Inhibitor Cocktail; Roche, Mannheim, Germany) on ice for 1 minute, with 2-second pulses every 2 seconds. The samples were then vortexed for 30 minutes at room temperature. Insoluble debris were pelleted by centrifugation at about 100,000×g (98,235×g) for 60 minutes at 4°C. The supernatant protein was quantified by the Bradford Assay Kit (BioRad, Hercules, CA) using bovine serum albumin as a standard. A total of 200 µg of each sample was reduced and denatured. The cysteines were then blocked as described in the 8-Plex iTRAQ protocol (Applied BioSystems, Foster City, CA). Each sample was digested with 20 µl of 0.25 µg/µl sequencing-grade modified trypsin solution (1∶20; Promega, Madison, WI) at 37°C, overnight. Samples were dried in a centrifugal vacuum concentrator, reconstituted with 30 µl dissolution buffer, and acidified with 0.1% formic acid to a pH of 2. The peptides were desalted with an Oasis HLB column (Waters, Milford, MA) and labeled with the iTRAQ tags (113–121 m/z): ESC (iTRAQ 113 and 117); EB6 (iTRAQ 114 and 118); EB12 (iTRAQ 115 and 119); and EB20 (iTRAQ 116 and 121). The labeled samples were then dried.

Off-line Strong Cation Exchange (SCX) Chromatography and On-line Nano-LC ESI-MS/MS Analysis

The peptides were separated, complexity was reduced, and all salts and urea were removed by pooling and injecting the acidified samples onto off-line strong cation exchange chromatography (SCX) columns. A total of 20 fractions were collected and dried by speed vacuum. Further separation was achieved by reverse-phase HPLC (Eksigent, Dublin, CA) that interfaced on-line to a QSTAR Elite mass spectrometer (Applied BioSystems). The mass spectrometer was operated in an information-dependent acquisition mode, whereby, following the interrogation of MS data (m/z 350–2000) using a 1-second survey scan, ions were selected for MS/MS analysis based on their intensity (>20 cpm) and charge state (+2, +3, and +4). Four product ion scans were set from each survey scan. Statistical evidence of differential expression of proteins was gained by performing two experiments for each biological replicate as described above (Figure 1).

[Figure omitted. See PDF.]

Figure 1. Experimental design of proteome analysis of hESCs using iTRAQ labeling.

Samples from undifferentiated hESCs and EBs at days 6, 12, and 20 were collected in three biological replicates. Similar amounts of proteins were digested into peptides using trypsin. Peptides were subsequently desalted and labeled with 8-plex iTRAQ reagents 113–121. Labeled peptides were pooled, fractionated into 20 SCX fractions, and then analyzed by reverse-phase LC-MS/MS.

https://doi.org/10.1371/journal.pone.0038532.g001

Data analysis and Interpretation

Relative abundances were quantitated, and peptides and proteins were identified with Protein-Pilot TM Software 2.0 (Applied BioSystems). Each MS/MS spectrum was searched for the Homo sapiens species in the UniProt database. The relative amount of a peptide in each sample was calculated by dividing the peak areas observed at 114.1, 115.1, and 116.1 m/z by those observed at 113.1 and 118.1 m/z, while those observed at 119.1 and 121.1 m/z were divided by that observed at 117.1 m/z. The peptides that did not have an iTRAQ modification were excluded. In addition, to avoid protein inference problem for proteins with high degree of sequence similarities, "Shared peptides", similar peptide sequence belongs to more than one protein reported in the results were excluded. The logarithm of each ratio was evaluated by a one-sample, unpaired t-test to determine statistical significance. Total significant proteins were clustered by the k-means clustering method with MATLAB version 7.3. The number of correct clusters was determined by measuring the average of intracluster and intercluster distances based on the similarity of a gene to the genes in its own cluster as compared to genes in other clusters [21], [22].

Western Blot Analysis

Samples of 50 µg of proteins from three biological replicates were separated by 12% SDS-PAGE electrophoresis at 120 V for 1 hour with a Mini-PROTEAN 3 electrophoresis cell (Bio-Rad). The proteins were transferred to a PVDF membrane (Amersham, Uppsala, Sweden) by semi-dry blotting (Bio-Rad) with Dunn carbonate transfer buffer (10 mM NaCHO3, 3 mM Na2CO3, 20% methanol). The membranes were blocked for 1.5 hours using a western blocker solution (Sigma, St. Louis, MO, W0138). Each membrane was incubated overnight at 4°C with one of these primary monoclonal antibodies: anti-ERP29 (1∶4000; Abcam); anti-NPM1 (1∶1000; Sigma); anti-HSC70 (1∶10000; Stressgen); anti-CALU (1∶4000; Santa Cruz Biotechnology, Santa Cruz, CA); or anti-STMN1 (1∶2000; Abcam). Next, the membranes were incubated for 2 hours at room temperature with the following peroxidase-conjugated secondary antibodies: anti-mouse (1∶180,000; Sigma, A9044); anti-rat (1∶160,000; Sigma, A5795), and anti-rabbit (1∶160,000; Sigma, A2074). Finally, the blots were visualized with ECL detection reagent (Sigma, CPS-1-120). The films were then scanned with a GS-800 densitometer (Bio-Rad), and quantitative analysis was performed with UVI bandmap software (UVItec, Cambridge, UK). The uniformity of the amounts of proteins loaded on the gels was investigated by staining the membranes with Fast Green (Sigma, F7252).

Results and Discussion

Cell Characterization

The hESCs grew as compact colonies with a high nuclear to cytoplasmic ratio and prominent nucleoli (Figure S1A and S1B). The colony showed a typical undifferentiated morphology with a distinct boundary, and each cell presents a compact morphology with a high nucleus to cytoplasmic ratio, that contains prominent nucleoli typical of undifferentiated hESCs. Moreover, the hESC line had a normal karyotype (46 XX) (Figure S1 C). To induce differentiation, hESCs were cultured as EBs (Figure S1 D-G). The high percentage of undifferentiated hESCs was confirmed by expression analysis of key hESC markers including Nanog, Oct-4, SSEA-4, Tra-1-60, Tra-1-81 (Figure S1 H and Table S1).

Proteome Analysis

We used an 8-plex iTRAQ system to analyze the proteomes of hESCs during proliferation and at different stages of differentiation in three biological replicates. For each biological replicate, the experiments were repeated twice. Protein lysates from ESCs and EBs on days 6, 12, and 20 after initiation of differentiation were used for iTRAQ labeling. These proteins represented the different stages of spontaneous differentiation. The proteins were labeled with eight different iTRAQ reagents (113–121) as follows: iTRAQ 113 and 117 for ESC; iTRAQ114 and118 for EB6; iTRAQ115 and 119 for EB12; and iTRAQ 116 and 121 for EB20. Labeled proteins were then pooled and fractionized by SCX. A total of 20 SCX fractions were analyzed by reverse-phase LC-MS/MS as described in Materials and Methods. MS/MS spectra were searched and quantitated.

A total of 2153 unique proteins were identified, of which 1648 proteins with quantitation ratios were considered for further analyses. By combining six replicates, 1032 nonredundant proteins were confidently identified and quantified by the criterion of unused protein score <1.3 (95% CI) per experiment. The complete list of all proteins identified from the three biologically independent replicates and the confidence scores is provided in Table S2.

Proteins with statistically significant changes during differentiation were identified by filtering according to these criteria: 1) they had to be present in at three replicates, including two biological replicates; 2) changes between stages had to be statistically significant (P<0.05); and 3) fold change had to be greater than 1.2. This approach allowed us to select 159 differentially expressed proteins for further analysis.

Out of 156 differentially expressed proteins, 81 were found to be downregulated and 79 upregulated in at least one EB compared to ESCs (Table S2 and S3). These proteins could be clustered into six different groups (Figure 2, Table S4), which represented downregulated proteins in clusters 1–3 and upregulated proteins in clusters 4–6. In cluster 1, the abundance of 37 proteins decreased significantly at all three stages. In cluster 2, 30 proteins were downregulated at EB20 (late stage of differentiation) and included such proteins as Lin28, the widely used marker of pluripotency. In cluster 3, the expression levels of 13 protein decreased at EB6 and EB12, however these did not change significantly at EB20. In cluster 4, 22 proteins were upregulated mainly in EB6 and EB12, with more pronounced upregulation in EB12. In cluster 5, 28 proteins significantly increased in abundance at EB20 but did not change at EB6 and EB12. In cluster 6, 29 proteins had upregulated expression in all three stages, but was more pronounced at EB20.

[Figure omitted. See PDF.]

Figure 2. K-mean clusters of differentially expressed proteins.

These proteins could be clustered into six different groups with clusters 1–3 representing downregulated proteins and clusters 4–6 representing upregulated proteins. In cluster 1, the abundance of 37 proteins decreased significantly at all three stages. In cluster 2, 30 proteins were downregulated at EB20, the late stage of differentiation. In cluster 3, the expression levels of 13 protein decreased at EB6 and EB12. In cluster 4, 22 proteins showed higher abundance in EB6 and EB12. In cluster 5, 28 proteins significantly increased in abundance at EB20 and in cluster 6, 29 proteins had higher abundance in all three stages that was more pronounced at EB20.

https://doi.org/10.1371/journal.pone.0038532.g002

Functional analysis of expression clusters with Panther software (www.pantherdb.org) revealed that proteins involved in nucleic acid binding, protein synthesis, and signaling, (particularly integrin signaling), were enriched in the downregulated protein clusters (clusters 1, 2, and 3), whereas the upregulated proteins in clusters 4, 5, and 6 were enriched in proteins involved in cytoskeleton structure (Figure 3). These proteins are listed in Table S5.

[Figure omitted. See PDF.]

Figure 3. The functional analysis of 6 k-mean clusters.

This list of proteins was employed to identify significantly activated pathways by comparing their functional annotations according to the PANTHER classification systems. Proteins involved in nucleic acid binding, protein synthesis, and signaling (particularly, proteins involved in integrin signaling) were enriched in the downregulated protein clusters (clusters 1, 2, and 3), whereas the upregulated proteins in clusters 4, 5, and 6 were enriched in proteins involved in cytoskeleton structure.

https://doi.org/10.1371/journal.pone.0038532.g003

iTRAQ Results Confirmed by Western Blot Analysis

The iTRAQ results were confirmed by Western blot analysis, which examined the expression levels of CALU, ERP31, NPM1, HSC70, and STMN1 (Figure 4). The levels of HSC70, ERP29, and NPM1 decreased during differentiation while those of CALU and STMN1 increased. Western blot analysis confirmed the results of the iTRAQ analysis, although the fold changes obtained by the two approaches differed slightly from each other (Figure 4).

[Figure omitted. See PDF.]

Figure 4. Western blot analysis of the total protein derived from Royan H5.

Fifty micrograms of protein from three biological replicates extracted from three independent replications of hESCs and Dif-ESCs of Royan H5 were subjected to SDS-PAGE followed by Western blotting. ESCs and Dif-ESCs were analyzed with antibodies against CALU, ERP31, NPM1, HSC70, and STMN1. The y-axis represents the area in Western blot for different stages (x-axis). Protein bands were quantified using UVI bandmap software. Error bars represent the standard deviation of three measurements.

https://doi.org/10.1371/journal.pone.0038532.g004

ESCs Express Relatively High Levels of Nucleic Acid Binding Proteins

Proteins downregulated during differentiation included the nucleic acid binding proteins DEAD box RNA helicase (DDX5), scaffold attachment factor B1 (SAFB1), HIST1H2AJ, HNRNPA2B1, HNRPA1, HNRPAB, nucleolin (NCL), RBMX, SSBP1, adenosine deaminase (ADAR1), HNRPDL, insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3), LIN28, PSPC1, SNRPD2, XRN2, and the widely used marker of pluripotency, LIN28. The downregulated protein IGF2BP3 is expressed mainly in embryonic development and in some tumors, and has been demonstrated to promote cell proliferation by inducing translation of IGF-II mRNA in K562 leukemia cells [15]. The downregulated protein NCL is a major nucleolar protein that plays a role in many pathways and functions [23]. Downregulation of NCL has been shown to be detrimental to the growth of ESCs and to increase the rate of apoptosis, which suggests its importance in maintaining the self-renewal of ESCs [24]. Downregulated ADAR1 is responsible for RNA editing by site-specific deamination of adenosines [25]. Knockdown of ADAR1 globally affects gene expression in hESCs and results in significantly increased RNA expression levels of genes involved in differentiation and developmental processes, including neurogenesis [26]. ADAR1 has been suggested to play a role in regulating hESC early differentiation [26]. Members of the heterogeneous nuclear ribonucleoprotein (HNRNP) family, which included HNRNPA2B1, HNRPA1, HNRPAB, HNRPDL, HNRNPI, and RBMX were downregulated, which was consistent with several studies [10], [11], [12].

Several identified proteins are multifunctional and contain numerous highly conserved functional domains. The SAFB1 protein can bind both DNA and RNA and is involved in RNA processing and stress response. It contains a transcriptional repression domain, can bind certain hormone receptors to repress their activity [27], and may also be involved in development, growth regulation, and reproduction [28]. DDX5 protein regulates aspects of RNA expression, including replication, transcription and splicing, and is required for cell growth [29]. It has been demonstrated to be a novel co-activator for Runx2, and can inhibit the osteogenic differentiation of mesenchymal progenitor cells [30].

EBs Increase Expression of Cytoskeleton-associated Proteins

The proteomes of EBs and ESCs exhibited differences in the expressions of cytoskeleton-pathway proteins. A total of 20 out of 26 cytoskeleton-associated proteins showed higher expression levels in EBs than ESCs. The altered proteins included CORO1C, GSN, MSN, ACTG1, CALD1, KRT18, KRT7, MYH9, MYL6, MYL9, PDLIM7, PDLIM1, TPM4, DSP, KRT19, KRT8, MAP1B, RDX, SNL, and VIL2. Cytoskeleton-associated proteins had downregulated expression in mouse, monkey, and human EB-mediated differentiation of ESCs [11], [12], [31]. The cell shape has been suggested to be a cue in the commitment process [32], and changes in cell shape may be transduced into a regulatory signal by several structures in the cell, including the actin cytoskeleton itself [33]. Mechanical tension can control the differentiation status of adult stromal stem cells through the actin filament complex [32].

Regulation of Calcium-binding Proteins

Calcium-binding proteins that included CANX, RCN2, MCP, ANXA5, GSN, CALU, MYL6, MYL9, RCN, S100A10, and ANXA1 were significantly regulated. Calcium-binding proteins, such as calreticulin, pyruvate dehydrogenase complex, and the translationally-controlled tumor protein are regulated during neural differentiation of mouse ESCs [34], implying that Ca2+ may play an important role in hESC differentiation. The Ca2+ ion is a highly versatile intracellular signal that regulates many different cellular functions, including fertilization, cell cycle, apoptosis, muscle contraction, vision, and memory.

Intracellular Ca2+ homeostasis is maintained through a complex interplay between intracellular stores such as the endoplasmic reticulum and extracellular Ca2+ that enters the cell through various transporters on the plasma membrane. Calcium regulates gene expression by modulating transcription factors [35] and mediates posttranslational modifications by modulating protein kinases, phosphatases, and Ca2+-sensitive adenylate cyclases [36]. Calcium modulates intercellular communication through gap junctions and triggers the terminal differentiation programs of cells. However, the details of Ca2+ homeostasis and signaling, and the mechanism by which Ca2+ regulates Ca2+-related proteins during ESC differentiation remain to be determined.

Ribosomal Proteins Downregulated

Ribosomal proteins, including RPL23, RPS24, RPL7, RPS7, RPLP2, and RPS19, were more abundant in undifferentiated ESCs than in differentiated cells. This finding has been supported by the observation that ESCs possess a surplus of free ribosomes as ribosomal subunits and single ribosomes that are recruited to actively translating polysomes during differentiation [37]. The protein synthesis capacity of ESCs has been suggested to be poised to allow rapid elevation of translation rate in response to differentiation signals [37].

Proteins Involved in Integrin Signaling were more Abundant in Undifferentiated ESCs

Extracellular matrix (ECM) signaling is predominantly transmitted via cell membrane receptors of the integrin family [38]. The spatially and temporally controlled engagement of different integrin receptors during embryogenesis demonstrates their roles during commitment and lineage determination of early embryogenesis [39].

Several proteins involved in integrin signaling including CD49B, FLNA, ACTN4, ITGB1, FLN2, and EDS4A were more abundant in ESCs than in EBs. An in vitro loss-of-function approach based on b1 integrin-deficient ESCs found that integrin-dependent mechanisms were involved in the regulation of Wnt-1 and BMP-4 expression [40]. Integrin signaling engineered in mouse ESCs demonstrated the critical role of simultaneous signaling from identified integrins in maintaining pluripotency [41].

[Figure omitted. See PDF.]

Table 1. Correlation between the expression patterns of differentially expressed proteins identified in this study and their corresponding mRNA at three different stages.

https://doi.org/10.1371/journal.pone.0038532.t001

The Correlation between mRNA and Protein Levels

Transcriptomics techniques (i.e., microarray) are powerful approaches to profile mRNA expression. The extent that changing mRNA expression patterns reflect corresponding changes in their cognate proteins is an important issue. We compared the expression patterns of proteins with the level of their corresponding mRNA previously analyzed by Fathi et al [11] were analyzed to determine the correspondence between protein and mRNA levels (Table 1 and Figure S2). Pearson’s method measured the correlations between 156 differentially expressed proteins and their corresponding mRNAs. The correlation between RNA and protein abundance levels was low (Figure S2). A trend of correlation between protein and mRNA changes was observed, but with many exceptions (Table 1).

Conclusion

We applied an iTRAQ-based quantitative proteomics approach to characterize changes in the proteome pattern of hESCs during differentiation. Our findings add insight to the understanding of the mechanisms involved in hESC proliferation and differentiation. Large-scale iTRAQ labeling methods are limited because combining labeled samples from different stages may result in peptides that have been derived from different stages dominating over those specifically found in only one or two stages. This limits their identification in regular data-dependent acquisition of highly abundant ions, such as several major transcription factors (Oct-4, Nanog, and SOX2) not identified by this study. Possibly, these proteins are highly expressed only in ESCs, and the peptides from these proteins become diluted upon mixing labeled peptides from several stages of differentiation. This work provides experimental evidence of the contributions of several candidate proteins and mechanisms in the differentiation of hESCs, but the function of these proteins in self-renewal and differentiation still needs to be precisely clarified.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Characterization of undifferentiated hESCs (Royan H5) and differentiated EBs at different time points. (A) Phase contrast photographs of a Royan H5 colony grown under feeder-free conditions and its high magnification photo. (B) Expression of alkaline phosphatase. (C) The karyotype of Royan H5. (D–G) EB formation by generating some clamps from undifferentiated hESCs at day 0 (D), day 6 (E), day 12 (F), and day 20. (G) After 12 days in suspension, EBs were plated on 0.1% gelatin-coated plates in the same medium to form a pool of spontaneously differentiated cells. The percentages of undifferentiated and differentiated hESCs are shown in (H) by the use of BD-FACS Caliber flow cytometry (Becton Dickinson). Data from three independent replicates were analyzed by WinMDI software (version 2.8). Bar = 500 µm.

https://doi.org/10.1371/journal.pone.0038532.s001

(DOCX)

Figure S2.

The heat map and Pearson correlation of relative mRNA and protein abundance for stages EB6/ESC, EB12/ESC, and EB20/ESC. The correlation was calculated for 184 mRNA, which paired 156 significantly changed proteins (there were more than one mRNA data for some proteins). The heat map is divided into 6 blocks, considering 6 different k-mean groups.

https://doi.org/10.1371/journal.pone.0038532.s002

(DOCX)

Table S1.

Flow cytometric analysis of three replicates regarding five hESC markers, Oct4, Nanog, SSEA-4, Tra 1-60 and Tra 1-81. Three replicates are showing similar patterns and stem cell related proteins were being down-regulated during EB formation.

https://doi.org/10.1371/journal.pone.0038532.s003

(DOCX)

Table S2.

The list of identified differentially expressed proteins in different stages of EB6, EB12 and EB20 compared to ESC.

https://doi.org/10.1371/journal.pone.0038532.s004

(DOCX)

Table S3.

Quantification details of different repeats (for all and regulated proteins). The quantification data for different stages and different replicates (shown in yellow highlighted columns) has been used for t-test analysis and resulted in identification of regulated proteins considering three criteria mentioned before. This table includes: N which presents the rank of a protein relative to all other proteins in the list of detected proteins in each repeat; Total ProtScore is an indicator of the total amount of evidence for a detected protein. The Total ProtScore is calculated using all of the peptides detected for the proteins, and does not show the confidence percentage of the identification of a protein. Unused ProtScore is an indicator of the protein confidence for a detected protein which is calculated from the peptide confidence for peptides from spectra that have not already been used by other proteins. %Cov represents the percentage of matching amino acids to at least one identified peptide having confidence greater than 0 divided by the total number of amino acids in the protein sequence. %Cov(50) and %Cov(95) are the percentage of matching amino acids to at least one identified peptide having confidence greater than or equal to 50 and 95, respectively, divided by the total number of amino acids in the protein sequence. Ratio is the average ratio for the protein in EBs compared to ESCs, which is corrected for experimental bias. The iTRAQ tags are iTRAQ 113 and 117 for ESC, iTRAQ 114 and 118 for EB6, iTRAQ 115 and 119 for EB12 and iTRAQ 116 and 121 EB20. P-value is a measure of the certainty that the average ratio differs from one (this p-value is for each repeat and is different from t-test analysis for all 6 replicates which has been done to detect regulated proteins with shown data in regulated proteins sheet). The error factor (EF) is a measure of the error in the average ratio, representing the 95% confidence interval of the average iTRAQ ratio as (ratio×EF) – (ratio/EF).

https://doi.org/10.1371/journal.pone.0038532.s005

(XLS)

Table S4.

List of the proteins depicted in the figure 2, with their expression ratio and p-values.

https://doi.org/10.1371/journal.pone.0038532.s006

(DOCX)

Table S5.

Classification of differentially expressed proteins according their protein classes. Table contains the class of each protein and also protein classes enriched in the every k-means clusters.

https://doi.org/10.1371/journal.pone.0038532.s007

(XLS)

Author Contributions

Conceived and designed the experiments: HB YKP GHS. Performed the experiments: MJ HJL MP AF. Analyzed the data: MJ SKJ SYC. Wrote the paper: YKP GHS.

Citation: Jadaliha M, Lee H-J, Pakzad M, Fathi A, Jeong S-K, Cho S-Y, et al. (2012) Quantitative Proteomic Analysis of Human Embryonic Stem Cell Differentiation by 8-Plex iTRAQ Labelling. PLoS ONE7(6): e38532. https://doi.org/10.1371/journal.pone.0038532

References

1. Chambers I, Silva J, Colby D, Nichols J, Nijmeijer B, et al. (2007) Nanog safeguards pluripotency and mediates germline development. Nature 450: 1230–1234.I. ChambersJ. SilvaD. ColbyJ. NicholsB. Nijmeijer2007Nanog safeguards pluripotency and mediates germline development.Nature45012301234

2. Chambers I, Smith A (2004) Self-renewal of teratocarcinoma and embryonic stem cells. Oncogene 23: 7150–7160.I. ChambersA. Smith2004Self-renewal of teratocarcinoma and embryonic stem cells.Oncogene2371507160

3. Molofsky AV, Pardal R, Iwashita T, Park IK, Clarke MF, et al. (2003) Bmi-1 dependence distinguishes neural stem cell self-renewal from progenitor proliferation. Nature 425: 962–967.AV MolofskyR. PardalT. IwashitaIK ParkMF Clarke2003Bmi-1 dependence distinguishes neural stem cell self-renewal from progenitor proliferation.Nature425962967

4. Nishino J, Kim I, Chada K, Morrison SJ (2008) Hmga2 promotes neural stem cell self-renewal in young but not old mice by reducing p16Ink4a and p19Arf Expression. Cell 135: 227–239.J. NishinoI. KimK. ChadaSJ Morrison2008Hmga2 promotes neural stem cell self-renewal in young but not old mice by reducing p16Ink4a and p19Arf Expression.Cell135227239

5. Chambers I, Colby D, Robertson M, Nichols J, Lee S, et al. (2003) Functional expression cloning of Nanog, a pluripotency sustaining factor in embryonic stem cells. Cell 113: 643–655.I. ChambersD. ColbyM. RobertsonJ. NicholsS. Lee2003Functional expression cloning of Nanog, a pluripotency sustaining factor in embryonic stem cells.Cell113643655

6. Boyer LA, Lee TI, Cole MF, Johnstone SE, Levine SS, et al. (2005) Core transcriptional regulatory circuitry in human embryonic stem cells. Cell 122: 947–956.LA BoyerTI LeeMF ColeSE JohnstoneSS Levine2005Core transcriptional regulatory circuitry in human embryonic stem cells.Cell122947956

7. Ivanova NB, Dimos JT, Schaniel C, Hackney JA, Moore KA, et al. (2002) A stem cell molecular signature. Science 298: 601–604.NB IvanovaJT DimosC. SchanielJA HackneyKA Moore2002A stem cell molecular signature.Science298601604

8. Bhattacharya B, Cai J, Luo Y, Miura T, Mejido J, et al. (2005) Comparison of the gene expression profile of undifferentiated human embryonic stem cell lines and differentiating embryoid bodies. BMC Dev Biol 5: 22.B. BhattacharyaJ. CaiY. LuoT. MiuraJ. Mejido2005Comparison of the gene expression profile of undifferentiated human embryonic stem cell lines and differentiating embryoid bodies.BMC Dev Biol522

9. Ramalho-Santos M, Yoon S, Matsuzaki Y, Mulligan RC, Melton DA (2002) "Stemness": transcriptional profiling of embryonic and adult stem cells. Science 298: 597–600.M. Ramalho-SantosS. YoonY. MatsuzakiRC MulliganDA Melton2002"Stemness": transcriptional profiling of embryonic and adult stem cells.Science298597600

10. Kurisaki A, Hamazaki TS, Okabayashi K, Iida T, Nishine T, et al. (2005) Chromatin-related proteins in pluripotent mouse embryonic stem cells are downregulated after removal of leukemia inhibitory factor. Biochem Biophys Res Commun 335: 667–675.A. KurisakiTS HamazakiK. OkabayashiT. IidaT. Nishine2005Chromatin-related proteins in pluripotent mouse embryonic stem cells are downregulated after removal of leukemia inhibitory factor.Biochem Biophys Res Commun335667675

11. Fathi A, Pakzad M, Taei A, Brink TC, Pirhaji L, et al. (2009) Comparative proteome and transcriptome analyses of embryonic stem cells during embryoid body-based differentiation. Proteomics 9: 4859–4870.A. FathiM. PakzadA. TaeiTC BrinkL. Pirhaji2009Comparative proteome and transcriptome analyses of embryonic stem cells during embryoid body-based differentiation.Proteomics948594870

12. Nasrabadi D, Rezaei Larijani M, Pirhaji L, Gourabi H, Shahverdi A, et al. (2009) Proteomic analysis of monkey embryonic stem cell during differentiation. J Proteome Res 8: 1527–1539.D. NasrabadiM. Rezaei LarijaniL. PirhajiH. GourabiA. Shahverdi2009Proteomic analysis of monkey embryonic stem cell during differentiation.J Proteome Res815271539

13. Ong SE, Mann M (2005) Mass spectrometry-based proteomics turns quantitative. Nat Chem Biol 1: 252–262.SE OngM. Mann2005Mass spectrometry-based proteomics turns quantitative.Nat Chem Biol1252262

14. Prokhorova TA, Rigbolt KT, Johansen PT, Henningsen J, Kratchmarova I, et al. (2009) Stable isotope labeling by amino acids in cell culture (SILAC) and quantitative comparison of the membrane proteomes of self-renewing and differentiating human embryonic stem cells. Mol Cell Proteomics 8: 959–970.TA ProkhorovaKT RigboltPT JohansenJ. HenningsenI. Kratchmarova2009Stable isotope labeling by amino acids in cell culture (SILAC) and quantitative comparison of the membrane proteomes of self-renewing and differentiating human embryonic stem cells.Mol Cell Proteomics8959970

15. Yocum AK, Gratsch TE, Leff N, Strahler JR, Hunter CL, et al. (2008) Coupled global and targeted proteomics of human embryonic stem cells during induced differentiation. Mol Cell Proteomics 7: 750–767.AK YocumTE GratschN. LeffJR StrahlerCL Hunter2008Coupled global and targeted proteomics of human embryonic stem cells during induced differentiation.Mol Cell Proteomics7750767

16. O'Brien RN, Shen Z, Tachikawa K, Lee PA, Briggs SP (2010) Quantitative proteome analysis of pluripotent cells by iTRAQ mass tagging reveals post-transcriptional regulation of proteins required for ES cell self-renewal. Mol Cell Proteomics 9: 2238–2251.RN O'BrienZ. ShenK. TachikawaPA LeeSP Briggs2010Quantitative proteome analysis of pluripotent cells by iTRAQ mass tagging reveals post-transcriptional regulation of proteins required for ES cell self-renewal.Mol Cell Proteomics922382251

17. Chaerkady R, Kerr CL, Marimuthu A, Kelkar DS, Kashyap MK, et al. (2009) Temporal analysis of neural differentiation using quantitative proteomics. J Proteome Res 8: 1315–1326.R. ChaerkadyCL KerrA. MarimuthuDS KelkarMK Kashyap2009Temporal analysis of neural differentiation using quantitative proteomics.J Proteome Res813151326

18. Baharvand H, Ashtiani SK, Taee A, Massumi M, Valojerdi MR, et al. (2006) Generation of new human embryonic stem cell lines with diploid and triploid karyotypes. Development, Growth and Differentiation 48: 117–128.H. BaharvandSK AshtianiA. TaeeM. MassumiMR Valojerdi2006Generation of new human embryonic stem cell lines with diploid and triploid karyotypes.Development, Growth and Differentiation48117128

19. Levenstein ME, Ludwig TE, Xu RH, Llanas RA, VanDenHeuvel-Kramer K, et al. (2006) Basic Fibroblast Growth Factor Support of Human Embryonic Stem Cell Self-Renewal. Stem Cells 24: 568.ME LevensteinTE LudwigRH XuRA LlanasK. VanDenHeuvel-Kramer2006Basic Fibroblast Growth Factor Support of Human Embryonic Stem Cell Self-Renewal.Stem Cells24568

20. Larijani MR, Seifinejad A, Pournasr B, Hajihoseini V, Hassani SN, et al. (2011) Long-term maintenance of undifferentiated human embryonic and induced pluripotent stem cells in suspension. Stem Cells and Development 20: 1911–1923.MR LarijaniA. SeifinejadB. PournasrV. HajihoseiniSN Hassani2011Long-term maintenance of undifferentiated human embryonic and induced pluripotent stem cells in suspension.Stem Cells and Development2019111923

21. Kaufman L, Rousseeuw PJ (1990) Finding groups in data. an introduction to cluster analysis. Wiley Series in Probability and Mathematical Statistics Applied Probability and Statistics, New York: Wiley, 1990. L. KaufmanPJ Rousseeuw1990Finding groups in data. an introduction to cluster analysis.Wiley Series in Probability and Mathematical Statistics Applied Probability and Statistics, New York: Wiley, 1990

22. Ray S, Turi RH (1999) Determination of number of clusters in k-means clustering and application in colour image segmentation. Proceedings of the 4th International Conference on Advances in Pattern Recognition and Digital Techniques (ICAPRDT'99), Calcutta, India. pp. 137–143.S. RayRH Turi1999Determination of number of clusters in k-means clustering and application in colour image segmentation.Proceedings of the 4th International Conference on Advances in Pattern Recognition and Digital Techniques (ICAPRDT'99), Calcutta, India137143

23. Srivastava M, Pollard HB (1999) Molecular dissection of nucleolin's role in growth and cell proliferation: new insights. FASEB J 13: 1911–1922.M. SrivastavaHB Pollard1999Molecular dissection of nucleolin's role in growth and cell proliferation: new insights.FASEB J1319111922

24. Li H, Wang B, Yang A, Lu R, Wang W, et al. (2009) Ly-1 antibody reactive clone is an important nucleolar protein for control of self-renewal and differentiation in embryonic stem cells. Stem Cells 27: 1244–1254.H. LiB. WangA. YangR. LuW. Wang2009Ly-1 antibody reactive clone is an important nucleolar protein for control of self-renewal and differentiation in embryonic stem cells.Stem Cells2712441254

25. Nishikura K (2006) Editor meets silencer: crosstalk between RNA editing and RNA interference. Nat Rev Mol Cell Biol 7: 919–931.K. Nishikura2006Editor meets silencer: crosstalk between RNA editing and RNA interference.Nat Rev Mol Cell Biol7919931

26. Osenberg S, Paz Yaacov N, Safran M, Moshkovitz S, Shtrichman R, et al. (2010) Alu sequences in undifferentiated human embryonic stem cells display high levels of A-to-I RNA editing. PLoS One 5: e11173.S. OsenbergN. Paz YaacovM. SafranS. MoshkovitzR. Shtrichman2010Alu sequences in undifferentiated human embryonic stem cells display high levels of A-to-I RNA editing.PLoS One5e11173

27. Oesterreich S (2003) Scaffold attachment factors SAFB1 and SAFB2: Innocent bystanders or critical players in breast tumorigenesis? J Cell Biochem 90: 653–661.S. Oesterreich2003Scaffold attachment factors SAFB1 and SAFB2: Innocent bystanders or critical players in breast tumorigenesis?J Cell Biochem90653661

28. Ivanova M, Dobrzycka KM, Jiang S, Michaelis K, Meyer R, et al. (2005) Scaffold attachment factor B1 functions in development, growth, and reproduction. Mol Cell Biol 25: 2995–3006.M. IvanovaKM DobrzyckaS. JiangK. MichaelisR. Meyer2005Scaffold attachment factor B1 functions in development, growth, and reproduction.Mol Cell Biol2529953006

29. Ford MJ, Anton IA, Lane DP (1988) Nuclear protein with sequence homology to translation initiation factor eIF-4A. Nature 332: 736–738.MJ FordIA AntonDP Lane1988Nuclear protein with sequence homology to translation initiation factor eIF-4A.Nature332736738

30. Jensen ED, Niu L, Caretti G, Nicol SM, Teplyuk N, et al. (2008) p68 (Ddx5) interacts with Runx2 and regulates osteoblast differentiation. J Cell Biochem 103: 1438–1451.ED JensenL. NiuG. CarettiSM NicolN. Teplyuk2008p68 (Ddx5) interacts with Runx2 and regulates osteoblast differentiation.J Cell Biochem10314381451

31. Baharvand H, Fathi A, Gourabi H, Mollamohammadi S, Salekdeh GH (2008) Identification of mouse embryonic stem cell-associated proteins. J Proteome Res 7: 412–423.H. BaharvandA. FathiH. GourabiS. MollamohammadiGH Salekdeh2008Identification of mouse embryonic stem cell-associated proteins.J Proteome Res7412423

32. McBeath R, Pirone DM, Nelson CM, Bhadriraju K, Chen CS (2004) Cell Shape, Cytoskeletal Tension, and RhoA Regulate Stem Cell Lineage Commitment. Developmental Cell 6: 483–495.R. McBeathDM PironeCM NelsonK. BhadrirajuCS Chen2004Cell Shape, Cytoskeletal Tension, and RhoA Regulate Stem Cell Lineage Commitment.Developmental Cell6483495

33. Hynes RO, Huang S, Chen CS, Ingber DE (1998) Control of Cyclin D1, p27Kip1, and Cell Cycle Progression in Human Capillary Endothelial Cells by Cell Shape and Cytoskeletal Tension. Molecular Biology of the Cell 9: 3179–3193.RO HynesS. HuangCS ChenDE Ingber1998Control of Cyclin D1, p27Kip1, and Cell Cycle Progression in Human Capillary Endothelial Cells by Cell Shape and Cytoskeletal Tension.Molecular Biology of the Cell931793193

34. Wang D, Gao L (2005) Proteomic analysis of neural differentiation of mouse embryonic stem cells. Proteomics 5: 4414–4426.D. WangL. Gao2005Proteomic analysis of neural differentiation of mouse embryonic stem cells.Proteomics544144426

35. Berridge MJ, Bootman MD, Roderick HL (2003) Calcium signalling: dynamics, homeostasis and remodelling. Nat Rev Mol Cell Biol 4: 517–529.MJ BerridgeMD BootmanHL Roderick2003Calcium signalling: dynamics, homeostasis and remodelling.Nat Rev Mol Cell Biol4517529

36. Yamaguchi T, Wallace DP, Magenheimer BS, Hempson SJ, Grantham JJ, et al. (2004) Calcium restriction allows cAMP activation of the B-Raf/ERK pathway, switching cells to a cAMP-dependent growth-stimulated phenotype. J Biol Chem 279: 40419–40430.T. YamaguchiDP WallaceBS MagenheimerSJ HempsonJJ Grantham2004Calcium restriction allows cAMP activation of the B-Raf/ERK pathway, switching cells to a cAMP-dependent growth-stimulated phenotype.J Biol Chem2794041940430

37. Sampath P, Pritchard DK, Pabon L, Reinecke H, Schwartz SM, et al. (2008) A hierarchical network controls protein translation during murine embryonic stem cell self-renewal and differentiation. Cell Stem Cell 2: 448–460.P. SampathDK PritchardL. PabonH. ReineckeSM Schwartz2008A hierarchical network controls protein translation during murine embryonic stem cell self-renewal and differentiation.Cell Stem Cell2448460

38. Hynes RO (1992) Integrins: versatility, modulation, and signaling in cell adhesion. Cell 69: 11–25.RO Hynes1992Integrins: versatility, modulation, and signaling in cell adhesion.Cell691125

39. Fassler R, Georges-Labouesse E, Hirsch E (1996) Genetic analyses of integrin function in mice. Curr Opin Cell Biol 8: 641–646.R. FasslerE. Georges-LabouesseE. Hirsch1996Genetic analyses of integrin function in mice.Curr Opin Cell Biol8641646

40. Czyz J, Wobus A (2001) Embryonic stem cell differentiation: the role of extracellular factors. Differentiation 68: 167–174.J. CzyzA. Wobus2001Embryonic stem cell differentiation: the role of extracellular factors.Differentiation68167174

41. Lee ST, Yun JI, Jo YS, Mochizuki M, van der Vlies AJ, et al. (2010) Engineering integrin signaling for promoting embryonic stem cell self-renewal in a precisely defined niche. Biomaterials 31: 1219–1226.ST LeeJI YunYS JoM. MochizukiAJ van der Vlies2010Engineering integrin signaling for promoting embryonic stem cell self-renewal in a precisely defined niche.Biomaterials3112191226