Abstract
In pulmonary sarcoidosis, CD4+ T-cells expressing T-cell receptor Vα2.3 accumulate in the lungs of HLA-DRB1*03+ patients. To investigate T-cell receptor-HLA-DRB1*03 interactions underlying recognition of hitherto unknown antigens, we performed detailed analyses of T-cell receptor expression on bronchoalveolar lavage fluid CD4+ T-cells from sarcoidosis patients.
Pulmonary sarcoidosis patients (n=43) underwent bronchoscopy with bronchoalveolar lavage. T-cell receptor α and β chains of CD4+ T-cells were analysed by flow cytometry, DNA-sequenced, and three-dimensional molecular models of T-cell receptor-HLA-DRB1*03 complexes generated.
Simultaneous expression of Vα2.3 with the Vβ22 chain was identified in the lungs of all HLA-DRB1*03+ patients. Accumulated Vα2.3/Vβ22-expressing T-cells were highly clonal, with identical or near-identical Vα2.3 chain sequences and inter-patient similarities in Vβ22 chain amino acid distribution. Molecular modelling revealed specific T-cell receptor-HLA-DRB1*03-peptide interactions, with a previously identified, sarcoidosis-associated vimentin peptide, (Vim)429–443 DSLPLVDTHSKRTLL, matching both the HLA peptide-binding cleft and distinct T-cell receptor features perfectly.
We demonstrate, for the first time, the accumulation of large clonal populations of specific Vα2.3/Vβ22 T-cell receptor-expressing CD4+ T-cells in the lungs of HLA-DRB1*03+ sarcoidosis patients. Several distinct contact points between Vα2.3/Vβ22 receptors and HLA-DRB1*03 molecules suggest presentation of prototypic vimentin-derived peptides.
Abstract
Clonal CD4+ lung T-cells associating with HLA-DRB1*03 molecules indicate specific antigens in pulmonary sarcoidosis http://ow.ly/UB81x
Introduction
CD4+ T-cells use clonally distributed αβ T-cell receptors (TCRs) to recognise antigenic peptides bound to HLA-DR, -DP and –DQ molecules on the surface of antigen-presenting cells. The random combination of variable (V), joining (J) and, for the β chain, diversity (D) gene segments and the similarly arbitrary addition and/or deletion of nucleotides at V–D–J junctions facilitate an extreme molecular diversity that underlies the unique antigen specificity of TCRs. Hypervariability primarily occurs in the complementarity-determining region (CDR) 3 loops, which constitute the main contact points for antigenic peptides presented by HLA molecules. V gene segments also contain germline-encoded CDR1 and CDR2 loops, which further mediate TCR recognition of peptide–HLA complexes by binding to the framework regions of HLA molecules, outside the peptide-binding groove [1].
In sarcoidosis, characteristic lung accumulations of T-cells expressing distinct TCR Vα or Vβ genes suggest the presence of specific antigens [2–4]. Lung-compartmentalised CD4+ T-cells are highly activated [5] and produce T-helper (Th) type 1 cytokines [6, 7]. More recently, an important role for Th17 [8] and regulatory T-cells Treg cells [9, 10] has also been implicated. We previously established that Vα2.3+ CD4+ T-cells always accumulate in the lungs of HLA-DRB1*0301+ and HLA-DRB3*0101+ [11] sarcoidosis patients with active disease, but not in patients with other inflammatory pulmonary diseases, such as allergic alveolitis or asthma, nor in healthy individuals [12, 13]. Earlier analyses of the TCR α chain of bronchoalveolar lavage (BAL) Vα2.3+ CD4+ T-cells indicated selection towards specific antigens [14]. Interestingly, HLA-DRB1*0301 and HLA-DRB3*0101 molecules are structurally and functionally similar and seem to present a highly similar peptide antigen repertoire, including identical epitopes [15]. Hence, it is plausible that sarcoidosis-specific antigen(s) could be presented by either HLA-DRB1*0301 or HLA-DRB3*0101 molecules, which may in turn be recognised by Vα2.3+ T-cells [11].
In order to investigate the TCR Vβ repertoire in BAL, we recently applied a panel of 21 TCR Vβ-specific antibodies, covering 70% of all Vβ gene segments, in sarcoidosis patients and healthy controls. Increased frequencies of T-cells expressing Vβ8, Vβ16, or Vβ22 were observed in the lungs of patients, and an association between lung-restricted Vβ22 expansions and HLA-DRB1*03+ was identified [16].
Here, we therefore performed a comprehensive analysis of the association between HLA-DRB1*03+ molecules and lung-accumulated Vα2.3+Vβ22+ CD4+ T-cells. Simultaneous expression of Vα2.3 and Vβ22 in BAL was assessed by flow cytometry, and such TCRs were further analysed by next-generation sequencing (NGS). We present a three-dimensional molecular model demonstrating hypothesised interactions between Vα2.3/Vβ22 TCRs and HLA-DRB1*03 molecules, in line with specific antigen presentation. Finally, by combining sequencing data and molecular modelling, we are able to propose a potential candidate antigen in the form of cytoskeletal protein vimentin, a finding that implicates autoimmune processes in sarcoidosis aetiology.
Materials and methods
Subjects, bronchoscopy and BAL
Bronchoscopy with BAL was performed as previously described [17] on 43 patients (16 females) with recent disease onset and a median age of 39 years (table 1). There was a preference for inclusion of patients with Löfgren's syndrome, as they are more often HLA-DRB1*03+, which associates strongly with lung accumulations of T-cells expressing the TCR Vα2.3 gene segment [13]. Löfgren's syndrome was defined in patients with an acute onset, usually with fever, chest radiographic findings with bilateral hilar lymphadenopathy, sometimes with pulmonary infiltrates, and with erythema nodosum and/or bilateral ankle arthritis. Thus, 26 patients had Löfgren's syndrome and 17 were classified as “non-Löfgren's syndrome” patients. All patients were HLA-typed (described in detail in the online data supplement) and diagnosed with sarcoidosis through typical clinical and radiographic manifestations, findings at bronchoscopy with BAL including an elevated CD4/CD8-ratio and positive biopsies (not required for diagnosing Löfgren's syndrome), in accordance with the criteria of the World Association of Sarcoidosis and other Granulomatous Disorders (WASOG) [18]. Written informed consent was obtained from all subjects and the Regional Ethical Review Board in Stockholm approved the study.
Clinical characteristics of the patient cohort with data shown for the entire patient group, as well as separately for HLA-DRB1*03 positive and negative patients
Sorting of Vα2.3+Vβ22+ CD4+ BAL T-cells and TCR α and β gene amplification
Phenotyping of Vα2.3+Vβ22+ T-cells by flow cytometry is described in the online data supplement. Expansions of lung T-cells expressing Vα2.3 and Vβ22 were defined based on their reference values of 3.5% [11] and 4.0% [16], respectively, in healthy individuals. The normal value of their co-expression was calculated to a median of 0.14% and the cut-off value for defining a Vα2.3+Vβ22+ T-cell expansion was set to 3×0.14=0.42%, thus approximated 0.5%. In a subset of patients with Vα2.3+Vβ22+ CD4+ BAL T-cell expansions, Vα2.3+Vβ22+, Vα2.3+Vβ22− or Vα2.3−Vβ22+ CD4+ BAL cells were isolated through fluorescence-activated cell sorting (FACS), using a MoFlo XDP cell sorter (Beckman Coulter; Brea, CA, USA). In a second subset of patients, we opted to use unsorted BAL cells as a means of avoiding cell distortion and mRNA damage that may occur during sorting. Due to the notable expansions of Vα2.3+Vβ22+ CD4+ T-cells, TCR Vα2.3 and Vβ22 chain mRNA should be present at high enough levels to be obtainable irrespective of sorting. The highly similar results showed both approaches (i.e. with or without prior FACS sorting) to be equally successful.
PCR amplification of Vα2.3 and Vβ22 genomic regions, respectively, was performed using a specific variable region oligonucleotide forward primer and a conserved constant region reverse primer for each TCR chain. Detailed mRNA extraction protocols, PCR conditions and primer information are provided in the online data supplement.
Statistical analysis
Relative numbers of Vα2.3+ and Vβ22+ BAL T-cells (expressed as percent of all CD4+ BAL T-cells) were calculated and compared between patient groups based on the presence or absence of distinct HLA-DRB1 and HLA-DRB3 alleles, respectively, using the two-tailed non-parametric Mann-Whitney U-test. A two-tailed Wilcoxon's signed rank test was applied for paired comparisons of Vα2.3+Vβ22+ versus Vα2.3−Vβ22− CD4+ T-cells. Analyses were performed using GraphPad Prism v.5.02 software (GraphPad Software Inc., La Jolla, CA, USA), and p<0.05 was considered significant.
DNA-sequencing and de novo transcript assembly
NGS was performed with the Illumina Nextera XT kit (Illumina, San Diego, CA, USA) according to protocol except for necessary modifications for adaption to the Agilent NGS workstation with MiSeq 2×250 sequencing (Agilent, Santa Clara, CA, USA). Additional details regarding the method and subsequent bioinformatics analysis is provided in the online data supplement. De novo transcript assembly and primary transcript analysis were performed using Trinity [19, 20].
Three-dimensional molecular modelling of ternary TCR/DRB1*0301/peptide complexes
Vα2.3/Vβ22 TCR sequences were retrieved from the IMGT (International ImMunoGeneTics Information System) database [21], and all modelling was performed manually using the COOT program [22]. A detailed description of the modelling procedure is provided in the online data supplement.
Results
TCRs of lung-accumulated T-cells in HLA-DRB1*03+ sarcoidosis patients
In HLA-DRB1*03+ patients, we verified previous findings of large populations of lung-accumulated CD4+ T-cells expressing the TCR Vα2.3 gene segment, constituting 33.3% (median) of all CD4+ BAL T-cells (figure 1a and table 2). In HLA-DRB1*03− patients, the corresponding median value was 5.2%, including six patients with close to 20% of BAL CD4+ T-cells expressing Vα2.3+ (figure 1a and table S2). These patients were found to be HLA-DRB1*03− but HLA-DRB3*01+ (tables S1 and S2). TCR Vβ22+ T-cells were also found to accumulate in the lungs of HLA-DRB1*03+ patients, with a median value of 9.2% of all BAL CD4+ T-cells, compared with 2.3% in HLA-DRB1*03− patients (table 2).
Vα2.3+Vβ22+ CD4+ T-cells accumulate in the lungs of HLA-DRB1*03+ patients with pulmonary sarcoidosis. a, b) Relative numbers of bronchoalveolar lavage (BAL) CD4+ T-cells expressing T-cell receptor (TCR) Vα2.3 (a) or TCR Vα2.3/Vβ22 (b) in DRB1*03+ (n=26) or DRB1*03− patients (n=17), respectively. ***: p<0.0001 (Mann–Whitney U test). c) Comparison of the percentage of CD4+ T-cells in BAL fluid that are Vα2.3+ Vβ22+ in DRB1*03+ patients (1–26, open bars) and in DRB1*03− patients (1–17, black bars), respectively, showing a significantly higher frequency of Vα2.3+ Vβ22+ CD4+ T-cells in HLA-DRB1*03+ than HLA-DRB1*03− patients. d) Representative contour plot of the Vα2.3+ Vβ22+ T-cell population in a DRB1*03+ sarcoidosis patient, as detected by multi-colour flow cytometry (Becton-Dickinson, San Jose, CA, USA). Cells are gated on CD3+ CD4+ lymphocytes, and populations were distributed as follows: Q1 20.8%, Q2 5.3%, Q3 71.0% and Q4 2.9%.
Summary of all bronchoalveolar lavage (BAL) CD4+ T-cells that express Vα2.3, Vβ22, and Vα2.3 together with Vβ22; the percentage of Vα2.3+ CD4+ BAL T-cells that express Vβ22, and the percentage of Vβ22+ CD4+ BAL T-cells that express Vα2.3 are also stated
Simultaneous expression of Vα2.3 and Vβ22 gene segments was found in 6.2% of all BAL CD4+ T-cells in HLA-DRB1*03+, compared with 0.3% in HLA-DRB1*03− patients (figure 1b and c, and table 2). No difference could be observed between men and women in terms of Vα2.3/Vβ22 expression. Moderate Vα2.3+Vβ22+ BAL T-cell expansions were found in some HLA-DRB1*03− patients (figure 1b and c); these patients were found to be DRB1*03− DRB3*01+. Relative percentages of Vα2.3+Vβ22+ CD4+ BAL T-cells in DRB1*03+, DRB1*03−/DRB3*01+, and DRB1*03− DRB3*01− patients, respectively, are presented in figure S1 and table S2. In HLA-DRB1*03+ patients, 21.0% of Vα2.3+ lung T-cells simultaneously expressed Vβ22, while 62.3% of Vβ22+ T-cells also expressed Vα2.3 (table 2 and figure S2). Representative flow cytometry staining plots of Vα2.3+Vβ22+ T-cells are shown in figure 1d and figure S2.
Vα2.3+Vβ22+ BAL T-cells differ with regard to phenotype
The lung-accumulated Vα2.3+Vβ22+ T-cell population expressed significantly higher levels of lymphocyte activation marker CD69 compared to Vα2.3−Vβ22− T-cells within the same individual, as assessed by flow cytometry (figure 2a). Furthermore, significantly fewer Vα2.3+Vβ22+ T-cells expressed CD27, which is lost upon prolonged activation and differentiation, compared to autologous Vα2.3−Vβ22− T-cells (figure 2b).
Vα2.3+Vβ22+ CD4+ T-cells have a higher activation and differentiation state than autologous Vα2.3− Vβ22− CD4+ T-cells. Expression of a) CD69 and b) CD27 in Vα2.3+Vβ22+ compared to Vα2.3−Vβ22− CD4+ bronchoalveolar lavage (BAL) T-cells from the same individual (n=12). Populations are expressed as percent of CD3+ CD4+ Vα2.3+Vβ22+ or Vα2.3−Vβ22− BAL lymphocytes. #: p=0.0049; ***: p=0.001 using Wilcoxon's signed rank test.
Vα2.3+Vβ22+ BAL T-cells are derived from a few dominant clones
TCR α and β chains were sequenced using specific primers for Vα2.3 and Vβ22 in combination with Cα and Cβ chains, respectively. In patients 1–4 and 8–10, lung-accumulated Vα2.3+Vβ22+ CD4+ T-cell populations were FACS-sorted before sequencing, while in patients 5–7 and patient 11, total BAL cells (containing Vα2.3+Vβ22+ T-cell expansions) were used. In patients 8–11, the mRNA concentration was too low to enable PCR amplification of the Vβ22 chain (table 3). In addition, in patients 2, 3, 12 and 13, FACS-sorted Vα2.3+Vβ22− T-cells and Vα2.3−Vβ22+ were sequenced (table S3).
Patients' T-cell receptor Vα2.3 and Vβ22 amino acid sequences determined by next-generation sequencing (NGS)
In all 13 cases, 1–4 dominant sequences were identified for the highly variable CDR3 region in any given patient, and in three patients only one dominant α chain sequence was observed (figure 3, table 3 and table S3). In six patients, identical or near-identical α chain CDR3 region sequences were identified, all expressing the same Jα gene segment 49*01 (table 3A and table S3). Specifically, the identical CDR3 amino acid sequence CVVNLAGNQFYF was observed in patients 5 and 7, and near-identical CDR3 amino acid sequences CVVNIAGNQFYF, CVVNRATGNQFYF, CVVNMAGNQFYF and CVVNPGTGNQFYF in patients 8, 10, 12 and 3, respectively (in the latter two cases from Vα2.3+Vβ22−-sorted BAL T-cells; table S3). Of note, the leucine residue (L) at CDR3 position 108 is conservatively replaced by an isoleucine (I) in patient 8 and by a methionine (M) in patient 12. Furthermore, the leucine residue at position 108 of the identical CDR3 amino acid sequences shared between patients 5 and 7 was encoded by two different nucleotide sequences (“CTA” in patient 5 and “CTT” in patient 7, respectively) (figure 3a and table S4A). Likewise, in patients 1, 2, 3, 4 and 6, the lysine residue at position 115 was coded for by AAG or AAA.
T-cell receptor (TCR) CDR3 region sequences and amino acid delineation of seven DRB1*03+ patients show inter-patient similarities in amino acid positioning and properties. TCR CDR3 region alignment was performed for the TCR α and β chains of the patients from whom sequencing data could be obtained for both chains. Patients 1–7 correspond to patients 1–7 in table 3. a) Vα2.3 CDR3 sequences and b) Vβ22 CDR3 sequences from the same patients are shown. Figures were generated with tools made available by IMGT (International ImMunoGeneTics Information System), and amino acid positions and chemical properties demarcated according to the IMGT definition.
For the β chains, we found 2–4 dominating clonal sequences in each individual patient, frequently with a preferred expression of Jβ gene segment 2–1*01 (found in four out of seven patient samples) (figure 3b and table 3). A positively charged arginine was found significantly more often at position 112 in the β chains of several patients (four out of 18 sequences, 22%) compared with published reference β chain sequences retrieved from IMGT/LIGM-DB [23] (one (1%) out of 73; p<0.001). In line with this finding, positively charged arginine or lysine amino acid residues at positions 112 or 113 appeared more frequently in our samples (eight (44%) out of 18 sequences, compared to nine (12%) out of 73 in reference β chain sequences; p<0.001) (figure 3b and figure S3). Finally, a negatively charged glutamic acid tended to appear at position 115 more often than in reference sequences (p=0.09). In addition, several of these identical amino acids were encoded by different nucleotide combinations (arginine at position 112–113 and glutamate at position 115, respectively), indicating an important role in specific interactions between these residues and HLA-DRB1*03+ molecules, as well as with HLA-restricted antigens (figure 3b and table S4B).
A three-dimensional molecular model of the Vα2.3+Vβ22+ TCR implicates specific residues for recognition of HLA-DRB1*03 molecules and a vimentin-derived peptide
No structural information has hitherto been provided for TCRs of the Vα2.3/Vβ22 class. An ensemble of three-dimensional molecular models of ternary TCR/HLA/peptide complexes was created based on the CDR loops derived from obtained sequencing results for patient TCRs. The proposed three-dimensional molecular model for the Vα2.3/Vβ22 TCR in complex with HLA-DRB1*03 (figure 4a) revealed that the glycine and asparagine residues Gly73 and Asn77, unique for HLA-DRB1*03, may play a key role through possible interactions with TCR residues Gln30 and Ser51 localised on the CDR1α and CDR2α loops, respectively (figure 4b). Moreover, the presence of a basic arginine (or lysine) residue at position 112 or 113 in close to half of the sequenced CDR3β loops enables interaction with negatively charged HLA-DRB1*03 residues Asp66 or Glu69 (figure 4c), or with a negatively charged or polar residue localised in the middle of the presented peptide. In addition to these observations, the model also indicated that Vβ22 residue Glu115, found in nearly all patients, may be important for selection of the Vα2.3 chain through specific interactions with residues Ser31 and Tyr50 localised on Vα2.3 CDR1α and CDR2α loops, respectively (figure 4b).
Conserved T-cell receptor (TCR) CDR3β residues interact with specific HLA and peptide residues, and preferentially select for Vα2.3 chains. a) Overall ribbon representation of a representative three-dimensional structure of a ternary Vα2.3/Vβ22 TCR/DRB1*03/peptide complex. TCR α and β chains are in yellow and pink, respectively. CDR1α, CDR2α, CDR3α, CDR1β, CDR2β and CDR3β regions are indicated as 1α, 2α, 3α, 1β, 2β and 3β, respectively. HLA α and β chains are in green and light blue, respectively. The peptide bound inside the HLA binding groove is in violet. b) HLA-DRB1*03 residues Gly73 and Asn77 could be responsible for specific selection of Vα2.3 TCR α chains. The HLA residue Gly73 allows for the CDR1α and CDR2α loops to come closer to the HLA heavy chain, possibly resulting in interactions between TCR residues Gln30 and Ser51 and HLA residues Asn77 and Asp76, respectively. Furthermore, a molecular model of the Vβ22 TCR CDR3β loop, found in patient 4, sequence 1 (CASSGGHGKGEQFF), indicates that TCR residue Glu115 (present in nine out of 18 sequences) is localised too far away from the peptide. Instead, this residue can form hydrogen bonds with residues Ser31 and Tyr50 localised on CDR1α and CDR2α, respectively, possibly selecting for a class of TCR α chain such as Vα2.3. c) The conformation of CDR3β loops depends on their length and amino acid compositions. However, the positively charged arginine residue Arg, found in a majority of Vβ22 TCR chains, at either position 112 or 113, is always localised at the top of the loop, independently of the CDR conformations. Residue Arg112/113 can thus interact either with conserved HLA class II β chain residues such as Glu69 or Asp66, or with a negatively charged residue within the middle part of the presented peptide. d) The three-dimensional model indicates that the vimentin-derived peptide Vim429–443 DSLPLVDTHSKRTLL can bind efficiently to HLA-DRB1*03 by making use of residues Leu433, Thr436, Ser438 and Thr441 as anchor positions, occupying pockets 1, 4, 6 and 9, respectively. Furthermore, a molecular model of the Vβ22 TCR CDR3β loop, found in patient 1, sequence 1 (CASSEQGRGETQYF), reveals that the TCR residue Glu114 (present in five out of 18 sequences) can form a salt bridge with a positively charged residue, such as the arginine localised at position 8 in Vim429–443 (Arg440). The same TCR model also indicates that the TCR residue Arg112 may interact with a negatively charged peptide residue such as Asp435 in Vim429–443.
Based on these findings and our previous discoveries regarding the vimentin-derived peptide Vim429–443 (DSLPLVDTHSKRTLL) [24], we hypothesised that vimentin could bind to HLA-DRB1*03 by using residues Leu433, Thr436, Ser438 and Thr441 as anchor positions and consequently occupying pockets 1, 4, 6 and 9 respectively. Other features, including peptide residues Asp435, His437 and Arg440 that can interact with TCR CDR3β residues Arg112 and Glu114, were also observed (figure 4d).
Discussion
Previous studies on TCR Vα2.3+ T-cells in sarcoidosis have not considered paired Vβ gene expression. In early attempts, PCR techniques were applied to address this question, but without any distinct findings [14]. Using new techniques, we initially found a higher-than-normal usage of Vβ22 in lung-accumulated T-cells from HLA-DRB1*03+ patients [16], and in the present report further establish that a substantial part of lung-accumulated Vα2.3+ T-cells also express Vβ22 in HLA-DRB1*03+ patients. Enlarged populations of T-cells expressing both Vα2.3 and Vβ22 were identified in HLA-DRB1*03+ patients, but not in HLA-DRB1*03− patients (with the exception of moderate increases in HLA-DRB3*01+ patients). Analyses of cell surface markers CD69 and CD27 suggested significantly higher levels of Vα2.3+Vβ22+ T-cells to be activated, and also to be differentiated to a higher degree compared to their Vα2.3−Vβ22− CD4+ T-cells counterparts, likely reflecting antigenic stimulation. Alternatively, the high CD69 expression reflects a tissue-resident memory T-cell phenotype, suggested to be of importance for human tissue-specific immune and inflammatory diseases [25].
Finally, molecular modelling pinpointed unique features of the HLA-DRB1*03 molecule and its interaction with specific residues of the Vα2.3+Vβ22+ TCR. This model also provided molecular insight into the preferential pairing of Vα2.3 with Vβ22, as well as with HLA-DRB1*03, and the ideal properties of this complex for accommodation of a potential antigenic peptide derived from cytoskeletal protein vimentin.
Vα2.3+Vβ22+ T-cells may have an enhanced capacity to recognise HLA-DRB1*03-epitope complexes, explaining their accumulation in the lungs of HLA-DRB1*03+ patients. In a recent study, we found that Vα2.3+ T-cells accumulate in BAL fluid but remain at normal levels in regional lymph nodes and blood, suggesting inhalation of a presumed antigen [26]. Given that 21% (median) of lung-accumulated Vα2.3+ T-cells in HLA-DRB1*03+ patients also express Vβ22, our previous studies focusing on Vα2.3+ T-cells alone might have overlooked differences in specificity and functional mechanisms, such as cytokine production and activation marker expression, between different Vα2.3+ T-cell populations. Considering the co-expression of Vβ22, antigen specificity and effector cell functionality can be studied in further detail, taking both the Vα and Vβ gene segments within these lung-accumulated T-cells into account.
The identification of a few dominant TCR Vα2.3 and Vβ22 sequences in each patient demonstrates a high degree of Vα2.3+Vβ22+ T-cell clonality in the lungs, strongly implicating specific antigen recognition in pulmonary sarcoidosis. These clones could be useful to further investigate candidate antigens, such as ESAT-6 [27] or mKatG [28], as well as to unravel molecular details underlying how such TCRs interact with HLA-DRB1*03-peptide complexes. The potential existence of specific antigens restricted to HLA-DRB1*03 is further supported by our findings of identical and near-identical TCR α chain sequences in different patients as well as of similar amino acids at distinct positions within the TCR β chains. Also, the fact that different nucleotide sequences code for the same amino acid residues indicates preferential selection at the protein level, as in the case of antigen recognition. Altogether, these findings strongly support our hypothesis that Vα2.3+Vβ22+ lung T-cells recognise and proliferate in response to sarcoidosis-associated and DRB1*0301-restricted antigens.
The inclusion of a three-dimensional molecular model for elucidating the possible arrangement of specific residues within the different CDR loops allowed us to establish common features of recognition in the Vα2.3/Vβ22 TCR/HLA-DRB1*03 molecular complex. It is well known that sarcoidosis, and especially variants with good prognosis such as Löfgren's syndrome, strongly associate with HLA-DRB1*0301 and HLA-DRB3*0101 alleles [29, 30]. Whether such a good prognosis depends on the HLA molecules themselves, linkage to other gene variants, or their connection to lung-accumulated Vα2.3+Vβ22+ T-cells that may be particularly efficient in eliminating antigens [31], remains to be elucidated. It is worth noting, however, that the more Vα2.3+ cells in BAL fluid, the better the prognosis [32].
Although the majority of residues that differ between HLA-DRB1*0301 and other HLA molecules are found within the peptide-binding cleft, alignment reveals that two specific residues, Gly73 and Asn77, localised on the HLA-molecule far from the presented epitope, are readily available for interaction with TCR residues. Interestingly, these positions are identical in HLA-DRB3*0101, which is also known to correlate with assembly of Vα2.3+ T-cells in the lungs [11], and to some degree with Vα2.3+Vβ22+ T-cells. We therefore hypothesise that Gly73 and Asn77 residues, in combination with specific, not yet identified, epitopes, could be responsible for attracting TCR Vα2.3+ T-cells to the lungs, since residues prominent among all the analysed CDR1α and CDR2α loops are closely positioned to both Gly73 and Asn77, as illustrated in figure 4a.
In addition, the repeated identical or near-identical CDR3 regions of Vα2.3-expressing T-cells can interact with the HLA-DRB1*03 molecule (figure 4). The charged amino acids found at distinct positions in a majority of Vβ22 CDR3 regions enable direct interaction with a peptide antigen harbouring charged amino acids at position 5 and 7, respectively. This suggests that the Vα2.3 gene segment allows for close association with HLA-DRB1*03 (or HLA-DRB3*01) molecules, while Vβ22 CDR3 regions directly recognise the actual antigenic peptide. Accordingly, the Vα2.3 gene segment may pair with other, still unidentified, Vβ segments to recognise other antigenic peptides in complex with HLA-DRB1*03 molecules; this may also explain the higher inter-patient similarity of Vα2.3 compared with Vβ22 CDR3 regions. In addition, since the charged Vβ22 amino acids at position 115 are located far away from the peptide, they may alternatively interact directly with a correspondingly charged amino acid within Vα2.3 CDR1 or CDR2 regions, and thus to some extent explain the preferential pairing of Vα2.3 with Vβ22.
Based on the strong TCR-HLA association presented here, we assessed the binding capacity of potential peptide antigens to the TCR-HLA complex. Interestingly, a previously identified vimentin-derived peptide, Vim429–443 (DSLPLVDTHSKRTLL), eluted from HLA-DR molecules of pooled BAL cells from HLA-DRB1*03+ sarcoidosis patients [24], matches perfectly into both the HLA peptide-binding pocket and the TCR Vβ22 CDR3 loop in the proposed molecular model. The same vimentin peptide has been found to stimulate peripheral blood T-cells of HLA-DRB1*03+ sarcoidosis patients [33].
The intracytoplasmic protein vimentin has been implicated as an auto-antigen in other diseases. Renal inflammation is a frequent and severe manifestation of systemic lupus erythematosus, with tubulointerstitial inflammation (TII) as a common variant. Vimentin was recently shown to be the main targeted auto-antigen in TII [34]. An altered and increased expression of vimentin locally in the kidneys, as well as an association between high titres of serum anti-vimentin autoantibodies and disease severity, was observed in patients with TII, implicating both systemic and local autoimmunity. As vimentin is a basic protein and can be secreted by activated macrophages, it may be abundant in inflammatory conditions [35]. Vimentin can bind to receptors (dectin 1) expressed on dendritic cells, macrophages and B cells, thereby loading and activating antigen-presenting cells [36]. Also, vimentin can be post-translationally modified and thereby become an auto-antigenic target, as shown in rheumatoid arthritis [37]. Whether sarcoidosis may be an autoimmune disease or not, or have an autoimmune component, is still not possible to conclude. While both associations with certain HLA types and reactivity against the self-protein vimentin are in line with the concept of autoimmunity, a seasonal distribution of disease onset [30], as well as the fact that sarcoidosis (especially Löfgren's syndrome) normally runs a self-limiting disease course, would argue against. Moreover, there are examples of transient autoimmunity following infection [38, 39], consistent with our previous finding that strong anti-vimentin T-cell responses were absent in patients after recovery [33].
The theoretical nature of the molecular model is a limitation to this study, and future functional analyses will aim to investigate any antigenic potential of vimentin. To this end, we plan to study B as well as T-cell responses against different forms of vimentin (recombinant protein, individual peptides, post-translationally modified peptides etc.). T-cell cytokine production and proliferation will be analysed using ELISpot and flow cytometry, and we will assess presence of any vimentin-specific Ig subclasses in both BAL fluid and serum. All patients will be HLA-typed, and to identify any associations of T and B cell responses with clinical features, a careful clinical characterisation, including monitoring of disease course and separate examination of specific subgroups e.g. Löfgren's syndrome patients, will have to be performed. However, we need also to consider the possibility that the ideal fit of vimentin into the peptide-binding cleft of HLA-DRB1*03 molecules reflects molecular mimicry, and that a “true” sarcoidosis antigen is in reality a foreign agent such as a microorganism.
In conclusion, the novel identification of lung-accumulated Vα2.3+Vβ22+ T-cell clones in strict association with HLA-DRB1*03 molecules, and inter-patient similarities between both TCR α and β chain sequences, strongly indicate recognition of a specific antigen in these patients. Furthermore, modelling of the TCR-HLA complex suggests a vimentin-derived peptide as a possible autoantigen in sarcoidosis given its potential for binding to the HLA-DRB1*03 molecule and interaction with TCR CDR3β residues. In our continuing search for a sarcoidosis-specific antigen, vimentin will be of particular interest for future studies.
Acknowledgements
The authors thank research nurses Heléne Blomqvist, Margitha Dahl and Gunnel de Forest as well as biomedical analyst Benita Dahlberg for skilful assistance with sample collections and preparations. Also, we would like to acknowledge support from Science for Life Laboratory, the National Genomics Infrastructure (NGI) and Uppmax for providing assistance in massively parallel sequencing and computational infrastructure.
Footnotes
Editorial comment in: Eur Respir J 2016; 47: 707–709 [DOI: 10.1183/13993003.01791-2015].
This article has supplementary material available from erj.ersjournals.com
Support statement: Financial support in the form of research grants was received from The Swedish Heart Lung Foundation, The Swedish Research Council, The Mats Kleberg Foundation, The King Oscar II Jubilee Foundation, The Stockholm County Council, The Swedish Association for Chest Physicians, The Center for Inflammatory Diseases, and Karolinska Institutet. Funding information for this article has been deposited with FundRef.
Conflict of interest: None declared.
- Received July 24, 2015.
- Accepted September 30, 2015.
- Copyright ©ERS 2016