Sequence variations in the collagen IX and XI genes are associated with degenerative lumbar spinal stenosis

A novel missense COL9A3 variant in a pedigree with multiple lumbar disc herniation

Trp3 allele in COL9A3 gene has been widely studied in populations with intervertebral disc disease. We identified a novel pathogenic variant in COL9A3 gene in a pedigree with multiple lumbar disc herniation (LDH). The proband was a 14-year-old boy who developed LDH at the L4/5 and L5/S1 spinal segments. His father, paternal aunt and grandfather were diagnosed with LDH at an age of 35, 30 and 23, respectively. By applying whole exome sequencing, a heterozygous missense variant (c.1150C > T, p.Arg384Trp) in COL9A3 was identified. According to the ACMG guidelines, this variant is predicted to be pathogenic. In addition, prediction tools found COL9A3 protein of this variant a reduced stability, some changed charge properties, and an altered spatial conformation. Findings expanded the mutational spectrum of LDH and contributed to the understanding of COL9A3 in the pathogenesis of LDH.

Biomarkers for intervertebral disc and associated back pain: From diagnosis to disease prognosis and personalized treatment

Biomarkers are commonly recognized as objective indicators of a medical state or clinical outcome and have been widely used as clinical and diagnostic tools and surrogate endpoints in many pathological conditions. In the context of intervertebral disc (IVD) and associated back pain, also known as degenerative disc disease (DDD), the use of biomarkers has been poorly explored. DDD is currently diagnosed using imaging techniques and subjective pain scales, limiting an objective association between DDD and pain levels, as well as an evaluation of disease progression. There is a need for objective and reliable measurements for DDD, pain and pathology progression. DDD predictors could also help clinicians in deciding on the optimal treatment for distinct patient groups. This review addresses the current candidate biomarkers in DDD, including imaging, genetic, metabolite and protein-based parameters, both at the tissue and systemic levels, that may become a major advance in the diagnosis and prognosis of the disease, as well as in the management of therapeutic approaches to DDD.

Disruption of Col9a2 expression leads to defects in osteochondral homeostasis and osteoarthritis-like phenotype in mice

Background/Objective

As one of the branched chains of Type IX collagen (Col9), Collagen IX alpha2 (Col9a2) has been reported to be associated with several orthopedic conditions. However, the relationship between Col9a2 and knee osteoarthritis (KOA) remains to be elucidated.

Methods

To probe the relationship between Col9a2 and KOA, we performed a systematic analysis of Col9a2-deficient (Col9a2-/-) mice using whole-mount skeletal staining, Micro-CT (μCT), biomechanics, histomorphometry, immunohistochemistry (IHC), immunofluorescence (IF) and Elisa.

Results

We found that the subchondral bone (SCB) in the knee joint of Col9a2-/- mice became sparse and deformed in the early stage, with altered bone morphometric parameters, reduced load-bearing capacity, dysfunctional bone homeostasis (decreased osteogenesis capacity and elevated bone resorption capacity), diminished cartilage proteoglycans and disrupted cartilage extracellular matrix (ECM) anabolism and catabolism compared with the Col9a2+/+ mice. In the late stage, the cartilage degeneration in Col9a2-/- mice were particularly pronounced compared to Col9a2+/+ mice, as evidenced by severe cartilage destruction and a marked reduction in cartilage thickness and area.

Conclusion

Overall, Col9a2 is essential for maintaining osteochondral homeostasis in the knee joint of mice, and the absence of this gene is accompanied by distinct sclerosis of the SCB and a reduction in load-bearing capacity; in the late stage, in the lack of SCB stress inhibition, excessive load is consistently exerted on the cartilage, ultimately leading to osteoarthritic-like articular cartilage damage.

Identification of Novel Genetic Markers for the Risk of Spinal Pathologies: A Genome-Wide Association Study of 2 Biobanks

BACKGROUND

Identifying genetic risk factors for spinal disorders may lead to knowledge regarding underlying molecular mechanisms and the development of new treatments.

METHODS

Cases of lumbar spondylolisthesis, spinal stenosis, degenerative disc disease, and pseudarthrosis after spinal fusion were identified from the UK Biobank. Controls were patients without the diagnosis. Whole-genome regressions were used to test for genetic variants potentially implicated in the occurrence of each phenotype. External validation was performed in FinnGen.

RESULTS

A total of 389,413 participants were identified from the UK Biobank. A locus on chromosome 2 spanning GFPT1, NFU1, AAK1, and LOC124906020 was implicated in lumbar spondylolisthesis. Two loci on chromosomes 2 and 12 spanning genes GFPT1, NFU1, and PDE3A were implicated in spinal stenosis. Three loci on chromosomes 6, 10, and 15 spanning genes CHST3, LOC102723493, and SMAD3 were implicated in degenerative disc disease. Finally, 2 novel loci on chromosomes 5 and 9, with the latter corresponding to the LOC105376270 gene, were implicated in pseudarthrosis. Some of these variants associated with spinal stenosis and degenerative disc disease were also replicated in FinnGen.

CONCLUSIONS

This study revealed nucleotide variations in select genetic loci that were potentially implicated in 4 different spinal pathologies, providing potential insights into the pathological mechanisms.

LEVEL OF EVIDENCE

Prognostic Level III. See Instructions for Authors for a complete description of levels of evidence.

Assessment of the COL9A1 Single Nucleotide Polymorphism with Severity of clubfoot in a paediatric population along with their biological mothers

Background

A frequent birth abnormality known as congenital talipes equinovarus (CTEV) results in severe abnormalities of one or both feet. The risk of CTEV has been shown to be significantly influenced by genetics. It is a typical congenital malformation. Diseases involving the articular cartilage and COL9A1 polymorphisms are linked.

Materials and method

This case parent dyad research was conducted in a premier care medical and educational facility by the Departments of Paediatric Orthopedics and Biochemistry. Children who had been diagnosed with clubfoot as well as their biological mothers took part in the study. 125 kids were able to join in the research when all the screening, inclusion, and exclusion criteria were met. Baseline demographic information was collected, including the child's age and sex, the kind of clubfoot, any family members who have the condition, and whether or not the mothers smoke or drink. Pirani score is used to clinically evaluate every case. Only one peripheral blood sample was taken from each patient, including their Biological Mothers.

Results

out of 125 children enrolled with biological mothers, Col9A1 SNP rs1135056 is substantially related.

Conclusion

Additionally, patients with the GG genotype for rs592121 have a higher chance of developing CTEV than those with other genotypes. In this investigation, we found possible associations between COL9A1 gene polymorphisms in the mother and offspring with the risk of CTEV. Our research may help us comprehend the genetic makeup of CTEV better and lay the groundwork for creative intervention strategies.

FokI polymorphism in the vitamin D receptor gene in patients with hip osteoarthritis: A case-control study

Objectives

This study aimed to examine whether there is an association between hip osteoarthritis and vitamin D receptor gene FokI polymorphisms.

Patients and methods

In this case-control study, a total of 162 volunteers (43 males, 119 females; mean age: 62.4±9.3; range, 50 to 80 years) were included between March 2011 and March 2012. The patient group included 80 individuals with a diagnosis of coxarthrosis. Eighty-two individuals with normal hip, low back, and sacroiliac joint examination were included in the control group. The American College of Rheumatology hip osteoarthritis classification criteria were used in the diagnosis of coxarthrosis. Analysis of FokI polymorphisms was performed with restriction fragment length polymorphisms.

Results

When the genotype and allele distributions of vitamin D receptor FokI polymorphisms were examined, no statistical difference was found in both groups.

Conclusion

No significant association between the hip osteoarthritis and FokI polymorphisms could be obtained.

Targeting Autophagy for Developing New Therapeutic Strategy in Intervertebral Disc Degeneration

Intervertebral disc degeneration (IVDD) is a prevalent cause of low back pain. IVDD is characterized by abnormal expression of extracellular matrix components such as collagen and aggrecan. In addition, it results in dysfunctional growth, senescence, and death of intervertebral cells. The biological pathways involved in the development and progression of IVDD are not fully understood. Therefore, a better understanding of the molecular mechanisms underlying IVDD could aid in the development of strategies for prevention and treatment. Autophagy is a cellular process that removes damaged proteins and dysfunctional organelles, and its dysfunction is linked to a variety of diseases, including IVDD and osteoarthritis. In this review, we describe recent research findings on the role of autophagy in IVDD pathogenesis and highlight autophagy-targeting molecules which can be exploited to treat IVDD. Many studies exhibit that autophagy protects against and postpones disc degeneration. Further research is needed to determine whether autophagy is required for cell integrity in intervertebral discs and to establish autophagy as a viable therapeutic target for IVDD.

Immuno-Modulatory Effects of Intervertebral Disc Cells

Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.

Prevalence and Definition of Multilevel Lumbar Developmental Spinal Stenosis

STUDY DESIGN

Cross-sectional study.

OBJECTIVES

To define multilevel lumbar developmental spinal stenosis (DSS) using a composite score model and to determine its prevalence.

METHODS

This was a cohort study of 2385 openly recruited subjects with lumbosacral (L1-S1) MRIs. All subjects with previous spinal surgery or spinal deformities were excluded. The anteroposterior (AP) vertebral canal diameter was measured by two independent observers. Any associations between level-specific vertebral canal diameter and subject body habitus were analysed with non-parametric tests. Three or more stenotic levels, equivalent to a composite score of 3 or more, were considered as multilevel DSS. The median values of these subjects' AP canal diameters were used to construct the multilevel DSS values. Receiver operating characteristic analysis was utilized to determine the ability of these cut-off values to screen for DSS by presenting their area under curve, sensitivity and specificity.

RESULTS

Subject body habitus was poorly correlated with AP vertebral canal diameter. Multilevel DSS was identified as L1<19 mm, L2<19 mm, L3<18 mm, L4<18 mm, L5<18 mm, S1<16 mm with 81%-96% sensitivity and 72%-91% specificity. The prevalence of multilevel DSS in this cohort was 7.3%.

CONCLUSIONS

Utilizing a large homogeneous cohort, the prevalence of multilevel DSS is determined. Our cut-offs provide high diagnostic accuracy. Patients with multiple levels that fulfil these criteria may be at-risk of spinal canal compressions at multiple sites.

LEVEL OF EVIDENCE

III.

Association between polymorphisms of collagen genes and susceptibility to intervertebral disc degeneration: a meta-analysis

Background

Collagens are important structural components of intervertebral disc. A number of studies have been performed for association between polymorphisms of collagen genes and risk of intervertebral disc degeneration (IVDD) but yielded inconsistent results. Here, we performed a meta-analysis to investigate the association of collagen IX alpha 2 (COL9A2) Trp2, collagen IX alpha 3 (COL9A3) Trp3, collagen I alpha 1 (COL1A1) Sp1 and collagen XI alpha 1 (COL11A1) C4603T polymorphisms with susceptibility to IVDD.

Method

Eligible studies were retrieved by searching MEDLINE, EMBASE, Web of Science prior to 31 March, 2021. Odds ratio (OR) and corresponding 95% confidence interval (CI) were calculated for association strength.

Results

A total of 28 eligible studies (31 datasets comprising 5497 cases and 5335 controls) were included. COL9A2 Trp2 carriers had an increased risk of IVDD than non-carriers in overall population (OR = 1.43, 95% CI 0.99–2.06, P = 0.058), which did not reach statistical significance. However, Trp2 carriers had 2.62-fold (95% CI 1.15–6.01, P = 0.022) risk than non-carriers in Caucasians. COL9A3 Trp3 was not associated with IVDD risk (OR = 1.28, 95% CI 0.81–2.02, P = 0.299). T allele and TT genotype of COL1A1 Sp1 (+ 1245G > T) were correlated with increased risk of IVDD. Significant associations were found between COL11A1 C4603T and IVDD risk under allelic (OR = 1.33, 95% CI 1.20–1.48), dominant (OR = 1.45, 95% CI 1.26–1.67), recessive (OR = 1.55, 95% CI 1.21–1.98) and homozygote model (OR = 1.81, 95% CI 1.40–2.34).

Conclusions

COL1A1 Sp1 and COL11A1 C4603T polymorphism are associated with IVDD risk while the predictive roles of collagen IX gene Trp2/3 need verification in more large-scale studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-021-02724-8.

High-Throughput Gene and Protein Analysis Revealed the Response of Disc Cells to Vitamin D, Depending on the VDR FokI Variants

Vitamin D showed a protective effect on intervertebral disc degeneration (IDD) although conflicting evidence is reported. An explanation could be due to the presence of the FokI functional variant in the vitamin D receptor (VDR), observed as associated with spine pathologies. The present study was aimed at investigating—through high-throughput gene and protein analysis—the response of human disc cells to vitamin D, depending on the VDR FokI variants. The presence of FokI VDR polymorphism was determined in disc cells from patients with discopathy. 1,25(OH)₂D₃ was administered to the cells with or without interleukin 1 beta (IL-1β). Microarray, protein arrays, and multiplex protein analysis were performed. In both FokI genotypes (FF and Ff), vitamin D upregulated metabolic genes of collagen. In FF cells, the hormone promoted the matrix proteins synthesis and a downregulation of enzymes involved in matrix catabolism, whereas Ff cells behaved oppositely. In FF cells, inflammation seems to hamper the synthetic activity mediated by vitamin D. Angiogenic markers were upregulated in FF cells, along with hypertrophic markers, some of them upregulated also in Ff cells after vitamin D treatment. Higher inflammatory protein modulation after vitamin D treatment was observed in inflammatory condition. These findings would help to clarify the clinical potential of vitamin D supplementation in patients affected by IDD.

Transcriptional profiling of intervertebral disc in a post-traumatic early degeneration organ culture model

INTRODUCTION

The goal of this study is to characterize transcriptome changes and gene regulation networks in an organ culture system that mimics early post-traumatic intervertebral disc (IVD) degeneration.

METHODS

To mimic a traumatic insult, bovine caudal IVDs underwent one strike loading. The control group was cultured under physiological loading. At 24 hours after one strike or physiological loading, RNA was extracted from nucleus pulposus (NP) and annulus fibrosus (AF) tissue. High throughput next generation RNA sequencing was performed to identify differentially expressed genes (DEGs) between the one strike loading group and the control group. Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes analyses were performed to analyze DEGs and pathways. Protein-protein interaction (PPI) network was analyzed with cytoscape software. DEGs were verified using qRT-PCR. Degenerated human IVD tissue was collected for immunofluorescence staining to verify the expression of DEGs in human disc tissue.

RESULTS

One strike loading resulted in significant gene expression changes compared with physiological loading. In total 253 DEGs were found in NP tissue and 208 DEGs in AF tissue. Many of the highly dysregulated genes have known functions in disc degeneration and extracellular matrix (ECM) homeostasis. ACTB, ACTG, PFN1, MYL12B in NP tissue and FGF1, SPP1 in AF tissue were verified by qRT-PCR and immunofluorescence imaging. The identified DEGs were involved in focal adhesion, ECM-receptor interaction, PI3K-AKT, and cytokine-cytokine receptor interaction pathways. Three clusters of PPI networks were identified. GO enrichment revealed that these DEGs were mainly involved in inflammatory response, the ECM and growth factor signaling and protein folding biological process.

CONCLUSION

Our study revealed different DEGs, pathways, biological process and PPI networks involved in post-traumatic IVD degeneration. These findings will advance the understanding of the pathogenesis of IVD degeneration, and help to identify novel biomarkers for the disease diagnosis.

The Role of Polymorphisms in Collagen-Encoding Genes in Intervertebral Disc Degeneration

(1) Background: The purpose of this review is to analyze domestic and foreign studies on the role of collagen-encoding genes polymorphism in the development of intervertebral discs (IVDs) degeneration in humans. (2) Methods: We have carried out a search for full-text articles published in e-Library, PubMed, Oxford Press, Clinical Case, Springer, Elsevier and Google Scholar databases. The search was carried out using keywords and their combinations. The search depth was 5 years (2016–2021). In addition, this review includes articles of historical interest. Despite an extensive search, it is possible that we might have missed some studies published in recent years. (3) Results: According to the data of genome-wide and associative genetic studies, the following candidate genes that play a role in the biology of IVDs and the genetic basis of the processes of collagen degeneration of the annulus fibrosus and nucleus pulposus of IVDs in humans are of the greatest interest to researchers: COL1A1, COL2A1, COL9A2, COL9A3, COL11A1 and COL11A2. In addition, the role of genes COL1A2, COL9A1 and others is being actively studied. (4) Conclusions: In our review, we summarized and systematized the available information on the role of genetic factors in IVD collagen fibers turnover and also focused on the functions of different types of collagen present in the IVD. Understanding the etiology of impaired collagen formation can allow doctors to prescribe pathogenetically-based treatment, achieving the most effective results.

Pedigree analysis of lumbar developmental spinal stenosis: Determination of potential inheritance patterns

Lumbar developmental spinal stenosis (DSS) refers to multilevel pre-existing narrowed spinal canals, which predispose to neural compromise. The objective of this study is to identify any inheritance pattern of DSS by utilizing pedigree charts. This was a case series of 13 families with a total of 80 subjects having magnetic resonance imaging (MRI) from L1 to S1. Cases (subjects with DSS) or controls (subjects without DSS) were identified by measuring their anteroposterior (AP) vertebral canal diameters. Multilevel model analyses were also performed to evaluate whether there is substantial clustering of observations within the families, and the effect of multilevel DSS. The intraclass correlation coefficient (ICC) and Akaike information criteria (AIC) were compared between models. Correlations between subject demographics and AP vertebral canal diameter were statistically insignificant at all levels. Only vertebral canal cross-sectional area, and axial and sagittal vertebral canal diameter were found to be statistically different between cases and controls at all levels (all p < .05). Both males and females were affected by DSS and there was no skipping of generation, which highly suggested DSS followed an autosomal dominant inheritance pattern. After accounting for multilevel DSS, there was a drop of more than 10 in AIC and some variances were also explained within families. This is the first study that suggests multilevel lumbar DSS to have an autosomal dominant inheritance pattern. Within families with a background of DSS, subjects had a smaller canal size, contributed by shortened axial and sagittal AP vertebral canal diameter, and smaller canal cross-sectional area.

Association of Single Nucleotide Polymorphism rs2228570 with Lumbar Disc Degeneration: A Case-Control Study and Meta-Analysis

Objective

To examine the association between single nucleotide polymorphisms (SNPs) rs2228570, rs731236, rs7975232, and rs1544410 and lumbar disc degeneration (LDD) predisposition.

Methods

A search strategy was carried out, and the data were extracted after being chosen by the inclusion and exclusion criteria. Pooled odds ratios and 95% confidence intervals were calculated to assess the association between the aforementioned SNPs and LDD under allelic, dominant, recessive, heterozygous, and homozygous genetic models. In addition, a case–control study involving 46 LDD cases and 45 controls was also performed in the analysis to verify the result.

Results

A total of 17 studies were included in this meta-analysis. The pooled results did not show any association between vitamin D receptor (VDR) gene polymorphisms and LDD. But, interestingly, in subgroup analysis, the rs2228570 polymorphism was associated with LDD under the allelic (OR = 0.70, 95% CI = 0.56–0.87, p = 0.002), recessive (OR = 0.60, 95% CI = 0.43–0.84, p = 0.003), and homozygous (OR = 0.47, 95% CI= 0.28–0.79, p = 0.004) genetic models in the Asian population. SNPs rs731236 and rs7975232 still did not show any obvious association. We obtained a similar result from the case–control study: rs2228570 had an obvious relationship with LDD under allelic and homozygous genetic models. At the same time, we found that rs2228570 was also associated with the degree of low back pain (visual analogue scale, VAS score) in LDD population.

Conclusion

SNP rs2228570 was significantly associated with LDD predisposition and the degree of low back pain in the Asian population.

Vitamin D receptor gene polymorphisms and risk of intervertebral disc degeneration: An updated meta-analysis based on 23 studies

BACKGROUND

Numerous studies have investigated the associations between Vitamin D receptor (VDR) gene polymorphisms and risk of intervertebral disc degeneration but the results remain controversial. This study aimed to drive a more precise estimation of association between VDR gene polymorphisms and risk of intervertebral disc degeneration.

METHODS

PubMed, EMBASE, Cochrane library, Web of Science and China Knowledge Resource Integrated Database for papers on VDR gene polymorphisms and risk of intervertebral disc degeneration were searched. The pooled odds ratios (ORs) with 95% confidence intervals (CIs) were used to assess the strength of association in the homozygote model, heterozygote model, dominant model, recessive model and an additive model.

RESULTS

Overall, 23 articles were included in the final meta-analysis. The subgroup analyses by ethnicity showed a significant association of VDR FokI mutation with disc degeneration risk in Caucasians (recessive model, OR with 95%CI 1.301, [1.041, 1.626]; additive model, OR with 95%CI 1.119, [1.006, 1.245]). The results of subgroup analyses by ethnicity showed a significant association of VDR TaqI mutation with disc degeneration risk in Asians but not in Caucasians. There was a significant association between VDR ApaI mutation and risk of disc degeneration and subgroup analyses by ethnicity showed a significant association in Caucasians and in Asians.

CONCLUSIONS

In summary, VDR FokI polymorphisms was associated with disc degeneration risk among Caucasians but not Asians, VDR TaqI polymorphisms was associated with disc degeneration risk among Asians but not Caucasians, while VDR ApaI polymorphism was associated with disc degeneration risk among Asians and Caucasians.

The identification of novel gene mutations for degenerative lumbar spinal stenosis using whole-exome sequencing in a Chinese cohort

Background

Degenerative lumbar spinal stenosis (DLSS) is a common lumbar disease that requires surgery. Previous studies have indicated that genetic mutations are implicated in DLSS. However, studies on specific gene mutations are scarce. Whole-exome sequencing (WES) is a valuable research tool that identifies disease-causing genes and could become an effective strategy to investigate DLSS pathogenesis.

Methods

From January 2016 to December 2017, we recruited 50 unrelated patients with symptoms consistent with DLSS and 25 unrelated healthy controls. We conducted WES and exome data analysis to identify susceptible genes. Allele mutations firstly identified potential DLSS variants in controls to the patients’ group. We conducted a site-based association analysis to identify pathogenic variants using PolyPhen2, SIFT, Mutation Taster, Combined Annotation Dependent Depletion, and Phenolyzer algorithms. Potential variants were further confirmed using manual curation and validated using Sanger sequencing.

Results

In this cohort, the major classification variant was missense_mutation, the major variant type was single nucleotide polymorphism (SNP), and the major single nucleotide variation was C > T. Multiple SNPs in 34 genes were identified when filtered allele mutations in controls to retain only patient mutations. Pathway enrichment analyses revealed that mutated genes were mainly enriched for immune response-related signaling pathways. Using the Novegene database, site-based associations revealed several novel variants, including HLA-DRB1, PARK2, ACTR8, AOAH, BCORL1, MKRN2, NRG4, NUP205 genes, etc., were DLSS related.

Conclusions

Our study revealed that deleterious mutations in several genes might contribute to DLSS etiology. By screening and confirming susceptibility genes using WES, we provided more information on disease pathogenesis. Further WES studies incorporating larger DLSS patient cohorts are required to comprehend the genetic landscape of DLSS pathophysiology fully.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-00981-4.

Vitamin D receptor gene polymorphisms and osteoarthritis: a meta-analysis

OBJECTIVE

To investigate the association of vitamin D receptor (VDR) gene polymorphisms with OA susceptibility.

METHODS

Meta-analyses were performed using allelic contrast, contrast of homozygotes, and recessive and dominant models to clarify the association between OA and VDR ApaI, BsmI, TaqI and FokI polymorphisms. Odds ratio (OR) and the corresponding 95% CI were obtained, and subgroup analyses were performed based on the ethnicity and OA sites.

RESULTS

A total of 18 studies with 2983 OA patients and 4177 controls were included in this meta-analysis. There were statistically significant associations in the spine between OA susceptibility and the VDR BsmI (B vs b: OR = 1.25, 95% CI: 1.03, 1.53, P = 0.026; BB vs bb: OR = 1.56, 95% CI: 1.02, 2.37, P = 0.038) and TaqI (T vs t: OR = 0.73, 95% CI: 0.54, 0.99, P = 0.044; TT vs Tt + tt: OR = 0.63, 95% CI: 0.42, 0.95, P = 0.028) polymorphisms, but not for the other polymorphisms. A statistically significant association was found between the VDR FokI polymorphism and OA susceptibility in the knee in the recessive model contrast (FF vs Ff + ff: OR = 0.63, 95% CI: 0.42, 0.95, P = 0.028), but this result was only pooled from one study. However, no significant associations were found between the VDR ApaI polymorphism and OA. Besides, ethnic stratification also indicated that there was no significant association between VDR gene polymorphism and OA in Caucasians or Asians.

CONCLUSION

Our meta-analysis suggests that the VDR BsmI and TaqI polymorphisms are associated with OA susceptibility in the spine. However, the VDR ApaI polymorphism is not a significant genetic risk factor for OA.

Directed Differentiation of Notochord-like and Nucleus Pulposus-like Cells Using Human Pluripotent Stem Cells

Intervertebral disc degeneration might be amenable to stem cell therapy, but the required cells are scarce. Here, we report the development of a protocol for directed in vitro differentiation of human pluripotent stem cells (hPSCs) into notochord-like and nucleus pulposus (NP)-like cells of the disc. The first step combines enhancement of ACTIVIN/NODAL and WNT and inhibition of BMP pathways. By day 5 of differentiation, hPSC-derived cells express notochordal cell characteristic genes. After activating the TGF-β pathway for an additional 15 days, qPCR, immunostaining, and transcriptome data show that a wide array of NP markers are expressed. Transcriptomically, the in vitro-derived cells become more like in vivo adolescent human NP cells, driven by a set of influential genes enriched with motifs bound by BRACHYURY and FOXA2, consistent with an NP cell-like identity. Transplantation of these NP-like cells attenuates fibrotic changes in a rat disc injury model of disc degeneration.

High Rate of Studies with Level 1 and 2 Evidence among the 100 Most Cited Articles in Lumbar Spinal Stenosis

BACKGROUND

 To date, no study has used bibliometric analysis to review the most influential articles in lumbar spinal stenosis. The objective of this study was to identify and analyze the characteristics and the level of evidence of the 100 most cited articles on lumbar spinal stenosis METHODS:  The Thomson Reuters Web of Science was accessed to find the 100 most cited articles on lumbar spinal stenosis. For each article, we recorded the number and density of citations, authors, country, journals and years, department, level of evidence, type of study, and if it was part of any multicenter studies.

RESULTS

 Until January 2017, the 100 most cited articles accumulated 11,136 citations (average: 259.05/y), ranging individually between 442 and 50 (average: 111.36). The first reference was published in 1974 in Clinical Orthopaedics and Related Research. Therapeutic studies (n = 40), the 1990s (n = 46), United States as country of origin (n = 51), Harvard University as institution (n = 16), Katz JN as author (n = 10), and Spine as journal (n = 48) have the hegemony. Many were multicenter (n = 42) and using level 2 evidence (n = 49). There is an inverse relationship between citation index and long-standing studies, maintenance of those most cited, and a temporary advance toward better levels of evidence.

CONCLUSION

 This bibliometric analysis reveals a good level of evidence in the published clinical series and includes 100 articles useful for the approach of lumbar spinal stenosis.

Integrated transcriptome and proteome analyses identify novel regulatory network of nucleus pulposus cells in intervertebral disc degeneration

Background

Degeneration of intervertebral disc is a major cause of lower back pain and neck pain. Studies have tried to unveil the regulatory network using either transcriptomic or proteomic analysis. However, neither have fully elucidated the exact mechanism of degeneration process. Since post-transcriptional regulation may affect gene expression by modulating the translational process of mRNA to protein product, a combined transcriptomic and proteomic study may provide more insight into the key regulatory network of Intervertebral disc degeneration.

Methods

In order to obtain the proteomic and transcriptomic data, we performed label-free proteome analysis on freshly isolated nucleus pulposus cells and obtained transcriptome profiling data from the Gene Expression Omnibus repository. To identify the key regulatory network of intervertebral disc degeneration in nucleus pulposus cells, we performed bioinformatic analyses and established a protein-RNA interacting network. To validate the candidate genes, we performed in vitro experimentation and immunochemistry labeling to identify their potential function during nucleus pulposus degeneration.

Results

The label-free proteome analysis identified altogether 656 proteins, and 503 of which were differentially expressed between nucleus pulposus cells from degenerated or normal disc cells. Using the existing nucleus pulposus transcriptomic profiling data, we integrated the proteomic and transcriptomic data of nucleus pulposus cells, and established a protein-RNA interacting network to show the combined regulatory network of intervertebral disc degeneration. In the network, we found 9 genes showed significant changes, and 6 of which (CHI3L1, KRT19, COL6A2, DPT, TNFAIP6 and COL11A2) showed concordant changes in both protein and mRNA level. Further functional analysis showed these candidates can significantly affect the degeneration of the nucleus pulposus cell when altering their expression.

Conclusions

This study is the first to use combined analysis of proteomic and transcriptomic profiling data to identify novel regulatory network of nucleus pulposus cells in intervertebral disc degeneration. Our established protein-RNA interacting network demonstrated novel regulatory mechanisms and key genes that may play vital roles in the pathogenesis of intervertebral disc degeneration.

Associations between vitamin D receptor gene polymorphisms and spinal degenerative disease: evidence from a meta-analysis based on 35 case-control studies

OBJECTIVE

Dozens of reports on the associations of vitamin D receptor (VDR) gene polymorphisms and susceptibility to spinal degenerative disease (SDD) were conducted with inconsistent findings. This study aimed to elucidate the associations through a meta-analysis approach.

METHODS

Databases of PubMed, EMBASE, Web of Science, CNKI, and Wanfang were searched until July 10, 2020. Study quality was evaluated by using Newcastle-Ottawa Scale (NOS). Odds ratios (ORs) and 95% confidence intervals (95%CIs) were calculated to evaluate the associations under allelic model (1 vs. 2), homozygous model (11 vs. 22), heterozygous model (12 vs. 22), dominant model (11 + 12 vs. 22), and recessive model (11 vs. 12 + 22).

RESULTS

A total of 5021 cases and 5746 controls from 35 studies were eligible to this meta-analysis. According to NOS, the included studies were in excellent quality. In the overall population, the pooled data indicated that ApaI was associated with a reduced SDD susceptibility (AA vs. Aa + aa, OR = 0.83, 95%CI 0.71 - 0.96, P = 0.010). But the association was not observed in FokI, TaqI, and BsmI polymorphisms. Subgroup analysis suggested that TaqI polymorphism was correlated to an elevated SDD risk in Asians (TT + Tt vs. tt, OR = 2.55, 95%CI 1.90 - 3.44, P < 0.001).

CONCLUSION

The present study indicates that ApaI polymorphism may contribute to a reduced risk to SDD in the overall population, and TaqI polymorphism confers an elevated susceptibility to SDD in Asians. While, BsmI and FokI polymorphisms appear to have no significant association with SDD.

Genetic risk factors for lumbar disc disease

AIM AND BACKGROUND

Lumbar disc degeneration (LDD) is thought to be multifactorial in origin. Very recently the focus has shifted to the involvement of a family of candidate genes in the pathogenesis of LDD. There is particular emphasis on the vitamin D receptor gene (VDR gene). The VDR polymorphisms FOK1, TAQ1, and APO1 have been variably associated with LDD.

OBJECTIVE

To evaluate the association between the FOK1/Taq1 genes and LDD.

MATERIALS AND METHODS

One hundred unrelated healthy (asymptomatic) individuals who presented for routine health checkup and 93 consecutive patients (43 males and 50 females) with no history of low back pain were enrolled in the study after informed consent was obtained. The MRI images of cases and controls were graded and peripheral blood samples were collected from all participants and sent for genetic analysis.

RESULTS

Individuals with the dominant genotype for Taq1 had a significantly higher association with LDD than those without it. There was no association between LDD and the Fok1 genotype.

CONCLUSION

Genetic predisposition is an important risk factor for LDD.

Genetic variant of COL11A2 gene is functionally associated with developmental dysplasia of the hip in Chinese Han population

OBJECTIVES

Developmental dysplasia of the hip (DDH) is a common skeletal disorder. This study was conducted to demonstrate the association between DDH and a polymorphism rs9277935 of COL11A2 gene.

RESULTS

A significant difference in genotype distribution in a recessive model (TT+GT vs. GG) between two groups (P=0.017) was demonstrated. Analysis in female patients showed significantly greater frequency of minor allele G(0.49 vs. 0.43, p=0.024) and significantly higher distribution of GG genotype (p=0.006). DDH patients were found to have significantly lower COL11A2 expression than controls. Moreover, DDH patients with rs9277935 genotype TT have a significantly increased expression of COL11A2 than those with genotype GG. COL11A2 demonstrated chondrogenic properties in vitro.

CONCLUSION

Polymorphism rs9277935 of gene COL11A2 is a functional variant regulating the expression and the chondrogenic properties of COL11A2 in DDH in Chinese Han population.

METHODS

A case-control candidate gene association study was conducted in 945 patients (350 radiologically confirmed DDH patients and 595 healthy controls). Difference of COL11A2 expression in hip joint tissue was compared between the patients and the controls. Allelic difference in Col11a2 expression by rs9277935 was assessed with luciferase activity. Chondrogenic effects of Col11a2 signaling on BMSCs were also determined in vitro.

Association of vitamin D receptor gene polymorphisms with disc degeneration

PURPOSE

Numerous candidate genes and single-nucleotide polymorphisms (SNPs) have been identified in the background of lumbar disc degeneration (LDD). However, in most of these underpowered studies, definitions of LDD are inconsistent; moreover, many of the findings have not been replicated and are contradictory. Our aim was to characterize LDD by well-defined phenotypes and possible endophenotypes and analyse the association between these and candidate vitamin D receptor (VDR) gene polymorphisms on a large (N = 1426) dataset.

METHODS

Seven candidate VDR SNPs were genotyped. Individual association, haplotype and gene-gene interaction analyses were performed. All degenerative endophenotypes were significantly associated with one or more candidate VDR gene variants.

RESULTS

Haplotype analyses confirmed the association between the 3'-end VDR variants (BsmI, ApaI, TaqI) and Modic changes as well as the relationship of 5'-end variants (Cdx2, A1012G) with endplate defects. We also found significant interactions between the 3'- and 5'-end regulatory regions and endplate defects. Based on our results, VDR and its gene variants are highly associated with specific degenerative LDD endophenotypes.

CONCLUSION

Understanding relationships between phenotype and gene variants is crucial for describing the pathways leading to the multifactorial, polygenic degeneration process and LDD-related conditions. These slides can be retrieved under Electronic Supplementary Material.

Vitamin D Receptor gene polymorphisms and plasma levels are associated with lumbar disc degeneration

The purpose of this study was to investigate the association of Vitamin D Receptor (VDR) gene polymorphisms and VDR levels with lumbar disc degeneration (LDD). TaqMan SNP Genotyping Assay was utilized to probe VDR gene polymorphisms including the FokI (rs2228570), ApaI (rs7975232) and TaqI (rs731236) in 454 patients with LDD and 485 controls. Enzyme-Linked Immunosorbent Assay (ELISA) was used to detect plasma VDR levels. The patients with LDD were divided into three subgroups (subgroup 1: lumbar disc herniation; subgroup 2: lumbar spinal stenosis; subgroup 3: lumbar spondylolisthesis) to further probe the association of plasma VDR levels and VDR gene polymorphisms and LDD. Moreover, immunohistochemistry (IHC) was implemented to evaluate VDR expression in lumbar degenerated disc and normal disc. Allele and genotype frequency of TaqI (rs731236) were significantly different in patients with LDD and controls (all P < 0.05). For TaqI polymorphism, the frequencies of T allele were significantly higher in the LDD patients compared with controls (OR = 1.319; 95%CI 1.091 to 1.595; P = 0.004, adjusted (OR = 1.319; 95%CI 1.091 to 1.595; P = 0.004, adjusted OR = 1.383; 95%CI 1.135 to 1.684; P = 0.016). Furthermore, the allele distribution showed a higher frequency of the T allele in the patients with lumbar disc herniation in subgroup 1 (OR = 1.384; 95% CI 1.105 to 1.732; P = 0.004, adjusted OR = 1.319; 95%CI 1.091 to 1.595; P = 0.016). Plasma VDR levels and VDR expression were significantly lower in patients with LDD compared with controls (all P < 0.05). Moreover, the TT genotype of TaqI polymorphism was significantly associated with lower plasma VDR levels in patients with LDD (P = 0.002). TaqI (rs731236) polymorphism was associated with a predisposition to LDD. Plasma VDR and VDR expression levels may be the marker for the occurrence and development of LDD.

Genetic variants in COL11A2 of lumbar disk degeneration among Chinese Han population

BACKGROUND

Lumbar disk disease (LDD) is a common musculoskeletal disorder. Several predisposing genetic and environmental risk factors have been established for symptomatic LDD.

METHODS

We conducted a case-control association study to investigate the role of the COL11A2 gene in LDD. Genotyping of 384 Chinese Han LDD patients and 384 Chinese Han controls was made for six single-nucleotide polymorphisms (SNPs) from COL11A2 by Agena Massarray. We evaluated these SNPs association with LDD using the chi-square test and genetic model analysis.

RESULTS

The strongest associations with LDD were observed for polymorphisms in rs2071025. Carriers of "A" allele had an increased risk of LDD (OR = 1.47, 95% CI = 1.20-1.80, p = 0.0002) as compared with the "G" allele in allele model. We found that rs2071025 were associated with LDD in female and male from the stratification analyses (p < 0.05). Genetic models showed that rs986522(C) significantly increased the risk of LDD in female; however, in males, we did not find significant associations between the rs986522 and LDD risk.

CONCLUSION

This study showed a genetic association with COL11A2 polymorphism in individuals with LDD. These data may provide novel insights into the pathogenesis of LDD, although further studies with larger numbers of participants worldwide are needed for validation of our conclusions.

Serum 25-hydroxyvitamin D, serum calcium and vitamin D receptor (VDR) polymorphisms in a selected population with lumbar disc herniation-A case control study

Background

Association of Vitamin D receptor (VDR) polymorphisms with lumbar disc herniation (LDH) have been identified in several ethnic groups globally. Despite abundant sunlight, vitamin D deficiency is reported in many tropical countries. As vitamin D is a key modulator for intestinal calcium absorption, low vitamin D could contribute to low serum calcium leading to abnormalities of skeletal homeostasis. Therefore, present study was aimed to study the association of serum 25-hydroxyvitamin D (25(OH)D), serum calcium and VDR polymorphisms in a selected Sri Lankan population.

Materials & methods

A case control study was conducted in 119 participants (cases = 51: controls = 68). Serum 25(OH)D levels were measured using ELISA. The VDR polymorphisms (Fok I and Taq I) were detected by polymerase chain reaction followed by restriction fragment length polymorphism.

Results

Findings indicated a significantly low (p = 0.000) 25(OH)D levels in cases (18.7±3.7 ng/mL) compared to controls(25.5±9.8 ng/mL) while 25(OH)D in both groups were below the reference range. Mean serum calcium levels in both groups were within normal reference range and was not significantly different among groups. Statistically significant association was not observed between VDR Fok I polymorphisms among cases and controls. Although Fok I polymorphism genotypes were in Hardy-Weinberg equilibrium (HWE), Taq I genotypes in controls violated HWE.

Conclusion

Present study confirms that insufficient serum 25(OH)D levels in cases have major contribution to LDH. VDR Fok I polymorphisms did not have any significant association with LDH in Sri Lankan ethnicity.

FokI as a genetic factor of intervertebral disc degeneration: A PRISMA-compliant systematic review of overlapping meta-analyses

The association of FokI (rs2228570), a polymorphism of the vitamin D receptor gene, with intervertebral disc degeneration (IDD) has been investigated in a multitude of studies. However, conflicting results of these studies led to emergence of several meta-analyses over the past few years. Despite the increased statistical power, these meta-analyses have failed to provide uniform and conclusive data on the relationship of FokI with IDD. The aim of this study was to present a comprehensive review based on the most up-to-date meta-analyses on the association of FokI with IDD. A comprehensive search of all major databases was conducted to identify meta-analyses investigating relation between FokI and IDD. No date or language restrictions were applied. The Jadad decision algorithm was utilized to evaluate included meta-analyses and identify the one providing the best evidence. A total of 7 meta-analyses (n = 2580 original patients), that included six to ten case control studies, analyzed the association of FokI polymorphism with IDD. The meta-analysis of the highest quality supported the notion that overall there is no statistically significant association between FokI polymorphism and IDD. However, the authors showed that Caucasians have a reduced risk of IDD and Hispanics have an increased risk of IDD in the dominant and dominant/homozygous/heterozygous models of FokI polymorphism. While currently there is no evidence of an association between FokI polymorphism and IDD in the general population, ethnic predisposition has been shown.

Genetic predictors of recovery in low back and lumbar radicular pain

Previous data suggest that persistent back pain may be associated with genetic variability. In this study, we assessed the correlation between 8 genetic polymorphisms (VDR, COL11, MMP1, MMP9, IL-1α, IL-1RN, OPRM1, COMT) and pain recovery in patients with low back pain (LBP) and lumbar radicular pain (LRP). In total, 296 patients with LBP or LRP were followed for 5 years. The patients underwent standardized clinical examination and completed pain and function questionnaires. Univariate linear regression associations with P values <0.1 were included in the multivariable analysis, adjusting for pain intensity at baseline, age, sex, smoking, body mass index, and LBP or LRP. Pain intensity at 5-year follow-up was associated with VDR rs731236 (B = -0.5, 95% confidence interval [CI] -0.9 to -0.1, P = 0.017), MMP9 rs17576 (B = 0.5, 95% CI 0.1-0.9, P = 0.022), and OPRM1 rs1799971 (B = -0.8, 95% CI -1.4 to -0.2, P = 0.006) in the univariate analyses. MMP9 rs17576 and OPRM1 rs1799971 remained significant (B = 0.4, 95% CI 0.05-0.8, P = 0.026 and B = -0.8, 95% CI -1.3 to -0.2, P = 0.007) in the multivariable model. Thus, the data demonstrated that the rare allele of MMP9 rs17576 was associated with poor pain recovery, whereas the rare allele of OPRM1 rs1799971 was associated with better pain recovery at 5-year follow-up in the LBP and LRP patients. In particular, the present study suggested that the OPRM1 rs179971 A>G in men was associated with better long-term pain recovery. In men, the OPRM1 rs1799971 explained 4.7% of the variance of pain intensity. We conclude that the MMP9 rs17576 and OPRM1 rs1799971 genotypes may affect 5-year recovery in patients with LBP and LRP.

Etiology of developmental spinal stenosis: A genome-wide association study

Our study aimed to identify possible single nucleotide polymorphisms (SNPs) via a genome-wide association study (GWAS) approach and a candidate gene platform that were associated with lumbar developmental spinal stenosis (DSS). Southern Chinese population-based study volunteers were assessed (age range: 18-55 years). DSS was defined as the anteroposterior bony spinal canal diameter on T1-weighted axial MRI of L1 to S1. Genotyping was performed using the Illumina HumanOmniZhongHua-8 BeadChip. Using the canal diameter as the quantitative trait, genomic statistical analyses was performed. A total of 469 subjects were recruited. The mean axial AP measurements noted were: L1: 21.8 mm, L2: 21.9 mm, L3: 22.4 mm, L4: 20.2 mm, L5: 19.6 mm, and S1: 17.3 mm. Q-Q plots of genome-wide associations found significant differences in L4 and L5 measurements. More significant SNPs were found on chromosomes 8, 11, and 18. Low-density lipoprotein receptor-related protein 5 on chromosome 11 was found to be an important functional gene in canal bony development via candidate gene approach. We found two clusters in the findings with one including the upper levels (L1-L4) and the other the lower levels (L5 and S1). This is the first GWAS addressing DSS. The presence of multiple SNPs suggests a multi-factorial origin of DSS. Further analyses noted region-specific genetic predisposition, delineating distinct upper to lower lumbar regions of DSS. With better understanding of the DSS phenotype and genetic markers, the at-risk population can be identified early, preventative measures can be initiated, lifestyle/activity modification can be implemented, and more novel and precision-based therapeutics can be developed. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1262-1268, 2018.

Collagen IX gene polymorphisms and lumbar disc degeneration: a systematic review and meta-analysis

Background

An increasing number of studies have investigated associations between collagen IX alpha 2 chain (COL9A2) and collagen IX alpha 3 chain (COL9A3) gene polymorphisms and the risk of lumbar disc degeneration (LDD). However, these studies have yielded contradictory results. The purpose of this meta-analysis is to investigate the association between the collagen IX gene polymorphisms (rs12077871, rs12722877, rs7533552 in COL9A2; rs61734651 in COL9A3) and LDD.

Methods

All relevant articles were collected from PubMed, Web of Science, and China National Knowledge Infrastructure (CNKI). The last electronic search was performed on September 1, 2017. The allele/genotype frequencies were extracted from each study. The odds ratio (OR) and 95% confidence interval (CI) were used to assess the strength of associations under the five comparison genetic models. Statistical analysis was performed by Review Manager (RevMan) 5.31 software.

Results

The meta-analysis of 10 case-control studies, including 2102 LDD cases and 2507 controls, indicated that COL9A2 gene (rs12077871, rs12722877, rs7533552) and COL9A3 gene (rs61734651) polymorphisms were not associated with LDD (rs12077871: T vs. C, OR = 1.85, 95% CI = 0.87–3.91, P = 0.11; rs12722877: G vs. C, OR = 0.83, 95% CI = 0.69–1.01, P = 0.06; rs7533552: G vs. A, OR = 1.11, 95% CI = 0.98–1.25, P = 0.09; rs61734651: T vs. C, OR = 1.57, 95% CI = 0.51–4.84, P = 0.43). The Egger text and the Begg funnel plot did not show any evidence of publication bias.

Conclusion

rs12077871, rs12722877, and rs7533552 variants in COL9A2 and rs61734651 variant in COL9A3 were not significantly associated with a predisposition to LDD. Large-scale and well-designed studies are needed to confirm this conclusion.

Genetic background of degenerative disc disease in the lumbar spine

This is a review paper on the topic of genetic background of degenerative disc diseases in the lumbar spine. Lumbar disc diseases (LDDs), such as lumbar disc degeneration and lumbar disc herniation, are the main cause of low back pain. There are a lot of studies that tried to identify the causes of LDDs. The causes have been categorized into environmental factors and genetic factors. Recent studies revealed that LDDs are mainly caused by genetic factors. Numerous studies have been carried out using the genetic approach for LDDs. The history of these studies is divided into three periods: (1) era of epidemiological research using familial background and twins, (2) era of genomic research using DNA polymorphisms to identify susceptible genes for LDDs, and (3) era of functional research to determine how the genes cause LDDs. This review article was undertaken to present the history of genetic approach to LDDs and to discuss the current issues and future perspectives.

Characteristics and determinants of clinical symptoms in radiographic lumbar spinal stenosis in a tertiary health care centre in sub-Saharan Africa

BACKGROUND

Lumbar spinal stenosis (LSS) refers to narrowing of the lumbar central spinal canal, lateral recess, and/or neuro-foramina. Radiographic LSS plays an important role in clinical LSS but is not solely accountable for the presence of symptoms. We sought to characterise clinical LSS and to determine factors associated with presence of symptoms of LSS in patients with radiographic LSS in a sub Saharan Africa setting.

METHODS

After prior ethical clearance, a case control study was done in a tertiary hospital in Douala-Cameroon, including 105 patients with radiographic LSS: 57 with symptoms of LSS (cases) and 58 with no symptoms (controls). Spinal stenosis was assessed using computed tomography (CT) scans. Data were analysed using SPSS version 23.

RESULTS

The mean age of our study participants was 53.4 ± 13.1 years. The mean age of onset of symptoms of LSS was 50.3 ± 11.6 years and the most common symptoms were Low back pain (100.0%), radicular symptoms (98.2%) and neurogenic claudication (98.2%). Obesity (p < 0.001) and a high waist circumference (p = 0.002) were significantly associated with presence of LSS symptoms in persons with radiographic LSS. After adjusting for body mass index, a positive family history of low back pain (p = 0.004), vertebra lesion at L2 (p = 0.034), L3 (p = 0.002), L4 (p = 0.025) and multiple (p = 0.008) levels, degenerative disc protrusion (p = 0.044), disc lesion at L3-L4 (p = 0.001), L4-L5 (p = 0.011) and multiple (p = 0.046) levels were significantly associated with presence of symptoms of LSS in persons with radiographic LSS.

CONCLUSION

Characteristics of clinical LSS have been described in this sub-Saharan Africa population. Obesity, a high waist circumference and a positive family history of low back pain are significantly associated with presence of symptoms of LSS in persons with radiographic LSS.

Lumbosacral stenosis in Labrador retriever military working dogs - an exomic exploratory study

Background

Canine lumbosacral stenosis is defined as narrowing of the caudal lumbar and/or sacral vertebral canal. A risk factor for neurologic problems in many large sized breeds, lumbosacral stenosis can also cause early retirement in Labrador retriever military working dogs. Though vital for conservative management of the condition, early detection is complicated by the ambiguous nature of clinical signs of lumbosacral stenosis in stoic and high-drive Labrador retriever military working dogs. Though clinical diagnoses of lumbosacral stenosis using CT imaging are standard, they are usually not performed unless dogs present with clinical symptoms. Understanding the underlying genomic mechanisms would be beneficial in developing early detection methods for lumbosacral stenosis, which could prevent premature retirement in working dogs. The exomes of 8 young Labrador retriever military working dogs (4 affected and 4 unaffected by lumbosacral stenosis, phenotypically selected by CT image analyses from 40 dogs with no reported clinical signs of the condition) were sequenced to identify and annotate exonic variants between dogs negative and positive for lumbosacral stenosis.

Results

Two-hundred and fifty-two variants were detected to be homozygous for the wild allele and either homozygous or heterozygous for the variant allele. Seventeen non-disruptive variants were detected that could affect protein effectiveness in 7 annotated (SCN1B, RGS9BP, ASXL3, TTR, LRRC16B, PTPRO, ZBBX) and 3 predicted genes (EEF1A1, DNAJA1, ZFX). No exonic variants were detected in any of the canine orthologues for human lumbar spinal stenosis candidate genes.

Conclusions

TTR (transthyretin) gene could be a possible candidate for lumbosacral stenosis in Labrador retrievers based on previous human studies that have reported an association between human lumbar spinal stenosis and transthyretin protein amyloidosis. Other genes identified with exonic variants in this study but with no known published association with lumbosacral stenosis and/or lumbar spinal stenosis could also be candidate genes for future canine lumbosacral stenosis studies but their roles remain currently unknown. Human lumbar spinal stenosis candidate genes also cannot be ruled out as lumbosacral stenosis candidate genes. More definitive genetic investigations of this condition are needed before any genetic test for lumbosacral stenosis in Labrador retriever can be developed.

Electronic supplementary material

The online version of this article (10.1186/s40575-017-0052-6) contains supplementary material, which is available to authorized users.

Genetic factors of cervical spondylotic myelopathy-a systemic review

BACKGROUND

Cervical spondylotic myelopathy (CSM) is a degenerative disorder of the neck. Recent studies have reported the roles of single nucleotide polymorphisms and abnormal gene expression in the etiology and development of CSM. However, a systemic review of these findings is currently unavailable.

METHODS

A systemic review of genetic factors of CSM was conducted through searching PubMed and EMbase databases. A total of 9 studies were included in this study, which included 8 genes: brain derived neurotrophic factor (BDNF), osteopontin (OPN), bone morphogenic protein (BMP) 4, collagen IX, vitamin D receptor (VDR), apolipoprotein E (ApoE), hypoxia-inducible factor α (HIF-1α), and cyclooxygenase 2 (COX-2).

RESULTS

The polymorphisms of 6 genes (OPN, BMP-4, collagen IX, VDR, HIF-1α) showed significant association with the susceptibility to or risk of CSM. The polymorphisms of 3 genes (BMP-4, ApoE4, HIF-1α) were significantly associated with the postoperative outcome. The polymorphism of BDNF, VDR, and expression of COX-2 were associated with the severity of disease.

CONCLUSION

This review demonstrates that 8 genes were associated with CSM although there is no repeated study. This review also suggests that large scale and high quality studies are needed to provide more reliable evidence for future evaluation.

Genetic analysis of the Vitamin D receptor start codon polymorphism (FokI) in cervical vertebra and lumbar spine pathologies: a meta-analysis

BACKGROUND

Vitamin D receptor (VDR) FokI polymorphism has been reported to influence the risk of spinal diseases. However, several studies suggest inconsistent results. Therefore, we performed this analysis to reveal the accurate relationship between VDR FokI polymorphism and spinal diseases.

MATERIALS AND METHODS

8 articles accord with the strict inclusion and exclusion criteria. 1116 cases and 1263 controls are entered into this analysis. The pooled odds ratios (ORs) and 95% confidence intervals (CI) are calculated to evaluate the association between VDR gene polymorphism and spinal diseases.

RESULT

The results suggest that allele F is a risk factor for spinal diseases and the difference is significant (F vs. f: OR = 1.151, 95% CI, 1.020-1.300). For the genotype analysis of VDR FokI, no statistical differences exist in the models of heterozygote comparison (Ff vs. ff), homozygote comparison (FF vs. ff) and dominant model (FF + Ff vs. ff) (p > 0.05). However, in recessive model (FF vs. Ff + ff), there is a significant association between VDR polymorphism and spinal diseases (OR = 1.209, 95% CI, 1.017-1.436). In subgroup analysis, the results show that allele F is a risk factor for spinal diseases in each estimation. In hospital-based subgroup, the significant differences exist in FF vs. ff and FF vs. Ff + FF models. In degenerative spine diseases group, the results are consistent with that of overall studies.

CONCLUSIONS

According to results of this meta-analysis, allele F is associated with the increased risk of spinal diseases. FF genotype may contribute to the susceptibility of spinal diseases. Therefore, VDR FokI polymorphism is related with spinal diseases.

Phenotyping of lumbosacral stenosis in Labrador retrievers using computed tomography

Deep phenotyping tools for characterizing preclinical morphological conditions are important for supporting genetic research studies. Objectives of this retrospective, cross-sectional, methods comparison study were to describe and compare qualitative and quantitative deep phenotypic characteristics of lumbosacral stenosis in Labrador retrievers using computed tomography (CT). Lumbosacral CT scans and medical records were retrieved from data archives at three veterinary hospitals. Using previously published qualitative CT diagnostic criteria, a board-certified veterinary radiologist assigned dogs as either lumbosacral stenosis positive or lumbosacral stenosis negative at six vertebral locations. A second observer independently measured vertebral canal area, vertebral fat area, and vertebral body area; and calculated ratios of vertebral canal area/vertebral body area and vertebral fat area/vertebral body area (fat area ratio) at all six locations. Twenty-five dogs were sampled (lumbosacral stenosis negative, 11 dogs; lumbosacral stenosis positive, 14 dogs). Of the six locations, cranial L6 was the most affected by lumbosacral stenosis (33%). Five of six dogs (83%) with clinical signs of lumbosacral pain were lumbosacral stenosis positive at two or more levels. All four quantitative variables were significantly smaller at the cranial aspects of the L6 and L7 vertebral foramina than at the caudal aspects (P < 0.0001). Fat area ratio was a significant predictor of lumbosacral stenosis positive status at all six locations with cranial L6 having the greatest predictive value (R² = 0.43) and range of predictive probability (25-90%). Findings from the current study supported the use of CT as a deep phenotyping tool for future research studies of lumbosacral stenosis in Labrador retrievers.

Identification of critical genes in nucleus pulposus cells isolated from degenerated intervertebral discs using bioinformatics analysis

Intervertebral disc (IVD) degeneration is a pathological process, which may lead to lower back pain. The present study aimed to investigate the pathogenesis of IVD degeneration. GSE42611 was downloaded from Gene Expression Omnibus, including 4 nucleus pulposus samples isolated from degenerated IVDs and 4 nucleus pulposus samples separated from normal IVDs. The differentially expressed genes (DEGs) between the degenerated and normal samples were screened using the limma package in R. Functional and pathway enrichment analyses were conducted separately for the upregulated and downregulated genes, using Database for Annotation, Visualization and Integrated Discovery software. In addition, protein-protein interaction (PPI) networks were constructed using the Search Tool for the Retrieval of Interacting Genes database and Cytoscape software. Finally, module analyses were conducted for the PPI networks using the MCODE plug-in in Cytoscape. A total of 558 DEGs were identified in the degenerated nucleus pulposus cells: 253 upregulated and 305 downregulated. Pathway enrichment analysis revealed that downregulated thrombospondin 1 (THBS1) was enriched in extracellular matrix-receptor interaction. Interleukin (IL)-6 in the PPI network for the upregulated genes and vascular endothelial growth factor A (VEGFA) in the PPI network for the downregulated genes had higher degrees. Additionally, four modules (µM1, µM2, µM3 and µM4) were identified from the PPI network for the upregulated genes. Four modules (dM1, dM2, dM3 and dM4) were identified from the PPI network for the downregulated genes. In the dM2 module, collagen genes and integrin subunit α4 (ITGA4) may interact with each other. Additionally, functional enrichment indicated that collagen genes were enriched in extracellular matrix organization. In conclusion, IL-6, VEGFA, THBS1, ITGA4 and collagen genes may contribute to the progression of IVD degeneration. These results suggested that the manipulation of these genes and their products may have potential as a novel therapeutic strategy for the treatment of patients with IVD.

Association between vitamin D receptor gene polymorphisms and intervertebral disc degeneration: A meta-analysis

BACKGROUND

Studies that have investigated the association between vitamin D receptor (VDR) gene polymorphisms and intervertebral disc degeneration (IDD) have yielded inconsistent results.

METHODS

To investigate the association between VDR gene polymorphisms and IDD, a systematic literature search for relevant published studies was performed on PubMed, Embase, Web of Science, Cochrane library, Wan-Fang, and CNKI databases. A random effects model was used for heterogeneous data; while a fixed effect model was used for homogenous data. Odds ratios (OR) and 95% confidence intervals (CI) were calculated to evaluate the strength of the association.

RESULTS

We observed no association between VDR FokI, TaqI-ApaI polymorphisms and IDD. However, on subgroup analysis by ethnicity, VDR FokI mutation was associated with a significantly lower risk for IDD [dominant model: OR = 0.78, 95% CI = 0.65-0.93; heterozygote model: OR = 0.76, 95% CI = 0.63-0.92; allele model: OR = 0.86, 95% CI = 0.75-0.98] among Caucasians.

CONCLUSION

These results suggest that the VDR FokI polymorphism may be associated with IDD among Caucasians. However, the association between VDR TaqI-ApaI polymorphisms and IDD in Asians is still not clear. Further well-designed studies are needed to arrive at a definitive conclusion.

PARAMETERS OF NUCLEAR MAGNETIC RESONANCE IN PATIENTS WITH CONGENITAL NARROWING OF THE LUMBAR SPINAL CANAL

ABSTRACT Objective: To compare the morphological parameters of magnetic resonance in patients with congenital narrowing of the lumbar spinal canal with patients with low back pain. Methods: A descriptive, retrospective, observational study was conducted with measurements in the axial and sagittal magnetic resonance sections of the vertebral body and canal of the lumbar spine of 64 patients with diagnosis of low back pain, which were compared with resonance images taken from 31 Mexican patients with congenital narrowing of the lumbar spinal canal. Results: The results show that patients with congenital narrowing of the lumbar spinal canal in the axial sections have a difference in diameters, being L2<13.9 mm, L3<13.3 mm, L4<12.9 mm, L5<13.1 mm, compared with controls L2<20.5 mm, L3<20.5 mm, L4<19.3 mm, L5<18.1 mm with p = 0.000. Conclusions: We found different measurements in the Mexican population compared to those found by similar studies. With the parameters obtained, it would be possible to make the proper diagnosis, surgical planning, and treatment.

Genome-Wide Meta-Analysis of Sciatica in Finnish Population

Sciatica or the sciatic syndrome is a common and often disabling low back disorder in the working-age population. It has a relatively high heritability but poorly understood molecular mechanisms. The Finnish population is a genetic isolate where small founder population and bottleneck events have led to enrichment of certain rare and low frequency variants. We performed here the first genome-wide association (GWAS) and meta-analysis of sciatica. The meta-analysis was conducted across two GWAS covering 291 Finnish sciatica cases and 3671 controls genotyped and imputed at 7.7 million autosomal variants. The most promising loci (p<1x10^-6) were replicated in 776 Finnish sciatica patients and 18,489 controls. We identified five intragenic variants, with relatively low frequencies, at two novel loci associated with sciatica at genome-wide significance. These included chr9:14344410:I (rs71321981) at 9p22.3 (NFIB gene; p = 1.30x10^-8, MAF = 0.08) and four variants at 15q21.2: rs145901849, rs80035109, rs190200374 and rs117458827 (MYO5A; p = 1.34x10^-8, MAF = 0.06; p = 2.32x10^-8, MAF = 0.07; p = 3.85x10^-8, MAF = 0.06; p = 4.78x10^-8, MAF = 0.07, respectively). The most significant association in the meta-analysis, a single base insertion rs71321981 within the regulatory region of the transcription factor NFIB, replicated in an independent Finnish population sample (p = 0.04). Despite identifying 15q21.2 as a promising locus, we were not able to replicate it. It was differentiated; the lead variants within 15q21.2 were more frequent in Finland (6–7%) than in other European populations (1–2%). Imputation accuracies of the three significantly associated variants (chr9:14344410:I, rs190200374, and rs80035109) were validated by genotyping. In summary, our results suggest a novel locus, 9p22.3 (NFIB), which may be involved in susceptibility to sciatica. In addition, another locus, 15q21.2, emerged as a promising one, but failed to replicate.

Vitamin D Receptor Gene, Matrix Metalloproteinase 3 Polymorphisms and the Risk of Intervertebral Disc Degeneration Susceptibility: Meta-Analysis

Several studies have evaluated the association between vitamin D receptor, matrix metalloproteinase 3 (MMP-3) polymorphisms and the risk of intervertebral disc degeneration susceptibility. The findings were inconsistent. This meta-analysis aimed to systematically assess the association between vitamin D receptor, MMP-3 polymorphisms and the risk of intervertebral disc degeneration susceptibility. A search of various databases was done covering all papers published until December 31th, 2014. Eight, 4, 3 studies were finally included that addressed the risk of intervertebral disc degeneration susceptibility and vitamin D receptor FokI (rs2228570), ApaI (rs7975232), and MMP-3 (rs731236) polymorphisms, respectively. FokI (f vs. F: summary odds ratio [OR], 1.13; 95% confidence interval [CI], 0.76–1.69; ff vs. FF: OR, 1.02; 95% CI, 0.59–1.77; ff vs. Ff/FF: OR, 1.05; 95% CI, 0.70–1.58), ApaI (a vs. A: OR, 0.73; 95% CI, 0.45–1.19; aa vs. AA: OR, 0.53; 95% CI, 0.22–1.25 p=0.14; aa vs. AA/Aa: OR, 0.69; 95% CI, 0.53–0.89) in the vitamin D receptor gene and MMP3 polymorphisms (5A vs. 6A: OR, 1.92; 95% CI, 0.77–4.80; 5A5A vs. 6A6A: OR, 2.17; 95% CI, 0.75–6.24; 5A5A vs. 5A6A/6A6A: OR, 1.58; 95% CI, 0.72–3.44) were not obviously associated with risk of intervertebral disc degeneration susceptibility. FokI, ApaI polymorphisms in the vitamin D receptor gene and MMP-3 polymorphism are not obvious risk factors for intervertebral disc degeneration susceptibility.

Vitamin D receptor gene polymorphisms and lumbar disc degeneration: a systematic review and meta-analysis

PURPOSE

To examine the association between Vitamin D receptor (VDR) gene polymorphisms and lumbar disc degeneration (LDD) predisposition.

METHODS

A comprehensive literature search was conducted to identify all the relevant studies. The allele/genotype frequencies were extracted from each study. We calculated the pooled odds ratios (ORs) and 95 % confidence intervals (CI) to assess the strength of the association between the VDR gene polymorphisms and LDD risk. Statistical analysis was performed using RevMan 5.31 software.

RESULTS

A total of 23 case-control studies (1835 cases and 1923 controls) were included in this systematic review. For the TaqI (rs731236), FokI (rs2228570) and ApaI (rs7975232) polymorphisms of VDR gene, nine studies, seven studies, and five studies, were eventually included in the meta-analysis, respectively. There was no evidence that the VDR gene polymorphisms (TaqI, FokI, ApaI) had significant associations with LDD risk.(for TaqI allelic comparison, OR = 1.07, 95 % CI 0.81-1.40, p = 0.64; for FokI allelic comparison, OR = 1.23, 95 % CI 0.83-1.82, p = 0.31; for ApaI allelic comparison, OR = 0.79, 95 % CI 0.55-1.14, p = 0.20). For stratified analyses by ethnicity and study design, no significant associations were found in Caucasian population and Asian population, as well as the population-based studies and hospital-based studies under all genetic models.

CONCLUSIONS

TaqI, FokI, and ApaI polymorphisms of VDR gene were not significantly associated with the predisposition of LDD. Large-scale and well-designed international studies are needed to further analyze this field.

BsmI, ApaI and TaqI Polymorphisms in the Vitamin D Receptor Gene (VDR) and Association with Lumbar Spine Pathologies: An Italian Case-Control Study

Three adjacent single nucleotide polymorphisms of the vitamin D receptor gene (VDR) BsmI (rs1544410), ApaI (rs7975232), and TaqI (rs731236) are commonly studied in several pathologies. We aimed to evaluate the distribution of VDR BsmI, ApaI, and TaqI allele, genotype, and haplotype frequencies in an Italian cohort of 266 patients with lumbar spine disorders assessed by Magnetic Resonance Imaging and 252 asymptomatic controls. The exposure to putative risk factors was evaluated by a questionnaire. Polymorphisms were detected by PCR-RFLP and TaqMan^® SNP Genotyping Assay. The results were statistically adjusted for the identified conventional risk factors. The three SNPs were in linkage disequilibrium. For all cases BbAaTT was a 3-fold risk factor OR = 3.38), whereas bbAATT (OR = 0.22), and bbaaTT (OR = 0.47) genotypes were found to be protective. Specifically, for patients affected by disc herniation only (n = 88) and all lumbar pathologies excluding stenosis and/or spondylolistesis (n = 215) B allele, Bb, Aa, and BbAaTT genotypes were risky, whereas b allele, bb, aa, and bbaaTT genotypes were protective. In patients affected by osteochondrosis with or without disc hernation (n = 50), T allele, Aa, and bbAaTT genotypes were risky, whereas t allele, AA, tt genotypes were protective. In patients affected by stenosis and/or spondylolistesis (n = 51) no significant associations were found. This is the first study showing an association of the three genetic VDR variants BsmI, ApaI, and TaqI and lumbar spine pathologies. Our study contributes to delineate genetic risk factors for specific subgroups of patients with lumbar spine pathologies highlighting the importance of haplotype analysis, and of detailed clinical evaluation of the patients for identification of genetic biomarkers.