|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





* Laboratory of Experimental Rheumatology, Robert S. Boas Center for Genomics and Human Genetics and
Graduate School of Molecular Medicine, North Shore-Long Island Jewish (LIJ) Research Institute, Manhasset, NY 11030;
Arthritis and Rheumatism Branch, National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS), National Institutes of Health (NIH), Bethesda, MD 20892;
Research Service, Veterans Affairs Medical Center and Department of Medicine/Rheumatology, University of Utah, Salt Lake City, UT 84132; and ¶ Department of Medicine, New York University School of Medicine, New York, NY
| Abstract |
|---|
|
|
|---|
mRNA were reduced. Although both DA.F344(Cia5) and DA.F344(Cia5a) rats were protected in CIA, the arthritis scores of DA.F344(Cia5) were significantly higher than those of DA.F344(Cia5a), suggesting the existence of a third locus where F344-derived alleles centromeric from Cia5a contribute to increased arthritis severity. The existence of the third locus was further supported by higher levels of autoantibodies against rat type II collagen in DA.F344(Cia5) congenics compared with DA.F344(Cia5a). Our results determined that Cia5 contains three major arthritis severity regulatory loci regulating central events in the pathogenesis of arthritis, and differentially influencing CIA and PIA. These loci are syntenic to regions on human chromosomes 17q and 5q implicated in the susceptibility to rheumatoid arthritis, suggesting that the identification of these genes will be relevant to human disease. | Introduction |
|---|
|
|
|---|
The identification of RA genes has been hampered by factors inherent to this complex trait such as variable penetrance, variable disease allele-associated relative risk, epistasis, genetic heterogeneity of human populations, and complex interactions between environmental and genetic factors (16). The study of rat models of arthritis in inbred strains that differ in their susceptibility to and severity of arthritis reduces several of these confounding variables (17, 18), thus facilitating the identification of loci regulating arthritis-related phenotypes, including disease severity (19, 20, 21, 22, 23, 24, 25, 26, 27, 28), chronicity (23, 29, 30), as well as specific disease features, such as the production of rheumatoid factors (31). Once identified, these rat arthritis genes will provide novel targets for focused candidate gene (or candidate pathway) analysis in RA case-control studies.
Previous genome-wide screens have identified several quantitative trait loci (QTL) regulating disease severity in various rat models of RA (19, 20, 21, 22, 23, 24, 25, 26, 27, 28). One of the identified QTLs, Cia5, was mapped to rat chromosome 10 in an intercross between the arthritis-susceptible DA and arthritis-resistant F344 rat strains studied for collagen-induced arthritis (CIA). Studies in congenic rats confirmed the arthritis-regulatory effect associated with Cia5 and suggested that this QTL could be accounted for by at least two distinct genes (32). Moreover, this QTL was involved in the regulation of arthritis severity in collagen-, pristane-, oil- (33), and adjuvant-induced arthritis (34). Cia5 colocalizes with QTLs involved in the regulation of other models of autoimmune diseases (35, 36, 37) and its human syntenic regions contain loci regulating different forms of autoimmune diseases (discussed in Ref.32), including RA (4, 38), suggesting that it harbors genes relevant not only to RA, but possibly to other diseases as well.
To localize and to narrow the intervals containing the two arthritis-regulatory genes, and to comprehensively characterize their regulatory effects in arthritis severity, joint histology, and synovial tissue cytokine gene expression, Cia5 subcongenics were generated and studied for their susceptibility to and severity of pristane-induced arthritis (PIA) and CIA, two well-established models of autoimmune erosive arthritis that share certain similarities in their pathogenesis and genetic regulation, but also have differences (23, 25, 33, 39, 40). We determined that Cia5 contains three different QTLs that regulate arthritis severity. In PIA, Cia5a and Cia5d also regulate synovial hyperplasia, pannus formation, angiogenesis, synovial inflammation, and bone and cartilage erosions, suggesting that these two genes control fundamental processes in the pathogenesis of arthritis. The novel QTL, Cia5d, is located in the more centromeric portion of Cia5, and regulates PIA only, while Cia5a regulates both PIA and CIA and was reduced to 10 Mb. A third and novel locus specific for CIA was identified in the region centromeric from Cia5a in which the arthritis severity favoring alleles come from F344. This third locus regulates levels of anti-collagen autoantibodies. The identification of the genes accounting for these three loci should increase our understanding of the processes regulating arthritis and point to novel candidate genes and pathways for studies in human disease.
| Materials and Methods |
|---|
|
|
|---|
Specific pathogen-free DA (DA/BklArb; arthritis-susceptible) and F344 (F344/Hsd, Harlan; arthritis-resistant) inbred rat strains were used in the breeding of the congenic and subcongenic strains. DA/BklArb were originally purchased from Bantin & Kingman, maintained at the Arthritis and Rheumatism Branch, NIAMS, NIH, and then transferred to the North Shore-LIJ Research Institute (NSLIJRI) (DA/BklArbNsi) and used as controls. All the experiments involving animals were reviewed and approved by the NSLIJRI Institutional Animal Care and Use Committee.
Construction of the genotype-guided Cia5 QTL-congenic and subcongenic lines
An 89.7-Mb interval, containing the 15 Mb (16.5 cM) two logarithm of odds support interval comprising Cia5, was introgressed from F344 into the DA rats through eight backcrosses (BC8) followed by at least five intercrosses, as previously described (32). DA.F344(Cia5a) subcongenics have been previously generated (32). The other subcongenic lines covering the Cia5 interval (see Fig. 1) were generated for the present study. DA.F344(Cia5) congenics were backcrossed twice with DA rats, and then screened for new recombinants within the Cia5 interval (see simple sequence length polymorphism, markers used on Fig. 1). Offspring (BC10) heterozygous at similar recombinant segments were brother-sister mated, and their offspring (BC10F1) genotyped to ensure homozygosity of the expected intervals. Homozygous subcongenics were used to expand the subcongenic lines. Experiments were done with offspring from second to fifth intercrosses (BC10F2-F5). Homozygous DA.F344(Cia5a) rats were backcrossed with DA animals to generate DA.F344(Cia5a)BC11 rats, heterozygous at the Cia5a interval. These rats were also studied in PIA.
|
Tail tips were excised from 3- to 4-wk-old rats and DNA was extracted with the DNeasy kit (Qiagen). PCR conditions have been previously reported and were set-up in 10-µl reactions (30, 32). Briefly, [32P]PCR products were resolved in acrylamide sequencing gels and fluorescent-based PCRs using a fluorescent-labeled forward primer were run in ABI 310 capillary genotypers/sequencers (Applied Biosystems). GENESCAN 3.1 software (Applied Biosystems) was used for data extraction and allele assignment of the fluorescent-labeled PCR products. All genotypes were manually checked by two readers and questionable readings rechecked or repeated. For marker details, see the Rat Genetic Database (
www.niams.nih.gov/rtbc/ratgbase/index.htm
) and the Rat Genome Database (
www.rgd.mcw.edu
).
Induction of PIA
Eight- to 12-wk-old rats received 150 µl of pristane (2,6,10,14-tetramethylpentadecane; Sigma-Aldrich) by intradermal injection (day 0) (33, 39). The dose was divided between two injection sites at the base of the tail.
Induction of CIA
Bovine type II collagen (BII; Chondrex) was dissolved overnight in 0.1 N acetic acid at 4°C (2 mg/ml) and emulsified with IFA (Difco) to a final concentration of 1 mg/ml. Eight- to 12-wk-old rats were injected intradermally at the base of the tail with 2 mg/kg of BII divided into six injection sites (day 0), and a booster injection of 100 µg of BII/IFA administered on day 7 (20). Serum was obtained on day 18 and stored at 80°C until tested.
Arthritis scoring
We used a previously described arthritis scoring system (19, 20, 21, 30) that evaluates individual joints and measures arthritis severity according to joint size as follows: 1) interphalangeal, metacarpophalangeal, and metatarsophalangeal joints in each one of the four lateral digits were scored: 0, no arthritis; 1, arthritis present; 2) wrist, midforepaw, ankle and midfoot joints were scored: 0, normal; 1, minimal swelling; 2, moderate swelling; 3, severe swelling; 4, severe swelling and non-weight bearing. The scores from all involved joints were added (maximum score per rat = 80). The same observer obtained the arthritis scores and weights on days 0, 14, 18, 21, 24, 28, and 31 following induction. The arthritis severity index (ASI), which is a measure of disease severity over time (area under the curve), was determined for each animal by summing the individual arthritis scores obtained over the course of the experiment. The ASI was used because it provides a more comprehensive evaluation of arthritis severity.
Histology and histologic scoring
At the end of the arthritis observation period (day 32), the right hind paw was fixed in 10% formaldehyde. Paws were then decalcified with a solution containing hydrochloric acid and 0.1 M EDTA (Cal-Ex; Fisher Scientific). Tissues were sectioned, embedded in paraffin, and slides prepared and stained with H&E and safranin-O. Slides were scored without knowledge of strain identity. We used a recently described comprehensive histologic scoring system (41). Briefly, tibiotalar, taluscalcaneal, and midfoot joints were histologically scored for the following parameters:
Synovial inflammation. Five high-power magnification fields (HMF) were scored for the percentage of infiltrating mononuclear cells as follows: 0, absent; 1, mild (110%); 2, moderate (1150%); 3, severe (51100%). The mean of the five HMF was used for analyses.
Synovial hyperplasia. 0, absent; 1, mild (510 layers); 2, moderate (1120 layers); 3, severe (>20 layers).
Extension of pannus formation based on the readers impression. 0, absent; 1, mild; 2, moderate; 3, severe.
Synovial fibrosis. 0, absent; 1, mild (110%); 2, moderate (1150%); 3, severe (51100%).
Synovial vascularity (angiogenesis). The number of vessels was counted in five HMF of synovial tissue and the mean used for analyses.
Cartilage erosion. Percentage of the cartilage surface that was eroded: 0, absent; 1, mild (110%); 2, moderate (1150%); 3, severe (51100%).
Cartilage degradation. Based on safranin-O staining of proteoglycans and described as the percentage of the cartilage that lost its staining: 0, none; 1, mild loss (110%); 2, moderate loss (1150%); 3, severe loss (51100%).
Bone erosion. 0, none; 1, minor erosion(s) observed only at HMF; 2, moderate erosion(s) observed at low magnification; 3, severe transcortical erosion(s).
Quantitative real-time PCR (RT-PCR)
Left ankle synovial tissue obtained after the completion of the arthritis observation period (day 32) was immediately frozen in liquid nitrogen. Tissues were subsequently homogenized and total RNA isolated with the RNeasy kit (Qiagen) and digested with DNase (Qiagen), according to the manufacturers protocol. Two micrograms of total RNA from each sample were used for cDNA synthesis (SuperScript II kit; Invitrogen Life Technologies). The quantitative real-time PCR method used was recently reported (41). Briefly, cDNA was optimized for relative gene expression by RT-PCR. The TaqMan 5' exonuclease assay (Applied Biosystems) was used for the quantification of the rat TNF-
and IL-1
genes. GAPDH was used as endogenous control. TaqMan primers and probes were designed using Primer Express software, version 1.5 (Applied Biosystems; Ref.41). Probes were labeled with FAM at the 5' end and TAMRA at the 3' end. Eurogentec qPCR mastermix reagents were used in an ABI 7700 Sequence Detection System (Applied Biosystems). The PCR mixture contained 1x mastermix, 500 nM of either TNF-
, IL-1
, or GAPDH forward and reverse primers, 200 nM TNF-
or IL-1
, or 100 nM GAPDH gene-specific TaqMan probes, and 5 µl of 1/20 diluted cDNA to a final volume of 25 µl. All samples were run in duplicates and the mean used for the analyses. Thermocycler conditions were 50°C for 2 min, 95°C for 10 min and 45 cycles of 95° C for 0.15 min and 60° C for 1 min. Data were analyzed using Sequence Detection System software version 1.9.1 (Applied Biosystems). Results were obtained as Ct (threshold cycle) values. The Ct value is inversely proportional to the starting template copy number. Relative expression of TNF-
and IL-1
in synovial tissues was adjusted for GAPDH in each sample (
Ct).
Ct values were compared with the Students t test. Fold-change differences in gene expression between DA, DA.F344(Cia5), and subcongenics were compared with the 2
Ct method (Ref.42 and ABI 7700 Sequence Detection System User Bulletin No. 2).
Measurement of anti-collagen Abs
Serum samples collected on day 18 after CIA induction were assayed for IgG Abs against BII and rat type II collagen (RII) using commercially available ELISA kits (Chondrex), according to the manufacturers instructions. The IgG Ab concentration in micrograms per milliliter was calculated using the equivalence ratio of 82 U/ml to 1 µg/ml.
Statistical analyses
Based on our previous observation of gender differences in the genetic regulation of CIA (32, 33), males and females were studied separately for their susceptibility to and severity of arthritis. Non-normally distributed data (medians) were compared with ANOVA on ranks with a pairwise multiple comparison procedure (Dunns method) for multiple groups, or with the Mann-Whitney rank sum test for two group comparisons. The Student t test was used to compare normally distributed data (gene expression). A p value of 0.05 or less was regarded as significant. All statistical analyses were done with SigmaStat 3.0 (SPSS).
| Results |
|---|
|
|
|---|
The introgression of the F344-derived Cia5-containing interval into the DA background, as in the DA.F344(Cia5) congenics (Fig. 1), was associated with a highly significant reduction of 97 and 98% in median ASI in male and female congenics, respectively, as compared with same sex DA rats (Fig. 2A) (median ASI: male DA = 127, male DA.F344(Cia5) = 4, p
0.001; female DA = 114, female DA.F344(Cia5) = 2, p
0.001). This protective effect could be noticed at day 14 (Fig. 2A) and persisted throughout the observation period, suggesting that the arthritis gene located within Cia5 regulates both early/induction and chronic/perpetuation stages during the course of disease pathogenesis.
|
0.001; female DA = 114, female DA.F344(Cia5a) = 12, p
0.001). To determine the mode of inheritance of Cia5a, DA.F344(Cia5a) rats were backcrossed with DA and males (BC11) heterozygous at the Cia5a interval and studied in PIA experiments. Cia5a heterozygous (het) rats developed similar ASI as those of Cia5a homozygous congenics (p = 0.064) and significantly lower ASI compared with DA (median male ASI: DA = 127, DA.F344(Cia5a)het = 14.5; p
0.001) (data not shown), suggesting that the DA-originated arthritis severity alleles operate in a recessive manner.
Both male and female DA.F344(Cia5e) subcongenics (60.57-Mb interval; Fig. 1) had significantly lower PIA scores compared with same sex DA (Figs. 2D and 3B), with reductions in median ASI of 69 and 76%, respectively (median ASI: male DA = 127, male DA.F344(Cia5e) = 39, p
0.001; female DA = 114, female DA.F344(Cia5e) = 27, p
0.001). Rats subcongenic for the most centromeric 41-Mb segment (47.3 Mb including the regions where recombinations took place) of Cia5, DA.F344(Cia5d) (Fig. 1), were also protected from PIA (Fig. 2C), with reductions in median ASI of 46% for males and 55% for females, compared with same sex DA controls (Fig. 3C) (Median ASI: male DA = 127, male DA.F344(Cia5d) = 69, p = 0.018; female DA = 114, female DA.F344(Cia5d) = 51, p = 0.024). These findings suggested that the same gene accounts for Cia5d and Cia5e, and that this gene is located within the interval shared by these two subcongenics.
|
0.01), DA.F344(Cia5a) (males and females, p
0.01), and DA.F344(Cia5e) (males, p = 0.015; females p = 0.007). DA.F344(Cia5a) and DA.F344(Cia5b) rats overlap part of their congenic interval (Fig. 1), and based on the phenotypic differences we were able to exclude that common region as the one containing the arthritis gene, thus narrowing the Cia5a interval to 10.1 Mb. Taken together, these data demonstrate that Cia5 contains two PIA-regulatory genes that influence disease expression both in males and females. One gene is located within the telomeric portion of the Cia5a interval that does not overlap Cia5b and the other is within the Cia5d and Cia5e shared centromeric interval. DA.F344(Cia5) congenics and DA.F344(Cia5a) subcongenics developed a significantly milder form of CIA
To determine whether the same intervals regulating PIA also regulated CIA, Cia5 congenics and subcongenics were studied in this model. Both DA.F344(Cia5) males and females tested for CIA were protected (Fig. 2F) and had lower median ASI (36 and 53%, respectively) than same sex DA rats (median ASI: male DA = 250, male DA.F344(Cia5) = 160, p
0.01; female DA = 330, female DA.F344(Cia5) = 154.5, p
0.01). DA.F344(Cia5a) males and females also displayed significantly lower median ASI than DA (92 and 98%, respectively) (Fig. 2G) (median ASI: male DA = 250, male DA.F344(Cia5a) = 20, p
0.001; female DA = 330, female DA.F344(Cia5a) = 5, p
0.001). Although the ASI of DA.F344(Cia5) congenics were significantly lower than those of DA, they were significantly higher than those of DA.F344(Cia5a) both in males and females (p
0.004), suggesting that F344-derived alleles at a novel locus centromeric from Cia5a but within Cia5 contribute to increased disease severity. The ASI of DA.F344(Cia5b), DA.F344(Cia5e), and DA.F344(Cia5d) subcongenics were not significantly different than same sex DA rats (Fig. 2, HJ). Taken together, these data show that Cia5 contains two distinct CIA-regulatory genes, one contained in the Cia5a interval also regulating PIA, and in which the arthritis severity favoring allele originates from DA. The other gene is contained within the Cia5e interval, centromeric from Cia5a, but unlike the Cia5d PIA-regulatory effect, the CIA regulatory allele originates from F344, suggesting that this is a different gene. DA.F344(Cia5d) rats were not protected in CIA. No gender-specific effect or differences were observed.
Histologic studies in PIA synovial tissues
Histologic findings did not differ between males and females within each strain. Therefore, data from both genders were combined for analyses. Joint histologic architecture was highly abnormal in DA rats with PIA (Table I; Fig. 4A), with pronounced synovial hyperplasia and pannus formation (Fig. 4C), increased number of synovial vessels (angiogenesis) (Fig. 4D), synovial infiltration with mononuclear cells (Fig. 4E), and extensive cartilage and bone erosive changes (Figs. 4, A and G), and cartilage loss of proteoglycans (Fig. 4G). However, both DA.F344(Cia5) (Fig. 4B), DA.F344(Cia5a), and DA.F344(Cia5d) (Fig. 4F) preserved a nearly normal joint histologic architecture (Table I). No synovial hyperplasia was observed in DA.F344(Cia5), DA.F344(Cia5a), and only mild focal synovial hyperplasia was observed in a few of the DA.F344(Cia5d) ankles (Fig. 4F). DA.F344(Cia5), DA.F344(Cia5a), and DA.F344(Cia5d) had a significantly lower number of synovial vessels and synovial-infiltrating mononuclear cells, and diminished cartilage and bone erosive changes (Fig. 4, B, F, and H). These histologic findings further support the observation of reduced clinical arthritis severity, and suggest that the Cia5a and Cia5d genes regulate processes central to the pathogenesis of synovial hyperplasia, pannus formation, and joint destruction in arthritis.
|
|
mRNA were significantly reduced in DA.F344(Cia5) and DA.F344(Cia5a), compared with DA rats
Quantitative RT-PCR analysis revealed a reduction of >50% in the levels of IL-1
in the synovial tissues of DA.F344(Cia5) and DA.F344(Cia5a) rats obtained at day 32, compared with DA (
Ct: DA = 8.09, Cia5 = 9.33 (p
0.01), Cia5a = 9.34 (p = 0.02), Cia5d = 9.15 (p = 0.14); the
Ct is inversely correlated with the number of copies of mRNA in the tissue) (Fig. 5). IL-1
levels in Cia5d also tended to be lower than DA, however that difference did not reach statistical significance. TNF-
expression showed a similar trend toward reduced expression in Cia5 congenics, but lower than 50% (
Ct: DA = 10.58, Cia5 = 11.43 (p = 0.09), Cia5a = 11.28 (p = 0.20), Cia5d = 11.13 (p = 0.46)). The mild differences in TNF-
expression between congenics and DA may be attributed to the use of day 32 tissues, a time point beyond the peak of synovial TNF-
expression.
|
High levels of IgG Abs against BII were detected in all strains immunized with BII/IFA, confirming appropriate immunization (Fig. 6A). There was good concordance between median levels of anti-BII and median levels of autoantibodies anti-RII in each strain. Levels of anti-BII and anti-RII were significantly lower in Cia5a congenics compared both with DA (p < 0.05) and with all the other tested congenics, particularly in the case of anti-RII (Fig. 6). Levels of anti-BII and anti-RII tended to be higher in Cia5 compared with all other strains, including DA. Levels of anti-BII and anti-RII Abs in DA rats were not significantly different from those in Cia5, Cia5b, Cia5d, and Cia5e. The concordance between arthritis severity scores and Ab levels in Cia5a suggests that the magnitude of the Ab response against the immunogen or the self-Ag may be involved in the regulation of disease severity. Additionally, the trend toward having higher levels of anti-collagen Abs in Cia5 rats, compared with both Cia5a and DA, further supports the concept outlined above that F344-derived alleles centromeric from Cia5a contribute to increasing disease severity. Interestingly, Cia5e subcongenics, which contain F344 alleles through the whole Cia5 interval centromeric from Cia5a, did not produce anti-RII levels similar to Cia5, suggesting that the interaction (epistasis) between F344 alleles at Cia5a and at Cia5e is required for the increased levels of anti-RII.
|
| Discussion |
|---|
|
|
|---|
To identify arthritis severity regulatory genes, we have conducted genetic analyses in rat intercrosses generated between arthritis-susceptible and -resistant strains studied in different models of autoimmune arthritis (19, 20, 21, 22, 30). The chromosome 10 QTL Cia5 was originally identified in a DA x F344 F2 intercross studied for CIA (20), and its arthritis regulatory effect confirmed in DA.F344(Cia5) congenics (32). Our previous study focused on clinical evaluation of arthritis severity and suggested the possibility that the Cia5 congenic interval contains at least two CIA and PIA regulatory loci, one of which is located within the Cia5a interval (33). In the present study, we report the generation of Cia5 subcongenic strains and the identification and localization of two novel arthritis-regulatory loci, Cia5d, a PIA-regulatory locus within the most centromeric 41.5-Mb interval (47 Mb including the recombination regions) of Cia5, and another locus contained within the Cia5e interval, centromeric from Cia5a, in which F344-derived alleles specifically increase disease severity in CIA. We also confirmed the arthritis-regulatory effect of Cia5a and reduced the interval within which it is located to 10 Mb.
In addition to evaluating clinical arthritis severity during a 31-day observation period, comprehensive histologic analyses that evaluated several different parameters and the quantification of synovial IL-1
and TNF-
mRNA were conducted. These analyses revealed that, in addition to reducing PIA clinical arthritis severity, the presence of F344 alleles at either Cia5a or Cia5d was associated with preservation of the joint architecture, significant reduction in synovial hyperplasia and pannus formation, reduction of the synovial inflammatory infiltration, and number of synovial vessels (angiogenesis). Additionally, Cia5a and Cia5d congenic strains had significant reduction in cartilage and bone erosive changes and both had reduced synovial levels of IL-1
message even when tissues were obtained at a later time point during disease course (day 32), compared with DA. These findings demonstrate that the regulatory effects of Cia5a or Cia5d persist beyond the early stages of disease and into the chronic/perpetuation stages. The synovial events affected by the presence of the Cia5a or Cia5d intervals have been previously shown to be at least in part interdependent. Specifically, angiogenesis is a critical event in the pathogenesis of autoimmune arthritis and in the development of synovial hyperplasia (47, 48). Similarly, synovial hyperplasia and pannus formation are typically associated with cartilage and bone erosions in arthritis (41, 49), and are partially dependent on the inflammatory mononuclear cellular infiltrate and the cytokines these infiltrating cells produce (50, 51, 52, 53). Moreover, several of the cytokines produced by the synovial tissues and infiltrating mononuclear cells have angiogenic properties (54). Therefore, while we cannot determine at this point which specific cellular and molecular events the Cia5a and Cia5d genes regulate, our observations demonstrate that these genes control fundamental events required for the development of synovial inflammation, arthritis, and articular damage. It is conceivable that these two genes will be key regulators of arthritis pathophysiologic events, similar to what has been described with NF-
B (50).
The presence of F344 alleles within the Cia5a interval was associated with significant reduction in arthritis scores both in PIA and CIA, suggesting that the arthritis gene contained within this interval regulates a cellular and molecular process critical for both models. Cia5a congenics had lower levels of anti-BII and autoantibodies to anti-RII, raising the possibility that this gene influences the arthritis phenotype via the regulation of B cell functions. Furthermore, our observation that the protective effects of the Cia5a F344 alleles were similar in DA.F344(Cia5a) heterozygotes and homozygotes indicates that the Cia5a DA allele acts recessively to confer increased arthritis severity. Knowledge of the mode of inheritance of the Cia5a gene will greatly facilitate the identification effort because heterozygous recombinant backcross offspring can be used in arthritis studies without having to further intercross recombinants to generate homozygosity. Cia5d congenic rats were only protected in PIA, and not in CIA, suggesting that this gene regulates a pathway critical for PIA and not for CIA. This observation is in agreement with the original localization of the Cia5 locus in the DA x F344 F2 genome-wide screen (20). However, a genome-wide screen done in a BB x BN F2 intercross studied for CIA identified a Cia5d-colocalizing locus, Cia16, further suggesting that this interval contains an arthritis-regulatory gene that is not limited to a single disease model. Moreover, one cannot completely exclude the possibility that Cia5d could interact with yet another locus in an epistatic manner in the regulation of CIA in DA x F344 congenics. Because Cia5d directly and significantly regulated PIA, this will be the model of choice in the process of narrowing down the critical region toward the identification of the specific gene.
Unexpectedly, Cia5 had significantly higher CIA arthritis severity scores than Cia5a, suggesting the possibility of a third regulatory locus located in the region centromeric from Cia5a (within Cia5e), and containing an F344-derived arthritis severity-favoring allele. This possibility was further supported by the higher levels of anti-BII and anti-RII observed in Cia5 congenics as compared with Cia5a. However, the arthritis severity and Ab levels in Cia5e were not different from DA, suggesting that the CIA arthritis-favoring effect associated with F344 alleles at Cia5e depend on the presence of F344 alleles at Cia5a (epistasis). Another possibility is that in the absence of F344 alleles at Cia5a the arthritis severity is so close to maximum that it leaves little room for the detection of any additional contributions to severity. To confirm and better resolve the location of this putative locus, Cia5e and Cia5d subcongenics will be bred with Cia5a to generate double congenics in an attempt to recreate the Cia5 phenotypes.
Our previous studies had suggested that Cia5a operated in a female-specific manner to regulate CIA severity (32). In the present study, both male and female DA.F344(Cia5a) congenics were significantly protected and developed a milder form of disease. Three possible explanations could account for these differences. First, the present study was conducted at the NSLIJRI, while the previous study was conducted at the NIH. Environmental differences between the two institutions, including the rat intestinal flora, could have influenced immune responses. Second, there were certain differences in reagents that could have affected the resulting immune responses. Specifically, the source of type II collagen was different, and while the same source of IFA was used, it is conceivable that the adjuvant effect could vary from the lot used at the NIH to the one used at the NSLIJRI. The male and female arthritis protection seen in the present study in association with the Cia5a interval in fact suggests a potentially broader relevance of our studies, because identification of the gene could benefit even more people than if it operated in a sex-specific manner.
IL-1
and TNF-
are two key proinflammatory cytokines involved in the pathogenesis of rodent arthritis and RA. In the rat synovium, levels of IL-1
and TNF-
protein correlate well with their respective amount of mRNA (55). Therefore, we used quantitative real-time PCR to study the changes in levels of these two cytokines in rats studied for PIA, and determined that IL-1
levels were reduced in DA.F344(Cia5), as well as in its subcongenic DA.F344(Cia5a) (and to a lesser degree Cia5d) rats compared with DA rats, suggesting that the genes accounting for these QTLs influence inflammatory processes upstream from this cytokine. Although DA.F344(Cia5) TNF-
mRNA levels tended to be lower than in DA, the difference did not reach statistical significance. One likely explanation is the timing of the synovial tissue collection. Specifically, we obtained synovial tissues after the completion of the arthritis observation period (day 32), a time point that was beyond the peak of the synovial levels of TNF-
in rats (55, 56). Interestingly, even at this late stage differences in cytokine expression could still be detected, underscoring the regulatory effect of Cia5a and Cia5d genes during later stages of the disease course.
The three arthritis regulatory intervals studied herein colocalize with loci regulating different forms of autoimmune arthritis in rats (19, 27, 28, 32, 33, 34, 57, 58). Additionally, Cia5d colocalizes with a murine arthritis locus Pgia7 (59), and most importantly, it contains SLC22A4 (human chromosome 5q23.3) (11), a gene recently implicated in the susceptibility to RA in Japan. Although candidate gene analyses are typically fruitless when dealing with a genomic interval as long as Cia5d, the rat homologue of SLC22A4, Octn1, was a compelling candidate. Preliminary cDNA sequencing analyses did not reveal any differences between the coding regions of DA and F344 Octn1 (M. Brenner and P. Gulko, unpublished observations). Additional candidate genes of interest in the Cia5d interval include IL-3, IL-4, IL-5, and TIM-1.
The Cia5a syntenic region on human chromosome 17q22-q25 contains two RA susceptibility loci recently identified in genome-wide scans (4, 38). Taken together, the colocalization of multiple arthritis regulatory loci and genes to the same, or syntenic, interval across different species, including humans, suggests that the findings from rat studies will be relevant to RA and may identify novel genes for focused case-control testing.
Cia5a and Cia5d also colocalize with several loci, or their syntenic regions, implicated in the regulation of rodent models of autoimmunity and inflammation (35, 36, 37, 60, 61, 62), or in the susceptibility to human autoimmune diseases (63, 64, 65, 66). These studies suggest that both Cia5a and Cia5d genes could be relevant not only to arthritis, but also to other autoimmune diseases. The concept of colocalization of non-MHC loci regulating different forms of autoimmune diseases has been previously proposed (67), and recently confirmed with the discovery that polymorphisms in CTLA4 (68) and PTPN22 (12, 69) regulate more than one type of autoimmune disease. We consider that Cia5a and Cia5d could become additional examples of such autoimmune regulatory genes common to more than one autoimmune disease.
In conclusion, the present study described two novel arthritis-regulatory locus on rat chromosome 10, a PIA-specific Cia5d and a CIA locus where F344-derived alleles within the Cia5e interval contribute to increased disease severity. Additionally, this study narrowed the previously reported locus Cia5a. Both Cia5a and Cia5d contain genes that regulate arthritis severity, synovial hyperplasia, pannus formation, synovial inflammation and angiogenesis, and joint damage, processes central to disease pathogenesis underscoring the potential relevance of the gene identification. We are in the process of reducing the critical regions containing Cia5a and Cia5d to
1 Mb for positional cloning of the specific genes. Once identified these genes will generate novel candidate genes or candidate pathways for focused RA case-control studies, as well as novel targets for the development of new therapies and tools for diagnosis and prognostication.
| Acknowledgments |
|---|
| Disclosures |
|---|
|
|
|---|
| Footnotes |
|---|
1 This work was funded by Grants R01-AR46213 and R01-AI54348 from the NIAMS and the National Institute of Allergy and Infectious Diseases (NIAID, NIH; to P.S.G.), by the Department of Veterans Affairs Biomedical Research Funds, and the Nora Eccles Treadwell Foundation (to M.M.G.). ![]()
2 Current address: Medical College of Ohio, Toledo, OH 43614. ![]()
3 Current address: Genomics and Genetic Branch, NIAMS, NIH, Bethesda, MD 20892. ![]()
4 Current address: MedImmune, Gaithersburg, MD 20878. ![]()
5 Address correspondence and reprint requests to Dr. Pércio S. Gulko, Laboratory of Experimental Rheumatology, Robert S. Boas Center for Genomics and Human Genetics, North Shore-Long Island Jewish Research Institute, 350 Community Drive, Room 139, Manhasset, NY, 11030. E-mail address: pgulko{at}nshs.edu ![]()
6 Abbreviations used in this paper: RA, rheumatoid arthritis; QTL, quantitative trait locus; CIA, collagen-induced arthritis; PIA, pristane-induced arthritis; BII, bovine type II collagen; RII, rat type II collagen; ASI, arthritis severity index; HMF, high-power magnification field; Ct, cycle threshold. ![]()
Received for publication December 8, 2004. Accepted for publication April 12, 2005.
| References |
|---|
|
|
|---|
BL as the second major histocompatibility complex-linked susceptibility locus for rheumatoid arthritis. Am. J. Hum. Genet. 72: 303-312.[Medline]
, interleukin-1 receptor antagonist, interleukin-4, and interleukin-10 gene polymorphisms: relationship to occurrence and severity of rheumatoid arthritis. Arthritis Rheum. 42: 1093-1100.[Medline]

CT) method. Methods 25: 402-408.[Medline]
v
3 antagonist. J. Clin. Invest. 103: 47-54.[Medline]
B activation provides the potential link between inflammation and hyperplasia in the arthritic joint. Proc. Natl. Acad. Sci. USA 95: 13859-13864.
1 (TGF-
1) induced neutrophil recruitment to synovial tissues: implications for TGF-
-driven synovial inflammation and hyperplasia. J. Exp. Med. 173: 1121-1132.This article has been cited by other articles:
![]() |
L Backdahl, J P Guo, M Jagodic, K Becanovic, B Ding, T Olsson, and J C Lorentzen Definition of arthritis candidate risk genes by combining rat linkage-mapping results with human case-control association data Ann Rheum Dis, December 1, 2009; 68(12): 1925 - 1932. [Abstract] [Full Text] [PDF] |
||||
![]() |
Q. Xiong, Y. Jiao, K. A. Hasty, J. M. Stuart, A. Postlethwaite, A. H. Kang, and W. Gu Genetic and Molecular Basis of Quantitative Trait Loci of Arthritis in Rat: Genes and Polymorphisms J. Immunol., July 15, 2008; 181(2): 859 - 864. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Brenner, T. Laragione, A. Mello, and P. S Gulko Cia25 on rat chromosome 12 regulates severity of autoimmune arthritis induced with pristane and with collagen Ann Rheum Dis, July 1, 2007; 66(7): 952 - 957. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Laragione, N. C. Yarlett, M. Brenner, A. Mello, B. Sherry, E. J. Miller, C. N. Metz, and P. S. Gulko The Arthritis Severity Quantitative Trait Loci Cia4 and Cia6 Regulate Neutrophil Migration into Inflammatory Sites and Levels of TNF-{alpha} and Nitric Oxide J. Immunol., February 15, 2007; 178(4): 2344 - 2351. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Joe Quest for arthritis-causative genetic factors in the rat Physiol Genomics, January 12, 2007; 27(1): 1 - 11. [Abstract] [Full Text] [PDF] |
||||
![]() |
M Brenner, C Braun, M Oster, and P S Gulko Thermal signature analysis as a novel method for evaluating inflammatory arthritis activity Ann Rheum Dis, March 1, 2006; 65(3): 306 - 311. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |