Luis Rodriguez-Rodriguez

Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2, 3, 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses—as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes—to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis

Using the Immunochip custom SNP array, which was designed for dense genotyping of 186 loci identified through genome-wide association studies (GWAS), we analyzed 11,475 individuals with rheumatoid arthritis (cases) of European ancestry and 15,870 controls for 129,464 markers. We combined these data in a meta-analysis with GWAS data from additional independent cases (n = 2,363) and controls (n = 17,872). We identified 14 new susceptibility loci, 9 of which were associated with rheumatoid arthritis overall and five of which were specifically associated with disease that was positive for anticitrullinated peptide antibodies, bringing the number of confirmed rheumatoid arthritis risk loci in individuals of European ancestry to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at 6 loci and identified association to low-frequency variants at 4 loci. Bioinformatic analyses generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.

A genetic study on C5-TRAF1 and progression of joint damage in rheumatoid arthritis

IntroductionThe severity of joint damage progression in rheumatoid arthritis (RA) is heritable. Several genetic variants have been identified, but together explain only part of the total genetic effect. Variants in Interleukin-6 (IL-6), Interleukin-10 (IL-10), C5-TRAF1, and Fc-receptor-like-3 (FCRL3) have been described to associate with radiographic progression, but results of different studies were incongruent. We aimed to clarify associations of these variants with radiographic progression by evaluating six independent cohorts.MethodsIn total 5,895 sets of radiographs of 2,493 RA-patients included in six different independent datasets from the Netherlands, Sweden, Spain and North-America were studied in relation to rs1800795 (IL-6), rs1800896 (IL-10), rs2900180 (C5-TRAF1) and rs7528684 (FCRL3). Associations were tested in the total RA-populations and in anti-citrullinated peptide antibodies (ACPA)-positive and ACPA-negative subgroups per cohort, followed by meta-analyses. Furthermore, the associated region C5-TRAF1 was fine-mapped in the ACPA-negative Dutch RA-patients.ResultsNo associations were found for rs1800795 (IL-6), rs1800896 (IL-10) and rs7528684 (FCRL3) in the total RA-population and after stratification for ACPA. Rs2900180 in C5-TRAF1 was associated with radiographic progression in the ACPA-negative population (P-value meta-analysis = 5.85 × 10−7); the minor allele was associated with more radiographic progression. Fine-mapping revealed a region of 66Kb that was associated; the lowest P-value was for rs7021880 in TRAF1. The P-value for rs7021880 in meta-analysis was 6.35 × 10−8. Previous studies indicate that the region of rs7021880 was associated with RNA expression of TRAF1 and C5.ConclusionVariants in IL-6, IL-10 and FCRL3 were not associated with radiographic progression. Rs2900180 in C5-TRAF1 and linked variants in a 66Kb region were associated with radiographic progression in ACPA-negative RA.

Identification of CSK as a systemic sclerosis genetic risk factor through Genome Wide Association Study follow-up

Systemic sclerosis (SSc) is complex autoimmune disease affecting the connective tissue; influenced by genetic and environmental components. Recently, we performed the first successful genome-wide association study (GWAS) of SSc. Here, we perform a large replication study to better dissect the genetic component of SSc. We selected 768 polymorphisms from the previous GWAS and genotyped them in seven replication cohorts from Europe. Overall significance was calculated for replicated significant SNPs by meta-analysis of the replication cohorts and replication-GWAS cohorts (3237 cases and 6097 controls). Six SNPs in regions not previously associated with SSc were selected for validation in another five independent cohorts, up to a total of 5270 SSc patients and 8326 controls. We found evidence for replication and overall genome-wide significance for one novel SSc genetic risk locus: CSK [P-value = 5.04 × 10(-12), odds ratio (OR) = 1.20]. Additionally, we found suggestive association in the loci PSD3 (P-value = 3.18 × 10(-7), OR = 1.36) and NFKB1 (P-value = 1.03 × 10(-6), OR = 1.14). Additionally, we strengthened the evidence for previously confirmed associations. This study significantly increases the number of known putative genetic risk factors for SSc, including the genes CSK, PSD3 and NFKB1, and further confirms six previously described ones.

High-density genotyping of immune loci in Koreans and Europeans identifies eight new rheumatoid arthritis risk loci.

OBJECTIVE A highly polygenic aetiology and high degree of allele-sharing between ancestries have been well elucidated in genetic studies of rheumatoid arthritis. Recently, the high-density genotyping array Immunochip for immune disease loci identified 14 new rheumatoid arthritis risk loci among individuals of European ancestry. Here, we aimed to identify new rheumatoid arthritis risk loci using Korean-specific Immunochip data. METHODS We analysed Korean rheumatoid arthritis case-control samples using the Immunochip and genome-wide association studies (GWAS) array to search for new risk alleles of rheumatoid arthritis with anticitrullinated peptide antibodies. To increase power, we performed a meta-analysis of Korean data with previously published European Immunochip and GWAS data for a total sample size of 9299 Korean and 45,790 European case-control samples. RESULTS We identified eight new rheumatoid arthritis susceptibility loci (TNFSF4, LBH, EOMES, ETS1-FLI1, COG6, RAD51B, UBASH3A and SYNGR1) that passed a genome-wide significance threshold (p<5×10(-8)), with evidence for three independent risk alleles at 1q25/TNFSF4. The risk alleles from the seven new loci except for the TNFSF4 locus (monomorphic in Koreans), together with risk alleles from previously established RA risk loci, exhibited a high correlation of effect sizes between ancestries. Further, we refined the number of single nucleotide polymorphisms (SNPs) that represent potentially causal variants through a trans-ethnic comparison of densely genotyped SNPs. CONCLUSIONS This study demonstrates the advantage of dense-mapping and trans-ancestral analysis for identification of potentially causal SNPs. In addition, our findings support the importance of T cells in the pathogenesis and the fact of frequent overlap of risk loci among diverse autoimmune diseases.

A genome-wide association study of rheumatoid arthritis without antibodies against citrullinated peptides

Introduction Rheumatoid arthritis (RA) patients can be classified based on presence or absence of anticitrullinated peptide antibodies (ACPA) in their serum. This heterogeneity among patients may reflect important biological differences underlying the disease process. To date, the majority of genetic studies have focused on the ACPA-positive group. Therefore, our goal was to analyse the genetic risk factors that contribute to ACPA-negative RA. Methods We performed a large-scale genome-wide association study (GWAS) in three Caucasian European cohorts comprising 1148 ACPA-negative RA patients and 6008 controls. All patients were screened using the Illumina Human Cyto-12 chip, and controls were genotyped using different genome-wide platforms. Population-independent analyses were carried out by means of logistic regression. Meta-analysis with previously published data was performed as follow-up for selected signals (reaching a total of 1922 ACPA-negative RA patients and 7087 controls). Imputation of classical HLA alleles, amino acid residues and single nucleotide polymorphisms was undertaken. Results The combined analysis of the studied cohorts resulted in identification of a peak of association in the HLA-region and several suggestive non-HLA associations. Meta-analysis with previous reports confirmed the association of the HLA region with this subset and an observed association in the CLYBL locus remained suggestive. The imputation and deep interrogation of the HLA region led to identification of a two amino acid model (HLA-B at position 9 and HLA-DRB1 at position 11) that accounted for the observed genome-wide associations in this region. Conclusions Our study shed light on the influence of the HLA region in ACPA-negative RA and identified a suggestive risk locus for this condition.

CCR5Δ32 variant and cardiovascular disease in patients with rheumatoid arthritis: a cohort study

IntroductionThe aim of our study was to analyze the influence of the CCR5Δ32 polymorphism in the risk of cardiovascular (CV) events and subclinical atherosclerosis among patients with rheumatoid arthritis (RA).MethodsA total of 645 patients fulfilling the American Rheumatism Association 1987 revised classification criteria for RA were studied. Patients were genotyped for the CCR5 rs333 polymorphism using predesigned TaqMan assays. Also, HLA DRB1 genotyping was performed using molecular-based methods. Carotid intima-media thickness, flow-mediated endothelium-dependent dilatation (FMD) and endothelium-independent vasodilatation, which were used as surrogate markers of subclinical atherosclerosis, were measured in a subgroup of patients with no clinical CV disease.ResultsA lower frequency of carriers of the CCR5Δ32 allele among patients with CV events (3.4% versus 11.3%, P = 0.025, odds ratio 0.28, 95% confidence interval (95% CI) 0.06 to 0.89) was observed. However, after adjusting for gender, age at time of RA diagnosis, and the presence of shared epitope, rheumatoid factor and classic CV risk factors in the Cox regression analysis, this reduction of CV events in CCR5Δ32 allele carriers was slightly outside the range of significance (P = 0.097; hazard ratio 0.37 (95% CI 0.12 to 1.19)). Carriers of the CCR5Δ32 deletion also showed higher FMD values than the remaining patients (CCR5/CCR5Δ32 patients: 7.03% ± 6.61% versus CCR5/CCR5 patients: 5.51% ± 4.66%). This difference was statistically significant when analysis of covariance was performed (P = 0.024).ConclusionsOur results show a potential influence of the CCR5Δ32 deletion on the risk of CV disease among patients with RA. This may be due to a protective effect of this allelic variant against the development of vascular endothelial dysfunction.

Cardiovascular risk assessment in patients with rheumatoid arthritis: The relevance of clinical, genetic and serological markers

Cardiovascular disease (CV) is the most common cause of premature mortality in patients with rheumatoid arthritis (RA). This is the result of an accelerated atherosclerotic process. Adequate CV risk stratification has special relevance in RA to identify patients at risk of CV disease. However, current CV risk screening and management strategies underestimate the actual CV risk in RA. Consequently, the search for additional tools that may help to identify those patients at high CV risk has become a key objective in the last years. In this regard, non-invasive surrogates, such as carotid ultrasonography, have been found to be excellent predictors of future CV events. In addition, several studies have revealed the relevance of a genetic component in the development of CV disease in RA patients. Besides an association with HLA-DRB1* shared epitope alleles other gene polymorphisms located inside and outside the HLA seem to influence the risk of cardiovascular disease in RA. Moreover, serum levels of some metabolic syndrome-related biomarkers, adipokines such as adiponectin and biomarkers of endothelial cell activation and inflammation such as Osteoprotegerin and Asymmetric dimethylarginine have recently been found useful for the prediction of CV disease in these patients. An update of the current knowledge on these potential markers, especially focused on new genetic and serological biomarkers is shown in this review.

Genetic Markers of Cardiovascular Disease in Rheumatoid Arthritis

Cardiovascular (CV) disease is the most common cause of premature mortality in patients with rheumatoid arthritis (RA). It is the result of an accelerated atherosclerotic process. Both RA and atherosclerosis are complex polygenic diseases. Besides traditional CV risk factors and chronic inflammation, a number of studies have confirmed the role of genetic factors in the development of the atherogenesis observed in RA. In this regard, besides a strong association between the HLA-DRB1∗04 shared epitope alleles and both endothelial dysfunction, an early step in the atherosclerotic process, and clinically evident CV disease, other polymorphisms belonging to genes implicated in inflammatory and metabolic pathways, located inside and outside the HLA region, such as the 308 variant (G > A, rs1800629) of the TNFA locus, the rs1801131 polymorphism (A > C; position + 1298) of the MTHFR locus, or a deletion of 32 base pairs on the CCR5 gene, seem to be associated with the risk of CV disease in patients with RA. Despite considerable effort to decipher the genetic basis of CV disease in RA, further studies are required to better establish the genetic influence in the increased risk of CV events observed in patients with RA.

Large-scale gene expression in bone marrow mesenchymal stem cells: a putative role for COL10A1 in osteoarthritis

Objectives To elucidate disease-specific molecular changes occurring in osteoarthritis (OA) by analysing the differential gene expression profiles of bone marrow mesenchymal stem cells (BM-MSCs) from patients with OA compared with those without OA. Methods Expression profiles of BM-MSCs from eight paired patients with OA and patients with hip fracture without signs of OA were compared by DNA microarray expression analysis and significant differences were evaluated by computational Gene Set Enrichment Analysis. To validate the involvement of COL10A1 as part of the most downregulated gene set in OA, three tagging single nucleotide polymorphisms were genotyped in 191 patients with OA and 283 controls. COL10A1 expression was also assessed by quantitative RT-PCR in additional subjects. Results Expression levels in 9% (1967) of the overall transcripts were significantly different (p<0.05) between MSCs from patients with OA and controls (532 genes reached twofold differences: 240 were upregulated and 292 were downregulated). Cell development and differentiation were the functional categories accounting for most genes with altered expression. Interestingly, several genes related to the Wnt/-catenin pathway and collagen genes were downregulated in MSCs from patients with OA. The collagen gene set was clearly downregulated in OA. Furthermore, the expression of COL10A1 was significantly reduced in patients with OA. A genetic association between the COL10A1 rs11965969 polymorphism and OA was also found. Conclusion COL10A1 downregulation seems to have a role in the establishment of a defective and/or unstable subchondral cartilage matrix in OA disease. It is proposed that OA may be linked to the intrinsic defective regenerative potential of BM-MSCs resulting from its reduced expression of fate commitment-related genes.