Candidate gene investigation of spinal degenerative osteoarthritis in Greek population

Bioinformatics analysis and experimental validation of key genes associated with lumbar disc degeneration and biomechanics

Background

Lumbar disc degeneration (LDD) is an important pathological basis for the development of degenerative diseases of the lumbar spine. Most clinical patients have low back pain as their main symptom. The deterioration of the biomechanical environment is an important cause of LDD. Although there is a large amount of basic research on LDD, there are fewer reports that correlate biomechanical mechanisms with basic research. Our research aims to identify 304 key genes involved in LDD due to biomechanical deterioration, using a bioinformatics approach. We focus on SMAD3, CAV1, SMAD7, TGFB1 as hub genes, and screen for 30 potential target drugs, offering novel insights into LDD pathology and treatment options.

Methods

The Gene Cards, GenCLip3, OMIM and Drugbank databases were explored to obtain genes associated with biomechanics and LDD, followed by making veen plots to obtain both co-expressed genes. GO enrichment analysis and KEGG pathway analysis of the co-expressed genes were obtained using the DAVID online platform and visualised via a free online website. Protein interaction networks (PPI) were obtained through the STRING platform and visualised through Cytoscape 3.9.0. These genes were predicted for downstream interaction networks using the STITCH platform. Then, the GSE56081 dataset was used to validate the key genes. RT-PCR was used to detect mRNA expression of core genes in the degenerated nucleus pulposus (NP) samples and western bolt was used for protein expression. Lastly, the obtained hub genes were searched in the drug database (DGIdb) to find relevant drug candidates.

Results

From the perspective of biomechanics-induced LDD, we obtained a total of 304 genes, the GO functional enrichment and KEGG pathway enrichment analysis showed that the functions of these genes are mostly related to inflammation and apoptosis. The PPI network was constructed and four Hub genes were obtained through the plug-in of Cytoscape software, namely SMAD3, CAV1, SMAD7 and TGFB1. The analysis of key genes revealed that biomechanical involvement in LDD may be related to the TGF-β signaling pathway. Validation of the GSE56081 dataset revealed that SMAD3 and TGFB1 were highly expressed in degenerating NP samples. RT-PCR results showed that the mRNA expression of SMAD3 and TGFB1 was significantly increased in the severe degeneration group; Western blot results also showed that the protein expression of TGFB1 and P-SMAD3 was significantly increased. In addition, we identified 30 potential drugs.

Conclusion

This study presented a new approach to investigate the correlation between biomechanical mechanisms and LDD. The deterioration of the biomechanical environment may cause LDD through the TGF-β signaling pathway. TGFB1 and SMAD3 are important core targets. The important genes, pathways and drugs obtained in this study provided a new basis and direction for the study, diagnosis and treatment of LDD.

Myofascial Pain as an Unseen Comorbidity in Osteoarthritis: A Scoping Review

OBJECTIVE

This review aimed to identify, summarize, and appraise the evidence supporting the coexistence of myofascial pain (MPS) and trigger points (MTrP) in osteoarthritis (OA), and the effectiveness of MTrPs treatments in OA-related pain and physical function outcomes.

METHODS

Three databases were searched from inception to June 2022. We included observational and experimental studies to fulfill our 2 study aims. Two independent reviewers conducted 2-phase screening procedures and risk of bias using checklist tools for cross-sectional, quasi-experimental, and randomized control trials. Patient characteristics, findings of active and latent MTrPs in relevant muscles, treatments, and pain and physical function outcomes were extracted from low-risk bias studies.

RESULTS

The literature search yielded 2898 articles, of which 6 observational and 7 experimental studies had a low bias risk and the data extracted. Active MTrPs in knee OA patients was more evident in the quadriceps and hamstring muscles than in healthy individuals. Dry needling on active MTrPs improved pain and physical function in the short term compared with sham treatment in hip OA patients. In knee OA, dry needling on latent or active MTrPs improved pain and functional outcomes compared with sham needling but did not result in better pain and physical outcomes when combined with a physical exercise program.

DISCUSSION

The presence of active versus latent MTrPs seems to be a more sensitive discriminating feature of OA given that latent is often present in OA and healthy individuals. Dry needling on active MTrPs improved pain and physical function in the short term compared with sham treatment in hip OA patients. However, the small sample size and the few number of studies limit any firm recommendation on the treatment.

REGISTRY

The study protocol was prospectively registered in Open Science Framework (https://doi.org/10.17605/OSF.IO/8DVU3).

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.

Association between gene polymorphisms of TGF-β and Smad3 and susceptibility to arthritis: a meta-analysis

OBJECTIVES

This meta-analysis was performed to investigate the associations between single-nucleotide polymorphisms (SNPs) in the TGF- β and Smad3 genes and arthritis.

METHODS

A meta-analysis was performed in STATA 14.0, with publication bias and meta-regression analysis. All types of arthritis were included, and subgroup analyses were performed to interpret variations among different types of arthritis.

RESULTS

Twenty-two qualified studieswere selected to analyze the pooled accuracy, and 4 SNP sites were involved. The analysis of the TGFB1 SNP rs1800470 showed an association with arthritis in allelic (P = 0.011), homozygous (P = 0.034) and recessive (P = 0.021) genetic models. The analysis of the TGFB1 SNP rs1800471 demonstrated a close association with rheumatoid arthritis (RA) in homozygous (P = 0.000, 95%) and recessive (P = 0.008) models. The analysis of the SMAD3 SNP rs12901499 revealed a close association with osteoarthritis (OA) in the allelic (P = 0.001) model.

CONCLUSION

This research showed that genetic variants of the TGF-β pathway impact arthritis. The polymorphisms rs1800470, rs1800471 and rs12901499 were correlated with a higher prevalence of arthritis.

The association between rs12901499 polymorphism in SMAD3 gene and risk of osteoarthritis: a meta-analysis

Purpose

This study was conducted to assess and synthesize the current evidence on the association between rs12901499 polymorphism in SMAD3 gene and risk of osteoarthritis (OA).

Materials and methods

Four electronic databases, including PubMed, Embase, ISI Web of Science, and CENTRAL were systematically searched for potential studies. Summary odds ratio and corresponding 95% CI were calculated to evaluate the association. Risk of bias was assessed through the Newcastle-Ottawa Scale. Subgroups and sensitivity analyses were performed using the RevMan 5.3 software. Publication bias was evaluated by Egger's and Begg's tests. Power analysis was conducted using the Power and Sample Size Calculation program.

Results

Eight case-control studies containing 5,625 patients with OA and 5,600 healthy controls were obtained for the meta-analysis. After excluding cohorts with inadequate power, the pooled data supported that G allele carriers of rs12901499 had a significantly increased risk of OA (odds ratio 1.31, 95% CI: 1.21 to 1.43, P<0.00001). When stratified by OA site and ethnicity, the association remained statistically significant.

Conclusion

The combined results evidently supported that rs12901499 polymorphism in SMAD3 gene is significantly associated with OA vulnerability across both Caucasian and Asian populations.

Association between SMAD3 gene rs12901499 polymorphism and knee osteoarthritis in a Chinese population

BACKGROUND

Many studies have revealed that transforming growth factor-beta (TGF-β) signals play important roles in maintaining normal status of articular cartilage in human osteoarthritis (OA). However, SMAD3 had inhibitory effect on TGF-β-induced chondrocyte maturation.

METHOD

To evaluate the association of SMAD3 genetic variants with the risk of knee OA, we conducted this hospital-based case-control study involving 350 knee patients with OA and 400 controls in a Chinese population. Genotyping was performed using a custom-by-design 48-Plex single-nucleotide polymorphism (SNP) Scan^™ Kit.

RESULTS

Our results indicate that the GG genotype of rs12901499 could decrease the risk of knee OA compared to AA genotype. However, stratified analyses by sex and age did not obtain positive findings with regard to the association between rs12901499 polymorphism and knee OA risk.

CONCLUSION

In conclusion, SMAD3 rs12901499 polymorphism may be involved in the development of knee OA. Larger studies with more diverse ethnic populations are needed to confirm these results.