Skip to main content
  • Research article
  • Open access
  • Published:

Development and validation of a nomogram predicting the risk of recurrent lumbar disk herniation within 6 months after percutaneous endoscopic lumbar discectomy

Abstract

Objective

To develop and validate a nomogram useful in predicting recurrent lumbar disk herniation (rLDH) within 6 months after percutaneous endoscopic lumbar discectomy (PELD).

Methods

Information on patients’ lumbar disk herniation (LDH) between January 2018 and May 2019 in addition to 26 other features was collected from the authors’ hospital. The least absolute shrinkage and selection operator (LASSO) method was used to select the most important risk factors. Moreover, a nomogram was used to build a prediction model using the risk factors selected from LASSO regression. The concordance index (C-index), the receiver operating characteristic (ROC) curve, and calibration curve were used to assess the performance of the model. Finally, clinical usefulness of the nomogram was analyzed using the decision curve and bootstrapping used for internal validation.

Results

Totally, 352 LDH patients were included into this study. Thirty-two patients had recurrence within 6 months while 320 showed no recurrence. Four potential factors, the course of disease, Pfirrmann grade, Modic change, and migration grade, were selected according to the LASSO regression model. Additionally, the C-index of the prediction nomogram was 0.813 (95% CI, 0.726-0.900) and the area under receiver operating characteristic curve (AUC) value was 0.798 while the interval bootstrapping validation C-index was 0.743. Hence, the nomogram might be a good predictive model.

Conclusion

Each variable, the course of disease, Pfirrmann grade, Modic change, and migration grade in the nomogram had a quantitatively corresponding risk score, which can be used in predicting the overall recurrence rate of rLDH within 6 months.

Introduction

Lumbar disk herniation (LDH) is a degenerative disease of the lumbar spine. The disease manifests the following primary signs and symptoms; radicular pain, sensory abnormalities, and weakness in the distribution of one or more lumbosacral nerve roots [1, 2]. Diagnosis of LDH mainly depends on clinical symptoms and imaging information. Presently, the treatment revolves around non-operative interventions, open discectomy, and minimally invasive surgery [3]. The beneficial outcomes of minimally invasive surgery have been confirmed. They include a decrease operative time, less blood loss, and quick return to normal routine [4].

Percutaneous endoscopic lumbar discectomy (PELD) is a common used minimally invasive procedure that has been proven to be effective in patients with LDH [5, 6]. The procedure has the advantage of requiring a smaller incision as well as conferring faster recovery, less damage to soft tissues, and fewer complications [7]. Clinically, PELD is further divided into two types of operations, namely, percutaneous endoscopic transforaminal discectomy (PETD) and percutaneous endoscopic interlaminar discectomy (PEID).

However, the rate of recurrent LDH (rLDH) has been reported to be between 5 and 15%, after the operation. Moreover, the duration between surgery and a recurrence ranges from a few days to a few years [8]. In the clinical setting, the definition of rLDH is the presence of herniated disk material at the same level with symptoms of oppression, but strict time interval is not necessary [9]. Reoperation within a short time presents financial constraints as well as physical and psychological pressure to patients. Moreover, many studies have highlighted the possible risk factors of rLDH after a minimally invasive or open surgery, including sex, age, height, weight, body mass index (BMI), smoking, occupation type, and Modic change [10, 11].

Therefore, given the many risk factors, an accurate prediction tool may be useful in forecasting rLDH. Generally, most studies focus on recurrence 6 months after the surgery. However, few have paid attention to the recurrence rates within the first 6 months after operation [12]. Similarly, this study explored the recurrence of LDH within 6 months of PELD. To the best of our knowledge, no study has provided an approach that could predict the probability of recurrence within 6 months after PELD. Moreover, inability to predict recurrence after surgery is an important cause of patient dissatisfaction. Consequently, the study set out to develop an accurate but simple method of predicting rLDH within 6 months after PELD, by assessing a group of possible risk factors.

Methods

Patients and risk factors

The data of all patients who had undergone PELD in the authors’ hospital between January 2018 and May 2019 was collected. The inclusion criteria were as follows: (1) patients who were diagnosed with LDH according to clinical manifestations and radiological characteristics and (2) patients who underwent PELD by the same team of surgeons. On the other hand, the exclusion criteria include (1) patients who had undergone an open surgery; (2) patients who lacked imaging data; (3) patients with recurrence of LDH whose primary surgeries were performed in other hospitals; (4) patients with recurrence of LDH whose primary operation was performed before 2018; and (5) patients who had recurrence after 6 months of operation. All the patients were informed consent and the study was approved by the Medical Ethics Committee of authors’ hospital.

Characteristics including demographics, radiological identifiable factors, and surgery-related information were collected from hospital medical records. Demographic and surgery-related information includes age, gender, height, weight, BMI, occupation, education, type of household registration, smoking, drinking, diabetes, hypertension, course of disease, type of surgery, and operative time. Radiological identifiable factors included Pfirrmann grade, level of the herniated disk, Modic change, herniated size, herniation direction, herniation location, migration grade, disk height, disk length, disk width, and disk size.

Statistical analysis

All the data is displayed in the Table 1. Statistical analysis was performed using the R software (Version 3.6.1). The R packages used include glmnet, rms, ROCR, and rmda.

Table 1 Characteristics of recurrence and non-recurrence patients

First, the least absolute shrinkage and selection operator (LASSO) method, a regression-based analysis method for penalizing the magnitude of the coefficients of prediction variables by imposing a constraint during parameter estimation, was used to select the most important risk factors to the rate of rLDH. This method could help reduce the number of variables and minimize the possibility of model overfitting [13].

Afterwards, the nomogram was used to build a predictive model using the risk factors obtained from LASSO. This nomogram is a logistic regression-based model, which can reduce statistical predictive models into a single numerical estimate of the probability of an event. The nomogram converts the regression coefficient of each covariate according to the formula, and visualizes the abstract results of logistic regression [14].

Thereafter, the concordance index (C-index), the receiver operating characteristic (ROC) curve, and calibration curve were used to assess the performance of the model. The C-index and the area under receiver operating characteristic curve (AUC) value measure discrimination with a range of 0.5 to 1.0. Therefore, the larger C-index and AUC value are, the more accurate the results were in distinguishing the subjects [15]. Furthermore, a calibration curve was plotted to assess the standardization of the nomogram.

Finally, bootstrapping was used (1000 bootstrap resamples) to calculate a corrected C-index, in order to verify the accuracy of the model [16].

Results

Characteristics of patients

A total of 352 LDH patients were enrolled in this study between January 2018 and May 2019. Thirty-two of them had rLDH within 6 months while 320 showed no recurrence of disease within the period. The patients include 230 males and 122 females within the range of 16-90 years of age. Among them, 17 showed recurrence within 1 month, 8 within 1 to 3 months after PELD, and 7 between 3 and 6 months after operation. All the details, including demographics, surgery-related information and radiological data of the two groups of patients were shown in Table 1. Disk height, length, width, and size were divided into two groups based on the median. Herniated size was classified as grade 1, 2, and 3 based on the Michigan State University (MSU) classification. In addition, the classification of Pfirrmann grade and Modic change were shown in Figs. 1 and 2. Figure 3 was the schematic representation of migration grade.

Fig. 1
figure 1

The classification of Pfirrmann grade. On T2-weighted image, grade I—the structure of the disk and disk height are normal, with a bright white signal intensity. Grade II—the structure of the disk is abnormal and the disk height is normal, with a white signal and a clear distinction between nucleus and anulus. Grade III—the structure of the disk is abnormal and the disk height is normal or slightly decreased, with an average gray signal intensity and an unclear distinction between nucleus and anulus. Grade IV—the structure of the disk is abnormal and the disk height is normal or moderately decreased, with a dark gray signal intensity and a lost distinction between nucleus and anulus. Grade V—the structure of the disk is abnormal and the disk space is collapsed. The distinction between nucleus and anulus is lost

Fig. 2
figure 2

The classification of Modic change. On T1-weighted and T2-weighted images, normal: isointense on T1w and T2w. Type I, hypointense on T1w and hyperintense on T2w. Type II, hyperintense on T1w and hyperintense or isointense on T2w. Type III, hypointense on T1w and T2w

Fig. 3
figure 3

Schematic representation of migration grade. The migration being smaller than the measured height of the posterior marginal disk space was described as a low-grade migration or as a high-grade migration

Of all the 26 risk factors, 4 were selected based on the LASSO regression method (Fig. 4a and b). These factors include the course of disease, Pfirrmann grade, Modic change, and migration grade.

Fig. 4
figure 4

The least absolute shrinkage and selection operator (LASSO) method of selecting risk factors. (a) The 26 feature LASSO coefficient profiles plot was produced against the log (lambda) sequence. (b) Using cross-validation via minimum criteria in the LASSO model, dotted vertical lines were drawn at the optimal values (4 factors)

Development and validation of the nomogram for rLDH prediction

The independent factors were incorporated and their predictive capability shown in the nomogram (Fig. 5). The calibration curve of the nomogram for the prediction of rLDH demonstrated good results in the 352 patients. However, the predictive capability of the nomogram decreased with an increase of recurrence risk. For instance, when the recurrence risk was less than about 0.3, the predictive capability was better than when the recurrence risk was more than 0.3 (Fig. 6). The C-index for the prediction nomogram was 0.813 (95% CI, 0.726-0.900) and the AUC value was 0.798 (Fig. 7) while the interval bootstrapping validation C-index was 0.743. All these indicated the nomogram model’s good predictive capability.

Fig. 5
figure 5

Nomogram to predict rLDH after PELD. Four simple factors selected from the LASSO analysis, including course of disease, Pfirrmann grade, Modic change, and migration grade were used to develop the nomogram

Fig. 6
figure 6

The calibration curves of the nomogram. The x-axis represents the predicted recurrence risk while the y-axis represents the actual probability. The diagonal line indicates an ideal prediction model. The solid line reflects the real performance of the nomogram. The closer the solid line gets to the diagonal one, the better the tool is at predicting recurrence

Fig. 7
figure 7

The ROC analysis for the nomogram

Discussion

Based on data from online public databases, a large number of nomograms have been used as predictive tools in tumor survival in recent years. However, the nomograms are scarcely used in other fields and no prognostic nomogram has been constructed for rLDH patients up to date. Nomogram is an important research method in the field of translational medicine, which converts simple data into clinical prediction models through mathematical modeling [17]. Therefore, this study was the first one to use a nomogram in predicting rLDH after PELD.

Based on factors such as surgeon preference, radiographic evidence, and severity of herniations, interventions vary between open and endoscopic discectomy. Over the past decade, extensive research has been done on the surgical treatment of LDH. Notably, the minimally invasive spinal surgery has been reported to have more advantages compared to the common open surgery. In this study, PELD, a common type of minimally invasive operation, capable of removing herniated disks, with a skin incision of only 7 mm, shorter hospitalization periods, and faster recovery was used [18]. However, despite the reported advantages of PELD, many surgeons still experience failure after surgery, with one of the common complications being early recurrence [19]. In this study, a good nomogram that could predict the recurrence rate as well as avoid some of the risk factors was developed.

According to previously published data, 26 possible risk factors were initially selected in the study. Additionally, such factors as sex, age, height, weight, BMI, smoking, occupation type, and Modic change were reported to be associated with rLDH [10, 11]. However, minor variations were reported in this study. Only four risk factors were selected based on the LASSO results, including the course of disease, Pfirrmann grade, Modic change, and migration grade. These factors could be easily obtained during routine clinical practice. The high C-index, the AUC value, and the calibration curve of the nomogram indicated both an accurate predictive and calibrative power of the method. Additionally, the great C-index in interval bootstrapping validation particularly confirmed that this easy-to-use nomogram could be widely applied in LDH patients.

Although several studies have conducted the multi-factor analysis of rLDH, each reported different results according to the risk factors. Similar to previous studies, Modic change was found to be an important risk factor of rLDH in this study [20]. Despite the importance of Modic change in evaluating the severity of disk degeneration [21], the etiology, and underlying mechanisms in relation to rLDH still needs further research. Moreover, the results of this study showed that the migration grade on rLDH was different from that reported in previous studies. For instance, a previous study reported that failure and recurrence rates were higher in the high-grade migration group [22]. However, our results showed that downward high-grade migration presented the least probability of rLDH. These differences might arise due to varying degrees of surgeons’ experience as well as technical factors [23].

Another index commonly used for evaluating degenerative disk disease is the Pfirrmann grade. In this study, the levels of herniated disks were all reported as Pfirrmann Grade III n. Previous studies reported that an increase in Pfirrmann grade resulted to cell apoptosis and a decrease in the moisture content of the lumbar disks. Consequently, this would result to a change in the intervertebral disk micro-nano environment [24]. Additionally, studies showed that a longer course of disease predicted slower recovery as well as reduced surgical efficacy [25]. Notably, patients with a disease course of more than 12 months before surgery were reported to less likely show optimal surgical outcomes compared to those whose disease course was less than 1 year [26].

This study showed that a simple four-risk-factor nomogram could easily be used in daily clinical work. However, a few limitations can be picked from this study. First, data was only collected for one and a half years and 32 patients showed recurrence within 6 months. This small sample size could potentially lead to a bias in the reported results. Second, only four risk factors were included in the analysis yet many others were able to cause rLDH, including personal income, postoperative care, and other conditions. Finally, only internal validation measures were used to ascertain the success of the nomogram. However, it would have been better to use data from external centers for a more comprehensive conclusion.

Conclusion

In summary, the study successfully constructed a nomogram capable of predicting rLDH after PELD using data from the authors’ hospital. By estimating the individual risk, the possibility of recurrence could be predicted. This allows patients to be well informed on such issues as recurrence, before undergoing an operation. Additionally, the model allows both doctors and patients to make important preopertative decisions such as whether to perform a PELD or a fusion surgery. Nonetheless, the nomogram still needs improvement which can be achieved by using data from more patients as well as applying external validation measures. Further studies are therefore required to overcome these shortcomings.

Abbreviations

LDH:

Lumbar disk herniation

rLDH:

Recurrent lumbar disk herniation

LASSO:

Least absolute shrinkage and selection operator

C-index:

Concordance index

ROC:

The receiver operating characteristic

AUC:

The area under receiver operating characteristic curve

PELD:

Percutaneous endoscopic lumbar discectomy

PEID:

Percutaneous endoscopic interlaminar discectomy

PETD:

Percutaneous endoscopic transforaminal discectomy

BMI:

Body mass index

CI:

Confidence interval

MSU:

Michigan State University

References

  1. Vroomen PC, de Krom MC, Wilmink JT, Kester AD, Knottnerus JA. Diagnostic value of history and physical examination in patients suspected of lumbosacral nerve root compression. J Neurol Neurosurg Psychiatry. 2002;72(5):630–4. https://doi.org/10.1136/jnnp.72.5.630.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Chen Z, Cao P, Zhou Z, Yuan Y, Jiao Y, Zheng Y. Overview: the role of Propionibacterium acnes in nonpyogenic intervertebral discs. Int Orthop. 2016;40(6):1291–8. https://doi.org/10.1007/s00264-016-3115-5.

    Article  PubMed  Google Scholar 

  3. Amin RM, Andrade NS, Neuman BJ. Lumbar disc herniation. Curr Rev Musculoskelet Med. 2017;10(4):507–16. https://doi.org/10.1007/s12178-017-9441-4.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Phan K, Xu J, Schultz K, Alvi MA, Lu VM, Kerezoudis P, et al. Full-endoscopic versus micro-endoscopic and open discectomy: a systematic review and meta-analysis of outcomes and complications. Clin Neurol Neurosurg. 2017;154:1–12. https://doi.org/10.1016/j.clineuro.2017.01.003.

  5. Maroon JC, Onik G, Vidovich DV. Percutaneous discectomy for lumbar disc herniation. Neurosurg Clin N Am. 1993;4(1):125–34. https://doi.org/10.1016/S1042-3680(18)30614-4.

    Article  CAS  PubMed  Google Scholar 

  6. Yong A, Sang-Ho L, Woo-Min P, Ho-Yeon L, Song-Woo S, Ho-Yeong K. Percutaneous endoscopic lumbar discectomy for recurrent disc herniation: surgical technique, outcome, and prognostic factors of 43 consecutive cases. Spine. 2004;29:E326–32.

    Article  Google Scholar 

  7. Eun SS, Lee SH, Sabal LA. Long-term follow-up results of percutaneous endoscopic lumbar discectomy. Pain Physician. 2016;19(8):E1161–6.

    PubMed  Google Scholar 

  8. Gaston P, Marshall RW. Survival analysis is a better estimate of recurrent disc herniation. J Bone Jnt Surg Brit Vol. 2003;85:535–7.

    Article  CAS  Google Scholar 

  9. Swartz KR, Trost GR. Recurrent lumbar disc herniation. Neurosurg Focus. 2003;15:E10.

    Article  PubMed  Google Scholar 

  10. Yaman ME, Kazanci A, Yaman ND, Bas F, Ayberk G. Factors that influence recurrent lumbar disc herniation. Hong Kong Med J. 2017;23(3):258–63. https://doi.org/10.12809/hkmj164852.

    Article  CAS  PubMed  Google Scholar 

  11. Park CH, Park ES, Lee SH, Lee KK, Kwon YK, Kang MS, et al. Risk factors for early recurrence after transforaminal endoscopic lumbar disc decompression. Pain Physician. 2019;22(2):E133–e138.

  12. Li Z, Yang H, Liu M, Lu M, Chu J, Hou S, et al. Clinical characteristics and risk factors of recurrent lumbar disk herniation: a retrospective analysis of three hundred twenty-one cases. Spine (Phila Pa 1976). 2018;43(21):1463–9. https://doi.org/10.1097/BRS.0000000000002655.

  13. Sauerbrei W, Royston P, Binder H. Selection of important variables and determination of functional form for continuous predictors in multivariable model building. Stat Med. 2007;26(30):5512–28. https://doi.org/10.1002/sim.3148.

    Article  PubMed  Google Scholar 

  14. Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol Off J Am Soc Clin Oncol. 2008;26:1364–70.

    Article  Google Scholar 

  15. Yin J, Yang Y, Ma K, Yang X, Lu T, Wang S, et al. Clinicopathological characteristics and prognosis of pulmonary pleomorphic carcinoma: a population-based retrospective study using SEER data. J Thor Dis. 2018;10(7):4262–73. https://doi.org/10.21037/jtd.2018.06.71.

  16. Pencina MJ, D’Agostino RB. Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimation. Stat Med. 2004;23(13):2109–23. https://doi.org/10.1002/sim.1802.

    Article  PubMed  Google Scholar 

  17. Mediouni M, Schlatterer DR, Madry H, Cucchiarini M, Rai B. A review of translational medicine. The future paradigm: how can we connect the orthopedic dots better? Curr Med Res Opin. 2018;34:1217–29.

    Article  PubMed  Google Scholar 

  18. Price JP, Dawson JM, Schwender JD, Schellhas KP. Clinical and radiologic comparison of minimally invasive surgery with traditional open transforaminal lumbar interbody fusion: a review of 452 patients from a single center. Clin Spine Surg. 2018;31(2):E121–e126. https://doi.org/10.1097/BSD.0000000000000581.

    Article  PubMed  Google Scholar 

  19. Kim HS, You JD, Ju CI. Predictive scoring and risk factors of early recurrence after percutaneous endoscopic lumbar discectomy. Biomed Res Int. 2019;2019:6492675.

    PubMed  PubMed Central  Google Scholar 

  20. Karchevsky M, Schweitzer ME, Carrino JA, Zoga A, Montgomery D, Parker L. Reactive endplate marrow changes: a systematic morphologic and epidemiologic evaluation. Skelet Radiol. 2005;34(3):125–9. https://doi.org/10.1007/s00256-004-0886-3.

    Article  Google Scholar 

  21. Modic MT, Steinberg PM, Ross JS, Masaryk TJ, Carter JR. Degenerative disk disease: assessment of changes in vertebral body marrow with MR imaging. Radiology. 1988;166(1):193–9. https://doi.org/10.1148/radiology.166.1.3336678.

    Article  CAS  PubMed  Google Scholar 

  22. Lee SH, Kang BU, Ahn Y, Choi G, Choi YG, Ahn KU, et al. Operative failure of percutaneous endoscopic lumbar discectomy: a radiologic analysis of 55 cases. Spine (Phila Pa 1976). 2006;31(10):E285–90. https://doi.org/10.1097/01.brs.0000216446.13205.7a.

  23. Yeung AT, Tsou PM. Posterolateral endoscopic excision for lumbar disc herniation: surgical technique, outcome, and complications in 307 consecutive cases. Spine (Phila Pa 1976). 2002;27(7):722–31. https://doi.org/10.1097/00007632-200204010-00009.

    Article  Google Scholar 

  24. Wei F, Zhong R, Zhou Z, Wang L, Pan X, Cui S, et al. In vivo experimental intervertebral disc degeneration induced by bleomycin in the rhesus monkey. BMC Musculoskelet Disord. 2014;15(1):340. https://doi.org/10.1186/1471-2474-15-340.

  25. Manchikanti L, Singh V, Falco FJ, Calodney AK, Onyewu O, Helm S 2nd, et al. An updated review of automated percutaneous mechanical lumbar discectomy for the contained herniated lumbar disc. Pain Physician. 2013;16:Se151–84.

  26. Jiang X, Zhou X, Xu N. Clinical effects of transforaminal and interlaminar percutaneous endoscopic discectomy for lumbar disc herniation: a retrospective study. Medicine. 2018;97(48):e13417. https://doi.org/10.1097/MD.0000000000013417.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This work was supported by the National Natural Science Foundation of China (81771348).

Author information

Authors and Affiliations

Authors

Contributions

Mengxian Jia, Yadong Sheng and Guoliang Chen performed the research study, analyzed the data and wrote this manuscript. Wenbin Zhang, Jiajin Lin, Sheng Lu, Fayi Li and Jinwei Ying collected data. Honglin Teng designed the research and provided guidance. All authors reviewed and approved the final manuscript.

Corresponding author

Correspondence to Honglin Teng.

Ethics declarations

Competing interests

None.

Additional information

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 licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jia, M., Sheng, Y., Chen, G. et al. Development and validation of a nomogram predicting the risk of recurrent lumbar disk herniation within 6 months after percutaneous endoscopic lumbar discectomy. J Orthop Surg Res 16, 274 (2021). https://doi.org/10.1186/s13018-021-02425-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-021-02425-2

Keywords