Interferon-α-mediated therapeutic resistance in early rheumatoid arthritis implicates epigenetic reprogramming.

OBJECTIVES: An interferon (IFN) gene signature (IGS) is present in approximately 50% of early, treatment naive rheumatoid arthritis (eRA) patients where it has been shown to negatively impact initial response to treatment. We wished to validate this effect and explore potential mechanisms of action. METHODS: In a multicentre inception cohort of eRA patients (n=191), we examined the whole blood IGS (MxA, IFI44L, OAS1, IFI6, ISG15) with reference to circulating IFN proteins, clinical outcomes and epigenetic influences on circulating CD19+ B and CD4+ T lymphocytes. RESULTS: We reproduced our previous findings demonstrating a raised baseline IGS. We additionally showed, for the first time, that the IGS in eRA reflects circulating IFN-α protein. Paired longitudinal analysis demonstrated a significant reduction between baseline and 6-month IGS and IFN-α levels (p<0.0001 for both). Despite this fall, a raised baseline IGS predicted worse 6-month clinical outcomes such as increased disease activity score (DAS-28, p=0.025) and lower likelihood of a good EULAR clinical response (p=0.034), which was independent of other conventional predictors of disease activity and clinical response. Molecular analysis of CD4+ T cells and CD19+ B cells demonstrated differentially methylated CPG sites and dysregulated expression of disease relevant genes, including PARP9, STAT1, and EPSTI1, associated with baseline IGS/IFNα levels. Differentially methylated CPG sites implicated altered transcription factor binding in B cells (GATA3, ETSI, NFATC2, EZH2) and T cells (p300, HIF1α). CONCLUSIONS: Our data suggest that, in eRA, IFN-α can cause a sustained, epigenetically mediated, pathogenic increase in lymphocyte activation and proliferation, and that the IGS is, therefore, a robust prognostic biomarker. Its persistent harmful effects provide a rationale for the initial therapeutic targeting of IFN-α in selected patients with eRA.

Cryptic splicing: common pathological mechanisms involved in male infertility and neuronal diseases.

High levels of transcription and alternative splicing are recognized hallmarks of gene expression in the testis and largely driven by cells in meiosis. Because of this, the male meiosis stage of the cell cycle is often viewed as having a relatively permissive environment for gene expression. In this review, we highlight recent findings that identify the RNA binding protein RBMXL2 as essential for male meiosis. RBMXL2 functions as a "guardian of the transcriptome" that protects against the use of aberrant (or "cryptic") splice sites that would disrupt gene expression. This newly discovered protective role during meiosis links with a wider field investigating mechanisms of cryptic splicing control that protect neurons from amyotrophic lateral sclerosis and Alzheimer's disease. We discuss how the mechanism repressing cryptic splicing patterns during meiosis evolved, and why it may be essential for sperm production and male fertility.

A Bayesian network approach incorporating imputation of missing data enables exploratory analysis of complex causal biological relationships.

Bayesian networks can be used to identify possible causal relationships between variables based on their conditional dependencies and independencies, which can be particularly useful in complex biological scenarios with many measured variables. Here we propose two improvements to an existing method for Bayesian network analysis, designed to increase the power to detect potential causal relationships between variables (including potentially a mixture of both discrete and continuous variables). Our first improvement relates to the treatment of missing data. When there is missing data, the standard approach is to remove every individual with any missing data before performing analysis. This can be wasteful and undesirable when there are many individuals with missing data, perhaps with only one or a few variables missing. This motivates the use of imputation. We present a new imputation method that uses a version of nearest neighbour imputation, whereby missing data from one individual is replaced with data from another individual, their nearest neighbour. For each individual with missing data, the subsets of variables to be used to select the nearest neighbour are chosen by sampling without replacement the complete data and estimating a best fit Bayesian network. We show that this approach leads to marked improvements in the recall and precision of directed edges in the final network identified, and we illustrate the approach through application to data from a recent study investigating the causal relationship between methylation and gene expression in early inflammatory arthritis patients. We also describe a second improvement in the form of a pseudo-Bayesian approach for upweighting certain network edges, which can be useful when there is prior evidence concerning their directions.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs.

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Histone ChIP-Seq identifies differential enhancer usage during chondrogenesis as critical for defining cell-type specificity.

Epigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterized the epigenome during the in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) was used to assess a range of N-terminal posttranscriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3, and H3K36me3) in both hMSCs and differentiated chondrocytes. Chromatin states were characterized using histone ChIP-seq and cis-regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis-related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de-methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this study with epigenomes of cell-types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell-types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers. Luciferase reporter assays confirmed that chondrocyte enhancers characterized in this study exhibited enhancer activity which may be modulated by DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross-talk between different epigenetic mechanisms during chondrocyte differentiation.

Lymphocyte DNA methylation mediates genetic risk at shared immune-mediated disease loci.

BACKGROUND: Defining regulatory mechanisms through which noncoding risk variants influence the cell-mediated pathogenesis of immune-mediated disease (IMD) has emerged as a priority in the post-genome-wide association study era. OBJECTIVES: With a focus on rheumatoid arthritis, we sought new insight into genetic mechanisms of adaptive immune dysregulation to help prioritize molecular pathways for targeting in this and related immune pathologies. METHODS: Whole-genome methylation and transcriptional data from isolated CD4+ T cells and B cells of more than 100 genotyped and phenotyped patients with inflammatory arthritis, all of whom were naive to immunomodulatory treatments, were obtained. Analysis integrated these comprehensive data with genome-wide association study findings across IMDs and other publicly available resources. RESULTS: We provide strong evidence that disease-associated DNA variants regulate cis-CpG methylation in CD4+ T and/or B cells at 37% RA loci. Using paired, cell-specific transcriptomic data and causal inference testing, we identify examples where site-specific DNA methylation in turn mediates gene expression, including FCRL3 in both cell types and ORMDL3/GSDMB, IL6ST/ANKRD55, and JAZF1 in CD4+ T cells. A number of genes regulated in this way highlight mechanisms common to RA and other IMDs including multiple sclerosis and asthma, in turn distinguishing them from osteoarthritis, a primarily degenerative disease. Finally, we corroborate the observed effects experimentally. CONCLUSIONS: Our observations highlight important mechanisms of genetic risk in RA and the wider context of immune dysregulation. They confirm the utility of DNA methylation profiling as a tool for causal gene prioritization and, potentially, therapeutic targeting in complex IMD.

DNA hypomethylation during MSC chondrogenesis occurs predominantly at enhancer regions.

Regulation of transcription occurs in a cell type specific manner orchestrated by epigenetic mechanisms including DNA methylation. Methylation changes may also play a key role in lineage specification during stem cell differentiation. To further our understanding of epigenetic regulation in chondrocytes we characterised the DNA methylation changes during chondrogenesis of mesenchymal stem cells (MSCs) by Infinium 450 K methylation array. Significant DNA hypomethylation was identified during chondrogenic differentiation including changes at many key cartilage gene loci. Integration with chondrogenesis gene expression data revealed an enrichment of significant CpGs in upregulated genes, while characterisation of significant CpG loci indicated their predominant localisation to enhancer regions. Comparison with methylation profiles of other tissues, including healthy and diseased adult cartilage, identified chondrocyte-specific regions of hypomethylation and the overlap with differentially methylated CpGs in osteoarthritis. Taken together we have associated DNA methylation levels with the chondrocyte phenotype. The consequences of which has potential to improve cartilage generation for tissue engineering purposes and also to provide context for observed methylation changes in cartilage diseases such as osteoarthritis.

Correlation of Infinium HumanMethylation450K and MethylationEPIC BeadChip arrays in cartilage.

DNA methylation of CpG sites is commonly measured using Illumina Infinium BeadChip platforms. The Infinium MethylationEPIC array has replaced the Infinium Methylation450K array. The two arrays use the same technology, with the EPIC array assaying almost double the number of sites than the 450K array. In this study, we compare DNA methylation values of shared CpGs of the same human cartilage samples assayed using both platforms. DNA methylation was measured in 21 human cartilage samples using the both 450K and EPIC arrays. Additional matched 450K and EPIC data in whole tumour and whole blood were downloaded from GEO GSE92580 and GSE86833, respectively. Data were processed using the Bioconductor package Minfi. DNA methylation of six CpG sites was validated for the same 21 cartilage samples by pyrosequencing. In cartilage samples, overall sample correlations of methylation values between arrays were high (Pearson's r > 0.96). However, 50.5% of CpG sites showed poor correlation (r < 0.2) between arrays. Sites with limited variance and with either very high or very low methylation levels in cartilage exhibited lower correlation values, corroborating prior studies in whole blood. Bisulphite pyrosequencing did not highlight one array as generating more accurate methylation values. For a specific CpG site, the array methylation correlation coefficient differed between cartilage, tumour, and whole blood, reflecting the difference in methylation variance between cell types. Researchers should be cautious when analysing methylation of CpG sites that show low methylation variance within the cell type of interest, regardless of the method used to assay methylation.

Expression analysis of the osteoarthritis genetic susceptibility mapping to the matrix Gla protein gene MGP.

BACKGROUND: Osteoarthritis (OA) is a common disease of older individuals that impacts detrimentally on the quality and the length of life. It is characterised by the painful loss of articular cartilage and is polygenic and multifactorial. Genome-wide association scans have highlighted over 90 osteoarthritis genetic signals, some of which reside within or close to highly plausible candidate genes. An example is an association to polymorphisms within and adjacent to the matrix Gla protein gene MGP. We set out to undertake a functional study of this gene. METHODS: Nucleic acid was extracted from cartilage, infrapatellar fat pad, synovium, trabecular bone, trapezium and peripheral whole blood from OA patients and also from mesenchymal stem cells (MSCs) subjected to chondrogenesis. Expression of MGP was measured by quantitative PCR (qPCR), RNA-sequencing and allelic expression imbalance (AEI) analysis. Matrix Gla protein was depleted in chondrocytes by knocking down MGP expression using RNA interference (RNAi) and the effect on a range of genes assessed by qPCR. RESULTS: MGP is expressed in joint tissues, blood and chondrocytes cultured from MSCs. There is a higher expression in diseased versus non-diseased cartilage. Polymorphisms that are associated with OA also correlate with the expression of MGP, with the OA risk-conferring allele showing significantly reduced expression in cartilage, fat pad and synovium but increased expression in blood. Depletion of Matrix Gla protein had a significant effect on the majority of genes tested, with an increased expression of catabolic genes that encode enzymes that degrade cartilage. CONCLUSIONS: MGP expression is subject to cis-acting regulators that correlate with the OA association signal. These are active in a range of joint tissues but have effects which are particularly strong in cartilage. An opposite effect is observed in blood, highlighting the context-specific nature of the regulation of this gene's expression. Recapitulation of the genetic deficit in cartilage chondrocytes is pro-catabolic.

Prioritization of PLEC and GRINA as Osteoarthritis Risk Genes Through the Identification and Characterization of Novel Methylation Quantitative Trait Loci.

OBJECTIVE: To identify methylation quantitative trait loci (mQTLs) correlating with osteoarthritis (OA) risk alleles and to undertake mechanistic characterization as a means of target gene prioritization. METHODS: We used genome-wide genotyping and cartilage DNA methylation array data in a discovery screen of novel OA risk loci. This was followed by methylation, gene expression analysis, and genotyping studies in additional cartilage samples, accompanied by in silico analyses. RESULTS: We identified 4 novel OA mQTLs. The most significant mQTL contained 9 CpG sites where methylation correlated with OA risk genotype, with 5 of the CpG sites having P values <1 × 10-10 . The 9 CpG sites reside in an interval of only 7.7 kb within the PLEC gene and form 2 distinct clusters. We were able to prioritize PLEC and the adjacent gene GRINA as independent targets of the OA risk. We identified PLEC and GRINA expression QTLs operating in cartilage, as well as methylation-expression QTLs operating on the 2 genes. GRINA and PLEC also demonstrated differential expression between OA hip and non-OA hip cartilage. CONCLUSION: PLEC encodes plectin, a cytoskeletal protein that maintains tissue integrity by regulating intracellular signaling in response to mechanical stimuli. GRINA encodes the ionotropic glutamate receptor TMBIM3 (transmembrane BAX inhibitor 1 motif-containing protein family member 3), which regulates cell survival. Based on our results, we hypothesize that in a joint predisposed to OA, expression of these genes alters in order to combat aberrant biomechanics, and that this is epigenetically regulated. However, carriage of the OA risk-conferring allele at this locus hinders this response and contributes to disease development.

Identification of a novel, methylation-dependent, RUNX2 regulatory region associated with osteoarthritis risk.

Osteoarthritis (OA) is a common, multifactorial and polygenic skeletal disease that, in its severest form, requires joint replacement surgery to restore mobility and to relieve chronic pain. Using tissues from the articulating joints of 260 patients with OA and a range of in vitro experiments, including CRISPR-Cas9, we have characterized an intergenic regulatory element. Here, genotype at an OA risk locus correlates with differential DNA methylation, with altered gene expression of both a transcriptional regulator (RUNX2), and a chromatin remodelling protein (SUPT3H). RUNX2 is a strong candidate for OA susceptibility, with its encoded protein being essential for skeletogenesis and healthy joint function. The OA risk locus includes single nucleotide polymorphisms (SNPs) located within and flanking the differentially methylated region (DMR). The OA association SNP, rs10948172, demonstrates particularly strong correlation with methylation, and two intergenic SNPs falling within the DMR (rs62435998 and rs62435999) demonstrate genetic and epigenetic effects on the regulatory activity of this region. We therefore posit that the OA signal mediates its effect by modulating the methylation of the regulatory element, which then impacts on gene expression, with RUNX2 being the principal target. Our study highlights the interplay between DNA methylation, OA genetic risk and the downstream regulation of genes critical to normal joint function.

Functional Characterization of the Osteoarthritis Genetic Risk Residing at ALDH1A2 Identifies rs12915901 as a Key Target Variant.

OBJECTIVE: To identify the functional single-nucleotide polymorphisms (SNPs) and mechanisms conferring increased risk of hand osteoarthritis (OA) at the ALDH1A2 locus, which is a retinoic acid regulatory gene. METHODS: Tissue samples from 247 patients with knee, hip, or hand OA who had undergone joint surgery were included. RNA-sequencing analysis was used to investigate differential expression of ALDH1A2 and other retinoic acid signaling pathway genes in cartilage. Expression of ALDH1A2 in joint tissues obtained from multiple sites was quantified using quantitative reverse transcription-polymerase chain reaction. Allelic expression imbalance (AEI) was measured by pyrosequencing. The consequences of ALDH1A2 depletion by RNA interference were assessed in primary human chondrocytes. In silico and in vitro analyses were used to pinpoint which, among 62 highly correlated SNPs, could account for the association at the locus. RESULTS: ALDH1A2 expression was observed across multiple joint tissue samples, including osteochondral tissue from the hand. The expression of ALDH1A2 and of several retinoic acid signaling genes was different in diseased cartilage compared to non-diseased cartilage, with ALDH1A2 showing lower levels in OA cartilage. Experimental depletion of ALDH1A2 resulted in changes in the expression levels of a number of chondrogenic markers, including SOX9. In addition, reduced expression of the OA risk-conferring allele was witnessed in a number of joint tissues, with the strongest effect in cartilage. The intronic SNP rs12915901 recapitulated the AEI observed in patient tissues, while the Ets transcription factors were identified as potential mediators of this effect. CONCLUSION: The ALDH1A2 locus seems to increase the risk of hand OA through decreased expression of ALDH1A2 in joint tissues, with the effect dependent on rs12915901. These findings indicate a mechanism that may now be targeted to modulate OA risk.

CD4+ and B Lymphocyte Expression Quantitative Traits at Rheumatoid Arthritis Risk Loci in Patients With Untreated Early Arthritis: Implications for Causal Gene Identification.

OBJECTIVE: Rheumatoid arthritis (RA) is a genetically complex disease of immune dysregulation. This study sought to gain further insight into the genetic risk mechanisms of RA by conducting an expression quantitative trait locus (eQTL) analysis of confirmed genetic risk loci in CD4+ T cells and B cells from carefully phenotyped patients with early arthritis who were naive to therapeutic immunomodulation. METHODS: RNA and DNA were isolated from purified B and/or CD4+ T cells obtained from the peripheral blood of 344 patients with early arthritis. Genotyping and global gene expression measurements were carried out using Illumina BeadChip microarrays. Variants in linkage disequilibrium (LD) with non-HLA RA single-nucleotide polymorphisms (defined as r2 ≥ 0.8) were analyzed, seeking evidence of cis- or trans-eQTLs according to whether the associated probes were or were not within 4 Mb of these LD blocks. RESULTS: Genes subject to cis-eQTL effects that were common to both CD4+ and B lymphocytes at RA risk loci were FADS1, FADS2, BLK, FCRL3, ORMDL3, PPIL3, and GSDMB. In contrast, those acting on METTL21B, JAZF1, IKZF3, and PADI4 were unique to CD4+ lymphocytes, with the latter candidate risk gene being identified for the first time in this cell subset. B lymphocyte-specific eQTLs for SYNGR1 and CD83 were also found. At the 8p23 BLK-FAM167A locus, adjacent genes were subject to eQTLs whose activity differed markedly between cell types; in particular, the FAM167A effect displayed striking B lymphocyte specificity. No trans-eQTLs approached experiment-wide significance, and linear modeling did not identify a significant influence of biologic covariates on cis-eQTL effect sizes. CONCLUSION: These findings further refine the understanding of candidate causal genes in RA pathogenesis, thus providing an important platform from which downstream functional studies, directed toward particular cell types, may be prioritized.

A consistent and potentially exploitable response during chondrogenesis of mesenchymal stem cells from osteoarthritis patients to the protein encoded by the susceptibility gene GDF5.

Osteoarthritis (OA) is a common joint disease characterised by the focal loss of the protective cartilage layer at the ends of the bones. It is painful, disabling, multifactorial and polygenic. The growth differentiation factor 5 gene GDF5 was one of the first reported OA susceptibility signals that showed consistent association to OA, with the transcript single nucleotide polymorphism (SNP) rs143383 demonstrating association in Asians and Europeans. The functional effect of the signal is reduced expression of the gene. The GDF5 protein is an extracellular matrix signalling molecule that is active during chondrogenesis and in mature chondrocytes. Due to the functional impact of the susceptibility, we previously assessed the effect of supplementing chondrocytes from OA patients with exogenous GDF5. Their response was highly discordant, precluding the application of GDF5 as a simple means of attenuating the genetic deficit. Since GDF5 is also active during development, we have now assessed the effect of exogenous GDF5 on bone marrow derived mesenchymal stem cells (MSCs) that are undergoing chondrogenesis during cartilage disc formation. MSCs from healthy donors and OA patients were studied and the effect of GDF5 was assessed by measuring the wet mass of the discs, by histological staining, and by monitoring the change in expression of anabolic, catabolic and hypertrophic protein-coding genes. The MSCs expressed the three principal GDF5 receptor genes and responded in a significantly anabolic manner (increase in wet mass, p = 0.0022; Bonferroni corrected p = 0.018) to a variant form of GDF5 that targets the most abundantly expressed receptor, BMPR-IA. GDF5 elicited significant (p < 0.05) changes in the expression of anabolic, catabolic and hypertrophic genes with several consistent effects in healthy donors and in OA patients. Our data implies that, unlike OA chondrocytes, OA MSCs do respond in a predictable, anabolic manner to GDF5, which could therefore provide a route to modulate the genetic deficit mediated by the rs143383 association signal.

Whole-genome sequencing identifies rare genotypes in COMP and CHADL associated with high risk of hip osteoarthritis.

We performed a genome-wide association study of total hip replacements, based on variants identified through whole-genome sequencing, which included 4,657 Icelandic patients and 207,514 population controls. We discovered two rare signals that strongly associate with osteoarthritis total hip replacement: a missense variant, c.1141G>C (p.Asp369His), in the COMP gene (allelic frequency = 0.026%, P = 4.0 × 10-12, odds ratio (OR) = 16.7) and a frameshift mutation, rs532464664 (p.Val330Glyfs*106), in the CHADL gene that associates through a recessive mode of inheritance (homozygote frequency = 0.15%, P = 4.5 × 10-18, OR = 7.71). On average, c.1141G>C heterozygotes and individuals homozygous for rs532464664 had their hip replacement operation 13.5 years and 4.9 years earlier than others (P = 0.0020 and P = 0.0026), respectively. We show that the full-length CHADL transcript is expressed in cartilage. Furthermore, the premature stop codon introduced by the CHADL frameshift mutation results in nonsense-mediated decay of the mutant transcripts.

The first international workshop on the epigenetics of osteoarthritis.

Osteoarthritis (OA) is a major clinical problem across the world, in part due to the lack of disease-modifying drugs resulting, to a significant degree, from our incomplete understanding of the underlying molecular mechanisms of the disease. Emerging evidence points to a role of epigenetics in the pathogenesis of OA, but research in this area is still in its early stages. In order to summarize current knowledge and to facilitate the potential coordination of future research activities, the first international workshop on the epigenetics of OA was held in Amsterdam in October 2015. Recent findings on DNA methylation and hydroxymethylation, histone modifications, noncoding RNAs, and other epigenetic mechanisms were presented and discussed. The workshop demonstrated the advantage of bringing together those working in this nascent field and highlights from the event are summarized in this report in the form of summaries from invited speakers and organizers.

Analysis of genetics and DNA methylation in osteoarthritis: What have we learnt about the disease?

Osteoarthritis (OA) is a chronic musculoskeletal disease characterised by the destruction of articular cartilage, synovial inflammation and bone remodelling. Disease aetiology is complex and highly heritable, with genetic variation estimated to contribute to 50% of OA occurrence. Epigenetic alterations, including DNA methylation changes, have also been implicated in OA pathophysiology. This review examines what genetic and DNA methylation studies have taught us about the genes and pathways involved in OA pathology. The influence of DNA methylation on the molecular mechanisms underlying OA genetic risk and the consequence of this interaction on disease susceptibility and penetrance are also discussed.

Functional Characterization of the Osteoarthritis Susceptibility Mapping to CHST11-A Bioinformatics and Molecular Study.

The single nucleotide polymorphism (SNP) rs835487 is associated with hip osteoarthritis (OA) at the genome-wide significance level and is located within CHST11, which codes for carbohydrate sulfotransferase 11. This enzyme post-translationally modifies proteoglycan prior to its deposition in the cartilage extracellular matrix. Using bioinformatics and experimental analyses, our aims were to characterise the rs835487 association signal and to identify the causal functional variant/s. Database searches revealed that rs835487 resides within a linkage disequilibrium (LD) block of only 2.7 kb and is in LD (r2 ≥ 0.8) with six other SNPs. These are all located within intron 2 of CHST11, in a region that has predicted enhancer activity and which shows a high degree of conservation in primates. Luciferase reporter assays revealed that of the seven SNPs, rs835487 and rs835488, which have a pairwise r2 of 0.962, are the top functional candidates; the haplotype composed of the OA-risk conferring G allele of rs835487 and the corresponding T allele of rs835488 (the G-T haplotype) demonstrated significantly different enhancer activity relative to the haplotype composed of the non-risk A allele of rs835487 and the corresponding C allele of rs835488 (the A-C haplotype) (p < 0.001). Electrophoretic mobility shift assays and supershifts identified several transcription factors that bind more strongly to the risk-conferring G and T alleles of the two SNPs, including SP1, SP3, YY1 and SUB1. CHST11 was found to be upregulated in OA versus non-OA cartilage (p < 0.001) and was expressed dynamically during chondrogenesis. Its expression in adult cartilage did not however correlate with rs835487 genotype. Our data demonstrate that the OA susceptibility is mediated by differential protein binding to the alleles of rs835487 and rs835488, which are located within an enhancer whose target may be CHST11 during chondrogenesis or an alternative gene.