Diminished cytokine-induced Jak/STAT signaling is associated with rheumatoid arthritis and disease activity.

Rheumatoid arthritis (RA) is a systemic and incurable autoimmune disease characterized by chronic inflammation in synovial lining of joints. To identify the signaling pathways involved in RA, its disease activity, and treatment response, we adapted a systems immunology approach to simultaneously quantify 42 signaling nodes in 21 immune cell subsets (e.g., IFNα→p-STAT5 in B cells) in peripheral blood mononuclear cells (PBMC) from 194 patients with longstanding RA (including 98 patients before and after treatment), and 41 healthy controls (HC). We found multiple differences between patients with RA compared to HC, predominantly in cytokine-induced Jak/STAT signaling in many immune cell subsets, suggesting pathways that may be associated with susceptibility to RA. We also found that high RA disease activity, compared to low disease activity, was associated with decreased (e.g., IFNα→p-STAT5, IL-10→p-STAT1) or increased (e.g., IL-6→STAT3) response to stimuli in multiple cell subsets. Finally, we compared signaling in patients with established, refractory RA before and six months after initiation of methotrexate (MTX) or TNF inhibitors (TNFi). We noted significant changes from pre-treatment to post-treatment in IFNα→p-STAT5 signaling and IL-10→p-STAT1 signaling in multiple cell subsets; these changes brought the aberrant RA signaling profiles toward those of HC. This large, comprehensive functional signaling pathway study provides novel insights into the pathogenesis of RA and shows the potential of quantification of cytokine-induced signaling as a biomarker of disease activity or treatment response.

Autoantigen can promote progression to a more aggressive TCL1 leukemia by selecting variants with enhanced B-cell receptor signaling.

(Auto)antigen engagement by the B-cell receptor (BCR) and possibly the sites where this occurs influence the outcome of chronic lymphocytic leukemia (CLL). To test if selection for autoreactivity leads to increased aggressiveness and if this selection plays out equally in primary and secondary tissues, we used T-cell leukemia (TCL)1 cells reactive with the autoantigen phosphatidylcholine (PtC). After repeated transfers of splenic lymphocytes from a single mouse with oligoclonal PtC-reactive cells, outgrowth of cells expressing a single IGHV-D-J rearrangement and superior PtC-binding and disease virulence occurred. In secondary tissues, increased PtC-binding correlated with enhanced BCR signaling and cell proliferation, whereas reduced signaling and division of cells from the same clone was documented in cells residing in the bone marrow, blood, and peritoneum, even though cells from the last site had highest surface membrane IgM density. Gene-expression analyses revealed reciprocal changes of genes involved in BCR-, CD40-, and PI3K-signaling between splenic and peritoneal cells. Our results suggest autoantigen-stimulated BCR signaling in secondary tissues promotes selection, expansion, and disease progression by activating pro-oncogenic signaling pathways, and that--outside secondary lymphoid tissues--clonal evolution is retarded by diminished BCR-signaling. This transferrable, antigenic-specific murine B-cell clone (TCL1-192) provides a platform to study the types and sites of antigen-BCR interactions and genetic alterations that result and may have relevance to patients.

Hypogalactosylation of serum N-glycans fails to predict clinical response to methotrexate and TNF inhibition in rheumatoid arthritis.

INTRODUCTION: Rheumatoid arthritis (RA) is associated with hypogalactosylation of immunoglobulin G (IgG). We examined whether a proxy measure for galactosylation of IgG N-glycans could predict response to therapy or was differentially affected by methotrexate (MTX) or TNF blockade. METHODS: Using a previously defined normal phase high-performance liquid chromatography approach, we ascertained the galactosylation status of whole serum N-glycans in two well-defined RA clinical cohorts: the Autoimmune Biomarkers Collaborative Network (n = 98) and Nested I (n = 64). The ratio of agalactosylated to monogalactosylated N-glycans in serum (sG0/G1) was determined before and during therapy with MTX or TNF inhibition and correlated with anticitrullinated peptide antibody (ACPA) status and clinical response as assessed by 28-joint Disease Activity Score utilizing C-reactive peptide and European League Against Rheumatism response criteria. RESULTS: RA patients from both cohorts exhibited elevation of sG0/G1 at baseline. Improvement in clinical scores correlated with a reduction in sG0/G1 (Spearman's ρ = 0.31 to 0.37; P < 0.05 for each cohort). However, pretreatment sG0/G1 was not predictive of clinical response. Changes in sG0/G1 were similar in the MTX and TNF inhibitor groups. Corrected for disease activity, ACPA positivity correlated with higher sG0/G1. CONCLUSIONS: Baseline serum N-glycan hypogalactosylation, an index previously correlated with hypogalactosylation of IgG N-glycans, did not distinguish patients with rheumatoid arthritis who were likely to experience a favorable clinical response to MTX or TNF blockade. Clinical improvement was associated with partial glycan normalization. ACPA-positive patients demonstrated enhanced N-glycan aberrancy compared with ACPA-negative patients.

Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor alpha therapy.

OBJECTIVE: Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy. METHODS: A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models. RESULTS: Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99). CONCLUSION: Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.

Interferon-regulated chemokines as biomarkers of systemic lupus erythematosus disease activity: a validation study.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by unpredictable flares of disease activity and irreversible damage to multiple organ systems. An earlier study showed that SLE patients carrying an interferon (IFN) gene expression signature in blood have elevated serum levels of IFN-regulated chemokines. These chemokines were associated with more-severe and active disease and showed promise as SLE disease activity biomarkers. This study was designed to validate IFN-regulated chemokines as biomarkers of SLE disease activity in 267 SLE patients followed up longitudinally. METHODS: To validate the potential utility of serum chemokine levels as biomarkers of disease activity, we measured serum levels of CXCL10 (IFNgamma-inducible 10-kd protein), CCL2 (monocyte chemotactic protein 1), and CCL19 (macrophage inflammatory protein 3beta) in an independent cohort of 267 SLE patients followed up longitudinally over 1 year (1,166 total clinic visits). RESULTS: Serum chemokine levels correlated with lupus activity at the current visit (P = 2 x 10(-10)), rising at the time of SLE flare (P = 2 x 10(-3)) and decreasing as disease remitted (P = 1 x 10(-3)); they also performed better than the currently available laboratory tests. Chemokine levels measured at a single baseline visit in patients with a Systemic Lupus Erythematosus Disease Activity Index of < or =4 were predictive of lupus flare over the ensuing year (P = 1 x 10(-4)). CONCLUSION: Monitoring serum chemokine levels in SLE may improve the assessment of current disease activity, the prediction of future disease flares, and the overall clinical decision-making.

Convergent Random Forest predictor: methodology for predicting drug response from genome-scale data applied to anti-TNF response.

Biomarker development for prediction of patient response to therapy is one of the goals of molecular profiling of human tissues. Due to the large number of transcripts, relatively limited number of samples, and high variability of data, identification of predictive biomarkers is a challenge for data analysis. Furthermore, many genes may be responsible for drug response differences, but often only a few are sufficient for accurate prediction. Here we present an analysis approach, the Convergent Random Forest (CRF) method, for the identification of highly predictive biomarkers. The aim is to select from genome-wide expression data a small number of non-redundant biomarkers that could be developed into a simple and robust diagnostic tool. Our method combines the Random Forest classifier and gene expression clustering to rank and select a small number of predictive genes. We evaluated the CRF approach by analyzing four different data sets. The first set contains transcript profiles of whole blood from rheumatoid arthritis patients, collected before anti-TNF treatment, and their subsequent response to the therapy. In this set, CRF identified 8 transcripts predicting response to therapy with 89% accuracy. We also applied the CRF to the analysis of three previously published expression data sets. For all sets, we have compared the CRF and recursive support vector machines (RSVM) approaches to feature selection and classification. In all cases the CRF selects much smaller number of features, five to eight genes, while achieving similar or better performance on both training and independent testing sets of data. For both methods performance estimates using cross-validation is similar to performance on independent samples. The method has been implemented in R and is available from the authors upon request: Jadwiga.Bienkowska@biogenidec.com.

Blood autoantibody and cytokine profiles predict response to anti-tumor necrosis factor therapy in rheumatoid arthritis.

INTRODUCTION: Anti-TNF therapies have revolutionized the treatment of rheumatoid arthritis (RA), a common systemic autoimmune disease involving destruction of the synovial joints. However, in the practice of rheumatology approximately one-third of patients demonstrate no clinical improvement in response to treatment with anti-TNF therapies, while another third demonstrate a partial response, and one-third an excellent and sustained response. Since no clinical or laboratory tests are available to predict response to anti-TNF therapies, great need exists for predictive biomarkers. METHODS: Here we present a multi-step proteomics approach using arthritis antigen arrays, a multiplex cytokine assay, and conventional ELISA, with the objective to identify a biomarker signature in three ethnically diverse cohorts of RA patients treated with the anti-TNF therapy etanercept. RESULTS: We identified a 24-biomarker signature that enabled prediction of a positive clinical response to etanercept in all three cohorts (positive predictive values 58 to 72%; negative predictive values 63 to 78%). CONCLUSIONS: We identified a multi-parameter protein biomarker that enables pretreatment classification and prediction of etanercept responders, and tested this biomarker using three independent cohorts of RA patients. Although further validation in prospective and larger cohorts is needed, our observations demonstrate that multiplex characterization of autoantibodies and cytokines provides clinical utility for predicting response to the anti-TNF therapy etanercept in RA patients.

Circulating immune complexes contain citrullinated fibrinogen in rheumatoid arthritis.

INTRODUCTION: There is increasing evidence that autoantibodies and immune complexes (ICs) contribute to synovitis in rheumatoid arthritis (RA), yet the autoantigens incorporated in ICs in RA remain incompletely characterised. METHODS: We used the C1q protein to capture ICs from plasma derived from human RA and control patients. Antibodies specific for immunoglobulin were used to detect ICs, and fibrinogen antibodies were used to detect fibrinogen-containing ICs. RA and control plasma were separated by liquid chromatography, and fractions then characterised by ELISA, immunoblotting and mass spectrometry. Immunohistochemical staining was performed on rheumatoid synovial tissue. RESULTS: C1q-immunoassays demonstrated increased levels of IgG (p = 0.01) and IgM (p = 0.0002) ICs in plasma derived from RA patients possessing anti-cyclic citrullinated peptide (CCP+) autoantibodies as compared with healthy controls. About one-half of the anti-CCP+ RA possessed circulating ICs containing fibrinogen (p = 0.0004). Fractionation of whole RA plasma revealed citrullinated fibrinogen in the high molecular weight fractions that contained ICs. Positive correlations were observed between fibrinogen-containing ICs and anti-citrullinated fibrinogen autoantibodies, anti-CCP antibody, rheumatoid factor and certain clinical characteristics. Immunohistochemical staining demonstrated co-localisation of fibrinogen, immunoglobulin and complement component C3 in RA pannus tissue. Mass spectrometry analysis of immune complexes immunoprecipitated from RA pannus tissue lysates demonstrated the presence of citrullinated fibrinogen. CONCLUSION: Circulating ICs containing citrullinated fibrinogen are present in one-half of anti-CCP+ RA patients, and these ICs co-localise with C3 in the rheumatoid synovium suggesting that they contribute to synovitis in a subset of RA patients.

Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis.

The prediction of response (or non-response) to anti-TNF treatment for rheumatoid arthritis (RA) is a pressing clinical problem. We conducted a genome-wide association study using the Illumina HapMap300 SNP chip on 89 RA patients prospectively followed after beginning anti-TNF therapy as part of Autoimmune Biomarkers Collaborative Network (ABCoN [Autoimmune Bio-markers Collaborative Network]) patient cohort. Response to therapy was determined by the change in Disease Activity Score (DAS28) observed after 14 wks. We used a two-part analysis that treated the change in DAS28 as a continuous trait and then incorporated it into a dichotomous trait of "good responder" and "nonresponder" by European League Against Rheumatism (EULAR) criteria. We corrected for multiple tests by permutation, and adjusted for potential population stratification using EIGENSTRAT. Multiple single nucleotide polymorphism (SNP) markers showed significant associations near or within loci including: the v-maf musculoaponeurotic fibrosarcoma oncogene homolog B (MAFB) gene on chromosome 20; the type I interferon gene IFNk on chromosome 9; and in a locus on chromosome 7 that includes the paraoxonase I (PON1) gene. An SNP in the IL10 promoter (rs1800896) that was previously reported as associated with anti-TNF response was weakly associated with response in this cohort. Replications of these results in independent and larger data sets clearly are required. We provide a reference list of candidate SNPs (P < 0.01) that can be investigated in future pharmacogenomic studies.

STAT4 and the risk of rheumatoid arthritis and systemic lupus erythematosus.

BACKGROUND: Rheumatoid arthritis is a chronic inflammatory disease with a substantial genetic component. Susceptibility to disease has been linked with a region on chromosome 2q. METHODS: We tested single-nucleotide polymorphisms (SNPs) in and around 13 candidate genes within the previously linked chromosome 2q region for association with rheumatoid arthritis. We then performed fine mapping of the STAT1-STAT4 region in a total of 1620 case patients with established rheumatoid arthritis and 2635 controls, all from North America. Implicated SNPs were further tested in an independent case-control series of 1529 patients with early rheumatoid arthritis and 881 controls, all from Sweden, and in a total of 1039 case patients and 1248 controls from three series of patients with systemic lupus erythematosus. RESULTS: A SNP haplotype in the third intron of STAT4 was associated with susceptibility to both rheumatoid arthritis and systemic lupus erythematosus. The minor alleles of the haplotype-defining SNPs were present in 27% of chromosomes of patients with established rheumatoid arthritis, as compared with 22% of those of controls (for the SNP rs7574865, P=2.81x10(-7); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.32). The association was replicated in Swedish patients with recent-onset rheumatoid arthritis (P=0.02) and matched controls. The haplotype marked by rs7574865 was strongly associated with lupus, being present on 31% of chromosomes of case patients and 22% of those of controls (P=1.87x10(-9); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.55). Homozygosity of the risk allele, as compared with absence of the allele, was associated with a more than doubled risk for lupus and a 60% increased risk for rheumatoid arthritis. CONCLUSIONS: A haplotype of STAT4 is associated with increased risk for both rheumatoid arthritis and systemic lupus erythematosus, suggesting a shared pathway for these illnesses.

Interaction between smoking, the shared epitope, and anti-cyclic citrullinated peptide: a mixed picture in three large North American rheumatoid arthritis cohorts.

OBJECTIVE: Recently, Swedish members of the Epidemiological Investigation of Rheumatoid Arthritis (EIRA) provided evidence that smoking may trigger RA-specific immune reactions to citrullinated protein in carriers of HLA-DR shared epitope alleles. In an effort to confirm this interaction between smoking and shared epitope alleles, we performed a case-only analysis of 3 North American RA cohorts. METHODS: A total of 2,476 white patients with RA were studied, 1,105 from the North American Rheumatoid Arthritis Consortium (NARAC) family collection, 753 from the National Inception Cohort of Rheumatoid Arthritis Patients (Inception Cohort), and 618 from the Study of New Onset Rheumatoid Arthritis (SONORA). All patients were HLA-DRB1 typed, and tested for anti-cyclic citrullinated peptide (anti-CCP) and rheumatoid factor. Information about smoking history was obtained by questionnaire. RESULTS: A significant association was found between smoking and the presence of anti-CCP in the NARAC and the Inception Cohort, but not in the SONORA. The shared epitope alleles consistently correlated with anti-CCP in all 3 populations. Using multiple logistic regression analyses, shared epitope alleles were still the most significant risk factor for anti-CCP positivity. Weak evidence of gene-environment interaction between smoking and shared epitope alleles for anti-CCP formation was found only in the NARAC. CONCLUSION: Unlike the EIRA data, we could not confirm a major gene-environment interaction for anti-CCP formation between shared epitope alleles and smoking in 3 North American RA cohorts. Our data indicate a need for further studies to address the full range of environmental factors other than smoking that may be associated with citrullination and RA.

A comparison of founder-only and all-pedigree-members genotype-expression association by regression analysis.

Genotype-expression association analysis using linear regression may produce different test results depending on whether founders only or all pedigreed members are used. This difference is not due to the correlation of samples within a pedigree, because linear mixed models have been applied to account for that correlation. We investigated the possibility that the difference is due to a dependence of expression levels on, among other things, the generation number in the pedigree. Indeed, of the 30 or so studied expression quantitative traits, several of them show significant dependence on the generation number. We propose to use all pedigree members in genotype-expression association analyses whenever the complete genotyping information is available.

Elevated serum levels of interferon-regulated chemokines are biomarkers for active human systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a serious systemic autoimmune disorder that affects multiple organ systems and is characterized by unpredictable flares of disease. Recent evidence indicates a role for type I interferon (IFN) in SLE pathogenesis; however, the downstream effects of IFN pathway activation are not well understood. Here we test the hypothesis that type I IFN-regulated proteins are present in the serum of SLE patients and correlate with disease activity. METHODS AND FINDINGS: We performed a comprehensive survey of the serologic proteome in human SLE and identified dysregulated levels of 30 cytokines, chemokines, growth factors, and soluble receptors. Particularly striking was the highly coordinated up-regulation of 12 inflammatory and/or homeostatic chemokines, molecules that direct the movement of leukocytes in the body. Most of the identified chemokines were inducible by type I IFN, and their levels correlated strongly with clinical and laboratory measures of disease activity. CONCLUSIONS: These data suggest that severely disrupted chemokine gradients may contribute to the systemic autoimmunity observed in human SLE. Furthermore, the levels of serum chemokines may serve as convenient biomarkers for disease activity in lupus.

PTPN22: setting thresholds for autoimmunity.

The 620W allelic variant of the intracellular tyrosine phosphatase, PTPN22, is associated with a number of different autoimmune disorders, and this provides direct evidence for common mechanisms underlying many of these diseases. The associated allele appears to influence thresholds for T cell receptor signaling, and a variety of disease models involving both central and peripheral tolerance can be proposed. However, given the fact that PTPN22 is expressed in a variety of immunologically relevant cell types, the precise mechanisms for these associations remain unclear. In general, the PTPN22 620W allele appears to play a role in autoimmune disorders that have a prominent humoral component, suggesting that further investigation of PTPN22 activity in B cells will be useful. From a genetic perspective, the data highlights the genetic heterogeneity underlying autoimmunity in different ethnic groups.

Gene expression profiling in human autoimmunity.

Human autoimmune diseases are well suited for the application of gene expression profiling. Sampling of blood cells and target tissues has already revealed many important pathways contributing to this spectrum of disorders, and many commonalities are emerging. For instance, clinically distinct diseases such as systemic lupus erythematosus, Sjögren's syndrome, dermatomyositis, and psoriasis all show evidence for dysregulation of the type I interferon pathway. These data suggest that autoimmune diseases will eventually be categorized at the level of gene expression. This work has led to advances in our understanding of disease pathogenesis and in the future promises to facilitate assessments of disease activity and improve targeting of therapies. Here, we review the literature on gene profiling in human autoimmune diseases and provide perspective on the current state of the art.

Overlapping probabilities of top ranking gene lists, hypergeometric distribution, and stringency of gene selection criterion.

When the same set of genes appear in two top ranking gene lists in two different studies, it is often of interest to estimate the probability for this being a chance event. This overlapping probability is well known to follow the hypergeometric distribution. Usually, the lengths of top-ranking gene lists are assumed to be fixed, by using a pre-set criterion on, e.g., p-value for the t-test. We investigate how overlapping probability changes with the gene selection criterion, or simply, with the length of the top-ranking gene lists. It is concluded that overlapping probability is indeed a function of the gene list length, and its statistical significance should be quoted in the context of gene selection criterion.

Regulation of anti-cyclic citrullinated peptide antibodies in rheumatoid arthritis: contrasting effects of HLA-DR3 and the shared epitope alleles.

OBJECTIVE: To examine the association between HLA-DRB1 alleles and the production of anti-cyclic citrullinated peptide (anti-CCP) and rheumatoid factor (RF) autoantibodies in patients with rheumatoid arthritis (RA). METHODS: We studied 1,723 Caucasian RA patients enrolled in the North American Rheumatoid Arthritis Consortium (NARAC) family cohort and the Study of New Onset Rheumatoid Arthritis (SONORA) cohort. All patients were tested for anti-CCP antibodies (by enzyme-linked immunosorbent assay), RF (by nephelometry), and HLA-DR genotype (by polymerase chain reaction and sequence-specific oligonucleotide hybridization). RESULTS: When controlled for the presence of RF, anti-CCP positivity was strongly associated with the HLA-DRB1 shared epitope (SE). In RF+ patients, the presence of the SE was very significantly associated with anti-CCP positivity, with an odds ratio (OR) of 5.8 and a 95% confidence interval (95% CI) of 4.1-8.3. This relationship was also seen in RF- patients (OR 3.1 [95% CI 1.8-5.3]). In contrast, RF positivity was not significantly associated with presence of the SE independently of anti-CCP antibodies. Strikingly, HLA-DRB1*03 was strongly associated with reduced anti-CCP titers, even after controlling for the presence of the SE and restricting the analysis to anti-CCP+ patients. HLA-DR3 was also associated with anti-CCP- RA in our population. CONCLUSION: The HLA-DRB1 SE is strongly associated with the production of anti-CCP antibodies, but not RF. In contrast, HLA-DR3 alleles are associated with anti-CCP- disease and with lower levels of anti-CCP antibodies, even when controlling for the SE. These data emphasize the complexity of the genetic effects of the major histocompatibility complex on the RA phenotype.