Dynamic proteomic profiling of human periodontal ligament stem cells during osteogenic differentiation
Stem Cell Research & Therapy volume 12, Article number: 98 (2021)
Human periodontal ligament stem cells (hPDLSCs) are ideal seed cells for periodontal regeneration. A greater understanding of the dynamic protein profiles during osteogenic differentiation contributed to the improvement of periodontal regeneration tissue engineering.
Tandem Mass Tag quantitative proteomics was utilized to reveal the temporal protein expression pattern during osteogenic differentiation of hPDLSCs on days 0, 3, 7 and 14. Differentially expressed proteins (DEPs) were clustered and functional annotated by Gene Ontology (GO) terms. Pathway enrichment analysis was performed based on the Kyoto Encyclopedia of Genes and Genomes database, followed by the predicted activation using Ingenuity Pathway Analysis software. Interaction networks of redox-sensitive signalling pathways and oxidative phosphorylation (OXPHOS) were conducted and the hub protein SOD2 was validated with western blotting.
A total of 1024 DEPs were identified and clustered in 5 distinctive clusters representing dynamic tendencies. The GO enrichment results indicated that proteins with different tendencies show different functions. Pathway enrichment analysis found that OXPHOS was significantly involved, which further predicted continuous activation. Redox-sensitive signalling pathways with dynamic activation status showed associations with OXPHOS to various degrees, especially the sirtuin signalling pathway. SOD2, an important component of the sirtuin pathway, displays a persistent increase during osteogenesis. Data are available via ProteomeXchange with identifier PXD020908.
This is the first in-depth dynamic proteomic analysis of osteogenic differentiation of hPDLSCs. It demonstrated a dynamic regulatory mechanism of hPDLSC osteogenesis and might provide a new perspective for research on periodontal regeneration.
Since they were first defined, human periodontal ligament stem cells (hPDLSCs) have been identified as the optimal seed cells for periodontal regeneration [1, 2]. They exhibit osteoblast-like characteristics and are capable of differentiating into cementoblasts or osteoblasts, showing great potential for the regeneration of functional periodontal tissues. Although efforts have been made in the past to promote the efficacy of osteogenesis of hPDLSCs, the clinical stem cell therapy effects are imprecise . The exact molecular mechanism and signalling pathways that regulate differentiation potency remain to be deciphered.
Previous studies have conducted high-throughput techniques to reveal gene expression alterations in hPDLSCs at the transcriptional level. Regrettably, they focused only on regulatory non-coding RNAs rather than functionally translational genes, or focused only on a single time point (day 7 or 14 of osteogenesis duration) and the final phenotype [4,5,6,7,8]. Osteogenesis is a dynamic cellular differentiation process that consists of distinct developmental phases including proliferation and matrix synthesis/maturation/mineralization stages, indicating a phasic pattern of osteogenesis-related protein synthesis . Zheng et al. investigated gene expression profiles altered continuously during the osteogenesis course of hPDLSCs and found some distinct expression patterns, such as constitutive downregulation of epigenetic regulators and dip patterns of osteoblast-associated genes . Although protein abundances somehow scale with mRNA levels, this association was limitative to an extent and mRNA levels should not be interpreted as the final output of gene expression [10, 11]. This addresses the need to investigate the execution of cell differentiation with proteomic techniques, which reveal the regulatory mechanism at the translational level.
To date, only two reports have examined the proteomic profile during the osteogenic differentiation of hPDLSCs. As early as 2009, Xiao et al. compared the protein profile between osteogenic-induced and non-induced hPDLSCs using gel-based mass spectrometry . They found that differentially expressed proteins were mainly related to the cytoskeleton, nuclear regulation, cell membrane binding, matrix synthesis, metabolic enzymes and signal transduction. Despite a new perspective on the osteogenic differentiation progress of hPDLSCs, the drawbacks of gel-based proteomic techniques, including low sensitivity, poor separation and resolution, undermine their application in accurate proteomic research . Our team recently carried out research to explore the role of integrin alpha 5 in the osteogenic differentiation of hPDLSCs using a label-based technique . This study involved protein profiling of late 14-day osteogenic differentiation to some extent, but it did not conduct an in-depth exploration of the dynamic protein expression pattern.
To address this question comprehensively, we conducted the first dynamic proteome analysis to reveal the temporal protein profile of hPDLSCs during the osteogenic differentiation. In the present study, we found that proteins involved in this process presented dynamic expression patterns and functioned distinctively. Oxidative phosphorylation (OXPHOS), the most significant pathway, was demonstrated to be continuously activated during the process. It was associated with several redox-sensitive signalling pathways (such as sirtuin and AMPK), which also showed altered activation/inhibition status. Furthermore, SOD2 might play a vital role in the osteogenesis of hPDLSCs.
Cell culture and osteogenic differentiation protocols
hPDLCs were obtained from three periodontally healthy donors (aged 20–25, molars or premolars) who requested tooth extraction for orthodontic treatment at the Department of Stomatology, Nanfang Hospital, Southern Medical University. This study was approved by the Ethics Committee of Nanfang Hospital, Southern Medical University (NFEC-2020-253), and informed consent was obtained from donors. The middle third of the periodontium was carefully scraped off with a surgical blade. After washing with PBS containing 100 U/mL penicillin/streptomycin, the periodontium was minced to approximately 1 mm2 per piece and scattered on the culture flask and maintained in low-glucose Dulbecco’s modified Eagle’s medium (Gibco, Invitrogen, NY, USA) with 10% foetal bovine serum (Gibco, Invitrogen, NY, USA) and 100 U/mL penicillin/streptomycin (Gibco, Invitrogen, NY, USA). After reaching 80% confluency, hPDLSCs were obtained utilizing the limiting dilution method, featuring clone formation capacity. Passages 4–6 of hPDLSCs were used in the subsequent experiments. We applied the chemically osteogenic differentiation protocol to induce osteogenesis of hPDLSCs. Briefly, hPDLSCs at a confluency of 80–90% were induced to differentiate using growth culture media supplemented with 10 mM beta-glycerophosphate (Sigma, St Louis, MO, USA), 50 μg/mL l-ascorbic acid-2-phosphate (WAKO, Osaka, Japan) and 100 nM dexamethasone (Sigma, St Louis, MO, USA). The osteogenic differentiation culture medium was changed every 3 days. Alkaline phosphatase (ALP) and alizarin red S (ARS) staining were performed to verify the establishment of the osteogenesis phenotype (details in supplementary materials).
Cell sample preparation and Tandem Mass Tag (TMT) proteomics analysis
The detailed materials and methods are described in supplementary materials.
Proteins in each differentiated time point were compared to any other time point, and those with fold change > 1.2 or < 0.83 and p value < 0.05 determined by one-way ANOVA followed by post hoc Tukey HSD test were defined as differentially expressed proteins (DEPs). Correlations among different time points were calculated in Prism using Pearson correlation. To address the temporal dynamics of osteogenic differentiation of hPDLSCs, we subjected the DEPs to the k-means clustering algorithm assigning the optimal cluster number based on the reckon of sum of the squared errors. Functional enrichment analyses were subsequently performed for each cluster. Gene Ontology (GO) annotations were retrieved from the GO database (http://www.geneontology.org/), and hypergeometric test was used to find significantly enriched GO terms based on GO::TermFinder . Pathway enrichment analysis was conducted in Metascape (http://metascape.org/gp/index.html) with the ontology source of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway using the p value set at 0.9 . Ingenuity Pathway Analysis (IPA) was used to predict the activation/inhibition scores of canonical pathways . The protein-protein interaction (PPI) network was built with the STRING online tool (https://string-db.org/) and visualized using Cytoscape. The minimum required score for PPI network construction was set to medium confidence (0.400) and the interaction derived from all provided sources including textmining, experiments, database, co-expression, neighbourhood, gene fusion and co-occurrence.
qRT-PCR and western blotting analysis
Total RNA and proteins of osteogenic hPDLSC samples at days 0, 3, 7 and 14 were extracted and tested. Please see details in supplementary materials.
GraphPad Prism was used to perform one-way ANOVA with post hoc Tukey HSD for comparisons among more than two groups. Significance was determined at p < 0.05. Data are presented as the mean ± SD. Band intensity in western blot images was quantified with ImageJ Software, and values are expressed as the means ± SD of at least three independent experiments. All experiments were performed at least in triplicate.
Summary of TMT proteomics analysis
To describe the proteome dynamics during osteogenic differentiation of hPDLSCs, we used the classical induction protocol to induce osteogenic-differentiated hPDLSCs (Fig. 1a). During the 14-day differentiation, the ALP activity continuously increased with a rapid increase at day 7, and the mineralized nodules were obvious at day 14 in the bright field as well as after ARS staining (Fig. 1b). Consistent with the morphological changes, osteogenesis-related proteins (RUNX2 and ALP) showed dynamic changes at the protein level (Fig. 1c). We then conducted TMT proteome analyses at days 0, 3, 7 and 14. A total of 7012 proteins with at least one unique peptide sequence and 1% false discovery rate confidence were identified, and the protein expression profiles at 4 time points are presented in Fig. 2a. We performed subsequent analysis on a set of 5806 proteins quantified across all time points with at least two replicates (variations of quantification data among biological samples were shown in Fig. S1). As shown in Fig. 2b, the correlation coefficients by Pearson correlation analysis were low between different time points, indicating that the protein profile among different osteogenesis phases was significantly distinct.
A total of 1024 DEPs were identified, and the number of DEPs across all time points is illustrated in Fig. 2c (details in Table S1). Proteins differed the most between the undifferentiated stage (D0) and terminal differentiated stage (D14). Among the 1024 DEPs, proteins associated with osteogenesis are listed in Table 1, including ALP and collagen proteins (COL1A1, COL1A2, COL3A1 and COL6A1) as well as some transduction signalling molecules (PDGFRA and TGFBR1). For FBN1, we have previously validated its involvement in osteogenic/odontongenic differentiation of dental pulp stem cells .
Dynamic expression and functional characterization of DEPs
To better monitor the temporal changes in protein expression, we subjected the 1024 DEPs to the k-means clustering algorithm. The optimal number of clusters was 5 (data not shown). Therefore, the quantitative temporal profiles partitioned into 5 clusters (Fig. 3a). The largest cluster was cluster 2, with 354 proteins whose expression continuously declined from day 0. The second largest cluster was cluster 4, with 346 proteins whose expression gradually increased during the whole differentiation course. Meanwhile, the expression of 80 proteins in cluster 5 changed similarly to that of proteins in cluster 4 but showed a more rapid increase tendency. Cluster 1 was composed of 143 proteins that peaked at day 3 and showed a small decrease or remained thereafter. Cluster 3 consisted of 101 proteins whose expression remained at the similarly level during the course.
Considering that proteins with different expression patterns may perform different functions, we tested the respective enrichment of GO terms. Generally, the GO enrichment results indicated that proteins with different tendencies show different functions. The results of GO enrichment are illustrated in the form of heatmaps (Fig. 3b–d, details in Table S2). Cluster 3, with proteins that changed little during the differentiation course, was not significantly enriched in any functional categories. Meanwhile, clusters 4 and 5 seemed to share similar terms to some extent.
Cluster 2, with proteins that continuously declined upon initiation of osteogenic differentiation, was enriched in DNA replication categories (Fig. 3b). This biological process (BP) term consisted of proteins involved in cell cycle progression, such as the cyclin-dependent kinase CDK9 and the minichromosome maintenance (MCM) protein family, which occupy a central role in DNA synthesis . Another commonly used proliferative marker PCNA and embryonic stem cell marker SSRP1  were also found in cluster 2. Cluster 1, where proteins increased to peak at day 3 and then declined slightly or maintained the expression level, was enriched in the extracellular matrix (ECM) and structure organization. The proteins involved included the osteogenic-related proteins annexin A2 (ANXA2), fibrillin-2 (FBN2) and fibulins (FBLN1 and FBLN5). For cluster 4, mitochondrial function, cellular respiratory process, OXPHOS, secretion function and catabolic process were obviously enriched. Meanwhile, the BP function of proteins from cluster 5 overlapped somewhat with that of cluster 4, but cluster 5 may also participate in cellular energy metabolism. Interestingly, the functions of ECM organization and structure were enriched in clusters 1, 2 and 5, representing different alteration tendencies. COL6A1, COL6A2, COL6A3 and COL18A1 belonged to cluster 5, while COL8A1 and COL11A1 were from cluster 1. COL1A1, COL1A2, COL3A1, COL5A1 and COL14A1 were from cluster 2.
When referring to molecular function GO analysis (Fig. 3c), functions distinguished more among different clusters. Cluster 1 describes the early osteogenesis phase and is involved in cell adhesion, molecular adaptors, structural molecules and ECM-related constituents. Cluster 2 described proteins that bind to diverse macromolecules, including RNA, DNA, growth factor and collagen. Clusters 4 and 5 both described proteins involved in the activity of oxidoreductase, catalytic and electron transfer. Cluster 4 also represented proteins associated with the activity of mitochondrial function, while proteins related to ECM structural constituent and nucleoside triphosphate adenylate kinase activity were for cluster 5.
In Cellular Component GO analysis (Fig. 3d), proteins from clusters 1 and 2 shared the same categories (extracellular organelle and cell junction/adhesion) to some degree, despite they representing different alteration tendencies. This was also the case between clusters 4 and 5, where proteins involved in mitochondrial constituents and secretory transporters presented an upregulation during the course of osteogenesis of hPDLSCs either in a moderate or rapid trend.
In summary, GO enrichment demonstrated that proteins associated with ECM organization/interaction, secretory function and mitochondrial activity played functions in a complex regulatory manner during osteogenesis of hPDLSCs.
OXPHOS was the most enriched pathway in KEGG analysis
A total of 1022 out of 1024 DEPs were mapped to the KEGG database using Metascape (http://metascape.org/gp/index.html), and their functional annotations were classified (Fig. 4a, b). The second largest pathway classification was “Metabolism”, annotated to 50 metabolic pathways occupying 21% (347) of 1022 DEPs. As shown in Fig. 4b, most of the DEPs (174) were clustered in the “Signal transduction” classification at level 2, annotated to 20 specific signalling pathways including osteogenesis-related pathways such as PI3K-Akt, HIF-1, mTOR, FoxO, TGF-β, AMPK, MAPK and Wnt signalling pathways (detailed enrichment results in Table S3). The top enriched pathways included metabolism-related pathways, biogenesis of proteins and amino acids, DNA replication, organelles, endocytosis and ECM interactions (Fig. 4c). It is worth noting that “oxidative phosphorylation” term was the most significant pathway.
Dynamic activation of redox-sensitive signalling pathways revealed by IPA
Apart from the pathways in which the DEPs participated, we also examined whether these pathways were suppressed or activated during the course. Here, we used IPA to identify activation or inhibition of the functional pathways. By comparing one time point with the previous one, we observed that various pathways showed altered activated/inhibited status. A total of 75 Ingenuity Canonical Pathways had valid data at least two time points (Table S4).
Undoubtedly, “oxidative phosphorylation” canonical pathway was also enriched in the IPA database with continuous activation in all three stages (Fig. 5a, b). The upregulations of 5 mitochondrial complexes were verified with western blotting (Fig. 5c). Noticeably, UQCRC2 (from complex III subunit), SDHB (from complex II subunit) and NDUFB8 (from complex I subunit) showed a remarkable increase during osteogenesis. Combined with the apparent effect of mitochondria revealed by GO analysis, we reasoned into the probable vital role of OXPHOS inside the mitochondria during the osteogenic differentiation of hPDLSCs. Increased mitochondrial activity and high levels of OXPHOS are inevitably accompanied by reactive oxygen species (ROS) by-products, which are thought to be primarily generated in mitochondria . ROS can regulate a diverse array of physiological processes via oxidation of signalling molecules [51, 52]. In addition, molecular signal transduction plays a crucial role in the osteogenic differentiation of stem cells . Thus, we further examined the time course changes in the activation profiles of the ROS-regulated signalling pathway, namely, redox-sensitive signalling pathways in IPA.
Based on a literature search, 6 redox-sensitive signalling pathways were closely related to osteogenesis. These pathways presented a dynamic activity profile (Fig. 5a and the detailed map in Fig. S2, S3, S4, S5, S6, S7). Wnt signalling was predicted to be non-activated at day 3 but activated at day 14. Sphingosine-1-phosphate (S1P) and sirtuin signalling pathways were predicted to be inhibited during the course, where the S1P signalling pathway was inhibited at days 3 and 14, while the sirtuin signalling pathway was suppressed at the early and mid-term course and showed no inhibition at the terminal stage. The ERK/MAPK signalling pathway represented activation status at days 3 and 14; however, the status of day 7 was not indicated in our data. Actin cytoskeleton signalling, a cytoplasmic structure-associated pathway, was predicted to be activated at the early osteogenesis stage but inhibited at the mid-term. Similarly, such dynamic changes could also be observed in the AMPK signalling pathway with activation at days 3 and 14 but inhibition at day 7.
From the above, redox-sensitive signalling transduction pathways might play roles in osteogenesis of hPDLSCs in a temporal manner, in line with the dynamic change in this cellular process.
The interactions between OXPHOS and redox-sensitive signalling pathways
To reveal the potential correlations between OXPHOS and redox-sensitive signalling pathways, a PPI network was utilized. As shown in Fig. S8, all six pathways were related to OXPHOS to varying degrees: 21 edges for sirtuin, 15 edges for AMPK, 7 edges for Actin cytoskeleton, 3 edges for S1P and 2 edges for both ERK/MAPK and Wnt/β-catenin signalling. The sirtuin signalling pathway not only correlated most but also shared the most components with OXPHOS. Illustrated with a detailed network (Fig. 6a), we found that SOD2, a main antioxidant, showed a high connection degree of 12 edges in the network. Additionally, among all eight antioxidants detected in our data, SOD2 showed the most remarkable upregulation within the course (Fig. 6b). Hence, we validated SOD2 mRNA and protein expression and found an upregulation starting at day 3 and a remarkable increase at day 14, consistent with the TMT data (Fig. 6c, d).
Osteogenic differentiation is a progressive physiological process, and related genes/proteins are expressed in a temporal and coordinating manner, which cannot be fully evaluated by limited morphological, cytochemical and biochemical testing. The present study was motivated by questions on the temporal dynamics of protein profiling during the osteogenic differentiation of hPDLSCs. Using label-based quantitative proteomics, which allowed the simultaneous profiling of approximately 7000 proteins of diverse time points, we first presented a large-scale dynamic proteomic profiling of this process. With clustering and functional annotation, we characterized the proteins from different expression patterns with distinct functions. Subsequently, we mapped activation profiles of canonical signalling pathways, which revealed the interactions between OXPHOS and redox-sensitive signalling pathways. Furthermore, SOD2, a continuously upregulated antioxidant enzyme, might coordinate the sirtuin signalling pathway and OXPHOS during the osteogenic differentiation process.
Clustering the proteins along with functional annotation helped us to identify the distinct functions of different protein classes. The quantitative temporal profiles were partitioned into 5 clusters. In the bone tissue lifespan, ECM is not a static but dynamic microenvironment with innate complexities, making it important in homeostasis and regeneration spatially and temporally [54, 55]. In our study, different types of ECM showed different expression patterns. Surprisingly, COL1A1 and COL1A2, along with other type I distribution-associated collagens, including COL3A1, COL5A1 and COL14A1 , continuously decreased. This was inconsistent with many other studies reporting increased COL1 mRNA expression during osteogenesis of MSCs. Nevertheless, a negative relationship was found between collagen mRNA levels and cell layer collagen protein accumulation in differentiated pre-osteoblasts . Collagen synthesis occurs early in the life of these cultures, whereas insoluble collagen deposition occurs later accompanied by a diminished rate of collagen synthesis. In the presence of ascorbic acid, which facilitates osteogenic differentiation by increasing collagen type 1 secretion , collagen type I protein was secreted and assembled into extracellular fibrils in a time-dependent manner while the intracellular collagen type I protein decreased . Along with the knowledge of the high intrinsic expression levels of COL1 mRNA in hPDLSCs [60,61,62], we suspected that the intrinsic intracellular collagen type I protein might be adequate for osteogenic differentiation and secretory function rather than protein synthesis playing a major role in this cellular differentiation process. However, further studies concerning the association between the intrinsic and secreted proteins during this process are needed.
COL6s, in contrast to COL1s, were increased throughout the course. In human osteoblast SaM-1 cells, the decline of COL6A1 at the mRNA and protein levels decreased COL1A1 mRNA in the early phase of mineralization . Another ECM FBLN1 was reported to be indispensable for the good organization of collagen fibre bundles but did not affect its expression in mouse PDLs . These results suggested the potential association between different types of ECM expressionally or functionally. The altered ECM protein pattern in our study coincided with the notion that the ECM is constantly undergoing remodelling. Nonetheless, to investigate the exact roles of ECM in the osteogenesis of hPDLSCs, methods for the detection of secretory proteins and ECM depositions should be conducted, not merely intracellular protein quantifications .
Compared to our results in PDLSCs, the ECM of MSC from the bone marrow presented an upregulation during osteogenesis, despite the detection at mRNA levels . It is a pity that so far, no study has investigated the dynamic protein expression patterns of BMSC during osteogenic differentiation. Comprehensive investigations should be carried out to answer the question whether these dynamic ECM expression patterns in MSCs were distinctive among origins.
Commonly, the mitochondria maintained a relatively low activity level in undifferentiated MSCs while activated in differentiated cells . Upon osteogenic induction of hBMSCs, mitochondrial biogenesis was intensified where the copy number of mitochondrial DNA, protein subunits of the respiratory enzymes and intracellular ATP content were increased . In our study, mitochondrial function was also activated, suggested by the upregulation of proteins functioned in mitochondrial activity and cellular energy-producing catabolic processes. Single-stranded DNA-binding protein, which is implicated in replication and maintenance of mitochondrial DNA , belonged to cluster 4 and continuously increased. The respiratory electron transport chain (ETC), located on the mitochondrial inner membrane, is composed of complexes I–IV of the respiratory chain and complex V, an ATP synthase, driving ATP synthesis . The ETC-related proteins all presented an increasing tendency in our study, indicating the high activity of the cellular respiration and the energy production. The well-known function of mitochondria is to provide energy for cell life activities through OXPHOS, which occurs in the mitochondrial cristae . This was also detected in our study, where the proteins from GO terms of OXPHOS, metabolic process and inner mitochondrial membrane complex tended to be upregulated.
The function of OXPHOS was further clarified in the KEGG pathway and IPA analyses. The high level of OXPHOS along with alterations in the distribution, number, mass and morphology of mitochondria are important characteristics during the osteogenic differentiation of MSCs , which corresponds with the increased protein expressions of OXPHOS complexes in western blotting validation. Commonly, stimulation of osteogenic differentiation seems to suppress adipogenic differentiation and this biological basis is increasingly implicated by the tight regulation by ROS, the by-products of OXPHOS. A study has revealed that the intracellular levels of ROS were dramatically reduced during the first 2 days of osteogenic induction but as osteogenic differentiation continued, gradual rebound ROS was noted and higher than that of undifferentiated hMSCs at the terminal stage of day 28 . However, in another study, the ROS level in BMSCs was significantly continuously increased during osteogenesis . Contradictories exist and the ROS level cannot be predicted by the activity of OXHPOS; hence, how the ROS level is altered during the osteogenic differentiation of hPDLSCs is urgently needed in future studies.
An increasing number of studies have paid more attention to the potential involvement of ROS in the signalling pathways [73, 74], i.e. redox-sensitive signalling pathways. Using IPA, we noticed that several redox-sensitive signalling pathways presented dynamic activities during the osteogenic differentiation of hPDLSCs as well as the connections with OXPHOS. Among these, the sirtuin signalling pathway was associated the most, in which SOD2 was one of the hub proteins. From the PPI, we noticed that numerous mitochondrial proteins belonged to both the pathways, indicating the intrinsic relationship between them. Sirtuins belong to the class III histone deacetylase (Hdac) family, named for their homology of the yeast protein silent information regulator (Sir)2 . Emerging studies have discussed the functions of different sirtuins in the field of osteogenesis during the last decade. SIRT1 mediated the enhancement of osteogenic differentiation by exogenous antioxidants [76, 77]. Upregulation of SIRT3 was found during the antioxidation by melatonin in osteogenesis along with strengthening the activity of SOD2 . Similarly, SIRT6 was also positively related to the osteogenic phenotype, in which downregulation of SIRT6 could inhibit osteogenesis of murine BMSCs, while overexpression could reverse it . However, downregulation of SIRT7 improved the osteogenic differentiation of hBMSCs by activating β-catenin . In our study, inhibited activity was found on day 3 and day 7, while non-activation was predicted in the terminal state. The function of sirtuins in osteogenesis and bone metabolism has just come into pictures, and the complexity and dual effects of sirtuins make it worthy of in-depth study. Although the sirtuin proteins did not show altered expression in our data, the downstream targets and other proteins participating in the sirtuin signalling pathway presented dynamic expression. FOXO3, a downstream protein of SIRT1, presented an upregulation at an early stage but maintained its expression after that in our study. A study found that in addition to the increasing expression of FOXO3a, overexpressed SIRT1 also deacetylated FOXO3a, thus leading to the restoration of osteogenesis , indicating that there might be a modulation of FOXO3 at the transcriptional and posttranscriptional levels in our study. SOD2 was the prominent antioxidant enzyme during hPDLSC osteogenesis and was also identified as an upregulation pattern in the osteogenesis of stem cells of other types . The levels of antioxidants are crucial for the correct lineage differentiation and osteogenesis potency of MSCs by modulating the antioxidative ability, redox-sensitive signalling pathway and mitochondrial function/activity . SOD2 was reported to be transcriptionally regulated by FOXO3 in the presence of dysregulated ROS , and there was a positive association between FOXO3 and SOD2 at the protein level via a Sirt1-mediated pathway during osteoblastic bone formation . The specific role of SOD2 and its exact relationship with FOXO3a in osteogenesis-related redox-sensitive pathways are worthy of further investigation.
Two signalling pathways showed dual dynamic activity in our study, namely, actin cytoskeleton signalling and AMPK signalling. The impaired migration ability of MSCs was associated with excessive cellular ROS which led to overpolymerization of the F-actin cytoskeleton . In the field of neuronal development and trafficking, the regulation of the cytoskeleton dynamics by differential ROS levels has been investigated . The dynamic actin cytoskeleton activity might be ROS-involved during the osteogenic differentiation of hPDLSCs. The AMPK signalling pathway, a widely studied pathway in cellular redox regulation, is activated by the increase in cellular ROS and conversely, it can also decrease ROS levels by increasing transcriptional regulation of antioxidant defences, including catalase and SOD2 via activation of PGC-1α and FOXO3 . A previous study  has demonstrated that AMPK activation was required at early differentiation but downregulation of AMPK at the mid-late period was also necessary. This was consistent with our findings despite the different durations of differentiation. Activation of AMPK at an early stage augments autophagy and leads to the initiation of differentiation, while accumulation of autophagosome components could impair differentiation, which could be downregulated by AMPK blockage. This suggests the necessity of dynamic signalling pathways as for their potential dual effects in different differentiated stages.
Other redox-associated signalling pathways, including Wnt/β-catenin signalling , S1P signalling  and the ERK/MAPK signalling pathway , also presented dynamic activity in our study and might regulate osteogenic differentiation in a temporal-dependent manner, where the activation status of the ERK/MAPK signalling pathway has been verified in a previous study .
At last, from our original data with respect to correlations among all samples, the intra-group correlation was not as high as expected (Fig. 2b). This might be attributed to the intrinsically donor-to-donor and intra-population heterogeneity, which often weakens the reproduction of experimental results in a certain degree [91, 92]. This suggested that we might need to expand the sample size in the future for further extrapolation of our conclusions. The sample size of the current study should be sufficient to draw a preliminary conclusion; however, we should take the sample size into account when it comes to primary mesenchymal stem cells.
Collectively, with the use of a large-scale quantitative proteomic technique, we first described the dynamic protein expression pattern during osteogenic differentiation of hPDLSCs. Biological functions differed among various clusters of proteins with different alteration tendencies. Pathway enrichment along with the activation prediction suggested that OXPHOS might play a vital role in the process, and the sirtuin signalling pathway was closely associated with OXPHOS. Furthermore, SOD2, associating the aforementioned two pathways, was confirmed to be continuously upregulated, representing the activation of the antioxidant system in the osteogenesis of hPDLSCs. Our comprehensive proteomics analysis demonstrated a dynamic regulatory mechanism during the osteogenic differentiation of hPDLSCs and might provide a new perspective for research on periodontal regeneration.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Adenosine 5′-monophosphate-activated protein kinase
Alizarin red S
Cyclin-dependent kinase 9
Collagen type XI alpha 1 chain
Collagen type XIV alpha 1 chain
Collagen type XVIII alpha 1 chain
Collagen type I alpha 1 chain
Collagen type I alpha 2 chain
Collagen type I alpha 3 chain
Collagen type V alpha 1 chain
Collagen type VI alpha 1 chain
Collagen type VI alpha 2 chain
Collagen type VI alpha 3 chain
Differentially expressed proteins
Extracellular regulated protein kinases
Electron transport chain
Forkhead box O
Hypoxia-inducible factor 1
Human periodontal ligament stem cells
Ingenuity Pathway Analysis
Kyoto Encyclopedia of Genes and Genomes
Mitogen-activated protein kinase
Mechanistic target of rapamycin kinase
Nicotinamide adenine dinucleotide
Proliferating cell nuclear antigen
Platelet-derived growth factor receptor alpha
Phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha
Reactive oxygen species
RUNX family transcription factor 2
Silent information regulator
Superoxide dismutase 2
Structure-specific recognition protein 1
Transforming growth factor
Transforming growth factor beta receptor 1
Tissue inhibitor of matrix metalloproteinase
Tandem Mass Tag
Seo BM, Miura M, Gronthos S, Bartold PM, Batouli S, Brahim J, et al. Investigation of multipotent postnatal stem cells from human periodontal ligament. Lancet. 2004;364:149–55.
Qiu W, Wu B-L, Fang F-C. Overview of noncoding RNAs involved in the osteogenic differentiation of periodontal ligament stem cells. World J Stem Cells. 2020;12:251–65.
Chen F-M, Gao L-N, Tian B-M, Zhang X-Y, Zhang Y-J, Dong G-Y, et al. Treatment of periodontal intrabony defects using autologous periodontal ligament stem cells: a randomized clinical trial. Stem Cell Res Ther. 2016;7:33.
Gu X, Li M, Jin Y, Liu D, Wei F. Identification and integrated analysis of differentially expressed lncRNAs and circRNAs reveal the potential ceRNA networks during PDLSC osteogenic differentiation. BMC Genet. 2017;18:100.
Hao Y, Ge Y, Li J, Hu Y, Wu B, Fang F. Identification of microRNAs by microarray analysis and prediction of target genes involved in osteogenic differentiation of human periodontal ligament stem cells. J Periodontol. 2017;88:1105–13.
Qu Q, Fang F, Wu B, Hu Y, Chen M, Deng Z, et al. Potential role of long non-coding RNA in osteogenic differentiation of human periodontal ligament stem cells. J Periodontol. 2016;87:e127–37.
Zheng Y, Li X, Huang Y, Jia L, Li W. The circular RNA landscape of periodontal ligament stem cells during osteogenesis. J Periodontol. 2017;88:906–14.
Zheng Y, Li X, Huang Y, Jia L, Li W. Time series clustering of mRNA and lncRNA expression during osteogenic differentiation of periodontal ligament stem cells. PeerJ. 2018;6:e5214.
Kulterer B, Friedl G, Jandrositz A, Sanchez-Cabo F, Prokesch A, Paar C, et al. Gene expression profiling of human mesenchymal stem cells derived from bone marrow during expansion and osteoblast differentiation. BMC Genomics. 2007;8:70.
Brent MR. Past roadblocks and new opportunities in transcription factor network mapping. Trends Genet. 2016;32:736–50.
Buccitelli C, Selbach M. mRNAs, proteins and the emerging principles of gene expression control. Nat Rev Genet. 2020;21:630–44.
Wu L, Wei X, Ling J, Liu L, Liu S, Li M, et al. Early osteogenic differential protein profile detected by proteomic analysis in human periodontal ligament cells. J Periodontal Res. 2009;44:645–56.
Meleady P. Two-dimensional gel electrophoresis and 2D-DIGE. Methods Mol Biol. 1664:2018, 3–4.
Wang H, Li J, Zhang X, Ning T, Ma D, Ge Y, et al. Priming integrin alpha 5 promotes the osteogenic differentiation of human periodontal ligament stem cells due to cytoskeleton and cell cycle changes. J Proteome. 2018;179:122–30.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20:3710–5.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Krämer A, Green J, Pollard J, Tugendreich S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics. 2014;30:523–30.
Koide H, Holmbeck K, Lui JC, Guo XC, Driggers P, Chu T, et al. Mice deficient in AKAP13 (BRX) are osteoporotic and have impaired osteogenesis. J Bone Miner Res. 2015;30:1887–95.
Genetos DC, Wong A, Weber TJ, Karin NJ, Yellowley CE. Impaired osteoblast differentiation in annexin A2- and -A5-deficient cells. PLoS One. 2014;9:e107482.
Nistala H, Lee-Arteaga S, Smaldone S, Siciliano G, Carta L, Ono RN, et al. Fibrillin-1 and -2 differentially modulate endogenous TGF-β and BMP bioavailability during bone formation. J Cell Biol. 2010;190:1107–21.
Arfat Y, Basra MAR, Shahzad M, Majeed K, Mahmood N, Munir H. miR-208a-3p suppresses osteoblast differentiation and inhibits bone formation by targeting ACVR1. Mol Ther Nucleic Acids. 2018;11:323–36.
Demetriou M, Binkert C, Sukhu B, Tenenbaum HC, Dennis JW. Fetuin/alpha2-HS glycoprotein is a transforming growth factor-beta type II receptor mimic and cytokine antagonist. J Biol Chem. 1996;271:12755–61.
Kawaki H, Kubota S, Suzuki A, Suzuki M, Kohsaka K, Hoshi K, et al. Differential roles of CCN family proteins during osteoblast differentiation: involvement of Smad and MAPK signaling pathways. Bone. 2011;49:975–89.
Marini JC, Forlino A, Cabral WA, Barnes AM, San Antonio JD, Milgrom S, et al. Consortium for osteogenesis imperfecta mutations in the helical domain of type I collagen: regions rich in lethal mutations align with collagen binding sites for integrins and proteoglycans. Hum Mutat. 2007;28:209–21.
Maehata Y, Takamizawa S, Ozawa S, Izukuri K, Kato Y, Sato S, et al. Type III collagen is essential for growth acceleration of human osteoblastic cells by ascorbic acid 2-phosphate, a long-acting vitamin C derivative. Matrix Biol. 2007;26:371–81.
Wang C, Gu W, Sun B, Zhang Y, Ji Y, Xu X, et al. CTHRC1 promotes osteogenic differentiation of periodontal ligament stem cells by regulating TAZ. J Mol Histol. 2017;48:311–9.
Xu Q, Wu N, Cui L, Wu Z, Qiu G, Filamin B. The next hotspot in skeletal research? J Genet Genomics. 2017;44:335–42.
Fahmy-Garcia S, Farrell E, Witte-Bouma J, Robbesom-van den Berge I, Suarez M, Mumcuoglu D, et al. Follistatin effects in migration, vascularization, and osteogenesis in vitro and bone repair in vivo. Front Bioeng Biotechnol. 2019;7:38.
Bivi N, Condon KW, Allen MR, Farlow N, Passeri G, Brun LR, et al. Cell autonomous requirement of connexin 43 for osteocyte survival: consequences for endocortical resorption and periosteal bone formation. J Bone Miner Res. 2012;27:374–89.
Woods S, Barter MJ, Elliott HR, McGillivray CM, Birch MA, Clark IM, et al. miR-324-5p is up regulated in end-stage osteoarthritis and regulates Indian Hedgehog signalling by differing mechanisms in human and mouse. Clin Exp Pharmacol Physiol. 2019;77:87–100.
Gazzerro E, Smerdel-Ramoya A, Zanotti S, Stadmeyer L, Durant D, Economides AN, et al. Conditional deletion of gremlin causes a transient increase in bone formation and bone mass. J Biol Chem. 2007;282:31549–57.
Karsdal MA, Andersen TA, Bonewald L, Christiansen C. Matrix metalloproteinases (MMPs) safeguard osteoblasts from apoptosis during transdifferentiation into osteocytes: MT1-MMP maintains osteocyte viability. DNA Cell Biol. 2004;23:155–65.
Fratzl-Zelman N, Bächinger HP, Vranka JA, Roschger P, Klaushofer K, Rauch F. Bone matrix hypermineralization in prolyl-3 hydroxylase 1 deficient mice. Bone. 2016;85:15–22.
Gruenwald K, Castagnola P, Besio R, Dimori M, Chen Y, Akel NS, et al. Sc65 is a novel endoplasmic reticulum protein that regulates bone mass homeostasis. J Bone Miner Res. 2014;29:666–75.
Moenning A, Jäger R, Egert A, Kress W, Wardelmann E, Schorle H. Sustained platelet-derived growth factor receptor alpha signaling in osteoblasts results in craniosynostosis by overactivating the phospholipase C-gamma pathway. Mol Cell Biol. 2009;29:881–91.
Christiansen HE, Schwarze U, Pyott SM, AlSwaid A, Al Balwi M, Alrasheed S, et al. Homozygosity for a missense mutation in SERPINH1, which encodes the collagen chaperone protein HSP47, results in severe recessive osteogenesis imperfecta. Am J Hum Genet. 2010;86:389–98.
Guo W, Fan Z, Wang S, Du J. ALK5 is essential for tooth germ differentiation during tooth development. Biotech Histochem. 2019;94:481–90.
Fu R, Lv WC, Xu Y, Gong MY, Chen XJ, Jiang N, et al. Endothelial ZEB1 promotes angiogenesis-dependent bone formation and reverses osteoporosis. Nat Commun. 2020;11:460.
Case N, Xie Z, Sen B, Styner M, Zou M, O'Conor C, et al. Mechanical activation of β-catenin regulates phenotype in adult murine marrow-derived mesenchymal stem cells. J Orthop Res. 2010;28:1531–8.
Wu F, Jiao J, Liu F, Yang Y, Zhang S, Fang Z, et al. Hypermethylation of Frizzled1 is associated with Wnt/β-catenin signaling inactivation in mesenchymal stem cells of patients with steroid-associated osteonecrosis. Exp Mol Med. 2019;51:1–9.
Govoni KE, Baylink DJ, Mohan S. The multi-functional role of insulin-like growth factor binding proteins in bone. Pediatr Nephrol. 2005;20:261–8.
Vas V, Kovács T, Körmendi S, Bródy A, Kudlik G, Szeder B, et al. Significance of the Tks4 scaffold protein in bone tissue homeostasis. Sci Rep. 2019;9:5781.
Izu Y, Ezura Y, Mizoguchi F, Kawamata A, Nakamoto T, Nakashima K, et al. Type VI collagen deficiency induces osteopenia with distortion of osteoblastic cell morphology. Tissue Cell. 2012;44:1–6.
Chen C, Akiyama K, Wang D, Xu X, Li B, Moshaverinia A, et al. mTOR inhibition rescues osteopenia in mice with systemic sclerosis. J Exp Med. 2015;212:73–91.
Drabek K, van de Peppel J, Eijken M, van Leeuwen JP. GPM6B regulates osteoblast function and induction of mineralization by controlling cytoskeleton and matrix vesicle release. J Bone Miner Res. 2011;26:2045–51.
Yuan H, Zhao H, Wang J, Zhang H, Hong L, Li H, et al. MicroRNA let-7c-5p promotes osteogenic differentiation of dental pulp stem cells by inhibiting lipopolysaccharide-induced inflammation via HMGA2/PI3K/Akt signal blockade. 2019;46:389–97.
Chen Z, Zhang K, Qiu W, Luo Y, Pan Y, Li J, et al. Genome-wide identification of long noncoding RNAs and their competing endogenous RNA networks involved in the odontogenic differentiation of human dental pulp stem cells. Stem Cell Res Ther. 2020;11:114.
Neves H, Kwok HF. In sickness and in health: the many roles of the minichromosome maintenance proteins. Biochim Biophys Acta Rev Cancer. 1868;2017:295–308.
Müller L, Gerighausen D, Farman M, Zeckzer D. Sierra platinum: a fast and robust peak-caller for replicated ChIP-seq experiments with visual quality-control and -steering. BMC Bioinformatics. 2016;17:377.
Atashi F, Modarressi A, Pepper MS. The role of reactive oxygen species in mesenchymal stem cell adipogenic and osteogenic differentiation: a review. Stem Cells Dev. 2015;24:1150–63.
Finkel T. Signal transduction by reactive oxygen species. J Cell Biol. 2011;194:7–15.
Rhee SG. Cell signaling. H2O2, a necessary evil for cell signaling. Science. 2006;312:1882–3.
Long F. Building strong bones: molecular regulation of the osteoblast lineage. Nat Rev Mol Cell Biol. 2011;13:27–38.
Salehi H, Collart-Dutilleul P-Y, Gergely C, Cuisinier FJG. Confocal Raman microscopy to monitor extracellular matrix during dental pulp stem cells differentiation. J Biomed Opt. 2015;20:076013.
Zhu JX, Sasano Y, Takahashi I, Mizoguchi I, Kagayama M. Temporal and spatial gene expression of major bone extracellular matrix molecules during embryonic mandibular osteogenesis in rats. Histochem J. 2001;33:25–35.
Kadler KE, Baldock C, Bella J, Boot-Handford RP. Collagens at a glance. J Cell Sci. 2007;120:1955.
Hong H-H, Pischon N, Santana RB, Palamakumbura AH, Chase HB, Gantz D, et al. A role for lysyl oxidase regulation in the control of normal collagen deposition in differentiating osteoblast cultures. J Cell Physiol. 2004;200:53–62.
Langenbach F, Handschel J. Effects of dexamethasone, ascorbic acid and β-glycerophosphate on the osteogenic differentiation of stem cells in vitro. Stem Cell Res Ther. 2013;4:117.
Lu Y, Kamel-El Sayed SA, Wang K, Tiede-Lewis LM, Grillo MA, Veno PA, et al. Live imaging of type I collagen assembly dynamics in osteoblasts stably expressing GFP and mCherry-tagged collagen constructs. J Bone Miner Res. 2018;33:1166–82.
Dong R, Du J, Wang L, Wang J, Ding G, Wang S, et al. Comparison of long noncoding RNA and mRNA expression profiles in mesenchymal stem cells derived from human periodontal ligament and bone marrow. Biomed Res Int. 2014;2014:317853.
Hakki SS, Kayis SA, Hakki EE, Bozkurt SB, Duruksu G, Unal ZS, et al. Comparison of mesenchymal stem cells isolated from pulp and periodontal ligament. J Periodontol. 2015;86:283–91.
Yang H, Gao L-N, An Y, Hu C-H, Jin F, Zhou J, et al. Comparison of mesenchymal stem cells derived from gingival tissue and periodontal ligament in different incubation conditions. Biomaterials. 2013;34:7033–47.
Harumiya S, Gibson MA, Koshihara Y. Antisense suppression of collagen VI synthesis results in reduced expression of collagen I in normal human osteoblast-like cells. Biosci Biotechnol Biochem. 2002;66:2743–7.
Ganburged G, Suda N, Saito M, Yamazaki Y, Isokawa K, Moriyama K. Dilated capillaries, disorganized collagen fibers and differential gene expression in periodontal ligaments of hypomorphic fibrillin-1 mice. Cell Tissue Res. 2010;341:381–95.
Choi Y-A, Lim J, Kim KM, Acharya B, Cho J-Y, Bae Y-C, et al. Secretome analysis of human BMSCs and identification of SMOC1 as an important ECM protein in osteoblast differentiation. J Proteome Res. 2010;9:2946–56.
Forni MF, Peloggia J, Trudeau K, Shirihai O, Kowaltowski AJ. Murine mesenchymal stem cell commitment to differentiation is regulated by mitochondrial dynamics. Stem Cells. 2016;34:743–55.
Chen C-T, Shih Y-RV, Kuo TK, Lee OK, Wei Y-H. Coordinated changes of mitochondrial biogenesis and antioxidant enzymes during osteogenic differentiation of human mesenchymal stem cells. Stem Cells. 2008;26:960–8.
Wanrooij S, Falkenberg M. The human mitochondrial replication fork in health and disease. Biochim Biophys Acta. 2010;1797:1378–88.
Senior AE, Nadanaciva S, Weber J. The molecular mechanism of ATP synthesis by F1F0-ATP synthase. Biochim Biophys Acta. 2002;1553:188–211.
Ernster L, Schatz G. Mitochondria: a historical review. J Cell Biol. 1981;91:227s–55s.
Li Q, Gao Z, Chen Y, Guan M-X. The role of mitochondria in osteogenic, adipogenic and chondrogenic differentiation of mesenchymal stem cells. Protein Cell. 2017;8:439–45.
Geng W, Shi H, Zhang X, Tan W, Cao Y, Mei R. Substance P enhances BMSC osteogenic differentiation via autophagic activation. Mol Med Report. 2019;20:664–70.
Anilkumar N, Sirker A, Shah AM. Redox sensitive signaling pathways in cardiac remodeling, hypertrophy and failure. Front Biosci (Landmark Ed). 2009;14:3168–87.
Pérez S, Taléns-Visconti R, Rius-Pérez S, Finamor I, Sastre J. Redox signaling in the gastrointestinal tract. Free Radic Biol Med. 2017;104:75–103.
Bradley EW, Carpio LR, van Wijnen AJ, McGee-Lawrence ME, Westendorf JJ. Histone deacetylases in bone development and skeletal disorders. Physiol Rev. 2015;95:1359–81.
Abed É, Delalandre A, Lajeunesse D. Beneficial effect of resveratrol on phenotypic features and activity of osteoarthritic osteoblasts. Arthritis Res Ther. 2017;19:151.
Huang R-X, Tao J. Nicotinamide mononucleotide attenuates glucocorticoid-induced osteogenic inhibition by regulating the SIRT1/PGC-1α signaling pathway. Mol Med Report. 2020;22:145–54.
Zhou W, Liu Y, Shen J, Yu B, Bai J, Lin J, et al. Melatonin increases bone mass around the prostheses of OVX rats by ameliorating mitochondrial oxidative stress via the SIRT3/SOD2 signaling pathway. Oxidative Med Cell Longev. 2019;2019:4019619.
Sun H, Wu Y, Fu D, Liu Y, Huang C. SIRT6 regulates osteogenic differentiation of rat bone marrow mesenchymal stem cells partially via suppressing the nuclear factor-κB signaling pathway. Stem Cells. 2014;32:1943–55.
Chen EEM, Zhang W, Ye CCY, Gao X, Jiang LLJ, Zhao TTF, et al. Knockdown of SIRT7 enhances the osteogenic differentiation of human bone marrow mesenchymal stem cells partly via activation of the Wnt/beta-catenin signaling pathway. Cell Death Dis. 2017;8:e3042.
Sun W, Qiao W, Zhou B, Hu Z, Yan Q, Wu J, et al. Overexpression of Sirt1 in mesenchymal stem cells protects against bone loss in mice by FOXO3a deacetylation and oxidative stress inhibition. Metabolism. 2018;88:61–71.
Chang S, Xia B, Wu Y, Peng T, Jiang Q, Li J, et al. Forkhead box O3 protects the heart against paraquat-induced aging-associated phenotypes by upregulating the expression of antioxidant enzymes. Aging Cell. 2019;18:e12990.
Shi D, Li X, Chen H, Che N, Zhou S, Lu Z, et al. High level of reactive oxygen species impaired mesenchymal stem cell migration via overpolymerization of F-actin cytoskeleton in systemic lupus erythematosus. Pathol Biol (Paris). 2014;62:382–90.
Wilson C, González-Billault C. Regulation of cytoskeletal dynamics by redox signaling and oxidative stress: implications for neuronal development and trafficking. Front Cell Neurosci. 2015;9:381.
Zhao X, Petursson F, Viollet B, Lotz M, Terkeltaub R, Liu-Bryan R. Peroxisome proliferator-activated receptor γ coactivator 1α and FoxO3A mediate chondroprotection by AMP-activated protein kinase. Arthritis Rheumatol. 2014;66:3073–82.
Xi G, Rosen CJ, Clemmons DR. IGF-I and IGFBP-2 stimulate AMPK activation and autophagy, which are required for osteoblast differentiation. Endocrinology. 2016;157:268–81.
Shares BH, Busch M, White N, Shum L, Eliseev RA. Active mitochondria support osteogenic differentiation by stimulating β-catenin acetylation. J Biol Chem. 2018;293:16019–27.
Golan K, Vagima Y, Ludin A, Itkin T, Cohen-Gur S, Kalinkovich A, et al. S1P promotes murine progenitor cell egress and mobilization via S1P1-mediated ROS signaling and SDF-1 release. Blood. 2012;119:2478–88.
Khalid S, Yamazaki H, Socorro M, Monier D, Beniash E, Napierala D. Reactive oxygen species (ROS) generation as an underlying mechanism of inorganic phosphate (P(i))-induced mineralization of osteogenic cells. Free Radic Biol Med. 2020;153:103–11.
Jaiswal RK, Jaiswal N, Bruder SP, Mbalaviele G, Marshak DR, Pittenger MF. Adult human mesenchymal stem cell differentiation to the osteogenic or adipogenic lineage is regulated by mitogen-activated protein kinase. J Biol Chem. 2000;275:9645–52.
Phinney DG. Functional heterogeneity of mesenchymal stem cells: implications for cell therapy. J Cell Biochem. 2012;113:2806–12.
Mo M, Wang S, Zhou Y, Li H, Wu Y. Mesenchymal stem cell subpopulations: phenotype, property and therapeutic potential. Cell Mol Life Sci. 2016;73:3311–21.
We thank Beijing igeneCode Biotech Co., Ltd. for their assistance in bioinformatic analysis. We also appreciated PhD. Dongqiang Zeng from the department of Oncology of Nanfang Hospital of Southern Medical University for the illustration of Fig. 4.
This research was supported by the grants from the National Natural Science Foundation of China (Grant No. 81870755, 81700943), Natural Science Foundation of Guangdong Province (Grant No. 2020A1515011455), China Postdoctoral Science Foundation (Grant No. 2019 M663009, 2018 M640804), Medical Scientific Research Foundation of Guangdong Province of China (Grant No. A2020090, A2017525) and President Foundation of Nanfang Hospital, Southern Medical University (Grant No. 2019B002).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Nanfang Hospital, Southern Medical University. Written informed consent was obtained from all adult patients.
Consent for publication
The authors have declared no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Boxplots of normalized protein ratio of all 12 samples.
Map of the “Sirtuin signalling pathway” canonical pathway in the Ingenuity Pathway Analysis database. Proteins labelled in red, yellow, green, blue and pink are proteins from clusters 1, 2, 3, 4 and 5 of differentially expressed proteins, respectively. Grey indicates proteins that are not differentially expressed.
Map of the “Actin Cytoskeleton signalling pathway” canonical pathway in the Ingenuity Pathway Analysis database. Proteins labelled in red, yellow, green, blue and pink are proteins from clusters 1, 2, 3, 4 and 5 of differentially expressed proteins, respectively. Grey indicates proteins that are not differentially expressed.
Map of the “AMPK signalling pathway” canonical pathway in the Ingenuity Pathway Analysis database. Proteins labelled in red, yellow, green, blue and pink are proteins from clusters 1, 2, 3, 4 and 5 of differentially expressed proteins, respectively. Color grey referrer to the proteins not differentially expressed.
Map of “ERK/MAPK signalling pathway” canonical pathway in the Ingenuity Pathway Analysis database. Proteins labelled in red, yellow, green, blue and pink are proteins from clusters 1, 2, 3, 4 and 5 of differentially expressed proteins, respectively. Grey indicates proteins that are not differentially expressed.
Map of the “Sphingosine-1-phosphate signalling pathway” canonical pathway in the Ingenuity Pathway Analysis database. Proteins labelled in red, yellow, green, blue and pink are proteins from clusters 1, 2, 3, 4 and 5 of differentially expressed proteins, respectively. Grey indicates proteins that are not differentially expressed.
Map of the “Wnt/β-catenin signalling pathway” canonical pathway in the Ingenuity Pathway Analysis database. Proteins labelled in red, yellow, green, blue and pink are proteins from clusters 1, 2, 3, 4 and 5 of differentially expressed proteins, respectively. Grey indicates proteins that are not differentially expressed.
Additional file 8: Figure S8. The protein-protein interaction network of oxidative phosphorylation and redox-sensitive signalling pathways. The interaction between oxidative phosphorylation and redox-sensitive signalling pathways (sirtuin signalling pathway, AMPK signalling, actin cytoskeleton signalling, Sphingosine-1-phosphate signalling, ERK/MAPK signalling and Wnt/β-catenin signalling pathways). The sirtuin signalling pathway was the most associated signalling pathway.
Sample set of TMT proteomic analysis. Three six-plex TMT sets were used for the four time point analysis (D0, D3, D7 and D14) with triplicate biological replication. A sample combined with all 12 sample was labeled with Label 130 and used as a reference sample for normalization between different sets. D0, 3, 7, 14 means samples undifferented, differentiated for 3, 7 and 14 days.
Details of differentially expressed proteins.
Details of GO enrichments of each cluster of differentially expressed proteins.
Details of KEGG pathway enrichments of differentially expressed proteins.
Details of activation prediction of canonical pathways from Ingenuity Pathway Analysis.
Supplemental materials and methods.
About this article
Cite this article
Li, J., Wang, Z., Huang, X. et al. Dynamic proteomic profiling of human periodontal ligament stem cells during osteogenic differentiation. Stem Cell Res Ther 12, 98 (2021). https://doi.org/10.1186/s13287-020-02123-6
- Human periodontal ligament stem cell
- Osteogenic differentiation
- Dynamic proteomics
- Oxidative phosphorylation