DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation.

BACKGROUND: Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions. RESULTS: The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression. CONCLUSIONS: In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.

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.

Regulation of microRNA-221, -222, -21 and -27 in articular cartilage subjected to abnormal compressive forces.

KEY POINTS: microRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression. miRs are involved in regulating cellular activities in response to mechanical loading in all physiological systems, although it is largely unknown whether this response differs with increasing magnitudes of load. miR-221, miR-222, miR-21-5p and miR-27a-5p were significantly increased in ex vivo cartilage explants subjected to increasing load magnitude and in in vivo joint cartilage exposed to abnormal loading. TIMP3 and CPEB3 are putative miR targets in chondrocytes Identification of mechanically regulated miRs that have potential to impact on tissue homeostasis provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology, in all physiological systems. ABSTRACT: MicroRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression and are involved in mechano-regulation of cellular activities in all physiological systems. It is unknown whether such epigenetic mechanisms are regulated in response to increasing magnitudes of load. The present study investigated mechano-regulation of miRs in articular cartilage subjected to 'physiological' and 'non-physiological' compressive loads in vitro as a model system and validated findings in an in vivo model of abnormal joint loading. Bovine full-depth articular cartilage explants were loaded to 2.5 MPa (physiological) or 7 MPa (non-physiological) (1 Hz, 15 min) and mechanically-regulated miRs identified using next generation sequencing and verified using a quantitative PCR. Downstream targets were verified using miR-specific mimics or inhibitors in conjunction with 3'-UTR luciferase activity assays. A subset of miRs were mechanically-regulated in ex vivo cartilage explants and in vivo joint cartilage. miR-221, miR-222, miR-21-5p and miR-27a-5p were increased and miR-483 levels decreased with increasing load magnitude. Tissue inhibitor of metalloproteinase 3 (TIMP3) and cytoplasmic polyadenylation element binding protein 3 (CPEB3) were identified as putative downstream targets. Our data confirm miR-221 and -222 mechano-regulation and demonstrates novel mechano-regulation of miR-21-5p and miR-27a-5p in ex vivo and in vivo cartilage loading models. TIMP3 and CPEB3 are putative miR targets in chondrocytes. Identification of specific miRs that are regulated by increasing load magnitude, as well as their potential to impact on tissue homeostasis, has direct relevance to other mechano-sensitive physiological systems and provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology.

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.

The long non-coding RNA ROCR contributes to SOX9 expression and chondrogenic differentiation of human mesenchymal stem cells.

Long non-coding RNAs (lncRNAs) are expressed in a highly tissue-specific manner and function in various aspects of cell biology, often as key regulators of gene expression. In this study, we established a role for lncRNAs in chondrocyte differentiation. Using RNA sequencing we identified a human articular chondrocyte repertoire of lncRNAs from normal hip cartilage donated by neck of femur fracture patients. Of particular interest are lncRNAs upstream of the master chondrocyte transcription factor SOX9 locus. SOX9 is an HMG-box transcription factor that plays an essential role in chondrocyte development by directing the expression of chondrocyte-specific genes. Two of these lncRNAs are upregulated during chondrogenic differentiation of mesenchymal stem cells (MSCs). Depletion of one of these lncRNAs, LOC102723505, which we termed ROCR (regulator of chondrogenesis RNA), by RNA interference disrupted MSC chondrogenesis, concomitant with reduced cartilage-specific gene expression and incomplete matrix component production, indicating an important role in chondrocyte biology. Specifically, SOX9 induction was significantly ablated in the absence of ROCR, and overexpression of SOX9 rescued the differentiation of MSCs into chondrocytes. Our work sheds further light on chondrocyte-specific SOX9 expression and highlights a novel method of chondrocyte gene regulation involving a lncRNA.

Predicting drug-free remission in rheumatoid arthritis: A prospective interventional cohort study.

BACKGROUND: Many patients with rheumatoid arthritis (RA) achieve disease remission with modern treatment strategies. However, having achieved this state, there are no tests that predict when withdrawal of therapy will result in drug-free remission rather than flare. We aimed to identify predictors of drug-free remission in RA. METHODS: The Biomarkers of Remission in Rheumatoid Arthritis (BioRRA) Study was a unique, prospective, interventional cohort study of complete and abrupt cessation of conventional synthetic disease-modifying anti-rheumatic drugs (DMARDs). Patients with RA of at least 12 months duration and in clinical and ultrasound remission discontinued DMARDs and were monitored for six months. The primary outcome was time-to-flare, defined as disease activity score in 28 joints with C-reactive protein (DAS28-CRP) ≥ 2.4. Baseline clinical and ultrasound measures, circulating inflammatory biomarkers, and peripheral CD4+ T cell gene expression were assessed for their ability to predict time-to-flare and flare/remission status by Cox regression and receiver-operating characteristic (ROC) analysis respectively. RESULTS: 23/44 (52%) eligible patients experienced an arthritis flare after a median (IQR) of 48 (31.5-86.5) days following DMARD cessation. A composite score incorporating five baseline variables (three transcripts [FAM102B, ENSG00000228010, ENSG00000227070], one cytokine [interleukin-27], one clinical [Boolean remission]) differentiated future flare from drug-free remission with an area under the ROC curve of 0.96 (95% CI 0.91-1.00), sensitivity 0.91 (0.78-1.00) and specificity 0.95 (0.84-1.00). CONCLUSION: We provide proof-of-concept evidence for predictors of drug-free remission in RA. If validated, these biomarkers could help to personalize immunosuppressant withdrawal: a therapy paradigm shift with ensuing patient and economic benefits.

Expression of STAT3-regulated genes in circulating CD4+ T cells discriminates rheumatoid arthritis independently of clinical parameters in early arthritis.

OBJECTIVES: Dysregulated signal transduction and activator of transcription-3 (STAT3) signalling in CD4+ T cells has been proposed as an early pathophysiological event in RA. We sought further evidence for this observation, and to determine its clinical relevance. METHODS: Microarray technology was used to measure gene expression in purified peripheral blood CD4+ T cells from treatment-naïve RA patients and disease controls newly recruited from an early arthritis clinic. Analysis focused on 12 previously proposed transcripts, and concurrent STAT3 pathway activation was determined in the same cells by flow cytometry. A pooled analysis of previous and current gene expression findings incorporated detailed clinical parameters and employed multivariate analysis. RESULTS: In an independent cohort of 161 patients, expression of 11 of 12 proposed signature genes differed significantly between RA patients and controls, robustly validating the earlier findings. Differential regulation was most pronounced for the STAT3 target genes PIM1, BCL3 and SOCS3 (>1.3-fold difference; P < 0.005), each of whose expression correlated strongly with paired intracellular phospho-STAT3. In a meta-analysis of 279 patients the same three genes accounted for the majority of the signature's ability to discriminate RA patients, which was found to be independent of age, joint involvement or acute phase response. CONCLUSION: The STAT3-mediated dysregulation of BCL3, SOCS3 and PIM1 in circulating CD4+ T cells is a discriminatory feature of early RA that occurs independently of acute phase response. The mechanistic and functional implications of this observation at a cellular level warrant clarification.

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.

Expression analysis of the osteoarthritis genetic susceptibility locus mapping to an intron of the MCF2L gene and marked by the polymorphism rs11842874.

BACKGROUND: Osteoarthritis (OA) is a painful, debilitating disease characterised by loss of articular cartilage with concurrent changes in other tissues of the synovial joint. Genetic association studies have shown that a number of common variants increase the risk of developing OA. Investigating their activity can uncover novel causal pathways and potentially highlight new treatment targets. One of the reported OA association signals is marked by the single nucleotide polymorphism (SNP) rs11842874 at chromosome 13q34. rs11842874 is positioned within a small linkage disequilibrium (LD) block within intron 4 of MCF2L, a gene encoding guanine-nucleotide exchange factor DBS. There are no non-synonymous SNPs that correlate with this association signal and we therefore set out to assess whether its effect on OA susceptibility is mediated by alteration of MCF2L expression. METHODS: Nucleic acid was extracted from cartilage, synovial membrane or infrapatellar fat pad tissues from OA patients. Expression of MCF2L was measured by quantitative PCR and RNA-sequencing whilst the presence of DBS was studied using immunohistochemistry. The functional effect of SNPs within the 13q34 locus was assessed using public databases and in vitro using luciferase reporter analysis. RESULTS: MCF2L gene and protein expression are detectable in joint tissues, with quantitative differences in the expression of the gene and in the transcript isoforms expressed between the tissues tested. There is an expression quantitative trait locus (eQTL) operating within synovial membrane tissue, with possession of the risk-conferring A allele of rs11842874 correlating with increased MCF2L expression. SNPs within the rs11842874 LD block reside within transcriptional regulatory elements and their direct analysis reveals that several show quantitative differences in regulatory activity at the allelic level. CONCLUSIONS: MCF2L is subject to a cis-acting eQTL in synovial membrane that correlates with the OA association signal. This signal contains several functional SNPs that could account for the susceptibility and which therefore merit further investigation. As far as we are aware, this is the first example of an OA susceptibility locus operating as an eQTL in synovial membrane tissue but not in cartilage.

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.

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.

IL-6 Mediated Transcriptional Programming of Naïve CD4+ T Cells in Early Rheumatoid Arthritis Drives Dysregulated Effector Function.

Objective: We have previously shown that increased circulating interleukin-6 (IL-6) results in enhanced CD4+ T cell signaling via signal transduction and activator of transcription-3 (STAT3) in early rheumatoid arthritis (RA). We tested the hypothesis that transcriptional "imprinting" of T-cells by this mechanism skews downstream effector responses, reinforcing immune dysregulation at a critical, but targetable, disease phase. Methods: We modeled naïve CD4+ T cell exposure to pathophysiological concentrations of IL-6 in vitro, assessing the dynamic transcriptional and functional consequences for downstream effector cells utilizing microarray and flow cytometry. Fresh blood from treatment-naïve early arthritis patients was phenotyped in parallel for comparison. Results: T cell sensitivity to IL-6 was most marked in the naïve subset, and related to gp130 rather than IL-6R expression. Exposure of healthy naïve CD4+ T cells to IL-6 induced the same STAT3 target genes as previously seen to discriminate RA patients from disease controls. After TCR stimulation IL-6 pre-exposed cells exhibited enhanced proliferative capacity, activation, and a propensity toward Th1 differentiation, compared to non-exposed cells. An entirely analogous phenotype was observed in early RA compared to control CD4+ T cells. Conclusions: Sustained IL-6 exposure at a critical point in the natural history of RA "primes" the adaptive immune system to respond aberrantly to TCR stimulation, potentiating disease induction with implications for the optimal timing of targeted therapy.

Severe type I interferonopathy and unrestrained interferon signaling due to a homozygous germline mutation in STAT2.

Excessive type I interferon (IFNα/ß) activity is implicated in a spectrum of human disease, yet its direct role remains to be conclusively proven. We investigated two siblings with severe early-onset autoinflammatory disease and an elevated IFN signature. Whole-exome sequencing revealed a shared homozygous missense Arg148Trp variant in STAT2, a transcription factor that functions exclusively downstream of innate IFNs. Cells bearing STAT2R148W in homozygosity (but not heterozygosity) were hypersensitive to IFNα/ß, which manifest as prolonged Janus kinase-signal transducers and activators of transcription (STAT) signaling and transcriptional activation. We show that this gain of IFN activity results from the failure of mutant STAT2R148W to interact with ubiquitin-specific protease 18, a key STAT2-dependent negative regulator of IFNα/ß signaling. These observations reveal an essential in vivo function of STAT2 in the regulation of human IFNα/ß signaling, providing concrete evidence of the serious pathological consequences of unrestrained IFNα/ß activity and supporting efforts to target this pathway therapeutically in IFN-associated disease.

Symptom-based stratification of patients with primary Sjögren's syndrome: multi-dimensional characterisation of international observational cohorts and reanalyses of randomised clinical trials.

BACKGROUND: Heterogeneity is a major obstacle to developing effective treatments for patients with primary Sjögren's syndrome. We aimed to develop a robust method for stratification, exploiting heterogeneity in patient-reported symptoms, and to relate these differences to pathobiology and therapeutic response. METHODS: We did hierarchical cluster analysis using five common symptoms associated with primary Sjögren's syndrome (pain, fatigue, dryness, anxiety, and depression), followed by multinomial logistic regression to identify subgroups in the UK Primary Sjögren's Syndrome Registry (UKPSSR). We assessed clinical and biological differences between these subgroups, including transcriptional differences in peripheral blood. Patients from two independent validation cohorts in Norway and France were used to confirm patient stratification. Data from two phase 3 clinical trials were similarly stratified to assess the differences between subgroups in treatment response to hydroxychloroquine and rituximab. FINDINGS: In the UKPSSR cohort (n=608), we identified four subgroups: Low symptom burden (LSB), high symptom burden (HSB), dryness dominant with fatigue (DDF), and pain dominant with fatigue (PDF). Significant differences in peripheral blood lymphocyte counts, anti-SSA and anti-SSB antibody positivity, as well as serum IgG, κ-free light chain, ß2-microglobulin, and CXCL13 concentrations were observed between these subgroups, along with differentially expressed transcriptomic modules in peripheral blood. Similar findings were observed in the independent validation cohorts (n=396). Reanalysis of trial data stratifying patients into these subgroups suggested a treatment effect with hydroxychloroquine in the HSB subgroup and with rituximab in the DDF subgroup compared with placebo. INTERPRETATION: Stratification on the basis of patient-reported symptoms of patients with primary Sjögren's syndrome revealed distinct pathobiological endotypes with distinct responses to immunomodulatory treatments. Our data have important implications for clinical management, trial design, and therapeutic development. Similar stratification approaches might be useful for patients with other chronic immune-mediated diseases. FUNDING: UK Medical Research Council, British Sjogren's Syndrome Association, French Ministry of Health, Arthritis Research UK, Foundation for Research in Rheumatology. VIDEO ABSTRACT.

DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation

BACKGROUND: Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions. RESULTS: The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression. CONCLUSIONS: In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.

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.