Skip to main content
Download PDF

Polymorphism rs143384 GDF5 reduces the risk of knee osteoarthritis development in obese individuals and increases the disease risk in non-obese population

Arthroplasty volume 6, Article number: 12 (2024) Cite this article

Abstract

Background

We investigated the effect of obesity on the association of genome-wide associative studies (GWAS)-significant genes with the risk of knee osteoarthritis (KOA).

Methods

All study participants (n = 1,100) were divided into 2 groups in terms of body mass index (BMI): BMI ≥ 30 (255 KOA patients and 167 controls) and BMI < 30 (245 KOA and 433 controls). The eight GWAS-significant KOA single nucleotide polymorphisms (SNP) of six candidate genes, such as LYPLAL1 (rs2820436, rs2820443), SBNO1 (rs1060105, rs56116847), WWP2 (rs34195470), NFAT5 (rs6499244), TGFA (rs3771501), GDF5 (rs143384), were genotyped. Logistic regression analysis (gPLINK online program) was used for SNPs associations study with the risk of developing KOA into 2 groups (BMI ≥ 30 and BMI < 30) separately. The functional effects of KOA risk loci were evaluated using in silico bioinformatic analysis.

Results

Multidirectional relationships of the rs143384 GDF5 with KOA in BMI-different groups were found: This SNP was KOA protective locus among individuals with BMI ≥ 30 (OR 0.41 [95%CI 0.20–0.94] recessive model) and was disorder risk locus among individuals with BMI < 30 (OR 1.32 [95%CI 1.05–1.65] allele model, OR 1.44 [95%CI 1.10–1.86] additive model, OR 1.67 [95%CI 1.10–2.52] dominant model). Polymorphism rs143384 GDF5 manifested its regulatory effects in relation to nine genes (GDF5, CPNE1, EDEM2, ERGIC3, GDF5OS, PROCR, RBM39, RPL36P4, UQCC1) in adipose tissue, which were involved in the regulation of pathways of apoptosis of striated muscle cells.

Conclusions

In summary, the effect of obesity on the association of the rs143384 GDF5 with KOA was shown: the “protective” value of this polymorphism in the BMI ≥ 30 group and the “risk” meaning in BMI < 30 cohort.

Introduction

Osteoarthritis (OA) is a whole-joint disease involving all joint tissues (cartilage, subchondral bone, synovial membrane, meniscus, and infrapatellar fat pad) [1]. Knee osteoarthritis (KOA) represents the most common joint disease with a systemic metabolic component [2]. KOA affects 16% of the population over the age of 15 worldwide, and in 2020 about 654.1 million people over the age of 40 suffered from this condition [3]. The prevalence of KOA is constantly on the rise, primarily due to the increasing average life expectancy, as well as higher rates of obesity among the population [4]. Across the globe, KOA is considered to be a significant public health problem that has serious social and economic consequences [5]. The main reasons for KOA patients to seek medical help are pain and loss of joint function [6]. Total knee replacement is currently the most used treatment option for end-stage KOA [7]. The total number of knee replacements is expected to grow to 3.48 million by 2030 [8]. The material costs associated with KOA account for about 0.5% of gross domestic product in developed countries [9].

A number of factors, such as age, female sex, obesity, genetics, joint injuries, vitamin D deficiency, etc., have been identified to be the leading risk factors for KOA [10, 11]. Among them, obesity and overweight are two major modifiable risk factors [12]. Traditionally, one of the KOA causes is believed to be the over-weight-related biomechanical load on the joint [13]. At the same time, it is known that adipokines and proinflammatory cytokines produced by systemic and local adipose tissues are involved in cartilage degradation, synovial membrane inflammation, and bone erosion [2]. In addition, adipose tissue in the knee joint (the infrapatellar and suprapatellar fat pads, and other small fat pads such as posterior knee fat pad and posterior suprapatellar fat pad [prefemoral]) can interact with neighboring tissues, thereby potentially affecting homeostasis of joint and leading to destructive processes in KOA due to pro-inflammatory mediators [2, 14, 15]. Therefore, KOA is presently considered to be a disease entity and is aggravated by a metabolic component associated with adipose tissue [13]. Obese or overweight people are three times more likely to develop KOA than individuals with normal body weight [16]. It is known that the progression of the disease is more often observed in obese/overweight KOA patients [17]. Obviously, BMI plays a substantial role in the predisposition to KOA, but the mechanisms (including genetic one) underlying this relationship remain unclear.

The KOA is a polygenic disease [18, 19]. Thanks to genome-wide association studies (GWAS), to date, more than 80 polymorphic loci associated with the development of KOA are known [20]. Several studies have examined the relationship of various genetic variants with BMI in KOA patients [21,22,23,24,25,26,27,28]. At the same time, despite numerous data indicating a significant relationship between BMI and the development/progression of KOA [12, 13, 17, 29], genetic studies revealing the role of individual GWAS-significant polymorphic loci in the disease formation in interaction with BMI are very limited [25, 28].

Therefore, in this study, we investigated the possible effect of obesity on the GWAS-significant genetic association with the risk of developing KOA.

Materials and methods

KOA patients and controls

The study was of “patient-control” design, and involved 1,100 subjects (500 patients with KOA and 600 individuals without KOA), who were divided into two groups in terms of BMI: group I, including individuals with a BMI ≥ 30 (255 KOA patients and 167 controls); and group II, consisting of subjects with a BMI < 30 (245 KOA patients and 433 controls). The anthropometric indicators (weight, height) were gathered by previously outlined standard methods [30]. BMI was computed by using the standard method (ratio between body weight (in kilograms) and height (in meters) in squared [kg/m2]) [31]. We used well-accepted BMI grading, i.e., < 18.5 (underweight), 18.5–24.9 (normal weight), 25.0–29.9 (overweight), and ≥ 30 (obese) [32]. The KOA patients for the study were selected by certified orthopedic-traumatologists over the period from February 2016 to December 2018, based on “Belgorod City Hospital No. 2” (Department of Orthopedics and Traumatology). The study was approved by the ethics committee of this hospital.

Several inclusion criteria were used in the formation of KOA and control cohorts:

  • (1) individuals of European origin, who were born and living in the Central region of Russia and were not related to each other [33,34,35]; (2) aged 40 years or older; (3) availability of informed consent to take part in the study; (4) the KOA group included patients with [36], (i) primary KOA of the knee joint, diagnosed against the American College of Rheumatology [37], (ii) KOA radiological stage by J. Kellgren-J. Lawrence (K/L) ≥ 2 [38], (iii) the presence of pain syndrome more than 40 points on the Visual Analog Scale (VAS) [39]; (5) the control group included subjects who did not have any pathology of the musculoskeletal system. Exclusion criteria were as follows: (1) the presence of severe hypertension, coronary heart disease, diabetes mellitus, renal-hepatic insufficiency, oncological diseases, systemic connective tissue diseases, joint injuries in the anamnesis, inflammatory joint diseases, congenital malformations of the musculoskeletal system, (2) refusal to participate in the study.

SNP selection criteria and genotyping

The eight GWAS-significant KOA SNPs of six candidate genes, such as LYPLAL1—lysophospholipase like 1 (rs2820436, rs2820443), SBNO1—strawberry notch homolog 1 (rs1060105, rs56116847), WWP2—WW domain containing E3 ubiquitin protein ligase 2 (rs34195470), NFAT5—nuclear factor of activated T cells 5 (rs6499244), TGFA—transforming growth factor alpha (rs3771501), GDF5—growth differentiation factor 5 (rs143384), were selected for the genetic study based on previously registered GWAS associations (P ≤ 5 × 10–8) of these loci with KOA in European populations [40,41,42,43,44] (Table S1, Supplementary Information) and existence of the functional value [45,46,47]. To determine the loci functionality, the HaploReg database was used [48] (Table S2, Supplementary Information).

Genomic DNA of participants was isolated from peripheral blood (The buffy coat containing leucocytes was used.) by using a well-established phenol–chloroform-ethanol extraction/concentration method based on a previously published laboratory protocol [49]. The purity and concentration of the isolated DNA samples were measured on a NanoDrop spectrophotometer [50]. DNA materials of KOA patients and controls (each PCR tablet contained DNA samples of patients and controls) were genotyped by real-time PCR on the CFX96-Real-Time PCR System (Bio-Rad Laboratories, Hercules, CA, USA) [51, 52] and by using specially-developed reagent kits (TestGen, Ulyanovsk, Russia). The sequences of oligonucleotide primers and probes used in SNP genotyping are presented in Table S3 (Supplementary Information). To control the quality of experimental data, ≈7%–10% of the randomly selected DNA specimen were re-genotyped [53, 54]. A virtually complete coincidence was achieved between the repeated genotyping results with the primary data (an error of no more than 1%).

Statistical and bioinformatic analysis

For all the considered loci in KOA patients and controls in the two study subgroups (BMI ≥ 30 and BMI < 30), we evaluated the correspondence of the observed genotype distribution to the expected one according to the Hardy–Weinberg pattern [55, 56]. The association between SNP and KOA was investigated in the two groups (BMI ≥ 30 and BMI < 30) separately by using the logistic regression (allelic/additive/dominant/recessive genetic models [57] were considered), with adjustments made for age, sex, BMI, occupation-related physical workload, hereditary burden, the presence of concomitant pathology of the cardiovascular, musculoskeletal systems, height and leisure time physical activity (Table 1). All calculations were carried out by employing the gPLINK software [58] and were subjected to calibration for multiple comparisons (a well-established permutation test was applied) [59, 60]. Finally, a Pperm. ≤ 0.025 was considered to be statistically significant (Bonferroni correction was introduced for the number of groups compared (n = 2)—with/without obesity) [61]. For individual SNPs, statistical power was estimated by utilizing Quanto (v.1.2.4) [62].

Table 1 Phenotypic characteristics of the study participants
Full size table

Functionality of KOA-associated loci (epigenetic; eQTL; sQTL; protein structure change (amino acid substitution) [63]) and SNPs strongly linked with them (parameter r2 ≥ 0.80 [64]) were estimated by using modern bioinformatic online resources (in silico procedures) [65,66,67]: (a) Blood eQTL browser [68], (b) PolyPhen-2 [69], (c) GeneMANIA [70], (d) HaploReg [48], (e) GTExproject [71], (f) SIFT [72].

Results

The main phenotypic parameters of KOA patients and KOA-free individuals in the two groups, grouped in terms of the presence/absence of obesity (BMI ≥ 30 and BMI < 30) are given in Table 1. It was found that the KOA patients with BMI ≥ 30, as well as those with BMI < 30, compared with their corresponding controls, had significantly higher BMI (P = 0.0001 and P < 1 × 10–6, respectively), hereditary burden (P = 0.0005 and P = 0.0005), incidences of cardiovascular (P = 0.0005 and P = 0.0005) and musculoskeletal diseases (P = 0.0006 and P = 0.0005) diseases. Among the KOA subjects (BMI ≥ 30 and BMI < 30), in comparison with their respective controls, the percentage of individuals with a high level of professional physical activity was significantly higher (1.89 times, P = 0.002, and 2.05 times, P = 0.0005) and the proportion of individuals with a low level of professional physical activity was significantly lower (2.42 times, P = 0.0005, and 1.87 times, P = 0.0005, respectively). Additionally, in KOA patients without obesity (BMI < 30), the percentage of individuals with low physical activity in their free time was significantly higher (1.25 times, P = 0.001) and the proportion of individuals with regular physical activity was significantly lower (3 times, P = 0.002), compared to the controls (Table 1). The above-mentioned environmental KOA risk/protective factors were included in the association analysis as covariates.

The statistical materials in Table S4 (Supplementary Information) (BMI < 30 cohort) and Table S5 (Supplementary Information) (BMI ≥ 30 subject) demonstrate that the distribution (observed/expected) of the studied SNPs followed the HWE law (the Bonferroni correction based on the number of examined loci was used (Pbonf. = 0.00625 [0.05/8]).

Multidirectional relationships of the rs143384 GDF5 with KOA in BMI-different groups were found: allele G of this SNP was a KOA protective genetic variant in individuals with BMI ≥ 30 (OR 0.41 [95%CI 0.20–0.94], P = 0.019, Pperm. = 0.020, power 87.23%, recessive model) and was a disease risk variant in subjects with BMI < 30 (OR 1.32 [95%CI 1.05–1.65], P = 0.016, Pperm. = 0.018, allele model; OR 1.44 [95%CI 1.10–1.86], P = 0.007, Pperm. = 0.009, power 89.33%, additive model; OR 1.67 [95%CI 1.10–2.52], P = 0.015, P perm. = 0.012, power 81.01%, dominant model) (Table 2).

Table 2 Associations of the studied gene polymorphisms with KOA in subjects with BMI < 30 and BMI ≥ 30
Full size table

Functionality of KOA-associated rs143384 GDF5 (in silico data)

The polymorphism rs143384 (located in the 5’-UTR region of the GDF5 gene) and 9 SNPs strongly linked to it exhibited various epigenetic effects (They are significant for the chromatin structure in the regions of potential promoters and enhancers, and affect the interaction of DNA with many transcription factors such as Ascl2, Foxa, TFE, Ets, Pitx2, SP2, LUN-1, EBF, Mxi1, Myf, Myc, NRSF, TAL1, YY1, Zfx, E2A, ELF1, etc.) (Table 3), including cell cultures of adipose (adipose derived mesenchymal stem cells, epigenomeID-E025/mesenchymal stem cells derived adipocyte cultured cells, epigenomeID-E023/nuclei of adipose, epigenome ID-E063) (Data were obtained from the Haploreg database [48]).

Table 3 Regulatory effects of the KOA-associated polymorphism rs143384 and SNPs in high LD (r2 ≥ 0.80) (data of the Haploreg database)
Full size table

The Blood eQTL browser showed that the minor allele G rs143384 is associated (PFDR = 0) with reduced mRNA level of UQCC (Z parameter -6.35) and CEP250 (-5.74) genes and a high production of EIF6 mRNA (11.29) in peripheral blood (Table S6, Supplementary Information). In addition, the involvement of the three loci (rs6060402, rs224329, rs224333) highly coupled with rs143384 in transcriptional regulation of the above three genes in peripheral blood was displayed in Table S7 (Supplementary Information).

Based on experimental data of GTEx portal, rs143384 GDF5 has been identified as a modulator of multiple genes expression (21 genes/more 30 organs) and alternative splicing (8 genes/above 20 organs), including eight genes in adipose tissue (expression quantitative locus/trait [eQTL]: CPNE1, EDEM2, GDF5, PROCR, RPL36P4, UQCC1; splicing quantitative locus/trait [sQTL]: RBM39, UQCC1, ERGIC3) (Tables S8 and S9, Supplementary Information). Remarkably, the G allele of the rs143384 locus was correlated with low expression/splicing of four/two genes (CPNE1, EDEM2, PROCR, UQCC1/ERGIC3, RBM39) in adipose tissue and high expression/splicing of two/one genes (GDF5, RPL36P4/UQCC1) in this tissue (Tables S8 and S9, Supplementary Information). Among nine high-linked SNPs, eight loci were eQTL (21 genes including seven genes in adipose tissue: CEP250, CPNE1, EDEM2, PROCR, RP4-614O4.13, RPL36P4, UQCC1) (Table S10, Supplementary Information) and sQTL (10 genes including five genes in adipose tissue: EIF6, ERGIC3, FER1L4, RBM39, UQCC1) (Table S11, Supplementary Information).

Overall, first of all, we found very pronounced rs143384 GDF5 functionality in relation to 26 genes in a variety of organs (more than thirty ones) which interactions due co-expression (the percentage contribution was the highest and amounted to 85.77%), physical interactions (12.16%) and co-localization (2.07%) (Fig. 1, GeneMANIA data) with the leading role of paired interactions such as LAP3–RBM39, NQO2–NQO1, TRPC4AP–MYH7B, BRD2–EPB41L1, DPM3–CEP250 (weight indicators 0.21–0.62) (Table S12, Supplementary Information). Secondly, considerable functionality of the rs143384 GDF5 in adipose tissue in relation to nine genes (CPNE1, EDEM2, ERGIC3, GDF5, GDF5OS, PROCR, RBM39, RPL36P4, UQCC1) was found with complete dominance (100%) of co-expression in their interactions (Fig. 2 and Table S13, GeneMANIA data) and involved above genes set in regulation of the pathways of apoptosis of striated muscle cells (PFDR = 0.004).

Fig. 1
figure 1

The interaction networks of the candidate genes associated with rs143384 (eQTL/sQTL/regulatory effects this SNP) inferred using GeneMANIA (http://genemania.org)

Full size image
Fig. 2
figure 2

The interaction networks of the candidate genes in adipose tissue associated with rs143384 (eQTL/sQTL/regulatory effects this SNP) inferred using GeneMANIA (http://genemania.org)

Full size image

Discussion

In the present study, the effect of obesity on the association of the rs143384 GDF5 with KOA was shown: allele G of this SNP was a KOA protective factor in individuals with BMI ≥ 30 (OR 0.41) and disease risk marker in individuals with BMI < 30 (OR 1.32–1.67). Polymorphism rs143384 GDF5 exerted its regulatory effects in relation to 9 genes in adipose tissue.

A multitude of literature data confirmed that high BMI and obesity are the leading risk factors for the development and progression of OA [29, 73, 74]. In the sample we studied, obesity was also a significant risk factor for KOA (OR = 6.73, P = 0.005). It is known that the key points in the pathogenesis of KOA in obesity are determined by excessive mechanical stress on the joint, chronic inflammation in adipose tissue and dyslipidemia, secretion of proinflammatory cytokines and adipokines by adipose tissue; cytokine secretion by infrapatellar adipose tissue [2]. Adipokines (leptin, resistin, etc.) and cytokines (TNFA, IL1, IL6) produced by adipose tissue (both local and systemic adipose tissue), in turn, can affect pathological processes in the tissues of the joint and bones, such as cartilage degradation, inflammatory processes in the synovial membrane, bone erosion [2, 75, 76]. It is worth noting that infrapatellar adipose tissue or Goff’s fat cushion plays important roles in the pathogenesis of KOA [2]. On the one hand, the damping role of this adipose tissue is known to be due to the damping of mechanical stress under load on the joint (A redistribution of adipose tissue takes place.) [77]. On the other hand, KOA is often accompanied by inflammation of infrapatellar adipose tissue with increased expression of inflammatory mediators such as IL-6, adipsin, visfatin and adiponectin, which also support inflammatory processes in other joint tissues [78, 79]. It has been shown that an increase in body mass index by 5 kg/m2 is associated with a rise of 35% in the risk of developing KOA [80]. It is known that obesity often leads to the progression of KOA [17] and can cause a more severe course of the disease [81]. A study by Vasilic-Brasnjevic et al. showed that obesity (BMI ≥ 30 kg/m2) was a risk factor for the development of stage 3–4 KOA in patients older than 50 years [81]. The presence of overweight, grade I and II obesity increased the risk of KOA by 2, 3.1 times and 4.7 times, respectively [73]. Takahashi et al. [17] demonstrated that 75% of obese KOA patients had disease progression (assessed on the Kellgren-Lawrence scale).

OA and obesity are diseases resulting from the interaction of multiple genetic and environmental factors [82,83,84,85] and sharing common pathophysiological mechanisms [29]. Obesity is considered a chronic inflammatory disease characterized by the production of cytokines and cytokine-like molecules (adipokines) that can affect various body tissues [86], including knee joint tissues. It is worth noting that the inflammatory reaction is also one of the main pathogenetic links of OA [2, 87]. The relationship between different genetic variants and BMI in KOA patients has been shown [21,22,23,24,25,26,27,28], among which there are a large number of genes associated with metabolic disorders (FTO, ADIPOQ, LEP, SREBP2) and other genes (GDF5, TGFB1, etc.). At the same time, the effect of overweight and obesity on the association of GWAS-significant loci (rs8044769 FTO and rs143383 GDF5) with KOA has been examined only in a small number of studies [25]. It should be noted that no significant association was found between the GWAS locus rs8044769 FTO and KOA in overweight and obese patients [25]. Conversely, Zhang et al. showed that the rs143383 GDF5 was associated with KOA both in subjects with BMI ≥ 24 (OR = 2.36–2.45) and in those with BMI < 24 (OR = 1.63–3.77). Thus, the allele T of the rs143383 GDF5 was a risk factor for KOA in both groups [28]. According to our data, the variant G rs143384 GDF5 was a KOA risk factor in individuals with BMI < 30 (OR = 1.32–1.67) and a protective factor against KOA in BMI ≥ 30 subjects (OR = 0.41). The association of GDF5 gene (rs143384) with KOA was established in four previously published GWAS [40, 42,43,44]. Two papers [40, 43] reported the association of the allelic variant A (rs143384) with KOA in Europeans (parameter OR = 1.10 was the same in both studies), and one study [44] showed the association in mixed samples of European and Asian origins (OR = 1.07). It is worth noting that, in these three GWAS, the allele A of the GDF5 locus (rs143384) is risky for the development of KOA. Another study [42] demonstrated that the allele G of rs143384 was a KOA protective factor in Caucasians (OR = 0.91). It should be noted that, in our work, the allele G of rs143384 also was of protective value for KOA in the BMI ≥ 30 group.

There are a number of studies demonstrating the relationship between the rs143384 GDF5 and various musculoskeletal pathologies of the lower extremities, including OA of other sites or body parts [88,89,90,91,92,93,94]. Some studies have identified associations of the rs143384 allele variant A with knee pain [90, 92, 94]. Other studies have shown the connection of this GDF5 gene locus with hip dysplasia [89, 91], OA of the hand [93], and congenital hip dislocation [88].

The relationship between rs143384 GDF5 and body weight, as well as various anthropometric indices (body fat distribution, waist-to-hip ratio, waist-hip index, etc.), which may be associated with overweight or obesity, was demonstrated in previously GWAS [95,96,97,98,99,100,101]. The G allele (rs143384) has been found to be associated with lower body fat distribution (leg fat ratio) (β = -0.031, P = 3 × 10–43) [98], waist-to-hip ratio adjusted for BMI (β = -0.035, P = 3 × 10–28) [99], waist-hip index (β = -0.031, P = 6 × 10–23) [99]; in turn, allele A (rs143384) was linked to a higher waist-to-hip ratio adjusted for BMI (β = 0.02, P = 2 × 10–27) [97]. On the contrary, other studies [95, 96] showed that the A allele of rs143384 was correlated with a lower hip circumference adjusted for BMI (β = -0.044, P = 1 × 10–31) [95], (β = -0.042, P = 3 × 10–7) [96]. Association of rs143384 GDF5 with body weight has been shown in two papers [96; 100], in which the G allele was associated with weight gain (β = 0.028, P = 3 × 10–57) in the mixed samples of Europeans and Asians [100], and the A allele had a link with weight loss (β = -0.041, P = 2 × 10–10) in Europeans [96]. Hübel et al. revealed that rs143384 GDF5 was associated with fat-free muscle mass (β = -0.390, P = 6 × 10–68) [101] and a study by Guilherme et al. found that the G rs143384 allele of the GDF5 gene was associated with a low BMI in Caucasians (P = 1.2 × 10–14) [102]. Thus, it should be mentioned that, on the one hand, the association of rs143384 GDF5 with various anthropometric characteristics was proven in several previously GWAS; on the other hand, there is inconsistencies among the results about the association of this allelic variant with the aforementioned characteristics (risk/protective effect on BMI/body fat distribution/waist-to-hip ratio of different allelic variants of rs143384) in various cohorts (populations). Our study also revealed a multidirectional nature of the association between the rs143384 GDF5 and KOA association (the risky nature in individuals with BMI < 30 and the protective role in the group with BMI ≥ 30).

Interestingly, this study (in silico materials) demonstrated that, in adipose tissue, rs143384 GDF5 had considerable functionality (expression; splicing; epigenetic) in relation to nine genes (GDF5, CPNE1, EDEM2, ERGIC3, GDF5OS, PROCR, RBM39, RPL36P4, UQCC1) involved in regulation of pathways of apoptosis of striated muscle cells. Moreover, the G allele rs143384 was associated with increased GDF5 gene expression. Premised on this, it can be assumed that in obese individuals with the G allele, the amount of the protein product of the GDF5 gene will be maximum (plenty of adipose tissue due to GDF5 production area and the presence of a highly productive allele G rs143384) and significantly exceed the level of GDF5 expression (GDF5 production) in obese individuals without the G allele (a lot of adipose tissue but the presence of a low-productive allele A rs143384). This may explain the protective value for KOA of the highly productive allele G rs143384 in obese individuals, established in our study. GDF5 (growth differentiation factor 5) is a member of the bone morphogenetic protein (BMP) gene family and the TGF-beta superfamily and plays an important role in skeletal development [103], inflammatory reactions, and tissue damage [104]. Overexpression of GDF5 in human mesenchymal stem cells leads to increased chondrogenesis in vitro [105]. In mice models of OA, high expression of GDF5 in the cartilage was detected during its recovery after unilateral destabilization of the medial meniscus [106]. Allelic variants A and G rs143384 exert an important modifying effect on the KOA-risk impact of other loci. It was revealed that the T allele rs143383 (It was associated with the OA risk), which is linked to rs143384 (r2 = 0.82), caused the reduced transcription of the GDF5 gene in chondrogenic cells [107,108,109]. It has been shown that the rs143384 locus is able to influence the “phenotypic effects” of rs143383 with respect to the GDF5 gene (These results were obtained on the model of the luciferase reporter assays of GDF5 promoter/5’-UTR constructs in the chondrogenic (CH8), adipogenic (SW872) and osteogenic (MG63) cell lines). The T allele of rs143383, which is risky for OA, causes a decrease in luciferase activity relative to the alternative allele C for it only in the presence of the A allele of rs143384 [110]. Increased expression of the GDF5 gene was observed in brown adipose tissue in obese mice [111]. The study by Yang et al. showed that systemic overexpression of GDF5 in adipocytes reduced non-alcoholic liver obesity caused by a high-fat diet in mice [112]. Pei et al. exhibited that GDF5 played an adipogenic role in the differentiation of 3T3-L1 preadipocytes [113]. Thus, the GDF5 gene is characterized by a pleotropic effect and, accordingly, affects not only KOA, but also the processes taking place in adipose tissue, which is just as important, if indirectly, for the pathophysiology of KOA. In general, as we can assume in obese individuals, the highly productive G allele rs143384 (which determines overexpression of GDF5) acts as a protective factor against KOA due to the apparent effects of high concentrations of GDF5 (increased chondrogenesis, etc.).

At the same time, we obtained data on the KOA-risk role of polymorphism rs143384 GDF5 (allele G) in non-obesity individuals. We speculate that this relationship may be based on the following mechanisms. Firstly, a significant disadvantage of expression of GDF5 in individuals with a low content of adipose tissue (little adipose tissue due to small source of production of GDF5) and consequently weak chondrogenic and adipogenic effects of GDF5 led to an increased risk of developing KOA. Secondly, in individuals with low fat mass, an increased risk of KOA development in the presence of the G allele rs143384 may be associated with other genes whose expression/splicing level is affected by this polymorphism (CPNE1, EDEM2, PROCR, UQCC1, RPL36P4/ERGIC3, RBM39). For instance, due to a significant “deficiency” of protein products of genes (CPNE1, EDEM2, PROCR, UQCC1), their expression can be extremely reduced in individuals carrying reduced fat mass and a low-productive G allele rs143384 (KOA risk factor in individuals without obesity), which may, as an important pathogenetic factor, significantly contribute to the development of KOA. So, CPNE1, encoding Copine1, a soluble calcium-dependent membrane-binding protein, affects the length of myotubes (knockdown of CPNE1 gene increases the length of myotubes) and works as a modifier of muscle mass in humans in vitro, though it is not definitively clear how alterations in myogenesis indicators in vitro relate to the hypertrophy/hyperplasia of fiber in vivo [114]. EDEM2 encodes an ER degradation enhancer, mannosidase alpha-like 2, involved in carbohydrate metabolism (EDEM 2 identifies misfolded endoplasmic reticulum glycoproteins and targets them for destruction), and its expression in the skeletal muscle tissue of geriatric vs. young adult animals (dogs) differed significantly, depending on the diet [115]. PROCR, encoding the endothelial protein C receptor), is a “key” regulator of the protein C pathway mediating the interaction between coagulation and pro-inflammatory/anti-inflammatory processes in vessels [116]. QCC1 encodes a trans-membrane protein ubiquinol-cytochrome-c reductase complex chaperone and is involved in the pathophysiology of OA [117]. Moreover, it is also associated with lipid metabolism (arm fat mass) [118]. However, it is important to emphasize that there is currently no definitive or flimsy evidence on this issue in the literature, and further epidemiological and experimental studies on this theme are needed.

The data obtained in the work on the genetic features of KOA in individuals with and without obesity is traumatologically and orthopedically of practical value and can help distinguish between individuals at risk for KOA development and clinically healthy population. Taking into account the presence/absence of obesity will allow for timely implementation of measures aimed at preventing the disease (for example, achieving weight loss in obese individuals with a genetic high-risk factor for KOA (allele A rs143384), etc.).

Conclusion

This study showed that obesity exerted an effect on the associations of the rs143384 GDF5 with the KOA risk. This polymorphism is of “protective” value in the BMI ≥ 30 subjects and a “risk” for the development of KOA in those with BMI < 30.

Availability of data and materials

The data generated in the present study are available from the corresponding author upon reasonable request.

Abbreviations

BMI:

Body Mass Index

GWAS:

Genome-Wide Association Studies

KOA:

Knee Osteoarthritis

SNP:

Single Nucleotide Polymorphism

References

  1. Poole AR. Osteoarthritis as a whole joint disease. HSS J. 2012;8(1):4–6. https://doi.org/10.1007/s11420-011-9248-6.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Chang J, Liao Z, Lu M, Meng T, Han W, Ding C. Systemic and local adipose tissue in knee osteoarthritis. Osteoarthritis Cartilage. 2018;26(7):864–71. https://doi.org/10.1016/j.joca.2018.03.004.004.

    Article  PubMed  CAS  Google Scholar 

  3. Cui A, Li H, Wang D, Zhong J, Chen Y, Lu H. Global, regional prevalence, incidence and risk factors of knee osteoarthritis in population-based studies. EClinicalMedicine. 2020;29–30:100587. https://doi.org/10.1016/j.eclinm.2020.100587.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Leyland KM, Gates LS, Sanchez-Santos MT, Nevitt MC, Felson D, Jones G, et al. Knee osteoarthritis and time-to all-cause mortality in six community-based cohorts: an international meta-analysis of individual participant-level data. Aging Clin Exp Res. 2021;33(3):529–45. https://doi.org/10.1007/s40520-020-01762-2.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Spector AL, Matsen E, Egede LE. Trends and racial/ethnic differences in health care spending stratified by gender among adults with arthritis in the United States 2011–2019. Int J Environ Res Public Health. 2022;19(15):9014. https://doi.org/10.3390/ijerph19159014.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Muoh O, Malemud CJ, Askare AD. Clinical significance and implications of genetic and genomic studies in patients with osteoarthritis. Adv Genomics Genet. 2014;4:193–206. https://doi.org/10.2147/AGG.S64284.

    Article  Google Scholar 

  7. Dong Y, Zhang P, Fan L. Recognition of factors of postoperative complications of knee osteoarthritis patients and comprehensive nursing intervention. Comput Math Methods Med. 2021;2021:1840613. https://doi.org/10.1155/2021/1840613.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Healy WL, Rana AJ, Iorio R. Hospital economics of primary total knee arthroplasty at a teaching hospital. Clin Orthop Relat Res. 2011;469(1):87–94. https://doi.org/10.1007/s11999-010-1486-2.

    Article  PubMed  Google Scholar 

  9. Puig-Junoy J, Ruiz ZA. Socio-economic costs of osteoarthritis: a systematic review of cost-of-illness studies. Semin Arthritis Rheum. 2015;44(5):531–41. https://doi.org/10.1016/j.semarthrit.2014.10.012.

    Article  PubMed  Google Scholar 

  10. Das SK, Farooqi A. Osteoarthritis. Best Pract Res Clin Rheumatol. 2008;22(4):657–75. https://doi.org/10.1016/j.berh.2008.07.002.

    Article  PubMed  CAS  Google Scholar 

  11. Michael JW, Schlüter-Brust KU, Eysel P. The epidemiology, etiology, diagnosis, and treatment of osteoarthritis of the knee. Dtsch Arztebl Int. 2010;107(9):152–62. https://doi.org/10.3238/arztebl.2010.0152.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Georgiev T, Angelov AK. Modifiable risk factors in knee osteoarthritis: treatment implications. Rheumatol Int. 2019;39(7):1145–57. https://doi.org/10.1007/s00296-019-04290-z.

    Article  PubMed  Google Scholar 

  13. Santangelo KS, Radakovich LB, Fouts J, Foster MT. Pathophysiology of obesity on knee joint homeostasis: contributions of the infrapatellar fat pad. Horm Mol Biol Clin Investig. 2016;26(2):97–108. https://doi.org/10.1515/hmbci-2015-0067.

    Article  PubMed  CAS  Google Scholar 

  14. Belluzzi E, El Hadi H, Granzotto M, Rossato M, Ramonda R, Macchi V, et al. Systemic and local adipose tissue in knee osteoarthritis. J Cell Physiol. 2017;232(8):1971–8. https://doi.org/10.1002/jcp.25716.

    Article  PubMed  CAS  Google Scholar 

  15. Belluzzi E, Stocco E, Pozzuoli A, Granzotto M, Porzionato A, Vettor R, et al. Contribution of infrapatellar fat pad and synovial membrane to knee osteoarthritis pain. Biomed Res Int. 2019;2019:6390182. https://doi.org/10.1155/2019/6390182.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Blagojevic M, Jinks C, Jeffery A, Jordan KP. Risk factors for onset of osteoarthritis of the knee in older adults: a systematic review and meta-analysis. Osteoarthritis Cartilage. 2010;18(1):24–33. https://doi.org/10.1016/j.joca.2009.08.010.

    Article  PubMed  CAS  Google Scholar 

  17. Takahashi A, Umehara J, Kamimura M, Aizawa T, Itoi E. Obesity is a risk factor for osteoarthritis progression and spontaneous osteoporosis is a risk for the development of spontaneous osteonecrosis in patients with medial meniscus posterior root tear. J Orthop Sci. 2021;26(5):844–9. https://doi.org/10.1016/j.jos.2020.09.001.

    Article  PubMed  Google Scholar 

  18. Panoutsopoulou K, Southam L, Elliott KS, Wrayner N, Zhai G, Beazley C, et al. Insights into the genetic architecture of osteoarthritis from stage 1 of the arcOGEN study. Ann Rheum Dis. 2011;70(5):864–7. https://doi.org/10.1136/ard.2010.141473.

    Article  PubMed  CAS  Google Scholar 

  19. Aubourg G, Rice SJ, Bruce-Wootton P, Loughlin J. Genetics of osteoarthritis. Osteoarthritis Cartilage. 2022;30(5):636–49. https://doi.org/10.1016/j.joca.2021.03.002.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. GWAS Catalog. https://www.ebi.ac.uk/gwas/search?query=osteoarthritis,%20knee. Accessed 18 March 2023.

  21. Qin J, Shi D, Dai J, Zhu L, Tsezou A, Jiang Q. Association of the leptin gene with knee osteoarthritis susceptibility in a Han Chinese population: a case-control study. J Hum Genet. 2010;55(10):704–6. https://doi.org/10.1038/jhg.2010.86.

    Article  PubMed  Google Scholar 

  22. Elliott KS, Chapman K, Day-Williams A, Panoutsopoulou K, Southam L, Lindgren CM, et al. Evaluation of the genetic overlap between osteoarthritis with body mass index and height using genome-wide association scan data. Ann Rheum Dis. 2013;72(6):935–41. https://doi.org/10.1136/annrheumdis-2012-202081.

    Article  PubMed  Google Scholar 

  23. Muthuri SG, Doherty S, Zhang W, Maciewicz RA, Muir KR, Doherty M. Gene-environment interaction between body mass index and transforming growth factor beta 1 (TGFβ1) gene in knee and hip osteoarthritis. Arthritis Res Ther. 2013;15(2):R52. https://doi.org/10.1186/ar4214.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wang Y, Chu M, Rong J, Xing B, Zhu L, Zhao Y, et al. No association of the single nucleotide polymorphism rs8044769 in the fat mass and obesity-associated gene with knee osteoarthritis risk and body mass index: a population-based study in China. Bone Joint Res. 2016;5(5):169–74. https://doi.org/10.1302/2046-3758.55.2000589.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Dai J, Ying P, Shi D, Hou H, Sun Y, Xu Z, et al. FTO variant is not associated with osteoarthritis in the Chinese Han population: replication study for a genome-wide association study identified risk loci. J Orthop Surg Res. 2018;13(1):65. https://doi.org/10.1186/s13018-018-0769-2.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Jiang L, Zhu X, Rong J, Xing B, Wang S, Liu A, et al. Obesity, osteoarthritis and genetic risk: the rs182052 polymorphism in the ADIPOQ gene is potentially associated with risk of knee osteoarthritis. Bone Joint Res. 2018;7(7):494–500. https://doi.org/10.1302/2046-3758.77.BJR-2017-0274.R1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Poornima S, Subramanyam K, Khan IA, Sumanlatha G, Hasan Q. Role of SREBP2 gene polymorphism on knee osteoarthritis in the South Indian hyderabad population: a hospital based study with G595C variant. J Orthop. 2019;16(3):293–7. https://doi.org/10.1016/j.jor.2019.05.001.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhang S, Wang J, Ji H, Jia H, Guan D. Interaction between GDF5 gene polymorphisms and environment factors increased the risk of knee osteoarthritis: a case-control study. Biosci Rep. 2019;39(2):BSR20182423. https://doi.org/10.1042/BSR20182423.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Kulkarni K, Karssiens T, Kumar V, Pandit H. Obesity and osteoarthritis. Maturitas. 2016;89:22–8. https://doi.org/10.1016/j.maturitas.2016.04.006.

    Article  PubMed  Google Scholar 

  30. Ryzhkov II, Borzilov EE, Churnosov MI, Ataman AV, Dedkov AA, Polonikov AV. Transforming growth factor beta 1 is a novel susceptibility gene for adolescent idiopathic scoliosis. Spine (Phila Pa 1976). 2013;38:E699–704. https://doi.org/10.1097/BRS.0b013e31828de9e1.

    Article  PubMed  Google Scholar 

  31. Pavlova N, Demin S, Churnosov M, Reshetnikov E, Aristova I, Churnosova M, et al. The modifying effect of obesity on the association of matrix metalloproteinase gene polymorphisms with breast cancer risk. Biomedicines. 2022;10(10):2617. https://doi.org/10.3390/biomedicines10102617.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Reshetnikov E, Ponomarenko I, Golovchenko O, Sorokina I, Batlutskaya I, Yakunchenko T, et al. The VNTR polymorphism of the endothelial nitric oxide synthase gene and blood pressure in women at the end of pregnancy. Taiwan J Obstet Gynecol. 2019;58(3):390–5. https://doi.org/10.1016/j.tjog.2018.11.035.

    Article  PubMed  Google Scholar 

  33. Golovchenko IO. Genetic determinants of sex hormone levels in endometriosis patients. Res Results Biomed. 2023;9:5–21. https://doi.org/10.18413/2658-6533-2023-9-1-0-1. In Russian.

    Article  Google Scholar 

  34. Reshetnikova Y, Churnosova M, Stepanov V, Bocharova A, Serebrova V, Trifonova E, et al. Maternal age at menarche gene polymorphisms are associated with offspring birth weight. Life (Basel). 2023;13(7):1525. https://doi.org/10.3390/life13071525.

    Article  ADS  PubMed  CAS  Google Scholar 

  35. Tikunova E, Ovtcharova V, Reshetnikov E, Dvornyk V, Polonikov A, Bushueva O, et al. Genes of tumor necrosis factors and their receptors and the primary open angle glaucoma in the population of Central Russia. Int J Ophthalmol. 2017;10(10):1490–4. https://doi.org/10.18240/ijo.2017.10.02.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Novakov V, Novakova O, Churnosova M, Sorokina I, Aristova I, Polonikov A, et al. Intergenic Interactions of SBNO1, NFAT5 and GLT8D1 determine the susceptibility to knee osteoarthritis among Europeans of Russia. Life (Basel). 2023;13(2):405. https://doi.org/10.3390/life13020405.

    Article  ADS  PubMed  CAS  Google Scholar 

  37. Altman R, Asch E, Bloch D, Bole G, Borenstein D, Brandt K, et al. Development of criteria for the classification and reporting of osteoarthritis. Classification of osteoarthritis of the knee. Diagnostic and Therapeutic Criteria Committee of the American Rheumatism Association. Arthritis Rheum. 1986;29(8):1039–49. https://doi.org/10.1002/art.1780290816.

    Article  PubMed  CAS  Google Scholar 

  38. Kellgren JH, Lawrence JS. Radiological assessment of osteoarthrosis. Ann Rheum Dis. 1957;16(4):494–502. https://doi.org/10.1136/ard.16.4.494.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Hawker GA, Mian S, Kendzerska T, French M. Measures of adult pain: Visual Analog Scale for Pain (VAS Pain), Numeric Rating Scale for Pain (NRS Pain), McGill Pain Questionnaire (MPQ), Short-Form McGill Pain Questionnaire (SF-MPQ), Chronic Pain Grade Scale (CPGS), Short Form-36 Bodily Pain Scale (SF-36 BPS), and Measure of Intermittent and Constant Osteoarthritis Pain (ICOAP). Arthritis Care Res (Hoboken). 2011;63(11):S240–52. https://doi.org/10.1002/acr.20543.

    Article  PubMed  Google Scholar 

  40. Styrkarsdottir U, Lund SH, Thorleifsson G, Zink F, Stefansson OA, Sigurdsson JK, et al. Meta-analysis of Icelandic and UK data sets identifies missense variants in SMO, IL11, COL11A1 and 13 more new loci associated with osteoarthritis. Nat Genet. 2018;50(12):1681–7. https://doi.org/10.1038/s41588-018-0247-0.

    Article  PubMed  CAS  Google Scholar 

  41. Zengini E, Hatzikotoulas K, Tachmazidou I, Steinberg J, Hartwig FP, Southam L, et al. Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis. Nat Genet. 2018;50(4):549–58. https://doi.org/10.1038/s41588-018-0079-y.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Styrkarsdottir U, Stefansson OA, Gunnarsdottir K, Thorleifsson G, Lund SH, Stefansdottir L, et al. GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures. Nat Commun. 2019;10(1):2054. https://doi.org/10.1038/s41467-019-09860-0.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  43. Tachmazidou I, Hatzikotoulas K, Southam L, Esparza-Gordillo J, Haberland V, Zheng J, et al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat Genet. 2019;51(2):230–6. https://doi.org/10.1038/s41588-018-0327-1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Boer CG, Hatzikotoulas K, Southam L, Stefánsdóttir L, Zhang Y, Coutinho de Almeida R, et al. Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell. 2021;184(18):4784-4818.e17. https://doi.org/10.1016/j.cell.2021.07.038.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Polonikov AV, Bushueva OY, Bulgakova IV, Freidin MB, Churnosov MI, Solodilova MA, et al. A comprehensive contribution of genes for aryl hydrocarbon receptor signaling pathway to hypertension susceptibility. Pharmacogenet Genomics. 2017;27(2):57–69. https://doi.org/10.1097/FPC.0000000000000261.

    Article  PubMed  CAS  Google Scholar 

  46. Ivanova T, Churnosova M, Abramova M, Plotnikov D, Ponomarenko I, Reshetnikov E, et al. Sex-specific features of the correlation between GWAS-noticeable polymorphisms and hypertension in Europeans of Russia. Int J Mol Sci. 2023;24(9):7799. https://doi.org/10.3390/ijms24097799.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Moskalenko M, Ponomarenko I, Reshetnikov E, Dvornyk V, Churnosov M. Polymorphisms of the matrix metalloproteinase genes are associated with essential hypertension in a Caucasian population of Central Russia. Sci Rep. 2021;11(1):5224. https://doi.org/10.1038/s41598-021-84645-4.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  48. Ward LD, Kellis M. HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016;44(D1):D877–81. https://doi.org/10.1093/nar/gkv1340.

    Article  PubMed  CAS  Google Scholar 

  49. Eliseeva N, Ponomarenko I, Reshetnikov E, Dvornyk V, Churnosov M. LOXL1 gene polymorphism candidates for exfoliation glaucoma are also associated with a risk for primary open-angle glaucoma in a Caucasian population from central Russia. Mol Vis. 2021;27:262–9.

    PubMed  PubMed Central  CAS  Google Scholar 

  50. Pavlova N, Demin S, Churnosov M, Reshetnikov E, Aristova I, Churnosova M, et al. Matrix metalloproteinase gene polymorphisms are associated with breast cancer in the Caucasian women of Russia. Int J Mol Sci. 2022;23(20):12638. https://doi.org/10.3390/ijms232012638.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Reshetnikov E, Zarudskaya O, Polonikov A, Bushueva O, Orlova V, Krikun E, et al. Genetic markers for inherited thrombophilia are associated with fetal growth retardation in the population of Central Russia. J Obstet Gynaecol Res. 2017;43(7):1139–44. https://doi.org/10.1111/jog.13329.

    Article  PubMed  CAS  Google Scholar 

  52. Moskalenko MI, Milanova SN, Ponomarenko IV, Polonikov AV, Churnosov MI. Study of associations of polymorphism of matrix metalloproteinases genes with the development of arterial hypertension in men. Kardiologiia. 2019;59(7S):31–9. https://doi.org/10.18087/cardio.2598. In Russian.

    Article  PubMed  CAS  Google Scholar 

  53. Golovchenko O, Abramova M, Ponomarenko I, Reshetnikov E, Aristova I, Polonikov A, et al. Functionally significant polymorphisms of ESR1and PGR and risk of intrauterine growth restriction in population of Central Russia. Eur J Obstet Gynecol Reprod Biol. 2020;253:52–7. https://doi.org/10.1016/j.ejogrb.2020.07.045.

    Article  PubMed  CAS  Google Scholar 

  54. Starikova D, Ponomarenko I, Reshetnikov E, Dvornyk V, Churnosov M. Novel data about association of the functionally significant polymorphisms of the MMP9 gene with exfoliation glaucoma in the Caucasian population of central Russia. Ophthalmic Res. 2021;64(3):458–64. https://doi.org/10.1159/000512507.

    Article  PubMed  CAS  Google Scholar 

  55. Bushueva O, Solodilova M, Churnosov M, Ivanov V, Polonikov A. The flavin-containing monooxygenase 3 gene and essential hypertension: the joint effect of polymorphism E158K and cigarette smoking on disease susceptibility. Int J Hypertens. 2014;2014:712169. https://doi.org/10.1155/2014/712169.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Minyaylo O, Ponomarenko I, Reshetnikov E, Dvornyk V, Churnosov M. Functionally significant polymorphisms of the MMP-9 gene are associated with peptic ulcer disease in the Caucasian population of Central Russia. Sci Rep. 2021;11(1):13515. https://doi.org/10.1038/s41598-021-92527-y.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  57. Golovchenko I, Aizikovich B, Golovchenko O, Reshetnikov E, Churnosova M, Aristova I, et al. Sex Hormone candidate gene polymorphisms are associated with endometriosis. Int J Mol Sci. 2022;23(22):13691. https://doi.org/10.3390/ijms232213691.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75. https://doi.org/10.1086/519795.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Che R, Jack JR, Motsinger-Reif AA, Brown CC. An adaptive permutation approach for genome-wide association study: evaluation and recommendations for use. BioData Min. 2014;7:9. https://doi.org/10.1186/1756-0381-7-9.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Ponomarenko I, Reshetnikov E, Polonikov A, Sorokina I, Yermachenko A, Dvornyk V, et al. Candidate genes for age at menarche are associated with endometrial hyperplasia. Gene. 2020;757:144933. https://doi.org/10.1016/j.gene.2020.144933.

    Article  PubMed  CAS  Google Scholar 

  61. Abramova M, Churnosova M, Efremova O, Aristova I, Reshetnikov E, Polonikov A, et al. Effects of pre-pregnancy overweight/obesity on the pattern of association of hypertension susceptibility genes with preeclampsia. Life (Basel). 2022;12(12):2018. https://doi.org/10.3390/life12122018.

    Article  ADS  PubMed  CAS  Google Scholar 

  62. Gauderman W, Morrison J. QUANTO 1.1: A computer program for power and sample size calculations genetic–epidemiology studies. 2006. http://hydra.usc.edu/gxe. Accessed 18 March 2023.

  63. Churnosov M, Abramova M, Reshetnikov E, Lyashenko IV, Efremova O, Churnosova M, et al. Polymorphisms of hypertension susceptibility genes as a risk factors of preeclampsia in the Caucasian population of central Russia. Placenta. 2022;129:51–61. https://doi.org/10.1016/j.placenta.2022.09.010.

    Article  PubMed  CAS  Google Scholar 

  64. Ivanova T, Churnosova M, Abramova M, Ponomarenko I, Reshetnikov E, Aristova I, et al. Risk Effects of rs1799945 Polymorphism of the HFE Gene and Intergenic Interactions of GWAS-significant loci for arterial hypertension in the Caucasian population of Central Russia. Int J Mol Sci. 2023;24(9):8309. https://doi.org/10.3390/ijms24098309.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Polonikov A, Kharchenko A, Bykanova M, Sirotina S, Ponomarenko I, Bocharova A, et al. Polymorphisms of CYP2C8, CYP2C9 and CYP2C19 and risk of coronary heart disease in Russian population. Gene. 2017;627:451–9. https://doi.org/10.1016/j.gene.2017.07.004.

    Article  PubMed  CAS  Google Scholar 

  66. Polonikov A, Rymarova L, Klyosova E, Volkova A, Azarova I, Bushueva O, et al. Matrix metalloproteinases as target genes for gene regulatory networks driving molecular and cellular pathways related to a multistep pathogenesis of cerebrovascular disease. J Cell Biochem. 2019;120(10):16467–82. https://doi.org/10.1002/jcb.28815.

    Article  PubMed  CAS  Google Scholar 

  67. Sirotina S, Ponomarenko I, Kharchenko A, Bykanova M, Bocharova A, Vagaytseva K, et al. A novel polymorphism in the promoter of the CYP4A11 gene is associated with susceptibility to coronary artery disease. Dis Markers. 2018;2018:5812802. https://doi.org/10.1155/2018/5812802.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. Westra HJ, Peters MJ, Esko T, Yaghootkar H, Schurmann C, Kettunen J, et al. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat Genet. 2013;45(10):1238–43. https://doi.org/10.1038/ng.2756.

  69. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;7:7–20. https://doi.org/10.1002/0471142905.hg0720s76.

    Article  Google Scholar 

  70. Franz M, Rodriguez H, Lopes C, Zuberi K, Montojo J, Bader GD, et al. GeneMANIA update 2018. Nucleic Acids Res. 2018;46(W1):W60–4. https://doi.org/10.1093/nar/gky311.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  71. GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science. 2020;369(6509):1318–30. https://doi.org/10.1126/science.aaz1776.

    Article  CAS  Google Scholar 

  72. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–81. https://doi.org/10.1038/nprot.2009.86.

    Article  PubMed  CAS  Google Scholar 

  73. Reyes C, Leyland KM, Peat G, Cooper C, Arden NK, Prieto-Alhambra D. Association between overweight and obesity and risk of clinically diagnosed knee, hip, and hand osteoarthritis: a population-based cohort study. Arthritis Rheumatol. 2016;68(8):1869–75. https://doi.org/10.1002/art.39707.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Khan B, Khan OY, Zehra S, Azhar A, Fatima S. Association between obesity and risk of knee osteoarthritis. Pak J Pharm Sci. 2020;33(1):295–8.

    PubMed  Google Scholar 

  75. Azamar-Llamas D, Hernández-Molina G, Ramos-Ávalos B, Furuzawa-Carballeda J. Adipokine contribution to the pathogenesis of osteoarthritis. Mediators Inflamm. 2017;2017:5468023. https://doi.org/10.1155/2017/5468023.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  76. Molnar V, Matisic V, Kodvanj I, Bjelica R, Jelec Z, Hudetz D, et al. Cytokines and chemokines involved in osteoarthritis pathogenesis. Int J Mol Sci. 2021;22(17):9208. https://doi.org/10.3390/ijms22179208.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Fontanella CG, Carniel EL, Frigo A, Macchi V, Porzionato A, Sarasin G, et al. Investigation of biomechanical response of Hoffa’s fat pad and comparative characterization. J Mech Behav Biomed Mater. 2017;67:1–9. https://doi.org/10.1016/j.jmbbm.2016.11.024.

    Article  PubMed  Google Scholar 

  78. Fontanella CG, Belluzzi E, Pozzuoli A, Favero M, Ruggieri P, Macchi V, et al. Mechanical behavior of infrapatellar fat pad of patients affected by osteoarthritis. J Biomech. 2022;131:110931. https://doi.org/10.1016/j.jbiomech.2021.110931.

    Article  PubMed  Google Scholar 

  79. Clockaerts S, Bastiaansen-Jenniskens YM, Runhaar J, Osch GJ, Offel JF, Verhaar JA, et al. The infrapatellar fat pad should be considered as an active osteoarthritic joint tissue: a narrative review. Osteoarthritis Cartilage. 2010;18:876–82. https://doi.org/10.1016/j.joca.2010.03.014.

    Article  PubMed  CAS  Google Scholar 

  80. Jiang L, Tian W, Wang Y, Rong J, Bao C, Liu Y, et al. Body mass index and susceptibility to knee osteoarthritis: a systematic review and meta-analysis. Joint Bone Spine. 2012;79(3):291–7. https://doi.org/10.1016/j.jbspin.2011.05.015.

    Article  PubMed  Google Scholar 

  81. Vasilic-Brasnjevic S, Marinkovic J, Vlajinac H, Vasiljevic N, Jakovljevic B, Nikic M, et al. Association of body mass index and waist circumference with severity of knee osteoarthritis. Acta Reumatol Port. 2016;41(3):226–31.

    PubMed  CAS  Google Scholar 

  82. Albuquerque D, Nobrega C, Manco L, Padez C. The contribution of genetics and environment to obesity. Br Med Bull. 2017;123(1):159–73. https://doi.org/10.1093/bmb/ldx022.

    Article  PubMed  Google Scholar 

  83. Snochowska A, Szmigielska P, Brzeziańska-Lasota E, Tomaszewski W. Genetic and epigenetic interactions in the etiopathogenesis of osteoarthritis. Selected molecular factors in OA etiopathogenesis. Ortop Traumatol Rehabil. 2017;19(3):227–37. https://doi.org/10.5604/15093492.1240791.

    Article  PubMed  Google Scholar 

  84. Stefik D, Vranic V, Ivkovic N, Abazovic D, Maric D, Vojvodic D, et al. An insight into osteoarthritis susceptibility: integration of immunological and genetic background. Bosn J Basic Med Sc. 2021;21(2):155–62. https://doi.org/10.17305/bjbms.2020.4735.

    Article  CAS  Google Scholar 

  85. Novakov VB, Novakova ON, Sorokina IN, Batlutskaya IV, Efremova OA, Orlova VS. Genetic markers of knee osteoarthritis in women of the Central Chernozem region of Russia. Res Results Biomed. 2023;9(2):191–205. https://doi.org/10.18413/2658-6533-2023-9-2-0-4.

    Article  Google Scholar 

  86. Coelho M, Oliveira T, Fernandes R. State of the art paper biochemistry of adipose tissue: an endocrine organ. Arch Med Sci. 2013;9(2):191–200. https://doi.org/10.5114/aoms.2013.33181.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  87. Nedunchezhiyan U, Varughese I, Sun AR, Wu X, Crawford R, Prasadam I. Obesity, inflammation, and immune system in osteoarthritis. Front Immunol. 2022;13:907750. https://doi.org/10.3389/fimmu.2022.907750.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  88. Rouault K, Scotet V, Autret S, Gaucher F, Dubrana F, Tanguy D, et al. Evidence of association between GDF5 polymorphisms and congenital dislocation of the hip in a Caucasian population. Osteoarthritis Cartilage. 2010;18(9):1144–9. https://doi.org/10.1016/j.joca.2010.05.018.

    Article  PubMed  CAS  Google Scholar 

  89. Hatzikotoulas K, Roposch A, Shah KM, Clark MJ, Bratherton S, Limbani V, et al. Genome-wide association study of developmental dysplasia of the hip identifies an association with GDF5. Commun Biol. 2018;1:56. https://doi.org/10.1038/s42003-018-0052-4.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  90. Meng W, Adams MJ, Palmer CAN, Shi J, Auton A, Ryan KA, et al. Genome-wide association study of knee pain identifies associations with GDF5 and COL27A1 in UK Biobank. Commun Biol. 2019;2:321. https://doi.org/10.1038/s42003-019-0568-2.

    Article  PubMed  PubMed Central  Google Scholar 

  91. Kenanidis E, Gkekas NK, Karasmani A, Anagnostis P, Christofilopoulos P, Tsiridis E. Genetic predisposition to developmental dysplasia of the hip. J Arthroplasty. 2020;35(1):291-300.e1. https://doi.org/10.1016/j.arth.2019.08.031.

    Article  PubMed  Google Scholar 

  92. Tsepilov YA, Freidin MB, Shadrina AS, Sharapov SZ, Elgaeva EE, Zundert JV, et al. Analysis of genetically independent phenotypes identifies shared genetic factors associated with chronic musculoskeletal pain conditions. Commun Biol. 2020;3(1):329. https://doi.org/10.1038/s42003-020-1051-9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  93. Boer CG, Yau MS, Rice SJ, Coutinho de Almeida R, Cheung K, Styrkarsdottir U, et al. Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene. Ann Rheum Dis. 2021;80(3):367–75. https://doi.org/10.1136/annrheumdis-2020-217834.

    Article  PubMed  CAS  Google Scholar 

  94. Yan S, Nie H, Bu G, Yuan W, Wang S. The effect of common variants in GDF5 gene on the susceptibility to chronic postsurgical pain. J Orthop Surg Res. 2021;16(1):420. https://doi.org/10.1186/s13018-021-02549-5.

    Article  PubMed  PubMed Central  Google Scholar 

  95. Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Mägi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518(7538):187–96. https://doi.org/10.1038/nature14132.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  96. Tachmazidou I, Süveges D, Min JL, Ritchie GRS, Steinberg J, Walter K, et al. Whole-genome sequencing coupled to imputation discovers genetic signals for anthropometric traits. Am J Hum Genet. 2017;100(6):865–84. https://doi.org/10.1016/j.ajhg.2017.04.014.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  97. Lotta LA, Wittemans LBL, Zuber V, Stewart ID, Sharp SJ, Luan J, et al. Association of genetic variants related to gluteofemoral vs abdominal fat distribution with type 2 diabetes, coronary disease, and cardiovascular risk factors. JAMA. 2018;320(24):2553–63. https://doi.org/10.1001/jama.2018.19329.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  98. Rask-Andersen M, Karlsson T, Ek WE, Johansson A. Genome-wide association study of body fat distribution identifies adiposity loci and sex-specific genetic effects. Nat Commun. 2019;10(1):339. https://doi.org/10.1038/s41467-018-08000-4.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  99. Christakoudi S, Evangelou E, Riboli E, Tsilidis KK. GWAS of allometric body-shape indices in UK Biobank identifies loci suggesting associations with morphogenesis, organogenesis, adrenal cell renewal and cancer. Sci Rep. 2021;11(1):10688. https://doi.org/10.1038/s41598-021-89176-6.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  100. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. 2021;53(10):1415–24. https://doi.org/10.1038/s41588-021-00931-x.

    Article  PubMed  CAS  Google Scholar 

  101. Hübel C, Gaspar HA, Coleman JRI, Finucane H, Purves KL, Hanscombe KB, et al. Genomics of body fat percentage may contribute to sex bias in anorexia nervosa. Am J Med Genet B Neuropsychiatr Genet. 2019;180(6):428–38. https://doi.org/10.1002/ajmg.b.32709.

    Article  PubMed  CAS  Google Scholar 

  102. Guilherme JPLF, Semenova EA, Larin AK, Yusupov RA, Generozov EV, Ahmetov II. Genomic predictors of brisk walking are associated with elite sprinter status. Genes. 2022;13(10):1710. https://doi.org/10.3390/genes13101710.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  103. Wu DD, Li GM, Jin W, Li Y, Zhang YP. Positive selection on the osteoarthritis-risk and decreased-height associated variants at the GDF5 gene in East Asians. PLoS One. 2012;7(8):e42553. https://doi.org/10.1371/journal.pone.0042553.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  104. Chen W, Meng J, Qian H, Deng Z, Chen S, Xu H, et al. A Study of IL-1β, MMP-3, TGF-β1, and GDF5 polymorphisms and their association with primary frozen shoulder in a Chinese Han population. Biomed Res Int. 2017;2017:3681645. https://doi.org/10.1155/2017/3681645.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  105. Coleman CM, Vaughan EE, Browe DC, Mooney E, Howard L, Barry F. Growth differentiation factor-5 enhances in vitro mesenchymal stromal cell chondrogenesis and hypertrophy. Stem Cells Dev. 2013;22(13):1968–76. https://doi.org/10.1089/scd.2012.0282.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  106. Kania K, Colella F, Riemen AHK, Wang H, Howard KA, Aigner T, et al. Regulation of Gdf5 expression in joint remodelling, repair and osteoarthritis. Sci Rep. 2020;10(1):157. https://doi.org/10.1038/s41598-019-57011-8.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  107. Miyamoto Y, Mabuchi A, Shi D, Kubo T, Takatori Y, Saito S, et al. A functional polymorphism in the 5’ UTR of GDF5 is associated with susceptibility to osteoarthritis. Nat Genet. 2007;39(4):529–33. https://doi.org/10.1038/2005.

    Article  PubMed  CAS  Google Scholar 

  108. Southam L, Rodriguez-Lopez J, Wilkins JM, Pombo-Suarez M, Snelling S, Gomez-Reino JJ, et al. An SNP in the 5’-UTR of GDF5 is associated with osteoarthritis susceptibility in Europeans and with in vivo differences in allelic expression in articular cartilage. Hum Mol Genet. 2007;16(18):2226–32. https://doi.org/10.1093/hmg/ddm174.

    Article  PubMed  CAS  Google Scholar 

  109. Chapman K, Takahashi A, Meulenbelt I, Watson C, Rodriguez-Lopez J, Egli R, et al. A meta-analysis of European and Asian cohorts reveals a global role of a functional SNP in the 5’ UTR of GDF5 with osteoarthritis susceptibility. Hum Mol Genet. 2008;17(10):1497–504. https://doi.org/10.1093/hmg/ddn038.

    Article  PubMed  CAS  Google Scholar 

  110. Egli RJ, Southam L, Wilkins JM, Lorenzen I, Pombo-Suarez M, Gonzalez A, et al. Functional analysis of the osteoarthritis susceptibility-associated GDF5 regulatory polymorphism. Arthritis Rheum. 2009;60(7):2055–64. https://doi.org/10.1002/art.24616.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  111. Hinoi E, Nakamura Y, Takada S, Fujita H, Iezaki T, Hashizume S, et al. Growth differentiation factor-5 promotes brown adipogenesis in systemic energy expenditure. Diabetes. 2014;63(1):162–75. https://doi.org/10.2337/db13-0808.

    Article  PubMed  CAS  Google Scholar 

  112. Yang Y, Zhang W, Wu X, Wu J, Sun C, Luo F, et al. Systemic overexpression of GDF5 in adipocytes but not hepatocytes alleviates high-fat diet-induced nonalcoholic fatty liver in mice. Can J Gastroenterol Hepatol. 2021;2021:8894685. https://doi.org/10.1155/2021/8894685.

    Article  PubMed  PubMed Central  Google Scholar 

  113. Pei Z, Yang Y, Kiess W, Sun C, Luo F. Dynamic profile and adipogenic role of growth differentiation factor 5 (GDF5) in the differentiation of 3T3-L1 preadipocytes. Arch Biochem Biophys. 2014;560:27–35. https://doi.org/10.1016/j.abb.2014.07.025.

    Article  PubMed  CAS  Google Scholar 

  114. Hernandez Cordero AI, Gonzales NM, Parker CC, Sokolof G, Vandenbergh DJ, Cheng R, et al. Genome-wide associations reveal human-mouse genetic convergence and modifiers of myogenesis, CPNE1 and STC2. Am J Hum Genet. 2019;105(6):1222–36. https://doi.org/10.1016/j.ajhg.2019.10.014.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  115. Middelbos IS, Vester BM, Karr-Lilienthal LK, Schook LB, Swanson KS. Age and diet affect gene expression profile in canine skeletal muscle. PLoS One. 2009;4(2):e4481. https://doi.org/10.1371/journal.pone.0004481.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  116. Stacey D, Chen L, Stanczyk PJ, Howson JMM, Mason AM, Burgess S, et al. Elucidating mechanisms of genetic cross-disease associations at the PROCR vascular disease locus. Nat Commun. 2022;13(1):1222. https://doi.org/10.1038/s41467-022-28729-3.

    Article  ADS  PubMed  PubMed Central  CAS  Google Scholar 

  117. Wang W, Ou Z, Peng J, Zhou Y, Wang N. A transcriptome-wide association study provides new insights into the etiology of osteoarthritis. Ann Transl Med. 2022;10(20):1116. https://doi.org/10.21037/atm-22-4471.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  118. Neville MJ, Wittemans LBL, Pinnick KE, Todorčević M, Kaksonen R, Pietiläinen KH, et al. Regional fat depot masses are influenced by protein-coding gene variants. PLoS One. 2019;14(5):e0217644. https://doi.org/10.1371/journal.pone.0217644.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

None.

Funding

This research received no external funding.

Author information

Authors and Affiliations

  1. Department of Medical Biological Disciplines, Belgorod State National Research University, Belgorod, 308015, Russia

    Vitaly Novakov, Olga Novakova, Maria Churnosova, Inna Aristova, Marina Ponomarenko, Yuliya Reshetnikova, Vladimir Churnosov, Inna Sorokina, Irina Ponomarenko, Olga Efremova, Valentina Orlova, Irina Batlutskaya, Alexey Polonikov, Evgeny Reshetnikov & Mikhail Churnosov

  2. Department of Biology, Medical Genetics and Ecology and Research Institute for Genetic and Molecular Epidemiology, Kursk State Medical University, Kursk, 305041, Russia

    Alexey Polonikov

Authors
  1. Vitaly Novakov
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Olga Novakova
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Maria Churnosova
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Inna Aristova
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Marina Ponomarenko
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yuliya Reshetnikova
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Vladimir Churnosov
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Inna Sorokina
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Irina Ponomarenko
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Olga Efremova
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Valentina Orlova
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Irina Batlutskaya
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Alexey Polonikov
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Evgeny Reshetnikov
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Mikhail Churnosov
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

A.P., O.E., and V.O.: conceptualization; V.N., M.P., I.S. and I.P.: data curation; O.N., I.A., Y.R., E.R. and I.B.: formal analysis; M.C. (Mikhail Churnosov): project administration; V.N. and M.C. (Maria Churnosova): writing—original draft; M.C. (Mikhail Churnosov), V.C., E.R. and I.A.: writing—review and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mikhail Churnosov.

Ethics declarations

Ethics approval and consent to participate

This study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Local Ethical Committee of the Belgorod State University (28 January 2016, No. 4).

Consent for publication

Informed consent was obtained from all subjects involved in this study.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

The literature data about associations of the studied polymorphisms of the candidate genes with ОА.

Additional file 2: Table S2.

The regulatory potential of the studied SNPs.

Additional file 3: Table S3.

Sequence of oligonucleotide primers and probes for SNP candidate genes for KОА.

Additional file 4: Table S4.

The allele and genotype frequencies of the studied SNPs in the KОА and control groups with BMI < 30.

Additional file 5: Table S5.

The allele and genotype frequencies of the studied SNPs in the KОА and control groups with BMI ≥ 30.

Additional file 6: Table S6.

Cis-eQTL values of the KOA-associated polymorphism rs143384 in blood.

Additional file 7: Table S7.

Cis-eQTL values of the SNPs in high LD (r2≥0.80) with KOA-associated polymorphism rs143384 in blood.

Additional file 8: Table S8.

eQTL values of the KOA-associated polymorphism rs143384.

Additional file 9: Table S9.

sQTL values of the KOA-associated polymorphisms rs143384.

Additional file 10: Table S10.

Effect of SNPs in high LD (r2 ≥ 0.80) with the KOA-associated polymorphisms rs143384 on gene expression level.

Additional file 11: Table S11.

Effect of SNPs in high LD (r2 ≥ 0.80) with the KOA-associated polymorphisms rs143384 on alternative splicing (cis-sQTL).

Additional file 12: Table S12.

Results of the gene-gene interaction analysis of the candidate genes associated with rs143384 (functionality in various tissue/organs.

Additional file 13: Table S13.

Results of the gene-gene interaction analysis of the candidate genes associated with rs143384 (functionality in adipose tissue).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Novakov, V., Novakova, O., Churnosova, M. et al. Polymorphism rs143384 GDF5 reduces the risk of knee osteoarthritis development in obese individuals and increases the disease risk in non-obese population. Arthroplasty 6, 12 (2024). https://doi.org/10.1186/s42836-023-00229-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42836-023-00229-9

Keywords

Arthroplasty

ISSN: 2524-7948

Contact us