Molecular profile of major growth factors in lumbar intervertebral disc herniation: Correlation with patient clinical and epidemiological characteristics
Abstract
The involvement of growth factors (GFs) in the pathogenesis of lumbar intervertebral disc (ID) herniation and the spontaneous resorption of herniated ID fragments remains only partially elucidated. A simultaneous assessment of the transcript levels of numerous GFs and their association with clinical and epidemiological profiles of human ID herniation would provide valuable insight into the biology and clinical course of the disease. In the present study, we examined simultaneously the transcript levels of vascular endothelial growth factor (VEGF), transforming growth factor β1 (TGF-β1), basic fibroblast growth factor 2 (bFGF2), platelet derived growth factor (PDGF) isoforms and receptors, epidermal growth factor (EGF) and insulin growth factor-1 (IGF-1) in herniated and control ID specimens and investigated their correlation with the clinicopathological profiles of patients suffering from symptomatic lumbar ID herniation. GF mRNA expression levels were determined by RT-qPCR in 63 surgical specimens from lumbar herniated discs and 10 control ID specimens. Multiple positive correlations were observed between the transcript levels of the GFs examined in the ID herniation group. VEGF mRNA expression was significantly increased in the protruding compared with the extruded discs. Intense and acute pain significantly upregulated the PDGF transcript levels. Significant negative correlations were observed between the patient body mass index and the transcript levels of VEGF and PDGF receptors. Our findings support the hypothesis of the involvement of GFs in the natural history of ID herniation. GFs synergistically act in herniated IDs. Increased VEGF expression possibly induces the neovascularization process in the earliest stages of ID herniation. PDGF-C and -D play a role in the acute phase of radiculopathy in a metabolic response for tissue healing. A molecular effect, in addition to the biomechanical effect of obesity in the pathogenesis of ID herniation is also implied.
Introduction
Intervertebral disc (ID) herniation is a common, yet poorly understood disorder. Behavioral, hereditary and environmental factors cumulatively affect the nutrition and homeostasis of the disc, ultimately resulting in its structural breakdown and in the occurrence of herniation. The incidence of symptomatic lumbar ID herniation in the United States is estimated at 1–2%, for which approximately 200,000 discectomies are performed annually (1). The direct costs for conservative or surgical treatment and he indirect costs deriving from the loss of productivity are enormous (2).
However, ID herniation is considered to present a favorable natural history. Sequential magnetic resonance imaging (MRI) studies have documented that the majority of disc herniations diminish in size (3). Saal et al reported that almost 80% of disc herniations treated non-operatively decreased by >50% in follow-up examinations (4). In addition, a positive correlation has been found between the regression of disc herniation and the remission of symptoms (3,4). Various molecules, such as pro-inflammatory cytokines, growth factors (GFs) and matrix metalloproteinases (MMPs) have been implicated in the mechanisms that regulate ID degradation and herniation, as well as the spontaneous resorption of herniated disc fragments (5).
GFs are polypeptides involved in the modulation of cell growth and differentiation that exert their biologic effects through endocrine, paracrine and autocrine pathways (6). Over 18 different GFs have been identified, including vascular endothelial growth factor (VEGF), transforming growth factor β1 (TGF-β1), basic fibroblast growth factor (bFGF)-2, platelet derived growth factor (PDGF), epidermal growth factor (EGF) and insulin growth factor-1 (IGF-1). The PDGF family in particular is formed by the PDGF-A, -B, -C and -D isoforms that are encoded by four genes located at different chromosomes and may form homodimers, as well as heterodimers. These molecules signal through two cell-surface receptors, PDGF-Rα and β (7).
To date, the expression of various GFs in the degenerate and herniated ID has been investigated in numerous studies. The increased protein expression of IGF-1, FGF-2, VEGF, PDGF and TGF-β has been found by means of immunohistochemistry (8–11). The transcript levels of TGF-β and VEGF have also been examined in human ID herniation specimens (12). Consequently, preliminary studies have examined the effect of administering different GFs on ID physiology and degeneration in vitro and in vivo (13,14). However, a simultaneous assessment of the molecular profile of numerous GFs in human herniated IDs has not been performed to date, at least to the best of our knowledge. In addition, clinical and epidemiological determinants, such as the extent of disc herniation, the duration and intensity of symptoms, age and employment may considerably affect the pathophysiology of ID degeneration and the potential for regression of herniated ID fragments. The association between the clinical, epidemiological and molecular profile of ID herniation would potentially provide valuable insight into the biology and clinical course of the disorder. Therefore, the purpose of the present study was to simultaneously examine the molecular profiles of numerous GFs by means of reverse transcriprtion-quantitative polymerase chain reaction (RT-qPCR), and to investigate potential correlations with the clinical and epidemiological characteristics of patients suffering from lumbar ID herniation.
Materials and methods
Tissue samples
ID tissue samples were collected from a total of 63 patients who underwent posterior open discectomy for lumbar ID herniation in our department during a recruitment period of one year (2008). The patients reported symptoms of radiculopathy prior to surgery. Patients suffering from spinal stenosis, spondylolisthesis, scoliosis or systemic inflammatory disorders were excluded from the study. After excision, an experienced pathologist, in addition to the clinicians present, examined the tissue samples and identified the pulposus sections in order to limit the percentage of end-plates or annulus. Samples were then immediately refrigerated at −80°C until RNA extraction. Ten cadaveric ID tissue samples were also obtained within 12 h of the patient succumbing to the disease, referring to various diseases other than ID herniation. The samples were used as controls of normal growth factor mRNA expression. Donors were aged between 20 and 46 years and were known to have a negative history of lower back pain, spinal trauma or systemic inflammatory disease. All the patients or relatives signed an informed consent form approved by the University Ethics Committee to participate in the present study. The study design was in accordance with the Declaration of Helsinki Guidelines.
Prior to surgery, the patients routinely underwent an MRI examination for the classification of ID herniation. The recommendations of the combined task force from the North American Spine Society, the American Society of Spine Radiology and the American Society of Neuroradiology based on the imaging characterization of Milette (15) were used to determine the type of ID herniation. Pre-operative pain intensity was measured subjectively using the visual analogue scale score (VAS). Patient employment was defined as heavy if it included manual labour with repetitive loading of the spine. A summary of patient and control group baseline data is presented in Table I. There were no significant differences in terms of age, gender, body mass index (BMI), employment and smoking habits between the control and ID herniation groups (minimum P=0.24).
RNA extraction and RT-qPCR
Total RNA was isolated from fresh tissue and homogenized with a power homogenizer using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). RNA concentration and purity were determined on a UV spectrophotometer (Hitachi Instruments Inc., San Jose, CA, USA) by absorbance measurements (260 nm absorbance and 260/280 nm absorbance ratio). RNA integrity was examined by 1% agarose gel electrophoresis and ethidium bromide staining.
Reverse transcription reactions for the preparation of first-strand cDNA from 2 µg of total RNA were performed using the AffinityScript™ Multi Temperature cDNA synthesis kit (Stratagene, La Jolla, CA, USA). Random hexamers were used as amplification primers. Quantitative PCR (qPCR) reactions were performed using the Mx3000P real-time PCR system with SYBR®-Green I Master Mix (both from Stratagene). Data were collected and analyzed using the Mx3000P real-time PCR software version 2.00, Build 215 Schema 60 (Stratagene). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal control to normalize the mRNA expression levels of the GFs examined. The primer pair sequences used are listed in Table II. PCR products were analyzed by electrophoresis on 2% agarose gels, stained with ethidium bromide and photographed on a UV light transilluminator. GF transcription levels were calculated using the following formula: normalized sample or control = (1 + EGF)−∆∆CqGF/(1 + EGAPDH)−∆∆CqGAPDH. Procedures were repeated with a cDNA template, synthesized three times from the same RNA. The mRNA levels of each sample for each gene tested represent the mean value of data acquired from three independent RT-qPCR experiments. The reproducibility of the RT-qPCR results for the same samples was 99%.
Statistical analysis
The one-sample Kolmogorov-Smirnov test was used to assess the normality of the distribution of the mRNA expression values for the genes studied. Accordingly, the mRNA expression levels of the GFs examined in the control and herniated groups were compared using parametric and non-parametric procedures, respectively. Differences in GF expression between herniated and control IDs and between subgroups of different clinical characteristics were examined using parametric (one-way and repeated measures ANOVA) and non-parametric (Kruskal-Wallis and Mann-Whitney) tests accordingly. Pairwise comparisons were performed to further explore the interaction. Categorical variables were compared with the use of χ2 test. The Spearman rank correlation test (non-parametric) was also employed to examine pair-wise correlations in growth factor mRNA expression in control and herniated discs. The significance level was set at P<0.05. Statistical calculations were performed using SPSS software, version 15 (SPSS, Inc., Chicago, IL, USA).
Results
GF mRNA levels in normal and herniated discs
The transcript levels of the GFs examined were not significant between the control and ID herniation groups (minimum P=0.26, Mann-Whitney test, Table III). With the exception of EGF that was not expressed in neither of the 2 groups, all the GFs examined were expressed in the two study groups. Specifically, the mRNA expression status (expressed or not) of the genes examined ranged from 40 to 100% in the controls and 35 to 100% in the ID herniation groups (Table III). PDGF-C and -D mRNA expression was observed in all the herniated IDs, while IGF was the least commonly expressed in the same group (34.9%).
Pair-wise GF mRNA co-expression analysis
Pair-wise correlations between the genes examined in the ID herniation groups and the controls are presented in Tables IV and andV,V, respectively (Spearman correlations). Compared with the controls, an increased number of positive correlations was observed in the ID herniation group between most of the genes examined.
Correlation of transcript levels with clinicopathological characteristics
Age and gender
Age did not significantly affect the mRNA expression levels of the genes examined (minimum P=0.07, Kruskal-Wallis test). Male patients exhibited increased mRNA expression levels of PDGF-C and -D compared with female patients (P=0.007 and 0.024, respectively, Mann-Whitney test) (absolute values not shown).
BMI
Significant negative correlations were observed between the patient BMI values and the transcript levels of VEGF, PDGFRA and PDGFRB (P=0.036, 0.037 and 0.02, respectively, Spearman correlation).
Duration of symptoms
The transcript levels of PDGF-C and -D were significantly affected by the duration of patient symptoms (P=0.02 and 0.05, respectively, Kruskal-Wallis test). Specifically, patients suffering from acute pain (<3 months) exhibited an increased mRNA expression of PDGF-C compared with those suffering for 3–10 months and for >10 months (P=0.04 and 0.01, respectively, Mann-Whitney test). Similarly, the PDGF-D expression levels were also significantly higher in patients suffering from acute (<3 months) compared with chronic pain (P=0.017, Mann-Whitney test, Fig. 1).
Pain intensity
The intensity of patient pain significantly affected the expression levels of PDGF-B, -C and -D (P=0.044, 0.024 and 0.05, Kruskal-Wallis test). Specifically, patients experiencing moderate pain (VAS 5–7) exhibited an increased mRNA expression of PDGF-B, -C and -D compared with those with low intensity pain (VAS 0–4) (P=0.01, 0.018 and 0.05, respectively, Mann-Whitney test). In addition, patients with intense pain (VAS 8–10) also exhibited a significantly increased expression of PDGF-C and -D compared with those with low intensity pain (P=0.026 and 0.03, respectively, Mann-Whitney test, Fig. 2).
Type of ID herniation
The transcript levels of VEGF significantly correlated with the type of ID herniation (P=0.023, Kruskal-Wallis test). Specifically, ID herniation samples with protrusion exhibited an increased mRNA expression of VEGF compared with the extruded IDs (P=0.019, Mann-Whitney test, Fig. 3).
Employment
The transcript levels of the genes examined did not differ significantly between patients with either light or heavy employment (minimum P=0.15, Mann-Whitney test) (absolute values not shown).
Smoking
The mRNA expression levels of the GFs examined did not differ significantly between smoking and non-smoking patients (minimum P=0.087, Mann-Whitney test). Accordingly, the amount of cigarette smoking did not significantly affect the transcript levels of any of the genes studied (minimum P=0.29, Kruskal-Wallis test) (absolute values not shown).
Discussion
The current study examined the molecular profiles of numerous GFs by means of RT-qPCR in surgically resected lumbar ID herniation specimens and matched controls. Correlations of GF transcript levels with the patient clinical and epidemiological characteristics were also explored. Multiple positive correlations between the GFs were identified in the group of herniated discs. The VEGF transcript levels were significantly higher in protruding compared with the extruded discs. The duration and intensity of patient symptoms significantly affected the PDGF transcript levels. BMI exhibited significant negative correlations with the transcript levels of VEGF and both PDGF-receptor isotypes (data not shown).
There have been limited reports in the literature examining the molecular profiles of GFs in ID herniation. In lumbar herniated IDs, Ahn et al detected TGF-β mRNA expression in 5 out of 10 specimens examined (16), whereas Kontinnen et al found a decreased mRNA expression of the same factor (12). Control groups were not included in these studies. Kontinnen et al also failed to identify the mRNA or protein expression of EGF, which they attributed to the reduced anabolic potential of the disc. The findings of the present study did not identify any differences in the mRNA expression of GFs in a substantial number of herniated and control ID specimens. By contrast, immunohistochemistry, has produced more variable results. Specchia et al reported stronger staining for IGF-1 in chondrocytes of herniated compared with control discs (10). Similarly, Tolonen et al observed significantly stronger staining for PDGF and TGF-β in human herniated ID tissue (8,9). Given that no differences in GF expression are evident at the mRNA level compared with the controls, these findings suggest that post-transcriptional modifications may occur in the setting of ID herniation.
The transcript levels of the GFs examined significantly correlated with each other in herniated IDs. In vitro experiments have demonstrated a synergistic effect of VEGF and TGF-β1 in the induction of angiogenesis (17). In addition, PDGF has been found to upregulate VEGF in mural cells, thereby further enhancing angiogenesis (18). IGF-I has also been shown to induce VEGF mRNA and protein expression in cultured colorectal carcinoma cells (19). Most of these GFs are also under the regulatory effect of inflammatory cytokines, such as interleukin-1 (IL-1) and tumor necrosis factor-α (TNF-α). Both TNF-α and IL-1 have been shown to interact with numerous GFs in human degenerated discs and trigger a sequence of macrophage and MMP activation that potentiates the spontaneous resorption of herniated ID fragments (20–22). Multiple positive correlations between numerous MMPs have also been demonstrated in human ID herniation (5). Our findings, therefore, confirm a synergistic activation of GFs in ID herniation as part of a metabolic cascade that ultimately promotes the regression process of herniated discs.
Our results demonstrated that the transcript levels of VEGF were significantly higher in protruding compared with extruded IDs. VEGF is a potent mitogen for vascular endothelial cells. It induces new vessel formation through cell proliferation and luminal formation and has also been found to increase blood vessel permeability (23). Using immunohistochemistry, Tolonen et al and Koike et al confirmed the expression of VEGF in extruded and sequestrated human IDs (9,24). A significant correlation was also observed between the expression of VEGF and the presence of vascular endothelial cells (24). Protruding ID herniation samples were not included in these studies. Ozaki et al reported the presence of small blood vessels in 73.4% of herniated ID tissues examined histopathologically, regardless of hernia type (25). They also suggested that neovascularization was present in the outermost annulus layers, even in protruding herniations. To the best of our knowledge, the present study is the first to confirm the expression of VEGF in protruding hernias at the mRNA level, thereby supporting the hypothesis that VEGF upregulation plays a key role even in the earliest stages of the neovascularization process that follows ID herniation.
PDGF is considered a key mediator in the tissue repair process. It stimulates mitogenicity and chemotaxis of neutrophils and macrophages in a variety of mesenchymal-derived cells, such as osteoblasts, chondrocytes and fibroblasts (7). It has been shown that PDGF exerts its effect on the repair process, not by altering its normal sequence, but by increasing its rate (26). It has been hypothesized that PDGF induces the production of extracellular matrix components and cellular proliferation in the ID in a similar mechanism as in the articular cartilage (9). It has also been shown to participate in the resorption process of herniated ID fragments by upregulating MMP expression (22). The findings of the present study indicated that the PDGF transcript levels were significantly higher in patients suffering from intense and acute pain. These findings support the role of PDGF as a key modulator of the repair process of damaged connective tissue, and indicate that it is primarily expressed in the acute phase of sciatica, during which it participates in a metabolic attempt to repair the herniated ID tissue.
Notably, PDGF isoforms -C and -D were the ones principally upregulated in patients with intense and acute pain. Among the PDGF isoforms, PDGF-BB is considered ‘universal’, mainly due to its ability to bind to both known forms of PDGF-receptor isotypes. In the musculoskeletal system, recombinant human PDGF-BB injection therapy has already been employed in experimental and clinical studies to promote fracture healing and foot and ankle fusion rates (27,28). It is considered important in upregulating cytokines that enhance osseous and soft-tissue healing and regeneration. However, novel information on the biological actions of PDGF-C and -D isoforms has recently been added. PDGF-D was found to strongly attract macrophages to its expression sites and to inhibit the vascular leakage that results from growth factor-induced angiogenesis. When co-expressed with VEGF-E, it promotes the in vitro stability of the newly generated vessels induced by VEGF-E alone, leading to enlarged and less leaky vessels (29). Both neovascularization and subsequent macrophage infiltration are critical stages of the acute repair response following ID herniation. Our findings therefore suggest a more favorable effect of PDGF-CC and -DD isoforms in the acute phase of disc herniation and potentially provide an alternative target for growth factor therapy research. PDGF-CC isoform is known to bind to PDGF-R A/A and A/B dimers, whereas PDGF-DD binds to B/B dimers. B-receptor expression was found to be low in connective tissue cells and to increase during inflammation (30). In the present study, the expression of each receptor isotype was detected in approximately 50% of the herniated specimens. However, it was not significant compared with the controls.
BMI exhibited significant negative correlations with the transcript levels of VEGF and both PDGF-receptor isotypes. To the best of our knowledge, there is a paucity of data available on patient BMI in studies evaluating the molecular mechanisms implicated in ID herniation. As a consequence, BMI has not been previously associated with the expression of GFs in ID herniation. In a recent histo-morphological study, Weiler et al correlated histological disc alterations in surgical lumbar ID specimens with patient BMI, suggesting an accelerated course of ID degeneration in obese individuals (31). Occupational and lifestyle factors, including obesity, have been implicated in the pathogenesis of ID herniation, mainly on a biomechanical basis focusing on increased loading and structural failure of the disc (32,33). However, a metabolic pathway may also exist. The findings of this study indicated that, in the obese population, a decrease in the expression of GFs directly associated with the acute repair process and regression of herniated IDs through peripheral neovascularization and chemotaxis. A previous study demonstrated the overexpression of TNF-α in the adipose tissue of obese individuals, as well as a consequent decrease with weight loss (34). In the ID, TNF-α has been shown to play a central signaling role in the metabolic cascade that regulates the natural history of degeneration and herniation by promoting the expression of other cytokines, MMPs and various GFs (16). Although a clear association cannot be established solely on the available evidence, these findings imply a molecular effet, in addition to the biomechanical effect of obesity in the pathogenesis and natural history of ID herniation.
In conclusion, the present study provides evidence of the simultaneous molecular profile of numerous GFs in herniated lumbar ID surgical specimens supporting the hypothesis of growth factor involvement in the natural history of ID herniation. GFs synergistically act in herniated IDs. Increased VEGF expression possibly induces the neovascularization process in the earliest stages of ID herniation. PDGF isoforms -C and -D potentially play a role in the acute phase of radiculopathy in a metabolic response for tissue healing. A molecular effect, in addition to the biomechanical effect of obesity in the pathogenesis and natural history of ID herniation is also suggested.