The arthritis severity locus Cia5a regulates the expression of inflammatory mediators including Syk pathway genes and proteases in pristane-induced arthritis

Background Cia5a is a locus on rat chromosome 10 that regulates disease severity and joint damage in two models of rheumatoid arthritis, collagen- and pristane-induced arthritis (PIA). In this study, we aimed to identify cellular and molecular processes regulated by Cia5a using microarray-based gene expression analysis of synovial tissues from MHC identical DA (severe erosive disease) and DA.F344(Cia5a) congenics (mild non-erosive disease) rats. Results Synovial tissues from six DA and eight DA.F344(Cia5a) rats were analyzed 21 days after the induction of PIA using the Illumina RatRef-12 BeadChip (21,922 genes) and selected data confirmed with qPCR. There was a significantly increased expression of pro-inflammatory mediators such as Il1b (5-fold), Il18 (3.9-fold), Cxcl1 (10-fold), Cxcl13 (7.5-fold) and Ccl7 (7.9-fold), and proteases like Mmp3 (23-fold), Mmp9 (32-fold), Mmp14 (4.4-fold) and cathepsins in synovial tissues from DA, with reciprocally reduced levels in congenics. mRNA levels of 47 members of the Spleen Tyrosine Kinase (Syk) pathway were significantly increased in DA synovial tissues compared with DA.F344(Cia5a), and included Syk (5.4-fold), Syk-activating receptors and interacting proteins, and genes regulated by Syk such as NFkB, and NAPDH oxidase complex genes. Nuclear receptors (NR) such as Rxrg, Pparg and Rev-erba were increased in the protected congenics, and so was the anti-inflammatory NR-target gene Scd1 (54-fold increase). Tnn (72-fold decrease) was the gene most significantly increased in DA. Conclusions Analyses of gene expression in synovial tissues revealed that the arthritis severity locus Cia5a regulates the expression of key mediators of inflammation and joint damage, as well as the expression of members of the Syk pathway. This expression pattern correlates with disease severity and joint damage and along with the gene accounting for Cia5a could become a useful biomarker to identify patients at increased risk for severe and erosive disease. The identification of the gene accounting for Cia5a has the potential to generate a new and important target for therapy and prognosis.


Background
Rheumatoid arthritis (RA) is a common, chronic and potentially debilitating form of autoimmune erosive arthritis. Advances in the understanding of RA pathogenesis have led to the development of new and better treatments [1][2][3]. Yet, sustained remission is still rarely achieved [4], and more effective therapies are needed.
The identification of genes implicated in the regulation of arthritis severity and articular damage has the potential to generate new and potentially better targets for therapies aimed at preserving joint architecture and function, and reducing the risk of developing joint deformities. Yet, little is known about those genes [5], and the large cohorts of RA patients used in genome-wide association studies for susceptibility were not designed to address disease severity and articular damage.
We have previously identified several disease severity and articular damage quantitative trait loci (QTL) in rat models of RA [6][7][8][9][10]. Using a combination of positional cloning and functional studies that include transcriptome analyses of synovial cells and synovial tissues we are beginning to understand the molecular processes regulating arthritis severity and joint damage in pristane-and collagen-induced arthritis (PIA and CIA) [10][11][12][13][14]. Similar strategies have been successfully used to identify other autoimmunity genes in rodent models [15,16].
Cia5a is a 20.6Mb QTL on rat chromosome 10 that regulates arthritis severity, cartilage and bone damage, synovial hyperplasia and inflammation in both PIA and CIA [9,10]. In the present study we used synovial tissues from arthritis-protected DA.F344(Cia5a) congenics and from arthritis-susceptible and MHC-identical DA rats in a microarray analysis of gene expression. We determined that the Cia5a locus regulates the expression of several genes central to RA pathogenesis and joint damage, such as cytokines Il1b and Il18, chemokines, proteases, mediators of the synthesis of reactive oxygen species and prostaglandins, and genes involved in Toll-like receptor signaling. Additionally, the expression of 47 members of the Syk kinase pathway genes, including NFκB genes were significantly regulated by the Cia5a locus. Furthermore, the presence of F344 alleles at the Cia5a interval was associated with increased expression of antiinflammatory genes, including nuclear receptors and Timp3, suggesting that the Cia5a locus contains a gene involved in maintaining an inflammation-free synovial tissue.

Results
DA.F344(Cia5a) congenics develop a mild form of PIA with a distinct pattern of gene expression compared with DA rats DA.F344(Cia5a) rats developed a significantly milder form of PIA compared with DA rats [median arthritis severity score (25-75 percentiles), DA=26.5 (17-36.9), DA.F344(Cia5a)=5.5 (3.6-7.2); p=0.002, Mann-Whitney test; Figure 1A and B]. 36% (7,925) of the genes in the RatRef-12 BeadChip were consistently expressed in synovial tissues. Nearly one-third of these genes (2,648) met the filtering criteria for differential expression (fold-difference ≥1.5 and p≤0.01). The presence of F344 alleles at the Cia5a interval, as in DA.F344(Cia5a) congenic rats, was associated with increased expression of 1,241 genes and reduced expression of 1,407 genes compared with DA. 134 genes had a ≥5-fold difference between strains ( Figure 1C). 46 genes had a ≥10-fold difference in expression, of which 19 were increased and 27 decreased in congenics, compared with DA (Tables 1 and 2).
Expression of pro-inflammatory genes, proteases (including matrix metalloproteases, MMPs) and adhesion molecules was significantly increased in DA and decreased in DA.F344(Cia5a) The 1,407 genes with increased expression in DA and reciprocally decreased expression in DA.F344(Cia5a) congenics included pro-inflammatory cytokines and chemokines implicated in arthritis pathogenesis such as Il1b (5.17-fold on microarray, and 2.46-fold on qPCR), Il18, Mif, Ccl2, Ccl7 and Cxcl13 (Table 3 and Additional file  1: Table S3 and Additional file 2: Table S4). Genes with significantly decreased expression in congenics also included those implicated in the development of cartilage and bone erosions such as MMPs (Mmp3 [24-fold], Mmp9 and Mmp14), and other proteases (cathepsins D, E, K and S) ( Table 3 and Figure 2). Interestingly, Syk (see below) has been shown to regulate the expression of differentially expressed MMPs such as Mmp3 [17] and Mmp9 [18], further suggesting a potential central role for Syk in arthritis and a Syk-regulatory effect of Cia5a. Components of the extracellular matrix (ECM; Cthrc1, Col12a1, Emilin1) also had reduced expression in congenics, and together with the levels of proteases suggested that there was reduced matrix turnover and reduced degradation, compared with arthritic DA rats (Table 3).
Adhesion molecules required for leukocyte migration into the synovium were increased in DA synovial tissues and decreased in DA.F344(Cia5a), including integrins Itga5, Itgam, Itgb2, Itgb7, and Cd44 (Table 3). Cadherin-11 (Cdh11), a FLS-specific gene required for cell-cell interactions and implicated in FLS invasion and synovial hyperplasia was also decreased in DA.F344(Cia5a) congenics, consistent with the non-hyperplastic synovial tissue previously described in this strain, as opposed to the highly hyperplastic synovial tissue seen in DA [10].
The gene with the most significantly increased expression in DA versus DA.F344(Cia5a) was Tnn (Tenascin N; Table 1 and Figure 2). Tnn has been implicated in osteogenesis and angiogenesis but not in arthritis or inflammation.
These results demonstrate that DA rats with PIA have increased synovial expression of many genes implicated in RA pathogenesis, further validating the molecular similarities between PIA and RA, and underscoring the potential relevance of both Cia5a in arthritis pathogenesis and the present study in discovering new key genes and pathways regulating arthritis.   (Table 4). Interestingly, genes that neutralize ROS (Cat, Sod1, Gss) went on the opposite direction with increased expression in congenics (Table 5).
Additionally, Syk and Vav1 expression levels correlated with the cumulative arthritis severity score (Pearson's correlation coefficient of 0.8 and p=0.0006 for both genes). Taken together, these observations suggest that the Cia5a QTL contains an arthritis gene that directly or indirectly regulates the expression of Syk pathway genes, providing a possible mechanistic explanation for this locus' effect on the regulation of disease severity.   Table 3) were expressed in increased levels in DA and decreased in congenics. Interestingly, and in line with these observations, the expression levels of negative regulators of TLR signaling such as Ptpn11 and Pparg was conversely increased in congenics, suggesting that the arthritis gene located within the Cia5a QTL might mediate the balance between activating and inhibitory signals implicated in TLR signaling. Genes implicated in the synthesis of prostaglandins and leukotrienes (Pla2g4a, Ptgs2/Cox2, Ptges) were also increased in DA (Table 3), while genes that counteract eicosanoid-mediated inflammation (Ptgis, Cyp2j3) were increased in congenics (Table 5).
Increased expression of anti-inflammatory genes, including nuclear receptors (NRs), in synovial tissues from DA.F344(Cia5a) Congenics Several genes with known anti-inflammatory and cytokinesuppressing properties were expressed in increased levels in DA.F344(Cia5a) synovial tissues, and reduced in DA. These Figure 2 qPCR validation of the microarray results. Genes expressed in increased levels in DA (six genes, black bars) and genes expressed in increased levels in DA.F344(Cia5a) congenics (three genes, grey bars) were selected for qPCR confirmation. Fold-differences were log-transformed. The same RNAs used in the microarray experiments were used for qPCR. ΔCt was used for statistical analyses; all genes had p≤0.04, t-test..   Syk gene is underlined. * t-test; §=All four chemokine and cytokine listed genes are NFkB-inducible targets. ** Ezrin is also involved in cytoskeleton regulation. ¶ Rac also regulates NFkB activity, cytoskeleton and cell proliferation. included the NRs Lxra, Pparg, Rev-erba, Rora, Thra, and Thrb (Table 5). Scd1 was the gene with the most significantly increased expression in DA.F344(Cia5a) congenics with a 55-fold difference compared with DA (Tables 2 and 5). Scd1 has been shown to reduce cytokine levels and to have anti-inflammatory activity [19,20]. We have previously reported that Scd1 is expressed in significantly reduced levels in synovial tissues from rats with severe arthritis, and increased in the synovial tissues of yet another arthritis-protected congenic strain [14].
Scd1 and some of the other genes up-regulated in DA. F344(Cia5a) synovial tissues such as Adipoq, Cidea, Cd36, Fabp4, Gpd1, Lpl, Lpin1, Mgst1, Plin, Pck1, Slc2a4 and Srebf1, are known to be inducible by NRs (Table 2  and Table 5). These observations suggest that NRs were not only expressed in increased levels but also had increased activity in synovial tissues from DA.F344 (Cia5a) compared with DA.
Genes located within the Cia5a interval have significantly different expression levels 75 of the 7,925 genes expressed by all samples were located within the Cia5a interval. 21 of these 75 had increased expression in DA synovial tissues, and 11 were increased in the congenics. 14 of these 32 differentially expressed genes had ≥2-fold-difference. Sphk1 and Sectm1b were the genes contained within the Cia5a interval with the most significantly increased expression in DA (7.58 and 7.61-fold, respectively), while Itgb4 and DIgr1 and were those with the most significantly increased expression in DA.F344(Cia5a) congenics (2.77 and 2.85-fold, respectively) ( Table 6). Additionally, four genes located within the Cia5a interval were expressed only or predominantly in DA synovium, while two other genes were expressed predominantly in DA.F344(Cia5a) ( Table 6). It is conceivable that these differences in expression levels of genes located within the Cia5a interval could be explained at least in part by a polymorphism in the 5' untranslated region (UTR) that affects a transcription factor binding site in cis, thus affecting transcription efficiency, or a 3' UTR polymorphism affecting mRNA stability.
Gene targets of microRNAs (miRNA) contained within the Cia5a interval were not differentially expressed The Cia5a interval contains six predicted miRNAs. We considered the possibility that polymorphisms in one of those six miRNAs could account for the Cia5a effect on gene expression and arthritis severity. In that case, such polymorphism would affect the miRNA activity on the transcription of its target genes. Therefore, we look for the differential expression of targets of all six predicted miRNAs located within the Cia5a interval. A list of target and non-target genes was generated for each of the six miRNAs, but no significant over-representation of targets was detected (Figure 4), suggesting that polymorphisms affecting the expression or function of the miRNAs contained within the Cia5a interval are less likely to explain the differences in gene expression identified in this study.
Analyses of cell type specific genes suggests synovial tissue cellularity differences between DA and DA.F344(Cia5a) Congenics 13 genes known to be specifically expressed by the cell types of interest to this study were used to compare DA and DA.F344(Cia5a) synovial tissues (Additional file 3: Table S2). The expression levels of those genes suggested increased numbers of FLS (consistent with synovial hyperplasia), macrophages, dendritic cells (DC), neutrophils and T cells in the synovial tissues of DA, compared with congenics (Additional file 3: Table S2). No gene specific for B cells, NK cells or Tregs were among those differentially expressed between the two strains, suggesting that the number of these cells in the synovial tissues of these two strains was not significantly different.

Discussion
Disease severity and articular damage are associated with increased risk for disability, joint deformities and reduced life expectancy in patients with RA [21][22][23]. Yet, little is known about the genes implicated in the regulation of disease severity and articular damage genes in RA, and these genes could be the most relevant targets for new therapies aimed at preserving the joint architecture and function.
We have previously identified Cia5a, a 20.6 Mb arthritis severity and joint damage regulatory locus, on rat chromosome 10 [10]. Cia5a co-localizes with other arthritis severity loci identified in other rodent models of arthritis such as oil-induced arthritis (Oia3) [24], and CIA in a DAxACI intercross (Cia27) [25]. There have been no genome-wide association or linkage studies of disease severity and joint damage in RA, and therefore, it is unknown whether the Cia5a syntenic region on human chromosome 17q22-q25 harbors a severity or joint damage arthritis regulatory gene. However, the human 17q22-q25 region contains a locus previously linked with RA susceptibility [26]. In the present study we analyzed synovial tissues from DA rats, which develop severe arthritis (PIA) with pronounced synovial hyperplasia and cartilage and bone destruction, and synovial tissues from the DA.F344(Cia5a) congenics, which develop mild and non-erosive disease. These two strains share the same MHC and are genetically identical except for the presence of F344 alleles at the Cia5a interval, underscoring the magnitude of the effect of this single locus on clinical disease, on histologic joint damage [10] and on gene expression (present study). DA.F344(Cia5a) congenics had significantly reduced expression of genes previously implicated in RA pathogenesis, RA severity and articular damage, including Il1b, Il18, Mif, Mmp3 and Mmp14. These and other similarities between DA rats and RA synovial tissues' gene expression, such as increased expression of chemokines, matrix proteins, adhesion molecules, mediators of innate immune responses, and others, underscore and further validate the potential clinical relevance of our model and discovery strategy.
We identified a new role for Cia5a on the regulation of the expression of members of the Syk pathway, where forty-seven genes directly or indirectly related to Syk activation were expressed in increased levels in DA, and significantly reduced levels in DA.F344(Cia5a) congenics. Syk is a tyrosine kinase that phosphorylates ITAM motifs in trans-membrane receptors or adaptors, and interacts with partners like Vav, PI3K and Slp76 [27]. Syk activation mediates signaling through several cell surface receptors, including those with significantly different levels in this study such C-lectin type receptors, Fcer1g, Fcgr2a, Trem2, Tyrobp, integrins, and the T-cell receptor (TCR) (Figure 3). Resident and infiltrating inflammatory cells in the RA synovial pannus, such as mast cells, macrophages, B and T cells, express these Syk-activating receptors. These resident cells and infiltrating cells have been implicated in arthritis pathogenesis and joint damage, raising the possibility that part of their effect may be mediated by Syk-activating receptors.
Analyses of cell-specific genes suggested reduced numbers of macrophages, dendritic cells, neutrophils and T cells in the synovial tissues of congenics compared with DA, which is in agreement with our previous histologic analyses and might explain part, but not all of the differences in expression of Syk genes. Additionally, DA.F344 (Cia5a) congenics had significantly lower levels of the FLS-specific gene Cdh11, compatible with the reduced synovial hyperplasia that we have previously described.
Syk pathway members regulate several cellular processes implicated in arthritis pathogenesis and articular damage, ranging from the production of reactive oxygen species, NFκB activation and the transcription of pro-inflammatory mediators such as Il1b and Ccl2, to the cell proliferation required for the development of synovial hyperplasia, and actin cytoskeleton rearrangements [27]. NFκB activity is Figure 4 Differentially expressed predicted targets of the miRNAs located within the Cia5a interval. Six miRNAs map to the Cia5a interval. The number of predicted target genes of each miRNA that were up-regulated and down-regulated in DA and DA.F344(Cia5a) congenics was not statistically significant (Chi-square test with the Yates correction). Numbers (percentages intra-parenthesis) of genes expressed in increased levels (dark grey), reduced levels (light grey) or unchanged (white) in DA compared with DA.F344(Cia5a) are shown. regulated by Syk and by several other pathways including TLRs and cytokine receptors [28]. The NFκB pathway has a central role in the production of pro-inflammatory cytokines such as IL-1β, IL-6 and TNFα, in the development of synovial hyperplasia and in disease severity [29][30][31]. Actin cytoskeleton rearrangements are also regulated by the Syk pathway [27], and are required for the migration of inflammatory cells into the synovial tissue, and for synovial cells and synovial tissue invasion and destruction of cartilage [13,32]. Therefore, our observations suggest that a gene located within the Cia5a interval is a new regulator of the expression of Syk pathway genes implicated in key processes in arthritis pathogenesis.
The precise mechanisms whereby Cia5a regulates the expression of Syk genes remain unclear, and might reflect differences in tissue cellularity, multiple cell-activating processes, or a polymorphism in transcription factor located within the Cia5a interval that affects transcription. Studies by our group of synovial tissues obtained from four different congenic strains yielded different results in gene expression (Brenner et al., manuscript in preparation) [12,14,33], suggesting that the Syk-regulatory effect of Cia5a is a specific observation, and not simply related cellularity differences or inflammation.
Syk has been recently implicated in arthritis pathogenesis and joint damage, and Syk-deficient mice are protected from autoantibody-induced erosive arthritis [34], and treatment with a SYK inhibitor significantly reduced disease severity and joint erosions and damage in collagen-induced arthritis [35]. Both the total and phosphorylated forms of SYK are expressed in increased levels in RA synovial tissues compared with osteoarthritis, and SYK inhibition reduced the expression of IL-6 and MMP-3 [17]. More importantly, the use of a SYK inhibitor significantly reduced disease activity in patients with RA [36], with 67%, 43% and 28% of patients achieving ACR20, ACR50 and ACR70, respectively, in a 3month double-blind and placebo-controlled study [37]. Therefore, it is conceivable that the Syk pathway genes differentially expressed in this study could help identify patients more likely to benefit from therapy with SYK inhibitors. Additionally, Syk is critical to TNFα-induced responses [38], raising the possibility that the Syk pathway 47-gene signature could be used to predict increased TNFα activity prior to choosing a biologic therapeutic agent. Additionally, the increased expression of Syk pathway genes could identify patients at increased risk to develop erosive disease and could become a prognostic tool. Lastly, the Cia5a gene itself has the potential to become a new target for therapies aimed at reducing articular damage via inhibition of Syk pathway genes, including processes downstream from Syk such as NFκB.
While several genes with pro-inflammatory, proteolytic, innate immunity and inflamasome-related activity were expressed in reduced levels in DA.F344(Cia5a) congenics, groups of genes with known anti-inflammatory properties were expressed in increased levels in congenics. These genes included Timp3, Ptpn11, antagonists of reactive oxygen species (Cat, Gss, Sod1) and nuclear receptors. Nuclear receptors such as Lxra, Pparg and Rora have been shown to interfere with NFκB and AP-1 activation [39][40][41], and to have anti-inflammatory and arthritis-suppressive properties [42][43][44][45]. Rxrg was another nuclear receptor expressed in significantly increased levels in DA.F344(Cia5a) congenics. While Rxrg itself has not been studied in the context of arthritis, it dimerizes with Lxra, Pparg, and with Vdr, and is required for their anti-inflammatory activity. Additionally, several nuclear receptor-inducible genes, including the inflammation-suppressor Scd1 [20] were expressed in increased levels in the synovial tissues of the congenics. These observations suggest that not only nuclear receptor levels were increased, but also their activity. We have recently identified a similar nuclear receptor expression signature in another arthritis-protective congenic strain, DA.ACI(Cia25) [14], suggesting that this effect is not specific to the Cia5a locus, but more broadly correlates with preservation of both a normal synovial environment and articular architecture.
The gene with the most significantly increased expression in DA compared with congenics was Tnn (Tenascin N). While little is known about this secreted extracellular matrix glycoprotein, it has been implicated in cancer-associated angiogenesis [46], and in integrindependent cancer motility [47]. Another member of the tenascin family, Tenascin C (Tnc), was recently shown to be an endogenous activator of TLR4, an inducer of IL-6 and TNFα, and was required for joint damage in arthritic mice [48], suggesting that Tnn could have a function similar to Tnc in arthritis.
Lastly, we considered the possibility that a polymorphism in the 5' UTR or 3'UTR region of the gene accounting for Cia5a could interfere with its transcription and/ or mRNA stability, respectively, leading to increased or reduced gene-specific mRNA levels. We looked for differentially expressed genes and genes preferentially expressed by only one of the strains and located within Cia5a as a clue to the above possibility. Thirty-eight genes met these criteria, and particularly the most significant sixteen genes are interesting candidates that will be studied in detail (Table 6).

Conclusion
In conclusion, in the present study we identified a pattern of gene expression regulated by Cia5a, which included several inflammatory mediators and 47 members of the Syk pathway. Levels of several mediators of arthritis pathogenesis, synovial hyperplasia and articular damage were also reduced in DA.F344(Cia5a) congenics, underscoring the importance of the gene accounting for this locus. Increased expression of nuclear receptors correlated with joint preservation, and a new potential mediator of inflammation, Tnn, was identified for the first time in synovial tissues. Our observations suggest that the gene accounting for Cia5a has the potential to become an important new target for therapies aimed at preserving joint architecture free of damage, and reducing inflammation. (Cia5a) congenic rats were generated as previously described [9]. Briefly, a 20.6Mb interval from chromosome 10 from the arthritis-resistant F344 strain was introgressed into arthritis-susceptible DA rats through genotype-guided breeding ( Figure 1A). This strategy selected for F344 alleles at the Cia5a interval while excluding donor genome contamination at other loci known to regulate arthritis [10,49]. Experiments were done with 8-12 week-old male rats homozygous at the congenic interval. All experiments were conducted under an Institutional Animal Care and Use Committee (IACUC)-approved protocol.

PIA and tissue collection
Male DA (n=6) and DA.F344(Cia5a) (n=8) congenic rats were anesthetized and injected intradermally with 150 μl of pristane (MP Bio, Solon, OH) divided into two injection sites at the base of the tail (day 0) [50]. Arthritis severity was assessed with a previously described 80-point scoring system [51]. Ankle synovial tissues were collected 21 days post-induction of arthritis.

RNA extraction
Total RNA was extracted from synovial tissues using the RNeasy Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's instructions and including a DNase treatment step. RNA was quantified and assessed for purity using the NanoDrop spectrophotometer (Rockland, DE). RNA integrity was verified with the BioAnalyzer 2100 (Agilent, Palo Alto, CA).

Microarray
All reagents and procedures were previously optimized for use with the Illumina Whole-Genome Expression platform [12]. Briefly, total RNA (200 ng) was amplified and biotinylated using the TotalPrep labeling kit (Ambion, Austin, TX). Each individual sample was hybridized to one individual array in the RatRef-12 Expression BeadChip (Illumina, San Diego, CA), which contains 22,522 probes covering 21,922 rat genes selected primarily from the NCBI RefSeq database (Release 16). Hybridization was done in Illumina IntelliHyb chambers, followed by washing and staining with Cy3streptavidin. The BeadChip was scanned on a highresolution Illumina BeadArray reader using a twochannel 0.8 μm resolution confocal laser scanner.
cDNA synthesis and quantitative real-time PCR (qPCR) expression analysis Differences in the expression of selected genes from the microarray analyses were validated with qPCR. The qPCR conditions have been described elsewhere [12]. Briefly, total RNA (200 ng) from each sample was used for cDNA synthesis using Superscript III (Invitrogen). Primers and qPCR probes were designed to target the same exons as the corresponding Illumina RatRef-12 Expression BeadChip probes (Additional file 4: Table S1). We used Universal ProbeLibrary (Roche, Indianapolis, IN) and Taqman (ABI, Applied Biosystems, Foster City, CA) probes labeled with FAM at the 5' end and TAMRA at 3' end. Reactions were prepared in duplicates with Eurogentec qPCR MasterMix (San Diego, CA), and run on an ABI Prism 7700 thermocycler using SDS software version 1.9.1 (ABI). Ct (threshold cycle) values were adjusted for Gapdh in each sample (ΔCt). Expression levels (ΔCt) were compared using the t-test and a p-value ≤0.01 was considered significant. Fold-differences were calculated with the 2 -ΔΔCt method [52].

MicroRNAs (miRNA)
We considered the possibility that polymorphisms in a miRNA located within the Cia5a interval could account for the Cia5a effect on gene expression and arthritis severity. In that case, such polymorphism would affect the miRNA activity on the transcription of its target genes. Therefore, we looked for miRNAs mapping to the Cia5a interval using the miRBase [53]. Target genes for the miRNAs contained within the Cia5a interval were predicted with TargetScan [54]. Enrichment for differentially expressed predicted targets of miRNAs located within Cia5a was calculated using the Chi-square test with the Yates correction.

Cellular subset gene expression signatures
Differences in tissue resident and infiltrating cell populations can affect the interpretation of gene expression analyses. We looked for cell-specific genes using the GNF Mouse GeneAtlas V3 (Affymetrix MOE430, GEO code GSE10246), a database containing gene expression information for 96 resting and stimulated mouse cell types and tissues [55,56], as well as the BioGPS website (www.biogps.org, Scripps Research Institute). The GNF Mouse GeneAtlas V3 did not include fibroblast-like synoviocytes (FLS) or regulatory T cells (Tregs). Therefore, additional non-redundant cell signature genes were obtained from the literature to represent FLS and Tregs [57][58][59][60][61][62][63][64][65]. We generated a list of genes specific for B cells, T cells, Treg cells, NK cells, FLS, dendritic cells, mast cells, macrophages and neutrophils. We next looked for those cell-specific genes within the list of genes differentially expressed between DA and DA.F344(Cia5a) congenics, as well as in the list of genes preferentially expressed, or only expressed in one strain and not in the other in order to gain insight into differences in cell populations.

Microarray analysis and statistics
Microarray fluorescence intensities were extracted using BeadStudio 2.0 (Illumina). Fluorescent intensities were background-subtracted and then normalized using the cubic spline algorithm. Normalized data were log 2 -transformed prior to all analyses. Probes consistently expressed in all arrays were included in the analyses. Genes with ≥1.5-fold difference in intensity between DA and DA.F344(Cia5a) and a t-test p-value ≤ 0.01 were considered differentially expressed and selected for pathway detection analyses using IPA 5.5.1 (Ingenuity Systems, Redwood City, CA), as well as public online databases (Ensembl, Genecards, Oncomine, BioGPS, Rat Genome Database) and literature search (Pubmed).
Strain-specific (genes only in one of the strains), or preferential strain (genes expressed in a higher percentage of rats of one strain, and in lower percentage of rats of the other strain) gene expression was determined with the Fisher's exact test.
Enrichment for biological functions and disease groups was determined with the IPA software and calculated using the Fisher's exact test with the Benjamini-Hochberg correction and a cutoff p-value of ≤0.05. Enrichment for differentially expressed genes within specific cell subsets, or genes located within the Cia5a interval was calculated using the Fisher's exact test. Non-normally distributed arthritis severity scores were compared using the Mann-Whitney test.
Funded by a Postdoctoral Fellowship Award from the New Jersey Chapter of the Arthritis Foundation to Dr. M. Brenner, and by the National Institutes of Health grants R01-AR46213, R01-AR052439 (NIAMS) and R01-AI54348 (NIAID) to Dr. P. Gulko.