Selective IL-27 production by intestinal regulatory T cells permits gut-specific regulation of TH17 cell immunity

Regulatory T cells (Treg cells) are instrumental in establishing immunological tolerance. However, the precise effector mechanisms by which Treg cells control a specific type of immune response in a given tissue remains unresolved. By simultaneously studying Treg cells from different tissue origins under systemic autoimmunity, in the present study we show that interleukin (IL)-27 is specifically produced by intestinal Treg cells to regulate helper T17 cell (TH17 cell) immunity. Selectively increased intestinal TH17 cell responses in mice with Treg cell-specific IL-27 ablation led to exacerbated intestinal inflammation and colitis-associated cancer, but also helped protect against enteric bacterial infection. Furthermore, single-cell transcriptomic analysis has identified a CD83+CD62Llo Treg cell subset that is distinct from previously characterized intestinal Treg cell populations as the main IL-27 producers. Collectively, our study uncovers a new Treg cell suppression mechanism crucial for controlling a specific type of immune response in a particular tissue and provides further mechanistic insights into tissue-specific Treg cell-mediated immune regulation.

Selective IL-27 production by intestinal regulatory T cells permits gut-specific regulation of Th17 immunity

Regulatory T (Treg) cells are instrumental in establishing immunological tolerance. However, the precise effector mechanisms by which Treg cells control a specific type of immune response in a given tissue remains unresolved. By simultaneously studying Treg cells from different tissue origins under systemic autoimmunity, here we show that IL-27 is specifically produced by intestinal Treg cells to regulate Th17 immunity. Selectively increased intestinal Th17 responses in mice with Treg cell-specific IL-27 ablation led to exacerbated intestinal inflammation and colitis-associated cancer, but also helped protect against enteric bacterial infection. Furthermore, single-cell transcriptomic analysis has identified a CD83+TCF1+ Treg cell subset that is distinct from previously characterized intestinal Treg cell populations as the main IL-27 producers. Collectively, our study uncovers a novel Treg cell suppression mechanism crucial for controlling a specific type of immune response in a particular tissue, and provides further mechanistic insights into tissue-specific Treg cell-mediated immune regulation.

Large-scale plasma proteomics can reveal distinct endotypes in chronic obstructive pulmonary disease and severe asthma

BACKGROUND

Chronic airway diseases including chronic obstructive pulmonary disease (COPD) and asthma are heterogenous in nature and endotypes within are underpinned by complex biology. This study aimed to investigate the utility of proteomic profiling of plasma combined with bioinformatic mining, and to define molecular endotypes and expand our knowledge of the underlying biology in chronic respiratory diseases.

METHODS

The plasma proteome was evaluated using an aptamer-based affinity proteomics platform (SOMAscan®), representing 1238 proteins in 34 subjects with stable COPD and 51 subjects with stable but severe asthma. For each disease, we evaluated a range of clinical/demographic characteristics including bronchodilator reversibility, blood eosinophilia levels, and smoking history. We applied modified bioinformatic approaches used in the evaluation of RNA transcriptomics.

RESULTS

Subjects with COPD and severe asthma were distinguished from each other by 365 different protein abundancies, with differential pathway networks and upstream modulators. Furthermore, molecular endotypes within each disease could be defined. The protein groups that defined these endotypes had both known and novel biology including groups significantly enriched in exosomal markers derived from immune/inflammatory cells. Finally, we observed associations to clinical characteristics that previously have been under-explored.

CONCLUSION

This investigational study evaluating the plasma proteome in clinically-phenotyped subjects with chronic airway diseases provides support that such a method can be used to define molecular endotypes and pathobiological mechanisms that underpins these endotypes. It provided new concepts about the complexity of molecular pathways that define these diseases. In the longer term, such information will help to refine treatment options for defined groups.

Searchlight: automated bulk RNA-seq exploration and visualisation using dynamically generated R scripts

BACKGROUND

Once bulk RNA-seq data has been processed, i.e. aligned and then expression and differential tables generated, there remains the essential process where the biology is explored, visualized and interpreted. Without the use of a visualisation and interpretation pipeline this step can be time consuming and laborious, and is often completed using R. Though commercial visualisation and interpretation pipelines are comprehensive, freely available pipelines are currently more limited.

RESULTS

Here we demonstrate Searchlight, a freely available bulk RNA-seq visualisation and interpretation pipeline. Searchlight provides: a comprehensive statistical and visual analysis, focusing on the global, pathway and single gene levels; compatibility with most differential experimental designs irrespective of organism or experimental complexity, via three workflows; reports; and support for downstream user modification of plots via user-friendly R-scripts and a Shiny app. We show that Searchlight offers greater automation than current best tools (VIPER and BioJupies). We demonstrate in a timed re-analysis study, that alongside a standard bulk RNA-seq processing pipeline, Searchlight can be used to complete bulk RNA-seq projects up to the point of manuscript quality figures, in under 3 h.

CONCLUSIONS

Compared to a manual R based analysis or current best freely available pipelines (VIPER and BioJupies), Searchlight can reduce the time and effort needed to complete bulk RNA-seq projects to manuscript level. Searchlight is suitable for bioinformaticians, service providers and bench scientists. https://github.com/Searchlight2/Searchlight2 .

Sensing of apoptotic cells through Axl causes lung basal cell proliferation in inflammatory diseases

Epithelial cell proliferation, division, and differentiation are critical for barrier repair following inflammation, but the initial trigger for this process is unknown. Here we define that sensing of apoptotic cells by the TAM receptor tyrosine kinase Axl is a critical indicator for tracheal basal cell expansion, cell cycle reentry, and symmetrical cell division. Furthermore, once the pool of tracheal basal cells has expanded, silencing of Axl is required for their differentiation. Genetic depletion of Axl triggers asymmetrical cell division, leading to epithelial differentiation and ciliated cell regeneration. This discovery has implications for conditions associated with epithelial barrier dysfunction, basal cell hyperplasia, and continued turnover of dying cells in patients with chronic inflammatory pulmonary diseases.

MEK inhibition drives anti-viral defence in RV but not RSV challenged human airway epithelial cells through AKT/p70S6K/4E-BP1 signalling

BACKGROUND

The airway epithelium is a major target tissue in respiratory infections, and its antiviral response is mainly orchestrated by the interferon regulatory factor-3 (IRF3), which subsequently induces type I (β) and III (λ) interferon (IFN) signalling. Dual specificity mitogen-activated protein kinase kinase (MEK) pathway contributes to epithelial defence, but its role in the regulation of IFN response in human primary airway epithelial cells (AECs) is not fully understood. Here, we studied the impact of a small-molecule inhibitor (MEKi) on the IFN response following challenge with two major respiratory viruses rhinovirus (RV2) and respiratory syncytial virus (RSVA2) and a TLR3 agonist, poly(I:C).

METHODS

The impact of MEKi on viral load and IFN response was evaluated in primary AECs with or without a neutralising antibody against IFN-β. Quantification of viral load was determined by live virus assay and absolute quantification using qRT-PCR. Secretion of cytokines was determined by AlphaLISA/ELISA and expression of interferon-stimulated genes (ISGs) was examined by qRT-PCR and immunoblotting. A poly(I:C) model was also used to further understand the molecular mechanism by which MEK controls IFN response. AlphaLISA, siRNA-interference, immunoblotting, and confocal microscopy was used to investigate the effect of MEKi on IRF3 activation and signalling. The impact of MEKi on ERK and AKT signalling was evaluated by immunoblotting and AlphaLISA.

RESULTS

Here, we report that pharmacological inhibition of MEK pathway augments IRF3-driven type I and III IFN response in primary human AECs. MEKi induced activation of PI3K-AKT pathway, which was associated with phosphorylation/inactivation of the translational repressor 4E-BP1 and activation of the protein synthesis regulator p70 S6 kinase, two critical translational effectors. Elevated IFN-β response due to MEKi was also attributed to decreased STAT3 activation, which consequently dampened expression of the transcriptional repressor of IFNB1 gene, PRDI-BF1. Augmented IFN response translated into inhibition of rhinovirus 2 replication in primary AECs but not respiratory syncytial virus A2.

CONCLUSIONS

Our findings unveil MEK as a key molecular mechanism by which rhinovirus dampens the epithelial cell's antiviral response. Our study provides a better understanding of the role of signalling pathways in shaping the antiviral response and suggests the use of MEK inhibitors in anti-viral therapy against RV.

Upper Airway Cell Transcriptomics Identify a Major New Immunological Phenotype with Strong Clinical Correlates in Young Children with Acute Wheezing

Asthma exacerbations are triggered by rhinovirus infections. We employed a systems biology approach to delineate upper-airway gene network patterns underlying asthma exacerbation phenotypes in children. Cluster analysis unveiled distinct IRF7^hi versus IRF7^lo molecular phenotypes, the former exhibiting robust upregulation of Th1/type I IFN responses and the latter an alternative signature marked by upregulation of cytokine and growth factor signaling and downregulation of IFN-γ. The two phenotypes also produced distinct clinical phenotypes. For IRF7^lo children, symptom duration prior to hospital presentation was more than twice as long from initial symptoms (p = 0.011) and nearly three times as long for cough (p < 0.001), the odds ratio of admission to hospital was increased more than 4-fold (p = 0.018), and time to recurrence was shorter (p = 0.015). In summary, our findings demonstrate that asthma exacerbations in children can be divided into IRF7^hi versus IRF7^lo phenotypes with associated differences in clinical phenotypes.

Pain in knee osteoarthritis is associated with variation in the neurokinin 1/substance P receptor (TACR1) gene

BACKGROUND

Substance P (SP) is a pain- and inflammation-related neuropeptide which preferentially binds to the neurokinin receptor 1 (NK₁ ). SP and NK₁ receptors have been implicated in joint pain, inflammation and damage in animal models and human studies of osteoarthritis (OA). The aim of this study was to test if genetic variation at the neurokinin 1 receptor gene (TACR1) is associated with pain in individuals with radiographic knee OA.

METHODS

Participants from the Genetics of OA and Lifestyle study were used for the discovery group (n = 1615). Genotype data for six SNPs selected to cover most variation in the TACR1 gene were used to test for an association with symptomatic OA. Replication analysis was performed using data from the Chingford 1000 Women Study, Hertfordshire Cohort Study, Tasmanian Older Adult Cohort Study and the Clearwater OA Study. In total, n = 1715 symptomatic OA and n = 735 asymptomatic OA individuals were analysed.

RESULTS

Out of six SNPs tested in the TACR1 gene, one (rs11688000) showed a nominally significant association with a decreased risk of symptomatic OA in the discovery cohort. This was then replicated in four additional cohorts. After adjusting for age, gender, body mass index and radiographic severity, the G (minor) allele at rs11688000 was associated with a decreased risk of symptomatic OA compared to asymptomatic OA cases (p = 9.90 × 10^-4 , OR = 0.79 95% 0.68-0.90 after meta-analysis).

CONCLUSIONS

This study supports a contribution from the TACR1 gene in human OA pain, supporting further investigation of this gene's function in OA.

SIGNIFICANCE

This study contributes to the knowledge of the genetics of painful osteoarthritis, a condition which affects millions of individuals worldwide. Specifically, a contribution from the TACR1 gene to modulating pain sensitivity in osteoarthritis is suggested.

Bone Area Provides a Responsive Outcome Measure for Bone Changes in Short-term Knee Osteoarthritis Studies

Objective.

To analyze the 3-D bone area from an osteoarthritis (OA) cohort demonstrating no change in cartilage thickness.

Methods.

Twenty-seven women with painful medial knee OA had magnetic resonance images at 0, 3, and 6 months. Images were analyzed using active appearance models.

Results.

At 3 and 6 months, the mean change in medial femoral bone area was 0.34% (95% CI 0.04–0.64) and 0.61% (95% CI 0.32–0.90), respectively. Forty-one percent of the subjects had progression greater than the smallest detectable difference at 6 months.

Conclusion.

In this small cohort at high risk of OA progression, bone area changed at 3 and 6 months when cartilage morphometric measures did not.

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10-8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.

Beer and wine consumption and risk of knee or hip osteoarthritis: a case control study

Introduction

The aim of this study was to investigate the association between alcoholic and non-alcoholic beverages and knee or hip osteoarthritis (OA).

Methods

We conducted a case–control study of Caucasian men and women aged 45 to 86 years of age from Nottingham, UK. Cases had clinically severe symptoms and radiographic knee or hip OA; controls had no symptoms and no radiographic knee or hip OA. Exposure information was sought using interview-based questionnaires and a semi-quantitative food frequency questionnaire to assess beverage consumption at ages 21 to 50 years. Odds ratios (ORs), adjusted ORs (aORs), 95% confidence intervals (CI) and P values were estimated using logistic regression models.

Results

A total of 1,001 knee OA, 993 hip OA and 933 control participants were included in the study. Increasing beer consumption was associated with an increasing risk of OA (P for trend ≤0.001). Compared to those who did not consume beer, aORs for people who consumed 20 or more servings of beer were 1.93 (95% CI 1.26 to 2.94) and 2.15 (95% CI 1.45 to 3.19) for knee OA and hip OA, respectively. In contrast, increasing levels of wine consumption were associated with decreased likelihood of knee OA (P for trend <0.001). Compared to those who did not consume wine, aOR for knee OA among those who consumed 4 to 6 glasses of wine per week and ≥7 glasses of wine per week was 0.55 (95% CI 0.34 to 0.87) and 0.48 (95% CI 0.29 to 0.80), respectively. No association was identified between non-alcoholic beverages and knee or hip OA.

Conclusions

Beer consumption appears to be a risk factor for knee and hip OA whereas consumption of wine has a negative association with knee OA. The mechanism behind these findings is speculative but warrants further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0534-4) contains supplementary material, which is available to authorized users.

The Axl receptor tyrosine kinase is a discriminator of macrophage function in the inflamed lung

Much of the biology surrounding macrophage functional specificity has arisen through examining inflammation-induced polarizing signals, but this also occurs in homeostasis, requiring tissue-specific environmental triggers that influence macrophage phenotype and function. The TAM receptor family of receptor tyrosine kinases (Tyro3, Axl and MerTK) mediates the non-inflammatory removal of apoptotic cells by phagocytes through the bridging phosphatidylserine-binding molecules growth arrest-specific 6 (Gas6) or Protein S. We show that one such TAM receptor (Axl) is exclusively expressed on mouse airway macrophages, but not interstitial macrophages and other lung leukocytes, under homeostatic conditions and is constitutively ligated to Gas6. Axl expression is potently induced by granulocyte-macrophage colony-stimulating factor expressed in the healthy and inflamed airway, and by type I interferon or Toll-like receptor-3 stimulation on human and mouse macrophages, indicating potential involvement of Axl in apoptotic cell removal under inflammatory conditions. Indeed, an absence of Axl does not cause sterile inflammation in health, but leads to exaggerated lung inflammatory disease upon influenza infection. These data imply that Axl allows specific identification of airway macrophages, and that its expression is critical for macrophage functional compartmentalization in the airspaces or lung interstitium. We propose that this may be a critical feature to prevent excessive inflammation because of secondary necrosis of apoptotic cells that have not been cleared by efferocytosis.

The effect of FTO variation on increased osteoarthritis risk is mediated through body mass index: a Mendelian randomisation study

OBJECTIVE

Variation in the fat mass and obesity-associated (FTO) gene influences susceptibility to obesity. A variant in the FTO gene has been implicated in genetic risk to osteoarthritis (OA). We examined the role of the FTO polymorphism rs8044769 in risk of knee and hip OA in cases and controls incorporating body mass index (BMI) information.

METHODS

5409 knee OA patients, 4355 hip OA patients and up to 5362 healthy controls from 7 independent cohorts from the UK and Australia were genotyped for rs8044769. The association of the FTO variant with OA was investigated in case/control analyses with and without BMI adjustment and in analyses matched for BMI category. A mendelian randomisation approach was employed using the FTO variant as the instrumental variable to evaluate the role of overweight on OA.

RESULTS

In the meta-analysis of all overweight (BMI≥25) samples versus normal-weight controls irrespective of OA status the association of rs8044769 with overweight is highly significant (OR[CIs] for allele G=1.14 [01.08 to 1.19], p=7.5×10(-7)). A significant association with knee OA is present in the analysis without BMI adjustment (OR[CIs]=1.08[1.02 to 1.14], p=0.009) but the signal fully attenuates after BMI adjustment (OR[CIs]=0.99[0.93 to 1.05], p=0.666). We observe no evidence for association in the BMI-matched meta-analyses. Using mendelian randomisation approaches we confirm the causal role of overweight on OA.

CONCLUSIONS

Our data highlight the contribution of genetic risk to overweight in defining risk to OA but the association is exclusively mediated by the effect on BMI. This is consistent with what is known of the biology of the FTO gene and supports the causative role of high BMI in OA.

Self-reported adult footwear and the risks of lower limb osteoarthritis: the GOAL case control study

Background

Biomechanical factors may play a role in osteoarthritis (OA) development and progression. Previous biomechanical studies have indicated that types of footwear may modulate forces across the knee joint, and high heeled womens’ shoes in particular are hypothesised to be detrimental to lower limb joint health. This analysis of data from a case control study investigated persistent users of different adult footwear for risks of knee and hip OA. Our underlying hypotheses were that high heeled, narrow heeled, and hard soled shoe types were putative risk factors for lower limb OA.

Methods

Data on footwear were initially obtained from participants during the Genetics of Osteoarthritis and Lifestyle (GOAL) hospital-based, case control study using standardised interview-delivered questionnaires. An additional questionnaire was later sent to GOAL study participants to verify findings and to further investigate specific shoe use per decade of life. Persistent users of footwear types (high or narrow heel; sole thickness or hardness) were identified from early adulthood. Participants were grouped into single sex knee OA, hip OA or control groups. Adjusted odds ratios (aOR) and 95% confidence interval (CI) were calculated.

Results

Univariate analysis of persistent users of women’s high heeled and narrow heeled shoes during early adulthood showed negative associations with knee OA and hip OA. After logistic regression, persistent narrow heel users were associated with less risk of OA (knee OA aOR 0.59, 95% CI 0.35 – 1.00 and hip aOR: 0.50, 95% CI 0.30 – 0.85), and other analyses were not statistically significant. Further analysis suggested that women with hip OA may have stopped wearing high and narrow heeled footwear to attenuate hip pain in early adulthood. Consistent associations between shoe soles and OA were not found.

Conclusions

In general, persistent users of high and narrow heeled shoes during early adulthood had a negative association with knee or hip OA. This does not necessarily imply a causal relationship, as changing footwear during early adulthood to modulate index joint pain may provide a possible explanation. Despite the findings of previous biomechanical studies of high heels, we did not find a positive association between women’s shoes and lower limb osteoarthritis.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2474-15-308) contains supplementary material, which is available to authorized users.

Evidence of changes to skeletal muscle contractile properties during the initiation of disease in the ageing guinea pig model of osteoarthritis

Background

Osteoarthritis (OA) is the most common joint disorder in the world and represents the leading cause of pain and disability in the elderly population. Advancing age remains the single greatest risk factor for OA. Several studies have characterised disease development in the guinea pig ageing model of OA in terms of its joint histopathology and inflammatory cytokine profile. However, the quadriceps muscle has yet to be studied in relation to age-related disease onset or early disease progression. Therefore, we examined whether the initiation of OA in the Dunkin Hartley guinea pig is associated with changes in the quadriceps skeletal muscle. Male Dunkin Hartley guinea pigs (N = 24) were group housed with free access to standard guinea pig chow and water. At 2, 3, 5 and 7 months of age, six animals were selected based on their proximity to the median weight of the cohort. OA severity was graded at each time point by the assessment of toluidine blue stained step coronal sections of the total knee joint. Serum CTX II was measured as a potential biomarker of OA severity. Myosin Heavy Chain (MHC) isoforms were determined by a validated real-time PCR assay. Oxidative and glycolytic potential was determined in quadriceps homogenates via the measurement of ICDH and LDH activity.

Results

Initiation of OA in the DH strain guinea pig occurred between 2 and 3 months of age and progressed until 7 months when the final analyses were conducted. Serum CTX II significantly decreased during this early period of OA initiation and levels were unrelated to the histopathological severity of knee OA at any of the time points assessed. MHC mRNA measurements revealed a significant elevation in MHC IIX mRNA (associated with fast-twitch skeletal muscle fibres) coincident with the initiation of OA at 3 months of age, with preliminary findings suggestive of a positive correlation to OA severity at this time point.

Conclusions

These preliminary findings suggest that disease initiation in the ageing guinea pig model of OA is not associated with overt quadriceps muscle atrophy but instead is coincident with altered expression of mRNAs associated with quadriceps skeletal muscle contractile properties (specifically fast-twitch MHC IIX).

Gene-environment interaction between body mass index and transforming growth factor beta 1 (TGFβ1) gene in knee and hip osteoarthritis

Introduction

The objective was to investigate potential gene-environment interaction between body mass index (BMI) and each of eight TGFβ1 polymorphisms in knee and hip osteoarthritis (OA).

Methods

We conducted a case-control study of Caucasian men and women aged 45 to 86 years from Nottingham, United Kingdom (Genetics of OA and Lifestyle (GOAL) study). Cases had clinically severe symptoms and radiographic knee or hip OA; controls had no symptoms and no radiographic knee/hip OA. We used logistic regression to investigate the association of TGFβ1 polymorphisms and OA when stratifying by BMI. Knee and hip OA were analyzed separately with adjustment for potential confounders. Additive and multiplicative interactions were examined.

Results

2,048 cases (1,042 knee OA, 1,006 hip OA) and 967 controls were studied. For hip OA, the highest risk was in overweight (BMI ≥25 kg/m²) individuals with the variant allele of single-nucleotide polymorphism (SNP) rs1800468 (odds ratio (OR) 2.21, 95% confidence interval (CI) 1.55, 3.15). Evaluation of gene-environment interaction indicated significant synergetic interaction (relative excess risk due to interaction (RERI) = 0.93, synergy index (SI) = 4.33) with an attributable proportion due to interaction (AP) of 42% (AP = 0.42; 95% CI 0.16, 0.68). Multiplicative interaction was also significant (OR for interaction (ORINT) = 2.27, P = 0.015). For knee OA, the highest risk was in overweight individuals with homozygous genotype 11 of SNP rs2278422 (OR = 6.95, P <0.001). In contrast, the variant allele indicated slightly lower risks (OR = 4.72, P <0.001), a significant antagonistic interaction (RERI = -2.66, SI = 0.59), AP = -0.56 (95%CI -0.94, -0.17) and a significant multiplicative interaction (ORINT = 0.47, P = 0.013).

Conclusion

TGFβ1 gene polymorphisms interact with being overweight to influence the risk of large joint OA.

The genetic contribution to severe post-traumatic osteoarthritis

Objective

to compare the combined role of genetic variants loci associated with risk of knee or hip osteoarthritis (OA) in post-traumatic (PT) and non-traumatic (NT) cases of clinically severe OA leading to total joint replacement.

Methods

A total of 1590 controls, 2168 total knee replacement (TKR) cases (33.2% PT) and 1567 total hip replacement (THR) cases (8.7% PT) from 2 UK cohorts were genotyped for 12 variants previously reported to be reproducibly associated with risk of knee or hip OA. A genetic risk score was generated and the association with PT and NT TKR and THR was assessed adjusting for covariates.

Results

For THR, each additional genetic risk variant conferred lower risk among PT cases (OR=1.07, 95% CI 0.96 to 1.19; p=0.24) than NT cases (OR 1.11, 95% CI 1.06 to 1.17; p=1.55×10⁻⁵). In contrast, for TKR, each risk variant conferred slightly higher risk among PT cases (OR 1.12, 95% CI 1.07 to 1.19; p=1.82×10⁻⁵) than among NT cases (OR 1.08, 95% CI 1.03 to 1.1; p=0.00063).

Conclusions

Based on the variants reported to date PT TKR cases have at least as high a genetic contribution as NT cases.

Temporal relationship between serum adipokines, biomarkers of bone and cartilage turnover, and cartilage volume loss in a population with clinical knee osteoarthritis

OBJECTIVE

The association of obesity with both hand and knee osteoarthritis (OA) is suggestive of a link between dysfunctional metabolism and joint integrity. Given the role of adipokines in mediating bone and cartilage homeostasis, we undertook this study to examine the relationship between adipokines and bone and cartilage biomarkers in a population of subjects with OA, and to determine whether adipokine levels predicted 2-year cartilage integrity.

METHODS

One hundred seventeen subjects underwent magnetic resonance imaging at baseline and at 2-year followup. Cartilage volume was assessed from these images. Serum adipokine levels were measured at baseline. Bone and cartilage biomarker levels were measured at baseline and at 2-year followup. Linear regression was used to examine the relationship between baseline levels of adipokines and adipokine receptors (leptin, soluble leptin receptor [sOB-Rb], resistin, and adiponectin) and changes in levels of bone biomarkers (osteocalcin, N-terminal type I procollagen propeptide [PINP], C-terminal crosslinking telopeptide of type I collagen, N-terminal crosslinking telopeptide of type I collagen, or C-terminal crosslinking telopeptide of type I collagen generated by matrix metalloproteinases), levels of cartilage biomarkers (cartilage oligomeric matrix protein, N-terminal type IIA procollagen propeptide [PIIANP], or C2C), cartilage defects score, and cartilage volume over 2 years.

RESULTS

Baseline leptin was associated with increased levels of bone formation biomarkers (osteocalcin and PINP) over 2 years, while sOB-Rb was associated with reduced levels of osteocalcin. Baseline sOB-Rb was associated with reduced levels of the cartilage formation biomarker PIIANP, an increased cartilage defects score, and increased cartilage volume loss over 2 years. All results were independent of age, sex, and body mass index.

CONCLUSION

The findings of this study support the concept that serum adipokines may provide a nonmechanical link between obesity and joint integrity (which may be mediated by bone and cartilage turnover) that subsequently results in changes to the cartilage defects score and cartilage volume loss. This may facilitate our understanding of the mechanisms by which obesity is involved in the pathogenesis of OA.

Degradation of myofibrils from rabbit, chicken and beef by cathepsin l and lysosomal lysates

The degradation of rabbit, chicken and beef myofibrils by cathepsin L or lysosomal lysates was studied by SDS-polyacrylamide-gel electrophoresis and electron microscopy (EM). Similar degradation patterns were observed for each myofibrillar preparation incubated with cathepsin L, except that myosin heavy chain and tropomyosin of beef were more susceptible than those of rabbit and chicken. Otherwise, troponin T, troponin in I and C-protein were rapidly degraded with slower degradation of titin, nebulin, myosin heavy chain, α-actinin, α-tropomyosin, actin and myosin light chains, LC1 and LC2. However, the component of 30 000 Mr was found to be further degraded to smaller peptides. Degradation at pH 5·5 (approximate post-mortem limit value) was faster than at pH 6·0 but slower than at pH 5·0. A number of new protein bands were identified (130 000, 120 000, 90 000, 85 000, 80 000, 31 000 and 30 000 Mr). The degradation patterns of rabbit myofibrils by rabbit muscle lysosomal lysates were similar to that of myofibrils incubated with purified cathepsin L except for the retention of the 30 000 Mr component and reduced degradation of actin, due presumably to the reduced amount or stability of cathepsin L in the crude enzyme preparations. Electron micrographs revealed that myofibrillar degradation by cathepsin L occurred preferentially at the Z-lines leading to removal of the Z-line proteins and fracturing of the myofibrils at these sites. Catheptic damage was seen to be most rapid in chicken myofibrils and least rapid in beef myofibrils consistent with the more rapid conditioning process in chicken.

Genetic contribution to radiographic severity in osteoarthritis of the knee

OBJECTIVE

Knee osteoarthritis (OA) has a significant genetic component. The authors have assessed the role of three variants reported to influence risk of knee OA with p<5×10-8 in determining patellofemoral and tibiofemoral Kellgren Lawrence (K/L) grade in knee OA cases.

METHODS

3474 knee OA cases with sky-line and weight-bearing antero-posterior x-rays of the knee were selected based on the presentation of K/L grade ≥2 at either the tibiofemoral or patellofemoral compartments for one or both knees. Patients belonging to three UK cohorts, were genotyped for rs143383, rs4730250 and rs11842874 mapping to the GDF5, COG5 and MCF2L genes, respectively. The association between tibiofemoral K/L grade and patellofemoral K/L grade was assessed after adjusting for age, gender and body mass index.

RESULTS

No significant association was found between the rs4730250 and radiographic severity. The rs11842874 mapping to MCF2L was found to be nominally significantly associated with patellofemoral K/L grade as a quantitative trait (p=0.027) but not as a binary trait. The GDF5 single nucleotide polymorphism rs143383 was associated with tibiofemoral K/L grade (β=0.05 (95% CI 0.02 to 0.08) p=0.0011).

CONCLUSIONS

Our data indicate that within individuals affected by radiographic knee OA, OAGDF5 has a modest but significant effect on radiographic severity after adjustment for the major risk factors.

A role for PACE4 in osteoarthritis pain: evidence from human genetic association and null mutant phenotype

OBJECTIVES

The aim of this study was to assess if genetic variation in the PACE4 (paired amino acid converting enzyme 4) gene Pcsk6 influences the risk for symptomatic knee osteoarthritis (OA).

METHODS

Ten PCSK6 single nucleotide polymorphisms were tested for association in a discovery cohort of radiographic knee OA (n=156 asymptomatic and 600 symptomatic cases). Meta-analysis of the minor allele at rs900414 was performed in three additional independent cohorts (total n=674 asymptomatic and 2068 symptomatic). Pcsk6 knockout mice and wild-type C57BL/6 mice were compared in a battery of algesiometric assays, including hypersensitivity in response to intraplantar substance P, pain behaviours in response to intrathecal substance P and pain behaviour in the abdominal constriction test.

RESULTS

In the discovery cohort of radiographic knee OA, an intronic single nucleotide polymorphism at rs900414 was significantly associated with symptomatic OA. Replication in three additional cohorts confirmed that the minor allele at rs900414 was consistently increased among asymptomatic compared to symptomatic radiographic knee OA cases in all four cohorts. A fixed-effects meta-analysis yielded an OR=1.35 (95% CI 1.17 to 1.56; p=4.3×10(-5) and no significant between-study heterogeneity). Studies in mice revealed that Pcsk6 knockout mice were significantly protected against pain in a battery of algesiometric assays.

CONCLUSIONS

These results suggest that a variant in PCSK6 is strongly associated with protection against pain in knee OA, offering some insight as to why, in the presence of the same structural damage, some individuals develop chronic pain and others are protected. Studies in Pcsk6 null mutant mice further implicate PACE4 in pain.

Large scale replication study of the association between HLA class II/BTNL2 variants and osteoarthritis of the knee in European-descent populations

Osteoarthritis (OA) is the most common form of arthritis and a major cause of disability. This study evaluates the association in Caucasian populations of two single nucleotide polymorphisms (SNPs) mapping to the Human Leukocyte Antigen (HLA) region and deriving from a genome wide association scan (GWAS) of knee OA in Japanese populations. The frequencies for rs10947262 were compared in 36,408 controls and 5,749 knee OA cases from European-descent populations. rs7775228 was tested in 32,823 controls and 1,837 knee OA cases of European descent. The risk (major) allele at rs10947262 in Caucasian samples was not significantly associated with an odds ratio (OR)  = 1.07 (95%CI 0.94 -1.21; p = 0.28). For rs7775228 the meta-analysis resulted in OR = 0.94 (95%CI 0.81-1.09; p = 0.42) for the allele associated with risk in the Japanese GWAS. In Japanese individuals these two SNPs are in strong linkage disequilibrium (LD) (r(2) = 0.86) with the HLA class II haplotype DRB1*1502 DQA1*0103 DQB1*0601 (frequency 8%). In Caucasian and Chinese samples, using imputed data, these SNPs appear not to be in LD with that haplotype (r(2)<0.07). The rs10947262 and rs7775228 variants are not associated with risk of knee OA in European descent populations and they do not appear tag the same HLA class II haplotype as they do in Japanese individuals.

The Ile585Val TRPV1 variant is involved in risk of painful knee osteoarthritis

Objective

To assess if a coding variant in the gene encoding transient receptor potential cation channel, subfamily V, member 1 (TRPV1) is associated with genetic risk of painful knee osteoarthritis (OA).

Methods

The Ile585Val TRPV1 variant encoded by rs8065080 was genotyped in 3270 cases of symptomatic knee OA, 1098 cases of asymptomatic knee OA and 3852 controls from seven cohorts from the UK, the USA and Australia. The genetic association between the low-pain genotype Ile–Ile and risk of symptomatic and asymptomatic knee OA was assessed.

Results

The TRPV1 585 Ile–Ile genotype, reported to be associated with lower thermal pain sensitivity, was associated with a lower risk of symptomatic knee OA in a comparison of symptomatic cases with healthy controls, with an odds ratio (OR) of 0.75 (95% CI 0.64 to 0.88; p=0.00039 by meta-analysis) after adjustment for age, sex and body mass index. No difference was seen between asymptomatic OA cases and controls (OR=1.02, 95% CI 0.82 to 1.27 p=0.86) but the Ile–Ile genotype was associated with lower risk of symptomatic versus asymptomatic knee OA adjusting for covariates and radiographic severity (OR=0.73, 95% CI 0.57 to 0.94 p=0.0136). TRPV1 expression in articular cartilage was increased by inflammatory cytokines (tumour necrosis factor α and interleukin 1). However, there were no differences in TRPV1 expression in healthy and arthritic synovial tissue.

Conclusions

A genotype involved in lower peripheral pain sensitivity is significantly associated with a decreased risk of painful knee OA. This indicates a role for the pro-nociceptive gene TRPV1 in genetic susceptibility to symptomatic knee OA, which may also be influenced by a role for this molecule in cartilage function.

Large-scale meta-analysis of interleukin-1 beta and interleukin-1 receptor antagonist polymorphisms on risk of radiographic hip and knee osteoarthritis and severity of knee osteoarthritis

OBJECTIVE

To clarify the role of common genetic variation in the Interleukin-1β (IL1B) and Interleukin-1R antagonist (IL1RN) genes on risk of knee and hip osteoarthritis (OA) and severity of knee OA by means of large-scale meta-analyses.

METHODS

We searched PubMed for articles assessing the role of IL1B and IL1RN polymorphisms/haplotypes on the risk of hip and/or knee OA. Novel data were included from eight unpublished studies. Meta-analyses were performed using fixed- and random-effects models with a total of 3595 hip OA and 5013 knee OA cases, and 6559 and 9132 controls respectively. The role of ILRN haplotypes on radiographic severity of knee OA was tested in 1918 cases with Kellgren-Lawrence (K/L) 1 or 2 compared to 199 cases with K/L 3 or 4.

RESULTS

The meta-analysis of six published studies retrieved from the literature search and eight unpublished studies showed no evidence of association between common genetic variation in the IL1B or IL1RN genes and risk of hip OA or knee OA (P>0.05 for rs16944, rs1143634, rs419598 and haplotype C-G-C (rs1143634, rs16944 and rs419598) previously implicated in risk of hip OA). The C-T-A haplotype formed by rs419598, rs315952 and rs9005, previously implicated in radiographic severity of knee OA, was associated with reduced severity of knee OA (odds ratio (OR)=0.71 95%CI 0.56-0.91; P=0.006, I(2)=74%), and achieved borderline statistical significance in a random-effects model (OR=0.61 95%CI 0.35-1.06 P=0.08).

CONCLUSION

Common genetic variation in the Interleukin-1 region is not associated with prevalence of hip or knee OA but our data suggest that IL1RN might have a role in severity of knee OA.

Lifetime body mass index, other anthropometric measures of obesity and risk of knee or hip osteoarthritis in the GOAL case-control study

OBJECTIVE

To examine the risk of large joint osteoarthritis (OA) in those becoming overweight during early adult life, and to assess the risks associated with high body mass index (BMI) and other anthropometric measures of obesity.

METHODS

BMI, waist and hip circumference were measured in the GOAL case-control study comprising hip OA cases (n=1007), knee OA cases (n=1042) and asymptomatic controls (n=1121). Retrospective estimates of lifetime weight, body shape and other risk factors were collected using an interview-lead questionnaire. Odds ratios (ORs), adjusted OR (aOR), 95% confidence intervals (CIs) and P values were calculated using logistic regression analysis.

RESULTS

BMI was associated with knee OA (aOR 2.68, 95% CI 2.33-3.09, P-trend<0.001) and hip OA (aOR 1.65, 95% CI 1.46-1.87, P-trend<0.001). Those who became overweight earlier in adulthood showed higher risks of lower limb OA (P-trend<0.001 for knee OA and hip OA). Self-reported body shape was also associated with knee OA and hip OA, following a similar pattern to current and life-course BMI measures. Waist:hip ratio (WHR) at time of examination did not associate with OA independently of BMI, except in women-only analysis. Waist circumference was associated with lower limb OA risk.

CONCLUSIONS

Becoming overweight earlier in adult life increased the risks of knee OA and hip OA. Different distribution patterns of adiposity may be related to OA risk in women.

Involvement of different risk factors in clinically severe large joint osteoarthritis according to the presence of hand interphalangeal nodes

OBJECTIVE

To quantify the differences in risk factors influencing total hip replacement (THR) and total knee replacement (TKR) based on the presence versus absence of multiple interphalangeal nodes in 2 or more rays of the fingers of each hand in patients with large joint osteoarthritis (OA).

METHODS

A group of 3,800 patients with large joint OA who underwent total joint replacement (1,201 of whom had the nodal phenotype) and 1,906 control subjects from 2 case-control studies and a population-based cohort in the UK were studied. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated for the risk of total joint replacement in association with age, sex, body mass index (BMI), height, and prevalence of the T allele in the GDF5 rs143383 polymorphism. ORs for total joint replacement were compared between cases of nodal OA and cases of non-nodal OA and between patients who underwent TKR and those who underwent THR.

RESULTS

Age, sex, and BMI had significantly higher ORs for an association with total joint replacement in nodal OA cases than in non-nodal OA cases. The GDF5 polymorphism was significantly associated with THR in cases of nodal OA, but not in cases of non-nodal OA, and increased height was a risk factor for THR in non-nodal OA cases only. Female sex was a protective risk factor for TKR in non-nodal OA cases (OR 0.60, 95% CI 0.52-0.70) but was predisposing for TKR in the nodal form of OA (OR 1.83, 95% CI 1.49-2.26). The nodal phenotype was associated with a significantly higher risk of undergoing both THR and TKR (OR 1.46, 95% CI 1.09-1.94) and also a significantly higher risk of bilateral TKR (OR 1.70, 95% CI 1.37-2.11), but, paradoxically, was associated with a lower risk of bilateral THR (OR 0.72, 95% CI 0.56-0.91).

CONCLUSION

Nodal and non-nodal forms of large joint OA have significantly different risk factors and outcomes, indicating a different etiology for the 2 forms of OA. With regard to the likelihood of undergoing THR, this appears to be, at least in part, genetically determined.

Role of the Nav1.7 R1150W amino acid change in susceptibility to symptomatic knee osteoarthritis and multiple regional pain

OBJECTIVE

To assess the genetic association of pain in patients with knee osteoarthritis (OA) and those with multiple regional pain with the R1150W variant in the α-subunit of the voltage-gated sodium channel Na(V)1.7.

METHODS

Knee OA patients from 2 UK cohorts (1,411 from the Genetics of Osteoarthritis and Lifestyle study and 267 from the Hertfordshire Cohort Study; 74% with symptomatic OA) with Western Ontario and McMaster Universities OA Index (WOMAC) pain scores were genotyped for rs6746030 (encoding the R1150W change). One hundred seventy-six knee OA patients (53% symptomatic) from the Clearwater Osteoarthritis Study were also tested. A total of 4,295 samples (both affected and unaffected OA) from all 3 studies with data on multiple regional pain were tested. Fixed-effects meta-analyses were carried out with the WOMAC, symptomatic OA (adjusting for radiographic severity), and multiple regional pain as outcomes.

RESULTS

No association with the WOMAC was seen in the UK cohorts. Overall, the meta-analysis of WOMAC yielded a summary statistic of β = 0.47 (95% confidence interval [95% CI] 0.04, 0.89; P = 0.030) for the variant allele. The meta-analysis of symptomatic versus asymptomatic OA did not demonstrate an association with rs6746030 (odds ratio [OR] 0.90 [95% CI 0.71, 1.15], P = 0.38). The meta-analysis of multiple regional pain resulted in a significant OR of 1.40 (95% CI 1.08, 1.80; P = 0.0085). No interstudy heterogeneity was seen for any of the analyses.

CONCLUSION

We find evidence that the R1150W amino acid change in the Na(V)1.7 α-chain is associated with multiple regional pain. This variant is confirmed to be involved in genetic susceptibility to pain, but it does not appear to have a major role in OA-specific pain.

Anatomically corresponded regional analysis of cartilage in asymptomatic and osteoarthritic knees by statistical shape modelling of the bone

Magnetic resonance imaging (MRI) is emerging as the method of choice for measuring cartilage loss in osteoarthritis (OA), but current methods of analysis are imperfect for therapeutic clinical trials. In this paper, we present and evaluate, in two multicenter multivendor studies, a new method for anatomically corresponded regional analysis of cartilage (ACRAC) that allows analysis of knee cartilage morphology in anatomically corresponding focal regions defined on the bone surface. In our first study, 3-D knee MR Images were obtained from 19 asymptomatic female volunteers, followed by segmentations of the bone and cartilage. Minimum description length (MDL) statistical shape models (SSMs) were constructed from the segmented bone surfaces, providing mean bone shapes and a dense set of anatomically corresponding positions on each individual bone, the accuracy of which were measured using repeat images from a subset of the volunteers. Cartilage thicknesses were measured at these locations along 3-D normals to the bone surfaces, yielding corresponded cartilage thickness maps. Functional subregions of the joint were defined on the mean bone shapes, and propagated, using the correspondences, to each individual. ACRAC improved reproducibility, particularly in the central, load bearing subregions of the joint, compared with measures of volume obtained directly from the segmented cartilage surfaces. In our second study, MR Images were obtained from 31 female patient-volunteers with knee OA at baseline and six months. We obtained manual segmentations of the cartilage, and automatic segmentations of the bone using active appearance models (AAMs) built from the bone SSMs of the first study. ACRAC enabled the detection of significant thickness loss in the central, load-bearing regions of the whole femur (-5.57% p = 0.01, annualized) and the medial condyle (-13.08% , p = 0.024 Bonferroni corrected, annualized). We conclude that statistical shape modelling of bone surfaces defines correspondences invariant to individual joint size or shape, providing focal measures of cartilage with improved reproducibility compared to whole compartment measures. It permits the identification of anatomically equivalent regions, and provides the ability to identify the main load-bearing regions of the joint, based on the imputed premorbid state. The method permitted detection of tiny morphological change in cartilage thickness over six months in a small study, and may be useful for OA disease analysis and treatment monitoring.

Self-reported knee and foot alignments in early adult life and risk of osteoarthritis

OBJECTIVE

To determine whether self-reported early adult life malalignment of knees or feet are risk factors for knee or hip osteoarthritis (OA).

METHODS

Participants in the Genetics of Osteoarthritis and Lifestyle case-control database were sent a questionnaire (n = 3,022) containing line-drawing instruments for self-reported knee and foot alignment at ages 20-29 years. Respondents were categorized as having straight, valgus, or varus knee, and straight, toe-in, or toe-out feet. Radiographic criteria were used to define current isolated knee or hip OA, combined knee and hip OA, or non-OA controls. Odds ratios (ORs), adjusted ORs, and 95% confidence intervals (95% CIs) were calculated and logistic regression was performed.

RESULTS

The response rate was 72%; 87.5% of responders (n = 1,901) completed the alignment questions. Increased risk of isolated knee OA occurred with early adult varus (adjusted OR 5.16, 95% CI 2.87-9.41) and valgus knees (adjusted OR 3.16, 95% CI 1.04-9.64). The positive association between knee OA and toe-in foot was explained by varus knee. There was an increased risk of combined knee/hip OA from varus (adjusted OR 4.52, 95% CI 2.39-8.53) and valgus knees (adjusted OR 3.07, 95% CI 0.99-9.54). Varus knee was associated with risk of medial tibiofemoral OA, whereas valgus knee was associated with risk of lateral tibiofemoral and lateral patellofemoral OA. Toe-out foot was associated with reduced medial patellofemoral OA. For knee OA, a multiplicative interaction analysis between occupational risks and varus/valgus yielded an OR(int) of 3.20 (95% CI 1.08-9.49).

CONCLUSION

Constitutional alignment of the leg in terms of varus or valgus knee or foot rotation may be a significant factor in determining development and distribution of knee but not hip OA.

Mild acetabular dysplasia and risk of osteoarthritis of the hip: a case-control study

OBJECTIVE

To determine whether mild variation in acetabular depth (AD) and shape is a risk factor for osteoarthritis (OA) of the hip.

METHODS

The unaffected contralateral hip of patients with unilateral hip OA was compared with hips of asymptomatic controls without hip OA, derived from the Nottingham Genetics Osteoarthritis and Lifestyle case-control study. Standardised anteroposterior x-rays of the pelvis were used to measure centre edge (CE) angle and AD. Cut-off points for narrow CE angle and shallow AD were calculated from the control group (mean -1.96 × SD). The relative risk of hip OA associated with each feature was estimated using OR and 95% CI and adjusted risks were calculated by logistic regression.

RESULTS

In controls, both the CE angle and the AD were lower in the left hip than in the right hip. The CE angle related to age in both hips, and AD of the right hip was lower in men than in women. The contralateral unaffected hip in patients with unilateral hip OA had a decreased CE angle and AD compared with controls, irrespective of side. The lowest tertile of the CE angle in contralateral hips was associated with an eightfold risk of OA (aOR 8.06, 95% CI 4.87 to 13.35) and the lowest tertile of AD was associated with a 2.5-fold risk of OA (aOR 2.53, 95% CI 1.28 to 5.00). Significant increases in the risk of OA were also found as the CE angle and AD decreased.

CONCLUSION

Constitutional mild acetabular dysplasia appears to increase the risk of hip OA.

A meta-analysis of interleukin-6 promoter polymorphisms on risk of hip and knee osteoarthritis

OBJECTIVE

Interleukin-6 is a pro-inflammatory cytokine involved in the pathogenesis of osteoarthritis (OA). We investigated the role of two single nucleotide polymorphisms (SNPs) mapping to the promoter of the IL-6 gene on genetic susceptibility to hip and knee OA.

METHODS

The -174G/C (rs1800795) and -597G/A (rs1800797) SNPs, implicated in the literature in risk of hip and hand OA, were genotyped in 2511 controls, 1101 hip OA cases and 1904 knee OA cases from four cohorts from the UK and Estonia. Data were analysed in conjuntion with published data on rs1800797 from the Genetics of OA and Lifestyle study (UK) on 791 controls, 1034 knee and 997 hip OA cases and rs1800795 data on 75 hip OA cases and 96 controls from Italy. Cases included both radiographic OA only and radiographic and symptomatic OA. Fixed and random-effects meta-analysis models were tested.

RESULTS

No significant association was found with hip OA or knee OA with either SNP nor with the haplotypes formed by them. For individual SNPs the smallest P-value for hip OA was observed using a random-effects model for rs1800795 OR(Gallele)=1.066 (95% CI 0.89-1.28) P<0.49, and significant heterogeneity between cohorts (I(2)=65%, P<0.034) was detected. For knee OA the smallest P-value was seen for rs1800797 OR(Aallele)=1.055 (95%CI 0.98-1.12) P<0.18, no significant heterogeneity was observed (I(2)=0%, P<0.68).

CONCLUSIONS

Our data do not support a role for the -174 and -597 IL-6 promoter polymorphisms in genetic susceptibility to knee or hip OA in Caucasian populations.

Effects of a metalloproteinase inhibitor on osteochondral angiogenesis, chondropathy and pain behavior in a rat model of osteoarthritis

OBJECTIVE

To investigate the effects of a matrix metalloproteinase (MMP) inhibitor on joint pathology and pain behavior in the rat meniscal transection (MNX) model of osteoarthritis (OA) and evaluate which aspects of structural disease modification contribute to symptom improvement.

METHODS

OA pathology was induced in male Lewis rats, by transecting the medial collateral ligament with (MNX) or without (SHAM) a full thickness cut through the meniscus. MNX animals were orally administered an equipotent MMP 2, 8, 9, 12, 13 inhibitor (0.25, 1 and 5 mg/kg/day) or vehicle from day 1. Chondropathy, osteophytosis, osteochondral vascularity were assessed from toluidine blue stained coronal sections of the total knee joint and weight-bearing asymmetry by incapacitance. Group differences were evaluated using 1-way analysis of variance (ANOVA) and associations as Spearman's correlation coefficients.

RESULTS

Treatment with the MMP inhibitor reduced weight-bearing asymmetry from day 14 onwards, and attenuated chondropathy (both P<0.05). Osteochondral vascularity was elevated in MNX compared with SHAM-operated animals (P<0.001) and reduced, dose-dependently, by MMP inhibitor treatment (r=-0.89, P<0.05). Reduced osteochondral vascularity and chondropathy were associated with the amelioration of weight-bearing asymmetry (both P<0.05).

CONCLUSION

Here we show that treatment with a MMP inhibitor reduces joint damage, osteochondral angiogenesis and behavioral evidence of pain. The association between osteochondral angiogenesis and pain behavior may be explained by perivascular nerve growth or stimulation of subchondral nerves following loss of osteochondral integrity. Our data suggest that targeting angiogenesis may have utility in the treatment of pain associated with structural damage in OA.

Measurement and visualisation of focal cartilage thickness change by MRI in a study of knee osteoarthritis using a novel image analysis tool

We describe the application of a novel analysis method that provides detailed maps of changes in cartilage thickness measured from MRI scans for individuals and cohorts of patients together with regional measures. A cohort of osteoarthritis patients was imaged using a 1.0 T MR scanner over a 36-month period. Hyaline cartilage was manually segmented from a three-dimensional (3D) spoiled gradient-echo sequence with fat suppression. Representative outlines of the bone surfaces of the distal femur and proximal tibia were automatically generated from T₂ weighted images using statistical models of the shape and appearance of the bones. Cartilage thickness was measured from a dense set of points representing the bony surface. The models of the bones provided a common frame of reference, relative to which change maps were generated and aggregated across the cohort and anatomically corresponding subregions of the joint to be identified. In the reproducibility arm involving six patients, the thickness of cartilage had coefficients of variation of 2.66% within the tibiofemoral joint and 2.94% within the medial femoral condyle region. In the 9 patients (6 female, 3 male) who completed the 36-month study, the most striking observation was that lack of change in global measures of cartilage thickness concealed substantial focal changes. Specifically, the cartilage thickness within the tibiofemoral joint decreased by 0.85% per annum (95% CI -2.13% to 0.45%) with the medial femoral condyle as the region with the most significant change, decreasing by 2.43% per annum (uncorrected 95% CI -4.31% to 0.51%).

Characterization of Fibrodysplasia in the Dog Following Inhibition of Metalloproteinases

Musculoskeletal side effects are a widely reported consequence of administration of particular matrix metalloproteinase inhibitors (MMPi) in clinical trials. We describe here histopathological findings during dog studies with a fairly selective MMPi AZM551248, that are consistent with these human clinical changes. They were characterized by a dose-and time-dependent formative connective tissue alteration we have termed “fibrodysplasia.” The most sensitive site was the subcuticular connective tissue, although musculoskeletal tissues were also extensively involved. In the subcutis, changes occurred initially around pre-existing blood vessels, but then more diffusely. There was proliferation of cells showing myofibroblast differentiation identified by elevated levels of alpha-smooth muscle actin, fibronectin, and transforming growth factor β, and the deposition of collagen type III with a lesser quantity of collagen type I. On longer-term administration at lower doses, there was evidence of active fibrodysplasia arising and resolving during the dosing period, resulting in the multifocal deposition of mature collagen. Although there was organ specificity, essentially identical changes occurred at multiple connective tissue sites. We conclude that MMPi-induced fibrodysplasia in animals and, by inference, musculoskeletal side effects in humans are potentially diffuse connective tissue disorders.

Attempt to replicate published genetic associations in a large, well-defined osteoarthritis case-control population (the GOAL study)

OBJECTIVE

Published studies have tested over 90 genes for association with osteoarthritis (OA), but few positives reported have been independently replicated. Using a new case-control study, our aim was to attempt the replication of findings from 12 genes reported to have significant genetic association with OA and to further examine the role of genetic variation in six of these genes.

METHODS

A case-control study was undertaken in Nottingham, UK. Hospital-referred index cases with symptomatic, radiographic OA (ROA) of the knee (n=1040) or hip (n=1004) were recruited. Asymptomatic controls (n=1123) were recruited from intravenous urography waiting lists and screened for radiographic hip and knee OA. Sixty-eight polymorphisms were genotyped in IL1A, IL1B, IL1RN, IL4R, IL6, COL2A1, ADAM12, ASPN, IGF1, TGFB1, ESR1 and VDR. Statistical analysis compared allele or genotype frequencies of these polymorphisms in all asymptomatic controls and the subset of controls without ROA vs all OA, knee OA and hip OA. The analyses were adjusted for age, gender and body mass index.

RESULTS

We were unable to replicate any of the published genetic associations investigated. Our extended exploratory analyses identified some associations between polymorphisms in TGFB1, IGF1 and IL1RN and OA; but the strength of evidence varied with the control group used.

CONCLUSION

Lack of replication is common and could be due to differences in study design, phenotype, populations examined or the occurrence of false positives in the initial study. Variants within TGFB1, IGF1 and IL1RN could have a role in OA susceptibility; however, replication of these findings is required in an independent study.

Evaluation of a magnetic resonance biomarker of osteoarthritis disease progression: doxycycline slows tibial cartilage loss in the Dunkin Hartley guinea pig

The objective was to assess the effect of doxycycline treatment on a magnetic resonance imaging (MRI) biomarker of cartilage volume loss, and on matrix metalloproteinase (MMP) activity in a guinea pig osteoarthritis model. Guinea pigs (9 months old) were dosed with vehicle or doxycycline, 0.6, 3.0 mg/kg/day for 66 days. Fat-suppressed 3D gradient-echo MRI of the left knee was acquired pre- and post dosing. Change in medial tibial plateau (MTP) cartilage volume (MT.VC) was determined using image analysis. At termination, MTP cartilage was removed from knees and proteolytic MMP activity determined using a fluorescent peptide substrate assay. Vehicle-treated animals lost 20.5% (95% CI mean 25.6-15.1) MT.VC. The doxycycline (0.6 mg/kg/day) group lost 8.6% (P < 0.05, 95% CI 20.6 to -5.3) whilst the 3.0 mg/kg/day group lost 10.0% (P < 0.05, 95% CI 13.9-6.0%). Endogenous levels of active MMPs were below limits of detection in all samples. However, doxycycline treatment ablated amino phenyl mercuric acid activated MMP-13 and MMP-8 levels, reduced MMP-9 levels by 65% and MMP-1 levels by 24%. Doxycycline treatment resulted in partial protection from MT.VC loss and was associated with complete reduction in MMP-13 and MMP-8, and partial reduction in MMP-9 activity. These data imply a role of MMPs in cartilage degeneration but incomplete protection suggests that additional doxycycline insensitive mechanisms are important in this model. The protective effect of doxycycline correlates with the clinical finding of lessened joint space narrowing, strengthens the utility of this animal model in identifying disease-modifying osteoarthritic drugs and supports the use of MRI biomarkers of cartilage loss.

Bone marrow lesions from osteoarthritis knees are characterized by sclerotic bone that is less well mineralized

Introduction

Although the presence of bone marrow lesions (BMLs) on magnetic resonance images is strongly associated with osteoarthritis progression and pain, the underlying pathology is not well established. The aim of the present study was to evaluate the architecture of subchondral bone in regions with and without BMLs from the same individual using bone histomorphometry.

Methods

Postmenopausal female subjects (n = 6, age 48 to 90 years) with predominantly medial compartment osteoarthritis and on a waiting list for total knee replacement were recruited. To identify the location of the BMLs, subjects had a magnetic resonance imaging scan performed on their study knee prior to total knee replacement using a GE 1.5 T scanner with a dedicated extremity coil. An axial map of the tibial plateau was made, delineating the precise location of the BML. After surgical removal of the tibial plateau, the BML was localized using the axial map from the magnetic resonance image and the lesion excised along with a comparably sized bone specimen adjacent to the BML and from the contralateral compartment without a BML. Cores were imaged via microcomputed tomography, and the bone volume fraction and tissue mineral density were calculated for each core. In addition, the thickness of the subchondral plate was measured, and the following quantitative metrics of trabecular structure were calculated for the subchondral trabecular bone in each core: trabecular number, thickness, and spacing, structure model index, connectivity density, and degree of anisotropy. We computed the mean and standard deviation for each parameter, and the unaffected bone from the medial tibial plateau and the bone from the lateral tibial plateau were compared with the affected BML region in the medial tibial plateau.

Results

Cores from the lesion area displayed increased bone volume fraction but reduced tissue mineral density. The samples from the subchondral trabecular lesion area exhibited increased trabecular thickness and were also markedly more plate-like than the bone in the other three locations, as evidenced by the lower value of the structural model index. Other differences in structure that were noted were increased trabecular spacing and a trend towards decreased trabecular number in the cores from the medial location as compared with the contralateral location.

Conclusions

Our preliminary data localize specific changes in bone mineralization, remodeling and defects within BMLs features that are adjacent to the subchondral plate. These BMLs appear to be sclerotic compared with unaffected regions from the same individual based on the increased bone volume fraction and increased trabecular thickness. The mineral density in these lesions, however, is reduced and may render this area to be mechanically compromised, and thus susceptible to attrition.

Nonspherical femoral head shape (pistol grip deformity), neck shaft angle, and risk of hip osteoarthritis: a case-control study

OBJECTIVE

To determine whether 2-dimensional measures of femoral head shape and angle are associated with hip osteoarthritis (OA).

METHODS

We compared cases with symptomatic radiographic hip OA with asymptomatic controls with no radiographic hip OA. On anteroposterior pelvis radiographs, we measured "pistol grip deformity" for each hip (visually categorized as nonspherical, indeterminate, or spherical), the femoral head-to-femoral neck ratio as an interval measure of femoral head shape, and the femoral neck shaft angle. The relative risk of hip OA associated with each feature was estimated using odds ratios (ORs) and 95% confidence intervals (95% CIs), adjusted for possible confounders using a logistic regression model.

RESULTS

Of 1,007 cases, 965 had definite radiographic hip OA; of 1,123 controls, 1,111 had no radiographic OA. The prevalence of pistol grip deformity in at least 1 hip was 3.61% in controls and 17.71% in cases (OR 6.95 [95% CI 4.64-10.41]), and the prevalence of abnormal femoral head-to-femoral neck ratio in at least 1 hip was 3.70% in controls and 24.27% in cases (OR 12.08 [95% CI 8.05-18.15]). The risk of hip OA increased as the femoral head-to-femoral neck ratio decreased (P for trend<0.001) and with each extreme of neck shaft angle (P<0.05). In cases with unilateral hip OA, the prevalence of abnormal femoral head-to-femoral neck ratio in the unaffected hip was 2 times greater than that in controls (OR 1.82 [95% CI 1.07-3.07]); in contrast, an abnormally low, but not abnormally high, neck shaft angle was more common in unaffected hips than in controls (OR 1.79 [95% CI 1.03-3.14]).

CONCLUSION

Our findings indicate that pistol grip deformity is associated with hip OA. The increased prevalence of pistol grip deformity and an abnormally low neck shaft angle in unaffected hips of cases with unilateral OA suggests that they are risk factors for development of hip OA. However, both a nonspherical head shape and an increase in neck shaft angle may occur as a consequence of OA.

Index to ring finger length ratio and the risk of osteoarthritis

OBJECTIVE

To determine the relationship between the index to ring finger (2D:4D) length ratio and the risk of knee and hip osteoarthritis (OA).

METHODS

We conducted a case-control study, in which cases with persistent symptoms and radiographic evidence of knee or hip OA were compared with controls with no symptoms and no radiographic evidence of knee or hip OA. Hand radiographs were visually classified as type 1 (index finger longer than the ring finger), type 2 (index finger equal to the ring finger), or type 3 (index finger shorter than the ring finger). The 2D:4D phalangeal and metacarpal length ratios were measured separately. The odds ratio (OR) and 95% confidence interval (95% CI) were calculated and adjusted for possible confounding factors using a logistic regression model.

RESULTS

Of 2,049 cases, 1,013 had radiographic evidence of knee OA and 995 had hip OA. Of 1,123 controls, 836 had no knee OA and 1,050 had no hip OA. The type 3 finger pattern was associated with knee OA (OR 1.94, 95% CI 1.54-2.44), and the risk was greater in women (OR 3.05, 95% CI 2.08-4.47) than in men (OR 1.45, 95% CI 1.08-1.95). There was a dose-response relationship between both 2D:4D phalangeal and metacarpal length ratios and the risk of knee OA. The risk of hip OA was inconsistent.

CONCLUSION

Compared with types 1 and 2, the type 3 "male" pattern 2D:4D length ratio is associated with OA, especially knee OA. The risk is independent of other major OA risk factors.

Radiographic assessment of the index to ring finger ratio (2D:4D) in adults

The smaller index to ring finger (2D:4D) ratio has been considered as a 'male finger pattern' and is associated with sporting ability and a number of conditions. However, the ratio may vary according to what is measured, the hand selected and the method used. This study aimed to determine: (1) which bones (phalanges, metacarpals or both) account for variation in the 2D:4D ratio; (2) whether the ratio shows right-left symmetry or relates to hand dominance; and (3) the correlation between visual classification and measured determinations of the ratio based on radiographs. Hand radiographs obtained as part of a large osteoarthritis genetic study were examined. Each hand was classified visually into three types according to the relative length of the index and ring finger: Type 1 (index longer than ring), Type 2 (index = ring) and Type 3 (index shorter than ring). For both index and ring fingers we measured (1) from base of proximal to tip of distal phalanx and (2) metacarpal length. Reproducibility of the classification and measurements were examined using kappa and intraclass correlation coefficient; symmetry between left and right hands was examined using Bland and Altman's agreement analysis; and correlation between visual classification and 2D:4D ratio data was analysed using the anova linearity test. Data were obtained from 3172 radiographs (1636 men, 1536 women; mean age 67 +/- 7.9 years, range 45-86 years). Prevalence of Type 3 hand was 61% in men and 37% in women (P < 0.001). Men had smaller 2D:4D ratios than women for phalanges (0.908 versus 0.922, P < 0.01), metacarpals (1.152 versus 1.157, P < 0.01) and the sum of phalanges plus metacarpals (1.005 versus 1.015, P < 0.01). The mean difference between right and left was -0.001 (95% limit of agreement -0.035, 0.032) for the phalangeal ratio and 0.003 (95% limit of agreement -0.051 to 0.057) for the metacarpal ratio. The 2D:4D ratio did not associate with handedness or age. There was a linear trend between the visual classification of hand type and the 2D:4D ratio data (P < 0.001). More technical difficulties (due to positioning, finger trauma, osteoarthritis) were encountered with the phalangeal ratio and visual categorization than with the metacarpal ratio: the latter could be measured in 98.7% of the study population. We concluded that measured 2D:4D ratios and visual categorization can be derived from hand radiographs. The phalanges and metacarpals both contribute to the variation in 2D:4D ratio with smaller ratios observed in men than in women. The ratio is symmetrical with only very small differences between right and left hands. Visual classification may be a useful simple tool for future epidemiological studies but is more prone to bias from positioning than direct measurement. If radiographs are used for this purpose, we recommend the metacarpal ratio with measurement of a single index hand or an average of both as it is least affected by bias from malpositioning, trauma or common joint disease.

Measurement of cartilage volumes in rheumatoid arthritis using MRI

MRI is a valuable imaging modality for assessment of the articular cartilage in rheumatoid arthritis (RA) and is potentially of use in monitoring disease progression and response to therapy. In this study, we investigated the sources of error in volume measurements obtained by segmentation of MR images of knee cartilage in patients with RA and followed cartilage volume in a group of RA patients for 12 months. 23 RA patient volunteers were recruited for knee imaging. Six subjects were imaged at baseline only, six were imaged at baseline and again within an hour in the same imaging session, six subjects were imaged at baseline and 7 days, and 17 subjects were imaged at baseline, 4+/-2 months and 12 months. Imaging was performed at 1.0 T using a three-dimensional spoiled gradient-echo sequence with fat-suppression. Manual image segmentation was performed once or twice on the lateral tibial, medial tibial, patellar and femoral compartment by either one or two segmenters. Coefficients of variation (CoV) for repeated volume measurement of total cartilage were 2.2% (same segmenter, same scan), 5.2% (different segmenter, same scan), 4.9% (same segmenter, different scan, same session), and 4.4% (same segmenter, different scan, different session). Over the 12 month duration of the study there was no significant change in total cartilage volume, nor were there significant changes in volume in any individual compartment. This measurement technique is reproducible, but any net change in cartilage volume over 1 year is very small.

Integrin alpha5beta1 and ADAM-17 interact in vitro and co-localize in migrating HeLa cells

Tumor necrosis factor (TNF) alpha-converting enzyme (TACE/ADAM-17) has diverse roles in the proteolytic processing of cell surface molecules and, due to its ability to process TNFalpha, is a validated therapeutic target for anti-inflammatory therapies. Unlike a number of other ADAM proteins, which interact with integrin receptors via their disintegrin domains, there is currently no evidence for an ADAM-17-integrin association. By analyzing the adhesion of a series of cell lines with recombinant fragments of the extracellular domain of ADAM-17, we now demonstrate a functional interaction between ADAM-17 and alpha(5)beta(1) integrin in a trans orientation. Because ADAM-17-mediated adhesion was sensitive to RGD peptides and EDTA, and the integrin-binding site within ADAM-17 was narrowed down to the disintegrin/cysteine-rich region, the two molecules appear to have a ligand-receptor relationship mediated by the alpha(5)beta(1) ligand binding pocket. Intriguingly, ADAM-17 and alpha(5)beta(1) were found to co-localize in both membrane ruffles and focal adhesions in HeLa cells. When confluent HeLa cell monolayers were wounded, ADAM-17 and alpha(5)beta(1) redistributed to the leading edge and co-localized, which is suggestive of a cis orientation. We postulate that the interaction of ADAM-17 with alpha(5)beta(1) may target or modulate its metalloproteolytic activity.

Magnetic resonance imaging measurement of knee cartilage volume in a multicentre study

OBJECTIVE

To investigate the variability between different high-field scanners in magnetic resonance imaging (MRI) measurement of knee cartilage volume in healthy female volunteers.

METHODS

Five volunteers had both knees scanned using three different MRI scanners. Cartilage volume in each compartment was measured from the images by image segmentation. The data were analysed using analysis of variance models.

RESULTS

The mean total cartilage volume of the 10 knees scanned at three different centres was 16.15, 16.40 and 15.63 ml for the Siemens, GE and Philips scanners respectively. Small systematic differences were seen in the total knee cartilage volume results.

CONCLUSIONS

Although there were small systematic differences in knee cartilage volume, the three MRI scanners gave broadly similar results.

Characterisation of the guinea pig model of osteoarthritis by in vivo three-dimensional magnetic resonance imaging

OBJECTIVE

To characterise longitudinal changes in joint integrity and cartilage volume in vivo in the guinea pig spontaneous osteoarthritis (OA) model by magnetic resonance imaging (MRI).

METHODS

Guinea pigs knee were imaged in vivo by high-resolution three-dimensional (3D) MRI between the ages of 3 and 12 months. Image analysis was performed to assess qualitative knee joint changes between 3 and 12 months (n=16) and quantitative volumetric changes of the medial tibial cartilage between 9 and 12 months (n=7). After imaging, animals were killed and knees were assessed macroscopically and histologically.

RESULTS

From 3 to 6 months qualitative observation by MRI and histopathology indicated localised cartilage swelling on the medial tibial plateau. At 6 months, bone cysts had developed in the epiphysis. At 9 months, we observed by MRI and histopathology, fragmentation of the medial tibial cartilage in areas not protected by the meniscus. Cartilage degeneration had intensified at 12 months with evidence of widespread loss of cartilage throughout the tibial plateau. Segmentation of the MR cartilage images showed a 36% loss of volume between 9 and 12 months.

CONCLUSIONS

We have achieved 3D image acquisition and segmentation of knee cartilage in a guinea pig model of chronic OA, which permits measurements previously only possible in man. High resolution and short acquisition time allowed qualitative longitudinal characterisation of the entire knee joint and enabled us to quantify for the first time longitudinal tibial cartilage volume loss associated with disease progression.

No loss of cartilage volume over three years in patients with knee osteoarthritis as assessed by magnetic resonance imaging

OBJECTIVE

Magnetic resonance imaging (MRI) has the potential to provide accurate quantification of structural changes in joint disease, with sensitivity to change, as it can provide direct visualization of the cartilage and bone. In this study, we investigated whether knee cartilage volume, as assessed by MRI, is sensitive to change over time in patients with osteoarthritis (OA).

DESIGN

Sixteen patient volunteers (10 male, six female) with established OA of the knee were entered into the study and demographic data recorded. At baseline, 12 months and 37+/-2 months, patients underwent simple measures of disease severity, as well as extended weight-bearing AP knee X-rays. In addition the patient's index knee was imaged using MR at 1.0 T using a 3-D spoiled gradient-echo sequence with fat-suppression, repetition time 50 ms, echo time 11 ms, flip-angle 40 degrees, sagittal slice thickness 1.56 mm and in-plane pixel resolution 0.55 mm. Manual image segmentation was performed on all knee cartilage compartments and the respective cartilage volumes determined.

RESULTS

Eleven of the original patients recruited completed the 3-year study. Radiographic features indicated that the majority had a spectrum of well-established OA at entry. The average decrease in medial tibiofemoral joint space width was 0.21+/-0.37 mm (mean+/-S.D.). Comparison of MR images at baseline and 37+/-2 months indicated little evidence of cartilage lesion shape or size change in any of the compartments. There was no significant MRI volume change in any of the knee cartilage compartments over the course of 1 year. The change in total knee cartilage volume, as measured by MRI, was a loss of only 1.6%, or 0.36+/-1.3 ml (mean+/-S.D.), over the 3 years.

CONCLUSIONS

The failure to identify loss of cartilage volume over 3 years in this cohort of patients with established knee OA using MRI challenges the face validity of this endpoint to assess structural changes in OA.