Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access November 2, 2021

Proteoglycan-depleted regions of annular injury promote nerve ingrowth in a rabbit disc degeneration model

  • Long Xin , Weixin Xu , Jian Wang , Fang Yu , Shunwu Fan , Xinwei Xu and Yang Yang EMAIL logo
From the journal Open Medicine

Abstract

Background

To assess the effects of proteoglycan-depleted regions of annular disruptions on nerve ingrowth in the injury site in vivo.

Methods

New Zealand white rabbits (n = 18) received annular injuries at L3/4, L4/5, and L5/6. The experimental discs were randomly assigned to four groups: (a) an annular defect was created; (b) an annular defect implanted with a poly lactic-co-glycolic acid (PLGA)/fibrin/PBS plug; (c) an annular defect implanted with a PLGA/fibrin/chondroitinase ABC (chABC) plug; and (d) an uninjured L2/3 disc (control). Disc degeneration was evaluated by radiography, MRI, histology, and analysis of the proteoglycan (PG) content. Immunohistochemical detection of nerve fibers and chondroitin sulfate (CS) was performed.

Results

The injured discs produced progressive and reliable disc degeneration. In the defective discs, the lamellated appearance of AF (Annulus fibrosus) was replaced by extensive fibrocartilaginous-like tissue formation outside the injured sites. In contrast, newly formed tissue was distributed along small fissures, and small blood vessels appeared in the outer part of the disrupted area in the PLGA/fibrin/PBS discs. More sprouting nerve fibers grew further into the depleted annulus regions in the PLGA/fibrin/chABC discs than in the control discs and those receiving PLGA/fibrin/PBS. In addition, the innervation scores of the PLGA/fibrin/chABC discs were significantly increased compared with those of the PLGA/fibrin/PBS discs and defected discs.

Conclusion

ChABC-based PLGA/fibrin gel showed promising results by achieving biointegration with native annulus tissue and providing a local source for the sustained release of active chABC. Disc-derived PG-mediated inhibition of nerve and blood vessel ingrowth was abrogated by chABC enzymatic deglycosylation in an annular-injured rabbit disc degeneration model.

1 Introduction

Low back pain is an extremely common socioeconomic problem in orthopedics. Although the pathophysiology of this condition remains uncertain, degeneration of the intervertebral disc (IVD) is believed to be a major cause of low back pain [1,2,3]. This degenerative change ultimately results in increased matrix degradation, proinflammatory cytokine expression, and inferior mechanical properties [4,5,6]. More recently, progressive neurovascular growth in annular fissures and sensitization of nerve fibers have been observed in the degenerated IVD [7,8]. Ingrowth of nociceptive nerve fibers deep within disrupted IVDs is believed to be a potential source of discogenic pain in clinical cases [9,10,11,12].

It has been shown that aggrecan, a proteoglycan (PG) found in the disc, may act as a barrier to inhibit nerve fiber growth in vivo and in vitro [13,14,15,16]. However, the precise inhibitory mechanism of aggrecan remains poorly understood. In previous work, we established a minimally invasive annulotomy-induced rabbit model and demonstrated the growth of orientated growth-permissive nerves into annulus defects during the degenerative process. In addition, loss of PG at the injury site is involved in IVD degeneration and facilitates the infiltration of blood vessels and nerve fibers into disrupted disc tissue [17,18].

Chondroitinase ABC (chABC ) is a well-studied bacteria-derived enzyme that specifically depolymerizes glycosaminoglycan (GAG) side chains, inducing mild disc degeneration and favoring axonal regeneration by reducing chondroitin sulfate proteoglycan (CSPG) at the injury site [19,20,21,22,23,24]. According to our previous work, PLGA/fibrin constructs as a delivery system exhibit good biocompatibility, safety, and biodegradability [18,25], and promote tissue repair without compromising the function of the incorporated proteins [26,27,28,29,30]. Therefore, as an ideal carrier of chABC, PLGA/fibrin constructs represent a simple and desirable alternative to achieve local delivery in a controlled manner. In this system, cleavage of GAGs via digestion overcomes CSPG-mediated inhibition at the injured site of the degenerated IVD.

In the present study, we hypothesized that PG-depleted regions of the annular disruptions could facilitate nerve and blood vessel ingrowth into the deeper parts of the degenerated disc in a rabbit animal model. Here fibrin-based polymeric carriers loaded with chABC were implanted into the injured annular site to provide sustained local delivery in vivo, and the aim of this study was thus to assess the effects of degrading CSPG on neurite ingrowth during IVD degeneration.

2 Methods and materials

2.1 Preparation of PLGA-coated fibrin gel/ChABC constructs

PLGA sponges were fabricated using a porogen leaching method as reported previously [31]. The PLGA/fibrin/chABC constructs were prepared as follows. Fibrinogen was isolated from fresh human plasma, and the final concentration used in all experiments was based on previously established protocols [25,32]. ChABC (2U, Sigma, St. Louis, MO) was dissolved in buffer with a final concentration of 0.10 U/mL according to the product’s instructions and as previously reported [19,33]. Thereafter, fibrinogen (20 mg/mL) and chABC (0.10 U/mL) solutions were homogenized and sterilized by filtering through syringe filters. PLGA plugs (1.8 mm diameter and 4 mm length) were immersed into the homogeneous fibrinogen/chABC solution (1 mL) under reduced pressure. Composite constructs were lyophilized, imaged using a SEM, and subsequently stored at −20°C until further use.

2.2 Evaluation of chABC release from PLGA/fibrin in vitro

To obtain release profiles for each PLGA/fibrin sample, a gel loaded with chABC was suspended in 100 μL of PBS (pH 7.4) in polypropylene tubes placed in a shaker bath (37°C) at 100 rpm. Then, 10 μL of the supernatant was collected from each PLGA/fibrin sample at predetermined time intervals (1, 2, 3, 4, 5, 7, 9, 12 and 14 days) and individually stored at −20°C. According to a previously reported protocol [34], the enzymatic activity of the solution was assessed by measuring the formation of unsaturated disaccharides via chondroitin sulfate (CS) degradation over time (Sigma enzymatic assay-EC 4.2.2.4), as per the manufacturer’s instructions. Finally, the percentage of sample absorbance relative to the standard was calculated, and the composite construct’s accumulated release kinetic curve was obtained.

2.3 Animal surgery

A total of 18 New Zealand rabbits (age 8.60 ± 1.25 months and weight 3.42 ± 0.16 kg) were supplied by the Laboratory Animal Center, Zhejiang Province. Our Institutional Animal Care and Use Committee approved all aspects of the study. Rabbits were randomly allocated to 1-, 3-, and 6-month survival groups (n = 6 in each group). The surgical implantation procedure was described previously [17]. Consecutive levels of rabbit IVD, including L3/4, L4/5, and L5/6, were exposed using an anterolateral retroperitoneal approach. Annular injuries were randomly allocated to three disc levels: (1) in the annular defects group: annular defects (diameter 1.8 mm and 4 mm depth) were created using a minitrephine as previously described [17,18]; (2) in the PLGA/fibrin/PBS group: annular defects were filled with a PLGA/fibrin gel plug loaded with phosphate buffered saline; (3) in the PLGA/fibrin/chABC group: annular defects were filled with a PLGA/fibrin gel plug loaded with chABC; and (4) in the intact group: the L2/3 disc served as the uninjured control. Finally, the wound was closed in layers. Following surgery, rabbits were permitted free cage activity and food and water ad libitum.

2.4 Radiographic analyses and magnetic resonance imaging (MRI)

X-rays were obtained under general anesthesia (sodium pentobarbital, 30 mg/kg) at 1, 3, and 6 months after surgery (n = 6 per time point). The comparison of the disc height index (DHI) was based on the ratio of the injured disc height to the sum of the height of the two adjacent vertebral bodies as described previously [35,36]. DHI% values are expressed as DHI normalized to the baseline preoperative measurement (postoperative DHI/baseline DHI × 100). In vivo MRI of the lumbar spine was performed using a 1.5-T MRI system (Signa, General Electric, Milwaukee, Wisconsin). Midsagittal T2-weighted images were obtained for analysis as previously described according to Pfirrmann’s classification scores [37] based on the changes in the degree and the area of the signal intensity. All measurements were performed using the picture archiving and communication system (PACS) routinely used in the local hospital.

2.5 Tissue harvesting

At each of the three time points (1, 3 and 6 months), 6 rabbits were euthanized by an intravenous sodium pentobarbital overdose for histological and immunohistological analyses. The experimental IVDs (L2/3, L3/4, L4/5, and L5/6) were removed from each lumbar spine under sterile conditions. The specimens were then dissected sagittally and divided into two symmetric portions. Using one half of each disc (with the defect side) was used for histological analysis. From the other half of the samples (with the contralateral side), the nucleus pulposus (NP) was bluntly separated from the AF and then snap-frozen in liquid nitrogen. Samples were then stored at −80°C in preparation for s-GAG analysis.

2.6 Sulfated-GAG content measurement

NP samples were isolated from discs at each level at 6 months after surgery (n = 6). The PG content was quantified using the 1,9-dimethylmethylene blue DMMB assay [38,39]. Briefly, each lyophilized sample was digested with 125 μg/mL of papain (Sangon Inc., Shanghai; PRC) in sterile PBS, 5 mM of EDTA, and 5 mM of cysteine·HCl at pH 6.8 and 60°C overnight. After complete digestion, 20 μL of the papain digest was added to 200 μL of DMMB reagent (Sigma-Aldrich, St. Louis, Missouri, USA), and absorbance was detected at 520 nm. A Quant-iT™ PicoGreen® dsDNA Assay Kit (Invitrogen™) was used to determine DNA content in the discs, according to the manufacturer’s protocol (Thermo Fisher Scientific Inc., Waltham, MA; USA). The total sGAG in the disc for each group was normalized according to the tested amount of DNA, and the sGAG/DNA ratio was measured and reported.

2.7 Histology and immunohistochemistry

The specimens were fixed in 10% formalin, decalcified in ethylenediamine tetraacetic acid (EDTA), and processed for paraffin sectioning. Tissue blocks were paraffin embedded and sectioned at a thickness of 5 μm. Sections of the IVD samples were stained with hematoxylin/eosin (HE) to observe degenerative changes or with safranin-O/fast green to determine the PG content. Alternatively, the sections were subjected to immunohistological analysis for the nerve marker Protein gene product 9.5 (PGP9.5). CS-56 immunostaining was used to identify intact CS. All the stained sections were analyzed under an optical microscope (Leica Microscope, Wetzlar, Germany). Briefly, the epitopes for PGP9.5 immunohistology were first subjected to heat-induced retrieval. The sections were then blocked with hydrogen peroxide and 25% of normal bovine serum albumin (BSA)/tris buffered saline (TBS) and incubated overnight at 4°C with the primary mouse monoclonal antibody against human PGP9.5 (diluted 1:80, Abcam, Cambridge, GB). For CS-56 staining, endogenous peroxidase activity in the sections was blocked with 0.5% of hydrogen peroxide for 1 h, and the sections were then washed three times in TBS. Then, the samples were blocked with 25% of normal BSA for 1 h (both at RT), followed by an overnight incubation at 4°C with a primary anti-CS-56 mouse antibody (diluted 1:100, Sigma, St. Louis, MO). The sections were then washed again and incubated with a biotinylated goat anti-mouse IgG antibody (1:200; Vector Laboratories) overnight. The sections were then processed with the avidin–biotin amplification method with conjugated peroxidase (Vectastain Elite ABC Kit; Vector) and visualized with diaminobenzidine (DAB; Sigma). The sections were then counterstained with hematoxylin (blue) and mounted with Aquatex for observation under light microscopy. The control immunoglobulins consistently yielded negative results. The ingrowth of immunoreactive nerve fibers in the specimens was scored using a previously described grading scale [18,40].

2.8 Statistical analysis

Results are expressed as mean values ± standard error of the mean. Statistical analysis was performed using SPSS 22.0 software (SPSS Inc., Chicago, IL, USA). Significant differences in the radiograph measurements were analyzed by repeated-measurement analysis of variance (ANOVA) and Fisher’s least significant difference (LSD) test. The effect of time after surgery was analyzed with the Kruskal–Wallis test. Mann–Whitney U tests were used to analyze the MRI score, innervation, and biochemical data. The significance level was set at p < 0.05.

  1. Ethics approval and consent to participate: All the procedures performed in this study involving animals were approved by the institutional review board and animal care committee of Tongde Hospital. The protocols were conducted in accordance with the Guidance for the Care and Use of Laboratory Animals, as formulated by the Ministry of Science and Technology of the People’s Republic of China and the “Principles of laboratory animal care” (NIH publication No. 86-23, revised 1985).

3 Results

3.1 PLGA/fibrin/ChABC morphology and ChABC release kinetics

Fabricated PLGA sponges filled with fibrin gel and the chABC coating of the PLGA/fibrin construct were assessed by SEM and are shown in Figure 1a and b, respectively. The SEM micrographs revealed interconnected micropores with a mean pore size of 350 µm. The profiles of chABC release from the PLGA/fibrin scaffold over 14 days are shown in Figure 1c. The curve of chABC from fibrin-based carriers showed a significant initial burst release during the first 5 days, and the chABC release continued steadily thereafter. Overall, the cumulative release rate of active chABC reached approximately 80% after 8 days.

Figure 1 
                  SEM micrographs of (a) the PLGA/fibrin gel sponge and (b) the composite sponge filled with fibrin gel and chABC at different magnifications; (c) in vitro cumulative release profile of chABC/PLGA/fibrin.
Figure 1

SEM micrographs of (a) the PLGA/fibrin gel sponge and (b) the composite sponge filled with fibrin gel and chABC at different magnifications; (c) in vitro cumulative release profile of chABC/PLGA/fibrin.

3.2 Radiographic and MRI assessment

Compared with the normal control group, the disc height of the injured groups showed a slow but progressive decrease that was sustained for up to 6 months. At 6 months, a narrowing disc space and bridging osteophytes were observed in the two scaffold groups (Figure 2a and b). The DHI% of the two scaffold groups and the AF defect group were significantly reduced compared with those of normal control groups at all postoperative periods (p < 0.05), whereas no significant difference was observed between the PLGA/fibrin/PBS and PLGA/fibrin/chABC groups (p > 0.05). Notably, DHI% in the defect group demonstrated a significant decrease at 1, 3, and 6 months after surgery (p < 0.01) (Figure 3).

Figure 2 
                  Representative lateral radiographs (a and b) and MRI (c and d) of the rabbit lumbar spine obtained 6 months after surgery. White arrows indicate that the injured disc received a different treatment.
Figure 2

Representative lateral radiographs (a and b) and MRI (c and d) of the rabbit lumbar spine obtained 6 months after surgery. White arrows indicate that the injured disc received a different treatment.

Figure 3 
                  Changes in disc height index after surgery. There was a slow, progressive decrease in the disc height of operated discs over the follow-up period. The DHI values in both scaffold groups and the AF defect group were significantly reduced compared with that of the control group after surgery (*,#
                     p < 0.05, vs the control group). DHI% in the defect group was significantly decreased at 1, 3, and 6 months after surgery (+
                     p < 0.01, vs the two scaffold groups).
Figure 3

Changes in disc height index after surgery. There was a slow, progressive decrease in the disc height of operated discs over the follow-up period. The DHI values in both scaffold groups and the AF defect group were significantly reduced compared with that of the control group after surgery (*,# p < 0.05, vs the control group). DHI% in the defect group was significantly decreased at 1, 3, and 6 months after surgery (+ p < 0.01, vs the two scaffold groups).

Serial MRI scans of rabbit lumbar spines revealed that the appearance of normal control discs remained relatively constant. Progressive decreases in the signal intensity in the NP area were apparent for each of the injured discs at all time points, as shown in Figure 2c and d. The MRI grades of the injured groups were progressively increased compared with those of the normal control groups at the postoperative time point (p < 0.01). Additionally, the MRI grades of the AF defect and PLGA/fibrin/chABC discs were significantly increased compared with those of the PLGA/fibrin/PBS or control discs after surgery (p < 0.05). However, no significant differences in the MRI grades were noted between the PLGA/fibrin/chABC and AF defect groups at 1 and 6 months after surgery (p > 0.05) (Figure 4).

Figure 4 
                  Changes in MRI grade after surgery. The MRI grades of the injured groups were progressively increased compared with those of the normal control groups at the postoperative time point (*p < 0.01, vs the control group). Both the MRI grades of AF defect and PLGA/fibrin/chABC discs were significantly increased compared with those of the PLGA/fibrin/PBS or control discs after surgery (#
                     p < 0.05, vs the control group or PLGA/fibrin/PBS group).
Figure 4

Changes in MRI grade after surgery. The MRI grades of the injured groups were progressively increased compared with those of the normal control groups at the postoperative time point (*p < 0.01, vs the control group). Both the MRI grades of AF defect and PLGA/fibrin/chABC discs were significantly increased compared with those of the PLGA/fibrin/PBS or control discs after surgery (# p < 0.05, vs the control group or PLGA/fibrin/PBS group).

3.3 Nuclear sulfated-GAG content

The s-GAG content of the three injured groups decreased significantly compared to that of the control group at different postoperative time points (p < 0.01). The s-GAG content of the PLGA/fibrin/chABC group exhibited a significant decrease compared with that of the PLGA/fibrin/PBS and AF defect groups after surgery (p < 0.05). The decrease in the sGAG/DNA ratio observed in the AF defect group was significantly more pronounced than that in the PLGA/fibrin/PBS group at 1 and 6 months after surgery (p < 0.05) (Figure 5).

Figure 5 
                  Changes in proteoglycan content of the NP at 6 months after surgery. The s-GAG content of the three injured groups decreased significantly compared to that of the control group at different postoperative time points (*p < 0.01, vs the control group). The PLGA/fibrin/chABC group exhibited a significant decrease in the sGAG/DNA ratio compared with the PLGA/fibrin/PBS or AF defect groups after surgery (#
                     p < 0.05 vs the PLGA/fibrin/PBS group or AF defect group). The decrease in the sGAG/DNA ratio observed in the AF defect group was significantly more pronounced than that in the PLGA/fibrin/PBS group at 1 and 6 months after surgery (+
                     p < 0.05).
Figure 5

Changes in proteoglycan content of the NP at 6 months after surgery. The s-GAG content of the three injured groups decreased significantly compared to that of the control group at different postoperative time points (*p < 0.01, vs the control group). The PLGA/fibrin/chABC group exhibited a significant decrease in the sGAG/DNA ratio compared with the PLGA/fibrin/PBS or AF defect groups after surgery (# p < 0.05 vs the PLGA/fibrin/PBS group or AF defect group). The decrease in the sGAG/DNA ratio observed in the AF defect group was significantly more pronounced than that in the PLGA/fibrin/PBS group at 1 and 6 months after surgery (+ p < 0.05).

3.4 Histologic assessment

At 6 months, HE staining demonstrated that the uninjured, healthy discs (control group) exhibited a characteristically well-organized intact AF with concentric lamellae. Upon safranin-O staining, uninjured discs exhibited a minimal disruption of the PG matrix within the annulus (Figure 6a and b). In the AF defect group, HE staining revealed the loss of the lamellated appearance of AF and the replacement of the annular defect by extensive fibrocartilaginous-like tissue that formed outside the injured sites. Some blood vessels and small fissures with a limited depth were typically distributed in the outer scar tissue and AF but did not extend to the inner AF (Figure 6c). Safranin-O staining indicated the presence of a PG-rich content in the fibrocartilaginous tissue (Figure 6d). In the PLGA/fibrin/PBS discs, HE staining revealed small residues of the PLGA scaffold with a naturally irregular form that were enfolded by newly formed tissue and well-integrated within the inner AF (Figure 6e). In addition, safranin-O staining revealed an open concave cavity with an increased loss of the PG content (Figure 6f). In contrast, advanced degeneration was observed in the PLGA/fibrin/chABC group, including NP fibrosis, disorganization of the AF, and new tissue clusters extending further into the deeper inner AF along the fissures. In addition, blood vessels and small fissures were commonly observed in the PLGA/fibrin/PBS group (Figure 6g). The PG content was markedly reduced in the injured region (Figure 6h).

Figure 6 
                  HE (a, c, e, and g) and safranin-O/fast green (b, d, f, and h) staining of the lesion site 6 months after surgery. The uninjured AF displayed a multilamellar structure rich in proteoglycans as demonstrated by strong safranin-O staining (a and b), loss of AF structural integrity, and replacement by extensive fibrocartilaginous-like tissue in the injured sites. Safranin-O staining indicated the presence of proteoglycan-rich content in the fibrocartilaginous tissue (c and d). Small irregular remnants of the PLGA scaffold were well integrated within the inner AF. Safranin-O staining revealed an open concave cavity with increased severe loss of the proteoglycan content (e and f). Advanced degeneration was observed in the PLGA/fibrin/chABC group, including NP fibrosis, disorganization of the AF, and reparative tissue extending further into the deeper inner AF along the fissures (g). The proteoglycan content was markedly reduced in the injured region (h). White arrows indicate areas of some blood vessels and small fissures.
Figure 6

HE (a, c, e, and g) and safranin-O/fast green (b, d, f, and h) staining of the lesion site 6 months after surgery. The uninjured AF displayed a multilamellar structure rich in proteoglycans as demonstrated by strong safranin-O staining (a and b), loss of AF structural integrity, and replacement by extensive fibrocartilaginous-like tissue in the injured sites. Safranin-O staining indicated the presence of proteoglycan-rich content in the fibrocartilaginous tissue (c and d). Small irregular remnants of the PLGA scaffold were well integrated within the inner AF. Safranin-O staining revealed an open concave cavity with increased severe loss of the proteoglycan content (e and f). Advanced degeneration was observed in the PLGA/fibrin/chABC group, including NP fibrosis, disorganization of the AF, and reparative tissue extending further into the deeper inner AF along the fissures (g). The proteoglycan content was markedly reduced in the injured region (h). White arrows indicate areas of some blood vessels and small fissures.

3.5 CS-56 immunolabeling for CS-GAG digestion

CS-56 immunolabeling in the lesion site indicated that CS was degraded following treatment with chABC. In the control group, the outer lamellar AF was weakly stained with the CS-56 antibody (Figure 7a). In the injured discs, scar tissue was formed on the surface of the injury site. A fuzzy CS-56-positive deposit was distributed in the scar tissue (Figure 7b). In the PLGA/fibrin/PBS group, small residues of the PLGA scaffold were surrounded by newly repaired tissue penetrating deeper toward the NP. CS-56-positive deposits were relatively negligible in the clusters of the newly regenerated tissue (Figure 7c). In contrast, in the chABC/PLGA/fibrin-treated groups, CS-56 immunolabeling was high in the region of the newly generated tissue at the lesion site, but intense reactivity was noted overall (Figure 7d).

Figure 7 
                  CS-56 immunohistochemical staining at the lesion site. (a) The outer lamellar AF was weakly stained with the CS-56 antibody in the control group. (b) In the injured discs, scar tissue formed on the surface of the injury site. A fuzzy CS-56-positive deposit was distributed in the scar tissue. (c) In the PLGA/Fibrin/PBS discs, newly formed tissue penetrated deeper toward the NP. CS-56-positive deposits were relatively negligible in the reparative tissue. (d) In the PLGA/Fibrin/chABC discs, CS-56-positive staining showed intense reactivity at the site of injury following treatment with chABC. White arrows indicate some regions of “intense” CS-56 reactivity. NP (nucleus pulposus).
Figure 7

CS-56 immunohistochemical staining at the lesion site. (a) The outer lamellar AF was weakly stained with the CS-56 antibody in the control group. (b) In the injured discs, scar tissue formed on the surface of the injury site. A fuzzy CS-56-positive deposit was distributed in the scar tissue. (c) In the PLGA/Fibrin/PBS discs, newly formed tissue penetrated deeper toward the NP. CS-56-positive deposits were relatively negligible in the reparative tissue. (d) In the PLGA/Fibrin/chABC discs, CS-56-positive staining showed intense reactivity at the site of injury following treatment with chABC. White arrows indicate some regions of “intense” CS-56 reactivity. NP (nucleus pulposus).

3.6 Innervation of the IVD

PGP9.5-positive nerves were sparsely noted in the lamellae and adjacent connective tissue of the outer AF (Figure 8a). In the AF defect group, small nerves in the vicinity of the outer AF were distributed along the fissures, but these nerves rarely invaded the deeper inner AF (Figure 8b). In the PLGA/fibrin/PBS group, nerve fibers reactive for PGP9.5 were identified within the newly generated tissue of the inner AF (Figure 8c). In contrast, more sprouting PGP9.5-immunoreative fibers in the PLGA/fibrin/chABC group were localized predominantly in the vicinity of the vascularized repair tissue and extended further into the inner AF (Figure 8d).

Figure 8 
                  PGP9.5 immunohistochemical staining at 6 months after surgery. In the intact AF, PGP9.5-positive nerves were distributed sparsely in the lamellae and adjacent connective tissue of the outer AF (a). Small nerves in the vicinity of the outer AF were distributed along the fissures but rarely invaded the deeper inner AF in the AF defect group (b). PGP9.5-positive nerve fibers were present within the newly generated tissue of inner AF (c). In the PLGA/fibrin/chABC discs, more sprouting nerve fibers were predominantly identified in the vicinity of vascularized repair tissue and extended further into the inner AF (d). White arrows indicate areas of PGP9.5-positive nerve fibers. Long white arrows indicate the radial direction (from nucleus to outer annulus). R, radial.
Figure 8

PGP9.5 immunohistochemical staining at 6 months after surgery. In the intact AF, PGP9.5-positive nerves were distributed sparsely in the lamellae and adjacent connective tissue of the outer AF (a). Small nerves in the vicinity of the outer AF were distributed along the fissures but rarely invaded the deeper inner AF in the AF defect group (b). PGP9.5-positive nerve fibers were present within the newly generated tissue of inner AF (c). In the PLGA/fibrin/chABC discs, more sprouting nerve fibers were predominantly identified in the vicinity of vascularized repair tissue and extended further into the inner AF (d). White arrows indicate areas of PGP9.5-positive nerve fibers. Long white arrows indicate the radial direction (from nucleus to outer annulus). R, radial.

The semiquantitative innervation scores were slightly increased in the three injured groups over the follow-up period. Notably, the innervation grade of the PLGA/fibrin/chABC group was significantly increased compared with that of the AF defect or PLGA/fibrin/PBS groups at 1 and 6 months after surgery (p < 0.05). The innervation score in the PLGA/fibrin/PBS group was significantly increased compared with that of the AF defect group at 1 and 6 months (p < 0.05), but no significant differences were noted among the different time points (p > 0.05; Figure 9).

Figure 9 
                  Changes in innervation scores at 6 months after surgery. The innervation score of the PLGA/fibrin/chABC group was significantly increased compared with that of the AF defect or PLGA/fibrin/PBS groups at 1 and 6 months after surgery (1 and 6 months, *p < 0.01). The innervation score of the PLGA/fibrin/PBS discs was significantly increased compared with that of the AF injured discs at 1 and 6 months (†
                     p < 0.05 vs the AF defect group).
Figure 9

Changes in innervation scores at 6 months after surgery. The innervation score of the PLGA/fibrin/chABC group was significantly increased compared with that of the AF defect or PLGA/fibrin/PBS groups at 1 and 6 months after surgery (1 and 6 months, *p < 0.01). The innervation score of the PLGA/fibrin/PBS discs was significantly increased compared with that of the AF injured discs at 1 and 6 months ( p < 0.05 vs the AF defect group).

4 Discussion

In our previous study [17,18,41], an AF injury model was established by creating a defect, which resulted in reproducible and degenerative MRI, radiograph, and histologic changes. Interestingly, degenerative changes occurred that were mainly characterized by regenerated tissue formation and the ingrowth of nerves and blood vessels into the disrupted disc tissue. More recent studies have demonstrated that aggrecan derived from both the AF and NP has inhibitory effects on nerve ingrowth into the IVD [13,18,42]. Thus, this study was conducted to evaluate the effects of degrading CS (CS)-PG on neurite ingrowth. Specifically, PLGA/fibrin gel with/without chABC was anchored in the annular defects during the degenerative process of the IVD. The results showed that fibrin-based PLGA loaded with chABC integrated well with the native annulus tissue and provided a local source for the sustained release of the active chABC. Importantly, more nerve fibers were orientated conductively to grow into the deeper PG-depleted regions of the injured annulus following chABC enzymatic deglycosylation. To the best of our knowledge, few reports have employed in vivo testing to assess the effects of PG-depleted sites in the annular wall on nerve ingrowth in a rabbit disc degeneration model.

It is well-known that chABC more specifically degrades the CS side chains of PGs [24,43]. Indeed, the use of a proteolytic enzyme for the induction of mild disc degeneration has been successfully achieved in both rabbits [44,45] and large sheep [19,33]. Furthermore, the reduction in PGs using a fibrin-based delivery system to achieve enzymatic digestion depends on its dosage and activity [25,46,47,48]. Based on the promising preliminary results, a dose of 0.1 U/mL of chABC incorporated into PLGA/fibrin gel was tested in the present study. Consistent with previous studies [23,34], the rate and duration of in vitro chABC release from the fibrin-based polymer were characterized by an initial burst (70–80%) that was sustained for up to 8 days. Thereafter, a small amount of chABC was still released from the PLGA/fibrin gel, achieving a delayed-release effect, as indicated in Figure 1c. In addition, CS-56 immunolabeling of the lesion site revealed that CS was locally degraded following treatment with chABC. Taken together, the results indicated that the administration of fibrin-based polymeric carriers at the injured site achieved sustained, localized release of chABC, and in vivo chABC activity was maintained after surgery, as demonstrated by the continuous significant decrease in the s-GAG content of the PLGA/fibrin/chABC discs compared with the two control discs. Conceivably, the present study offers a faithful analysis of the actual process of GAG reduction and its activity because it is very unlikely that the in vitro release kinetics would accurately reflect the more relevant in vivo release kinetics at the lesion site.

We observed PGP9.5-immunoreactive fibers in the defect group exclusively in the superficial area of scar tissue, and similar findings have been described in previous studies [18]. However, increased nerve ingrowth into PG-depleted discs sealed with PLGA/fibrin/chABC yielded a significantly increased innervation score compared to PLGA/fibrin/PBS or empty injured discs. In addition, this finding indicated that advanced degeneration with small fissures and vascularization occurs significantly more often and penetrates deeper into the PLGA/fibrin/chABC discs in the early postoperative period (Figure 6g). Eventually, the degenerative area appeared increasingly innervated in the chABC-treated discs following the reduction in PG. Furthermore, recent studies have shown that nerves and blood vessels grow into those regions of human lumbar discs that are damaged and depleted of PGs, supporting the present results from rabbit discs [12,15,49]. Several factors potentially influence nerve ingrowth into the injured IVD: (i) more extensive PG reduction in the depleted annulus regions is conducive to neural and vascular ingrowth following the chABC enzymatic deglycosylation; (ii) loss of structural integrity results in the early leakage of NP with the depletion of PGs, thus increasing the attraction of ingrowing nerves; (iii) the disruption of the tight collagen network that entraps the PGs eventually produces more new tissue (e.g., periannular innervated and vascularized granulation tissue) ingrowth along the local AF deficiency; and (iv) ingrowth is potentially inhibited by scar tissue due to the reduced PG loss compared with the newly regenerated tissue in the AF injured track. A further limitation of using enzyme degradation in an in vivo model involves the possibility that local undegraded islands of CSPGs may exist and are resorbed by cells synthesizing and replenishing the lost GAGs. Thus, these areas do not completely mimic a natural ‘degenerative’ process. Further studies are necessary to evaluate neural ingrowth in the degenerated IVD upon application of a control enzyme (e.g., hyaluronidase or matrix metalloproteinases) as a treatment option in a large animal model with longer follow-up times and dose-related deglycosylation.

5 Conclusion

Our results indicate that the chABC-based PLGA/fibrin gel showed promising results, achieving biointegration with native annulus tissue and providing a local source for the sustained release of active chABC. Disc PG-mediated inhibition of nerve and blood vessel ingrowth was abrogated by deglycosylation in an annular-injured rabbit disc degeneration model. Moreover, PG reduction may play a direct role in nerve and blood vessel ingrowth into degenerated IVDs. Intact GAGs, such as CS, secreted by disc cells are potential candidates that could be useful to reduce neural and vascular ingrowth associated with discogenic pain in degenerated IVDs.

Abbreviations

AF

annulus fibrosus

ANOVA

analysis of variance

BSA

bovine serum albumin

chABC

chondroitinase ABCDMMB

CS

chondroitin sulfate

DHI

disc height index

DMMB

1,9-dimethylmethylene blue

EDTA

ethylenediamine tetraacetic acid

GAG

glycosaminoglycan

HE

hematoxylin/eosin

IVD

intervertebral disc

LSD

least significant difference

MRI

magnetic resonance imaging

NP

nucleus pulposus

PBS

phosphate-buffered saline

PLGA

poly(lactic-co-glycolic acid)

PG

proteoglycan

PGP9.5

protein gene product 9.5

PACS

picture archiving and communication system

TBS

tris-buffered saline

Acknowledgement

The authors gratefully acknowledge the financial support by Zhejiang Provincial Health commission of Foundation (No. 2020362988; 2020PY003; 2019KY050) and Zhejiang Provincial of Administration of Traditional Chinese Medicine of Foundation (2019ZA026; 2019ZB027 and 2019ZB029).

  1. Funding information: Zhejiang Provincial Health commission of Foundation (No. 2020362988; 2020PY003; 2019KY050) and Zhejiang Provincial of Administration of Traditional Chinese Medicine of Foundation (2019ZA026; 2019ZB027 and 2019ZB029).

  2. Author contributions: All the authors were equally involved in conception and design of the study. Long X, Jian W and Weixin X performed the surgical treatment and animal care. Fang Y and Xinwei X performed the acquisition, analysis and interpretation of the data. Long X, Fang Y, Shunwu F and Yang Y completed the drafting of the manuscript and revised it critically for important intellectual content. All the authors read and approved the final manuscript.

  3. Conflict of interest: The authors declare that they have no conflicts of interest.

  4. Consent for publication: Data from individuals are not applicable in the study. This manuscript has never been published. We have a registered preprint online, found in a website accessible from: “https://www.researchsquare.com/article/rs19538/v1,” DOI: 10.21203/rs.3.rs19538/v1.

  5. Data availability statement: The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

[1] Patrick N, Emanski E, Knaub MA. Acute and chronic low back pain. Med Clin North Am. 2014;98(4):777–89, xii.10.1016/j.mcna.2014.03.005Search in Google Scholar

[2] Hoy D, Brooks P, Blyth F, Buchbinder R. The epidemiology of low back pain. Best Pract Res Clin Rheumatol. 2010;24(6):769–81.10.1016/j.berh.2010.10.002Search in Google Scholar

[3] Hughes SP, Freemont AJ, Hukins DW, McGregor AH, Roberts S. The pathogenesis of degeneration of the intervertebral disc and emerging therapies in the management of back pain. J Bone Jt Surg Br Volume. 2012;94(10):1298–304.10.1302/0301-620X.94B10.28986Search in Google Scholar

[4] Adams MA, Roughley PJ. What is intervertebral disc degeneration, and what causes it? Spine. 2006;31(18):2151–61.10.1097/01.brs.0000231761.73859.2cSearch in Google Scholar

[5] Colombier P, Clouet J, Hamel O, Lescaudron L, Guicheux J. The lumbar intervertebral disc: from embryonic development to degeneration. Joint, Bone, SpineRev du Rhumatisme. 2014;81(2):125–9.10.1016/j.jbspin.2013.07.012Search in Google Scholar

[6] Freemont AJ. The cellular pathobiology of the degenerate intervertebral disc and discogenic back pain. Rheumatol (Oxford, Engl). 2009;48(1):5–10.10.1093/rheumatology/ken396Search in Google Scholar

[7] Stefanakis M, Al-Abbasi M, Harding I, Pollintine P, Dolan P, Tarlton J, et al. Annulus fissures are mechanically and chemically conducive to the ingrowth of nerves and blood vessels. Spine. 2012;37(22):1883–91.10.1097/BRS.0b013e318263ba59Search in Google Scholar

[8] Aoki Y, Nakajima A, Ohtori S, Takahashi H, Watanabe F, Sonobe M, et al. Increase of nerve growth factor levels in the human herniated intervertebral disc: can annular rupture trigger discogenic back pain? Arthritis Res Ther. 2014;16(4):R159.10.1186/ar4674Search in Google Scholar

[9] Freemont AJ, Peacock TE, Goupille P, Hoyland JA, O’Brien J, Jayson MI. Nerve ingrowth into diseased intervertebral disc in chronic back pain. Lancet (London, Engl). 1997;350(9072):178–81.10.1016/S0140-6736(97)02135-1Search in Google Scholar

[10] Ohtori S, Inoue G, Miyagi M, Takahashi K. Pathomechanisms of discogenic low back pain in humans and animal models. Spine J: Off J North Am Spine Soc. 2015;15(6):1347–55.10.1016/j.spinee.2013.07.490Search in Google Scholar PubMed

[11] Rea WJ, Kapur S, Mutagi H. Intervertebral disc as a source of pain. Continuing Edu Anaes Crit Care Pain Adv. 2012;12(6):279–82.10.1093/bjaceaccp/mks028Search in Google Scholar

[12] Binch AL, Cole AA, Breakwell LM, Michael AL, Chiverton N, Creemers LB, et al. Nerves are more abundant than blood vessels in the degenerate human intervertebral disc. Arthritis Res Ther. 2015;17:370.10.1186/s13075-015-0889-6Search in Google Scholar PubMed PubMed Central

[13] Johnson WE, Caterson B, Eisenstein SM, Roberts S. Human intervertebral disc aggrecan inhibits endothelial cell adhesion and cell migration in vitro. Spine. 2005;30(10):1139–47.10.1097/01.brs.0000162624.95262.73Search in Google Scholar PubMed

[14] Johnson WE, Caterson B, Eisenstein SM, Hynds DL, Snow DM, Roberts S. Human intervertebral disc aggrecan inhibits nerve growth in vitro. Arthritis Rheumatism. 2002;46(10):2658–64.10.1002/art.10585Search in Google Scholar PubMed

[15] Purmessur D, Cornejo MC, Cho SK, Roughley PJ, Linhardt RJ, Hecht AC, et al. Intact glycosaminoglycans from intervertebral disc-derived notochordal cell-conditioned media inhibit neurite growth while maintaining neuronal cell viability. Spine J: Off J North Am Spine Soc. 2015;15(5):1060–9.10.1016/j.spinee.2015.02.003Search in Google Scholar PubMed PubMed Central

[16] Melrose J, Roberts S, Smith S, Menage J, Ghosh P. Increased nerve and blood vessel ingrowth associated with proteoglycan depletion in an ovine anular lesion model of experimental disc degeneration. Spine. 2002;27(12):1278–85.10.1097/00007632-200206150-00007Search in Google Scholar PubMed

[17] Xin L, Zhang C, Zhong F, Fan S, Wang W, Wang Z. Minimal invasive annulotomy for induction of disc degeneration and implantation of poly(lactic-co-glycolic acid) (PLGA) plugs for annular repair in a rabbit model. Eur J Med Res. 2016;21(1):7.10.1186/s40001-016-0202-4Search in Google Scholar PubMed PubMed Central

[18] Xin L, Xu W, Yu L, Fan S, Wang W, Yu F, et al. Effects of annulus defects and implantation of poly(lactic-co-glycolic acid) (PLGA)/fibrin gel scaffolds on nerves ingrowth in a rabbit model of annular injury disc degeneration. J Ortho Surg Res. 2017;12(1):73.10.1186/s13018-017-0572-5Search in Google Scholar PubMed PubMed Central

[19] Gullbrand SE, Malhotra NR, Schaer TP, Zawacki Z, Martin JT, Bendigo JR, et al. A large animal model that recapitulates the spectrum of human intervertebral disc degeneration. Osteoarthr Cartil. 2017;25(1):146–56.10.1016/j.joca.2016.08.006Search in Google Scholar PubMed PubMed Central

[20] Shen Y, Tenney AP, Busch SA, Horn KP, Cuascut FX, Liu K, et al. PTPsigma is a receptor for chondroitin sulfate proteoglycan, an inhibitor of neural regeneration. Sci (N York, NY). 2009;326(5952):592–6.10.1126/science.1178310Search in Google Scholar PubMed PubMed Central

[21] Jefferson SC, Tester NJ, Howland DR. Chondroitinase ABC promotes recovery of adaptive limb movements and enhances axonal growth caudal to a spinal hemisection. J Neurosci. 2011;31(15):5710–20.10.1523/JNEUROSCI.4459-10.2011Search in Google Scholar PubMed PubMed Central

[22] Zukor K, He Z. Regenerative medicine: drawing breath after spinal injury. Nature. 2011;475(7355):177–8.10.1038/475178aSearch in Google Scholar PubMed

[23] Lee JC, Min HJ, Lee S, Seong SC, Lee MC. Effect of chondroitinase ABC on adhesion and behavior of synovial membrane-derived mesenchymal stem cells in rabbit partial-thickness chondral defects. J Ortho Res. 2013;31(8):1293–301.10.1002/jor.22353Search in Google Scholar PubMed

[24] Muir E, De Winter F, Verhaagen J, Fawcett J. Recent advances in the therapeutic uses of chondroitinase ABC. Exp Neurol. 2019;321:113032.10.1016/j.expneurol.2019.113032Search in Google Scholar PubMed

[25] Wang W, Li B, Yang J, Xin L, Li Y, Yin H, et al. The restoration of full-thickness cartilage defects with BMSCs and TGF-beta 1 loaded PLGA/fibrin gel constructs. Biomaterials. 2010;31(34):8964–73.10.1016/j.biomaterials.2010.08.018Search in Google Scholar PubMed

[26] Li XK, Cai SX, Liu B, Xu ZL, Dai XZ, Ma KW, et al. Characteristics of PLGA-gelatin complex as potential artificial nerve scaffold. Colloids Surf B Biointerfaces. 2007;57(2):198–203.10.1016/j.colsurfb.2007.02.010Search in Google Scholar PubMed

[27] Xiong Y, Zeng YS, Zeng CG, Du BL, He LM, Quan DP, et al. Synaptic transmission of neural stem cells seeded in 3-dimensional PLGA scaffolds. Biomaterials. 2009;30(22):3711–22.10.1016/j.biomaterials.2009.03.046Search in Google Scholar PubMed

[28] Lee DY, Choi BH, Park JH, Zhu SJ, Kim BY, Huh JY, et al. Nerve regeneration with the use of a poly(l-lactide-co-glycolic acid)-coated collagen tube filled with collagen gel. J Cranio-Maxillo-Facial Surg. 2006;34(1):50–6.10.1016/j.jcms.2005.07.011Search in Google Scholar PubMed

[29] Han J, Lazarovici P, Pomerantz C, Chen X, Wei Y, Lelkes PI. Co-electrospun blends of PLGA, gelatin, and elastin as potential nonthrombogenic scaffolds for vascular tissue engineering. Biomacromolecules. 2011;12(2):399–408.10.1021/bm101149rSearch in Google Scholar PubMed

[30] Pan Z, Ding J. Poly(lactide-co-glycolide) porous scaffolds for tissue engineering and regenerative medicine. Interface Focus. 2012;2(3):366–77.10.1098/rsfs.2011.0123Search in Google Scholar PubMed PubMed Central

[31] Zhou QGY, Gao C. Microstructure and mechanical properties of poly(l-lactide) scaffolds fabricated by gelatin particle leaching method. J Appl Polym Sci. 2005;98(3):1373–9.10.1002/app.22289Search in Google Scholar

[32] Dresdale A, Rose EA, Jeevanandam V, Reemtsma K, Bowman FO, Malm JR. Preparation of fibrin glue from single-donor fresh-frozen plasma. Surgery. 1985;97(6):750–5.Search in Google Scholar

[33] Hoogendoorn RJ, Wuisman PI, Smit TH, Everts VE, Helder MN. Experimental intervertebral disc degeneration induced by chondroitinase ABC in the goat. Spine. 2007;32(17):1816–25.10.1097/BRS.0b013e31811ebac5Search in Google Scholar

[34] Hyatt AJ, Wang D, Kwok JC, Fawcett JW, Martin KR. Controlled release of chondroitinase ABC from fibrin gel reduces the level of inhibitory glycosaminoglycan chains in lesioned spinal cord. J Controlled Rel. 2010;147(1):24–9.10.1016/j.jconrel.2010.06.026Search in Google Scholar

[35] Masuda K, Aota Y, Muehleman C, Imai Y, Okuma M, Thonar EJ, et al. A novel rabbit model of mild, reproducible disc degeneration by an anulus needle puncture: correlation between the degree of disc injury and radiological and histological appearances of disc degeneration. Spine. 2005;30(1):5–14.10.1097/01.brs.0000148152.04401.20Search in Google Scholar

[36] Masuda K, Imai Y, Okuma M, Muehleman C, Nakagawa K, Akeda K, et al. Osteogenic protein-1 injection into a degenerated disc induces the restoration of disc height and structural changes in the rabbit anular puncture model. Spine. 2006;31(7):742–54.10.1097/01.brs.0000206358.66412.7bSearch in Google Scholar

[37] Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine. 2001;26(17):1873–8.10.1097/00007632-200109010-00011Search in Google Scholar

[38] Farndale RW, Buttle DJ, Barrett AJ. Improved quantitation and discrimination of sulphated glycosaminoglycans by use of dimethylmethylene blue. Biochim et Biophys Acta. 1986;883(2):173–7.10.1016/0304-4165(86)90306-5Search in Google Scholar

[39] Lao L, Tan H, Wang Y, Gao C. Chitosan modified poly(L-lactide) microspheres as cell microcarriers for cartilage tissue engineering. Colloids Surf B Biointerfaces. 2008;66(2):218–25.10.1016/j.colsurfb.2008.06.014Search in Google Scholar PubMed

[40] Aoki Y, Akeda K, An H, Muehleman C, Takahashi K, Moriya H, et al. Nerve fiber ingrowth into scar tissue formed following nucleus pulposus extrusion in the rabbit anular-puncture disc degeneration model: effects of depth of puncture. Spine. 2006;31(21):E774–80.10.1097/01.brs.0000238681.71537.41Search in Google Scholar PubMed

[41] Xin L, Han GC, Zhao FD, Zhao X, Li G, Fan SW. In vivo study of innervation of degenerative intervertebral discs in rabbit anular-injury model. Zhejiang da xue xue bao Yi xue ban = J Zhejiang Univ Med Sci. 2009;38(5):485–92.Search in Google Scholar

[42] Johnson WE, Sivan S, Wright KT, Eisenstein SM, Maroudas A, Roberts S. Human intervertebral disc cells promote nerve growth over substrata of human intervertebral disc aggrecan. Spine. 2006;31(11):1187–93.10.1097/01.brs.0000217669.04903.61Search in Google Scholar PubMed

[43] Katsuhiko Ishibashi HI. Chemonucleolysis with chondroitin sulfate ABC endolyase as a novel minimally invasive treatment for patients with lumbar intervertebral disc herniation. J Spine Surg. 2019;5:S115–21.10.21037/jss.2019.04.24Search in Google Scholar

[44] Norcross JP, Lester GE, Weinhold P, Dahners LE. An in vivo model of degenerative disc disease. J Ortho Res. 2003;21(1):183–8.10.1016/S0736-0266(02)00098-0Search in Google Scholar

[45] Imai Y, Okuma M, An HS, Nakagawa K, Yamada M, Muehleman C, et al. Restoration of disc height loss by recombinant human osteogenic protein-1 injection into intervertebral discs undergoing degeneration induced by an intradiscal injection of chondroitinase ABC. Spine. 2007;32(11):1197–205.10.1097/BRS.0b013e3180574d26Search in Google Scholar PubMed

[46] Varden LJ, Nguyen DT, Michalek AJ. Slow depressurization following intradiscal injection leads to injectate leakage in a large animal model. JOR Spine. 2019;2(3):e1061.10.1002/jsp2.1061Search in Google Scholar PubMed PubMed Central

[47] Fusellier M, Colombier P, Lesoeur J, Samy Y, Madec S, Gauthier O, et al. Longitudinal comparison of enzyme-and laser-treated intervertebral disc by MRI, X-ray, and histological analyses reveals discrepancies in the progression of disc degeneration: a rabbit study. BioMed Res Int. 2016;2016:1–12.10.1155/2016/5498271Search in Google Scholar PubMed PubMed Central

[48] O’Connell GD, Nims RJ, Green J, Cigan AD, Ateshian GA, Hung CT. Time and dose-dependent effects of chondroitinase ABC on growth of engineered cartilage. Eur Cell Mater. 2014;27:312–20.10.22203/eCM.v027a22Search in Google Scholar

[49] Lama P, Le Maitre CL, Harding IJ, Dolan P, Adams MA. Nerves and blood vessels in degenerated intervertebral discs are confined to physically disrupted tissue. J Anat. 2018;233(1):86–97.10.1111/joa.12817Search in Google Scholar PubMed PubMed Central

Received: 2021-03-03
Revised: 2021-07-30
Accepted: 2021-08-23
Published Online: 2021-11-02

© 2021 Long Xin et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 8.12.2023 from https://www.degruyter.com/document/doi/10.1515/med-2021-0363/html
Scroll to top button