- Research article
- Open Access
Bioinformatic analysis and experimental identification of blood biomarkers for chronic nonunion
Journal of Orthopaedic Surgery and Research volume 15, Article number: 208 (2020)
Incomplete fracture healing may lead to chronic nonunion; thus, determining fracture healing is the primary issue in the clinical treatment. However, there are no validated early diagnostic biomarkers for assessing chronic nonunion. In this study, bioinformatics analysis combined with an experimental verification strategy was used to identify blood biomarkers for chronic nonunion.
First, differentially expressed genes in chronic nonunion were identified by microarray data analysis. Second, Dipsaci Radix (DR), a traditional Chinese medicine for fracture treatment, was used to screen the drug target genes. Third, the drug-disease network was determined, and biomarker genes were obtained. Finally, the potential blood biomarkers were verified by ELISA and qPCR methods.
Fifty-five patients with open long bone fractures (39 healed and 16 nonunion) were enrolled in this study, and urgent surgical debridement and the severity of soft tissue injury had a significant effect on the prognosis of fracture. After the systems pharmacology analysis, six genes, including QPCT, CA1, LDHB, MMP9, UGCG, and HCAR2, were chosen for experimental validation. We found that all six genes in peripheral blood mononuclear cells (PBMCs) and serum were differentially expressed after injury, and five genes (QPCT, CA1, MMP9, UGCG, and HCAR2) were significantly lower in nonunion patients. Further, CA1, MMP9, and QPCT were markedly increased after DR treatment.
CA1, MMP9, and QPCT are biomarkers of nonunion patients and DR treatment targets. However, HCAR2 and UGCG are biomarkers of nonunion patients but not DR treatment targets. Therefore, our findings may provide valuable information for nonunion diagnosis and DR treatment.
ISRCTN, ISRCTN13271153. Registered 05 April 2020—Retrospectively registered.
Fracture healing is a complex process and is dependent on multiple factors . Limited fracture healing can occur in 5–10% of fracture cases , resulting in chronic nonunion and functional disability, which can have a devastating impact on the patient’s quality of life . The determination of nonunion is essential for fracture diagnosis and subsequent treatment, but there is a lack of objective tools to assess fracture healing, making nonunion as an uncertain outcome . Therefore, establishing a new diagnosis method is vital for chronic nonunion diagnosis and treatment. Undoubtedly, the circulating biomarkers of bone fracture healing are gaining popularity as possible early predictors of chronic nonunion . However, there are currently no valid biomarkers for chronic nonunion diagnosis in the blood.
In fracture healing, monocytes are involved in angiogenesis  and differentiate into osteoclasts . Many bone formation-related proteins were changed in human peripheral blood mononuclear cells (PBMCs) after bone fracture , which indicated that the altered genes of PBMCs might be early predictors of chronic nonunion. To investigate this possibility, the differentially expressed genes were analyzed using microarray datasets, and potential biomarkers were predicted by a systems pharmacology approach.
Dipsaci Radix (DR) is derived from Dipsacus asperoides, a classical, traditional Chinese medicine with a long history of safe use for the treatment of bone fractures [9, 10]. DR can also be used for the treatment of lower back pain and traumatic hematoma . Moreover, growing evidence indicates that DR could improve osteoporosis by increasing bone density and bone strength [12, 13]. In addition, many studies have confirmed that DR can effectively reestablish the dynamic balance between bone formation and bone resorption . However, no research has examined the potential of DR in the treatment of chronic nonunion.
In this study, bioinformatics analysis methods, including microarray data analysis and an integrated systems pharmacology approach, were used to predict the potential biomarkers of chronic nonunion and mechanisms of DR treatment. First, the differentially expressed chronic nonunion genes in PBMCs were identified using microarray datasets. Second, the potential active ingredients of DR were used to screen the possible target genes. Third, the DR target genes that were also potential biomarkers of chronic nonunion were obtained. Finally, the predicted biomarkers of chronic nonunion were verified by experimental methods, such as enzyme-linked immunosorbent assay (ELISA) and quantitative polymerase chain reaction (qPCR), with blood samples from healed and nonunion patients. Our study provides a more specific and effective way to investigate blood biomarkers for chronic nonunion and provides new insight into the mechanisms of DR in the treatment of chronic nonunion.
Materials and methods
Microarray data and difference analysis
Microarray datasets GSE93138, GSE93213, and GSE93215 were downloaded from the Gene Expression Omnibus (GEO) database  and collected using the GPL6244 platform (Affymetrix Human Gene 1.0 ST Array). These microarray datasets form the study “Biomarker Identification in Fracture Healing.” This study aimed to identify potential peripheral blood biomarkers for normal healers, slow healers, and nonunion patients , which is consistent with the design of our study. Difference analysis was performed by R script using the limma (linear models for microarray analysis) R package, with p < 0.05 and |logFC| > 1 as cutoff values for screening differentially expressed genes (DEGs). The DEGs of each comparison group are shown as volcano plots. The interaction of the DEG sets was obtained by Venn diagram.
DR ingredient profiling and drug target gene identification
The bioactive ingredients of DR were collected from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database . We set drug-likeness (DL) ≥ 0.18 and oral bioavailability (OB) ≥ 30% as the threshold for bioactive ingredients. The structures of the ingredients were downloaded from the PubChem database . Afterward, the target genes corresponding to the compounds were screened from the PharmMapper database  and the Swiss Target Prediction database . The details of the data are shown in a previous study .
The STRING database  was used to analyze the protein-protein interactions (PPIs) of the DEGs of healed and nonunion patients, and the hub genes were counted by R script. The interactions of the DEGs and DR target genes were obtained by Venn diagram. Based on the potential target genes, the Cytoscape software (version: 3.7.2) was used to construct the drug-disease network. The GeneCards database was used to analyze the function of the potential biomarkers.
Gene ontology enrichment analysis
The GO enrichment analysis was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) , which is an integrated online biological knowledge base and analytical tool. In our study, the target genes were mapped into DAVID to identify the biological processes, molecular function, and cellular components of the predicted target genes involved.
Patients and treatment
Patients admitted to the Department of Orthopedics of Xuanwu Hospital from August 2018 to July 2019 were enrolled in this study. Participants were skeletally mature, and all were diagnosed with only one open fracture of a long bone, including the humerus, radius/ulna, femur, and tibia/fibula. Exclusion criteria included patients who had a severe head injury, renal insufficiency, liver disease, systemic inflammation (CRP > 0.5 mg/dL), and osteoporosis. Patients whose initial surgical debridement occurred within 8 h and those whose surgical debridement occurred after 8 h were treated with internal fixation according to the type of fracture by well-trained orthopedic surgeons. Ultimately, 55 participants were enrolled in the study and set as the acute injury (AI) group. Then, the AI group was categorized into healed and nonunion groups according to the amendment from the Food and Drug Administration (FDA) on the diagnostic criteria for nonunion. The FDA defines nonunion as a fractured bone that has not completely healed within 9 months of the initial injury and without signs of healing for at least 3 months . In this study, we chose 9 months as the time point to define nonunion. A total of 39 healed patients and 16 nonunion patients were selected in this study, and the clinical characteristics of the patients are shown in Table 1. Moreover, 30 healthy participants were included as the healthy control (HC) group, and no participant in the HC group experienced any fracture. Some of the bone fracture patients were taking DR (21 of the healed patients and 7 of the nonunion patients), and all the patients were followed up every 3 months until 12 months. All procedures performed in the study were approved by the Ethics Committee of Xuanwu Hospital, Capital Medical University (Beijing, China), and informed consent was obtained from all individuals involved in this study or their guardians.
Serum collection and ELISA
Whole blood samples were taken from all participants and centrifuged for 5 min at 4000×g, and serum was isolated and stored at − 80 °C. The serum concentrations of the biomarkers were measured using separate ELISA kits (Mlbio, Shanghai, China) according to the manufacturer’s instructions.
PBMC isolation and qPCR analysis
Five milliliters of peripheral blood was taken from all participants. Ficoll-Paque density gradient centrifugation was used to isolate the PBMCs, and the cells were washed with 4 °C phosphate-buffered saline (PBS). Then, the RNeasy Mini Kit (Qiagen, Germany) was used to extract the total RNA, which was reverse transcribed by the ThermoScript™ RT-PCR System (Invitrogen, Carlsbad, CA, USA). mRNA expression was detected by a Talent qPCR kit (TIANGEN, Beijing, China). The primers used in this study are shown in Table 2.
Statistical analysis was performed by the GraphPad Prism software version 7.0 (CA, USA). All data are displayed as the mean ± SEM. The comparisons between the two groups were analyzed using two-tailed unpaired Student’s t tests or one-way or two-way ANOVA. Categorical data were analyzed using the chi-squared (χ2) test. Values of p < 0.05 were considered significant.
The clinical characteristics of patients with open fractures of long bones
The clinical characteristics of the nonunion group and the healed group are shown in Table 1. There were no significant differences in age, gender, or DR treatment (p > 0.05). Eleven of nonunion patients were smokers, compared to 19 healed patients (p > 0.05). The majority of fractures occurred in healthy patients, and 5 nonunion patients had comorbidities, compared to 9 healed patients (p > 0.05). Most fractures were the result of a traffic accident, and 13 of the nonunion fractures were caused by traffic injury, compared to 27 of the healed fractures (p > 0.05). The larger proportion of injuries resulted from multi-trauma, and 11 of the nonunion fractures were from multi-trauma, compared to 23 of the healed fractures (p > 0.05). The severity of the soft tissue injury was assessed using the Gustilo classification, and 5 nonunion fractures were Gustilo grade III, compared to 3 healed fractures (p < 0.05). The initial surgical debridement times of 6 nonunion and 27 healed fractures were within 8 h, and 10 nonunion and 12 healed fractures received the initial surgical debridement after 8 h (p < 0.05). The results above indicate that urgent surgical debridement and the severity of soft tissue injury have a significant effect on the nonunion of open long bone fractures.
Identification of differentially expressed genes in the PBMCs of healed and nonunion patients
Three PBMC gene expression datasets (GSE93138, GSE93213, and GSE93215) of healed and nonunion patients were downloaded from the GEO database. Gene comparison analysis was performed on these two groups, and there were a total of 258 differentially expressed genes between healed and nonunion patients (Fig. 1a). To explore the relationship between these 258 genes, protein-protein interaction (PPI) analysis was performed using the STRING database. A network was generated, and hub genes were analyzed using the Cytoscape software. The top 20 hub genes in the network were marked in red and yellow (the deeper shades of red indicated greater connections of nodes), and other genes connected to hub genes were presented as blue nodes (Fig. 1b). A bar plot of the number of hub gene links is shown in Fig. 1b, and we found that hub genes, such as RPS27A, RPS17, and RBX1, may play a critical role in the biological activity of chronic nonunion.
Prediction of gene function using enrichment analysis for the differentially expressed genes
To gain a comprehensive understanding of the differentially expressed genes, we used DAVID to perform GO enrichment analysis. A variety of GO terms were enriched, including 43 biological processes, 21 cellular components, and 17 molecular functions. The top 10 GO terms are shown in Fig. 2. The biological processes (such as innate immune response and translational initiation, Fig. 2a), molecular function (such as GTPase activity and protein binding, Fig. 2b), and cellular components (such as extracellular exosome and cytosol, Fig. 2c) may be involved in the biological activity in chronic nonunion.
Prediction of nonunion-related biomarkers by a systems pharmacology approach
In this study, we chose Dipsaci Radix (DR), a classical medicine for the treatment of bone fractures, to predict biomarkers for nonunion patients. Ten potential active compounds were retrieved from the TCMSP database, and a total of 443 target genes were obtained using the PharmMapper and Swiss Target Prediction databases. Finally, the interaction of drug target genes and nonunion-related genes was determined using a Venn diagram. As shown in Fig. 3a, we obtained six interacting genes. Moreover, we constructed an interactive network to connect DR, DR compounds, target genes, and nonunion (Fig. 3b). In our study, we found that five active ingredients of DR (Japonine, Gentisin, Gauloside A, Sylvestroside III, and 3,5-Di-O-caffeoylquinic acid) and six target genes (QPCT, CA1, LDHB, MMP9, UGCG, and HCAR2) were related to nonunion. To better understand the function of the potential biomarkers, the GeneCards database was used to find acute phase symptom-related genes. As shown in Table S1, all six biomarkers were related to the acute phase. Five biomarkers were related to multiorgan failure (CA1, LDHB, MMP9, UGCG, and HCAR2) and inflammation (QPCT, CA1, LDHB, MMP9, and HCAR2), four biomarkers were related to the shock (CA1, LDHB, MMP9, and UGCG), and two were related to acidosis (CA1 and LDHB). Based on the systems pharmacology analysis of these predictions and the results of the network analysis, the six potential biomarkers of nonunion were validated using ELISA and qPCR methods.
Validation of potential biomarkers in serum by ELISA
ELISA was performed to test the six potential biomarkers in serum samples from healthy control (HC) and acute injury (AI) patients. As shown in Fig. 4, the serum concentrations of all six potential biomarkers were increased considerably in the AI group (p < 0.05). The results showed that after bone fracture, the six genes were all significantly increased, which may be closely related to the fracture healing process. However, the serum concentrations of all six potential biomarkers were not significantly different between fracture-only and multi-trauma patients (Fig. S1). To further observe the changes in the six biomarkers at different times, we divided the AI patients into four groups: healed, nonunion, healed+DR, and nonunion+DR, and then detected the potential biomarkers at the 0, 3-, 6-, 9-, and 12-month time points. We found that the serum concentrations of CA1, MMP9, QPCT, HCAR2, and UGCG were considerably decreased in the nonunion groups (p < 0.01, Fig. 5a–e), but there were no significant changes in LDHB (Fig. 5f). DR treatment significantly increased serum concentrations of CA1, MMP9, and QPCT in both the healed and nonunion groups (p < 0.05, Fig. 5a–c), but the HCAR2, UGCG, and LDHB concentrations were not changed (Fig. 5d–f). The six potential biomarkers were all increased after acute injury and gradually decreased after 6 months. The results above indicated that CA1, MMP9, QPCT, HCAR2, and UGCG could be biomarkers to identify nonunion patients, and the treatment of DR may target CA1, MMP9, and QPCT; accelerate fracture healing; and minimize delayed healing and nonunion.
Further validation of potential biomarkers in PBMCs by qPCR analysis
To further validate the CA1, MMP9, and QPCT in nonunion patients, we performed qPCR to identify the changes in the three potential biomarkers in PBMCs from different groups of patients. As shown in Fig. 6, the PBMC mRNA expression levels of CA1, MMP9, and QPCT were significantly changed in the nonunion groups (p < 0.05). DR treatment significantly increased CA1, MMP9, and QPCT, PBMC mRNA expression levels in both the healed and nonunion groups (p < 0.05, Fig. 6a–c), consistent with the ELISA results. The results confirmed that CA1, MMP9, and QPCT could be biomarkers to identify nonunion patients and DR treatment targets.
Fracture healing is a continuous process that greatly depends on the location and type of fracture, the choice of treatment, and other factors related to the host and injury . No standard diagnostic criteria for chronic nonunion has been shown to delay the requirement for treatment, which costs significantly more than those uncomplicated fracture healing . Aside from radiological and clinical examination, serologic markers show promise in predicting the status and quality of fracture healing . The ideal biomarker of bone healing must have the characteristics of accuracy, high sensitivity, specificity, and rapidity, and be inexpensive, and able to predict the outcome of bone nonunion . However, due to the lack of relevant clinical evidence, it is difficult to determine which biomarkers can be reliably used for clinical follow-up and prediction of chronic nonunion .
In our study, we predicted the potential blood biomarkers of nonunion through a series of bioinformatics analyses and combined with subsequent experimental verification. Based on this strategy, we predicted 258 differentially expressed genes in PBMCs between healed and nonunion patients. To further understand the function of these genes, PPI analysis and GO enrichment analysis were conducted. We found that hub genes, such as RPS27A, RPS17, and RBX1, may play a critical role in the biological activity in nonunion patients. The innate immune response, GTPase activity, extracellular exosomes, and cytosol may be involved in the altered biological activity in nonunion patients.
After analysis of these 258 genes by a systems pharmacology approach using the DR drug; six genes, including QPCT, CA1, LDHB, MMP9, UGCG, and HCAR2, were chosen for experimental validation by ELISA and qPCR methods. We verified that all six genes were differentially expressed after acute injury, and five genes (QPCT, CA1, MMP9, UGCG, and HCAR2) were significantly lower in nonunion patients. Meanwhile, CA1, MMP9, and QPCT were markedly increased in patients using DR treatment, suggesting that the proteins might serve as potential blood biomarkers for nonunion and may be potential DR treatment targets. It is worth noting that the six potential biomarkers were related to the acute phase, shock, acidosis, multiorgan failure, and inflammation, but the specific changes of the biomarkers in these manifestations require further study.
Matrix metalloproteases (MMPs) are a family of 23 zinc-dependent proteolytic enzymes that can cleave the extracellular matrix (ECM). MMPs play an essential role in tissue regeneration and bone remodeling processes . Previous studies have shown that MMP9 contributes to bone healing processes, which are essential during fracture repair [30, 31]. Furthermore, MMP9 can be considered an early marker of tissue healing; a high level of MMP9 is indicative of the beginning of the bone remodeling process . In this study, we predicted and verified that the level of MMP9 in serum and PBMCs was significantly lower in the nonunion group than in the healed group. Moreover, DR treatment can promote MMP9 expression, which indicates that MMP9 is not only a biomarker of nonunion patients but also a treatment target of DR.
Carbonic anhydrase I (CA1) is a member of the carbonic anhydrase (CA) family, which catalyzes the reversible hydration and dehydration reactions of CO2/H2CO3 . A study has shown that CA1 stimulates calcium deposition and cell calcification, which is an essential step for new bone formation . Moreover, CA1 recruitment of CO32− ions is also regarded as an essential part of bone fracture healing . In the present study, CA1 was significantly changed in the nonunion group compared with the healed group and can also be increased by DR treatment, which indicated that CA1 was a potential biomarker and treatment target of nonunion patients.
Pituitary glutaminyl cyclotransferase (QPCT), also known as pituitary glutaminyl cyclase, can convert active forms of gonadotropin-releasing hormone (GnRH) peptides to protected forms . A previous study reported that the QPCT gene could affect bone mineral density (BMD) among postmenopausal Japanese women , which was also verified in the Chinese population , and together, these findings indicated that the QPCT gene is one of the osteoporosis susceptibility genes. However, the function of QPCT in bone formation has not been reported. In the present study, QPCT was significantly increased after bone fracture and DR treatment but much lower in the nonunion group than in the healed group, which indicated that QPCT has a function in promoting the remodeling process. Further studies that focus on the roles of QPCT in the diagnosis and treatment of nonunion patients are essential and of great interest.
Hydroxycarboxylic acid receptor 2 (HCAR2) is an endogenous ketone produced by fatty acid oxidation in liver mitochondria during carbohydrate deficiency . HCAR2 has anti-inflammatory and antioxidative properties on immune and epithelial cells , but the roles of HCAR2 in bone fracture healing are still unknown. For the first time, we have demonstrated that HCAR2 significantly increases after bone fracture and has the potential to identify the nonunion but it is not the target of DR treatment.
UDP-glucose ceramide glycosyltransferase (UGCG) is the only enzyme responsible for the de novo production of glucosylceramide (GlcCer), which is essential for proper cell function . LDHB is one isoenzyme of human lactate dehydrogenase (LDH), which is a crucial glycolytic enzyme  and is distributed in different tissues . Neither UGCG nor LDHB is well-known in the bone repair process. This is the first study presenting a relationship between LDHB, UGCG, and bone fracture. We found that both LDHB and UGCG were increased in bone fracture patients, but only UGCG was significantly lower in the nonunion group, with no significant changes during DR treatment. The results above indicated that UGCG might predict nonunion, but neither UGCG nor LDHB was DR targets.
Our study provides valuable information to investigate blood biomarkers for chronic nonunion and the potential mechanisms of DR in the treatment of bone fracture. Despite the large patient cohort, only a small number of patients could be included in the study due to a lack of objective tools to assess nonunion. Although we attempted to increase the accuracy of the results, there was inevitable interference from other factors, and the predicted genes need to be validated on large-scale blood samples further.
This bioinformatics analysis combined with the experimental verification strategy provides five potential blood biomarkers for nonunion patients and three DR treatment targets and is the first study to use such an approach for predicting nonunion blood biomarkers. Early diagnosis of chronic nonunion will help clinicians take timely countermeasures to improve bone healing, which will result in better clinical management of patients. In addition, further prospective clinical studies will evaluate the predictive power of these biomarkers for the prognosis of bone fractures. Changes in the potential biomarkers also need to be further studied in nonunion of other fractures in addition to long bone fractures.
Availability of data and materials
The datasets used and/or analyzed during the current study are available either online or from the corresponding author on reasonable request.
Peripheral blood mononuclear cells
Enzyme-linked immunosorbent assay
Quantitative polymerase chain reaction
Food and Drug Administration (FDA)
Gene Expression Omnibus
Differentially expressed genes
Traditional Chinese Medicine Systems Pharmacology
Database for Annotation, Visualization, and Integrated Discovery
Carbonic anhydrase I
Pituitary glutaminyl cyclotransferase
Bone mineral density
Hydroxycarboxylic acid receptor 2
UDP-glucose ceramide glycosyltransferase
Squyer ER, Dikos GD, Kaehr DM, Maar DC, Crichlow RJ. Early prediction of tibial and femoral fracture healing: are we reliable? Injury. 2016;47(12):2805–8. https://doi.org/10.1016/j.injury.2016.10.036.
Einhorn TA. Enhancement of fracture-healing. J Bone Joint Surg Am. 1995;77(6):940–56. https://doi.org/10.2106/00004623-199506000-00016.
Brinker MR, Hanus BD, Sen M, O'Connor DP. The devastating effects of tibial nonunion on health-related quality of life. J Bone Joint Surg Am. 2013;95(24):2170–6. https://doi.org/10.2106/JBJS.L.00803.
Morshed S, Corrales L, Genant H, Miclau T 3rd. Outcome assessment in clinical trials of fracture-healing. J Bone Joint Surg Am. 2008;90(Suppl 1):62–7. https://doi.org/10.2106/JBJS.G.01556.
Cox G, Einhorn TA, Tzioupis C, Giannoudis PV. Bone-turnover markers in fracture healing. J Bone Joint Surg Br. 2010;92(3):329–34. https://doi.org/10.1302/0301-620X.92B3.22787.
Urbich C, Heeschen C, Aicher A, Dernbach E, Zeiher AM, Dimmeler S. Relevance of monocytic features for neovascularization capacity of circulating endothelial progenitor cells. Circulation. 2003;108(20):2511–6. https://doi.org/10.1161/01.CIR.0000096483.29777.50.
Zhao Y, Glesne D, Huberman E. A human peripheral blood monocyte-derived subset acts as pluripotent stem cells. Proc Natl Acad Sci U S A. 2003;100(5):2426–31. https://doi.org/10.1073/pnas.0536882100.
Zhang Z, Shively JE. Generation of novel bone forming cells (monoosteophils) from the cathelicidin-derived peptide LL-37 treated monocytes. PloS one. 2010;5(11):e13985. https://doi.org/10.1371/journal.pone.0013985.
He JB, Chen MH, Lin DK. New insights into the tonifying kidney-yin herbs and formulas for the treatment of osteoporosis. Arch Osteoporos. 2017;12(1):14. https://doi.org/10.1007/s11657-016-0301-4.
Wang L, Li Y, Guo Y, Ma R, Fu M, Niu J, Gao S, Zhang D. Herba Epimedii: an ancient Chinese herbal medicine in the prevention and treatment of osteoporosis. Curr Pharm Des. 2016;22(3):328–49. https://doi.org/10.2174/1381612822666151112145907.
Hung TM, Na M, Thuong PT, Su ND, Sok D, Song KS, Seong YH, Bae K. Antioxidant activity of caffeoyl quinic acid derivatives from the roots of Dipsacus asper Wall. Journal of ethnopharmacology. 2006;108(2):188–92. https://doi.org/10.1016/j.jep.2006.04.029.
Niu YB, Li YH, Kong XH, Zhang R, Sun Y, Li Q, Li C, Liu L, Wang J, Mei QB. The beneficial effect of Radix Dipsaci total saponins on bone metabolism in vitro and in vivo and the possible mechanisms of action. Osteoporos Int. 2012;23(11):2649–60. https://doi.org/10.1007/s00198-012-1932-y.
Wong RW, Rabie AB, Hagg EU. The effect of crude extract from Radix Dipsaci on bone in mice. Phytotherapy research : PTR. 2007;21(6):596–8. https://doi.org/10.1002/ptr.2126.
Clough E, Barrett T. The Gene Expression Omnibus Database. Methods Mol Biol. 2016;1418:93–110. https://doi.org/10.1007/978-1-4939-3578-9_5.
McKinley TO, Lisboa FA, Horan AD, Gaski GE, Mehta S. Precision medicine applications to manage multiply injured patients with orthopaedic trauma. J Orthop Trauma. 2019;33(Suppl 6):S25–9. https://doi.org/10.1097/BOT.0000000000001468.
Ru J, Li P, Wang J, Zhou W, Li B, Huang C, Li P, Guo Z, Tao W, Yang Y, Xu X, Li Y, Wang Y, Yang L. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13. https://doi.org/10.1186/1758-2946-6-13.
Kim S. Getting the most out of PubChem for virtual screening. Expert Opin Drug Discov. 2016;11(9):843–55. https://doi.org/10.1080/17460441.2016.1216967.
Wang X, Shen Y, Wang S, Li S, Zhang W, Liu X, Lai L, Pei J, Li H. PharmMapper 2017 update: a web server for potential drug target identification with a comprehensive target pharmacophore database. Nucleic Acids Res. 2017;45(W1):W356–60. https://doi.org/10.1093/nar/gkx374.
Gfeller D, Grosdidier A, Wirth M, Daina A, Michielin O, Zoete V. SwissTargetPrediction: a web server for target prediction of bioactive small molecules. Nucleic Acids Res. 2014 42 (Web Server issue):W32-38. https://doi.org/10.1093/nar/gku293.
Zhang W, Xue K, Gao Y, Huai Y, Wang W, Miao Z, Dang K, Jiang S, Qian A. Systems pharmacology dissection of action mechanisms of Dipsaci Radix for osteoporosis. Life sciences. 2019;235:116820. https://doi.org/10.1016/j.lfs.2019.116820.
von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31(1):258–61. https://doi.org/10.1093/nar/gkg034.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):P3.
Hak DJ, Fitzpatrick D, Bishop JA, Marsh JL, Tilp S, Schnettler R, Simpson H, Alt V. Delayed union and nonunions: epidemiology, clinical issues, and financial aspects. Injury. 2014;45(Suppl 2):S3–7. https://doi.org/10.1016/j.injury.2014.04.002.
Cunningham BP, Brazina S, Morshed S, Miclau T 3rd. Fracture healing: a review of clinical, imaging and laboratory diagnostic options. Injury. 2017;48(Suppl 1):S69–75. https://doi.org/10.1016/j.injury.2017.04.020.
Antonova E, Le TK, Burge R, Mershon J. Tibia shaft fractures: costly burden of nonunions. BMC Musculoskelet Disord. 2013;14:42. https://doi.org/10.1186/1471-2474-14-42.
Chaverri D, Vives J. Toward the clinical use of circulating biomarkers predictive of bone union. Biomark Med. 2017;11(12):1125–33. https://doi.org/10.2217/bmm-2017-0180.
Pountos I, Georgouli T, Pneumaticos S, Giannoudis PV. Fracture non-union: can biomarkers predict outcome? Injury. 2013;44(12):1725–32. https://doi.org/10.1016/j.injury.2013.09.009.
Cheng A, Krishnan L, Pradhan P, Weinstock LD, Wood LB, Roy K, Guldberg RE. Impaired bone healing following treatment of established nonunion correlates with serum cytokine expression. J Orthop Res. 2019;37(2):299–307. https://doi.org/10.1002/jor.24186.
Galliera E, Randelli P, Dogliotti G, Dozio E, Colombini A, Lombardi G, Cabitza P, Corsi MM. Matrix metalloproteases MMP-2 and MMP-9: are they early biomarkers of bone remodelling and healing after arthroscopic acromioplasty? Injury. 2010;41(11):1204–7. https://doi.org/10.1016/j.injury.2010.09.024.
Colnot C, Thompson Z, Miclau T, Werb Z, Helms JA. Altered fracture repair in the absence of MMP9. Development (Cambridge, England). 2003;130(17):4123–33. https://doi.org/10.1242/dev.00559.
Vu TH, Shipley JM, Bergers G, Berger JE, Helms JA, Hanahan D, Shapiro SD, Senior RM, Werb Z. MMP-9/gelatinase B is a key regulator of growth plate angiogenesis and apoptosis of hypertrophic chondrocytes. Cell. 1998;93(3):411–22. https://doi.org/10.1016/s0092-8674(00)81169-1.
Haubruck P, Heller R, Apitz P, Kammerer A, Alamouti A, Daniel V, Schmidmaier G, Moghaddam A. Evaluation of matrix metalloproteases as early biomarkers for bone regeneration during the applied Masquelet therapy for non-unions. Injury. 2018;49(10):1732–8. https://doi.org/10.1016/j.injury.2018.07.015.
Supuran CT. Carbonic anhydrases--an overview. Curr Pharm Des. 2008;14(7):603–14. https://doi.org/10.2174/138161208783877884.
Mirjafari P, Asghari K, Mahinpey N. Investigating the application of enzyme carbonic anhydrase for CO2 sequestration purposes. Ind Eng Chem Res. 2007;46(3):921–6. https://doi.org/10.1021/ie060287u.
Chang X, Zheng Y, Yang Q, Wang L, Pan J, Xia Y, Yan X, Han J. Carbonic anhydrase I (CA1) is involved in the process of bone formation and is susceptible to ankylosing spondylitis. Arthritis research & therapy. 2012;14(4):R176. https://doi.org/10.1186/ar3929.
Fischer WH, Spiess J. Identification of a mammalian glutaminyl cyclase converting glutaminyl into pyroglutamyl peptides. Proc Natl Acad Sci U S A. 1987;84(11):3628–32. https://doi.org/10.1073/pnas.84.11.3628.
Ezura Y, Kajita M, Ishida R, Yoshida S, Yoshida H, Suzuki T, Hosoi T, Inoue S, Shiraki M, Orimo H, Emi M. Association of multiple nucleotide variations in the pituitary glutaminyl cyclase gene (QPCT) with low radial BMD in adult women. Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 2004;19(8):1296–301. https://doi.org/10.1359/JBMR.040324.
Huang QY, Kung AW. The association of common polymorphisms in the QPCT gene with bone mineral density in the Chinese population. J Hum Genet. 2007;52(9):757–62. https://doi.org/10.1007/s10038-007-0178-6.
Dedkova EN, Blatter LA. Role of beta-hydroxybutyrate, its polymer poly-beta-hydroxybutyrate and inorganic polyphosphate in mammalian health and disease. Front Physiol. 2014;5:260. https://doi.org/10.3389/fphys.2014.00260.
Schaub A, Futterer A, Pfeffer K. PUMA-G, an IFN-gamma-inducible gene in macrophages is a novel member of the seven transmembrane spanning receptor superfamily. Eur J Immunol. 2001;31(12):3714–25. https://doi.org/10.1002/1521-4141(200112)31:12<3714::aid-immu3714>3.0.co;2-1.
Ichikawa S, Sakiyama H, Suzuki G, Hidari KI, Hirabayashi Y. Expression cloning of a cDNA for human ceramide glucosyltransferase that catalyzes the first glycosylation step of glycosphingolipid synthesis. Proc Natl Acad Sci U S A. 1996;93(10):4638–43. https://doi.org/10.1073/pnas.93.10.4638.
Song T, Gan W, Chen J, Huang L, Yin H, He T, Huang H, Hu X. Antibodies against Clonorchis sinensis LDH could cross-react with LDHB localizing on the plasma membrane of human hepatocarcinoma cell SMMC-7721 and induce apoptosis. Parasitology research. 2016;115(4):1595–603. https://doi.org/10.1007/s00436-015-4895-z.
Chen R, Zhou X, Yu Z, Liu J, Huang G. Low expression of LDHB correlates with unfavorable survival in hepatocellular carcinoma: strobe-compliant article. Medicine (Baltimore). 2015;94(39):e1583. https://doi.org/10.1097/MD.0000000000001583.
The authors would like to thank the patients and their families who participated in this study.
This work was supported by grants from the Xuanwu Hospital Funding Project (No. XWJK-2019033). The funding body had no effect on either the design of the study and collection, analysis, or interpretation of data nor in writing the manuscript.
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. All procedures performed in studies were approved by the Ethics Committee of Xuanwu Hospital, Capital Medical University (Beijing, China).
Consent for publication
Informed consent was obtained from all individual participants included in the study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1. Validation of the six potential biomarkers between fracture only and multi-trauma patients by ELISA. Concentration of CA1 (A), MMP9 (B), QPCT (C), HCAR2 (D), UGCG (E), and LDHB (F) in serum samples from the fracture only (n = 16 in healed group, n = 5 in nonunion group) and multi-trauma patients (n = 23 in healed group, n = 11 in nonunion group). The error bars represent means ± SEM. There has no significant difference between fracture only and multi-trauma patients. Table S1. Functional analysis of the six potential blood biomarkers.
About this article
Cite this article
Wu, J., Liu, L., Hu, H. et al. Bioinformatic analysis and experimental identification of blood biomarkers for chronic nonunion. J Orthop Surg Res 15, 208 (2020). https://doi.org/10.1186/s13018-020-01735-1
- Dipsaci Radix (DR)