Transcriptome Analysis Describing New Immunity and Defense Genes in Peripheral Blood Mononuclear Cells of Rheumatoid Arthritis Patients

Abstract

Background

Large-scale gene expression profiling of peripheral blood mononuclear cells from Rheumatoid Arthritis (RA) patients could provide a molecular description that reflects the contribution of diverse cellular responses associated with this disease. The aim of our study was to identify peripheral blood gene expression profiles for RA patients, using Illumina technology, to gain insights into RA molecular mechanisms.

Methodology/Principal Findings

The Illumina Human-6v2 Expression BeadChips were used for a complete genome-wide transcript profiling of peripheral blood mononuclear cells (PBMCs) from 18 RA patients and 15 controls. Differential analysis per gene was performed with one-way analysis of variance (ANOVA) and P values were adjusted to control the False Discovery Rate (FDR<5%). Genes differentially expressed at significant level between patients and controls were analyzed using Gene Ontology (GO) in the PANTHER database to identify biological processes. A differentially expression of 339 Reference Sequence genes (238 down-regulated and 101 up-regulated) between the two groups was observed. We identified a remarkably elevated expression of a spectrum of genes involved in Immunity and Defense in PBMCs of RA patients compared to controls. This result is confirmed by GO analysis, suggesting that these genes could be activated systemically in RA. No significant down-regulated ontology groups were found. Microarray data were validated by real time PCR in a set of nine genes showing a high degree of correlation.

Conclusions/Significance

Our study highlighted several new genes that could contribute in the identification of innovative clinical biomarkers for diagnostic procedures and therapeutic interventions.

About the Authors:

Vitor Hugo Teixeira

Contributed equally to this work with: Vitor Hugo Teixeira, Robert Olaso

* E-mail: [email protected]

Affiliations GenHotel-EA3886, Evry University - Paris 7 University Medical School, AutoCure European Consortium member, Evry, France, Faculty of Medicine, University of Coimbra, Coimbra, Portugal

Robert Olaso

Contributed equally to this work with: Vitor Hugo Teixeira, Robert Olaso

Affiliation: CEA - IG – CNG, Functional genomic platform, Department of Translational Research, Evry, France

Marie-Laure Martin-Magniette

Affiliations UMR - AgroParisTech/INRA, MIA 518, Paris, France, URGV- UMR - INRA 1165, CNRS 8114, Evry University, Evry, France

Sandra Lasbleiz

Affiliation: Lariboisière Hospital, APHP, Paris, France

Laurent Jacq

Affiliation: GenHotel-EA3886, Evry University - Paris 7 University Medical School, AutoCure European Consortium member, Evry, France

Catarina Resende Oliveira

Affiliation: Center for Neurosciences and Cell Biology, Faculty of Medicine, University of Coimbra, Coimbra, Portugal

Pascal Hilliquin

Affiliation: Sud Francilien Hospital Center, Corbeil-Essonnes, France

Ivo Gut

Affiliation: CEA - IG - CNG - Department of Translational Research, Evry, France

François Cornelis

Affiliations GenHotel-EA3886, Evry University - Paris 7 University Medical School, AutoCure European Consortium member, Evry, France, Lariboisière Hospital, APHP, Paris, France, Sud Francilien Hospital Center, Corbeil-Essonnes, France, Clinical Genetics Unit, Lariboisière Hospital, APHP, Paris, France

Elisabeth Petit-Teixeira

Affiliation: GenHotel-EA3886, Evry University - Paris 7 University Medical School, AutoCure European Consortium member, Evry, France

Introduction

Rheumatoid Arthritis (RA) is an autoimmune disease characterized by chronic and persistent joint synovial tissue inflammation associated with the destruction of affected joints [1]. The multifactorial nature of the disease provides a high RA heterogeneity with specific combinations of a genetic background and environmental factors that influence the susceptibility, severity and outcome of the disease [2]. The RA heterogeneity is demonstrated by the presence of distinct autoantibody specificities, like rheumatoid factor (RF) and anti-cyclic citrullinated peptide antibodies (ACPA) in the serum [3], [4], the differential responsiveness to treatment [5], [6], and the variability in clinical presentation. In addition, several gene expression profiling studies of synovial tissues and peripheral blood mononuclear cells (PBMCs) from RA patients showed marked variation in gene expression profiles that allowed to identify distinct molecular disease mechanisms involved in RA pathology [7], [8]. The relative contribution of the different mechanisms may vary among patients and in different stages of disease. Thus, the broad goals of expression profiling in RA are to (i) improve our understanding of the pathogenic mechanisms underlying RA, (ii) identify new drugs targets, (iii) assess activity of the disease, (iv) predict future outcomes, such as responsiveness therapy, overall disease severity, and organ specific risk and (v) develop new diagnostic tests [9].

PBMCs gene expression profiling allows both pathogenetic and pathophysiological processes identification as demonstrated in several types of diseases: cancer [10], asthma [11], systemic lupus erythematosus (SLE) [12], cardiovascular diseases [13] and psychiatric disorders [14]. Pathogenetic processes are primarily associated with the cause of a disease. Then, microarrays could lead to the identification of abnormal genes and gene activities that may not be only limited to PBMCs, but could occur in cells of pathological tissue as well. In contrast, pathophysiological changes in lymphocytic gene expression are considered an essentially normal reaction of the immune system to a pathological stimulus. Therefore, pathophysiological gene profiles may be shared in a variety of diseases, whereas pathogenetic gene expression is expected to be disease specific [15]. The differences in expression profiles provide opportunities to stratify RA patients based on molecular criteria that may require different treatment strategies.

Considering several comparative studies, Illumina and other microrray technologies have similar performances [16]–[18]. However, these studies showed that each approach was able to detect specific genes, meaning an increase in knowledge by each platform. To complete previous studies on RA with other microarrays [7], we applied Illumina large-scale gene expression profiling in PBMCs of RA patients to potentially gain insights into molecular mechanism of this disease. We identified new genes involved in different functional Immunity and Defense related mechanisms as pro-inflammation, anti-microbial activity, cellular stress and immunomodulatory functions in Rheumatoid Arthritis.

Materials and Methods

Study population

The study and all protocols presented here were approved by the Ethics Committees of Bicêtre and Saint Louis Hospitals (Paris, France) and all study participants provided written informed consent. All RA patients satisfied the revised criteria of the American College of Rheumatology [19] according to the rheumatologist in charge of the patient. A rheumatologist university fellow reviewed all clinical data. Characteristics of the 18 RA French Caucasian Patients are reported in Table 1. Among the control group consisted of 15 RA French Caucasian healthy individuals, 11 were females (mean ± Standard Deviation (SD) age at enrolment 56.9 ± 6.6). In all comparisons mentioned, the groups were age and sex-matched.

[Figure omitted. See PDF.]

Table 1. Clinical and demographic characteristics of the RA patients.

https://doi.org/10.1371/journal.pone.0006803.t001

Isolation of total RNA

Peripheral blood (PB) was drawn in PAXgene RNA isolation tubes (PreAnalytix) from 18 patients and 15 controls. Total RNA was isolated from PBMCs using the PAXgene RNA isolation kit (PreAnalytix). Total RNA yield (ng) was determined spectrophotometrically using the NanoDrop ND-1000 (Wilmington). Total RNA profiles were recorded using a Bioanalyzer 2100 (Agilent). RNA integrity numbers were determined and the mean value was 8.07+/−0.51 SD and a Coefficient of Variation (CV) of 6.4%.

Probe synthesis, hybridization and detection

cRNA was synthesized, amplified and purified using the Illumina TotalPrep RNA Amplification Kit (Ambion Inc.) following manufacturer recommendations. Briefly, 200 ng of RNA was reverse transcribed. After second strand synthesis, the cDNA was transcribed in vitro and cRNA labelled with biotin-16-UTP. Labelled probe hybridization to Illumina BeadChips human-6v2 was carried out using Illumina's BeadChip 6v2 protocol. These beadchips contain 48,701 unique 50-mer oligonucleotides in total, with hybridization to each probe assessed at ∼30 different beads on average. 22,403 probes (46%) are targeted at Reference Sequence (RefSeq) [20] transcripts and the remaining 26,298 probes (54%) are for other transcripts, generally less well characterized (including predicted transcripts).

Beadchips were scanned on the Illumina BeadArray 500GX Reader using Illumina BeadScan image data acquisition software (version 2.3.0.13). Illumina BeadStudio software (version 1.5.0.34) was used for preliminary data analysis. To assess quality metrics of each run, several quality control procedures were implemented. Total RNA control samples were analyzed with each run. The Illumina BeadStudio software was used to view control summary reports, scatter plots of the total RNA control results from different days and scatter plots of daily run samples. The scatter plots compared control against control or sample against sample and calculated a correlation coefficient (Figure 1). Viewing the scatter plots determined whether controls across different days varied in quality, indicating a reduction in assay performance, and highlighted those samples that were of lower quality. The control summary report is generated by the BeadStudio software, which evaluates the performance of the built-in controls of the BeadChips across particular runs. This allows the user to look for variations in signal intensity, hybridization signal, background signal and the background to noise ratio for all samples analyzed in that run. Data are expressed as log2 ratios of fluorescence intensities of the experimental and the common reference sample. The Illumina data were then normalized using the ‘normalize quantiles’ function in the BeadStudio Software.

[Figure omitted. See PDF.]

Figure 1. Scatter plots representation of signal intensities.

Scatter plot of technical replicates: (A) same reference RNA undergoing two different hybridization or (B) same reference RNA from two different labeling runs. Typical scatter plots of data obtained from two patients (C) or a patient and a control (D). Pearson correlation coefficient is indicated in each scatter plot.

https://doi.org/10.1371/journal.pone.0006803.g001

All microarray data reported in this study is described in accordance with MIAME guidelines and have been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) public repository, and they are accessible through GEO accession (GSE15573).

Real-Time PCR

Total RNA was reverse transcribed using Superscript III and oligo(dT) primers (Invitrogen) according to the manufacturer's instructions. Real-time quantitative PCR was carried out using the SYBR-green master mix (Applied BioSystems) in an Mx 3005P thermocycler (Agilent). PCR conditions were 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec and, 60°C for 1 min. At the end of the amplification reaction, melting curve analyses were performed to confirm the specificity as well as the integrity of the PCR products by the presence of a single peak. Gene-specific primers were designed inside or nearby the microarray sequence targeted, using Primer Express Software (PE Applied Biosystems). Primers sequences of all genes analyzed are provided in Table S1. Absence of cross contamination and primer dimer was checked on genomic DNA and water. From a list of 8 housekeeping genes (HKG), we chose HMBS and ALDOA, which meet the criteria of less variation between samples and compatible expression level with the studied genes. The geometric mean of HKG expression was used to normalize the expression of genes of interest [21]. Standard curves were generated from assays made with serial dilutions of reference cDNA to calculate PCR efficiencies (100%+/−15%, with r2 > = 0.997). Ct samples were transformed into quantity values using the formula (1+Efficiency)Ct. Only means of triplicate with a CV of less than 10% were analyzed. Inter-plate variation was below 8%.

Statistical analysis

Statistical analysis on microarray data was performed using one-way analysis of variance (ANOVA) per gene where the normalized signals are explained by the patient status. One contrast was built to determine an expression difference between controls and RA patients [22]. Since the number of individuals is large, the residual variance was used to calculate the statistic test [23]. The raw P values were adjusted by the Benjamini-Hochberg procedure, which controls the False Discovery Rate (FDR) [24]. For the contrast, a gene is considered differentially expressed if the Benjamini-Hochberg-corrected P value is less than 0.05.

Genes that were expressed at significantly different levels between patients and controls were analyzed by supervised hierarchical clustering (uncentered correlation, complete linkage) [25] to visualize the correlation of co-expressed genes in Treeview (available at http://rana.lbl.gov/EisenSoftware.htm).

For an interpretation of the biological processes that are represented by the genes that show a significantly different level of expression in RA patients compared to the controls, we applied Gene Ontology analysis in the PANTHER database at http://www.pantherdb.org (Applied Biosystems) [26]. PANTHER uses the binomial statistics tool to compare our gene list to a reference list (NCBI: Homo sapiens genes) determining the statistically significant over- or under- representation of PANTHER biological process [27]. After, for each biological process in PANTHER, the genes associated with that term are evaluated according to the likelihood that their fold changes were drawn randomly from the overall distribution of fold changes. The Mann-Whitney U Test (Wilcoxon Rank-Sum Test) is used to determine the P value that, say, if specific biological process genes have random fold changes relative to the overall list of values that was input. A significant P value indicates that the distribution (fold change) for this category is non-random and different from the overall distribution [28]. In both statistical tests, processes with a P value <0.05 were considered significant after Bonferroni correction which was applied to adjust for multiple comparisons. Correlations between two set of data were measured using Pearson coefficient.

Results

Gene expression profiling in PMBCs of RA patients

Genome-wide transcriptional profiles of PBMCs from 18 RA patients and 15 age and sex-matched controls were measured on microarrays that contain 48,701 unique 50-mer oligonucleotides in total, with a mean ∼30 hybridizations per sequence. Data were analyzed using ANOVA. Using this test with a FDR of 5% we identified 380 transcripts with significant expression. The proportion of detected transcripts was substantially higher among RefSeq genes (91%) than non-RefSeq genes (9%), reflecting the greater degree of knowledge and certainty about the existence of RefSeq transcripts. Four genes, represented more than once in this list, were averaged from sequences with the same Unigene identifier. Significant difference in expression level between the two groups was observed for 339 RefSeq genes. Among them, 238 were downregulated (Table S2) and 101 were upregulated (Table S3). The significant gene expression differences between RA patients and controls were visualized in a cluster diagram (Figure 2).

[Figure omitted. See PDF.]

Figure 2. Cluster diagram of the expression of 339 significantly expressed genes in 18 RA patients and 15 controls.

Genes are organized by hierarchical clustering based on overall similarity in expression patterns. Red represents relative expression greater than the median expression level across all samples, and green represents an expression level lower than the median. Black indicates intermediate expression.

https://doi.org/10.1371/journal.pone.0006803.g002

Genes upregulated in RA

To categorize the up-regulated 101 genes into functional biological groups we used the PANTHER database (described in Methods section). We observed an elevated expression of a spectrum of genes involved in Immunity and Defense, nucleoside, nucleotide and nucleic acid metabolism, signal transduction, protein metabolism and modification, mRNA transcription, transport and developmental processes in the peripheral blood of RA patients compared to controls. Compared to a NCBI Homo sapiens reference list, the differentially up-regulated genes list revealed three biological processes significantly over represented (P<0.05) (Figure 3A). Then, the fold change of the genes associated to each biological process was compared to the overall distribution of fold changes. Immunity and Defense was the only significant functional biological process after Bonferroni correction (P = 0.03) (Figure 3B).

[Figure omitted. See PDF.]

Figure 3. Gene ontology analysis of the most representative biological processes represented by the 101 up-regulated genes in RA patients.

Gene Ontology (GO) analysis in the PANTHER database was applied for an interpretation of the biological processes that are represented by the genes showing a higher significantly different expression level in RA patients compared to the controls. A and B - Biological processes with a P value <0.05 (red bars) were considered significant after Bonferroni correction. The GO analysis P-value was plotted on the y axis versus biological processes on the x axis. A - The binomial statistics tool was used to compare our gene list to a reference list (NCBI: Homo sapiens genes) determining the statistically significant (P value) over- or under- representation of PANTHER biological process. B - The Mann-Whitney U Test was used to determine significant P values which indicate that the distribution (fold change) for each biological process is non-random and different from the overall distribution. OP - Oxidative Phosphorylation; ET - Electron Transport; ID - Immunity and Defense; NNNAM - Nucleoside, Nucleotide and Nucleic Acid Metabolism; T – Transport.

https://doi.org/10.1371/journal.pone.0006803.g003

This cluster of genes involved in Immunity and Defense process contains the S100 family proteins S100 calcium-binding protein A8 (S100A8), S1000A9 and S100A12, the orosomucoid family proteins ORM1 and ORM2, as well as other inflammatory mediators like lymphocyte antigen 96 (LY96), cathelicidin antimicrobial peptide (CAMP), thioredoxin (TXN), allograft inflammatory factor 1 (AIF1), nuclear factor of activated T-cells 5 (NFAT5), F5 (coagulation factor V), SLC11A1 (solute carrier family 11, member 1) and PGLYRP1 (peptidoglycan recognition protein 1). These genes are characterized in Table 2.

[Figure omitted. See PDF.]

Table 2. Differentially up-regulated transcripts linked to Immunity and Defense biological process in RA patients.

https://doi.org/10.1371/journal.pone.0006803.t002

Genes downregulated in RA

Compared to the same NCBI Homo sapiens reference list, the mainstream of the genes that showed a lower expression in RA patients are linked to different biological processes such as nucleoside, nucleotide and nucleic acid metabolism, mRNA transcription and regulation, cell cycle, intracellular protein traffic. A small number of down-regulated genes were also involved in oncogenesis like runt-related transcription factor 3 (RUNX3), SMAD family member 7 (SMAD7), PHD finger protein 17 (PHF17) as well as interleukin-1 receptor-associated kinase 2 (IRAK2), interleukin 2 receptor, beta (IL2RB), CD96 and SH2B adaptor protein 3 (SH2B3) related to immune functions. Therefore, in PANTHER database classification we found four significant functional biological processes in our 238 down-regulated genes list. However, no significant biological process was found (after Bonferroni correction) when we evaluate all the down-regulated genes associated to biological processes, according with their fold changes drawn randomly from the overall distribution of fold changes (data not shown).

Real time PCR validation

In all samples, we confirmed the expression of five up-regulated genes (LY96, S100A12, ORM2, ORM1, RPL31) and four down regulated genes (IL2RB, DNMT1, RUNX3 and IRF1) in RA patients, by real-time PCR. From microarray and real time PCR data, we calculated the RA patients/controls ratio for each genes expression. The qPCR expression data of the nine genes showed a high correlation with the microarray expression data (r = 0.937) (Figure 4).

[Figure omitted. See PDF.]

Figure 4. Correlation between microarray and real time PCR data.

The scatter plot compares mean expression of RA patients/controls ratio for nine genes. Each point represents the RA patients/controls ratio from the microarray (y axis) and real time PCR (x axis). Pearson correlation coefficient is indicated in the scatter plot.

https://doi.org/10.1371/journal.pone.0006803.g004

Discussion

Microarray technology has been used to discriminate differences in gene expression profiles in tissues and PBMCs. Both synovial tissue and PBMCs have been used to evaluate differences in the gene expression profiles in RA [7], [8]. If expression-based profiling is to be of practical importance, sample accessibility becomes crucial. In this context, peripheral mononuclear cells are key sentinels of host defence, being used to identify novel disease mediators, disease variants and treatment responses [7], [29], [30].

Transcriptome studies using Illumina and other technologies showed that each approach was able to detect specific genes, meaning an increase in knowledge by each platform [16]–[18]. To complete previous studies on RA with Affymetrix or double colour microarrays, we decided to use Illumina technology. Our study did not confirm a specific expression for the genes regulated by interferon type I, as described in a RA large-scale expression profiling [8] (two colours analysis protocol using Stanford University microarrays) (data not shown). Our analysis revealed only one significantly increased biological mechanism: Immunity and Defense. This process was already highlighted by other studies in RA, as several genes that we described (S100A8, S1000A9, S100A12 and AIF1). Additionally, our study identified new genes like LY96/MD-2, NFAT5, TXN, CAMP/LL37, ORM1, ORM2, SLC11A1, PGLYRP1 and F5. These genes are involved in different functional Immunity and Defense related mechanisms as pro-inflammation, anti-microbial activity, oxidative and osmotic cellular stress and immunomodulatory functions (Figure 5).

[Figure omitted. See PDF.]

Figure 5. Biological functions of Immunity and Defense genes highlighted in our Rheumatoid Arthritis study.

The genes differentially expressed in Immunity and Defense process are stratified in four functional related mechanisms: pro-inflammation, anti-microbial activity, oxidative and osmotic cellular stress and immunomodulatory functions.

https://doi.org/10.1371/journal.pone.0006803.g005

The S100 calcium-binding proteins (S100A8, S100A9, and S100A12) are characterized by strong prevalence in cells of myeloid origin. Activated phagocytes expressing S100A8 and S100A9 proteins are among the first cells infiltrating inflammatory lesions in the synovium [31]. The expression of S100A8 and S100A9 was found to be strongest at the cartilage-pannus junction, which is the prime site of cartilage destruction and bone erosion in arthritis [32]. S100A12 is strongly expressed in inflamed synovial tissue, whereas it is nearly undetectable in synovia of control subjects or patients after successful treatment [33]. Further studies showed increased levels of S100A8/S100A9 and S100A12 concentrations in synovial fluid (SF) and serum in RA patients [34]–[36]. In addition, S100A9, S100A8 and S100A12 levels were associated with body mass index, presence of ACPA and RF and presence of ACPA, respectively [37]. Furthermore, several expression profiling studies in PBMCs showed a highly expression of S100A8, S100A9 and S100A12 in RA patients compared to controls [8], [38], [39]. Therefore, pro-inflammatory S100 proteins are attractive therapeutic targets for immune interventions in the treatment of RA.

Human allograft inflammatory factor-1 (AIF-1) is a Ca2+-binding EF-hand protein encoded within the major histocompatibility complex class (MHC) III region of chromosome 6. AIF-1 is produced by macrophages and lymphocytes, and its synthesis is mediated by several cytokines, such as interferon-γ (IFN-γ) [40]. Kimura et al. have demonstrated that AIF-1 is expressed in synovial and mononuclear cells in RA synovial tissue and increases proliferation of cultured synoviocytes and IL-6 production by these synoviocytes and PBMCs [41]. A recent study, showed an increased expression of AIF1 mRNA in RA patients PMBCs compared with controls as well as in synovial macrophages in the lining layer of all the inflamed RA synovial membranes compared with non-inflamed OA controls [42]. AIF-1 plays therefore an important role in the pathogenesis of RA by affecting key processes such as the activation of synovial cell proliferation and the inflammatory cytokine cascade including IL-6 in joints and may represent a new molecular target in RA therapy.

The LY96 (MD-2) acts as an extracellular adaptor protein in the activation of TLR4 by the lipopolysaccharide (LPS) of Gram negative bacteria, a well-known inducer of the innate immune response [43]. Knockout studies in mice have demonstrated that MD-2 is indispensable for LPS responses [44]. Multiple evidences points to the potential importance of TLR signalling in RA pathogenesis by means of the presence of different TLR ligands and functional TLR receptors in inflamed joints from patients with RA [45]. TLR4 mRNA is highly expressed in the synovium at early stages of RA as well as at later stages of the disease. In vitro stimulation of the RA synovial fibroblasts with TLR4 ligand (LPS) produces a wide range of proinflammatory cytokines, chemokines, and tissue destructive enzymes [46]. Recently, the potential value of TLR4 or signals derived from this receptor as therapeutic targets has become clearer - mainly as a result of studies in animal models of joint inflammation but also in human RA [47]. Interestingly, in our study, mRNA TLR4 was up-regulated in RA patients compared to controls but not significantly. Thus, the TLR4-MD2-LPS complex could be involved in the activation of synovial fibroblasts and contribute to the development of synovial inflammation and joint destruction.

NFAT5, the primordial member of the NFAT family, is expressed by almost all cells and is activated in response to osmotic stress. In lymphocytes, NFAT5 controls the osmotic stress-induced expression of several cytokines, including tumour-necrosis factor (TNF) and lymphotoxin-β. NFAT5-deficient mice have impaired T-cell function under hyperosmotic conditions and decreased cellularity of the thymus and spleen [48]. In RA, NFAT5 mRNA is expressed in proliferating RA-synovial fibroblasts (RA-SF) but not in nonproliferating RA-SF. Furthermore, NFAT5 mRNA is expressed in RA synovium - but not in normal individuals - as well as at sites of bone destruction. NFAT5 could be then related not only with proliferation but also with the activation and invasion of RA-SF in vivo [49].

Oxidative stress to essential cell components caused by oxygen free radicals is generally considered as a serious mechanism in RA pathogenicity [50]. Thioredoxin (TXN), a cellular reducing catalyst induced by oxidative stress, is involved in the stimulation of the DNA-binding activity of NF-kb transcription factor [51]. Increased cytokine production driven by NF-κb can enhance expression of vascular adhesion molecules that attract leucocytes into the joint, as well as matrix metalloproteinases that help to degrade the extracellular matrix [50]. TXN concentrations were found significantly elevated in SF and serum of RA patients. The positive correlation between the SF TXN and the serum C Reactive Protein in the absence of a high concentration of SF TNF-α may indicate that TXN is involved in the prolongation and persistence of the RA inflammation, because high concentrations of TXN could stimulate NF-kB activation in the presence of the otherwise insufficient concentration of TNF-α [52]. Thus, TXN monitoring in RA patients could provide useful information regarding the extent of oxidative stress. Furthermore, truncated thioredoxin (Txn80) stimulates monocytes/macrophages to induce IL-12, implying that it is involved in immune inflammatory reactions directing Th1 immunity and IFN-γ production [53]. Moreover, Kim et al. have demonstrated that human TXN is a novel target gene induced by IFN-γ [54]. RA is a Th1-driven disease which has IFN-γ as characteristic Th1 cytokine and subsequently the TXN mRNA up-regulation in PBMCs could suggest an involvement in RA IFN-γ pathway.

CAMP (LL37) is an antimicrobial peptide which has a broad of antimicrobial activity. LL-37 has the potential to participate in the innate immune response both by killing bacteria and by recruiting a cellular immune response [55]. Gilliet and Lande recently found that LL37, overexpressed in psoriatic skin, is the key mediator of plasmacytoid dendritic cells (pDCs) activation in psoriasis. LL37 converts nonstimulatory self-DNA into a potent trigger of pDCs to produce IFN. pDCs respond to self-DNA if coupled with an antimicrobial peptide, suggesting that modified self-DNA drives autoimmunity in psoriasis by activating TLR9 [56]. In RA, recent studies have described pDCs as perpetuators of synovial inflammation and modulators of B cell responses in the synovial tissue [57]. The mRNA up-regulation of the LL37 gene found in RA patients comparing to controls in our study could be due to an immunological response to infectious agents such as bacteria and viruses. Moreover the LL37 over-expression could be implicated in a pDCs-dependent mechanism involved in the perpetuation of RA inflammation through the abolition of self-tolerance and subsequent emergence of self-reactive lymphocytes.

Human alpha-1-acid glycoprotein (AGP) - also called orosomucoid - is a 37-kDa molecule consisting of a heavily glycosylated single polypeptide chain. Alpha-1-acid glycoprotein 1 (AGP1) and alpha-1-acid glycoprotein 2 (AGP2), coded by ORM1 and ORM2 genes respectively, are positive acute phase proteins [58], [59]. AGPs plasma concentration may increase several fold during acute phase reactions such as inflammation or chronic disease. AGPs have a strong immunomodulatory function [60]. It was shown that under pathological conditions not only the total concentration of AGPs but its glycosylation pattern may be altered [59]. Smith and colleagues have demonstrated that the AGPs populations in the serum and synovial fluid of RA patients are distinct in terms of glycosylation pattern. This discovery has direct functional significance since only the serum AGPs population is capable of blocking leucocyte adhesion [61]. Furthermore, Haston et al. have shown that AGPs can influence MMP-13 activity. It is hypothesized that AGPs may form part of a negative feedback mechanism which is inadequate to prevent disease progression in rheumatoid arthritis. These processes may exacerbate the increased turnover of collagen characteristic of the disease [62]. These results suggest an interesting role for AGP in RA pathogenesis.

SLC11A1 (formerly called NRAMP1) is a gene that is important in macrophage-mediated natural resistance to a variety of intracellular pathogens. Exogenous and endogenous agents that mediate inflammation by activating the macrophage can cause NRAMP1 translocation to the membrane of the phagolysosome, where it serves as a cation transporter. The significant increase in iron deposition observed in the synovial membrane of rheumatoid arthritis patients, and foam cells in atherosclerotic lesions, could be attributable to NRAMP1 [63]. In human RA synovium NRAMP1 was detected in macrophages and neutrophils in the linning and subinitimal zone, as well as in inflammatory infiltrates, but was absent in fibroblasts [64]. NRAMP1 has also pleiotropic effects on macrophage function, including upregulation of chemokine/cytokine gene, TNFα, IL-1β, inducible nitric oxide synthase (iNOS), MHC expression as well as tumoricidal and antimicrobial activity. These effects are involved in resistance to infection and may also be involved in induction and maintenance of autoimmune disease [65]. TNFα and IL-1β play important roles in inflammation and tissue destruction of RA [66], [67]. Bacterial or viral infection may play a role in triggering the development of RA and TNFα and iNOS are key players in enhanced antimicrobial activity of activated macrophages [65].

PGRPs are innate immunity proteins, recognizing bacterial peptidoglycan, and acting in antibacterial immunity. PGRP-1 seems to belong to the innate immune arm of effectors molecules, such as antimicrobial peptides and C-type lectins, among others. This protein was shown to be almost exclusively present as a soluble protein in the granules of polymorphonuclear leucocytes (PMN) [68]. Saha et al. examined the immunomodulating activities of the PGRPs in a peptidoglycan-induced arthritis mice model. They showed that a systemic injection of peptidoglycan (PG) or muramyldipeptide (MDP) induces an acute arthritis of the joints of the feet in BalbC mice. Peptidoglycan-induced arthritis PGLYRP-1−/− mice, had a MDP-induced activation of proinflammatory genes than WT mice. Moreover, PGLYRP-1−/− mice have longer-lasting MDP-induced arthritis than WT mice. The anti-inflammatory function of PGLYRP-1 manifests itself only in the later stages of MDP-induced arthritis, which is consistent with the local release of PGLYRP-1 from PMN granules after PMNs' arrival into the mice foot [69]. These data point to that PGRP-1 could have a specialized but nevertheless significant role in signalling events like arthritis in mammals.

The multi-step coagulation complex system is activated by tissue factor (TF), which is exposed to blood. In this process, factor 5 (FV) is cleaved and activates factor Va (FVa). After several proteins interactions, the prothrombinase complex (FXa–FVa) converts prothrombin to thrombin which is generated in a large amount. Thrombin activates FV and FVIII and platelets converting fibrinogen to a fibrin clot [70]. Accumulation of fibrin in the RA synovium exceeds that in control tissue by a wide margin and represents one of the most striking pathologic features of rheumatoid synovitis. For some time, this fibrin deposition has been considered to be a serious contributor to permanent damage by maintaining a vicious circle of inflammation [71], [72]. In the extravascular coagulation at the arthritic synovial joint sequential activity of factor Xa (in the presence of cofactor Va) and of thrombin leads to fibrin deposition in the joint [73]. Thus, mRNA F5 up-regulation in RA patients could enhance the production of factor V that under cleavage produces factor Va. This Va increase could lead to the augmentation of thrombin and consequently fibrin in RA joints.

In conclusion, our study highlighted several new genes (LY96/MD-2, NFAT5, TXN, CAMP/LL37, ORM1, ORM2, SLC11A1, PGLYRP1 and F5) in PBMCs of RA patients that could contribute in the identification of innovative clinical biomarkers for diagnostic procedures and therapeutic interventions. Nevertheless, comparative analysis with another disease involving an inflammatory process could clarify the relation between the expression profiling and the pathophysiological processes specifically involved in Rheumatoid Arthritis.z

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Primers sequences of the nine genes analyzed by real-time PCR

https://doi.org/10.1371/journal.pone.0006803.s001

(0.04 MB DOC)

Table S2.

List of 238 downregulated genes differentially expressed between RA patients and controls

https://doi.org/10.1371/journal.pone.0006803.s002

(0.28 MB DOC)

Table S3.

List of 101 upregulated genes differentially expressed between RA patients and controls

https://doi.org/10.1371/journal.pone.0006803.s003

(0.14 MB DOC)

Acknowledgments

We thank Jeanne-Antide Perrier and Magali de Surville for microarray and real time PCR technical assistance. We are grateful to the RA patients and individual controls for their participation in this study.

Author Contributions

Conceived and designed the experiments: VHT EPT. Performed the experiments: VHT RO. Analyzed the data: VHT RO MLMM EPT. Contributed reagents/materials/analysis tools: VHT RO SL LJ PH IG FC EPT. Wrote the paper: VHT RO MLMM CRO IG EPT.

Citation: Teixeira VH, Olaso R, Martin-Magniette M-L, Lasbleiz S, Jacq L, Oliveira CR, et al. (2009) Transcriptome Analysis Describing New Immunity and Defense Genes in Peripheral Blood Mononuclear Cells of Rheumatoid Arthritis Patients. PLoS ONE4(8): e6803. https://doi.org/10.1371/journal.pone.0006803

References

1. Seldin MF, Amos CI, Ward R, Gregersen PK (1999) The genetics revolution and the assault on rheumatoid arthritis (review). Arthritis Rheum 42: 1071–1079.MF SeldinCI AmosR. WardPK Gregersen1999The genetics revolution and the assault on rheumatoid arthritis (review).Arthritis Rheum4210711079

2. Firestein GS (2003) Evolving concepts of rheumatoid arthritis [review]. Nature 423: 356–361.GS Firestein2003Evolving concepts of rheumatoid arthritis [review].Nature423356361

3. Zendman AJ, van Venrooij WJ, Pruijn GJ (2006) Use and significance of anti-CCP autoantibodies in rheumatoid arthritis. Rheumatology (Oxford) 45: 20–25.AJ ZendmanWJ van VenrooijGJ Pruijn2006Use and significance of anti-CCP autoantibodies in rheumatoid arthritis.Rheumatology (Oxford)452025

4. Hueber W, Kidd BA, Tomooka BH, Lee BJ, Bruce B, et al. (2005) Antigen microarray profiling of autoantibodies in rheumatoid arthritis. Arthritis Rheum 52: 2645–2655.W. HueberBA KiddBH TomookaBJ LeeB. Bruce2005Antigen microarray profiling of autoantibodies in rheumatoid arthritis.Arthritis Rheum5226452655

5. Lipsky PE, van der Heijde DM, St Clair EW, Furst DE, Breedveld FC, et al. (2000) Infliximab and methotrexate in the treatment of rheumatoid arthritis. Anti-Tumor Necrosis Factor Trial in Rheumatoid Arthritis with Concomitant Therapy Study Group. N Engl J Med 343: 1594–1602.PE LipskyDM van der HeijdeEW St ClairDE FurstFC Breedveld2000Infliximab and methotrexate in the treatment of rheumatoid arthritis. Anti-Tumor Necrosis Factor Trial in Rheumatoid Arthritis with Concomitant Therapy Study Group.N Engl J Med34315941602

6. Edwards JC, Szczepanski L, Szechinski J, Filipowicz-Sosnowska A, Emery P, et al. (2004) Efficacy of B-cell-targeted therapy with rituximab in patients with rheumatoid arthritis. N Engl J Med 350: 2572–2581.JC EdwardsL. SzczepanskiJ. SzechinskiA. Filipowicz-SosnowskaP. Emery2004Efficacy of B-cell-targeted therapy with rituximab in patients with rheumatoid arthritis.N Engl J Med35025722581

7. Toonen EJ, Barrera P, Radstake TR, van Riel PL, Scheffer H, et al. (2008) Gene expression profiling in rheumatoid arthritis: current concepts and future directions. Ann Rheum Dis 67(12): 1663–1669.EJ ToonenP. BarreraTR RadstakePL van RielH. Scheffer2008Gene expression profiling in rheumatoid arthritis: current concepts and future directions.Ann Rheum Dis67(12)16631669

8. van der Pouw Kraan TC, Wijbrandts CA, van Baarsen LG, Voskuyl AE, Rustenburg F, et al. (2007) Rheumatoid arthritis subtypes identified by genomic profiling of peripheral blood cells: assignment of a type I interferon signature in a subpopulation of patients. Ann Rheum Dis 66(8): 1008–1014.TC van der Pouw KraanCA WijbrandtsLG van BaarsenAE VoskuylF. Rustenburg2007Rheumatoid arthritis subtypes identified by genomic profiling of peripheral blood cells: assignment of a type I interferon signature in a subpopulation of patients.Ann Rheum Dis66(8)10081014

9. Baechler EC, Batliwalla FM, Reed AM, Peterson EJ, Gaffney PM, et al. (2006) Gene expression profiling in human autoimmunity. Immunol Rev 210: 120–137.EC BaechlerFM BatliwallaAM ReedEJ PetersonPM Gaffney2006Gene expression profiling in human autoimmunity.Immunol Rev210120137

10. Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, et al. (2000) Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature 403(6769): 503–511.AA AlizadehMB EisenRE DavisC. MaIS Lossos2000Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling.Nature403(6769)503511

11. Brutsche MH, Brutsche IC, Wood P, Brass A, Morrison N, et al. (2001) Apoptosis signals in atopy and asthma measured with cDNA arrays. Clin Exp Immunol 123(2): 181–187.MH BrutscheIC BrutscheP. WoodA. BrassN. Morrison2001Apoptosis signals in atopy and asthma measured with cDNA arrays.Clin Exp Immunol123(2)181187

12. Mandel M, Achiron A (2006) Gene expression studies in systemic lupus erythematosus. Lupus 15(7): 451–456.M. MandelA. Achiron2006Gene expression studies in systemic lupus erythematosus.Lupus15(7)451456

13. Henriksen PA, Kotelevtsev Y (2002) Application of gene expression profiling to cardiovascular disease. Cardiovasc Res 54(1): 16–24.PA HenriksenY. Kotelevtsev2002Application of gene expression profiling to cardiovascular disease.Cardiovasc Res54(1)1624

14. Colangelo V, Schurr J, Ball MJ, Pelaez RP, Bazan NG, et al. (2002) Gene expression profiling of 12633 genes in Alzheimer hippocampal CA1: transcription and neurotrophic factor down-regulation and up-regulation of apoptotic and pro-inflammatory signaling. J Neurosci Res 70(3): 462–473.V. ColangeloJ. SchurrMJ BallRP PelaezNG Bazan2002Gene expression profiling of 12633 genes in Alzheimer hippocampal CA1: transcription and neurotrophic factor down-regulation and up-regulation of apoptotic and pro-inflammatory signaling.J Neurosci Res70(3)462473

15. Gladkevich A, Nelemans SA, Kauffman HF, Korf J (2005) Microarray profiling of lymphocytes in internal diseases with an altered immune response: potential and methodology. Mediators Inflamm 14: 317–330.A. GladkevichSA NelemansHF KauffmanJ. Korf2005Microarray profiling of lymphocytes in internal diseases with an altered immune response: potential and methodology.Mediators Inflamm14317330

16. MAQC Consortium (2006) The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol 24(9): 1151–161.MAQC Consortium2006The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements.Nat Biotechnol24(9)1151161

17. Barnes M, Freudenberg J, Thompson S, Aronow B, Pavlidis P (2005) Experimental comparison and cross-validation of the Affymetrix and Illumina gene expression analysis platforms. Nucleic Acids Res 33(18): 5914–5923.M. BarnesJ. FreudenbergS. ThompsonB. AronowP. Pavlidis2005Experimental comparison and cross-validation of the Affymetrix and Illumina gene expression analysis platforms.Nucleic Acids Res33(18)59145923

18. Maouche S, Poirier O, Godefroy T, Olaso R, Gut I, et al. (2008) Performance comparison of two microarray platforms to assess differential gene expression in human monocyte and macrophage cells. BMC Genomics 9: 302.S. MaoucheO. PoirierT. GodefroyR. OlasoI. Gut2008Performance comparison of two microarray platforms to assess differential gene expression in human monocyte and macrophage cells.BMC Genomics9302

19. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, et al. (1988) The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 31: 315–324.FC ArnettSM EdworthyDA BlochDJ McShaneJF Fries1988The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis.Arthritis Rheum31315324

20. Pruitt KD, Tatusova T, Maglott DR (2007) NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res 35: D61–65.KD PruittT. TatusovaDR Maglott2007NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins.Nucleic Acids Res35D6165

21. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3(7): J. VandesompeleK. De PreterF. PattynB. PoppeN. Van Roy2002Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.Genome Biol3(7)

22. Gentleman R, Carey V (2002) Bioconductor. R. GentlemanV. Carey2002Bioconductor.RNews, 2, 1116. RNews, 2, 1116.

23. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, et al. (2003) Exploration, Normalization, and Summaries of High Density Oligonucleotide Array Probe Level Data. Biostatistics 4(2): 249–264.RA IrizarryB. HobbsF. CollinYD Beazer-BarclayKJ Antonellis2003Exploration, Normalization, and Summaries of High Density Oligonucleotide Array Probe Level Data.Biostatistics4(2)249264

24. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing, J R Statistics Soc 57: 289–300.Y. BenjaminiY. Hochberg1995Controlling the false discovery rate: a practical and powerful approach to multiple testing,J R Statistics Soc57289300

25. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 95: 14863–14868.MB EisenPT SpellmanPO BrownD. Botstein1998Cluster analysis and display of genome-wide expression patterns.Proc Natl Acad Sci U S A951486314868

26. Mi H, Lazareva-Ulitsky B, Loo R, Kejariwal A, Vandergriff J, et al. (2005) The PANTHER database of protein families, subfamilies, functions and pathways. Nucleic Acids Res 33: D284–288.H. MiB. Lazareva-UlitskyR. LooA. KejariwalJ. Vandergriff2005The PANTHER database of protein families, subfamilies, functions and pathways.Nucleic Acids Res33D284288

27. Cho RJ, Campbell MJ (2000) Transcription, genomes, function. Trends in Genetics 16: 409–415.RJ ChoMJ Campbell2000Transcription, genomes, function.Trends in Genetics16409415

28. Clark AG, Glanowski S, Nielsen R, Thomas PD, Kejariwal A, et al. (2003) Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios. Science 302(5652): 1960–1963.AG ClarkS. GlanowskiR. NielsenPD ThomasA. Kejariwal2003Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios.Science302(5652)19601963

29. Shaffer AL, Wright G, Yang L, Powell J, Ngo V, et al. (2006) A library of gene expression signatures to illuminate normal and pathological lymphoid biology. Immunol Rev 210: 67–85.AL ShafferG. WrightL. YangJ. PowellV. Ngo2006A library of gene expression signatures to illuminate normal and pathological lymphoid biology.Immunol Rev2106785

30. Yamagata T, Benoist C, Mathis DA (2006) Shared gene-expression signature in innate-like lymphocytes. Immunol Rev 210: 52–66.T. YamagataC. BenoistDA Mathis2006Shared gene-expression signature in innate-like lymphocytes.Immunol Rev2105266

31. Foell D, Roth J (2004) Proinflammatory S100 proteins in arthritis and autoimmune disease. Arthritis Rheum 50(12): 3762–3771.D. FoellJ. Roth2004Proinflammatory S100 proteins in arthritis and autoimmune disease.Arthritis Rheum50(12)37623771

32. Youssef P, Roth J, Frosch M, Costello P, Fitzgerald O, et al. (1999) Expression of myeloid related proteins (MRP) 8 and 14 and the MRP8/14 heterodimer in rheumatoid arthritis synovial membrane. J Rheumatol 26: 2523–2528.P. YoussefJ. RothM. FroschP. CostelloO. Fitzgerald1999Expression of myeloid related proteins (MRP) 8 and 14 and the MRP8/14 heterodimer in rheumatoid arthritis synovial membrane.J Rheumatol2625232528

33. Foell D, Kane D, Bresnihan B, Vogl T, Nacken W, et al. (2003) Expression of the pro-inflammatory protein S100A12 (EN-RAGE) in rheumatoid and psoriatic arthritis. Rheumatology (Oxford) 42: 1383–1389.D. FoellD. KaneB. BresnihanT. VoglW. Nacken2003Expression of the pro-inflammatory protein S100A12 (EN-RAGE) in rheumatoid and psoriatic arthritis.Rheumatology (Oxford)4213831389

34. Frosch M, Strey A, Vogl T, Wulffraat NM, Kuis W, et al. (2000) Myeloid-related proteins 8 and 14 are specifically secreted during interaction of phagocytes and activated endothelium and are useful markers for monitoring disease activity in pauciarticular-onset juvenile rheumatoid arthritis. Arthritis Rheum 43: 628–637.M. FroschA. StreyT. VoglNM WulffraatW. Kuis2000Myeloid-related proteins 8 and 14 are specifically secreted during interaction of phagocytes and activated endothelium and are useful markers for monitoring disease activity in pauciarticular-onset juvenile rheumatoid arthritis.Arthritis Rheum43628637

35. Kane D, Roth J, Frosch M, Vogl T, Bresnihan B, et al. (2003) Increased perivascular synovial membrane expression of myeloid related proteins in psoriatic arthritis. Arthritis Rheum 48: 1676–1685.D. KaneJ. RothM. FroschT. VoglB. Bresnihan2003Increased perivascular synovial membrane expression of myeloid related proteins in psoriatic arthritis.Arthritis Rheum4816761685

36. Rouleau P, Vandal K, Ryckman C, Poubelle PE, Boivin A, et al. (2003) The calcium-binding protein S100A12 induces neutrophil adhesion, migration, and release from bone marrow in mouse at concentrations similar to those found in human inflammatory arthritis. Clin Immunol 107: 46–54.P. RouleauK. VandalC. RyckmanPE PoubelleA. Boivin2003The calcium-binding protein S100A12 induces neutrophil adhesion, migration, and release from bone marrow in mouse at concentrations similar to those found in human inflammatory arthritis.Clin Immunol1074654

37. Chen YS, Yan W, Geczy CL, Brown MA, Thomas R (2009) Serum levels of soluble receptor for advanced glycation end products and of S100 proteins are associated with inflammatory, autoantibody, and classical risk markers of joint and vascular damage in rheumatoid arthritis. Arthritis Res Ther 11: R39.YS ChenW. YanCL GeczyMA BrownR. Thomas2009Serum levels of soluble receptor for advanced glycation end products and of S100 proteins are associated with inflammatory, autoantibody, and classical risk markers of joint and vascular damage in rheumatoid arthritis.Arthritis Res Ther11R39

38. Bovin LF, Rieneck K, Workman C, Nielsen H, Sørensen SF, et al. (2004) Blood cell gene expression profiling in rheumatoid arthritis. Discriminative genes and effect of rheumatoid factor. Immunol Lett 93: 217–226.LF BovinK. RieneckC. WorkmanH. NielsenSF Sørensen2004Blood cell gene expression profiling in rheumatoid arthritis. Discriminative genes and effect of rheumatoid factor.Immunol Lett93217226

39. Batliwalla FM, Baechler EC, Xiao X, Li W, Balasubramanian S, et al. (2005) Peripheral blood gene expression profiling in rheumatoid arthritis. Genes Immun 6: 388–397.FM BatliwallaEC BaechlerX. XiaoW. LiS. Balasubramanian2005Peripheral blood gene expression profiling in rheumatoid arthritis.Genes Immun6388397

40. Utans U, Arceci RJ, Yamashita Y, Russell ME (1995) Cloning and characterization of allograft inflammatory factor-1: a novel macrophage factor identified in rat cardiac allografts with chronic rejection. J Clin Invest 95: 2954–2962.U. UtansRJ ArceciY. YamashitaME Russell1995Cloning and characterization of allograft inflammatory factor-1: a novel macrophage factor identified in rat cardiac allografts with chronic rejection.J Clin Invest9529542962

41. Kimura M, Kawahito Y, Obayashi H, Ohta M, Hara H, et al. (2007) A critical role for allograft inflammatory factor-1 in the pathogenesis of rheumatoid arthritis. J Immunol 178(5): 3316–3322.M. KimuraY. KawahitoH. ObayashiM. OhtaH. Hara2007A critical role for allograft inflammatory factor-1 in the pathogenesis of rheumatoid arthritis.J Immunol178(5)33163322

42. Harney SM, Vilariño-Güell C, Adamopoulos IE, Sims AM, Lawrence RW, et al. (2008) Fine mapping of the MHC Class III region demonstrates association of AIF1 and rheumatoid arthritis. Rheumatology (Oxford) 47(12): 1761–1767.SM HarneyC. Vilariño-GüellIE AdamopoulosAM SimsRW Lawrence2008Fine mapping of the MHC Class III region demonstrates association of AIF1 and rheumatoid arthritis.Rheumatology (Oxford)47(12)17611767

43. Viriyakosol S, Tobias PS, Kitchens RL, Kirkland TN (2001) MD-2 binds to bacterial lipopolysaccharide. J Biol Chem 276: 38044–38051.S. ViriyakosolPS TobiasRL KitchensTN Kirkland2001MD-2 binds to bacterial lipopolysaccharide.J Biol Chem2763804438051

44. Nagai Y, Akashi S, Nagafuku M, Ogata M, Iwakura Y, et al. (2002) Essential role of MD-2 in LPS responsiveness and TLR4 distribution. Nat Immunol 3(7): 667–672.Y. NagaiS. AkashiM. NagafukuM. OgataY. Iwakura2002Essential role of MD-2 in LPS responsiveness and TLR4 distribution.Nat Immunol3(7)667672

45. Brentano F, Kyburz D, Schorr O, Gay R, Gay S (2005) The role of Toll-like receptor signalling in the pathogenesis of arthritis. Cell Immunol 233: 90–96.F. BrentanoD. KyburzO. SchorrR. GayS. Gay2005The role of Toll-like receptor signalling in the pathogenesis of arthritis.Cell Immunol2339096

46. Ospelt C, Brentano F, Rengel Y, Stanczyk J, Kolling C, et al. (2008) Overexpression of toll-like receptors 3 and 4 in synovial tissue from patients with early rheumatoid arthritis: toll-like receptor expression in early and longstanding arthritis. Arthritis Rheum 58(12): 3684–3692.C. OspeltF. BrentanoY. RengelJ. StanczykC. Kolling2008Overexpression of toll-like receptors 3 and 4 in synovial tissue from patients with early rheumatoid arthritis: toll-like receptor expression in early and longstanding arthritis.Arthritis Rheum58(12)36843692

47. O'Neill LA (2008) Primer: Toll-like receptor signaling pathways—what do rheumatologists need to know? Nat Clin Pract Rheumatol 4(6): 319–327.LA O'Neill2008Primer: Toll-like receptor signaling pathways—what do rheumatologists need to know?Nat Clin Pract Rheumatol4(6)319327

48. Macian F (2005) NFAT proteins: key regulators of T-cell development and function. Nat Rev Immunol 5(6): 472–484.F. Macian2005NFAT proteins: key regulators of T-cell development and function.Nat Rev Immunol5(6)472484

49. Masuda K, Masuda R, Neidhart M, Simmen BR, Michel BA, et al. (2002) Molecular profile of synovial fibroblasts in rheumatoid arthritis depends on the stage of proliferation. Arthritis Res Ther 4(5): R8.K. MasudaR. MasudaM. NeidhartBR SimmenBA Michel2002Molecular profile of synovial fibroblasts in rheumatoid arthritis depends on the stage of proliferation.Arthritis Res Ther4(5)R8

50. Filippin LI, Vercelino R, Marroni NP, Xavier RM (2008) Redox signalling and the inflammatory response in rheumatoid arthritis. Clin Exp Immunol 152(3): 415–422.LI FilippinR. VercelinoNP MarroniRM Xavier2008Redox signalling and the inflammatory response in rheumatoid arthritis.Clin Exp Immunol152(3)415422

51. Matthews JR, Wakasugi N, Virelizier JL, Yodoi J, Hay RT (1992) Thioredoxin regulates the DNA binding activity of NF-kB by reduction of a disulfide bond involving cysteine 62. Nucleic Acids Res 20(15): 3821–3830.JR MatthewsN. WakasugiJL VirelizierJ. YodoiRT Hay1992Thioredoxin regulates the DNA binding activity of NF-kB by reduction of a disulfide bond involving cysteine 62.Nucleic Acids Res20(15)38213830

52. Yoshida S, Katoh T, Tetsuka T, Uno K, Matsui N, et al. (1999) Involvement of thioredoxin in rheumatoid arthritis: its costimulatory roles in the TNF-alpha-induced production of IL-6 and IL-8 from cultured synovial fibroblasts. J Immunol 163(1): 351–358.S. YoshidaT. KatohT. TetsukaK. UnoN. Matsui1999Involvement of thioredoxin in rheumatoid arthritis: its costimulatory roles in the TNF-alpha-induced production of IL-6 and IL-8 from cultured synovial fibroblasts.J Immunol163(1)351358

53. Pekkari K, Avila-Cariño J, Bengtsson A, Gurunath R, Scheynius A, et al. (2001) Truncated thioredoxin(Trx80) induces production of interleukin-12 and enhances CD14 expression in human monocytes. Blood 97: 3184–3190.K. PekkariJ. Avila-CariñoA. BengtssonR. GurunathA. Scheynius2001Truncated thioredoxin(Trx80) induces production of interleukin-12 and enhances CD14 expression in human monocytes.Blood9731843190

54. Kim SH, Oh J, Choi JY, Jang JY, Kang MW, et al. (2008) Identification of human thioredoxin as a novel IFN-gamma-induced factor: mechanism of induction and its role in cytokine production. BMC Immunol 9: 64.SH KimJ. OhJY ChoiJY JangMW Kang2008Identification of human thioredoxin as a novel IFN-gamma-induced factor: mechanism of induction and its role in cytokine production.BMC Immunol964

55. Gallo RL, Murakami M, Ohtake T, Zaiou M (2002) Biology and clinical relevance of naturally occurring antimicrobial peptides. J Allergy Clin Immunol 110(6): 823–831.RL GalloM. MurakamiT. OhtakeM. Zaiou2002Biology and clinical relevance of naturally occurring antimicrobial peptides.J Allergy Clin Immunol110(6)823831

56. Gilliet M, Lande R (2008) Antimicrobial peptides and self-DNA in autoimmune skin inflammation. Curr Opin Immunol 20(4): 401–407.M. GillietR. Lande2008Antimicrobial peptides and self-DNA in autoimmune skin inflammation.Curr Opin Immunol20(4)401407

57. Lebre MC, Tak PP (2009) Dendritic cells in rheumatoid arthritis: Which subset should be used as a tool to induce tolerance? Hum Immunol. In press. MC LebrePP Tak2009Dendritic cells in rheumatoid arthritis: Which subset should be used as a tool to induce tolerance?Hum Immunol. In press

58. Yuasa I, Umetsu K, Vogt U, Nakamura H, Nanba E, et al. (1997) Human orosomucoid polymorphism: molecular basis of the three common ORM1 alleles, ORM1*F1, ORM1*F2, and ORM1*S. Hum Genet 99: -393–398.I. YuasaK. UmetsuU. VogtH. NakamuraE. Nanba1997Human orosomucoid polymorphism: molecular basis of the three common ORM1 alleles, ORM1*F1, ORM1*F2, and ORM1*S.Hum Genet99-393398

59. Ceciliani F, Pocacqua V (2007) The acute phase protein alpha1-acid glycoprotein: a model for altered glycosylation during diseases. Curr Protein Pept Sci 8: 91–108.F. CecilianiV. Pocacqua2007The acute phase protein alpha1-acid glycoprotein: a model for altered glycosylation during diseases.Curr Protein Pept Sci891108

60. Hochepied T, Berger FG, Baumann H, Libert C (2003) Alpha(1)-acid glycoprotein: an acute phase protein with inflammatory and immunomodulating properties. Cytokine Growth Factor Rev 14(1): 25–34.T. HochepiedFG BergerH. BaumannC. Libert2003Alpha(1)-acid glycoprotein: an acute phase protein with inflammatory and immunomodulating properties.Cytokine Growth Factor Rev14(1)2534

61. Smith KD, Pollqcchi A, Field M, Watson J (2002) The heterogeneity of the glycosylation of the alpha-1-acid glycoprotein between the sera and synovial fluid in Rheumatoid Arthritis. Biomed Chromatogr 16: 261–266.KD SmithA. PollqcchiM. FieldJ. Watson2002The heterogeneity of the glycosylation of the alpha-1-acid glycoprotein between the sera and synovial fluid in Rheumatoid Arthritis.Biomed Chromatogr16261266

62. Haston JL, FitzGerald O, Kane D, Smith KD (2003) The influence of alpha1-acid glycoprotein on collagenase-3 activity in early rheumatoid arthritis. Biomed Chromatogr 17(6): 361–364.JL HastonO. FitzGeraldD. KaneKD Smith2003The influence of alpha1-acid glycoprotein on collagenase-3 activity in early rheumatoid arthritis.Biomed Chromatogr17(6)361364

63. Awomoyi AA (2007) The human solute carrier family 11 member 1 protein (SLC11A1): linking infections, autoimmunity and cancer? FEMS Immunol Med Microbiol 49(3): 324–329.AA Awomoyi2007The human solute carrier family 11 member 1 protein (SLC11A1): linking infections, autoimmunity and cancer?FEMS Immunol Med Microbiol49(3)324329

64. Telfer JF, Brock JH (2002) Expression of ferritin, transferrin receptor, and non-specific resistance associated macrophage proteins 1 and 2 (Nramp1 and Nramp2) in the human rheumatoid synovium. Ann Rheum Dis 61(8): 741–744.JF TelferJH Brock2002Expression of ferritin, transferrin receptor, and non-specific resistance associated macrophage proteins 1 and 2 (Nramp1 and Nramp2) in the human rheumatoid synovium.Ann Rheum Dis61(8)741744

65. Blackwell JM (1996) Structure and function of the natural-resistance associated macrophage protein (Nramp1), a candidate protein for infectious and autoimmune disease susceptibility. Mol Med Today 2: 205–211.JM Blackwell1996Structure and function of the natural-resistance associated macrophage protein (Nramp1), a candidate protein for infectious and autoimmune disease susceptibility.Mol Med Today2205211

66. Chu CQ, Field M, Feldmann M, Maini RN (1991) Localization of tumor necrosis factor alpha in synovial tissues and at the cartilage–pannus junction in patients with rheumatoid arthritis. Arthritis Rheum 34: 1125–1132.CQ ChuM. FieldM. FeldmannRN Maini1991Localization of tumor necrosis factor alpha in synovial tissues and at the cartilage–pannus junction in patients with rheumatoid arthritis.Arthritis Rheum3411251132

67. Miossec P, Dinarello CA, Ziff M (1986) Interleukin-1 lymphocyte chemotactic activity in rheumatoid arthritis synovial fluid. Arthritis Rheum 29: 461–470.P. MiossecCA DinarelloM. Ziff1986Interleukin-1 lymphocyte chemotactic activity in rheumatoid arthritis synovial fluid.Arthritis Rheum29461470

68. Boneca IG (2009) Mammalian PGRPs in the spotlight. Cell Host Microbe 5(2): 109–111.IG Boneca2009Mammalian PGRPs in the spotlight.Cell Host Microbe5(2)109111

69. Saha S, Qi J, Wang S, Wang M, Li X, et al. (2009) PGLYRP-2 and Nod2 are both required for peptidoglycan-induced arthritis and local inflammation. Cell Host Microbe 5: 137–150.S. SahaJ. QiS. WangM. WangX. Li2009PGLYRP-2 and Nod2 are both required for peptidoglycan-induced arthritis and local inflammation.Cell Host Microbe5137150

70. Dahlbäck B, Villoutreix BO (2005) The anticoagulant protein C pathway. FEBS Lett 579(15): 3310–3316.B. DahlbäckBO Villoutreix2005The anticoagulant protein C pathway.FEBS Lett579(15)33103316

71. Barnhart MI, Riddle JM, Bluhm GB, Quintana C (1967) Fibrin promotion and lysis in arthritic joints. Ann Rheum Dis 26(3): 206–18.MI BarnhartJM RiddleGB BluhmC. Quintana1967Fibrin promotion and lysis in arthritic joints.Ann Rheum Dis26(3)20618

72. Weinberg JB, Pippen AMM, Greenberg CS (1991) Extravascular fibrin formation and dissolution in synovial tissue of patients with osteoarthritis and rheumatoid arthritis. Arthritis Rheum 34: 996–1005.JB WeinbergAMM PippenCS Greenberg1991Extravascular fibrin formation and dissolution in synovial tissue of patients with osteoarthritis and rheumatoid arthritis.Arthritis Rheum349961005

73. Busso N, Hamilton JA (2002) Extravascular coagulation and the plasminogen activator/plasmin system in rheumatoid arthritis. Arthritis Rheum 46(9): 2268–2279.N. BussoJA Hamilton2002Extravascular coagulation and the plasminogen activator/plasmin system in rheumatoid arthritis.Arthritis Rheum46(9)22682279