Kornelis S M van der Geest1, Wayel H Abdulahad1, Gerda Horst1, Pedro G Lorencetti1, Johan Bijzet1, Suzanne Arends1, Marieke van der Heiden2, Anne-Marie Buisman3, Bart-Jan Kroesen4, Elisabeth Brouwer1, Annemieke M H Boots1. 1. Department of Rheumatology and Clinical Immunology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands. 2. Department of Rheumatology and Clinical Immunology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands; Center for Immunology of Infectious Diseases and Vaccines, National Institute for Public Health and the Environment, Bilthoven, The Netherlands. 3. Center for Immunology of Infectious Diseases and Vaccines, National Institute for Public Health and the Environment, Bilthoven, The Netherlands. 4. Department Laboratory Medicine, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands.
Abstract
Measuring changes of the T cell receptor (TCR) repertoire is important to many fields of medicine. Flow cytometry is a popular technique to study the TCR repertoire, as it quickly provides insight into the TCR-Vβ usage among well-defined populations of T cells. However, the interpretation of the flow cytometric data remains difficult, and subtle TCR repertoire changes may go undetected. Here, we introduce a novel means for analyzing the flow cytometric data on TCR-Vβ usage. By applying economic statistics, we calculated the Gini-TCR skewing index from the flow cytometric TCR-Vβ analysis. The Gini-TCR skewing index, which is a direct measure of TCR-Vβ distribution among T cells, allowed us to track subtle changes of the TCR repertoire among distinct populations of T cells. Application of the Gini-TCR skewing index to the flow cytometric TCR-Vβ analysis will greatly help to gain better understanding of the TCR repertoire in health and disease.
Measuring changes of the T cell receptor (TCR) repertoire is important to many fields of medicine. Flow cytometry is a popular technique to study the TCR repertoire, as it quickly provides insight into the TCR-Vβ usage among well-defined populations of T cells. However, the interpretation of the flow cytometric data remains difficult, and subtle TCR repertoire changes may go undetected. Here, we introduce a novel means for analyzing the flow cytometric data on TCR-Vβ usage. By applying economic statistics, we calculated the Gini-TCR skewing index from the flow cytometric TCR-Vβ analysis. The Gini-TCR skewing index, which is a direct measure of TCR-Vβ distribution among T cells, allowed us to track subtle changes of the TCR repertoire among distinct populations of T cells. Application of the Gini-TCR skewing index to the flow cytometric TCR-Vβ analysis will greatly help to gain better understanding of the TCR repertoire in health and disease.
A broad T cell receptor (TCR) repertoire is considered critical for optimal immunity against a wide variety of antigens [1,2]. Contraction of the TCR repertoire is associated with poor vaccine responses in aged individuals [3], progression to AIDS in HIV-infectedpatients [4], and poor survival in cancerpatients [5]. Assessment of the T cell receptor (TCR) repertoire is therefore relevant to researchers in many fields of medicine.Several techniques are currently used to study the TCR repertoire [6]. Moderate to high resolution assessment of the TCR repertoire is provided by PCR-based methods, such as TCR spectratyping and sequencing. As these methods are relatively labor-intensive and preferably require cell-sorting of highly pure T cell populations, many researchers turn to flow cytometry [6]. Flow cytometry quickly measures the proportional TCR-Vβ usage in multiple T cell subsets on a per-cell basis, without the need for cell-sorting. Although this assessment can provide helpful information, an accurate and reliable way for analyzing the flow cytometric data on the TCR repertoire is currently lacking.Here, we introduce economic statistics to improve the analysis of flow cytometric data on TCR-Vβ usage. We noticed that the distribution of TCR-Vβ families among T cells resembles the distribution of income among people (Fig 1A–1C). Economists typically study the distribution of income by constructing Lorenz curves and calculating the Gini index. The Gini index, with scores ranging from 0 to 100, is a direct measure of income distribution [7-9]. By applying the Gini index to the flow cytometric TCR-Vβ analysis, we could directly measure the distribution of 24 TCR-Vβ families among multiple, well-defined T cell subsets. In this context, low Gini index values indicated equal distribution of TCR-Vβ families (i.e. broad repertoire), whereas high values pointed to unequal distribution of TCR-Vβ families (i.e. repertoire skewing). Although the Gini index has been used in TCR sequencing studies [10,11], we here demonstrate that the Gini index, hence referred to as the Gini-TCR skewing index, also aids the analysis of flow cytometric data on TCR-Vβ usage. Importantly, the Gini-TCR skewing index allowed us to detect subtle changes of the TCR repertoire among multiple, well-defined T cell subpopulations.
Fig 1
Schematic overview showing the relation between T cell receptor (TCR) Vβ diversity, distribution and percentages.
(A) Schematic drawing illustrating broad versus contracted TCR-Vβ diversity. (B) Schematic drawing illustrating distribution and proportional usage of TCR Vβ families when TCR-Vβ diversity is broad or contracted. (C) Schematic drawing illustrating distribution of income among people.
Schematic overview showing the relation between T cell receptor (TCR) Vβ diversity, distribution and percentages.
(A) Schematic drawing illustrating broad versus contracted TCR-Vβ diversity. (B) Schematic drawing illustrating distribution and proportional usage of TCR Vβ families when TCR-Vβ diversity is broad or contracted. (C) Schematic drawing illustrating distribution of income among people.
Methods
Subjects
Heparinized blood samples were acquired from 27 healthy volunteers. Eight men and 19 women were included (age range 22–81). Samples of 5 children (age 9 years) undergoing DTaP-IPV vaccination were also collected. Written informed consent was obtained from adult volunteers or from parents on behalf of the children. The study was approved by the institutional review board of the UMCG (METc2012/375) and the Central Committee on Research Involving Human Subjects in the Netherlands (CCMO; ISRCTN64117538). None of the study participants had an overt history of infection, cancer or auto-immune disease.
Flow cytometry
Whole blood samples (150 μL) or isolated peripheral blood mononuclear cells were incubated with fluorochrome-conjugated antibodies (Table 1 and S1 Fig) for 45 minutes at room temperature. Subsequently, whole blood samples were lysed with 2 mL of 1x Lysing solution (BD Biosciences) for 10 minutes at room temperature. Finally, samples were washed twice with phosphate buffered saline containing 1% bovine serum albumin. Samples were measured immediately on a LSR-II flow cytometer (BD Biosciences). The flow cytometric data was obtained with FACS Diva (BD Biosciences) and analyzed with Kaluza Software (Beckman Coulter). Fluorescence-minus-one (FMO) controls served as negative controls for the flow cytometric staining.
Table 1
Monoclonal antibodies.
Antibody specificity
Fluorochrome
Clone
Manufacturer
Vβ1
PE+FITC
BL37.2
Beckman Coulter
Vβ2
PE+FITC
MPB2D5
Beckman Coulter
Vβ3
FITC
CH92
Beckman Coulter
Vβ4
PE+FITC
WJF24
Beckman Coulter
Vβ5.1
PE+FITC
IMMU157
Beckman Coulter
Vβ5.2
PE
36213
Beckman Coulter
Vβ5.3
PE
3D11
Beckman Coulter
Vβ7.1
PE+FITC
ZOE
Beckman Coulter
Vβ7.2
FITC
ZIZOU4
Beckman Coulter
Vβ8
FITC
56C5.2
Beckman Coulter
Vβ9
PE
FIN9
Beckman Coulter
Vβ11
PE
C21
Beckman Coulter
Vβ12
FITC
VER2.32
Beckman Coulter
Vβ13.1
PE
IMMU222
Beckman Coulter
Vβ13.2
PE
H132
Beckman Coulter
Vβ13.6
PE+FITC
JU74.3
Beckman Coulter
Vβ14
FITC
CAS1.1.3
Beckman Coulter
Vβ16
FITC
TAMAYA1.2
Beckman Coulter
Vβ17
PE+FITC
E17.5F3
Beckman Coulter
Vβ18
PE
BA62.6
Beckman Coulter
Vβ20
FITC
ELL1.4
Beckman Coulter
Vβ21.3
FITC
IG125
Beckman Coulter
Vβ22
PE+FITC
IMMU546
Beckman Coulter
Vβ23
PE
AF23
Beckman Coulter
CD28
APC
CD28.2
Beckman Coulter
CD31
Pacific blue
5.6E
Beckman Coulter
CD4
APC-AF700
13B8.2
Beckman Coulter
CD45RA
APC-AF750
2H4
Beckman Coulter
CD8
Krome orange
B9.11
Beckman Coulter
TCRγδ
PC5.5
IMMU510
Beckman Coulter
CCR7
PE-Cy7
3D12
BD Biosciences
Overview of fluorochrome-conjugated monoclonal antibodies applied in the study.
Overview of fluorochrome-conjugated monoclonal antibodies applied in the study.
Calculation of Gini-TCR skewing index
Gini-TCR skewing index values were calculated for T cell subsets of individual blood donors with use of the Gini index, which is generally used to measure income distribution [7,8]. A Microsoft Excel file allowing automatic calculation of Gini-TCR skewing index from percentages of 24 TCR-Vβ families is provided in the Supporting Information (S1 File). A simplified overview to calculate the Gini-TCR skewing index is provided (S2 Fig). Briefly, proportions of all 24 TCR-Vβ families within a T cell subset were arranged from small to large. To construct Lorentz curves, the cumulative proportions of all TCR-Vβ family proportions were normalized to a total of 100%. Lorenz curves for individual blood donors were plotted with the cumulative percentage of the 24 TCR-Vβ families studied on the x-axis, and the cumulative proportion of CD4 or CD8 T cells that were covered by these 24 TCR-Vβ families on the y-axis. Next, the surface under the Lorenz curve was subtracted from the surface under the equidistribution line (i.e. line of perfectly equal distribution). To calculate the Gini-TCR skewing index, the surface between the equidistribution line and the Lorenz curve was divided by the total surface under the equidistribution line (100 * 100 / 0.5 = 5000) and then multiplied by 100. An example of the complete calculation of a Gini-TCR skewing index is shown in Table 2.
Table 2
Example showing calculation of Gini-TCR skewing index from TCR-Vβ usage of an individual blood donor.
#
Vβ usage (%)
Ranked Vβ usage (%)
Ranked Vβ usage (%)—Corrected to 100%
Cumulative TCR Vβ usage (%)
Surface under Lorenz curve
1
Vβ1
3.23
Vβ23
0.58
0.78
0.78
A
3.25
2
Vβ2
10.84
Vβ11
0.77
1.03
1.81
B
5.40
3
Vβ3
7.62
Vβ7.2
0.78
1.05
2.86
C
9.74
4
Vβ4
2.68
Vβ5.2
0.92
1.24
4.10
D
14.49
5
Vβ5.1
5.65
Vβ5.3
0.92
1.24
5.33
E
19.64
6
Vβ5.2
0.92
Vβ16
1.18
1.58
6.92
F
25.52
7
Vβ5.3
0.92
Vβ20
1.48
1.99
8.90
G
32.96
8
Vβ7.1
1.97
Vβ13.6
1.56
2.10
11.00
H
41.47
9
Vβ7.2
0.78
Vβ21.3
1.88
2.52
13.52
I
51.09
10
Vβ8
3.84
Vβ7.1
1.97
2.65
16.17
J
61.86
11
Vβ9
4.72
Vβ13.2
2.13
2.86
19.03
K
73.33
12
Vβ11
0.77
Vβ12
2.22
2.98
22.01
L
85.50
13
Vβ12
2.22
Vβ4
2.68
3.60
25.61
M
99.21
14
Vβ13.1
3.34
Vβ1
3.23
4.34
29.95
N
115.75
15
Vβ13.2
2.13
Vβ18
3.33
4.47
34.42
O
134.10
16
Vβ13.6
1.56
Vβ13.1
3.34
4.49
38.91
P
152.77
17
Vβ14
3.41
Vβ14
3.41
4.58
43.49
Q
171.65
18
Vβ16
1.18
Vβ22
3.69
4.96
48.44
R
191.52
19
Vβ17
5.72
Vβ8
3.84
5.16
53.60
S
212.59
20
Vβ18
3.33
Vβ9
4.72
6.34
59.94
T
236.54
21
Vβ20
1.48
Vβ5.1
5.65
7.59
67.53
U
265.55
22
Vβ21.3
1.88
Vβ17
5.72
7.68
75.21
V
297.36
23
Vβ22
3.69
Vβ3
7.62
10.23
85.44
W
334.69
24
Vβ23
0.58
Vβ2
10.84
14.56
100.00
X
386.34
Total
74.46
74.46
100.00
3022.32
First TCR-Vβ families are ranked based on their proportional size. Then values are normalized to 100%. Subsequently, the surface under the Lorentz curve is calculated. In this blood donor, the surface between Lorenz curve and equidistribution line = 5000–3022.32 = 1977.68. Gini-TCR skewing index = 1977.68/5000*100 = 39.55.
First TCR-Vβ families are ranked based on their proportional size. Then values are normalized to 100%. Subsequently, the surface under the Lorentz curve is calculated. In this blood donor, the surface between Lorenz curve and equidistribution line = 5000–3022.32 = 1977.68. Gini-TCR skewing index = 1977.68/5000*100 = 39.55.
Hierarchical clustering analysis
Gini-TCR skewing index values of all T cell subsets from individual donors were plotted in a heat map using Microsoft Excel 2007. Colors were appointed based on absolute Gini-TCR skewing index values. Subsequently, hierarchical clustering analysis of T cell subsets and blood donors was performed in Genesis 1.7.6 (Genomics and Bioinformatics Graz) using Euclidean distance and compete linkage clustering.
CMV ELISA
CMV-specific IgG was determined in serum samples using an in-house ELISA. 96-well ELISA plates (Greiner) were coated with lysates of CMV-infected fibroblasts overnight. Lysates of non-infected fibroblasts were used as negative controls. Following the coating, dilutions of serum samples were incubated for 1 hour. Goat anti-human IgG was added and incubated for 1 hour. Samples were incubated with phosphatase for 15 minutes and the reaction was stopped with NaOH. The plates were scanned on a Versamax reader (Molecular Devices). A pool of sera from 3 CMV-seropositive individuals with known concentrations of CMV-specific IgG was used to quantify levels of CMV-specific IgG in the tested samples.
Vaccination
Nine year old children in good health were vaccinated with a combination vaccine against diphtheria, tetanus, acellular pertussis and inactivated poliovirus (DTaP-IPV). Peripheral blood mononuclear cells were isolated by density centrifugation before vaccination, as well as 28 days and 1 year after vaccination. Samples were viably stored at -135°C with RPMI + 10% FCS and 10% DMSO. To perform flow cytometric staining, the samples were thawed in a 37°C water bath and washed twice with RPMI + 10% FCS.
Statistics
Gini-TCR skewing index values between two T cell subsets within individuals were compared with the Wilcoxon signed rank test. In case multiple T cell subsets were compared, the Wilcoxon signed rank test was preceded by the Friedman test. Gini-TCR skewing index values in CMV seropositive and seronegative individuals were compared with the Mann Whitney U test. Two-sided p values < 0.05 were considered statistically significant. Statistical analysis was performed with IBM SPSS Statistics 20 (SPSS, Chicago, IL, USA).
Results
Calculation of the Gini-TCR skewing index
To demonstrate that flow cytometry and economic statistics can be combined, we started by calculating Gini-TCR skewing index values for CD4 and CD8 T cells from healthy individuals. Therefore, we first measured the proportional usage of 24 TCR-Vβ families in CD4 and CD8 T cells with a commercially available panel of monoclonal antibodies (S1 Fig). In accordance with previous reports, the antibody panel recognized TCR-Vβ families on 66 ± 1% (mean ± SEM) of CD4 T cells and 60 ± 2% of CD8 T cells [12,13]. At this stage, it was difficult to conclude whether the distribution of TCR-Vβ families differed between CD4 and CD8 T cells (Fig 2A). Next, the proportional TCR-Vβ usage of CD4 and CD8 T cells from individual blood donors was plotted in Lorenz curves with the cumulative percentage of the 24 TCR-Vβ families studied on the x-axis, and the cumulative proportion of CD4 or CD8 T cells that was covered by these 24 TCR-Vβ families on the y-axis (Fig 2B and S2 Fig). A Lorenz curve far from the equidistribution line refers to unequal distribution and skewing of the TCR repertoire. Finally, Gini-TCR skewing index values were calculated for CD4 and CD8 T cells of individual blood donors (S2 Fig). In some individuals, Gini-TCR skewing index values were markedly higher for CD8 T cells than CD4 T cells, whereas the values were similar in others (Fig 2C). This difference could be explained by cytomegalovirus (CMV) infection, as only CD8 T cells of CMV-seropositive individuals demonstrated increased Gini-TCR skewing index values (Fig 2D). Furthermore, we noted that variation of Gini-TCR skewing index values in CMV-seronegative could be explained by age, as we observed a substantial, aging-associated increase in Gini-TCR skewing index values among CD8 T cells of CMV-seronegative individuals (Fig 2E). In contrast, the TCR-Vβ repertoire of CD8 T cells was already skewed in young CMV-seropositive donors and did not modulate further with age. In addition, Gini-TCR skewing index values only modestly increased with age among CD4 T cells, irrespective of CMV-serostatus. With these experiments, we demonstrated that flow cytometry can be combined with economic statistics and confirmed previous findings that CMV-infection and aging are associated with profound skewing of the TCR repertoire in CD8 T cells [14-16].
Fig 2
T cell receptor (TCR) Vβ distribution analysis in CD4 and CD8 T cells.
(A) Traditional graph indicating proportional TCR-Vβ usage in CD4 and CD8 T cells of 27 healthy individuals. Bars and error bars represent median and interquartile ranges. (B) Lorenz curve of TCR-Vβ distribution in CD4 and CD8 T cells of a representative healthy individual. (C) Gini-TCR skewing index values in CD4 and CD8 T cells of 27 healthy individuals. Lines represent medians. (D) Gini-TCR skewing index values in CD4 and CD8 T cells of CMV-seropositive (CMV+) and seronegative (CMV-) individuals. Lines represent medians. (E) Relation between age and Gini-TCR skewing index values in CD4 and CD8 T cells from 27 healthy CMV+ and CMV- individuals. Statistical significance was tested with the Wilcoxon signed rank test or Mann Whitney U test and is indicated as * p<0.05, ** p<0.01, *** p<0.001.
T cell receptor (TCR) Vβ distribution analysis in CD4 and CD8 T cells.
(A) Traditional graph indicating proportional TCR-Vβ usage in CD4 and CD8 T cells of 27 healthy individuals. Bars and error bars represent median and interquartile ranges. (B) Lorenz curve of TCR-Vβ distribution in CD4 and CD8 T cells of a representative healthy individual. (C) Gini-TCR skewing index values in CD4 and CD8 T cells of 27 healthy individuals. Lines represent medians. (D) Gini-TCR skewing index values in CD4 and CD8 T cells of CMV-seropositive (CMV+) and seronegative (CMV-) individuals. Lines represent medians. (E) Relation between age and Gini-TCR skewing index values in CD4 and CD8 T cells from 27 healthy CMV+ and CMV- individuals. Statistical significance was tested with the Wilcoxon signed rank test or Mann Whitney U test and is indicated as * p<0.05, ** p<0.01, *** p<0.001.
Gini-TCR skewing index in T cell differentiation subsets
To further validate our method, we determined Gini-TCR skewing index values in well-defined T cell differentiation subsets. T cells were therefore divided into naive (TNaive), central memory (TCM), effector memory (TEM) and terminally differentiated (TTD) T cells based on CD45RA and CCR7 expression [17] as shown in S3 Fig TCR sequencing and spectratyping studies have shown that naive T cells demonstrate the broadest TCR repertoire, whereas the repertoire becomes more contracted in memory T cells [18-20]. We indeed observed that CD4 TNaive cells show significantly lower Gini-TCR skewing index values than the three CD4 memory cells subsets (Fig 3A). Among the memory cells, CD4 TCM cells demonstrated significantly less TCR skewing than CD4 TEM and TTD cells, as previously found with TCR sequencing [18]. Further delineation of CD4 TNaive cells based on CD31 expression (S3 Fig) confirmed the previous observation that CD31 negative T cells (central naive) have a slightly more skewed TCR repertoire than CD31 positive T cells (thymic naive) [21], as demonstrated in Fig 3B. Analysis of TCR-Vβ distribution in CD8 differentiation subsets (S3 Fig) resulted in similar findings as observed in CD4 T cells, although the TCR repertoire skewing was more pronounced in the CD8 TEM and TTD populations (Fig 3C). Taken together, calculation of the Gini-TCR skewing index in naive and memory T cell subsets provided similar results as previously reported with TCR sequencing and spectratyping [18-21].
Fig 3
T cell receptor (TCR) Vβ distribution in CD4 and CD8 T cell subsets.
(A) Gini-TCR skewing index values in naive (TNaive; CD45RA+CCR7+), central memory (TCM; CD45RA-CCR7+), effector memory (TEM; CD45RA-CCR7-) and terminally differentiated (TTD; CD45RA+CCR7-) CD4 T cells of 27 healthy individuals. (B) Gini-TCR skewing index values in CD31+ and CD31- naive CD4 T cells of 27 healthy individuals. (C) Gini-TCR skewing index values in CD8 TNaive, TCM, TEM and TTD cells of 27 healthy individuals. Lines represent medians. Statistical significance was tested with the Wilcoxon signed rank test and is indicated as * p<0.05, ** p<0.01, *** p<0.001. In case of multiple testing, the Wilcoxon signed rank test was preceded by the Friedman test. (D) Heat map based on absolute Gini-TCR skewing index values in CD4 and CD8 T cell differentiation subsets of healthy donors in which CMV status was determined. Age is shown between ‘()’. The dendrogram was based on hierarchical clustering analysis with use of Euclidean distance and complete linkage clustering.
T cell receptor (TCR) Vβ distribution in CD4 and CD8 T cell subsets.
(A) Gini-TCR skewing index values in naive (TNaive; CD45RA+CCR7+), central memory (TCM; CD45RA-CCR7+), effector memory (TEM; CD45RA-CCR7-) and terminally differentiated (TTD; CD45RA+CCR7-) CD4 T cells of 27 healthy individuals. (B) Gini-TCR skewing index values in CD31+ and CD31- naive CD4 T cells of 27 healthy individuals. (C) Gini-TCR skewing index values in CD8 TNaive, TCM, TEM and TTD cells of 27 healthy individuals. Lines represent medians. Statistical significance was tested with the Wilcoxon signed rank test and is indicated as * p<0.05, ** p<0.01, *** p<0.001. In case of multiple testing, the Wilcoxon signed rank test was preceded by the Friedman test. (D) Heat map based on absolute Gini-TCR skewing index values in CD4 and CD8 T cell differentiation subsets of healthy donors in which CMV status was determined. Age is shown between ‘()’. The dendrogram was based on hierarchical clustering analysis with use of Euclidean distance and complete linkage clustering.
Hierarchical clustering analysis of Gini-TCR skewing index values
As the distribution of the TCR repertoire was simultaneously measured in a large number of distinct CD4 and CD8 T cell subsets, we next determined to what extent TCR skewing in one subset is linked to TCR skewing in another subset. Hierarchical clustering analysis readily revealed that CD8 TEM and CD8 TTD cells formed a separate cluster from all other CD4 and CD8 T cell subsets (Fig 3D). This finding indicates that TCR skewing in CD8 TEM and TTD cells occurs relatively independently from TCR skewing in the other subsets. An additional hierarchical clustering analysis showed that Gini-TCR skewing index values of CMV-seropositive and CMV-seronegative individuals largely formed two separate clusters. This finding further underscores the impact that CMV has on the TCR repertoire in humans [14].
Tracking Gini-TCR skewing index changes following an antigenic challenge
Finally, we determined whether the Gini-TCR skewing index allowed us to detect TCR repertoire changes following an antigenic challenge. Therefore, we studied the TCR repertoire of CD4 and CD8 TNaive and TEM cells in 9-year old children receiving a combination vaccine against diphtheria, tetanus, acellular pertussis and inactivated poliovirus (DTaP-IPV). We hypothesized that DTaP-IPV vaccination will recruit naive T cells towards the effector memory T cell pool and thereby cause skewing of the naive TCR repertoire. In accordance with our expectations, we observed clear skewing of the TCR repertoire among CD8 TNaive cells at 28 days after vaccination (Fig 4A, left panel). In parallel, we observed a slight decrease in Gini-TCR skewing index values, i.e. broadening of the TCR repertoire, among CD8 TEM cells at 28 days (Fig 4A, right panel). The latter finding was only temporary, as Gini-TCR skewing index values of CD8 TEM cells returned to baseline levels at 1 year after vaccination. In contrast, Gini-TCR skewing index values remained elevated among CD8 TNaive cells after 1 year follow-up. These TCR repertoire changes were associated with a stable numerical decrease in CD8 TNaive cells (Fig 4B, left panel) and increase in CD8 TEM cells (Fig 4B, right panel). In essence, we observed comparable modulations within the CD4 T cell compartment (Fig 4C and 4D), albeit substantially less pronounced. These findings indeed implied that naive T cells are recruited towards the effector memory pool upon DTaP-IPV vaccination. Taken together, our results show that the Gini-TCR skewing index provides stable values over time, but is sensitive enough to detect TCR repertoire changes following an antigenic challenge.
Fig 4
Analysis of naive and effector memory CD4 and CD8 T cells before and after DTaP-IPV vaccination.
Five healthy children received a combination vaccine against diphtheria, tetanus, acellular pertussis and inactivated poliovirus (DTaP-IPV). Blood samples were collected before vaccination, and 28 days and 1 year after vaccination. (A) Gini-TCR skewing index values and (B) proportions of naive (TNaive) and effector memory (TEM) CD8 T cells before and after vaccination. (C) Gini-TCR skewing index values and (D) proportions of CD4 TNaive and TEM cells before and after vaccination. Boxes show median with interquartile range, whiskers are drawn according to the Tukey method.
Analysis of naive and effector memory CD4 and CD8 T cells before and after DTaP-IPV vaccination.
Five healthy children received a combination vaccine against diphtheria, tetanus, acellular pertussis and inactivated poliovirus (DTaP-IPV). Blood samples were collected before vaccination, and 28 days and 1 year after vaccination. (A) Gini-TCR skewing index values and (B) proportions of naive (TNaive) and effector memory (TEM) CD8 T cells before and after vaccination. (C) Gini-TCR skewing index values and (D) proportions of CD4 TNaive and TEM cells before and after vaccination. Boxes show median with interquartile range, whiskers are drawn according to the Tukey method.
Discussion
We here report a novel application of the Gini index for analyzing flow cytometric data on the TCR-Vβ repertoire. Flow cytometry is widely used to investigate the humanTCR repertoire, as it quickly provides insight into the TCR repertoire of many different T cell subsets without the need for laborious cell-sorting. Until now however, the analysis and interpretation of the flow cytometric data on TCR-Vβ usage was difficult and subtle TCR repertoire changes were easily missed [6]. In the current study, we show that economic statistics, which have previously been applied in TCR sequencing studies [10,11], allowed us to quantify the flow cytometric data on proportional TCR-Vβ usage into direct measures of TCR-Vβ distribution: the Gini-TCR skewing index. Using the Gini-TCR skewing index, we were capable of detecting subtle changes of the TCR repertoire among multiple, well-defined populations of T cells.We validated the Gini-TCR skewing index by comparing our results to previous findings obtained with TCR sequencing and spectratyping [14,18-21]. We confirmed that chronic CMV infection and aging are associated with substantial skewing of the TCR repertoire in CD8 T cells [14-16]. Moreover, we found that naive and memory T cell populations display differences in TCR repertoire skewing. Naive T cells are known to have the broadest TCR repertoire, whereas the repertoire becomes skewed in memory T cells [18-20]. By calculating the Gini-TCR skewing index, we could confirm that naive T cells have a less skewed TCR-Vβ repertoire than memory T cells. In addition, we verified that CD31+ thymic naive CD4+ T cells have a broader TCR-Vβ repertoire than CD31- central naive CD4+ T cells [21]. Furthermore, we confirmed that central memory CD4 and CD8 T cells demonstrate less TCR repertoire skewing than effector memory and terminally differentiated memory T cells [18]. Interestingly, we could even detect changes of the TCR-Vβ repertoire following an antigenic challenge (i.e. DTaP-IPV vaccination). Our findings therefore confirm that distinct populations of T cells are characterized by distinct skewing of the TCR repertoire, and that these differences can be detected with the Gini-TCR skewing index. Furthermore, our data underscore that future TCR repertoire studies should take into account that age, antigenic challenges and the presence of chronic viral infections (i.e. CMV) have a profound impact on the TCR repertoire.In summary, we introduce a novel and rapid approach for analyzing flow cytometric data on TCR-Vβ usage. By applying economic statistics, we could calculate the Gini-TCR skewing index and readily track subtle changes of the TCR-Vβ repertoire among a large number of T cell subsets. Our study therefore reveals interesting synergy between economic statistics and the flow cytometric analysis of the TCR repertoire. We expect that application of the Gini-TCR skewing index to the flow cytometric TCR-Vβ analysis will contribute to a better understanding of the TCR repertoire in health and disease. The development of novel monoclonal antibodies for additional TCR-Vβ families could make the flow cytometric assessment of the TCR repertoire also more suitable for use in clinical diagnostics [4,5,22].
Flow cytometric analysis of T cell receptor (TCR) Vβ usage.
(A) Representative flow cytometric plot showing gating of CD4 T cells. Additional negative gating for CD8 T cells and γδ T cells was performed. (B) Representative flow cytometric plot showing gating of CD8 T cells. Additional negative gating for CD4 T cells and γδ T cells was performed. (C) Representative flow cytometric plots showing staining for 24 TCR-Vβ families in 8 separate tubes. Anti-TCR-Vβ antibodies were labeled with PE, FITC and PE+FITC.(DOCX)Click here for additional data file.
Simplified calculation of Gini-TCR skewing index based on the Gini index.
A step-by-step explanation shows how to calculate Gini-TCR skewing index values. For clarity purposes, only 5 hypothetical TCR-Vβ families are shown instead of 24 TCR-Vβ families.(DOCX)Click here for additional data file.
Flow cytometric analysis of CD4 and CD8 T cell differentiation subsets.
(A) Representative flow cytometric staining for CD45RA and CCR7 to identify naïve, central memory (CM), effector memory (EM) and terminally differentiated (TD) CD4 and CD8 T cells. The flow cytometric gates were based on fluorescence-minus-one (FMO) controls. (B) Representative flow cytometric staining for CD31 in CD4 T cells and FMO control for CD31.(DOCX)Click here for additional data file.
Microsoft Excel file allowing automatic calculation of the Gini-TCR skewing index.
Percentages of the 24 TCR-Vβ families can be copied into the depicted section (colored cells). Following arrangement of values from low to high, the Gini-TCR skewing index will be calculated (in blue).(XLSX)Click here for additional data file.
Authors: Qian Qi; Yi Liu; Yong Cheng; Jacob Glanville; David Zhang; Ji-Yeun Lee; Richard A Olshen; Cornelia M Weyand; Scott D Boyd; Jörg J Goronzy Journal: Proc Natl Acad Sci U S A Date: 2014-08-25 Impact factor: 11.205
Authors: G Gorochov; A U Neumann; A Kereveur; C Parizot; T Li; C Katlama; M Karmochkine; G Raguin; B Autran; P Debré Journal: Nat Med Date: 1998-02 Impact factor: 53.440
Authors: Paul L Klarenbeek; Paul P Tak; Barbera D C van Schaik; Aeilko H Zwinderman; Marja E Jakobs; Zhuoli Zhang; Antoine H C van Kampen; René A W van Lier; Frank Baas; Niek de Vries Journal: Immunol Lett Date: 2010-07-16 Impact factor: 3.685
Authors: Olga V Britanova; Ekaterina V Putintseva; Mikhail Shugay; Ekaterina M Merzlyak; Maria A Turchaninova; Dmitriy B Staroverov; Dmitriy A Bolotin; Sergey Lukyanov; Ekaterina A Bogdanova; Ilgar Z Mamedov; Yuriy B Lebedev; Dmitriy M Chudakov Journal: J Immunol Date: 2014-02-07 Impact factor: 5.422
Authors: Marita Bosticardo; Francesca Pala; Enrica Calzoni; Ottavia M Delmonte; Kerry Dobbs; Cameron L Gardner; Nicolo' Sacchetti; Tomoki Kawai; Elizabeth K Garabedian; Debbie Draper; Jenna R E Bergerson; Suk See DeRavin; Alexandra F Freeman; Tayfun Güngör; Nicholas Hartog; Steven M Holland; Donald B Kohn; Harry L Malech; Mary Louise Markert; Katja G Weinacht; Anna Villa; Christopher S Seet; Amelie Montel-Hagen; Gay M Crooks; Luigi D Notarangelo Journal: Blood Adv Date: 2020-06-23
Authors: Marieke van der Heiden; Menno C van Zelm; Sophinus J W Bartol; Lia G H de Rond; Guy A M Berbers; Annemieke M H Boots; Anne-Marie Buisman Journal: Sci Rep Date: 2016-05-31 Impact factor: 4.379
Authors: Ling Huang; Michiel G H Betjes; Mariska Klepper; Anton W Langerak; Carla C Baan; Nicolle H R Litjens Journal: Front Immunol Date: 2017-12-15 Impact factor: 7.561
Authors: Kostas Patas; Anne Willing; Cüneyt Demiralay; Jan Broder Engler; Andreea Lupu; Caren Ramien; Tobias Schäfer; Christian Gach; Laura Stumm; Kenneth Chan; Marissa Vignali; Petra C Arck; Manuel A Friese; Ole Pless; Klaus Wiedemann; Agorastos Agorastos; Stefan M Gold Journal: Front Immunol Date: 2018-02-20 Impact factor: 7.561
Authors: Kornelis S M van der Geest; Qi Wang; Thijs M H Eijsvogels; Hans J P Koenen; Irma Joosten; Elisabeth Brouwer; Maria T E Hopman; Joannes F M Jacobs; Annemieke M H Boots Journal: Immun Ageing Date: 2017-02-22 Impact factor: 6.400
Authors: Rosanne D Reitsema; Annemieke M H Boots; Kornelis S M van der Geest; Maria Sandovici; Peter Heeringa; Elisabeth Brouwer Journal: Front Immunol Date: 2021-03-18 Impact factor: 7.561
Authors: Ling Huang; Anton W Langerak; Ingrid L M Wolvers-Tettero; Ruud W J Meijers; Carla C Baan; Nicolle H R Litjens; Michiel G H Betjes Journal: Immun Ageing Date: 2015-12-14 Impact factor: 6.400
Authors: Kornelis S M van der Geest; Bart-Jan Kroesen; Gerda Horst; Wayel H Abdulahad; Elisabeth Brouwer; Annemieke M H Boots Journal: Front Immunol Date: 2018-04-19 Impact factor: 7.561
Authors: Rob R H van den Brom; Kornelis S M van der Geest; Elisabeth Brouwer; Geke A P Hospers; Annemieke M H Boots Journal: Cancer Immunol Immunother Date: 2018-03-15 Impact factor: 6.968