Association of GDF5, SMAD3 and RUNX2 polymorphisms with temporomandibular joint osteoarthritis in female Han Chinese

Correlation of growth differentiation factor-5 + 104T>C polymorphism with the risk of knee, hand, and hip osteoarthritis: a case-control study and meta-analysis based on 47 case-control studies

Osteoarthritis (OA) arises from a intricate interplay of genetic and environmental factors. Numerous studies have explored the link between the growth differentiation factor 5 (GDF-5) +104T>C polymorphism and OA risk, but the findings have been inconclusive. We carried out a case-control study with 704 OA cases and 418 healthy controls. Furthermore, we conducted a meta-analysis by thoroughly searching the literature for relevant studies published until 1 September, 2023. The combined odds ratio and 95% confidence intervals were used to assess the correlation's strength. A total of 47 independent case-control studies, including 17,602 OA cases and 30,947 controls, were analyzed. Of these, 31 studies (11,176 cases, 16,724 controls) focused on knee OA, 8 studies (3,973 cases, 8,055 controls) examined hip OA, and 6 studies (2244 cases, 5965 controls) investigated hand OA. Overall, our findings suggest that the GDF-5 + 104T>C polymorphism has a protectibe role in development of OA in global scale. Subgroup analyses by ethnicity indicated that this genetic variation provides protection against OA in Caucasian, Asian, and African populations. Further subgroup analysis based on the type of OA showed a decreased risk of knee and hand OA associated with this variation, but not for hip OA. Our combined data indicates that the GDF-5 + 104T>C polymorphism offers protection against the development of OA in general, as well as knee and hand OA. Nevertheless, there was no correlation found between this polymorphism and the development of hip OA.

Transcriptional profiles of TGF-β superfamily members in the lumbar DRGs and the effects of activins A and C on inflammatory pain in rats

Signaling by the transforming growth factor (TGF)-β superfamily is necessary for proper neural development and is involved in pain processing under both physiological and pathological conditions. Sensory neurons that reside in the dorsal root ganglia (DRGs) initially begin to perceive noxious signaling from their innervating peripheral target tissues and further convey pain signaling to the central nervous system. However, the transcriptional profile of the TGF-β superfamily members in DRGs during chronic inflammatory pain remains elusive. We developed a custom microarray to screen for transcriptional changes in members of the TGF-β superfamily in lumbar DRGs of rats with chronic inflammatory pain and found that the transcription of the TGF-β superfamily members tends to be downregulated. Among them, signaling of the activin/inhibin and bone morphogenetic protein/growth and differentiation factor (BMP/GDF) families dramatically decreased. In addition, peripherally pre-local administration of activins A and C worsened formalin-induced acute inflammatory pain, whereas activin C, but not activin A, improved formalin-induced persistent inflammatory pain by inhibiting the activation of astrocytes. This is the first report of the TGF-β superfamily transcriptional profiles in lumbar DRGs under chronic inflammatory pain conditions, in which transcriptional changes in cytokines or pathway components were found to contribute to, or be involved in, inflammatory pain processing. Our data will provide more targets for pain research.

Transcriptional Profiling of TGF-β Superfamily Members in Lumbar DRGs of Rats Following Sciatic Nerve Axotomy and Activin C Inhibits Neuropathic Pain

BACKGROUND

Neuroinflammation and cytokines play critical roles in neuropathic pain and axon degeneration/regeneration. Cytokines of transforming growth factor-β superfamily have implications in pain and injured nerve repair processing. However, the transcriptional profiles of the transforming growth factor-β superfamily members in dorsal root ganglia under neuropathic pain and axon degeneration/regeneration conditions remain elusive.

OBJECTIVE

We aimed to plot the transcriptional profiles of transforming growth factor-β superfamily components in lumbar dorsal root ganglia of sciatic nerve-axotomized rats and to further verify the profiles by testing the analgesic effect of activin C, a representative cytokine, on neuropathic pain.

METHODS

Adult male rats were axotomized in sciatic nerves, and lumbar dorsal root ganglia were isolated for total RNA extraction or section. A custom microarray was developed and employed to plot the gene expression profiles of transforming growth factor-β superfamily components. Realtime RT-PCR was used to confirm changes in the expression of activin/inhibin family genes, and then in situ hybridization was performed to determine the cellular locations of inhibin α, activin βC, BMP-5 and GDF-9 mRNAs. The rat spared nerve injury model was performed, and a pain test was employed to determine the effect of activin C on neuropathic pain.

RESULTS

The expression of transforming growth factor-β superfamily cytokines and their signaling, including some receptors and signaling adaptors, were robustly upregulated. Activin βC subunit mRNAs were expressed in the small-diameter dorsal root ganglion neurons and upregulated after axotomy. Single intrathecal injection of activin C inhibited neuropathic pain in spared nerve injury model.

CONCLUSION

This is the first report to investigate the transcriptional profiles of members of transforming growth factor-β superfamily in axotomized dorsal root ganglia. The distinct cytokine profiles observed here might provide clues toward further study of the role of transforming growth factor-β superfamily in the pathogenesis of neuropathic pain and axon degeneration/regeneration after peripheral nerve injury.

Genetic overlap between temporomandibular disorders and primary headaches: A systematic review

Primary headache disorders (PHD), specifically migraine, are strongly associated with temporomandibular disorders (TMD), sharing some patterns of orofacial pain. Both disorders have significant genetic contributions already studied. PRISMA guidelines were followed to conduct this systematic review, which comprehensively summarize and discuss the genetic overlap between TMD and PHD to aid future research in potential therapy targets. This review included eight original articles published between 2015 and 2020, written in English and related to either TMD and/or PHD. The genes simultaneously assessed in PHD and TMD studies were COMT, MTHFR, and ESR1. COMT was proved to play a critical role in TMD pathogenesis, as all studies have concluded about its impact on the occurrence of the disease, although no association with PHD was found. No proof on the impact of MTHFR gene regulation on either TMD or PHD was found. The most robust results are concerning the ESR1 gene, which is present in the genetic profile of both clinical conditions. This novel systematic review highlights not only the need for a clear understanding of the role of ESR1 and COMT genes in pain pathogenesis, but it also evaluates their potential as a promising therapeutic target to treat both pathologies.

Effect of IL-6 -174G/C and -572G/C variants on susceptibility to osteoarthritis in Turkish population

Osteoarthritis (OA) is a complex disorder characterized by degenerative articular cartilage in which inflammatory mechanisms play a major role in the pathogenesis. Interleukin-6 (IL6), a multifunctional cytokine, can trigger osteoclast differentiation and bone resorption. Our purpose in this study was to evaluate the association of IL-6 -174 G/C (rs1800795) and -572 G/C (rs1800796) variants with the susceptibility to OA. One hundred fifty OA patients and 150 healthy individuals were enrolled in the study. Polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) was used for genotyping the IL-6 gene variants. The results of analyses were evaluated for statistical significance. The pain intensity was assessed using the Visual Analogue Scale (VAS). There was a statistically significant difference in the genotype and allele frequencies of the IL-6 -174 G/C variant between patients with OA and control groups (p = 0.001, p = 0.002, respectively). IL-6 -174 G/C GG genotype and G allele were more prevalent in patients with OA. We found that the IL-6 -572 G/C variant was not different between patients and controls in either genotype distribution and allele frequency. IL-6 174 G/C and -572 G/C loci GG-GG combined genotype was significantly higher in OA patients (p = 0.00). Our study suggests that there was a strong association between the IL-6 -174 G/C variant and OA in the Turkish population. Further studies on populations of different ethnic background are necessary to prove the association of IL-6 variants with OA.

Joint disease-specificity at the regulatory base-pair level

Given the pleiotropic nature of coding sequences and that many loci exhibit multiple disease associations, it is within non-coding sequence that disease-specificity likely exists. Here, we focus on joint disorders, finding among replicated loci, that GDF5 exhibits over twenty distinct associations, and we identify causal variants for two of its strongest associations, hip dysplasia and knee osteoarthritis. By mapping regulatory regions in joint chondrocytes, we pinpoint two variants (rs4911178; rs6060369), on the same risk haplotype, which reside in anatomical site-specific enhancers. We show that both variants have clinical relevance, impacting disease by altering morphology. By modeling each variant in humanized mice, we observe joint-specific response, correlating with GDF5 expression. Thus, we uncouple separate regulatory variants on a common risk haplotype that cause joint-specific disease. By broadening our perspective, we finally find that patterns of modularity at GDF5 are also found at over three-quarters of loci with multiple GWAS disease associations.

Frequency of Growth Differentiation Factor 5 rs143383 and asporin D-repeat polymorphisms in patients with hand and knee osteoarthritis in Kurdistan province, Iran

AIM

Osteoarthritis (OA) is the most common chronic joint disorder, resulting from the breakdown of joint cartilage. It occurs in the knees, hands, and hips, leading to pain, stiffness, inflammation, and swelling.

METHODS

In this study, 100 hand and knee OA patients, meeting the American College of Rheumatology criteria were included in the case group, and 100 healthy individuals were allocated to the control group. Blood samples were collected from the participants. After DNA extraction, genotyping was carried out for GDF5 rs143383 C/T polymorphism by allele-specific polymerase chain reaction (ASPCR) and for D-repeat alleles of asporin (ASPN) by conventional PCR assay.

RESULTS

The results showed that the frequency of the D14 allele of ASPN was significantly higher than other alleles in the case group (P = .0001). Also, the frequency of the D14 allele among women was significantly higher than in men (P = .004). Moreover, the frequency of the TT allele in GDF5 rs143383 C/T polymorphism was significantly higher than the CC and CT alleles in the case group, compared with the control group (P = .001). A significant difference was found between the TT allele and other alleles in female and male patients compared with the control group (P = .02).

CONCLUSIONS

The D14 allele of the ASPN gene and TT allele of the GDF5 gene (rs143383 + 104T/C) are associated with hand and knee OA in the Kurdish population, indicating that these alleles could be risk factors for OA, at least in our populations.

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.

Detection of insertions/deletions (InDels) within the goat Runx2 gene and their association with litter size and growth traits

Runt-related transcription factor 2 (Runx2) is characterized by its critical functions in osteoblastic and ovulatory processes. The goal of this study was to explore the insertion/deletion (indel) variants of this gene and to evaluate their association with productive traits. Herein, a 12 bp and 6 bp insertion within the Runx2 gene was uncovered in Shaanbei white cashmere goats (SBWC; n = 1200). Chi-square analysis revealed that the 12 bp insertion was related to litter size (p < 0.01). Further association analysis also found this insertion was significantly associated with litter size (p = 1.1E-5). Interestingly, this insertion was also significantly associated with chest circumference (p = 0.018). Additionally, the 6 bp insertion was associated with body length (p = 0.003), chest width (p = 0.011), and chest circumference (p = 0.005). Furthermore, diplotype associations also uncovered that the combined genotypes of these two indels also significantly affected litter size and growth traits (p < 0.05). These findings suggested that these two insertions within the Runx2 gene were significantly associated with reproduction and growth traits, which would make them beneficial functional DNA markers that can be used in goat breeding.

Genome-wide association reveals contribution of MRAS to painful temporomandibular disorder in males

Painful temporomandibular disorders (TMDs) are the leading cause of chronic orofacial pain, but its underlying molecular mechanisms remain obscure. Although many environmental factors have been associated with higher risk of developing painful TMD, family and twin studies support a heritable genetic component as well. We performed a genome-wide association study assuming an additive genetic model of TMD in a discovery cohort of 999 cases and 2031 TMD-free controls from the Orofacial Pain: Prospective Evaluation and Risk Assessment (OPPERA) study. Using logistic models adjusted for sex, age, enrollment site, and race, we identified 3 distinct loci that were significant in combined or sex-segregated analyses. A single-nucleotide polymorphism on chromosome 3 (rs13078961) was significantly associated with TMD in males only (odds ratio = 2.9, 95% confidence interval: 2.02-4.27, P = 2.2 × 10). This association was nominally replicated in a meta-analysis of 7 independent orofacial pain cohorts including 160,194 participants (odds ratio = 1.16, 95% confidence interval: 1.0-1.35, P = 2.3 × 10). Functional analysis in human dorsal root ganglia and blood indicated this variant is an expression quantitative trait locus, with the minor allele associated with decreased expression of the nearby muscle RAS oncogene homolog (MRAS) gene (beta = -0.51, P = 2.43 × 10). Male mice, but not female mice, with a null mutation of Mras displayed persistent mechanical allodynia in a model of inflammatory pain. Genetic and behavioral evidence support a novel mechanism by which genetically determined MRAS expression moderates the resiliency to chronic pain. This effect is male-specific and may contribute to the lower rates of painful TMD in men.

Association between GDF5 rs143383 genetic polymorphism and musculoskeletal degenerative diseases susceptibility: a meta-analysis

BACKGROUND

Several studies have assessed the association between GDF5 rs143383 polymorphism and the susceptibility of musculoskeletal degenerative diseases, such as intervertebral disc degeneration (IDD) and osteoarthritis (OA), but the results are inconsistent. The aim of our study was to evaluate the association between them comprehensively.

METHODS

A systematical search was conducted on PubMed, Scopus, Web of Science (WOS), Embase, and the Cochrane Library databases updated to April 20, 2018. Eligible studies about polymorphisms in GDF5 gene and risk of IDD or OA were included. Pooled odds ratios (ORs) and 95% confidence intervals (95% CIs) were utilized.

RESULTS

Fifteen studies with a total of 5915 cases and 12,252 controls were finally included in our study. Meta-analysis of GDF5 rs143383 polymorphism was statistically associated with increased risk of musculoskeletal degenerative diseases under each genetic model (allele model: OR = 1.32, 95% CI 1.19-1.48, P = 0.000; homozygote model: OR = 1.80, 95%CI 1.49-2.16, P = 0.000; heterozygote model: OR = 1.37, 95%CI 1.21-1.55, P = 0.000; dominant model: OR = 1.56, 95%CI 1.39-1.75, P = 0.000; recessive model: OR = 1.39, 95%CI 1.20-1.60, P = 0.000). Stratified analyses based on disease type showed a significant association between the GDF5 rs143383 polymorphism and increased risk of IDD and OA under all genetic models studied. When stratified with ethnicity, pooled outcomes revealed that this polymorphism was significantly related with increased risk of musculoskeletal degenerative diseases in both Asian and Caucasian populations under all genetic models studied.

CONCLUSIONS

The present study suggested that GDF5 rs143383 polymorphism was significantly associated with susceptibility to musculoskeletal degenerative diseases.

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.

SMAD3 gene rs12901499 polymorphism increased the risk of osteoarthritis

A growing body of evidence suggested that smad family member 3 gene rs12901499 polymorphism was associated with the risk of osteoarthritis. However, the results of previous studies were conflicting. In the present study, we assessed whether smad family member 3 gene rs12901499 polymorphism was associated with the risk of osteoarthritis by the meta-analysis. We searched in the databases of PubMed, Embase, and CNKI. Pooled odds ratios and 95% confidence intervals were calculated. Seven papers involving 11 studies (5344 cases and 9080 controls) analyzed the association between smad family member 3 gene rs12901499 polymorphism and osteoarthritis risk. This meta-analysis confirmed that smad family member 3 gene rs12901499 polymorphism increased the risk of osteoarthritis. Stratification analysis of ethnicity found that rs12901499 polymorphism increased the risk of osteoarthritis among both Asians and Caucasians [G vs A: Asians, OR and 95%CI, 1.34(1.07, 1.69), P=0.012; Caucasians, OR and 95%CI, 1.21(1.13, 1.29), P<0.001]. In addition, subgroup analysis by type of osteoarthritis revealed that smad family member 3 gene rs12901499 polymorphism was correlated with the increased risk of hip osteoarthritis, but not associated with knee osteoarthritis. Sensitivity analysis did not draw different findings. In conclusion, this meta-analysis indicates that smad family member 3 gene rs12901499 polymorphism increased the risk of osteoarthritis.

Combined bioinformatics analysis reveals gene expression and DNA methylation patterns in osteoarthritis

Osteoarthritis (OA) is a common type of arthritis, which may cause pain and disability. Alterations in gene expression and DNA methylation have been proven to be associated with the development of OA. The aim of the present study was to identify potential therapeutic targets and associated processes for OA via the combined analysis of gene expression and DNA methylation datasets. The gene expression and DNA methylation profiles were obtained from the Gene Expression Omnibus, and differentially expressed genes (DEGs) and differentially methylated sites (DMSs) were identified in the present study, using R programming software. The enriched functions of DEGs and DMSs were obtained via the Database for Annotation, Visualization and Integrated Discovery. Finally, cross analysis of DEGs and DMSs was performed to identify genes that exhibited differential expression and methylation simultaneously. The protein‑protein interaction (PPI) network of overlaps between DEGs and DMSs was obtained using the Human Protein Reference Database; the topological properties of PPI network overlaps were additionally obtained. Hub genes in the PPI network were further confirmed via reverse transcription‑quantitative polymerase chain reaction (RT‑qPCR). The results of the present study revealed that the majority of DEGs and DMSs were upregulated and hypomethylated in patients with OA, respectively. DEGs and DMSs were primarily involved in inflammatory, immune and gene expression regulation‑associated processes and pathways. Cross analysis revealed 30 genes that exhibited differential expression and methylation in OA simultaneously. Topological analysis of the PPI network revealed that numerous genes, including G protein subunit α1 (GNAI1), runt related transcription factor 2 (RUNX2) and integrin subunit β2 (ITGB2), may be involved in the development of OA. Additionally, RT‑qPCR analysis of GNAI1, RUNX2 and ITGB2 provided further confirmation. Numerous known and novel therapeutic targets were obtained via network analysis. The results of the present study may be beneficial for the diagnosis and treatment of OA.

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.

Intra-articular therapy with recombinant human GDF5 arrests disease progression and stimulates cartilage repair in the rat medial meniscus transection (MMT) model of osteoarthritis

OBJECTIVE

Investigation of osteoarthritis (OA) risk alleles suggests that reduced levels of growth and differentiation factor-5 (GDF5) may be a precipitating factor in OA. We hypothesized that intra-articular recombinant human GDF5 (rhGDF5) supplementation to the OA joint may alter disease progression.

METHODS

A rat medial meniscus transection (MMT) joint instability OA model was used. Animals received either one intra-articular injection, or two or three bi-weekly intra-articular injections of either 30 μg or 100 μg of rhGDF5 beginning on day 21 post surgery after structural pathology had been established. Nine weeks after MMT surgery, joints were processed for histological analysis following staining with toluidine blue. Control groups received intra-articular vehicle injections, comprising a glycine-buffered trehalose solution. OA changes in the joint were evaluated using histopathological end points that were collected by a pathologist who was blinded to treatment.

RESULTS

Intra-articular rhGDF5 supplementation reduced cartilage lesions on the medial tibial plateau in a dose-dependent manner when administered therapeutically to intercept OA disease progression. A single 100 μg rhGDF5 injection on day 21 slowed disease progression at day 63. A similar effect was achieved with two bi-weekly injections of 30 μg. Two bi-weekly injections of 100 μg or three bi-weekly injections of 30 μg stopped progression of cartilage lesions. Importantly, three biweekly injections of 100 μg rhGDF5 stimulated significant cartilage repair.

CONCLUSIONS

Intra-articular rhGDF5 supplementation can prevent and even reverse OA disease progression in the rat MMT OA model. Collectively, these results support rhGDF5 supplementation as an intra-articular disease modifying OA therapy.

Association between GDF5 +104T/C polymorphism and knee osteoarthritis in Caucasian and Asian populations: a meta-analysis based on case-control studies

Background

Osteoarthritis (OA) is a degenerative joint disease with a complex genetic background. Variants in growth differentiation factor 5 (GDF5) have been reported to be associated with rheumatoid arthritis (RA) in several ethnic populations. The present study aimed to assess the association between the GDF5 +104T/C polymorphism and the susceptibility of the knee to OA through a meta-analysis of available case-control studies.

Methods

The PubMed and Science Direct citation databases were used to search electronic literature in order to identify studies published between January 2007 and July 2016 that evaluated the association between the GDF5 +104T/C polymorphism and the susceptibility of the knee to OA. Different genetic models were used to assess the pooled and stratified data.

Results

A positive association was found in all pooled studies (OR = 0.808, 95 % CI = 0.754–0.866, p < 0.001). Regarding genotypes, significant associations were found using a dominant model (OR = 0.777, 95 % CI = 0.708–0.852, p < 0.001), a recessive model (OR = 0.723, 95%CI = 0.623–0.839, p < 0.001), and an additive model (CC vs TT OR = 0.648, 95 % CI = 0.552–0.760, p < 0.001; CC vs CT OR = 0.801, 95 % CI = 0.685–0.936, p = 0.005). Meta-analysis data were stratified by ethnicity, and the GDF5 C allele was found to be positively associated with OA of the knee in both Caucasians and Asians, as were the GDF5 TC and CC genotypes. In addition, using an additive model, the CC genotype was found to be significantly associated with OA of the knee in both Caucasians and Asians when comparing CC vs TT genotypes, but not in Caucasians when comparing TT vs CT genotypes.

Conclusions

Meta-analysis results indicated that the GDF5 +104T/C polymorphism is a protective factor for OA among Caucasian and Asian populations.