scDrugPrio: a framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

BACKGROUND: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. METHODS: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. RESULTS: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. CONCLUSIONS: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package ( https://github.com/SDTC-CPMed/scDrugPrio ).

Seven-chain adaptive immune receptor repertoire analysis in rheumatoid arthritis reveals novel features associated with disease and clinically relevant phenotypes.

BACKGROUND: In rheumatoid arthritis (RA), the activation of T and B cell clones specific for self-antigens leads to the chronic inflammation of the synovium. Here, we perform an in-depth quantitative analysis of the seven chains that comprise the adaptive immune receptor repertoire (AIRR) in RA. RESULTS: In comparison to controls, we show that RA patients have multiple and strong differences in the B cell receptor repertoire including reduced diversity as well as altered isotype, chain, and segment frequencies. We demonstrate that therapeutic tumor necrosis factor inhibition partially restores this alteration but find a profound difference in the underlying biochemical reactivities between responders and non-responders. Combining the AIRR with HLA typing, we identify the specific T cell receptor repertoire associated with disease risk variants. Integrating these features, we further develop a molecular classifier that shows the utility of the AIRR as a diagnostic tool. CONCLUSIONS: Simultaneous sequencing of the seven chains of the human AIRR reveals novel features associated with the disease and clinically relevant phenotypes, including response to therapy. These findings show the unique potential of AIRR to address precision medicine in immune-related diseases.

Mechanisms by which the cystic fibrosis transmembrane conductance regulator may influence SARS-CoV-2 infection and COVID-19 disease severity.

Patients with cystic fibrosis (CF) exhibit pronounced respiratory damage and were initially considered among those at highest risk for serious harm from SARS-CoV-2 infection. Numerous clinical studies have subsequently reported that individuals with CF in North America and Europe-while susceptible to severe COVID-19-are often spared from the highest levels of virus-associated mortality. To understand features that might influence COVID-19 among patients with cystic fibrosis, we studied relationships between SARS-CoV-2 and the gene responsible for CF (i.e., the cystic fibrosis transmembrane conductance regulator, CFTR). In contrast to previous reports, we found no association between CFTR carrier status (mutation heterozygosity) and more severe COVID-19 clinical outcomes. We did observe an unexpected trend toward higher mortality among control individuals compared with silent carriers of the common F508del CFTR variant-a finding that will require further study. We next performed experiments to test the influence of homozygous CFTR deficiency on viral propagation and showed that SARS-CoV-2 production in primary airway cells was not altered by the absence of functional CFTR using two independent protocols. On the contrary, experiments performed in vitro strongly indicated that virus proliferation depended on features of the mucosal fluid layer known to be disrupted by absent CFTR in patients with CF, including both low pH and increased viscosity. These results point to the acidic, viscous, and mucus-obstructed airways in patients with cystic fibrosis as unfavorable for the establishment of coronaviral infection. Our findings provide new and important information concerning relationships between the CF clinical phenotype and severity of COVID-19.

Longitudinal analysis of blood DNA methylation identifies mechanisms of response to tumor necrosis factor inhibitor therapy in rheumatoid arthritis.

BACKGROUND: Rheumatoid arthritis (RA) is a chronic, immune-mediated inflammatory disease of the joints that has been associated with variation in the peripheral blood methylome. In this study, we aim to identify epigenetic variation that is associated with the response to tumor necrosis factor inhibitor (TNFi) therapy. METHODS: Peripheral blood genome-wide DNA methylation profiles were analyzed in a discovery cohort of 62 RA patients at baseline and at week 12 of TNFi therapy. DNA methylation of individual CpG sites and enrichment of biological pathways were evaluated for their association with drug response. Using a novel cell deconvolution approach, altered DNA methylation associated with TNFi response was also tested in the six main immune cell types in blood. Validation of the results was performed in an independent longitudinal cohort of 60 RA patients. FINDINGS: Treatment with TNFi was associated with significant longitudinal peripheral blood methylation changes in biological pathways related to RA (FDR<0.05). 139 biological functions were modified by therapy, with methylation levels changing systematically towards a signature similar to that of healthy controls. Differences in the methylation profile of T cell activation and differentiation, GTPase-mediated signaling, and actin filament organization pathways were associated with the clinical response to therapy. Cell type deconvolution analysis identified CpG sites in CD4+T, NK, neutrophils and monocytes that were significantly associated with the response to TNFi. INTERPRETATION: Our results show that treatment with TNFi restores homeostatic blood methylation in RA. The clinical response to TNFi is associated to methylation variation in specific biological pathways, and it involves cells from both the innate and adaptive immune systems. FUNDING: The Instituto de Salud Carlos III.

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.

Synovial Immunohistological Biomarkers of the Classification of Undifferentiated Arthritis Evolving to Rheumatoid or Psoriatic Arthritis.

Background: Undifferentiated arthritis (UA) is defined as an inflammatory arthritis that does not fulfill criteria for a definite diagnosis. Delay in reaching a specific diagnostic and therapy may lead to impaired functional outcomes. Our aim was to identify synovial biomarkers associated with definitive diagnostic classification in patients with UA. Methods: DMARD-naïve UA patients with available initial synovial tissue (ST) and a final diagnosis of rheumatoid arthritis (RA) or psoriatic arthritis (PsA) during follow-up were included and compared with patients with well-defined disease (RA or PsA). Clinical, arthroscopic, and pathological data were compared between groups. Pathology included quantitative immunohistochemical (IHC) analysis of cell types and human interferon-regulated MxA. Principal component analysis (PCA) was performed to extract disease patterns. Results: One hundred and five patients were included: 31 patients with DMARD-naïve UA (19 evolving to RA and 12 to PsA during a median follow up of 7 years), 39 with established RA, and 35 with established PsA. ST from the UA group showed higher macrophage density compared with the established RA and PsA groups. Patients with UA evolving to RA (UA-RA) showed higher MxA expression and CD3+ T-cell density compared with established RA. UA patients evolving to PsA (UA-PsA) showed increased vascularity and lining synovial fibroblast density compared with established PsA. Synovitis of UA-PsA patients showed more mast cells and lining fibroblasts compared with UA-RA. No between-group differences in local or systemic inflammation markers were found. Conclusions: Our results show differences in the cellular composition of UA synovium compared with RA and PsA, with higher density of the cellular infiltrate in the UA groups. Initial expression of the interferon inducible gene MxA could be a biomarker of progression to RA, while higher mast cell and fibroblastic density may be associated with PsA progression.

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.

Food groups associated with immune-mediated inflammatory diseases: a Mendelian randomization and disease severity study.

BACKGROUND/OBJECTIVES: Immune-mediated inflammatory diseases (IMIDs) are prevalent diseases. There is, however, a lack of understanding of the link between diet and IMIDs, how much dietary patterns vary between them and if there are food groups associated with a worsening of the disease. SUBJECTS/METHODS: To answer these questions we analyzed a nation-wide cohort of n = 11,308 patients from six prevalent IMIDs and 2050 healthy controls. We compared their weekly intake of the major food categories, and used a Mendelian randomization approach to determine which dietary changes are caused by disease. Within each IMID, we analyzed the association between food frequency and disease severity. RESULTS: After quality control, n = 11,230 recruited individuals were used in this study. We found that diet is profoundly altered in all IMIDs: at least three food categories are significantly altered in each disease (P < 0.05). Inflammatory bowel diseases showed the largest differences compared to controls (n ≥ 8 categories, P < 0.05). Mendelian randomization analysis supported that some of these dietary changes, like vegetable reduction in Crohn's Disease (P = 2.5 × 10-10, OR(95% CI) = 0.73(0.65, 0.80)), are caused by the disease. Except for Psoriatic Arthritis and Systemic Lupus Erythematosus, we have found ≥2 food groups significantly associated with disease severity in the other IMIDs (P < 0.05). CONCLUSIONS: This cross-disease study demonstrates that prevalent IMIDs are associated to a significant change in the normal dietary patterns. This variation is highly disease-specific and, in some cases, it is caused by the disease itself. Severity in IMIDs is also associated with specific food groups. The results of this study underscore the importance of studying diet in IMIDs.

Lower peripheral helper T cell levels in the synovium are associated with a better response to anti-TNF therapy in rheumatoid arthritis.

BACKGROUND: The mechanisms by which only some rheumatoid arthritis (RA) patients respond favorably to TNF blockade are still poorly characterized. The goal of this study was to identify biological features that explain this differential response using a multilevel transcriptome analysis of the synovial membrane. METHODS: Synovial samples from 11 patients on anti-TNF therapy were obtained by arthroscopy at baseline and week 20. Analysis of the synovial transcriptome was performed at the gene, pathway, and cell-type levels. Newly characterized pathogenic cell types in RA, peripheral helper T cells (TPH), and CD34-THY1+ fibroblasts were estimated using a cell-type deconvolution approach. TPH association was validated using immunofluorescence. External validation was performed on an independent dataset. RESULTS: After multiple-test correction, 16 and 4 genes were differentially expressed at baseline and week 20, respectively. At the pathway level, 86 and 17 biological processes were significantly enriched at baseline and week 20, respectively. Longitudinal expression changes were associated with a drastic decrease of innate immune activity (P < 5e-30), and an activation of the bone and cartilage regeneration processes (P < 5e-10). Cell-type deconvolution revealed a significant association between low TPH cells at baseline and a better response (P = 0.026). Lower TPH cells were maintained in good responders up to week 20 (P = 0.032). Immunofluorescent analyses confirmed the accuracy of the cell-type estimation (r2 = 0.58, P = 0.005) and an association with response. TPH association with anti-TNF response was validated in an independent sample of RA patients (P = 0.0040). CONCLUSIONS: A lower abundance in the synovial membrane of the pathogenic T cell type newly associated with RA, peripheral helper T lymphocyte, is associated with a good response to anti-TNF therapy. Major changes in the myeloid cell compartment were also observed in response to therapy. The results of this study could help develop more effective therapies aimed at treating the pathogenic mechanisms in RA that are currently not well targeted by anti-TNF agents.

Transcriptomic identification of TMIGD1 and its relationship with the ileal epithelial cell differentiation in Crohn's disease.

Crohn's disease (CD) is a complex and multifactorial illness. There are still considerable gaps in our knowledge regarding its pathophysiology. A transcriptomic approach could shed some light on little-known biological alterations of the disease. We therefore aimed to explore the ileal transcriptome to gain knowledge about CD. We performed whole transcriptome gene expression analysis on ileocecal resections from CD patients and inflammatory bowel disease-free controls, as well as on a CD-independent cohort to replicate selected results. Normalized data were hierarchically clustered, and gene ontology and the molecular network were studied. Cell cultures and molecular methods were used for further evaluations. Genome-wide expression data analysis identified a robust transmembrane immunoglobulin domain-containing 1 (TMIGD1) gene underexpression in CD tissue, which was even more marked in inflamed ileum, and which was replicated in the validation cohort. Immunofluorescence showed TMIGD1 to be located in the apical microvilli of well-differentiated enterocytes but not in intestinal crypt. This apical TMIGD1 was lower in the noninflamed tissue and almost disappeared in the inflamed mucosa of surgical resections. In vitro studies showed hypoxic-dependent TMIGD1 decreased its expression in enterocyte-like cells. The gene enrichment analysis linked TMIGD1 with cell recovery and tissue remodeling in CD settings, involving guanylate cyclase activities. Transcriptomics may be useful for finding new targets that facilitate studies of the CD pathology. This is how TMIGD1 was identified in CD patients, which was related to multiciliate ileal epithelial cell differentiation.NEW & NOTEWORTHY This is a single-center translational research study that aimed to look for key targets involved in Crohn's disease and define molecular pathways through different functional analysis strategies. With this approach, we have identified and described a novel target, the almost unknown TMIGD1 gene, which may be key in the recovery of injured mucosa involving intestinal epithelial cell differentiation.

Sampling time-dependent artifacts in single-cell genomics studies.

Robust protocols and automation now enable large-scale single-cell RNA and ATAC sequencing experiments and their application on biobank and clinical cohorts. However, technical biases introduced during sample acquisition can hinder solid, reproducible results, and a systematic benchmarking is required before entering large-scale data production. Here, we report the existence and extent of gene expression and chromatin accessibility artifacts introduced during sampling and identify experimental and computational solutions for their prevention.

A cross-disease meta-GWAS identifies four new susceptibility loci shared between systemic sclerosis and Crohn's disease.

Genome-wide association studies (GWASs) have identified a number of genetic risk loci associated with systemic sclerosis (SSc) and Crohn's disease (CD), some of which confer susceptibility to both diseases. In order to identify new risk loci shared between these two immune-mediated disorders, we performed a cross-disease meta-analysis including GWAS data from 5,734 SSc patients, 4,588 CD patients and 14,568 controls of European origin. We identified 4 new loci shared between SSc and CD, IL12RB2, IRF1/SLC22A5, STAT3 and an intergenic locus at 6p21.31. Pleiotropic variants within these loci showed opposite allelic effects in the two analysed diseases and all of them showed a significant effect on gene expression. In addition, an enrichment in the IL-12 family and type I interferon signaling pathways was observed among the set of SSc-CD common genetic risk loci. In conclusion, through the first cross-disease meta-analysis of SSc and CD, we identified genetic variants with pleiotropic effects on two clinically distinct immune-mediated disorders. The fact that all these pleiotropic SNPs have opposite allelic effects in SSc and CD reveals the complexity of the molecular mechanisms by which polymorphisms affect diseases.

Functional rare variants influence the clinical response to anti-TNF therapy in Crohn's disease.

BACKGROUND: The effect of low-frequency functional variation on anti-tumor necrosis factor alpha (TNF) response in Crohn's disease (CD) patients remains unexplored. The objective of this study was to investigate the impact of functional rare variants in clinical response to anti-TNF therapy in CD. METHODS: CD anti-TNF naïve patients starting anti-TNF treatment due to active disease [Crohn's Disease Activity Index (CDAI > 150)] were included. The whole genome was sequenced using the Illumina Hiseq4000 platform. Clinical response was defined as a CDAI score <150 at week 14 of anti-TNF treatment. Low-frequency variants were annotated and classified according to their damaging potential. The whole genome of CD patients was screened to identify homozygous loss-of-function (LoF) variants. The TNF signaling pathway was tested for overabundance of damaging variants using the SKAT-O method. Functional implication of the associated rare variation was evaluated using cell-type epigenetic enrichment analyses. RESULTS: A total of 41 consecutive CD patients were included; 3250 functional rare variants were identified (2682 damaging and 568 LoF variants). Two homozygous LoF mutations were found in HLA-B and HLA-DRB1 genes associated with lack of response and remission, respectively. Genome-wide LoF variants were enriched in epigenetic marks specific for the gastrointestinal tissue (colon, p = 4.11e-4; duodenum, p = 0.011). The burden of damaging variation in the TNF signaling pathway was associated with response to anti-TNF therapy (p = 0.016); damaging variants were enriched in epigenetic marks from CD8+ (p = 6.01e-4) and CD4+ (p = 0.032) T cells. CONCLUSIONS: Functional rare variants are involved in the response to anti-TNF therapy in CD. Cell-type enrichment analysis suggests that the gut mucosa and CD8+ T cells are the main mediators of this response.

A Combined Transcriptomic and Genomic Analysis Identifies a Gene Signature Associated With the Response to Anti-TNF Therapy in Rheumatoid Arthritis.

Background: Rheumatoid arthritis (RA) is the most frequent autoimmune disease involving the joints. Although anti-TNF therapies have proven effective in the management of RA, approximately one third of patients do not show a significant clinical response. The objective of this study was to identify new genetic variation associated with the clinical response to anti-TNF therapy in RA. Methods: We performed a sequential multi-omic analysis integrating different sources of molecular information. First, we extracted the RNA from synovial biopsies of 11 RA patients starting anti-TNF therapy to identify gene coexpression modules (GCMs) in the RA synovium. Second, we analyzed the transcriptomic association between each GCM and the clinical response to anti-TNF therapy. The clinical response was determined at week 14 using the EULAR criteria. Third, we analyzed the association between the GCMs and anti-TNF response at the genetic level. For this objective, we used genome-wide data from a cohort of 348 anti-TNF treated patients from Spain. The GCMs that were significantly associated with the anti-TNF response were then tested for validation in an independent cohort of 2,706 anti-TNF treated patients. Finally, the functional implication of the validated GCMs was evaluated via pathway and cell type epigenetic enrichment analyses. Results: A total of 149 GCMs were identified in the RA synovium. From these, 13 GCMs were found to be significantly associated with anti-TNF response (P < 0.05). At the genetic level, we detected two of the 13 GCMs to be significantly associated with the response to adalimumab (P = 0.0015) and infliximab (P = 0.021) in the Spain cohort. Using the independent cohort of RA patients, we replicated the association of the GCM associated with the response to adalimumab (P = 0.0019). The validated module was found to be significantly enriched for genes involved in the nucleotide metabolism (P = 2.41e-5) and epigenetic marks from immune cells, including CD4+ regulatory T cells (P = 0.041). Conclusions: These findings show the existence of a drug-specific genetic basis for anti-TNF response, thereby supporting treatment stratification in the search for response biomarkers in RA.

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.

Genome-wide pathway analysis identifies VEGF pathway association with oral ulceration in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a genetically complex rheumatic disease characterized by heterogeneous clinical manifestations of unknown etiology. Recent studies have suggested the existence of a genetic basis for SLE heterogeneity. The objective of the present study was to identify new genetic variation associated with the clinically relevant phenotypes in SLE. METHODS: A two-stage pathway-based approach was used to identify the genetic variation associated with the main clinical phenotypes in SLE. In the discovery stage, 482 SLE patients were genotyped using Illumina Human Quad610 microarrays. Association between 798 reference genetic pathways from the Molecular Signatures Database and 11 SLE phenotypes was tested using the set-based method implemented in PLINK software. Pathways significantly associated after multiple test correction were subsequently tested for replication in an independent cohort of 425 SLE patients. Using an in silico approach, we analyzed the functional effects of common SLE therapies on the replicated genetic pathways. The association of known SLE risk variants with the development of the clinical phenotypes was also analyzed. RESULTS: In the discovery stage, we found a significant association between the vascular endothelial growth factor (VEGF) pathway and oral ulceration (P value for false discovery rate (P FDR) < 0.05), and between the negative regulation signaling pathway of retinoic acid inducible gene-I/melanoma differentiation associated gene 5 and the production of antinuclear antibodies (P FDR < 0.05). In the replication stage, we validated the association between the VEGF pathway and oral ulceration. Therapies commonly used to treat mucocutaneous phenotypes in SLE were found to strongly influence VEGF pathway gene expression (P = 4.60e-4 to 5.38e-14). Analysis of known SLE risk loci identified a strong association between PTPN22 and the risk of hematologic disorder and with the development of antinuclear antibodies. CONCLUSIONS: The present study has identified VEGF genetic pathway association with the risk of oral ulceration in SLE. New therapies targeting the VEGF pathway could be more effective in reducing the severity of this phenotype. These findings represent a first step towards the understanding of the genetic basis of phenotype heterogeneity in SLE.

A functional variant of TLR10 modifies the activity of NFkB and may help predict a worse prognosis in patients with rheumatoid arthritis.

BACKGROUND: Toll-like receptor (TLR) family members are key players in inflammation. TLR10 has been poorly studied in chronic inflammatory disorders, and its clinical relevance in rheumatoid arthritis (RA) is as yet unknown. We aimed at identifying TLR10 variants within all coding regions of the gene in patients with RA as well as studying their functional and clinical significance. METHODS: TLR10 gene variants were studied by performing sequencing of 66 patients with RA and 30 control subjects. A selected variant, I473T, was then analyzed in 1654 patients and 1702 healthy control subjects. The capacity of this TLR10 variant to modify the transcriptional activity of nuclear factor kappa-light-chain-enhancer of activated B cells (NFkB) was determined by using a luciferase reporter assay and analyzing the expression of NFkB target genes by quantitative polymerase chain reaction. Differences between groups were analyzed by using the Mann-Whitney U test and the unpaired two-tailed Student's t test. RESULTS: We detected ten missense variants in the TLR10 gene and focused on the I473T substitution based on allele frequencies and the predicted functional impact. I473T variant is not associated with susceptibility to RA, but it significantly correlates with erosive disease in patients seropositive for antibodies to citrullinated protein antigens (p = 0.017 in the total cohort and p = 0.0049 in female patients) and with a lower response to infliximab treatment as measured by the change in Disease Activity Score in 28 joints (p = 0.012) and by the European League Against Rheumatism criteria (p = 0.049). Functional studies showed that TLR10 reduced activation of the NFkB inflammatory pathway in hematopoietic cells, whereas the I473T variant lacked this inhibitory capacity. Consistently, after exposure to infliximab, cells expressing the I437T variant showed higher NFkB activity than cells carrying wild-type TLR10. CONCLUSIONS: A TLR10 allelic variant, I473T, has impaired NFkB inhibitory activity and is highly associated with disease severity and low response to infliximab in patients with RA.

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 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.