The role of single nucleotide polymorphisms of IL-1A -889C&gt;T (rs1800587), TNF-A -238G&gt;A (rs361525), and VDR TaqI (rs731236) on susceptibility to herniated nucleus pulposus: a systematic review and meta-analysis

Azharuddin Azharuddin; Muhammad Ilmawan; Jonny Karunia Fajar; Marhami Fahriani; Sukamto S. Mamada; Helnida Anggun Maliga; Firzan Nainu; Kuldeep Dhama; Harapan Harapan; Rahadyan Magetsari

doi:10.12688/f1000research.53235.3

Home Browse The role of single nucleotide polymorphisms of IL-1A -889C>T (rs1800587),...

ALL Metrics

Views

Downloads

Get PDF

Get XML

Export

▬

✚

Systematic Review

Revised

The role of single nucleotide polymorphisms of IL-1A -889C>T (rs1800587), TNF-A -238G>A (rs361525), and VDR TaqI (rs731236) on susceptibility to herniated nucleus pulposus: a systematic review and meta-analysis

[version 3; peer review: 2 approved]

Previous title: The role of single nucleotide polymorphisms of IL-1A -889C>T (rs1800587), TNF-A -238G>A (rs361525), and VDR TaqI (rs731236) on susceptibility to herniated nucleus pulposus

Azharuddin Azharuddin ^1,2, Muhammad Ilmawan³, Jonny Karunia Fajar^4,5, [...] Marhami Fahriani⁴, Sukamto S. Mamada⁶, Helnida Anggun Maliga³, Firzan Nainu⁶, Kuldeep Dhama⁷, Harapan Harapan⁴, Rahadyan Magetsari^8,9

Azharuddin Azharuddin ^1,2, Muhammad Ilmawan³, [...] Jonny Karunia Fajar^4,5, Marhami Fahriani⁴, Sukamto S. Mamada⁶, Helnida Anggun Maliga³, Firzan Nainu⁶, Kuldeep Dhama⁷, Harapan Harapan⁴, Rahadyan Magetsari^8,9

PUBLISHED 23 Aug 2021

Author details Author details

¹ Department of Orthopedic and Traumatology, School of Medicine, Universitas Syiah Kuala, Banda Aceh, Aceh, 23111, Indonesia
² Department of Orthopedic and Traumatology, Dr. Zainoel Abidin Hospital, Banda Aceh, Aceh, 24415, Indonesia
³ Faculty of Medicine, Universitas Brawijaya, Malang, East Java, 65117, Indonesia
⁴ Medical Research Unit, School of Medicine, Universitas Syiah Kuala, Banda Aceh, Aceh, 23111, Indonesia
⁵ Brawijaya Internal Medicine Research Center, Department of Internal Medicine, Faculty of Medicine, Universitas Brawijaya, Malang, East Java, 65145, Indonesia
⁶ Faculty of Pharmacy, Hasanuddin University, Tamalanrea, Makassar, South Sulawesi, 90245, Indonesia
⁷ Division of Pathology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
⁸ Orthopaedic and Traumatology Division, Department of Surgery, Faculty of Medicine, Universitas Gadjah Mada, Yogyakarta, DKI Yogyakarta, 55281, Indonesia
⁹ Orthopaedic and Traumatology Division, Department of Surgery, Dr. Sardjito Hospital, Yogyakarta, DKI Yogyakarta, 55281, Indonesia

Azharuddin Azharuddin
Roles: Conceptualization, Data Curation, Funding Acquisition, Investigation, Methodology, Project Administration, Resources, Validation, Writing – Original Draft Preparation, Writing – Review & Editing

Muhammad Ilmawan
Roles: Data Curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – Original Draft Preparation

Jonny Karunia Fajar
Roles: Data Curation, Formal Analysis, Investigation, Methodology, Resources, Software, Supervision, Validation

Marhami Fahriani
Roles: Project Administration, Validation, Writing – Original Draft Preparation, Writing – Review & Editing

Sukamto S. Mamada
Roles: Investigation, Validation, Writing – Original Draft Preparation, Writing – Review & Editing

Helnida Anggun Maliga
Roles: Data Curation, Investigation, Methodology, Software, Validation, Writing – Original Draft Preparation

Firzan Nainu
Roles: Investigation, Methodology, Validation, Writing – Original Draft Preparation

Kuldeep Dhama
Roles: Project Administration, Resources, Supervision, Validation, Writing – Review & Editing

Harapan Harapan
Roles: Conceptualization, Investigation, Methodology, Resources, Supervision, Validation, Writing – Original Draft Preparation, Writing – Review & Editing

Rahadyan Magetsari
Roles: Conceptualization, Supervision, Validation, Writing – Review & Editing

OPEN PEER REVIEW

REVIEWER STATUS

This article is included in the Cell & Molecular Biology gateway.

Abstract

Background: The pathogenesis of herniated nucleus pulposus (HNP) is complex and may involve the wide variety of gene polymorphism. However, the reports from the existing studies are inconclusive. The objective of this study was to determine the role of single nucleotide polymorphisms (SNPs) in interleukin 1 alpha (IL-1A), tumor necrosis factor-alpha (TNF-A), and vitamin D receptor (VDR) genes on the susceptibility to herniated nucleus pulposus (HNP).
Methods: Four databases (PubMed, Embase, Cochrane, and Web of Science) were searched as of April 1 ^st, 2021. Authors, publication year, targeted genes, genotype and allele frequency in each case and control groups were collected. Newcastle-Ottawa scale was used to evaluate the publication quality. The pooled estimates of association of IL-1A -889C>T (rs1800587), TNF-A -238G>A (rs361525), and VDR TaqI (rs731236) and susceptibility to HNP were assessed using Z test.
Results: We screened 3,067 unique studies for eligibility and three, two and nine case-control studies on IL-1A -889C>T, TNF-A -238G>A, and VDR TaqI were included, respectively, in our meta-analysis. The studies consisting 369 HNP cases and 433 controls for IL-1A -889C>T, 252 cases and 259 controls for TNF-A -238G>A and 1130 cases and 2096 controls for VDR TaqI. Our pooled estimates indicated that there was no significant association of those SNPs with the susceptibility to HNP in any genotype, dominant model, recessive model, or allele comparations.
Conclusion: Although individual studies suggested the important role of gene expression dysregulation associated with SNPs in IL-1A, TNF-A, and VDR, our data indicated that IL-1A -889C>T, TNF-A -238G>A, and VDR TaqI had weak association with HNP susceptibility in both genotypes and allele distributions. However, since heterogeneity was identified among studies included in this meta-analysis, further meta-analysis with a larger population and subgroup analysis on specific population are warranted to support this finding.

Keywords

Spinal disc herniation, SNP, IL-1A, TNF-A, VDR

Corresponding author: Azharuddin Azharuddin

Competing interests: No competing interests were disclosed.

Grant information: The author(s) declared that no grants were involved in supporting this work.

Copyright: © 2021 Azharuddin A et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite: Azharuddin A, Ilmawan M, Fajar JK et al. The role of single nucleotide polymorphisms of IL-1A -889C>T (rs1800587), TNF-A -238G>A (rs361525), and VDR TaqI (rs731236) on susceptibility to herniated nucleus pulposus: a systematic review and meta-analysis [version 3; peer review: 2 approved]. F1000Research 2021, 10:419 (https://doi.org/10.12688/f1000research.53235.3) First published: 25 May 2021, 10:419 (https://doi.org/10.12688/f1000research.53235.1) Latest published: 23 Aug 2021, 10:419 (https://doi.org/10.12688/f1000research.53235.3)

Revised Amendments from Version 2

We have added additional discussion about compression of our finding with previous studies based on recommendation of the reviewer. We added: The role of gene polymorphism of IL-1A (rs1800587) was rarely reported. Studies found that IL-1A (rs1800587) has no association with rheumatoid arthritis and systemic sclerosis.^{52, 53} Another study also found no association between TNF-A gene variants and the risk of bone-joint and spinal tuberculosis.⁵⁴ Furthermore, a study assessing the role of VDR gene polymorphism in the case of osteoarthritis also failed to show the potential contribution of VDR gene variant in the development of osteoarthritis.⁵⁵ These evidences suggested that, while theoretically the gene variants of IL-A, TNF-A, and VDR may contribute to the development of bone diseases, the evidence reveals otherwise.

See the authors' detailed response to the review by Ismail Hadisoebroto Dilogo
See the authors' detailed response to the review by Mahir Gachabayov

Introduction

Herniated nucleus pulposus (HNP) or disc herniation is the most common spinal degenerative disease associated with lower back pain and radicular pain of the lower extremities due to nerve compression.¹ HNP is also the most common cause of persistent sciatic pain due to displacement of the nucleus pulposus beyond the intravertebral disc space.² The most prevalent HNP locations are between the L4 and L5 vertebrae and between the L5 and S1 vertebrae³ whilst the highest incidence is observed amongst people aged 30-50 years old.⁴ Diabetes,⁵ smoking,^5,6 obesity,^5-7 type of occupation,⁸ age,⁷ and gender^4,9 have all been associated with a high risk of developing disc degenerative diseases. However, it has been suggested that genetic factors also play a vital role in susceptibility to disc degenerative diseases. A study showed that individuals aged younger than 30 years who have a family history of disc herniation have a 14.5 times higher risk of developing disc protrusion than individuals who have no family history.⁸ Family history is also attributed to a 5.1 times higher risk of disc herniation in people aged between 30-50 years old.⁸ The Twin Spine Study found that heredity substantially influences disc degeneration by 43-77%.^10,11

The intervertebral disc (IVD) consists of two different components: the nucleus pulposus (NP) and the annulus fibrosus (AF),¹² where proteoglycans (mostly found in NP) acts as an internal semi-fluid mass and collagen (mostly found in AF) acts as a laminar fibrous container.¹³ Genes encoding components of IVD such as collagens I,¹⁴ collagens IX,¹⁵ collagens XI,¹⁶ aggrecan,¹⁷ cartilage intermediate layer protein (CILP),¹⁸ and vitamin D receptor (VDR)¹⁹ have previously been studied to determine susceptibility to lumbar disc diseases. Other factors such as increased production of extracellular matrix-degrading enzymes (encoded by matrix metalloproteinase 3 gene (MMP-3) and MMP-9²⁰ and increased expression of inflammatory cytokines such as interleukin-1 alpha (IL-1A), IL-18,²¹ IL-6, and tumor necrosis factor-alpha (TNF-α)²²) are commonly found in disc degeneration. Excessive synthesis, secretion, and biological activity of these inflammatory mediators are associated with tissue destruction and are therefore commonly found in inflammatory disorders including disc degeneration.²³

One of the mechanisms that alters the production of protein mediators in the human body are single-nucleotide polymorphisms (SNPs). These genetic variations, single nucleotide changes at specific positions in a gene, may influence gene expression and hence associate to particular disease. A three-fold increase in susceptibility of disc degeneration was observed in individuals with a TT genotype compared to those without the allele (CC genotype) on SNP IL-1A -889C>T (rs1800587).²⁴ People with minor allele of IL-1A -889C>T (T allele) also had a 2.4-fold increased risk of disc bulges²⁴ and a 2.5-fold increased risk of endplate modic change.²⁵ A study in an Iranian population found that among nine SNPs on pro-inflammatory cytokine genes (IL-1, IL-6 and TNF-A), no association to IVD degeneration was found except for two SNPs in the TNF-A gene (TNF-A−308 G/A and TNF-A −238 G/A).¹⁹ TNF-α plays important role in the pathophysiology of HNP such as upregulating the activity and the gene expression of MMP, stimulating other cytokines such as IL-1, IL-6, and IL-8, stimulating cell migration, altering endothelial permeability, and decreasing the synthesis of collagen and proteoglycan.²⁶ A study reported that G allele and GG genotype of TNF-A 238G>A (rs361525) were 2.51 times and 2.98 times, respectively, more prevalent in patients with HNP compared to healthy controls.¹⁹

Several roles of VDR such as regulating chondrocyte proliferation and differentiation, bone mineralization and remodeling, and matrix production have previously been demonstrated.²⁷ VDR's role in spinal degenerative disorder has been studied in Italian,²⁸ Turkish,²⁹ and Southern European populations.³⁰ A study in a Chinese population suggested that subjects with the t allele of VDR TaqI (rs731236) had a 2.61 times higher risk to have degenerative disc disease.³¹ Moreover, individuals aged younger than 40 years who had the t allele were almost six times more likely to develop disc degeneration and 7.17 times more likely to develop disc bulge compared to those without the t allele.³¹ However, studies in Danish³² and Mexican populations³³ contradict previous results suggesting a role of VDR TaqI in disc degenerative disease. This conflicting role of SNPs inIL-1A, TNF-A and VDR on HNP therefore needs to be further evaluated. This study sought to determine the association of IL-1A -889C>T (rs1800587), TNF-A 238G>A (rs361525), and VDR TaqI (rs731236) in susceptibility to HNP.

Methods

Study design and protocol

A systematic review and meta-analysis were conducted to assess the association of three SNPs, IL-1A (rs1800587), TNF-A (rs361525), and VDR (rs731236), on susceptibility to HNP. The outcome variable of this study was the risk or susceptivity to have HNP while the response variables were the SNPs in three genes: IL-1A (rs1800587), TNF-A (rs361525), and VDR (rs731236). We searched databases for relevant studies, then extracted and analyzed data from those studies to achieve the pooled odds ratios (ORs) and 95% confidence interval (95%CI) using a random or fixed effect model depending on the data. This study was conducted in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guideline.^34,35 The protocol of this study has been registered in PROSPERO (reg. number CRD42021249187).

Literature search strategy

The literature searches were conducted on PubMed, Embase, Cochrane, and Web of Science. The searches were conducted using the keywords:‘degenerative disc disease’ AND (‘IL-1A’ OR ‘rs1800587’ OR ‘-889C>T’) OR (‘TNF-A’ OR ‘rs361525’ OR ‘-238G>A’) OR (‘VDR’ OR ‘rs731236’ OR ‘TaqI’) AND ‘gene polymorphism’, including all results up to April 1^st, 2021. The keywords were adapted from Medical Subject Heading (MeSH). Additional studies were also retrieved from the references of relevant papers. If two or more studies with the same study data were identified, the most recent study was used. The processes were conducted by three independent authors (JKF, MI, HAM).

Study eligibility

To be eligible for the meta-analysis, a study had to meet all the inclusion criteria below: (1) the study design should be case-control, cross-sectional, or cohort design; (2) the study should evaluate the association of IL-1 (rs1800587), TNF-A (rs361525), or VDR (rs731236) on HNP and have case and control groups; and (3) studies should present genotype frequency or minor allele frequency (MAF). All studies with duplicate records, poor quality or which had deviation from Hardy-Weinberg Equilibrium (HWE) were excluded.³⁶

Data extraction

Important information from the studies such as first author name, year of publication, names of targeted gene and the SNP, genotype frequency, or MAF from case and control groups were collected. The allele frequency and MAF were recalculated using Mendel’s law. Data extraction processes were conducted by three independent authors (JKF, MI, HAM) and consensus established together with senior authors (AA, HH) if discrepancies were found.

Quality assessment

The quality of the included studies was evaluated using Newcastle-Ottawa Score (NOS)³⁷ by three independent authors (JKF, MI, HAM). This evaluation was conducted to ensure the quality of three fundamental methodological parameters of the studies: patient selection (four points), comparability of the groups (two points), and ascertainment of exposure (three points); NOS ranged from 0 to 9. Each study was then categorized based on the NOS: (1) good quality (NOS ≥ 7); (2) moderate quality (NOS ≥ 5); or (3) poor quality (NOS < 5). Consensus was established if discrepancies were found.

Covariates and sub-group analysis

The outcome measure in our study was the incidence of HNP while the predictor covariates were the gene polymorphisms of IL-1A (rs1800587), TNF-A (rs361525), and VDR (rs731236). All genetic models were applied to describe the role of each gene variant in the pathogenesis of HNP. For IL-1A (rs1800587), the allele models were C vs. T and T vs. C; and the genotype models were CC vs. CT + TT, CT vs. CC + TT, and TT vs. CC + CT. For TNF-A (rs361525), the allele models were G vs. A and A vs. G; and the genotype models were GG vs. GA + AA, GA vs. GG + AA, and AA vs. GG + GA. For VDR (rs731236), the allele models were T vs. C and C vs. T; and the genotype models were TT vs. TC + CC, TC vs. TT + CC, and CC vs. TT + TC.

Statistical analysis

To assess the association of IL-1A (rs1800587), TNF-A (rs361525), and VDR (rs731236) on HNP, a Z-test was employed. The Egger test was used to evaluate the publication bias and a p < 0.05 indicated the possibility of publication bias in each calculated result. The Q test was used to evaluate the heterogeneity and decide between random and fixed-effect models for OR calculation. If heterogeneity was indicated (p-value less than 0.10), the random effect model was used; otherwise, the fixed-effect model was used. All analyses were performed using ‘meta’,³⁸ ’ metafor’,³⁹ and ‘dmetar’⁴⁰ packages in R version 4.0.4.⁴¹

Results

Study eligibility results

The literature searches yielded 3,199 articles of which 3,067 references were retained after removing duplicates. Screening of the titles and abstracts excluded 2,965 articles as they did not meet the inclusion criteria. After a further screening of full text, an additional 90 studies were excluded due to lack of relevance (n = 90), incomplete data (n = 5), and HWE deviation (n = 6) (Figure 1). 12 studies were included in the meta-analysis: three studies for IL-1A (rs1800587),^19,42,43 two studies for TNF-A (rs361525),^19,44 and nine studies for VDR (rs731236).^28,42,45-51 The summary of studies included in the meta-analysis is presented in Table 1.

Figure 1. Flowchart of the result of literature searches according to the preferred reporting items of systematic reviews and meta-analyses (PRISMA).

Table 1. Study characteristics of IL-1A -889C>T, TNF-A 238G>A, and VDR TaqI (rs731236) on HNP.

Gene	SNP code	Author	Year	Country	Study design	NOS	Study groups
							Case								Control
							Genotype				Allele			HWE	Genotype				Allele			HWE
							AA	AB	BB	n	A	B	n	HWE	AA	AB	BB	n	A	B	n
IL1A	-889C>T						CC	CT	TT		C	T			CC	CT	TT		C	T
		Serrano et al⁴²	2014	Mexico	Case-control	8	51	45	4	100	147	53	200	2.41	55	35	10	100	145	55	200	1.49
		Chen et al⁴³	2018	China	Case-control	8	87	78	28	193	252	134	386	2.27	102	81	14	197	285	109	394	0.15
		Abdollahzade et al¹⁹	2018	Iran	Case-control	8	33	33	10	76	99	53	152	0.15	62	62	12	136	186	86	272	0.40
TNF-A	238G>A						GG	GA	AA		G	A		GG	GA	AA		G	A
		Aparicio et al⁴⁴	2011	Spain	Case-control	7	45	5	0	50	95	5	100	0.14	113	9	0	122	235	9	244	0.18
		Abdollahzade et al¹⁹	2018	Iran	Case-control	8	61	15	0	76	137	15	152	0.91	79	57	1	137	215	59	274	7.32
VDR	TaqI						TT	TC	CC		T	C			TT	TC	CC		T	C
		Chen et al⁴⁵	2012	China	Case-control	8	79	2	0	81	160	2	162	0.01	86	14	1	101	186	16	202	0.25
		Cheung et al⁴⁶	2012	China	Case-control	9	92	15	1	108	199	17	216	0.19	103	3	0	106	209	3	212	0.02
		Colombini et al²⁸	2016	Italy	Case-control	8	114	117	35	266	345	187	532	0.33	106	109	37	252	321	183	504	1.06
		Eser et al⁴⁷	2010	Turkey	Case-control	8	65	67	17	150	198	102	300	0.00	67	66	16	150	201	99	300	0.00
		Li et al⁴⁸	2018	China	Case-control	8	114	6	0	114	234	6	240	0.08	109	11	0	120	229	11	240	0.28
		Oishi et al⁴⁹	2003	Japan	Case-control	8	31	8	0	39	70	8	78	0.51	16	5	0	21	37	5	42	0.38
		Omair et al⁵⁰	2012	Norway	Case-control	9	53	70	23	147	176	116	292	0.00	61	92	35	188	214	163	376	0.00
		Serrano et al⁴²	2014	Mexico	Case-control	8	69	27	4	69	165	35	200	0.42	62	35	3	100	159	41	200	0.54
		Yuan et al⁵¹	2010	China	Case-control	8	156	22	0	156	334	22	356	0.77	256	28	0	284	540	28	568	0.76

Distribution of allele and genotype frequency of IL-1A -889C>T, TNF-A 238G>A and VDR TaqI

Our data indicated that the TT genotype of IL-1A -889C>T was 1.37 more frequent in HNP patients than in controls, while the distribution of alleles and other genotypes were similar between patients and healthy controls. CT genotype and T allele of TNF-A 238G>A were both 1.6 times more frequent in healthy controls than in HNP cases. No difference in the distribution of alleles or genotypes between HNP and controls was observed in VDR TaqI (Table 1).

Association between alleles and genotypes of IL-1A, TNF-A, and VDR polymorphism and HNP

Our pooled estimates suggested that no IL-1A -889C>T genotypes were associated with the risk of HNPs with CC vs. CT+TT (OR: 0.82, 95%CI: 0.62, 1.09), CT vs. CC+TT (OR: 0.07; 95%CI: 0.81, 1.42), and TT vs. CT+CC (OR: 1.20; 95%CI: 0.13, 11.37) (Table 2 and Figure 2). The pooled data also suggested that allele frequency of IL-1A -889C>T had no significant association with the susceptibility to HNP with OR: 0.83; 95%CI: 0.67, 1.02 for C allele compared to T allele.

Figure 2. The forest plot of the association of IL-1A-889C>T and HNP.

(A) CC vs CT+TT (OR: 0.82; 95%CI: 0.62, 1.09); p-value0.170, p-Het 0.867, and p-Egger 0.193, (B) CT vs CC+TT (OR1.07; CI95%: 0.81, 1.42); p-value 0.637; p-Het 0.373; p-Egger 0.675, (C) TT vs CT+CC (OR: 1.20; CI95%: 0.13, 11.37); p-value 0.764, p-Het 0.039, p-Egger 0.177, (D) C vs T (OR0.83; 95%CI: 0.67, 1.02); p-value 0.081; p-Het 0.369; p-Egger 0.221), and (E) T vs C (OR: 1.21; CI95%: 0.98, 1.50); p-value 0.081; p-Het 0.369; p-Egger 0.221.

Table 2. Associations of genotypes and alleles ofIL-1A -889C>T, TNF-A 238G>A, and VDR TaqI (rs731236) on HNP.

Gene	Allele/genotype model	Number of studies	Model	OR (CI 95%)	p-value	p-Het	p-Egger
IL-1A	CC vs. CT+TT	3	Fixed	0.82 (0.62, 1.09)	0.170	0.867	0.193
	CT vs. CC+TT	3	Fixed	1.07 (0.81, 1.42)	0.637	0.373	0.675
	TT vs. CT+CC	3	Random	1.20 (0.13, 11.37)	0.764	0.039	0.177
	C vs. T	3	Fixed	0.83 (0.67, 1.02)	0.081	0.369	0.221
	T vs. C	3	Fixed	1.21 (0.98, 1.50)	0.081	0.369	0.221
TNF-A	GG vs. GA+AA	2	Random	1.60 (0.00,12882.74)	0.628	0.034	0.889
	GA vs. GG+AA	2	Random	0.63 (0.00, 4211.31)	0.629	0.038	0.864
	AA vs. GG+GA	2	Random	1.34 (0.16, 6712.96)	0.955	0.587	<0.001
	G vs. A	2	Random	1.60 (0.00, 12882.74)	0.628	0.034	0.744
	A vs. G	2	Random	0.67 (0.00, 1555.67)	0.629	0.056	0.744
VDR	TT vs. TC+CC	9	Random	2.65 (0.60, 11.85)	0.172	<0.001	0.986
	TC vs. TT+CC	9	Random	1.01 (0.55, 1.85)	0.964	0.041	0.791
	CC vs. TC+TT	9	Fixed	0.94 (0.69, 1.28)	0.688	0.970	0.3871
	T vs. C	9	Random	1.06 (0.91, 1.22)	0.775	0.023	0.874
	C vs. T	9	Random	0.92 (0.50, 1.72)	0.775	0.023	0.874

Pooled estimates for allele and genotype distribution of the TNF-A 238G>A also had no significant association with the risk for HNP. No association was observed between genotype models and the risk of HNP: GG vs. GA+AA (OR: 1.60; 95%CI: 0.00, 12882.74), GA vs. GG+AA (OR: 0.63; 95%CI: 0.00, 4211.31), and AA vs. GG+GA (OR: 1.34; 95%CI: 0.16, 6712.96) (Table 2 and Figure 3). Distribution of the allele also had no strong association with HNP susceptibility.

Figure 3. The forest plot of the association of TNF-A 238G>A and HNP.

(A) GG vs GA+AA (OR1.60; CI95%: 0.00, 12882.74); p-value 0.628, p-Het 0.034, p-Egger 0.889, (B) GA vs GG+AA (OR0.63; CI95%: 0.00, 4211.31); p-value 0.629, p-Het 0.038, p-Egger 0.864, (C) AA vs GG+GA (OR1.05; CI95%: 0.16, 6712.96); p-value 0.955, p-Het 0.587, p-Egger <0.001, (D) G vs A (OR1.49; 95%CI: 0.00, 3473.60); p-value 0.629, p-Het 0.034, p-Egger 0.744, (E) A vs G (OR: 0.67; 95%CI:0.00, 1555.67); p-value 0.629, p-Het 0.056, p-Egger 0.744.

Our estimates for genotypes of VDR TaqI (rs731236) suggested that none of the genotypes were associated with susceptibility to degenerative disc disease HNP with OR: 2.65; 95%CI: 0.60, 11.85 for TT vs. TC+CC, OR: 1.01; 95%CI: 0.55, 1.85 for TC vs. TT+CC and OR: 0.94; 95%CI: 0.69, 1.28 (Table 2 and Figure 4). None of the alleles of VDR TaqI (rs731236) were associated with HNP; people with T allele had OD 1.06 with 95%CI: 0.91, 1.22 for HNP.

Figure 4. The forest plot of the association of VDR TaqI (rs731236) and HNP.

(A) TT vs TC+CC (OR: 2.65; CI95%: 0.60, 11.85); p-value 0.172, p-Het <0.001, p-Egger 0.986, (B) TC vs TT+CC (OR1.01; CI95%: 0.55, 1.85); p-value 0.964, p-Het 0.041, p-Egger 0.791, (C) CC vs TC+TT (OR 0.94; CI95%: 0.69, 1.28); p-value 0.688, p-Het 0.970, p-Egger 0.387, (D) T vs C (OR: 1.06; CI95%: 0.91, 1.22); p-value 0.775, p-Het 0.023, p-Egger 0.874), (E) C vs T (OR: 0.92; 95%CI: 0.50, 1.72); p-value 0.775, p-Het 0.023, p-Egger 0.874.

Discussion

Our present failed to clarify the role of IL-1A (rs1800587), TNF-A (rs361525), and VDR (rs731236) on the pathogenesis of HNP. To the best of our knowledge, our current study is the first study providing the holistic gene polymorphism in the case of HNP. Therefore, the comprehensive comparison in the context of methodological quality between our study and previous studies was unable to perform. The role of gene polymorphism of IL-1A (rs1800587) was rarely reported. Studies found that IL-1A (rs1800587) has no association with rheumatoid arthritis and systemic sclerosis.⁵²^,⁵³ Another study also found no association between TNF-A gene variants and the risk of bone-joint and spinal tuberculosis.⁵⁴ Furthermore, a study assessing the role of VDR gene polymorphism in the case of osteoarthritis also failed to show the potential contribution of VDR gene variant in the development of osteoarthritis.⁵⁵ These evidences suggested that, while theoretically the gene variants of IL-A, TNF-A, and VDR may contribute to the development of bone diseases, the evidence reveals otherwise. For the negative findings of our study, several possible reasons might be proposed. First, HNP is a complex disease and is cause by multiple factors. Thus, no single factor such as a single SNP is responsible for the whole pathogenesis. Second, large variations in the number of samples or allele frequencies among studies in our meta-analysis also contribute to the findings. This probably relates to differences in populations where study data were collected. Finally, the small number of samples significantly influenced the results of our meta-analysis. Therefore, studies with larger sample sizes and sub-analyses for different populations such as Asian, Caucasian, and other populations are warranted whenever more data are available.

HNP occurs when the central part of the intervertebral disc, the nucleus pulposus, herniates through the surrounding part of the disc, the annulus fibrosus. The damage to the annulus fibrosus resulting in the herniated disc may be associated with factors such as gender, age, certain activities such as lifting of weights and carrying, and being overweight.⁵⁶ For example, the degeneration of disc organization could occur during aging as the regulation of the extracellular matrix (ECM), a major component of the disc, is damaged during the aging process.⁵⁶ The herniated disc or HNP may occur as a result of several pathological mechanisms. Those mechanisms ultimately cause imbalances in disc composition that are directly linked to the quality of the ECM.⁵⁶ Therefore, the balance between the ECM and its degrading enzymes, such as matrix metalloproteases (MMPs), seems to be the key to maintaining normal disc function. Interestingly, of several mechanical pathways, the regulation of the ECM and the MMPs could also be determined by the occurrence of single nucleotide polymorphisms in genes responsible for ECM regulation as explained below.⁵⁷

Two major structural proteins that are important in the matrix structure of the disc are collagen and proteoglycan.⁵⁶ The main proteoglycan found in the normal intervertebral disc is aggrecan.⁵⁸ Although both annulus fibrosus and nucleus pulposus are mainly composed of water, proteoglycan, and collagen, the level of those contents differ between the structures. The annulus fibrosus consists of 70% water, 15% collagen, and 5% proteoglycan, while the composition in the nucleus pulposus is 77% water, 4% collagen, and 14% proteoglycan.⁵⁶ Any event causing disturbances in those ratios and/or in synthesis and/or degradation of those structural proteins could lead to herniated disc problems. For example, hypoxic and acidic conditions could repress the synthesis of the matrix leading to the dysfunctionalities of the cells of the disc.^59,60

The role of inflammation in predisposing the disc to damage has also been revealed. Specifically, proinflammatory cytokines are found to play a specific role in herniated and degenerated discs.^58,61-63 A controlled immunohistochemical study observed the accumulation of inflammatory cells, mainly macrophages, in herniated disc cells indicating the role of proinflammatory cytokines in the disease.⁶² For example, IL-1β could induce the annulus fibrosus to generate inflammatory factors leading to the impairment of proteoglycan aggregation.⁶³ Another cytokine, TNF-α, is also involved in the development of intervertebral disc problems.⁶⁴ However, it seems that its effect is less significant than IL-1.⁶⁵ This finding may be related to its relatively lower expression compared to IL-1 in the normal and healthy disc.^56,66

Penetration of those inflammatory cells or proteins could be caused by matrix loss. In normal conditions, aggrecan should prevent the penetration of various compounds, especially serum proteins and cytokines.⁶⁷ Therefore, in addition to its pivotal role in maintaining sufficient hydration to the disc, proteoglycan loss could stimulate the movement of cytokines towards the disc activating the inflammation cascade.⁵⁶

One of the mechanisms by which the proinflammatory cytokines, such as IL-1α and TNF-α, generate problems in the intervertebral disc is related to their effect on inducing MMP production.^68-71 The exaggerated activity of MMPs causes excessive degradation of collagen and proteoglycan.⁵⁸ Another mechanism is associated with the activity of the cytokines in inhibiting tissue inhibitors of MMPs (TIMPs) which are responsible for terminating the action of MMPs.⁵⁸ Taken together, those actions ultimately impair disc functionality.

As the normal intervertebral disc is relatively avascular and aneural,⁵⁶ the nutritional supply to the disc depends on the ability of the nutrients to diffuse from the closest vascularized structure outside the disc which are the vertebral bodies.⁷² The nutrients then penetrate the cartilaginous endplate and finally reach the annulus fibrosus and nucleus pulposus.^72,73 Accordingly, calcification of the endplate would diminish diffusion of vital nutrients, leading to the death of the disc cells.⁷⁴ Therefore, VDR plays a critical role as this ligand-dependent transcription factor is involved in regulating calcium homeostasis and bone mineralization in the body, including in the intervertebral disc.^75,76 It has been known that genetic polymorphisms occurring in genes encoding VDR are associated with intervertebral disc problems,⁷⁶ including herniated disc.

Our current study had several limitations. First, the pertinent confounding factors that might affect the final findings of our study were not included in analysis, such as mechanical and behavioral factors and the levels of proteoglycans. Second, since the included articles in our current study were non-randomized controlled trials (RCTs), the final findings might have the higher risk of bias. Therefore, the up-coming meta-analysis was expected to involve only RCTs studies. Third, the limited reports on the context of gene polymorphism in HNP had made our study included the limited number of studies. Therefore, our findings should be interpreted with caution.

In conclusion, our results suggest that well-regulated IL-1A, TNF-A, and VDR are important for normal intervertebral discs and that dysregulation of these could negatively affect the intervertebral discs. Some individual studies found that SNPs in IL-1A (rs1800587), TNF-A (rs361525), and VDR (rs731236) were associated with the susceptibility to HNP, however, our meta-analysis suggested that the effects are not robust.

Data availability

Underlying data

All data underlying the results are available as part of the article and no additional source data are required.

Reporting guidelines

Figshare: PRISMA checklist for ‘The role of single nucleotide polymorphisms of IL-1A -889C>T (rs1800587), TNF-A -238G>A (rs361525), and VDR TaqI (rs731236) on susceptivity of herniated nucleus pulposus’, https://doi.org/10.6084/m9.figshare.14479233.⁷⁷

References

1. Kim DH, Lee SH, Kim KT, et al.: Association of interleukin-1 receptor antagonist gene polymorphism with response to conservative treatment of lumbar herniated nucleus pulposus. Spine. 2010 Jul 15; 35(16): 1527–31. PubMed Abstract | Publisher Full Text
2. De Cicco FL, Willhuber GOC: Nucleus pulposus herniation. StatPearls [Internet]. 2020.
3. Jordan J, Konstantinou K, O'Dowd J: Herniated lumbar disc. BMJ Clin Evid. 2009; 2009: 1118. eng. PubMed Abstract | Free Full Text
4. Fjeld OR, Grøvle L, Helgeland J, et al.: Complications, reoperations, readmissions, and length of hospital stay in 34 639 surgical cases of lumbar disc herniation. Bone Joint J. 2019; 101-B(4): 470–7. PubMed Abstract | Publisher Full Text
5. Jakoi AM, Pannu G, D'Oro A, et al.: The clinical correlations between diabetes, cigarette smoking and obesity on intervertebral degenerative disc disease of the lumbar spine. Asian Spine J. 2017; 11(3): 337–47. PubMed Abstract | Publisher Full Text | Free Full Text
6. Lener S, Wipplinger C, Hartmann S, et al.: The impact of obesity and smoking on young individuals suffering from lumbar disc herniation: a retrospective analysis of 97 cases. Neurosurg Rev. 2020; 43(5): 1297–303. PubMed Abstract | Publisher Full Text | Free Full Text
7. Teraguchi M, Yoshimura N, Hashizume H, et al.: Prevalence and distribution of intervertebral disc degeneration over the entire spine in a population-based cohort: the Wakayama Spine Study. Osteoarthritis Cartilage. 2014; 22(1): 104–10. PubMed Abstract | Publisher Full Text
8. Y-g Z, Sun Z, Zhang Z, et al.: Risk factors for lumbar intervertebral disc herniation in Chinese population: A case-control study. Spine. 2009; 34(25). PubMed Abstract | Publisher Full Text
9. Wang YXJ, Griffith JF, Ma HT, et al.: Relationship between gender, bone mineral density, and disc degeneration in the lumbar spine: a study in elderly subjects using an eight-level MRI-based disc degeneration grading system. Osteoporosis Int. 2011; 22(1): 91–6. PubMed Abstract | Publisher Full Text
10. Battié MC, Videman T, Kaprio J, et al.: The Twin Spine Study: Contributions to a changing view of disc degeneration. Spine J. 2009; 9(1): 47–59. PubMed Abstract | Publisher Full Text
11. Battié MC, Videman T, Gibbons LE, et al.: Determinants of lumbar disc degeneration: A study relating lifetime exposures and magnetic resonance imaging findings in identical twins. Spine. 1995; 20(24): 2601–12. Publisher Full Text
12. Niosi CA, Oxland TR: Degenerative mechanics of the lumbar spine. Spine J. 2004; 4(6): S202–S8. PubMed Abstract | Publisher Full Text
13. Yoon ST, Patel NM: Molecular therapy of the intervertebral disc. Eur Spine J. 2006; 15 Suppl 3(Suppl 3): S379–S88. PubMed Abstract | Publisher Full Text | Free Full Text
14. Pluijm SMF, van Essen HW, Bravenboer N, et al.: Collagen type I α1 Sp1 polymorphism, osteoporosis, and intervertebral disc degeneration in older men and women. Ann Rheum Dis. 2004; 63(1): 71–7. PubMed Abstract | Publisher Full Text | Free Full Text
15. Annunen S, Paassilta P, Lohiniva J, et al.: An allele of COL9A2 associated with intervertebral disc disease. Science. 1999; 285(5426): 409. PubMed Abstract | Publisher Full Text
16. Mio F, Chiba K, Hirose Y, et al.: A functional polymorphism in COL11A1 which encodes the alpha 1 chain of type XI collagen, is associated with susceptibility to lumbar disc herniation. Am J Hum Genet. 2007; 81(6): 1271–7. PubMed Abstract | Publisher Full Text | Free Full Text
17. Kawaguchi Y, Osada R, Kanamori M, et al.: Association between an aggrecan gene polymorphism and lumbar disc degeneration. Spine. 1999; 24(23): 2456–60. PubMed Abstract | Publisher Full Text
18. Seki S, Kawaguchi Y, Chiba K, et al.: A functional SNP in CILP, encoding cartilage intermediate layer protein, is associated with susceptibility to lumbar disc disease. Nat Genet. 2005; 37(6): 607–12. PubMed Abstract | Publisher Full Text
19. Abdollahzade S, Hanaei S, Sadr M, et al.: Significant association of TNF-α, but not other pro-inflammatory cytokines, single nucleotide polymorphisms with intervertebral disc degeneration in Iranian population. Clin Neurol Neurosurg. 2018; 173: 77–83. PubMed Abstract | Publisher Full Text
20. Takahashi M, Haro H, Wakabayashi Y, et al.: The association of degeneration of the intervertebral disc with 5a/6a polymorphism in the promoter of the human matrix metalloproteinase-3 gene. J Bone Joint Surg Br. 2001; 83-B(4): 491–5. PubMed Abstract | Publisher Full Text
21. Videman T, Saarela J, Kaprio J, et al.: Associations of 25 structural, degradative, and inflammatory candidate genes with lumbar disc desiccation, bulging, and height narrowing. Arthritis Rheum. 2009; 60(2): 470–81. PubMed Abstract | Publisher Full Text
22. Yang Q, Liu Y, Guan Y, et al.: Vitamin D Receptor gene polymorphisms and plasma levels are associated with lumbar disc degeneration. Sci Rep. 2019; 9(1): 7829. PubMed Abstract | Publisher Full Text | Free Full Text
23. Kalb S, Martirosyan NL, Kalani MYS, et al.: Genetics of the degenerated intervertebral disc. World Neurosurg. 2012; 77(3): 491–501. PubMed Abstract | Publisher Full Text
24. Solovieva S, Kouhia S, Leino-Arjas P, et al.: Interleukin 1 polymorphisms and intervertebral disc degeneration. Epidemiology. 2004; 15(5): 626–33. PubMed Abstract | Publisher Full Text
25. Karppinen J, Solovieva S, Luoma K, et al.: Modic changes and interleukin 1 gene locus polymorphisms in occupational cohort of middle-aged men. Eur Spine J. 2009; 18(12): 1963–70. PubMed Abstract | Publisher Full Text | Free Full Text
26. Weiler C, Nerlich AG, Bachmeier BE, et al.: Expression and distribution of tumor necrosis factor alpha in human lumbar intervertebral discs: a study in surgical specimen and autopsy controls. Spine. 2005; 30(1): 44–53. PubMed Abstract | Publisher Full Text
27. Zwerina K, Baum W, Axmann R, et al.: Vitamin D receptor regulates TNF-mediated arthritis. Ann Rheum Dis. 2011; 70(6): 1122–9. PubMed Abstract | Publisher Full Text
28. Colombini A, Brayda-Bruno M, Lombardi G, et al.: BsmI, ApaI and TaqI polymorphisms in the vitamin D receptor gene (VDR) and association with lumbar spine pathologies: an Italian case-control study. PloS one. 2016; 11(5): e0155004. PubMed Abstract | Publisher Full Text | Free Full Text
29. Hekimoglu B, Eraltan H, Sari H, et al.: Investigation of the vitamin D receptor (VDR) gene polymorphisms in lumbar disc herniation in Turkish patients. Experimed. 2020; 9: 93–8. Publisher Full Text
30. Toktaş ZO, Ekşi MŞ, Yılmaz B, et al.: Association of collagen I, IX and vitamin D receptor gene polymorphisms with radiological severity of intervertebral disc degeneration in Southern European Ancestor. Eur Spine J. 2015; 24(11): 2432–41. PubMed Abstract | Publisher Full Text
31. Cheung KMC, Chan D, Karppinen J, et al.: Association of the Taq I allele in vitamin D receptor with degenerative disc disease and disc bulge in a Chinese population. Spine. 2006; 31(10): 1143–8. PubMed Abstract | Publisher Full Text
32. Eskola PJ, Kjaer P, Daavittila IM, et al.: Genetic risk factors of disc degeneration among 12-14-year-old Danish children: a population study. Int J Mol Epidemiol Genet. 2010; 1(2): 158–65. PubMed Abstract | Free Full Text
33. Cervin Serrano S, González Villareal D, Aguilar-Medina M, et al.: Genetic polymorphisms of interleukin-1 alpha and the vitamin D receptor in Mexican Mestizo patients with intervertebral disc degeneration. Int J Genomics. 2014; 2014: 302568. PubMed Abstract | Publisher Full Text | Free Full Text
34. Moher D, Liberati A, Tetzlaff J, et al.: Preferred reporting items for systematic reviews and meta-analyses: The PRISMA statement. Plos Med. 2009 Jul; 6(7): e1000097. PubMed Abstract | Publisher Full Text | Free Full Text
35. Page MJ, Moher D, Bossuyt PM, et al.: PRISMA 2020 explanation and elaboration: updated guidance and exemplars for reporting systematic reviews. BMJ. 2021: 372. PubMed Abstract | Publisher Full Text | Free Full Text
36. Rodriguez S, Gaunt TR, Day IN: Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies. Am J Epidemiol. 2009; 169(4): 505–14. PubMed Abstract | Publisher Full Text | Free Full Text
37. Ilmawan M, Satriawan AG, Purnamasari Y, et al.: The association between glutathione S-transferase P1 105IIe>Val gene polymorphism and the risk of bladder cancer: A meta-analysis. Türkiye Klinikleri Tip Bilimleri Dergisi. 2020; 40(2): 250–9. Publisher Full Text
38. Balduzzi S, Rücker G, Schwarzer G: How to perform a meta-analysis with R: a practical tutorial. Evid Based Ment Health. 2019; 22(4): 153–60. PubMed Abstract | Publisher Full Text
39. Viechtbauer W: Conducting meta-analyses in R with the metafor package. J Stat Softw. 2010; 36(3): 1–48. Publisher Full Text
40. Harrer M, Cuijpers P, Furukawa T, et al.: dmetar: Companion R Package For The Guide 'Doing Meta-Analysis in R'. R package version 00. 2019; 9000.
41. R Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria; 2020. Reference Source
42. Cervin Serrano S, González Villareal D, Aguilar-Medina M, et al.: Genetic polymorphisms of interleukin-1 alpha and the vitamin d receptor in mexican mestizo patients with intervertebral disc degeneration. Int J Genomics. 2014; 2014: 302568. PubMed Abstract | Publisher Full Text | Free Full Text
43. Chen Y, Ma H, Bi D, et al.: Association of interleukin 1 gene polymorphism with intervertebral disc degeneration risk in the Chinese Han population. Biosci Rep. 2018; 38(4). PubMed Abstract | Publisher Full Text | Free Full Text
44. Aparicio JP, Bances IF, Fernández EL-A, et al.: The IL-1 β (+ 3953 T/C) gene polymorphism associates to symptomatic lumbar disc herniation. Eur Spine J. 2011; 20(3): 383–9. PubMed Abstract | Publisher Full Text | Free Full Text
45. Chen W, Li G, Sun H, et al.: Association of vitamin D receptor gene polymorphism in Han people with lumbar degenerative disc disease. Afr J Pharm Pharmacol. 2012; 6(16): 1211–5. Publisher Full Text
46. Cheung KM, Chan D, Karppinen J, et al.: Association of the Taq I allele in vitamin D receptor with degenerative disc disease and disc bulge in a Chinese population. Spine. 2006; 31(10): 1143–8. PubMed Abstract | Publisher Full Text
47. Eser B, Cora T, Eser O, et al.: Association of the polymorphisms of vitamin D receptor and aggrecan genes with degenerative disc disease. Genet Test Mol Biomark. 2010; 14(3): 313–7. PubMed Abstract | Publisher Full Text
48. Li L, Ni D, Zhu F, et al.: No association between VDR gene polymorphisms and lumbar disc herniation in a Chinese population. Int J Clin Exp Med. 2018; 11: 1009–14.
49. Oishi Y, Shimizu K, Katoh T, et al.: Lack of association between lumbar disc degeneration and osteophyte formation in elderly Japanese women with back pain. Bone. 2003; 32(4): 405–11. PubMed Abstract | Publisher Full Text
50. Omair A, Lie BA, Reikeras O, et al.: An association study of interleukin 18 receptor genes (IL18R1 and IL18RAP) in lumbar disc degeneration. Open Orthop J. 2012; 6: 164. PubMed Abstract | Publisher Full Text | Free Full Text
51. Yuan H-Y, Tang Y, Liang Y-X, et al.: Matrix metalloproteinase-3 and vitamin d receptor genetic polymorphisms, and their interactions with occupational exposure in lumbar disc degeneration. J Occup Health. 2010; 52(1): 23–30. PubMed Abstract | Publisher Full Text
52. Su H, Rei N, Zhang L, et al.: Meta-analyses of IL1A polymorphisms and the risk of several autoimmune diseases published in databases. PLoS One. 2018; 13: e0198693.
53. Wang H, Wang Y: Association between interleukin-1 alpha rs1800587 polymorphism and rheumatoid arthritis: A meta-analysis. Sci Prog. 2021; 104: 368504211013647.
54. Appelkvist EL: Dolichol biosynthesis in rat liver peroxisomes. Acta Chem Scand B. 1987; 41: 73–5.
55. Liu H, He H, Li S, et al.: Vitamin D receptor gene polymorphisms and risk of osteoarthritis: a meta-analysis. Exp Biol Med (Maywood) 2014; 239: 559–67.
56. Raj PP: Intervertebral disc: anatomy-physiology-pathophysiology-treatment. Pain Pract. 2008 Jan-Feb; 8(1): 18–44. PubMed Abstract | Publisher Full Text
57. Cervin Serrano S, González Villareal D, Aguilar-Medina M, et al.: Genetic polymorphisms of interleukin-1 alpha and the vitamin d receptor in mexican mestizo patients with intervertebral disc degeneration. Int J Genomics. 2014; 2014: 302568. PubMed Abstract | Publisher Full Text | Free Full Text
58. Mayer JE, Iatridis JC, Chan D, et al.: Genetic polymorphisms associated with intervertebral disc degeneration. Spine J. 2013 Mar; 13(3): 299–317. PubMed Abstract | Publisher Full Text | Free Full Text
59. Ishihara H, Urban JP: Effects of low oxygen concentrations and metabolic inhibitors on proteoglycan and protein synthesis rates in the intervertebral disc. J Orthop Res. 1999 Nov; 17(6): 829–35. PubMed Abstract | Publisher Full Text
60. Ohshima H, Urban JP: The effect of lactate and pH on proteoglycan and protein synthesis rates in the intervertebral disc. Spine. 1992 Sep; 17(9): 1079–82. PubMed Abstract | Publisher Full Text
61. Akyol S, Eraslan BS, Etyemez H, et al.: Catabolic cytokine expressions in patients with degenerative disc disease. Turk Neurosurg. 2010 Oct; 20(4): 492–9. PubMed Abstract | Publisher Full Text
62. Grönblad M, Virri J, Tolonen J, et al.: A controlled immunohistochemical study of inflammatory cells in disc herniation tissue. Spine (Phila Pa 1976). 1994 Dec 15; 19(24): 2744–51. PubMed Abstract | Publisher Full Text
63. Rannou F, Corvol MT, Hudry C, et al.: Sensitivity of anulus fibrosus cells to interleukin 1 beta. Comparison with articular chondrocytes. Spine (Phila Pa 1976). 2000 Jan; 25(1): 17–23. PubMed Abstract | Publisher Full Text
64. Hong J, Yan J, Chen J, et al.: Identification of key potential targets for TNF-α/TNFR1-related intervertebral disc degeneration by bioinformatics analysis. Connect Tissue Res. 2020 Aug 26: 1–11. PubMed Abstract | Publisher Full Text
65. Le Maitre CL, Hoyland JA, Freemont AJ: Catabolic cytokine expression in degenerate and herniated human intervertebral discs: IL-1beta and TNFalpha expression profile. Arthritis Res Ther. 2007; 9(4): R77. PubMed Abstract | Publisher Full Text | Free Full Text
66. Le Maitre CL, Freemont AJ, Hoyland JA: A preliminary in vitro study into the use of IL-1Ra gene therapy for the inhibition of intervertebral disc degeneration. Int J Exp Pathol. 2006 Feb; 87(1): 17–28. PubMed Abstract | Publisher Full Text | Free Full Text
67. Maroudas A: Biophysical chemistry of cartilaginous tissues with special reference to solute and fluid transport. Biorheology. 1975 Jun; 12(3–4): 233–48. PubMed Abstract | Publisher Full Text
68. El-Omar EM, Carrington M, Chow WH, et al.: Interleukin-1 polymorphisms associated with increased risk of gastric cancer. Nature. 2000 Mar 23; 404(6776): 398–402. PubMed Abstract | Publisher Full Text
69. Nee LE, McMorrow T, Campbell E, et al.: TNF-alpha and IL-1beta-mediated regulation of MMP-9 and TIMP-1 in renal proximal tubular cells. Kidney Int. 2004 Oct; 66(4): 1376–86. PubMed Abstract | Publisher Full Text
70. Shingu M, Nagai Y, Isayama T, et al.: The effects of cytokines on metalloproteinase inhibitors (TIMP) and collagenase production by human chondrocytes and TIMP production by synovial cells and endothelial cells. Clin Exp Immunol. 1993 Oct; 94(1): 145–9. PubMed Abstract | Publisher Full Text | Free Full Text
71. Wisithphrom K, Murray PE, Windsor LJ: Interleukin-1 alpha alters the expression of matrix metalloproteinases and collagen degradation by pulp fibroblasts. J Endod. 2006 Mar; 32(3): 186–92. PubMed Abstract | Publisher Full Text
72. Fournier DE, Kiser PK, Shoemaker JK, et al.: Vascularization of the human intervertebral disc: A scoping review. JOR Spine. 2020; 3(4): e1123. PubMed Abstract | Publisher Full Text | Free Full Text
73. Urban JP, Smith S, Fairbank JC: Nutrition of the intervertebral disc. Spine. 2004 Dec 1; 29(23): 2700–9. PubMed Abstract | Publisher Full Text
74. Roberts S, Urban JP, Evans H, et al.: Transport properties of the human cartilage endplate in relation to its composition and calcification. Spine. 1996 Feb 15; 21(4): 415–20. PubMed Abstract | Publisher Full Text
75. Kongsbak M, Levring T, Geisler C, et al.: The Vitamin D Receptor and T Cell Function. Front Immunol. 2013 2013-June-18; 4(148). PubMed Abstract | Publisher Full Text | Free Full Text
76. Videman T, Leppävuori J, Kaprio J, et al.: Intragenic polymorphisms of the vitamin D receptor gene associated with intervertebral disc degeneration. Spine. 1998 Dec 1; 23(23): 2477–85. PubMed Abstract | Publisher Full Text
77. Harapan H: The role of single nucleotide polymorphisms of IL-1A -889C>T (rs1800587), TNF-A -238G>A (rs361525), and VDR TaqI (rs731236) on susceptivity of herniated nucleus pulposus. figshare. Figure. 2020. Publisher Full Text

Comments on this article Comments (0)

Version 3

VERSION 3 PUBLISHED 25 May 2021