Association study of candidate genes for the prevalence and progression of knee osteoarthritis

A novel biomarker of MMP-cleaved cartilage intermediate layer protein-1 is elevated in patients with rheumatoid arthritis, ankylosing spondylitis and osteoarthritis

Rheumatic joints have an altered cartilage turnover. Cartilage intermediate layer protein 1 (CILP-1) is secreted from articular chondrocytes and deposited into the cartilage extracellular matrix. We developed an immunoassay targeting a Matrix Metalloproteinase (MMP)-generated neo-epitope of CILP-1, named CILP-M. Human articular cartilage was cleaved with proteolytic enzymes and CILP-M levels were measured. We also quantified CILP-M in two studies from patients with rheumatoid arthritis (RA), ankylosing spondylitis (AS) and osteoarthritis (OA) and explored the monitoring and prognostic potential of CILP-M in TNF-α inhibitory treatment and modified Stoke AS Spine Score (mSASSS) progression. CILP-M was generated by MMP-1, -8 and -12. In the discovery study, CILP-M was significantly higher in patients with RA, AS and OA than healthy donors (p < 0.01, p < 0.001, p < 0.05) with an area under the curve (AUC) between the diseased groups and healthy donors > 0.95 (p < 0.001). In the validation study, patients with RA and AS had significantly higher CILP-M levels than healthy controls (p < 0.001) and AUC > 0.90 (p < 0.001). Patients with AS treated with TNF- α inhibitory treatment in the validation study had significantly lower CILP-M levels after treatment (p = 0.004). CILP-M may provide useful insights into cartilage degradation processes in rheumatic diseases.

Mechanotransduction pathways in articular chondrocytes and the emerging role of estrogen receptor-α

In the synovial joint, mechanical force creates an important signal that influences chondrocyte behavior. The conversion of mechanical signals into biochemical cues relies on different elements in mechanotransduction pathways and culminates in changes in chondrocyte phenotype and extracellular matrix composition/structure. Recently, several mechanosensors, the first responders to mechanical force, have been discovered. However, we still have limited knowledge about the downstream molecules that enact alterations in the gene expression profile during mechanotransduction signaling. Recently, estrogen receptor α (ERα) has been shown to modulate the chondrocyte response to mechanical loading through a ligand-independent mechanism, in line with previous research showing that ERα exerts important mechanotransduction effects on other cell types, such as osteoblasts. In consideration of these recent discoveries, the goal of this review is to position ERα into the mechanotransduction pathways known to date. Specifically, we first summarize our most recent understanding of the mechanotransduction pathways in chondrocytes on the basis of three categories of actors, namely mechanosensors, mechanotransducers, and mechanoimpactors. Then, the specific roles played by ERα in mediating the chondrocyte response to mechanical loading are discussed, and the potential interactions of ERα with other molecules in mechanotransduction pathways are explored. Finally, we propose several future research directions that may advance our understanding of the roles played by ERα in mediating biomechanical cues under physiological and pathological conditions.

Expression of CILP-2 and DDR2 and ultrastructural changes in the articular cartilage of patients with knee osteoarthritis undergoing total knee arthroplasty: a pilot morphological study

The aim of the study was to correlate the immunohistochemical expression of cartilage intermediate layer protein 2 (CILP-2) and discoidin domain receptor 2 (DDR2), and the ultrastructural changes in the cartilage with the degree of articular cartilage damage in osteoarthritis (OA) patients. Cartilage samples were obtained from twenty patients aged from 46 to 68 years undergoing total knee arthroplasty. In each patient, medial and lateral tibial plateau samples were analysed applying OARSI histopathology grading. Positive correlation was noted between the extent of CILP-2 staining intensity and OARSI grades. Abundant staining for CILP-2 was found in the superficial and middle layers and in the pericellular matrix (PCM) of the deep zone. Transmission electron microscopy studies demonstrated strong damage of chondrocytes, the organelles were often diminished or focally aggregated. As a characteristic finding, PCM was frequently expanded, which may reflect a pathogenic step in OA progression. In conclusion, CILP-2 may potentially be a relevant marker of OA progression as its expression correlated better with cartilage damage than the known marker of articular cartilage damage, DDR2.

The association between ADAM12 gene polymorphisms and osteoarthritis: an updated meta-analysis

BACKGROUND

Osteoarthritis of the knee is an irreversible disease that causes great pain, and genetic factors play an important role in its occurrence and development. There have been many studies on the correlation between ADAM12 polymorphisms and genetic susceptibility to osteoarthritis, but the results remain inconclusive.

METHODS

Papers from PubMed, Web of Science, EMbase, Springer, SCOPUS, Google Scholar and other databases were systematically retrieved with a cut-off of January 2022. All case-control studies on ADAM12 rs3740199, rs1871054, rs1044122, and rs1278279 polymorphisms and osteoarthritis were searched. Fixed or random effects models were used for pooled analysis with OR values and 95% confidence intervals (CI), and publication bias was assessed. In addition, the false-positive reporting probability test was used to assess the confidence of a statistically significant association.

RESULTS

Eleven articles were included, which included 3332 patients with osteoarthritis and 5108 healthy controls. Meta-analysis showed that the rs1871054 polymorphism of ADAM12 was associated with osteoarthritis in dominant, recessive, allelic, and homozygote genetic models [C vs. T: OR = 1.34 95% CI (1.05, 1.71), P < 0.001]. Our subgroup analysis revealed an association between the ADAM12 polymorphism rs1871054 in Asians and osteoarthritis [C vs. T: OR = 1.61, 95% CI (1.25, 2.08), P < 0.001], albeit this was only for three studies. In addition, the ADAM12 polymorphism rs1871054 is associated with osteoarthritis in patients younger than 60 years of age [C vs. T: OR = 1.39, 95% CI (1.01, 1.92), P = 0.289]; however, the ADAM12 gene rs3740199, rs1044122, and rs1278279 site polymorphisms were not significantly. Furthermore, when assessing the confidence of the positive results, the positive results were found to be credible (except for Age < 60).

CONCLUSION

Polymorphism at the rs1871054 site of ADAM12 is associated with genetic susceptibility to osteoarthritis, but rs3740199, rs1044122, and rs1278279 site polymorphisms are not.

Polymorphism of LYPLAL1 and TGFA Genes Associated With Progression of Knee Osteoarthritis in Residents Central Chernozem Region of Russia

Background. Кnee osteoarthritis (OA) is a multifactorial disease in which genetic factors play an important role. The share of the hereditary component in the development of OA, according to various literature sources, ranges from 40 to 65%. Кnee OA is a progressive disease that leads to a decrease in the quality of life and disability. The study aimed to evaluate the role of polymorphic markers of candidate genes rs2820436 and rs2820443 LYPLAL1, rs3771501 TGFA, rs11177 GNL3, rs6976 GLT8D1 in the progression of knee OA in the population of the Central Chernozem Region of Russia. Methods. The study was performed in a case-control design on a sample of 500 patients with knee OA. Case patients with III-IV stages of the disease according to KellgrenLawrence (n = 325), control (individuals who do not have the analyzed sign III-IV stages of the disease) patients with stage II (n = 175). Genotyping of five single nucleotide polymorphisms (SNPs) of candidate genes was performed using the polymerase chain reaction method for DNA synthesis. The study of the associations of the studied polymorphic loci, the calculation of haplotype frequencies and the analysis of their relationship with the progression of knee OA was carried out by the method of logistic regression in the program PLINK v 2.050. Results. Significant associations with the progression of OA of the knee were established for allelic variant A rs2820436 of LYPLAL1 gene according to allelic (OR = 1.48, p = 0.010, pperm = 0.012), additive (OR = 1.58, p = 0.009, pperm = 0.010), dominant (OR = 1.61, p = 0.024, pperm = 0.030) genetic models and A/A genotype of the same polymorphism (OR = 2.53, p = 0.041). The genotypes C/C rs2820436 LYPLAL1 (OR = 0.67, p = 0.043), A/G rs3771501 TGFA (OR = 0.67, p = 0.042) have a protective role in the progression of the disease. It was found that the frequency of the AC haplotype of haploblock rs2820436-rs2820443 in the group of patients with III-IV stages of the disease was significantly higher than in patients with stage II (OR = 1.83, p = 0.002, pperm = 0.002). The identified molecular genetic markers rs2820436 and rs2820443 of LYPLAL1 gene, rs3771501 of TGFA gene are associated both with the risk of developing OA according to previous genome-wide studies and, according to our data, are associated with the progression of knee OA. Conclusions. Genetic risk factors for the development of knee OA of III-IV radiological stages are allelic variant A and genotype A/A rs2820436 of LYPLAL1 gene, haplotype AC of haploblock rs2820436-rs2820443 in the population of the Central Chernozem Region of Russia. Genotypes C/C rs2820436 of LYPLAL1 gene and A/G rs3771501 of TGFA gene have a protective value in the progression of this disease.

Functional biomaterials for osteoarthritis treatment: From research to application

Osteoarthritis (OA) is a common disease that endangers millions of middle-aged and elderly people worldwide. Researchers from different fields have made great efforts and achieved remarkable progress in the pathogenesis and treatment of OA. However, there is still no cure for OA. In this review, we discuss the pathogenesis of OA and summarize the current clinical therapies. Moreover, we introduce various natural and synthetic biomaterials for drug release, cartilage transplantation, and joint lubricant during the OA treatment. We also present our perspectives and insights on OA treatment in the future. We hope that this review will foster communication and collaboration among biological, clinical, and biomaterial researchers, paving the way for OA therapeutic breakthroughs.

Association between CILP and IL-1α polymorphisms and phenotype-dependent intervertebral disc degeneration susceptibility: A meta-analysis

Background: The relationship between CILP (1184T>C) and IL-1α(+889C/T) polymorphisms and intervertebral disc degeneration (IDD) have been explored in several studies but the results were conflicting. The aim of the study was to evaluate and synthesize the currently available data on the association between CILP (1184T>C) and IL-1α(+889C/T) polymorphisms and susceptibility of phenotype-dependent radiologic IDD (RIDD) and symptomatic intervertebral disk herniation (SIDH). Methods: A computerized literature search was in PubMed, Cochrane Library, Embase, China National Knowledge Infrastructure database, and Web of Science. The pooled results were presented as odds ratios (ORs) with 95% confidence intervals (CIs). Moreover, the false-positive report probability (FPRP) test and trial sequential analysis (TSA) were applied to estimate the significant results. Results: Our evidence demonstrated that IL-1α(+889C/T) was significant associated with RIDD (allele model: OR = 1.34, 95%CI 1.03-1.74, p = 0.029) and SIDH (allele model: OR = 1.28, 95% CI 1.03-1.60, p = 0.028). However, the results were not noteworthy under the FPRP test and TSA analysis. Additionally, CILP (1184T>C) polymorphism was significantly associated with RIDD with adequate evidence (allele model: OR = 1.27, 95% CI 1.09-1.48, p = 0.002) instead of SIDH. Conclusion: The current meta-analysis illustrated firm evidence that CILP (1184T>C) polymorphism was significantly associated with the susceptibility of RIDD. However, the significant associations between IL-1α(+889C/T) and RIDD and SIDH were less credible. Thus, more multi-center studies with diverse populations were required to verify the results.

Serum Concentrations of Cartilage Intermediate Layer Protein 2 Were Higher in Overweight and Obese Subjects

Background

Cartilage intermediate layer protein 2 (CILP2) is associated with a variety of plasma lipoproteins and lipid traits. However, the correlation between CILP2 and obesity remains unknown. The aim of this study was to investigate the relationship between circulating CILP2 levels and obesity based on body mass index (BMI).

Methods

A total of 252 subjects were divided into three groups: normal weight (n = 124), overweight (n = 94), and obese (n = 34). Metabolic parameters were measured in a fasting state. Serum CILP2 concentration was tested by enzyme-linked immunosorbent assay. Multivariate linear regression analysis was used to explore the relationship between CILP2 and obesity. We also conducted bioinformatics analysis to further explore the genes and signaling pathways related to CILP2.

Results

The concentrations of serum CILP2 in the overweight and obese groups were significantly higher than that in the normal weight group. In multiple linear regression analysis, BMI was positively correlated with CILP2 concentration after controlling gender and age. Being overweight and obese were independently correlated with CILP2 concentration after adjusting for gender, age, SBP, DBP, FBG, 2-hour OGTT blood glucose (2h-BG), fasting blood insulin (FIns), TG, TC, HDL-C, LDL-C, and FFA. Bioinformatics analysis showed that the genes related to CILP2 are primarily associated with lipid metabolism and insulin resistance.

Conclusion

We speculate that CILP2 may attribute to metabolic disorders in obesity.

Protein N-glycosylation aberrations and glycoproteomic network alterations in osteoarthritis and osteoarthritis with type 2 diabetes

Whether the relationship between type 2 diabetes mellitus (T2DM) and osteoarthritis (OA) can be solely attributed to the shared risk factors, such as obesity, remains controversial. Several studies have revealed the critical role of abnormal glycosylation in the pathogenesis of OA and T2DM. Therefore, we speculate that T2DM may contribute to the pathogenesis of OA through the intrinsic mechanisms of N-glycosylation aberrations. Using N-glycoproteomics, we compared the changes in N-glycosylated protein abundance in cartilage samples from patients with OA without and with T2DM (DM-OA), and from patients with traumatic joint injury (NC) as controls. We identified 847 N-glycosylation sites corresponding to 729 peptides fragments from 374 proteins. The number of N-glycosylated proteins in the DM-OA group tended to decrease compared with that in the OA and NC groups. We identified 22 upregulated and 1 down-regulated N-glycosylated peptides in the OA group compared to the NC group, while only fibronectin 1 (FN1) at position N1007, cartilage intermediate layer protein 1 (CILP) at N346, and collagen type VI alpha 1 chain (COL6A1) at N804, were also identified in the DM-OA group. Compared to the OA group, the downregulation of secreted protein acidic and rich in cysteine (SPARC) at N116, collagen type VI alpha 1 chain (COL6A2) at N785, and asporin (ASPN) at N282, and the upregulation of complement component C8 alpha chain (C8α) at N437, were the most remarkable alterations in the DM-OA group. The differentially expressed N-glycosylated proteins between the OA and DM-OA groups were mainly located extracellularly and enriched in the KEGG pathways involving PI3K/Akt signaling, focal adhesion, and ECM-receptor interaction. Their predicted protein–protein interactions were also depicted. We were thus able to show the general characteristics of N-glycosylation aberrations in OA and DM-OA. Moreover, the upregulated glycosylated complement C8α in the DM-OA group might augment membrane attack complex activity, thereby exacerbating cartilage destruction. Although further confirmation is required, our hypothesis proposes a possible explanation for the deduction that T2DM is an independent risk factor for OA.

Interleukin-1 genetic polymorphisms in knee osteoarthritis: What do we know? A meta-analysis and systematic review

PURPOSE

Interleukin-1 is the main proinflammatory cytokine in osteoarthritis (OA). Several single-nucleotide polymorphisms (SNPs) within the IL-1 gene cluster (IL-1β, IL-1R1, and IL-1RN) have been determined, but their associations with knee OA remain poorly understood. The present study aimed to identify the associations between IL-1 SNPs and knee osteoarthritis.

METHODS

This meta-analysis and systematic review included all comparative studies published in the MEDLINE/PubMed, Embase, Google Scholar, and Cochrane Library databases. We performed a systematic search to identify relevant studies on the evaluation of the correlation between the IL-1 gene and knee OA published up to February 2020 that met the eligibility criteria. Nine studies on a total of 2256 knees with OA and 3527 healthy knees met the eligibility criteria. Results associated with IL-1A, IL-1B, IL-1R1, and IL-1RN SNPs were extracted and compared between knees with OA and healthy knees. Methodological quality was assessed using the Newcastle-Ottawa scale (NOS). All studies with fair or good quality were included.

RESULTS

The meta-analysis showed that the risk of knee OA is decreased by the IL-1RN*1 and IL-1RN*1/*1 genotypes and increased by the IL-1RN*2 and I-L1RN*1/*2 genotypes. The systematic review revealed only two studies associating the IL-1RN allele, none associating the IL-1B polymorphism, and only one study associating IL-1A and IL-1R1 polymorphisms with knee OA.

CONCLUSIONS

Several IL-1RN alleles and genotypes play a role in knee OA but other genetic variations in the IL-1 region were still conflicting in its association with knee OA.

Osteoarthritis: Insights Offered by the Study of Bone Mass Genetics

PURPOSE OF REVIEW

This paper reviews how bone genetics has contributed to our understanding of the pathogenesis of osteoarthritis. As well as identifying specific genetic mechanisms involved in osteoporosis which also contribute to osteoarthritis, we review whether bone mineral density (BMD) plays a causal role in OA development.

RECENT FINDINGS

We examined whether those genetically predisposed to elevated BMD are at increased risk of developing OA, using our high bone mass (HBM) cohort. HBM individuals were found to have a greater prevalence of OA compared with family controls and greater development of radiographic features of OA over 8 years, with predominantly osteophytic OA. Initial Mendelian randomisation analysis provided additional support for a causal effect of increased BMD on increased OA risk. In contrast, more recent investigation estimates this relationship to be bi-directional. However, both these findings could be explained instead by shared biological pathways. Pathways which contribute to BMD appear to play an important role in OA development, likely reflecting shared common mechanisms as opposed to a causal effect of raised BMD on OA. Studies in HBM individuals suggest this reflects an important role of mechanisms involved in bone formation in OA development; however further work is required to establish whether the same applies to more common forms of OA within the general population.

CILP2 overexpression correlates with tumor progression and poor prognosis in patients with colorectal cancer in The Cancer Genome Atlas (TCGA) study

Background

Genetic alterations play an important role in the progression of colorectal cancer (CRC). Identifying new biomarkers to assess the prognosis of patients with CRC is critical. Cartilage intermediate layer protein 2 (CILP2) gene, screened from TCGA database by bioinformatics, may be closely related to the progression of CRC. CILP2 was barely reported with clinical features of tumors.

Materials and methods

Clinical information and RNA-seq data were derived from TCGA colorectal carcinoma cohort. CILP2 expression at mRNA level was estimated by bioinformatical analysis of TCGA cases. Tissue microarray (TMA) was constructed containing paraffin-embedded 64 pairs of CRC and matched adjacent normal tissues. The expression at the protein level was detected in 64 pairs of CRC and matched adjacent normal tissues by immunohistochemical analysis. CILP2 expression level and its clinical value were estimated by bioinformatical analysis with linear and logistic regression. Survival analysis was performed between high and low groups of CILP2 expression by Cox regression analysis, and the P value was calculated by the log-rank test. The Kaplan-Meier curves were tested by the log-rank test.

Results

CILP2 was statistically significantly higher expressed in the CRC tissues when compared with paired adjacent normal tissues in TCGA cohort (P < 0.001) and in the TMA cohort (P = 0.001). Also, CILP2 high expression was strongly correlated with T3/4 stage (P = 0.001), N1/2/3 stage (P = 0.005), M1 stage (P = 0.048), and higher clinical stage (UICC 2010 stage) (P < 0.001) in TCGA cohort, and also positively associated with T3/4 stage (P = 0.022) and higher clinical stage (UICC 2010 stage) (P = 0.03) in TMA cohort. Furthermore, CILP2 overexpression predicted poor prognosis and could be an independent prognostic factor (P = 0.003).

Conclusion

We revealed that CILP2 is associated with advanced stages and could play a role as an independent predictor of poor survival in CRC.

Identification of TGFβ signatures in six murine models mimicking different osteoarthritis clinical phenotypes

OBJECTIVE

TGFβ is a key player in cartilage homeostasis and OA pathology. However, few data are available on the role of TGFβ signalling in the different OA phenotypes. Here, we analysed the TGFβ pathway by transcriptomic analysis in six mouse models of OA.

METHOD

We have brought together seven expert laboratories in OA pathophysiology and, used inter-laboratories standard operating procedures and quality controls to increase experimental reproducibility and decrease bias. As none of the available OA models covers the complexity and heterogeneity of the human disease, we used six different murine models of knee OA: from post-traumatic/mechanical models (meniscectomy (MNX), MNX and hypergravity (HG-MNX), MNX and high fat diet (HF-MNX), MNX and seipin knock-out (SP-MNX)) to aging-related OA and inflammatory OA (collagenase-induced OA (CIOA)). Four controls (MNX-sham, young, SP-sham, CIOA-sham) were added. OsteoArthritis Research Society International (OARSI)-based scoring of femoral condyles and ribonucleic acid (RNA) extraction from tibial plateau samples were done by single operators as well as the transcriptomic analysis of the TGFβ family pathway by Custom TaqMan® Array Microfluidic Cards.

RESULTS

The transcriptomic analysis revealed specific gene signatures in each of the six models; however, no gene was deregulated in all six OA models. Of interest, we found that the combinatorial Gdf5-Cd36-Ltbp4 signature might discriminate distinct subgroups of OA: Cd36 upregulation is a hallmark of MNX-related OA while Gdf5 and Ltbp4 upregulation is related to MNX-induced OA and CIOA.

CONCLUSION

These findings stress the OA animal model heterogeneity and the need of caution when extrapolating results from one model to another.

MFN2 contributes to metabolic disorders and inflammation in the aging of rat chondrocytes and osteoarthritis

OBJECTIVE

Metabolic disorders and inflammation of chondrocytes are major pathological changes in aging cells and osteoarthritis (OA). Recent studies demonstrated age-related mitochondrial dysfunction may be a key contributing factor in the development of OA. Mitofusin 2 (MFN2) is a key regulator of mitochondrial fusion, cell metabolism, autophagy and apoptosis. This study was performed to ascertain whether MFN2 was involved in the aging of chondrocytes and OA.

METHODS

Metabolic measurements were taken in rat chondrocytes between different ages (3-week, 5-month, 12-month). MFN2 activity was detected in both human and rat chondrocytes during aging and OA. Then, knockdown of MFN2 with small interfering RNA (siRNA) was performed to confirm whether MFN2 contributes to metabolic changes. Lentiviruses were used to establish MFN2-overexpression/knockdown OA models both in vivo and in vitro to confirm whether MFN2 contributes to OA progress. Further, regulatory mechanism of MFN2 was assessed and interaction between MFN2 and PARKIN was performed.

RESULTS

A metabolic shift to mitochondrial respiration was confirmed in rat chondrocytes during aging. MFN2 expression was elevated in both human and rat chondrocytes during aging and OA. Knockdown of MFN2 with siRNA reversed the age-related metabolic changes in rat chondrocytes. Overexpression of MFN2 exacerbated inflammation and OA progress, while knockdown of MFN2 ameliorated inflammation and OA progress. Further, MFN2 could be ubiquitinated by PARKIN, declined PARKIN expression during aging and OA might result in elevated MFN2 expression.

CONCLUSIONS

Elevated MFN2 contributes to metabolic changes and inflammation during aging of rat chondrocytes and osteoarthritis.

Genetics in Cartilage Lesions: Basic Science and Therapy Approaches

Cartilage lesions have a multifactorial nature, and genetic factors are their strongest determinants. As biochemical and genetic studies have dramatically progressed over the past decade, the molecular basis of cartilage pathologies has become clearer. Several homeostasis abnormalities within cartilaginous tissue have been found, including various structural changes, differential gene expression patterns, as well as altered epigenetic regulation. However, the efficient treatment of cartilage pathologies represents a substantial challenge. Understanding the complex genetic background pertaining to cartilage pathologies is useful primarily in the context of seeking new pathways leading to disease progression as well as in developing new targeted therapies. A technology utilizing gene transfer to deliver therapeutic genes to the site of injury is quickly becoming an emerging approach in cartilage renewal. The goal of this work is to provide an overview of the genetic basis of chondral lesions and the different approaches of the most recent systems exploiting therapeutic gene transfer in cartilage repair. The integration of tissue engineering with viral gene vectors is a novel and active area of research. However, despite promising preclinical data, this therapeutic concept needs to be supported by the growing body of clinical trials.

Mayaro Virus Infects Human Chondrocytes and Induces the Expression of Arthritis-Related Genes Associated with Joint Degradation

Mayaro virus (MAYV) is an emerging arthritogenic alphavirus belonging to the Togaviridae family. Infection leads to a dengue-like illness accompanied by severe polyarthralgia. However, the molecular and cellular mechanisms of arthritis as a result of MAYV infection remain poorly understood. In the present study, we assess the susceptibility of human chondrocytes (HC), fibroblast-like synoviocytes and osteoblasts that are the major cell types involved in osteoarthritis, to infection with MAYV. We show that these cells are highly permissive to MAYV infection and that viral RNA copy number and viral titers increase over time in infected cells. Knowing that HC are the primary cells in articular cartilage and are essential for maintaining the cartilaginous matrix, gene expression studies were conducted in MAYV-infected primary HC using polymerase chain reaction (PCR) arrays. The infection of the latter cells resulted in an induction in the expression of several matrix metalloproteinases (MMP) including MMP1, MMP7, MMP8, MMP10, MMP13, MMP14 and MMP15 which could be involved in the destruction of articular cartilage. Infected HC were also found to express significantly increased levels of various IFN-stimulated genes and arthritogenic mediators such as TNF-α and IL-6. In conclusion, MAYV-infected primary HC overexpress arthritis-related genes, which may contribute to joint degradation and pathogenesis.

Association of OPG gene polymorphisms with the risk of knee osteoarthritis among Chinese people

BACKGROUND

Osteoarthritis (OA) is usually recognized to have a genetic factor, and in our study, we performed a case-control study to analyze the association between 14 single nucleotide polymorphisms (SNPs) in OPG and the risk of knee OA in a Chinese Han population.

METHODS

Fourteen OPG SNPs were assayed using MassARRAY in 393 patients clinically and radiographically diagnosed with knee OA and in 500 controls. Allelic and genotypic frequencies were compared between the groups. Logistic regression adjusting for age and gender was used to estimate risk associations between specific genotypes and knee OA by computing odds ratios (ORs) and 95% confidence intervals (95% CIs).

RESULTS

We found that the minor alleles of six SNPs in OPG were associated with an increased or decreased risk of knee OA in the allelic model analysis. In the genetic model analysis, we found that rs1905786, rs1032128, rs3134058, rs11573828, rs11573849, rs3134056, and rs1564861 were associated with an increased or decreased risk of knee OA before adjusted by sex and age. And after adjustment, three SNPs (rs1485286, rs1905786, and rs1032128) were identified to have a negative effect on knee OA.

CONCLUSION

Our results verify that genetic variants of OPG contribute to knee OA susceptibility in the population of northern China. These genetic associations may identify individuals at a particularly high risk of developing knee OA.

The Importance of the Knee Joint Meniscal Fibrocartilages as Stabilizing Weight Bearing Structures Providing Global Protection to Human Knee-Joint Tissues

The aim of this study was to review aspects of the pathobiology of the meniscus in health and disease and show how degeneration of the meniscus can contribute to deleterious changes in other knee joint components. The menisci, distinctive semilunar weight bearing fibrocartilages, provide knee joint stability, co-ordinating functional contributions from articular cartilage, ligaments/tendons, synovium, subchondral bone and infra-patellar fat pad during knee joint articulation. The meniscus contains metabolically active cell populations responsive to growth factors, chemokines and inflammatory cytokines such as interleukin-1 and tumour necrosis factor-alpha, resulting in the synthesis of matrix metalloproteases and A Disintegrin and Metalloprotease with ThromboSpondin type 1 repeats (ADAMTS)-4 and 5 which can degrade structural glycoproteins and proteoglycans leading to function-limiting changes in meniscal and other knee joint tissues. Such degradative changes are hall-marks of osteoarthritis (OA). No drugs are currently approved that change the natural course of OA and translate to long-term, clinically relevant benefits. For any pharmaceutical therapeutic intervention in OA to be effective, disease modifying drugs will have to be developed which actively modulate the many different cell types present in the knee to provide a global therapeutic. Many individual and combinatorial approaches are being developed to treat or replace degenerate menisci using 3D printing, bioscaffolds and hydrogel delivery systems for therapeutic drugs, growth factors and replacement progenitor cell populations recognising the central role the menisci play in knee joint health.

In vivo identification of Bmp2-correlation networks during fracture healing by means of a limb-specific conditional inactivation of Bmp2

Bmp2 is known to play an essential role in the initiation of fracture healing via periosteal activation. Specifically, activation and subsequent differentiation of periosteal progenitor cells requires Bmp2 signaling for activation of the osteo-chondrogenic pathway. Here, we explored the interactive transcriptional gene-gene interplays between Bmp2 and 150 known candidate genes during fracture repair. We constructed the interactive Bmp2 signaling pathways in vivo, by comparing gene expression levels prior and 24 h post femur fracture, in presence (wild type) and in absence of Bmp2 (Bmp2c/c;Prx1::cre limb-specific conditional knockout). Twenty-six differentially expressed genes (pre- vs. post-fracture), which demonstrated high correlations within each experimental condition, were used to construct the co-expression networks. Topological dynamic shifts across different co-expression networks characterized the 26 differentially expressed genes as non-redundant focal linking hubs, redundant connecting hubs, periphery genes, or non-existent. Top-ranked up- or down-regulated genes were identified and discussed. Protein-protein interactions in public databases support our findings. Thus, the co-expression networks from this study can be used for future experimental hypotheses.

A Sparse Regression Method for Group-Wise Feature Selection with False Discovery Rate Control

The method of Sorted L-One Penalized Estimation, or SLOPE, is a sparse regression method recently introduced by Bogdan et. al. [1] . It can be used to identify significant predictor variables in a linear model that may have more unknown parameters than observations. When the correlations between predictor variables are small, the SLOPE method is shown to successfully control the false discovery rate (the expected proportion of the irrelevant among all selected predictors) at a user specified level. However, the requirement for nearly uncorrelated predictors is too restrictive for genomic data, as demonstrated in our recent study [2] by an application of SLOPE to realistic simulated DNA sequence data. A possible solution is to divide the predictor variables into nearly uncorrelated groups, and to modify the procedure to select entire groups with an overall significant group effect, rather than individual predictors. Following this motivation, we extend SLOPE in the spirit of Group LASSO to Group SLOPE, a method that can handle group structures between the predictor variables, which are ubiquitous in real genomic data. Our theoretical results show that Group SLOPE controls the group-wise false discovery rate (gFDR), when groups are orthogonal to each other. For use in non-orthogonal settings, we propose two types of Monte Carlo based heuristics, which lead to gFDR control with Group SLOPE in simulations based on real SNP data. As an illustration of the merits of this method, an application of Group SLOPE to a dataset from the Framingham Heart Study results in the identification of some known DNA sequence regions associated with bone health, as well as some new candidate regions. The novel methods are implemented in the R package grpSLOPEMC , which is publicly available at https://github.com/agisga/grpSLOPEMC.

Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage

Aim

Osteoarthritis (OA) is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control the expression of genes and are likely to regulate the OA transcriptome. We performed integrative genomic analyses to define methylation-gene expression relationships in osteoarthritic cartilage.

Patients and Methods

Genome-wide DNA methylation profiling of articular cartilage from five patients with OA of the knee and five healthy controls was conducted using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, California). Other independent genome-wide mRNA expression profiles of articular cartilage from three patients with OA and three healthy controls were obtained from the Gene Expression Omnibus (GEO) database. Integrative pathway enrichment analysis of DNA methylation and mRNA expression profiles was performed using integrated analysis of cross-platform microarray and pathway software. Gene ontology (GO) analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID).

Results

We identified 1265 differentially methylated genes, of which 145 are associated with significant changes in gene expression, such as DLX5, NCOR2 and AXIN2 (all p-values of both DNA methylation and mRNA expression < 0.05). Pathway enrichment analysis identified 26 OA-associated pathways, such as mitogen-activated protein kinase (MAPK) signalling pathway (p = 6.25 × 10-4), phosphatidylinositol (PI) signalling system (p = 4.38 × 10-3), hypoxia-inducible factor 1 (HIF-1) signalling pathway (p = 8.63 × 10-3 pantothenate and coenzyme A (CoA) biosynthesis (p = 0.017), ErbB signalling pathway (p = 0.024), inositol phosphate (IP) metabolism (p = 0.025), and calcium signalling pathway (p = 0.032).

Conclusion

We identified a group of genes and biological pathwayswhich were significantly different in both DNA methylation and mRNA expression profiles between patients with OA and controls. These results may provide new clues for clarifying the mechanisms involved in the development of OA.

Cite this article: A. He, Y. Ning, Y. Wen, Y. Cai, K. Xu, Y. Cai, J. Han, L. Liu, Y. Du, X. Liang, P. Li, Q. Fan, J. Hao, X. Wang, X. Guo, T. Ma, F. Zhang. Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage. Bone Joint Res 2018;7:343–350. DOI: 10.1302/2046-3758.75.BJR-2017-0284.R1.

Functional polymorphisms in asporin and CILP together with joint loading predispose to hand osteoarthritis

Background

Osteoarthritis (OA) is the most common degenerative joint disease afflicting people in the Western world and has a strong genetic influence. The aim of this study was to examine the association of two known functional polymorphisms in the TGF-β inhibiting genes, asporin (ASPN) and cartilage intermediate layer protein (CILP), with hand OA and potential gene-occupational hand loading interaction.

Results

Statistically significant interaction of the CILP rs2073711 T and ASPN D15 alleles with hand OA was observed (OR = 2.48, 95% CI 1.27–4.85, p = 0.008) in a Finnish hand OA cohort of 543 women (aged 45–63). When stratified by variation in working tasks, low variation of working tasks increased the risk further (OR = 3.00, 95% CI 1.35–6.66, p = 0.007). Based on the analysis of ASPN and CILP protein-coding regions, functional studies were performed with one observed variant, rs41278695 in the ASPN gene. Analyses showed that bone morphogenetic protein 2 (BMP2) mediated expression of aggrecan (Agc1) and type II collagen (Col2a1) was significantly suppressed (p = 0.011 and p = 0.023, respectively) in a murine chondrocytic cell line (ATDC5) with cells stably expressing ASPN rs41278695.

Conclusions

The carriage of either ASPN D15 or CILP rs2073711 TT is associated with increased risk of symmetrical hand OA, particularly in individuals with low variation in work tasks. ASPN rs41278695 SNP had an effect on Agc1 and Col2a1 gene expression when induced with BMP-2 suggesting an effect on the cartilage extracellular matrix composition.

TGF-β Family Signaling in Connective Tissue and Skeletal Diseases

The transforming growth factor β (TGF-β) family of signaling molecules, which includes TGF-βs, activins, inhibins, and numerous bone morphogenetic proteins (BMPs) and growth and differentiation factors (GDFs), has important functions in all cells and tissues, including soft connective tissues and the skeleton. Specific TGF-β family members play different roles in these tissues, and their activities are often balanced with those of other TGF-β family members and by interactions with other signaling pathways. Perturbations in TGF-β family pathways are associated with numerous human diseases with prominent involvement of the skeletal and cardiovascular systems. This review focuses on the role of this family of signaling molecules in the pathologies of connective tissues that manifest in rare genetic syndromes (e.g., syndromic presentations of thoracic aortic aneurysm), as well as in more common disorders (e.g., osteoarthritis and osteoporosis). Many of these diseases are caused by or result in pathological alterations of the complex relationship between the TGF-β family of signaling mediators and the extracellular matrix in connective tissues.

Opportunities and challenges for the discovery and validation of proteomic biomarkers for common arthritic diseases

Osteoarthritis (OA) and rheumatoid arthritis (RA) are most prevalent among all the rheumatic diseases, and currently, there are no reliable biochemical measures for early diagnosis or for predicting who is likely to progress. Early diagnosis is important for making decisions on treatment options and for better management of patients. This narrative review highlights the first-generation biomarkers identified over the last two decades and focuses on the discovery and validation of candidate OA biomarkers from recent mass-spectrometry-based proteomic studies for diagnosis and monitoring disease outcomes in human. It discusses the challenges and opportunities for discovery of novel biomarkers and progress in the development of techniques for measuring biomarkers, and provides directions for future discovery and validation of biomarkers for OA and rheumatoid arthritis.

Association of ADAM12 gene polymorphisms with knee osteoarthritis susceptibility

Previous studies that evaluated the association between a disintegrin and metalloprotease 12 (ADAM12) gene polymorphisms and knee osteoarthritis (KOA) have given controversial and indefinite results. Therefore, we performed a meta-analysis to confirm this correlation. We searched the PubMed, Embase, and SinoMed databases for all papers published up to April 11, 2017. Overall, five different studies, totaling 2,353 cases and 3,668 controls, were retrieved on the basis of the search criteria for KOA susceptibility related to four polymorphisms (rs3740199, rs1278279, rs1871054, and rs1044122) in the ADAM12 gene. Odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the strength of this association. Publication bias was assessed using Egger's and Begg's tests. The rs3740199 G/C polymorphism was found to be associated with increased KOA risk in men (e.g., CG versus GG: OR = 1.44, 95% CI = 1.02–2.04, P = 0.040), but not in the overall analysis and in analyses of other subgroups. Significantly increased associations were also found for the rs1871054 polymorphism (e.g., C versus T allele: OR = 1.85, 95% CI = 1.49–2.30, P < 0.001). However, there were no associations for the rs1278279 and rs1044122 polymorphisms. Furthermore, no obvious evidence of publication bias was detected. Our study indicated that the rs1871054 polymorphism of ADAM12 was significantly associated with increased KOA risk.

Integrative epigenomics, transcriptomics and proteomics of patient chondrocytes reveal genes and pathways involved in osteoarthritis

Osteoarthritis (OA) is a common disease characterized by cartilage degeneration and joint remodeling. The underlying molecular changes underpinning disease progression are incompletely understood. We investigated genes and pathways that mark OA progression in isolated primary chondrocytes taken from paired intact versus degraded articular cartilage samples across 38 patients undergoing joint replacement surgery (discovery cohort: 12 knee OA, replication cohorts: 17 knee OA, 9 hip OA patients). We combined genome-wide DNA methylation, RNA sequencing, and quantitative proteomics data. We identified 49 genes differentially regulated between intact and degraded cartilage in at least two -omics levels, 16 of which have not previously been implicated in OA progression. Integrated pathway analysis implicated the involvement of extracellular matrix degradation, collagen catabolism and angiogenesis in disease progression. Using independent replication datasets, we showed that the direction of change is consistent for over 90% of differentially expressed genes and differentially methylated CpG probes. AQP1, COL1A1 and CLEC3B were significantly differentially regulated across all three -omics levels, confirming their differential expression in human disease. Through integration of genome-wide methylation, gene and protein expression data in human primary chondrocytes, we identified consistent molecular players in OA progression that replicated across independent datasets and that have translational potential.

Association of ADAM 12 polymorphisms with the risk of knee osteoarthritis: meta-analysis of 5048 cases and 6848 controls

Several studies have suggested the association between ADAM 12 polymorphisms and the risk of osteoarthritis (OA), but the results remained controversial. Therefore, we designed a meta-analysis to systematically evaluate the association on this issue. A literature search for eligible studies was conducted in PubMed, Web of Science and Google Scholar databases. The association between ADAM 12 polymorphisms and knee OA risk was calculated by odds ratios (ORs) and 95% confidence intervals (CIs). Study heterogeneity, sensitivity and publication bias analyses were also conducted. Ten articles covering 5048 cases and 6848 controls met our criteria for the final analysis. We found that the rs1871054 was significantly associated with the risk of knee OA (allele model OR 1.72, 95% CI 1.43-2.07, P < 0.001; additive model: OR 2.06, 95% CI 1.19-3.56, P = 0.010; dominant model: OR 2.45, 95% CI 1.85-3.25, P < 0.001; recessive model: OR 1.54, 95% CI 1.13-2.10, P = 0.007). rs1044122 was significantly associated with knee OA susceptibility in recessive model (OR 1.45, 95% CI 1.03-2.04, P = 0.031). For rs3740199 and rs1278279, no significant associations with knee OA were found. In the stratified analysis by gender, significant association was identified with the risk of knee OA for rs3740199 in men in allele model (OR 2.41, 95% CI 1.51-3.84, P < 0.001), dominant model (OR 2.68, 95% CI 1.17-6.14, P = 0.02) and recessive model (OR 3.51, 95% CI 1.68-7.36, P = 0.001), but not for additive model (OR 1.30, 95% CI 0.81-2.08, P = 0.28). This meta-analysis suggests that the ADAM 12 genetic polymorphisms rs1871054 and rs1044122 might be associated with risk of knee OA; rs3740199 might be associated with risk of knee OA in men. Further well-designed and large scale studies are warranted to validate these associations.

Increased DNA Methylation and Reduced Expression of Transcription Factors in Human Osteoarthritis Cartilage

OBJECTIVE

To analyze the methylome of normal and osteoarthritic (OA) knee articular cartilage and to determine the role of DNA methylation in the regulation of gene expression in vitro.

METHODS

DNA was isolated from human normal (n = 11) and OA (n = 12) knee articular cartilage and analyzed using the Infinium HumanMethylation450 BeadChip array. To integrate methylation and transcription, RNA sequencing was performed on normal and OA cartilage and validated by quantitative polymerase chain reaction. Functional validation was performed in the human TC28 cell line and primary chondrocytes that were treated with the DNA methylation inhibitor 5-aza-2'-deoxycytidine (5-aza-dC).

RESULTS

DNA methylation profiling revealed 929 differentially methylated sites between normal and OA cartilage, comprising a total of 500 individual genes. Among these, 45 transcription factors that harbored differentially methylated sites were identified. Integrative analysis and subsequent validation showed a subset of 6 transcription factors that were significantly hypermethylated and down-regulated in OA cartilage (ATOH8, MAFF, NCOR2, TBX4, ZBTB16, and ZHX2). Upon 5-aza-dC treatment, TC28 cells showed a significant increase in gene expression for all 6 transcription factors. In primary chondrocytes, ATOH8 and TBX4 were increased after 5-aza-dC treatment.

CONCLUSION

Our findings reveal that normal and OA knee articular cartilage have significantly different methylomes. The identification of a subset of epigenetically regulated transcription factors with reduced expression in OA may represent an important mechanism to explain changes in the chondrocyte transcriptome and function during OA pathogenesis.

ADAMTS and ADAM metalloproteinases in osteoarthritis - looking beyond the 'usual suspects'

INTRODUCTION

Matrix metalloproteinases (MMPs) and 'aggrecanase' a disintegrin and metalloproteinase with thrombospondin motifs (ADAMTSs) are well established to play key roles in osteoarthritis (OA) through degradation of extracellular matrix (ECM) type II collagen and aggrecan, and are thus potential targets for development of OA therapies.

OBJECTIVE

This paper aims to provide a comprehensive review of the expression and potential roles of other, lesser-known ADAMTSs and related adamalysins (or a disintegrin and metalloproteinases (ADAMs)) in cartilage, with a view to identifying potentially protective or homeostatic metalloproteinases in the joint and informing consequent selective inhibitor design.

DESIGN

A comprehensive literature search was performed using PubMed terms 'osteoarthritis' and 'ADAMTS' or 'ADAM'.

RESULTS

Several ADAMTSs and ADAMs were identified as having reportedly increased expression in OA. These include enzymes likely to play roles in cartilage matrix anabolism (e.g., the procollagen N-proteinases ADAMTS-2, ADAMTS-3 and ADAMTS-14), chondrocyte differentiation and proliferation (e.g., ADAM9, ADAM10, ADAM12), as well as enzymes contributing to cartilage catabolism (e.g., Cartilage oligomeric protein (COMP)-degrading ADAMTS-7 and ADAMTS-12).

CONCLUSIONS

In addition to the well-characterised MMPs, ADAMTS-4 and ADAMTS-5, many other ADAMTSs and ADAMs are expressed in cartilage and several show significantly altered expression in OA. Studies aimed at elucidating the pathophysiological roles of these enzymes in cartilage will contribute to our understanding of OA pathogenesis and enable design of targeted inhibitors that effectively target metalloproteinase-mediated cartilage degradation while sparing cartilage repair pathways.

The association of ADAM12 polymorphism with osteoarthritis susceptibility: a meta-analysis

BACKGROUND

The pathology of osteoarthritis (OA) is partly attributed to genetic factors; however, the role of ADAM12 polymorphism is still controversial. It is necessary to perform a meta-analysis to investigate this possible correlation.

METHODS

Case-control studies on the association between OA susceptibility and ADAM12 polymorphism were comprehensively collected by searching PubMed, Embase, and Web of Science. Odds ratios (ORs) and 95% confidence intervals (CIs) were pooled to evaluate OA risk that was possibly conferred by ADAM12 variant. The analyses were performed not only among general population but also in male and female groups.

RESULTS

A total of 8 studies with 10 populations were finally included in this meta-analysis. In the general population, 4 comparisons were carried out (C allele vs G allele, CC vs GG, GC + CC vs GG, and CC vs GC + GG) and found that ADAM12 rs3740199 polymorphism was not associated with increased OA vulnerability. On the other hand, the analyses stratified by gender made 5 comparisons (C allele vs G allele, CC vs GG, GC vs GG, GC + CC vs GG, and CC vs GC + GG). It was shown that rs3740199 polymorphism (GC + CC vs GG) was a risk factor for OA among male patients (OR =1.45, 95% CI =1.04-2.01). Sensitivity analysis indicated that it was an unstable outcome. No correlation was identified in women. Neither heterogeneity nor publication bias was detected in the analyses mentioned above.

CONCLUSION

ADAM12 rs3740199 polymorphism is likely to be associated with OA susceptibility among male patients, other than the general population. More studies are needed to confirm this observation. The mechanism by which ADAM12 variant plays a role in OA pathogenesis is also warranted and important for interpreting this possible correlation.

Genetic effects of rs3740199 polymorphism in ADAM12 gene on knee osteoarthritis: a meta-analysis

BACKGROUND

Knee osteoarthritis (OA) is a complex arthritic condition in which genetic factors play an important role. ADAM12 gene is one of the recognized candidate genes although the results are conflicting. To derive a more precise estimation of the association between rs3740199 polymorphism in ADAM12 gene and risk of knee OA, we performed a meta-analysis based on six related studies, including a total of 2185 cases and 3716 controls.

METHODS

A comprehensive search was performed to identify related studies up to April 14, 2017. We used odds ratios (ORs) with 95% confidence intervals (CIs) to assess the strength of the association. Different genetic models were used to assess the pooled and stratified data.

RESULTS

Overall, no significant association was found in all genetic models (C vs. G, OR = 0.983, 95% CI = 0.910-1.061; CC vs. GG, OR = 1.033, 95% CI = 0.851-1.255; CG vs. GG, OR = 1.030, 95% CI = 0.877-1.209; CC/CG vs. GG, OR = 1.031, 95% CI = 0.886-1.201; CC vs. CG/GG, OR = 1.017, 95% CI = 0.868-1.190). When stratified by ethnicity, no significant association was found.

CONCLUSIONS

This meta-analysis suggested that the rs3740199 polymorphism does not contribute to the development of knee OA. Additional well-designed large studies are required to confirm these findings in different populations.

The genetics behind osteoarthritis: Asian focus

Osteoarthritis (OA), a degenerative joint disease, is the most common form of arthritis in the elderly. Problems arising from the condition are not only health, physical, and psychological, but also economical and social. OA is a complex disease causes by environmental factors (obesity, female gender, advancing age) and genetics. A strong genetic contribution to primary OA has been observed in several studies. Asian population is undoubtedly genetically different from European, whom most studies have been performed. This review systematically compares results of molecules involved with OA from studies performed on Asian and Caucasian populations. The clear differences between Asian and Caucasian populations may be from the sex- and ethnic-specific nature of the disease, as well as from the inclusion criteria of the studies. Global collaboration is highly significant in understanding genetic underlying OA. Knowledge gained from genetic study is important to the development of prevention and therapeutic intervention, and to identify individuals at risk of developing severe or progressive OA.

Association Between Cartilage Intermediate Layer Protein and Degeneration of Intervertebral Disc: A Meta-analysis

STUDY DESIGN

Meta-analysis to collect all eligible studies for investigating whether cartilage intermediate layer protein (CILP) gene rs2073711 polymorphism is associated with susceptibility to intervertebral disc degeneration (IVD).

OBJECTIVE

To investigate the association between CILP gene rs2073711 and IVD.

SUMMARY OF BACKGROUND DATA

CILP is one of the limited genes that is being suggested to be associated with IVD. However, previous association study results of CILP and IVD are conflicted.

METHODS

The PubMed, ISI Web of Science, Science Direct Online, and China National Knowledge Infrastructure database were searched for relevant articles till November 1, 2015. The odds ratio (OR) and related 95% confidence interval (CI) were drawn out and merged to assess the strength of association by STATA version 12.0.

RESULTS

A total of 1551 IVD cases and 1793 controls from the 5 studies were used in this study, comprising four Asian populations and two European populations. Meta-analysis found that the CILP 2073711 polymorphism was significantly associated with IVD risk. The estimated OR of risk allele T of rs2073711 was 1.36. (95% CI = 1.18-1.55, P < 0.001). No significant heterogeneity and publication bias were observed in this study.

CONCLUSION

Our results confirm the positive association between CILP and IVD, providing novel clues for clarifying the role of CILP in the development of IVD.

LEVEL OF EVIDENCE

2.

A bivariate genome-wide association study identifies ADAM12 as a novel susceptibility gene for Kashin-Beck disease

Kashin-Beck disease (KBD) is a chronic osteoarthropathy, which manifests as joint deformities and growth retardation. Only a few genetic studies of growth retardation associated with the KBD have been carried out by now. In this study, we conducted a two-stage bivariate genome-wide association study (BGWAS) of the KBD using joint deformities and body height as study phenotypes, totally involving 2,417 study subjects. Articular cartilage specimens from 8 subjects were collected for immunohistochemistry. In the BGWAS, ADAM12 gene achieved the most significant association (rs1278300 p-value = 9.25 × 10⁻⁹) with the KBD. Replication study observed significant association signal at rs1278300 (p-value = 0.007) and rs1710287 (p-value = 0.002) of ADAM12 after Bonferroni correction. Immunohistochemistry revealed significantly decreased expression level of ADAM12 protein in the KBD articular cartilage (average positive chondrocyte rate = 47.59 ± 7.79%) compared to healthy articular cartilage (average positive chondrocyte rate = 64.73 ± 5.05%). Our results suggest that ADAM12 gene is a novel susceptibility gene underlying both joint destruction and growth retardation of the KBD.

Association of rs1137101 polymorphism in LEPR and susceptibility to knee osteoarthritis in a Northwest Chinese Han population

Background

Osteoarthritis (OA) is a complex arthritic condition in which the genetic factor plays a crucial role. A single nucleotide polymorphism (SNP), rs1137101 (Gln223Arg) of leptin receptor (LEPR) gene has been demonstrated to be associated with susceptibility to knee OA. However, this association in Chinese Han population has never been examined. The present study aimed to investigate whether Gln223Arg was related to knee OA susceptibility in a Northwest Chinese population with Han ethnicity.

Methods

Gln223Arg polymorphisms were genotyped in 587 patients with confirmed knee OA and in 628 age- and sex-matched healthy controls using polymerase chain reaction-restriction fragment length polymorphism analysis. Besides, LEPR genotypes were verified by direct DNA sequencing analysis on PCR products.

Results

The genotype and allele frequencies in LEPR SNP rs1137101 were significantly different between cases and control groups (chi-square = 6.52, P = 0.038 for genotype and chi-square = 5.06, P = 0.024 for allele frequencies; respectively). Rs1137101 was correlated with knee OA in the dominant genetic model (GG + GA versus AA) (P = 0.016) and a higher G allele frequency existed (P = 0.024) among all patients with knee OA and controls. On stratification analysis, the genotype GG and G allele were associated with susceptibility to knee OA in females, both young (≤65 years) and old groups (>65 years) and patients with mild knee OA.

Conclusions

Our finding suggested that the genetic variant of LEPR gene rs1137101 is independently related to knee OA susceptibility in Northwest Chinese population with Han ethnicity and may serve as a potential biomarker to determine risk of knee OA.

Association between ADAM12 polymorphism and knee osteoarthritis in Thai population

BACKGROUND

Osteoarthritis (OA), a common degenerative joint disorder in the elderly, is characterized by the destruction of articular cartilage, bony outgrowths at joint margins, and synovitis. The objective of this study was to evaluate whether there is an association between the ADAM12 (rs3740199) polymorphism and susceptibility to knee OA in a Thai population.

METHODS

Genomic deoxyribonucleic acid (DNA) was isolated from 200 Thai knee OA patients and 200 healthy controls. High resolution melting analysis was used to detect ADAM12 polymorphisms. The melt profile of all DNA samples was generated on the CFX96™ real-time polymerase chain reaction system and analyzed by Precision Melt Analysis™ software. The genotype distributions and allele frequencies of ADAM12 were compared between groups using the StatCalc program.

RESULTS

The significant associations were shown from the C allele (OR=2.10, 95% CI=1.16-3.79, P=0.008) and the CC genotype (OR=4.28, 95% CI=1.21-15.72, P=0.01) in male knee OA patients. No significant association was observed in female patients.

CONCLUSION

The rs3740199 in ADAM12 was associated with knee OA susceptibility in Thai male patients, and individuals with the CC genotype carried the highest risk when compared with the GG and GC genotypes.

CLINICAL RELEVANCE

The rs3740199 polymorphism of the ADAM12 gene can potentially be used to determine genetically high-risk subgroup of knee osteoarthritis and to better understand the pathogenesis of knee osteoarthritis.

ADAM33 and ADAM12 genetic polymorphisms and their expression in Egyptian children with asthma

BACKGROUND

The ADAM family is involved in some pathologic processes, such as inflammation and asthma.

OBJECTIVES

To assess the association between ADAM33 and ADAM12 single-nucleotide polymorphisms (SNPs) with asthma risk and severity and to investigate the effect of ADAM33 and ADAM12 polymorphisms on expression of these proteases in sputum.

METHODS

Two SNPs of the ADAM33 gene, F+1 (rs511898) G/A and ST+4 (rs44707) A/C, and 2 SNPs of the ADAM12 gene, rs3740199 and rs1871054, were analyzed in 400 asthma cases and 200 controls aged 3 to 14 years using the polymerase chain reaction-restriction fragment length polymorphism method. Messenger RNA expression profile of ADAM33 and ADAM12 proteases in sputum from studied groups was determined by reverse transcription polymerase chain reaction.

RESULTS

ADAM33 F+1 homozygous mutant genotype (AA) and ST+4 heterozygous and homozygous mutant genotype (AC and CC) and mutant alleles of both polymorphisms were significantly associated with asthma risk and severity in moderate and severe subgroups. Patients with the ADAM12 (rs3740199) CC genotype were at increased risk for moderate and severe asthma. Messenger RNA levels of ADAM12 were significantly increased in asthmatic children compared with controls, whereas we were not able to detect the expression of ADAM33 in the sputum of the groups studied. The ADAM12 expression was significantly higher in homozygous CC (variant type) compared with homozygous GG (wild type) of both ADAM12 rs3740199 and rs1871054 in the asthmatic group.

CONCLUSION

Our analysis suggests a likely role for ADAM33 and ADAM12 in the development of asthma in Egyptian children. Furthermore, ADAM12 polymorphisms may affect ADAM12 expression in asthma.

What Are the Prognostic Factors for Radiographic Progression of Knee Osteoarthritis? A Meta-analysis

BACKGROUND

A previous systematic review on prognostic factors for knee osteoarthritis (OA) progression showed associations for generalized OA and hyaluronic acid levels. Knee pain, radiographic severity, sex, quadriceps strength, knee injury, and regular sport activities were not associated. It has been a decade since the literature search of that review and many studies have been performed since then investigating prognostic factors for radiographic knee OA progression.

QUESTIONS/PURPOSES

The purpose of this study is to provide an updated systematic review of available evidence regarding prognostic factors for radiographic knee OA progression.

METHODS

We searched for observational studies in MEDLINE and EMBASE. Key words were: knee, osteoarthritis (or arthritis, or arthrosis, or degenerative joint disease), progression (or prognosis, or precipitate, or predictive), and case-control (or cohort, or longitudinal, or follow-up). Studies fulfilling the inclusion criteria were assessed for methodologic quality according to established criteria for reviews on prognostic factors in musculoskeletal disorders. Data were extracted and results were pooled if possible or summarized according to a best-evidence synthesis. A total of 1912 additional articles were identified; 43 met our inclusion criteria. The previous review contained 36 articles, thus providing a new total of 79 articles. Seventy-two of the included articles were scored high quality, the remaining seven were low quality.

RESULTS

The pooled odds ratio (OR) of two determinants showed associations with knee OA progression: baseline knee pain (OR, 2.38 [95% CI, 1.74-3.27) and Heberden nodes (OR, 2.66 [95% CI, 1.46-8.84]). Our best-evidence synthesis showed strong evidence that varus alignment, serum hyaluronic acid, and tumor necrosis factor-α are associated with knee OA progression. There is strong evidence that sex, former knee injury, quadriceps strength, smoking, running, and regular performance of sports are not associated with knee OA progression. Evidence for the majority of determined associations, however, was limited, conflicting, or inconclusive.

CONCLUSIONS

Baseline knee pain, presence of Heberden nodes, varus alignment, and high levels of serum markers hyaluronic acid and tumor necrosis factor-α predict knee OA progression. Sex, knee injury, and quadriceps strength, among others, did not predict knee OA progression. Large variation remains in definitions of knee OA and knee OA progression. Clinical studies should use more consistent definitions of these factors to facilitate data pooling by future meta-analyses.

Analysis of single nucleotide polymorphisms within ADAM12 and risk of knee osteoarthritis in a Chinese Han population

Objective. Osteoarthritis (OA) is a complex arthritic condition in which the genetic factor plays a major role. One of the candidate genes of is the ADAM12 gene, but no consistency has been reached till now. This study aims to investigate the potential role of four single nucleotide polymorphisms (SNPs) of the ADAM12 gene in susceptibility to knee OA and its progression in Chinese Han population. Methods. The rs1278279, rs3740199, rs1044122, and rs1871054 polymorphisms were genotyped and compared in a population based cohort consisting of 164 OA subjects and 200 age- and gender-matched controls. Results. The SNP rs1871054 was found with increased risk of OA susceptibility in comparing the genotype frequencies between the case and control groups no matter for which model of comparison (allele level, dominant model, recessive model, and extreme genotype model). Additionally, the SNP rs1871054 was found associated with increased OA severity according to the K/L grade. Conclusion. In summary, we have identified that the rs1871054 variant within the ADAM12 gene is a risk factor for increased osteoarthritis susceptibility and severity.

Improved prediction of knee osteoarthritis progression by genetic polymorphisms: the Arthrotest Study

OBJECTIVE

The aim of this study was to develop a genetic prognostic tool to predict radiographic progression towards severe disease in primary knee OA (KOA) patients.

METHODS

This investigation was a cross-sectional, retrospective, multicentric association study in 595 Spanish KOA patients. Caucasian patients aged ≥40 years at the time of diagnosis of primary KOA of Kellgren-Lawrence grade 2 or 3 were included. Patients who progressed to Kellgren-Lawrence score 4 or who were referred for total knee replacement within 8 years after diagnosis were classified as progressors to severe disease. Clinical variables of the initial stages of the disease (gender, BMI, age at diagnosis, OA in the contralateral knee, and OA in other joints) were registered as potential predictors. Single nucleotide polymorphisms and clinical variables with an association of P < 0.05 were included in the multivariate analysis using forward logistic regression.

RESULTS

A total of 23 single nucleotide polymorphisms and the time of primary KOA diagnosis were significantly associated with KOA severe progression in the exploratory cohort (n = 220; P < 0.05). The predictive accuracy of the clinical variables was limited: area under the curve (AUC) = 0.66. When genetic variables were added to the clinical model (full model), the prediction of KOA progression was significantly improved (AUC = 0.82). Combining only genetic variables (rs2073508, rs10845493, rs2206593, rs10519263, rs874692, rs7342880, rs780094 and rs12009), a predictive model with good accuracy was also obtained (AUC = 0.78). The predictive ability for KOA progression of the full model was confirmed on the replication cohort (two-sample Z-test; n = 62; P = 0.190).

CONCLUSION

An accurate prognostic tool to predict primary KOA progression has been developed based on genetic and clinical information from OA patients.

Quantitative proteomic profiling of human articular cartilage degradation in osteoarthritis

Osteoarthritis (OA) is the most common rheumatic pathology and is characterized primarily by articular cartilage degradation. Despite its high prevalence, there is no effective therapy to slow disease progression or regenerate the damaged tissue. Therefore, new diagnostic and monitoring tests for OA are urgently needed, which would also promote the development of alternative therapeutic strategies. In the present study, we have performed an iTRAQ-based quantitative proteomic analysis of secretomes from healthy human articular cartilage explants, comparing their protein profile to those from unwounded (early disease) and wounded (advanced disease) zones of osteoarthritic tissue. This strategy allowed us to identify a panel of 76 proteins that are distinctively released by the diseased tissue. Clustering analysis allowed the classification of proteins according to their different profile of release from cartilage. Among these proteins, the altered release of osteoprotegerin (decreased in OA) and periostin (increased in OA), both involved in bone remodelling processes, was verified in further analyses. Moreover, periostin was also increased in the synovial fluid of OA patients. Altogether, the present work provides a novel insight into the mechanisms of human cartilage degradation and a number of new cartilage-characteristic proteins with possible biomarker value for early diagnosis and prognosis of OA.

Associations of polymorphisms in the bone morphogenetic protein-2 gene with risk and prognosis of osteosarcoma in a Chinese population

Osteosarcoma is the most common type of bone cancer in adolescence. Bone morphogenetic protein-2 (BMP-2) plays important roles in the development of bone and cartilage and in inhibiting the tumorigenicity of cancer stem cells in human osteosarcoma cell line. The aim of this study was to examine whether polymorphisms in the BMP2 gene are associated with osteosarcoma risk and prognosis in Chinese population. Five single nucleotide polymorphisms (SNP) in the BMP2 gene were genotyped in a case-control study, including 203 osteosarcoma patients and 406 cancer-free controls. We found that rs3178250 TT genotype was associated with significant increased osteosarcoma risk (age-adjusted odds ratio (OR) = 2.06, 95% confidence intervals (CI) of 1.23-3.45) compared with CC genotype. Subjects carrying the AA genotype of rs1005464 had significant decreased cancer risk (age-adjusted OR = 0.44, 95% CI of 0.23-0.85) compared with those carrying the GG genotype. Haplotype analysis also showed that carriers of the G-T-T-G and A-T-T-G haplotypes (rs235764-rs3178250-rs235768-rs1005464) had significant increased risks of osteosarcoma (age-adjusted OR = 1.85, 95% CI of 1.28-2.66 and age-adjusted OR = 1.51, 95% CI of 1.06-2.16) compared with the G-C-T-A haplotype carriers. Besides, rs1005464 was an independent prognostic factor for osteosarcoma patients (GA vs. GG: age-adjusted hazard radio (HR) = 0.60, 95% CI of 0.36-0.99). Our data suggest that genetic mutations in the BMP2 gene are associated with osteosarcoma risk and prognosis in a Chinese population.

Mitochondrial DNA (mtDNA) haplogroups influence the progression of knee osteoarthritis. Data from the Osteoarthritis Initiative (OAI)

OBJECTIVE

To evaluate the influence of the mtDNA haplogroups on knee osteoarthritis progression in Osteoarthritis Initiative (OAI) participants through longitudinal data from radiographs and magnetic resonance imaging (MRI).

METHODS

Four-year knee osteoarthritis progression was analyzed as increase in Kellgren and Lawrence (KL) grade, in addition to increase in OARSI atlas grade for joint space narrowing (JSN), osteophytes and subchondral sclerosis in the tibia medial compartment of 891 Caucasian individuals from the progression subcohort. The influence of the haplogroups on the rate of structural progression was also assessed as the four-year change in minimum joint space width (mJSW in millimetres) in both knees of (n = 216) patients with baseline unilateral medial-tibiofemoral JSN. Quantitative cartilage measures from longitudinal MRI data were those related to cartilage thickness and volume with a 24 month follow-up period (n = 381).

RESULTS

During the four-year follow-up period, knee OA patients with the haplogroup T showed the lowest increase in KL grade (Hazard Risk [HR] = 0.499; 95% Confidence Interval [CI]: 0.261-0.819; p<0.05) as well as the lowest cumulative probability of progression for JSN (HR = 0.547; 95% CI: 0.280-0.900; p<0.05), osteophytes (HR = 0.573; 95% CI: 0.304-0.893; p<0.05) and subchondral sclerosis (HR = 0.549; 95% CI: 0.295-0.884; p<0.05). They also showed the lowest decline in mJSW (standardized response means (SRM) = -0.39; p = 0.037) in those knees without baseline medial JSN (no-JSN knees). Normalized cartilage volume loss was significantly lower in patients carrying the haplogroup T at medial tibia femoral (SRM = -0.33; p = 0.023) and central medial femoral (SRM = -0.27; p = 0.031) compartments. Cartilage thickness loss was significantly lower in carriers of haplogroup T at central medial tibia-femoral (SRM = -0.42; p = 0.011), medial tibia femoral (SRM = -0.32; p = 0.018), medial tibia anterior (SRM = +0.31; p = 0.013) and central medial femoral (SRM = -0.19; p = 0.013) compartments.

CONCLUSIONS

Mitochondrial genome seems to play a role in the progression of knee osteoarthritis. mtDNA variation could improve identification of patients predisposed to faster or severe progression of the disease.

Associations between vitamin D receptor gene polymorphisms and osteoarthritis: an updated meta-analysis

OBJECTIVE

Vitamin D receptor (VDR) gene polymorphisms may be associated with the risk of OA, however, evidence for this is controversial. This meta-analysis aims to confirm whether VDR gene polymorphisms are associated with OA.

METHODS

Meta-analyses on the association between OA and VDR ApaI, BsmI, TaqI and FokI polymorphisms were conducted using allele and homozygote contrast and contrasts in the recessive and dominant models. Stratification analyses by different demographic regions (Europe vs Asian) were also performed and pooled odds ratios (ORs) were obtained using the random effects model if the results were heterogeneous.

RESULTS

A total of 13 relevant studies involving OA patients (n = 2104) and controls (n = 2939) were included in the analysis. There were significant associations between VDR ApaI polymorphisms and OA in the Asian population (A vs a: OR= 1.16, 95% CI 1.02, 1.32, P = 0.025; AA vs Aa/aa: OR= 1.36, 95% CI 1.04, 1.77, P = 0.025; AA vs aa: OR= 1.35, 95% CI 1.00, 1.80, P = 0.047), but not in the whole population. There was also a statistically significant association between FokI polymorphism and OA (FF vs Ff/ff: OR= 0.65, 95% CI 0.44, 0.95, P = 0.024); however, this result was derived from only two studies. No significant associations were found between VDR TaqI and BsmI polymorphisms and OA.

CONCLUSION

There are modest but statistically significant associations between VDR ApaI polymorphisms and the susceptibility of OA in the Asian population.

Mitochondrial DNA haplogroups modulate the radiographic progression of Spanish patients with osteoarthritis

Not all patients with osteoarthritis (OA) show the same disease progression, as some of them remain relatively stable over time, while others progress to severe structural deterioration of the joint. In this sense, the main goal of both genetic and protein biomarkers in OA is to predict not only the risk of OA at an earlier stage of the disease but also which OA patients are more likely to progress to severe disease. Taking into account the incidence of the mitochondria and the mtDNA haplogroups in the pathogenesis of OA, the main objective of this work was to evaluate the incidence of the mtDNA haplogroups in the radiographic progression of the OA disease in a well-characterized follow-up cohort of Spanish patients. DNA from 281 OA patients from Hospital Universitario A Coruña was isolated to determine the European mtDNA haplogroups. Knee or hip radiographs from all affected joints were obtained at two time points with at least 36 months apart. Radiographs were evaluated using the Kellgren/Lawrence (K/L) scale; radiographic OA progression was defined as any radiographic worsening of the K/L joint score. Statistical analyses included Kaplan-Meier survival curves and Cox regression models. Patients belonging to the cluster TJ showed a slower radiographic OA progression than patients in the cluster KU (p = 0.036). Moreover, patients carrying the most common mtDNA haplogroup H are more apt to require total joint replacement surgery than non-H patients (p = 0.049). The inherited mitochondrial variants influence the radiographic progression of OA and could be considered among the genetic variants taken into account when the radiographic progression of OA is analyzed.

Evaluation of the effectiveness and tolerability of controlled-release diclofenac-potassium versus immediate-release diclofenac-potassium in the treatment of knee osteoarthritis

Abstract.

BACKGROUND

A controlled-release (CR) form of diclofenac-potassium has been developed, which delivers 100 mg over the course of 24 hours. This formulation is administered QD and provides steady plasmatic levels of the drug.

OBJECTIVE

The aim of this study was to compare the effectiveness and tolerability of CR diclofenac-potassium versus the immediate-release (IR) formulation, when used for treatment of pain in patients with knee osteoarthritis.

METHODS

This prospective, randomized, double-blind, comparative, multi-center, parallel-group study was conducted in male and female patients who had been previously diagnosed with knee osteoarthritis. Inclusion criteria included knee joint pain and ≥3 of the following: age >50 years, morning rigidity lasting <30 minutes, crackling in the joint, pain with applied pressure to the bones, bone hypertrophy, absence of articular heat, and a radiology status of I to III on the Kellgren-Lawrence scale. Patients were randomly divided into 1 of 2 equal-sized groups: 1 group received diclofenac-potassium IR 50 mg BID for 30 days and 1 group received diclofenac-potassium CR 100 mg QD for 30 days. Patients were assessed at baseline and again at 15 and 30 days after initiation of treatment with a physical examination, pain measurement via 100-mm visual analog scale (VAS), and Western Ontario McMaster (WOMAC) osteoarthritis index questionnaire. Adverse events (AEs) were assessed by direct interrogation, hematology controls, blood chemistry, hepatic tests, coagulation tests, and urine tests performed on patients before treatment initiation and on day 30.

RESULTS

Sixty-five patients were screened and 62 patients (mean [SD]age, 61.8 [8.9] years; mean [SD] weight, 71.3 [12.4] kg; female sex, 55 [88.9%]) were included in the study; each study group had 31 patients. After 30 days, both products were equally effective in relieving pain, as measured by VAS (IR, 17.3 vs CR, 21.6; P = NS), and changes in the WOMAC score (IR, 14.5 vs CR, 19.2; P = NS). Significantly more patients in the IR group reported feeling better after 30 days than in the CR group (94% vs 76%; P = 0.002) and, according to the physician's opinion, significantly more patients treated with diclofenac-potassium IR felt better (97% vs 83%; P = 0.03). Significantly more patients in the IR group required rescue medication than those in the CR group (36% vs 26%; P = 0.03). In the CR group, 7 patients experienced AEs: 6 were gastrointestinal (ie, pyrosis, epigastralgia, dyspepsia) and 1 patient experienced increased arterial pressure. One patient from this group discontinued treatment due to a lack of efficacy. In the IR group, 6 patients experienced AEs (ie, tachycardia, epigastralgia, and pyrosis). One patient discontinued because of AEs, and 3 withdrew due to a lack of efficacy.

CONCLUSION

Based on the results from this small study in a Venezuelan population, both IR and CR formulations of diclofenac-potassium have similar effectiveness and tolerability profiles.