# IMGT/HighV-QUEST Statistical Significance of IMGT Clonotype (AA) Diversity per Gene for Standardized Comparisons of Next Generation Sequencing Immunoprofiles of Immunoglobulins and T Cell Receptors

Safa Aouinti, Dhafer Malouche, Véronique Giudicelli, Sofia Kossida, Marie-Paule Lefranc

► To cite this version:

Safa Aouinti, Dhafer Malouche, Véronique Giudicelli, Sofia Kossida, Marie-Paule Lefranc. IMGT/HighV-QUEST Statistical Significance of IMGT Clonotype (AA) Diversity per Gene for Standardized Comparisons of Next Generation Sequencing Immunoprofiles of Immunoglobulins and T Cell Receptors. PLoS ONE, 2015, 10 (11), pp.e142353. 10.1371/journal.pone.0142353. [hal-01228428](https://hal.science/hal-01228428)

HAL Id: hal-01228428

https://hal.science/hal-01228428

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# IMGT/HighV-QUEST Statistical Significance of IMGT Clonotype (AA) Diversity per Gene for Standardized Comparisons of Next Generation Sequencing Immunoprofiles of Immunoglobulins and T Cell Receptors

Safa Aouinti¹,²,³, Dhafer Malouche²,³, Véronique Giudicelli¹, Sofia Kossida¹*, Marie-Paule Lefranc¹*

1 IMGT, the international ImMunoGeneTics information system, Laboratoire d'ImmunoGénétique Moléculaire LIGM, Institut de Génétique Humaine IGH, UPR CNRS 1142, Montpellier University, Montpellier cedex 5, France, 2 Higher School of Statistics and Information Analysis, University of Carthage, Tunis, Tunisia, 3 National School of Engineers of Tunis, Information Technology and Communications Department, Laboratory U2S, University of Tunis El-Manar, Tunis, Tunisia

* sofia.kossida@igh.cnrs.fr (SK); marie-paule.lefranc@igh.cnrs.fr (M-PL)

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Aouinti S, Malouche D, Giudicelli V, Kossida S, Lefranc M-P (2015) IMGT/HighV-QUEST Statistical Significance of IMGT Clonotype (AA) Diversity per Gene for Standardized Comparisons of Next Generation Sequencing Immunoprofiles of Immunoglobulins and T Cell Receptors. PLoS ONE 10(11): e0142353. doi:10.1371/journal.pone.0142353

**Editor:** Rachel Louise Allen, University of London, St George's, UNITED KINGDOM

**Received:** September 15, 2015

**Accepted:** September 18, 2015

**Published:** November 5, 2015

**Copyright:** © 2015 Aouinti et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Sequencing data is available in the NCBI Sequence Read Archive under the accession code SRX326382.

**Funding:** The authors received no specific funding for this work. Authors Safa Aouinti, Véronique Giudicelli, Sofia Kossida and Marie-Paule Lefranc are affiliated to IMGT, the international ImMunoGeneTics information system (Montpellier University and CNRS). IMGT, the international ImMunoGeneTics information system (Montpellier University and CNRS) provided support in the form of salaries for

## Abstract

The adaptive immune responses of humans and of other jawed vertebrate species (gnasthostomata) are characterized by the B and T cells and their specific antigen receptors, the immunoglobulins (IG) or antibodies and the T cell receptors (TR) (up to 2.10¹² different IG and TR per individual). IMGT, the international ImMunoGeneTics information system (http://www.imgt.org), was created in 1989 by Marie-Paule Lefranc (Montpellier University and CNRS) to manage the huge and complex diversity of these antigen receptors. IMGT built on IMGT-ONTOLOGY concepts of identification (keywords), description (labels), classification (gene and allele nomenclature) and numerotation (IMGT unique numbering), is at the origin of immunoinformatics, a science at the interface between immunogenetics and bioinformatics. IMGT/HighV-QUEST, the first web portal, and so far the only one, for the next generation sequencing (NGS) analysis of IG and TR, is the paradigm for immune repertoire standardized outputs and immunoprofiles of the adaptive immune responses. It provides the identification of the variable (V), diversity (D) and joining (J) genes and alleles, analysis of the V-(D)-J junction and complementarity determining region 3 (CDR3) and the characterization of the 'IMGT clonotype (AA)' (AA for amino acid) diversity and expression. IMGT/HighV-QUEST compares outputs of different batches, up to one million nucleotide sequences for the statistical module. These high throughput IG and TR repertoire immunoprofiles are of prime importance in vaccination, cancer, infectious diseases, autoimmunity and lymphoproliferative disorders, however their comparative statistical analysis still remains a challenge. We present a standardized statistical procedure to analyze IMGT/HighV-QUEST outputs for the evaluation of the significance of the IMGT clonotype (AA) diversity differences in proportions, per gene of a given group, between NGS IG and TR
authors SA, VG, SK and M-PL, but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the "author contributions" section.

**Competing Interests:** The authors have declared that no competing interests exist.

repertoire immunoprofiles. The procedure is generic and suitable for evaluating significance of the IMGT clonal type (AA) diversity and expression per gene, and for any IG and TR immunoprofiles of any species.

## Introduction

IMGT, the international ImMunoGeneTics information system (http://www.imgt.org) [1], was created in 1989 by Marie-Paule Lefranc, Laboratoire d'ImmunoGénétique Moléculaire LIGM (Montpellier University and CNRS) at Montpellier, France, in order to standardize and to manage the complexity and the diversity of immunogenetics data. IMGT, built on IMGT-ONTOLOGY [2], is at the origin of immunoinformatics [3], a science at the interface between immunogenetics and bioinformatics.

The adaptive immune response was acquired by jawed vertebrates (or gnathostomata) more than 450 million years ago and is found in all extant jawed vertebrate species from fishes to humans [3]. The potential antigen receptor repertoire of each individual is estimated to comprise about 10¹² different immunoglobulins (IG) or antibodies [4] and 10¹² different T cell receptors (TR) [5] per individual. This huge diversity is created by combinatorial and junctional diversity (together with somatic hypermutations for IG) and the limiting factor is only the number of B and T cells that an organism is genetically programmed to produce [3].

IG are made of two identical heavy (H) chains and two identical light (L) (kappa or lambda) chains, encoded by genes located in three major loci: the IG heavy (IGH) locus, IG kappa (IGK) locus and IG lambda (IGL) locus [4, 6]. TR are made of two chains, alpha and beta, or gamma and delta, encoded by genes located in four major loci: the TR alpha (TRA), TR beta (TRB), TR gamma (TRG) and TR delta (TRD) [5, 7] (see IMGT Repertoire, http://www.imgt.org/IMGTindex/locus.html and IMGT/GENE-DB [8]).

There are four IG or TR gene types: variable (V), diversity (D) (only for IGH, TRB and TRD), joining (J) and constant (C) genes, which define 24 IG and TR groups (e.g., IGHV, IGHD, IGHJ, ..., TRBV, TRBD, TRBJ, ...) [2, 3] (http://www.imgt.org/IMGTindex/group.html). The V, D, J and C genes contribute to the IG and TR chain synthesis [3-5]. The variable domain at the N-terminal end of each IG or TR chain results from a V-(D)-J rearrangement whereas the remaining of the chain, or constant region, is encoded by a C gene [3-5].

The analysis of the immune antigen receptor (IG and TR) repertoires has greatly benefited from the next generation sequencing (NGS) technologies. The vast amount of generated data necessitated the development of novel methods and analysis tools. IMGT/HighV-QUEST [9], a high throughput version of IMGT/V-QUEST [10-14] was implemented by IMGT in October 2010 and is the first reference NGS web portal for IG and TR. IMGT/HighV-QUEST analyzes up to 1,000,000 IG and TR sequences from NGS high throughput and deep sequencing [9, 15, 16] and compares outputs of different batches, up to one million nucleotide sequences for the statistical module. The analysis is based on the IMGT-ONTOLOGY concepts of identification, description, classification and numerotation [2].

IMGT/HighV-QUEST [9, 15] uses the same algorithm as IMGT/V-QUEST [10-14]. It identifies the variable (V), diversity (D) and joining (J) genes and alleles, describes the V-REGION mutations, determines the hot spot positions in the closest germline V gene and allele, incorporates IMGT/JunctionAnalysis for the analysis of the V-J and V-D-J junctions [17, 18] and IMGT/Automat [19, 20] for a full V-J-REGION and V-D-J-REGION annotation.
As a novel NGS functionality, IMGT/HighV-QUEST identifies clonotypes (for IG or TR) which are defined as 'IMGT clonotype (AA)' (AA for amino acid) and are characterized by a unique V-(D)-J rearrangement (IMGT genes and alleles determined at the nucleotide level), conserved CDR3-IMGT anchors (cysteine C 104, tryptophan W 118 or phenylalanine F 118), and a unique CDR3-IMGT AA junction sequence [16]. Each 'IMGT clonotype (AA)' is characterized by a selected unique representative sequence. For the first time for NGS antigen receptor data analysis, the IMGT standardized approach allows a clear distinction between the clonotype diversity (numbers of IMGT clonotypes (AA) per V, D or J gene), and the clonotype expression (numbers of sequences assigned, unambiguously, to a given IMGT clonotype (AA) per V, D or J gene) [16], making the IMGT/HighV-QUEST NGS immunoprofiles available for statistical analysis.

In this work, we analyze IMGT/HighV-QUEST NGS outputs on IMGT clonotypes (AA) diversity per gene between NGS sets with the aim of defining a standardized IMGT/HighV-QUEST statistical procedure for the evaluation of the significance of the IMGT clonotype (AA) diversity differences in proportions between high throughput repertoire immunoprofiles.

## Materials and Methods

### IMGT/HighV-QUEST outputs

IMGT/HighV-QUEST outputs are from the analysis performed in [16] which aimed to minimize experimental biases. The analysis of T cell receptors results were selected as they do not undergo somatic mutations as opposed to immunoglobulins. The choice of the 454 Genome Sequencer FLX (GSFLX) Titanium (Roche) technology allowed to obtain sequences long enough to cover the full V-DOMAIN [3]. Amplicon libraries of the corresponding TRB V-D-J transcripts were prepared using anchored 5' rapid amplification of cDNA ends (RACE) [21, 22] polymerase chain reaction (PCR) to exclude amplification bias. The quality of the results was confirmed by the identification of transcripts of pseudogenes (previously not known as being transcribed) and by the strikingly similar pattern of distribution of TRBV genes and alleles of the common IMGT clonotypes (AA) between sets of a given T cell population at different time points [16]. Sequencing data is available in the NCBI Sequence Read Archive under the accession code SRX326382 [16].

Eight from the twelve sets analyzed in [16] were selected for this study. They correspond to two T cell populations (CD4⁻ and CD4⁺) at four time points (pre-vaccination (Pre), and day 3 (d3), day 8 (d8) and day 26 (d26) post-vaccination) of a single individual vaccinated against H1N1 influenza virus. The sets are identified by barcodes or multiplex identifiers (MID) (MID1, MID2, MID4, MID5, MID7, MID8, MID10 and MID11) (Table 1) which allow to

<table><caption>Table 1. Number (nb) of in-frame *Homo sapiens* TRB IMGT clonotypes (AA) per group in CD4⁻ versus (⁻) CD4⁺ populations at four time points.</caption><thead><tr><th rowspan="2">Compared sets* CD4<sup>-</sup>(set i) ~ CD4<sup>+</sup>(set j)</th><th rowspan="2">Time points</th><th colspan="3">Nb of IMGT clonotypes (AA)** per group in CD4<sup>-</sup>(n<sub>i</sub>) ~ CD4<sup>+</sup>(n<sub>j</sub>)</th></tr><tr><th>TRBV</th><th>TRBD</th><th>TRBJ</th></tr></thead><tbody><tr><td>CD4<sup>-</sup>(MID1) ~ CD4<sup>+</sup>(MID2)</td><td>Pre</td><td>2234 ~ 2523</td><td>2124 ~ 2412</td><td>2234 ~ 2523</td></tr><tr><td>CD4<sup>-</sup>(MID4) ~ CD4<sup>+</sup>(MID5)</td><td>d3</td><td>1389 ~ 2859</td><td>1323 ~ 2751</td><td>1389 ~ 2859</td></tr><tr><td>CD4<sup>-</sup>(MID7) ~ CD4<sup>+</sup>(MID8)</td><td>d8</td><td>1309 ~ 1966</td><td>1242 ~ 1886</td><td>1309 ~ 1966</td></tr><tr><td>CD4<sup>-</sup>(MID10) ~ CD4<sup>+</sup>(MID11)</td><td>d26</td><td>1924 ~ 2875</td><td>1853 ~ 2769</td><td>1924 ~ 2875</td></tr></tbody></table>

<footnote>* Sets identified by barcodes or multiplex identifiers (MID).</footnote>

<footnote>** IMGT clonotypes (AA) are defined by a unique V-(D)-J rearrangement and characterized by their V gene and J gene [16] (therefore identical $n_i$ and $n_j$ per TRBV and TRBJ group at a given time point). The lower numbers of IMGT clonotypes (AA) for TRBD group $n_i$ ~ $n_j$ are due to the difficulty of identifying the TRBD genes and alleles, owing to their short length and high % of identity [5]. Six IMGT clonotypes (AA) with an outlier CDR3-IMGT length (> 60 AA or < 4 AA) were removed: two from MID1 (63 AA and 65 AA), two from MID2 and MID5 (2 AA) and two from MID2 and MID10 (3 AA).</footnote>
differentiate between the two T cell populations and the four time points. The reads were run with IMGT/HighV-QUEST program version 1.1.3, IMGT/V-QUEST program version 3.2.31 and IMGT/V-QUEST reference directory release 201338-1 (http://www.imgt.org) [9, 15, 16].

## IMGT/HighV-QUEST statistical procedure

**Differences in proportions.** The purpose was to compare, between two sets, the differences in proportions of IMGT clonotypes (AA) per gene of a given group (clonotype diversity). For this analysis, only IMGT clonotypes (AA) which are in-frame (no frameshift in the junction) identified by IMGT/HighV-QUEST are selected. Statistical analyses applied to data sets were performed using R software (www.r-project.org; version 3.2.0).

Let *m* be the number of genes analyzed in a given group of a locus, for example, the *Homo sapiens* TRB locus [8], 54 genes for the TRBV group, 2 genes for the TRBD group and 13 genes for the TRBJ group. Let each gene of a given group, identified by its IMGT gene name (e.g., TRBV2, TRBJ1-1...) [5, 7] be indexed by *k* (*k* = 1, ..., *m*).

Sets to be compared corresponding to T cell populations CD4⁺ and CD4⁻ were indexed by {i, j} (*i* ≠ *j*). The numbers (nb) of IMGT clonotypes (AA) per group (TRBV, TRBD or TRBJ) in CD4⁺ versus CD4⁻ populations (*n*ᵢ ~ *n*ⱼ) at four time points (Pre, d3, d8 and d26) are shown in Table 1.

Let $(X_r)_{r=1..n_i}^k$ and $(Y_s)_{s=1..n_j}^k$ be two random variables defined as follows:

$$ (X_r)_{r=1..n_i}^k = \begin{cases} 1, & \text{if the } r^{\text{th}} \text{ IMGT clonotype (AA) in the set } i \\ & \text{with the gene } k \text{ of a given group} \\ 0, & \text{else} \end{cases} $$

and

$$ (Y_s)_{s=1..n_j}^k = \begin{cases} 1, & \text{if the } s^{\text{th}} \text{ IMGT clonotype (AA) in the set } j \\ & \text{with the gene } k \text{ of a given group} \\ 0, & \text{else} \end{cases} $$

where $n_i$ (resp., $n_j$) is the number of IMGT clonotypes (AA) in the set *i* (resp., set *j*) for a given group (Table 1).

($X_r)_{r=1..n_i}^k$ (resp., $\hat{Y}_s)_{s=1..n_j}^k$) were considered to be independent and identically distributed (i. i.d.) Bernoulli random variables of parameter $p_i^k$ (resp., $\hat{p}_j^k$). The parameters $p_i^k$ and $\hat{p}_j^k$ vary in $\Theta = [0, 1]^2$ and represent the probabilities of finding at least one IMGT clonotype (AA) with the gene *k* in the set *i* and the set *j*, respectively. $\hat{p}_i^k$ and $\hat{p}_j^k$ are the estimators of $p_i^k$ and $\hat{p}_j^k$, respectively, defined as follows:

$$ \begin{cases} \hat{p}_i^k = \frac{\text{Number of IMGT clonotypes (AA) in the set } i \text{ with the gene } k \text{ of a given group}}{\text{Number of IMGT clonotypes (AA) in the set } i \text{ for a given group}} \\ = \frac{1}{n_i} \sum_{r=1}^{n_i} X_r \\ \hat{p}_j^k = \frac{\text{Number of IMGT clonotypes (AA) in the set } j \text{ with the gene } k \text{ of a given group}}{\text{Number of IMGT clonotypes (AA) in the set } j \text{ for a given group}} \\ = \frac{1}{n_j} \sum_{s=1}^{n_j} Y_s \end{cases} $$
For example, the number of IMGT clotypes (AA) with the gene k TRBV5-1 of the TRBV group is 244 in the set i MID1 and 420 in the set j MID2 (S1 Table). Knowing that the number of IMGT clotypes (AA) is 2234 in the set i MID1 and 2523 in the set j MID2 (Table 1), the estimators of $p_{\text{MID1}}^{\text{TRBV5-1}}$ and $p_{\text{MID2}}^{\text{TRBV5-1}}$ are as follows:

$$ \left\{ \begin{array}{l} \hat{p}_{\text{MID1}}^{\text{TRBV5-1}} = \frac{244}{2234} = 0.109221 \\ \hat{p}_{\text{MID2}}^{\text{TRBV5-1}} = \frac{420}{2523} = 0.166468 \end{array} \right. $$

The difference in proportions $\hat{p}_i^k - \hat{p}_j^k$ for TRBV5-1 between MID1 and MID2 is -0.057247 (S1 Table).

Approximate 95% confidence intervals (CI) [23, 24] were calculated for reference to provide a magnitude of the true difference in proportions $p_i^k - p_j^k$ as follows:

$$ (\hat{p}_i^k - \hat{p}_j^k) \pm z(1 - \frac{\alpha}{2}) \cdot \sigma \hat{p}_i^k - \hat{p}_j^k \quad (1) $$

where $z(1 - \frac{\alpha}{2})$ is the $1 - \frac{\alpha}{2}$ percentile of a standard normal distribution (assumption of the distribution proven in Step 2) (so, at $\alpha = 5\%$, $1 - \frac{\alpha}{2} = 0.975$ and $z = 1.96$) and where $\sigma \hat{p}_i^k - \hat{p}_j^k$ is the standard deviation (SD) (square root of the variance) estimated as follows:

$$ \sigma \hat{p}_i^k - \hat{p}_j^k = \sqrt{\frac{\hat{p}_i^k(1 - \hat{p}_i^k)}{n_i} + \frac{\hat{p}_j^k(1 - \hat{p}_j^k)}{n_j}} \quad (2) $$

For the example of the TRBV5-1 gene, the standard deviation is calculated as follows:

$$ \sigma \hat{p}_{\text{MID1}}^{\text{TRBV5-1}} - \hat{p}_{\text{MID2}}^{\text{TRBV5-1}} = \sqrt{\frac{0.109221(1-0.109221)}{2234} + \frac{0.166468(1-0.166468)}{2523}} = 0.009927 \text{ and the 95\% CI for the difference in proportions } \hat{p}_{\text{MID1}}^{\text{TRBV5-1}} - \hat{p}_{\text{MID2}}^{\text{TRBV5-1}} \text{ is given by: } -0.057247 \pm 1.96(0.009927) = [-0.076704, -0.037790] \text{ (S1 Table).} $$

**Significance of the differences in proportions.** The hypothesis to be tested is to determine whether the difference between two proportions is significant. In order to evaluate the significance of the differences in proportions, a two-tailed statistical test was performed by following these steps:

**Step 1.** Set up the null hypothesis $H_0^k$ (there is no difference between $p_i^k$ and $p_j^k$) and its alternative $H_1^k$ defined as follows:

$$ \left\{ \begin{array}{l} H_0^k : p_i^k - p_j^k = 0 \\ H_1^k : p_i^k - p_j^k \neq 0 \end{array} \right. $$

**Step 2.** In this study, the sample size is large enough to expect a standard normal distribution and to choose the z-score test statistic ($z = \frac{\text{score} - \mu(\text{mean})}{\sigma(\text{standard deviation SD})}$) which defines the number of standard deviations an element (here, the difference in proportions $\hat{p}_i^k - \hat{p}_j^k$) is away from the mean (here, the hypothesized difference in proportions $p_i^k - p_j^k$), with

$$ z^k = \frac{(\hat{p}_i^k - \hat{p}_j^k) - \mu_{\hat{p}_i^k - \hat{p}_j^k}}{\sigma_{\hat{p}_i^k - \hat{p}_j^k}} = \frac{(\hat{p}_i^k - \hat{p}_j^k) - (p_i^k - p_j^k)}{\sigma_{\hat{p}_i^k - \hat{p}_j^k}} \quad (3) $$

where $p_i^k - p_j^k$ replaces $\mu_{\hat{p}_i^k - \hat{p}_j^k}$ in the classical z-score formula.
To reject or accept the null hypothesis a significance level of 5% ($\alpha$ = 0.05) was chosen. This is a two-tailed test since a value on either side of the standard normal distribution ($z$ > 1.96 or $z$ < 1.96) would cause one to reject the null hypothesis.

Note that in the case of null or too small proportions (where $n_i \hat{p}_i^k < 5$, $n_i(1 - \hat{p}_i^k) < 5$ and $n_j \hat{p}_j^k < 5$, $n_j(1 - \hat{p}_j^k) < 5$), the z-score results are not considered (e.g., the case of the TRBV5-4 and TRBV7-7 genes having a null proportion in the set MID1).

Under the null hypothesis $H_0^k$, the difference $p_i^k - p_j^k$ is null and the z-score formula (3) becomes as follows:

$$z^k = \frac{\hat{p}_i^k - \hat{p}_j^k}{\sigma_{\hat{p}_i^k - \hat{p}_j^k}} \sim \mathcal{N}(0, 1) \quad (4)$$

Replacing in Eq (4) the denominator $\sigma_{\hat{p}_i^k - \hat{p}_j^k}$ Eq (2) by the pooled estimate of the standard deviation $\sigma_{\hat{p}_ij^k}$ of the weighted pooled proportion $\hat{p}_{ij}^k = \frac{n_i \hat{p}_i^k + n_j \hat{p}_j^k}{n_i + n_j}$ [25, 26], the z-score formula becomes as follows:

$$z^k = \frac{\hat{p}_i^k - \hat{p}_j^k}{\sqrt{\hat{p}_{ij}^k(1 - \hat{p}_{ij}^k)\left(\frac{1}{n_i} + \frac{1}{n_j}\right)}} \sim \mathcal{N}(0, 1) \quad (5)$$

The pooled estimation of the standard deviation is valid only if the null hypothesis is true and it is used in order to reduce the variability of the difference between the two tested sample proportions. Thus, the weighted pooled proportion allows the control of this variability by assigning more weight to the larger sample size.

For the example of TRBV5-1 gene, if the null hypothesis $H_0^k : p_{\text{TRBV5-1}}^{\text{MID1}} = p_{\text{TRBV5-1}}^{\text{MID2}}$ is true, the weighted pooled proportion $\hat{p}_{ij}^k$ is $\hat{p}_{\text{TRBV5-1}}^{\text{MID1MID2}} = \frac{(2234 \times 0.109221) + (2523 \times 0.166468)}{2234 + 2523} = 0.139583$, the pooled estimate of the standard deviation $\sigma_{\hat{p}_{ij}^k}$ is $\sigma_{\hat{p}_{\text{TRBV5-1}}^{\text{MID1MID2}}} = \sqrt{0.139583(1 - 0.139583)(\frac{1}{2234} + \frac{1}{2523})} = 0.010067$ and the test statistic is $z_{\text{TRBV5-1}} = \frac{-0.057247}{0.010067} = -5.686$ (S2 Table).

At the significance level of 5% ($\alpha$ = 0.05) level and two-tailed test, the null hypothesis $H_0^{\text{TRBV5-1}}$ is rejected since -5.686 (more than 5 standard deviations from zero) is much less than -1.96 and falls in the $\frac{\alpha}{2} = 0.025$ left zone of the standard normal distribution.

**Step 3.** Each z-score for each difference in proportions leads to a p-value obtained from computed tables. The p-value is the probability $\mathbb{P}$ of observing a value Z greater than the absolute value of the z-score observed under the null hypothesis $H_0$. The p-value allows, as the critical value (step 2), to reject or not the null hypothesis, at a defined rejection level (usually $\alpha$ = 0.01 or 0.05 or 0.1) (the critical value approach (step 2) and the p-value approach should always lead to the same conclusion). It also allows to evaluate the significance of the deviation from the null hypothesis. The p-value calculation is based on the alternative hypothesis $H_1^k$ as follows: $2\mathbb{P}(Z \ge |z|)$ for $H_1^k : p_i^k \ne p_j^k$ (i.e., $\mathbb{P}(Z \ge z)$ for $H_1^k : p_i^k > p_j^k$ and $\mathbb{P}(Z \le z)$ for $H_1^k : p_i^k < p_j^k$).

For the example of the TRBV5-1 gene, the computed p-value $2\mathbb{P}(Z \ge |5.686|)$ is 1.3E-08. Since this p-value is much less than the rejection level 0.05, the null hypothesis is rejected and the difference in proportions of the TRBV5-1 gene between MID1 and MID2 is declared to be significant at $\alpha$ = 0.05.
**Error rates control and multiple testing procedures.** The intuitive idea is to test *p*-value obtained for each gene, individually from each *z*-score test at a chosen *α* level (usually *α* = 0.05, implying that it is acceptable to have a 5% probability of incorrectly rejecting the null hypothesis). However, when multiple tests of hypotheses are conducted simultaneously, more than 5% of them are very likely to be statistically significant purely by chance; for example, when *m* tests (here, number of genes) are performed, the experimentwise significance level will be given by 1 – (1 – *α*)<sup>*m*</sup> ≈ *ma* when *α* is small and will represent the global Type I error rate (i.e., risk to incorrectly reject a true null hypothesis *H*<sub>0</sub> (false positive), here, for instance, the difference in proportions of IMGT clonotypes (AA) with the gene *k* between two sets would be declared as significant while this is not the case). For that reason, an adjustment of the *p*-values must be made through a multiple testing procedure.

There are two basic and commonly used strategies for dealing with multiple hypotheses testing by controlling different error rates which incorporate the number of performed tests: the family-wise error rate (FWER) defined by Tukey [27] and the false discovery rate (FDR) of Benjamini & Hochberg [28]. The FWER is the probability to make one or more Type I errors or false positives (i.e., when the statistical test leads to reject a null hypothesis $H_0^k$ when it is true) among all the *m* hypotheses tested, whereas the FDR is the expected proportion of Type I errors among rejected hypotheses. The FDR procedures are generally more powerful (i.e., they provide a high probability that the null hypothesis $H_0^k$ is correctly rejected by the statistical test when the alternative hypothesis $H_1^k$ is true) than the FWER procedures, considered as more conservative [28–30]. Therefore, the FDR procedures are generally applied in an explanatory context or in the case of large scale multiple testing problems. The FWER procedures are more appropriate in confirmatory or regulatory setting due to their strict control of risk of Type I errors. In some situations, the FWER is likely to be a secure measure of controlling the error rates especially when the test is falsely declared as significant. For these reasons, the choice of the adequate error measure is heavily influenced by the scientific question asked.

In this study, the following multiple testing procedures were applied: Bonferroni [31], Holm [32], Hochberg [33], ŠidákSS and ŠidákSD [34, 35], controlling the FWER, as well as Benjamini & Hochberg (BH) [28] and Benjamini & Yekutieli (BY) [36], controlling the FDR (Table 2). Each procedure belongs to one of the two main categories of algorithms for generating multiple testing adjustments: single-step and stepwise procedures (step-down or step-up algorithms). In single-step procedures, a single criterion is used to evaluate the importance of all test statistics or corresponding *p*-values without taking into consideration their order or their number, namely, a critical value is used to assess each null hypothesis $H_0^k$ independently of the results obtained by testing all the other hypotheses [37].

Stepwise procedures take into account results obtained for the previous hypotheses. They are improved in terms of power while maintaining the control of the Type I error rate. Tests

<table><caption>**Table 2. Properties of the multiple testing procedures.**</caption><thead><tr><td>Procedures</td><td>Type of control</td><td>Algorithm structure</td><td>Dependence of p-values under H<sub>0</sub></td><td>Properties</td></tr></thead><tbody><tr><td>Bonferonni [31]</td><td>FWER</td><td>Single-step</td><td>Ignorance</td><td>The most conservative</td></tr><tr><td>Šidák (SS) [34, 35]</td><td>FWER</td><td>Single-step</td><td>Independence</td><td>Less conservative than Bonferroni</td></tr><tr><td>Holm [32]</td><td>FWER</td><td>Step-down</td><td>Ignorance</td><td>Less conservative than Bonferroni</td></tr><tr><td>Šidák (SD) [34, 35]</td><td>FWER</td><td>Step-down</td><td>Dependence</td><td>Similar to Holm</td></tr><tr><td>Hochberg [33]</td><td>FWER</td><td>Step-up</td><td>Independence</td><td>Step-up of Holm</td></tr><tr><td>Benjamini & Hochberg (BH) [28]</td><td>FDR</td><td>Step-up</td><td>Independence</td><td>The least conservative</td></tr><tr><td>Benjamini & Yekutieli (BY) [36]</td><td>FDR</td><td>Step-up</td><td>Ignorance</td><td>More conservative than BH</td></tr></tbody></table>

doi:10.1371/journal.pone.0142353.t002
are considered successively either in descending order of their corresponding *p*-values (step-down) or in ascending order (step-up). For step-down procedures, null hypotheses having the largest test statistic (smallest unadjusted *p*-values) are regarded successively. As soon as one of the null hypotheses does not get rejected, the step-down procedure stops and no more hypotheses are rejected. For step-up procedures, null hypotheses having the least test statistic (largest unadjusted *p*-values) are regarded successively. As soon as one of the null hypotheses is rejected, the step-up procedure stops and all residual more significant hypotheses are rejected [34].

All procedures mentioned above take into account the total number *m* of tested hypotheses. For this reason and in order to refine the results by reducing the error risks, TRBV genes not found in the IMGT clonotypes (AA) of the two compared sets (i.e., the difference in proportions is null) were excluded when applying multiple testing procedures.

# Results

## IMGT clonotypes (AA) for 10,000 per group

As a first step, in order to permit a visual comparative analysis of IMGT clonotype (AA) diversity for each gene per group (TRBV, TRBD or TRBJ) between two sets, the numbers of IMGT clonotypes (AA) for a given gene obtained from the IMGT/HighV-QUEST outputs were normalized for 10,000 IMGT clonotypes (AA) per group and represented by bar graphs (Fig 1) for the CD4⁻ and CD4⁺ populations (juxtaposed colored top and bottom bars, respectively). The four panels correspond to the four time points (Pre, d3, d8, d26).

Differences in proportions and corresponding 95% CI were calculated for the IMGT clonotypes (AA) with a gene *k* for a given group (TRBV, TRBD or TRBJ), between the two T cell populations (CD4⁻ and CD4⁺) at four time points (Pre, d3, d8, d26). S1 Table reports, as examples, the values for the differences in proportions found significant in this study (discussed in next paragraphs).

## Significance using multiple testing procedures

In order to evaluate the significance of the differences in proportions found in the two T cell populations CD4⁻ and CD4⁺, at the four time points, the multiple testing procedures cited in Material & Methods (Table 2) were applied. Unadjusted *p*-values and adjusted *p*-values from the multiple testing procedures were calculated. In Table 3 are shown the smallest unadjusted *p*-values (rawp) (*< 0.05*) and the corresponding adjusted values after application of the FWER (Bonferroni, Holm, Hochberg, ŠidákSS and ŠidákSD) or FDR (BH and BY) controlling procedures. Adjusted *p*-values below the rejection level 0.05 under a given procedure are highlighted in bold and correspond to significant tests (Table 3).

The importance of adjusted *p*-values can be illustrated with the example of the set comparison CD4⁻(MID1) ~ CD4⁺(MID2) (Pre) (Table 3), where 17 differences in proportions of IMGT clonotypes (AA) are declared to be significant before controlling the error rate (unadjusted *p*-values (rawp) column), whereas only six are retained by all adjusting procedures (*p*-values < 0.05 and highlighted in bold). Differences in proportions for four additional genes are retained by BH (FDR controlling procedure). Similarly, in the set comparison CD4⁻(MID4) ~ CD4⁺(MID5) at d3, 24 differences in proportions are declared to be significant before adjustment of *p*-values, whereas only 10 are retained by all multiple testing procedures (Table 3). Differences in proportions for one additional gene (TRBV15) is similarly retained by Holm, Hochberg, ŠidákSD, BH and BY procedures, whereas more genes are retained by BH and BY (one and four, respectively) in line with the FDR procedures being more powerful than the FWER procedures (Table 3).
<figure><img src="image_4.png" /><figcaption>**Fig 1. Normalized bar graphs of the proportions.** The normalized bar graphs of the proportions for *Homo sapiens* TRB IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ) are shown between two T cell populations (CD4⁺ and CD4⁻) at four time points (Pre, d3, d8 and d26). IMGT clonotypes (AA) proportions are normalized for 10,000 IMGT clonotypes (AA) per group. Juxtaposed colored bars in each panel correspond to CD4⁺ (top) and CD4⁻ (bottom).</figcaption></figure>

doi:10.1371/journal.pone.0142353.g001
<table><caption>**Table 3. Unadjusted and adjusted *p*-values from multiple testing procedures.** Unadjusted *p*-values (rawp) and adjusted *p*-values from multiple testing procedures were calculated, for the significance of the differences in proportions for *Homo sapiens* TRB IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ), between two T cell populations (CD4⁺ and CD4⁻) at four time points (Pre, d3, d8 and d26).</caption><thead><tr><th rowspan="2">Compared sets [time points]</th><th rowspan="2">Unadjusted p-values (rawp)</th><th colspan="7">Adjusted p-values from multiple testing procedures</th></tr><tr><th>Bonferroni</th><th>Holm</th><th>Hochberg</th><th>ŠidákSS</th><th>ŠidákSD</th><th>BH</th><th>BY</th></tr></thead><tbody><tr><td colspan="9"><strong>CD4<sup>+</sup>(MID1) ~ CD4<sup>+</sup>(MID2) [Pre]</strong></td></tr><tr><td>Homsap TRBV5-1 F (F)</td><td>1.3E-08</td><td>7.93E-07</td><td>7.93E-07</td><td>7.93E-07</td><td>7.93E-07</td><td>7.93E-07</td><td>7.93E-07</td><td>3.72E-06</td></tr><tr><td>Homsap TRBV27 F</td><td>4.81E-07</td><td>2.93E-05</td><td>2.89E-05</td><td>2.89E-05</td><td>2.93E-05</td><td>2.89E-05</td><td>1.47E-05</td><td>6.89E-05</td></tr><tr><td>Homsap TRBJ2-1 F</td><td>1.2E-05</td><td>0.00073</td><td>0.00071</td><td>0.00071</td><td>0.00073</td><td>0.00071</td><td>0.00024</td><td>0.00114</td></tr><tr><td>Homsap TRBV7-2 F (F)</td><td>7E-05</td><td>0.00447</td><td>0.00425</td><td>0.00425</td><td>0.00446</td><td>0.00424</td><td>0.00112</td><td>0.00525</td></tr><tr><td>Homsap TRBV7-9 F (F)</td><td>0.00011</td><td>0.00697</td><td>0.00652</td><td>0.00652</td><td>0.00695</td><td>0.00649</td><td>0.00127</td><td>0.00596</td></tr><tr><td>Homsap TRBV20-1 F (F)</td><td>0.00012</td><td>0.00761</td><td>0.00699</td><td>0.00699</td><td>0.00758</td><td>0.00696</td><td>0.00127</td><td>0.00596</td></tr><tr><td>Homsap TRBV7-6 F (F)</td><td>0.00243</td><td>0.14798</td><td>0.13343</td><td>0.13343</td><td>0.13771</td><td>0.12505</td><td>0.02114</td><td>0.09928</td></tr><tr><td>Homsap TRBJ2-5 F</td><td>0.00304</td><td>0.18524</td><td>0.16398</td><td>0.16398</td><td>0.16933</td><td>0.15145</td><td>0.02315</td><td>0.10874</td></tr><tr><td>Homsap TRBV25-1 F</td><td>0.00439</td><td>0.26766</td><td>0.23256</td><td>0.23256</td><td>0.23528</td><td>0.2079</td><td>0.02974</td><td>0.13967</td></tr><tr><td>Homsap TRBJ2-7 F ORF</td><td>0.00597</td><td>0.36418</td><td>0.31045</td><td>0.31045</td><td>0.30599</td><td>0.26756</td><td>0.03642</td><td>0.17103</td></tr><tr><td>Homsap TRBV4-1 F (F)</td><td>0.00940</td><td>0.57319</td><td>0.47922</td><td>0.47922</td><td>0.4378</td><td>0.38214</td><td>0.05008</td><td>0.23518</td></tr><tr><td>Homsap TRBV6-4 F</td><td>0.01023</td><td>0.62418</td><td>0.51163</td><td>0.51163</td><td>0.46602</td><td>0.40206</td><td>0.05008</td><td>0.23518</td></tr><tr><td>Homsap TRBV11-2 F [F] (F)</td><td>0.01067</td><td>0.65102</td><td>0.52295</td><td>0.52295</td><td>0.48031</td><td>0.40889</td><td>0.05008</td><td>0.23518</td></tr><tr><td>Homsap TRBJ1-5 F</td><td>0.01314</td><td>0.80141</td><td>0.63061</td><td>0.62678</td><td>0.55368</td><td>0.46996</td><td>0.05423</td><td>0.25469</td></tr><tr><td>Homsap TRBV12-4 F (F)</td><td>0.01334</td><td>0.81348</td><td>0.63061</td><td>0.62678</td><td>0.55911</td><td>0.46996</td><td>0.05423</td><td>0.25469</td></tr><tr><td>Homsap TRBV9 F (F)</td><td>0.02552</td><td>1</td><td>1</td><td>0.97616</td><td>0.79345</td><td>0.69559</td><td>0.09731</td><td>0.45700</td></tr><tr><td>Homsap TRBV12-3 F</td><td>0.04187</td><td>1</td><td>1</td><td>0.97616</td><td>0.92639</td><td>0.85407</td><td>0.15023</td><td>0.70553</td></tr><tr><td colspan="9"><strong>CD4<sup>+</sup>(MID4) ~ CD4<sup>+</sup>(MID5) [d3]</strong></td></tr><tr><td>Homsap TRBV5-1 F (F)</td><td>1.00E-15</td><td>5.70E-14</td><td>5.70E-14</td><td>5.70E-14</td><td>5.90E-14</td><td>5.90E-14</td><td>5.70E-14</td><td>2.64E-13</td></tr><tr><td>Homsap TRBV27 F</td><td>8.72E-13</td><td>5.14E-11</td><td>5.06E-11</td><td>5.06E-11</td><td>5.14E-11</td><td>5.06E-11</td><td>2.03E-11</td><td>9.49E-11</td></tr><tr><td>Homsap TRBV7-2 F (F)</td><td>1.03E-12</td><td>6.10E-11</td><td>5.90E-11</td><td>5.90E-11</td><td>6.10E-11</td><td>5.90E-11</td><td>2.03E-11</td><td>9.49E-11</td></tr><tr><td>Homsap TRBJ2-7 F ORF</td><td>1.32E-08</td><td>7.76E-07</td><td>7.37E-07</td><td>7.37E-07</td><td>7.76E-07</td><td>7.37E-07</td><td>1.94E-07</td><td>9.05E-07</td></tr><tr><td>Homsap TRBV7-9 F (F)</td><td>5.45E-08</td><td>3.22E-06</td><td>3.00E-06</td><td>3.00E-06</td><td>3.22E-06</td><td>3.00E-06</td><td>6.43E-07</td><td>3.00E-06</td></tr><tr><td>Homsap TRBV4-1 F (F)</td><td>7.79E-08</td><td>4.60E-06</td><td>4.21E-06</td><td>4.21E-06</td><td>4.60E-06</td><td>4.21E-06</td><td>7.66E-07</td><td>3.57E-06</td></tr><tr><td>Homsap TRBJ2-5 F</td><td>9.18E-08</td><td>5.41E-06</td><td>4.86E-06</td><td>4.86E-06</td><td>5.41E-06</td><td>4.86E-06</td><td>7.73E-07</td><td>3.61E-06</td></tr><tr><td>Homsap TRBV4-3 F (F)</td><td>8.52E-05</td><td>0.00503</td><td>0.00443</td><td>0.00443</td><td>0.00502</td><td>0.00442</td><td>0.00063</td><td>0.00293</td></tr><tr><td>Homsap TRBJ2-1 F</td><td>0.00041</td><td>0.02424</td><td>0.02095</td><td>0.02095</td><td>0.02395</td><td>0.02074</td><td>0.00269</td><td>0.01256</td></tr><tr><td>Homsap TRBV20-1 F (F)</td><td>0.00061</td><td>0.03602</td><td>0.03052</td><td>0.03052</td><td>0.03539</td><td>0.03007</td><td>0.0036</td><td>0.0168</td></tr><tr><td>Homsap TRBV15 F (F)</td><td>0.00088</td><td>0.05199</td><td>0.04317</td><td>0.04317</td><td>0.05068</td><td>0.04227</td><td>0.00473</td><td>0.02204</td></tr><tr><td>Homsap TRBV7-6 F (F)</td><td>0.00190</td><td>0.11211</td><td>0.09121</td><td>0.09121</td><td>0.10615</td><td>0.08725</td><td>0.00934</td><td>0.04356</td></tr><tr><td>Homsap TRBV12-4 F (F)</td><td>0.00288</td><td>0.17004</td><td>0.13545</td><td>0.13545</td><td>0.15657</td><td>0.12685</td><td>0.01308</td><td>0.06099</td></tr><tr><td>Homsap TRBV24-1 F</td><td>0.00476</td><td>0.28063</td><td>0.21880</td><td>0.21880</td><td>0.24520</td><td>0.19694</td><td>0.02005</td><td>0.09348</td></tr><tr><td>Homsap TRBJ1-5 F</td><td>0.01167</td><td>0.68846</td><td>0.52509</td><td>0.52509</td><td>0.49968</td><td>0.41032</td><td>0.04501</td><td>0.20991</td></tr><tr><td>Homsap TRBJ2-2 F</td><td>0.01221</td><td>0.72022</td><td>0.53711</td><td>0.53711</td><td>0.51551</td><td>0.41749</td><td>0.04501</td><td>0.20991</td></tr><tr><td>Homsap TRBV9 F (F)</td><td>0.01517</td><td>0.89520</td><td>0.65244</td><td>0.65244</td><td>0.59427</td><td>0.48182</td><td>0.05266</td><td>0.24556</td></tr><tr><td>Homsap TRBV18 F</td><td>0.01726</td><td>1</td><td>0.72497</td><td>0.72497</td><td>0.64203</td><td>0.51872</td><td>0.05658</td><td>0.26384</td></tr><tr><td>Homsap TRBV25-1 F</td><td>0.01860</td><td>1</td><td>0.76243</td><td>0.76243</td><td>0.66961</td><td>0.53681</td><td>0.05775</td><td>0.26928</td></tr><tr><td>Homsap TRBJ1-3 F</td><td>0.02756</td><td>1</td><td>1</td><td>0.98127</td><td>0.80770</td><td>0.67299</td><td>0.08129</td><td>0.37909</td></tr><tr><td>Homsap TRBV14 F (F)</td><td>0.03018</td><td>1</td><td>1</td><td>0.98127</td><td>0.83607</td><td>0.69739</td><td>0.0848</td><td>0.39546</td></tr><tr><td>Homsap TRBD2 F</td><td>0.04008</td><td>1</td><td>1</td><td>0.98127</td><td>0.91049</td><td>0.78868</td><td>0.10281</td><td>0.47944</td></tr><tr><td>Homsap TRBD1 F</td><td>0.04008</td><td>1</td><td>1</td><td>0.98127</td><td>0.91049</td><td>0.78868</td><td>0.10281</td><td>0.47944</td></tr><tr><td>Homsap TRBV12-3 F</td><td>0.04966</td><td>1</td><td>1</td><td>0.98127</td><td>0.95046</td><td>0.84016</td><td>0.12207</td><td>0.56926</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3.** (Continued)</caption><thead><tr><th rowspan="2">Compared sets [time points]</th><th rowspan="2">Unadjusted p-values (rawp)</th><th colspan="7">Adjusted p-values from multiple testing procedures</th></tr><tr><th>Bonferroni</th><th>Holm</th><th>Hochberg</th><th>ŠidákSS</th><th>ŠidákSD</th><th>BH</th><th>BY</th></tr></thead><tbody><tr><td colspan="9"><strong>CD4'(MID7) ~ CD4'(MID8) [d8]</strong></td></tr><tr><td>Homsap TRBV27 F</td><td>1.36E-13</td><td>8.04E-12</td><td>8.04E-12</td><td>8.04E-12</td><td>8.04E-12</td><td>8.04E-12</td><td>8.04E-12</td><td>3.75E-11</td></tr><tr><td>Homsap TRBJ2-7 F ORF</td><td>6.45E-10</td><td>3.80E-08</td><td>3.74E-08</td><td>3.74E-08</td><td>3.80E-08</td><td>3.74E-08</td><td>1.42E-08</td><td>6.64E-08</td></tr><tr><td>Homsap TRBV5-1 F (F)</td><td>7.24E-10</td><td>4.27E-08</td><td>4.13E-08</td><td>4.13E-08</td><td>4.27E-08</td><td>4.13E-08</td><td>1.42E-08</td><td>6.64E-08</td></tr><tr><td>Homsap TRBV4-1 F (F)</td><td>6.44E-08</td><td>3.80E-06</td><td>3.60E-06</td><td>3.60E-06</td><td>3.80E-06</td><td>3.60E-06</td><td>9.49E-07</td><td>4.43E-06</td></tr><tr><td>Homsap TRBV7-2 F (F)</td><td>8.27E-07</td><td>4.88E-05</td><td>4.55E-05</td><td>4.55E-05</td><td>4.88E-05</td><td>4.55E-05</td><td>9.76E-06</td><td>4.55E-05</td></tr><tr><td>Homsap TRBJ2-5 F</td><td>1.08E-05</td><td>0.00064</td><td>0.00058</td><td>0.00058</td><td>0.00064</td><td>0.00058</td><td>0.00011</td><td>5.00E-04</td></tr><tr><td>Homsap TRBJ2-1 F</td><td>0.00024</td><td>0.0142</td><td>0.01275</td><td>0.01275</td><td>0.0141</td><td>0.01267</td><td>0.00203</td><td>0.00946</td></tr><tr><td>Homsap TRBV7-9 F (F)</td><td>5.00E-04</td><td>0.02948</td><td>0.02598</td><td>0.02598</td><td>0.02906</td><td>0.02565</td><td>0.00369</td><td>0.01718</td></tr><tr><td>Homsap TRBJ2-3 F</td><td>0.00150</td><td>0.08825</td><td>0.07628</td><td>0.07628</td><td>0.08453</td><td>0.0735</td><td>0.00981</td><td>0.04573</td></tr><tr><td>Homsap TRBJ1-1 F</td><td>0.00212</td><td>0.1249</td><td>0.10585</td><td>0.10585</td><td>0.11754</td><td>0.10054</td><td>0.01249</td><td>0.05825</td></tr><tr><td>Homsap TRBJ2-2 F</td><td>0.00471</td><td>0.27763</td><td>0.23057</td><td>0.23057</td><td>0.24292</td><td>0.20635</td><td>0.02524</td><td>0.11769</td></tr><tr><td>Homsap TRBV25-1 F</td><td>0.00914</td><td>0.53897</td><td>0.43848</td><td>0.43848</td><td>0.41809</td><td>0.35628</td><td>0.04491</td><td>0.20944</td></tr><tr><td>Homsap TRBV4-3 F (F)</td><td>0.01085</td><td>0.64019</td><td>0.50998</td><td>0.50998</td><td>0.47465</td><td>0.40116</td><td>0.04925</td><td>0.22964</td></tr><tr><td>Homsap TRBV5-4 F (F)</td><td>0.01423</td><td>0.83981</td><td>0.65477</td><td>0.64053</td><td>0.5708</td><td>0.48287</td><td>0.05599</td><td>0.26108</td></tr><tr><td>Homsap TRBV7-7 F (F)</td><td>0.01423</td><td>0.83981</td><td>0.65477</td><td>0.64053</td><td>0.5708</td><td>0.48287</td><td>0.05599</td><td>0.26108</td></tr><tr><td>Homsap TRBV9 F (F)</td><td>0.01776</td><td>1</td><td>0.78139</td><td>0.78139</td><td>0.65257</td><td>0.54544</td><td>0.06346</td><td>0.29591</td></tr><tr><td>Homsap TRBV20-1 F (F)</td><td>0.01828</td><td>1</td><td>0.78622</td><td>0.78622</td><td>0.66336</td><td>0.54774</td><td>0.06346</td><td>0.29591</td></tr><tr><td>Homsap TRBV12-3 F</td><td>0.02632</td><td>1</td><td>1</td><td>0.87592</td><td>0.79274</td><td>0.67382</td><td>0.08628</td><td>0.40233</td></tr><tr><td>Homsap TRBV18 F</td><td>0.04074</td><td>1</td><td>1</td><td>0.87592</td><td>0.91403</td><td>0.81825</td><td>0.1265</td><td>0.58987</td></tr><tr><td colspan="9"><strong>CD4'(MID10) ~ CD4'(MID11) [d26]</strong></td></tr><tr><td>Homsap TRBV27 F</td><td>2.14E-26</td><td>1.31E-24</td><td>1.31E-24</td><td>1.31E-24</td><td>1.27E-24</td><td>1.23E-24</td><td>1.31E-24</td><td>6.14E-24</td></tr><tr><td>Homsap TRBJ2-7 F ORF</td><td>4.54E-17</td><td>2.77E-15</td><td>2.72E-15</td><td>2.72E-15</td><td>2.78E-15</td><td>2.68E-15</td><td>1.38E-15</td><td>6.50E-15</td></tr><tr><td>Homsap TRBV7-2 F (F)</td><td>1.70E-14</td><td>1.06E-12</td><td>1.03E-12</td><td>1.03E-12</td><td>1.06E-12</td><td>1.03E-12</td><td>3.54E-13</td><td>1.66E-12</td></tr><tr><td>Homsap TRBV5-1 F (F)</td><td>9.10E-14</td><td>5.54E-12</td><td>5.26E-12</td><td>5.26E-12</td><td>5.53E-12</td><td>5.26E-12</td><td>1.38E-12</td><td>6.50E-12</td></tr><tr><td>Homsap TRBV4-1 F (F)</td><td>2.41E-10</td><td>1.47E-08</td><td>1.38E-08</td><td>1.38E-08</td><td>1.47E-08</td><td>1.38E-08</td><td>2.94E-09</td><td>1.38E-08</td></tr><tr><td>Homsap TRBV4-3 F (F)</td><td>1.66E-07</td><td>1.01E-05</td><td>9.31E-06</td><td>9.31E-06</td><td>1.01E-05</td><td>9.31E-06</td><td>1.69E-06</td><td>7.94E-06</td></tr><tr><td>Homsap TRBJ2-1 F</td><td>3.15E-07</td><td>1.92E-05</td><td>1.73E-05</td><td>1.73E-05</td><td>1.92E-05</td><td>1.73E-05</td><td>2.74E-06</td><td>1.29E-05</td></tr><tr><td>Homsap TRBJ2-5 F</td><td>5.49E-07</td><td>3.35E-05</td><td>2.97E-05</td><td>2.97E-05</td><td>3.35E-05</td><td>2.97E-05</td><td>4.19E-06</td><td>1.97E-05</td></tr><tr><td>Homsap TRBV20-1 F (F)</td><td>6.34E-07</td><td>3.87E-05</td><td>3.36E-05</td><td>3.36E-05</td><td>3.87E-05</td><td>3.36E-05</td><td>4.30E-06</td><td>2.02E-05</td></tr><tr><td>Homsap TRBV6-2 F (P)</td><td>1.17E-06</td><td>7.15E-05</td><td>6.09E-05</td><td>6.09E-05</td><td>7.15E-05</td><td>6.09E-05</td><td>6.84E-06</td><td>3.21E-05</td></tr><tr><td>Homsap TRBV25-1 F</td><td>1.23E-06</td><td>7.53E-05</td><td>6.29E-05</td><td>6.29E-05</td><td>7.53E-05</td><td>6.29E-05</td><td>6.84E-06</td><td>3.21E-05</td></tr><tr><td>Homsap TRBJ2-3 F</td><td>0.00005</td><td>0.00289</td><td>0.00237</td><td>0.00237</td><td>0.00288</td><td>0.00236</td><td>0.00024</td><td>0.00113</td></tr><tr><td>Homsap TRBV12-4 F (F)</td><td>0.00011</td><td>0.00667</td><td>0.00536</td><td>0.00536</td><td>0.00665</td><td>0.00534</td><td>0.00051</td><td>0.00241</td></tr><tr><td>Homsap TRBV18 F</td><td>0.00048</td><td>0.02900</td><td>0.02282</td><td>0.02282</td><td>0.02859</td><td>0.02256</td><td>0.00207</td><td>0.00973</td></tr><tr><td>Homsap TRBJ1-4 F</td><td>0.00080</td><td>0.04853</td><td>0.03739</td><td>0.03739</td><td>0.04739</td><td>0.03672</td><td>0.00324</td><td>0.01519</td></tr><tr><td>Homsap TRBV6-1 F</td><td>0.00154</td><td>0.09386</td><td>0.07078</td><td>0.07078</td><td>0.08965</td><td>0.06838</td><td>0.00587</td><td>0.02755</td></tr><tr><td>Homsap TRBV7-9 F (F)</td><td>0.00301</td><td>0.18336</td><td>0.13527</td><td>0.13320</td><td>0.16777</td><td>0.12670</td><td>0.01026</td><td>0.04818</td></tr><tr><td>Homsap TRBV7-6 F (F)</td><td>0.00303</td><td>0.18467</td><td>0.13527</td><td>0.13320</td><td>0.16885</td><td>0.12670</td><td>0.01026</td><td>0.04818</td></tr><tr><td>Homsap TRBV13 F (F)</td><td>0.00666</td><td>0.40628</td><td>0.28640</td><td>0.28640</td><td>0.33478</td><td>0.24975</td><td>0.02138</td><td>0.10042</td></tr><tr><td>Homsap TRBV12-3 F</td><td>0.01313</td><td>0.80117</td><td>0.55162</td><td>0.55162</td><td>0.55357</td><td>0.42609</td><td>0.04006</td><td>0.18812</td></tr><tr><td>Homsap TRBV29-1F (F)</td><td>0.01522</td><td>0.92828</td><td>0.62393</td><td>0.62393</td><td>0.60758</td><td>0.46673</td><td>0.04420</td><td>0.20759</td></tr><tr><td>Homsap TRBJ1-2 F</td><td>0.02597</td><td>1.00000</td><td>1.00000</td><td>0.98697</td><td>0.79914</td><td>0.65096</td><td>0.07201</td><td>0.33818</td></tr><tr><td>Homsap TRBV11-2 F [F] (F)</td><td>0.02979</td><td>1.00000</td><td>1.00000</td><td>0.98697</td><td>0.84195</td><td>0.69256</td><td>0.07901</td><td>0.37105</td></tr><tr><td>Homsap TRBV6-4 F</td><td>0.03543</td><td>1.00000</td><td>1.00000</td><td>0.98697</td><td>0.88923</td><td>0.74607</td><td>0.08710</td><td>0.40904</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3.** (Continued)</caption><thead><tr><th rowspan="2">Compared sets [time points]</th><th rowspan="2">Unadjusted p-values (rawp)</th><th colspan="7">Adjusted p-values from multiple testing procedures</th></tr><tr><th>Bonferroni</th><th>Holm</th><th>Hochberg</th><th>ŠidákSS</th><th>ŠidákSD</th><th>BH</th><th>BY</th></tr></thead><tbody><tr><td>Homsap TRBV5-4 F (F)</td><td>0.03570</td><td>1.00000</td><td>1.00000</td><td>0.98697</td><td>0.89109</td><td>0.74607</td><td>0.08710</td><td>0.40904</td></tr><tr><td>Homsap TRBV9 F (F)</td><td>0.04172</td><td>1.00000</td><td>1.00000</td><td>0.98697</td><td>0.92569</td><td>0.78435</td><td>0.09505</td><td>0.44637</td></tr><tr><td>Homsap TRBJ1-3 F</td><td>0.04207</td><td>1.00000</td><td>1.00000</td><td>0.98697</td><td>0.92733</td><td>0.78435</td><td>0.09505</td><td>0.44637</td></tr></tbody></table>

Unadjusted *p*-values significant at 5% (< 0.05) are reported. Adjusted *p*-values significant at 5% (< 0.05) under the given controlling procedure are highlighted in bold.

The functionality is defined for each allele of a gene according to the rules described in http://www.imgt.org/IMGTScientificChart/SequenceDescription/IMGTfunctionality.html and can be functional, ORF or pseudogene.

- F (Functional) if the coding region has an open reading frame without stop codon, and if there is no described defect in the splicing sites, recombination signals and/or regulatory elements.

- ORF (Open Reading Frame) if the coding region has an open reading frame.

- P (Pseudogene) if the coding region has stop codon(s) and/or frameshift mutation(s).

Note that 6 pseudogenes of the TRBV group (TRBV1, TRBV3-2, TRBV12-1, TRBV12-2, TRBV21-1 and TRBV26), that cannot give productive TR chains, were excluded from the study.

In the table the functionality of all known alleles of each gene is reported and shown as follows:

- between parentheses, (F) and (P), when the accession number refers to rearranged genomic DNA or cDNA and the corresponding germline gene has not yet been isolated.

- between brackets, [F] and [P], when the accession number refers to genomic DNA, but not known as being germline or rearranged.

doi:10.1371/journal.pone.0142353.t003

## Visualization plots

In order to visualize the statistical analysis results obtained using the multiple testing procedures, two types of plots (line graphs and scatter plots, Fig 2) were generated for each comparison using the Bioconductor (http://www.bioconductor.org) R package multtest [38].

The line graphs (Fig 2, left panel) represent the number of rejected hypotheses (number of significant values in Table 3) against the Type I error rate (α-level). Dotted lines represent unadjusted *p*-values (rawp) whereas colored lines represent adjusted *p*-values of the seven procedures. These graphs allow to visualize the number of rejected null hypotheses for a chosen α-level under a given procedure. Thus, by fixing the smallest adjusted *p*-values obtained by the procedures at a Type I error rate (α-level = 0.05) (vertical line in Fig 2), the numbers of rejected null hypotheses (numbers of adjusted *p*-values in bold in Table 3), for example with the BH procedure (red line), are 10 (Pre), 16 (d3), 13 (d8) and 21 (d26), whereas with the Bonferroni procedure (purple line), they are 6 (Pre), 10 (d3), 8 (d8) and 15 (d26). The line graphs also allow to identify, for a selected number of rejected null hypotheses, the α-level which is required to get that number for a given procedure.

The scatter plots (Fig 2, right panel) display negative decimal logarithms (-log₁₀) of unadjusted *p*-values (black symbols) and adjusted *p*-values obtained by each multiple testing procedure (colored symbols) against the test statistics z-scores on abscissa (values in S2 Table).

As the *p*-values range from 0 to 1, their negative decimal logarithms are used in the plots to provide a larger visual representation. Thus, a *p*-value of 0.05 corresponds to -log₁₀(0.05) ≈ 1.3 (horizontal line in Fig 2). Therefore, significant differences in proportions of IMGT clonotypes (AA) with a gene *k* at a Type I error rate α-level = 0.05 for the different procedures correspond to those represented with a -log₁₀(*p*-value) > 1.3 and with a z-score less than -1.96 (< -1.96) for negative differences or with a z-score greater than 1.96 (> 1.96) for positive differences. Differences are negative if the proportion of IMGT clonotypes (AA) is greater in the second set
<figure><img src="image_5.png" /><figcaption>**Fig 2. Multiple testing procedures visualization plots.** Multiple testing procedures visualization plots are displayed for comparison of the differences in proportions for IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ), between two T cell populations (CD4⁺ and CD4⁻) at four time points (Pre, d3, d8 and d26). The following procedures: Bonferroni, Holm, Hochberg, ŠidákSS and ŠidákSD, BH and BY were applied. Left panel (A, C, E, G): Line graphs showing the number of rejected null hypotheses against the Type I error rate. Dotted lines represent unadjusted *p*-values (rawp) whereas colored lines represent adjusted *p*-values of the seven procedures. A vertical line corresponds to a Type I error rate (α-level) at 0.05 (significance level of 5%). Right panel (B, D, F, H): Negative decimal logarithms (-log₁₀) of unadjusted *p*-values (black symbols) and adjusted *p*-values (colored symbols) against the test statistic z-scores. Two areas in scatter plots (top left and top right) correspond to significant differences in proportions and they are delimited at a significance level of 5% (0.05) by -log₁₀(*p*-values) > 1.3 (horizontal line) and by z-scores (< -1.96 for negative differences or > 1.96 for positive differences) (vertical line).</figcaption></figure>

doi:10.1371/journal.pone.0142353.g002
(CD4+) than in the first set (CD4⁻) ($p_{\{CD4⁻\}}^k < p_{\{CD4⁺\}}^k$) and, alternatively, differences are posi-
tive if the proportion of the IMGT clonotypes (AA) is greater in the first set (CD4⁻) than in the
second set (CD4⁺) ($p_{\{CD4⁻\}}^k > p_{\{CD4⁺\}}^k$).

The genes with the highest significant differences in proportion of IMGT clonotypes (AA) between the two T cell populations CD4⁻ and CD4⁺ are in the top left ($p_{\{CD4⁻\}}^k < p_{\{CD4⁺\}}^k$) and top right ($p_{\{CD4⁻\}}^k > p_{\{CD4⁺\}}^k$) of each scatter plot (Fig 2, right panel) and can be identified in S2 Table. For examples, the three genes identified with highly significant negative difference in proportions for CD4⁻(MID1) ~ CD4⁺(MID2) (Pre) are: TRBV5-1 (z-score = -5.686, -log₁₀(-rawp) = 7.886, -log₁₀(adjusted p-values) = 6.101 (FWER procedure and BH) and 5.429 (BY)), TRBV7-2 (z-score = -3.965, -log₁₀(rawp) = 4.135, log₁₀(adjusted p-values) from 2.350 to 2.373 (FWER), 2.952 (BH) and 2.280 (BY)) and TRBV20-1 (z-score = -3.837, -log₁₀(rawp) = 3.904, -log₁₀(adjusted p-values) from 2.119 to 2.157 (FWER), 2.897 (BH) and 2.225 (BY)) (S2 Table and top left of Fig 2B). The three genes identified with the highly significant positive differences for CD4⁻(MID1) ~ CD4⁺(MID2) (Pre) are: TRBV27 (z-score = 5.034, -log₁₀(rawp) = 6.318, -log₁₀(adjusted p-values) from 4.532 to 4.540 (FWER procedures), 4.833 (BH) and 4.162 (BY)), TRBJ2-1 (z-score = 4.378, -log₁₀(rawp) = 4.922, log₁₀(adjusted p-values) from 3.136 to 3.151 (FWER), 3.613 (BH) and 2.942 (BY)) and TRBV7-9 (z-score = 3.858, -log₁₀(rawp) = 3.942, -log₁₀(adjusted p-values) from 2.157 to 2.187 (FWER), 2.897 (BH) and 2.225 (BY)) (S2 Table and top right of Fig 2B).

Difference in proportions graph

The difference in proportions graph (Fig 3) shows the differences in proportions with signifi-
cance and confidence intervals (CI) of IMGT clonotypes (AA) with a gene *k* of a given group
(TRBV, TRBD or TRBJ) between the two T cell populations (CD4⁻ and CD4⁺) at the four time
points (Pre, d3, d8 and d26). The point in the middle of each 95% CI bar corresponds to the
value of the difference in proportions of IMGT clonotypes (AA) with a gene *k* found in two dif-
ferent sets. Differences in proportions and 95% CI lower and upper bounds are reported in S1
Table for all cases with unadjusted p-values < 0.05 of Table 3. Negative differences in propor-
tions of the IMGT clonotypes (AA) (proportion greater in the second set (CD4⁺) than in the
first set (CD4⁻)) are displayed on the left of the vertical reference line (abscissa equal to 0).
Alternatively, positive differences in proportions of the IMGT clonotype (AA) (proportion
greater in the first set (CD4⁻) than in the second set (CD4⁺)) are displayed on the right of the
vertical line. CI bars are colored based on the interpretation test before adjustment (legend in
Fig 3): red for non-significant difference (i.e., -1.96 < z-score < 1.96) or light blue for signifi-
cant difference (i.e., z-score < -1.96 or z-score > 1.96). If a CI bar does not contain 0 (i.e., does
not cross the vertical reference line), this is a clear indication that there is a statistically signifi-
cant difference. More the point in the middle of a CI bar is far from the vertical reference line,
more the significance of the difference in proportions is important. Lower levels of variance
yield shorter CI bar widths, and therefore more accurate estimates.

Considering the results of multiple testing procedures, the CI bars declared as significant by
the seven adjustment procedures (Table 3) are colored in dark blue. If the significance is vali-
dated by two or more multiple testing procedures, CI bars are colored in pink. If the signifi-
cance is validated only by the BH procedure, CI bars are colored in green.

As examples, the negative differences in proportions of the IMGT clonotypes (AA) of the TRBV5-1 (-0.057247), TRBV7-2 (-0.030195) and TRBV20-1 (-0.044506) genes (S1 Table) between sets CD4⁻ and CD4⁺ at Pre are displayed on the left of the vertical reference line in panel Pre(MID1) ~ Pre(MID2) (Fig 3) ($p_{\{CD4⁻\}}^k < p_{\{CD4⁺\}}^k$). The positive difference in
CD4⁻~CD4⁺

<figure><img src="image_6.png" /><figcaption>**Fig 3. Differences in proportions graph.** The difference in proportions graph with significance and confidence interval (CI) bars are shown for *Homo sapiens* TRB IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ) between two T cell populations (CD4⁻ and CD4⁺) at four time points (Pre, d3, d8 and d26). Negative differences in proportions ($p_{(CD4⁻)} - p_{(CD4⁺)} < 0$) are shown on the left of the vertical line (abscissa equal to 0) and positive differences ($p_{(CD4⁻)} - p_{(CD4⁺)} > 0$) are shown on the right of the vertical line. CI bars colors indicated in the legend correspond to the test interpretation before adjustment of *p*-values (rawp) (significant differences validated by the seven procedures (All_p): dark blue, by two or more multiple testing procedures (Min_2p): pink, and only by BH (Only_BH): green).</figcaption></figure>

doi:10.1371/journal.pone.0142353.g003
proportions of IMGT clonotypes (AA) of the TRBV7-9 (+0.018678), TRBV27 (+0.031817) and TRBJ2-1 (+0.048132) are displayed on the right of the vertical line ($p_{\{CD4^-\}}^k > p_{\{CD4^+\}}^k$). The CI bars of these 6 genes colored in dark blue indicate that the significance was validated by all procedures (i.e., adjusted *p*-values < 0.05, Table 3). In the same panel, the CI bars of TRBV7-6, TRBV25-1, TRBJ2-5 and TRBJ2-7 colored in green indicate that the significance is validated by only the BH procedure (Table 3). The CI bars colored in pink, present in the other panels, indicate that the significance of the difference in proportions is validated by two or more procedures, for example, TRBV7-6 gene at d3 validated by the BH (*P* = 0.009) and the BY procedures (*P* = 0.043) (Table 3, d3) and (Fig 3: ‘d3’ panel).

## Synthesis graph

To facilitate the comparison with the experimental results, a synthesis graph (Fig 4) is proposed that combines the display, in parallel panels, of the normalized bar graphs of the proportions and the differences in proportions with significance and confidence intervals (CI bars), described in the previous sections. A few examples are given to illustrate the interpretation and comparisons (Figs 4 and 5).

The important differences in proportions of IMGT clonotypes (AA) of the TRBV5-1, TRBV7-2, TRBV20-1 and TRBJ2-5 genes with higher clonal type diversity in CD4⁺ compared to CD4⁻ at the four time points (Pre, d3, d8, d26) was validated by all multiple testing procedures at an α-level = 0.05 (Table 3) (Fig 4) (dark blue CI bars) with the exception of TRBV20-1 at d8 (only validated before adjustment: unadjusted *p*-value = 0.01, light blue CI bar) and of TRBJ2-5 (Pre) where the significance was validated only by the BH method (*P* = 0.023, Table 3) and Fig 4 (green CI bar in ‘Pre’ panel).

The higher clonal type diversity of the TRBJ2-7 gene in CD4⁻ compared to CD4⁺ observed before vaccination (Pre), at d3, d8 and d26 post-vaccination was validated by all multiple testing procedures (*P* < 10⁻⁶, Table 3) (Fig 4, dark blue CI bar), whereas before vaccination (Pre) it is validated only by the BH procedure (*P* = 0.036, Table 3).

The slight differences in proportions of the IMGT clonotypes (AA) of the TRBV4-1 and TRBV4-3 genes were declared as significant at d3, d8 and d26 post-vaccination by all adjustment methods, except for TRBV4-3 at day 8, when the significance was validated only by the BH procedure. Slight significant differences in clonal type diversity were also observed at the four time points for TRBV7-9 but only declared significant at Pre, d3, d8 by all adjustment methods and at d26 by the BH and BY procedures. For TRBV12-4 gene, the significance of the clonal type diversity was validated by all adjusting procedure at d26, but only by the BH procedure at d3. For TRBJ2-2 gene at d3 and d8 post-vaccination, differences in proportions were declared as significant only by the BH procedure.

A potentially important finding by comparing CD4⁻ and CD4⁺ populations is the fact that the vaccination did not change the sign of significant differences in proportions found before vaccination. In others words, if the proportion of IMGT clonotypes (AA) having a specific gene *k* in CD4⁻ exceeds that in CD4⁺ at baseline (i.e., $p_{\{CD4^-\}} - p_{\{CD4^+\}} > 0$), this positive sign remains unchanged until d26. The repertoire characteristics in terms of IMGT clonal type diversity (AA) per gene is therefore kept at any given time point and the differences in proportion-sobserved are maintained between the two T cell populations CD4⁻ and CD4⁺.

## Discussion

In this study, we define a standardized IMGT/HighV-QUEST procedure for the evaluation of the significance of the IMGT clonotypes (AA) diversity differences in proportions between NGS immunoprofiles. We statistically analyzed the diversity of IMGT clonotypes (AA) of each
CD4⁻~CD4⁺ [Pre,d3]

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Synthesis graph (part 1: Pre, d3). The synthesis graph is displayed for Homo sapiens TRB IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ) between two T cell populations (CD4⁻ and CD4⁺) at two time points (Pre and d3), with for each time point, two panels. Left panel: normalized bar graph (IMGT clonotypes (AA) proportions normalized for 10,000 IMGT clonotypes (AA) per group), with juxtaposed colored bars corresponding to CD4⁻ (top) and CD4⁺ (bottom). Right panel: difference in proportions graph with significance and confidence interval (CI) bars. CI bar colors correspond to the test interpretation before adjustment of *p*-values (rawp) (non-significant: red, significant: light blue) and after adjustment by the multiple testing procedures (significant differences validated by the seven procedures (All_p): dark blue, by two or more multiple testing procedures (Min_2p): pink, and only by BH (Only_BH): green).</figcaption></figure>

doi:10.1371/journal.pone.0142353.g004
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Synthesis graph (part 2: d8, d26). The synthesis graph is displayed for *Homo sapiens* TRB IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ) between two T cell populations (CD4⁻ and CD4⁺) at two time points (d8 and d26), with for each time point, two panels. Left panel: normalized bar graph (IMGT clonotypes (AA) proportions normalized for 10,000 IMGT clonotypes (AA) per group), with juxtaposed colored bars corresponding to CD4⁻ (top) and CD4⁺ (bottom). Right panel: difference in proportions graph with significance and confidence interval (CI) bars. CI bar colors correspond to the test interpretation before adjustment of *p*-values (rawp) (non-significant: red, significant: light blue) and after adjustment by the multiple testing procedures (significant differences validated by the seven procedures (All_p): dark blue, by two or more multiple testing procedures (Min_2p): pink, and only by BH (Only_BH): green).</figcaption></figure>

doi:10.1371/journal.pone.0142353.g005
gene of a given group (TRBV, TRBD or TRBJ) between two T cells populations (CD4⁻ and CD4⁺) at four time points pre-vaccination (Pre) and d3, d8 and d26 post-vaccination of a single individual vaccinated against H1N1 influenza virus. We first applied a test statistic z-score using a weighted pooled proportion to compare the IMGT clonotypes (AA) diversity of each gene of a given group between sets.

To control introduced errors when testing many null hypotheses, we applied FWER (Bonferroni, Holm, Hochberg, ŠiddákSS and Šidák SD) and FDR (BH and BY) controlling procedures. We observed that FWER procedures, and especially Bonferroni, Holm and Hochberg, gave very similar adjusted *p*-values (i.e., same number of rejected null hypotheses). This is most likely due to the independence structure of the unadjusted *p*-values [30, 34]. We also noted that the FDR (BY and especially BH) procedures gave more rejected null hypotheses confirming that they are more powerful and less conservative than the FWER procedures.

These observations are important for the choice of the controlling procedures for multiple significance testing as shown in this study for the evaluation of the statistical significance of IMGT clonotypes (AA) diversity per gene. FWER controlling procedures (Bonferroni, Holm, Hochberg, ŠidákSS and ŠidákSD) have been criticized for genomics studies as being very stringent methods by avoiding Type I errors and therefore leading to loss of information in some treated cases with many Type II errors (β-level) (i.e., when the statistical test of the hypotheses leads to accept a null hypothesis $H_0^k$ when it is false) [30]. Bonferroni, the most popular method, is also the most conservative one according to several simulation studies [39, 40], however, it is valid even if the test statistics are dependent. Šidák (ŠidákSS), a single-step procedure based on the Bonferroni method, is used in the case of independently distributed unadjusted *p*-values and it does not ensure a good control for random distributions of the test statistics [34]. Šidák (ŠidákSD), a step-down procedure is developed to improve the conservative results of single-step methods. Holm, a step-down procedure, is also considered as less conservative than the Bonferroni method and it makes no hypotheses on the dependence of the unadjusted *p*-values. Hochberg, a step-up modification of the Holm procedure, is more powerful than the Bonferroni correction and makes assumptions on the dependence structure of the unadjusted *p*-values. Our study confirms a similarity in the results found by using the Bonferroni, Holm, Hochberg, ŠiddákSS and Šidák SD procedures in terms of the number of rejected null hypotheses.

The improvements of other procedures for corrections of multiple comparison tests following the Bonferroni procedure focus, in general, on the dependence or ignorance of test statistics under null hypotheses. Despite the assumption of the independence between IMGT clonotypes (AA) with a gene *k* of a given group in this study (i.e., independence of unadjusted *p*-values), strong consideration will be given to results validated by all procedures to keep a Type I error control as reliable as possible.

Our study confirms that the control of the FDR is more powerful than the control of the FWER, with the BH method being even less conservative than the BY procedure as previously shown [28, 30, 34, 39]. For this reason, a special consideration should be given to results obtained by the BH procedure.

IMGT/HighV-QUEST represents a major breakthrough for the NGS analysis and the comparison of the antigen receptor V domain repertoires and immunoprofiles of the adaptive immune response [9, 15]. Indeed, IMGT/HighV-QUEST provides standardized and high-quality outputs based on the IMGT-ONTOLOGY concepts of identification, description, classification and numerotation [2]. IMGT/HighV-QUEST comprises more than 500 features per sequence allowing a detailed characterization of the antigen receptor (IG and TR) V domain repertoires. Comparison between batches up to 1 million sequences can be performed online.
The analysis of these repertoires are invaluable for the comparison of immune profiles in protective (vaccination, infections, cancers) or pathogenic (autoimmunity, lymphoproliferative disorders) immune responses [41].

One main feature of immunoprofiles is the comparison of the IMGT clonotypes (AA) diversity and expression. The IMGT/HighV-QUEST statistical procedure implemented for evaluating the significance of the IMGT clonotypes (AA) adds a further major step in the repertoire analysis. Based on the IMGT-ONTOLOGY concepts, the analysis is generic. It is applicable to the IMGT clonotypes (AA) diversity of the four TR loci (TRA, TRB, TRG, TRD) [5] and of the three IG loci (IGH, IGK, IGL) [4]. IMGT/HighV-QUEST outputs provide the number of sequences per IMGT clonotype (AA) (IMGT clonotype expression) based on the same standardized criteria and IMGT-ONTOLOGY concepts [16]. The procedure is therefore suitable for evaluating significance of the IMGT clonotype (AA) diversity as well as the significance of IMGT clonotype (AA) expression changes of a gene for a given group. More generally, this procedure can also be used for statistical comparison of IMGT clonotypes (AA) of different B or T cell populations in order to detect significant changes in immunoprofiles observed in normal or pathological situations, in one individual or between individuals.

## Supporting Information

**S1 Table.** Number (nb) of *Homo sapiens* TRB IMGT clonotypes (AA) and normalized nb for 10,000. The number (nb) of *Homo sapiens* TRB IMGT clonotypes (AA) and normalized nb for 10,000 per gene and per group (TRBV, TRBD or TRBJ), differences in proportions and corresponding 95% confidence intervals (CI) are given between two T cell populations (CD4⁻ and CD4⁺) at four time points (Pre, d3, d8 and d26).
(PDF)

**S2 Table.** Test statistics (z-scores) and negative decimal logarithms (-log₁₀). The test statistics (z-scores) and negative decimal logarithms (-log₁₀) for unadjusted (rawp) and adjusted p-values for *Homo sapiens* TRB IMGT clonotypes (AA) with a gene of a given group (TRBV, TRBD or TRBJ) are given between two T cell populations (CD4⁻ and CD4⁺) at four time points (Pre, d3, d8 and d26).
(PDF)

## Acknowledgments

We thank Patrice Duroux and Arthur Lavoie for IMGT/HighV-QUEST. We are grateful to Gérard Lefranc for helpful comments. IMGT is Academic Institutional Member of the International Medical Informatics Association (IMIA) and of the Global Alliance for the Genomics and Health (GA4GH).

## Author Contributions

Conceived and designed the experiments: DM M-PL. Performed the experiments: SA. Analyzed the data: SA. Wrote the paper: SA M-PL. Supervised the project: DM VG SK M-PL.

## References

1. Lefranc M-P, Giudicelli V, Duroux P, Jabado-Michaloud J, Folch G, Aouinti S, et al. IMGT, the international ImMunoGeneTics information system 25 years on. Nucleic Acids Research. 2015; 43 part 11: e413–422. doi: 10.1093/nar/gku1056

2. Giudicelli V, Lefranc M-P. IMGT-ONTOLOGY 2012. Frontiers in genetics. 2012; 3:79. doi: 10.3389/fgene.2012.00079 PMID: 22654892
3. Lefranc M-P. Immunoglobulin (IG) and T cell receptor (TR) genes: IMGT and the birth and rise of immunoinformatics. Frontiers in Immunology. 2014; 5(22). doi: 10.3389/fimmu.2014.00022 PMID: 24600447

4. Lefranc M-P, Lefranc G. Immunoglobulin FactsBook. Academic Press; 2001. pp.1-458.

5. Lefranc M-P, Lefranc G. The T cell receptor FactsBook. Academic Press; 2001. pp.1-398.

6. Lefranc M-P. Nomenclature of the human immunoglobulin genes. In: Coligan JE, Bierer BE, Margulies DE, Shevach EM, Strober W, editors. Current Protocols in Immunology. John Wiley and Sons, Hoboken N.J. 2000. pp. A.1P.1–A.1P.37.

7. Lefranc M-P. Nomenclature of the human T cell receptor genes. In: Coligan JE, Bierer BE, Margulies DE, Shevach EM, Strober W, editors. Current Protocols in Immunology. John Wiley and Sons, Hoboken NJ. 2000. pp. A.1O.1–A.1O.23.

8. Giudicelli V, Chaume D, Lefranc M-P. IMGT/GENE-DB: a comprehensive database for human and mouse immunoglobulin and T cell receptor genes. Nucleic Acids Research. 2005; 33:D256–D261. doi: 10.1093/nar/gki010 PMID: 15608191

9. Alamyar E, Giudicelli V, Li S, Duroux P, Lefranc M-P. IMGT/High V-QUEST: The IMGT web portal for immunoglobulin (IG) or antibody and T cell receptor (TR) analysis from NGS high throughput and deep sequencing. Immunome Research. 2012; 8:26. Available: www.imgt.org/PDF/ImmunomeRes/8_1-15_2012.pdf

10. Giudicelli V, Chaume D, Lefranc M-P. IMGT/V-QUEST, an integrated software for immunoglobulin and T cell receptor V-J and V-D-J rearrangement analysis. Nucleic Acids Research. 2004; 32:W435–40.

11. Giudicelli V, Lefranc M-P. Interactive IMGT on-line tools for the analysis of immunoglobulin and T cell receptor repertoires. In: Veskler BA, editor. New Research on Immunology. Nova Science Publishers Inc.; New York, NY, USA. 2005. pp. 77–105.

12. Brochet X, Lefranc M-P, Giudicelli V. IMGT/V-QUEST: the highly customized and integrated system for IG and TR standardized V-J and V-D-J sequence analysis. Nucleic acids research. 2008; 36(suppl 2): W503–508.

13. Giudicelli V, Lefranc M-P. IMGT standardized analysis of immunoglobulin rearranged sequences. In: Ghia P, Rosenquist R, Davi F, editor. Immunoglobulin Gene Analysis in Chronic Lymphocytic Leukemia. Milan, Italy: Wolters Kluwer Health Italy Ltd; 2009. pp. 33–52.

14. Giudicelli V, Brochet X, Lefranc M-P. IMGT/V-QUEST: IMGT standardized analysis of the immunoglobulin (IG) and T cell receptor (TR) nucleotide sequences. Cold Spring Harbor Protocols. 2011; 6: 695–715.

15. Alamyar E, Duroux P, Lefranc M-P, Giudicelli V. IMGT tools for the nucleotide analysis of immunoglobulin IG and T cell receptor (TR) V-(D)-J repertoires, polymorphisms, and IG mutations: IMGT/V-QUEST and IMGT/HighV-QUEST for NGS. Methods in Molecular Biology. 2012; 882:569–604. doi: 10.1007/978-1-61779-842-9_32 PMID: 22665256

16. Li S, Lefranc M-P, Miles JJ, Alamyar E, Giudicelli V, Duroux P, et al. IMGT/HighV-QUEST paradigm for T cell receptor IMGT clonotype diversity and next generation repertoire immunoprofiling. Nature Communications. 2013; 4:2333. doi: 10.1038/ncomms3333 PMID: 23995877

17. Yousfi Monod M, Giudicelli V, Chaume D, Lefranc M-P. IMGT/JunctionAnalysis: the first tool for the analysis of the immunoglobulin and T cell receptor complex V-J and V-D-J JUNCTIONs. Bioinformatics (Oxford, England). 2004; 20 Suppl 1:i379–385. doi: 10.1093/bioinformatics/bth945

18. Giudicelli V, Lefranc M-P. IMGT/junctionanalysis: IMGT standardized analysis of the V-J and V-D-J junctions of the rearranged immunoglobulins (IG) and T cell receptors (TR). Cold Spring Harb Protoc. 2011; 6:716–725.

19. Giudicelli V, Protat C, Lefranc M-P. The IMGT strategy for the automatic annotation of IG and TR cDNA sequences: IMGT/Automat. In: 'Proceedings of the European Conference on Computational Biology (ECCB 2003)', Data and Knowledge Bases, Poster DKB_31. Institut National de Recherche en Informatique et en Automatique. 2003. pp. 103–104.

20. Giudicelli V, Chaume D, Jabado-Michaloud J, Lefranc M-P. Immunogenetics sequence annotation: the strategy of IMGT based on IMGT-ONTOLOGY. Studies in health technology and informatics. 2005; 116:3–8. PMID: 16160227

21. Frohman MA, Dush MK, Martin GR. Rapid production of full-length cDNAs from rare transcripts: amplification using a single gene-specific oligonucleotide primer. Proc Nati Acad Sci USA. 1988; 85(23):8998–9002. doi: 10.1073/pnas.85.23.8998

22. Freeman JD, Warren RL, Webb JR, Nelson BH, Holt RA. Profiling the T-cell receptor beta-chain repertoire by massively parallel sequencing. Genome research. 2009; 19(10):1817–1824. doi: 10.1101/gr.092924.109 PMID: 19541912

23. Pieres AM, Amado C. Interval estimators for a binomial proportion: comparison of twenty methods. REVSTAT—Statistical Journal. 2008; 6:165–197. Available: ine-Inx04.ine.pt/revstat/pdf/rs080204.pdf
24. Zhang H, Andres Gutierrez Rojas H, Cepeda Cuervo E. Confidence and credibility intervals for the difference of two proportions. Revista Colombiana de Estadistica. 2010; 33(1):63-88. Available: www.scielo.org.co/pdf/rce/v33n1/v33n1a05.pdf

25. LeBlanc DC. Statistics: concepts and applications for science. Jones and Bartlett; 2004.

26. De Muth JE. Basic statistics and pharmaceutical statistical applications. 2nd ed. Pharmacy Education Series. Taylor & Francis; 2006.

27. Tukey JW. The problem of multiple comparisons. In: Mimeographed Notes. Princeton, NJ: Princeton University; 1953.

28. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological). 1995; 57(1):289-300.

29. Zhong S, Tian L, Li C, Storch KF, Wong WH. Comparative analysis of gene sets in the Gene Ontology space under the multiple hypothesis testing framework. In: Proc IEEE Comput Syst Bioinform Conf. 2004:425-435. PMID: 16448035

30. Goeman JJ, Solari A. Tutorial in biostatistics: multiple hypothesis testing in genomics. Statistics in Medicine. 2012. pp.1-27.

31. Bonferroni CE. Teoria statistica delle classi e calcolo delle probabilità. Pubblicazioni del R Istituto Superiore di Scienze Economiche e Commerciali di Firenze. 1936; 8:3-62.

32. Holm S. A simple sequentially rejective multiple test procedure. Scandinavian journal of statistics. 1979; 6:65-70. Available: www.ime.usp.br/~abe/lista/pdf4R8xPVzCnX.pdf

33. Hochberg Y. A sharper Bonferroni procedure for multiple tests of significance. Biometrika. 1988; 75 (4):800-802.

34. Dudoit S, van der Laan MJ. Multiple testing procedures with application to genomics. Springer Series in Statistics; 2008.

35. Šidák Z. Rectangular confidence regions for the means of multivariate normal distributions. Journal of the American Statistical Association. 1967; 62(318):626-633. doi: 10.1080/01621459.1967.10482935

36. Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Annals of statistics. 2001; 29(4):1165-1188.

37. Dudoit S, Shaffer J, Boldrick J. Multiple hypothesis testing in microarray experiments. Statistical Science. 2003; 18(1):71-103.

38. Pollard KS, Dudoit S, van der Laan MJ. Multiple testing procedures: R multtest package and applications to genomics, in Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Springer; 2005.

39. Narum SR. Beyond Bonferroni: Less conservative analyses for conservation genetics. Conservation Genetics. 2006; 7(5):783-787. doi: 10.1007/s10592-006-9189-7

40. Korn EL, Troendle JF, McShane LM, Simon R. Controlling the number of false discoveries: application to high-dimensional genomic data. Journal of Statistical Planning and Inference. 2004; 124(2):379-398. doi: 10.1016/S0378-3758(03)00211-8

41. Giudicelli V, Duroux P, Lavoie A, Aouinti S, Lefranc M-P, Kossida S. From IMGT-ONTOLOGY to IMGT/HighV-QUEST for NGS Immunoglobulin (IG) and T cell Receptor (TR) Repertoires in Autoimmune and Infectious Diseases. Autoimmun Infec Dis. 2015. 1(1):