Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2018.8517 mmr-17-04-5238 Articles Intermittent cyclic mechanical tension altered the microRNA expression profile of human cartilage endplate chondrocytes FengChencheng * LiuMinghan * FanXin YangMinghui LiuHuan ZhouYue Department of Orthopedics, Xinqiao Hospital, Third Military Medical University, Chongqing 400037, P.R. China Correspondence to: Professor Huan Liu or Professor Yue Zhou, Department of Orthopedics, Xinqiao Hospital, Third Military Medical University, 183 Xinqiao Main Street, Shapingba, Chongqing 400037, P.R. China, E-mail: 20016040@163.com, E-mail: happyzhou@vip.163.com

Contributed equally

 042018 30012018 17 4 5238 5246 14112016 03012018 Copyright: © Feng et al. 2018 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Previous studies have identified the association between cartilage endplate (CEP) degeneration and abnormal mechanical loading. Several studies have reported that intermittent cyclic mechanical tension (ICMT) regulates CEP degeneration via various biological processes and signaling pathways. However, the functions of microRNAs in regulating the cellular responses of CEP chondrocytes to ICMT remain to be elucidated. The current study determined the differentially expressed microRNAs in human CEP chondrocytes exposed to ICMT using microarray analysis. A total 21 significantly upregulated and 62 downregulated miRNAs were identified compared with the control. The findings were subsequently partially validated by reverse transcription-quantitative polymerase chain reaction. Potential target genes of the significantly differentially expressed miRNAs were predicted using bioinformatics analysis and were used for Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The present study revealed that the significantly differentially expressed microRNAs were involved in various signaling pathways and biological processes that are crucial to regulating the responses of CEP chondrocytes to ICMT. The current study provided a global view of microRNA expression in CEP chondrocytes under mechanical stimulation, suggesting that microRNAs are important for regulating the mechanical response of CEP chondrocytes. Additionally, it provided a novel insight into the association between mechanical stress and the establishment and progression of intervertebral disc degeneration.

 microRNA cyclic mechanical tension cartilage endplate chondrocytes cartilage endplate calcification intervertebral disc degeneration Introduction

Intervertebral disc degeneration (IDD) is a major contributor of low back pain (LBP) and has been recognized to be associated with various etiological factors, including aging, smoking, genetic predisposition, trauma, infection and abnormal mechanical loading (1–7). IDD is a disc cell-mediated pathological process (8). The primary cellular events in the process of IDD have been identified as disorder of extracellular matrix (ECM) metabolism in disc cells, aberrant production of matrix proteases, inflammatory cytokines and chemokines, cell senescence, apoptosis and cell death (9–12). It is of note that these cellular events have been identified to be strongly associated with a reduced nutrient supply in intervertebral disc (IVD) (13,14).

The nutrient supply of disc cells primarily depends on the diffusion of nutrient solutes from blood vessels through cartilage endplate (CEP) (15,16). However, the reduced permeability of CEP due to calcification impedes the availability of nutrients to disc cells and suppresses the clearance of metabolites. Therefore, the microenvironment of the disc becomes harsh, impairing the viability and activity of disc cells and consequently contributing to IDD (13,17). The causes of CEP calcification remain to be elucidated. The underlying molecular mechanism is not well understood. Previous studies have determined the association between CEP calcification and abnormal mechanical loading of the spine (18,19). Excessive mechanical loading has been identified to lead to the damage of chondrocytes and ECM in CEP (8). CEP chondrocytes, the basic units of CEP, have an important function in sensing and responding the mechanical stress loaded on CEP. Therefore, investigating the responses of CEP chondrocytes to mechanical loading may provide a novel insight into the mechanism of CEP calcification and the pathogenesis of IDD.

Cyclic mechanical tension (CMT) generated by the Flexcell system in vitro was widely used to simulate the mechanical strain which CEP chondrocytes undergo in vivo (20–24). The effects of CMT on the calcification of CEP chondrocytes have been discussed intensively in previous studies (20–26). It varies with the types of CMT: Intermittent CMT (ICMT, stimulation for several h each day in consecutive days) significantly upregulated the expression of calcification-associated genes, such as collagen type I, type X, osteocalcin and osteopontin and downregulated the expression of anti-calcification genes, including progressive ankylosis (ANK) gene, extracellular nucleotide phosphatase/phosphodiesterase 1 (ENPP1) and transforming growth factor-β1 (TGF-β1) leading to ICMT-induced calcification of CEP chondrocytes (21,22). Conversely, continuous CMT (CCMT, continuous stimulation for several days) upregulated the expression of progressive ankylosis (ANK) and TGF-β1 and protected CEP chondrocytes from calcification (20,23). Additionally, CMT has been identified to regulate the cartilage matrix metabolism, autophagy and cytoskeleton arrangement in CEP chondrocytes (22,24–26). In summary, the response of CEP chondrocytes to CMT was complicated, and involved in complex biological processes and signaling pathways. However, the regulatory mechanism underlying the response of CEP chondrocytes to CMT remains to be elucidated.

MicroRNAs (miRNAs) are members of the family of small non-coding RNAs and are 20–22 nucleotides in length. They bind to their target mRNAs in order to trigger degradation or to suppress their translation (27). MiRNAs have been determined to regulate various biological processes, such as cell proliferation, differentiation, apoptosis, senescence and development (27). Additionally, the involvement of miRNAs in the development and progression of various diseases has been discussed in previous studies, including carcinoma, neurodegenerative disorders, inflammatory diseases and reproductive disorders (28–34). The importance of miRNAs in musculoskeletal disorders, such as IDD, bone tumors, osteoarthritis and osteonecrosis of the femoral head, have been previously reported (35–40). However, the miRNA expression profile of CEP chondrocytes under CMT stimulation remains to be elucidated. Elucidating the altered miRNA expression profile of CEP chondrocytes under CMT stimulation helps us understand the mechanism underlying the regulatory effects of CMT on the biological processes and signaling pathways in CEP chondrocytes.

The present study aimed to determine the differentially expressed miRNAs of CEP chondrocytes under ICMT stimulation. CEP chondrocytes were isolated from human CEP specimens. Subsequently, ICMT stimulation was applied and the total RNA of CEP chondrocytes was extracted for miRNA microarray. Bioinformatics analysis was performed to determine the significantly differentially expressed miRNAs and their target genes. The signaling pathways and biological processes in which these target genes were involved were predicted. The findings of the present study identified the miRNA expression profile associated with the biological responses of CEP under ICMT stimulation, which may aid in identifying the mechanism underlying the mechanical force-induced changes of CEP. Additionally, the present study improved the clinical understanding of the pathogenesis of IDD.

 Materials and methods <sec> <title>Ethical statement

All experimental procedures for human IVD samples were approved by the Ethical Committee of Xinqiao Hospital (Chongqing, China) and were in consistent with the Declaration of Helsinki. Signed informed consent was obtained from the participating patients.

 Human CEP specimens

CEP specimens used in this study were obtained from patients who underwent posterior discectomy and a fusion procedure for IDD at the orthopedics department of Xinqiao Hospital between September 2015 and November 2015. The degenerative changes of CEP were evaluated using the method previously described by Thompson et al (41). In order to reduce the bias resulting from age and different CEP degenerative degrees, the present study selected the patients between the ages of 40 and 50. Additionally, the modic types of the selected patients were type 0, normal. The specimens from 6 patients were used in the present study. Detailed information about the 6 patients and samples is presented in Table I.

 CEP chondrocyte culture

CEP chondrocytes were isolated from IVD specimens according to the protocol described in previous studies (42,43). IVD specimens were placed in a sterilized culture dish and washed with phosphate buffered saline (PBS). CEP tissues were carefully separated from IVD specimens under a dissecting microscope. Subsequently, CEP tissues were minced into 1-mm3 tissue blocks followed by digestion with Dulbecco's modified Eagle's medium (DMEM)/F-12 medium (Invitrogen; Thermo Fisher Scientific, Inc., Waltham, MA, USA) containing 0.2% type II collagenase (Sigma-Aldrich; Merck Millipore, Darmstadt, Germany) for 2 h at 37°C. The cell suspension was filtered with a 70 µm cell strainer and subsequently centrifuged for 10 min at 100 × g and then re-suspended in DMEM/F12 containing 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin (Invitrogen; Thermo Fisher Scientific, Inc.). The isolated CEP chondrocytes were cultured at 37°C and 5% CO2. The medium was replaced twice a week. When the cells reached about 80% confluence, they were sub-cultivated. The cells at the second passage were used in the experiments.

 ICMT treatment

Following trypsinization, CEP chondrocytes were seeded on 6-well flexible silicone rubber BioFlex plates coated with collagen type I (Flexcell International Corporation, Burlington, NC, USA) at a density of 1.0×105 cells/well. When the cells reached 80–90% confluence, ICMT at 0.5 Hz sinusoidal curve at of 10% elongation was applied for 4 h/day in 14 consecutive days using a FX-5000T Flexcell Tension Plus system (Flexcell International Corporation). The culture medium was replaced every 2 days. The cells were cultured in a humidified atmosphere at 37°C and 5% CO2 ICMT treatment. The morphology of the cells was observed and imaged using a phase contrast microscope (Olympus Corporation, Tokyo, Japan). The cells were harvested immediately after following ICMT treatment. The CEP chondrocytes from the same patients without ICMT stimulation were used as the control.

 RNA isolation and microRNA microarray

CEPs isolated from 3 patients (case no. 1, 2 and 3) were treated with ICMT and then termed the ICMT samples (n=3). The CEPs from the same patient without ICMT treatment were used as the control samples (n=3). Total RNA of all samples (n=6) lysed with 1 ml TRIzol reagent (Takara Bio Inc., Otsu, Japan) were sent to Kangchen Biotech Co., Ltd. (Shanghai, China) for microRNA microarray analysis. RNA quantity and quality were determined by a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Inc.). All samples met the quality control standards. RNA was labeled using a miRCURY Array Power Labeling kit (Exiqon, Woburn, MA, USA) and then was hybridized to a one-color miRCURY Array (Exiqon, Woburn). The slides were scanned with the Axon GenePix 4000B microarray scanner (Molecular Devices, LLC, Sunnyvale, CA, USA). After image scanning, GenePix Pro version 6.0 was used to read the raw intensity of the image. The intensity of green signal was calculated after background subtraction and four replicated spots of each probe on the same slide have been averaged. The median normalization method was used to obtain the ‘Normalized Data’, as follows: Normalized Data=(Foreground-Background)/median, the median is 50 percent quantile of microRNA intensity that is >30 in all samples after background correction. The difference in the normalized data of miRNAs between the ICMT group (n=3) and the control group (n=3) was statistically tested using a paired sample Student's t-test. The fold-change of the upregulated miRNAs was calculated using the following equation: Fold-change=ICMT normalized data/Control normalized data. The fold-change of the downregulated miRNAs was calculated as following: Fold-change=−(Control normalized data/the normalized data of the ICMT groups). The differentially expressed miRNAs were identified according to their fold-change (fold-change ≥2 or ≤-2). The significantly differentially expressed miRNAs were selected according to the following thresholds P<0.05 and fold-change ≥2 or ≤-2.

 Reversed transcription-quantitative polymerase chain reaction (RT-qPCR)

A total of 10 significantly differentially expressed miRNAs with relatively high fold-changes (fold-change ≥4 or ≤-4) were selected for validation using the RT-qPCR method, including 5 upregulated miRNAs (hsa-let-7a, hsa-miR-29c, hsa-miR-142, hsa-miR-181a, hsa-miR-19b) and 5 downregulated miRNAs (hsa-miR-548q, hsa-miR-637, hsa-miR-614, hsa-miR-581, hsv1-miR-H14). RNA from the CEPs of the patients (case no. 1, 2, 3, 4, 5 and 6) was reverse transcribed using a microRNA First Strand cDNA Synthesis kit (Sangon Biotech Co., Ltd., Shanghai, China) according to the manufacturer's protocol. The reaction was performed at 37°C for 60 min followed by 85°C for 5 min. qPCR was performed in triplicate on a StepOnePlus Real-Time PCR system (Applied Biosystems; Thermo Fisher Scientific, Inc.) with SYBR Premix Ex Taq™ II (Takara Bio) according to the protocol provided by the manufacturer. The following thermocycling conditions were used for the PCR: 95°C for 30 sec; 40 cycles of 95°C for 5 sec and 60°C for 34 sec. PCR products were subjected to melting curve analysis. Relative expression was calculated using the 2−ΔΔCq method (44). U6 was used as the internal reference. Mean Cq values were normalized to U6. The forward primers of the miRNAs investigated in the present study were listed in Table II. The reverse primers used for the miRNAs were universal PCR primer R.

 Bioinformatics analysis

The present study used three popular databases MiRBase (www.mirbase.org) (45), TargetScan (www.targetscan.org) (46) and MiRanda (www.microrna.org) (47) to predict the target genes of the significantly differentially expressed miRNAs. Genes that overlapped all three databases were selected as the target genes for further function annotation analysis. The Gene Ontology analysis (GO; www.geneontology.org) was used to annotate the functions of the target genes. The Kyoto Encyclopedia of Genes and Genomes pathway analysis (KEGG; www.genome.jp/kegg) was performed to identify the enriched signaling pathways by the target genes. P<0.05 was considered to indicate statistically significant difference.

 Statistical analysis

Data were presented as the mean ± standard error of the mean. Data from the RT-qPCR were statistically tested using Kruskal-Wallis non-parametric analysis and Mann-Whitney U post-hoc tests as previously described (48). The remaining data were analyzed using SPSS version 22.0 (IBM Corporation, Armonk, NY, USA) and GraphPad Prism version 6 (GraphPad Software Inc., La Jolla, CA, USA). P<0.05 was considered to indicate a statistically significant difference.

 Results <sec> <title>Morphological changes of CEP chondrocytes subjected to ICMT

The morphology of CEP chondrocytes subjected to ICMT changed from a polygonal morphology into a spindle morphology. Additionally, ICMT induced an orderly alignment of CEP chondrocytes with a certain direction (Fig. 1).

 Differentially expressed miRNAs of CEP chondrocytes induced by ICMT

All differentially expressed miRNAs were filtered according to their fold-changes (fold-change ≥2 or ≤-2). There were 86 upregulated miRNAs and 188 downregulated miRNAs compared with the control. From these, 83 miRNAs (21 upregulated and 62 downregulated) were statistically significant (P<0.05; Tables III and IV). To validate the microarray results, 10 significantly differentially expressed miRNAs with high fold-changes (fold-change ≥4 or ≤-4), including 5 upregulated miRNAs (hsa-let-7a, hsa-miR-29c, hsa-miR-142, hsa-miR-181a, hsa-miR-19b) and 5 downregulated miRNAs (hsa-miR-548q, hsa-miR-637, hsa-miR-614, hsa-miR-581, hsv1-miR-H14) were selected as representative miRNAs for RT-qPCR assays. With the exception of has-miR-614, these miRNAs were significantly regulated by ICMT, which was consistent with the microarray results (Fig. 2).

 KEGG pathway analysis

The target genes of the 83 significantly differentially expressed miRNAs were identified based on the overlapping of 3 databases (MiRBase, MiRanda and TargetScan). A total of 1,836 genes were identified to be targeted by the significantly regulated miRNAs. The KEGG pathways analysis of these target genes revealed that the genes targeted by significantly differentially expressed miRNAs were enriched in various signaling pathways. The top 20 signaling pathways regulated by the significantly differentially expressed miRNAs are presented in Fig. 3, including the mitogen-activated protein kinase (MAPK) and Wnt signaling pathways, focal adhesion, regulation of actin cytoskeleton, GnRH signaling pathway, cell cycle, ErbB signaling pathway, insulin signaling pathway, TGF-β signaling pathway and adherens junction.

 GO analysis

The significantly differentially expressed miRNAs were involved in different biological processes, including regulation of transcription, signal transduction, development, cell adhesion, cell cycle, negative regulation of cell proliferation and chromatin modification (Fig. 4). The top 20 biological processes enriched by the target genes of significantly differentially expressed miRNAs are presented in Fig. 4.

 Discussion

Previous studies have reported that ICMT markedly induced the calcification of CEP chondrocytes (21,22). Furthermore, CMT has been identified to regulate the matrix metabolism, autophagy and cytoskeleton arrangement of CEP chondrocytes through MAPK and Wnt/β-catenin pathways (22,24–26). Briefly, CMT is involved in the regulation of complex biological processes and signaling pathways of CEP chondrocytes. It is of note that CEP calcification is a major pathological characteristic of degenerative discs, as it impedes the availability of nutrients to disc cells and suppresses the clearance of metabolites (13,17). The presence of disc cell autophagy and the disturbance of the balance between matrix anabolism and catabolism in human degenerative disc cells has also been identified as causal factor of IDD (49–52). Therefore, mechanical tension is a widely recognized as a primary contributor of IDD. However, the mechanism underlying the effect of ICMT on the viability and function of CEP chondrocytes remains to be fully elucidated. The present study determined that 83 miRNAs were significantly differentially expressed in human CEP chondrocytes in response to ICMT. Additionally, GO and KEGG pathway analyses revealed that the differentially expressed miRNAs were widely involved in the regulation of various cellular processes and signaling pathways of human CEP chondrocytes, suggesting that they are crucial intermediators of the effects of ICMT on the viability and function of CEP chondrocytes. More importantly, the present study suggests that miRNA dysregulation has an important function in the establishment and progression of IDD through mediation of the effect of mechanical loadings on the structural and functional homeostasis of IVDs.

The target genes of the significantly differentially expressed miRNAs were enriched in various biological processes according to the GO analysis. It is of note that several GO terms were associated with regulation of cell transcription, including regulation of transcription (DNA-dependent), positive regulation of transcription from RNA polymerase II promoter, negative regulation of transcription from RNA polymerase II promoter and chromatin modification. The present findings suggest that there is a regulatory effect of ICMT on the transcriptome of human CEPs. This is consistent with the findings of previous studies (21,24). ICMT has been identified to regulate the transcription of ECM-associated genes, including TGF-β gene and the genes associated with calcification, including collagen type I, type X, osteocalcin, osteopontin, ANK and ENPP1. Subsequently, the CEP chondrocytes undergo a phenotype shift from a cartilage-specific to calcification-associated, which also promotes IDD progression. Based on the current findings, miRNAs are crucial regulators of CEP chondrocyte transcription under ICMT stimulation. Therefore, it is possible that this ICMT-induced phenotype shift is mediated by miRNAs, indicating that miRNA dysregulation contributes to the initiation and progression of CEP calcification. The effect of miRNAs on the transcriptome of human CEP chondrocytes should be discussed in further studies. In addition, there were some GO terms identified to be associated with cell signal transduction, including signal transduction, protein amino acid phosphorylation, modification-dependent protein catabolism and protein amino acid dephosphorylation. Considering the essential function of signal transduction in the regulation of cell viability and function, the current findings indicate that miRNAs affect the biological responses of CEP chondrocytes to mechanical stimulation by regulation of signal transduction pathways in CEP chondrocytes. However, the function of miRNAs in regulating mechanical force-induced signal transduction remains unclear and should be investigated in future research.

The KEGG pathway analysis was used to annotate the functional roles of significantly differentially expressed miRNAs in regulating the biological responses of CEP chondrocytes to ICMT. When exposed to ICMT, human CEP chondrocytes changed from a randomly distributed polygonal morphology into a spindle morphology with a specific directional alignment, which was consistent with the findings of previous studies (20,26). Multiple mechanotransduction pathways, including focal adhesion pathway, action cytoskeleton regulation pathway and adherens junction pathway were previously identified to be responsible for these changes (53–55). The findings of the present study predicted that these signaling pathways were regulated by the significantly differentially expressed miRNAs under ICMT stimulation, suggesting that miRNAs may regulate the migration and morphology of CEP cells via these signaling pathways.

The proliferation of CEP chondrocytes is crucial for the maintenance of the number of functional cells in CEP, which is regulated by cell cycle progression. In the current study, some GO terms and KEGG pathways enriched by the target genes of miRNAs were associated with cell cycle and negative regulation of cell proliferation, suggesting that ICMT regulates cell cycle progression to affect the proliferation of CEP chondrocytes via miRNAs. ICMT has been previously reported to promote the proliferation of CEP chondrocytes (26). However, the mechanism underlying this effect had not been previously investigated. The present study investigated miRNA regulation in order to elucidate the effect of mechanical stimulation on CEP chondrocyte proliferation.

TGF-β1 is an established regulator of cell proliferation, differentiation and ECM metabolism. Previous studies have demonstrated the function of TGF-β1 in the ICMT-induced calcification of CEP chondrocytes (21,22). TGF-β1 induced the expression of ANK and ENPP1 to slow down the process of CEP calcification. However, ICMT reduced the expression of TGF-β1 and in turn reduced the expression of ANK and ENPP1 leading to calcification of the CEP cells (21,22). Additionally, TGF-β1 has been demonstrated to promote the ECM anabolism of disc cells (56). Previous studies have determined that the downregulation of collagen type II and aggrecan was accompanied by the downregulation of TGF-regulationnion of collagen type he ECM anabolis (21,26). It is of note that the present study identified an association between the significantly differentially expressed miRNAs and the TGF-β signaling pathway, indicating the roles of differentially expressed miRNAs in the pathogenesis of IDD through regulation of the matrix metabolism and calcification of CEP chondrocytes via the TGF-β signaling pathway.

The activation of Wnt pathway has been identified in human degenerative CEP specimens (26). ICMT has been identified to activate the Wnt signaling pathway in order to induce the loss of the chondrogenic phenotype of human CEP chondrocytes (26). Furthermore, the activation of the Wnt signaling pathway has been found to promote IDD by inducing nucleus pulposus (NP) cell senescence and upregulating the expression of matrix metalloproteinase-3 (MMP-3), MMP-7, MMP-9 and MMP-10 in NP cells (48). Therefore, the Wnt signaling pathway is involved in the occurrence and development of IDD. Additionally, the KEGG pathway analysis performed in the present study suggested that the Wnt signaling pathway is associated with the significantly differentially expressed miRNAs under ICMT stimulation. The function of the Wnt pathway in regulating biological responses of CEP chondrocytes under ICMT provided novel insights into the involvement of miRNA dysregulation in the pathogenesis of IDD.

The MAPK signaling pathway has been previously revealed to regulate the expression of ANK and be responsible for mechanotransduction (57–59). CMT has been previously reported to activate the p38-MAPK signaling pathway to regulate the expression of ANK in CEP chondrocytes (23), indicating the crucial role of the MAPK signaling pathway in the regulation of ICMT-induced calcification of CEP chondrocytes. According the findings of the current study, the MAPK signaling pathway may be regulated by the differentially expressed miRNA under ICMT stimulation. The mechanism underlying the regulatory effects of miRNAs on the MAPK signaling pathway will be discussed in future studies.

There are two major limitations to the current study. First one is that the sample size of miRNA microarray analysis was relatively small. Therefore, more RT-qPCR validation is required to confirm the effect of ICMT on the expression of miRNAs in human CEP chondrocytes. Secondly, further studies using dysregulated miRNAs are necessary to elucidate the effects of miRNAs on the biological functions of CEP chondrocytes in detail.

In conclusion, the present study identified the differentially expressed miRNAs of human CEP chondrocytes under ICMT stimulation. The bioinformatics analyses suggest that miRNAs are essential for the regulation of the biological responses of CEP chondrocytes to ICMT, particularly in terms of the ICMT-induced calcification of CEP chondrocytes. The present study improved the current understanding of the involvement of miRNAs in the pathogenesis of CEP calcification, suggesting the possible function of miRNA dysregulation in the pathogenesis of IDD. Determining the detailed regulatory functions of these differentially expressed miRNAs contributes to identifying novel approaches to protect CEP cells from calcification under abnormal mechanical stress and subsequently suppressing the process of IDD.

 Acknowledgements

The current study was supported by the National Natural Science Foundation of China (grant nos. 81271982, 81472076 and 81572186).

 Abbreviations IVD

intervertebral disc

 IDD

intervertebral disc degeneration

 NP

nucleus pulposus

 AF

annulus fibrosus

 ICMT

intermittent cyclic mechanical tension

 ECM

extracellular matrix

 ANK

progressive ankylosis

 ENPP1

extracellular nucleotide phosphatase/phosphodiesterase 1

 TGF-β1

transforming growth factor-β1

 CMT

cyclic mechanical tension

 LBP

lower back pain

 CEP

cartilage endplate

 miRNA

microRNA

 MMP

matrix metalloproteinase

 FBS

fetal bovine serum

 PBS

phosphate-buffered saline

 References RobertsSEvansHTrivediJMenageJHistology and pathology of the human intervertebral discJ Bone Joint Surg Am88Suppl 2S10S14200610.2106/00004623-200604002-00003 KanayamaMTogawaDTakahashiCTeraiTHashimotoTCross-sectional magnetic resonance imaging study of lumbar disc degeneration in 200 healthy individualsJ Neurosurg Spine11501507200910.3171/2009.5.SPINE0867519929349 CheungKMKarppinenJChanDHoDWSongYQShamPCheahKSLeongJCLukKDPrevalence and pattern of lumbar magnetic resonance imaging changes in a population study of one thousand forty-three individualsSpine (Phila Pa 1976)34934940200910.1097/BRS.0b013e3181a01b3f19532001 BattieMCVidemanTKaprioJGibbonsLEGillKManninenHSaarelaJPeltonenLThe twin spine study: Contributions to a changing view of disc degenerationSpine J94759200910.1016/j.spinee.2008.11.01119111259 VergroesenPPKingmaIEmanuelKSHoogendoornRJWeltingTJvan RoyenBJvan DieënJHSmitTHMechanics and biology in intervertebral disc degeneration: A vicious circleOsteoarthritis Cartilage2310571070201510.1016/j.joca.2015.03.02825827971 WangDNastoLARoughleyPLemeASHoughtonAMUsasASowaGLeeJNiedernhoferLShapiroSSpine degeneration in a murine model of chronic human tobacco smokersOsteoarthritis Cartilage20896905201210.1016/j.joca.2012.04.010225314583389285 StirlingAWorthingtonTRafiqMLambertPAElliottTSAssociation between sciatica and Propionibacterium acnesLancet35720242025200110.1016/S0140-6736(00)05109-611438138 AdamsMARoughleyPJWhat is intervertebral disc degeneration and what causes it?Spine (Phila Pa 1976)3121512161200610.1097/01.brs.0000231761.73859.2c16915105 RisbudMVShapiroIMRole of cytokines in intervertebral disc degeneration: Pain and disc contentNat Rev Rheumatol104456201410.1038/nrrheum.2013.16024166242 VoNVHartmanRAYurubeTJacobsLJSowaGAKangJDExpression and regulation of metalloproteinases and their inhibitors in intervertebral disc aging and degenerationSpine J13331341201310.1016/j.spinee.2012.02.027233694953637842 JiangLZhangXZhengXRuANiXWuYTianNHuangYXueEWangXXuHApoptosis, senescence and autophagy in rat nucleus pulposus cells: Implications for diabetic intervertebral disc degenerationJ Orthop Res31692702201310.1002/jor.2228923238821 GruberHEIngramJANortonHJHanleyENHanleyENJrSenescence in cells of the aging and degenerating intervertebral disc: Immunolocalization of senescence-associated beta-galactosidase in human and sand rat discsSpine (Phila Pa 1976)32321327200710.1097/01.brs.0000253960.57051.de17268263 HuangYCUrbanJPLukKDIntervertebral disc regeneration: Do nutrients lead the way?Nat Rev Rheumatol10561566201410.1038/nrrheum.2014.9124914695 VoNVHartmanRAPatilPRRisbudMVKletsasDIatridisJCHoylandJALe MaitreCLSowaGAKangJDMolecular mechanisms of biological aging in intervertebral discsJ Orthop Res3412891306201610.1002/jor.23195268902034988945 UrbanJPSmithSFairbankJCNutrition of the intervertebral discSpine (Phila Pa 1976)2927002709200410.1097/01.brs.0000146499.97948.5215564919 GrunhagenTShirazi-AdlAFairbankJCUrbanJPIntervertebral disk nutrition: A review of factors influencing concentrations of nutrients and metabolitesOrthop Clin North Am42465477201110.1016/j.ocl.2011.07.01021944584 BibbySRUrbanJPEffect of nutrient deprivation on the viability of intervertebral disc cellsEur Spine J13695701200410.1007/s00586-003-0616-x150485603454063 BianQLiangQQWanCHouWLiCGZhaoYJLuSShiQWangYJProlonged upright posture induces calcified hypertrophy in the cartilage end plate in rat lumbar spineSpine (Phila Pa 1976)3620112020201110.1097/BRS.0b013e3181ffde9221304433 PengBHouSShiQJiaLThe relationship between cartilage end-plate calcification and disc degeneration: An experimental studyChin Med J (Engl)114308312200111780320 XuHGHuCJWangHLiuPYangXMZhangYWangLTEffects of mechanical strain on ANK, ENPP1 and TGF-β1 expression in rat endplate chondrocytes in vitroMol Med Rep4831835201121674130 XuHGZhangXHWangHLiuPWangLTZuoCJTongWXZhangXLIntermiåttent cyclic mechanical tension-induced calcification and downregulation of ankh gene expression of end plate chondrocytesSpine (Phila Pa 1976)3711921197201210.1097/BRS.0b013e318244d98922695244 XuHGLiZRWangHLiuPWangLTZuo CJ TongWXZhangXLIntermittent cyclic mechanical tension-induced down-regulation of ectonucleotide pyrophosphatase phosphodiesterase 1 gene expression is mainly dependent on TGF-β1 in end-plate chondrocytesOrthop Surg54045201310.1111/os.1202823420746 XuHZhangXWangHZhangYShiYZhangXContinuous cyclic mechanical tension increases ank expression in endplate chondrocytes through the TGF-β1 and p38 pathwayEur J Histochem57e28201310.4081/ejh.2013.e28240852773794359 XuHGYuYFZhengQZhangWWangCDZhaoXYTongWXWangHLiuPZhangXLAutophagy protects end plate chondrocytes from intermittent cyclic mechanical tension induced calcificationBone66232239201410.1016/j.bone.2014.06.01824970040 XuHGMaMMZhengQShenXWangHZhangSFXuJJWangCDZhangXLP120-catenin protects endplate chondrocytes from intermittent cyclic mechanical tension induced degeneration by inhibiting the expression of RhoA/ROCK-1 signaling pathwaySpine (Phila Pa 1976)4112611271201610.1097/BRS.000000000000153226913467 XuHGZhengQSongJXLiJWangHLiuPWangJWangCDZhangXLIntermittent cyclic mechanical tension promotes endplate cartilage degeneration via canonical Wnt signaling pathway and E-cadherin/β-catenin complex cross-talkOsteoarthritis Cartilage24158168201610.1016/j.joca.2015.07.01926247612 AmbrosVmicroRNAs: Tiny regulators with great potentialCell107823826200110.1016/S0092-8674(01)00616-X11779458 CroceCMCauses and consequences of microRNA dysregulation in cancerNat rev Genet10704714200910.1038/nrg2634197631533467096 TeagueEMPrintCGHullMLThe role of microRNAs in endometriosis and associated reproductive conditionsHum Reprod Update16142165201010.1093/humupd/dmp03419773286 AlexandrovPNDuaPLukiwWJUp-Regulation of miRNA-146a in progressive, age-related inflammatory neurodegenerative disorders of the human CNSFront Neurol5181201410.3389/fneur.2014.00181253248234179622 LvZShiQHuangWXingCHaoYFengXYangYZhangAKongQYukiNWangYMicroRNA expression profiling in Guillain-Barre syndromeJ Neuroimmunol3011215201610.1016/j.jneuroim.2016.10.01427836180 MenYFanYShenYLuLKallenANThe Steroidogenic acute regulatory protein (StAR) is regulated by the H19/let-7 AxisEndocrinology15840240920165413078 LuLKatsarosDCanutoEMBigliaNRischHAYuHLIN-28B/let-7a/IGF-II axis molecular subtypes are associated with epithelial ovarian cancer prognosisGynecol Oncol141121127201610.1016/j.ygyno.2015.12.03526751131 LuLKatsarosDRischHACanutoEMBigliaNYuHMicroRNA let-7a modifies the effect of self-renewal gene HIWI on patient survival of epithelial ovarian cancerMol Carcinog55357365201610.1002/mc.2228525630839 Ohrt-NissenSDøssingKBRossingMLajerCVikesåJNielsenFCFriis-HansenLDahlBCharacterization of miRNA expression in human degenerative lumbar disksConnect Tissue Res54197203201310.3109/03008207.2013.78159423586579 ZhaoBYuQLiHGuoXHeXCharacterization of microRNA expression profiles in patients with intervertebral disc degenerationInt J Mol Med334350201410.3892/ijmm.2013.154324173697 YuanHFVon RoemelingCGaoHDZhangJGuoCAYanZQAnalysis of altered microRNA expression profile in the reparative interface of the femoral head with osteonecrosisExp Mol Pathol98158163201510.1016/j.yexmp.2015.01.00225612520 MiyakiSAsaharaHMacro view of microRNA function in osteoarthritisNat Rev Rheumatol8543552201210.1038/nrrheum.2012.128228902453572197 WangXQianWWuZBianYWengXPreliminary screening of differentially expressed circulating microRNAs in patients with steroidinduced osteonecrosis of the femoral headMol Med Rep1031183124201410.3892/mmr.2014.266025323050 LiuXCheLXieYKHuQJMaCJPeiYJWuZGLiuZHFanLYWangHQNoncoding RNAs in human intervertebral disc degeneration: An integrated microarray studyGenom Data58081201510.1016/j.gdata.2015.05.027264842304583642 ThompsonKJDagherAPEckelTSClarkMReinigJWModic changes on MR images as studied with provocative diskography: Clinical relevance-a retrospective study of 2457 disksRadiology250849855200910.1148/radiol.250308047419244050 XiongCJHuangBZhouYCunYPLiuLTWangJLiCQPanYWangHMacrophage migration inhibitory factor inhibits the migration of cartilage end plate-derived stem cells by reacting with CD74PLoS One7e43984201210.1371/journal.pone.0043984229528373428348 LiuLTHuangBLiCQZhuangYWangJZhouYCharacteristics of stem cells derived from the degenerated human intervertebral disc cartilage endplatePLoS One6e26285201110.1371/journal.pone.0026285220288473196539 LivakKJSchmittgenTDAnalysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) methodMethods25402408200110.1006/meth.2001.126211846609 Griffiths-JonesSGrocockRJvan DongenSBatemanAEnrightAJmiRBase: microRNA sequences, targets and gene nomenclatureNucleic Acids Res34D140D144200610.1093/nar/gkj11216381832 LewisBPShihIHJones-RhoadesMWBartelDPBurgeCBPrediction of mammalian microRNA targetsCell115787798200310.1016/S0092-8674(03)01018-314697198 JohnBEnrightAJAravinATuschlTSanderCMarksDSHuman MicroRNA targetsPLoS Biol2e363200410.1371/journal.pbio.002036315502875521178 HiyamaASakaiDRisbudMVTanakaMAraiFAbeKMochidaJEnhancement of intervertebral disc cell senescence by WNT/β-catenin signaling-induced matrix metalloproteinase expressionArthritis Rheum6230363047201010.1002/art.27599205335443622204 YeWXuKHuangDLiangAPengYZhuWLiCAge-related increases of macroautophagy and chaperone-mediated autophagy in rat nucleus pulposusConnect Tissue Res52472478201110.3109/03008207.2011.56433621591930 JiangWZhangXHaoJShenJFangJDongWWangDZhangXShuiWLuoYSIRT1 protects against apoptosis by promoting autophagy in degenerative human disc nucleus pulposus cellsSci Rep47456201410.1038/srep07456255038524264007 GruberHEHoelscherGLIngramJABetheaSHanleyENJr.Autophagy in the degenerating human intervertebral disc: In vivo molecular and morphological evidence and induction of autophagy in cultured annulus cells exposed to proinflammatory cytokines-implications for disc degenerationSpine (Phila Pa 1976)40773782201510.1097/BRS.000000000000086526091153 KeplerCKPonnappanRKTannouryCARisbudMVAndersonDGThe molecular basis of intervertebral disc degenerationSpine J13318330201310.1016/j.spinee.2012.12.00323537454 EnglWArasiBYapLLThieryJPViasnoffVActin dynamics modulate mechanosensitive immobilization of E-cadherin at adherens junctionsNat Cell Biol16587594201410.1038/ncb297324859003 XuYWangQLiYGanYLiPLiSZhouYZhouQCyclic tensile strain induces tenogenic differentiation of tendon-derived stem cells in bioreactor cultureBiomed Res Int2015790804201510.1155/2015/790804262299624502284 TamielloCBuskermolenABCBaaijensFPTBroersJLVBoutenCVCHeading in the Right Direction: Understanding Cellular Orientation Responses to Complex Biophysical EnvironmentsCell Mol Bioeng91237201610.1007/s12195-015-0422-726900408 WalshAJBradfordDSLotzJCIn vivo growth factor treatment of degenerated intervertebral discsSpine (Phila Pa 1976)29156163200410.1097/01.BRS.0000107231.67854.9F14722406 DingLHeyingENicholsonNStroudNJHomandbergGABuckwalterJAGuoDMartinJAMechanical impact induces cartilage degradation via mitogen activated protein kinasesOsteoarthritis Cartilage1815091517201010.1016/j.joca.2010.08.014208131943013628 TetsunagaTNishidaKFurumatsuTNaruseKHirohataSYoshidaASaitoTOzakiTRegulation of mechanical stress-induced MMP-13 and ADAMTS-5 expression by RUNX-2 transcriptional factor in SW1353 chondrocyte-like cellsOsteoarthritis Cartilage19222232201110.1016/j.joca.2010.11.00421094261 CailottoFBianchiASebillaudSVenkatesanNMoulinDJouzeauJYNetterPInorganic pyrophosphate generation by transforming growth factor-beta-1 is mainly dependent on ANK induction by Ras/Raf-1/extracellular signal-regulated kinase pathways in chondrocytesArthritis Res Ther9R122200710.1186/ar2330180348742246241

Morphology of CEP chondrocytes with ICMT stimulation. Magnification, ×200. CEP chondrocytes without ICMT stimulation had a polygonal morphology with a random distribution. CEP chondrocytes exposed to ICMT had a spindle morphology, additionally they had orderly alignment with a certain direction. CEP, cartilage endplate; ICMT, intermittent cyclic tension.

Validation of the representative microRNAs by reverse transcription-quantitative polymerase chain reaction. (A) Upregulated and (B) downregulated microRNAs. *P<0.05 vs. control. ICMT, intermittent cyclic tension; miR, microRNA.

Kyoto Encyclopedia of Genes and Genomes pathway analysis based on the genes targeted by the significantly differentially expressed microRNAs. Top 20 signaling pathways enriched by the target genes of the significantly differentially expressed microRNAs. MAPK, mitogen-activated protein kinase.

Gene Ontology analysis based on the genes targeted by the significantly differentially expressed microRNAs. Top 20 biological processes enriched by the target genes of the significantly differentially expressed microRNAs.

Patient information and the usage of cartilage endplate specimens.

Case no. Age Gender Level Modic type Validation type
1 43 Female L4/L5 0 Microarray
2 45 Male L3/L4 0 Microarray
3 49 Male L5/S1 0 Microarray
4 47 Male L5/S1 0 PCR
5 44 Female L5/S1 0 PCR
6 48 Male L5/S1 0 PCR

Forward primers used in reverse transcription-quantitative polymerase chain reaction analysis.

microRNA Forward
hsa-miR-548q GGCCGCTGGTGCAAAAGTAA
hsa-miR-637 ACTGGGGGCTTTCGGGC
hsa-miR-614 GAACGCCTGTTCTTGCCAG
hsv1-miR-H14 AGTCGCACTCGTCCCTGG
hsa-miR-581 GGTCTTGTGTTCTCTAGATCAG
hsa-miR-142-3p GTGTAGTGTTTCCTACTTTATGGA
hsa-let-7a-5p GCTGAGGTAGTAGGTTGTATAG
hsa-miR-181a-5p AACATTCAACGCTGTCGGTGA
hsa-miR-19b-3p GGCTGTGCAAATCCATGCAAA
hsa-miR-29c-3p GCTAGCACCATTTGAAATCGGT

Significantly upregulated microRNAs in cartilage endplate chondrocytes under cyclic mechanical tension.

ID Name Fold-change P-value
42,865 hsa-miR-181a-5p 15.84 3.29×10−2
10,947 hsa-miR-142-3p 14.65 2.29×10−2
11,041 hsa-miR-29c-3p 10.53 2.00×10−2
10,998 hsa-miR-19b-3p 10.07 3.11×10−2
11,013 hsa-miR-181a-3p   9.22 4.65×10−3
4610 hsa-miR-126-3p   8.47 1.94×10−2
9938 hsa-let-7i-5p   6.00 3.82×10−2
145,693 hsa-miR-92a-3p   5.19 4.85×10−2
10,990 hsa-miR-196a-5p   5.04 8.43×10−3
17,506 hsa-miR-24-3p   4.90 1.11×10−2
32,884 hsa-miR-342-3p   4.76 2.17×10−2
147,162 hsa-let-7a-5p   4.44 3.52×10−2
19,588 hsa-miR-17-3p   3.63 2.25×10−2
146,049 hsa-miR-28-5p   3.15 2.17×10−2
33,902 hsa-miR-128-3p   2.88 1.30×10−2
10,925 hsa-miR-10b-5p   2.64 3.42×10−2
17,935 hsa-miR-101-5p   2.52 2.21×10−2
42,783 hsa-miR-197-3p   2.47 9.19×10−2
42,829 hsa-miR-127-3p   2.25 1.96×10−2
168,925 hsa-miR-1273g-3p   2.18 1.72×10−2
169,385 hsa-miR-4500   2.02 3.12×10−2

Significantly downregulated microRNAs in cartilage endplate chondrocytes under cyclic mechanical tension.

ID Name Fold-change P-value
145,750 hsa-miR-614 −15.33 1.40×10−2
147,815 hsv1-miR-H14-5p   −7.84 1.47×10−2
14,962 hsa-miR-581   −5.81 5.64×10−2
146,148 hsa-miR-548q   −5.60 4.05×10−2
169,118 hsa-miR-5009-3p   −5.55 4.19×10−2
148,677 hsa-miR-637   −5.29 1.54×10−2
147,623 hsa-miR-4304   −5.17 1.33×10−2
168,895 hsa-miR-548ag   −4.40 1.58×10−2
46,606 hsa-miR-1288-3p   −4.15 3.34×10−2
17,402 ebv-miR-BART2-5p   −4.12 3.10×10−2
168,667 hsa-miR-4999-3p   −3.87 1.65×10−2
45,775 hsa-miR-1279   −3.84 2.98×10−2
148,140 hsa-miR-181d-3p   −3.84 1.75×10−2
169,131 hsa-miR-4724-3p   −3.82 3.26×10−2
46,235 hsa-miR-524-5p   −3.79 8.34×10−3
42,493 hsa-miR-892b   −3.58 1.87×10−2
169,281 hsa-miR-4752   −3.27 4.72×10−2
169,084 hsa-miR-5708   −3.17 3.51×10−3
147,866 hsa-miR-3134   −3.12 4.62×10−2
42,786 hsa-miR-188-3p   −3.11 2.14×10−2
11,154 hsa-miR-517c-3p   −3.07 4.48×10−2
148,474 hsa-miR-3622a-5p   −3.00 2.71×10−2
169,318 hsa-miR-4649-3p   −2.95 3.40×10−2
169,164 hsa-miR-4705   −2.92 2.64×10−2
42,773 ebv-miR-BART17-3p   −2.90 4.87×10−2
17,851 hsa-miR-200c-5p   −2.88 3.22×10−2
168,698 hsa-miR-3127-3p   −2.73 4.35×10−2
168,886 hsa-miR-4745-5p   −2.71 2.72×10−2
148,215 hsa-miR-3591-3p   −2.71 3.65×10−2
147,979 hsa-miR-3150a-3p   −2.68 1.33×10−2
46,634 hsa-miR-1281   −2.65 4.57×10−2
148,102 hsa-miR-3619-5p   −2.60 4.67×10−2
147,977 hsa-miR-3190-3p   −2.48 1.02×10−2
46,800 hsa-miR-1224-3p   −2.40 3.47×10−4
168,745 hsa-miR-4667-3p   −2.40 3.02×10−2
146,006 hsa-miR-670-5p   −2.40 3.03×10−2
169,080 hsa-miR-4684-5p   −2.36 2.83×10−2
146,042 hsv1-miR-H8-3p   −2.33 4.63×10−2
42,613 ebv-miR-BART19-5p   −2.33 3.63×10−2
46,363 hsa-miR-1272   −2.32 2.82×10−2
46,556 hsa-miR-623   −2.28 2.30×10−2
46,705 hsa-miR-548k   −2.28 5.75×10−3
17,336 hsa-miR-618   −2.28 2.53×10−3
146,169 mcv-miR-M1-3p   −2.22 4.99×10−2
46,210 hsa-miR-1249-3p   −2.20 2.25×10−2
21,498 hsa-miR-654-3p   −2.19 3.65×10−2
169,166 hsa-miR-4683   −2.17 1.80×10−2
17,306 ebv-miR-BART12   −2.15 2.81×10−2
42,584 hsa-miR-432-3p   −2.15 2.51×10−3
168,904 hsa-miR-4473   −2.12 5.84×10−3
17,299 hcmv-miR-UL22A-3p   −2.09 4.27×10−2
148,063 hsa-miR-3713   −2.09 2.82×10−2
169,038 hsa-miR-488-3p   −2.09 3.31×10−3
168,603 hsa-miR-4664-5p   −2.07 2.58×10−3
148,398 hsa-miR-3908   −2.07 1.93×10−3
147,837 hsa-miR-3119   −2.07 4.05×10−2
168,850 hsa-miR-3191-5p   −2.07 2.24×10−2
169,185 hsa-miR-5187-3p   −2.07 2.71×10−3
46,408 hsa-miR-1322   −2.06 1.04×10−2
42,457 hsa-miR-323a-5p   −2.05 1.68×10−2
168,841 hsa-miR-5588-3p   −2.04 9.37×10−3
148,643 hsa-miR-642a-5p   −2.04 3.38×10−2