Processing speed mediates the relationship between DDR1 and psychosocial functioning in euthymic patients with bipolar disorder presenting psychotic symptoms.

The DDR1 locus is associated with the diagnosis of schizophrenia and with processing speed in patients with schizophrenia and first-episode psychosis. Here, we investigated whether DDR1 variants are associated with bipolar disorder (BD) features. First, we performed a case‒control association study comparing DDR1 variants between patients with BD and healthy controls. Second, we performed linear regression analyses to assess the associations of DDR1 variants with neurocognitive domains and psychosocial functioning. Third, we conducted a mediation analysis to explore whether neurocognitive impairment mediated the association between DDR1 variants and psychosocial functioning in patients with BD. Finally, we studied the association between DDR1 variants and white matter microstructure. We did not find any statistically significant associations in the case‒control association study; however, we found that the combined genotypes rs1264323AA-rs2267641AC/CC were associated with worse neurocognitive performance in patients with BD with psychotic symptoms. In addition, the combined genotypes rs1264323AA-rs2267641AC/CC were associated with worse psychosocial functioning through processing speed. We did not find correlations between white matter microstructure abnormalities and the neurocognitive domains associated with the combined genotypes rs1264323AA-rs2267641AC/CC. Overall, the results suggest that DDR1 may be a marker of worse neurocognitive performance and psychosocial functioning in patients with BD, specifically those with psychotic symptoms.

Prevalence of non-coeliac gluten sensitivity and assessment of the response to gluten-free diet in a cohort of patients with fibromyalgia.

OBJECTIVES: To establish the prevalence of non-coeliac gluten sensitivity (NCGS) in a cohort of fibromyalgia patients and to evaluate their clinical response to a six-week gluten-free diet (GFD), the improvement in their symptoms, the percentage of diet responders who did not fulfil the diagnostic criteria for NCGS and the baseline characteristics that were associated with diet response and diagnostic criteria fulfilment. METHODS: Uncontrolled prospective experimental study in a cohort of patients with fibromyalgia from a specialized hospital unit. The percentage of patients that fulfilled the Salerno Experts' Criteria, that responded to GFD, that improved their symptomatology and baseline characteristics associated with GFD response and diagnostic criteria fulfilment was analysed. RESULTS: In total, 142 patients were selected and a NCGS prevalence of 5.6% was observed. A total of 21.8% responded to GFD due to their improvement in intestinal symptoms. In total, 74.2% of the responders did not fulfil the Salerno Experts' Criteria. The presence of diarrhoea and intraepithelial lymphocytosis and lower levels of anxiety were predictive factors of GFD response. No predictive factors of NCGS criteria fulfilment were found due to the low number of discriminators between gluten and placebo. CONCLUSIONS: A NCGS prevalence similar to that estimated in the general population was found. A GFD cannot be systematically recommended to all patients with fibromyalgia, although it could be evaluated in those with diarrhoea or intraepithelial lymphocytosis to evaluate if there are improvements in their intestinal symptoms.

Association of class II HLA alleles with susceptibility to develop immune-mediated diseases in Paraguayan patients.

Genetic and nongenetic factors are involved in the pathogenesis of immune-mediated inflammatory diseases (IMIDs). The best-known genetic factor for susceptibility to IMIDs is the human leukocyte antigen (HLA). The aim of the present study was to evaluate the association of HLA class II genes with the risk of systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and systemic sclerosis (SSc) in the Paraguayan population. We included 254 patients with IMIDs (101 SLE, 103 RA, and 50 SSc) and 50 healthy controls. The haplotypes of five genes corresponding to HLA class II genes and their relationship to the IMIDs studied were determined. Note that 84.6% were women, with a mean age of 43.4 ± 14 years. Among the associated HLA alleles, we found the previously identified risk factors in other populations like HLA-DRB1*03:01 and HLA-DRB1*14:02 for RA, as well as new ones not previously identified, such as DPA1*02:01 for SLE and, DB1*02:01 for RA and SSc. In the genetic association analysis, already known associations have been replicated, and unpublished associations have been identified in Paraguayan patients with IMIDs. This is the first genetic association study in Paraguayan patients with IMIDs.

Trajectories in early rheumatoid arthritis related fatigue over 10 years: results from the ESPOIR cohort.

OBJECTIVES: In a cohort of early rheumatoid arthritis (RA) patients, we aimed to determine and characterise fatigue trajectories over 10 years of follow-up and identify predictors of trajectory membership. METHODS: We selected patients fulfilling the 2010 ACR/EULAR criteria for RA included in the ESPOIR cohort. We used a cluster analysis to obtain fatigue (assessed by fatigue visual analogue scale) trajectories over the course of 10 years from enrolment. Chi-square tests or ANOVA were performed to evaluate differences of baseline variables between fatigue trajectories. Using a multinomial logistic regression we were able to identify predictors of trajectory membership. RESULTS: We analysed 598 patients with mean disease duration at enrolment of 26.2±40.9 days. Cluster analysis revealed 3 trajectories: high (18%), moderate (52%) and low fatigue (30%). Compared to patients with moderate or low fatigue trajectory, patients with high fatigue trajectory were predominantly women and reported significantly higher duration and intensity of morning stiffness, HAQ score, tender joints count, levels of pain, number of awakenings due to arthritis, frequency of fibromyalgic RA, levels of physician and patient global assessment, more frequent sleep problems, and increased psychological distress. Female patients with pain, psychological distress and presence of sicca symptoms had a higher risk of being in the high trajectory group. CONCLUSIONS: These findings suggest that levels of fatigue are rather stable over time in each trajectory. Baseline clinical measures and baseline patient-reported measures of functional status better distinguished the three fatigue trajectories. We did not find any differences between trajectories in baseline laboratory measures. Inflammatory activity was not a predictor of being in the high trajectory fatigue group.

The ovarian reserve as measured by the anti-Müllerian hormone is not diminished in patients with rheumatoid arthritis compared to the healthy population.

OBJECTIVES: Rheumatoid arthritis (RA) is the most prevalent chronic inflammatory arthritis, affecting 0.5-1% worldwide population and predominates in females. Altered fertility has been reported due to a decrease in ovarian reserve secondary to sustained inflammation. The anti-Müllerian Hormone (AMH) is currently the most reliable biomarker of ovarian reserve. However, few and contradictory studies have been reported to analyse the relationship between fertility in RA female patients and AMH. The aim of present study is to determine the AMH serum concentrations in a long-standing RA patient group and control group. We also sought to determine the correlation between AMH serum levels and disease activity measured by different parameters and the effect of biological DMARDs. METHODS: Serum AMH levels were measured in 60 women with long-standing RA aged 20-50 y.o. and compared to 59 healthy women. AMH was assessed by an electrochemiluminescence immunoassay method (ECLIA, Roche Diagnostics) and a large data set of clinical and molecular data was annotated. Demographic parameters, RA disease activity measured by DAS28 score and inflammatory biomarkers such as ESR, CRP, lymphocyte CD4+, CD8+, NK cells, IL-10 and IL-6 were determined. A comprehensive gynaecological self-administered questionnaire was given. Serum AMH levels were age-correlated. Differences between groups were calculated using Student's t-test or Mann-Whitney U-test for continuous variables and Fisher's exact test for categorical variables. Multivariate analysis was conducted by the partial correlation coefficient. Linear regression analysis was performed to study the effect of different variables on proportional AMH change. p-values <0.005 were considered significant. RESULTS: The median age was similar in AR and control groups (37.4±6.23 vs. 37.3±6.27 p=0.937). Mean disease duration was 8.37±5.36 years. The number of previous treatments was <3 in 71.7% of patients and ≥3 in 28.3%. Disease activity measured by DAS28 was 2.89± 1.54. The age-adjusted mean serum concentration of AMH was 1.27 ng/ml [IQR 0.42; 2.24] in RA patients and 1.31 ng/ml [IQR 0.46; 3.09] in controls (p=0.608). Neither disease activity (p=0.862), nor current or previous bDMARDs treatments (p=0.871) were associated with AMH levels. However, a negative linear correlation was observed between AMH and IL-10 levels (p=0.033). CONCLUSIONS: Our study shows that ovarian reserve determined by AMH serum levels is not reduced in rheumatoid arthritis patients compared with healthy controls. In our series, AMH levels were not affected by disease activity, however, a significant correlation was observed between AMH and IL-10 levels. These results support the role of cytokines profile in the female reproductive system and will focus further investigations in this critical area, mainly once biological DMARDs have been recommended in RA pregnant patients.

Dietary habits in patients with fibromyalgia: a cross-sectional study.

OBJECTIVES: To the scarce information on dietary habits in fibromyalgia (FM), it is added that there are no comparative studies with other rheumatic diseases. The objective of this study was to characterise the dietary habits of patients with FM by comparing, for the first time, with healthy controls (HC) and rheumatoid arthritis (RA). METHODS: This cross-sectional, observational study was based on data obtained from the Dietfibrom project for FM and from the IMID Consortium for RA and HC. All participants completed a food frequency questionnaire evaluating their weekly dietary intake of main food groups. The three cohorts were compared using a multiple logistic regression model adjusted for age, sex, and body mass index. RESULTS: After quality control, n=287 FM, n=1,983 HC and n=1,942 RA patients were analysed. We found that FM had a profound impact in the diet compared to HC, reducing the consumption of dairy (OR=0.32, p<0.0001), bread and/or whole grain cereals (OR=0.59, p=0.0006), fresh fruit (OR=0.66, P=0.008), and fish (OR=0.64, p=0.002). These same four food groups were also significantly reduced in FM patients in comparison to RA patients (p<0.0005 in all cases). Additionally, a lower consumption of pasta, rice and/or potatoes was also observed in FM compared to RA (OR=0.72, p=0.028). CONCLUSIONS: The present cross-sectional study shows that FM is associated to a significant change in the normal dietary patterns. These results underscore the importance of diet in this prevalent disease and are a warning of the potential long-range effects of a deficient nutritional status.

Interactions between rheumatoid arthritis antibodies are associated with the response to anti-tumor necrosis factor therapy.

BACKGROUND: Blocking of the Tumor Necrosis Factor (TNF) activity is a successful therapeutic approach for 50-60% of rheumatoid arthritis (RA) patients. However, there are yet no biomarkers to stratify patients for anti-TNF therapy. Rheumatoid factor (RF) and anti-cyclic-citrullinated antibodies (anti-CCP) have been evaluated as biomarkers of response but the results have shown limited consistency. Anti-carbamylated protein (anti-CarP) and anti-peptidylarginine deiminase type 4 (anti-PAD4) antibodies have been much less studied. Despite being linked to common immune processes, the interaction between these markers has not been evaluated yet. Our aim was to analyze the interaction between these four antibodies in relation to the response to anti-TNF therapy. METHODS: For this objective, a prospective cohort of n = 80 RA patients starting anti-TNF therapy was recruited. Serum determinations at baseline were performed for RF, anti-CCP, anti-CarP and anti-PAD4 antibodies using enzyme-linked immunosorbent assays (ELISA). The clinical response to anti-TNF therapy was determined at week 12 using the change in DAS28 score. Association was performed using multivariate linear regression adjusting for baseline DAS28, sex and age. RESULTS: The interaction between pairs of antibodies was tested by the addition of an interaction term. We found two highly significant antibody interactions associated with treatment response: anti-CarP with anti-PAD4 (p = 0.0062), and anti-CCP with RF (p = 0.00068). The latter antibody interaction was replicated in an independent retrospective cohort of RA patients (n = 199, p = 0.04). CONCLUSIONS: The results of this study suggest that antibody interaction effects are important factors in the response to anti-TNF therapy in RA.

A questionnaire-based study on contraceptive practice in patients with rheumatic disease found no significant difference in age-matched healthy controls.

OBJECTIVE: Birth control is crucial in preventing unplanned pregnancy. The study analyzed contraceptive practice in women and men with rheumatic disease. METHODS: A questionnaire-based study investigated the actual contraceptive practices in patients of reproductive age from three European countries and compared them to age-matched healthy women and men. Associations between patient characteristics and contraception behavior were analyzed by association analysis. RESULTS: No significant difference in the frequency of contraception use was found in 133 rheumatic patients compared to 122 healthy controls. The main reason for not using contraception was lack of partner or the wish to become pregnant, whereas the current use of contraception was predominantly to limit family size in general or at this stage of life. Both patients and controls preferred barrier methods (48% and 45%, respectively) followed by hormonal contraceptives (31% and 38%, respectively). Characteristics associated with less use of contraception in patients were living single, having no children, and for being religious, whereas gender and education had no influence. Treatment with teratogenic drugs was no major patient concern, and 13 of 30 female patients using methotrexate, mycophenolate mofetil, or leflunomide did not practice birth control. CONCLUSION: Patients used contraception less frequently than healthy individuals, and the main reason for use was to limit family size. Contraception should be an integral part of counseling patients of fertile age, since the patient-preferred methods in case of active disease or therapy with teratogenic drugs were unreliable for the prevention of pregnancy.

Epigenome-wide association study of rheumatoid arthritis identifies differentially methylated loci in B cells.

Epigenetic regulation of immune cell types could be critical for the development and maintenance of autoimmune diseases like rheumatoid arthritis (RA). B cells are highly relevant in RA, since patients express autoantibodies and depleting this cell type is a successful therapeutic approach. Epigenetic variation, such as DNA methylation, may mediate the pathogenic activity of B cells. In this study, we performed an epigenome-wide association study (EWAS) for RA with three different replication cohorts, to identify disease-specific alterations in DNA methylation in B cells. CpG methylation in isolated B lymphocytes was assayed on the Illumina HumanMethylation450 BeadChip in a discovery cohort of RA patients (N = 50) and controls (N = 75). Differential methylation was observed in 64 CpG sites (q < 0.05). Six biological pathways were also differentially methylated in RA B cells. Analysis in an independent cohort of patients (N = 15) and controls (N = 15) validated the association of 10 CpG sites located on 8 genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3 and TNFRSF13C, and 2 intergenic regions. Differential methylation at the CBL signaling pathway was replicated. Using an additional case-control cohort (N = 24), the association between RA risk and CpGs cg18972751 at CD1C (P = 2.26 × 10-9) and cg03055671 at TNFSF10 (P = 1.67 × 10-8) genes was further validated. Differential methylation at genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3, TNFRSF13C and intergenic region chr10p12.31 was replicated in a cohort of systemic lupus erythematosus (SLE) patients (N = 47) and controls (N = 56). Our results highlight genes that may drive the pathogenic activity of B cells in RA and suggest shared methylation patterns with SLE.

Genetic association between CD96 locus and immunogenicity to anti-TNF therapy in Crohn's disease.

The production of antibodies to anti-tumor necrosis factor alpha (TNF) agents is one of the main causes of treatment failure in Crohn's disease (CD). To date, however, the contribution of genetics to anti-TNF immunogenicity in CD is still unknown. The objective of the present study was to identify genetic variation associated with anti-TNF immunogenicity in CD. We performed a two-stage genome-wide association study in a cohort of 96 and 123 adalimumab-treated patients, respectively. In the discovery stage, we identified a genome-wide significant association between the CD96 locus and the production of antibodies to anti-TNF treatment (P = 1.88e-09). This association was validated in the replication stage (P < 0.05). The risk allele for anti-TNF immunogenicity was found to be also associated with a lack of response to anti-TNF therapy (P = 0.019). These findings represent an important step toward the understanding of the immunogenicity-based mechanisms that underlie anti-TNF response in CD.

Pro- and anti-inflammatory eicosanoids in psoriatic arthritis.

INTRODUCTION: Eicosanoids are biological lipids that serve as both activators and suppressors of inflammation. Eicosanoid pathways are implicated in synovitis and joint destruction in inflammatory arthritis, yet they might also have a protective function, underscoring the need for a comprehensive understanding of how eicosanoid pathways might be imbalanced. Until recently, sensitive and scalable methods for detecting and quantifying a high number of eicosanoids have not been available. OBJECTIVE: Here, we intend to describe a detailed eicosanoid profiling in patients with psoriatic arthritis (PsA) and evaluate correlations with parameters of disease activity. METHODS: Forty-one patients with PsA, all of whom satisfied the CASPAR classification criteria for PsA, were studied. Outcomes reflecting the activity of peripheral arthritis as well as skin psoriasis, Disease Activity Score (DAS)28, Clinical Disease Index (CDAI) and Body Surface Area (BSA) were assessed. Serum eicosanoids were determined by LC-MS, and the correlation between metabolite levels and disease scores was evaluated. RESULTS: Sixty-six eicosanoids were identified by reverse-phase LC/MS. Certain eicosanoids species including several pro-inflammatory eicosanoids such as PGE2, HXB3 or 6,15-dk,dh,PGF1a correlated with joint disease score. Several eicosapentaenoic acid (EPA)-derived eicosanoids, which associate with anti-inflammatory properties, such as 11-HEPE, 12-HEPE and 15-HEPE, correlated with DAS28 (Disease Activity Score) and CDAI (Clinical Disease Activity Index) as well. Of interest, resolvin D1, a DHA-derived anti-inflammatory eicosanoid, was down-regulated in patients with high disease activity. CONCLUSION: Both pro- and anti-inflammatory eicosanoids were associated with joint disease score, potentially representing pathways of harm as well as benefit. Further studies are needed to determine whether these eicosanoid species might also play a role in the pathogenesis of joint inflammation in PsA.

aRNApipe: a balanced, efficient and distributed pipeline for processing RNA-seq data in high-performance computing environments.

SUMMARY: The wide range of RNA-seq applications and their high-computational needs require the development of pipelines orchestrating the entire workflow and optimizing usage of available computational resources. We present aRNApipe, a project-oriented pipeline for processing of RNA-seq data in high-performance cluster environments. aRNApipe is highly modular and can be easily migrated to any high-performance computing (HPC) environment. The current applications included in aRNApipe combine the essential RNA-seq primary analyses, including quality control metrics, transcript alignment, count generation, transcript fusion identification, alternative splicing and sequence variant calling. aRNApipe is project-oriented and dynamic so users can easily update analyses to include or exclude samples or enable additional processing modules. Workflow parameters are easily set using a single configuration file that provides centralized tracking of all analytical processes. Finally, aRNApipe incorporates interactive web reports for sample tracking and a tool for managing the genome assemblies available to perform an analysis. AVAILABILITY AND DOCUMENTATION: https://github.com/HudsonAlpha/aRNAPipe ; DOI: 10.5281/zenodo.202950. CONTACT: rmyers@hudsonalpha.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Intravenous administration of expanded allogeneic adipose-derived mesenchymal stem cells in refractory rheumatoid arthritis (Cx611): results of a multicentre, dose escalation, randomised, single-blind, placebo-controlled phase Ib/IIa clinical trial.

OBJECTIVES: To evaluate the safety and tolerability of the intravenous administration of Cx611, a preparation of allogeneic expanded adipose-derived stem cells (eASCs), in patients with refractory rheumatoid arthritis (RA), as well as to obtain preliminary clinical efficacy data in this population. METHODS: It is a multicentre, dose escalation, randomised, single-blind (double-blind for efficacy), placebo-controlled, phase Ib/IIa clinical trial. Patients with active refractory RA (failure to at least two biologicals) were randomised to receive three intravenous infusions of Cx611: 1 million/kg (cohort A), 2 million/kg (cohort B), 4 million/kg (cohort C) or placebo, on days 1, 8 and 15, and they were followed for therapy assessment for 24 weeks. RESULTS: Fifty-three patients were treated (20 in cohort A, 20 in cohort B, 6 in cohort C and 7 in placebo group). A total of 141 adverse events (AEs) were reported. Seventeen patients from the group A (85%), 15 from the group B (75%), 6 from the group C (100%) and 4 from the placebo group (57%) experienced at least one AE.Eight AEs from 6 patients were grade 3 in intensity (severe), 5 in cohort A (lacunar infarction, diarrhoea, tendon rupture, rheumatoid nodule and arthritis), 2 in cohort B (sciatica and RA) and 1 in the placebo group (asthenia). Only one of the grade 3 AEs was serious (the lacunar infarction). American College of Rheumatology 20 responses for cohorts A, B, C and placebo were 45%, 20%, 33% and 29%, respectively, at month 1, and 25%, 15%, 17% and 0%, respectively, at month 3. CONCLUSIONS: The intravenous infusion of Cx611 was in general well tolerated, without evidence of dose-related toxicity at the dose range and time period studied. In addition, a trend for clinical efficacy was observed. These data, in our opinion, justify further investigation of this innovative therapy in patients with RA. TRIAL REGISTRATION NUMBERS: EudraCT: 2010-021602-37; NCT01663116; Results.

A genome-wide association study identifies a novel locus at 6q22.1 associated with ulcerative colitis.

The genetic analysis of ulcerative colitis (UC) has provided new insights into the etiology of this prevalent inflammatory bowel disease. However, most of the heritability of UC (>70%) has still not been characterized. To identify new risk loci for UC we have performed the first genome-wide association study (GWAS) in a Southern European population and undertaken a meta-analysis study combining the newly genotyped 825 UC patients and 1525 healthy controls from Spain with the six previously published GWAS comprising 6687 cases and 19 718 controls from Northern-European ancestry. We identified a novel locus with genome-wide significance at 6q22.1 [rs2858829, P = 8.97 × 10(-9), odds ratio (OR) (95% confidence interval, CI] = 1.12 (1.08-1.16)] that was validated with genotype data from a replication cohort of the same Southern European ancestry consisting in 1073 cases and 1279 controls [combined P = 7.59 × 10(-10), OR (95% CI) = 1.12 (1.08-1.16)]. Furthermore, we confirmed the association of 33 reported associations with UC and we nominally validated the GWAS results of nine new risk loci (P < 0.05, same direction of effect). SNP rs2858829 lies in an intergenic region and is a strong cis-eQTL for FAM26F gene, a gene that is shown to be selectively upregulated in UC colonic mucosa with active inflammation. Our results provide new insight into the genetic risk background of UC, confirming that there is a genetic risk component that differentiates from Crohn's Disease, the other major form of inflammatory bowel disease.

Identification of risk loci for Crohn's disease phenotypes using a genome-wide association study.

BACKGROUND & AIMS: Crohn's disease is a highly heterogeneous inflammatory bowel disease comprising multiple clinical phenotypes. Genome-wide association studies (GWASs) have associated a large number of loci with disease risk but have not associated any specific genetic variants with clinical phenotypes. We performed a GWAS of clinical phenotypes in Crohn's disease. METHODS: We genotyped 576,818 single-nucleotide polymorphisms in a well-characterized cohort of 1090 Crohn's disease patients of European ancestry. We assessed their association with 17 phenotypes of Crohn's disease (based on disease location, disease behavior, disease course, age at onset, and extraintestinal manifestations). A total of 57 markers with strong associations to Crohn's disease phenotypes (P < 2 × 10(-4)) were subsequently analyzed in an independent replication cohort of 1296 patients of European ancestry. RESULTS: We replicated the association of 4 loci with different Crohn's disease phenotypes. Variants in MAGI1, CLCA2, 2q24.1, and LY75 loci were associated with a complicated stricturing disease course (Pcombined = 2.01 × 10(-8)), disease location (Pcombined = 1.3 × 10(-6)), mild disease course (Pcombined = 5.94 × 10(-7)), and erythema nodosum (Pcombined = 2.27 × 10(-6)), respectively. CONCLUSIONS: In a GWAS, we associated 4 loci with clinical phenotypes of Crohn's disease. These findings indicate a genetic basis for the clinical heterogeneity observed for this inflammatory bowel disease.

Urine metabolome profiling of immune-mediated inflammatory diseases.

BACKGROUND: Immune-mediated inflammatory diseases (IMIDs) are a group of complex and prevalent diseases where disease diagnostic and activity monitoring is highly challenging. The determination of the metabolite profiles of biological samples is becoming a powerful approach to identify new biomarkers of clinical utility. In order to identify new metabolite biomarkers of diagnosis and disease activity, we have performed the first large-scale profiling of the urine metabolome of the six most prevalent IMIDs: rheumatoid arthritis, psoriatic arthritis, psoriasis, systemic lupus erythematosus, Crohn's disease, and ulcerative colitis. METHODS: Using nuclear magnetic resonance, we analyzed the urine metabolome in a discovery cohort of 1210 patients and 100 controls. Within each IMID, two patient subgroups were recruited representing extreme disease activity (very high vs. very low). Metabolite association analysis with disease diagnosis and disease activity was performed using multivariate linear regression in order to control for the effects of clinical, epidemiological, or technical variability. After multiple test correction, the most significant metabolite biomarkers were validated in an independent cohort of 1200 patients and 200 controls. RESULTS: In the discovery cohort, we identified 28 significant associations between urine metabolite levels and disease diagnosis and three significant metabolite associations with disease activity (P FDR < 0.05). Using the validation cohort, we validated 26 of the diagnostic associations and all three metabolite associations with disease activity (P FDR < 0.05). Combining all diagnostic biomarkers using multivariate classifiers we obtained a good disease prediction accuracy in all IMIDs and particularly high in inflammatory bowel diseases. Several of the associated metabolites were found to be commonly altered in multiple IMIDs, some of which can be considered as hub biomarkers. The analysis of the metabolic reactions connecting the IMID-associated metabolites showed an over-representation of citric acid cycle, phenylalanine, and glycine-serine metabolism pathways. CONCLUSIONS: This study shows that urine is a source of biomarkers of clinical utility in IMIDs. We have found that IMIDs show similar metabolic changes, particularly between clinically similar diseases and we have found, for the first time, the presence of hub metabolites. These findings represent an important step in the development of more efficient and less invasive diagnostic and disease monitoring methods in IMIDs.

A genome-wide association study identifies SLC8A3 as a susceptibility locus for ACPA-positive rheumatoid arthritis.

OBJECTIVE: RA patients with serum ACPA have a strong and specific genetic background. The objective of the study was to identify new susceptibility genes for ACPA-positive RA using a genome-wide association approach. METHODS: A total of 924 ACPA-positive RA patients with joint damage in hands and/or feet, and 1524 healthy controls were genotyped in 582 591 single-nucleotide polymorphisms (SNPs) in the discovery phase. In the validation phase, the most significant SNPs in the genome-wide association study representing new candidate loci for RA were tested in an independent cohort of 863 ACPA-positive patients with joint damage and 1152 healthy controls. All individuals from the discovery and validation cohorts were Caucasian and of Southern European ancestry. RESULTS: In the discovery phase, 60 loci not previously associated with RA risk showed evidence for association at P < 5×10(-4) and were tested for replication in the validation cohort. A total of 12 loci were replicated at the nominal level (P < 0.05, same direction of effect as in the discovery phase). When combining the discovery and validation cohorts, an intronic SNP in the Solute Carrier family 8 gene (SLC8A3) was found to be associated with ACPA-positive RA at a genome-wide level of significance RA [odds ratio (95% CI): 1.42 (1.25, 1.6), Pcombined = 3.19×10(-8)]. CONCLUSIONS: SLC8A3 was identified as a new risk locus for ACPA-positive RA. This study demonstrates the advantage of analysing relevant subsets of RA patients to identify new genetic risk variants.

A deletion at ADAMTS9-MAGI1 locus is associated with psoriatic arthritis risk.

OBJECTIVE: Copy number variants (CNVs) have been associated with the risk to develop multiple autoimmune diseases. Our objective was to identify CNVs associated with the risk to develop psoriatic arthritis (PsA) using a genome-wide analysis approach. METHODS: A total of 835 patients with PsA and 1498 healthy controls were genotyped for CNVs using the Illumina HumanHap610 BeadChip genotyping platform. Genomic CNVs were characterised using CNstream analysis software and analysed for association using the χ(2) test. The most significant genomic CNV associations with PsA risk were independently tested in a validation sample of 1133 patients with PsA and 1831 healthy controls. In order to test for the specificity of the variants with PsA aetiology, we also analysed the association to a cohort of 822 patients with purely cutaneous psoriasis (PsC). RESULTS: A total of 165 common CNVs were identified in the genome-wide analysis. We found a highly significant association of an intergenic deletion between ADAMTS9 and MAGI1 genes on chromosome 3p14.1 (p=0.00014). Using the independent patient and control cohort, we validated the association between ADAMTS9-MAGI1 deletion and PsA risk (p=0.032). Using next-generation sequencing, we characterised the 26 kb associated deletion. Finally, analysing the PsC cohort we found a lower frequency of the deletion compared with the PsA cohort (p=0.0088) and a similar frequency to that of healthy controls (p>0.3). CONCLUSIONS: The present genome-wide scan for CNVs associated with PsA risk has identified a new deletion associated with disease risk and which is also differential from PsC risk.

Risk variants for psoriasis vulgaris in a large case-control collection and association with clinical subphenotypes.

Recent genome-wide association studies (GWASs) have identified >20 new loci associated with the susceptibility to psoriasis vulgaris (PsV) risk. We investigated the association of PsV and its main clinical subphenotypes with 32 loci having previous genome-wide evidence of association with PsV (P < 5e-8) or strong GWAS evidence (P < 5e-5 in discovery and P < 0.05 in replication sample) in a large cohort of PsV patients (n = 2005) and controls (n = 1497). We provide the first independent replication for COG6 (P = 0.00079) and SERPINB8 (P = 0.048) loci with PsV. In those patients having developed psoriatic arthritis (n = 955), we found, for the first time, a strong association with IFIH1 (P = 0.013). Analyses of clinically relevant PsV subtypes yielded a significant association of severity of cutaneous disease with variation at LCE3D locus (P = 0.0005) in PsV and nail involvement with IL1RN in purely cutaneous psoriasis (PsC, P = 0.007). In an exploratory analysis of epistasis, we replicated the previously described HLA-C-ERAP1 interaction with PsC. Our findings show that common genetic variants associated with a complex phenotype like PsV influence different subphenotypes of high clinical relevance.