Systemic sclerosis (SSc) is an autoimmune disease characterized by vascular damage and fibrosis of the skin and internal organs. Because activated and oligoclonally expanded CD8+ T cells can be detected in peripheral blood and lungs of SSc patients, effector memory CD8+ T cells may play a critical role for organ involvement in SSc; however, the pathogenic functions of effector memory CD8+ T cells remain incompletely understood. In this study, we performed DNA microarray analysis of the sort-purified effector memory CD8+ T cells from SSc patients and healthy controls, and showed that the expression of genes related to immune response and cell adhesion, including CD226 (also known as DNAX accessory molecule-1 [DNAM-1]), was significantly altered. Moreover, detailed analysis of CD226 revealed that CD226highCD8+ T cells were increased in SSc patients (mean, 50.7%) compared with healthy controls (32.9%) and were appreciably associated with the severity of skin sclerosis and interstitial lung disease. Furthermore, CD226+CD8+ T cells produced higher amount of various cytokines than CD226 ones, and CD226highCD8+ T cells from SSc patients showed upregulated IL-13 production and positive correlation with the cytotoxic capacity of CD8+ T cells against HUVECs. Finally, the neutralization of CD226 in CD8+ T cells impaired costimulation, cytokine productions, and cytolysis against HUVECs. These findings indicate that upregulated CD226 expression on CD8+ T cells reflects disease severity and is involved in SSc pathogenesis via the production of various cytokines, including profibrotic IL-13 and endothelial cell injury, and that CD226 may be a useful target in the treatment of SSc.

Systemic sclerosis (SSc) is an autoimmune disease characterized by vascular damage, inflammation, immune dysfunction, and fibrosis of the skin and internal organs (13). Endothelial vascular damage is an early event in SSc (1, 4) and leads to activation of the immune response with release of various cytokines, autoantibody production, fibroblast activation, and enhanced collagen synthesis (5). In this process, various types of immune cells, such as T cells, B cells, and macrophages, play important roles.

Previous studies have reported that CD8+ T cells infiltrated the skin (6) and lung (7, 8) in SSc patients and that CD8+ T cells in lung and peripheral blood showed an activated phenotype (7, 9) and Ag-driven oligoclonal expansion (8, 10, 11). These findings indicate that effector or memory CD8+ T cells, or both, are involved in SSc pathogenesis. Recent reports have shown that effector CD8+ T cells were directly involved in modulating dermal fibrosis in SSc via dysregulated IL-13 production (6, 12). Memory T cells contain two populations of central memory and effector memory cells characterized by distinct homing potentials and effector functions (13). Effector memory CD8+ T cells are memory cells that have lost the lymph node homing receptor CCR7 and CD62L expression and display characteristic sets of chemokine receptors and adhesion molecules required for homing to inflamed tissues (14). This suggests that effector memory CD8+ T cells are critical for organ involvement in SSc; however, the pathogenic functions of effector memory CD8+ T cells remain incompletely understood.

Comprehensive gene expression analysis by cDNA microarray is one of the most useful tools for understanding disease pathogenesis. PBMCs constitute the most commonly used samples for microarray analysis, because they are easily obtained and are potentially representative of multiorgan conditions (15). However, microarray analyses of PBMCs have the possibility to miss slight but important changes in minor populations because of their heterogeneity (16). Therefore, it is important to prepare well-purified subpopulations of PBMCs prior to performing microarray analysis.

In this study, we first purified effector memory CD8+ T cells by cell sorting and subsequently performed cDNA microarray analysis of the sort-purified effector memory CD8+ T cells from SSc patients and healthy controls (HCs). Based on the results of the microarray analysis, we investigated the characteristic functions of CD226+CD8+ T cells in SSc pathogenesis.

We studied 50 Japanese patients (9 for microarray analysis and 41 for FACS analysis) who were treated for SSc at the Kyushu University hospital and 31 HCs (5 and 26, respectively). We included SSc patients who fulfilled the 1980 classification criteria of the American College of Rheumatology for SSc (17) and had no other autoimmune diseases. This study was approved by the ethics committee of our institution, and the principles of the Helsinki Declaration were followed throughout the study. Informed consent was obtained from all participants.

We obtained the information from the medical records of the patients, including demographic data, clinical manifestations, laboratory findings, and medications. SSc patients were classified as having limited cutaneous or diffuse cutaneous disease according to the criteria of LeRoy et al. (18). The disease duration was calculated from the time of onset of the first non-Raynaud phenomenon. Skin sclerosis was evaluated with the modified Rodnan skin thickness score (MRSS) (19) and graded as mild (MRSS, 1–14) or moderate-to-severe (MRSS, 15–39) (20). SSc-related interstitial lung disease (ILD) was classified as extensive or limited disease based on combined evaluation with chest high-resolution computed tomography (HRCT) and pulmonary functional tests (21). In brief, extensive ILD was defined as HRCT extent >30%, or HRCT extent ranging 10–30% and forced vital capacity (FVC) < 70%; limited ILD as HRCT extent <10% or HRCT extent ranging 10–30% and FVC ≥ 70%. FVC was expressed as a percentage of predicted values according to the prediction equations of the Japanese Respiratory Society (22). Serum anti-topoisomerase I Ab and anticentromere Ab levels were measured using fluorescence enzyme immunoassay (Phadia AB, Uppsala, Sweden).

PBMCs were separated from the heparinized fresh blood using centrifugation on Lymphocyte Separation Medium gradients (MP Biomedicals, Santa Ana, CA). PBMCs obtained were stained with mAb mixtures in PBS containing 2% heat-inactivated FBS (Thermo Fisher Scientific, Waltham, MA) at 4°C for 30 min in the dark and sorted using the BD FACSAria flow cytometer (BD Biosciences, San Jose, CA). The >97% purity of each population was confirmed.

For the microarray analysis, PE-conjugated anti-CD62L (Dreg56); allophycocyanin-conjugated anti-CD3 (UCTH1); allophycocyanin-Cy7-conjugated anti-CD8 (SK1); PE-Cy5-conjugated anti-CD4 (RPA-T4), anti-CD19 (HIB19), anti-CD20 (2H7), and anti-CD56 (B159; all from BD Biosciences); PC5-conjugated anti-CD14 (RMO52; Beckman Coulter, Brea, CA); FITC-conjugated anti-CD45RO (UCHL1; BioLegend, San Diego, CA); and PE-Cy7-conjugated anti-CD45RA (HI100; eBioscience, San Diego, CA) were used. Effector memory CD8+ T cells (CD3+CD8+CD45RACD45RO+CD62L cells) were suspended in TRIzol reagent (Invitrogen, Carlsbad, CA) for RNA isolation.

For the analysis of the cytokine profiles of CD8+ T cells with or without CD226, FITC-conjugated anti-CD3 (UCTH1), allophycocyanin-conjugated anti-CD8 (RPA-T8; both from BD Biosciences), and PE-conjugated anti-CD226 (DX11; Miltenyi Biotec, Bergisch Gladbach, Germany) were used. The sorted CD226+CD8+, CD226highCD8+, CD226lowCD8+, or CD226CD8+ T cells were suspended in serum-free AIM-V medium (Invitrogen).

Total RNA was extracted from the sort-purified effector memory CD8+ T cells (5000 cells) using TRIzol, and biotinylated cRNA was synthesized with two rounds of amplification steps using the MessageAmpII aRNA Amplification Kit and Illumina TotalPrep RNA Amplification Kit (both from Applied Biosystems, Foster City, CA). Next, 750 ng of cRNA from the sample was hybridized to the HumanHT-12 v4 Expression BeadChip (Illumina, San Diego, CA). After staining and washing, the BeadChip was scanned with an Illumina Bead Array reader, and the data were compiled using Bead Studio software (Illumina).

The raw signal intensities of all samples were log2-transformed and normalized using a quantile algorithm with a “preprocessCore” library package (23) of Bioconductor software (24). The genes flagged as present (detection p < 0.01) in all samples were selected and used for further analysis. The linear models for microarray analysis (limma) package (25) of the Bioconductor software was used to identify the differentially expressed genes. Functional annotation clustering analysis was performed with Database for Annotation, Visualization, and Integrated Discovery (DAVID) using default parameters on its Web site (http://david.abcc.ncifcrf.gov/home.jsp) (26, 27). The heat map was generated using MeV software (28). Hierarchical clustering was performed using the Pearson correlation metric and average linkage clustering.

Microarray data have been deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GSE63903; http://www.ncbi.nlm.nih.gov/geo).

Freshly isolated PBMCs were stained with the following direct conjugated mAbs: FITC-conjugated anti-CD3 (UCTH1), anti-CD45RO (UCHL1), or anti-CD62L (DREG-56; Beckman Coulter); PE-conjugated anti-CD226 (DX11); PerCP-conjugated anti-CD8 (RPA-T8) or anti-CD19 (HIB19; both from BioLegend); PE-Cy5–conjugated anti-CD45RO (UCHL1); PE-Cy7–conjugated anti-CD3 (SK7); allophycocyanin-conjugated anti-CD4 (RPA-T4), anti-CD62L (DREG-56), anti-CD56 (B159; all five from BD Biosciences), or anti-TCR γ/δ (B1; BioLegend); Alexa Fluor 647-conjugated anti-chemoattractant receptor-homologous molecule on Th2 cells (CRTH2; clone BM16; BioLegend); and allophycocyanin-Cy7–conjugated anti-CD8 (SK1). Isotype control Abs (BD Biosciences) were used to determine the level of background staining. Samples were acquired with the BD FACSCalibur (BD Biosciences) or FACSAria flow cytometer. Data analysis was performed using FlowJo software (Tree Star, Ashland, OR).

The sorted CD226+CD8+ or CD226CD8+ T cells (2 × 105 cells) were stimulated with a combination of 50 ng/ml PMA (Sigma, St Louis, MO) and 1 μg/ml ionomycin (Sigma) for 6 h in serum-free AIM V medium. After centrifugation, the supernatants were collected and stored at −20°C until the analysis. The protein levels of IFN-γ, TNF-α, IL-4, IL-5, IL-10, IL-13, and granzyme B were measured using cytometric beads array (CBA; BD Biosciences) on the BD FACSCalibur flow cytometer according to the manufacturer’s instructions.

Heparinized whole blood (200 μl), diluted 1:1 with sterile RPMI 1640 (Wako, Osaka, Japan) medium, was stimulated with PMA (50 ng/ml) and ionomycin (1 μg/ml) for 5 h in the presence of GolgiStop (BD Biosciences). The activated cultures were treated with 4 ml RBC Lysis Solution (Miltenyi Biotec) for 10 min. After centrifugation, the cells were suspended in PBS containing 2% FBS and stained with FITC-conjugated anti-CD3 (UCTH1), PerCP-conjugated anti-CD8 (RPA-T8), and PE-conjugated anti-CD226 (DX11). The cells were subsequently fixed and permeabilized using BD fixation–permeabilization solution (BD Biosciences) and stained with allophycocyanin-conjugated anti–IL-13 (JES10-5A2; BD Biosciences). The samples were analyzed using the BD FACSCalibur flow cytometer.

The sorted CD226highCD8+, CD226lowCD8+, or CD226CD8+ T cells (4 × 104 cells) were stimulated with PMA (50 ng/ml) and ionomycin (1 μg/ml) for 4 h in serum-free AIM V medium. The activated cells were harvested using TRIzol reagent, total RNA was extracted, and cDNA was synthesized with QuantiTect Reverse Transcription Kit (Qiagen, Hilden, Germany). Quantitative PCR was performed using the StepOnePlus Real-Time PCR System (Applied Biosystems) with TaqMan Gene Expression Master Mix (Applied Biosystems). The following primers were purchased from Applied Biosystems: GAPDH (Hs02758991_g1), IL-13 (Hs00174379_m1), and GATA binding protein 3 (GATA3; Hs00231122_m1). Relative gene expression quantification was calculated according to the comparative cycle threshold method using GAPDH as an endogenous control.

CD8+ T cells were negatively isolated from PBMCs using magnetic beads (CD8+ T Cell Isolation kit; Miltenyi Biotec). The purity of the CD8+ T cells was >95%, as determined by flow cytometry. Freshly isolated cells were cultured in complete RPMI 1640 medium supplemented with 10% FBS, 50 IU/ml penicillin, and 20 μg/ml streptomycin (both from Life Technologies, Grand Island, NY).

HUVECs were purchased from Lonza (Basel, Switzerland) and cultured in fully supplemented endothelial growth medium (EGM-2; Lonza) at 37°C in a 5% CO2 humidified atmosphere. HUVECs were used up in passage 4–7.

HUVECs (1 × 104 cells), which expressed high levels of CD155 (one of the CD226 ligands), were plated in 96-well plates and grown overnight to confluent monolayers in EGM-2 medium. After removing the medium, the cells were washed and subsequently cocultured with purified CD8+ T cells (1 × 105 cells) in RPMI 1640 medium supplemented with 10% FBS in the presence of Anti-Biotin MACSiBead Particles loaded with CD3-Biotin (anti-CD3–coated beads; both from Miltenyi Biotec), and anti-CD226 blocking Ab (DX11; 10 μg/ml) or isotype control Ab (both from BD Biosciences). In this assay, CD8+ T cells were activated with anti-CD3–coated beads and costimulated through CD226/CD155 interactions. After 72 h of incubation, the supernatants were collected, and the protein levels of IL-13 and IFN-γ were measured using the CBA system.

The cytotoxicity of CD8+ T cells against HUVECs was evaluated using flow cytometry (29). In brief, HUVECs were labeled with 0.25 μM CFSE (Sigma) and subsequently cocultured with CD8+ T cells at E:T ratios of 10:1 (CD8+ T cells, 1 × 105 cells; HUVECs, 1 × 104 cells) and 40:1 (CD8+ T cells, 2 × 105 cells; HUVECs, 5 × 103 cells). Six hours later, the cells were collected and labeled with 7-aminoactinomycin D (7-AAD; BD Biosciences), and then the percentage of dead HUVECs (CFSE+7-AAD+ cells) was analyzed using the BD FACSCalibur flow cytometer. The percentage of specific lysis was calculated using the formula: (Cell death in sample − Spontaneous cell death) / (100 − Spontaneous cell death). In some experiments, CD8+ T cells were activated using anti-CD3–coated beads; anti-CD226 blocking Ab (DX11) or isotype control Ab was added.

The differences between two groups were analyzed using Student or Welch t test according to their distributions. Multiple comparisons were analyzed using Tukey-Kramer honestly significant difference test. Comparison between the same individual was performed using a paired t test. The relations between two continuous variables were analyzed using the Spearman rank correlation. All tests were two-tailed, and p values < 0.05 were considered significant. All analyses were performed using JMP statistical software (SAS Institute, Cary, NC).

To reveal the immune dysfunction of effector memory CD8+ T cells in SSc patients, we performed comprehensive gene expression analysis by cDNA microarray of the sort-purified effector memory CD8+ T cells from SSc patients (n = 9; Table I) and HCs (n = 5). We obtained 55 significantly altered genes (limma p < 0.05 and absolute fold change, > 2) between the samples of effector memory CD8+ T cells from SSc patients and those from HCs (Fig. 1A). We next completed DAVID functional annotation clustering analysis to explore the functional roles of the differentially expressed genes and found two significantly enriched functional clusters (defined as an Enrichment Score > 1.3; Supplemental Table I). One cluster contained eight genes associated with an immune response: CCL3-like 3 (CCL3L3); CCL4-like 1 (CCL4L1); IFN-γ; IFN-induced protein 44-like (IFI44L); phospholipase C, γ 2 (PLCG2); MHC class II DR α (HLA-DRA); IL-32; and complement factor D (adipsin) (CFD). The other cluster contained seven genes associated with cell adhesion: CCL4L1; G protein-coupled receptor 56 (GPR56); integrin, α M (complement component 3 receptor 3 subunit; ITGAM); IL-32; adhesion molecule, interacts with CXADR Ag 1 (AMICA1); ninjurin 1 (NINJ1); and fasciculation and elongation protein ζ 1 (zygin I) (FEZ1). Because the results indicate that the functions of immune response and cell adhesion were affected in SSc patients and may be critical in SSc pathogenesis, we also focused on the genes related to both clusters. In the 61 genes selected, the genes with a relatively high fold change were CCL4L1; CD300a molecule (CD300A); CD226 molecule (CD226); VCAM1; and spondin 2, extracellular matrix protein (SPON2; Fig. 1B).

Table I.
Clinical characteristics of SSc patients
FeatureMicroarray (n = 9)FACS (n = 41)p Values
Age, mean ± SD (y) 53.6 ± 12.3 57.9 ± 11.4 0.32 
Female, n (%) 9 (100) 37 (90) 
Disease duration, mean ± SD (y) 3.6 ± 3.3 7.0 ± 4.9 0.06 
Diffuse cutaneous SSc, n (%) 6 (67) 27 (66) 
MRSS, mean ± SD 16.4 ± 11.1 10.0 ± 6.7 0.19 
ILD, n (%) 7 (78) 23 (56) 0.28 
FVC, % predicted, mean ± SD 88.5 ± 10.3 85.9 ± 17.1 0.58 
DLCO, % predicted, mean ± SD 50.0 ± 12.6 66.4 ± 21.0 0.10 
Antinuclear Ab, n (%) 8 (89) 37 (90) 
Anti-topoisomerase I Ab, n (%) 5 (56) 15 (37) 0.45 
Anticentromere Ab, n (%) 1 (11) 8 (20) 0.66 
Low-dose corticosteroids, n (%) 6 (67) 20 (49) 0.47 
Immunosuppressive drug use, n (%) 2 (22) 9 (22) 
FeatureMicroarray (n = 9)FACS (n = 41)p Values
Age, mean ± SD (y) 53.6 ± 12.3 57.9 ± 11.4 0.32 
Female, n (%) 9 (100) 37 (90) 
Disease duration, mean ± SD (y) 3.6 ± 3.3 7.0 ± 4.9 0.06 
Diffuse cutaneous SSc, n (%) 6 (67) 27 (66) 
MRSS, mean ± SD 16.4 ± 11.1 10.0 ± 6.7 0.19 
ILD, n (%) 7 (78) 23 (56) 0.28 
FVC, % predicted, mean ± SD 88.5 ± 10.3 85.9 ± 17.1 0.58 
DLCO, % predicted, mean ± SD 50.0 ± 12.6 66.4 ± 21.0 0.10 
Antinuclear Ab, n (%) 8 (89) 37 (90) 
Anti-topoisomerase I Ab, n (%) 5 (56) 15 (37) 0.45 
Anticentromere Ab, n (%) 1 (11) 8 (20) 0.66 
Low-dose corticosteroids, n (%) 6 (67) 20 (49) 0.47 
Immunosuppressive drug use, n (%) 2 (22) 9 (22) 

DLCO, diffusing capacity for carbon monoxide; FVC, forced vital capacity.

FIGURE 1.

Gene expression profiles of effector memory CD8+ T cells from 9 SSc patients and 5 HCs. (A) The heat map generated by hierarchical clustering shows 55 significantly altered genes (limma p < 0.05; absolute fold change [FC] > 2) between the samples from SSc patients and those from HCs. Enriched functional clusters of these differentially expressed genes are shown in detailed in Supplemental Table I. *Genes annotated to cluster 1 (including immune response). Genes annotated to cluster 2 (including cell adhesion). (B) Top five upregulated or downregulated genes related to both clusters.

FIGURE 1.

Gene expression profiles of effector memory CD8+ T cells from 9 SSc patients and 5 HCs. (A) The heat map generated by hierarchical clustering shows 55 significantly altered genes (limma p < 0.05; absolute fold change [FC] > 2) between the samples from SSc patients and those from HCs. Enriched functional clusters of these differentially expressed genes are shown in detailed in Supplemental Table I. *Genes annotated to cluster 1 (including immune response). Genes annotated to cluster 2 (including cell adhesion). (B) Top five upregulated or downregulated genes related to both clusters.

Close modal

In this study, we attempted to validate CD226 expression at a protein level and to investigate the characteristic functions of CD226 on CD8+ T cells in SSc pathogenesis, because CD226 has been reported as a genetic risk factor for the susceptibility to SSc in genome-wide association studies (30).

To validate the results of the microarray analysis, we investigated CD226 expression by flow cytometry on CD8+ T cells from 41 SSc patients (mean age, 57.8 y; female, 37) and 23 HCs (mean age, 55.7 y; female, 19). No significant differences were found between SSc patients and HCs in terms of age and sex. As described in Table I, the clinical characteristics of the SSc patients investigated were generally similar between two analysis groups.

CD226 was constitutively expressed on CD8+ T cells (Fig. 2A). CD226 expression had no association with age, sex, or medication (data not shown). The proportion of CD226+CD8+ T cells was significantly higher in SSc patients than in HCs (Fig. 2A), whereas CD226 expression on CD4+ T cells, CD19+ B cells, NK cells, and γδ cells was almost the same between two groups (Fig. 2B–E). For further validation, we also analyzed CD226 expression on effector memory CD8+ T cells and found that CD226-expressing cells were increased in SSc patients compared with HCs, as expected from the microarray results (Fig. 2F).

FIGURE 2.

CD226 expression on CD8+ T cells is upregulated in SSc patients compared with HCs. CD226 expression on lymphocyte subsets from 41 SSc patients and 23 HCs was analyzed using flow cytometry. (A) Representative histograms show the percentage of CD226+ cells on gated CD3+CD8+ cells from SSc patients and HCs (left panels). The proportion of CD226+CD8+ T cells was compared between SSc patients and HCs (right panel). (BF) The frequency of CD226-expressing cells in (B) CD3+CD4+ T cells, (C) CD19+ B cells, (D) γδ cells, (E) NK cells, and (F) effector memory CD8+ T cells (CD8+CD45RO+CD62L+ cells) was compared between SSc patients and HCs. (G) Representative dot plots show CD45RO and CD62L expression on CD226CD8+, CD226lowCD8+, and CD226highCD8+ T cells. Numbers indicate the proportion of each population. (H) The percentage of CD226highCD8+ T cells was compared between SSc patients and HCs. Each data point represents a single subject. Horizontal lines show the mean. CM, central memory T cells; EM, effector memory T cells; EMRA, terminal differentiated effector memory T cells; naive, naive T cells.

FIGURE 2.

CD226 expression on CD8+ T cells is upregulated in SSc patients compared with HCs. CD226 expression on lymphocyte subsets from 41 SSc patients and 23 HCs was analyzed using flow cytometry. (A) Representative histograms show the percentage of CD226+ cells on gated CD3+CD8+ cells from SSc patients and HCs (left panels). The proportion of CD226+CD8+ T cells was compared between SSc patients and HCs (right panel). (BF) The frequency of CD226-expressing cells in (B) CD3+CD4+ T cells, (C) CD19+ B cells, (D) γδ cells, (E) NK cells, and (F) effector memory CD8+ T cells (CD8+CD45RO+CD62L+ cells) was compared between SSc patients and HCs. (G) Representative dot plots show CD45RO and CD62L expression on CD226CD8+, CD226lowCD8+, and CD226highCD8+ T cells. Numbers indicate the proportion of each population. (H) The percentage of CD226highCD8+ T cells was compared between SSc patients and HCs. Each data point represents a single subject. Horizontal lines show the mean. CM, central memory T cells; EM, effector memory T cells; EMRA, terminal differentiated effector memory T cells; naive, naive T cells.

Close modal

As shown in Fig. 2A and 2G, the CD226+CD8+ T cells were subdivided into two subpopulations on the basis of the intensity of expression of CD226 and differentiation markers (CD45RO and CD62L). CD8+ T cells with low intensity of CD226 (CD226lowCD8+ T cells) mainly contained naive cells; in contrast, CD8+ T cells with high (CD226highCD8+ T cells) were heterogeneous and included memory and effector cells. In this line, further analysis with attention to CD226highCD8+ T cells, which may be more pathogenic, revealed that the percentage of the cells was markedly greater in SSc patients than in HCs (Fig. 2H).

To study the roles of upregulated CD226 expression on CD8+ T cells in SSc patients, we investigated the association between CD226highCD8+ T cells and clinical manifestations of SSc. The percentage of CD226highCD8+ T cells had no association with disease duration (data not shown). Concerning the subgroups of SSc, the expansion of CD226highCD8+ T cells was significantly larger in patients with diffuse cutaneous SSc than in those with limited cutaneous SSc (Fig. 3A). SSc patients with anti-topoisomerase I Abs also showed a significantly higher proportion of CD226highCD8+ T cells than those with anticentromere Abs did (Fig. 3B). The percentage of CD226highCD8+ T cells was significantly greater in patients with moderate to severe skin sclerosis (MRSS, 15–39) than in those with mild (MRSS, 1–14; Fig. 3C). When we studied the association of CD226 expression on CD8+ T cells with ILD, which is one of the major organ involvement in SSc, the frequency of CD226highCD8+ T cells was higher in patients with extensive ILD than in those with limited ILD and those without ILD (Fig. 3D). These results indicate that CD226highCD8+ T cells are associated with the severity of skin and lung involvement.

FIGURE 3.

CD226highCD8+ T cells are increased in patients with diffuse cutaneous SSc (dcSSc) and are associated with the severity of skin and lung involvement. (A) The percentage of CD226highCD8+ T cells was compared between dcSSc patients (n = 27) and patients with limited cutaneous SSc (lcSSc; n = 14). (B) The proportion of CD226highCD8+ T cells was compared between patients with anti-topoisomerase I Abs (ATA; n = 15) and those with anticentromere Abs (ACA; n = 8). (C) The percentage of CD226highCD8+ T cells was compared between patients with mild skin sclerosis (MRSS, 1–14; n = 31) and those with moderate to severe (MRSS, 15–39; n = 10). (D) The frequency of CD226highCD8+ T cells was compared among SSc patients without ILD (ILD−; n = 18), those with limited ILD (n = 18), and those with extensive ILD (n = 5). Each data point represents a single subject. Horizontal lines show the mean.

FIGURE 3.

CD226highCD8+ T cells are increased in patients with diffuse cutaneous SSc (dcSSc) and are associated with the severity of skin and lung involvement. (A) The percentage of CD226highCD8+ T cells was compared between dcSSc patients (n = 27) and patients with limited cutaneous SSc (lcSSc; n = 14). (B) The proportion of CD226highCD8+ T cells was compared between patients with anti-topoisomerase I Abs (ATA; n = 15) and those with anticentromere Abs (ACA; n = 8). (C) The percentage of CD226highCD8+ T cells was compared between patients with mild skin sclerosis (MRSS, 1–14; n = 31) and those with moderate to severe (MRSS, 15–39; n = 10). (D) The frequency of CD226highCD8+ T cells was compared among SSc patients without ILD (ILD−; n = 18), those with limited ILD (n = 18), and those with extensive ILD (n = 5). Each data point represents a single subject. Horizontal lines show the mean.

Close modal

To investigate the roles of the elevated CD226 expression level of CD8+ T cells in SSc pathogenesis, we examined the functional properties of CD226+CD8+ T cells. We first compared cytokine production after in vitro stimulation with PMA/ionomycin between CD226+CD8+ and CD226CD8+ T cells using the CBA system. We confirmed that cytokine production was not affected by anti-CD226 Ab staining for cell sorting (data not shown). In HCs, the sort-purified CD226+CD8+ T cells showed higher levels of production of various cytokines (IFN-γ, IL-4, IL-5, and IL-13) than CD226CD8+ T cells did (Fig. 4A). Interestingly, IL-4, IL-5, and IL-13 were almost exclusively produced by CD226+CD8+ T cells, although these type 2 cytokines were not produced as much as type 1 cytokines (Fig. 4A). We next analyzed cell surface expression of CRTH2, the most reliable marker for the cells producing these type 2 cytokines (Tc2) (31), on CD8+ T cells and found that CRTH2+ cells were enriched in CD226highCD8+ T cells (Fig. 4B).

FIGURE 4.

IL-13 production of CD226highCD8+ T cells is higher than that of CD226lowCD8+ and CD226CD8+ T cells. (A) The sort-purified CD226CD8+ or CD226+CD8+ T cells were stimulated with PMA/ionomycin for 6 h. Cytokine production was measured in the supernatants using the CBA system. Bars are the mean ± SEM of 7 HCs. (B) Histograms show CD226 expression on CRTH2+CD3+CD8+ and CD3+CD8+ T cells. Data shown are representative of three HCs and three SSc patients. (C) Whole blood was stimulated with PMA/ionomycin for 5 h in the presence of GolgiStop. After RBC lysis, the cells were stained and analyzed using flow cytometry. Representative contour plots show IL-13 positive cells in CD226CD3+CD8+, CD226lowCD3+CD8+, and CD226highCD3+CD8+ T cells (left panel). Percentages indicate the proportion of each population. The proportion of IL-13 positive cells in each population was compared between 10 SSc patients and 12 HCs (right panel). (D) The sort-purified CD226CD8+, CD226lowCD8+, or CD226highCD8+ T cells were stimulated with PMA/ionomycin for 4 h. The mRNA expression levels of GATA3 and IL-13 were assessed by quantitative real-time PCR. Bars are the mean ± SEM of three HCs. (E) The purified CD8+ T cells were activated with anti-CD3–coated beads and cocultured with HUVECs (for costimulation through CD226/CD155 interactions) for 72 h in the presence of anti-CD226 blocking Ab or isotype control Ab. IL-13 and IFN-γ production was measured using the CBA system. Results were normalized to their respective isotype control. Bars are the mean ± SEM of four HCs.

FIGURE 4.

IL-13 production of CD226highCD8+ T cells is higher than that of CD226lowCD8+ and CD226CD8+ T cells. (A) The sort-purified CD226CD8+ or CD226+CD8+ T cells were stimulated with PMA/ionomycin for 6 h. Cytokine production was measured in the supernatants using the CBA system. Bars are the mean ± SEM of 7 HCs. (B) Histograms show CD226 expression on CRTH2+CD3+CD8+ and CD3+CD8+ T cells. Data shown are representative of three HCs and three SSc patients. (C) Whole blood was stimulated with PMA/ionomycin for 5 h in the presence of GolgiStop. After RBC lysis, the cells were stained and analyzed using flow cytometry. Representative contour plots show IL-13 positive cells in CD226CD3+CD8+, CD226lowCD3+CD8+, and CD226highCD3+CD8+ T cells (left panel). Percentages indicate the proportion of each population. The proportion of IL-13 positive cells in each population was compared between 10 SSc patients and 12 HCs (right panel). (D) The sort-purified CD226CD8+, CD226lowCD8+, or CD226highCD8+ T cells were stimulated with PMA/ionomycin for 4 h. The mRNA expression levels of GATA3 and IL-13 were assessed by quantitative real-time PCR. Bars are the mean ± SEM of three HCs. (E) The purified CD8+ T cells were activated with anti-CD3–coated beads and cocultured with HUVECs (for costimulation through CD226/CD155 interactions) for 72 h in the presence of anti-CD226 blocking Ab or isotype control Ab. IL-13 and IFN-γ production was measured using the CBA system. Results were normalized to their respective isotype control. Bars are the mean ± SEM of four HCs.

Close modal

We further performed a detailed analysis of IL-13, because it was reported as a profibrotic cytokine (32) and it had the highest levels among these type 2 cytokines. FACS analysis of cytokine production from CD8+ T cells also revealed that CD226high subpopulation was a major source of IL-13 in SSc patients and HCs (Fig. 4C). IL-13–producing CD226highCD8+ T cells were significantly increased in SSc patients compared with HCs (Fig. 4C). We then assessed whether CD226 was involved in the IL-13 production of CD8+ T cells. The mRNA expression analysis revealed that GATA3 (master transcription factor of Tc2) and IL-13 expression after in vitro stimulation with PMA/ionomycin were markedly higher in CD226highCD8+ T cells than in CD226lowCD8+ and CD226CD8+ T cells (Fig. 4D), suggesting that CD226highCD8+ T cells had more potent capacity to differentiate into Tc2 and produce IL-13. Next, to clarify the importance of costimulatory signals through CD226, we cocultured purified CD8+ T cells with HUVECs, which expressed high levels of CD155 (one of the CD226 ligands), in the presence or absence of anti-CD226 blocking Ab. In this assay, the blocking costimulatory signals through CD226 significantly suppressed effector functions and impaired the IL-13 and IFN-γ production of CD8+ T cells, measured using the CBA system (Fig. 4E). These findings suggest that CD226highCD8+ T cells are involved in SSc pathogenesis by producing abundant cytokines including IL-13 via CD226/CD155 costimulatory signaling.

To evaluate the effect of CD226 expression on the cytotoxic capacity of CD8+ T cells, we analyzed granzyme B production after in vitro stimulation with PMA/ionomycin by the CBA system. In HCs, CD226+CD8+ T cells showed higher levels of granzyme B production than CD226CD8+ T cells did, indicating that CD226+CD8+ T cells were highly cytotoxic populations (Fig. 5A).

FIGURE 5.

CD226 is involved in the upregulated cytotoxic capacity of CD8+ T cells in SSc patients. (A) The sort-purified CD226CD8+ or CD226+CD8+ T cells were stimulated with PMA/ionomycin for 6 h. Granzyme B production was measured in the supernatants using the CBA system. Bars are the mean ± SEM of seven HCs. (B) Purified CD8+ T cells were cocultured with CFSE-labeled HUVECs for 6 h at the E:T ratio of 10:1 or 40:1. Dead HUVECs were defined as CFSE+7-AAD+ cells by flow cytometry. The cytotoxic capacity of CD8+ T cells against HUVECs was compared between SSc patients and HCs. (C) Correlation between the cytotoxic capacity and the percentage of CD226highCD8+ T cells in SSc patients and HCs (left, E:T = 10:1; right, 40:1). (D) Purified CD8+ T cells were activated with anti-CD3–coated beads and cocultured with CFSE-labeled HUVECs for 6 h in the presence of anti-CD226 blocking Ab or isotype control Ab. Each data point represents a single subject (two SSc patients and three HCs). Horizontal lines show the mean.

FIGURE 5.

CD226 is involved in the upregulated cytotoxic capacity of CD8+ T cells in SSc patients. (A) The sort-purified CD226CD8+ or CD226+CD8+ T cells were stimulated with PMA/ionomycin for 6 h. Granzyme B production was measured in the supernatants using the CBA system. Bars are the mean ± SEM of seven HCs. (B) Purified CD8+ T cells were cocultured with CFSE-labeled HUVECs for 6 h at the E:T ratio of 10:1 or 40:1. Dead HUVECs were defined as CFSE+7-AAD+ cells by flow cytometry. The cytotoxic capacity of CD8+ T cells against HUVECs was compared between SSc patients and HCs. (C) Correlation between the cytotoxic capacity and the percentage of CD226highCD8+ T cells in SSc patients and HCs (left, E:T = 10:1; right, 40:1). (D) Purified CD8+ T cells were activated with anti-CD3–coated beads and cocultured with CFSE-labeled HUVECs for 6 h in the presence of anti-CD226 blocking Ab or isotype control Ab. Each data point represents a single subject (two SSc patients and three HCs). Horizontal lines show the mean.

Close modal

To confirm this phenotype, we performed a CD8+ T cell–mediated cytotoxicity assay. Because endothelial cell injury was considered as one of the earliest events of SSc (4), we selected HUVECs as the target of the CD8+ T cell–mediated cytotoxicity assay. The cytotoxic ability of CD8+ T cells was significantly upregulated in SSc patients compared with HCs at both E:T ratios tested (Fig. 5B). Furthermore, the cytotoxic capacity was positively correlated with the percentage of CD226highCD8+ T cells, and this correlation was remarkable in SSc patients (Fig. 5C). We next performed this assay in the presence or absence of anti-CD226 blocking Ab to establish whether CD226 was involved in HUVEC cytolysis. The cytotoxicity of anti-CD3–activated CD8+ T cells against HUVECs was significantly inhibited by CD226 neutralization (Fig. 5D). These data indicate that CD226 is involved in the upregulated cytotoxic capacity of CD8+ T cells in SSc patients.

In this study, we showed that genes related to immune response and cell adhesion were differently expressed in effector memory CD8+ T cells and that the expression of CD226, which is associated with both functions, is higher on CD8+ T cells from SSc patients than from HCs. We also demonstrated that CD226highCD8+ T cells are associated with skin and lung involvements and may be involved in SSc pathogenesis by upregulated production of cytokines, such as IL-13, and by cytotoxic capacity against endothelial cells.

CD226 is a member of the Ig-superfamily of receptors and is constitutively expressed on the cell surface of T cells, NK cells, monocytes, platelets, and a subset of B cells (33). It is also known as a costimulatory and adhesion molecule (33, 34) and mediates activation signals for cytotoxicity by CD8+ T cells and NK cells via its binding to CD155, CD112, or both on target cells (35, 36). The important function of CD226 is well known as the immune surveillance of tumors (37). Concerning autoimmune diseases such as SSc, the nonsynonymous rs763361 polymorphism in CD226 has been identified as a genetic risk factor in the Western white population (30, 38). More recently, Avouac et al. (39) reported that mice deficient for CD226 were protected from bleomycin-induced dermal fibrosis in animal models of SSc. In humans, however, the roles of CD226 in SSc pathogenesis remain unknown.

In this study, we showed that CD226+CD8+ T cells, especially CD226highCD8+ T cells, were expanded in SSc patients compared with HCs. The expansion was striking in patients with diffuse cutaneous SSc and in those with ILD. Furthermore, the expanded CD226highCD8+ T cells reflected the severity of skin and lung involvement. Recent reports about analysis of CD226 polymorphism have shown that the CD226 rs763361 T allele was associated with the diffuse cutaneous SSc subtype (30) and SSc-related fibrosing alveolitis subsets (30, 38). Although the level of CD226 expression was not influenced by the CD226 rs763361 genotype (30), this CD226 polymorphism substituted a glycine to a serine residue at amino acid position 307 in the intracellular domain and has a possibility to affect downstream signaling and cell functions (40). These findings suggest that the CD226 polymorphism is also associated with the characteristic functions of CD226+CD8+ T cells.

IL-13 is one of cytokines that mediate tissue fibrosis (32, 41) and plays a pivotal function in SSc pathogenesis (42, 43). In this study, we found that CD226+CD8+ T cells produced a large amount of cytokines including IL-13 and that the IL-13 production of CD226highCD8+ T cells was upregulated in SSc patients compared with HCs. The roles of CD226 in the functions of CD8+ T cells were well understood about type 1 cytokine production and cytotoxicity (33), but not about type 2 cytokine production. In the current study, CD226+CD8+ T cells produced higher amount of type 2 cytokines (IL-4, IL-5, and IL-13) as well as type 1 cytokines (IFN-γ and TNF-α) than CD226-negative ones. Furthermore, CD226highCD8+ T cells had potent capacity to differentiate into Tc2, and the blockage of CD226 signaling inhibited IL-13 and IFN-γ production. These results indicate that CD226 is an important costimulator of CD8+ T cells in producing not only IFN-γ, but also IL-13. Recent studies about CD4+ T cells revealed that CD226 was involved in Th1 differentiation (44, 45), that IL-13 production did not correlate with CD226 expression (46), and that the blockage of CD226 signaling did not affect IL-13 production (46). Although this discrepancy might reflect the difference between CD8+ T cells and CD4+ T cells, further study is required to investigate the detailed mechanism of the IL-13 production of CD8+ T cells.

Endothelial cell apoptosis could be a primary event in the pathogenesis of SSc (4). NK cells (47) and γδ cells (48) were reported as the mediators of the cytotoxicity; however, it remains unclear whether CD8+ T cells are involved in endothelial cell injury. On this point, we observed that the cytotoxic capacity of CD8+ T cells against endothelial cells was upregulated in SSc patients. In addition, the cytotoxic capacity was correlated with CD226 expression on CD8+ T cells and partially suppressed by blocking CD226. In CD226-deficient mouse, activated CD226−/− T cells showed impaired cytotoxic activity because of dysfunctional cell–cell contact required for the efficient killing of target cells (49), and CD8+ T cells required CD226 for costimulation when recognizing Ag presented by nonprofessional APCs (50). Our results and these findings indicate that interaction between CD226 on CD8+ T cells and its ligand CD155 on endothelial cells is critical for both cell–cell contact and costimulatory signals. Upregulated interaction of CD8+ T cells and endothelial cells was supported by the recent report that increased expression of CD226 on perivascular inflammatory cells (e.g., T cells) in the lesional skin of SSc patients (39).

The limitations of this study include small sample size, cross-sectional design, and race bias. Further investigation with larger sample size and various populations is needed to confirm the results of this study.

In conclusion, upregulated CD226 expression on CD8+ T cells is associated with skin and lung involvement, reflects disease severity, and may be involved in SSc pathogenesis via the production of cytokines such as profibrotic IL-13 and endothelial cell injury. Because blockage of CD226 impairs the IL-13 production and cytotoxic capacity of CD8+ T cells, CD226 may be a useful target in the treatment of SSc.

This work was supported by grants from the Ministry of Health, Labour, and Welfare (H24-Jitsuyoka-Kokusai-004) and the Ministry of Education, Culture, Sports, Science and Technology of Japan (Japan Society for the Promotion of Science, KAKENHI Grants 24591465 and 15K09527).

The microarray data presented in this article have been submitted to the National Center for Biotechnology Information Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) under accession number GSE63903.

The online version of this article contains supplemental material.

Abbreviations used in this article:

7-AAD

7-aminoactinomycin D

CBA

cytometric beads array

CRTH2

chemoattractant receptor-homologous molecule on Th2 cells

DAVID

Database for Annotation, Visualization, and Integrated Discovery

FVC

forced vital capacity

GATA3

GATA binding protein 3

HC

healthy control

HRCT

high-resolution computed tomography

ILD

interstitial lung disease

limma

linear models for microarray analysis

MRSS

modified Rodnan skin thickness score

SSc

systemic sclerosis.

1
Gabrielli
A.
,
Avvedimento
E. V.
,
Krieg
T.
.
2009
.
Scleroderma.
N. Engl. J. Med.
360
:
1989
2003
.
2
Gu
Y. S.
,
Kong
J.
,
Cheema
G. S.
,
Keen
C. L.
,
Wick
G.
,
Gershwin
M. E.
.
2008
.
The immunobiology of systemic sclerosis.
Semin. Arthritis Rheum.
38
:
132
160
.
3
Abraham
D. J.
,
Varga
J.
.
2005
.
Scleroderma: from cell and molecular mechanisms to disease models.
Trends Immunol.
26
:
587
595
.
4
Sgonc
R.
,
Gruschwitz
M. S.
,
Dietrich
H.
,
Recheis
H.
,
Gershwin
M. E.
,
Wick
G.
.
1996
.
Endothelial cell apoptosis is a primary pathogenetic event underlying skin lesions in avian and human scleroderma.
J. Clin. Invest.
98
:
785
792
.
5
Baraut
J.
,
Michel
L.
,
Verrecchia
F.
,
Farge
D.
.
2010
.
Relationship between cytokine profiles and clinical outcomes in patients with systemic sclerosis.
Autoimmun. Rev.
10
:
65
73
.
6
Fuschiotti
P.
,
Larregina
A. T.
,
Ho
J.
,
Feghali-Bostwick
C.
,
Medsger
T. A.
 Jr.
2013
.
Interleukin-13-producing CD8+ T cells mediate dermal fibrosis in patients with systemic sclerosis.
Arthritis Rheum.
65
:
236
246
.
7
Luzina
I. G.
,
Atamas
S. P.
,
Wise
R.
,
Wigley
F. M.
,
Choi
J.
,
Xiao
H. Q.
,
White
B.
.
2003
.
Occurrence of an activated, profibrotic pattern of gene expression in lung CD8+ T cells from scleroderma patients.
Arthritis Rheum.
48
:
2262
2274
.
8
Lympany
P. A.
,
Southcott
A. M.
,
Welsh
K. I.
,
Black
C. M.
,
Boylston
A. W.
,
du Bois
R. M.
.
1999
.
T-cell receptor gene usage in patients with fibrosing alveolitis and control subjects.
Eur. J. Clin. Invest.
29
:
173
181
.
9
Gustafsson
R.
,
Tötterman
T. H.
,
Klareskog
L.
,
Hällgren
R.
.
1990
.
Increase in activated T cells and reduction in suppressor inducer T cells in systemic sclerosis.
Ann. Rheum. Dis.
49
:
40
45
.
10
Tiev
K. P.
,
Abriol
J.
,
Burland
M. C.
,
Antonelli
D.
,
Klatzmann
D.
,
Cabane
J.
,
Boyer
O.
.
2005
.
T cell repertoire in patients with stable scleroderma.
Clin. Exp. Immunol.
139
:
348
354
.
11
Sakkas
L. I.
,
Platsoucas
C. D.
.
2004
.
Is systemic sclerosis an antigen-driven T cell disease?
Arthritis Rheum.
50
:
1721
1733
.
12
Fuschiotti
P.
,
Medsger
T. A.
 Jr.
,
Morel
P. A.
.
2009
.
Effector CD8+ T cells in systemic sclerosis patients produce abnormally high levels of interleukin-13 associated with increased skin fibrosis.
Arthritis Rheum.
60
:
1119
1128
.
13
Sallusto
F.
,
Lenig
D.
,
Förster
R.
,
Lipp
M.
,
Lanzavecchia
A.
.
1999
.
Two subsets of memory T lymphocytes with distinct homing potentials and effector functions.
Nature
401
:
708
712
.
14
Sallusto
F.
,
Geginat
J.
,
Lanzavecchia
A.
.
2004
.
Central memory and effector memory T cell subsets: function, generation, and maintenance.
Annu. Rev. Immunol.
22
:
745
763
.
15
Oertelt
S.
,
Selmi
C.
,
Invernizzi
P.
,
Podda
M.
,
Gershwin
M. E.
.
2005
.
Genes and goals: an approach to microarray analysis in autoimmunity.
Autoimmun. Rev.
4
:
414
422
.
16
Lyons
P. A.
,
Koukoulaki
M.
,
Hatton
A.
,
Doggett
K.
,
Woffendin
H. B.
,
Chaudhry
A. N.
,
Smith
K. G. C.
.
2007
.
Microarray analysis of human leucocyte subsets: the advantages of positive selection and rapid purification.
BMC Genomics
8
:
64
.
17
Masi
A. T.
,
Rodnan
G. P.
,
Medsger
T. A.
.
1980
.
Preliminary criteria for the classification of systemic sclerosis (scleroderma). Subcommittee for scleroderma criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee.
Arthritis Rheum.
23
:
581
590
.
18
LeRoy
E. C.
,
Black
C.
,
Fleischmajer
R.
,
Jablonska
S.
,
Krieg
T.
,
Medsger
T. A.
 Jr.
,
Rowell
N.
,
Wollheim
F.
.
1988
.
Scleroderma (systemic sclerosis): classification, subsets and pathogenesis.
J. Rheumatol.
15
:
202
205
.
19
Clements
P. J.
,
Lachenbruch
P. A.
,
Seibold
J. R.
,
Zee
B.
,
Steen
V. D.
,
Brennan
P.
,
Silman
A. J.
,
Allegar
N.
,
Varga
J.
,
Massa
M.
, et al
.
1993
.
Skin thickness score in systemic sclerosis: an assessment of interobserver variability in 3 independent studies.
J. Rheumatol.
20
:
1892
1896
.
20
Medsger
T. A.
 Jr.
,
Bombardieri
S.
,
Czirjak
L.
,
Scorza
R.
,
Della Rossa
A.
,
Bencivelli
W.
.
2003
.
Assessment of disease severity and prognosis.
Clin. Exp. Rheumatol.
21
(
3
,
Suppl 29
)
S42
S46
.
21
Goh
N. S. L.
,
Desai
S. R.
,
Veeraraghavan
S.
,
Hansell
D. M.
,
Copley
S. J.
,
Maher
T. M.
,
Corte
T. J.
,
Sander
C. R.
,
Ratoff
J.
,
Devaraj
A.
, et al
.
2008
.
Interstitial lung disease in systemic sclerosis: a simple staging system.
Am. J. Respir. Crit. Care Med.
177
:
1248
1254
.
22
Sasaki
E.
,
Nakamura
M.
,
Kida
K.
.
2001
.
Reference values for spirogram and blood gas analysis in Japanese non-smoking healthy adults.
J. Jpn. Respir. Soc.
39
:
383
399
.
23
Bolstad
B. M.
,
Irizarry
R. A.
,
Astrand
M.
,
Speed
T. P.
.
2003
.
A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.
Bioinformatics
19
:
185
193
.
24
Gentleman
R. C.
,
Carey
V. J.
,
Bates
D. M.
,
Bolstad
B.
,
Dettling
M.
,
Dudoit
S.
,
Ellis
B.
,
Gautier
L.
,
Ge
Y.
,
Gentry
J.
, et al
.
2004
.
Bioconductor: open software development for computational biology and bioinformatics.
Genome Biol.
5
:
R80
.
25
Smyth
G. K.
2005
.
Limma: linear models for microarray data
. In
Bioinformatics and computational biology solutions using R and Bioconductor.
Gentleman
R.
,
Carey
V. J.
,
Huber
W.
,
Irizarry
R. A.
,
Dudoit
S.
, eds.
Springer
,
New York, NY
, p.
397
420
.
26
Huang
D. W.
,
Sherman
B. T.
,
Tan
Q.
,
Collins
J. R.
,
Alvord
W. G.
,
Roayaei
J.
,
Stephens
R.
,
Baseler
M. W.
,
Lane
H. C.
,
Lempicki
R. A.
.
2007
.
The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists.
Genome Biol.
8
:
R183
.
27
Huang
W.
,
Sherman
B. T.
,
Lempicki
R. A.
.
2009
.
Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.
Nat. Protoc.
4
:
44
57
.
28
Saeed
A. I.
,
Sharov
V.
,
White
J.
,
Li
J.
,
Liang
W.
,
Bhagabati
N.
,
Braisted
J.
,
Klapa
M.
,
Currier
T.
,
Thiagarajan
M.
, et al
.
2003
.
TM4: a free, open-source system for microarray data management and analysis.
Biotechniques
34
:
374
378
.
29
Lecoeur
H.
,
Février
M.
,
Garcia
S.
,
Rivière
Y.
,
Gougeon
M. L.
.
2001
.
A novel flow cytometric assay for quantitation and multiparametric characterization of cell-mediated cytotoxicity.
J. Immunol. Methods
253
:
177
187
.
30
Dieudé
P.
,
Guedj
M.
,
Truchetet
M. E.
,
Wipff
J.
,
Revillod
L.
,
Riemekasten
G.
,
Matucci-Cerinic
M.
,
Melchers
I.
,
Hachulla
E.
,
Airo
P.
, et al
.
2011
.
Association of the CD226 Ser(307) variant with systemic sclerosis: evidence of a contribution of costimulation pathways in systemic sclerosis pathogenesis.
Arthritis Rheum.
63
:
1097
1105
.
31
Cosmi
L.
,
Annunziato
F.
,
Galli
M.
,
Maggi
R.
,
Nagata
K.
,
Romagnani
S.
.
2000
.
CRTH2 is the most reliable marker for the detection of circulating human type 2 Th and type 2 T cytotoxic cells in health and disease.
Eur. J. Immunol.
30
:
2972
2979
.
32
Wynn
T. A.
2003
.
IL-13 effector functions.
Annu. Rev. Immunol.
21
:
425
456
.
33
Shibuya
A.
,
Campbell
D.
,
Hannum
C.
,
Yssel
H.
,
Franz-Bacon
K.
,
McClanahan
T.
,
Kitamura
T.
,
Nicholl
J.
,
Sutherland
G. R.
,
Lanier
L. L.
,
Phillips
J. H.
.
1996
.
DNAM-1, a novel adhesion molecule involved in the cytolytic function of T lymphocytes.
Immunity
4
:
573
581
.
34
Shibuya
A.
,
Tahara-Hanaoka
S.
,
Shibuya
K.
.
2005
.
DNAM-1 (CD226): A Two-Sword Fencer for Innate and Adaptive Immunity.
Curr. Med. Chem.
4
:
53
58
.
35
Bottino
C.
,
Castriconi
R.
,
Pende
D.
,
Rivera
P.
,
Nanni
M.
,
Carnemolla
B.
,
Cantoni
C.
,
Grassi
J.
,
Marcenaro
S.
,
Reymond
N.
, et al
.
2003
.
Identification of PVR (CD155) and Nectin-2 (CD112) as cell surface ligands for the human DNAM-1 (CD226) activating molecule.
J. Exp. Med.
198
:
557
567
.
36
Tahara-Hanaoka
S.
,
Shibuya
K.
,
Onoda
Y.
,
Zhang
H.
,
Yamazaki
S.
,
Miyamoto
A.
,
Honda
S.
,
Lanier
L. L.
,
Shibuya
A.
.
2004
.
Functional characterization of DNAM-1 (CD226) interaction with its ligands PVR (CD155) and nectin-2 (PRR-2/CD112).
Int. Immunol.
16
:
533
538
.
37
Chan
C. J.
,
Andrews
D. M.
,
Smyth
M. J.
.
2012
.
Receptors that interact with nectin and nectin-like proteins in the immunosurveillance and immunotherapy of cancer.
Curr. Opin. Immunol.
24
:
246
251
.
38
Bossini-Castillo
L.
,
Simeon
C. P.
,
Beretta
L.
,
Broen
J. C.
,
Vonk
M. C.
,
Ríos-Fernández
R.
,
Espinosa
G.
,
Carreira
P.
,
Camps
M. T.
,
Castillo
M. J.
, et al
Spanish Scleroderma Group
.
2012
.
A multicenter study confirms CD226 gene association with systemic sclerosis-related pulmonary fibrosis.
Arthritis Res. Ther.
14
:
R85
.
39
Avouac
J.
,
Elhai
M.
,
Tomcik
M.
,
Ruiz
B.
,
Friese
M.
,
Piedavent
M.
,
Colonna
M.
,
Bernhardt
G.
,
Kahan
A.
,
Chiocchia
G.
, et al
.
2013
.
Critical role of the adhesion receptor DNAX accessory molecule-1 (DNAM-1) in the development of inflammation-driven dermal fibrosis in a mouse model of systemic sclerosis.
Ann. Rheum. Dis.
72
:
1089
1098
.
40
Maier
L. M.
,
Hafler
D. A.
.
2009
.
Autoimmunity risk alleles in costimulation pathways.
Immunol. Rev.
229
:
322
336
.
41
Wynn
T. A.
2008
.
Cellular and molecular mechanisms of fibrosis.
J. Pathol.
214
:
199
210
.
42
Fuschiotti
P.
2011
.
Role of IL-13 in systemic sclerosis.
Cytokine
56
:
544
549
.
43
Greenblatt
M. B.
,
Aliprantis
A. O.
.
2013
.
The immune pathogenesis of scleroderma: context is everything.
Curr. Rheumatol. Rep.
15
:
297
.
44
Shibuya
K.
,
Shirakawa
J.
,
Kameyama
T.
,
Honda
S.
,
Tahara-Hanaoka
S.
,
Miyamoto
A.
,
Onodera
M.
,
Sumida
T.
,
Nakauchi
H.
,
Miyoshi
H.
,
Shibuya
A.
.
2003
.
CD226 (DNAM-1) is involved in lymphocyte function-associated antigen 1 costimulatory signal for naive T cell differentiation and proliferation.
J. Exp. Med.
198
:
1829
1839
.
45
Dardalhon
V.
,
Schubart
A. S.
,
Reddy
J.
,
Meyers
J. H.
,
Monney
L.
,
Sabatos
C. A.
,
Ahuja
R.
,
Nguyen
K.
,
Freeman
G. J.
,
Greenfield
E. A.
, et al
.
2005
.
CD226 is specifically expressed on the surface of Th1 cells and regulates their expansion and effector functions.
J. Immunol.
175
:
1558
1565
.
46
Lozano
E.
,
Joller
N.
,
Cao
Y.
,
Kuchroo
V. K.
,
Hafler
D. A.
.
2013
.
The CD226/CD155 interaction regulates the proinflammatory (Th1/Th17)/anti-inflammatory (Th2) balance in humans.
J. Immunol.
191
:
3673
3680
.
47
Sgonc
R.
,
Gruschwitz
M. S.
,
Boeck
G.
,
Sepp
N.
,
Gruber
J.
,
Wick
G.
.
2000
.
Endothelial cell apoptosis in systemic sclerosis is induced by antibody-dependent cell-mediated cytotoxicity via CD95.
Arthritis Rheum.
43
:
2550
2562
.
48
Kahaleh
M. B.
,
Fan
P. S.
,
Otsuka
T.
.
1999
.
Gammadelta receptor bearing T cells in scleroderma: enhanced interaction with vascular endothelial cells in vitro.
Clin. Immunol.
91
:
188
195
.
49
Ramsbottom
K. M.
,
Hawkins
E. D.
,
Shimoni
R.
,
McGrath
M.
,
Chan
C. J.
,
Russell
S. M.
,
Smyth
M. J.
,
Oliaro
J.
.
2014
.
Cutting edge: DNAX accessory molecule 1-deficient CD8+ T cells display immunological synapse defects that impair antitumor immunity.
J. Immunol.
192
:
553
557
.
50
Gilfillan
S.
,
Chan
C. J.
,
Cella
M.
,
Haynes
N. M.
,
Rapaport
A. S.
,
Boles
K. S.
,
Andrews
D. M.
,
Smyth
M. J.
,
Colonna
M.
.
2008
.
DNAM-1 promotes activation of cytotoxic lymphocytes by nonprofessional antigen-presenting cells and tumors.
J. Exp. Med.
205
:
2965
2973
.

The authors have no financial conflicts of interest.

Supplementary data