# A TCRβ Repertoire Signature Can Predict Experimental Cerebral Malaria

Encarnita Mariotti-Ferrandiz, Hang-Phuong Pham, Sophie Dulauroy, Olivier

Gorgette, David Klatzmann, Pierre-André Cazenave, Sylviane Pied, Adrien

Six

► To cite this version:

Encarnita Mariotti-Ferrandiz, Hang-Phuong Pham, Sophie Dulauroy, Olivier Gorgette, David Klatzmann, et al.. A TCRβ Repertoire Signature Can Predict Experimental Cerebral Malaria. PLoS ONE, 2016, 11 (2), pp.e0147871. 10.1371/journal.pone.0147871 . hal-01284361

HAL Id: hal-01284361

https://hal.sorbonne-universite.fr/hal-01284361

Submitted on 7 Mar 2016

**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
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# A TCRβ Repertoire Signature Can Predict Experimental Cerebral Malaria

Encarnita Mariotti-Ferrandiz¹,²,³,⁴,⁵✉, Hang-Phuong Pham⁴,⁵✉✉a, Sophie Dulauroy¹,²,³✉b, Olivier Gorgette¹,²,³✉c, David Klatzmann⁴,⁵, Pierre-André Cazenave¹,²,³✉d, Sylviane Pied¹,²,³✉d, Adrien Six¹,²,³,⁴,⁵*

1 Sorbonne Universités, UPMC Univ Paris 06, URA 1961 CNRS, F-75005, Paris, France, 2 CNRS, URA 1961 CNRS, F-75005, Paris, France, 3 Institut Pasteur, Immunophysiopathologie infectieuse, F-75015, Paris, France, 4 Sorbonne Universités, UPMC Univ Paris 06, UMR 7211, UPMC Immunology-Immunopathology-Immunotherapy, F-75013, Paris, France, 5 Inserm, U959, Immunology-Immunopathology-Immunotherapy, F-75013, Paris, France

✉ These authors contributed equally to this work.

✉a Current address: ILTOO Pharma, iPEPS—ICM Hôpital Pitié Salpêtrière, F-75013, Paris, France

✉b Current address: Développement des tissus lymphoïdes, CNRS URA 1961, Institut Pasteur, F-75015, Paris, France

✉c Current address: Plate-forme de microscopie ultrastructurale, Institut Pasteur, F-75015, Paris, France

✉d Current address: CIIL, CNRS 8204-INSERM 1019-Université Lille 2, Institut Pasteur de Lille, F-59000, Lille, France

* adrien.six@upmc.fr

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

OPEN ACCESS

**Citation:** Mariotti-Ferrandiz E, Pham H-P, Dulauroy S, Gorgette O, Klatzmann D, Cazenave P-A, et al. (2016) A TCRβ Repertoire Signature Can Predict Experimental Cerebral Malaria. PLoS ONE 11(2): e0147871. doi:10.1371/journal.pone.0147871

**Editor:** Leonardo Jose de Moura Carvalho, Oswaldo Cruz Institute (IOC-Fiocruz), BRAZIL

**Received:** November 4, 2015

**Accepted:** December 4, 2015

**Published:** February 4, 2016

**Copyright:** © 2016 Mariotti-Ferrandiz 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:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work has been funded by Assistance Publique-Hôpitaux de Paris, Université Pierre and Marie Curie (Paris VI), LabEx Transimmunum (ANR-11-IDEX-0004-02), LabEx Parafrap (ANR-11-LABX-0024) and LIA CNRS DBT SIGID ("Systems Immunology and genetics of infectious diseases"). EMF & HPP have been supported by a doctoral fellowship from Ministère de la Recherche et de la Technologie, as well as by Fondation de France and Institut Pasteur for EMF.

## Abstract

Cerebral Malaria (CM) is associated with a pathogenic T cell response. Mice infected by *P. berghei* ANKA clone 1.49 (PbA) developing CM (CM⁺) present an altered PBL TCR repertoire, partly due to recurrently expanded T cell clones, as compared to non-infected and CM⁻ infected mice. To analyse the relationship between repertoire alteration and CM, we performed a kinetic analysis of the TRBV repertoire during the course of the infection until CM-related death in PbA-infected mice. The repertoires of PBL, splenocytes and brain lymphocytes were compared between infected and non-infected mice using a high-throughput CDR3 spectratyping method. We observed a modification of the whole TCR repertoire in the spleen and blood of infected mice, from the fifth and the sixth day post-infection, respectively, while only three TRBV were significantly perturbed in the brain of infected mice. Using multivariate analysis and statistical modelling, we identified a unique TCRβ signature discriminating CM⁺ from CTR mice, enriched during the course of the infection in the spleen and the blood and predicting CM onset. These results highlight a dynamic modification and compartmentalization of the TCR diversity during the course of PbA infection, and provide a novel method to identify disease-associated TCRβ signature as diagnostic and prognostic biomarkers.

## Introduction

Cerebral malaria (CM) represents a global health disease caused by *Plasmodium falciparum* infection. Despite efforts made in controlling *Plasmodium* infection spreading in the last decade, its burdens remains extensive, endemically accounting for 30% of the 627,000 infection-related deaths as estimated in 2012 [1]. CM is obviously associated with neurological
**Competing Interests:** The authors have declared that no competing interests exist.

features caused by the binding onto endothelial cells (EC) of *P. falciparum* parasitized red blood cells (pRBC) leading to their sequestration in the brain microvessels [2,3]. Intravascular leukocytes and platelets infiltration has been also observed in brains of Malawian children dead of CM [4]. Mouse models of experimental CM (ECM), greatly contribute to the description of the immune response in cerebral malaria, confirmed the major role of T lymphocytes in the neuropathogenesis [5]. Particularly, this emphasizes the major implication of both CD4⁺ and CD8⁺ Tαβ cells in the development of ECM [6–11]. Indeed, we and others observed the sequestration of Tαβ lymphocytes in the brain of mice developing CM (CM⁺) [12,13]. However, evidence regarding the natural antigenic specificities of these infiltrating T cells is still poor. Studies using recombinant parasite constitutively expressing the OVA peptide showed that pathogenic CD8⁺ T cells reaching the brain are specific for this exogenous peptide [14], confirming that “parasite-specific” CD8⁺ T cells are induced during infection. Very recently, two groups identified PbA epitopes recognized by different CD8⁺ T cells, each of which being characterized by different T cell receptors (TCR), which in turn reach and might damage the brain [15–17]. Although, most of those epitopes are associated with an enrichment of specific CD8⁺ T cells in both spleen and brain of PbA infected mice, none of them can protect mice from ECM outcome. Given that *Plasmodium* is characterized by a differential pattern of protein expression through his life-cycle and a high diversity of molecules, including antigen, superantigen and mitogen [18–21], it is conceivable that ECM outcome is the results of a synergic action of these several molecules leading to inappropriate responses that, in turn, scramble or divert the protective appropriate response. As a consequence, T cell repertoire might be profoundly altered in contrast with a more classical restricted clonal response.

In fact, we previously showed that blood TCRβ repertoire of CM⁺ mice is greatly perturbed compared to healthy mice and also to infected mice without cerebral symptoms. This perturbation is partly due to recurrently expanded T cell clones [22]. However, it remains unclear whether those modifications are the cause or the consequence of the disease. In order to address the quality of lymphocyte responses during the course of experimental CM infection, we described their antigen-specific receptor diversity, produced by somatic DNA rearrangements of V, (D) and J segments later spliced to C segments [23], using CDR3 spectratyping and the ISEApeaks strategies [24–26] on blood and spleen lymphocytes, from day 3 post-infection (p-i) until the ECM-related death of PbA-infected mice. Additionally, we characterized the whole brain Tαβ cell repertoire in naïve and PbA-infected mice. Using a microarray-derived analysis and prediction modelling, we looked for TCRβ peak signatures. Our results showed that splenic and blood TCRβ repertoires are progressively and broadly modified simultaneously with disease development with spleen modifications appearing before blood modifications. Importantly, we explored the whole TCRαβ repertoire in mouse brain and we showed that, although peculiar in naïve mice, there are few but major modifications following the infection, suggesting a particular response in the brain of PbA-infected mice. Finally, we identified a list of TCRβ peaks forming a signature associated with ECM development and appearing during the course of the infection. Altogether, these data strongly support the idea that T cell diversity as a whole must be taken into account for drug and vaccine development against infectious diseases.

# Material and Methods

## Mice and Parasites

Seven-week-old B10.D2 mice were purchased from Harlan UK Limited one week before infection and housed in filter-topped cages under specific pathogen-free conditions in the Institut Pasteur animal facilities of the Immunulogy Department at Institut Pasteur (Paris, France),
under a 12h:12h light:dark cyle, and received the same food and water provided by the animal
facility staff. Food and water bottle were changed daily and made available 24h a day for all the
animals of the same cage during the whole experiment timeline. Up to five sex-matched mice
were grouped per cage and followed during the experiment timeline. The clone 1.49L of *Plas-
modium berghei* ANKA was kindly given by Dr. Walliker (Institute of Genetics, Edinburg, UK)
and is maintained in our laboratory on C57BL/6J female mice. This clone induces in mice a
neurological syndrome partly mimicking the one of human CM. Blood stages of the parasite
were cryopreserved in liquid nitrogen as stabilates in Alserver's solution containing 10%
glycerol.

Infection

Mice infection was induced in 75 mice by intraperitoneal injection of 10⁶ parasitized Red Blood Cells (pRBC). 35 mice were sacrificed at day 3 (n = 5), 4 (n = 10), 5 (n = 10) and 6 (n = 10) post-infection, constituting four groups for the kinetic study. A total of 40 mice were sacrificed at the stage of cerebral malaria later referred to as CM+. Animals under continuous human observation were classified as CM+ when decreased body temperature and ataxia or hemiplegia or paraplegia or convulsions were recorded. Animals may have experienced the suffering associated with CM+ signs, though not excessive since they were sacrificed whenever diagnosed as such. To prevent excessive suffering, we implemented a shift of 2 people every 8 hours (3 groups of 2 people per 24h) to ensure human presence, animal observation and handling during the whole experiment duration. To minimize animal suffering, animals where under continuous observation 24h per day during the whole experiment duration by 2-people shifts every 8h, immediately sacrificed whenever diagnosed with CM+ signs and handled according to recommended regulations. At the same time, 24 B10.D2 mice which received either PBS or uninfected C57BL/6 Red Blood Cells (RBC) were sacrificed to constitute the control group (CTR). Parasitemia was measured on Giemsa-stained thin blood smears from day 4 p-i for all infected mice. All individuals included in this study displayed parasite positive staining from day 4 until ECM development, with 1 to 20% of parasitized RBCs (data not shown).

Ethics Statements

Institut Pasteur animal facility was accredited by the French Ministry of Agriculture to perform experiments on live mice, in application of the French (Decree 87-848 issued on 19/10/1987) and European (Directive 86/609/CEE) regulations on care and protection of Laboratory Animals. All animal experiments were approved and conducted in accordance with the Institut Pasteur Biosafety Committee (Paris) and performed in compliance with the NIH Animal Welfare Assurance #A5476-01 issued on 02/07/2007. All efforts were made to minimize animal suffering, mice euthanasia was performed using Carbon Dioxide flow in chamber for 3 min.

Cell Preparation

For each group (CTR, CM+) and subgroups (day3, 4, 5 and 6 p-i), blood and spleen were
removed. Blood was obtained on heparin by retroorbital or intracardiac punction. Mononu-
clear cells were isolated on Ficoll-Hypaque gradient (Pharmacia, France). Spleen was removed
and cells suspended in 3% FCS-PBS. Brain was harvested from CTR, day 6 p-i and CM+ ani-
mals and maintained for thirty minutes in Hepes medium containing 0.05% of Collagenase
(Sigma). Cells were then isolated on 30% Percoll gradient in DMEM medium and suspended
in 3% FCS-PBS. Residual RBCs were removed from all the samples by hypotonic shock using
ammonium chloride (ACK) lysis buffer for 1 to 2 minutes at room temperature. Cell
preparations were then washed twice with 3% FCS-PBS. Lymphoid cells were counted using
Malassez cell in presence of eosin to exclude dead cells.

TCRB Repertoire

In order to avoid bias, samples were treated randomly. Depending on the sample, total RNA
was extracted from 250,000 to 8.10⁶ mononuclear cells using the TRI REAGENT kit (Molecu-
lar Research Center, Cincinnati, Ohio). 20 µg of glycogen (Roche, Meylan, France) was used to
ensure optimal precipitation of RNA and pellet visualization. cDNA was synthesised for 1h30
at 42°C using retrotranscriptase for AM virus (Roche-Boehringer) and an inhibitor of RNAse
(RNAsin, Promega) was added to avoid RNA degradation. PCR were done using cDNA corre-
sponding to 250,000 mononuclear cells. Protocols for TCR TRBV-TRBC CDR3 spectratyping
have already been explained elsewhere in detail [22,24]. TRBC and TRBV primer sequences
were as described earlier [22]. As TRBV12-3 [27], TRBV24 [28] and TRBV21 [29] are not
functional in B10.D2 mice, they were not analysed. TRBV16 was excluded from this study
afterwards due to technical problem with the corresponding primer. PCR products were loaded
on a 36-well ABI373 or 96-well ABI377 automated sequencer (Applied Biosystems, Foster city,
CA) and separated according to their nucleotide length forming a profile of peaks for each
primer combination, spaced by 3 nucleotides as expected for in-frame transcripts. Each peak
corresponds to a CDR3 length. The Immunoscope software [24] was used to obtain peak area
and nucleotide length and CDR3 profile displays from sequencer raw data. IMGT nomencla-
ture has been used for TRBV genes [30].

CDR3 Spectratype Analysis

We used the ISEApeaks software package (2000–2002 Institut Pasteur, Paris, France) to extract, smooth, manage and analyse the data [25,26]. For each CDR3 length profile, the peak distribution is calculated as the percentage of each peak, obtained by dividing its area by the total area of all peaks within the profile. Briefly, for each TRBV-TRBC combination analysed, a reference repertoire is computed as the average CDR3 peak distribution of the samples belonging from the control group. Then, the distance between the peak distribution of each sample and that of the reference repertoire is calculated for each TRBV-TRBC combination. This distance, named DBV-BC, is the perturbation index reflecting the perturbation of the repertoire against a reference repertoire [31,32]. Since perturbation score distributions are often skewed, we log-transformed them for statistical analysis. Then, for each sample of each experimental group, the average of the DBV-BC (µDBV-BC) perturbations of all the TRBV-TRBC combinations is calculated, reflecting the global perturbation of the repertoire. DBV-BC perturbations indices range from 0 (identical profiles) to 100 (completely different profiles). To deal with possible missing values while keeping statistical power for unbiased conclusions, we applied multiple imputation as proposed by Rubin et al. [33]: (1). Missing data were imputed randomly 1000 times to produce 1000 complete datasets using predictive mean matching (ppm) method [34,35]; (2). On each of 1000 complete datasets, linear regressions were used to model the relationship between perturbation scores and day post-infection for each TRBV; (3). Estimated regression coefficient (or slope) and its standard error were pooled from 1000 analyses to give a final result; (4). A test is performed to assess whether the regression coefficient is significantly different from 0.

Multivariate methods were used to analyse data such as Principal Component Analysis (PCA) and unsupervised clustering methods of perturbation scores to study globally the modification of repertoire between organs and after infection and explore the underlying structure of data. We defined global perturbation index as the average perturbation score across all
TRBV (μμDBV-BC). ANOVA test was used to compare the global perturbation index within days post-infection and organ. To control the false discovery rate in multiple testing, we adjusted p-values using Benjamini-Hochberg's method [36]. Statistical significant level was fixed at α = 0.05 (type I error: the probability to reject the null hypothesis when it is true).

TCRβ Peak Signature Discovery

Gene set enrichment analysis, a microarray-based method, allows to determine whether an a
priori gene set is significantly enriched in one of two biological conditions. For each organ,
TCRβ peaks were first ranked according to t-statistic based on pooled coefficient computed for
differential expression between uninfected and CM+ mice by multiple imputation process as
described above. For each signature, GSEA provided normalized enrichment score, p-value
and q-value for multiple testing problem (details on how GSEA works are reported in Subra-
manian et al. [37]). TCRβ peak sets or clusters were formed using pvclust, an unsupervised
bootstrap-based version of hierarchical classification tree [38], which assesses p-values indicat-
ing how strong are the cluster. Significantly enriched clusters identified by GSEA are called sig-
natures. To test the predictive power of significantly enriched signatures, random forest
models were trained on our datasets and cross-validated on an independent dataset (24). Miss-
ing datasets were imputed 100 times to produce 100 complete datasets using the above imputa-
tion process. A random forest model was trained on each of the 100 complete datasets with 100
trees each (intermediate models) and combined within a single model (final model). This
model was tested on a test dataset. For cross-validation, random forest model was trained on
the test dataset and tested on the training dataset. Prediction accuracy was set as the good pre-
diction rate based on the confusion matrix (S1 Fig).

Statistical Analysis Tools

Statistical analysis was performed using R platform v3.0.2 with the following packages in addition: mice 2.17 [39], ade4 1.5-2, survival 2.37-4, pvclust 1.2-2, randomForest 4.6-7 [40]. These tools are available at CRAN repository http://www.r-project.org. GSEA 2.0.13 software was installed from http://www.broadinstitute.org/gsea.

Results

Compartmentalized Repertoire Diversity during PbA Infection

Our previous study showed that blood Tαβ cell repertoire is highly perturbed in PbA-infected B10.D2 mice developing ECM (CM+) [22]. We aimed at determining whether these modifications occur in the early days before the disease and could be used as a signature of the development of ECM in mice. For this purpose, we analysed the TCRβ diversity of B10.D2 mice from day 3 until ECM-related death. As described in Material and Methods and in Fig 1A, we infected 75 mice with 1.10⁶ PbA pRBC, among which 35 were sacrificed at day 3 (5 mice), 4 (10 mice), 5 (10 mice) and 6 (10 mice) post-injection (p-i) constituting four groups for the kinetic study and 40 mice were killed when they developed ECM (as described in material and methods), constituting the CM+ group. The control uninfected group was composed by 24 mice half of which received PBS or 1.10⁶ non-parasitized RBCs. Spleen, blood and brain were harvested as described in Material and Methods, leucocytes were isolated and following RNA extraction and cDNA synthesis, we performed the combined CDR3 spectratyping technique and ISE/Apeaks strategy in order to evaluate the perturbation of the repertoire. Perturbation score was calculated against the average repertoire of the spleen control uninfected group as the polyclonal reference repertoire for all the groups (see Material and Methods and [26,31]).
<figure><img src="image_4.png" /><figcaption>**Fig 1. Kinetic analysis of the TCR TRBV-TRBC repertoire during the course of PbA infection.** (A) Experimental procedure showing the preparation of samples from both infected and non-infected mice, as well as the parallel analysis of mice during the course of the infection ("kinetic" samples) and at the time of ECM onset ("CM+" samples). Days post-infection (days p-i) indicate the time at which animals of the corresponding groups are sacrificed. Organs harvested for each group are indicated. (B) Modification of the TRBV-TRBC repertoire in spleen, blood and brain of B10.D2 mice during the course of *Plasmodium berghei* ANKA infection. Average DBV-BC perturbations across all TRBVs and individuals in each group (µµDBV-BC) in the spleen (black), the blood (red) and the brain (green) are shown for the control uninfected (CTR), day 3 p-i (d3), day 4 p-i (d4), day 5 p-i (d5), day 6 p-i (d6) and CM+ groups. DBV-BC perturbations were computed with ISEApeaks using CTR Spleen as the reference group.</figcaption></figure>

doi:10.1371/journal.pone.0147871.g001

Perturbation score distributions were log-transformed and multiple imputation was applied to correct for bias due to missing data. Descriptive analyses were performed on the average of 1000 complete datasets. Fig 1B represents the average of global perturbation score (µµDBV-BC) of splenocytes, PBLs and brain repertoires during the course of infection. The perturbation increased in all organs compared to the organ related CTR groups and was different between organs (two-way ANOVA p-values showed significant perturbation between organs (p<0.0001) and between day p-i (p<0.0001)). These results suggest a compartmentalization of the response to PbA infection. We compared, for each compartment, the global perturbation scores of infected groups against the CTR group. In order to explore globally the repertoire modification across all TRBV, we performed a principal component analysis (PCA) on perturbation scores computed on spleen and blood repertoire (Fig 2). PCA plots show a separation of day 5, day 6 p-i and CM+ groups (negative PCA scores) from CTR, day 3 and day 4 p-i groups (positive PCA scores) in the spleen (Fig 2A). In the blood, day 6 p-i was separated from the other groups (Fig 2B). These observations were confirmed by t-tests showing that the perturbation index is significantly different from day 5 (p<0.0001) in the spleen (Fig 2C) and day 6 (p = 0.0003) in the blood (Fig 2D), until the development of ECM (p<0.0001 for both compartments) when compared to CTR group. Finally, the perturbation index is higher in CM+ blood compared to CM+ spleen (Fig 2C and 2D), as shown previously [22]. These results indicate that during infection, the T cell repertoire is modified early in the spleen compared to the blood, suggesting a dynamic with time of the modifications between both compartments.
<figure><img src="image_5.png" /><figcaption>**Fig 2. Differential kinetics of TCR TRBV-TRBC repertoire perturbation in the spleen and blood of infected mice.** (A-B) Kinetic representation of global perturbation scores across all TRBVs in the spleen (A) and the blood (B) using principal component analysis (PCA). Progressive modification of the repertoire is diagrammed by the shift of day p-i-related groups from the right to the left on the first PCA component (PC1). Colors correspond to analyzed groups. (C-D) Mean DBV-BC perturbations across all TRBVs (µDBV-BC) in the spleen (C) and the blood (D). DBV-BC were computed as in Fig 1B. Black dots represent individual mouse global perturbation scores. Red dots represent the average global perturbation score for each group. Statistical comparisons were performed using the two-way ANOVA test for the difference between organs and between infected groups. Tests were significant for organs (p<0.0001) and day post infection (p<0.0001) at α = 0.05.</figcaption></figure>

doi:10.1371/journal.pone.0147871.g002

Progressive but Massive Modifications of Spleen and Blood Repertoires
during PbA Infection

We further analysed the perturbation for each TRBV-TRBC combinations accounting for a
total of 2079 variables. As previously, we used a multiple imputation process (see Material and
Methods). 1000 complete datasets were generated and linear regressions were fitted for each
TRBV-TRBC combination to compare perturbation score in each day p-i to uninfected group.
Regression coefficients estimated over complete datasets were pooled to produce a unique coef-
ficient and its confidence interval. Using such information, t-tests were performed to compare
each infected group to the CTR in the corresponding organ. Results are shown in the Table 1.
Except for TRBV14, all other TRBVs are deeply altered. In the spleen, significant alterations
were observed at day 5 p-i for 8 out of 20 TRBVs to which 7 other TRBVs were added at day 6
p-i. TRBV31 is transiently perturbed at day 5. In blood, significant perturbations appeared at
day 6 for five TRBVs followed by an alteration of the overall repertoire, except for TRBV13-2,
TRBV14. Interestingly, before the development of ECM, some TRBVs perturbed in the blood
are not perturbed in the spleen, suggesting that the kinetic of the modifications is proper to
each compartment, the blood containing circulating cells from the whole body. Altogether,
these results show that TCR diversity is highly modified not only in CM+ mice but also during
the course of PbA infection, supporting the hypothesis that the alteration of the TCRβ reper-
toire is an upstream and rather early process with regards to disease development, and thus
probably a cause of ECM onset. In addition, the alteration of several TRBVs indicates multiple
immunogenic sources.
<table><caption>Table 1. Estimated regression coefficient and 95% confidence interval of log perturbation values between infected and control groups in spleen and blood.</caption><thead><tr><th>TRBV</th><th>Group</th><th>D3</th><th>D4</th><th>D5</th><th>D6</th><th>CM+</th></tr></thead><tbody><tr><td rowspan="3">1</td><td rowspan="2">Spleen</td><td>0.07<br/>[-0.46;0.6]</td><td>0.25<br/>[-0.14;0.64]</td><td>0.52<br/>[0.15;0.9]*</td><td>0.66<br/>[0.27;1.05]**</td><td>1.24<br/>[0.94;1.54]***</td></tr><tr><td>-0.06<br/>[-0.57;0.45]</td><td>-0.21<br/>[-0.62;0.21]</td><td>-0.02<br/>[-0.43;0.39]</td><td>0.32<br/>[-0.11;0.75]</td><td>1.16<br/>[0.84;1.47]***</td></tr><tr><td>Blood</td><td>0.2<br/>[-0.22;0.61]</td><td>0.3<br/>[-0.01;0.6]</td><td>0.23<br/>[-0.06;0.53]</td><td>0.37<br/>[0.07;0.68]*</td><td>0.61<br/>[0.37;0.84]***</td></tr><tr><td rowspan="3">2</td><td rowspan="2">Spleen</td><td>0.38<br/>[-0.04;0.8]</td><td>0.15<br/>[-0.18;0.47]</td><td>0.36<br/>[0.03;0.7]</td><td>0.47<br/>[0.11;0.83]*</td><td>0.77<br/>[0.5;1.04]***</td></tr><tr><td>-0.05<br/>[-0.48;0.38]</td><td>-0.18<br/>[-0.16;0.53]</td><td>-0.29<br/>[-0.05;0.64]</td><td>0.54<br/>[0.18;0.9]*</td><td>1.05<br/>[0.78;1.32]***</td></tr><tr><td>Blood</td><td>0.15<br/>[-0.31;0.6]</td><td>0.3<br/>[-0.06;0.66]</td><td>0.1<br/>[-0.26;0.46]</td><td>0.75<br/>[-0.08;0.68]</td><td>0.52<br/>[0.45;1.01]***</td></tr><tr><td rowspan="3">3</td><td rowspan="2">Spleen</td><td>0.09<br/>[-0.37;0.56]</td><td>0.1<br/>[-0.25;0.44]</td><td>0.54<br/>[0.21;0.87]*</td><td>0.72<br/>[0.38;1.07]***</td><td>1.02<br/>[0.76;1.29]***</td></tr><tr><td>-0.05<br/>[-0.48;0.38]</td><td>0.18<br/>[-0.16;0.53]</td><td>0.29<br/>[-0.05;0.64]</td><td>0.54<br/>[0.18;0.9]*</td><td>1.05<br/>[0.78;1.32]***</td></tr><tr><td>Blood</td><td>0.15<br/>[-0.31;0.6]</td><td>0.4<br/>[0.01;0.79]</td><td>0.2<br/>[-0.19;0.59]</td><td>0.75<br/>[0.36;1.14]**</td><td>0.52<br/>[0.22;0.82]**</td></tr><tr><td rowspan="3">4</td><td rowspan="2">Spleen</td><td>0.15<br/>[-0.38;0.68]</td><td>0.4<br/>[0.01;0.79]</td><td>0.2<br/>[-0.19;0.59]</td><td>0.3<br/>[-0.08;0.68]</td><td>0.73<br/>[0.45;1.01]***</td></tr><tr><td>-0.05<br/>[-0.48;0.38]</td><td>0.18<br/>[-0.16;0.53]</td><td>0.2<br/>[-0.19;0.59]</td><td>0.54<br/>[0.18;0.9]*</td><td>0.52<br/>[0.45;1.01]***</td></tr><tr><td>Blood</td><td>0.15<br/>[-0.31;0.6]</td><td>0.3<br/>[-0.06;0.66]</td><td>0.1<br/>[-0.26;0.46]</td><td>0.3<br/>[-0.08;0.68]</td><td>0.52<br/>[0.45;1.01]***</td></tr><tr><td rowspan="3">5</td><td rowspan="2">Spleen</td><td>0.16<br/>[-0.27;0.59]</td><td>0<br/>[-0.32;0.31]</td><td>0.37<br/>[0.07;0.68]*</td><td>0.56<br/>[0.25;0.88]**</td><td>0.68<br/>[0.44;0.93]***</td></tr><tr><td>-0.3<br/>[-0.67;0.07]</td><td>-0.13<br/>[-0.43;0.17]</td><td>-0.41<br/>[-0.71;-0.12]</td><td>0.01<br/>[-0.31;0.33]</td><td>0.59<br/>[0.35;0.82]***</td></tr><tr><td>Blood</td><td>0.11<br/>[-0.29;0.52]</td><td>0.06<br/>[-0.29;0.4]</td><td>0.56<br/>[0.21;0.9]*</td><td>0.31<br/>[-0.05;0.67]</td><td>0.14<br/>[-0.14;0.41]</td></tr><tr><td rowspan="3">12-1</td><td rowspan="2">Spleen</td><td>-0.41<br/>[-0.9;0.09]</td><td>-0.15<br/>[-0.51;0.21]</td><td>0.03<br/>[-0.32;0.38]</td><td>0.31<br/>[-0.05;0.67]</td><td>0.59<br/>[0.35;0.82]***</td></tr><tr><td>-0.03<br/>[-0.6;0.54]</td><td>0.1<br/>[-0.37;0.57]</td><td>0.49<br/>[0.03;0.96]</td><td>0.46<br/>[0.03;0.89]</td><td>0.76<br/>[0.38;1.14]***</td></tr><tr><td>Blood</td><td>0.11<br/>[-0.29;0.52]</td><td>0.06<br/>[-0.29;0.4]</td><td>0.56<br/>[0.21;0.9]*</td><td>0.32<br/>[-0.07;0.7]</td><td>0.5<br/>[0.23;0.77]***</td></tr><tr><td rowspan="3">12-2</td><td rowspan="2">Spleen</td><td>-0.18<br/>[-0.72;0.37]</td><td>-0.31<br/>[-0.78;0.17]</td><td>0.43<br/>[0.01;0.85]</td><td>0.45<br/>[0.03;0.88]</td><td>0.56<br/>[0.21;0.91]**</td></tr><tr><td>-0.03<br/>[-0.6;0.54]</td><td>0.1<br/>[-0.37;0.57]</td><td>0.49<br/>[0.03;0.96]</td><td>0.46<br/>[0.03;0.89]</td><td>0.76<br/>[0.38;1.14]***</td></tr><tr><td>Blood</td><td>0.18<br/>[-0.72;0.37]</td><td>-0.31<br/>[-0.78;0.17]</td><td>0.43<br/>[0.01;0.85]</td><td>0.45<br/>[0.03;0.88]</td><td>0.56<br/>[0.21;0.91]**</td></tr><tr><td rowspan="3">13-1</td><td rowspan="2">Spleen</td><td>-0.22<br/>[-0.65;0.22]</td><td>-0.37<br/>[-0.69;-0.05]</td><td>-0.16<br/>[-0.47;0.15]</td><td>0.22<br/>[-0.1;0.54]</td><td>0.28<br/>[0.03;0.52]*</td></tr><tr><td>-0.25<br/>[-0.71;0.21]</td><td>-0.01<br/>[-0.37;0.35]</td><td>0.21<br/>[-0.18;0.59]</td><td>0.03<br/>[-0.36;0.42]</td><td>0.61<br/>[0.32;0.9]***</td></tr><tr><td>Blood</td><td>0.12<br/>[-0.38;0.62]</td><td>0.38<br/>[0.02;0.75]</td><td>0.27<br/>[-0.09;0.62]</td><td>0.59<br/>[0.22;0.96]**</td><td>0.41<br/>[0.12;0.69]**</td></tr><tr><td rowspan="3">13-2</td><td rowspan="2">Spleen</td><td>0.22<br/>[-0.19;0.63]</td><td>0.07<br/>[-0.25;0.39]</td><td>0.26<br/>[-0.07;0.6]</td><td>0.31<br/>[-0.05;0.67]</td><td>0.19<br/>[-0.07;0.45]</td></tr><tr><td>-0.25<br/>[-0.71;0.21]</td><td>-0.01<br/>[-0.37;0.35]</td><td>0.21<br/>[-0.18;0.59]</td><td>0.03<br/>[-0.36;0.42]</td><td>0.61<br/>[0.32;0.9]***</td></tr><tr><td>Blood</td><td>0.12<br/>[-0.38;0.62]</td><td>0.38<br/>[0.02;0.75]</td><td>0.27<br/>[-0.09;0.62]</td><td>0.59<br/>[0.22;0.96]**</td><td>0.41<br/>[0.12;0.69]**</td></tr><tr><td rowspan="3">13-3</td><td rowspan="2">Spleen</td><td>0<br/>[-0.4;0.4]</td><td>0.05<br/>[-0.24;0.35]</td><td>0.15<br/>[-0.14;0.43]</td><td>0.4<br/>[0.1;0.69]*</td><td>0.49<br/>[0.27;0.72]***</td></tr><tr><td>-0.04<br/>[-0.4;0.4]</td><td>0.24<br/>[-0.24;0.35]</td><td>0.09<br/>[-0.25;0.42]</td><td>0.37<br/>[0.01;0.73]</td><td>0.43<br/>[0.17;0.69]**</td></tr><tr><td>Blood</td><td>0<br/>[-0.4;0.4]</td><td>0.02<br/>[-0.24;0.35]</td><td>0.09<br/>[-0.25;0.42]</td><td>0.37<br/>[0.01;0.73]</td><td>0.49<br/>[0.27;0.72]***</td></tr><tr><td rowspan="3">14</td><td rowspan="2">Spleen</td><td>-0.13<br/>[-0.45;0.19]</td><td>0<br/>[-0.24;0.23]</td><td>0.17<br/>[-0.07;0.42]</td><td>0.07<br/>[-0.17;0.3]</td><td>0.05<br/>[-0.14;0.23]</td></tr><tr><td>-0.04<br/>[-0.35;0.26]</td><td>-0.08<br/>[-0.32;0.17]</td><td>-0.1<br/>[-0.35;0.15]</td><td>0.03<br/>[-0.23;0.3]</td><td>-0.06<br/>[-0.26;0.14]</td></tr><tr><td>Blood</td><td>0<br/>[-0.4;0.4]</td><td>0.16<br/>[-0.24;0.56]</td><td>0.17<br/>[-0.19;0.54]</td><td>0.26<br/>[-0.13;0.66]</td><td>0.73<br/>[0.43;1.02]***</td></tr><tr><td>15</td><td>Spleen</td><td>-0.17<br/>[-0.69;0.35]</td><td>0.16<br/>[-0.24;0.56]</td><td>0.17<br/>[-0.19;0.54]</td><td>0.26<br/>[-0.13;0.66]</td><td>0.73<br/>[0.43;1.02]***</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>TRBV</th><th>Group</th><th>D3</th><th>D4</th><th>D5</th><th>D6</th><th>CM+</th></tr></thead><tbody><tr><td></td><td>Blood</td><td>-0.18<br>[-0.67;0.31]</td><td>-0.26<br>[-0.64;0.12]</td><td>0.17<br>[-0.23;0.56]</td><td>0.58<br>[0.15;1.01]*</td><td>0.61<br>[0.29;0.92]***</td></tr><tr><td rowspan="2">16</td><td>Spleen</td><td>0.09<br>[-0.24;0.42]</td><td>-0.02<br>[-0.27;0.23]</td><td>0.12<br>[-0.12;0.36]</td><td>0.05<br>[-0.19;0.3]</td><td>0.11<br>[-0.07;0.3]</td></tr><tr><td>Blood</td><td>-0.27<br>[-0.58;0.03]</td><td>-0.17<br>[-0.41;0.08]</td><td>-0.29<br>[-0.55;-0.04]</td><td>-0.07<br>[-0.32;0.19]</td><td>0.11<br>[-0.08;0.3]</td></tr><tr><td rowspan="2">17</td><td>Spleen</td><td>0.4<br>[-0.13;0.93]</td><td>0.34<br>[-0.09;0.76]</td><td>0.24<br>[-0.14;0.61]</td><td>0.19<br>[-0.19;0.57]</td><td>0.5<br>[0.19;0.81]**</td></tr><tr><td>Blood</td><td>0.2<br>[-0.29;0.68]</td><td>-0.06<br>[-0.45;0.33]</td><td>0.07<br>[-0.35;0.49]</td><td>0.07<br>[-0.34;0.48]</td><td>0.36<br>[0.04;0.68]*</td></tr><tr><td rowspan="2">19</td><td>Spleen</td><td>0.2<br>[-0.27;0.67]</td><td>0.19<br>[-0.15;0.54]</td><td>0.49<br>[0.14;0.83]*</td><td>0.64<br>[0.3;0.98]**</td><td>0.55<br>[0.29;0.81]***</td></tr><tr><td>Blood</td><td>-0.02<br>[-0.45;0.41]</td><td>0.37<br>[0.04;0.71]</td><td>0.17<br>[-0.18;0.52]</td><td>0.7<br>[0.33;1.06]**</td><td>0.56<br>[0.29;0.83]***</td></tr><tr><td rowspan="2">20</td><td>Spleen</td><td>-0.2<br>[-0.69;0.29]</td><td>0.17<br>[-0.19;0.53]</td><td>0.44<br>[0.09;0.79]*</td><td>0.63<br>[0.26;0.99]**</td><td>0.44<br>[0.16;0.73]**</td></tr><tr><td>Blood</td><td>0.38<br>[-0.16;0.91]</td><td>0.14<br>[-0.28;0.56]</td><td>-0.12<br>[-0.58;0.33]</td><td>0.1<br>[-0.37;0.56]</td><td>0.63<br>[0.28;0.98]**</td></tr><tr><td rowspan="2">23</td><td>Spleen</td><td>0.24<br>[-0.21;0.68]</td><td>0.38<br>[0.05;0.71]</td><td>0.43<br>[0.12;0.74]*</td><td>0.55<br>[0.22;0.87]**</td><td>0.61<br>[0.36;0.85]***</td></tr><tr><td>Blood</td><td>-0.02<br>[-0.5;0.45]</td><td>0.06<br>[-0.31;0.43]</td><td>0<br>[-0.38;0.38]</td><td>0.35<br>[-0.04;0.75]</td><td>0.63<br>[0.32;0.95]***</td></tr><tr><td rowspan="2">26</td><td>Spleen</td><td>0.15<br>[-0.29;0.59]</td><td>0.03<br>[-0.3;0.36]</td><td>0.75<br>[0.44;1.07]***</td><td>1.05<br>[0.72;1.38]***</td><td>0.75<br>[0.5;1]***</td></tr><tr><td>Blood</td><td>0.3<br>[-0.19;0.8]</td><td>-0.25<br>[-0.63;0.13]</td><td>0.03<br>[-0.39;0.44]</td><td>0.7<br>[0.29;1.11]*</td><td>1.24<br>[0.92;1.55]***</td></tr><tr><td rowspan="2">29</td><td>Spleen</td><td>-0.34<br>[-0.85;0.18]</td><td>-0.22<br>[-0.6;0.15]</td><td>0.09<br>[-0.28;0.45]</td><td>0.5<br>[0.12;0.87]*</td><td>0.56<br>[0.27;0.85]***</td></tr><tr><td>Blood</td><td>0.04<br>[-0.41;0.48]</td><td>-0.11<br>[-0.43;0.2]</td><td>0.17<br>[-0.15;0.5]</td><td>0.49<br>[0.12;0.85]*</td><td>0.74<br>[0.48;0.99]***</td></tr><tr><td rowspan="2">30</td><td>Spleen</td><td>0.33<br>[-0.18;0.83]</td><td>0.16<br>[-0.22;0.53]</td><td>0.37<br>[0.01;0.73]</td><td>0.71<br>[0.34;1.09]**</td><td>0.28<br>[-0.01;0.58]</td></tr><tr><td>Blood</td><td>0.21<br>[-0.3;0.71]</td><td>-0.06<br>[-0.46;0.35]</td><td>0.25<br>[-0.14;0.64]</td><td>0.42<br>[0;0.85]</td><td>0.36<br>[0;0.71]</td></tr><tr><td rowspan="2">31</td><td>Spleen</td><td>-0.09<br>[-0.55;0.37]</td><td>0.08<br>[-0.26;0.41]</td><td>0.51<br>[0.18;0.84]*</td><td>0.34<br>[-0.01;0.68]</td><td>0.41<br>[0.16;0.67]**</td></tr><tr><td>Blood</td><td>0.3<br>[-0.19;0.79]</td><td>0.02<br>[-0.37;0.4]</td><td>0.11<br>[-0.3;0.51]</td><td>0.03<br>[-0.38;0.44]</td><td>0.49<br>[0.18;0.8]**</td></tr></tbody></table>

doi:10.1371/journal.pone.0147871.t001

## Targeted T Cell Repertoire Modifications in the Brain of PbA-Infected Mice

Next, we wanted to characterize the diversity of the brain T cell population repertoire. The perturbation index was calculated as previously using CTR Spleen as the reference repertoire.

<footnote>*** adjusted p<0.001</footnote>

<footnote>** adjusted p<0.01</footnote>

<footnote>* adjusted p<0.05</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig 3.** TRBV-TRBC repertoire differences in spleen, blood and brain of naïve B10.D2 mice. (A) Mean DBV-BC perturbations across all TRBVs (µDBV-BC) in the spleen (left), blood (center) and brain (right) of CTR uninfected mice are shown for each individual. DBV-BC perturbations were computed with ISEApeaks using CTR Spleen as the reference group. Statistical comparisons were performed using pairwise t-test with correction for FDR between groups depicting significant p-value (p<0.0001) for each comparison. (B) PCA on DBV-BC log perturbation scores separating repertoires of the spleen, blood and brain on the x axis (PC1).</figcaption></figure>

doi:10.1371/journal.pone.0147871.g003

First, we focused our attention on the CTR group. Only a limited number of T cells (2 to 3.10⁴ cells) are found in the brain of naïve mice [12]. In order to characterize the whole repertoire, we pooled six brains of naïve mice per control sample to ensure enough material for full TCRβ repertoire analysis and retain the representative range of daily circulating brain T cell population [41]. We performed pairwise comparisons by t-test to confront perturbation scores of CTR brain against CTR spleen and CTR blood and found that perturbation scores of CTR groups are significantly different from each other (Figs 3A and 2B, adjusted p<0.0001 for each comparison). This result revealed a clear compartmentalization of the T cell repertoire already at steady-state and a higher inter-individual variability of blood TRBV.

We therefore analysed the global perturbation score of the brain T cell repertoire of CTR and infected mice (at day 6 p-i as well as CM⁺ mice). Globally, day 6 and CM⁺ brain repertoires are significantly different compared to CTR Spleen (p<0.0001) moreover, a significant perturbation between day 6 p-i, but not CM⁺, and uninfected brain repertoire was also detected (p = 0.0238) (Fig 4A). The high inter-individual variability of CM⁺ brain sample repertoires

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Modification of the TRBV-TRBC repertoire in the brain of B10.D2 mice during the course of PbA infection. (A) Mean DBV-BC perturbations across all TRBVs (µDBV-BC) in the brain of CTR uninfected (CTR), day 6 p-i (Day 6) and CM⁺ mice are shown for each individual. DBV-BC were computed as in Fig 1B. Black dots represent individual mouse scores. Lines represent the average score for each group. Average DBV-BC of CTR spleen is indicated as the reference (Ctr_sp). (B) Survival curve of B10.D2 mice infected with 10⁶ PbA-PRBC. All mice developed ECM symptoms. 75% of the mice died between day 6 and day 8. (C) Regression slope and 95% confidence band of parasitemia (%) over day post infection. The infection has no effect on the % of parasitemia (no significant slope). (D) DBV-BC perturbations scores of five TRBV. Statistical tests were performed as described in Fig 2.</figcaption></figure>

doi:10.1371/journal.pone.0147871.g004
might explain these inconsistent results. To ensure that this variability was not due to the large
timeframe of ECM onset as shown in Fig 4B, we subdivided the CM⁺ group into two groups:
CM⁺ mice identified between day 6 and 8 p-i and CM⁺ identified between day 10 and 14, given
that all the CM⁺ mice died after similar clinical manifestations. Global perturbation scores cal-
culated for both groups showed no significant differences (data not shown). These results are
consistent with the comparable level of parasitemia observed for mice developing ECM
between day 7 and 14 p-i (Fig 4C, not significant regression coefficient; 0 is included in the
confidence interval and p = 0.2701). This again strongly supports the hypothesis that impor-
tant TCR repertoire perturbations are a cause of ECM development. However, another non-
exclusive explanation would be that few TRBVs are perturbed in the brain during infection.
Since all TRBVs in the brain are significantly perturbed at day 6 p-i and CM⁺ compared to the
uninfected spleen group (data not shown), we therefore analysed the perturbation per TRBV in
the brain against the uninfected brain group (Fig 4D). t-tests on pooled coefficients against
uninfected mice showed that three TRBVs are significantly perturbed in infected CM⁺ mice
compared to CTR brain: TRBV1 (adjusted p = 0.0044), TRBV29 (adjusted p = 0.0237) and
TRBV3 (adjusted p = 0.0237). TRBV13-3 was significantly perturbed at day 6 p-i (adjusted
p = 0.007) but not in CM⁺ brain repertoires.

Altogether, our data show that, although there is a massive perturbation of the whole spleen
and blood repertoire of infected mice, only three TRBV are perturbed in the brain of infected
mice compared to uninfected mice.

A TCRβ Repertoire Signature Predicts Cerebral Malaria

Furthermore, we tested the impact of ECM onset on TRBV peak composition. In other words, we looked for CM⁺ signatures as a combination of individual peaks taken separately from several TRBV that could predict the CM⁺ onset. The peak database was tested in GSEA software with pre-ranked option to determine whether our *a priori* defined peak sets are significantly enriched in the CM⁺ group. Peak sets are generated separately for each organ by pvclust and assembled within a peak set database. We retained clusters having at least 95% of confidence. In the presence of missing data, peak values are imputed using the multiple imputation process (see Material and Methods). Peaks were then ranked according to the pooled t-statistic value comparing uninfected samples to CM⁺ samples. GSEA reports (Fig 5A) identified a peak set that was significantly enriched in the CM⁺ group in the spleen (NES = 2.1744, q-value = 0), in the blood (NES = 2.3300, q-value = 0) and in the brain (NES = 2.2873, q-value = 0). This peak set contains peaks from TRBV1, TRBV3, TRBV4, TRBV5, TRBV12-2, TRBV13-2, TRBV19, TRBV26 and TRBV29 (S2 Fig) suggesting a rather restricted pathogenic T cell response. In order to assess the predictive power of this signature, we applied the process described in the methodology section separately for spleen and blood. For validation, we used previous data published in [22] as a test set to predict ("AC set"). As a training set, we used our current data from spleen and blood ("EMF set"). Due to a high number of missing data on the "EMF set", we applied the multiple imputation process for peak data. Random forest models were trained on 20 CTR and 23 CM⁺ in spleen and on 20 CTR and 21 CM⁺ in blood and tested on 6 CTR and 9 CM⁺ in spleen and on 7 CTR and 13 CM⁺ in blood. Accuracy was assessed on train and test sets for each step. Random forest models showed very good accuracy > 90% of good prediction (Table 2). Altogether, the identified peak set allows to discriminate CM⁺ from CTR mice regardless of the organs, as validated on two independent datasets. We then asked whether the combinations of most discriminant peaks were different between organs, Random forest models providing score indexes to rank variables according to their importance. In the spleen, the most
A

Spleen

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

Blood

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

Brain

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

B

<figure><img src="image_11.png" /><figcaption>**Fig 5.** A unique TRBV signature discriminates CTR from CM⁺ spleen, blood and brain repertoires. (A) Panels of GSEA report for the peak set BL_PVCLUST_37 showing significant enrichment for CM⁺ condition compared to uninfected group in spleen (left) and blood (middle) and in the brain (right). The enrichment figures show ranked peaks according to the pooled t-statistic across imputed dataset (bottom). Peak positions are indicated on the ranked list (middle). The enrichment score (ES) is the maximum of the running sum (top). (B) Normalized GSEA enrichment scores post-infection growth curves in three compartments for the peak set BL_PVCLUST_37. In comparison to the CTR groups, the set is significantly enriched in day 5 post-infection in spleen (solid line), in day 6 in blood (dashed line) and in ECM state in brain (dotted line).</figcaption></figure>

doi:10.1371/journal.pone.0147871.g005

<table><caption>**Table 2.** Accuracy of random forest models.</caption><thead><tr><th colspan="4">Train on EMF, test on AC</th><th colspan="2">Train on AC, test on EMF</th></tr><tr><th>Prediction accuracy</th><th>Intermediate model Imputed train data</th><th>Final model Imputed train data</th><th>Final model test data</th><th>Final model test data</th><th>Final model train data</th></tr></thead><tbody><tr><td>Spleen</td><td>97,00%</td><td>100,00%</td><td>100,00%</td><td>93,30%</td><td>92,70%</td></tr><tr><td>Blood</td><td>94,60%</td><td>99,56%</td><td>90,00%</td><td>95,00%</td><td>90,24%</td></tr></tbody></table>

doi:10.1371/journal.pone.0147871.t002
discriminant peaks are TRBV26_008, TRBV1_009, TRBV3_010 and TRBV5_010 whereas in the blood, the most discriminant peaks are TRBV1_009, TRBV19_006, TRBV26_008, TRBV5_010.

To determine whether this signature can predict ECM outcome, we performed GSEA between each “kinetic” group and the CTR group in three compartments separately (due to small sample size, day 3 and 4 were pooled). During the course of infection, the peak set is significantly enriched at day 5 p-i in the spleen (NES = 1.9833; q-value = 0.0049), at day 6 p-i in the blood (NES = 1.9387; q-value = 0.0023) and in CM⁺ mice for the brain. Until day 6 p-i, NES is higher in spleen than in blood (Fig 5B), and then, it seems to converge to an equivalent value for the three compartments. These results are similar to those obtained from the analyses of global perturbation scores above. These findings reveal that the compartmentalized TCRβ repertoire changes can be detected at the level of TCR peaks. Moreover, the progression of NES in the spleen and in the blood underlines the diagnostic value of this TCRβ signature.

## Discussion

Cerebral malaria is a complex situation in which T cells are “necessary” to the development of the disease. Indeed, in the absence of T cells, ECM-susceptible mice are protected [8,42]. It is well described that CD8 T cells display a pathogenic role [12,13], mainly through their cytotoxic activity as shown by the protection associated with the perforin deficiency [43]. However, although up to seven target epitopes have been identified [15–17], none of them could definitely explain ECM development in mice [17]. A possible explanation is that as shown by Poh et al., part of these PbA epitopes are also shared by non-ECM inducing *Plasmodium* strains. An alternative explanation is that ECM occurs under the synergic/progressive presentation of several epitopes, including the CD8⁺ epitopes identified by the different studies but also still unknown CD4⁺ epitopes, given the critical role of CD4⁺ T cells in ECM development [8–10].

Here, we report a comprehensive analysis of the TCRβ repertoire during the course of PbA infection, including the analysis of brain T cell repertoires. Our results confirmed our previous observations of a massive perturbation of the repertoire associated with ECM onset, particularly in the blood. Furthermore, we showed that the modifications of the repertoire appear progressively in the spleen and in the blood, with TRBVs being differentially affected over time till all get perturbed. Strikingly, the comparison of brain repertoires between healthy and CM⁺ mice revealed a perturbation limited to three TRBVs out of 23. These three TRBVs are perturbed in the spleen and the blood, following a similar kinetic. These results indicate a targeted alteration of the brain TCRβ diversity, compared to the blood and the spleen. With regards to the important increase of cell numbers in the brain of CM⁺ animals, this suggests a selective recruitment of particular T cell clones in the brain directed against some particular parasitic or self-antigens.

As this is the first study showing the global repertoire of the brain T cells in naïve mice, we wanted to ensure that the CTR brain repertoire observed are not the result of blood contamination. In the case of CM⁺ brain samples, this can be excluded since most TRBV are perturbed in the blood when they are not in the brain. Surprisingly, no difference is observed between CTR blood and CTR brain repertoires. In a preliminary analysis, we observed that the blood repertoire is heterogeneous in CTR mice (data not shown). Since our perturbation index is a distance between each sample and an average reference repertoire, the DBV-BC values can differ depending on the reference repertoire used. Thus, we calculated the perturbation using CTR Blood as reference and observed that the brain T cell repertoire is indeed different compared to blood repertoire in naïve mice. This reflects the physiological compartmentalization of T cell diversity in naïve animals and highlights the importance of the microenvironment in shaping the TCR repertoire.
Finally, we could identify a TCRβ CDR3 length signature that discriminates CTR from CM⁺ repertoires in spleen, blood and brain. Using multivariate statistical modelling, we showed that this signature allows to discriminate the CTR and CM⁺ repertoires from spleen and blood obtained in a previous study, supporting the robustness and reliability of the approach. Thus, two phenomena seem to be involved in the disease. On the one hand, an alteration of the whole repertoire in spleen and blood could be associated with parasite mitogens or superantigens. Indeed, a superantigenic activity has been shown in PbA-infected C57BL/6 directing Vβ8.1 (TRBV13-3) expressing cells [20,21]. Additionally, a recent study on malaria-susceptible West-African children showed that the memory B cell repertoire is diverse, suggesting a response to several *Plasmodium* epitopes [44]. On the other hand, the identification of a unique and organ-independent TCRβ signature associated with CM⁺ onset suggests a response to specific antigenic peptides, possibly targeting antigens expressed in the cerebral compartment. The current hypothesis suggests that pathogenic T cells should be activated upon presentation of parasitic peptides [5]. Under these conditions, T cells could be selectively activated and attracted to the brain where they would play their pathogenic role, or have an autoimmune activity against brain molecules, through molecular mimicry, leading to brain damage. The recent identification of parasite-derived CD8 T cell epitopes supports this hypothesis [15-17]. Moreover, in those studies, CD8 T cells from CM⁺ mouse splenocytes also respond to the identified epitopes, again in line with our observations. Although the mouse strains used in our and the discussed studies are different, it can be noted that the TRBV signature we identified comprises the same TRBV genes as those found in response to the epitopes tested by Poh et al. and Howland et al. [15,17]. Moreover, we have shown that human CM is associated with an autoimmune blood B-cell repertoire, directed against a human brain protein [45], suggesting that *Plasmodium* can indirectly induce an autoimmunity-related process. These hypotheses are not exclusive and the activation site of pathogenic T cells remains to be elucidated. Our data suggest that the observed compartmentalization of the T cell repertoire reflects the selective migration of activated T cells from the spleen to the brain.

The complexity of the observed modifications is consistent with *Plasmodium* infection complexity. So far, therapies are oriented toward the development of vaccines targeting specific dominant antigens expressed by *Plasmodium*, and unfortunately all failed in inducing an efficient response against the parasite, as it has been confirmed in ECM by Poh et al. [17]. Finally, we identify a CM⁺ signature significantly enriched in the spleen and the blood during the course of the infection, strongly supporting the interest of diagnosis value of TCR repertoire studies. This study supports the effort made by several laboratory to analyse in parallel the modifications of the TCR repertoire diversity and the parasite antigenic variability to develop accurate and efficient therapies.

**Supporting Information**

S1 Fig. Random Forest prediction modelling approach.
(TIF)

S2 Fig. TCRβ signature peaks.
(TIF)

**Acknowledgments**

This work has been funded by Assistance Publique-Hôpitaux de Paris, Université Pierre and Marie Curie (Paris VI), LabEx Transimmunom (ANR-11-IDEX-0004-02), LabEx Parafrap (ANR-11-LABX-0024) and LIA CNRS DBT SIGID ("Systems Immunology and genetics of
infectious diseases"). EMF & HPP have been supported by a doctoral fellowship from Ministère
de la Recherche et de la Technologie, as well as by Fondation de France and Institut Pasteur for
EMF.

Author Contributions

Conceived and designed the experiments: AS SP PAC EMF. Performed the experiments: EMF
SD OG. Analyzed the data: EMF HPP AS. Contributed reagents/materials/analysis tools: SA SP
PAC DK. Wrote the paper: EMF HPP SA SP PAC DK.

References

1. World Health Organization. World malaria report 2013. [S.I.]; 2014.

2. MacPherson GG, Warrell MJ, White NJ, Looareesuwan S, Warrell DA. A quantitative ultrastructural analysis of parasitized erythrocyte sequestration. Am J Pathol. 1985; 119: 385-401. PMID: 3893148

3. Berendt AR, Turner GDH, Newbold CI. Cerebral malaria: the sequestration hypothesis. Parasitol Today. 1994; 10: 412-412. PMID: 15275553

4. Grau GE, Mackenzie CD, Carr RA, Redard M, Pizzolato G, Allasia C, et al. Platelet accumulation in brain microvessels in fatal pediatric cerebral malaria. J Infect Dis. 2003; 187: 461-466. PMID: 12552430

5. Renia L, Potter SM, Mauduit M, Rosa DS, Kayibanda M, Deschemin J-C, et al. Pathogenic T cells in cerebral malaria. Int J Parasitol. 2006; 36: 547-554. PMID: 16600241

6. Curfs JH, Schetters TP, Hermsen CC, Jerusalem CR, Van Zon AA, Eling WM. Immunological aspects of cerebral lesions in murine malaria. Clin Exp Immunol. 1989; 75: 136. PMID: 2649283

7. Finley RW, Mackey LJ, Lambert PH. Virulent P. berghei malaria: prolonged survival and decreased cerebral pathology in cell-dependent nude mice. J Immunol. 1982; 129: 2213-2218. PMID: 6749988

8. Grau GE, Piguet PF, Engers HD, Louis JA, Vassalli P, Lambert PH. L3T4+ T lymphocytes play a major role in the pathogenesis of murine cerebral malaria. J Immunol. 1986; 137: 2348-2354. PMID: 3093572

9. Hermsen CC, Van De Wiel T, Mommers E, Sauerwein RW, Eling W. Depletion of CD4+ or CD8+ T-cells prevents Plasmodium berghei induced cerebral malaria in end-stage disease. Parasitology. 1997; 114: 7-12. PMID: 9011069

10. Yanez DM, Manning DD, Cooley AJ, Weidanz WP, Van Der Heyde HC. Participation of lymphocyte subpopulations in the pathogenesis of experimental murine cerebral malaria. J Immunol. 1996; 157: 1620-1624. PMID: 8759747

11. Boubou MI, Collette A, Voegtlé D, Mazier D, Cazenave PA, Pied S. T cell response in malaria pathogenesis: selective increase in T cells carrying the TCR V(beta)8 during experimental cerebral malaria. Int Immunol. 1999; 11: 1553-62. PMID: 10464176

12. Belnoue E, Kayibanda M, Vigario AM, Deschemin J-C, Rooijen NV, Viguier M, et al. On the Pathogenic Role of Brain-Sequestered {alpha}{beta} CD8+ T Cells in Experimental Cerebral Malaria. J Immunol. 2002; 169: 6369-6375. PMID: 12444144

13. Bagot S, Nogueira F, Collette A, do Rosario V, Lemonier F, Cazenave P-A, et al. Comparative Study of Brain CD8+ T Cells Induced by Sporozoites and Those Induced by Blood-Stage Plasmodium berghei ANKA Involved in the Development of Cerebral Malaria. Infect Immun. 2004; 72: 2817-2826. doi: 10.1128/IAI.72.5.2817-2826.2004 PMID: 15102792

14. Miyakoda M, Kimura D, Yuda M, Chinzei Y, Shibata Y, Honma K, et al. Malaria-Specific and Nonspecific Activation of CD8+ T Cells during Blood Stage of J Immunol. 2008.

15. Howland SW, Poh CM, Gun SY, Claser C, Malleret B, Shastri N, et al. Brain microvessel cross-presentation is a hallmark of experimental cerebral malaria: Antigen-specific CD8+ T cells in cerebral malaria. EMBO Mol Med. 2013; 5: 984-999. doi: 10.1002/emmm.201202273

16. Lau LS, Fernandez-Ruiz D, Mollard V, Sturm A, Neller MA, Cozijnsen A, et al. CD8+ T Cells from a Novel T Cell Receptor Transgenic Mouse Induce Liver-Stage Immunity That Can Be Boosted by Blood-Stage Infection in Rodent Malaria. Mota MM, editor. PLoS Pathog. 2014; 10: e1004135. doi: 10.1371/journal.ppat.1004135 PMID: 24854165

17. Poh CM, Howland SW, Grotenbreg GM, Renia L. Damage to the Blood-Brain Barrier during Experimental Cerebral Malaria Results from Synergistic Effects of CD8+ T Cells with Different Specificities. Infect Immun. 2014; 82: 4854-4864. doi: 10.1128/IAI.02180-14 PMID: 25156726
18. Ballet JJ, Druilhe P, Querleux MA, Schmitt C, Agrapart M. Parasite-derived mitogenic activity for human T cells in Plasmodium falciparum continuous cultures. Infect Immun. 1981; 33: 758–762. PMID: 7026445

19. Gabrielsen AA, Jensen JB. Mitogenic activity of extracts from continuous cultures of Plasmodium falciparum. Am J Trop Med Hyg. 1982; 31: 441–448. PMID: 6211106

20. Pied S, Voegtlé D, Marussig M, Renia L, Miltgen F, Mazier D, et al. Evidence for superantigenic activity during murine malaria infection. Int Immunol. 1997; 9: 17–25. PMID: 9043944

21. Gorgette O, Existe A, Boubou MI, Bagot S, Guenet J-L, Mazier D, et al. Deletion of T Cells Bearing the V8.1 T-Cell Receptor following Mouse Mammary Tumor Virus 7 Integration Confers Resistance to Murine Cerebral Malaria. Infect Immun. 2002; 70: 3701–3706. doi: 10.1128/IAI.70.7.3701–3706.2002 PMID: 12065512

22. Collette A, Bagot S, Ferrandiz ME, Cazenave P-A, Six A, Pied S. A profound alteration of blood TCRB repertoire allows prediction of cerebral malaria. J Immunol. 2004; 173: 4568–4575. PMID: 15383590

23. Davis MM, Bjorkman PJ. T-cell antigen receptor genes and T-cell recognition. Nature. 1988; 334: 395–402. PMID: 3043226

24. Pannetier C, Cochet M, Darche S, Casrouge A, Zöller M, Kourilsky P. The sizes of the CDR3 hypervariable regions of the murine T-cell receptor beta chains vary as a function of the recombined germ-line segments. Proc Natl Acad Sci U S A. 1993; 90: 4319–4323. PMID: 8483950

25. Collette A, Six A. ISEApeaks: an Excel platform for GeneScan and Immunoscope data retrieval, management and analysis. Bioinforma Oxf Engl. 2002; 18: 329–30.

26. Collette A, Cazenave P-A, Pied S, Six A. New methods and software tools for high throughput CDR3 spectratyping. Application to T lymphocyte repertoire modifications during experimental malaria. J Immunol Methods. 2003; 278: 105–116. doi: 10.1016/S0022-1759(03)00225-4 PMID: 12957400

27. Chou HS, Anderson SJ, Louie MC, Godambe SA, Pozzi MR, Behlke MA, et al. Tandem linkage and unusual RNA splicing of the T-cell receptor beta-chain variable-region genes. Proc Natl Acad Sci. 1987; 84: 1992–1996. PMID: 3470773

28. Wade T, Bill J, Marrack P, Palmer E, Kappler JW. Molecular basis for the nonexpression of V{beta}17 in some strains of mice. J Immunol. 1988; 141: 2165–2167. PMID: 2844897

29. Louie MC, Nelson CA, Loh DY. Identification and characterization of new murine T cell receptor beta chain variable region (V beta) genes. J Exp Med. 1989; 170: 1987–1998. PMID: 2531193

30. Lefranc M-P. WHO-IUIS Nomenclature Subcommittee for immunoglobulins and T cell receptors report. Immunogenetics. 2007; 59: 899–902. doi: 10.1007/s00251-007-0260-4 PMID: 18046549

31. Gorochov G, Neumann AU, Kereveur A, Parizot C, Li T, Katlama C, et al. Perturbation of CD4+ and CD8+ T-cell repertoires during progression to AIDS and regulation of the CD4+ repertoire during antiviral therapy. Nat Med. 1998; 4: 215–221. PMID: 9461196

32. Han M, Harrison L, Kehn P, Stevenson K, Currier JR, Robinson MA. Invariant or highly conserved TCR alpha are expressed on double-negative (CD31CD42CD82) and CD81 T Cells. J Immunol. 1999; 163: 301–311. PMID: 10384129

33. Rubin DB. Multiple Imputation for Nonresponse in Surveys. J. Wiley & Sons. New-York; 1987.

34. Heitjan DF, Little RJA. Multiple imputation for the Fatal Accident Reporting System. Appl Stat. 1991; 40: 1329.

35. Schenker N, Taylor JMG. Partially parametric techniques for multiple imputation. Comput Stat Data Anal. 1996; 22: 425–446.

36. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995; 57: 289–300.

37. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102: 15545–15550. PMID: 16199517

38. Shimodaira H. Approximately unbiased tests of regions using multistep-multiscale bootstrap resampling. Ann Stat. 2004; 32: 2616–2641.

39. Van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate imputation by chained equations in R. J Stat Softw. 2011; 45: 1–67.

40. Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002; 2: 18–22.

41. Engelhardt B, Ransohoff RM. The ins and outs of T-lymphocyte trafficking to the CNS: anatomical sites and molecular mechanisms. Trends Immunol. 2005; 26: 485–495. PMID: 16039904

42. Grau GE, Gretener D, Lambert PH. Prevention of murine cerebral malaria by low-dose cyclosporin A. Immunology. 1987; 61: 521. PMID: 3327806
43. Nitcheu J, Bonduelle O, Combadiere C, Tefit M, Seilhean D, Mazier D, et al. Perforin-dependent brain-infiltrating cytotoxic CD8+ T lymphocytes mediate experimental cerebral malaria pathogenesis. J Immunol. 2003; 170: 2221–2228. PMID: 12574396

44. Zinocker S, Schindler CE, Skinner J, Rogosch T, Waisberg M, Schickel J-N, et al. The V Gene Repertoires of Classical and Atypical Memory B Cells in Malaria-Susceptible West African Children. J Immunol. 2015; 194: 929–939. doi: 10.4049/jimmunol.1402168 PMID: 25556245

45. Guiyedi V, Chanseaud Y, Fesel C, Snounou G, Rousselle J-C, Lim P, et al. Self-Reactivities to the Non-Erythroid Alpha Spectrin Correlate with Cerebral Malaria in Gabonese Children. Rodrigues M, editor. PLoS ONE. 2007; 2: e389. doi: 10.1371/journal.pone.0000389 PMID: 17460756