|
|
||||||||

* Medical Inflammation Research, Lund University, Lund, Sweden; and
Department of Laboratory Medicine, Section of Clinical Chemistry, Lund University, Malmö University Hospital, Malmö, Sweden
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Rheumatoid arthritis (RA)4 is one such common and complex disease. RA is a chronic inflammatory disease affecting
1% of the human population. RA has a clear genetic predisposition with heritability estimated to
60% (2). The genetic influence is to a significant degree associated with genes in the MHC. Efforts have been made to identify non-MHC susceptibility genes in human RA (exemplified by Refs. 3 , 4).
A well-established strategy to explore the genetic causative factors and molecular mechanism for common disease is to use inbred animal models. There are a number of different animal models for RA (reviewed in Ref. 5). The most widely used model is collagen-induced arthritis (CIA) (reviewed in Ref. 6), which meets most of the criteria for diagnosis of RA. Both humoral and cellular immune mechanisms are involved in the induction and mediation of CIA. Various genetic linkage studies have until today located >30 genomic loci potentially harboring genes predisposing to CIA and other models for human RA in mice and rats. Subsequent isolation of the loci in congenic strains has provided fundamentally new knowledge on the genetic control of these diseases. The strategy has successfully identified two important genes regulating development of arthritis in mice, the MHC class II gene Aq controlling immune response (7) and the Ncf1 gene controlling oxidative burst (8, 9). Both of these genes are associated with T cell activation and arthritis severity. However, in many cases, the identification of the actual disease causing or modifying gene within the susceptibility loci has turned out to be more difficult than expected. The phenotype determined by a susceptibility locus often depends on the genetic background in which the locus is introduced, suggesting that genetic modifiers interact with the locus and thereby have a role in determining the functional expression of genetic variation (10).
Using animal models, the new strategy of partial advanced intercross (PAI) have recently been developed to investigate genetic interactions by repetitively crossing two congenic strains, both made on the same genetic background (11, 12). This strategy was used in the present study to investigate the interaction between the two arthritis susceptibility loci Cia2 and Cia9 in mice. We have previously identified the Cia2 and Cia9 loci as the main genetic factors influencing the arthritis development in a F2 intercross between the arthritis susceptible strain C57BL/10.Q (B10.Q) and the resistant NOD.Q strain (13). The Cia2 locus associated with CIA was also identified in (DBA/1 x SWR/J) F2 crosses (14, 15) and associated with arthritis in the KXB/N arthritis model in a cross between B6 and NOD (16). The Cia2 locus harbors a strong candidate gene for arthritis susceptibility, the complement component C5-encoding gene Hc. Both the NOD and the SWR/J strains carry a 2-bp deletion in an exon near the 5' end of the Hc gene resulting in frame-shift and termination resulting in C5 protein deficiency (17, 18). The complement C5 plays multiple roles in the process of inflammation. Upon activation of C5, C5a fragment is released and acts as chemoattractant for neutrophils, monocytes, macrophages, and eosinophils, induces production of cytokines and other proinflammatory mediators, and enhances vascular permeability (19, 20). The other product of C5 activation, C5b, initiates assembly of membrane attack complex (MAC), which can further contribute to inflammation due to its ability to lyse cells and to induce drastic changes in cells when present at sublytic concentrations. The involvement of C5 in the pathogenic mechanisms of arthritis has been demonstrated in several reports (21, 22, 23, 24). It has even been suggested that C5 sufficiency is an absolute requirement for CIA to develop (14, 22). The hypothesis that the Cia2 locus was almost exclusively mediated by a deficient C5 gene was, however, not compatible with some of the reported findings. First, mice with NOD-derived Cia2 could develop arthritis, although at low incidence rate, and second, mice with Cia2 heterozygosity, with undisturbed C5 levels in the serum, had a dramatically reduced arthritis severity (13). In addition, the Cia2 loci seemed to influence the functional expression of the Cia9 locus on chromosome 1, indicating an interactive influence. To resolve these questions, we made an extensive backcross experiment to investigate the possibility that there are other dominant loci that could interact. Furthermore, we established congenic strains of the Cia2 and Cia9 loci and made a PAI between them to pinpoint the interactive effects. Finally, we analyzed the role of Cia2 in modulating complement functions to provide an explanation for a possible C5-mediated effect.
| Materials and Methods |
|---|
|
|
|---|
Mice were bred and kept in the animal house of Section for Medical Inflammation Research, Lund University. The mice were kept in a climate-controlled environment with 12-h light/dark cycles and fed with standard rodent chow and water ad libitum (as defined at
www.inflam.lu.se
). The congenic strain for the Cia2, B10.Q.NOD-Cia2, was constructed by the speed congenic technique (25) as described previously (13). For this study, the congenic strain had been backcrossed to B10.Q strain for additional 5 generations, yielding in total 10 generations of backcrossing. The congenic fragment from NOD.Q was
54-Mb long, ranging from D2Mit116 to D2Mit91. The congenic strain for the Cia9 locus (B10.Q.NOD-Cia9) was constructed using the same technique as the Cia2 congenic (25). The congenic fragment from NOD.Q was
150-Mb long, ranging from D1Mit235 to D1Mit17 (see Fig. 6).
|
Genotyping
Microsatellite markers were purchased from Interactiva Biotechnology. The order of the markers was based on the map available from the Jackson laboratory (
www.jax.org
).
Backcross.
A full genomic scan was performed on 395 backcross animals using 162 microsatellite markers. The complete list of markers is available on
www.inflam.lu.se
.
All markers were assayed by PCR as follows: Genomic DNA (15 ng) was amplified in a final volume of 9 µl containing dNTP (200 µM), MgCl2 (1.5 mM), primer (1.5 pmol of each), and TaqGold DNA polymerase (0.2 U) (PerkinElmer Biosystems). The forward primer was labeled with a fluorescent dye. Amplification conditions were as follows: 95°C for 1.1 min, followed by 33 cycles of 95°C for 30 s, 55°C for 75 s, 72°C for 75 s, and a final extension at 72°C for 7 min. The reactions were performed using an ABI 877 Integrated Thermal cycler (Applied Biosystems), followed by an automatic pooling of the products. The PCR products were size-fractionated on 4% polyacrylamide gels, on an ABI 377 sequencer, and the size of the fragments were determined using the Genescan software version 3.1 (Applied Biosystems) with TAMRA GS-350 or 500 as internal size standard.
Partial advanced intercross.
DNA was prepared from tail biopsies by a quick alkaline lysis protocol (26). Briefly, tail biopsies were incubated in 50 mM NaOH for 2 h at 95°C, vortexed, and neutralized in 1 M Tris-HCl (pH 8.0). After centrifugation, the supernatant was used for PCR. Genotypes from 502 PAI mice were analyzed from 11 markers spanning 46.0189.3 Mb on chromosome 1 and 6 markers from 12.0 to 38.2 Mb on chromosome 2 (see
www.inflam.lu.se
). Genotypes from 162 microsatellite markers were analyzed.
PCR was performed with 10 ng of DNA in a reaction volume of 10 µl containing the following: 10 mM Tris-HCl (pH 9.0), 50 mM KCl, 1.5 mM MgCl2, 0.3 µM forward and reverse primer (MWG-Biotec and Applied Biosystems), 100 µM dNTP (Amersham Biosciences), and 0.25 U Taq (Amersham Biosciences). Forward primers were labeled with fluorescent dyes. PCR was performed in a thermal cycler (MJ Research) under the following amplification conditions: denaturation at 94°C for 2.5 min, annealing at 56°C for 45 s, polymerization at 72°C for 1 min, followed by 30 cycles of 94°C for 30 s, 56°C for 45 s, and 72°C for 1 min. The final cycle ended by elongation at 72°C for 6 min. PCR products were analyzed on a MegaBACE 1000 (Amersham Biosciences), according to the manufacturers protocol.
Induction and evaluation of CIA
At the age of 1012 wk, mice were immunized at the base of the tail with 100 µg of rat collagen type II (CII) emulsified 1:1 in CFA (Difco) in a volume of 100 µl. At 35 days postimmunization, the mice were given a boost injection of 50 µg of CII emulsified in IFA (Difco) in 100 µl. Arthritis was evaluated blindly using a scoring system, as described previously (27), based on the number of inflamed joints in each paw. Briefly, each toe and knuckle was given a score of 1 (total 10/limb), and inflamed wrist of ankle was given a score of 5, giving a maximum score per limb of 15 and a maximal score per mouse of 60. Severity of arthritis was evaluated as maximum disease score and as accumulated score (area under the curve), generated by calculating the sum of the disease scores for each mouse, respectively.
Measurement of anti-CII IgG levels
Mice were bled at day 35 after induction of CIA, and the level of mouse anti-CII IgG Abs (µg/ml) were measured in the sera, using ELISA as previously described (28).
Detection of complement factor 5 in sera
The levels of complement factor 5 in sera were detected using the ELISA as described previously (29). Sera (from day 35 after induction) were added to 96-well plates (Costar, Corning) coated with a chicken anti-mouse C5 Ab. To detect C5 levels, a biotinylated chicken anti-mouse Ab was used followed by a streptavidin alkaline phosphatase (Jackson ImmunoResearch Laboratories) incubation. Paranitrophenol was used as a chromogenic substrate, and the absorbance was determined in a Titertek multiscan spectrophotometer.
Abs for flow cytometry
For staining of surface epitopes, the following mAbs were used: anti-CD4 (H129.19) and anti-CD8 (53.6.7), purchased from BD Pharmingen. Anti-B220 (RA3-6B2), anti-CD4 (CT-CD4), and anti-CD8 (CT-CD8) (all three from Caltag Laboratories), anti-CD11B (M1/70), anti-Fc
RII/III (2.4.G2.), and anti-CD3
(145-2C11) were purified from culture supernatants by affinity chromatography on protein G-Sepharose and conjugated.
Tissue preparation and in vitro flow cytometry
At the end of the experiment spleens were removed, and cell suspension was prepared in PBS by mechanically mashing the tissue. Cells were washed and resuspended in staining buffer containing 0.5% BSA (Sigma-Aldrich) and 0.01% NaN3 in PBS. A total of 12 x 106 cells was stained for various surface epitopes for 20 min at 4°C. An anti-Fc
RII/III Ab was used to inhibit Ab binding to the Fc
R. A typical forward and side scatter gate for activated lymphocytes was set to exclude dead cells and aggregates. Typically, 104 events from the lymphocyte gate were collected and analyzed using a FACSort (BD Biosciences) and BD Biosciences software. Quadrant and histogram statistics were placed on the basis of the staining of the negative controls. Less than 0.5% positively stained cells was not regarded as significant.
Complement assay
Activity of the alternative pathway was studied using hemolytic assay for measurement of the alternative pathway including activity of C5. Rabbit erythrocytes were washed twice in Mg-EGTA buffer (2.5 mM veronal buffer (pH 7.3), 70 mM NaCl, 140 mM glucose, 0.1% gelatin, 7 mM MgCl2, and 10 mM EGTA) and resuspended at 0.5 x 109 cells/ml. Mouse serum (5 µl) was incubated with 25 µl of erythrocyte suspension in total volume of 100 µl of Mg-EGTA buffer. After 4 h of incubation at 37°C, the degree of lysis was determined by spectrophotometric measurement of the amount of released hemoglobin at 405 nm. For each sample, an absorbance provided by serum itself was subtracted, and the values are shown as Z, which equals the natural negative logarithm of 1 percentage of lysis.
Statistical and linkage analysis
Linkage analysis was performed using R/qtl QTL mapping software (30) in R (R Development Core Team, 2004). The genome scan was performed using scanone (imputation model) (31) under the assumption that sex influences the arthritis development as an interactive covariate. The chromosome X was omitted. Arthritis severity was analyzed as maximum arthritic score and accumulted arthritic score (area under the curve). QTL analysis was calculated under the assumption that gender influenced the arthritis development as an interactive covariate. To investigate gene interactions, two-dimensional genome scan with a two-QTL model was performed. In the analysis, LOD scores are calculated testing the full model, in the presence of covariates, to each of two alternatives: The joint LOD (LODjnt (joint)) score compares the full model to the null model, and the epistasis LOD score (LODint) compares the full model to the additive model. The full model is as follows: (y = m + b[q1] + b[q2] + b[q1 x q2] + Ag + Zd[q1] + Zd[q2] + Zd[q1 x q2] + e), where q1 and q2 are the unknown genotypes at two locations, A is a matrix of covariates, and Z is a matrix of QTL-interacting covariates. The null model is as y = m + Ag + e and the additional model is as y = m + b[q1] + b[q2] + Ag + Zd[q1] + Zd[q2] + e. Empirical significance levels (70, 95, and 99%) were established by permutation tests (1000 permutations) (32) for each phenotype.
The marker map covered the NOD.Q x B10.Q backcross genome at an
20 cM intermarker distance.
Other statistical analyses were performed using Mann-Whitney U test for comparisons between groups. Differences were regarded as significant if the p value was <0.05.
| Results |
|---|
|
|
|---|
We have previously identified Cia9 and Cia2 as loci regulating arthritis severity in a F2 intercross between the B10.Q and NOD.Q strains (13). Because we were surprised that we did not find any other loci, in particular loci earlier identified to be associated with diabetes using crosses with the same genetic backgrounds, we designed a powerful backcross experiment. We produced 395 (B10.Q x NOD.Q) x B10.Q backcross mice in an effort to find additional modifying loci and to define the relative importance of Cia2 and Cia9. As expected, the arthritis developed variably from none to a severe acute disease that in most affected animals resulted in joint deformation. Genotypes from 162 markers genome-wide were analyzed, and a genetic linkage analysis was performed. Again, we found that the Cia2 and Cia9 were the strongest loci affecting the arthritis development, as in the earlier F2 intercross (Table I). To our surprise, only one additional locus was identified, a locus in the centromeric part of chromosome 10 (7.0 cM), that influenced arthritis severity (LOD = 2.94) (Table I). No other loci were detected in the present backcross. We conclude that in crosses between the NOD and C57BL/10 genomes, the polygenic segregation is mainly dependent on only two dominant loci, the Cia9 locus promoting arthritis and Cia2 protecting when originating from the NOD strain.
|
We could not identify any loci linked to immune response against CII measured as levels of anti-CII Abs in the present backcross (data not shown). Neither did we detect any difference in anti-CII Ab production between B10.Q.NOD-Cia2 congenic mice and wild-type B10.Q (data not shown). These data confirmed our previous finding of no loci linkage to anti-collagen Ab production in the F2 intercross between B10.Q and NOD.Q (13). Instead, to further dissect the mechanism of arthritis in the backcross, we characterized the mice regarding distribution of immune and inflammatory cells in blood. Frequency distribution was determined by flow cytometry on spleen cells taken at the end point of the experiment. The frequency of the different cell types were used as quantitative traits in linkage analysis to identify loci determining the cell distributions. We identified a locus on chromosome 12 involved in the distribution of both T and B cells (Table II). The NOD allele at D12Mit76 reduced the frequency of B cells and thereby the ratio of B cells:T cells (data not shown). The chromosome 12 locus also influenced the ratio of CD4-positive to CD8-positive T cells by the NOD allele increasing the ratio. The Cia2 locus influenced the frequency of CD11b-positive cells (Table II), and the higher expression of CD11B-positive cells correlated with arthritis development (p = 0.011) (Fig. 1). Moreover, the arthritis-suggestive locus on chromosome 10 was associated with a decreased frequency of CD11B-positive cells in blood. However, the Cia9 locus did not influence the frequency or distribution of any of the investigated cells.
|
|
In the backcross, the linkage to Cia2 was significantly associated with reduction of arthritis development. If the mice were homozygous B10.Q at the Cia9 locus, the effect of B10.Q/NOD.Q heterozygosity at Cia2, mediated almost complete inhibition of arthritis development (thick line, Fig. 2). However, introduction of one NOD allele at Cia9 partially overcame the inhibitory effect of Cia2, enabling the development of arthritis (thin line, Fig. 2). The interaction between the two loci was highly significant (LOD score of 16.1 (LODjnt)).
|
To confirm the linkage to Cia2 and Cia9, the fragments were introgressed on B10.Q, and subsequently the congenic strains were used to establish a PAI cohort following an earlier outlined strategy (11) (Fig. 3). Using this strategy, we obtained mice with combinations of the genotypes of the two loci to determine differences in arthritis incidence and severity, without influences from the rest of the genome. The dominant protective effect of Cia2 was clearly confirmed in the PAI, seen by almost complete inhibition of arthritis in mice carrying one or two NOD.Q alleles at the Cia2 locus (Table III and Fig. 4B,
and
![]()
). The Cia9 by itself was influencing the disease by increasing the arthritis severity in the mice homozygous and heterozygous for NOD alleles at the peak marker (D1Mit270) of Cia9, located close to the FcR region (Fig. 4B, and
) as compared with B10.Q wild-type mice (Fig. 4B,
). Interestingly, the inheritance was different at the more centromeric part of the congenic region on chromosome 1 (at marker D1Mit14) where heterozygosity enhanced the disease, whereas homozygosity of the NOD allele in fact protects against arthritis (Fig. 4A). This argues for the presence of at least two genes in the Cia9 fragment.
|
|
|
Using the B10.Q.NOD-Cia2 congenic mice, we could confirm the observation from the backcross that the Cia2 locus completely suppresses arthritis development in a dominant fashion (Fig. 5), because none of the homozygous or heterozygous Cia2 congenic mice developed CIA. Within this fragment the NOD allele of C5 (Hc) has a mutation prohibiting C5 expression, and because C5 has been suggested to be of crucial importance for arthritis susceptibility, it is the major candidate gene. The levels of C5 in blood were reduced in mice heterozygous for the NOD allele of C5, although not statistically significantly lower than in wild-type B10.Q mice (C5 sufficient) (Fig. 6). We then investigated effect of Cia2 on functional effects of complement. We analyzed the complement activity using functional hemolytic assay, dependent on the presence of C5, after arthritis induction but before disease onset as well as during active inflammation (days 0, 35, and 70). As expected, we could demonstrate no complement C5-dependent lysis in the homozygous Cia2 congenic mice (Fig. 7). Interestingly, at the time of induction (day 0) as well as day 35 after disease induction, the effect on C5-dependent lysis was significantly less pronounced in the heterozygous Cia2 congenics compared with the wild-type B10.Q (Fig. 7). This finding supports a dominant protective effect of the mutated C5 gene because this led to a decreased functional complement activity. However, it does not exclude additional modifying effects by other polymorphic genes within the Cia2 fragment.
|
|
| Discussion |
|---|
|
|
|---|
Our results indicate that the Cia2 locus in the Cia2-congenic strain inhibits arthritis in a dominant fashion, and we show that the function of C5 is disturbed in the heterozygous animals. Although the levels of C5 in serum are not significantly reduced in heterozygous mice compared with wild type, there are detectable differences using a functional hemolytic assay. We conclude that the Hc gene plays the major role in the Cia2 loci. However, it is important to emphasize that a major role for C5 in the Cia2 still does not exclude a role for additional modifier genes. Such genes could only be determined through reducing the fragment in smaller pieces in subcongenic strains.
The role of complement in the development of arthritis is perhaps underestimated as a pathogenic factor that may differ between arthritis and diabetes. Complement activation can be observed in joint spaces of patients with RA, and it correlates with severity of the disease (35). Levels of C5a are high in synovial fluid of patients with RA (36). Further support for in vivo role of complement in arthritis is provided by the fact that mice lacking various complement components such as C5, C3, or factor B due to genetic manipulations do not develop arthritis (23, 37). Furthermore, monoclonal anti-C5 Abs are able to prevent CIA in mice (24) and are now subjected to clinical test in patients with RA. In addition to blocking of induction of the disease, these Abs block progression and partially reverse joint destruction in established disease. Although complement is one of the players in the development of arthritis, there are many additional ones. It is important to delineate the role of each complement component to rationally design complement inhibitors to be used in the clinic. Although inhibitors of C3-convertase (enzymatic complex activating C3 into C3b and C3a) may protect against arthritis, they are only moderately efficient (21) or may cause more side effects than therapy aiming at inhibition of C5. For example, opsonization with C3b is the main mechanism by which complement contributes to phagocytosis of pathogenic microorganisms. The hypothesis that C5a and/or MAC are main contributors to the disease is further supported by beneficial effects of oral C5aR antagonist in treatment of arthritis in rat (38) as well as increased arthritis severity in mice lacking the MAC inhibitor, CD59 (39). Interestingly, in the PAI experiment, the level of complement activity was up-regulated at day 70, which is consistent with the fact that several complement factors are acute phase proteins up-regulated upon inflammation. Although our study shows the genetic complexity of development of arthritis, it seems clear that C5 deficiency is involved, explaining a large part of the arthritis variance and confirming C5 as an important target for anti-inflammatory therapy of arthritis.
| Acknowledgments |
|---|
| Disclosures |
|---|
|
|
|---|
| Footnotes |
|---|
1 This work was supported by grants from the Swedish Research Council, Swedish Foundation for Strategic Research (INGVAR), Greta and Johan Kocks Research Foundation, King Gustav Vs 80th Anniversary Foundation, Anna Greta Crafoord (Crafoord) Foundation, Kock Foundation, Österlund Foundation, Swedish Association against Rheumatism, Swedish Medical Research Council, Strategic Research Foundation, and the European Union projects MUGEN and AUTOCURE. ![]()
2 Current address: Cartela AB, Box 709, Lund 22007, Sweden. ![]()
3 Address correspondence and reprint requests to Dr. Rikard Holmdahl, Medical Inflammation Research, BMC, I11, Lund University, Lund 22184, Sweden. E-mail address: rikard.holmdahl{at}med.lu.se ![]()
4 Abbreviations used in this paper: RA, rheumatoid arthritis; CIA, collagen-induced arthritis; PAI, partial advanced intercross; MAC, membrane attack complex; CII, collagen type II; LOD, logarithm of the odds; jnt, joint. ![]()
Received for publication January 6, 2006. Accepted for publication June 22, 2006.
| References |
|---|
|
|
|---|
RIIb and is not associated with loci controlling diabetes. Eur. J. Immunol. 31: 1847-1856. [Medline]
-chain transgenic mice. Arthritis Rheum. 41: 256-262. [Medline]
haplotype and complement C5 play no significant role for the resistance to collagen induced arthritis in the SWR mouse. Immunology 73: 191-196. [Medline]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] |
||||
![]() |
A M Blom, K S Nandakumar, and R Holmdahl C4b-binding protein (C4BP) inhibits development of experimental arthritis in mice Ann Rheum Dis, January 1, 2009; 68(1): 136 - 142. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Boross, P. L. van Lent, J. Martin-Ramirez, J. van der Kaa, M. H. C. M. Mulder, J. W. C. Claassens, W. B. van den Berg, V. L. Arandhara, and J. S. Verbeek Destructive Arthritis in the Absence of Both Fc{gamma}RI and Fc{gamma}RIII J. Immunol., April 1, 2008; 180(7): 5083 - 5091. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |