Osteoarthritis year 2013 in review: genetics and genomics

Progress in genetic research has delivered important highlights in the last year. One of the widest impact is the publication of the Encyclopedia of DNA Elements (ENCODE) project showing the impressive complexity of the human genome and providing information useful for all areas of genetics. More specific of osteoarthritis (OA) has been the incorporation of DOT1-like, histone H3 methyltransferase (DOT1L) to the list of 11 OA loci with genome-wide significant association, the demonstration of significant overlap between OA genetics and height or body mass index (BMI) genetics, and the tentative prioritization of HMG-box transcription factor 1 (HBP1) in the 7q22 locus based on functional analysis. In addition, the first large scale analysis of DNA methylation has found modest differences between OA and normal cartilage, but has identified a subgroup of OA patients with a very differentiated phenotype. The role of DNA methylation in regulation of NOS2, SOX9, MMP13 and IL1B has been further clarified. MicroRNA expression studies in turn have shown some replication of differences between OA and control cartilage from previous profiling studies and have identified potential regulators of TGF beta signaling and of IL1 beta effects. In addition, non-coding RNAs showed promising results as serum biomarkers of cartilage damage. Gene expression microarray studies have found important differences between studies of hip or knee OA that reinforce the idea of joint specificity in OA. Expression differences between articular cartilage and other types of cartilage highlighted the WNT pathway whose regulation is proposed as critical for maintaining the articular cartilage phenotype. Many of these results need confirmation but they signal the exciting progress that is taking place in all areas of OA genetics, indicate questions requiring more study and augur further interesting discoveries. (C) 2013 Osteoarthritis Research Society International

Genetic distance as an alternative to physical distance for definition of gene units in association studies

Background: Some association studies, as the implemented in VEGAS, ALIGATOR, i-GSEA4GWAS, GSA-SNP and other software tools, use genes as the unit of analysis. These genes include the coding sequence plus flanking sequences. Polymorphisms in the flanking sequences are of interest because they involve cis-regulatory elements or they inform on untyped genetic variants trough linkage disequilibrium. Gene extensions have customarily been defined as ±50 Kb. This approach is not fully satisfactory because genetic relationships between neighbouring sequences are a function of genetic distances, which are only poorly replaced by physical distances. Results: Standardized recombination rates (SRR) from the deCODE recombination map were used as units of genetic distances. We searched for a SRR producing flanking sequences near the ±50 Kb offset that has been common in previous studies. A SRR≥2 was selected because it led to gene extensions with median length=45.3 Kb and the simplicity of an integer value. As expected, boundaries of the genes defined with the ±50 Kb and with the SRR≥2 rules were rarely concordant. The impact of these differences was illustrated with the interpretation of top association signals from two large studies including many hits and their detailed analysis based in different criteria. The definition based in genetic distance was more concordant with the results of these studies than the based in physical distance. In the analysis of 18 top disease associated loci form the first study, the SRR≥2 genes led to a fully concordant interpretation in 17 loci; the ±50 Kb genes only in 6. Interpretation of the 43 putative functional genes of the second study based in the SRR≥2 definition only missed 4 of the genes, whereas the based in the ±50 Kb definition missed 10 genes. Conclusions: A gene definition based on genetic distance led to results more concordant with expert detailed analyses than the commonly used based in physical distance. The genome coordinates for each gene are provided to maintain a simple use of the new definitions

Particular association of clinical and genetic features with autoimmunity to citrullinated α-enolase in rheumatoid arthritis.

OBJECTIVE: To confirm that the presence of anti-citrullinated alpha-enolase peptide 1 (anti-CEP-1) antibodies identifies a subgroup of patients with rheumatoid arthritis (RA). METHODS: DNA and serum samples were obtained from 451 patients with RA and 279 healthy control subjects, all of whom were of Spanish ancestry. Antibodies to cyclic citrullinated peptide (CCP) and CEP-1 were measured by enzyme-linked immunosorbent assay. HLA-DRB1 and the R620W single-nucleotide polymorphism of PTPN22 were genotyped. RESULTS: Anti-CEP-1 and anti-CCP antibodies were observed in 26.8% and 71.2% of the patients with RA, respectively. Most of the patients (86.6%) with anti-CEP-1 antibodies also had anti-CCP antibodies. Erosive arthritis, rheumatoid factor (RF) positivity, and the presence of the HLA shared epitope (especially the DRB1*04 alleles) were disproportionately associated with the group of patients with both antibodies. In addition, evidence of a significant interaction between the shared epitope and the risk allele of PTPN22 was observed only in these patients. In contrast, the association with these clinical and genetic features was weaker in patients with anti-CCP antibodies but lacking anti-CEP-1 antibodies. These results were obtained in patients in whom the prevalence of RA risk factors differed from that in other previously studied patients. CONCLUSION: We observed that autoimmunity against citrullinated alpha-enolase may identify a subset of patients with a higher frequency of joint erosions and RF positivity. In addition, we confirmed the disproportionately large effect of the susceptibility alleles of HLA-DRB1 and their interaction with PTPN22 in this subset of patients. These results extend, confirm, and generalize the evidence supporting the specificity of the anti-CEP-1 antibody-positive subgroup of patients with RA among anti-CCP antibody-positive patients with RA

Ultrasonographic assessment of enthesitis in HLA-B27 positive patients with rheumatoid arthritis, a matched case-only study

Introduction HLA-B27 has a modifier effect on the phenotype of multiple diseases, both associated and non-associated with it. Among these effects, an increased frequency of clinical enthesitis in patients with Rheumatoid Arthritis (RA) has been reported but never explored again. We aimed to replicate this study with a sensitive and quantitative assessment of enthesitis by using standardized ultrasonography (US). Methods The Madrid Sonography Enthesitis Index (MASEI) was applied to the US assessment of 41 HLA-B27 positive and 41 matched HLA-B27 negative patients with longstanding RA. Clinical characteristics including explorations aimed to evaluate spondyloarthrtitis and laboratory tests were also done. Results A significant degree of abnormalities in the entheses of the patients with RA were found, but the MASEI values, and each of its components including the Doppler signal, were similar in HLA-B27 positive and negative patients. An increase of the MASEI scores with age was identified. Differences in two clinical features were found: a lower prevalence of rheumatoid factor and a more common story of low back pain in the HLA-B27 positive patients than in the negative. The latter was accompanied by radiographic sacroiliitis in two HLA-B27 positive patients. No other differences were detected. Conclusion We have found that HLA-B27 positive patients with RA do not have more enthesitis as assessed with US than the patients lacking this HLA allele. However, HLA-B27 could be shaping the RA phenotype towards RF seronegativity and axial involvement.The study was supported by grants 10CSA918040PR from the Xunta de Galicia (http://www.sergas.e/MostrarContidos_N3_T01.aspx?IdPaxina=10142) and PI08/0744 of the Instituto de Salud Carlos III (http://www.isciii.es/) that are partially financed by the European Regional Development Fund of the European UnionS

Assessment of global DNA methylation in peripheral blood cell subpopulations of early rheumatoid arthritis before and after methotrexate

INTRODUCTION: DNA methylation is an epigenetic mechanism regulating gene expression that has been insufficiently studied in the blood of rheumatoid arthritis (RA) patients, as only T cells and total peripheral blood mononuclear cells (PBMCs) from patients with established RA have been studied and with conflicting results. METHOD: Five major blood cell subpopulations: T, B and NK cells, monocytes, and polymorphonuclear leukocytes, were isolated from 19 early RA patients and 17 healthy controls. Patient samples were taken before and 1 month after the start of treatment with methotrexate (MTX). Analysis included DNA methylation with high-performance liquid chromatography-electrospray ionization-tandem mass spectrometry-selected reaction monitoring (HPLC-ESI-MS/MS-SRM) and expression levels of seven methylation-specific enzymes by quantitative polymerase chain reaction (qPCR). RESULTS: Disease-modifying anti-rheumatic drug (DMARD)-naive early RA patients showed global DNA hypomethylation in T cells and monocytes, together with a lower expression of DNA methyltrasnferase 1 (DNMT1), the maintenance DNA methyltransferase, which was also decreased in B cells. Furthermore, significantly increased expression of ten-eleven translocation1 (TET1), TET2 and TET3, enzymes involved in demethylation, was found in monocytes and of TET2 in T cells. There was also modest decreased expression of DNMT3A in B cells and of growth arrest and DNA-damage-inducible protein 45A (GADD45A) in T and B cells. Treatment with MTX reverted hypomethylation in T cells and monocytes, which were no longer different from controls, and increased global methylation in B cells. In addition, DNMT1 and DNMT3A showed a trend to reversion of their decreased expression. CONCLUSIONS: Our results confirm global DNA hypomethylation in patients with RA with specificity for some blood cell subpopulations and their reversal with methotrexate treatment. These changes are accompanied by parallel changes in the levels of enzymes involved in methylation, suggesting the possibility of regulation at this level

Non-HLA genes PTPN22, CDK6 and PADI4 are associated with specific autoantibodies in HLA-defined subgroups of rheumatoid arthritis

Introduction: Genetic susceptibility to complex diseases has been intensively studied during the last decade, yet only signals with small effect have been found leaving open the possibility that subgroups within complex traits show stronger association signals. In rheumatoid arthritis (RA), autoantibody production serves as a helpful discriminator in genetic studies and today anti-citrullinated cyclic peptide (anti-CCP) antibody positivity is employed for diagnosis of disease. The HLA-DRB1 locus is known as the most important genetic contributor for the risk of RA, but is not sufficient to drive autoimmunity and additional genetic and environmental factors are involved. Hence, we addressed the association of previously discovered RA loci with disease-specific autoantibody responses in RA patients stratified by HLA-DRB1*04. Methods: We investigated 2178 patients from three RA cohorts from Sweden and Spain for 41 genetic variants and four autoantibodies, including the generic anti-CCP as well as specific responses towards citrullinated peptides from vimentin, alpha-enolase and type II collagen. Results: Our data demonstrated different genetic associations of autoantibody-positive disease subgroups in relation to the presence of DRB1*04. Two specific subgroups of autoantibody-positive RA were identified. The SNP in PTPN22 was associated with presence of anti-citrullinated enolase peptide antibodies in carriers of HLA-DRB1*04 (Cochran-Mantel-Haenszel test P = 0.0001, P corrected <0.05), whereas SNPs in CDK6 and PADI4 were associated with anti-CCP status in DRB1*04 negative patients (Cochran-Mantel-Haenszel test P = 0.0004, P corrected <0.05 for both markers). Additionally we see allelic correlation with autoantibody titers for PTPN22 SNP rs2476601 and anti-citrullinated enolase peptide antibodies in carriers of HLA-DRB1*04 (Mann Whitney test P = 0.02) and between CDK6 SNP rs42041 and anti-CCP in non-carriers of HLA-DRB1*04 (Mann Whitney test P = 0.02). Conclusion: These data point to alternative pathways for disease development in clinically similar RA subgroups and suggest an approach for study of genetic complexity of disease with strong contribution of HLA

Replication of recently identified systemic lupus erythematosus genetic associations: a case–control study

Introduction We aimed to replicate association of newly identified systemic lupus erythematosus (SLE) loci. Methods We selected the most associated SNP in 10 SLE loci. These 10 SNPs were analysed in 1,579 patients with SLE and 1,726 controls of European origin by single-base extension. Comparison of allele frequencies between cases and controls was done with the Mantel–Haenszel approach to account for heterogeneity between sample collections. Results A previously controversial association with a SNP in the TYK2 gene was replicated (odds ratio (OR) = 0.79, P = 2.5 × 10-5), as well as association with the X chromosome MECP2 gene (OR = 1.26, P = 0.00085 in women), which had only been reported in a single study, and association with four other loci, 1q25.1 (OR = 0.81, P = 0.0001), PXK (OR = 1.19, P = 0.0038), BANK1 (OR = 0.83, P = 0.006) and KIAA1542 (OR = 0.84, P = 0.001), which have been identified in a genome-wide association study, but not found in any other study. All these replications showed the same disease-associated allele as originally reported. No association was found with the LY9 SNP, which had been reported in a single study. Conclusions Our results confirm nine SLE loci. For six of them, TYK2, MECP2, 1q25.1, PXK, BANK1 and KIAA1542, this replication is important. The other three loci, ITGAM, STAT4 and C8orf13-BLK, were already clearly confirmed. Our results also suggest that MECP2 association has no influence in the sex bias of SLE, contrary to what has been proposed. In addition, none of the other associations seems important in this respectThe present work was supported by Fondo de Investigacion Sanitaria of the Instituto de Salud Carlos III (Spain), grants 04/1651 and 06/0620 that are partially financed by the Fondo Europeo de Desarrollo Regional program of the European Union, by grants from the Xunta de Galicia, and by BMBF KN Rheuma grant C2.12 (to TW)S

Lack of replication of higher genetic risk load in men than in women with systemic lupus erythematosus

Introduction: We aimed to replicate a recent study which showed higher genetic risk load at 15 loci in men than in women with systemic lupus erythematosus (SLE). This difference was very significant, and it was interpreted as indicating that men require more genetic susceptibility than women to develop SLE. Methods: Nineteen SLE-associated loci (thirteen of which are shared with the previous study) were analyzed in 1,457 SLE patients and 1,728 healthy controls of European ancestry. Genetic risk load was calculated as sex-specific sum genetic risk scores (GRS(s)). Results: Our results did not replicate those of the previous study at either the level of individual loci or the global level of GRS(s). GRS(s) were larger in women than in men (4.20 ± 1.07 in women vs. 3.27 ± 0.98 in men). This very significant difference (P < 10(-16)) was more dependent on the six new loci not included in the previous study (59% of the difference) than on the thirteen loci that are shared (the remaining 41%). However, the 13 shared loci also showed a higher genetic risk load in women than in men in our study (P = 6.6 × 10(-7)), suggesting that heterogeneity of participants, in addition to different loci, contributed to the opposite results. Conclusion: Our results show the lack of a clear trend toward higher genetic risk in one of the sexes for the analyzed SLE loci. They also highlight several limitations of assessments of genetic risk load, including the possibility of ascertainment bias with loci discovered in studies that have included mainly women