Toll-like receptor 3 increases antigen-presenting cell responses to a pro-apoptotic stimulus, yet does not contribute to systemic lupus erythematosus genetic susceptibility.

OBJECTIVES: TLR3 mediates skin solar injury by binding nuclear material released from apoptotic keratinocytes, resulting in the production of pro-inflammatory cytokines. Because the TLR3 gene is located in 4q35, a known systemic lupus erythematosus (SLE) susceptibility locus, we wondered whether TLR3 single nucleotide polymorphisms (SNPs) were associated with inflammatory mechanisms relevant to the development of SLE, and disease susceptibility. METHODS: Functional assays were carried out in TLR3-transfected HEK293 cells and in monocyte-derived dendritic cells (moDCs). TLR3 and IFNß immunofluorescence studies were performed in skin samples from 7 SLE patients and 3 controls. We performed a SNP association study in a discovery cohort of 153 patients and 105 controls, followed by a confirmation study in an independent cohort of 1,380 patients and 2,104 controls. RESULTS: TLR3 and IFNß are overexpressed in SLE skin lesions. TLR3 overexpression in HEK293 cells amplifies their sensitivity to a pro-apoptotic stimulus. Taking advantage of a naturally occurring polymorphic TLR3 variant (rs3775291) that weakly versus strongly responds to poly I:C stimulation, we found that TLR3 is associated with amplified apoptotic responses, production of the Ro/SSA autoantigen and increased maturation of myeloid-derived dendritic cells (moDC) after exposure to UV irradiation. However, TLR3 SNPs are not associated with susceptibility to SLE in a large population of patients and controls. CONCLUSIONS: TLR3 is overexpressed in SLE skin lesions and amplifies apoptotic and inflammatory responses to UV-irradiation in antigen-presenting cells in vitro. However, TLR3 SNPs do not impact susceptibility to the development of the disease.

Evaluation of a clinical pharmacogenetics model to predict methotrexate response in patients with rheumatoid arthritis.

Variability of response to treatment hinders successful management of rheumatoid arthritis (RA). Consequently, a clinical pharmacogenetics model for predicting response to methotrexate (CP-MTX) has been previously proposed that includes four clinical variables (disease activity, sex, the presence of rheumatoid factor and smoking status) and four SNPs (rs2236225, rs17602729, rs1127354, and rs2372536) in genes of the folate pathway. It showed good performance, but failed to attract attention, likely, in relation with lack of clear clinical benefit. Here, we have revised the value of the CP-MTX model directly addressing its clinical benefit by focusing on the expected benefit-cost of the predictions. In addition, our study included a much larger number of RA patients (n = 720) in MTX monotherapy than previous studies. Benefit of CP-MTX prediction was defined as the patients that would have received combination therapy as first treatment because they were correctly predicted as non-responders to MTX monotherapy. In contrast, cost of CP-MTX prediction was defined as the responder patients that were wrongly predicted as non-responders. Application of CP-MTX predictions to our patients showed a good benefit-cost relationship, with half of the 66.7% non-responders to MTX monotherapy rightly directed to alternative treatments (a benefit of 33.3%) at the cost of 8.5% wrongly predicted non-responders. These benefits-costs were consistent with reanalysis of the previously published studies. Therefore, predictions of CP-MTX showed a good benefit-cost relationship for informing MTX prescription.

A combined large-scale meta-analysis identifies COG6 as a novel shared risk locus for rheumatoid arthritis and systemic lupus erythematosus.

OBJECTIVES: During the last years, genome-wide association studies (GWASs) have identified a number of common genetic risk factors for rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE). However, the genetic overlap between these two immune-mediated diseases has not been thoroughly examined so far. The aim of the present study was to identify additional risk loci shared between RA and SLE. METHODS: We performed a large-scale meta-analysis of GWAS data from RA (3911 cases and 4083 controls) and SLE (2237 cases and 6315 controls). The top-associated polymorphisms in the discovery phase were selected for replication in additional datasets comprising 13 641 RA cases and 31 921 controls and 1957 patients with SLE and 4588 controls. RESULTS: The rs9603612 genetic variant, located nearby the COG6 gene, an established susceptibility locus for RA, reached genome-wide significance in the combined analysis including both discovery and replication sets (p value=2.95E-13). In silico expression quantitative trait locus analysis revealed that the associated polymorphism acts as a regulatory variant influencing COG6 expression. Moreover, protein-protein interaction and gene ontology enrichment analyses suggested the existence of overlap with specific biological processes, specially the type I interferon signalling pathway. Finally, genetic correlation and polygenic risk score analyses showed cross-phenotype associations between RA and SLE. CONCLUSIONS: In conclusion, we have identified a new risk locus shared between RA and SLE through a meta-analysis including GWAS datasets of both diseases. This study represents the first comprehensive large-scale analysis on the genetic overlap between these two complex disorders.

Replication study of polymorphisms associated with response to methotrexate in patients with rheumatoid arthritis.

About 70 genetic studies have already addressed the need of biomarkers to predict the response of patients with rheumatoid arthritis (RA) to methotrexate (MTX) treatment. However, no genetic biomarker has yet been sufficiently validated. Here, we aimed to replicate a selection of 25 SNPs in the largest collection of patients up to date, which consisted of 915 patients treated with MTX. The change in disease activity (measured as ΔDAS28) from baseline was considered the primary outcome. In addition, response according to widely used criteria (EULAR) was taken as secondary outcome. We considered consistency between outcomes, P values accounting for the number of SNPs, and independence from potential confounders for interpretation of the results. Only the rs1801394 SNP in MTRR fulfilled the high association standards. Its minor allele was associated with less improvement than the major allele according to ΔDAS28 (p = 0.0016), and EULAR response (p = 0.004), with independence of sex, age, baseline DAS28, smoking, seropositivity, concomitant corticosteroid use or previous treatments. In addition, previous evidence suggests the association of this SNP with response to MTX in another autoimmune disease, juvenile idiopathic arthritis, and with high intracellular folate levels, which could contribute to poor response.

Corrigendum: Simplified Assay for Epigenetic Age Estimation in Whole Blood of Adults.

[This corrects the article on p. 126 in vol. 7, PMID: 27471517.].

Anti-carbamylated protein autoantibodies associated with mortality in Spanish rheumatoid arthritis patients.

Patients with rheumatoid arthritis (RA) have an increased mortality rate that is associated with the presence of RA-specific autoantibodies in many studies. However, the relative role of rheumatoid factor (RF), anti-CCP antibodies and the most recently established RA-autoantibodies, directed against carbamylated proteins (anti-CarP antibodies), is unclear. Here, we have assessed the role of these three antibodies in 331 patients with established RA recruited from 2001 to 2009 and followed until November 2015. During this time, 124 patients died (37.5%). This death rate corresponds to a mortality rate 1.53 (95% CI 1.26 to 1.80) folds the observed in the reference population. We used for analysis of all-cause mortality the Cox proportional hazard regression model with adjustment for age, sex and smoking. It showed a trend for association with increased mortality of each of the three RA autoantibodies in antibody-specific analysis (hazards ratio (HR) from 1.37 to 1.79), but only the HR of the anti-CarP antibodies was significant (HR = 1.79, 95% CI 1.23 to 2.61, p = 0.002). In addition, the multivariate analysis that included all autoantibodies showed a marked decrease in the HR of RF and of anti-CCP antibodies, whereas the HR of anti-CarP remained significant. This increase was specific of respiratory system causes of death (HR = 3.19, 95% CI 1.52 to 6.69, p = 0.002). Therefore, our results suggest a specific relation of anti-CarP antibodies with the increased mortality in RA, and drive attention to their possible connection with respiratory diseases.

Simplified Assay for Epigenetic Age Estimation in Whole Blood of Adults.

Biological age is not always concordant with chronological age and the departures are of interest for understanding how diseases and environmental insults affect tissue function, organismal health, and life expectancy. The best-known biological age biomarker is telomere length, but there are more accurate biomarkers as the recently developed based in epigenetic, transcriptomic, or biochemical changes. The most accurate are the epigenetic biomarkers based on specific changes in DNA methylation referred as DNA methylation age measures (DmAM). Here, we have developed and validated a new DmAM that addresses some limitations of the previously available. The new DmAM includes the study in whole blood (WB) of 8 CpG sites selected as the most informative on a training set of 390 healthy subjects. The 8 CpG DmAM showed better accuracy than other DmAM based in few CpG in an independent validation set of 335 subjects. Results were not significantly influenced by sex, smoking, or variation in blood cell subpopulations. In addition, the new 8 CpG DmAM was amenable to study in a single multiplex reaction done with methylation-sensitive single-nucleotide primer extension (MS-SNuPE), a methodology based on commercially available reagents and run in capillary electrophoresis sequencers. In this way, the high cost of DNA methylation microarrays or of a pyrosequencer, which are needed for alternative DmAM, was avoided. Performance of the DmAM with MS-SNuPE was assessed in a set of 557 donors, showing high call rate (>97%), low CV (<3.3%) and high accuracy (Mean Absolute Deviation = 6.07 years). Therefore, the 8 CpG DmAM is a feasible and accurate tool for assessing the epigenetic component of biological age in blood of adults.

Specific premature epigenetic aging of cartilage in osteoarthritis.

Osteoarthritis (OA) is a disease affecting multiple tissues of the joints in the elderly, but most notably articular cartilage. Premature biological aging has been described in this tissue and in blood cells, suggesting a systemic component of premature aging in the pathogenesis of OA. Here, we have explored epigenetic aging in OA at the local (cartilage and bone) and systemic (blood) levels. Two DNA methylation age-measures (DmAM) were used: the multi-tissue age estimator for cartilage and bone; and a blood-specific biomarker for blood. Differences in DmAM between OA patients and controls showed an accelerated aging of 3.7 years in articular cartilage (95% CI = 1.1 to 6.3, P = 0.008) of OA patients. By contrast, no difference in epigenetic aging was observed in bone (0.04 years; 95% CI = -1.8 to 1.9, P = 0.3) and in blood (-0.6 years; 95% CI = -1.5 to 0.3, P = 0.2) between OA patients and controls. Therefore, premature epigenetic aging according to DNA methylation changes was specific of OA cartilage, adding further evidence and insight on premature aging of cartilage as a component of OA pathogenesis that reflects damage and vulnerability.

The genetic legacy of the pre-colonial period in contemporary Bolivians.

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.