Introduction
Mechanical stress enhances bone metabolism and periodontal tissue remodeling.1 During orthodontic tooth movement (OTM), bone remodeling is initiated via the periodontal ligament.2 As the main mesenchymal stem cells (MSCs) in the periodontal ligament, periodontal ligament stem cells (PDLSCs) play an important role in mechanical signal transduction. Currently, a consensus has been reached that cyclic mechanical tension is a strong driver of the differentiation of PDLSCs into the osteoblast lineage.24 Mechanical tension activates calcium channels,5 which activate the ERK1/2 and P38 MAPK pathways through integrin-FAK or protein kinase (PKC)-SR signaling6 and induce the phosphorylation of Runt-associated transcription factor 2 (Runx2),7 promoting osteogenic precursor cell synthesis and the transcription of mineralizable proteins.8 At present, the TGF-, BMP, MAPK, Notch, Wnt, Hedgehog, FGF, and Hippo signaling pathways have been found to be involved in this process.
Force parameters (including magnitude, frequency, and duration) are crucial for well-regulated tissue remodeling. However, numerous in vitro studies performed to date show enormous heterogeneity in tensile force parameters.9 In different studies, cyclic tension was applied with magnitudes ranging from 1% to 24%, frequencies ranging from 0.1 Hz to 1.0 Hz, and stimuli duration ranging from 1 hour to 6 days, thereby reducing comparability between different studies.9 To establish strategies to optimize tensile force parameters, it is of particular importance to understand how different tensile force parameters affect the osteogenic differentiation of PDLSCs.
The effects of different tensile force magnitudes and durations have been investigated in some studies. Among the magnitudes, a magnitude of 10% generally led to a lower level of inflammation and a higher level of osteogenesis,10 whereas a magnitude of 12% was found to correlate well with strain conditions at the mid-root under physiological loading conditions11,12 and to induce optimal effects in both the proliferation and osteogenesis of PDLSCs.13 Cyclic tension alone at 3000 strain significantly enhanced SATB Homeobox 2 (Satb2) after 3 h of loading and significantly upregulated Runx2 after 6 h.14 The synthesis of BMP9 increased under 6-h continuously applied cyclic tension.15 In addition, 12% cyclic tensile force gradually upregulated the expression of Runx2, alkaline phosphatase (ALP), and osteocalcin (OCN) with force durations of 6 h, 12 h, and 24 h, respectively.16,17 The protein level of osterix increased stepwise following 3 h, 6 h, 12 h, and 24 h of exposure to tensile strain.14 Recently, temporal gene expression patterns were delineated.17
Tensile frequency varies largely among different studies. The ROCKTAZ pathway and its interaction with Cbf1 were found to be essential for the cyclic tension (12% elongation, 0.1 Hz)-induced osteogenic differentiation in PDLSCs.18 Cyclic tension (10% elongation, 0.5 Hz) stimulated the osteogenic differentiation of PDLSCs by inhibiting miR-129-5p expression and activating the BMP2/Smad pathway.17 LncRNAs-miRNAs-mRNAs networks in PDLSCs were depicted under cyclic tension (10% elongation, 1.0 Hz).19 However, there have been rare studies examining the impact of different cyclic tensile frequencies on osteogenesis of PDLSCs and the expression of relevant genes thus far. The low-magnitude high-frequency (LMHF) vibration approach was excluded because it is used to simulate a masticatory force, while cyclic tension is used to simulate an orthodontic force, and the two methods of force application are completely different.2,20 Previous animal studies on long bone distraction osteogenesis have shown that loading frequency affects the osteogenic response of bone tissue.21 The mechano-regulation of trabecular bone adaptation is logarithmically dependent on the loading frequency.22 Therefore, we hypothesized that tensile frequency would affect osteogenesis of PDLSCs, in which some tensile frequency-sensitive genes may play an important role. To test our hypotheses, human PDLSCs were subjected to cyclic mechanical tension at different frequencies of 0.10.7 Hz to examine the osteoblastic differentiation of PDLSCs, and high-throughput sequencing was performed to characterize the frequency-course expression patterns of mRNA during the osteogenic differentiation of PDLSCs. This study aimed to investigate the effects of tensile frequency on the osteogenic differentiation of PDLSCs as well as the relevant molecular mechanisms.
Healthy periodontal ligament tissues were scraped from the middle third of tooth roots, which were extracted for orthodontic reasons, with informed consents. All donors were aged from 14 to 16 years and had no systemic or oral diseases. The periodontal ligament tissues were cut into small pieces and enzymatically digested for 40 min at 37C with collagenase I (3 mg/mL, Sigma-Aldrich, St. Louis, MO, USA) and dispase II (4 mg/mL; Sigma-Aldrich). The cells were seeded in 25 cm2 flasks (Falcon, BD Biosciences, Franklin Lakes, NJ, USA) with -minimal essential medium (-MEM; Sigma-Aldrich) supplemented with 10% fetal bovine serum (FBS, Gibco, Life Technologies Co., Grand Island, NY, USA) and antibiotics (100 U/mL penicillin and 100 g/mL streptomycin, Hyclone, Logan, Utah, USA), and incubated in a humidified atmosphere (37C, 5% CO2). The medium was changed every 3 days. After reaching 80% of confluence, the cells were detached with 0.25% trypsin/EDTA (Gibco, Life Technologies Co., Grand Island, NY, USA), and single-cell suspensions were cloned with the limiting-dilution method to purify the stem cells.17 Cell clusters from the colony were trypsinized and serially sub-cultured.
The third-passage PDLSCs were sub-cultured into six-well plates until confluent. The culture medium was then removed, and cells were fixed with 4% formaldehyde (Zhonghuihecai, Xian, CN) for 20 min, permeabilized with 0.3% Triton X-100 (Zhonghuihecai) for 5 min, and incubated with primary antibodies (anti-pan-cytokeratin, 1:300, Abcam, Cambridge, MA, USA; anti-vimentin, 1:500, Abcam; anti-STRO-1, 1:200, Abcam; anti-CD146, 1:200, Abcam) overnight at 4C. The cells were then washed with PBS and incubated with CY3/FITC-conjugated secondary antibodies (1:500, Zhuangzhi, Xian, CN) in darkness for 30 min, and then washed with PBS. Finally, the nuclei were counterstained with 4, 6-diamidino-2-phenylindole (DAPI, Zhuangzhi), and fluorescent images were captured with a fluorescence microscope (Olympus, Japan).
For osteogenic and adipogenic differentiation, PDLSCs were seeded into six-well plates at a density of 2105 cells/well, and after reaching 80% confluence, the medium was replaced with an osteogenic or adipogenic inductive medium (Osteogenesis or Adipogenesis Differentiation Kit, Cyagen, USA). Seven days after osteogenic induction, the cells were stained with ALP (ALP staining kit, Solarbio, CN), and 14 days after osteogenic incubation, the cells were stained with Alizarin Red (ARS staining kit, Cyagen, USA). After 21 days of adipogenic incubation, the cells were stained with Oil Red O (Cyagen, USA).
Flexcell FX-5000T Tension Plus System (Flexcell International Corporation, Hillsborough, NC, USA) was used to mimic the tensile force exerted on PDLSCs during OTM, according to previous studies.23 Cyclic tensile loading experiments were performed on the fourth-passage PDLSCs from four different healthy donors in triplicate. PDLSCs were seeded onto six-well type I collagen (COL-I)-coated silicone culture plates (Flexcell International Corporation) at a density of 2105 cells/well. Upon reaching 80% confluence, the cells were serum-starved overnight, and the medium was changed to osteogenic medium (Cyagen, USA). Cyclic tensile force (12% bottom membrane elongation) was applied to different plates at different frequencies of 0.1 Hz, 0.5 Hz, and 0.7 Hz. Control cells were cultured under identical culture condition but without mechanical stimulation.
After 10 h of cyclic tensile force, the cells were collected. Proteins were isolated, electrophoretically separated, and immunoblotted as previously described.23 Briefly, PDLSCs were lysed with RIPA buffer containing 1% phenylmethanesulfonyl fluoride (PMSF, proteinase inhibitor, Zhonghuihecai) and 1% phosphatase inhibitor (Zhonghuihecai). After centrifugation, the supernatant was collected and measured quantitatively using a BCA Protein Assay Kit (Absin, Shanghai, CN). Total protein from cell lysates (20 g/lane) was separated by SDS-PAGE gels (Beyotime, Hangzhou, CN) and then transferred onto a polyvinylidene difluoride (PVDF) membrane (EMD Millipore, Billerica, MA, USA). After blocking with 5% skimmed milk in tris-buffered saline tween-20 (TBST) for 2 h at room temperature, the membranes were incubated overnight at 4C with primary antibodies (runt-related transcription factor 2 (Runx2), 1:500, ImmunoWay, USA; COL-I, 1:1000, Proteintech, USA; Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 1:10000, Proteintech, USA), followed by incubation with horseradish peroxidase (HRP)-conjugated secondary antibody (1:2000, Proteintech, USA) for 2 h at room temperature. The protein expression was visualized using ChemiDocTM XRS+ (Bio-Rad Laboratories, Inc., Hercules, CA, USA) with an enhanced chemiluminescence (ECL) kit (Millipore, Billerica, MA, USA). GAPDH was used as an internal control for normalization.
The fourth-passage PDLSCs from three donors were used for RNA sequencing after tension loading. Total RNA was extracted from the four groups of cells (normal PDLSCs and PDLSCs tensioned at frequencies of 0.1 Hz, 0.5 Hz, and 0.7 Hz for 6 h) using the Trizol (Sigma-Aldrich), according to the manufacturers protocols. After digestion with DNase, rRNA were depleted using a Ribo-Zero magnetic kit, and sequencing libraries were constructed as previously described.24 The sequencing of the cDNA library was carried out by Gene Denovo Biotechnology Co. (Guangzhou, China). The gene expression level was evaluated by reads per kilobase transcriptome per million mapped reads (RPKM). Requirements for filtering differentially expressed genes (DEGs) were as follows: (1) |log2 (fold-change)| 1; (2) p value < 0.05. DESeq2 (differential gene expression analysis based on the negative binomial distribution)25 was used to calculate p values and adjusted p (adj. p) values. Heatmaps and volcano plots analyses were used to visualize these DEGs using the Complex Heatmap package and ggplots2 package of R software. The online tool Venny 2.1 (http://bioinfogp.cnb.csic.es/tools/venny/index.html) was applied to identify the common DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used for annotation visualization and integrated discovery. Raw data of the performed RNA sequencing (RNA-seq) were recorded in the SRA database with the SRA accession: PRJNA665587.
Gene expression pattern analysis is used to cluster genes of similar expression patterns for multiple samples in a tensile frequency order. To examine the expression pattern of DEGs in different frequencies, the expression data of each sample (in the order of treatment) were normalized to 0, log2 (v1/v0), and log2 (v2/v0), and then clustered using the Short Time-series Expression Miner software (STEM).26 The parameters were set as follows: (1) maximum unit change in model profiles between frequency points was 1; (2) maximum output profiles number was 20 (similar profiles will be merged; and (3) minimum ratio of fold change of DEGs was no less than 2.0. The clustered profiles with p value < 0.05 were considered as significant profiles. Then, the DEGs in all profiles or in each profile were subjected to GO and KEGG pathway enrichment analysis. STRING (https://string-db.org/) was used for PPI network analysis. The DEGs were mapped onto the PPI network with a minimum interaction score of 0.4. Cytoscape v3.7.1 software was used to visualize the PPI network. Gene network clustering analysis was performed to identify the key PPI network modules, using the MCODE and cytoHuba app from the Cytoscape software suite. Adj. p value < 0.05 was set as the significance threshold.
Total RNA was extracted using Trizol (Sigma-Aldrich) according to the manufacturers protocols. Quantitative real-time PCR (RTqPCR) was performed in triplicate using the Power SYBR Green PCR Mastermix (Applied Biosystems, Foster City, CA, USA). Sequences of the primers used are shown in Table S1. The mean expression values were calculated relative to GAPDH, which was used as an internal control for normalization.
For differential gene expression analysis, a likelihood ratio test was used assuming an underlying zero-inflated negative binomial distribution.25,27 FDR-corrected p-values were calculated using the Benjamini-Hochberg procedure.28 K-means cluster method29 was used for frequency cluster expression pattern analysis by the STEM software. P and Q values of GO/KEGG functional analysis and frequency cluster analysis were calculated based on a hypergeometric distribution.30 Values of the relative protein and mRNA expression were expressed as mean SD within each group. One-way ANOVA followed by SNK post hoc tests was used in the Western blotting assay and RTqPCR assay. The significance threshold was set at 0.05.
Immunofluorescent staining revealed that the isolated cells were positive for CD146, vimentin, and STRO-1 (Figure 1AC, respectively) but negative for pan-cytokeratin (Figure 1D), confirming that the cells were mesenchymal stem cells of mesodermal origin. The osteogenic potential of PDLSCs was determined by positive ALP staining after 7 days of osteogenic induction (Figure 1E and F) and red mineralized matrix nodules in Alizarin Red staining after 14 days of osteogenic induction (Figure 1G and H). The presence of red lipid droplets in Oil Red staining after 21 days of adipogenic induction (Figure 1I) indicated the adipogenic differentiation of PDLSCs.
Figure 1 Identification of PDLSCs. The Immunofluorescence showed that the cells were positive for CD146 (A, red), vimentin (B, green), and STRO-1 (C, red) but negative for pan-cytokeratin (D). Scale bar = 50 m. ALP staining was positive both visually (E, black) and under the microscope (F, black, Scale bar = 100 m), after 7 days of osteogenesis induction culture. After 14 days of osteogenesis induction, ARS staining was observed to be positive by the naked eye (G, red), and mineralized nodules were obvious under the microscope (H, red, Scale bar = 100 m). After 21 days of adipogenic induction, oil red O-positive lipid clusters were observed microscopically (I, red, Scale bar = 100 m).
PDLSCs were exposed to 12% cyclic tension, which contributed to the cellular reorientation, including an irregular arrangement at the central region and a parallel arrangement at the peripheral region of the plate (Figure 2A). As shown in Figure 2B and C, after consecutive 10-h cyclic tensile force loading, the protein expression levels of Runx2 and COLI increased with increasing tensile frequency from 0.1 Hz to 0.7 Hz and were remarkably higher than that in the group without tension application (p < 0.05). The result suggested that mechanical tension (12% deformation) upregulated the osteogenesis of PDLSCs in a frequency-dependent manner. Higher frequencies of cyclic tension were associated with higher osteogenic differentiation of PDLSCs.
Figure 2 Cyclic tension promoted the osteogenic differentiation of PDLSCs. Under cyclic equibiaxial tension, PDLSCs reoriented in parallel alignment at the peripheral region of the plate, while in random orientation at the center of the plate (A, Scale bar = 100 m). Western blotting was used to detect protein levels of osteogenesis-related genes, COL-I and Runx2, at different tensile frequencies (B and C). **p < 0.01, ***p < 0.001, vs control group.
It has been previously shown that cyclic mechanical tensile stress can improve osteogenesis of PDLSCs, and that consecutive 46 h of tension can significantly upregulate the mRNA expression of osteogenesis-related genes.16,17 To gain insight into the molecular mechanisms by which mechanical tension stimulates osteoblast differentiation of PDLSCs, total RNA was extracted from PDLSCs to conduct RNA-seq after 6 h of cyclic tension. The mRNA expression profiles of the PDLSCs at different tensile frequencies (0.1 Hz, 0.5 Hz, and 0.7 Hz) were detected. Comparative expression analyses were performed according to the different frequencies of the tensile stress (0.1 Hz vs control, 0.5 Hz vs control, and 0.7 Hz vs control). Heat maps of the top 40 DEGs (Figure 3AC) and volcano plots (Figure 3DF) were depicted. In total, 50 mRNAs were upregulated, and 261 mRNAs were downregulated at 0.1 Hz. At 0.5 Hz, 656 mRNAs were upregulated, and 1474 mRNAs were downregulated. At 0.7Hz, 139 mRNAs were upregulated, and 194 mRNAs were downregulated. A Venny analysis (Figure 3G and H) showed that 78 genes were simultaneously upregulated and 118 were simultaneously downregulated among the 0.1 Hz, 0.5 Hz, and 0.7 Hz groups. The GO analysis (Figure 4A) demonstrated that changes in biological processes (BPs) were mainly enriched in metabolic process, response to stimulus, biological regulation, signaling, and localization. Changes in Cellular Components (CCs) were mainly enriched in organelle, membrane, macromolecular complex, and membrane-enclosed lumen. Moreover, binding, catalytic activity, and nucleic acid binding transcription factor activity emerged as the highest-ranked Molecular Function (MF) groups. As shown in Figure 4B, DNA replication, cell cycle, and the TNF signaling pathway were significantly enriched in the KEGG pathway. Within the primary category Environmental Information Processing, Signal transduction, and Signaling molecules and interaction were strongly enriched (Figure 4C).
Figure 3 Identification of DEGs among different frequencies. (AC) Heatmaps of the top 40 DEGs between 0.1 Hz/0.5 Hz/0.7 Hz and static culture, respectively. Red rectangles represent high expression, and blue rectangles represent low expression. (DF) Volcano plot of DEGs between 0.1 Hz/0.5 Hz/0.7 Hz and static culture, respectively. The red plots represent upregulated genes, the blue plots represent downregulated genes, and the black plots represent nonsignificant genes. (G) Venn diagram of upregulated DEGs among 0.1 Hz, 0.5 Hz, and 0.7 Hz. (H) Venn diagram of downregulated DEGs among all the three frequencies.
Figure 4 Functional enrichment analysis of all DEGs. (A) GO enrichment analysis of all DEGs among different frequencies. (B) Top 20 pathways of the KEGG enrichment analysis of all DEGs among different frequencies, with the KEGG pathway annotation (C). The screening criteria for significance were p value < 0.05.
The sequencing data were normalized to the control, and trend analyses of DEGs were identified using STEM. In Figure 5A, within the 20 model profiles, eight mRNA trend profiles were statistically significant. The profile number assigned by STEM was on the top left corner of each profile box, p value was on the bottom left, and the number of the cardinality of each cluster was on the top right corner. As shown in Table 1, among different profiles, the top-ranked KEGG pathways were mainly in the metabolic pathways, PI3K-Akt signaling pathway, cytokine-cytokine receptor interaction, and MAPK signaling pathway. A continuous downregulation pattern was found in profile 0 (Figure 5B), in which the high-ranked BPs, CCs, and MFs in GO enrichment (Figure 5C) were similar to those in Figure 4A, and inflammatory pathways such as arachidonic acid metabolism, peroxisome, and cytokinecytokine receptor interaction were strongly enriched (Figure 5D).
Table 1 Top 10 of KEGG Enrichment Among Different Profiles
Figure 5 Frequency series clustering analysis on expression profiles of mRNAs by STEM. (A) Within the 20 model profiles, eight mRNA trend profiles were statistically significant. The number at the upper-left corner of each profile box was the profile number assigned by STEM, the number on the bottom left was the p value, and the number on the top-right corner was the number of genes within each cluster. (B) Persistently downregulated genes along frequency were clustered in profile 0. (C) GO enrichment of profile 0. (D) Top 10 pathways of the KEGG enrichment of profile 0.
The interactions of 194 DEGs in profile 0 were analyzed using the STRING online database, and the PPI network was obtained using the Cytoscape software (Figure 6A). The MCODE plugin was then used to investigate the key PPI network modules, and one key module with four genes (EYA1, SIX5, SALL1, FRAS1) was identified (Figure 6B). The cytoHubba plugin was then used to analyze hub genes with maximum correlation criterion (MCC)/Degree, and genes with the top 10 scores were respectively identified. The intersection (EYA1, SALL1) of hub genes according to the above three methods were selected for further RTqPCR validation (Figure 6C). The results of RTqPCR (Figure 6D and E) showed that the mRNA expression of EYA1 and SALL1 decreased with increasing frequency from 0.1 Hz to 0.7 Hz, which were highly consistent with our high-throughput sequencing.
Figure 6 Identification and validation of mechanofreqency-sensitive hub genes. (A) The interaction network between proteins coded by the DEGs in profile 0. The nodes represent genes, and the edges represent links between genes. Blue represents downregulated genes. (B) The highest scoring module was extracted by MCODE. (C) The intersection was obtained among modules measured by MCODE, and the top 10 highly connected genes were identified using MCC and Degree in cytoHubba. (D and E) Validation of the expression of the two intersection genes, EYA1 and SALL1, using RTqPCR. *p < 0.05, vs control group.
OTM is based on remodeling processes in the periodontal ligament and the alveolar bone. PDLSCs play an important role in mechano-transduction and in promoting periodontal tissue regeneration in OTM.31,32 It is evident that cyclic tensile force regulates the osteogenic differentiation of PDLSCs.33 A complex network of signaling molecules regulates the osteoblastic differentiation of PDLSCs under cyclic tension.5,17,34 The heterogeneity of mechanical force parameters (duration, magnitude, frequency, and others)9 led to the heterogeneity of the osteogenic phenotype and gene regulation in different studies. This study focused on the effect of tensile frequency on the osteogenic differentiation of PDLSCs and attempted to elucidate a potential mechanism.
In the current study, we successfully isolated and characterized human PDSLCs. The cells at the 4th to 6th passage were used, whose phenotype was generally believed to be maintained.35 Mechanical tension was applied using the Flexcell tension system, which has been widely used in PDLSCs studies.2,20 We observed that the cellar reorientation after force loading was similar to that of the previous study,23 and this could be attributed to the mechano-responsive stress fibersfocal adhesion system.36 The present study showed that cyclic mechanical tension (magnitude: 12% deformation, duration: 10 h) in the range of 0.10.7 Hz promoted osteogenic marker genes including Runx2 and COL-I in PDLSCs, and their protein expressions increased with the increasing tensile frequency. Runx2 is an osteogenic lineage commitment specific transcription factor, which binds to the specific cis-acting elements of osteoblasts to promote the transcription and translation of OCN, osteopontin (OPN), bone sialoprotein (BSP), and COL-I.8 COL-I, which acts as a template onto which minerals are deposited to form bone matrix,37 is the major constituent of extracellular matrix in the periodontal ligament and bone, and is confirmed to be essential for osteogenesis in response to tension during OTM.38 Consistently, upregulation of RUNX2 and COL-I in response to tension was reported in most studies.9 To the best of our knowledge, the present study is the first to reveal the frequency dependence during cyclic tension in enhancing the expression of osteogenic markers within the first 10 h of cyclic tension application, which may inform a new method of accelerating OTM.
Using RNA-seq, we observed that mRNAs in strained PDLSCs were mainly enriched in response to the stimulus process, signal transduction, and relative pathways such as mismatch repair, TNF signaling pathway, and FOXO signaling pathway, which were associated with cell survival and differentiation as well as immune and inflammatory responses.39,40 The STEM platform was also used to investigate how gene expression profiles change with tensile frequency during the osteoblast differentiation of PDLSCs under cyclic tension. Eight trend profiles were noted as significant. Genes in these profiles were mainly enriched in the metabolic pathways, PI3K-Akt signaling pathway, cytokine-cytokine receptor interaction, and MAPK signaling pathway. The PI3K-Akt signaling pathway has been reported to be involved in the mechanical force-induced osteoblast differentiation of PDLSCs.41,42 The MAPK signaling pathway also have been found to participate in the mechano-transduction of PDLSCs.43,44
Genes in profile 0 showed a continuous downward trend from 0.1 Hz to 0.7 Hz, an inverse trend of osteogenic genes, and were mainly enriched for pathways related to an inflammatory response, such as arachidonic acid metabolism and cytokinecytokine receptor interaction. Cyclooxygenase-2 (COX2) and prostaglandin E synthase (PTGES) participate in the arachidonic acid metabolism pathway. TNF Receptor Superfamily Member 14 (TNFRSF14) and Interleukin 9 (IL9) are involved in the cytokinecytokine receptor interaction. PTGES is induced by inflammatory mediators.45 COX2 is involved in the synthesis of prostaglandin E2 (PGE2), which is a potent pro-inflammatory cytokine, and participates in bone resorption.2 The expression of COX2 and PGE2 after mechanical stimulation was previously reported to be correlated with force duration and force magnitude,9 and it showed a negative correlation with tensile frequency in the present study. TNFRSF14 is a membrane-bound receptor leading to the induction of proinflammatory genes by activating the NF-B pathway.46 IL-9 plays a role in regulating inflammatory immunity, and it has demonstrated pro-inflammatory activity in several mouse models of inflammation.47 The response of periodontal ligament to mechanical stress generated by OTM is known as an aseptic transitory inflammatory process, which is regulated by various cytokines and chemokines.48 Proinflammatory cytokines activate matrix metalloproteinases (MMPs), degrade the ECM, and inhibit the expression of COL-I.49 Increased osteogenesis is usually accompanied by lower levels of inflammatory cytokines and chemokines.10,50,51 Accordingly, in the present study, with increasing frequency, the osteogenic commitment increased, and the suppression of pro-inflammatory genes and the relative inflammatory response pathway were observed.
Furthermore, through PPI network screening, we identified two candidate genes, EYA1 and SALL1, which were specifically sensitized to tensile frequency. The result was validated by RTqPCR, which confirmed the decreased expression of EYA1 and SALL1 with increasing frequency of tension stimulation. EYA1 is a conserved critical regulator of organ-specific stem cells.52 SALL1 is also considered a stem cell marker.53 The osteoblastic differentiation of PDLSCs increased with increasing tensile frequency; thus, stemness and related genes correspondingly reduced. In view of the high-throughput sequencing and validation after 6 h of tensile force exposure, further studies over a longer period are needed. Whether overexpression of EYA1 and SALL1 would reverse the frequency-dependent trend of the osteogenic differentiation of PDLSCs also deserves further study.
In the present study, the role of tensile frequency on osteogenic commitment of PDLSCs were identified, and the mRNA transcriptomes of PDLSCs during the osteogenic differentiation under cyclic tension with different frequencies were delineated. Frequency series clustering were defined using STEM, and tensile frequency-sensitive genes were identified. This study extends the knowledge about the role of tensile frequency in cyclic tension induced PDLSCs osteogenesis.
The osteoblastic differentiation of PDLSCs under mechanical tensile force is frequency dependent. EYA1 and SALL1 were identified as potential important tensile frequency-sensitive genes, which may contribute to the cyclic tension-induced osteogenic differentiation of PDLSCs in a frequency-dependent manner.
Raw data of the performed RNA sequencing (RNA-seq) have been recorded in the SRA database with the SRA accession: PRJNA665587 (https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA665587). Other data in this study are available from the corresponding author Xi Chen upon request.
The study was performed in accordance with the principles stated in the Declaration of Helsinki and approved by the Medical Ethics Committee of the First Affiliated Hospital of Medical College of Xi an Jiaotong University (No: XJTU1AF2019LSK-078). Informed consent was obtained from all donors and their legal guardians involved in the study. Written informed consent was also obtained from the donors and their legal guardians to publish this paper.
All authors made a significant contribution to the work reported, whether in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising, or critically reviewing the article; gave final approval of the version to be published; agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.
This research was funded by the Key Research and Development project of Shaanxi Province under Grant 2018SF-037.
The authors declare that they have no competing interests.
1. Wang H, Sun W, Ma J, Pan Y, Wang L, Zhang W. Polycystin-1 mediates mechanical strain-induced osteoblastic mechanoresponses via potentiation of intracellular calcium and Akt/-catenin pathway. PLoS One. 2014;9(3):e91730. doi:10.1371/journal.pone.0091730
2. Li M, Zhang C, Yang Y. Effects of mechanical forces on osteogenesis and osteoclastogenesis in human periodontal ligament fibroblasts: a systematic review of in vitro studies. Bone Joint Res. 2019;8(1):1931. doi:10.1302/2046-3758.81.Bjr-2018-0060.R1
3. Chang M, Lin H, Fu H, Wang B, Han G, Fan M. MicroRNA-195-5p regulates osteogenic differentiation of periodontal ligament cells under mechanical loading. J Cell Physiol. 2017;232(12):37623774. doi:10.1002/jcp.25856
4. Yang Y, Wang BK, Chang ML, Wan ZQ, Han GL. Cyclic stretch enhances osteogenic differentiation of human periodontal ligament cells via YAP activation. Biomed Res Int. 2018;2018:2174824. doi:10.1155/2018/2174824
5. Jin SS, He DQ, Wang Y, et al. Mechanical force modulates periodontal ligament stem cell characteristics during bone remodelling via TRPV4. Cell Prolif. 2020;53(10):e12912. doi:10.1111/cpr.12912
6. Katz S, Boland R, Santilln G. Modulation of ERK 1/2 and p38 MAPK signaling pathways by ATP in osteoblasts: involvement of mechanical stress-activated calcium influx, PKC and Src activation. Int J Biochem Cell Biol. 2006;38(12):20822091. doi:10.1016/j.biocel.2006.05.018
7. Ren D, Wei F, Hu L, Yang S, Wang C, Yuan X. Phosphorylation of Runx2, induced by cyclic mechanical tension via ERK1/2 pathway, contributes to osteodifferentiation of human periodontal ligament fibroblasts. J Cell Physiol. 2015;230(10):24262436. doi:10.1002/jcp.24972
8. Xu J, Li Z, Hou Y, Fang W. Potential mechanisms underlying the Runx2 induced osteogenesis of bone marrow mesenchymal stem cells. Am J Transl Res. 2015;7(12):25272535.
9. Sun C, Janjic Rankovic M, Folwaczny M, Otto S, Wichelhaus A, Baumert U. Effect of tension on human periodontal ligament cells: systematic review and network analysis. Front Bioeng Biotechnol. 2021;9:695053. doi:10.3389/fbioe.2021.695053
10. Sun C, Janjic Rankovic M, Folwaczny M, et al. Effect of different parameters of in vitro static tensile strain on human periodontal ligament cells simulating the tension side of orthodontic tooth movement. Int J Mol Sci. 2022;23(3). doi:10.3390/ijms23031525
11. Natali AN, Pavan PG, Scarpa C. Numerical analysis of tooth mobility: formulation of a non-linear constitutive law for the periodontal ligament. Dent Mater. 2004;20(7):623629. doi:10.1016/j.dental.2003.08.003
12. Pinkerton MN, Wescott DC, Gaffey BJ, Beggs KT, Milne TJ, Meikle MC. Cultured human periodontal ligament cells constitutively express multiple osteotropic cytokines and growth factors, several of which are responsive to mechanical deformation. J Periodontal Res. 2008;43(3):343351. doi:10.1111/j.1600-0765.2007.01040.x
13. Liu J, Li Q, Liu S, et al. Periodontal ligament stem cells in the periodontitis microenvironment are sensitive to static mechanical strain. Stem Cells Int. 2017;2017:1380851. doi:10.1155/2017/1380851
14. Tang N, Zhao Z, Zhang L, et al. Up-regulated osteogenic transcription factors during early response of human periodontal ligament stem cells to cyclic tensile strain. Arch Med Sci. 2012;8(3):422430. doi:10.5114/aoms.2012.28810
15. Tantilertanant Y, Niyompanich J, Everts V, Supaphol P, Pavasant P, Sanchavanakit N. Cyclic tensile force stimulates BMP9 synthesis and in vitro mineralization by human periodontal ligament cells. J Cell Physiol. 2019;234(4):45284539. doi:10.1002/jcp.27257
16. Shen T, Qiu L, Chang H, et al. Cyclic tension promotes osteogenic differentiation in human periodontal ligament stem cells. Int J Clin Exp Pathol. 2014;7(11):78727880.
17. Wu X, Li Y, Cao Z, Xie Y, Fu C, Chen H. Mechanism of cyclic tensile stress in osteogenic differentiation of human periodontal ligament stem cells. Calcif Tissue Int. 2021;108(5):640653. doi:10.1007/s00223-020-00789-x
18. Wang Y, Hu B, Hu R, et al. TAZ contributes to osteogenic differentiation of periodontal ligament cells under tensile stress. J Periodontal Res. 2020;55(1):152160. doi:10.1111/jre.12698
19. Wang H, Feng C, Li M, Zhang Z, Liu J, Wei F. Analysis of lncRNAs-miRNAs-mRNAs networks in periodontal ligament stem cells under mechanical force. Oral Dis. 2021;27(2):325337. doi:10.1111/odi.13530
20. Yang L, Yang Y, Wang S, Li Y, Zhao Z. In vitro mechanical loading models for periodontal ligament cells: from two-dimensional to three-dimensional models. Arch Oral Biol. 2015;60(3):416424. doi:10.1016/j.archoralbio.2014.11.012
21. Hsieh YF, Turner CH. Effects of loading frequency on mechanically induced bone formation. J Bone Miner Res. 2001;16(5):918924. doi:10.1359/jbmr.2001.16.5.918
22. Scheuren AC, Vallaster P, Kuhn GA, et al. Mechano-regulation of trabecular bone adaptation is controlled by the local in vivo environment and logarithmically dependent on loading frequency. Front Bioeng Biotechnol. 2020;8:566346. doi:10.3389/fbioe.2020.566346
23. Yang SY, Wei FL, Hu LH, Wang CL. PERK-eIF2-ATF4 pathway mediated by endoplasmic reticulum stress response is involved in osteodifferentiation of human periodontal ligament cells under cyclic mechanical force. Cell Signal. 2016;28(8):880886. doi:10.1016/j.cellsig.2016.04.003
24. Borodina T, Adjaye J, Sultan M. A strand-specific library preparation protocol for RNA sequencing. Methods Enzymol. 2011;500:7998. doi:10.1016/b978-0-12-385118-5.00005-0
25. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. doi:10.1186/s13059-014-0550-8
26. Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinform. 2006;7:191. doi:10.1186/1471-2105-7-191
27. Brulois K, Rajaraman A, Szade A, et al. A molecular map of murine lymph node blood vascular endothelium at single cell resolution. Nat Commun. 2020;11(1):3798. doi:10.1038/s41467-020-17291-5
28. Stanton BZ, Hodges C, Calarco JP, et al. Smarca4 ATPase mutations disrupt direct eviction of PRC1 from chromatin. Nat Genet. 2017;49(2):282288. doi:10.1038/ng.3735
29. Beauchaine TP, Beauchaine RJ. A comparison of maximum covariance and K-means cluster analysis in classifying cases into known taxon groups. Psychol Methods. 2002;7(2):245261. doi:10.1037/1082-989x.7.2.245
30. Iacono G, Dubos A, Mziane H, et al. Increased H3K9 methylation and impaired expression of Protocadherins are associated with the cognitive dysfunctions of the Kleefstra syndrome. Nucleic Acids Res. 2018;46(10):49504965. doi:10.1093/nar/gky196
31. Tantilertanant Y, Niyompanich J, Everts V, Supaphol P, Pavasant P, Sanchavanakit N. Cyclic tensile force-upregulated IL6 increases MMP3 expression by human periodontal ligament cells. Arch Oral Biol. 2019;107:104495. doi:10.1016/j.archoralbio.2019.104495
32. Symmank J, Zimmermann S, Goldschmitt J, et al. Mechanically-induced GDF15 secretion by periodontal ligament fibroblasts regulates osteogenic transcription. Sci Rep. 2019;9(1):11516. doi:10.1038/s41598-019-47639-x
33. Yu N, Prodanov L, Te Riet J, et al. Regulation of periodontal ligament cell behavior by cyclic mechanical loading and substrate nanotexture. J Periodontol. 2013;84(10):15041513. doi:10.1902/jop.2012.120513
34. Wei FL, Wang JH, Ding G, et al. Mechanical force-induced specific MicroRNA expression in human periodontal ligament stem cells. Cells Tissues Organs. 2014;199(56):353363. doi:10.1159/000369613
35. Jnsson D, Nebel D, Bratthall G, Nilsson BO. The human periodontal ligament cell: a fibroblast-like cell acting as an immune cell. J Periodontal Res. 2011;46(2):153157. doi:10.1111/j.1600-0765.2010.01331.x
36. Livne A, Bouchbinder E, Geiger B. Cell reorientation under cyclic stretching. Nat Commun. 2014;5:3938. doi:10.1038/ncomms4938
37. Rather HA, Jhala D, Vasita R. Dual functional approaches for osteogenesis coupled angiogenesis in bone tissue engineering. Mater Sci Eng C Mater Biol Appl. 2019;103:109761. doi:10.1016/j.msec.2019.109761
38. Jacobs C, Grimm S, Ziebart T, Walter C, Wehrbein H. Osteogenic differentiation of periodontal fibroblasts is dependent on the strength of mechanical strain. Arch Oral Biol. 2013;58(7):896904. doi:10.1016/j.archoralbio.2013.01.009
39. Chen G, Goeddel DV. TNF-R1 signaling: a beautiful pathway. Science. 2002;296(5573):16341635. doi:10.1126/science.1071924
40. Arden KC. FoxO: linking new signaling pathways. Mol Cell. 2004;14(4):416418. doi:10.1016/s1097-2765(04)00213-8
41. Qi L, Zhang Y. The microRNA 132 regulates fluid shear stress-induced differentiation in periodontal ligament cells through mTOR signaling pathway. Cell Physiol Biochem. 2014;33(2):433445. doi:10.1159/000358624
42. Jiang N, He D, Ma Y, et al. Force-induced autophagy in periodontal ligament stem cells modulates M1 macrophage polarization via AKT signaling. Front Cell Dev Biol. 2021;9:666631. doi:10.3389/fcell.2021.666631
43. Ziegler N, Alonso A, Steinberg T, et al. Mechano-transduction in periodontal ligament cells identifies activated states of MAP-kinases p42/44 and p38-stress kinase as a mechanism for MMP-13 expression. BMC Cell Biol. 2010;11:10. doi:10.1186/1471-2121-11-10
44. Tang M, Peng Z, Mai Z, et al. Fluid shear stress stimulates osteogenic differentiation of human periodontal ligament cells via the extracellular signal-regulated kinase 1/2 and p38 mitogen-activated protein kinase signaling pathways. J Periodontol. 2014;85(12):18061813. doi:10.1902/jop.2014.140244
45. Blair PJ, Hwang SJ, Shonnard MC, et al. The role of prostaglandins in disrupted gastric motor activity associated with type 2 diabetes. Diabetes. 2019;68(3):637647. doi:10.2337/db18-1064
46. Shui JW, Steinberg MW, Kronenberg M. Regulation of inflammation, autoimmunity, and infection immunity by HVEM-BTLA signaling. J Leukoc Biol. 2011;89(4):517523. doi:10.1189/jlb.0910528
47. Goswami R, Kaplan MH. A brief history of IL-9. J Immunol. 2011;186(6):32833288. doi:10.4049/jimmunol.1003049
48. Lee SI, Park KH, Kim SJ, Kang YG, Lee YM, Kim EC. Mechanical stress-activated immune response genes via Sirtuin 1 expression in human periodontal ligament cells. Clin Exp Immunol. 2012;168(1):113124. doi:10.1111/j.1365-2249.2011.04549.x
49. Sun C, Liu F, Cen S, et al. Tensile strength suppresses the osteogenesis of periodontal ligament cells in inflammatory microenvironments. Mol Med Rep. 2017;16(1):666672. doi:10.3892/mmr.2017.6644
50. Loi F, Crdova LA, Pajarinen J, Lin TH, Yao Z, Goodman SB. Inflammation, fracture and bone repair. Bone. 2016;86:119130. doi:10.1016/j.bone.2016.02.020
51. Orapiriyakul W, Tsimbouri MP, Childs P, et al. Nanovibrational stimulation of mesenchymal stem cells induces therapeutic reactive oxygen species and inflammation for three-dimensional bone tissue engineering. ACS Nano. 2020;14(8):1002710044. doi:10.1021/acsnano.0c03130
52. Sun Y, Kaneko S, Li XK, The LX. PI3K/Akt signal hyperactivates Eya1 via the SUMOylation pathway. Oncogene. 2015;34(19):25272537. doi:10.1038/onc.2014.179
53. Griesche N, Luttmann W, Luttmann A, Stammermann T, Geiger H, Baer PC. A simple modification of the separation method reduces heterogeneity of adipose-derived stem cells. Cells Tissues Organs. 2010;192(2):106115. doi:10.1159/000289586
Read more:
Effect of tensile frequency on osteogenic differentiation | IJGM - Dove Medical Press
