Sequence variant at 8q24.21 associates with sciatica caused by lumbar disc herniation

Transcriptomic imputation identifies tissue-specific genes associated with cervical myelopathy

BACKGROUND CONTEXT

Degenerative cervical myelopathy (DCM) is a progressive spinal condition that can lead to severe neurological dysfunction. Despite its degenerative pathophysiology, family history has shown to be a largely important factor in incidence and progression, suggesting that inherent genetic predisposition may play a role in pathophysiology.

PURPOSE

To determine the tissue-specific, functional genetic basis of hereditary predisposition to cervical myelopathy.

STUDY DESIGN

Retrospective case-control study using patient genetics and matched EHR from the Mount Sinai BioMe Biobank.

METHODS

In a large, diverse, urban biobank of 32,031 individuals, with 558 individuals with cervical myopathy, we applied transcriptomic imputation to identify genetically regulated gene expression signatures associated with DCM. We performed drug-repurposing analysis using the CMAP database to identify candidate therapeutic interventions to reverse the cervical myelopathy-associated gene signature.

RESULTS

We identified 16 genes significantly associated with DCM across 5 different tissues, suggesting tissue-specific manifestations of inherited genetic risk (upregulated: HES6, PI16, TMEM183A, BDH2, LINC00937, CLEC4D, USP43, SPATA1; downregulated: TTC12, CDK5, PAFAH1B2, RCSD1, KLHL29, PTPRG, RP11-620J15.3, C1RL). Drug repurposing identified 22 compounds with the potential to reverse the DCM-associated signature, suggesting points of therapeutic intervention.

CONCLUSIONS

The inherited genetic risk for cervical myelopathy is functionally associated with genes involved in tissue-specific nociceptive and proliferative processes. These signatures may be reversed by candidate therapeutics with nociceptive, calcium channel modulating, and antiproliferative effects.

CLINICAL SIGNIFICANCE

Understanding the genetic basis of DCM provides critical insights into the hereditary factors contributing to the disease, allowing for more personalized and targeted therapeutic approaches. The identification of candidate drugs through transcriptomic imputation and drug repurposing analysis offers potential new treatments that could significantly improve patient outcomes and quality of life by addressing the underlying genetic mechanisms of DCM.

Incidence of and risk factors for lumbar disc herniation with radiculopathy in adults: a systematic review

BACKGROUND

Lumbar disc herniation (LDH) with radiculopathy is associated with greater pain, disability, healthcare use, and costs compared with nonspecific low back pain. Reliable information about its incidence and risk factors were lacking.

QUESTIONS

(1) What is the incidence of lumbar disc herniation (LDH) with radiculopathy in adults? (2) What are the risk factors for LDH with radiculopathy in adults?

METHODS

Systematic review. We searched five electronic databases from 1970 to September 2023. Eligible cohort and case-control studies were identified and independently assessed for risk of bias. A qualitative best evidence synthesis of low and moderate risk of bias studies was conducted.

RESULTS

We critically reviewed 87 studies and synthesised data from 59 (68%) studies; 12 were of low and 47 of moderate risk of bias. The lower and upper bound limits of the 95% CIs of annual incidence estimates ranged from 0.3 to 2.7 per 1000 persons for surgical case definitions, from 0.04 to 1.5 per 1,000 persons for hospital-based case definitions, and from 0.1 to 298.3 per 1,000 persons for clinical case definitions. Factors associated with the development of LDH with radiculopathy included middle-age (30-50 years), smoking, higher BMI, presence of cardiovascular risk factors (in women), and greater cumulative occupational lumbar load by forward bending postures and manual materials handling, with effect sizes ranging from ranging from 1.1 (1.0-1.3) to 3.7 (2.3-6.0).

CONCLUSIONS

Incidence of LDH varies in different populations and according to case definition. Risk factors include individual, behavioural, and work-related variables. Our findings support the need to develop standardised case definitions that validly classify the clinical spectrum of LDH and for future low risk of bias studies examining causal relationships for LDH with radiculopathy in adults.

Genome-wide meta-analysis conducted in three large biobanks expands the genetic landscape of lumbar disc herniations

Given that lumbar disc herniation (LDH) is a prevalent spinal condition that causes significant individual suffering and societal costs, the genetic basis of LDH has received relatively little research. Our aim is to increase understanding of the genetic factors influencing LDH. We perform a genome-wide association analysis (GWAS) of LDH in the FinnGen project and in Estonian and UK biobanks, followed by a genome-wide meta-analysis to combine the results. In the meta-analysis, we identify 41 loci that have not been associated with LDH in prior studies on top of the 23 known risk loci. We detect LDH-associated loci in the vicinity of genes related to inflammation, disc-related structures, and synaptic transmission. Overall, our research contributes to a deeper understanding of the genetic factors behind LDH, potentially paving the way for the development of new therapeutics, prevention methods, and treatments for symptomatic LDH in the future.

Insight into the structure, function and the tumor suppression effect of gasdermin E

Gasdermin E (GSDME), which is also known as DFNA5, was first identified as a deafness-related gene that is expressed in cochlear hair cells, and mutation of this gene causes autosomal dominant neurogenic hearing loss. Later studies revealed that GSDME is mostly expressed in the kidney, placenta, muscle and brain cells, but it is expressed at low levels in tumor cells. The GSDME gene encodes the GSDME protein, which is a member of the gasdermin (GSDM) family and has been shown to participate in the induction of apoptosis and pyroptosis. The current literature suggests that Caspase-3 and Granzyme B (Gzm B) can cleave GSDME to generate the active N-terminal fragment (GSDME-NT), which integrates with the cell membrane and forms pores in this membrane to induce pyroptosis. Furthermore, GSDME also forms pores in mitochondrial membranes to release apoptosis factors, such as cytochrome c (Cyt c) and high-temperature requirement protein A2 (HtrA2/Omi), and subsequently activates the intrinsic apoptosis pathway. In recent years, GSDME has been shown to exert tumor-suppressive effects, suggesting that it has potential therapeutic effects on tumors. In this review, we introduce the structure and function of GSDME and the mechanism by which it induces cell death, and we discuss its tumor suppressive effect.

Percutaneous endoscopic lumbar discectomy combined with platelet-rich plasma injection for lumbar disc herniation: analysis of clinical and imaging outcomes

OBJECTIVE

To evaluate the clinical efficacy and imaging outcomes of percutaneous endoscopic lumbar discectomy (PELD) combined with platelet-rich plasma (PRP) for the treatment of lumbar disc herniation (LDH).

METHODS

A total of 155 patients with LDH between January 2020 and June 2022 were retrospective analyzed, of which 75 underwent PELD with PRP and 80 underwent PELD only. Clinical functional scores and imaging data were compared. Clinical functional scores included visual analog scale of leg pain (VAS-LP) and back pain (VAS-BP), Japanese Orthopedic Association score (JOA), Oswestry Disability Index (ODI) and modified MacNab criteria. Imaging data included disc height index (DHI), spinal cross-sectional area (SCSA), disc protrusion size (DPZ), and ratio value of disc grey scales (RVG).

RESULTS

Both groups showed clinical improvement, and VAS-LP, VAS-BP, JOA and ODI were significantly improved in the PRP group compared with the control group at 3, 6 and 12 months postoperatively (P < 0.05). At the last follow-up, the differences in SCSA, DPZ and RVG between the two groups were statistically significant (P < 0.05), with the PRP group being superior to the control group. The excellent and good rates of the modified Macnab criteria in the PRP group and control group were 93.3% and 90%, respectively, with no statistically significant difference (P > 0.05). No serious complications occurred during the follow-up period.

CONCLUSION

PELD combined with PRP is a safe and effective method for treating patients with LDH. PRP injection was beneficial for delaying disc degeneration and promoting disc remodeling.

Multi-Trait Exome-Wide Association Study of Back Pain-Related Phenotypes

Back pain (BP) is a major contributor to disability worldwide, with heritability estimated at 40-60%. However, less than half of the heritability is explained by common genetic variants identified by genome-wide association studies. More powerful methods and rare and ultra-rare variant analysis may offer additional insight. This study utilized exome sequencing data from the UK Biobank to perform a multi-trait gene-based association analysis of three BP-related phenotypes: chronic back pain, dorsalgia, and intervertebral disc disorder. We identified the SLC13A1 gene as a contributor to chronic back pain via loss-of-function (LoF) and missense variants. This gene has been previously detected in two studies. A multi-trait approach uncovered the novel FSCN3 gene and its impact on back pain through LoF variants. This gene deserves attention because it is only the second gene shown to have an effect on back pain due to LoF variants and represents a promising drug target for back pain therapy.

A systematic review of genome-wide association studies for pain, nociception, neuropathy, and pain treatment responses

ABSTRACT

Pain is the leading cause of disability worldwide, imposing an enormous burden on personal health and society. Pain is a multifactorial and multidimensional problem. Currently, there is (some) evidence that genetic factors could partially explain individual susceptibility to pain and interpersonal differences in pain treatment response. To better understand the underlying genetic mechanisms of pain, we systematically reviewed and summarized genome-wide association studies (GWASes) investigating the associations between genetic variants and pain/pain-related phenotypes in humans. We reviewed 57 full-text articles and identified 30 loci reported in more than 1 study. To check whether genes described in this review are associated with (other) pain phenotypes, we searched 2 pain genetic databases, Human Pain Genetics Database and Mouse Pain Genetics Database. Six GWAS-identified genes/loci were also reported in those databases, mainly involved in neurological functions and inflammation. These findings demonstrate an important contribution of genetic factors to the risk of pain and pain-related phenotypes. However, replication studies with consistent phenotype definitions and sufficient statistical power are required to validate these pain-associated genes further. Our review also highlights the need for bioinformatic tools to elucidate the function of identified genes/loci. We believe that a better understanding of the genetic background of pain will shed light on the underlying biological mechanisms of pain and benefit patients by improving the clinical management of pain.

Genetics of Intervertebral Disc Degeneration

PURPOSE OF REVIEW

Intervertebral disc degeneration is a contributor to chronic back pain. While a part of the natural aging process, early or rapid intervertebral disc degeneration is highly heritable. In this review, we summarize recent progress towards unraveling the genetics associated with this degenerative process.

RECENT FINDINGS

Use of large cohorts of patient data to conduct genome-wide association studies (GWAS) for intervertebral disc disease, and to lesser extent for aspects of this process, such as disc height, has resulted in a large increase in our understanding of the genetic etiology. Genetic correlation suggests that intervertebral disc disease is pleiotropic with risk factors for other diseases such as osteoporosis. The use of Mendelian Randomization is slowly establishing what are the causal relationships between intervertebral disc disease and factors previously correlated with this disease. The results from these human genetic studies highlight the complex nature of this disease and have the potential to lead to improved clinical management of intervertebral disc disease. Much additional work should now be focused on characterizing the causative relationship various co-morbid conditions have with intervertebral disc degeneration and on finding interventions to slow or halt this disease.

Macrophages and Intervertebral Disc Degeneration

The intervertebral disc (IVD) aids in motion and acts to absorb energy transmitted to the spine. With little inherent regenerative capacity, degeneration of the intervertebral disc results in intervertebral disc disease, which contributes to low back pain and significant disability in many individuals. Increasing evidence suggests that IVD degeneration is a disease of the whole joint that is associated with significant inflammation. Moreover, studies show elevated macrophage accumulation within the IVD with increasing levels of disease severity; however, we still need to understand the roles, be they causative or consequential, of macrophages during the degenerative process. In this narrative review, we discuss hallmarks of IVD degeneration, showcase evidence of macrophage involvement during disc degeneration, and explore burgeoning research aimed at understanding the molecular pathways regulating macrophage functions during intervertebral disc degeneration.

LINC00511/hsa-miR-573 axis-mediated high expression of Gasdermin C associates with dismal prognosis and tumor immune infiltration of breast cancer

Breast cancer (BC) is considered the second commonest human carcinoma and the most incident and mortal in the female population. Despite promising treatments for breast cancer, mortality rates of metastatic disease remain high. Gasdermin C (GSDMC) is an affiliate of the gasdermin (GSDM) family, which is involved in the process of pyroptosis. Pyroptosis is implicated in tumorigenesis, but the role of GSDMC in cancer cells is yet to be fully elucidated. In this study, we investigated the role and mechanism of GSDMC in breast cancer. We conducted a pan-cancer analysis of the expression and prognosis of GSDMC utilizing multidimensional data from The Cancer Genome Atlas (TCGA). We investigated GSDMC expression levels in 15 BC tissues and matched adjacent normal tissues by immunohistochemistry (IHC). Further verification was performed in the Gene Expression Omnibus (GEO) database. We discovered that elevated GSDMC expression was considerably linked to a worse prognosis in breast invasive carcinoma (BRCA). Next, we identified noncoding RNAs (ncRNAs) which contributing to higher expression of GSDMC by a series of expression, survival, and correlation analysis. We finally identified LINC00511/hsa-miR-573 axis to be the most promising ncRNA-associated pathways that account for GSDMC in BRCA. Furthermore, we demonstrated the significant correlations between GSDMC expression and immune infiltrates, immune checkpoints, and immune markers in BRCA. This study illustrated that ncRNAs-mediated upregulation of GSDMC linked to dismal prognosis and also exhibited a correlation with tumor immune cell infiltration in BRCA. It is anticipated to offer novel ideas for the link between pyroptosis and tumor immunotherapy.

CircRNA expression profile and potential role of hsa_circ_0040039 in intervertebral disc degeneration

PURPOSE

Circular RNAs (circRNAs) play an critical role in the pathological processes associated with IDD. However, the potential roles of circRNAs in IDD remain largely unclear. Here, we identify the circRNAs expression profiles and elucidate the potential role of candidate circRNAs in the pathogenesis of intervertebral disc degeneration (IDD) through microarray data and bioinformatics analyses.

METHODS

We obtained the datasets of microarrays (GSE67566 and GSE116726) from the Gene Expression Omnibus database. The differentially expressed circRNAs and miRNAs were identified using the Limma R package. The target miRNAs and target genes of the candidate circRNAs were predicted using an online tool. Functional enrichment analyses of the target genes were performed using the clusterProfiler R package. A protein-protein interaction (PPI) network was constructed using STRING.

RESULTS

A total of 104 differentially expressed circRNAs were identified between the IDD and the control groups, including 41 upregulated circRNAs and 63 downregulated circRNAs (cutoff criteria (|log2 fold change| > 2, P < .05)). Hsa_circ_0040039, which was the most upregulated circRNA (log2 fold change = 2.95), was selected for further analysis. The regulatory circRNA-miRNA-mRNA network comprised hsa_circ_0040039, 2 target miRNAs (hsa-miR-424-5p and hsa-miR-15b-5p), and 77 target genes. Functional enrichment analysis showed that the 77 promising target genes are mainly enriched in the ubiquitin proteasome system and Wnt signaling pathway. Further, the PPI network showed that the top 3 hub genes are BRTC, SIAH1, and UBE2V1.

CONCLUSIONS

A total of 104 differentially expressed circRNAs were identified between the IDD and control groups. Hsa_circ_0040039 may serve as a sponge of hsa-miR-424-5p and hsa-miR-15b-5p, to regulate the expression of downstream genes (such as BRTC, SIAH1, and UBE2V1); thus, it may be involved in IDD-associated pathological processes via the Wnt/β-catenin signaling pathway. Further studies are required to confirm the potential roles of hsa_circ_0040039 in IDD.

Association between MIR31HG polymorphisms and the risk of Lumbar disc herniation in Chinese Han population

Lumbar disc herniation (LDH) is a common spinal disease that endangers human health. Genetic factors play a vital role in the progression of LDH. This study aimed to explore the relationship of the MIR31HG polymorphism with LDH risk in the Chinese population. Seven candidate SNPs on MIR31HG in 504 patients with LDH and 503 healthy people were genotyped by Agena MassARRAY platform. Logistic regression was used to calculate the relationship between MIR31HG polymorphism and LDH risk under different genetic models. Multi-factor dimensionality reduction (MDR) analysis was performed to evaluate the SNP-SNP interaction. We found that rs10965059 was significantly associated with a decreased risk of LDH under the dominant (OR = 0.46, 95% CI: 0.34-0.62, P < 0.001), log-additive (OR = 0.59, 95%CI: 0.45-0.76, P < 0.001), and codominant (OR = 0.40, 95%CI: 0.29-0.55, P < 0.001) models in the overall analysis. In the subgroup analyses of age, male, and complications, we found that rs10965059 was associated with a reduced risk of LDH. However, there was no significant correlation between MiR-31HG polymorphisms and risk of LDH in females. In addition, the three SNPs (rs72703442-rs2025327-rs55683539) was mapped to a 26kb LD block with D' >0.96, suggesting a significant linkage disequilibrium presence among each pair SNPs. MDR analysis showed that the best single-locus and multi-locus models for the prediction of LDH risk were rs10965059 and seven-locus models, respectively, and both of them increased LDH risk. Our results shown that in the Chinese Han population, the MIR31HG polymorphism rs10965059 was involved in a risk to symptomatic LDH, which provides a scientific basis for early screening, prevention, diagnosis and treatment of local LDH high-risk populations.

Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology

Back pain is a common and debilitating disorder with largely unknown underlying biology. Here we report a genome-wide association study of back pain using diagnoses assigned in clinical practice; dorsalgia (119,100 cases, 909,847 controls) and intervertebral disc disorder (IDD) (58,854 cases, 922,958 controls). We identify 41 variants at 33 loci. The most significant association (ORIDD = 0.92, P = 1.6 × 10-39; ORdorsalgia = 0.92, P = 7.2 × 10-15) is with a 3'UTR variant (rs1871452-T) in CHST3, encoding a sulfotransferase enzyme expressed in intervertebral discs. The largest effects on IDD are conferred by rare (MAF = 0.07 - 0.32%) loss-of-function (LoF) variants in SLC13A1, encoding a sodium-sulfate co-transporter (LoF burden OR = 1.44, P = 3.1 × 10-11); variants that also associate with reduced serum sulfate. Genes implicated by this study are involved in cartilage and bone biology, as well as neurological and inflammatory processes.

Low Back Pain With Persistent Radiculopathy; the Clinical Role of Genetic Variants in the Genes SOX5, CCDC26/GSDMC and DCC

In a recently published genome-wide association study (GWAS) chronic back pain was associated with three loci; SOX5, CCDC26/GSDMC and DCC. This GWAS was based on a heterogeneous sample of back pain disorders, and it is unknown whether these loci are of clinical relevance for low back pain (LBP) with persistent radiculopathy. Thus, we examine if LBP with radiculopathy 12 months after an acute episode of LBP with radiculopathy is associated with the selected single nucleotide polymorphisms (SNPs); SOX5 rs34616559, CCDC26/GSDMC rs7833174 and DCC rs4384683. In this prospective cohort study, subjects admitted to a secondary health care institution due to an acute episode of LBP with radiculopathy, reported back pain, leg pain, and Oswestry Disability Index (ODI), were genotyped and followed up at 12 months (n = 338). Kruskal-Wallis H test showed no association between the SNPs and back pain, leg pain or ODI. In conclusion, LBP with radiculopathy 12 months after an acute episode of LBP with radiculopathy, is not associated with the selected SNPs; SOX5 rs34616559, CCDC26/GSDMC rs7833174 and DCC rs4384683. This absent or weak association suggests that the SNPs previously associated with chronic back pain are not useful as prognostic biomarkers for LBP with persistent radiculopathy.

3D assessment of intervertebral disc degeneration in zebrafish identifies changes in bone density that prime disc disease

Back pain is a common condition with a high social impact and represents a global health burden. Intervertebral disc disease (IVDD) is one of the major causes of back pain; no therapeutics are currently available to reverse this disease. The impact of bone mineral density (BMD) on IVDD has been controversial, with some studies suggesting osteoporosis as causative for IVDD and others suggesting it as protective for IVDD. Functional studies to evaluate the influence of genetic components of BMD in IVDD could highlight opportunities for drug development and repurposing. By taking a holistic 3D approach, we established an aging zebrafish model for spontaneous IVDD. Increased BMD in aging, detected by automated computational analysis, is caused by bone deformities at the endplates. However, aged zebrafish spines showed changes in bone morphology, microstructure, mineral heterogeneity, and increased fragility that resembled osteoporosis. Elements of the discs recapitulated IVDD symptoms found in humans: the intervertebral ligament (equivalent to the annulus fibrosus) showed disorganized collagen fibers and herniation, while the disc center (nucleus pulposus equivalent) showed dehydration and cellular abnormalities. We manipulated BMD in young zebrafish by mutating sp7 and cathepsin K, leading to low and high BMD, respectively. Remarkably, we detected IVDD in both groups, demonstrating that low BMD does not protect against IVDD, and we found a strong correlation between high BMD and IVDD. Deep learning was applied to high-resolution synchrotron µCT image data to analyze osteocyte 3D lacunar distribution and morphology, revealing a role of sp7 in controlling the osteocyte lacunar 3D profile. Our findings suggest potential avenues through which bone quality can be targeted to identify beneficial therapeutics for IVDD.

Genome-wide association studies of low back pain and lumbar spinal disorders using electronic health record data identify a locus associated with lumbar spinal stenosis

ABSTRACT

Identifying genetic risk factors for lumbar spine disorders may lead to knowledge regarding underlying mechanisms and the development of new treatments. We conducted a genome-wide association study involving 100,811 participants with genotypes and longitudinal electronic health record data from the Electronic Medical Records and Genomics Network and Geisinger Health. Cases and controls were defined using validated algorithms and clinical diagnostic codes. Electronic health record-defined phenotypes included low back pain requiring healthcare utilization (LBP-HC), lumbosacral radicular syndrome (LSRS), and lumbar spinal stenosis (LSS). Genome-wide association study used logistic regression with additive genetic effects adjusting for age, sex, site-specific factors, and ancestry (principal components). A fixed-effect inverse-variance weighted meta-analysis was conducted. Genetic variants of genome-wide significance (P < 5 × 10-8) were carried forward for replication in an independent sample from UK Biobank. Phenotype prevalence was 48.8% for LBP-HC, 19.8% for LSRS, and 7.9% for LSS. No variants were significantly associated with LBP-HC. One locus was associated with LSRS (lead variant rs146153280:C>G, odds ratio [OR] = 1.17 for G, P = 2.1 × 10-9), but was not replicated. Another locus on chromosome 2 spanning GFPT1, NFU1, and AAK1 was associated with LSS (lead variant rs13427243:G>A, OR = 1.10 for A, P = 4.3 × 10-8) and replicated in UK Biobank (OR = 1.11, P = 5.4 × 10-5). This was the first genome-wide association study meta-analysis of lumbar spinal disorders using electronic health record data. We identified 2 novel associations with LSRS and LSS; the latter was replicated in an independent sample.

Pyroptosis, metabolism, and tumor immune microenvironment

In response to a wide range of stimulations, host cells activate pyroptosis, a kind of inflammatory cell death which is provoked by the cytosolic sensing of danger signals and pathogen infection. In manipulating the cleavage of gasdermins (GSDMs), researchers have found that GSDM proteins serve as the real executors and the deterministic players in fate decisions of pyroptotic cells. Whether inflammatory characteristics induced by pyroptosis could cause damage the host or improve immune activity is largely dependent on the context, timing, and response degree. Here, we systematically review current points involved in regulatory mechanisms and the multidimensional roles of pyroptosis in several metabolic diseases and the tumor microenvironment. Targeting pyroptosis may reveal potential therapeutic avenues.

Sources of Individual Differences in Pain

Pain is an immense clinical and societal challenge, and the key to understanding and treating it is variability. Robust interindividual differences are consistently observed in pain sensitivity, susceptibility to developing painful disorders, and response to analgesic manipulations. This review examines the causes of this variability, including both organismic and environmental sources. Chronic pain development is a textbook example of a gene-environment interaction, requiring both chance initiating events (e.g., trauma, infection) and more immutable risk factors. The focus is on genetic factors, since twin studies have determined that a plurality of the variance likely derives from inherited genetic variants, but sex, age, ethnicity, personality variables, and environmental factors are also considered.

Genome-wide Association Analysis Across 16,956 Patients Identifies a Novel Genetic Association Between BMP6, NIPAL1, CNGA1 and Spondylosis

STUDY DESIGN

A case-control genome-wide association study (GWAS) on spondylosis.

OBJECTIVE

Leveraging Geisinger's MyCode initiative's multimodal dataset, we aimed to identify genetic associations with degenerative spine disease.

SUMMARY OF BACKGROUND DATA

Degenerative spine conditions are a leading cause of global disability; however, the genetic underpinnings of these conditions remain under-investigated. Previous studies using candidate-gene approach suggest a genetic risk for degenerative spine conditions, but large-scale GWASs are lacking.

METHODS

We identified 4434 patients with a diagnosis of spondylosis using ICD diagnosis codes with genotype data available. We identified a population-based control of 12,522 patients who did not have any diagnosis for osteoarthritis. A linear-mix, additive genetic model was employed to perform the genetic association tests adjusting for age, sex, and genetic principal components to account for the population structure and relatedness. Gene-based association tests were performed and heritability and genetic correlations with other traits were investigated.

RESULTS

We identified a genome-wide significant locus at rs12190551 (odds ratio = 1.034, 95% confidence interval 1.022-1.046, P = 8.5 × 10-9, minor allele frequency = 36.9%) located in the intron of BMP6. Additionally, NIPAL1 and CNGA1 achieved Bonferroni significance in the gene-based association tests. The estimated heritability was 7.19%. Furthermore, significant genetic correlations with pain, depression, lumbar spine bone mineral density, and osteoarthritis were identified.

CONCLUSION

We demonstrated the use of a massive database of genotypes combined with electronic health record data to identify a novel and significant association spondylosis. We also identified significant genetic correlations with pain, depression, bone mineral density, and osteoarthritis, suggesting shared genetic etiology and molecular pathways with these phenotypes.Level of Evidence: N/A.

Pharmacological inhibition of mTORC1 activity protects against inflammation-induced apoptosis of nucleus pulposus cells

Lumbar disc herniation is a common disease characterized by the degeneration of intervertebral discs (IVDs), accompanied by imbalance of metabolic and inflammatory homeostasis. Current studies establish that IVD degeneration is induced by increased apoptosis of nucleus pulposus (NP) cells. However, the underlying mechanisms of NP cell survival/apoptosis are not well elucidated. Here, we reveal a novel mechanism by which mTORC1 signaling controls NP cell survival through regulating metabolic homeostasis. We demonstrated that hyperactivated mTORC1 activity induced by inflammatory cytokines engenders the apoptosis of NP cells, whereas pharmacological inhibition of mTORC1 activity promotes NP cell survival. Using an integrative approach spanning metabolomics and biochemical approaches, we showed that mTORC1 activation enhanced glucose metabolism and lactic acid production, and therefore caused NP cell apoptosis. Our study identified mTORC1 in NP cells as a novel target for IVD degeneration, and provided potential strategies for clinical intervention of lumbar disc herniation.

Two GWAS-identified variants are associated with lumbar spinal stenosis and Gasdermin-C expression in Chinese population

The aim of this study is to investigate the expression levels of genome-wide association studies (GWAS)-identified variants near Gasdermin-C (GSDMC) and its association with lumbar disc degeneration (LDD) in a Chinese population. In accordance with previously reported findings, our study involved the top 4 variants; rs6651255, rs7833174, rs4130415, and rs7816342. A total of 800 participants, 400 LDD patients and 400 controls were involved in the study. The LDD patients were divided into two mutually exclusive subgroups: subgroup 1: lumbar disc herniation; subgroup 2: lumbar spinal stenosis. Genotyping were performed using TaqMan assay, and Enzyme-Linked Immunosorbent Assay (ELISA) used to measure the plasma GSDMC levels, while quantitative reverse-transcription (qRT)-PCR and immunohistochemistry (IHC) were used to evaluate the GSDMC expression levels. Among the studied variants, there were no statistically significant differences in allelic and genotypic frequencies between LDD patients and their controls (all P > 0.05). However, the subgroup analysis revealed a significant association between rs6651255 and rs7833174 in patients with lumbar spinal stenosis (subgroup 2). Furthermore, the max-statistic test revealed that the inheritance models of two variants of lumbar spinal stenosis were represented by the recessive model. The plasma and mRNA expression levels of GSDMC were significantly higher in patients with lumbar spinal stenosis compared with the control group (P < 0.05). Furthermore, the CC genotypes of rs6651255 and rs7833174 were significantly associated with increased plasma expression levels of GSDMC in patients with lumbar spinal stenosis (P < 0.01). Two GWAS-identified variants (rs6651255 and rs7833174) near GSDMC were associated with a predisposition to lumbar spinal stenosis. GSDMC protein and mRNA expression levels may have prognostic qualities as biomarkers for the existence, occurrence or development of lumbar spinal stenosis.

Can Implementation of Genetics and Pharmacogenomics Improve Treatment of Chronic Low Back Pain?

Etiology of back pain is multifactorial and not completely understood, and for the majority of people who suffer from chronic low back pain (cLBP), the precise cause cannot be determined. We know that back pain is somewhat heritable, chronic pain more so than acute. The aim of this review is to compile the genes identified by numerous genetic association studies of chronic pain conditions, focusing on cLBP specifically. Higher-order neurologic processes involved in pain maintenance and generation may explain genetic contributions and functional predisposition to formation of cLBP that does not involve spine pathology. Several genes have been identified in genetic association studies of cLBP and roughly, these genes could be grouped into several categories, coding for: receptors, enzymes, cytokines and related molecules, and transcription factors. Treatment of cLBP should be multimodal. In this review, we discuss how an individual's genotype could affect their response to therapy, as well as how genetic polymorphisms in CYP450 and other enzymes are crucial for affecting the metabolic profile of drugs used for the treatment of cLBP. Implementation of gene-focused pharmacotherapy has the potential to deliver select, more efficacious drugs and avoid unnecessary, polypharmacy-related adverse events in many painful conditions, including cLBP.

Sequence variation at 8q24.21 and risk of back pain

Back pain (BP) is a common condition of major social importance and poorly understood pathogenesis. Intervertebral lumbar disc degeneration in all its guises is one of the major biological risk factors for BP. Previously, we identified the locus at 8q24.21 associated with chronic BP, which has been found elsewhere associated with sciatica after surgery for lumbar disc herniation. In the current study we used co-localisation methods to identify the gene most likely to harbor the causal variant. We show that the same functional variant at the 8q24.21 locus is responsible for both lumbar disc degeneration and BP, and we also studied the effects of this locus on related phenotypes. Our results link the locus to intervertebral disc and bone mineral density, but not to anthropometric measurements, thus corroborating the epidemiological evidence. Moreover, the same functional variant at the locus is more likely to affect the expression of the nearby FAM49B gene, rather than the GSDMC gene, which was previously proposed as a causative one for BP.

Insight into the genetic architecture of back pain and its risk factors from a study of 509,000 individuals

Back pain (BP) is a common condition of major social importance and poorly understood pathogenesis. Combining data from the UK Biobank and CHARGE consortium cohorts allowed us to perform a very large genome-wide association study (total N = 509,070) and examine the genetic correlation and pleiotropy between BP and its clinical and psychosocial risk factors. We identified and replicated 3 BP-associated loci, including one novel region implicating SPOCK2/CHST3 genes. We provide evidence for pleiotropic effects of genetic factors underlying BP, height, and intervertebral disk problems. We also identified independent genetic correlations between BP and depression symptoms, neuroticism, sleep disturbance, overweight, and smoking. A significant enrichment for genes involved in the central nervous system and skeletal tissue development was observed. The study of pleiotropy and genetic correlations, supported by the pathway analysis, suggests at least 2 strong molecular axes of BP genesis, one related to structural/anatomical factors such as intervertebral disk problems and anthropometrics, and another related to the psychological component of pain perception and pain processing. These findings corroborate with the current biopsychosocial model as a paradigm for BP. Overall, the results demonstrate BP to have an extremely complex genetic architecture that overlaps with the genetic predisposition to its biopsychosocial risk factors. The work sheds light on pathways of relevance in the prevention and management of low BP.

Association of GSDMC polymorphisms with lumbar disc herniation among Chinese Han population

Lumbar disc herniation (LDH) is a relatively common spinal disease, but its pathogenesis is still unknown. Numerous studies have shown that LDH is closely correlated with inflammation, and it has been found to be related to some single nucleotide polymorphisms (SNPs). Our purpose is to explore the correlation between gene polymorphisms of GSDMC and LDH risk, which is of great significance for the study of the pathogenesis of LDH. DNA was extracted from 508 LDH patients and 508 controls. We select SNPs with minor allele frequency >5% in GSDMC gene from 1,000 genome project (http://www.internationalgenome.org/). Then, genotyping was performed using Agena MassARRAY. We used unconditional logistic regression analysis to calculate odds ratios (ORs) and 95% confidence intervals (CIs). The haplotype construction and analysis in GSDMC were applied to detect the association. We identified that rs77681114 in the GSDMC gene was significantly associated with a decreased risk of LDH in the alleles model (OR = 0.81, 95% CI = 0.66-0.99, p = .049) and the log-additive model (OR = 0.81, 95% CI = 0.65-0.99, p = .049) adjusted by age and gender. The haplotype "AG" constructed by rs77681114 and rs4285452 (OR = 1.24, 95% CI = 1.01-1.53, p = .039) was associated with increased risk of LDH. After age and gender stratification, rs77681114 protected LDH risk at age 49 or older in allelic model (p = .010), co-dominant model (p = .006), dominant model (p = .029), recessive model (p = .011) and log-additive model (p = .005). Rs77681114 had protective effect on female LDH risk in both co-dominant models (p = .033) and recessive models (p = .043). These studies indicated that genetic polymorphisms of GSDMC can relatively reduce the risk of LDH.

Multiple variants in collagen genes are associated with the susceptibility to lumbar disc herniation in the Chinese population

PURPOSE

Lumbar disc herniation (LDH) is an important cause of back pain and sciatica, but its aetiology is not fully understood. Single-nucleotide polymorphisms (SNPs) in specific collagen genes are known to increase the risk of lumbar disc degeneration. We performed a case-control study among the Chinese Han population to investigate whether genetic variations in collagen genes were associated with the risk of LDH or not.

METHODS

We genotyped SNPs selected from 1000 Genome Projects using Agena MassARRAY technology. Three hundred and eighty-four LDH cases were compared with 384 controls of similar age, using the odds ratio and 95% confidence interval to calculate the susceptibility in several genetic models.

RESULTS

Our results revealed that subjects with the rs6122316-C variant of the COL9A3 gene had a higher likelihood of LDH than subjects with the allele T variant in both the codominant and recessive models. In addition, after gender stratification analysis, we found significant associations between rs16970089 and rs740024 and LDH risk in females. Age stratification analysis illustrated that rs16970089 and rs6122316 were also correlated with LDH risk in people over 50 years. The smoking stratification illustrated that rs2071358 and rs740024 had an increased association with LDH risk in smokers. And after drinking stratification, we also observed the significance between rs740024 and LDH risk.

CONCLUSIONS

Variants in genes for COL1A1, COL9A3 and COL2A1 significantly influence the risk of LDH. Large and well-designed studies are needed to confirm and explain these conclusions. These slides can be retrieved under Electronic Supplementary Material.

Investigation of the STOX1 polymorphism on lumbar disc herniation

Background

Lumbar disc herniation (LDH) is a common musculoskeletal disorder affliction and associated with several genes polymorphism. Storkhead box 1 (STOX1) gene is a transcriptional factor related with several signaling pathways including inflammatory pathway. However, little is known about single‐nucleotide polymorphisms (SNPs) of STOX1 associated with LDH risk.

Methods

We conducted a case–control study among 508 LDH cases and well‐matched 508 controls, and six candidate SNPs in STOX1 were genotyped by Agena MassARRAY. Chi‐squared test, genetic model, and haploview analysis were used to evaluate associations. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated by unconditional logistic regression.

Results

In the allelic model analysis, we found the minor allele “T” of rs7903209 and “A” of rs4472827 were associated with an increased risk of LDH (p = .029, p = .016). Furthermore, in the genotype model analysis, rs7903209 polymorphism was associated with the increased susceptibility of LDH based on dominant (p = .033) and additive model (p = .024); and rs4472827 variant was found to play a harmful role in the LDH risk based on genotype (p = .014), dominant (p = .012), and additive model (p = .015). In the haplotype analysis, the haplotype “GT” in block (rs10998461 and rs10998468) decreased LDH risk (OR = 0.7, 95% CI = 0.52–0.93, p = .016). Functional assessment indicated that rs7903209 and rs4472827 polymorphisms may influence the expression of STOX1.

Conclusion

Our results provide evidence for polymorphisms of rs7903209 and rs4472827 in STOX1 associated with LDH risk in Chinese Han population.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

Intervertebral Disc Biology: Genetic Basis of Disc Degeneration

PURPOSE OF REVIEW

This review aims to highlight recent advances in understanding the genetic basis of intervertebral disc degeneration (IDD).

RECENT FINDINGS

It has been known for some time that IDD is highly heritable. Recent studies, and in particular the availability of agnostic techniques such as genome-wide association studies, have identified new variants in a variety of genes which contribute to the risk of IDD and to back pain.

SUMMARY

A variety of genetic variants are involved in IDD. Some are shared with variants predisposing to back pain, but few have been identified reliably in either phenotype. Further research is required to explain fully the high heritability and how the genetic variants influence cell biology to lead to IDD.

Trans-Ethnic Polygenic Analysis Supports Genetic Overlaps of Lumbar Disc Degeneration With Height, Body Mass Index, and Bone Mineral Density

Lumbar disc degeneration (LDD) is age-related break-down in the fibrocartilaginous joints between lumbar vertebrae. It is a major cause of low back pain and is conventionally assessed by magnetic resonance imaging (MRI). Like most other complex traits, LDD is likely polygenic and influenced by both genetic and environmental factors. However, genome-wide association studies (GWASs) of LDD have uncovered few susceptibility loci due to the limited sample size. Previous epidemiology studies of LDD also reported multiple heritable risk factors, including height, body mass index (BMI), bone mineral density (BMD), lipid levels, etc. Genetics can help elucidate causality between traits and suggest loci with pleiotropic effects. One such approach is polygenic score (PGS) which summarizes the effect of multiple variants by the summation of alleles weighted by estimated effects from GWAS. To investigate genetic overlaps of LDD and related heritable risk factors, we calculated the PGS of height, BMI, BMD and lipid levels in a Chinese population-based cohort with spine MRI examination and a Japanese case-control cohort of lumbar disc herniation (LDH) requiring surgery. Because most large-scale GWASs were done in European populations, PGS of corresponding traits were created using weights from European GWASs. We calibrated their prediction performance in independent Chinese samples, then tested associations with MRI-derived LDD scores and LDH affection status. The PGS of height, BMI, BMD and lipid levels were strongly associated with respective phenotypes in Chinese, but phenotype variances explained were lower than in Europeans which would reduce the power to detect genetic overlaps. Despite of this, the PGS of BMI and lumbar spine BMD were significantly associated with LDD scores; and the PGS of height was associated with the increased the liability of LDH. Furthermore, linkage disequilibrium score regression suggested that, osteoarthritis, another degenerative disorder that shares common features with LDD, also showed genetic correlations with height, BMI and BMD. The findings suggest a common key contribution of biomechanical stress to the pathogenesis of LDD and will direct the future search for pleiotropic genes.