- Research article
- Open access
- Published:
IL17RC affects the predisposition to thoracic ossification of the posterior longitudinal ligament
Journal of Orthopaedic Surgery and Research volumeĀ 14, ArticleĀ number:Ā 210 (2019)
Abstract
Background
Thoracic ossification of the posterior longitudinal ligament (T-OPLL) can cause thoracic spinal stenosis, which results in intractable myelopathy and radiculopathy. The etiology of T-OPLL is unknown and the condition is difficult to treat surgically. Whole-genome sequencing identified a genetic variant at rs199772854 of the interleukin 17 receptor C (IL17RC) gene as a potentially pathogenic locus associated with T-OPLL. We aimed to determine whether the rs199772854A site mutation causes abnormal expression of the IL17RC in Han Chinese patients with T-OPLL and predict the possible pathogenic mechanisms of T-OPLL. Analyses were performed to determine whether IL17RC is involved in the pathogenicity of T-OPLL.
Methods
Peripheral blood and OPLL tissue were collected from a total of 72 patients with T-OPLL disease (36 patients carrying the rs199772854A site mutation in IL17RC and 36 wild-type patients). The expression of IL17RC was analyzed by enzyme-linked immunosorbent assay, reverse transcription-quantitative polymerase chain reaction, immunohistochemistry, and Western blotting.
Results
rs199772854A mutation resulted in markedly increased IL17RC gene expression levels in peripheral blood samples and the OPLL tissue obtained following clinical surgery (Pā<ā0.05).
Conclusions
The results suggest that the rs199772854A site mutation of IL17RC can significantly increase the expression of IL17RC. The IL17RC gene rs199772854A site polymorphism is a potential pathogenic mutation in T-OPLL disease, which may be associated with the occurrence of T-OPLL.
Background
Thoracic ossification of the posterior longitudinal ligament (T-OPLL) is characterized by pathological heterotopic ossification of this region. T-OPLL is one of the common factors that cause thoracic spinal stenosis, which results in intractable myelopathy and radiculopathy. Early diagnosis of the disease is difficult, as the majority of patients develop symptoms only when the ossified ligament severely compresses the spinal cord. T-OPLL causes a much higher rate of disability than cervical OPLL (C-OPLL). T-OPLL has a marked ethnic predilection, as this disease occurs more frequently in Japanese and Chinese individuals. The prevalence of T-OPLL in individuals of Japanese ethnicity is 1.6ā1.9%, and the mean age of onset is >ā40Ā years old [1, 2].
Conservative treatments are typically ineffective for T-OPLL, with surgery as the only effective treatment. However, due to the unique anatomy and pathophysiological factors associated with for T-OPLL, postoperative cerebrospinal fluid leakage, paraplegia, infection, and other complications are common, occurring in 9.6ā40.8% of patients that receive surgical treatment [3,4,5], and the pathogenesis of T-OPLL remains unclear. Most studies suggest that C-OPLL is a āgeneticā disease [6,7,8,9,10,11,12,13,14,15,16], the thoracic spine experiences less local biological stress than the cervical spine; thus, we speculate that genetic factors may contribute to the development of T-OPLL. Our previous whole-genome sequencing and candidate gene-association studies demonstrated that the presence of the rs199772854A single-nucleotide polymorphism (SNP) in the interleukin 17 receptor C (IL17RC) gene is potentially associated with T-OPLL susceptibility [17, 18]. Therefore, we hypothesize that IL17RC might be involved in the formation of OPLL of the thoracic spine.
The IL17RC gene is located at chromosome region 3p25.3-3p24.1 and encodes a type I transmembrane protein [19]. IL17RC is a proinflammatory receptor with a crucial role in the development of osteoblasts and accelerates osteoblast differentiation [20]. As OPLL promotes bone formation in ligament tissue, OPLL is associated with increased bone mineral density.
The present study aimed to determine whether the rs199772854A site mutation causes abnormal expression of the IL17RC gene in patients with T-OPLL among a Han Chinese population and to determine whether IL17RC is involved in the pathogenicity of T-OPLL.
Materials and methods
Disease criteria and patients
This retrospective study protocol was approved by the ethical committee for human subjects of the Peking University Third Hospital (Beijing, China). Informed consent was provided by all participating individuals. Unrelated northern Chinese Han patients with T-OPLL carrying the rs199772854A site mutation in IL17RC and unrelated northern Chinese Han patients with T-OPLL carrying the wild-type rs199772854C site were enrolled in this study between May 2014 and July 2018. Diagnosis of T-OPLL was performed by orthopedic spine specialists based on clinical symptoms and computed tomography (CT) of the thoracic spine. The appearance of T-OPLL observed in CT was classified as segmental, continuous, mixed, or local disease type [21]. Neurological status was evaluated by the Japanese Orthopedic Association (JOA) score for thoracic myelopathy (maximum 11 points). The posterior longitudinal ligament specimens of the thoracic spine in patients with T-OPLL were collected during circumferential decompression surgery (Fig.Ā 1).
Plasma IL17RC enzyme-linked immunosorbent assay
Plasma collection and storage were performed using standard methods. Plasma IL17RC levels were quantified using commercially available enzyme-linked immunosorbent assay (ELISA) kits (Trust Specialty Zeal, Inc., San Francisco, CA, USA). All samples were assayed according to the manufacturerās instructions and were run in duplicate. The optical density of each well was determined using a microplate reader at 450Ā nm. No interference and no cross reactivity were expected based on the manufacturerās instructions.
Reverse transcription-quantitative polymerase chain reaction
Total RNA was purified from blood using the SK1321 RNA Blood Mini Kit (Sangon Biotech Co., Ltd., Shanghai, China). A one-column DNase digest (Sangon Biotech Co., Ltd.) was performed before the clean-up step to eliminate residual genomic DNA. cDNA was synthesized from total RNA (2Ā Ī¼g) using a RevertAid Premium Reverse Transcriptase kit (Thermo Fisher Scientific, Inc., Waltham, MA, USA). Relative qPCR was applied to quantify the mRNAs levels of IL17RC using SYBR Green Real-Time PCR master mix on the LightCycler480 Real-Time System (Roche Diagnostics, Basel, Switzerland). All experiments were performed in triplicate and normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Details of the primer sequences are listed in TableĀ 1.
Hematoxylin-eosin staining and immunohistochemistry analysis
Serial 5-mm-thick sections were prepared from paraffin-embedded thoracic spine specimens for staining. Hematoxylin-eosin staining was performed in an autostainer machine (Leica Microsystems GmbH, Mannheim, Germany) using standard procedures. Sections for immunohistochemical (IHC) staining were deparaffinized using xylene and dehydrated in serially graded ethanol solutions. The sections were washed in distilled water, treated with a 0.3% H2O2 solution dissolved in absolute methanol at 20Ā Ā°C for 15Ā min, and then rinsed with PBS (pH 7.4). Antigen retrieval was performed using a high temperature and high pressure method. The sections were incubated with primary polyclonal mouse anti-human IL17RC antibody (1:200; ab69673; Abcam, Cambridge, MA, USA) at 4Ā Ā°C overnight in a humidified chamber. Sections were washed with PBS three times for 5Ā min each wash and then incubated with horseradish peroxidase-conjugated goat anti-mouse IgG in a humidified chamber for 30Ā min at room temperature. Sections were rinsed with PBS (pH 7.4), and antibody binding was visualized by incubation with a diaminobenzidine (DAB) solution (ZLI-9017; OriGene Technologies, Inc., Beijing, China). Sections were washed in water to remove excess DAB and counterstained with hematoxylin to visualize nuclei. Negative control sections were incubated with PBS instead of the primary antibodies under the same conditions.
Western blot analysis
Tissue lysates were obtained using ice-cold RIPA lysis buffer (Beyotime Institute of Biotechnology, Haimen, China) containing 100Ā mM PMSF as a protease inhibitor. Total protein (100Ā Ī¼g) was separated in a Bis-Tris polyacrylamide gel and transferred onto a nitrocellulose membrane. The membrane was then incubated in 1% bovine serum albumin containing primary rabbit anti-human polyclonal antibodies at 4Ā Ā°C overnight. Following incubation with horseradish peroxidase-conjugated goat anti-rabbit antibody at room for 1Ā h, proteins were detected using electrochemiluminescence (EMD Millipore, Billerica, MA, USA). The following primary and secondary antibodies were used: anti-IL17RC (1:1,000; ab69673; Abcam) and goat anti-mouse antibody (1:2,500; CW0102M; Beijing Kangwei Century Biotechnology Co., Ltd., Beijing, China). The blots were detected using a Kodak film developer (Fujifilm, Tokyo, Japan). Protein levels were quantified by densitometry analysis using Image-Pro Plus 6.0 software (Media Cybernetics, Inc., Rockville, MD, USA). Beta-actin was used as the endogenous control.
Statistical analysis
All statistical analyses were performed using SPSS v17.0 software (SPSS, Inc., Chicago, IL, USA). Descriptive data for continuous variables are presented as the meanāĀ±āstandard deviations. Studentās t test was used to compare age, gender, and JOA score differences between T-OPLL patients with different IL17RC gene variants. The differences in T-OPLL subtypes between patients with or without IL17RC gene mutation were applied using one-way analysis of variance with post hoc Fisherās test. Pāā¤ā0.05 was considered to be statistically significant.
Results
Genotype-phenotype analysis
A total of 36 unrelated northern Chinese Han patients with T-OPLL carrying the rs199772854A site mutation in IL17RC (16 men, mean age 56.11āĀ±ā11.08Ā years; 20 women, mean age 53.25āĀ±ā6.65Ā years) and 36 unrelated northern Chinese Han patients with T-OPLL carrying the wild-type rs199772854C site (14 men, mean age 58.43āĀ±ā8.89Ā years; 22 women, mean age 53.45āĀ±ā6.51Ā years) were enrolled in this study. Phenotype-genotype associations were analyzed among the T-OPLL patients with or without the rs199772854A mutation (nā=ā36 per group; TableĀ 2). No differences were found between these two groups in terms of sex, age, and JOA score at diagnosis. Additionally, radiological analysis of T-OPLL morphology revealed that the mutation-positive patients and mutation-negative patients exhibited no difference in the disease type classification (segmental, continuous, mixed, or local).
Analysis of IL17RC levels in the blood of patients with T-OPLL
The plasma concentration of IL17RC was shown in Fig.Ā 2; plasma IL17RC concentration was significantly higher (~ā3-fold higher) in T-OPLL patients with rs199772854A mutation (9.91āĀ±ā1.91Ā Ī¼g/l) compared with T-OPLL patients carrying the wild-type rs199772854C variant (3.34āĀ±ā0.69āĪ¼g/l; PĀ <ā0.01). Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis was performed using RNA extracted from peripheral blood cells (Fig.Ā 3) and demonstrated that IL17RC mRNA levels were ~ā6-fold higher in T-OPLL patients carrying the rs199772854A mutation than T-OPLL patients with the wild-type rs199772854C site (PĀ <ā0.001). Compared with wild-type T-OPLL patients, the rs199772854A mutation significantly increased IL17RC gene expression, suggesting that this is a potential pathogenic locus that alters IL17RC gene expression in cells.
IL17RC protein expression in T-OPLL patients with rs199772854A mutation
Hematoxylin-eosin staining, shown in Fig.Ā 4, revealed that there were more chondrocytes in the ossification zone (open arrows) and zigzag tidal traces formed by calcification were visible near the chondrocytes (solid arrows) in T-OPLL patients. IHC revealed that IL17RC protein was positively expressed in the ossified areas of T-OPLL patients with rs199772854A mutation, but no nuclear reactivity was observed in the normal fibers. However, there was weakly positive IL17RC in T-OPLL patients carrying the wild-type rs199772854C variant (Fig.Ā 5). Western blot (WB) analysis revealed that the expression of IL17RC protein was significantly higher in T-OPLL patients with the IL17RC gene rs199772854A mutation than T-OPLL patients carrying the wild-type rs199772854C variant (Fig.Ā 6).
Discussion
T-OPLL poses a major challenge for spinal surgeons. Although the exact cause of the disease is still unclear, it is generally believed that genetic factors have an important role in the development of the disease. Some scholars believe that the accumulation of harmful missense mutations in the human genome creates the genetic basis for various complex diseases [22]. Whole genome sequencing identified a Cā>āA mutation at the rs199772854 locus of the IL17RC gene. Disease progression is also affected by the gene expression in peripheral blood cells, with several genes reported to exhibit higher expression in the peripheral blood of OPLL patients compared to healthy controls [9, 23].
In this study, peripheral blood of T-OPLL patients with or without an rs199772854 mutation was collected and analyzed to assess the role of the mutated gene locus. The results of ELISA and qPCR analysis demonstrated that the expression of IL17RC was significantly higher in the peripheral blood of T-OPLL patients carrying the rs199772854A mutation compared to patients without the mutation. IHC and WB were performed to determine IL17RC expression in patients with T-OPLL, and the results demonstrated that the expression of IL17RC protein was significantly higher in the T-OPLL patients with the mutation compared to those with the wild-type. Therefore, the rs199772854A site mutation may lead to overexpression of IL17RC. We also showed that there was no difference in the OPLL classification and JOA score between T-OPLL patients with and without the rs199772854A mutation. The possible association between the rs199772854A site mutation in IL17RC and the severity of the T-OPLL phenotypes requires larger-scale studies in the future.
The rs199772854 site of the IL17RC gene is located in the promoter region. An SNP in the promoter region can increase or reduce gene transcriptional activity by altering the binding efficiency of transcription factors to various sequence elements, thus interfering with the gene expression process and potentially leading to disease occurrence. However, the mechanism and role of IL17RC gene regulation during ectopic osteogenesis remains unclear. IL17RC may accelerate bone metabolism via the transforming growth factor-Ī² signaling pathway [24, 25]. The majority of studies suggest that IL17RC has a major role in disease pathogenesis through its function in the IL17 signaling axis. Some studies have shown that the IL17 axis affects bone formation and remodeling, and can also protect bone mass if bone loss occurs due to infection or hormonal imbalance [26, 27]. In addition, the IL17 axis can induce osteoblastic differentiation of bone marrow-derived mesenchymal stem cells [28]. It has become increasingly clear that inflammation-mediated imbalance in bone is a major feature of various bone diseases [29]; this imbalance is caused by increases in various cytokines in the inflammatory tissue. As an inflammatory factor, IL17 RC are expressed on various cells, such as osteoblasts, chondrocytes, and fibroblasts [30], and compression force can induces the expression of IL-17RC in osteoblast-like cells [31, 32]. In addition, the differentiation of osteoclast precursors into osteoclasts is suppressed by IL17RC [30]; IL17RC produced in response to compressive force may suppress osteoclastogenesis through the expression of OPG [31]. Additionally, IL17RC stimulates the secretion of other factors such as IL-6, tumor necrosis factor (TNF)-Ī± and IL-1Ī² in osteoclasts, further aggravating the inflammation [33].
The results of this study demonstrated that the expression of IL17RC in T-OPLL patients carrying an rs199772854A mutation was significantly higher in peripheral blood and tissues than in patients without the mutation. It is suggested that the rs199772854A site mutation can lead to overexpression of IL17RC and may induce pathological OPLL. However, this study lacks in-depth research and discussion on the mechanism by which IL17RC facilitates T-OPLL. In addition, due to the prevalence of T-OPLL disease is very rare, this studyās sample size is small. The association between heterotopic ossification and IL17RC requires further research.
Conclusions
In conclusion, the findings of this study suggest that the rs199772854A site mutation can lead to overexpression of IL17RC and is a potential pathogenic mutation associated with T-OPLL. The results provide a potential basis for the pathogenesis of T-OPLL and the pathogenic role of IL17RC in T-OPLL disease.
Further genetic studies with more participants are required to validate these findings, and in appropriate model systems.
Abbreviations
- C-OPLL:
-
Cervical ossification of the posterior longitudinal ligament
- CT:
-
Computed tomography
- ELISA:
-
Enzyme-linked immunosorbent assay
- IHC:
-
Immunohistochemical
- IL17RC:
-
Interleukin 17 receptor C
- JOA:
-
Japanese Orthopedic Association
- RT-qPCR:
-
Reverse transcription-quantitative polymerase chain reaction
- T-OPLL:
-
Thoracic ossification of the posterior longitudinal ligament
- WB:
-
Western blot
References
Fujimori T, Watabe T, Iwamoto Y, Hamada S, Iwasaki M, Oda T. Prevalence, concomitance, and distribution of ossification of the spinal ligaments. Spine. 2016;41:1668ā76.
Mori K, Imai S, Kasahara T, Nishizawa K, Mimura T, Matsusue Y. Prevalence, distribution, and morphology of thoracic ossification of the posterior longitudinal ligament in Japanese: results of CT-based cross-sectional study. Spine (Phila Pa 1976). 2014;39(5):394ā9.
Yu LJ, Li WJ, Guo SG, Zhao Y. Transforaminal thoracic interbody fusion: treatment of thoracic myelopathy caused by anterior compression. Orthopade. 2018;47(12):986ā92.
Imagama S, Ando K, Takeuchi K, et al. Perioperative complications after surgery for thoracic ossification of posterior longitudinal ligament- nationwide multicenter prospective study. Spine (Phila Pa 1976). 2018;43(23):E1389āE1397.
Hu P, Yu M, Liu X, Liu Z, Jiang L. A circumferential decompression-based surgical strategy for multilevel ossification of thoracic posterior longitudinal ligament. Spine J. 2015;15:2484ā92.
Xu C, Chen Y, Zhang H, et al. Integrated microRNA-mRNA analyses reveal OPLL specific microRNA regulatory network using high-throughput sequencing. Sci Rep. 2016;6:21580.
Yang H, Shi L, Shi G, et al. Connexin 43 affects osteogenic differentiation of the posterior longitudinal ligament cells via regulation of ERK activity by stabilizing Runx2 in ossification. Cell Physiol Biochem. 2016;38:237ā47.
Nakajima M, Kou I, Ohashi H, Ikegawa S. Identification and functional characterization of RSPO2 as a susceptibility gene for ossification of the posterior longitudinal ligament of the spine. Am J Hum Genet. 2016;99:202ā7.
Chen X, Guo J, Cai T, et al. Targeted next-generation sequencing reveals multiple deleterious variants in OPLL-associated genes. Sci Rep. 2016;6:26962.
Ikegawa S. Genetics of ossification of the posterior longitudinal ligament of the spine: a mini review. J Bone Meta. 2014;21:127ā32.
Guo Q, Lv S, Wu S, Tian X, Li Z. Association between single nucleotide polymorphism of IL15RA gene with susceptibility to ossification of the posterior longitudinal ligament of the spine. J Orthop Surg Res. 2014;9:103ā7.
Wei W, He HL, Chen CY, et al. Whole exome sequencing implicates PTCH1 and COL17A1 genes in ossification of the posterior longitudinal ligament of the cervical spine in Chinese patients. Genet Mol Res. 2014;13:1794ā804.
Ikegawa S. Genomic study of ossification of the posterior longitudinal ligament of the spine. Proc Jpn Acad Ser B. 2014;90:405ā12.
Nakajima M, Takahashi A, Tsuji T, et al. A genome-wide association study identifies susceptibility loci for ossification of the posterior longitudinal ligament of the spine. Nat Genet. 2014;46:1012ā6.
Ren Y, Feng J, Liu ZZ, Wan H, Li JH, Lin X. A new haplotype in BMP4 implicated in ossification of the posterior longitudinal ligament (OPLL) in a Chinese population. J Orthop Res. 2012;30(5):748ā56.
He Z, Zhu H, Ding L, Xiao H, Chen D, Xue F. Association of NPP1 polymorphism with postoperative progression of ossification of the posterior longitudinal ligament in Chinese patients. Genet Mol Res. 2013;12:4648ā55.
Wang P, Liu X, Zhu B, et al. Identification of susceptibility loci for thoracic ossification of the posterior longitudinal ligament by whole-genome sequencing. Mol Med Rep. 2018;17:2557ā64.
Wang P, Liu X, Zhu B, et al. Association of IL17RC and COL6A1 genetic polymorphisms with susceptibility to ossification of the thoracic posterior longitudinal ligament in Chinese patients. J Orthop Surg Res. 2018;13:109.
Ho AW, Gaffen SL. IL-17RC: a partner in IL-17 signaling and beyond. Semin Immunopathol. 2010;32:33ā42.
Sato K, Suematsu A, Okamoto K, et al. Th17 functions as an osteoclastogenic helper T cell subset that links T cell activation and bone destruction. J Exp Med. 2006;203:2673ā82.
Tsuyama N. Ossification of the posterior longitudinal ligament of the spine. Clin Orthop Relat Res. 1984;184:71ā84.
Kryukov GV, Pennacchio LA, Sunyaev SR. Most rare missense alleles are deleterious in humans: implications for complex disease and association studies. Am J Hum Genet. 2007;80:727ā39.
Niu C, Lin S, Yuan L, et al. Correlation of blood bone turnover biomarkers and Wnt signaling antagonists with AS, DISH, OPLL, and OYL. BMC Musculoskel Dis. 2017;18:s12817ā91.
Wang Y, Xing F, Ye S, et al. Jagged-1 signaling suppresses the IL-6 and TGF-beta treatment-induced Th17 cell differentiation via the reduction of RORgammat/IL-17A/IL-17F/IL-23a/IL-12rb1. Sci Rep. 2015;5:8234.
Mohammad KS, Chen CG, Balooch G, et al. Pharmacologic inhibition of the TGF-beta type I receptor kinase has anabolic and anti-catabolic effects on bone. PLoS One. 2009;4:e5275.
Lin D, Li L, Sun Y, Wang W, et al. IL-17 regulates the expressions of RANKL and OPG in human periodontal ligament cells via TRAF6/TBK1-JNK/NF-kappaB pathways. Immunology. 2014;144:472ā85.
Kotake S, Udagawa N, Takahashi N, et al. IL-17 in synovial fluids from patients with rheumatoid arthritis is a potent stimulator of osteoclastogenesis. J Clin Invest. 1999;103:1345ā52.
Huang H, Kim HJ, Chang EJ, et al. IL-17 stimulates the proliferation and differentiation of human mesenchymal stem cells: implications for bone remodeling. Cell Death Differ. 2009;16:1332ā43.
Croes M, ƶner FC, van Neerven D, et al. Proinflammatory T cells and IL-17 stimulate osteoblast differentiation. BONE. 2016;84:262ā70.
Kitami S, Tanaka H, Kawato T, et al. IL-17A suppresses the expression of bone resorption-related proteinases and osteoclast differentiation via IL-17RA or IL-17RC receptors in RAW264.7 cells. Biochimie. 2010;92:398ā404.
Zhang F, Wang CL, Koyama Y, et al. Compressive force stimulates the gene expression of IL-17s and their receptors in MC3T3-E1 cells. Connect Tissue Res. 2010;51:359ā69.
Rossini M, Viapiana O, Adami S, et al. Focal bone involvement in inflammatory arthritis: the role of IL17. Rheumato Int. 2016;36:469ā82.
Du Y, Tong Y, Mei W, et al. A truncated IL-17RC peptide ameliorates synovitis and bone destruction of arthritic mice. Adv Healthc Mater. 2016;5:2911ā21.
Acknowledgements
Not applicable.
Funding
This work was supported by the National Natural Science Foundation of China (81672201 and 81472041).
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Author information
Authors and Affiliations
Contributions
PW, XL, CK, and SL conceived of and designed the experiments. PW, YM, LY, and CL performed the experiments. ZT, LY, XL, and GH analyzed the data. PW wrote the paper. SL. revised the manuscript. All authors reviewed the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study protocol was approved by the Ethics Committee for Human Subjects of the Peking University Third Hospital. Informed consent was provided by all participating individuals.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisherās Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Wang, P., Liu, X., Liu, X. et al. IL17RC affects the predisposition to thoracic ossification of the posterior longitudinal ligament. J Orthop Surg Res 14, 210 (2019). https://doi.org/10.1186/s13018-019-1253-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13018-019-1253-3