A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip

Genomic Determinants of Knee Joint Biomechanics: An Exploration into the Molecular Basis of Locomotor Function, a Narrative Review

In recent years, the nexus between genetics and biomechanics has garnered significant attention, elucidating the role of genomic determinants in shaping the biomechanical attributes of human joints, specifically the knee. This review seeks to provide a comprehensive exploration of the molecular basis underlying knee joint locomotor function. Leveraging advancements in genomic sequencing, we identified specific genetic markers and polymorphisms tied to key biomechanical features of the knee, such as ligament elasticity, meniscal resilience, and cartilage health. Particular attention was devoted to collagen genes like COL1A1 and COL5A1 and their influence on ligamentous strength and injury susceptibility. We further investigated the genetic underpinnings of knee osteoarthritis onset and progression, as well as the potential for personalized rehabilitation strategies tailored to an individual's genetic profile. We reviewed the impact of genetic factors on knee biomechanics and highlighted the importance of personalized orthopedic interventions. The results hold significant implications for injury prevention, treatment optimization, and the future of regenerative medicine, targeting not only knee joint health but joint health in general.

Genetic Risk Factors for End-Stage Hip Osteoarthritis Treated With Total Hip Arthroplasty: A Genome-wide Association Study

BACKGROUND

Although a genetic component to hip osteoarthritis (OA) has been described, focused evaluation of the genetic components of end-stage disease is limited. We present a genomewide association study for patients undergoing total hip arthroplasty (THA) to characterize the genetic risk factors associated with end-stage hip osteoarthritis (ESHO), defined as utilization of the procedure.

METHODS

Patients who underwent primary THA for hip OA were identified in a national patient data repository using administrative codes. Fifteen thousand three hundred and fifty-five patients with ESHO and 374,193 control patients were identified. Whole genome regression of genotypic data for patients who underwent primary THA for hip OA corrected for age, sex, and body mass index (BMI) was performed. Multivariate logistic regression models were used to evaluate the composite genetic risk from the identified genetic variants.

RESULTS

There were 13 significant genes identified. Composite genetic factors resulted in an odds ratio 1.04 for ESHO (P < .001). The effect of genetics was lower than that of age (Odds Ratio (OR): 2.38; P < .001) and BMI (1.81; P < .001).

CONCLUSION

Multiple genetic variants, including 5 novel loci, were associated with end-stage hip OA treated with primary THA. Age and BMI were associated with greater odds of developing end-stage disease when compared to genetic factors.

Maternal pre-pregnancy BMI, offspring epigenome-wide DNA methylation, and childhood obesity: findings from the Boston Birth Cohort

BACKGROUND

Maternal pre-pregnancy obesity is an established risk factor for childhood obesity. Investigating epigenetic alterations induced by maternal obesity during fetal development could gain mechanistic insight into the developmental origins of childhood obesity. While obesity disproportionately affects underrepresented racial and ethnic mothers and children in the USA, few studies investigated the role of prenatal epigenetic programming in intergenerational obesity of these high-risk populations.

METHODS

This study included 903 mother-child pairs from the Boston Birth Cohort, a predominantly urban, low-income minority birth cohort. Mother-infant dyads were enrolled at birth and the children were followed prospectively to age 18 years. Infinium Methylation EPIC BeadChip was used to measure epigenome-wide methylation level of cord blood. We performed an epigenome-wide association study of maternal pre-pregnancy body mass index (BMI) and cord blood DNA methylation (DNAm). To quantify the degree to which cord blood DNAm mediates the maternal BMI-childhood obesity, we further investigated whether maternal BMI-associated DNAm sites impact birthweight or childhood overweight or obesity (OWO) from age 1 to age 18 and performed corresponding mediation analyses.

RESULTS

The study sample contained 52.8% maternal pre-pregnancy OWO and 63.2% offspring OWO at age 1-18 years. Maternal BMI was associated with cord blood DNAm at 8 CpG sites (genome-wide false discovery rate [FDR] < 0.05). After accounting for the possible interplay of maternal BMI and smoking, 481 CpG sites were discovered for association with maternal BMI. Among them 123 CpGs were associated with childhood OWO, ranging from 42% decrease to 87% increase in OWO risk for each SD increase in DNAm. A total of 14 identified CpG sites showed a significant mediation effect on the maternal BMI-child OWO association (FDR < 0.05), with mediating proportion ranging from 3.99% to 25.21%. Several of these 14 CpGs were mapped to genes in association with energy balance and metabolism (AKAP7) and adulthood metabolic syndrome (CAMK2B).

CONCLUSIONS

This prospective birth cohort study in a high-risk yet understudied US population found that maternal pre-pregnancy OWO significantly altered DNAm in newborn cord blood and provided suggestive evidence of epigenetic involvement in the intergenerational risk of obesity.

Novel genetic loci associated with osteoarthritis in multi-ancestry analyses in the Million Veteran Program and UK Biobank

Osteoarthritis is a common progressive joint disease. As no effective medical interventions are available, osteoarthritis often progresses to the end stage, in which only surgical options such as total joint replacement are available. A more thorough understanding of genetic influences of osteoarthritis is essential to develop targeted personalized approaches to treatment, ideally long before the end stage is reached. To date, there have been no large multiancestry genetic studies of osteoarthritis. Here, we leveraged the unique resources of 484,374 participants in the Million Veteran Program and UK Biobank to address this gap. Analyses included participants of European, African, Asian and Hispanic descent. We discovered osteoarthritis-associated genetic variation at 10 loci and replicated findings from previous osteoarthritis studies. We also present evidence that some osteoarthritis-associated regions are robust to population ancestry. Drug repurposing analyses revealed enrichment of targets of several medication classes and provide potential insight into the etiology of beneficial effects of antiepileptics on osteoarthritis pain.

Wear patterns in knee OA correlate with native limb geometry

Background: To date, the amount of cartilage loss is graded by means of discrete scoring systems on artificially divided regions of interest (ROI). However, optimal statistical comparison between and within populations requires anatomically standardized cartilage thickness assessment. Providing anatomical standardization relying on non-rigid registration, we aim to compare morphotypes of a healthy control cohort and virtual reconstructed twins of end-stage knee OA subjects to assess the shape-related knee OA risk and to evaluate possible correlations between phenotype and location of cartilage loss.

Methods: Out of an anonymized dataset provided by the Medacta company (Medacta International SA, Castel S. Pietro, CH), 798 end-stage knee OA cases were extracted. Cartilage wear patterns were observed by computing joint space width. The three-dimensional joint space width data was translated into a two-dimensional pixel image, which served as the input for a principal polynomial autoencoder developed for non-linear encoding of wear patterns. Virtual healthy twin reconstruction enabled the investigation of the morphology-related risk for OA requiring joint arthroplasty.

Results: The polynomial autoencoder revealed 4 dominant, orthogonal components, accounting for 94% of variance in the latent feature space. This could be interpreted as medial (54.8%), bicompartmental (25.2%) and lateral (9.1%) wear. Medial wear was subdivided into anteromedial (11.3%) and posteromedial (10.4%) wear. Pre-diseased limb geometry had a positive predictive value of 0.80 in the prediction of OA incidence (r 0.58, p &lt; 0.001).

Conclusion: An innovative methodological workflow is presented to correlate cartilage wear patterns with knee joint phenotype and to assess the distinct knee OA risk based on pre-diseased lower limb morphology. Confirming previous research, both alignment and joint geometry are of importance in knee OA disease onset and progression.

Genetics of osteoarthritis

Osteoarthritis genetics has been transformed in the past decade through the application of large-scale genome-wide association scans. So far, over 100 polymorphic DNA variants have been associated with this common and complex disease. These genetic risk variants account for over 20% of osteoarthritis heritability and the vast majority map to non-protein coding regions of the genome where they are presumed to act by regulating the expression of target genes. Statistical fine mapping, in silico analyses of genomics data, and laboratory-based functional studies have enabled the identification of some of these targets, which encode proteins with diverse roles, including extracellular signaling molecules, intracellular enzymes, transcription factors, and cytoskeletal proteins. A large number of the risk variants correlate with epigenetic factors, in particular cartilage DNA methylation changes in cis, implying that epigenetics may be a conduit through which genetic effects on gene expression are mediated. Some of the variants also appear to have been selected as humans adapted to bipedalism, suggesting that a proportion of osteoarthritis genetic susceptibility results from antagonistic pleiotropy, with risk variants having a positive role in joint formation but a negative role in the long-term health of the joint. Although data from an osteoarthritis genetic study has not yet directly led to a novel treatment, some of the osteoarthritis associated genes code for proteins that have available therapeutics. Genetic investigations are therefore revealing fascinating fundamental insights into osteoarthritis and can expose options for translational intervention.

Associations between adipokines gene polymorphisms and knee osteoarthritis: a meta-analysis

Background

Adipokines gene polymorphisms are speculated to be associated with the risk of knee osteoarthritis (OA), but evidence remains conflicting. This study therefore aimed to examine whether associations exist between adipokines gene polymorphisms and knee OA by considering the evidence collected from eligible studies through a meta-analysis.

Methods

A systematic search was performed on PubMed, Embase, Web of Science, China National Knowledge Infrastructure (CNKI), and Wanfang up to March 31, 2020. Meta-analysis was carried out by focusing on the associations between adipokines gene polymorphisms and knee OA with the allele model, dominant model, and recessive model.

Results

The present meta-analysis included 5 eligible studies for ADIPOQ rs1501299 with 1,021 cases and 1,097 controls, 3 eligible studies for ADIPOQ rs2241766 with 549 cases and 544 controls, 3 eligible studies for LEPR rs1137101 with 808 cases and 856 controls, 2 eligible studies for VISFATIN rs4730153 with 339 cases and 680 controls and 2 eligible studies for VISFATIN rs16872158 with 339 cases and 680 controls. Significant association was observed between LEPR rs1137101 and knee OA in the overall population (recessive: OR = 0.40, 95% CI 0.21–0.79). Limited data revealed that associations may exist between ADIPOQ rs2241766 and knee OA in Asians (dominant: OR = 1.35, 95% CI 1.03–1.78), between VISFATIN rs4730153 and knee OA in Asians (allele: OR = 0.58, 95% CI 0.41–0.83; dominant: OR = 0.57, 95% CI 0.39–0.83), and between VISFATIN rs16872158 and knee OA in Asians (allele: OR = 1.84, 95% CI 1.26–2.68; dominant: OR = 1.94, 95% CI 1.31–2.89).

Conclusions

Adipokines gene polymorphisms may be associated with knee OA. The association was observed in LEPR rs1137101 in the present study. In addition, limited data revealed that associations may also exist in ADIPOQ rs2241766, VISFATIN rs4730153 and VISFATIN rs16872158.

Prospero registration

CRD42020187664.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05111-4.

Inhibition of nuclear receptor RORα attenuates cartilage damage in osteoarthritis by modulating IL-6/STAT3 pathway

Osteoarthritis (OA) is characterized by cartilage destruction, chronic inflammation, and local pain. Evidence showed that retinoic acid receptor-related orphan receptor-α (RORα) is crucial in cartilage development and OA pathogenesis. Here, we investigated the role and molecular mechanism of RORα, an important member of the nuclear receptor family, in regulating the development of OA pathologic features. Investigation into clinical cartilage specimens showed that RORα expression level is positively correlated with the severity of OA and cartilage damage. In an in vivo OA model induced by anterior crucial ligament transaction, intra-articular injection of si-Rora adenovirus reversed the cartilage damage. The expression of cartilage matrix components type II collagen and aggrecan were elevated upon RORα blockade. RNA-seq data suggested that the IL-6/STAT3 pathway is significantly downregulated, manifesting the reduced expression level of both IL-6 and phosphorylated STAT3. RORα exerted its effect on IL-6/STAT3 signaling in two different ways, including interaction with STAT3 and IL-6 promoter. Taken together, our findings indicated the pivotal role of the RORα/IL-6/STAT3 axis in OA progression and confirmed that RORα blockade improved the matrix catabolism in OA chondrocytes. These results may provide a potential treatment target in OA therapy.

The Genetic Epidemiology of Joint Shape and the Development of Osteoarthritis

Congruent, low-friction relative movement between the articulating elements of a synovial joint is an essential pre-requisite for sustained, efficient, function. Where disorders of joint formation or maintenance exist, mechanical overloading and osteoarthritis (OA) follow. The heritable component of OA accounts for ~ 50% of susceptible risk. Although almost 100 genetic risk loci for OA have now been identified, and the epidemiological relationship between joint development, joint shape and osteoarthritis is well established, we still have only a limited understanding of the contribution that genetic variation makes to joint shape and how this modulates OA risk. In this article, a brief overview of synovial joint development and its genetic regulation is followed by a review of current knowledge on the genetic epidemiology of established joint shape disorders and common shape variation. A summary of current genetic epidemiology of OA is also given, together with current evidence on the genetic overlap between shape variation and OA. Finally, the established genetic risk loci for both joint shape and osteoarthritis are discussed.

Mean femoral cartilage thickness is higher in athletes as compared with sedentary individuals

PURPOSE

It is unclear how high-intensity physical activity (HIPA) affects the knee joint, specifically the femoral cartilage (FC). Therefore, the aims of this study were to evaluate FC thickness via ultrasound among elite athletes involved in different types of HIPA, and to determine whether there is a correlation between serum cartilage oligomeric matrix protein (sCOMP) and rectus femoris (RF) thickness.

METHOD

A total of 132 male individuals participated in this study and were assigned to two groups, the sedentary (n = 43, 23.9 ± 3.7) and athlete groups (n = 89, 22.7 ± 4.6), which did not significantly differ in age. The athletes were elite and performed HIPA during sports such as volleyball (n = 20), soccer (n = 21), basketball (n = 28), and weightlifting (n = 20). RF thickness and three (mid-point) measurements were obtained for each knee. The mean FC thickness for each knee was defined as the sum of the medial, lateral condyles, and intercondylar areas. Blood samples for sCOMP analyses were also obtained.

RESULTS

All the measurements of the FC of both knees were significantly higher in the athletes than in the sedentary individuals (p < 0.001 and p = 0.001). The mean right and left FC values were also higher in the athletes (p < 0.001). Multiple linear regression analysis showed that participation in sporting activities was a significant predictor associated with the right and left mean FC thickness (p < 0.001 for both). No significant differences in the sCOMP levels were found between the two groups.

CONCLUSION

It was found that the mean FC was higher among athletes than among sedentary individuals. As a result, it is suggested that sports' participation is an independent factor associated with the right and left mean FC thickness.

LEVEL OF EVIDENCE

III.

Genetic and Epigenetic Interplay Within a COLGALT2 Enhancer Associated With Osteoarthritis

OBJECTIVE

The osteoarthritis (OA)-associated single-nucleotide polymorphism (SNP) rs11583641 is located in COLGALT2, encoding a posttranslational modifier of collagen. In cartilage, the SNP genotype correlates with DNA methylation in a putative enhancer. This study was undertaken to characterize the mechanistic relationship between rs11583641, the putative enhancer, and COLGALT2 expression using cartilage samples from human patients and a chondrocyte cell model.

METHODS

Nucleic acids were extracted from articular cartilage samples obtained from patients with OA (n = 137). Samples were genotyped, and DNA methylation was quantified at 12 CpGs using pyrosequencing. The putative enhancer was deleted in Tc28a2 chondrocytes using clustered regularly interspaced short palindromic repeat/Cas9, and the impact on nearby gene expression was determined using real-time quantitative polymerase chain reaction. Targeted modulation of the epigenome using catalytically dead Cas9 (dCas9) constructs fused to DNA methyltransferase 3a or ten-eleven translocase 1 allowed for the investigation of a causal relationship between DNA methylation and enhancer activity.

RESULTS

The genotype at rs11583641 correlated with DNA methylation at 3 CpGs, and the presence of the OA risk allele, C, corresponded to reduced levels of methylation. Deletion of the enhancer resulted in a 2.7-fold reduction in COLGALT2 expression. Targeted methylation and demethylation of the CpGs had antagonistic effects on COLGALT2 expression. An allelic imbalance in the expression of COLGALT2 was identified in the cartilage from patients with OA, with relative overexpression of the OA risk allele. Allelic expression ratios correlated with DNA methylation at 4 CpGs.

CONCLUSION

COLGALT2 is a target of OA genetic risk at this locus. The genotype at rs11583641 impacts DNA methylation in a gene enhancer, which, in turn, modulates COLGALT2 expression. COLGALT2 encodes an enzyme that initiates posttranslational glycosylation of collagens and is therefore a compelling OA susceptibility target.

NCOA3 identified as a new candidate to explain autosomal dominant progressive hearing loss

Hearing loss is a frequent sensory impairment in humans and genetic factors account for an elevated fraction of the cases. We have investigated a large family of five generations, with 15 reported individuals presenting non-syndromic, sensorineural, bilateral and progressive hearing loss, segregating as an autosomal dominant condition. Linkage analysis, using SNP-array and selected microsatellites, identified a region of near 13 cM in chromosome 20 as the best candidate to harbour the causative mutation. After exome sequencing and filtering of variants, only one predicted deleterious variant in the NCOA3 gene (NM_181659, c.2810C > G; p.Ser937Cys) fit in with our linkage data. RT-PCR, immunostaining and in situ hybridization showed expression of ncoa3 in the inner ear of mice and zebrafish. We generated a stable homozygous zebrafish mutant line using the CRISPR/Cas9 system. ncoa3-/- did not display any major morphological abnormalities in the ear, however, anterior macular hair cells showed altered orientation. Surprisingly, chondrocytes forming the ear cartilage showed abnormal behaviour in ncoa3-/-, detaching from their location, invading the ear canal and blocking the cristae. Adult mutants displayed accumulation of denser material wrapping the otoliths of ncoa3-/- and increased bone mineral density. Altered zebrafish swimming behaviour corroborates a potential role of ncoa3 in hearing loss. In conclusion, we identified a potential candidate gene to explain hereditary hearing loss, and our functional analyses suggest subtle and abnormal skeletal behaviour as mechanisms involved in the pathogenesis of progressive sensory function impairment.

Accelerating functional gene discovery in osteoarthritis

Osteoarthritis causes debilitating pain and disability, resulting in a considerable socioeconomic burden, yet no drugs are available that prevent disease onset or progression. Here, we develop, validate and use rapid-throughput imaging techniques to identify abnormal joint phenotypes in randomly selected mutant mice generated by the International Knockout Mouse Consortium. We identify 14 genes with functional involvement in osteoarthritis pathogenesis, including the homeobox gene Pitx1, and functionally characterize 6 candidate human osteoarthritis genes in mouse models. We demonstrate sensitivity of the methods by identifying age-related degenerative joint damage in wild-type mice. Finally, we phenotype previously generated mutant mice with an osteoarthritis-associated polymorphism in the Dio2 gene by CRISPR/Cas9 genome editing and demonstrate a protective role in disease onset with public health implications. We hope this expanding resource of mutant mice will accelerate functional gene discovery in osteoarthritis and offer drug discovery opportunities for this common, incapacitating chronic disease.

Widespread regulation of gene expression by glucocorticoids in chondrocytes from patients with osteoarthritis as determined by RNA-Seq

Background

Intra-articular glucocorticoid (GC) injections are widely used as a symptomatic treatment for osteoarthritis (OA). However, there are also concerns about their potentially harmful effects, and their detailed effects on chondrocyte phenotype remain poorly understood.

Methods

We studied the effects of dexamethasone on gene expression in OA chondrocytes with RNA-Seq. Chondrocytes were isolated from the cartilage from OA patients undergoing knee replacement surgery and cultured with or without dexamethasone for 24 h. Total RNA was isolated and sequenced, and functional analysis was performed against the Gene Ontology (GO) database. Results for selected genes were confirmed with RT-PCR. We also investigated genes linked to OA in recent genome-wide expression analysis (GWEA) studies.

Results

Dexamethasone increased the expression of 480 and reduced that of 755 genes with a fold change (FC) 2.0 or greater. Several genes associated with inflammation and cartilage anabolism/catabolism as well as lipid and carbohydrate metabolism were among the most strongly affected genes. In the GO analysis, genes involved in the extracellular matrix organization, cell proliferation and adhesion, inflammation, and collagen synthesis were enriched among the significantly affected genes. In network analysis, NGF, PI3KR1, and VCAM1 were identified as central genes among those most strongly affected by dexamethasone.

Conclusions

This is the first study investigating the genome-wide effects of GCs on the gene expression in OA chondrocytes. In addition to clear anti-inflammatory and anticatabolic effects, GCs affect lipid and glucose metabolism in chondrocytes, an observation that might be particularly important in the metabolic phenotype of OA.

Gene-based association tests using GWAS summary statistics

MOTIVATION

A huge number of genome-wide association studies (GWAS) summary statistics freely available in databases provide a new material for gene-based association analysis aimed at identifying rare genetic variants. Only a few of the many popular gene-based methods developed for individual genotype and phenotype data are adapted for the practical use of the GWAS summary statistics as input.

RESULTS

We analytically prove and numerically illustrate that all popular powerful methods developed for gene-based association analysis of individual phenotype and genotype data can be modified to utilize GWAS summary statistics. We have modified and implemented all of the popular methods, including burden and kernel machine-based tests, multiple and functional linear regression, principal components analysis and others, in the R package sumFREGAT. Using real summary statistics for coronary artery disease, we show that the new package is able to detect genes not found by the existing packages.

AVAILABILITY AND IMPLEMENTATION

The R package sumFREGAT is freely and publicly available at: https://CRAN.R-project.org/package=sumFREGAT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolutionary Selection and Constraint on Human Knee Chondrocyte Regulation Impacts Osteoarthritis Risk

During human evolution, the knee adapted to the biomechanical demands of bipedalism by altering chondrocyte developmental programs. This adaptive process was likely not without deleterious consequences to health. Today, osteoarthritis occurs in 250 million people, with risk variants enriched in non-coding sequences near chondrocyte genes, loci that likely became optimized during knee evolution. We explore this relationship by epigenetically profiling joint chondrocytes, revealing ancient selection and recent constraint and drift on knee regulatory elements, which also overlap osteoarthritis variants that contribute to disease heritability by tending to modify constrained functional sequence. We propose a model whereby genetic violations to regulatory constraint, tolerated during knee development, lead to adult pathology. In support, we discover a causal enhancer variant (rs6060369) present in billions of people at a risk locus (GDF5-UQCC1), showing how it impacts mouse knee-shape and osteoarthritis. Overall, our methods link an evolutionarily novel aspect of human anatomy to its pathogenesis.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics

Although osteoarthritis (OA) aetiology is complex, genetic, genomic and epigenetic studies published within the last decade have advanced our understanding of the molecular processes underlying this common musculoskeletal disease. The purpose of this narrative review is to highlight the key research articles within the OA genetics, genomics and epigenetics fields that were published between April 2018 and April 2019. The review focuses on the identification of new OA genetic risk loci, genomics techniques that have been used for the first time in human cartilage and new publicly available databases, and datasets that will aid OA functional studies. Fifty-six new OA susceptibility loci were identified by two large scale genome wide association study meta-analyses, increasing the number of genome-wide significant risk loci to 90. OA risk variants are enriched near genes involved in skeletal development and morphology, and show genetic overlap with height, hip shape, bone area and developmental dysplasia of the hip. Several functional studies of OA loci were published, including a genome-wide analysis of genetic variation on cartilage gene expression. A specialised data portal for exploring cross-species skeletal transcriptomic datasets has been developed, and the first use of cartilage single cell RNAseq analysis reported. This year also saw the systematic identification of all microRNAs, long non-coding RNAs and circular RNAs expressed in human OA cartilage. Putative transcriptional regulatory regions have been mapped in human chondrocytes genome-wide, providing a dataset that will facilitate the prioritisation and characterisation of OA genetic and epigenetic loci.

Differences between race and sex in measures of hip morphology: a population-based comparative study

OBJECTIVE

This paper aims to (i) identify differences in measures of hip morphology between four racial groups using anteroposterior (AP) hip x-rays, and (ii) examine whether these differences vary by sex.

METHODS

912 hip x-rays (456 individuals) from four racial groups (European Caucasians, American Caucasians, African Americans and Chinese) were obtained. Males and females (45-75 years) with no radiographic hip OA (Kellgren and Lawrence < Grade 2 or Croft < Grade 1) were included. Eleven features of hip joint morphology were analysed. Linear regression with generalised estimating equations (GEE) was used to determine race and sex differences in hip morphology. Post-hoc Bonferroni procedure was used to adjust for multiple comparisons.

RESULTS

The final analysis included 875 hips. Chinese hips showed significant differences for the majority of measures to other racial groups. Chinese were characterised by more shallow and narrow acetabular sockets, reduced femoral head coverage, smaller femoral head diameter, and a lesser angle of alignment between the femoral neck and shaft. Variation was found between other racial groups, but with few statistically significant differences. The average of lateral centre edge angle, minimum neck width and neck length differed between race and sex (p-value for interaction < 0.05).

CONCLUSIONS

Significant differences were found in measures of morphology between Chinese hips compared to African Americans or Caucasian groups; these may explain variation in hip OA prevalence rates between these groups and the lower rate of hip OA in Chinese. Sex differences were also identified, which may further explain male-female prevalence differences for OA.

Inhibition of 4-aminobutyrate aminotransferase protects against injury-induced osteoarthritis in mice

Recently we demonstrated that ablation of the DNA methyltransferase enzyme, Dnmt3b, resulted in catabolism and progression of osteoarthritis (OA) in murine articular cartilage through a mechanism involving increased mitochondrial respiration. In this study, we identify 4-aminobutyrate aminotransferase (Abat) as a downstream target of Dnmt3b. Abat is an enzyme that metabolizes γ-aminobutyric acid to succinate, a key intermediate in the tricarboxylic acid cycle. We show that Dnmt3b binds to the Abat promoter, increases methylation of a conserved CpG sequence just upstream of the transcriptional start site, and inhibits Abat expression. Dnmt3b deletion in articular chondrocytes results in reduced methylation of the CpG sequence in the Abat promoter, which subsequently increases expression of Abat. Increased Abat expression in chondrocytes leads to enhanced mitochondrial respiration and elevated expression of catabolic genes. Overexpression of Abat in murine knee joints via lentiviral injection results in accelerated cartilage degradation following surgical induction of OA. In contrast, lentiviral-based knockdown of Abat attenuates the expression of IL-1β-induced catabolic genes in primary murine articular chondrocytes in vitro and also protects against murine articular cartilage degradation in vivo. Strikingly, treatment with the FDA-approved small-molecule Abat inhibitor, vigabatrin, significantly prevents the development of injury-induced OA in mice. In summary, these studies establish Abat as an important new target for therapies to prevent OA.

Inhibition of cyclooxygenase-2 activity in subchondral bone modifies a subtype of osteoarthritis

Osteoarthritis (OA) causes the destruction of joints. Its pathogenesis is still under investigation, and there is no effective disease-modifying therapy. Here, we report that elevated cyclooxygenase-2 (COX-2) expression in the osteocytes of subchondral bone causes both spontaneous OA and rheumatoid arthritis (RA). The knockout of COX-2 in osteocytes or treatment with a COX-2 inhibitor effectively rescues the structure of subchondral bone and attenuates cartilage degeneration in spontaneous OA (STR/Ort) mice and tumor necrosis factor-α transgenic RA mice. Thus, elevated COX-2 expression in subchondral bone induces both OA-associated and RA-associated joint cartilage degeneration. The inhibition of COX-2 expression can potentially modify joint destruction in patients with arthritis.

Cartilage and Bone Destruction in Arthritis: Pathogenesis and Treatment Strategy: A Literature Review

Arthritis is inflammation of the joints accompanied by osteochondral destruction. It can take many forms, including osteoarthritis, rheumatoid arthritis, and psoriatic arthritis. These diseases share one commonality—osteochondral destruction based on inflammation. The background includes a close interaction between osseous tissues and immune cells through various inflammatory cytokines. However, the tissues and cytokines that play major roles are different in each disease, and as a result, the mechanism of osteochondral destruction also differs. In recent years, there have been many findings regarding not only extracellular signaling pathways but also intracellular signaling pathways. In particular, we anticipate that the intracellular signals of osteoclasts, which play a central role in bone destruction, will become novel therapeutic targets. In this review, we have summarized the pathology of arthritis and the latest findings on the mechanism of osteochondral destruction, as well as present and future therapeutic strategies for these targets.

Genome Engineering for Osteoarthritis: From Designer Cells to Disease-Modifying Drugs

BACKGROUND

Osteoarthritis (OA) is a highly prevalent degenerative joint disease involving joint cartilage and its surrounding tissues. OA is the leading cause of pain and disability worldwide. At present, there are no disease-modifying OA drugs, and the primary therapies include exercise and nonsteroidal anti-inflammatory drugs until total joint replacement at the end-stage of the disease.

METHODS

In this review, we summarized the current state of knowledge in genetic and epigenetic associations and risk factors for OA and their potential diagnostic and therapeutic applications.

RESULTS

Genome-wide association studies and analysis of epigenetic modifications (such as miRNA expression, DNA methylation and histone modifications) conducted across various populations support the notion that there is a genetic basis for certain subsets of OA pathogenesis.

CONCLUSION

With recent advances in the development of genome editing technologies such as the CRISPR-Cas9 system, these genetic and epigenetic alternations in OA can be used as platforms from which potential biomarkers for the diagnosis, prognosis, drug response, and development of potential personalized therapeutic targets for OA can be approached. Furthermore, genome editing has allowed the development of "designer" cells, whereby the receptors, gene regulatory networks, or transgenes can be modified as a basis for new cell-based therapies.

Functional testing of thousands of osteoarthritis-associated variants for regulatory activity

To date, genome-wide association studies have implicated at least 35 loci in osteoarthritis but, due to linkage disequilibrium, the specific variants underlying these associations and the mechanisms by which they contribute to disease risk have yet to be pinpointed. Here, we functionally test 1,605 single nucleotide variants associated with osteoarthritis for regulatory activity using a massively parallel reporter assay. We identify six single nucleotide polymorphisms (SNPs) with differential regulatory activity between the major and minor alleles. We show that the most significant SNP, rs4730222, exhibits differential nuclear protein binding in electrophoretic mobility shift assays and drives increased expression of an alternative isoform of HBP1 in a heterozygote chondrosarcoma cell line, in a CRISPR-edited osteosarcoma cell line, and in chondrocytes derived from osteoarthritis patients. This study provides a framework for prioritization of GWAS variants and highlights a role of HBP1 and Wnt signaling in osteoarthritis pathogenesis.

Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage

OBJECTIVE

Multiple single-nucleotide polymorphisms (SNPs) conferring susceptibility to osteoarthritis (OA) mark imbalanced expression of positional genes in articular cartilage, reflected by unequally expressed alleles among heterozygotes (allelic imbalance [AI]). We undertook this study to explore the articular cartilage transcriptome from OA patients for AI events to identify putative disease-driving genetic variation.

METHODS

AI was assessed in 42 preserved and 5 lesioned OA cartilage samples (from the Research Arthritis and Articular Cartilage study) for which RNA sequencing data were available. The count fraction of the alternative alleles among the alternative and reference alleles together (φ) was determined for heterozygous individuals. A meta-analysis was performed to generate a meta-φ and P value for each SNP with a false discovery rate (FDR) correction for multiple comparisons. To further validate AI events, we explored them as a function of multiple additional OA features.

RESULTS

We observed a total of 2,070 SNPs that consistently marked AI of 1,031 unique genes in articular cartilage. Of these genes, 49 were found to be significantly differentially expressed (fold change <0.5 or >2, FDR <0.05) between preserved and paired lesioned cartilage, and 18 had previously been reported to confer susceptibility to OA and/or related phenotypes. Moreover, we identified notable highly significant AI SNPs in the CRLF1, WWP2, and RPS3 genes that were related to multiple OA features.

CONCLUSION

We present a framework and resulting data set for researchers in the OA research field to probe for disease-relevant genetic variation that affects gene expression in pivotal disease-affected tissue. This likely includes putative novel compelling OA risk genes such as CRLF1, WWP2, and RPS3.

Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data

Osteoarthritis is the most common musculoskeletal disease and the leading cause of disability globally. Here, we perform a genome-wide association study for osteoarthritis (77,052 cases and 378,169 controls), analysing 4 phenotypes: knee osteoarthritis, hip osteoarthritis, knee and/or hip osteoarthritis, and any osteoarthritis. We discover 64 signals, 52 of them novel, more than doubling the number of established disease loci. Six signals fine map to a single variant. We identify putative effector genes by integrating eQTL colocalization, fine-mapping, human rare disease, animal model, and osteoarthritis tissue expression data. We find enrichment for genes underlying monogenic forms of bone development diseases, and for the collagen formation and extracellular matrix organisation biological pathways. Ten of the likely effector genes, including TGFB1, FGF18, CTSK and IL11 have therapeutics approved or in clinical trials, with mechanisms of action supportive of evaluation for efficacy in osteoarthritis.

Elucidating the environmental risk factors for rheumatic diseases: An umbrella review of meta-analyses

AIMS

Although rheumatic diseases constitute a leading cause of disability, the environmental risk factors for these diseases are not clarified. In the present study, we aim to systematically appraise the epidemiological credibility of the environmental risk factors for rheumatic diseases.

METHODS

We systematically searched PubMed to capture meta-analyses of observational studies on environmental risk factors for the most prevalent rheumatic diseases. For each association, we estimated the summary effect size estimate, the 95% confidence and prediction intervals, and the I² metric. We further examined the presence of small-study effects and excess significance bias.

RESULTS

Overall, we identified 30 eligible papers describing 42 associations. Thirty-three associations were statistically significant at P < 0.05, whereas 13 of them were statistically significant at P < 1 × 10^-6 . Thirty-two associations had large or very large between-study heterogeneity. In 12 associations, evidence of small-study effects and/or excess significance bias was found. Six risk factors (nine associations) presented convincing or highly suggestive evidence of association: smoking and pack-years of smoking for rheumatoid arthritis; BMI (per 5 kg/m² increase) for gout and hip osteoarthritis; alcohol consumption for gout; BMI (overweight vs lean, obese vs lean), knee injury and participation in heavy work for knee osteoarthritis.

CONCLUSION

Our umbrella review indicated that a narrow range of risk factors has been examined for rheumatic diseases. Current evidence strongly supports that smoking, obesity, alcohol consumption, knee injury, and work activities are associated with risk for at least one rheumatic disease.

Mitochondrial DNA variants and pulmonary function in older persons

BACKGROUND

We provide the first examination of mitochondrial DNA (mtDNA) variants and pulmonary function in older persons.

METHODS

Cross-sectional associations between mtDNA variants and pulmonary function were evaluated as a combined p-values meta-analysis, using data from two independent cohorts of older persons. The latter included white and black participants, aged ≥70 years, from the Lifestyle Interventions and Independence for Elders study (LIFE) (N = 1247) and the Health, Aging and Body Composition study (Health ABC) (N = 731), respectively. Pulmonary function included the forced expiratory volume in one-second as a Z-score (FEV1z) and the maximal inspiratory pressure (MIP) in cm of water.

RESULTS

In black participants, significant associations were found between mtDNA variants and MIP: m.7146A > G, COI (p = 3E-5); m.7389 T > C, COI (p = 2E-4); m.15301G > A, CYB (p = 9E-5); m.16265A > G, HV1 (p = 9E-5); meta-analytical p-values <0.0002. Importantly, these mtDNA variants were unique to black participants and were not present in white participants. Moreover, in black participants, aggregate genetic effects on MIP were observed across mutations in oxidative phosphorylation complex IV (p = 0.004), complex V (p = 0.0007), and hypervariable (p = 0.003) regions. The individual and aggregate variant results were significant after adjustment for multiple comparisons. Otherwise, no significant associations were detected for MIP in whites or for FEV1z in whites or blacks.

CONCLUSIONS

We have shown that mtDNA variants of African origin are cross-sectionally associated with MIP, a measure of respiratory muscle strength. Thus, our results establish the rationale for longitudinal studies to evaluate whether mtDNA variants of African origin identify those at risk of subsequently developing a respiratory muscle impairment (lower MIP values).

Meta-analysis of Icelandic and UK data sets identifies missense variants in SMO, IL11, COL11A1 and 13 more new loci associated with osteoarthritis

Osteoarthritis has a highly negative impact on quality of life because of the associated pain and loss of joint function. Here we describe the largest meta-analysis so far of osteoarthritis of the hip and the knee in samples from Iceland and the UK Biobank (including 17,151 hip osteoarthritis patients, 23,877 knee osteoarthritis patients, and more than 562,000 controls). We found 23 independent associations at 22 loci in the additive meta-analyses, of which 16 of the loci were novel: 12 for hip and 4 for knee osteoarthritis. Two associations are between rare or low-frequency missense variants and hip osteoarthritis, affecting the genes SMO (rs143083812, frequency 0.11%, odds ratio (OR) = 2.8, P = 7.9 × 10^-12, p.Arg173Cys) and IL11 (rs4252548, frequency 2.08%, OR = 1.30, P = 2.1 × 10^-11, p.Arg112His). A common missense variant in the COL11A1 gene also associates with hip osteoarthritis (rs3753841, frequency 61%, P = 5.2 × 10^-10, OR = 1.08, p.Pro1284Leu). In addition, using a recessive model, we confirm an association between hip osteoarthritis and a variant of CHADL¹ (rs117018441, P = 1.8 × 10^-25, OR = 5.9). Furthermore, we observe a complex relationship between height and risk of osteoarthritis.

Chromatin accessibility landscape of articular knee cartilage reveals aberrant enhancer regulation in osteoarthritis

Osteoarthritis (OA) is a common joint disorder with increasing impact in an aging society. While genetic and transcriptomic analyses have revealed some genes and non-coding loci associated to OA, the pathogenesis remains incompletely understood. Chromatin profiling, which provides insight into gene regulation, has not been reported in OA mainly due to technical difficulties. Here, we employed Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) to map the accessible chromatin landscape in articular knee cartilage of OA patients. We identified 109,215 accessible chromatin regions for cartilages, of which 71% were annotated as enhancers. By overlaying them with genetic and DNA methylation data, we have determined potential OA-relevant enhancers and their putative target genes. Furthermore, through integration with RNA-seq data, we characterized genes that are altered both at epigenomic and transcriptomic levels in OA. These genes are enriched in pathways regulating ossification and mesenchymal stem cell (MSC) differentiation. Consistently, the differentially accessible regions in OA are enriched for MSC-specific enhancers and motifs of transcription factor families involved in osteoblast differentiation. In conclusion, we demonstrate how direct chromatin profiling of clinical tissues can provide comprehensive epigenetic information for a disease and suggest candidate genes and enhancers of translational potential.

Meta-analysis identifies mitochondrial DNA sequence variants associated with walking speed

Declines in walking speed are associated with a variety of poor health outcomes including disability, comorbidity, and mortality. While genetic factors are putative contributors to variability in walking, few genetic loci have been identified for this trait. We examined the role of mitochondrial genomic variation on walking speed by sequencing the entire mitochondrial DNA (mtDNA). Data were meta-analyzed from 1758 Lifestyle Interventions and Independence for Elders (LIFE) Study and replication data from 730 Health, Aging, and Body Composition (HABC) Study participants with baseline walking speed information. Participants were 69+ years old of diverse racial backgrounds (African, European, and other race/ethnic groups) and had a wide range of mean walking speeds [4-6 m (0.78-1.09 m/s) and 400 m (0.83-1.24 m/s)]. Meta-analysis across studies and racial groups showed that m.12705C>T, ND5 variant was significantly associated (p < 0.0001) with walking speed at both short and long distances. Replication and meta-analysis also identified statistically significant walking speed associations (p < 0.0001) between the m.5460.G>A, ND2 and m.309C>CT, HV2 variants at short and long distances, respectively. All results remained statistically significant after multiple comparisons adjustment for 499 mtDNA variants. The m.12705C>T variant can be traced to the beginnings of human global migration and that cells carrying this variant display altered tRNA expression. Significant pooled effects related to stopping during the long-distance walk test were observed across OXPHOS complexes I (p = 0.0017) and III (p = 0.0048). These results suggest that mtDNA-encoded variants are associated with differences in walking speed among older adults, potentially identifying those at risk of developing mobility impairments.

Single Nucleotide Polymorphism in the COL11A2 Gene Associated with Heat Pain Sensitivity in Knee Osteoarthritis

Pain is one of the most prominent symptoms of osteoarthritis. However, there is often discordance between the pain experienced by individuals with osteoarthritis and the degree of articular pathology. This suggests that individual differences, including genetic variability in the central processing of nociceptive stimuli, may impact the presentation of osteoarthritis. Here, we show that the single nucleotide polymorphism rs16868943 in the collagen gene COL11A2 is significantly associated with lowered heat pain tolerance on the arm in participants with knee osteoarthritis (P = 1.21 × 10⁻⁶, P = 0.0053 after Bonferroni correction, beta = −3.42). A total of 161 knee osteoarthritis participants were included and evaluated for heat, punctate and pressure pain sensitivity of the affected knee and the ipsilateral arm. Each participant was genotyped for 4392 single nucleotide polymorphisms in genes implicated in pain perception, inflammation and mood and tested for association with pain sensitivity. The minor A allele of single nucleotide polymorphism rs16868943 was significantly associated with lower arm heat pain tolerance after correction for age, gender, race, and study site. This single nucleotide polymorphism was also nominally associated with other measures of heat pain sensitivity, including lowered knee heat pain tolerance (P = 1.14 × 10⁻⁵, P = 0.05 after Bonferroni correction), lowered arm heat pain threshold (P = 0.0039, uncorrected) and lowered knee heat pain threshold (P = 0.003, uncorrected). Addition of genotypes from 91 participants without knee pain produced a significant interaction between knee osteoarthritis status and the rs16868943 single nucleotide polymorphism in heat pain tolerance (P = 1.71 × 10⁻⁵), such that rs16868943 was not associated with heat pain tolerance in participants without knee pain (P = 0.12, beta = 1.3). This is the first study to show genetic association with heat pain tolerance in individuals with osteoarthritis. The association is specific to participants who have already developed knee osteoarthritis, suggesting that the COL11A2 gene, which has previously been associated with familial osteoarthritis, may play a role in pain sensitization after the development of osteoarthritis.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

The Involvement of MicroRNAs in Osteoarthritis and Recent Developments: A Narrative Review

Background:

Osteoarthritis (OA) is the most common chronic joint disease and it may progressively cause disability and compromise quality of life. Lately, the role of miRNAs in the pathogenesis of OA has drawn a lot of attention. miRNAs are small, single-stranded, non-coding molecules of RNA which regulate gene expression at post-transcriptional level. The dysregulation of the expression of several miRNAs affects pathways involved in OA pathogenesis.

Objective:

The purpose of this article is to review the literature on the involvement of miRNAs in the pathogenesis of OA and the implications on its diagnosis and treatment.

Materials and Methods:

An extensive electronic literature search was conducted by two researchers from January 2008 to August 2017. Titles and abstracts of papers were screened by the authors for further inclusion in the present work. Finally, full texts of the selected articles were retrieved.

Results:

Abnormally expressed miRNAs enhance the production of cartilage degrading enzymes, inhibit the expression of cartilage matrix components, increase the production of proinflammatory cytokines, facilitate chondrocyte apoptosis, suppress autophagy in chondrocytes and are involved in pain-related pathways. miRNAs are also incorporated in extra-cellular membranous vesicles such as exosomes and participate in the intercellular communication in osteoarthritic joints.

Conclusion:

Ongoing research on miRNAs has potential implications in the diagnosis and treatment of OA. Their different levels in peripheral blood and synovial fluid between OA patients and healthy population makes them candidates for being used as biomarkers of the disease, while targeting miRNAs may be a novel therapeutic strategy in OA.

What did we learn from 'omics' studies in osteoarthritis

PURPOSE OF REVIEW

'Omics' technologies developed for the massive analysis of the major biologically relevant molecules (genes, proteins, metabolites) have been applied to the study of osteoarthritis (OA) for more than a decade.

RECENT FINDINGS

'Omics' studies have undoubtedly contributed to increase the knowledge on pathogenic processes related with OA and have provided hundreds to thousands of molecules that might have a putative biomarker utility for this disease.

SUMMARY

This review describes the most recent 'omics' studies in OA research, their conclusions, and discuss those remaining challenges. Still many validation studies must be performed in large and well-characterized cohorts for the translation of the findings from 'omics' strategies to clinical applications. The development of tools for the intelligent integration of 'omics' data with clinical and imaging information is also mandatory to take full profit of the work that has been already performed.

Association between Fas/FasL gene polymorphism and musculoskeletal degenerative diseases: a meta-analysis

Background

It was reported that Fas (rs1800682, rs2234767) and FasL (rs5030772, rs763110) gene polymorphism might be related to the risk of musculoskeletal degenerative diseases (MSDD), such as osteoarthritis (OA), intervertebral disc degeneration (IVDD) and rheumatoid arthritis (RA). However, data from different studies was inconsistent. Here we aim to elaborately summarize and explore the association between the Fas (rs1800682, rs2234767) and FasL (rs5030772, rs763110) and MSDD.

Methods

Literatures were selected from PubMed, Web of Science, Embase, Scopus and Medline in English and VIP, SinoMed, Wanfang and the China National Knowledge Infrastructure (CNKI) in Chinese up to August 21, 2017. All the researches included are case-control studies about human. We calculated the pooled odds ratios (ORs) with 95% confidence intervals (95% CI) to evaluate the strengths of the associations of Fas (rs1800682, rs2234767) and FasL (rs5030772, rs763110) polymorphisms with MSDD risk.

Results

Eleven eligible studies for rs1800682 with 1930 cases and 1720 controls, 6 eligible studies for rs2234767 with 1794 cases and 1909 controls, 3 eligible studies for rs5030772 with 367 cases and 313 controls and 8 eligible studies for rs763110 with 2010 cases and 2105 controls were included in this analysis. The results showed that the G allele of Fas (rs1800682) is associated with an increased risk of IVDD in homozygote and recessive models. The G allele of Fas (rs2234767) is linked to a decreased risk of RA but an enhanced risk of OA in allele and recessive models. In addition, the T allele of FasL (rs763110) is correlated with a reduced risk of IVDD in all of models. However, no relationship was found between FasL (rs5030772) and these three types of MSDD in any models.

Conclusions

Fas (rs1800682) and FasL (rs763110) polymorphism were associated with the risk of IVDD and Fas (rs2234767) was correlated to the susceptibility of OA and RA. Fas (rs1800682) and Fas (rs2234767) are more likely to be associated with MSDD for Chinese people. FasL (rs763110) is related to the progression of MSDD for both Caucasoid and Chinese race groups. But FasL (rs5030772) might not be associated with any types of MSDD or any race groups statistically.

Electronic supplementary material

The online version of this article (10.1186/s12891-018-2057-z) contains supplementary material, which is available to authorized users.

Association between EN1 rs4144782 and susceptibility of knee osteoarthritis: A case-control study

Osteoarthritis (OA) is a complex disease that affects the whole joint, resulting from the combined influence of biomechanical factors and genetic factors. The heritable component for primary OA accounts for about 60% of variation in population liability to the disease. So far, genome-wide association studies (GWAS) and candidate gene studies have established many OA-related loci. However, these findings account for only a rather small fraction of the genetic component. To further reveal the genetic architecture of OA, we conducted this case-control study to explore the association of locus EN1 rs4144782 and knee OA susceptibility in a Chinese population. EN1 rs4144782 was significantly associated with increased risk of knee OA (OR=1.26; 95% CI: 1.05-1.50, P value=0.012). In dominant model, compared with carriers of GG genotype, those with AG or AA genotype have an 1.44-fold increased risk of OA (OR: 1.44; 95% CI: 1.10-1.88; P value=0.008). Subgroup analyses didn't change the results materially. This should be the first association study of EN1 locus on risk of OA, and our finding suggested that the EN1 rs4144782 might contribute to the susceptibility of knee OA.

The association of AKNA gene polymorphisms with knee osteoarthritis suggests the relevance of this immune response regulator in the disease genetic susceptibility

Recent studies have identified AKNA as a potential susceptibility gene for several inflammatory diseases. Here, we aimed to assess the potential association of AKNA polymorphisms with knee osteoarthritis (KOA) susceptibility in a Mexican population, following STREGA recommendations. From a DNA bank of 181 KOA patients and 140 healthy controls, two AKNA SNPs were genotyped using TaqMan probes. The association between KOA susceptibility and AKNA polymorphisms genotypes was evaluated by multivariated logistic regression analysis. Information regarding patients' inflammatory biomarkers levels was obtained and their association with AKNA polymorphisms genotypes was assessed by lineal regression. We found a positive association with the recessive inheritance model of both AKNA polymorphisms (A/A genotype for both) and KOA susceptibility adjusting by age, body mass index (BMI), gender and place of birth (OR = 2.48, 95% CI 1.09-5.65 for rs10817595 polymorphism; and OR = 4.96; 95% CI 2.421-10.2 for rs3748176 polymorphism). Additionally these associations were also seen after stratifying patients by KOA severity and age. Furthermore the total leukocyte count was positively associated with rs10817595 AKNA polymorphism (β = 1.39; 95% CI 0.44-2.34) adjusting by age, BMI, gender, place of birth and disease severity. We suggest that regulatory and coding polymorphisms of the inflammatory modulator gene AKNA can influence the development of KOA. Further structural and functional studies might reveal the role of AKNA in OA and other rheumatic diseases.

Genome-wide association and functional studies identify a role for matrix Gla protein in osteoarthritis of the hand

OBJECTIVE

Osteoarthritis (OA) is the most common form of arthritis and the leading cause of disability in the elderly. Of all the joints, genetic predisposition is strongest for OA of the hand; however, only few genetic risk loci for hand OA have been identified. Our aim was to identify novel genes associated with hand OA and examine the underlying mechanism.

METHODS

We performed a genome-wide association study of a quantitative measure of hand OA in 12 784 individuals (discovery: 8743, replication: 4011). Genome-wide significant signals were followed up by analysing gene and allele-specific expression in a RNA sequencing dataset (n=96) of human articular cartilage.

RESULTS

We found two significantly associated loci in the discovery set: at chr12 (p=3.5 × 10-10) near the matrix Gla protein (MGP) gene and at chr12 (p=6.1×10-9) near the CCDC91 gene. The DNA variant near the MGP gene was validated in three additional studies, which resulted in a highly significant association between the MGP variant and hand OA (rs4764133, Betameta=0.83, Pmeta=1.8*10-15). This variant is high linkage disequilibrium with a coding variant in MGP, a vitamin K-dependent inhibitor of cartilage calcification. Using RNA sequencing data from human primary cartilage tissue (n=96), we observed that the MGP RNA expression of the hand OA risk allele was significantly lowercompared with the MGP RNA expression of the reference allele (40.7%, p<5*10-16).

CONCLUSIONS

Our results indicate that the association between the MGP variant and increased risk for hand OA is caused by a lower expression of MGP, which may increase the burden of hand OA by decreased inhibition of cartilage calcification.

Genetic Determinants of Radiographic Knee Osteoarthritis in African Americans

OBJECTIVE

The etiology of knee osteoarthritis (OA), the most common form of arthritis, is complex and may differ by race or ethnicity. In recent years, genetic studies have identified many genetic variants associated with OA, but nearly all the studies were conducted in European whites and Asian Americans. Few studies have focused on the genetics of knee OA in African Americans.

METHODS

We performed a genome-wide association study of radiographic knee OA in 1217 African Americans from 2 North American cohort studies: 590 subjects from the Johnston County Osteoarthritis Project and 627 subjects from the Osteoarthritis Initiative. Analyses were conducted in each cohort separately and combined in an inverse variance fixed effects metaanalysis, which were then included in pathway analyses. We additionally tested 12 single-nucleotide polymorphisms robustly associated with OA in European white populations for association in African Americans.

RESULTS

We identified a genome-wide significant variant in LINC01006 (minor allele frequency 12%; p = 4.11 × 10-9) that is less common in European white populations (minor allele frequency < 3%). Five other independent loci reached suggestive significance (p < 1 × 10-6). In pathway analyses, dorsal/ventral neural tube patterning and iron ion transport pathways were significantly associated with knee OA in African Americans (false discovery rate < 0.05). We found no evidence that previously reported OA susceptibility variants in European whites were associated with knee OA in African Americans.

CONCLUSION

These results highlight differences in the genetic architecture of knee OA between African American and European whites. This finding underscores the need to include more diverse populations in OA genetics studies.

Clinical significance of Matrilin-3 gene polymorphism in Egyptian patients with primary knee osteoarthritis

OBJECTIVE

Osteoarthritis (OA) is a multifactorial, degenerative, and inflammatory disorder of joints causing damage of the articular cartilage, formation of osteophytes, and eburination of the subchondral bone. Matrilin-3 (MATN-3) is a non-collagenous oligomeric extracellular matrix protein (ECM), which is the smallest member of the matrilin family. This study was conducted to identify the potential association and clinical significance of MATN-3 rs8176070 (SNP6) polymorphism in a series of Egyptian patients with primary knee OA.

MATERIAL AND METHODS

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) was used to determine genotypes of MATN-3 SNP6 for 50 primary knee OA patients in addition to 50 healthy subjects of the same sex and age range. Full history was obtained from OA patients, followed by clinical examination, together with clinical assessment of the severity of knee OA using Lequesne Algofunctional Index score and radiological grading using the Kellgren-Lawrence grade scale (KL).

RESULTS

With regard to genotypes of MATN-3 gene SNP6 (rs8176070), a statistically significant difference between OA patients and healthy control subjects was found for the B\b genotype and b allele (p=0.046 and 0.042 respectively), with the prevalence being higher in OA patients with a high risk to develop OA (Odds Ratio [OR]=2.250, 95% CI=1.011-5.008). Patients with the B\b genotype had worse clinical and radiological findings than those with B\B and b\b genotypes.

CONCLUSION

The investigated polymorphism in the MATN-3 gene may reflect the risk and severity of knee OA in Egyptian patients, particularly with the B\b genotype.

Systems approaches in osteoarthritis: Identifying routes to novel diagnostic and therapeutic strategies

Systems orientated research offers the possibility of identifying novel therapeutic targets and relevant diagnostic markers for complex diseases such as osteoarthritis. This review demonstrates that the osteoarthritis research community has been slow to incorporate systems orientated approaches into research studies, although a number of key studies reveal novel insights into the regulatory mechanisms that contribute both to joint tissue homeostasis and its dysfunction. The review introduces both top-down and bottom-up approaches employed in the study of osteoarthritis. A holistic and multiscale approach, where clinical measurements may predict dysregulation and progression of joint degeneration, should be a key objective in future research. The review concludes with suggestions for further research and emerging trends not least of which is the coupled development of diagnostic tests and therapeutics as part of a concerted effort by the osteoarthritis research community to meet clinical needs. © 2017 The Authors. Journal of Orthopaedic Research Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 35:1573-1588, 2017.

Identification of key genes associated with the effect of osmotic stimuli on intervertebral discs using microarray analysis

The present study aimed to explore the effect of osmotic stimuli on intervertebral discs (IVDs) using microarray analysis. Gene expression dataset GSE1648 was downloaded from the Gene Expression Omnibus database. There were 11 IVD cell samples in this dataset, which included 4 hyperosmotic stimuli samples, 3 hypoosmotic stimuli samples and 4 isosmotic stimuli samples. The differentially expressed genes (DEGs) in hyperosmotic or hypoosmotic IVD cells (designated DEGs-hyper or DEGs-hypo) were identified, compared with isosmotic cells, using the limma package of R software. The Database for Annotation, Visualization and Integrated Discovery was used to perform a Gene Ontology (GO) term enrichment analysis for the DEG sets. Protein-protein interaction (PPI) network and microRNA (miRNA) gene-regulatory network data for the DEG sets were obtained using the Human Protein Reference Database (HPRD) and the TargetScan database, respectively, and these networks were constructed and visualized using Cytoscape software. There was a total of 43 DEGs in DEGs-hyper and 9 in DEGs-hypo. Analysis of DEGs-hyper revealed that 41 GO terms were significantly enriched. In total, 376 pairs and 382 nodes were involved in the PPI network, and 1,314 miRNA-gene pairs and 422 nodes were contained in the miRNA-gene-regulated network. The results of the present study indicated that potential target genes (including NCOA3, SOS1, XPO1, ZBTB18, EFNB2 and SOBP) may be involved in the effect of osmotic stimuli on IVD, and the biological processes of apoptosis and cell death may be associated with the effect of high osmolality on IVD disease. The potential targets identified in the present study are more reliable than those identified by previous studies.

Genome-Wide Association Study of Radiographic Knee Osteoarthritis in North American Caucasians

OBJECTIVE

A major barrier to genetic studies of osteoarthritis (OA) is the need to obtain large numbers of individuals with standardized radiographic evaluations for OA. To address this gap, we performed a genome-wide association study (GWAS) of radiographically defined tibiofemoral knee OA in 3,898 cases and 3,168 controls from 4 well-characterized North American cohorts, and we performed replication analysis of previously reported OA loci.

METHODS

We performed meta-analysis using a 2-stage design. Stage 1 (discovery) consisted of a GWAS meta-analysis of radiographic knee OA carried out in the Osteoarthritis Initiative and the Johnston County Osteoarthritis Project. Knee OA was defined as definite osteophytes and possible joint space narrowing or total joint replacement in one or both knees. Stage 2 (validation) was performed in the Multicenter Osteoarthritis Study and the Genetics of Osteoarthritis study. We genotyped lead meta-analysis variants (P ≤ 1 × 10^-4 ) from stage 1 and tested the association between these variants and knee OA. We then combined results from all cohorts in a meta-analysis.

RESULTS

Lead variants from stage 1, representing 49 unique loci, were analyzed in stage 2; none met genome-wide significance in the combined analysis of stage 1 and stage 2. We validated 1 locus (rs4867568 near LSP1P3) with nominal significance (P < 0.05), which was also our top finding in the combined meta-analysis (odds ratio [OR] 0.84 [95% confidence interval (95% CI) 0.79-0.91], P = 3.02 × 10^-6 ). We observed nominally significant associations (P < 0.05) with 3 previously reported OA loci: rs143383 in GDF5 (OR 1.12 [95% CI 1.04-1.21], P = 2.13 × 10^-3 ), rs835487 in CHST11 (OR 0.93 [95% CI 0.85-0.99], P = 0.03), and rs8044769 in FTO (OR 1.10 [95% CI 1.03-1.19], P = 6.13 × 10^-3 ).

CONCLUSION

These findings provide suggestive evidence of a novel knee OA locus and confirm previously reported associations in GDF5, CHST11, and FTO.

The rs4238326 polymorphism in ALDH1A2 gene potentially associated with non-post traumatic knee osteoarthritis susceptibility: a two-stage population-based study

OBJECTIVE

A recent genome-wide association study reported significant associations of genetic variants within the ALDH1A2 gene with osteoarthritis (OA) of the hand in European populations. However, these findings have not been well generalized to other joints, or to other populations.

METHODS

We performed a two-stage population-based case-control study including 196 non-post traumatic knee OA cases and 442 controls in the first stage and independent 143 non-post traumatic knee OA cases and 238 controls in the second stage in a Chinese population by genotyping eight tagging polymorphisms in ALDH1A2.

RESULTS

In the first stage, the single nucleotide polymorphism (SNP) rs4238326 was found to be potentially associated with knee OA risk (additive model: odds ratio [OR] = 0.70; 95% confidence interval [95% CI] = 0.49-1.01; P = 0.055), which was further confirmed in the second stage with similar effect (additive model: OR = 0.60; 95% CI = 0.38-0.95; P = 0.029). After combining the two stages, we found that the variant C allele of rs4238326 was probably associated with decreased risk of knee OA (additive model: OR = 0.65; 95% CI = 0.49-0.86; P = 0.003). Furthermore, interaction analyses showed that rs4238326 interacted multiplicatively with age to contribute to knee OA risk (interaction P = 0.041).

CONCLUSIONS

These findings indicate that the SNP rs4238326 in ALDH1A2 gene may potentially modify individual susceptibility to knee OA in the Chinese population. Beyond that, further studies are warranted to validate and extend our findings, and future functional studies are required to clarify the possible mechanisms.

Radiographic endophenotyping in hip osteoarthritis improves the precision of genetic association analysis

OBJECTIVE

Osteoarthritis (OA) has a strong genetic component but the success of previous genome-wide association studies (GWAS) has been restricted due to insufficient sample sizes and phenotype heterogeneity. Our aim was to examine the effect of clinically relevant endophenotyping according to site of maximal joint space narrowing (maxJSN) and bone remodelling response on GWAS signal detection in hip OA.

METHODS

A stratified GWAS meta-analysis was conducted in 2118 radiographically defined hip OA cases and 6500 population-based controls. Signals were followed up by analysing differential expression of proximal genes for bone remodelling endophenotypes in 33 pairs of macroscopically intact and OA-affected cartilage.

RESULTS

We report suggestive evidence (p<5×10-6) of association at 6 variants with OA endophenotypes that would have been missed by using presence of hip OA as the disease end point. For example, in the analysis of hip OA cases with superior maxJSN versus cases with non-superior maxJSN we detected association with a variant in the LRCH1 gene (rs754106, p=1.49×10-7, OR (95% CIs) 0.70 (0.61 to 0.80)). In the comparison of hypertrophic with non-hypertrophic OA the most significant variant was located between STT3B and GADL1 (rs6766414, p=3.13×10-6, OR (95% CIs) 1.45 (1.24 to 1.69)). Both of these associations were fully attenuated in non-stratified analyses of all hip OA cases versus population controls (p>0.05). STT3B was significantly upregulated in OA-affected versus intact cartilage, particularly in the analysis of hypertrophic and normotrophic compared with atrophic bone remodelling pattern (p=4.2×10-4).

CONCLUSIONS

Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.

DNA methyltransferase 3b regulates articular cartilage homeostasis by altering metabolism

Osteoarthritis (OA) is the most common form of arthritis worldwide. It is a complex disease affecting the whole joint but is generally characterized by progressive degradation of articular cartilage. Recent genome-wide association screens have implicated distinct DNA methylation signatures in OA patients. We show that the de novo DNA methyltransferase (Dnmt) 3b, but not Dnmt3a, is present in healthy murine and human articular chondrocytes and its expression decreases in OA mouse models and in chondrocytes from human OA patients. Targeted deletion of Dnmt3b in murine articular chondrocytes results in an early-onset and progressive postnatal OA-like pathology. RNA-Seq and methylC-Seq analyses of Dnmt3b loss-of-function chondrocytes show that cellular metabolic processes are affected. Specifically, TCA metabolites and mitochondrial respiration are elevated. Importantly, a chondroprotective effect was found following Dnmt3b gain of function in murine articular chondrocytes in vitro and in vivo. This study shows that Dnmt3b plays a significant role in regulating postnatal articular cartilage homeostasis. Cellular pathways regulated by Dnmt3b in chondrocytes may provide novel targets for therapeutic approaches to treat OA.

The genetics revolution in rheumatology: large scale genomic arrays and genetic mapping

Susceptibility to rheumatic diseases, such as osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, systemic lupus erythematosus, juvenile idiopathic arthritis and psoriatic arthritis, includes a large genetic component. Understanding how an individual's genetic background influences disease onset and outcome can lead to a better understanding of disease biology, improved diagnosis and treatment, and, ultimately, to disease prevention or cure. The past decade has seen great progress in the identification of genetic variants that influence the risk of rheumatic diseases. The challenging task of unravelling the function of these variants is ongoing. In this Review, the major insights from genetic studies, gained from advances in technology, bioinformatics and study design, are discussed in the context of rheumatic disease. In addition, pivotal genetic studies in the main rheumatic diseases are highlighted, with insights into how these studies have changed the way we view these conditions in terms of disease overlap, pathways of disease and potential new therapeutic targets. Finally, the limitations of genetic studies, gaps in our knowledge and ways in which current genetic knowledge can be fully translated into clinical benefit are examined.