Sex differences in risk and heritability estimates on primary knee osteoarthritis leading to total knee arthroplasty: a nationwide population based follow up study in Danish twins

Shared genetic factors between osteoarthritis and cardiovascular disease may underlie common etiology

Osteoarthritis is one of the most common musculoskeletal diseases and increases the risk of severe cardiovascular disease, like heart attack and stroke. In some individuals, osteoarthritis and cardiovascular disease will co-occur. This co-occurrence might be due to shared risk factors, for example high age, lifestyle factors and/or a shared genetic liability for the two diseases. Here, we show that the correlation between osteoarthritis and cardiovascular disease can be explained by shared genetic factors, independent of high age and body weight, and also likely independent of lifestyle factors, like smoking and physical activity level. Findings suggest that genetic factors that are shared for osteoarthritis and cardiovascular disease may contribute to both diseases. Thus, the prevailing idea that osteoarthritis is predominantly a risk factor for cardiovascular disease is challenged. Our findings imply that the current diagnostic boundaries between these diseases may need to be re-evaluated.

Phenotype and energy metabolism differ between osteoarthritic chondrocytes from male compared to female patients: Implications for sexual dimorphism in osteoarthritis development?

OBJECTIVES

The prevalence and severity of knee osteoarthritis (OA) are greater in females than males. The purpose of this study was to determine whether there is an underlying difference in the biology of OA chondrocytes between males and females.

METHODS

Chondrocytes were obtained following knee arthroplasty from male and female patients with primary OA. Phenotype marker expression, glucose and fat consumption, and rates of glycolysis and oxidative phosphorylation were compared between females and males. RNAi was used to determine the consequences of differential expression of Sry-box transcription factor 9 (SOX9) and PGC1α between males and females.

RESULTS

OA chondrocytes from male donors showed elevated ribonucleic acid (RNA) and protein levels of SOX9, elevated COL2A1 protein synthesis, higher glucose consumption, and higher usage of glycolysis compared to females. OA chondrocytes from females had higher PGC1α protein levels, higher fat consumption, and higher oxidative energy metabolism than males. Knockdown of SOX9 reduced expression of COL2A1 to a greater extent in male OA chondrocytes than females whereas knockdown of PGC1α reduced COL2A1 expression in females but not males. Expression of ACAN and the glycolytic enzyme PGK1 was also reduced in males but not females following SOX9 knockdown.

CONCLUSIONS

OA chondrocyte phenotype and energy metabolism differ between males and females. Our results indicate transcriptional control of COL2A1 differs between the two. Differences in chondrocyte biology between males and females imply the underlying mechanisms involved in OA may also differ, highlighting the need to consider sex and gender when investigating pathogenesis and potential treatments for OA.

Genetic Influence on Osteoarthritis Versus Other Rheumatic Diseases

OBJECTIVE

We aimed to compare the genetic contribution to osteoarthritis (OA) versus other rheumatic/musculoskeletal diseases (RMDs) in the same population and to explore the role for any shared genetics between OA and other RMDs.

METHODS

In 59,970 Swedish twins aged 35 years or older, we estimated the heritability (the variance explained by genetic factors) of OA in peripheral joints, back and neck pain, shoulder pain (adhesive capsulitis, impingement syndrome, etc), rheumatoid arthritis, spondyloarthritis (SpA) and psoriatic arthritis, myalgia, and osteoporosis diagnosed in specialist and inpatient care. We also studied how much covariance between OA and each of the RMDs could be explained by genetics by studying phenotypic correlations in bivariate classical twin models.

RESULTS

Any-site OA and hip OA (50% and 64%) were among the most heritable RMDs (as compared with 23% for fibromyalgia [lowest] and 63% for SpA [highest]). The highest phenotypic correlations were between OA (any joint site) and shoulder pain in the same individual (r = 0.33, 95% confidence interval 0.31-0.35), of which 70% (95% confidence interval 52-88) could be explained by shared genetics. The phenotypic correlation between OA and back/neck pain was r = 0.25, with 25% to 75% explained by genetics. Phenotypic correlations between OA and each of the other RMDs were lower (r ~ 0.1 to r ~ 0.2), with inconclusive sources of variation.

CONCLUSION

OA has relatively large heritability as compared with other RMDs. The coexistence of OA and shoulder pain, as well as back pain, was common and could often be explained by genetic factors. Findings imply similar etiologies of OA and several pain conditions.

Genetic and environmental influence on resting state networks in young male and female adults: a cartographer mapping study

We propose a unique, minimal assumption, approach based on variance analyses (compared with standard approaches) to investigate genetic influence on individual differences on the functional connectivity of the brain using 65 monozygotic and 65 dizygotic healthy young adult twin pairs' low-frequency oscillation resting state functional Magnetic Resonance Imaging (fMRI) data from the Human Connectome Project. Overall, we found high number of genetically-influenced functional (GIF) connections involving posterior to posterior brain regions (occipital/temporal/parietal) implicated in low-level processes such as vision, perception, motion, categorization, dorsal/ventral stream visuospatial, and long-term memory processes, as well as high number across midline brain regions (cingulate) implicated in attentional processes, and emotional responses to pain. We found low number of GIF connections involving anterior to anterior/posterior brain regions (frontofrontal > frontoparietal, frontotemporal, frontooccipital) implicated in high-level processes such as working memory, reasoning, emotional judgment, language, and action planning. We found very low number of GIF connections involving subcortical/noncortical networks such as basal ganglia, thalamus, brainstem, and cerebellum. In terms of sex-specific individual differences, individual differences in males were more genetically influenced while individual differences in females were more environmentally influenced in terms of the interplay of interactions of Task positive networks (brain regions involved in various task-oriented processes and attending to and interacting with environment), extended Default Mode Network (a central brain hub for various processes such as internal monitoring, rumination, and evaluation of self and others), primary sensorimotor systems (vision, audition, somatosensory, and motor systems), and subcortical/noncortical networks. There were >8.5-19.1 times more GIF connections in males than females. These preliminary (young adult cohort-specific) findings suggest that individual differences in the resting state brain may be more genetically influenced in males and more environmentally influenced in females; furthermore, standard approaches may suggest that it is more substantially nonadditive genetics, rather than additive genetics, which contribute to the differences in sex-specific individual differences based on this young adult (male and female) specific cohort. Finally, considering the preliminary cohort-specific results, based on standard approaches, environmental influences on individual differences may be substantially greater than that of genetics, for either sex, frontally and brain-wide. [Correction added on 10 May 2023, after first online publication: added: functional Magnetic Resonance Imaging. Added: individual differences in, twice. Added statement between furthermore … based on standard approaches.].

Do associations with hand OA vary by knee osteoarthritis phenotype? Cross-sectional data from the Multicenter Osteoarthritis Study

Objective

Osteoarthritis (OA) is highly heterogeneous and has both biomechanical and systemic components that may not have the same etiology. We therefore aimed to identify specific knee OA phenotypes that may be more strongly associated with hand OA to refine the criteria used to define multi-joint OA.

Design

We assessed data from the Multicenter Osteoarthritis Study (MOST). We ascertained hand OA from bilateral hand photographs; scores for each joint row were summed to yield an aggregate hand OA score. Knee OA was ascertained from bilateral posteroanterior knee radiographs read for Kellgren-Lawrence grade and individual radiographic features. We tested associations between hand and knee OA with phenotypes including symptomatic OA, hyper- and atrophic knee OA, and one excluding post-traumatic OA. Associations between hand and knee OA were assessed with logistic regression, adjusted for age.

Results

We studied 2493 participants with hand and knee OA measures. Median age was 63 years with 57% women. 55% had an aggregate hand OA score ≥2; frequency of knee OA phenotypes ranged from 8% to 34%. The age-adjusted odds ratio (OR) was 1.14 (95% confidence interval (CI) ​= ​1.04-1.26) for knee OA per standard deviation of the hand OA aggregate score. Hand OA associations with symptomatic knee OA and knee OA excluding post-traumatic knee OA were OR ​= ​1.16 (95% CI ​= ​1.03-1.31) and OR ​= ​1.21 (95% CI ​= ​1.08-1.35), respectively. No other knee OA phenotype reached statistical significance.

Conclusions

Age-adjusted associations between hand and knee OA were modest and were largely similar across knee OA phenotypes.

Does lower educational attainment increase the risk of osteoarthritis surgery? a Swedish twin study

BACKGROUND

Previous studies have reported an inverse association between educational attainment and different osteoarthritis (OA) outcomes. However, none of the previous studies have accounted for potential confounding by early-life environment and genetics. Thus, we aimed to examine the association between educational attainment and knee and hip OA surgery using twin data.

METHODS

From the Swedish Twin Registry (STR), we identified dizygotic (DZ) and monozygotic (MZ) twins. All twins in the STR aged 35 to 64 years were followed from January the 1^st 1987 or the date they turned 35 years until OA surgery, relocation outside Sweden, death or the end of 2016 (18,784 DZ and 8,657 MZ complete twin pairs). Associations between educational attainment and knee and hip OA surgery were estimated in models matched on twin pairs, using Weibull within-between (WB) shared frailty model.

RESULTS

For knee OA surgery, the analysis matched on MZ twins yielded a within-estimate hazard ratio (HR) per 3 years of education, of 1.06 (95% CI: 0.81, 1.32), suggesting no association between the outcome and the individual´s education. Rather, there seemed to be a so called familial effect of education, with a between-pair estimate of HR = 0.71 (95% CI: 0.41, 1.01). For hip OA surgery, the within- and between-pair estimates for MZ twins were 0.92 (95% CI: 0.69, 1.14) and 1.15 (95% CI: 0.87, 1.42), respectively.

CONCLUSIONS

Our results suggest that the inverse associations between education and knee/hip OA surgery observed in cohort studies are potentially confounded by unobserved familial factors like genetics and/or early life exposures.

Effects of medications on incidence and risk of knee and hip joint replacement in patients with osteoarthritis: a systematic review and meta-analysis

Background

This systematic review and meta-analysis aimed to investigate the incidence and risk of knee and hip replacement in patients with osteoarthritis (OA) treated with different medications.

Methods

OVID MEDLINE, OVID EMBASE, Cochrane and Web of Science electronic databases were searched from inception to May 4th, 2022. Clinical trials, including randomized controlled trials, cohort studies and case–control studies, were selected. The meta-analysis effect size was estimated using either incidence with 95% confidence intervals (CIs) or odds ratio (OR)/relative risk (RR) with 95% CIs. The risk of bias and heterogeneity among studies were assessed and analyzed.

Results

Forty studies were included, involving 6,041,254 participants. The incidence of joint replacement in patients with OA varied according to the study design and treatments. The incidence of knee arthroplasty varied from 0 to 70.88%, while the incidence of hip arthroplasty varied from 11.71 to 96.43%. Compared to non-users, bisphosphonate users had a reduced risk of knee replacement (RR = 0.71, 95% CI: 0.66–0.77; adjusted hazard ratio [aHR] = 0.76, 95% CI: 0.70–0.83). Compared to intra-articular corticosteroid users, hyaluronic acid (HA) users had a higher risk of knee arthroplasty (RR = 1.76, 95% CI: 1.38–2.25). No publication bias was observed.

Conclusions

Bisphosphonate treatment is associated with a reduced risk of knee replacement. More studies are needed to validate our results due to the limited number of eligible studies and high heterogeneity among studies.

Characterization of Genetic Risk of End-Stage Knee Osteoarthritis Treated with Total Knee Arthroplasty: A Genome-Wide Association Study

BACKGROUND

End-stage knee osteoarthritis (OA) is a highly debilitating disease for which total knee arthroplasty (TKA) serves as an effective treatment option. Although a genetic component to OA in general has been described, evaluation of the genetic contribution to end-stage OA of the knee is limited. To this end, we present a genome-wide association study involving patients undergoing TKA for primary knee OA to characterize the genetic features of severe disease on a population level.

METHODS

Individuals with the diagnosis of knee OA who underwent primary TKA were identified in the U.K. Biobank using administrative codes. The U.K. Biobank is a data repository containing prospectively collected clinical and genomic data for >500,000 patients. A genome-wide association analysis was performed using the REGENIE software package. Logistic regression was also used to compare the total genetic risk between subgroups stratified by age and body mass index (BMI).

RESULTS

A total of 16,032 patients with end-stage knee OA who underwent primary TKA were identified. Seven genetic loci were found to be significantly associated with end-stage knee OA. The odds ratio (OR) for developing end-stage knee OA attributable to genetics was 1.12 (95% confidence interval [CI], 1.10 to 1.14), which was lower than the OR associated with BMI (OR = 1.81; 95% CI, 1.78 to 1.83) and age (OR = 2.38; 95% CI, 2.32 to 2.45). The magnitude of the OR for developing end-stage knee OA attributable to genetics was greater in patients <60 years old than in patients ≥60 years old (p = 0.002).

CONCLUSIONS

This population-level genome-wide association study of end-stage knee OA treated with primary TKA was notable for identifying multiple significant genetic variants. These loci involve genes responsible for cartilage development, cartilage homeostasis, cell signaling, and metabolism. Age and BMI appear to have a greater impact on the risk of developing end-stage disease compared with genetic factors. The genetic contribution to the development of severe disease is greater in younger patients.

LEVEL OF EVIDENCE

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

The genetic contribution to hand osteoarthritis

OBJECTIVE

To estimate the genetic contribution to doctor-diagnosed hand osteoarthritis (OA).

METHODS

Using data from the Swedish Twin Registry and National Patient Register, we conducted a 20-year population-based longitudinal cohort study including 59,970 twins aged 35 years or older. We studied inpatient and outpatient doctor-diagnosed hand OA using ICD-10 codes from 1997 until 2016, including both the distal/proximal interphalangeal (DIP/PIP) joints and/or the first carpometacarpal (CMC-1) joints. We calculated intra-pair correlation, estimated the heritability (i.e., the percentage variation in hand OA that can be explained by genetic factors) as well as a genetic risk.

RESULTS

Among 59,970 included persons, 936 had a hand OA diagnosis registered during the study period. The heritabilities of hand OA (any joint), CMC-1 OA and DIP/PIP OA were ∼87%, 86% and 48%, respectively, yet the two latter should be interpreted with care due to low numbers. Hand OA in any joint in both twins in a pair occurred more frequently in identical twins (54/554 = 9.7%, intra-pair correlation = 0.54, 95% CI = 0.44-0.63) than in fraternal twins (18/1,246 = 1.4%, intra-pair correlation = 0.10, 95% CI = -0.01-0.22). Identical twins who were diagnosed with hand OA in any joint had a far higher risk than fraternal twins with hand OA to also have their co-twin diagnosed with hand OA in any joint (Hazard Ratio = 6.98, 95% CI = 3.08-15.45).

CONCLUSION

The genetic contribution to hand OA is high and likely varying between 48% and 87%. Potential differential heritability by hand OA phenotypes should be further explored.

Body Mass Index Differentially Moderates Heritability of Total Joint Replacement Due to Hip and Knee Osteoarthritis: A Cohort Study of 29,893 Swedish Twin Pairs

BACKGROUND

Osteoarthritis and obesity are diseases with high prevalence, and they share common etiologies. We investigated the sex-specific genetic susceptibility to hip and knee osteoarthritis necessitating total joint replacement (TJR), and how body mass index (BMI) moderated the heritability of these osteoarthritis phenotypes.

METHODS

We linked 29,893 twin pairs with information on BMI in the Swedish Twin Registry with the Swedish National Patient Register to identify twins who underwent primary TJR of the hip or knee combined with a concomitant diagnosis of primary osteoarthritis of these joints. Structural equation modeling was used to calculate the heritability of hip and knee osteoarthritis treated with TJR, with estimates adjusted for the first available BMI, birth year, and sex. We also investigated how heritability varied with BMI treated as a continuous variable.

RESULTS

Similar heritability estimates for hip replacement (0.65 [95% confidence interval (CI), 0.59 to 0.70]) and knee replacement (0.57 [95% CI, 0.50 to 0.64]) were found. Heritability decreased with higher BMI in both sexes for hip replacement and in men for knee replacement. In contrast, heritability for knee replacement increased with higher BMI in women; the estimate was 0.37 (90% likelihood interval [LI], 0.25 to 0.49) for a BMI of 20 kg/m2 and 0.87 (90% LI, 0.68 to 0.94) for a BMI of 35 kg/m2.

CONCLUSIONS

In our population, heritability explained, on average, about half of the susceptibility to undergo primary TJR of the hip or knee with the indication of primary osteoarthritis, but it varied with BMI and sex. We demonstrated substantial heritability for knee replacement in obese women.

LEVEL OF EVIDENCE

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

Current perspectives on the role of biomechanical loading and genetics in development of disc degeneration and low back pain; a narrative review

Degenerative changes in the disc have long been of interest; they are thought to be strongly associated with low back pain and caused by inappropriate loading or through injury. However, independent of the magnitude of occupational spinal loading, twin studies find that the heritability of lumbar disc degeneration is 34-74%. This finding has led to intensive searches for susceptibility genes; some genes associated with disc degeneration have been identified, though all with small effects on the degenerative process. The complex nature of degenerative changes suggests that many different genes are involved, and that interactions with environmental factors are influential in progression of degeneration. Low back pain itself also appears heritable (30-46%). The most important clinical question though, is not how discs degenerate but is disc degeneration related to low back pain. Imaging studies find many people with degenerate discs or even with discs showing pathological features such as herniations, are asymptomatic. However results are obscured by the lack of consistent definitions of the phenotypes of disc degeneration and of low back pain. Epidemiological studies could help disentangle these complex relationships, but they will only be successful once consistent classifications and phenotypes of both disc degeneration and low back pain are developed.

Nature vs nurture in knee osteoarthritis - the importance of age, sex and body mass index

OBJECTIVE

(1) To estimate the life-time genetic contribution for knee osteoarthritis (OA) surgery and (2) to explore any differences in the genetic contribution across age, sex and body mass index (BMI).

METHODS

We studied the sex-specific genetic contribution to knee OA surgery in a prospective cohort study of 62,490 twins aged 35 years or older with a follow-up period of up to 47 years (10,092 identical and 21,153 non-identical twin pairs, 54% women). To study interactions with age, we graphed the heritabilities over the lifespan for men and women. We also studied the sex-specific heritability across strata of the median BMI to explore any interactions with BMI.

RESULTS

The overall heritability of knee OA surgery was 0.53 (95% confidence intervals [CI] = 0.31-0.75), with higher heritability among women (H² = 0.80 (95% CI = 0.73-0.87)) than men (H² = 0.39 (95% CI = 0.10-0.69)). For men, the heritability started to rise after age 68. The genetic contribution was particularly low in men above median BMI (H²_{≥23.7 kg/m2} = 0.08, 95% CI = -0.32-0.48). For women, the heritability was consistently high from age 50 to death, independently of BMI (H²_{≥22.5 kg/m2} = 0.77, 95% CI = 0.66-0.87).

CONCLUSION

There is a higher and more consistent genetic contribution for knee OA surgery in women than men. In men the genetic contribution was relatively low and varied with age and BMI.

Bone shape mediates the relationship between sex and incident knee osteoarthritis

Background

Knee bone shape differs between men and women and the incidence of knee osteoarthritis (OA) is higher in women than in men. Therefore, the purpose of the present study was to determine whether the observed difference in the incidence of knee radiographic OA (ROA) between men and women is mediated by bone shape.

Methods

We randomly sampled 304 knees from the OAI with incident ROA (i.e., development of Kellgren/Lawrence grade ≥ 2 by month 48) and 304 knees without incident ROA. We characterized distal femur and proximal tibia shape on baseline radiographs using Statistical Shape Modeling. If a specific bone shape was associated with the risk of incident ROA, marginal structural models were generated to assess the mediation effect of that bone shape on the relation of sex and risk of incident knee ROA adjusting for baseline covariates.

Results

Case and control participants were similar by age, sex and race, but case knees were from higher body mass index (BMI) participants (29.4 vs. 27.0; p < 0.001). Women had 49% increased odds of incident knee ROA compared with men (adjusted odds ratio (OR) = 1.49, 95% Confidence Interval (C.I.): 1.04, 2.12). There was an inconsistent mediation effect for tibial mode 2 between sex and incident knee ROA, with an indirect effect OR of 0.96 (95% C.I.: 0.91–1.00) and a direct effect OR of 1.56 (95% C.I.: 1.08–2.27), suggesting a protective effect for this mode. Similar findings were also observed for the mediation effect of tibia mode 10 and femur mode 4. These shape modes primarily involved differences in the angular relation of the heads to the shafts of the femur and tibia.

Conclusions

Distal femur and proximal tibia bone shapes partially and inconsistently mediated the relationship between sex and incident knee OA. Women had a higher risk of incident ROA, and specific bone shapes modestly protected them from even higher risk of ROA. The clinical significance of these findings warrant further investigation.

Twin studies for the prognosis, prevention and treatment of musculoskeletal conditions

BACKGROUND

Musculoskeletal conditions are highly prevalent in our ageing society and are therefore incurring substantial increases in population levels of years lived with disability (YLD). An evidence-based approach to the prognosis, prevention, and treatment of those disorders can allow an overall improvement in the quality of life of patients, while also softening the burden on national health care systems.

METHODS

In this Masterclass article, we provide an overview of the most relevant twin study designs, their advantages, limitations and major contributions to the investigation of traits related to the domain of musculoskeletal physical therapy.

CONCLUSIONS

Twin studies can be an important scientific tool to address issues related to musculoskeletal conditions. They allow researchers to understand how genes and environment combine to influence human health and disease. Twin registries and international collaboration through existing networks can provide resources for achieving large sample sizes and access to expertise in study design and analysis of twin data.

Genetic factors contribute more to hip than knee surgery due to osteoarthritis - a population-based twin registry study of joint arthroplasty

OBJECTIVE

To explore and quantify the relative strengths of the genetic contribution vs the contribution of modifiable environmental factors to severe osteoarthritis (OA) having progressed to total joint arthroplasty.

DESIGN

Incident data from the Norwegian Arthroplasty Registry were linked with the Norwegian Twin Registry on the National ID-number in 2014 in a population-based prospective cohort study of same-sex twins born 1915-60 (53.4% females). Education level and height/weight were self-reported and Body Mass Index (BMI) calculated. The total follow-up time was 27 years for hip arthroplasty (1987-2014, 424,914 person-years) and 20 years for knee arthroplasty (1994-2014, 306,207 person-years). We estimated concordances and the genetic contribution to arthroplasty due to OA in separate analyses for the hip and knee joint.

RESULTS

The population comprised N = 9058 twin pairs (N = 3803 monozygotic (MZ), N = 5226 dizygotic (DZ)). In total, 73% (95% confidence intervals (CI) = 66-78%) and 45% (95% CI = 30-58%) of the respective variation in hip and knee arthroplasty could be explained by genetic factors. Zygosity (as a proxy for genetic factors) was associated with hip arthroplasty concordance over time when adjusted for sex, age, education and BMI (HR = 2.98, 95% CI = 1.90-4.67 for MZ compared to DZ twins). Knee arthroplasty was to a greater extent dependent on BMI when adjusted for zygosity and the other covariates (HR = 1.15, 95% CI = 1.02-1.29).

CONCLUSION

Hip arthroplasty was strongly influenced by genetic factors whereas knee arthroplasty to a greater extent depended on a high BMI. The study may imply there is a greater potential for preventing progression of knee OA to arthroplasty in comparison with hip OA.