# Platelets Alter Gene Expression Profile in Human Brain Endothelial Cells in an In Vitro Model of Cerebral Malaria

Mathieu Barbier, Dorothée Faille, Béatrice B. Loriod, Julien Textoris, Claire Camus, Denis Puthier, Laurence Flori, Samuel Crocodile Wassmer, Genevieve Victorero, Marie-Christine Alessi, et al.

► To cite this version:

Mathieu Barbier, Dorothée Faille, Béatrice B. Loriod, Julien Textoris, Claire Camus, et al. Platelets Alter Gene Expression Profile in Human Brain Endothelial Cells in an In Vitro Model of Cerebral Malaria. PLoS ONE, 2011, 6 (5), pp.e19651. 10.1371/journal.pone.0019651. hal-01743259

HAL Id: hal-01743259

https://amu.hal.science/hal-01743259v1

Submitted on 26 Mar 2018

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 CC BY 4.0 - Attribution - International License
# Platelets Alter Gene Expression Profile in Human Brain Endothelial Cells in an *In Vitro* Model of Cerebral Malaria

Mathieu Barbier¹, Dorothée Faille³,⁴,⁶, Béatrice Loriod², Julien Textoris², Claire Camus¹,², Denis Puthier², Laurence Flori⁵, Samuel Crocodile Wassmer⁷, Geneviève Victorero², Marie-Christine Alessi⁴, Thierry Fusai³, Catherine Nguyen², Georges E. Grau⁶, Pascal Rihet¹,²*

¹Laboratoire de Pharmacogenétique des Maladies Parasitaires-EA 864, IFR 48, Faculté de Pharmacie, Aix-Marseille Université, Marseille, France, ²INSERM U928-TAGC, Aix-Marseille Université, IFR137, Marseille, France, ³Unité de Recherche en Biologie et en Épidémiologie Parasitaires, UMR 6236-URMITE-IMTSSA, Institut de Recherche Biomédicale des Armées-Antenne Marseille, Marseille, France, ⁴INSERM, UMR 626, Faculté de Médecine, Aix-Marseille Université, Marseille, France, ⁵Laboratoire de Génétique Animale et Biologie Intégrative, Domaine de Vilvert, INRA AgroParisTech, Jouy-en-Josas, France, ⁶Vascular Immunology Unit, Department of Pathology, University of Sydney, Camperdown, Australia, ⁷Department of Medical Parasitology, New York University School of Medicine, New York, New York, United States of America

## Abstract

Platelet adhesion to the brain microvasculature has been associated with cerebral malaria (CM) in humans, suggesting that platelets play a role in the pathogenesis of this syndrome. *In vitro* co-cultures have shown that platelets can act as a bridge between *Plasmodium falciparum*-infected red blood cells (pRBC) and human brain microvascular endothelial cells (HBEC) and potentiate HBEC apoptosis. Using cDNA microarray technology, we analyzed transcriptional changes of HBEC in response to platelets in the presence or the absence of tumor necrosis factor (TNF) and pRBC, which have been reported to alter gene expression in endothelial cells. Using a rigorous statistical approach with multiple test corrections, we showed a significant effect of platelets on gene expression in HBEC. We also detected a strong effect of TNF, whereas there was no transcriptional change induced specifically by pRBC. Nevertheless, a global ANOVA and a two-way ANOVA suggested that pRBC acted in interaction with platelets and TNF to alter gene expression in HBEC. The expression of selected genes was validated by RT-qPCR. The analysis of gene functional annotation indicated that platelets induce the expression of genes involved in inflammation and apoptosis, such as genes involved in chemokine-, TREM1-, cytokine-, IL10-, TGFβ-, death-receptor-, and apoptosis-signaling. Overall, our results support the hypothesis that platelets play a pathogenic role in CM.

**Citation:** Barbier M, Faille D, Loriod B, Textoris J, Camus C, et al. (2011) Platelets Alter Gene Expression Profile in Human Brain Endothelial Cells in an *In Vitro* Model of Cerebral Malaria. PLoS ONE 6(5): e19651. doi:10.1371/journal.pone.0019651

**Editor:** Rick Edward Paul, Institut Pasteur, France

**Received** December 21, 2010; **Accepted** April 12, 2011; **Published** May 16, 2011

**Copyright:** © 2011 Barbier 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.

**Funding:** MB was supported by a studentship from the French Ministry of Research and Technology. The study was supported by INSERM, by Aix-Marseille University, by NHMRC, by Australian Research Council and by the AL Kerr Bequest, Sydney Medical School. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

## Introduction

Malaria remains a major problem of public health in tropical developing countries worldwide. Close to 230 million people (World Health Organization) are infected annually with *Plasmodium falciparum*. About 1–2% out of clinical attacks of malaria are complicated by severe manifestations leading to death. The major syndromes of severe malaria are severe anaemia, severe respiratory distress, and cerebral malaria (CM) [1]. Severe malaria predominantly affects children under 5 years of age, non-immune adults, and pregnant women. More generally, the outcome of infection depends on numerous factors, including host and parasite genetics, host age, rates of infection, and host immune status.

CM is a complex neurological syndrome, which accounts for a large proportion of deaths due to malaria. CM is characterized by the sequestration of *P. falciparum*-infected red blood cells (pRBC), leukocytes, and platelets in brain microvessels. It has been proposed that this accumulation of pRBC interferes with microcirculation, and leads to ischemia [2]. In addition, many

observations provide evidence that leukocytes and inflammatory cytokines, such as tumor necrosis factor (TNF), are involved in malaria pathogenesis [2,3]. TNF produced by monocytes activates endothelial cells, thereby increasing the expression of cell adhesion molecules, such as P-selectin, ICAM-1 and VCAM-1. Increased expression of endothelial cell adhesion molecules contributes to enhanced pRBC sequestration, leukocyte and platelets adhesion.

There is a growing body of evidence that platelets (PL) are involved in CM pathogenesis: (i) thrombocytopenia is marked during *P. falciparum* and experimental CM; (ii) PL adhesion to the microvessels has been associated with CM in humans [4], and has been shown to be involved in mouse CM pathogenesis [4,5]; (iii) *in vitro* co-cultures showed that PL act as a bridge between pRBC and TNF-stimulated endothelium, and contribute to brain endothelial alterations [6,7].

These studies prompted us to analyze gene expression in endothelial cells co-incubated with PL in an *in vitro* model of CM. Microarray analyses showed that pro-inflammatory cytokines, including TNF, induce distinct gene expression programs in

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

<footnote>* E-mail: rihet@tagc.univ-mrs.fr</footnote>

<footnote>‡ Current address: UMR 745, Faculté des Sciences Pharmaceutiques et Biologiques, Université Paris Descartes, Paris, France</footnote>
different endothelial cell types [8,9]; this indicates that the endothelial cell type should be carefully chosen for an *in vitro* model of CM. Some studies have investigated the direct effect of pRBC alone on the human umbilical vein endothelial cell (HUVEC) transcriptome, and revealed changes in the expression of genes involved in tight junction, inflammatory response, and apoptosis [10,11,12], whereas co-cultures with human lung microvascular endothelial cells (HLEC) showed that pRBC induces apoptosis through a caspase-dependent pathway [13]. Recently, a genome wide transcriptional profiling of human brain microvascular endothelial cells (HBEC) after exposure with pRBC revealed up-regulated transcripts involved in the immune response, apoptosis, and NF-κB activation cascade [14]. Nonetheless, the influence of PL on endothelial cell gene expression has not been evaluated in those models. We therefore investigated transcriptional changes in HBEC in a tri-partite co-culture model including TNF, pRBC, and PL. This model is a well established *in vitro* model of CM, which has been validated for a number of interactions, as recently reviewed [15]. Here, we report that PL induced transcriptional changes in HBEC related to canonical pathways involved in inflammation or apoptosis.

# Materials and Methods

## Human brain endothelial cells

HBEC-5i were derived by Dorovonis-Zis and colleagues from small fragments of human cortex obtained from patients who died of various causes, devoid of any pathological abnormalities [16]. These cells were immortalized and characterized as described elsewhere [7]. HBEC-5i were grown to confluence in DMEM/F12 medium containing 15 mmol/L HEPES and L-glutamine pH 7.4 (Gibco Life Technologies, Carlsbad USA) supplemented with 10% (v/v) of fetal bovine serum. Conservation of the phenotype of this cell line was checked over passages and cell line was tested negative for mycoplasma contamination every 5 passages.

## Platelets

Venous blood was obtained by venipuncture from one healthy volunteer and anticoagulated with 3.2 g/L buffered sodium citrate. The volunteer had not taken any drugs for at least 1 month. The healthy volunteer gave a written informed consent. Experiments were performed after clearance from human ethics committee of the University of Sydney, which approved the study. Platelet-rich plasma was prepared by centrifugation at 200 g for 15 minutes at room temperature then PL were pelleted by centrifugation of the PL-rich plasma for 6 minutes at 2000 g and washed in HEPES buffer warmed at 37°C (0.137 M NaCl, 2.68 mM KCl, 1 mM MgCl₂, 1 mM CaCl₂, 5 mM HEPES, 0.1% glucose, pH = 6.8). PL were then adjusted to 1.10⁸/mL. PL functions were tested by agonist-induced platelet aggregation and PL activation was assessed by flow cytometric measurement of P-selectin in PL-rich plasma and washed PL suspension. PL of the volunteer had a normal reactivity compared to a reference population of 35 healthy controls, and washing steps did not alter either platelet reactivity or activation state.

## Parasites

The IPPAM strain of *P. falciparum* (gift from C. Behr, Institut Pasteur, Paris, France) was obtained by panning on endothelial cells expressing either GPIV (CD36) or ICAM-1 using the method described by Fried and Duffy [17]. Parasites were maintained in continuous culture at 2% hematocrit using type O+ human red blood cells (RBC) as described elsewhere [18]. Uninfected normal red blood cells (NRBC) used as controls were cultured the same

way for at least 2 weeks before experiments. Late trophozoite-stage pRBC preparations were enriched to 80 to 85% by gelatin flotation with Plasmion (Fresenius Kabi France, Couvier France) as described elsewhere [19]. NRBC or pRBC were adjusted to 2×10⁷/mL in the co-culture medium composed of RPMI 1640 with the pH adjusted to 6.8 before use.

## Human brain endothelial cells-platelets-parasitized red blood cells co-cultures

HBEC were grown to confluence in 6-well tissue culture plates and were either left unstimulated or incubated overnight with 10 ng/mL recombinant TNF (PeproTech, London UK). HBEC were then washed with PBS and incubated with PL, pRBC or NRBC according to the experimental conditions with the following sequence: HBEC were first incubated with 2×10⁸ PL per well as a first step, washed three times in PBS to remove unbound PL and 4×10⁷ pRBC or NRBC per well were then added as a second step. First and second incubation steps were carried out at 37°C for 90 min. Co-cultures without TNF, PL, pRBC, or NRBC were also performed. Table 1 summarizes all the experimental conditions (n = 24). Each experimental condition was done in triplicate.

## RNA isolation

After co-incubation, HBEC were washed three times with PBS and either harvested immediately (t=0 h) or further incubated for 5 h (t=5 h) and then harvested. Total RNA extraction of HBEC was carried out after harvesting using TRIzol (Gibco Life Technologies, Carlsbad USA) according to the manufacturer's instructions. RNA concentration was determined by reading the absorbance at 260/280 nm. The quality of RNA and the absence of any DNA contamination were checked on a 2100 Bioanalyzer (Agilent Technologies, Santa Clara USA). Each mRNA sample extracted from an individual sample was run on a single microarray.

## cDNA Microarray Hybridization

Hybridizations were carried out for each triplicate of the 24 conditions. All microarray procedures were done at our microarray facilities. Human cDNA microarray were designed and prepared as previously described [20]. PCR amplification was performed as previously described [21], and PCR products were spotted onto nylon membranes (Hybond-N+, Amersham, Pharmacia Biotech, UK) with a GMS-427 arrayer (Affymetrix, Santa Clara USA). 9216 clones were spotted; they corresponded to 8780 genes and 436 controls. Production, sensitivity and reproducibility

<table><caption>Table 1. Summary of the 24 experimental conditionsª.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="5">Platelets (PL)</th></tr><tr><th>-</th><th>-</th><th>TNF (T)</th><th>-</th><th>TNF (T)</th></tr></thead><tbody><tr><td rowspan="3">0 h<sup>b</sup></td><td>-</td><td>0 h</td><td>0h_T</td><td>0h_PL</td><td>0h_T_PL</td><td></td></tr><tr><td>NRBC (NR)</td><td>0hNR</td><td>0h_TNR</td><td>0hNR_PL</td><td>0hTNR_PL</td><td></td></tr><tr><td>pRBC (PR)</td><td>0h_PR</td><td>0h_T_PR</td><td>0h_PR_PL</td><td>0hT_PR_PL</td><td></td></tr><tr><td rowspan="3">5 h<sup>b</sup></td><td>-</td><td>5 h</td><td>5h_T</td><td>5h_PL</td><td>5hT_PL</td><td></td></tr><tr><td>NRBC (NR)</td><td>5hNR</td><td>5h_TNR</td><td>5hNR_PL</td><td>5hTNR_PL</td><td></td></tr><tr><td>pRBC (PR)</td><td>5h_PR</td><td>5h_T_PR</td><td>5h_PR_PL</td><td>5hT_PR_PL</td><td></td></tr></tbody></table>

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

<footnote>ªLabels used on Figure 1. Each experimental condition was done in triplicate.</footnote>

<footnote>ᵇTime after co-incubation.</footnote>

<footnote>doi:10.1371/journal.pone.0019651.t001</footnote>
A

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

B

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

<figure><img src="image_6.png" /></figure>
**Figure 1.** Hierarchical clustering of HBEC responses to NRBC, pRBC, TNF, and platelets. A) Hierarchical clustering of 72 samples was done on the basis of the expression of 469 genes that were selected by using a one-way ANOVA. Each row represents a gene, and each column an experimental sample. Red and green indicate expression levels above and below the median for each gene, respectively. Similarities in gene expression are represented by dendogrammes of samples and genes. B) Dendogramme of samples. 24 experimental conditions and 3 samples per condition are represented. For each condition, the presence (black) or the absence (white) of factors is specified. The factors are NRBC (−/+), pRBC (−/+), co-culture time (0 h and 5 h), PL (−/+), and TNF (−/+). doi:10.1371/journal.pone.0019651.g001

of the microarrays have been previously described [22]. Microarray were hybridized with ³³P-labelled probes, first with an oligonucleotide sequence common to all spotted PCR products (5'TCACACAGGAAACAGCTATGAC3') called vector probe to measure the amount of PCR products spotted onto the microarray. Then, after stripping, microarrays were hybridized with complex probes made from 5 µg of retrotranscribed total RNA. Probe preparations, hybridizations and washes were carried out as described previously [22]. Microarrays were exposed on a radioimager during 48 h and scanned using a Fuji BAS5000 machine (Fuji, France). Hybridization signals were quantified using BZscan 2 software [23] locally developed on the TAGC platforms (http://tagc.univ-mrs.fr/welcome/). All data are MIAME compliant and have been loaded into ArrayExpress database (http://www.ebi.ac.uk/microarray-as/ae/). The ArrayExpress accession number of this experiment is E-MEXP-1989.

## Microarray data Analysis

BZscan first applied a criterion of spot quality to select analyzable spots on the microarrays. 4502 genes passed through this selection for all conditions and were analyzed. Since the hybridization signal depends on the amount of PCR products spotted onto the microarray, the signal obtained with the vector probe was used to normalize the signal obtained with complex probes [24]. A specific R library, 'nylonmagic' (currently under development, script available upon request), that uses the 'S4' system of formal classes and methods was used to process and normalize nylon microarray

data [25]. After loading, data obtained with vector probe and complex probes were both corrected for neighborhood effects and local background as described [23]. Quantile normalization was applied to vector probe data (V) and complex probe data (C) to correct for global intensity and dispersion. In a third step, the correction by the vector signal was made by calculating a C/V ratio before log transformation (base 2). The basic approach relies on the assumption that the complex and vector signal intensities are related by a constant factor. However, as reported by others, in many cases, log(C/V) still depends on log(V) [24]. To correct for this bias, we used the method proposed by Talby *et al* [24]. For any sample, $C_k$ and $V_k$ denote the complex and vector probe signals, and k represents the array element. The curve $M_k = \log(C_k/V_k)$ versus $\log(V_k)$ was fitted using a robust scatter-plot smoother (lowess regression). The normalization adjustment was $M'_k = M_k - s(\log(V_k))$, where $s(\log(V_k))$ was the lowess fit to the M versus $\log(V)$. In the last step, data were centered relative to the median for each gene.

Unsupervised hierarchical clustering was used to investigate relationships between samples and relationships between genes, and to evaluate the effect of the different sources of variability (TNF, PL, RBC, time-points, and replicates). It was applied to median-centred data, using the Cluster and TreeView programs (average linkage clustering using Euclidian distance as similarity metric).

Statistical analysis was performed using the TIGR MeV (MultiExperiment Viewer) v4.1 software (http://www.tm4.org/mev.html) and the GeneANOVA programs [26]. One-way analysis of variance was applied to identify genes differentially expressed between the treatment groups. Paired t-tests and paired significant microarray analysis method (SAM) were used to evaluate the influence of one particular variability factor (TNF, PL, RBC, and time). The members of each pair differed from each other only with respect to the factor studied. Therefore, the paired design allowed isolating the contribution of the factor of interest. Empirical *P*-values were calculated after 10000 permutations. Moreover, we applied the adaptive false discovery rate (FDR) procedure to control the FDR for the paired t-tests and the paired SAM, respectively [27,28]. The FDR of a test is the expected proportion of false positive results among the declared significant ones. We generally used an FDR of 5%.

In addition, an ANOVA based on the microarray analysis was used to analyze simultaneously the effect of TNF, PL, RBC, and time on gene expression. A global ANOVA model gives an estimation of the contribution of each factor in the total variation of the whole data set; gene was considered a factor in this model, and interaction terms were included [29]; no multiple test correction was performed. Besides, a local ANOVA allows the determination of each contribution for each gene: it gives an estimation of the variation due to the factor studied, and the significance of the estimate. We applied the adaptive false discovery rate (FDR) procedure to control the FDR for the local ANOVA [27,28]. We generally used an FDR of 5%. In addition, a two-way ANOVA was performed to analyse for each gene the interaction RBC-by-TNF and the interaction RBC-by-PL. Empirical *P* values were calculated after permutations, and multiple testing was controlled by using the adaptive FDR.

<table><caption>**Table 2.** Multi-way analysis of variance for HBEC gene expression data.</caption><thead><tr><th>Factor</th><th>Sum of square</th><th>DF<sup>a</sup></th><th>Mean square</th><th>F</th><th>P-value</th></tr></thead><tbody><tr><td>gene</td><td>849219.05</td><td>4501</td><td>188.67</td><td>274.11</td><td>&lt;0.00001</td></tr><tr><td>RBC</td><td>0.37</td><td>2</td><td>0.18</td><td>0.27</td><td>-</td></tr><tr><td>TNF</td><td>0.17</td><td>1</td><td>0.17</td><td>0.24</td><td>-</td></tr><tr><td>PL</td><td>4.37</td><td>1</td><td>4.37</td><td>6.56</td><td>0.01013</td></tr><tr><td>time</td><td>0.1</td><td>1</td><td>0.1</td><td>0.15</td><td>-</td></tr><tr><td>gene × RBC</td><td>6341.37</td><td>9002</td><td>0.7</td><td>1.02</td><td>0.0416</td></tr><tr><td>gene × TNF</td><td>6497.82</td><td>4501</td><td>1.44</td><td>2.1</td><td>&lt;0.00001</td></tr><tr><td>gene × PL</td><td>4644.08</td><td>4501</td><td>1.03</td><td>1.5</td><td>&lt;0.00001</td></tr><tr><td>gene × time</td><td>5355.23</td><td>4501</td><td>1.19</td><td>1.73</td><td>&lt;0.00001</td></tr><tr><td>gene × time × TNF</td><td>3250.75</td><td>4501</td><td>0.72</td><td>1.09</td><td>0.00005</td></tr><tr><td>gene × time × PL</td><td>3264.17</td><td>4501</td><td>0.73</td><td>1.09</td><td>0.00002</td></tr><tr><td>gene × TNF × PL</td><td>3219.91</td><td>4501</td><td>0.72</td><td>1.07</td><td>0.00029</td></tr><tr><td>gene × time × RBC × TNF</td><td>6177.54</td><td>9002</td><td>0.69</td><td>1.03</td><td>0.0214</td></tr><tr><td>gene × time × RBC × PL</td><td>6256.12</td><td>9002</td><td>0.69</td><td>1.04</td><td>0.0026</td></tr><tr><td>Residual<sup>b</sup></td><td>176809.89</td><td>265625</td><td>0.67</td><td>-</td><td>-</td></tr><tr><td>Total</td><td>1070570.59</td><td>324143</td><td>3.3</td><td>-</td><td>-</td></tr></tbody></table>

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

<footnote><sup>a</sup>Degree of freedom.</footnote>

<footnote><sup>b</sup>The residual corresponds to the variance that was not explained by the model. doi:10.1371/journal.pone.0019651.t002</footnote>
A

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

B

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

<figure><img src="image_10.png" /></figure>
**Figure 2. Gene expression variance analysis.** Each gene was plotted as a point. The abscissa is the variation due to the factor normalized by the total variation of the gene, and the ordinate is the logarithm of the *P*-value. The nominal *P*-value that corresponded to a FDR of 5% is stated. The results are shown for the TNF (A) and PL (B) factors.
doi:10.1371/journal.pone.0019651.g002

Gene annotation of all 4502 analyzable genes/expressed sequence tags (ESTs) was performed using DAVID database (http://david.abcc.ncifcrf.gov/) [30]. This program, which allows a biological interpretation of gene clusters on the basis of gene ontology (GO) terms and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways (http://www.genome.jp/kegg/), were used to assess whether specific biological pathways were overrepresented among the differentially expressed genes and within specific gene clusters. A score based on Fisher's exact test reflected the probability that the prevalence of a particular term within a cluster was due to chance alone, given the prevalence of that term in the population of all genes under study. The FDR was used to account for multiple testing [27].

The Ingenuity Pathways Analysis (IPA) program permitted the determination of significant networks, top functions and canonical pathways associated with the differentially expressed genes (IPA 5.0, Ingenuity Systems Inc., USA). The IPA program searches the Ingenuity Pathway Knowledge Base (IPKB) for interactions (known from the literature) between the studied gene set and all other genes contained in IPKB and generates a series of networks. The genes selected by the IPA network analysis are called focus genes. Fisher's exact test was used to assign statistical significance, and each network's score was displayed as the −log (*P*-value). Significant networks with a score greater than 3 were selected (*P*<0.001) and top functions associated with these networks are presented. Canonical pathway analysis identified the most significant pathways, from the IPA library of canonical pathways. The association between the input data set and the canonical pathway was scored as the ratio of the number of genes from the data set that map to the pathway divided by the total number of genes that map to the canonical pathway; their significance was assessed by using a *P*-value calculated by using Fisher's exact test.

## Quantitative Real-Time PCR

RNA of HBEC used for Real-Time-quantitative Polymerase Chain Reaction (RT-qPCR) was obtained from another co-culture experiment. cDNAs were prepared from 2 µg of total RNA using random hexamers pd(N)₆ (GE Healthcare, UK) and M-MLV reverse transcriptase (Invitrogen Life Technologies, Carlsbad USA). RT-qPCR assays were performed on AB-7300 Real-Time PCR system using SYBR PCR Master Mix (Applied Biosystems Life Technologies, Carlsbad USA). Forward and reverse primers used are shown in Table S1. They were designed manually to have primers spanning exon-exon junctions to avoid the amplification of DNA contaminant. PCR were performed using the following conditions: 10 min at 95°C, 15 sec at 95°C and 1 min at 60°C during 40 cycles. Melting curve analysis was performed to check the specificity of the PCR and to verify the amplification efficiency. Relative quantification was calculated with the 2⁻⁰⁻⁰ΔΔC(t) method: ΔC(t) represents the difference between the threshold cycle of the gene of interest and the threshold cycle of the internal control gene (*ACTB*). Co-cultures were incubated during 5 h, and each tested condition was done twice.

## Results

We used an unsupervised hierarchical clustering method that grouped genes on the vertical axis and samples on the horizontal axis, on the basis of similarity in their expression profiles. The

similarities are summarized in a dendogramme in which the pattern and length of the branches reflect the relatedness of the samples. The discrimination between experimental conditions was unsatisfactory on this basis (data not shown). We assumed that the samples were incorrectly classified because of confounding genes, which were not differentially expressed between groups. Therefore, we filtered the genes by using a one-way ANOVA to visualize the effect of experimental conditions. We generated empirical *P*-value for each gene, and selected the genes for which *P*<0.05. We defined 469 informative genes, and we performed unsupervised hierarchical clustering of samples on the basis of the expression of selected genes. Figures 1A and 1B show the classification of samples obtained.

We further performed a multi-way ANOVA over the whole data set to estimate the contribution of each factor. Such estimates could be performed because we used a fully crossed factorial design, meaning that each state of one factor was found in combination with each state of other factors (Table 1). So, there were 72 measurements for each of the 4502 genes, three per combination of four factors (RBC, TNF, PL and time). Table 2 presents the results of a global ANOVA over the whole data set. It shows that the major source of variation was the gene factor, and that there was a highly significant interaction between the gene factor on the one hand, and time, TNF or PL on the other hand (*P*<0.0001), whereas the interaction between the gene factor and RBC appeared to be weaker (*P* = 0.04). In addition, the analysis revealed a significant interaction between gene, PL and TNF (*P* = 0.00029). Interestingly, RBC was included in high-order interaction terms, namely gene-by-time-by-RBC-by-TNF (*P* = 0.02), and gene-by-time-by-RBC-by-PL (*P* = 0.002).

Moreover, we performed an ANOVA for each gene to identify genes showing significant variation in their expression pattern. We assessed the effect of TNF, PL, RBC, and time on gene expression; we did not include interaction terms at this step. Figure 2A and 2B shows a graphical representation of the proportion of gene expression variance due to TNF and PL and the logarithm of *P*-value for each gene. Up to 90% of the variance was explained by TNF or PL for some genes (Figure 2A and 2B), whereas up to 50% of the variance was explained by time (Figure S1A). RBC explained less than 20% of gene expression variance (Figure S1B). Furthermore, we applied a FDR of 5% to account for multiple testing. We identified 39, 143, and 109 genes, the expression of which was significantly altered by PL (Table 3), TNF (Table S2), and time (Table S3), respectively, whereas there was no gene influenced by the factor RBC at this significance level.

We applied additional statistical analyses to identify genes, the expression of which was altered by one particular factor. We used the paired t-test and paired SAM method on the full data set (n = 4502), and applied a FDR of 5%. Paired t-tests yielded lists of 107 and 32 genes that showed a differential expression between experimental conditions with and without TNF and PL, respectively. Hierarchical clusterings of TNF- and PL-informative genes yielded with the t-test are presented in Figures 3 and 4 respectively. Seventy-one genes were also differentially expressed when comparing gene expression at 0 h and 5 h. Table 3 shows the genes transcriptionally modified by the platelets (PL) as determined by the three different methods. Comparison of the methods of analysis indicates that SAM identified the largest number of differentially expressed genes (54 genes), whereas

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 3.** List of genes transcriptionally altered by platelets.</caption><thead><tr><th>Gene</th><th>SAM<sup>a</sup></th><th>t-test</th><th>GeneAnova</th><th>% variance<sup>b</sup></th></tr></thead><tbody><tr><td>ADAM21</td><td>+</td><td>+</td><td>+</td><td>65.7</td></tr><tr><td>ADAM28</td><td>+</td><td>+</td><td>+</td><td>23.9</td></tr><tr><td>ANGPTL4</td><td>+</td><td>+</td><td>+</td><td>31.8</td></tr><tr><td>BCL3</td><td>+</td><td>+</td><td>+</td><td>24.4</td></tr><tr><td>BIRC2</td><td></td><td></td><td>+</td><td>17.4</td></tr><tr><td>BSCL2</td><td>+</td><td>+</td><td></td><td></td></tr><tr><td>CCL2</td><td>+</td><td>+</td><td>+</td><td>17.5</td></tr><tr><td>CCL4</td><td></td><td></td><td></td><td></td></tr><tr><td>CCL7</td><td></td><td></td><td></td><td></td></tr><tr><td>CCR7</td><td>+</td><td>+</td><td>+</td><td>19.5</td></tr><tr><td>CD83</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>CIRH1A</td><td>+</td><td></td><td>+</td><td>17.4</td></tr><tr><td>COPS7B</td><td></td><td></td><td>+</td><td>16.9</td></tr><tr><td>CSEN</td><td>+</td><td>(+)</td><td>+</td><td>17.1</td></tr><tr><td>CTGF</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>CXCL1</td><td>+</td><td>+</td><td>+</td><td>21</td></tr><tr><td>CXCL3</td><td>+</td><td>+</td><td>+</td><td>49.6</td></tr><tr><td>D2S448</td><td>+</td><td>+</td><td>+</td><td>16.6</td></tr><tr><td>DDIT4</td><td>+</td><td></td><td>+</td><td>18.3</td></tr><tr><td>DNASE1L3</td><td>+</td><td>+</td><td>+</td><td>18.2</td></tr><tr><td>DUSP5</td><td>+</td><td></td><td>+</td><td>8.7</td></tr><tr><td>EHD1</td><td>+</td><td>+</td><td>+</td><td>40.6</td></tr><tr><td>EPHA2</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>FAM19A4</td><td>+</td><td>+</td><td>+</td><td>24.5</td></tr><tr><td>FCER1G</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>FOSB</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>GLB1</td><td>+</td><td>+</td><td>+</td><td>17.4</td></tr><tr><td>H1FX</td><td></td><td></td><td></td><td></td></tr><tr><td>HBEGF</td><td>+</td><td>+</td><td>+</td><td>24.5</td></tr><tr><td>ID3</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>IER3</td><td>+</td><td>+</td><td>+</td><td>81.7</td></tr><tr><td>IL11</td><td>+</td><td>+</td><td>+</td><td>23.2</td></tr><tr><td>INHBA</td><td>+</td><td>+</td><td>+</td><td>52.2</td></tr><tr><td>INHBB</td><td>+</td><td>+</td><td>+</td><td>33.7</td></tr><tr><td>IRF1</td><td>+</td><td>+</td><td>+</td><td>30.6</td></tr><tr><td>KCNK3</td><td>+</td><td>+</td><td>+</td><td>40.2</td></tr><tr><td>MSC</td><td>+</td><td>+</td><td>+</td><td>39.6</td></tr><tr><td>MUC4</td><td>+</td><td>(+)</td><td>+</td><td>25.9</td></tr><tr><td>NFKBIA</td><td>+</td><td>+</td><td>+</td><td>66.7</td></tr><tr><td>ODC1</td><td></td><td></td><td></td><td></td></tr><tr><td>PAK1IP1</td><td>+</td><td>(+)</td><td>+</td><td>15.4</td></tr><tr><td>PDE1A</td><td>+</td><td>+</td><td>+</td><td>34.2</td></tr><tr><td>PDE7B</td><td></td><td></td><td></td><td></td></tr><tr><td>PGCP</td><td>+</td><td>+</td><td>+</td><td>18.3</td></tr><tr><td>PLK2</td><td>+</td><td>+</td><td>+</td><td>34</td></tr><tr><td>PSFL</td><td></td><td></td><td>+</td><td>17.5</td></tr><tr><td>PTK6</td><td>+</td><td>+</td><td>+</td><td>32.3</td></tr><tr><td>RAVER1</td><td>+</td><td>(+)</td><td></td><td></td></tr><tr><td>SACS</td><td>+</td><td>+</td><td></td><td></td></tr><tr><td>SLAMF8</td><td>+</td><td>+</td><td>+</td><td>71.4</td></tr></tbody></table>

<table><caption>**Table 3.** Cont.</caption><thead><tr><th>Gene</th><th>SAM<sup>a</sup></th><th>t-test</th><th>GeneAnova</th><th>% variance<sup>b</sup></th></tr></thead><tbody><tr><td>SORL1</td><td>+</td><td>(+)</td><td>+</td><td>20.2</td></tr><tr><td>STK23</td><td>+</td><td>+</td><td>+</td><td>44.3</td></tr><tr><td>TCEAL8</td><td></td><td></td><td></td><td></td></tr><tr><td>TNNI3</td><td>+</td><td>+</td><td></td><td></td></tr><tr><td>TPD52L1</td><td>+</td><td>+</td><td>+</td><td>17.8</td></tr><tr><td>TRAF4</td><td></td><td></td><td></td><td></td></tr><tr><td>UAP1</td><td></td><td></td><td></td><td></td></tr><tr><td>UCK2</td><td></td><td></td><td></td><td></td></tr></tbody></table>

+: informative genes with a FDR 5%.

(+): informative genes with a t-test 5% FDR 10%.

aSignificance Analysis of Microarray.

bpercentage of variance due to PL from multi-way ANOVA analysis.

doi:10.1371/journal.pone.0019651.t003

ANOVA and paired t-tests identified 39 and 32 genes respectively. When gene lists obtained with the t test were compared to those obtained with other methods, t-test exhibited 91% agreement (29 out of 32) with ANOVA, and 100% agreement with SAM (32 out of 32). Similar results were obtained when the genes transcriptionally modified by TNF and time were determined by the different methods (supplementary Tables S2 and S3, respectively). Overall, the different methods confirmed that the expression of a significant number of genes was influenced by TNF, platelets or time.

In contrast, there was no gene showing differential expression when comparing experimental conditions with NRBC and those with pRBC with a FDR of 5%. It should be stressed, however, that the SAM method identified 11 genes when we applied a FDR of 35%. Interestingly, TARBP2, ALS2CL, P4HA1, KLK7, POR, RANBP9, Q9H693, and SCARB1 that were identified by SAM were also differentially expressed as determined by ANOVA with a nominal threshold of *P*<0.01. Besides, we assumed that the effect of pRBC *per se* was sub-optimal, and that it could be detectable when accounting for putative interaction with TNF and/or PL. We found evidence of an interaction between gene, pRBC, TNF and time (*P*=0.02) and between gene, pRBC, PL and time (*P*=0.002) by using the global GeneANOVA (Table 2). This indicates that the effect of the interactions between pRBC and TNF and between pRBC and PL depended on time. Therefore, we separately analyzed gene expression at time 0 h and 5 h, performed a two-way ANOVA, and applied a FDR of 5% to identify genes involved. There was a significant interaction between pRBC and TNF at time 5 h for 4 genes (KTN1, RRBPI, GLIPR1, CORO2B) and a significant interaction between pRBC and PL for 4 genes at time 0 h (ADAMTS16, RREB1, EFNB3, PLXDC1), and for 3 genes at time 5 h (ASB12, AP3M2, LAMP1). This suggests that pRBC act in interaction with PL and TNF to induce a transcriptional response in HBEC.

The Gene Ontology (GO) terms and KEGG annotations of genes included in the different lists were analyzed by using DAVID tools. Table 4 shows the results of the functional annotation analysis. There was an over-representation of GO terms, such as “immune response”, “cell adhesion” and “apoptosis” in the list of genes regulated by TNF. The analysis of GO terms of genes regulated by PL showed an over-representation of “apoptosis” and “cytokine activity”. The analysis of GO terms of the 469 genes identified by the one-way ANOVA yielded similar results (data not shown).

<figure><img src="image_12.png" /></figure>
<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>
**Figure 3. Hierarchical clustering of TNF-induced genes.** Hundred and seven TNF-induced genes were selected on the basis of a Welch t-test and a FDR of 5%. The samples obtained from HBEC stimulated with TNF were compared with those from HBEC incubated without TNF. The presence (black) and the absence (white) of TNF are shown. Genes marked with an asterisk were validated by RT-qPCR.
doi:10.1371/journal.pone.0019651.g003

We used IPA software to visualize differentially expressed genes in networks of interacting genes and to identify the biological functions that were the most significant in our data set. We observed that 56 of the 58 genes differentially expressed between experimental conditions with and without PL could be assigned to one among five related networks at the 0.001 significance level (Table S4). The analysis of the canonical pathways indicated that HBEC has altered activities of several functional pathways, including cytokine-, chemokine-, death receptor-, apoptosis-, erythropoietin-, and TGFβ- signaling (Table S5). Figure 5 shows the most significant network, the top functions of which were haematological system development and immune system. We identified, in this network, several canonical pathways significantly represented; these included several cytokine- and chemokine- signaling pathways, but also glucocorticoid receptor-, PPAR- and TREM1-signaling pathways.

We further analyzed the expression of 11 selected genes by RT-qPCR to confirm their gene expression pattern. Transcript levels for each of the 11 genes were quantified by using SYBR Green technology. Figure 6A shows differences in their expression i) between endothelial cells incubated with TNF and those incubated without TNF and ii) between endothelial cells cocultured with PL and those incubated without PL. Microarray and RT-qPCR results were consistent, except for *CCR7*. Most importantly, the positive or negative changes of expression remained consistent. As expected, there was a difference of sensitivity between microarray and RT-qPCR methods. RT-qPCR results of selected genes ranged from 1.5-fold to 36-fold regulation (Figure 6A), whereas microarray results ranged from 1.5-fold to 10-fold regulation for the same genes (Figure 6B). The fold changes in transcript abundance of *CCR7* strongly differed between RT-qPCR and

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

**Figure 4. Hierarchical clustering of platelet-induced genes.** A Welch t-test and a FDR of 5% led to identify 32 PL-induced genes. The samples obtained from HBEC co-cultured with PL were compared with those from HBEC incubated without PL. The presence (black) and the absence (white) of PL are shown. Genes marked with an asterisk were validated by RT-qPCR.
doi:10.1371/journal.pone.0019651.g004

<figure><img src="image_16.png" /></figure>
<table><caption>**Table 4.** Functional annotation clustering of genes differentially expressed in the presence of TNF, platelets, and between 0 h and 5 h of incubation.</caption><thead><tr><th></th><th>Functional annotation (GO terms)</th><th>% of genes</th><th>P-value</th></tr></thead><tbody><tr><td rowspan="6">TNF-regulated genes<sup>a</sup>:<br>164 genes annotated/185 genes</td><td>immune response</td><td>16</td><td>6×10<sup>-6b</sup></td></tr><tr><td>plasma membrane</td><td>29</td><td>9×10<sup>-4b</sup></td></tr><tr><td>localization</td><td>29</td><td>7×10<sup>-3</sup></td></tr><tr><td>cell adhesion</td><td>10</td><td>3×10<sup>-2</sup></td></tr><tr><td>regulation of signal transduction</td><td>8</td><td>4×10<sup>-2</sup></td></tr><tr><td>apoptosis</td><td>11</td><td>5×10<sup>-2</sup></td></tr><tr><td rowspan="3">PL-regulated genes<sup>a</sup>:<br>53 genes annotated/58 genes</td><td>protein kinase cascade</td><td>7</td><td>5×10<sup>-2</sup></td></tr><tr><td>cytokine activity</td><td>15</td><td>2×10<sup>-5b</sup></td></tr><tr><td>apoptosis</td><td>25</td><td>2×10<sup>-4b</sup></td></tr><tr><td rowspan="4">Time-regulated genes<sup>a</sup>:<br>136 genes annotated/158 genes</td><td>I-kappaB kinase/NF-kappaB cascade</td><td>8</td><td>2×10<sup>-2</sup></td></tr><tr><td>B-cell activation</td><td>6</td><td>3×10<sup>-2</sup></td></tr><tr><td>heparin binding</td><td>6</td><td>3×10<sup>-2</sup></td></tr><tr><td>protein kinase cascade</td><td>10</td><td>1×10<sup>-3b</sup></td></tr><tr><td></td><td>response to biotic stimulus</td><td>6</td><td>2×10<sup>-2</sup></td></tr><tr><td></td><td>intracellular</td><td>74</td><td>3×10<sup>-2</sup></td></tr><tr><td></td><td>enzyme linked receptor protein signaling</td><td>6</td><td>3×10<sup>-2</sup></td></tr><tr><td></td><td>cell maturation</td><td>3</td><td>3×10<sup>-2</sup></td></tr><tr><td></td><td>JNK cascade</td><td>3</td><td>4×10<sup>-2</sup></td></tr></tbody></table>

microarray techniques. Transcript levels for *CCR7* were the lowest ones among those for the 10 selected genes by using RT-qPCR, and very low signal intensities were measured for *CCR7* by using the microarray technology; this indicates that differences in gene expression are more difficult to quantify by using microarray technology for low transcript levels.

Discussion

Here we have investigated the respective effects of TNF, PL and pRBC on human brain endothelial cells by profiling gene expression from *in vitro* co-cultures. Unsupervised hierarchical clustering of samples indicated an effect of a short exposure to PL, besides a strong effect of a prolonged exposure to TNF. We applied stringent statistical analyses to identify genes differentially expressed between the experimental conditions, and we searched for over-representation of functional annotations to identify pathways involved.

As expected, TNF regulated the expression of genes involved in inflammation, such as *IL32*, *IL11*, *CCL18*, *CXCL10*, *CCL2*, and *CCL5*, and apoptosis, such as *SERPINB2*, *TNFSF13B* and *TNFSF10*. Similar results have been previously obtained with non-cerebral endothelial cells [8,31,32,33]. Overall, our results confirm that endothelial cells activated by TNF contribute to endothelial cell apoptosis and leukocyte chemotaxis and activation. Interestingly, we detected a significant down-regulation of *LRPAP1* and *LRP1* in TNF-stimulated HBEC. LRP1 that binds to LRPAP1 has been shown to be involved in the clearance of the β-amyloid protein [34], which accumulates in the brain of patients with Alzheimer's disease [35]. Since the β-amyloid protein has been

also detected in brains during human and mouse CM [36,37], our results suggest that TNF inhibits the capacity of HBEC to eliminate the β-amyloid protein from brain. This further supports the hypothesis that CM and Alzheimer's disease share some common mechanisms of pathogenesis [37].

Here we show that platelets dramatically modulate the expression of genes involved in inflammation and apoptosis. The analysis of canonical pathways revealed the effect of PL on cytokine-, chemokine-, TGFβ-, death-receptor-, apoptosis-, erythropoietin-, and TREM1-signaling. These observations are in line with evidence that i) TGFβ participates in HBEC apoptosis induced by PL [38], ii) erythropoietin promotes endothelial cell integrity and angiogenesis [39] and has been associated with resistance against CM in humans and mice [40,41], iii) TREM-1, which is expressed on myeloid and endothelial cells, has been shown to be involved in the innate inflammatory response and sepsis [42,43], and its ligand is expressed on human PL [44]. All together, these data support the hypothesis that PL induce HBEC gene expression changes affecting the outcome of malaria infection. It should be stressed that we used platelets from one single donor, and that this is a limitation in our approach. However, those platelets did not show any abnormalities in term of reactivity and probably induced a typical response of HBEC.

It is likely that both TNF and PL can contribute to apoptosis and inflammation, suggesting that TNF and platelets work in synergy, contributing to the same cell response. Indeed the expression of some genes, including *ADAM28*, *DNASEIL3*, and *CCL2*, appeared to be regulated both by TNF and platelets, and the expression of other genes (*GPR83*, *IER3*, *RAB4*, and *HLA-DRB1*) was altered by an interaction between TNF and PL. This

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

<footnote><sup>a</sup>Informative genes identified by either paired t-test, paired SAM or ANOVA analyses.</footnote>

<footnote><sup>b</sup>FDR<5%.</footnote>

<footnote>doi:10.1371/journal.pone.0019651.t004</footnote>
<figure><img src="image_18.png" /><figcaption>**Figure 5.** A gene network of platelet-regulated genes. This network is the most significative network among five significative networks obtained using the 58 PL-regulated genes with IPA (supplementary data). Twenty genes of this network are differentially expressed in HBEC in presence of PL (colored in green for the down-regulated gene and red for up-regulated genes). The molecules hedged in orange belong to the main canonical pathways linked to this network (Acute phase response-, B cell receptor-, Chemokine-, Erythropoietin-, Glucocorticoid receptor-, IL6-, IL8-, IL10-, IL12-, PPAR-, RAR-, TREM-1-signaling, IRF, LXR/RXR, PPARa/RXRa, RAR activation, Role of pattern recognition receptor, Role of PKR in IFN induction, Hepatic cholestasis, Hepatic fibrosis). doi:10.1371/journal.pone.0019651.g005</figcaption></figure>

interaction underlines the role of TNF in the induction of PL adherence to cerebrovascular endothelium, notably *via* the overexpression of ICAM-1 on EC surface [45] and confirms that PL potentiate the pro-apoptotic effect of TNF on HBEC [7]. Interestingly, TNF and PL also induced the expression of genes with an anti-inflammatory and/or anti-apoptotic activity, such as *IL11* and *INHBA*. It should be stressed that we analyzed gene expression immediately and 5 hours after co-cultures, whereas apoptosis was detected in endothelial cells only 48 h after platelets co-incubation [7]. Chakravorty *et al* recently proposed that the balance between positive and negative regulation will determine endothelial pathology during malaria infection [12]. Our result highlights that a pro-apoptotic stimulus also induces regulatory genes, and fits well with the observation of a transient dysfunction of the blood brain barrier in human severe malaria [46].

When compared to previous reports [10,12,13,14], we found that the likelihood of pRBC-induced transcriptional changes were much lower than those triggered by PL or TNF, i.e. required a FDR of 35% instead of 5%. Such discrepancies may reflect

differences in endothelial cell lines, in the durations of exposure to pRBC, and in parasite strains. For instance, we co-cultured HBEC and *P. falciparum* from the IPPAM strain with a polymorphic adhesion phenotype, equally CD36- and ICAM-1-binding, whereas others incubated endothelial cells (HUVEC, HLEC, HBEC...) with *P. falciparum* ItG and ITO4-A4 (ICAM-1 binding) or 3D7 and FCR3 (CD36-binding) strains. In addition, we co-cultured HBEC and pRBC for a short time to detect early events, whereas others detected transcriptional changes after a prolonged co-culture, which may favour the effect of pRBC on endothelial cells [10,12,13,14]. We assume that our system was not optimal to detect the effect of pRBC *per se*, but might be more relevant physiologically to analyse the early effect of adherent pRBC. Besides, we cannot exclude that we missed some genes of interest that were not spotted on the microarrays. Interestingly, Viebig *et al.*, Pino *et al.*, Tripathi *et al.*, and Chakravorty *et al.* collectively provided evidence that pRBC represent a pro-inflammatory (cell adhesion molecules, chemotaxis and chemokine) and a pro-apoptotic stimulus for endothelial cells [10,12,13,14]. This suggests

<figure><img src="image_19.png" /></figure>
A

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

B

<figure><img src="image_21.png" /><figcaption>**Figure 6.** Comparison of the expression of selected genes assessed by RT-qPCR and by microarray technology. Five TNF-induced genes and 6 PL-induced genes were selected from the list of genes identified through Welch t-tests and multiple test corrections. For 2 genes (CCL2 and IL11), gene expression was regulated by TNF and PL. HBEC were co-cultured with NRBC in the presence or the absence of either TNF or PL to measure the expression of TNF- or PL-specific genes, respectively. A) Gene expression was assessed by RT-qPCR. Data were analyzed by the 2⁻⁰ΔΔC(t) method with ACTB as a control gene. The n-fold value represents the mean of TNF- and PL-induced expression levels compared with control conditions. Error bar represents SD between duplicates. B) Gene expression was assessed by microarray technology. The n-fold value was calculated on the basis of normalized data when comparing the levels of TNF- and PL-induced expression with those of control conditions. Error bar represents SD between triplicates.
doi:10.1371/journal.pone.0019651.g006</figcaption></figure>

that pRBC together with TNF and PL contribute to similar endothelial responses. Using a global ANOVA, we observed a time-dependent interaction between pRBC and TNF and between pRBC and PL, and we further identified some genes by using a two-way ANOVA, whose expression was influenced by these high-order interaction terms. Although the power may be insufficient to identify such genes in our analysis, our results suggest that pRBC together with TNF or PL contribute to endothelial responses.

In conclusion, this study shows that both TNF and PL dramatically alter gene expression levels in HBEC, and suggests that pRBC act in interaction with TNF and PL to induce transcriptional changes in HBEC. It also indicates that TNF and PL contribute to HBEC apoptosis. This study further supports the pathogenic role of platelets, and points out several altered canonical pathways, such as chemokine-, TGFβ-, death-receptor, erythropoietin-, and TREM1-signaling. We speculate that these observations can be used to develop novel strategies to prevent CM.

**Supporting Information**

**Figure S1** Gene expression variance analysis for time (A) and RBC (B) factors. Each gene was plotted as a point. The abscissa is the variation due to the factor normalized by the total variation of the gene, and the ordinate is the logarithm of the P-value. (DOC)

**Table S1** Primers used for the RT-qPCR. (DOC)

**Table S2** List of genes transcriptionally altered by TNF. (DOC)

**Table S3** List of genes transcriptionally altered by time. (DOC)

**Table S4** IPA Network analysis based on the 58 platelet-regulated genes. (DOC)

<figure><img src="image_22.png" /></figure>
**Table S5** IPA Canonical Pathways analysis based on the 58 platelet-regulated genes.
(DOC)

**Acknowledgments**

We thank Albin Fontaine for technical support. We acknowledge the technical support of Marseille-Nice Genopole.

**References**

1. Marsh K, Snow RW (1999) Malaria transmission and morbidity. *Parassitologia* 41: 241–246.

2. van der Heyde HC, Nolan J, Combes V, Gramaglia I, Grau GE (2006) A unified hypothesis for the genesis of cerebral malaria: sequestration, inflammation and hemostasis leading to microcirculatory dysfunction. *Trends Parasitol* 22: 503–508.

3. Flori L, Delahaye NF, Iraqi FA, Hernandez-Valladares M, Fumoux F, et al. (2005) TNF as a malaria candidate gene: polymorphism-screening and family-based association analysis of mild malaria attack and parasitemia in Burkina Faso. *Genes Immun* 6: 472–480.

4. Grau GE, Mackenzie CD, Carr RA, Redard M, Pizzolato G, et al. (2003) Platelet accumulation in brain microvessels in fatal pediatric cerebral malaria. *J Infect Dis* 187: 461–466.

5. Sun G, Chang WL, Li J, Berney SM, Kimpel D, et al. (2003) Inhibition of platelet adherence to brain microvasculature protects against severe Plasmodium berghei malaria. *Infect Immun* 71: 6553–6561.

6. Wassmer SC, Lepolard C, Traore B, Pouvelle B, Gysin J, et al. (2004) Platelets reorient Plasmodium falciparum-infected erythrocyte cytoadhesion to activated endothelial cells. *J Infect Dis* 189: 180–189.

7. Wassmer SC, Combes V, Candal FJ, Juhan-Vague I, Grau GE (2006) Platelets potentiate brain endothelial alterations induced by Plasmodium falciparum. *Infect Immun* 74: 645–653.

8. Sana TR, Janatpour MJ, Sathe M, McEvoy LM, McClanahan TK (2005) Microarray analysis of primary endothelial cells challenged with different inflammatory and immune cytokines. *Cytokine* 29: 256–269.

9. Viemann D, Goebeler M, Schmid S, Nordhues U, Klimmek K, et al. (2006) TNF induces distinct gene expression programs in microvascular and macrovascular human endothelial cells. *J Leukoc Biol* 80: 174–185.

10. Viebig NK, Wulbrand U, Forster R, Andrews KT, Lanzer M, et al. (2005) Direct activation of human endothelial cells by Plasmodium falciparum-infected erythrocytes. *Infect Immun* 73: 3271–3277.

11. Susomboon P, Maneerat Y, Dekumyoy P, Kalambaheti T, Iwagami M, et al. (2006) Down-regulation of tight junction mRNAs in human endothelial cells co-cultured with Plasmodium falciparum-infected erythrocytes. *Parasitol Int* 55: 107–112.

12. Chakravorty SJ, Carret C, Nash GB, Ivens A, Szestak T, et al. (2007) Altered phenotype and gene transcription in endothelial cells, induced by Plasmodium falciparum-infected red blood cells: Pathogenic or protective? *Int J Parasitol*.

13. Pino P, Vouldoukis I, Kolb JP, Mahmoudi N, Desportes-Livage I, et al. (2003) Plasmodium falciparum-infected erythrocyte adhesion induces caspase activation and apoptosis in human endothelial cells. *J Infect Dis* 187: 1283–1290.

14. Tripathi AK, Sha W, Shulaev V, Stins MF, Sullivan DJ Jr. (2009) Plasmodium falciparum-infected erythrocytes induce NF-kappaB regulated inflammatory pathways in human cerebral endothelium. *Blood* 114: 4243–4252.

15. Combes V, El-Assaad F, Faille D, Jambou R, Hunt NH, et al. (2010) Microvesiculation and cell interactions at the brain-endothelial interface in cerebral malaria pathogenesis. *Prog Neurobiol* 91: 140–151.

16. Dorovini-Zis K, Prameya R, Bowman PD (1991) Culture and characterization of microvascular endothelial cells derived from human brain. *Lab Invest* 64: 425–436.

17. Fried M, Duffy PE (1996) Adherence of Plasmodium falciparum to chondroitin sulfate A in the human placenta. *Science* 272: 1502–1504.

18. Trager W, Jensen JB (1976) Human malaria parasites in continuous culture. *Science* 193: 673–675.

19. Heidrich HG, Mrema JE, Vander Jagt DL, Reyes P, Rieckmann KH (1982) Isolation of intracellular parasites (Plasmodium falciparum) from culture using free-flow electrophoresis: separation of the free parasites according to stages. *J Parasitol* 68: 443–450.

20. Bertucci F, Pages C, Finetti P, Rochaix P, Lamant L, et al. (2007) Gene expression profiling of human melanoma cell lines with distinct metastatic potential identifies new progression markers. *Anticancer Res* 27: 3441–3449.

21. Puthier D, Joly F, Irla M, Saade M, Victorero G, et al. (2004) A general survey of thymocyte differentiation by transcriptional analysis of knockout mouse models. *J Immunol* 173: 6109–6118.

22. Bertucci F, Bernard K, Loriod B, Chang YC, Granjeaud S, et al. (1999) Sensitivity issues in DNA array-based expression measurements and performance of nylon microarrays for small samples. *Hum Mol Genet* 8: 1715–1722.

**Author Contributions**

Wrote the paper: MB DF PR. Conceived the study and participated in its design and coordination: MB DF GEG PR. Provided material and analytical tools: M-CA TF SCW CN. Performed the co-cultures: DF. Performed RNA extraction: MB. Responsible for the design and production of the microarray and participated in the acquisition of microarray data: BL GV JT DP. Performed cDNA hybridization and acquisition of microarray data: MB CC. Performed data analysis: MB LF PR. All authors contributed to and approved the final manuscript.

23. Lopez F, Rougemont J, Loriod B, Bourgeois A, Loi L, et al. (2004) Feature extraction and signal processing for nylon DNA microarrays. *BMC Genomics* 5: 38.

24. Talby L, Chambost H, Roubaud MC, N'Guyen C, Milili M, et al. (2006) The chemosensitivity to therapy of childhood early B acute lymphoblastic leukemia could be determined by the combined expression of CD34, SPI-B and BCR genes. *Leuk Res* 30: 665–676.

25. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, et al. (2004) Bioconductor: open software development for computational biology and bioinformatics. *Genome Biol* 5: R80.

26. Didier G, Brezellec P, Remy E, Henaut A (2002) GeneANOVA-gene expression analysis of variance. *Bioinformatics* 18: 490–491.

27. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate: a practical and powerful approach to multiple testing. *J Royal Stat Soc Ser B* 57: 289–300.

28. Reiner A, Yekutieli D, Benjamini Y (2003) Identifying differentially expressed genes using false discovery rate controlling procedures. *Bioinformatics* 19: 368–375.

29. Kerr MK, Martin M, Churchill GA (2000) Analysis of variance for gene expression microarray data. *J Comput Biol* 7: 819–837.

30. Dennis G, Jr., Sherman BT, Hosack DA, Yang J, Gao W, et al. (2003) *DAVID*: Database for Annotation, Visualization, and Integrated Discovery. *Genome Biol* 4: P3.

31. Zoellner H, Wojta J, Gallicchio M, McGrath K, Hamilton JA (1993) Cytokine regulation of the synthesis of plasminogen activator inhibitor-2 by human vascular endothelial cells. Comparison with plasminogen activator inhibitor-1 synthesis. *Thromb Haemost* 69: 135–140.

32. Hillyer P, Mordelet E, Flynn G, Male D (2003) Chemokines, chemokine receptors and adhesion molecules on different human endothelium: discriminating the tissue-specific functions that affect leucocyte migration. *Clin Exp Immunol* 134: 431–441.

33. Li JH, Kirkles-Smith NC, McNiff JM, Pober JS (2003) TRAIL induces apoptosis and inflammatory gene expression in human endothelial cells. *J Immunol* 171: 1526–1533.

34. Shibata M, Yamada S, Kumar SR, Calero M, Bading J, et al. (2000) Clearance of Alzheimer's amyloid-ss(1-40) peptide from brain by LDL receptor-related protein-1 at the blood-brain barrier. *J Clin Invest* 106: 1489–1499.

35. Donahue JE, Flaherty SL, Johanson CE, Duncan JA, 3rd, Silverberg GD, et al. (2006) RAGE, LRP-1, and amyloid-beta protein in Alzheimer's disease. *Acta Neuropathol* 112: 405–415.

36. Medana IM, Day NP, Hien TT, Mai NT, Bethell D, et al. (2002) Axonal injury in cerebral malaria. *Am J Pathol* 160: 655–666.

37. Delahaye NF, Coltel N, Puthier D, Barbier M, Benech P, et al. (2007) Gene expression analysis reveals early changes in several molecular pathways in cerebral malaria-susceptible mice versus cerebral malaria-resistant mice. *BMC Genomics* 8: 452.

38. Wassmer SC, de Souza JB, Frere C, Candal FJ, Juhan-Vague I, et al. (2006) TGF-beta1 released from activated platelets can induce TNF-stimulated human brain endothelium apoptosis: a new mechanism for microvascular lesion during cerebral malaria. *J Immunol* 176: 1180–1184.

39. Chong ZZ, Kang JQ, Maiese K (2002) Angiogenesis and plasticity: role of erythropoietin in vascular systems. *J Hematother Stem Cell Res* 11: 863–871.

40. Casals-Pascual C, Idro R, Gicheru N, Gwer S, Kitsao B, et al. (2008) High levels of erythropoietin are associated with protection against neurological sequelae in African children with cerebral malaria. *Proc Natl Acad Sci U S A* 105: 2634–2639.

41. Kaiser K, Texier A, Ferrandiz J, Buguet A, Meiller A, et al. (2006) Recombinant human erythropoietin prevents the death of mice during cerebral malaria. *J Infect Dis* 193: 987–995.

42. Chen LC, Laskin JD, Gordon MK, Laskin DL (2008) Regulation of TREM expression in hepatic macrophages and endothelial cells during acute bouchouxemia. *Exp Mol Pathol* 84: 145–155.

43. Haselmayer P, Grosse-Hovest L, von Landenberg P, Schild H, Radsak MP (2007) TREM-1 ligand expression on platelets enhances neutrophil activation. *Blood* 110: 1029–1035.

<figure><img src="image_23.png" /></figure>
45. Lou J, Donati YR, Juillard P, Giroud C, Vesin C, et al. (1997) Platelets play an important role in TNF-induced microvascular endothelial cell pathology. Am J Pathol 151: 1397-1405.

46. Brown H, Hien TT, Day N, Mai NT, Chuong LV, et al. (1999) Evidence of blood-brain barrier dysfunction in human cerebral malaria. Neuropathol Appl Neurobiol 25: 331-340.

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