Depletion of regulatory T cells leads to an exacerbation of delayedtype hypersensitivity arthritis in C57BL/6 mice that can be counteracted by IL-17 blockade

(2016). Depletion of regulatory T cells leads to an exacerbation of delayedtype hypersensitivity arthritis in C57BL/6 mice that can be counteracted by IL-17 blockade. Disease Models & Mechanisms , 9 , 427-440. https://doi.org/10.1242/dmm.022905


INTRODUCTION
Rodent models of arthritis have been extensively used in the process of elucidating the pathogenesis of rheumatoid arthritis (RA), a chronic inflammatory disease characterised by chronic joint inflammation and bone erosion, and they are instrumental in the development of new therapeutic strategies (Asquith et al., 2009).Delayed-type hypersensitivity arthritis (DTHA) is an arthritis model in C57BL/6 (B6) mice that is characterised by synchronised onset, 100% penetrance and low variation.DTHA is induced by modifying a protein-antigen-induced delayed-type hypersensitivity response in the paw by administration of an anti-type-II-collagen antibody cocktail (anti-CII) between immunisation and challenge.The result is severe inflammatory arthritis affecting one, predefined, hind paw (Atkinson et al., 2012).Notably, the anti-CII dose required for DTHA is five-to eight-times lower than required for induction of collagen-antibody-induced arthritis (CAIA) in B6 mice (Hutamekalin et al., 2009).
Owing to their well-established role in controlling inflammation in many animal models, we hypothesise that regulatory T cells (T regs ) play a role in the resolution of inflammation in DTHA.It has been shown in mouse models of arthritis that deficiency or depletion of T regs exacerbates disease, and increasing numbers of T regs can reduce disease activity (Gizinski and Fox, 2014;Kelchtermans et al., 2009;Miyara et al., 2014;Noack and Miossec, 2014;Zaiss et al., 2010).The development of a more sustained disease phenotype in the absence of T regs would let us study disease drivers unchecked by this immunoregulatory cell subset and to identify which disease driver mechanisms are suppressed by T regs in vivo in experimental arthritis.The purpose of the present study was therefore to investigate the mechanisms of self-limiting disease in DTHA through selective depletion of T regs after disease induction.T reg depletion studies in experimental arthritis models are scant and have mainly used the anti-CD25 approach.Anti-CD25 treatment in collagen-induced arthritis (CIA) accelerates disease (Kelchtermans et al., 2005;Morgan et al., 2003).Administration of anti-CD25 prior to induction exacerbates glucose-6-phospate isomerase (G6PI)-induced arthritis (Frey et al., 2010) and antigeninduced arthritis (AIA) (Frey et al., 2005).However, using anti-CD25 antibodies for depletion of T regs also targets effector T cells (T eff ).The Foxp3-DTR-eGFP mouse allows selective depletion of T regs without affecting T eff .This mouse expresses a fusion of a diphtheria toxin receptor (DTR) and enhanced green fluorescent protein (eGFP) under the control of the forkhead box protein 3 (Foxp3) gene locus, which means that T regs can be selectively depleted by administration of diphtheria toxin (DT) (Feuerer et al., 2009;Lahl et al., 2007).T reg depletion in Foxp3-DTR-eGFP mice on a DBA/1 background exacerbated G6PI-induced arthritis (Irmler et al., 2014).However, to our knowledge, this is the first study to address, in-depth, the use of the Foxp3-DTRe-GFP mouse for T reg depletion in experimental arthritis and the first to use the Foxp3-DTR-eGFP mouse on a B6 background in experimental arthritis.
RA has recently been associated with changes in the gut microbiota (Scher et al., 2013;Zhang et al., 2015).In spontaneous mouse models of autoimmune arthritis, joint inflammation is attenuated under germ-free conditions, but colonisation of the gut with commensal microbes is sufficient to elicit joint inflammation comparable to that observed in conventional mice (Abdollahi-Roodsaz et al., 2008;Wu et al., 2010).In these models, the colonisation of the gut resulted in a perturbed T reg /T eff balance, which was associated with disease onset and progression.The interplay between T regs and the gut microbiota modulates T reg abundance and function and might thereby also affect onset and progression of arthritis.Therefore, we also analysed the fecal microbiota following DTHA induction alone and in conjunction with T reg depletion.
In the present study, we found that depleting T regs after onset of DTHA led to an exacerbation of arthritis.Inflammatory cell infiltration, osteoclast activation and bone erosion were increased in T reg -depleted mice.T reg depletion also increased levels of antimutated-citrullinated-vimentin (MCV) antibodies, indicating an increase in autoimmunity associated with increased protein citrullination.Production of a wide range of cytokines and chemokines was increased, including IL-17, a cytokine involved in the pathogenesis of both RA and murine experimental arthritis.We show that treating T reg -depleted mice with anti-IL-17 monoclonal antibody (mAb) rescues them from exacerbation of disease through a reduction of circulating neutrophils and a reduction in IL-6 and receptor-activator of nuclear factor κB ligand (RANKL).

RESULTS
A highly activated and proliferating subset of T regs is found in the lymph node draining the arthritic paw early after disease onset We investigated the dynamics of T regs in DTHA by analysing the popliteal lymph node draining the arthritic paw (dPLN), the popliteal lymph node draining the control paw (ndPLN), paw infiltrate and blood at different time points during disease by flow cytometry (Fig. 1A).T regs were defined as live CD45 + TCRβ + CD4 + CD25 + FoxP3 + cells and T eff cells by inverting the CD25 + FoxP3 + gate within the CD45 + TCRβ + CD4 + gate (full gating strategy shown in Fig. S1A).T regs were most abundant in the dPLN on day 2 and 4 after DTHA induction and were still elevated compared to day 1 and naïve mice on day 7.The increase in T reg number in the dPLN was also reflected in blood (Fig. 1B).A dramatic increase in activated (CD44 + ) and proliferating (Ki67 + ) T regs was seen in the dPLN, but not in the ndPLN or in blood (Fig. 1B, bottom graphs).Analysis of Foxp3 mRNA showed a significantly higher Foxp3 expression in paws on day 3 and 8 compared to before arthritis induction (Fig. 1C).The same was seen with Ctla4 mRNA, which codes for the surface marker cytotoxic T-lymphocyte-associated protein 4 (CTLA4) found on T regs and memory T cells, but to a greater extent on T regs (Jago et al., 2004) (Fig. 1C).No increase in T reg numbers in blood or popliteal lymph nodes was observed during the immunisation phase of DTHA (from naïve to day 0), so we also analysed the inguinal lymph nodes (ILNs), which drain the immunisation sites, in naïve mice and immediately prior to methylated bovine serum albumin (mBSA) challenge.In these lymph nodes a dramatic increase in the frequency of both T reg and T eff cells was observed (Fig. 1D), but this increase was not reflected in the number of circulating T regs or T eff on day 0 (Fig. 1B).Interestingly, a subpopulation of CD25 hi FoxP3 hi T regs was found in the dPLN but not in the ndPLN or blood (Fig. 1E,F).These CD25 hi FoxP3 hi T regs were found to possess an activated phenotype and were actively proliferating already from day 1 post-DTHA induction (Fig. 1G).We also analysed cells isolated from the arthritic paw on day 2 and 4 after DTHA induction and found the fraction of T regs of total T cells to increase from 12.17±1.09%to 14.70±0.63%from day 2 to 4 (mean±s.e.m., P=0.079, Student's t-test).Analysis of immune cells isolated from non-inflamed paws is not feasible owing to the very low number of cells that it is possible to harvest.We also found subpopulations of CD25 hi FoxP3 hi and Ki67 + T regs in the paw.Among paw infiltrating T cells on day 2 and day 4, we found 33.72±4.16%and 23.04±1.55% of T regs to be CD25 hi FoxP3 hi (mean±s.e.m., P=0.044, Student's t-test) and 82.16±1.3% and 86.35±1.2% of T regs to be Ki67 + (mean±s.e.m., P=0.034, Student's t-test), respectively.In summary, we observed increased numbers of activated and proliferating T regs in the dPLN and paw early after mBSA challenge in DTHA.

Depletion of T regs after arthritis induction exacerbates inflammation and leads to a rapid increase in arthritic disease activity
In order to assess the importance of T regs in regulating DTHA severity, we selectively depleted T regs by administration of DT at 24 and 48 h after arthritis induction to Foxp3-DTR-eGFP mice (FoxP3-DTR + ) and littermate controls (FoxP3-DTR − ).This resulted in depletion of all eGFP + cells in FoxP3-DTR + mice as confirmed by flow cytometry (Fig. S2).However, eGFP + cells eventually repopulate the T reg niche after transient DT treatment, as shown in Fig. S2C.Depletion of T regs rapidly led to increased paw and ankle swelling, which persisted up to at least day 10 (Fig. 2A,B).In T reg -depleted mice, levels of the acute-phase protein serum amyloid P component (SAP) were elevated in serum on day 10 compared to control mice (Fig. 2C), as was C-terminal telopeptide of type I collagen (CTX-I; Fig. 2D), whereas tartrateresistant acid phosphatase (TRAP) serum levels were not increased (Fig. 2D).Histopathological analysis on day 10 showed increased infiltration of macrophages into the T reg -depleted paws (Fig. 2E,F).Histological scoring revealed a tendency towards increased bone erosion in T reg -depleted mice, whereas there was no difference in cartilage degradation and synovitis (Fig. 2G).However, inflammation and new bone formation at extra-articular sites were increased in T reg -depleted mice, and so was the histopathological sum score (Fig. 2H).Interestingly, we observed a correlation (r=0.5949,P=0.0092, n=18) between extra-articular inflammation and new bone formation (Fig. 2I), indicating that new bone formation in DTHA is an inflammation-associated process.

Depletion of T regs leads to an increase in biomarkers of arthritis and inflammatory cells in the affected paw and in circulation
Next, we investigated the mechanisms of the increased arthritis response following T reg depletion by sampling on day 4 and 7 post-arthritis induction.We found increased paw swelling in T reg -depleted mice (Fig. 3A), and increased levels of SAP, IL-6 and matrix metalloproteinase 3 (MMP3) on both day 4 and 7 in T regdepleted mice (Fig. 3B).Levels of CTX-I were increased on day 7 but not day 4. Neutrophil chemoattractant granulocyte colony stimulating factor (G-CSF) was increased on day 4, but not day 7 (Fig. 3B).Interestingly, we detected production of anti-MCV in DTHA, and, following T reg depletion, anti-MCV antibody levels were increased in serum (Fig. 3B).In order to investigate whether T regs play a role in regulating bone resorption in DTHA, we measured levels of RANKL, osteoprotegerin (OPG) and TRAP in whole paw tissue.We found that levels of both RANKL and TRAP increased in the arthritic paws on day 4 and 7, and the RANKL/OPG ratio increased after T reg depletion (Fig. 3C).We also investigated the cellular players in the increased inflammatory response after T reg depletion.Systemic neutrophilia and increased numbers of neutrophils in inflamed paws and dPLNs were observed (A) DTHA was induced by immunising mice with methylated bovine serum albumin (mBSA) emulsified in complete Freund's adjuvant (CFA) on day −7, followed by an intravenous (i.v.) dose of an anti-type-II-collagen antibody cocktail (anti-CII) on day −3.On day 0 the mice were challenged with mBSA in the right hind paw and with vehicle in the left hind paw.Δpaw swelling on day x was calculated by subtracting the baseline measurement made on day 0 prior to challenge from the measurement made on day x.n=5-35; mean±s.d.shown.(B) T regs were defined as live CD45 + TCRβ + CD4 + CD25 + FoxP3 + cells, and effector T cells (T eff ) by inverting the CD25 + FoxP3 + gate within the CD45 + TCRβ + CD4 + gate.T regs were further gated on the activation marker CD44 and the proliferation marker Ki67.Total numbers of T regs , T eff , CD44 + T regs and Ki67 + T regs are shown in blood (cells/ml), the popliteal lymph node draining the arthritic paw (dPLN) and the popliteal lymph node draining the control paw (ndPLN).Mean±s.e.m. is shown, n=5.(Fig. 4A,B).Increased numbers of macrophages in the inflamed paws and dPLNs were also observed (Fig. 4B and data not shown).Numbers of CD4 + T cells in the paw were increased compared to controls on day 7 but not day 4 (Fig. 4B), but there was no increase in CD4 + T cells in the dPLN (Fig. 4C).However, CD4 + T cells in the dPLN were found to be more activated (Fig. 4D).Taken together, these data show that depletion of T regs increases systemic inflammation, increases the production of autoantibodies, increases bone resorption, increases the numbers of neutrophils and macrophages in the inflamed paw and its draining lymph node, and enhances T-cell activation.
Concomitant depletion of T regs and blockade of IL-10 signalling increased paw and ankle swelling to an even larger degree than T reg depletion or IL-10 blockade alone We found that concomitant depletion of T regs and blockade of IL-10 signalling increased paw and ankle swelling to an even larger degree than T reg depletion or IL-10 blockade alone (Fig. S3A-C).Active bone remodelling on day 10 post-arthritis induction, measured through serum CTX-I, was only detectable in mice with both depleted T regs and blocked IL-10R signalling (Fig. S3D).

Depletion of T regs increases bacterial diversity transiently without affecting the overall major taxonomic composition of fecal microbiota
Fecal samples from the distal colon were collected before DTHA induction on day 0, and on day 4 and 9 following DTHA induction with or without concurrent depletion of T regs .The composition of the microbiota was investigated by sequencing of the V4 hypervariable region of the 16S rRNA gene.Principal coordinate analysis (PCoA) at the level of operational taxonomic units (OTUs) revealed no significant changes in overall clustering following DTHA induction (Fig. S4).However, we found that the bacterial composition of the T reg -depleted mice was more diverse on day 4 compared to the non-depleted mice (Fig. 5A), with a tendency towards an increase in alpha diversity for the T reg -depleted mice from day 0 to 4 (P=0.141),whereas beta diversity did not differ significantly (Fig. 5B).Interestingly, alpha diversity correlated inversely with ankle swelling for the T reg -depleted group at day 4 (Spearman's ρ=−0.964;P<0.01; n=7).Analysis at the phylum level revealed no significant differences between the two groups, or between days within each group (Fig. 5C).We detected no significant changes in the relative abundance of genera (Fig. 5D) between groups, or between days within each group, although the abundance of Lactobacillus tended to be higher on day 4 compared to day 0 and 9 (P=0.084and 0.12 for FoxP3-DTR + and P=0.229 and 0.305 for FoxP3-DTR − , respectively).Correlations between taxa abundance and ankle swelling, paw swelling, TRAP, and CTX-I, respectively, did not show high correlation coefficients or patterns of correlations (data not shown).Taken together, these data show that DTHA induction alone or in conjunction with the depletion of T regs induce only minor and transient compositional changes in the fecal bacterial microbiota.Thus, contrasting with results obtained by analysis of RA patients (Zhang et al., 2015), no significant arthritis-associated changes in the gut microbiota were detected in the mouse models.

Depletion of T regs leads to increased production of inflammatory mediators locally in the arthritic paws
To further elucidate the mechanisms of increased inflammation following T reg depletion in DTHA, we measured levels of a range of cytokines and chemokines in whole paw tissue.The cytokines IL1β, IL-17, IL-12( p70), IL-6 and IL-10 were all increased in paw tissue on both day 4 and 7 compared with controls (Fig. 6).IFNγ was increased only on day 4 and TNFα only on day 7 (Fig. 6).The chemokines CCL2 (MCP-1), CXCL10 (IP-10), CXCL5 (LIX), CCL5 (RANTES) and CXCL9 (MIG) were all increased in paw tissue on both day 4 and 7 compared with controls (Fig. 6).CXCL2 (MIP-2), CCL3 (MIP-1α) and CXCL1 (KC) were increased only on day 7 and tended to increase on day 4 (Fig. 6).Levels of G-CSF were increased in paw tissue on day 7 and tended to increase on day 4 (Fig. 6).Granulocyte-macrophage colony-stimulating factor (GM-CSF) was increased in paw tissue both on day 4 and 7 after DTHA induction (Fig. 6).Taken together, these data show a profile of inflammatory mediators that both favours and is a result of Th17 T-cell polarisation, T-cell infiltration and activation, and neutrophil and macrophage infiltration and activation.
Treatment with anti-IL-17 mAb rescues mice from exacerbation of disease caused by depletion of T regs Based on our findings so far, we hypothesised that IL-17 could play a part in the increased disease activity observed following T reg depletion.To test this, we set up a study where T reg -depleted mice were also treated with anti-IL-17 mAb.We found that anti-IL-17 treatment reduced paw swelling in mice both with and without intact T reg compartments (Fig. 7A).We also found that anti-IL-17 treatment of T reg -depleted mice reduced paw swelling to levels seen in isotype-treated, non-T reg -depleted mice, thus rescuing them from the exacerbated disease caused by the depletion of T regs (Fig. 7A).Starting anti-IL-17 treatment at immunisation (Fig. S5) did not reduce paw swelling to a greater degree than anti-IL-17 treatment begun at the time of arthritis induction (Fig. 7A).Anti-IL-17 treatment did not have an effect on the production of anti-MCV (Fig. 7B).There was a tendency towards reduced myeloperoxidase (MPO) levels in anti-IL-17-treated mice both with and without intact T reg compartments, compared with isotype-control-treated animals (Fig. 7C).In addition, neutrophil numbers in circulation were reduced after IL-17 blockade in mice both with and without an intact T reg compartment (Fig. 7D).We found that blockade of IL-17 reduced levels of IL-6, CXCL1, CXCL5 and G-CSF in arthritic paws of mice with intact T reg compartments and tended to reduce IL-6 and CXCL1 in mice depleted of T regs (Fig. 7E).IL-17 blockade reduced the levels of IL-6 and CXCL1 in mice that received T reg depletion+anti-IL-17 to levels seen in non-T reg -depleted+isotypetreated mice (Fig. 7E).Importantly, anti-IL-17 treatment reduced levels of RANKL in paw tissue in mice both with and without intact T reg compartments (Fig. 7F).In line with this, there was a tendency towards a reduction of TRAP levels after IL-17 blockade.T reg depletion reduced the RANKL/OPG ratio and IL-17 blockade increased the RANKL/OPG ratio, indicating that T reg depletion favours a more resorptive state, whereas IL-17 blockade reverses the effect of T reg depletion on bone erosion (Fig. 7F).In summary, we find that IL-17 blockade rescues mice from the increased disease activity observed after T reg depletion, and does so through a reduction in neutrophil numbers and local IL-6, CXCL1 and RANKL levels.

DISCUSSION
In RA, dysregulation of the T reg compartment has been reported, and T reg numbers and levels of IL-17 expression correlate with disease severity (Noack and Miossec, 2014;Zaiss et al., 2010).Different approaches to tipping the T eff /T reg balance towards T reg dominance have been put forward as potential therapeutic strategies in RA and other inflammatory diseases; however, none have yet been fully successful (Miyara et al., 2014).Knowledge of the roles played by T regs in experimental arthritis models is instrumental for the development of safe and effective therapeutic strategies involving T regs in humans.In this study, we investigated the kinetics of T regs in DTHA and characterised the exacerbated inflammatory response that resulted from selective depletion of T regs .To our knowledge, this is the first study to characterise the mechanisms of exacerbated disease in an arthritis model following depletion of T regs using Foxp3-DTR-eGFP mice.Transient DT treatment, although sufficient to give the desired biological effect in many cases, can also leave a residual number of eGFP + FoxP3 + cells, and these cells eventually repopulate the T reg niche (Lahl et al., 2007;Lahl and Sparwasser, 2011).However, because we were interested in examining the effects of T reg depletion on inflammation in the early stages of DTHA, our assessment was that this mouse was suitable for our purpose, despite its limitations.
T reg numbers in DTHA expanded rapidly and these cells were found to be highly activated and proliferating in lymphoid and paw tissue.This indicates that T regs are important in the early stages of DTHA.However, the fact that we observed severe inflammation in the affected paw despite the expansion of T regs could be the result of a dysregulation of the Th17/T reg balance, because increased Th17 activity is also seen in the early stages of DTHA (Atkinson et al., 2012).In support of this, we also show in the present study that neutralisation of IL-17 reduced severity of disease and could prevent the detrimental effects of T reg depletion in the model.
Selective T reg depletion increased arthritis severity, and the observed correlation between paw inflammation and new bone formation indicated that new bone formation could be linked to inflammation in DTHA.In support of this, we have previously shown that suppression of inflammation in DTHA with anti-TNFα antibodies (Atkinson et al., 2012) or anti-C5aR antibodies (Atkinson et al., 2015) could reduce bone formation.Inflammation has also been linked to bone formation in the CIA, proteoglycan-induced arthritis (PGIA) and SKG inflammatory arthritis models (Denninger et al., 2015;Schett et al., 2009;Haynes et al., 2012;Keller et al., 2012), and in ankylosing spondylitis (Maksymowych et al., 2013;Wanders et al., 2005).However, T reg depletion also led to increased bone resorption in DTHA in the present study.The presence of these two opposing phenotypes can seem counterintuitive, but the bone formation and bone erosion take place at different areas of the paw (Atkinson et al., 2012), and simultaneous bone erosion and bone formation has also been observed in the CIA model in a recent study (Denninger et al., 2015).The net effect, however, is a negative association between inflammation and bone mass, as previously described (Gilbert et al., 2013;Proulx et al., 2007;Walsh et al., 2005).
Depletion of T regs led to increased CD4 + T-cell activation and to increased recruitment of CD4 + T cells to the site of inflammation.This shows that the control of T eff activation in the dPLN and T eff numbers in the paw is one of the mechanisms of T reg -mediated suppression of DTHA.We also observed increased levels of anti-MCV, demonstrating an augmented immune response to a citrullinated self-antigen.Anti-MCV antibodies and other anticitrullinated protein antibodies (ACPAs) are used as prognostic biomarkers of disease and joint erosion in RA, and ACPA levels correlate with serum CTX-I in RA patients (Gudmann et al., 2015;Jilani and Mackworth-Young, 2015).The fact that we observe a break of self-tolerance to MCV in DTHA and that anti-MCV levels increase with increased disease severity adds to the strengths of DTHA as a model of RA.However, the current study did not address which, if any, role anti-MCV plays in the pathogenesis of DTHA.
Previous studies have shown that T regs are able to limit the accumulation of Th17 cells at inflammatory sites and draining lymph nodes (Lohr et al., 2006).Increased IL-17 production was observed in the T reg -depleted mice in the present study.This is supported by the observed increase in production of IL-6, IL-23 and IL-1β, which are cytokines that favour Th17 polarisation (Lubberts et al., 2005;  S1.MCP, monocyte chemoattractant protein; CCL, CC chemokine ligand; MIP, macrophage inflammatory protein; CXCL, CXC chemokine ligand; IP, interferon gamma-induced protein; MIG, monokine induced by gamma interferon; RANTES, regulated on activation, normal T cell expressed and secreted; G-CSF, granulocyte colony-stimulating factor; GM-CSF, granulocyte-macrophage colony-stimulating factor.Results presented in this figure are representative of two independent experiments.S1.This experiment was carried out once.(F) Levels of RANKL, OPG and TRAP were measured in whole paw homogenate supernatants from T reg -depleted mice (blue) and non-depleted mice (black) treated with anti-IL-17 mAb or isotype control, using ELISA.Mean±s.e.m. shown, n=10.ns, not significant, *P≤0.05.Dotted line represents mean OPG level in n=2 naïve mice.CCL, CC chemokine ligand; MIP, macrophage inflammatory protein; CXCL, CXC chemokine ligand; IP, interferon gamma-induced protein; G-CSF, granulocyte colony-stimulating factor; GM-CSF, granulocyte-macrophage colony-stimulating factor; MIG, monokine induced by gamma interferon; OPG, osteoprotegerin; RANKL, receptoractivator of nuclear factor κB ligand; RANTES, regulated on activation, normal T cell expressed and secreted; TRAP, tartrate-resistant acid phosphatase.Sonderegger et al., 2008;Sutton et al., 2006).IL-17 is an important cytokine in experimental arthritis when T cells play a part in the pathogenesis (Atkinson et al., 2012;Lubberts et al., 2005).Blockade of IL-17 has previously been shown to ameliorate experimental arthritis and is also a promising therapeutic avenue in RA (Lubberts et al., 2001(Lubberts et al., , 2004;;Pöllinger et al., 2011;van den Berg and McInnes, 2013).We found that treatment with anti-IL-17 could attenuate disease in DTHA, and the results indicate that IL-17 is more important for DTHA pathogenesis after the antigen challenge.The increased IL-17 production after T reg depletion could also contribute to the increased bone erosion observed, and results from other studies support this.IL-17 induces osteoclastogenesis from human monocytes in the absence of RANKL (Yago et al., 2009) and upregulates RANKL expression (Gizinski and Fox, 2014), and blocking IL-17 prevents joint destruction in CIA through an increased RANKL/OPG ratio (Lubberts et al., 2004).Our results thus suggest that the increase in bone erosion could be driven by an increase in the RANKL/OPG ratio driven by increased IL-17 activity.
We observed increased IL-10 production in arthritic paws after T reg depletion, which suggests that the contribution of IL-10 by T regs to immunoregulation in DTHA is negligible.In a study in which T regs were depleted using anti-CD25 mAb, we found that concomitant blockade of IL-10 signalling increased paw and ankle swelling to an even larger degree than T reg depletion or IL-10 blockade alone.In the absence of T regs , regulatory B cells and their release of IL-10 could represent a compensatory regulatory pathway.B-cell-derived IL-10 has been shown to regulate inflammation in other mouse models (Fillatreau et al., 2002;Mauri et al., 2003;Mizoguchi et al., 1997).In experimental autoimmune encephalitis (EAE), components from Mycobacterium tuberculosis in the complete Freund's adjuvant (CFA) that is used to induce disease provide Toll-like receptor (TLR) agonists that can trigger regulatory functions of B cells that are required for the resolution of EAE (Lampropoulou et al., 2008).
Because B cells are not required for the induction of DTHA (Tanaka et al., 2007) and because CFA containing M. tuberculosis is also used for induction, it could be speculated that B cells indeed play an immunoregulatory role in DTHA and that depletion of B cells on top of T reg depletion would further exacerbate DTHA.
We also observed an increased number of neutrophils in paws, dPLNs and blood as well as increased production of granulopoietic factors and neutrophil chemoattractants after T reg depletion.This could in part be driven by the increased IL-17 production, because IL-17 is known to stimulate neutrophil recruitment through induction of G-CSF, GM-CSF and chemoattractants, including CXCL1 (Hamilton, 2008;Kolls and Lindén, 2004;Park et al., 2005;Schwarzenberger et al., 2000).Increased levels of CXCL1 and CXCL2 at the site of inflammation have also previously been shown to result from T reg depletion in mice (Richards et al., 2010).Indeed, blockade of IL-17 reduced neutrophil numbers in T reg -depleted mice and reduced levels of CXCL1.IL-6 production can be induced by IL-17 (Kolls and Lindén, 2004) and, in line with this, we observed a major decrease of IL-6 production following IL-17 blockade.IL-6 plays an important role in neutrophil trafficking (Fielding et al., 2008;Lally et al., 2005), so the observed decrease in blood neutrophil numbers after IL-17 blockade could be a result of the decreased IL-6 levels.Taken together, these data support a regulatory role for IL-17 in neutrophil recruitment through IL-6 and CXCL1.
Recent findings have suggested an association between the gut microbiota and RA (Scher et al., 2013;Zhang et al., 2015), and interactions between the gut microbiota and immune function, including T regs , neutrophils and IL-17, are well established (Maynard et al., 2012).However, in our model, depletion of T regs did not significantly change the bacterial composition, indicating that immune-driven changes in the gut microbiota are unlikely to play an important role in disease progression in our model.We observed a tendency towards higher alpha diversity in T reg -depleted mice and that alpha diversity correlated inversely with ankle size at day 4. Interestingly, we observed a tendency towards an increased abundance of Lactobacillus on day 4, which also correlated with maximal ankle swelling.Lactobacillus has been shown to be overrepresented in RA patients compared to healthy subjects (Liu et al., 2013;Zhang et al., 2015), and the abundance of Lactobacillus was also shown to correlate with disease activity in RA patients.Furthermore, mono-colonisation with Lactobacillus bifidus was shown to be sufficient for inducing joint inflammation in a mouse model of spontaneous arthritis via TLR2-and TLR4-dependent signalling affecting the balance between T reg , Th17 and Th1 in the gut (Abdollahi-Roodsaz et al., 2008).When comparing the microbial changes observed in laboratory animals to humans, it is important to consider the marked environmental differences.Moreover, in DTHA, severe disease is observed after just a few days, whereas RA develops over several years.Although the composition of the microbiota can change within days, the short duration of the model might not be sufficient to markedly affect the composition of the gut microbiota in a manner that reflects the situation in human disease pathogenesis.Nevertheless, the findings suggest a link between Lactobacillus and joint inflammation both in rodent models and RA patients.

Conclusion
Taken together, our findings show that T regs are important for the containment of inflammation in DTHA, and that T reg depletion in DTHA acts, at least partly, through increasing IL-17 and RANKL activity and through increasing numbers of neutrophils in blood, dPLNs and paws.We show that IL-17 drives an increase in inflammation in the absence of T regs : administration of anti-IL-17 mAb just prior to T reg depletion rescues mice from the exacerbation of disease associated with T reg depletion through a reduction in IL-6 and CXCL1 production, reduced neutrophil numbers and a reduction in the RANKL/OPG ratio.Thus, an imbalance in the T reg and Th17 relationship could be a driver of inflammation in DTHA.To our knowledge, this is the first study using the Foxp3-DTR-eGFP mouse on a C57BL/6 (B6) background for T reg depletion in an arthritis model, and we here demonstrate the usefulness of the approach and build upon the already existing body of work in this new model.Many genetically modified mice are bred on a B6 background, and can thus be combined with the Foxp3-DTR-eGFP C57BL/6 mouse, making possible the study of T reg depletion in concert with the absence/overexpression of proteins relevant for disease pathways.However, the limitations of the Foxp3-DTR-eGFP mouse generated by BAC transgenesis must be kept in mind when choosing this model system.Given the relevance of T regs in RA and the possibility of developing T reg -mediated therapies for RA, using this approach for selective depletion of T regs could yield important results to help understand the diverse roles of T regs in RA and other inflammatory diseases.

Mice
Female C57BL/6J mice were purchased from Taconic, Ry, Denmark.Female Foxp3-DTR-eGFP mice were generated as described in Lahl et al. (2007), and bred and maintained at Forschungseinrichtungen für Experimentelle Medizin, Charité Universitätsmedizin Berlin or purchased from Jackson Laboratories (Bar Harbor, ME) on a NOD background (generated as described in Feuerer et al., 2009) and backcrossed for at least eight generations onto the C57BL/6 background at Taconic, Ry, Denmark.Foxp3-DTR-eGFP mice express a fusion protein consisting of the diphtheria toxin (DT) receptor (DTR) and enhanced green fluorescent protein (eGFP) under the control of the Foxp3 gene locus inserted by BAC transgenesis.This construct allows selective depletion of T regs by DT injection.Animals were housed in a facility with a 12-h light/dark cycle and with free access to water and standard rodent chow (Altromin ® ).All animal experiments performed in Germany were in accordance with the guideline 2010/63/EU from the European Union and the European Convention for the protection of vertebrate animals used for experimental and other scientific purposes.Animal protocols were approved by the ethics committee and the Berlin state authorities (LAGeSo registration # G0331/08).All animal experiments performed in Denmark were conducted according to Danish legislation and have been approved by the Danish Animal Inspectorate and the Novo Nordisk ethical review board.

Depletion of T regs and IL-10R blockade
For depletion of T regs , mice were dosed with 200 µg rat anti-mouse CD25 mAb (clone PC.61, BioXcell, West Lebanon, NH) or rat IgG1 antitrinitrophenol (TNP) isotype control (Novo Nordisk A/S, Måløv, Denmark) 11 and 8 days prior to immunisation.For blockade of IL-10R, mice were dosed with 250 µg rat anti-mouse IL-10R mAb (clone 1B1.3A,BioXcell, West Lebanon, NH) or rat IgG1 anti-TNP isotype control (Novo Nordisk A/S, Måløv, Denmark) three times weekly from the time of immunisation.Depletion of T regs was confirmed by flow cytometry on whole blood samples drawn 48 h after the second dose of anti-CD25.Mouse regulatory T cell staining kit (eBioscience, San Diego, CA) was used.Antibodies used were CD45-PerCP, CD4-APC and CD25-FITC all from BD Biosciences (NJ), and FoxP3-PE (eBioscience, San Diego, CA).

Anti-IL-17 mAb treatment
Mice were treated with 400 µg rat anti-mouse IL-17 mAb (clone 17F3, BioXcell, West Lebanon, NH) or rat IgG1 isotype control (clone MOPC-21, BioXcell, West Lebanon, NH) in 200 µl PBS on the day of arthritis induction and with 200 µg rat anti-mouse IL-17 mAb or isotype control in PBS after 48 h and every 48 h until the experiment was terminated on day 7.For results presented in Fig. S4, mice were dosed with 200 µg rat antimouse IL-17 mAb (clone 17F3, BioXcell, West Lebanon, NH) or rat IgG1 isotype control (clone MOPC-21, BioXcell, West Lebanon, NH) in 200 µl PBS from the time of immunisation and every 48 h until study termination.

Histopathology
Paws were processed and stained with hematoxylin and eosin (H&E), Safranin O and for tartrate-resistant acid phosphatase (TRAP) as previously described (Atkinson et al., 2012).TRAP stains osteoclasts red, and Safranin O stains cartilage red.The intensity of Safranin O staining is directly proportional to the proteoglycan content in cartilage.Pathological changes in the paws were assessed on H&E-, TRAP-and Safranin-O-stained sections.The extra-articular infiltration of inflammatory cells (assessed on a scale of 0-3) and arthritic changes were assessed separately.Arthritic changes were assessed on metatarsal and tarsal joints as synovitis, cartilage destruction and bone erosion, and scored separately on a 0-3 scale.For each of the three parameters of arthritic changes, an average between the two joint areas was calculated.In addition, new bone formation and extra-articular infiltration overall in the paw was scored on a 0-3 scale.New bone formation was defined as the presence of osteophytes, as described in Atkinson et al. (2012).The histology sum score was calculated by adding the five scores (extra-articular infiltration, synovitis, cartilage destruction, bone erosion and bone formation), whereas the extra-articular infiltration score is left out in the arthritis score.The person who performed the evaluation was blinded to the experimental setup.

Immunohistochemistry and digitalised image analysis
Immunohistochemical (IHC) detection of macrophages and neutrophils in the paws was performed as previously described (Atkinson et al., 2015).The sections were all digitally scanned and studied using a NanoZoomer Digital Pathology Virtual Slide Viewer (Hamamatsu Photonic, Shizuoka, Japan).Automated digital image analyses of the infiltrating macrophages (F4/80 + ) and neutrophils (Ly6B.2+ ) in the paws were performed using the Visiopharm Integrator System (VIS; version 4.2.2.0, Visiopharm, Hoersholm, Denmark).On individual digital images of the arthritic paw a region-of-interest (ROI) was automatically defined of the entire paw.The bone marrow was outlined manually and excluded from analysis, as was hair follicles and artefacts.Next, an analysis was run inside the ROI to detect the brown DAB (3-3′-diamino-benzidine-tetrahydrochloride) staining of the target cells.The results are given as percentage of tissue area positive for F4/ 80 or Ly6B.2 of the entire paw area.

Multiplex analysis of inflammatory markers in paw homogenate
Hind paws were sampled at selected times after DTHA induction and each paw was placed in 1.25 ml of an ice-cold custom-made homogenization buffer containing a solution of 200 mM NaCl, 5 mM EDTA, 10 mM Tris, 10% glycerin, 1 mM phenylmethylsulfonyl fluoride (PMSF), 1 mg/ml leupeptin and 28 mg/ml aprotinin with a pH value of 7.4 (Ampliqon, Skovlunde, Denmark).The paws were homogenized by using a T25 Ultraturrax homogeniser (IKA, Staufen, Germany) followed by centrifugation at 10,000 g for 15 min.The supernatants were decanted and centrifuged once more at 10,000 g for 15 min.The final supernatants were analysed undiluted for levels of IL-1β, TNFα, IFNγ, IL-6, IL-10, IL-17, IL-12( p40), MIP-2 (CXCL2), IP-10 (CXCL10), LIX (CXCL5), RANTES (CCL5), MIG (CXCL9), MIP-1α (CCL3), KC (CXCL1), MCP-1 (CCL2), G-CSF and GM-CSF using bead-based Luminex ® xMAP ® technology with Milliplex kits from Millipore (Billerica, MA) according to the manufacturer's instructions.For statistical analysis, any values below the detection limit were set to the detection limit for the analyte in question and any values above the detection limit were set to the upper detection limit for the analyte in question.

Enzyme-linked immunosorbent assays (ELISAs)
Levels of serum amyloid P component (SAP) were measured in serum using kits from Genway (San Diego, CA) according to the manufacturer's instructions.Levels of receptor-activator of nuclear factor kappa B (RANKL), osteoprotegerin (OPG) and IL-23 were measured in paw homogenate supernatants using kits from R&D Systems (Minneapolis, MN) according to the manufacturer's instructions.Levels of CTX-I and TRAP were measured in serum and paw homogenate supernatants using ELISA kits from Immunodiagnostic Systems (Boldon, UK) according to the manufacturer's instructions.Levels of MMP3, IL-6 and G-CSF were measured in serum using kits from R&D Systems (Minneapolis, MN) according to the manufacturer's instructions.Antibodies to mutated citrullinated vimentin (anti-MCV) were measured using kits from Orgentec Diagnostika GmbH (Mainz, Germany) with the modification that the detection antibody supplied in the kit was substituted with goat antimouse IgG-HRP (Invitrogen, Carlsbad, CA) diluted 1:3000 and the substrate solution supplied was substituted with substrate reagents from R&D Systems (Minneapolis, MN).

mRNA deep sequencing
Arthritic hind paws were removed at the hairline and homogenised in Qiagen RLT buffer (Qiagen, Germantown, MD) with 1% 2betamercaptoethanol (Sigma-Aldrich, St Louis, MO) and stored at −80°C.The Ambion Magmax prep protocol was used to extract total RNA.RNA quantity was assessed using the Nanodrop and RNA quality was assessed using the Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA).RNAseq libraries were prepared using the Illumina TruSeq Sample Prep Kit (Illumina, Inc., San Diego, CA).120 ng of total RNA was used for input.Barcode adapters were added to samples in such a way as to allow pooling of samples in flow cell lanes in a randomized pattern relative to sample annotation and source.Samples were sequenced using an Illumina HiSeq 2000 (Illumina, Inc., San Diego, CA) at a multiplexing level sufficient to generation 10-million to 25-million reads per sample.Following generation of sequence, reads were aligned to the mouse genome version NCBI m37 using TopHat (http://tophat.cbcb.umd.edu/).Quality control of sequencing data was accomplished by using the ShortRead package in R (http://www.r-project.org/).Aligned reads were mapped to Ensembl transcript models and converted to reads per kilobase of transcript per million mapped reads (RPKM) values using Cufflinks (http:// cufflinks.cbcb.umd.edu/).

Bacterial 16s rDNA sequencing
Bacterial DNA from fecal samples from the distal colon was extracted using the NucleoSpin soil kit (Macherey-Nagel, Düren, Germany) following the manufacturer's protocol.Yield and integrity of the DNA were assessed by Nanodrop and agarose gel electrophoresis, respectively.16S rDNA amplification and library generation were performed as previously described (Holm et al., 2015).PCR products were purified using Agencourt AMPure XP beads (Beckman Coulter, Brea, CA) and normalized to 16-21 ng/μl per sample.Subsequently, samples were pooled (5 μl of each sample) and quantified using a Qubit dsDNA HS Assay Kit and a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA).14 pM library and 0.7 pM PhiX Control v3 (Illumina, San Diego, CA) were sequenced using an Illumina MiSeq V2 PE500 cartridge (500 cycles) on an Illumina MiSeq.Sequences were analysed using QIIME v1.7.0 with default settings including quality-based sequence trimming, removal of primers, assembly of paired-end sequences and chimera checking (Caporaso et al., 2010).De novo operational taxonomic unit (OTU)-picking was performed by UCLUST (Edgar, 2010) utilizing 97% sequence similarity.OTUs were assigned against the Greengenes database v11_2 (DeSantis et al., 2006) using the RDP-classifier (Wang et al., 2007) with an 80% confidence threshold.Analyses were performed in R v3.2.0 using the metagenomeSeq (Paulson et al., 2013), PhyloSeq (McMurdie andHolmes, 2013) and Vegan (Oksanen et al., 2011) packages.Data was filtered for low-abundant OTUs by removing OTUs present in fewer than three of the 42 samples and with a relative abundance across all samples ≤0.5%.Analyses on filtered data were performed with an average of 20,731±7655 (s.d.) sequences per sample compared to an initial 24,974±9382 (s.d.) sequences before filtering.Alpha diversity (Shannon index) was calculated on unfiltered data.Read counts were normalized using metagenomeSeq (Paulson et al., 2013), which utilizes a cumulative-sum scaling where raw counts are divided by the cumulative sum of counts up to a particular quantile.Principal coordinate analysis (PCoA) using Bray-Curtis dissimilarity indices, analysis of differential abundance of taxa and beta diversity analysis using Whittaker's species turnover were performed on filtered and normalized data.

Statistics
Statistical analyses were conducted using GraphPad Prism software version 6.04.Non-parametric data or non-normal parametric data were analysed using the Mann-Whitney U-test, and parametric data were analysed using a two-sided unpaired Student's t-test or one-way ANOVA.For statistical analysis of the histology score data, a two-sided unpaired Student's t-test with Welch's correction was used.For bacterial 16S rDNA, statistical analyses were conducted using R differences in abundance of phyla, alpha diversity measures and beta diversity between and within the FoxP3-DTR − and FoxP3-DTR + groups over the course of the study using the Wilcoxon Rank-Sum and Wilcoxon Signed-Rank test, respectively, with Benjamini-Hochberg P-value adjustment.Adonis test of significance was performed using the calculated Bray-Curtis distance matrix, to test for overall differences between the FoxP3-DTR − and FoxP3-DTR + groups and within each group over the course of the study.Correlation analyses between taxonomic abundances, alpha diversity, serum markers, and paw and ankle size were performed in R using Spearman's rank correlation test with Benjamini-Hochberg P-value adjustment.Differences between groups and correlation coefficients were considered significant when P≤0.05 and levels of significance were assigned as *P≤0.05,**P≤0.01 and ***P≤0.001.In Figs 6 and 7E P-values were corrected for mass significance using the False Discovery Rate method.The raw and corrected P-values are shown in Table S1.

Fig. 1 .
Fig. 1.Dynamics and phenotypic characterisation of regulatory T cells (T regs ) during delayed-type hypersensitivity arthritis (DTHA).(A) DTHA was induced by immunising mice with methylated bovine serum albumin (mBSA) emulsified in complete Freund's adjuvant (CFA) on day −7, followed by an intravenous (i.v.) dose of an anti-type-II-collagen antibody cocktail (anti-CII) on day −3.On day 0 the mice were challenged with mBSA in the right hind paw and with vehicle in the left hind paw.Δpaw swelling on day x was calculated by subtracting the baseline measurement made on day 0 prior to challenge from the measurement made on day x.n=5-35; mean±s.d.shown.(B) T regs were defined as live CD45 + TCRβ + CD4 + CD25 + FoxP3 + cells, and effector T cells (T eff ) by inverting the CD25 + FoxP3 + gate within the CD45 + TCRβ + CD4 + gate.T regs were further gated on the activation marker CD44 and the proliferation marker Ki67.Total numbers of T regs , T eff , CD44 + T regs and Ki67 + T regs are shown in blood (cells/ml), the popliteal lymph node draining the arthritic paw (dPLN) and the popliteal lymph node draining the control paw (ndPLN).Mean±s.e.m. is shown, n=5.(C) Foxp3 and Ctla4 mRNA in arthritic paws measured by mRNA deep sequencing.RPKM, reads per kilobase of transcript per million mapped reads.Mean shown, n=4.**P≤0.01,***P≤0.001,difference to levels on day 0, one-way ANOVA.(D) Total cell counts in the inguinal lymph nodes.T regs were defined as live CD45 + TCRβ + CD4 + CD25 + FoxP3 + cells and T eff by inverting the CD25 + FoxP3 + gate within the CD45 + TCRβ + CD4 + gate, n=6, mean±s.e.m. shown.Student's t-test, ***P≤0.001.(E) Representative flow cytometry plots from day 2 post-DTHA induction showing CD25 hi FoxP3 hi T reg populations (box) in dPLN, ndPLN and blood.Cells were defined as live CD45 + TCRβ + CD4 + before gating on CD25 and FoxP3.The top panel shows FoxP3 staining and the bottom panel staining with IgG2a isotype control antibody.(F) Percentage of CD25 hi FoxP3 hi T regs in dPLN, ndPLN and blood after induction of DTHA and in naïve mice.Mean±s.e.m. shown, n=3-6.Difference to percentages in naïve mice analysed using one-way ANOVA, **P≤0.01,***P≤0.001.(G) Percentage of CD25 hi FoxP3 hi T regs in dPLN expressing CD44 (left) and Ki67 (right); mean±s.d.shown, n=3-6.All experiments presented in this figure were carried out once.Full flow cytometry gating strategy is shown in Fig. S1A.
Fig. 1.Dynamics and phenotypic characterisation of regulatory T cells (T regs ) during delayed-type hypersensitivity arthritis (DTHA).(A) DTHA was induced by immunising mice with methylated bovine serum albumin (mBSA) emulsified in complete Freund's adjuvant (CFA) on day −7, followed by an intravenous (i.v.) dose of an anti-type-II-collagen antibody cocktail (anti-CII) on day −3.On day 0 the mice were challenged with mBSA in the right hind paw and with vehicle in the left hind paw.Δpaw swelling on day x was calculated by subtracting the baseline measurement made on day 0 prior to challenge from the measurement made on day x.n=5-35; mean±s.d.shown.(B) T regs were defined as live CD45 + TCRβ + CD4 + CD25 + FoxP3 + cells, and effector T cells (T eff ) by inverting the CD25 + FoxP3 + gate within the CD45 + TCRβ + CD4 + gate.T regs were further gated on the activation marker CD44 and the proliferation marker Ki67.Total numbers of T regs , T eff , CD44 + T regs and Ki67 + T regs are shown in blood (cells/ml), the popliteal lymph node draining the arthritic paw (dPLN) and the popliteal lymph node draining the control paw (ndPLN).Mean±s.e.m. is shown, n=5.(C) Foxp3 and Ctla4 mRNA in arthritic paws measured by mRNA deep sequencing.RPKM, reads per kilobase of transcript per million mapped reads.Mean shown, n=4.**P≤0.01,***P≤0.001,difference to levels on day 0, one-way ANOVA.(D) Total cell counts in the inguinal lymph nodes.T regs were defined as live CD45 + TCRβ + CD4 + CD25 + FoxP3 + cells and T eff by inverting the CD25 + FoxP3 + gate within the CD45 + TCRβ + CD4 + gate, n=6, mean±s.e.m. shown.Student's t-test, ***P≤0.001.(E) Representative flow cytometry plots from day 2 post-DTHA induction showing CD25 hi FoxP3 hi T reg populations (box) in dPLN, ndPLN and blood.Cells were defined as live CD45 + TCRβ + CD4 + before gating on CD25 and FoxP3.The top panel shows FoxP3 staining and the bottom panel staining with IgG2a isotype control antibody.(F) Percentage of CD25 hi FoxP3 hi T regs in dPLN, ndPLN and blood after induction of DTHA and in naïve mice.Mean±s.e.m. shown, n=3-6.Difference to percentages in naïve mice analysed using one-way ANOVA, **P≤0.01,***P≤0.001.(G) Percentage of CD25 hi FoxP3 hi T regs in dPLN expressing CD44 (left) and Ki67 (right); mean±s.d.shown, n=3-6.All experiments presented in this figure were carried out once.Full flow cytometry gating strategy is shown in Fig. S1A.

Fig. 5 .
Fig. 5. T reg depletion in DTHA has limited effect on the overall composition of the microbiota in the distal colon.Fecal microbiota was analysed for bacterial 16S rRNA composition.(A) Alpha diversity based on Shannon index of unfiltered data for fecal samples.Median alpha diversity with upper and lower quantiles shown, n=7.*P≤0.05.(B) Beta diversity based on Whittaker's species turnover within each group.Median beta diversity with upper and lower quantiles indicated, n=7.(C) Taxa summary plot at phylum level showing changes in the microbiota composition of the fecal samples.Mean relative abundance shown, n=7.(D) Taxa summary plot at phylum level showing changes in the microbiota composition of the fecal samples.Mean relative abundance shown, n=7.Except for beta diversity, Wilcoxon rank-sum for tests between groups and Wilcoxon signed-rank for tests within groups.This experiment was carried out once.

Fig. 6 .
Fig. 6.Depletion of T regs augments production of inflammatory mediators locally in the affected paws.Whole paw homogenate supernatants were analysed for protein levels ( pg/g tissue) of a range of inflammatory markers using multiplex analysis.Mean±s.e.m. shown, n=10.Student's t-test, *P≤0.05;**P≤0.01;***P≤0.001.P-values were corrected for mass significance using the False Discovery Rate method.The raw and corrected P-values are shown in TableS1.MCP, monocyte chemoattractant protein; CCL, CC chemokine ligand; MIP, macrophage inflammatory protein; CXCL, CXC chemokine ligand; IP, interferon gamma-induced protein; MIG, monokine induced by gamma interferon; RANTES, regulated on activation, normal T cell expressed and secreted; G-CSF, granulocyte colony-stimulating factor; GM-CSF, granulocyte-macrophage colony-stimulating factor.Results presented in this figure are representative of two independent experiments.