# Genome-Wide Expression Profiling Deciphers Host Responses Altered during Dengue Shock Syndrome and Reveals the Role of Innate Immunity in Severe Dengue

Stephanie Devignot, Cedric Sapet, Veasna Duong, Aurélie Bergon, Pascal Rihet, Sivuth Ong, Patrich T. Lorn, North Chroeung, Sina Ngeav, Hugues J. Tolou, et al.

► To cite this version:

Stephanie Devignot, Cedric Sapet, Veasna Duong, Aurélie Bergon, Pascal Rihet, et al. Genome-Wide Expression Profiling Deciphers Host Responses Altered during Dengue Shock Syndrome and Reveals the Role of Innate Immunity in Severe Dengue. PLoS ONE, 2010, 5 (7), e11671. doi:10.1371/journal.pone.0011671

HAL Id: hal-01595911

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

Submitted on 17 Sep 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
# Genome-Wide Expression Profiling Deciphers Host Responses Altered during Dengue Shock Syndrome and Reveals the Role of Innate Immunity in Severe Dengue

Stéphanie Devignot¹, Cédric Sapet¹,², Veasna Duong²,³, Aurélie Bergon³, Pascal Rihet³, Sivuth Ong², Patrich T. Lorn⁴, North Chroeung⁴, Sina Ngeav⁴, Hugues J. Tolou¹, Philippe Buchy², Patricia Couissinier-Paris¹*

1. French Army Biomedical Research Institute (Institut de recherche biomédicale des armées, IRBA), Antenne de Marseille-IMTSSA; Unité de Virologie, Marseille, France,

2. Institut Pasteur in Cambodia, Department of Virology, Phnom Penh, Cambodia, 3 TAGC-INSERM U928, Marseille, France, 4 Kampong Cham Provincial Hospital, Kampong Cham, Cambodia

## Abstract

**Background:** Deciphering host responses contributing to dengue shock syndrome (DSS), the life-threatening form of acute viral dengue infections, is required to improve both the differential prognosis and the treatments provided to DSS patients, a challenge for clinicians.

**Methodology/Principal Findings:** Based on a prospective study, we analyzed the genome-wide expression profiles of whole blood cells from 48 matched Cambodian children: 19 progressed to DSS while 16 and 13 presented respectively classical dengue fever (DF) or dengue hemorrhagic fever grades I/II (DHF). Using multi-way analysis of variance (ANOVA) and adjustment of p-values to control the False Discovery Rate (FDR<10%), we identified a signature of 2959 genes differentiating DSS patients from both DF and DHF, and showed a strong association of this DSS-gene signature with the dengue disease phenotype. Using a combined approach to analyse the molecular patterns associated with the DSS-gene signature, we provide an integrative overview of the transcriptional responses altered in DSS children. In particular, we show that the transcriptome of DSS children blood cells is characterized by a decreased abundance of transcripts related to T and NK lymphocyte responses and by an increased abundance of anti-inflammatory and repair/remodeling transcripts. We also show that unexpected pro-inflammatory gene patterns at the interface between innate immunity, inflammation and host lipid metabolism, known to play pathogenic roles in acute and chronic inflammatory diseases associated with systemic vascular dysfunction, are transcriptionally active in the blood cells of DSS children.

**Conclusions/Significance:** We provide a global while non exhaustive overview of the molecular mechanisms altered in of DSS children and suggest how they may interact to lead to final vascular homeostasis breakdown. We suggest that some mechanisms identified should be considered putative therapeutic targets or biomarkers of progression to DSS.

**Citation:** Devignot S, Sapet C, Duong V, Bergon A, Rihet P, et al. (2010) Genome-Wide Expression Profiling Deciphers Host Responses Altered during Dengue Shock Syndrome and Reveals the Role of Innate Immunity in Severe Dengue. PLoS ONE 5(7): e11671. doi:10.1371/journal.pone.0011671

**Editor:** Patricia T. Bozza, Fundação Oswaldo Cruz, Brazil

**Received** January 20, 2010; **Accepted** June 22, 2010; **Published** July 20, 2010

**Copyright:** © 2010 Devignot 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:** This work was supported by a clinical research program PRC 2007-13 from the French Army Medical Health Service (Service de Santé des Armées Françaises). The funders had no role in the 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

Acute dengue virus infections are a major public health problem for many tropical and sub-tropical countries and an increasing risk for the worldwide population [1]. Symptomatic infections occur under a spectrum of diseases ranging from classical dengue fever (DF) to the most severe life-threatening dengue shock syndrome (DSS), a leading cause of childhood hospitalisation and death in endemic countries with limited health resources [1,2].

DSS is regarded as a vascular disease involving a complex interplay between virus, whole blood cells and microvascular territories [3,4], and thought to result largely from an aberrant host response to infection. As for other major systemic diseases, a

detrimental cytokine storm is thought to be central to the systemic microcirculatory failure and massive plasma leakage leading to cardiovascular decompensation characterizing DSS [5]. However, controversies exist regarding the nature of pathogenic host immune responses supporting this life-threatening syndrome [6-8]. Indeed, reactivation of cross-reactive memory T lymphocytes and increased infection of monocytes mediated by cross-reactive antibodies acquired during previous infections by distinct dengue virus serotypes, are the main hypothetic mechanisms proposed to explain the putative cytokine storm leading to plasma leakage [5,9]. However, those hypothesis fail to explain the occurrence of DSS in patients having primary dengue infection and their relevance to the pathophysiology of DSS disease is discussed [8,10].

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

<footnote>* E-mail: parisp@imtssa.fr</footnote>

<footnote>¶ These authors contributed equally to this work.</footnote>
Efforts to identify soluble biomarkers of severe dengue differentiating uncomplicated dengue infections from severe ones has led to the identification of a diversity of cytokines, chemokines, endothelial agonists or soluble endothelial molecules [11–18]. However, discrepancies in definition of dengue severity, variability in patients cohorts characteristics, as well as in techniques and markers investigated, have impaired the identification of reliable sets of DSS biomarkers and the possibility to get a global overview of biological markers altered during DSS.

Understanding the molecular basis of DSS and identifying relevant DSS biomarkers thus remains a major challenge [5,6]. Indeed, DSS occurs by the end of the acute infection in only a fraction of dengue-infected patients and current severity criteria, based on the 1997 World Health Organization (WHO) classification of dengue severity, fail to predict a significant proportion of patients who progress to life-threatening DSS [19–21]. Attempting to decipher molecular mechanisms underlying DSS by analyzing circulating whole blood cell genome-wide expression profiles is a relevant approach regarding the study of other systemic inflammatory syndromes, where a cognate cross-talk between endothelial vascular cells and blood cells occurs [22–24]. Whole blood represent a highly informative while complex cellular sample, that may reflect host pathophysiological responses ongoing at the time of blood sampling [22]. Furthermore, whole blood cells are easy to collect and store during field studies on large cohorts, reducing samples volumes required and limiting technical bias due to cell purification. However, due to the high cellular complexity of whole blood cells samples, whole gene expression patterns should be carefully analyzed and deciphered to allow returning to an integrative view of the molecular mechanisms altered during the pathophysiological process studied [25].

Such a bench-to-bedside medical research has gained more and more interest in the recent years. Indeed, it allowed improving the understanding of pathophysiological processes underlying systemic critical illnesses such as sterile and non sterile systemic inflammatory responses syndromes (SIRS), allowing the identification of relevant disease biomarkers and of new putative therapeutic targets [22–24,26].

Genome-wide expression studies aimed at deciphering molecular responses altered in the whole blood cells of adults [27] and children DSS patients [28] have been implemented recently by colleagues. They reported a decreased IFN type I-induced response and a benign transcriptional response at the time of cardiovascular decompensation [27,28], but failed in identifying biological pathways relevant to DSS pathophysiology and particularly inflammatory ones that could sustain microvascular dysfunction [28].

We report here the results of a prospective study comparing the whole blood genome-wide expression profiles of 48 matched Cambodian children recruited during the huge 2007 dengue outbreak who presented with classical dengue fever (DF), dengue hemorrhagic fever grades I/II (DHF) or dengue shock syndrome (DSS), according to the 1997 WHO classification of dengue severity [29]. Based on careful study design and statistic treatment of microarrays data, we identified a large and highly relevant gene signature of DSS never reported before, that discriminates DSS children from paediatric patients with DF or DHF grades I/II, who did not present severe clinical complications. Using an integrative analysis of the gene patterns altered in DSS children, we deciphered part of the complex interactive molecular processes occurring during DSS, highlighting similarities between DSS and other major inflammatory processes. Finally, we identified unexpected pro-inflammatory innate immune responses activated in the whole blood cells of DSS children that may play a major

role in DSS pathophysiology. The implications of present findings to the improvement of DSS prognosis and treatment are discussed.

# Materials and Methods

## Ethics statement

The global study and all protocols presented here were approved by the national Cambodian ethical committee. Written informed consent was obtained from the legal guardians of each child. To ensure strict anonymity regarding the patients, samples were encoded as PLxxx (Plasma Leakage).

## Patients and clinical data

Inclusion criteria retained were: age (1 to 15 years old); positive diagnosis of acute dengue infection assessed by different methods; absence of known chronic inflammatory disease or ongoing acute co-infection at the time of inclusion.

An eligible cohort of 83 dengue-infected children hospitalised at the Kampong Cham provincial hospital, Cambodia, was prospectively enrolled from July to September 2007 during the huge 2007 dengue outbreak in Cambodia, characterized by a high number of DSS cases.

Children diagnosed with acute dengue infections were classified at admission as classical dengue fever (DF), dengue hemorrhagic fever (DHF) or dengue shock syndrome (DSS) based on the 1997 WHO criteria [29]. Clinical and biological follow-up was done daily for each hospitalised patient. DSS patients were admitted to hospital intensive care unit where they received appropriate fluid resuscitation and were monitored for vital parameters. Children who required blood transfusion were not included in the study.

To increase the probability to identify gene signatures specific of DSS, we chose to include only symptomatic dengue-infected classified DF, DHF and DSS, but no healthy or non-dengue children in the present study. This is based on the rationale that comparing DF, DHF and DSS patients together should improve the probability to identify a DSS-specific gene signature, while including an external non dengue control group should increase the probability to identify a general dengue-related signature but should be less powerful at identifying a signature of severe dengue disease.

DF, DHF and DSS patients whole blood samples selected for the present study corresponded to comparable duration of illness after onset of fever: all were collected within a window of time comprised between 3 days and 7 days after onset of fever (being considered day 0). For most DSS patients, this generally corresponded to the day of cardiovascular decompensation (shock) or the day after, except for 3 (PL017, PL033, PL047) and 2 DSS (PL005, PL101) for whom blood was collected respectively 2 and 3 days after onset of shock.

Patients' samples selected for the present study were also carefully matched for age, gender, viral serotype (when identified) and immunological status (primary or secondary, according to reference assays described in diagnosis methods) towards dengue infection.

Diagnosis assays carried out as described thereafter, indicated that about 90% of all dengue-infected children had secondary infection.

## Dengue diagnosis and immunological status

All diagnosis assays were carried out at the Institut Pasteur in Cambodia, the National Reference Center for arboviral diseases in Cambodia. IgM capture ELISA and Hemagglutination-inhibition were performed on paired sera collected at admission and at discharge, and systematically tested for both dengue and Japanese

<figure><img src="image_3.png" /></figure>
Encephalitis virus, another flavivirus endemic in Cambodia, as described previously [30]. Virus isolation was carried out on earliest serum samples by inoculating permissive C6/36 and VERO E6 cells, followed by serotype-specific immunofluorescence [30]. Viral RNA was detected in specimens collected at early stage of the disease using a nested RT-PCR [31]. Primary or anamnetic "secondary" antibody response indicating previous infections by dengue viruses was determined from paired serum samples by hemagglutination-inhibition assay. Interpretation followed WHO recommendations [29].

## Blood sample preparation

Whole blood samples (2.5 ml) were collected on PAXgene™ Tubes (PreAnalytiX™) further stored at -80°C, before being sent to France in dry ice. Extraction of series of 24 matched samples (DF, DHF and DSS) was done using PAXgene™ Blood RNA kits (PreAnalytiX™) rapidly after collection. Purified total RNAs kept at -80°C were processed for hybridization on genome-wide DNA microarrays within one month.

## cRNAs preparation and microarrays hybridization

All RNAs were checked for integrity using the 2100 BioAnalyzer (Agilent Technologies) and quantified using a ND-1000 spectrophotometer (NanoDrop Technologies). Cyanine-3-labeled cRNA was generated from 0.3 µg of RNA using the One-Color Low RNA Input Linear Amplification kit (Agilent) according to the manufacturer's instructions, followed by purification on RNAeasy column (QIAGEN). All amplified cRNAs were checked for dye incorporation, cRNA yield and amplification profile. Only those fitting all quality criteria were fragmented for further hybridization on microarrays. Samples from DF, DHF and DSS patients were then carefully matched and hybridized onto Agilent Whole Human Genome (4×44K) Oligo Microarrays (G4112F). Microarrays were scanned using an Agilent DNA microarray scanner G2505B.

## Microarray data analysis

All microarray data is MIAME compliant and the raw and normalized data have been deposited in the MIAME compliant database Gene Expression Omnibus [32] (GEO Series accession number GSE17924 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc = GSE17924).

Individual microarray quality was evaluated based on QC report, pair-wise MA-plots, and box plots. Intra-array normalization of raw signals from the 48 microarrays was done using Feature Extraction software 9.1.3.1 (Agilent). Microarrays normalized data were further exported into the Limma package [33], for inter-array normalization using the quantile method [34].

Statistical analysis were performed using the TIGR MeV (MultiExperiment Viewer) v 4.4 software (http://www.tm4.org/mev.html), and the GeneANOVA program [35]. Multi-way ANOVA model was implemented: first, to identify differentially regulated genes when accounting for the multiple sources of variation in the microrray experiment; second, to evaluate the effect of the main variable, disease phenotype, relatively to that of other putative confounding variables such age, gender, duration of illness or microarray technical variability (independent extractions or hybridizations). Local ANOVA further determined the contribution of each covariate on the expression level of each gene. Multiple test correction was further carried out using the false discovery rate (FDR) method [36]. Cluster [37] and Tree View softwares [38] were used for unsupervised hierarchical clustering. Iterative SVM (Support Vector Machine) method

associated with leave-one-out cross-validation [39] was used to assess the robustness of DSS-gene signature.

## Real-time PCR validation of genes over and under-expressed in DSS patients

Briefly, total RNA extracted from whole blood samples was reverse-transcribed using the High Capacity cDNA RT kit (Applied Biosystems Inc) and random primers. Real-time PCR were carried out using the FastStart Universal Probe Master (ROX) (Roche) and real-time PCR primers designed using the Universal Probe Library (UPL) Assay Design Center (Roche). Amplification products were run on an ABI-PRISM 7900HT (Applied Biosystems). Cycle threshold Ct values were automatically calculated and value obtained for each gene amplified was normalized by subtracting the Ct corresponding to amplification of the HPRT1 gene (ΔCt) for the same sample. Correlation between ΔCt values obtained by real-time PCR and corresponding expression values from microarrays was estimated using Spearman correlation coefficient.

## Comprehensive overview of functional patterns altered during DSS

Bio-informatics-based analysis using the demonstration version 7.1 of Ingenuity Pathway Analysis software (IPA; Ingenuity® Systems, www.ingenuity.com) associated with manual and literature-based analysis was carried out to identify the most relevant functional processes associated with the identified DSS gene signature. This was done by combining most informative canonical pathways identified using IPA, genes having the strongest association with the disease phenotype based on ANOVA analysis, and similarities to molecular patterns altered in other systemic inflammatory processes associated with endothelial dysfunction.

## Results

### Patient characteristics

To identify gene patterns specifically altered in DSS patients, we compared three groups of carefully matched paediatric patients representing the main clinical forms of symptomatic dengue infections DF (n = 16), DHF (n = 13) and DSS (n = 19), according to the 1997 WHO classification criteria of dengue severity [29]. Altogether, DF, DHF and DSS represent different subtypes of the disease phenotype variable, further considered in this study. The clinical characteristics and values of haematological parameters are presented in table 1 (median values from each patient group) and table S1 (individual values from each of the 48 patients included). Supportive treatment provided to DSS patients are mentioned. As indicated, DSS children had significant lower relative neutrophil counts (median values / DF: 3900; DHF: 3950; DSS: 2500; p-value: 0.03; Kruskal Wallis test).

### Unsupervised hierarchical clustering discriminates DSS children from DF/DHF ones, revealing a DSS-gene signature

Since microarray data analysis can be affected by a number of bias [40], we put a particular care on study design and analysed data from the 48 normalized microarrays using multi-way analysis of variance (ANOVA) [41]. Indeed ANOVA evaluates the statistical probability (p-value), for each individual gene, that a difference in expression between the three patient groups could have been observed by chance. This allows revealing genes that show even small but highly significant changes in expression, regarding the studied phenotype.

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 1.** Clinical and biological characteristics of DF, DHF and DSS patient groups at the time of hospital admission.</caption><thead><tr><th></th><th>DF (n=16)</th><th>DHF (n=13)</th><th>DSS (n=19)</th></tr></thead><tbody><tr><td><strong>Patients characteristics</strong></td><td></td><td></td><td></td></tr><tr><td>gender, male n (%)</td><td>7 (43%)</td><td>4 (31%)</td><td>7 (37%)</td></tr><tr><td>age, median (IQR), years</td><td>8 (4-9)</td><td>7 (5-8)</td><td>8 (7-9)</td></tr><tr><td>weight, median (IQR), kg</td><td>18 (13-20)</td><td>15 (14-18)</td><td>19 (15-23)</td></tr><tr><td>hospital admission, median (IQR), day after onset of fever (D0)</td><td>2 (1-3)</td><td>2 (2-3)</td><td>4 (3-4)</td></tr><tr><td><strong>Dengue status</strong></td><td></td><td></td><td></td></tr><tr><td>viral serotype, n (DENV-1/DENV-2/DENV-3/DENV-4/unknown)</td><td>4/2/8/1/1</td><td>1/1/10/1/0</td><td>1/1/10/0/7</td></tr><tr><td>immunological status, secondary infections, n (%)</td><td>14 (88%)</td><td>12 (92%)</td><td>18 (95%)</td></tr><tr><td><strong>Clinical manifestations</strong></td><td></td><td></td><td></td></tr><tr><td>tourniquet test (pos/neg/not done) (%)</td><td>56%/44%/0%</td><td>54%/38%/8%</td><td>37%/32%/31%</td></tr><tr><td>hepatomegaly, n (%)</td><td>3 (19%)</td><td>6 (46%)</td><td>17 (89%)</td></tr><tr><td>gastro-intestinal bleeding, gingivorragy, hematemesis, melena, n (%)</td><td>0</td><td>1 (8%)</td><td>6 (32%)</td></tr><tr><td><strong>Blood pressure</strong></td><td></td><td></td><td></td></tr><tr><td>heart frequency, median (IQR), pulse per minute</td><td>113 (100-124) (n=14<sup>a</sup>)</td><td>120 (112-120)</td><td>Not perceptible (n=15<sup>a</sup>)</td></tr><tr><td>pulse pressure, median (IQR), mm Hg</td><td>40 (30-45)</td><td>30 (30-40)</td><td>15 (10-20) (n=15<sup>a</sup>)</td></tr><tr><td><strong>Haematological parameters</strong></td><td></td><td></td><td></td></tr><tr><td>thrombocytopenia (platelet count &lt;100000/mm<sup>3</sup>), %</td><td>15% (n=14<sup>a</sup>)</td><td>55% (n=11<sup>a</sup>)</td><td>94% (n=17<sup>a</sup>)</td></tr><tr><td>hematocrit, median (IQR), %</td><td>36.5 (35-39) (n=14<sup>a</sup>)</td><td>39.75 (38-42) (n=12<sup>a</sup>)</td><td>42.5 (38-45) (n=18<sup>a</sup>)</td></tr><tr><td>hemoconcentration (hematocrit &gt;20%), n (%)</td><td>1 (6%)</td><td>3 (23%)</td><td>17 (89%)</td></tr><tr><td>white blood cells, median (IQR), number/mm<sup>3</sup></td><td>6600 (5500-9900) (n=13<sup>a</sup>)</td><td>6450 (6200-7400) (n=10<sup>a</sup>)</td><td>6900 (4800-6900) (n=17<sup>a</sup>)</td></tr><tr><td>neutrophils, median (IQR), number/mm<sup>3</sup></td><td>3900 (2900-7600) (n=13<sup>a</sup>)</td><td>3950 (3500-4200) (n=10<sup>a</sup>)</td><td>2500 (2200-3800)(n=17<sup>a</sup>)</td></tr><tr><td>lymphocytes, median (IQR), number/mm<sup>3</sup></td><td>1600 (1400-2100) (n=13<sup>a</sup>)</td><td>1850 (1500-1900) (n=10<sup>a</sup>)</td><td>2200 (1500-3200) (n=17<sup>a</sup>)</td></tr><tr><td><strong>Supportive medical care</strong></td><td></td><td></td><td></td></tr><tr><td>oxygen supplementation, n (%)</td><td>0</td><td>0</td><td>15 (79%)</td></tr><tr><td>perfusion of colloid (dextran 40), n (%)</td><td>0</td><td>0</td><td>14 (74%)</td></tr><tr><td>perfusion of human plasma, n (%)</td><td>0</td><td>0</td><td>8 (42%)</td></tr></tbody></table>

DENV, dengue virus; DF, dengue fever; DHF, dengue hemorrhagic fever; DSS, dengue shock syndrome; IQR, interquartile range; n, number.

<sup>a</sup>n=x: with x: number of patients for which the data is available.

doi:10.1371/journal.pone.0011671.t001

Based on ANOVA analysis, lists of genes differentially expressed between DF, DHF and DSS groups were generated using different false discovery rate (FDR) ranging from 0.05 up to 10%. Indeed, low FDR provide more stringent statistical filter while they reduce the number and thus the enrichment of genes differentially expressed. At the opposite, higher FDR, while statistically accepting a higher number of false positive genes, also provide larger and enriched gene lists that should be more informative when searching to identify molecular pathways. Based on this rationale, we chose to work using the gene list generated at FDR 10 after we verified by a different statistical method currently used for the analysis of microarrays data, SAM (Significant Analysis of Microarray) [33], that most significant genes were commonly found by the two types of analyses (data not shown). The gene list generated at FDR10 included 2959 genes differentially expressed between DF, DHF and DSS patients groups (Table S2).

The biological relevance of those differentially expressed genes was assessed using local ANOVA that allows evaluating the contribution of the main variable, disease phenotype, and that of other putative confounding variables related to patients (age, gender, day of blood sampling, viral serotype) and to technical steps (effect of independent RNA extractions, amplifications and hybridization) on variations of expression levels of those 2959 genes. This confirmed that the disease

phenotype strongly influenced the variations of expression of the 2959 genes differentially expressed between the three patient groups, reinforcing the biological significance of this set of genes (Table S2).

Unsupervised hierarchical clustering based on the 2959 gene signature identified was then applied to the 48 children gene expression profiles. This allows clustering the patients whose gene expression profiles are the more similar independently of their disease phenotype subtype. As a result, the 48 patients expression profiles were organized in two major subsets (Figure 1): subset 1 (first dendrogram branch) includes both DF and DHF patients without distinction; subset 2 (second dendrogram branch) encompasses a sub-group (2a) of DF and DHF patients, and a distinct sub-group (2b) including 17 out of the 19 DSS patients, whatever they received or not plasma infusion, revealing a “DSS-gene signature” common to most DSS patients. Some few patients clustered however in unexpected subsets: two DSS patients (PL005, PL101) had gene expression profiles closer to those of the DF/DHF 2a subset, while one DF patient (PL064) and three DHF patients (PL037, PL058, PL070) gene expression profiles clustered within the DSS 2b subset.

We confirmed the robustness of the DSS-gene signature using the iterative Support Vector Machine (SVM) classifier learning method [39], which reclassified all the 19 DSS patients together.

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 1.** Unsupervised hierarchical clustering of whole blood cells expression profiles from the 48 dengue-infected children.</figcaption></figure>

The clustering is based on the 2959 gene list (3515 clones, detailed in Table S2) discriminating dengue fever (DF), dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS) patients. Each row represents a single transcript and each column represents a patient's sample. Color scale indicates the range of gene expression: black indicates median expression level, red greater expression, green lower expression. The 2 patient subsets identified are indicated. PLxxx, code relative to one patient. Black star: DSS patient sampled 3 days after shock. Orange star: Patients who received perfusion of human plasma before collection of blood samples.

doi:10.1371/journal.pone.0011671.g001

To validate microarray data, we carried out real-time RT-PCR focusing on nine genes strongly associated with the DSS-gene signature, using 15 patients samples (five from each disease phenotype subtype: DF, DHF and DSS). Results obtained strongly correlate microarray data (Figure S1).

## DSS-gene signature analysis identifies a diversity of genes and canonical molecular pathways related to immunity, inflammation and host metabolism

Filtering genes from those having the highest to the lowest statistical association with the disease phenotype variable (Table S2), relying on results from multi-way and local ANOVA revealed that the individual genes having the strongest association with the DSS phenotype subtype are, for a large part, related to innate immunity, inflammation and host lipid metabolism, a finding confirmed when the whole 2959 genes of the DSS-gene signature were processed through the knowledge-based IPA software.

Indeed, IPA analysis identified that 163 canonical pathways were significantly associated with those genes (data not shown) with a large proportion of immune-related pathways in the first top 30 (Figure 2). In particular, several under-expressed but partially redundant signaling canonical pathways related to T lymphocyte activation were identified, of which the T cell receptor (TCR) signaling pathway (Figure 3), which has the strongest association with the DSS-gene signature. Interestingly, a number of metabolic pathways, and particularly of lipid-signaling pathways, were significantly represented among the 163 DSS-related canonical pathways.

When comparing our results to those of colleagues who reported gene or protein signatures associated with DSS, we identified some transcripts encoding proteins considered putative markers of severe dengue. This includes non exhaustively the acute phase pentraxin-related protein PTX3 [15], the anti-inflammatory IL-10 [11] or the pro-inflammatory IL-18 [12] cytokine transcripts that have increased abundance in the DSS-gene signature, while having intermediate to low statistical association with the disease phenotype variable according to the multi-way ANOVA (Table S2).

IFN type I-related transcripts, of which abundance was shown to be decreased in DSS patients by others [27,28,42], represented only a limited number of genes associated with the DSS-gene signature. This was confirmed by the IPA analysis that did not identify IFN type I-related pathways among those strongly associated with the DSS-gene signature (Figure 2).

## DSS is associated with impaired expression of T and NK cell-related genes but increased expression of anti-inflammatory and repair/remodeling transcriptional responses

Integrative analysis of the most significant individual genes and canonical pathways extended the finding that a large and diverse set of genes related to T but also to NK lymphocyte activity is

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 2.** Top 30 canonical pathways identified from the DSS-gene signature using Ingenuity Pathway Analysis software. The significance of the association between data set and canonical pathway was estimated by the p-value (Fischer's exact test; left axis) and the ratio (right axis) of genes that maps to each canonical pathway.
doi:10.1371/journal.pone.0011671.g002</figcaption></figure>

under-expressed in DSS patients compared to DF and DHS counterparts (Table 2; non exhaustive list; individual p-values available in Table S2; and Figure 3). Those genes are critical to a number of T and NK-cell functions, including T and NK cell differentiation, receptor signaling, activation and proliferation, cytotoxic functions or recruitment of lymphocytes to peripheral tissues. Since lymphocyte counts did not differ between the DF, DHS and DSS children (p = 0.428; Kruskal Wallis test), we searched whether genes encoding factors regulating negatively T and NK functions were over-expressed in the DSS-gene signature. We identified that the two genes having the strongest association with the disease phenotype variable, encode two major immunomodulatory factors, the microsomal prostaglandin E synthase (PTGES/Agilent clone number A-24-P478940) and the complement regulatory protein CR1g/VSIG4, considered potent negative regulators of T and NK cell responses [43–45]. The decreased abundance of *Nfkb* signal transduction-related transcripts (Table 2), already reported in DSS patients by others [46], might be related to impaired expression of T and NK cell-related genes.

Our analysis also revealed that DSS whole blood cells from children over-expressed an enriched pattern of anti-inflammatory and repair/tissue remodeling genes (Table 3; non exhaustive list; individual p-values available in Table S2). Over-expressed anti-inflammatory genes identified encode molecules with diverse functions: the anti-inflammatory cytokine IL-10, a putative marker of severe dengue [11], serine proteases and metalloprotease inhibitors, IL-1β cytokine decoy receptor, free heme scavenger molecules, or complement-regulating receptors. Repair and remodeling genes over-expressed in the DSS-gene signature also

encode a diversity of molecules: the MMP-9, a matrix metalloprotease with key role in tissue remodeling and a candidate to dengue plasma leakage [47], the extracellular matrix molecules fibronectin, versican and collagens, the angiogenin and VEGF [17] endothelial agonists as well as the arginase 1 repair enzyme, which competes with the endothelial NOS (NOS3) for L-arginine biodisponibility [48].

Thus, DSS children whole blood cells have a global decreased abundance of T and NK cell-related transcripts but an increased abundance of anti-inflammatory and repair/remodeling transcripts at the time of cardiovascular decompensation.

**Pro-inflammatory innate defense and host-lipid
metabolism-related transcriptional responses are
activated in DSS children**

When searching for pro-inflammatory gene patterns that may be relevant to DSS pathophysiology and particularly to systemic inflammation and vascular dysfunction, we identified three major pro-inflammatory gene patterns. Interestingly, all are related to innate defense and host lipid metabolism, and considered major pathogenic mechanisms in other systemic inflammatory diseases.

As shown in table 4 (non exhaustive list; individual p-values available in Table S2), the first one is defined by a set of over-expressed genes strongly associated with the disease phenotype subtype DSS that encode highly pro-inflammatory microbicidal peptides and enzymes. This pattern includes non exhaustively the alpha defensins DEFA1, DEFA3 and DEFA4, the cathelicidin (CAMP) and lactoferrin (LTF) peptides, the neutrophil enzymes

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 3.** T Cell Receptor Signaling canonical pathway from Ingenuity Pathway Analysis. Genes in green and red are respectively under- and over-expressed in the DSS-gene signature. Genes in white are other genes present in the canonical pathway but absent from the DSS-gene signature. DSS: Dengue Shock Syndrome.
doi:10.1371/journal.pone.0011671.g003</figcaption></figure>

myeloperoxidase (MPO), neutrophil RNASE2, RNASE3, cathepsin G and neutrophil elastase (ELANE). Transcripts encoding the potent pro-inflammatory calgranulins proteins S100A8/9 and S100A12, characteristic of granulocyte neutrophil activity [51] and involved in a diversity of inflammatory diseases [56], as well as the granulocyte-related metalloprotease MMP8 are also over-expressed. Increased abundance of those transcripts cannot be explained by increase in granulocyte count, since DSS patients have lower relative granulocyte counts than DF and DHF counterparts (median values/DF: 3900; DHF: 3950; DSS: 2500; p-value: 0.03; Kruskal Wallis test), thus reflecting more likely cellular activation. Altogether those results show that a transcriptional pattern of innate defense genes is activated in the whole blood of DSS children.

The second pro-inflammatory gene pattern identified is typical
of altered homeostasis of cholesterol in monocytes/macrophages

that characterises inflammatory lipid-laden monocytes/macro-
phages (lipid-laden Mo/Mac), a subtype of foam cells initiating
vascular lesions in metabolic inflammatory diseases [57–59]
(Table 5; non exhaustive list; individual p-values available in
Table S2). Since the *PPARG* gene, which encodes a nuclear-
lipid receptor involved in lipid signaling and lipid-homeostasis
in inflammatory lipid-laden Mo/Mac [65], has a very strong
association with the dengue disease phenotype, we searched
whether other genes involved in cholesterol homeostasis in Mo/
Mac had altered expression in the DSS-gene signature.
Remarkably, we found a large lipid-laden Mo/Mac-related gene
expression pattern characterized in particular by a higher
abundance of transcripts encoding the key scavenger receptors
of modified low density lipoproteins OLR-1, CD36 and MSR1,
but a decreased abundance of transcripts encoding critical
cholesterol transporters such as the NPC1 [66] or the ABCA1-

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 2.** T lymphocytes and NK cells-related genes present in the DSS-gene signature.</caption><thead><tr><th>Function</th><th>Genes</th><th>P-value</th><th>Var.<sup>a</sup> (%)</th></tr></thead><tbody><tr><td>Th1 differentiation</td><td>RUNX3, STAT4, TBX21</td><td>&lt;0.00001 to 0.00242</td><td>25 to 42</td></tr><tr><td>Th2 differentiation</td><td>GATA3, STAT5A</td><td>0.00003 to 0.00225</td><td>17 to 32</td></tr><tr><td>Cytotoxic T lymphocyte functions</td><td>CTSW, PRF1</td><td>0.00005 to 0.00231</td><td>21 to 33</td></tr><tr><td>T lymphocyte activation</td><td>IL2RB, IL2RG</td><td>0.00014 to 0.00039</td><td>29 to 35</td></tr><tr><td>Cooperation with antigen-presenting cells</td><td>CD40LG</td><td>0.00105</td><td>21</td></tr><tr><td>Recruitment and interaction of T lymphocytes with endothelium</td><td>ITGAL, XCL1, XCL2</td><td>&lt;0.00001 to 0.00214</td><td>20 to 33</td></tr><tr><td>Inhibitory NK cell receptors</td><td>KLRD1</td><td>0.00001</td><td>31</td></tr><tr><td>Activating NK cell receptors</td><td>NCR1, NCR3, CD160</td><td>&lt;0.00001 to 0.00069</td><td>28 to 39</td></tr><tr><td>Cytotoxic molecules</td><td>GZMM</td><td>&lt;0.00001</td><td>32</td></tr><tr><td>Receptors for NK cells homing to peripheral tissues</td><td>S1PR5</td><td>&lt;0.00001</td><td>48</td></tr><tr><td>Differentiation factors of NK cells</td><td>FLT3LG, <b>IL15</b>, <b>IL17C</b>, KITLG</td><td>0.00088 to 0.00774</td><td>13 to 21</td></tr><tr><td>Suppression of T lymphocytes and NK cells response</td><td><b>PTGES</b>, <b>VSIG4</b></td><td>&lt;0.00001</td><td>60 to 63</td></tr><tr><td>NFKB-related genes</td><td><b>IRAK3</b>, TNIK, RELA, NFKBIB, TRAF1, TRAF2, TRAF6</td><td>0.00001 to 0.00506</td><td>9 to 31</td></tr></tbody></table>

HUGO gene names are indicated. When genes were represented by several clones on the microarray, p-value and variance medians were calculated. Genes in regular and bold are respectively under- and over-expressed in dengue shock syndrome patients.

like ABCA10 [86], which regulate the efflux of modified cholesterol from Mo/Mac. Other lipid-laden cells-related genes also have altered expression in the DSS-gene signature. In particular, the *PPARA* gene that negatively regulates the formation of lipid-laden Mo/Mac [87] has decreased abundance in DSS patients. At the opposite, transcripts encoding the chitinase 1, a marker of pro-inflammatory lipid-laden Mo/Mac [73], and the FABP4, SOCS6, RETN and IRS2 proteins involved in lipid-laden Mo/Mac-induced insulin-resistance and compensatory response [68–72], have all increased abundance, also strongly supporting a biological signature of foam cells. Interestingly, the *PCSK9* transcript, which encodes a secreted

protein that decreases the recycling of LDL to the liver by inducing the degradation of liver LDL receptors [74], is also over-expressed in the DSS signature and highly associated with the disease phenotype.

Thus, a gene expression pattern similar to that characterizing lipid-laden monocytes, is activated in the whole blood cells of DSS children at the time of cardiovascular decompensation.

The third pro-inflammatory gene pattern associated with the DSS-gene signature is characteristic of the metabolic pro-inflammatory arachidonic-acid pathway, one of the lipid metabolic pathways identified through IPA. As shown in table 5, the gene encoding the upstream cytosolic phospholipase

<table><caption>**Table 3.** Anti-inflammatory, tissue remodeling and repair genes present in the DSS-gene signature.</caption><thead><tr><th>Function</th><th>Gene Symbol</th><th>P-value</th><th>Var.<sup>a</sup> (%)</th></tr></thead><tbody><tr><td>Anti-inflammatory genes</td><td></td><td></td><td></td></tr><tr><td>immunoregulatory molecules</td><td><b>IL10</b></td><td>0.00430</td><td>20</td></tr><tr><td>anti-proteases</td><td><b>SERPINB2</b>, <b>SERPINB8</b>, <b>SERPINB10</b>, <b>SLPI</b></td><td>&lt;0.00001 to 0.00081</td><td>19 to 49</td></tr><tr><td>metalloproteinase inhibitor</td><td><b>TIMP1</b></td><td>0.00183</td><td>19</td></tr><tr><td>decoy receptor</td><td><b>IL1R2</b></td><td>0.00077</td><td>30</td></tr><tr><td>free-heme scavenger molecules</td><td><b>CD163</b>, <b>HP</b>, <b>HMOX1</b></td><td>&lt;0.00001 to 0.00064</td><td>26 to 46</td></tr><tr><td>complement regulatory molecules</td><td><b>CD55</b>, <b>VSIG4</b></td><td>&lt;0.00001 to 0.00096</td><td>24 to 60</td></tr><tr><td>Tissue remodeling and repair genes</td><td></td><td></td><td></td></tr><tr><td>metallopeptidase</td><td><b>MMP9</b></td><td>0.00001</td><td>33</td></tr><tr><td>extracellular matrix components</td><td><b>COL1A2<sup>b</sup></b>, <b>COL8A2<sup>b</sup></b>, <b>COL14A1<sup>b</sup></b>, <b>COL17A1<sup>b</sup></b>, <b>FN1<sup>b</sup></b>, <b>SDC1<sup>b</sup></b>, <b>VCAN<sup>b</sup></b></td><td>&lt;0.00001 to 0.00309</td><td>18 to 34</td></tr><tr><td>pro-angiogenic factors</td><td><b>ANG</b>, <b>VEGFA</b></td><td>0.00004 to 0.00236</td><td>25 to 30</td></tr><tr><td>others</td><td><b>ARG1</b>, <b>NOS3</b></td><td>&lt;0.00001 to 0.00054</td><td>18 to 44</td></tr></tbody></table>

HUGO gene names are indicated. When genes were represented by several clones on the microarray, p-value and variance medians were calculated. Genes in regular and bold are respectively under- and over-expressed in dengue shock syndrome patients.

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

<footnote><sup>a</sup>percentage of variance associated to disease phenotype.</footnote>

<footnote>doi:10.1371/journal.pone.0011671.t002</footnote>

<footnote><sup>a</sup>percentage of variance associated to disease phenotype.</footnote>

<footnote><sup>b</sup>Danger-associated molecular pattern (DAMP) activity.</footnote>

<footnote>doi:10.1371/journal.pone.0011671.t003</footnote>
<table><caption>**Table 4.** Pro-inflammatory innate immunity-related genes present in the DSS-gene signature.</caption><thead><tr><th>Function</th><th>Gene Symbol</th><th>P-value</th><th>Var.<sup>a</sup> (%)</th><th>Main cellular origin</th><th>Ref</th></tr></thead><tbody><tr><td rowspan="3">microbicidal peptides</td><td>DEFA1<sup>b</sup>, DEFA3<sup>b</sup>, DEFA4<sup>b</sup></td><td>&lt;0.00001 to 0.00007</td><td>0.25 to 0.44</td><td>PMN neutro, EpC</td><td>[49,50]</td></tr><tr><td>CAMP<sup>b</sup></td><td>&lt;0.00001</td><td>0.34</td><td>PMN neutro, Mo, mast cells, EpC</td><td>[49,50]</td></tr><tr><td>LTF<sup>b</sup></td><td>&lt;0.00001</td><td>0.41</td><td>PMN neutro, inflamed EpC</td><td>[50]</td></tr><tr><td rowspan="2">calgranulin proteins</td><td>S100A8<sup>b</sup>, S100A9<sup>b</sup></td><td>&lt;0.00001 to 0.00014</td><td>0.18 to 0.38</td><td>PMN neutro, Mo/Mac</td><td>[51,52]</td></tr><tr><td>S100A12<sup>b</sup></td><td>&lt;0.00001</td><td>0.33</td><td>PMN neutro</td><td>[51]</td></tr><tr><td rowspan="7">granulocyte enzymes</td><td>RNASE2<sup>b</sup></td><td>0.00017</td><td>0.25</td><td>Mo/Mac, Eo, EpC, PMN neutro</td><td>[53]</td></tr><tr><td>MPO<sup>b</sup></td><td>0.00024</td><td>0.25</td><td>PMN neutro, Mo, subtypes of tissue Mac</td><td>[50]</td></tr><tr><td>RNASE3<sup>b</sup></td><td>&lt;0.00001</td><td>0.29</td><td>Eo, Mo, PMN neutro</td><td>[54]</td></tr><tr><td>MMP8</td><td>&lt;0.00001</td><td>0.49</td><td>PMN neutro</td><td>[50]</td></tr><tr><td>CTSG</td><td>&lt;0.00001</td><td>0.36</td><td>PMN neutro</td><td>[50]</td></tr><tr><td>ELANE</td><td>&lt;0.00001</td><td>0.39</td><td>PMN neutro</td><td>[50]</td></tr><tr><td>IL-18</td><td>0.00052</td><td>0.21</td><td>Kupffer cells, activated Mac, Mo, DC, EpC</td><td>[55]</td></tr><tr><td>pro-inflammatory cytokines and related molecules</td><td>IL18BP</td><td>0.00710</td><td>0.20</td><td>T cells, peripheral blood leukocytes, EC</td><td>[55]</td></tr></tbody></table>

HUGO gene names are indicated. When genes were represented by several clones on the microarray, p-value and variance medians were calculated. Genes in regular and bold are respectively under- and over-expressed in dengue shock syndrome patients. DC, dendritic cell; EC, endothelial cell; Eo, eosinophil; EpC, epithelial cell; Mac, macrophage, Mo, monocyte; PMN neutro, polymorphonuclear neutrophil; RAGE, receptor for advanced glycation end products.

PLA2G4A, which is the initial rate-limiting enzyme that cleaves membrane phospholipids [78], is over-expressed. Similarly, most downstream key enzymes from the COX-2 and 5-LOX sub-pathways involved in the final synthesis, conversion and transport of inflammatory eicosanoids lipid mediators, are over-expressed. In particular, the transcript encoding the inducible microsomal prostaglandin E synthase PTGES that catalyzes the conversion of prostaglandin PGH2 to PGE2 in the COX-2 sub-pathway, and thought to play a pathogenic role in a number of inflammatory processes [88] is significantly increased and has the highest statistical association with the disease phenotype (62% of gene variance explained by the disease phenotype according to multi-way ANOVA). At the opposite the *PTGDS* transcript, which encodes the anti-inflammatory prostaglandin D2 synthase, has decreased abundance, a finding already reported in metabolic inflammatory processes [89]. Increased abundance of the transcript encoding the LTA4H enzyme that converts the LTA4 leukotriene to LTB4 reflects the activation of the 5-LOX sub-pathway. Finally, transcripts encoding the oxidative enzymes ALOX15B lipoxygenase [84] and cytochrome P450 epoxygenase family members [85], involved in the arachidonic acid metabolic pathway, are also significantly increased in the DSS-gene signature, also reflecting activation of those sub-pathways during DSS.

Thus, a transcriptional signature related to the lipid-related metabolic arachidonic acid pathway is activated in the whole blood cells of DSS children at the time of cardiovascular decompensation.

Discussion

Numerous studies have addressed the pathophysiology of DSS, the more frequent and severe complication of dengue infections. Despite important findings, only partial understanding of the cellular and molecular processes that may support this

life-threatening syndrome has been obtained, and we still lack a comprehensive overview of the complete figure of alterations that contribute to – or reflect – the setting-up of the shock syndrome. This could allow the improvement of patients’ management and treatment, a major challenge for clinicians. We designed a study aimed at analysing the quasi-global transcriptome of whole blood cells from dengue paediatric patients, looking at every modification that could make sense to the understanding of the pathogenic process. The capacity of such an exhaustive approach to identify relevant host responses, of which unsuspected pathways has been demonstrated in other systemic inflammatory syndromes such as human sepsis or post-trauma sterile SIRS [23–25]. We compared the transcriptome of blood cells from DSS paediatric patients at time of shock to those of children classified as DF or DHF grades I/II [29] matched for important variables such as age, gender, immune status towards dengue infection (primary or secondary infection) and time of disease evolution after onset of fever. Our study has produced significant results, further discussed in the context of DSS pathophysiology.

First, we identify a transcriptional signature of the DSS, differentiating DSS from the other forms of dengue infection and characterizing DSS as a unique and specific entity. Giving particular attention to study design and statistical analysis, we identify a large and robust gene expression profile of 2959 genes that discriminates DSS paediatric patients from other dengue patients, DF or DHF, who did not progress to shock, whatever the supportive treatment they received. Importantly, DSS children clustered together whatever they were considered as having primary or secondary dengue infection, while secondary infections represented the majority of DF, DHF and DSS children recruited (see table S1), as expected in hyper-endemic areas. The robustness of the DSS-associated gene signature was established by showing that the disease phenotype variable significantly affected expression levels of all the genes identified (multi-way ANOVA) and

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

<footnote><sup>a</sup>percentage of variance associated to disease phenotype.</footnote>

<footnote><sup>b</sup>Danger-associated molecular pattern (DAMP) activity.</footnote>

<footnote>doi:10.1371/journal.pone.0011671.t004</footnote>
<table><caption>**Table 5.** Pro-inflammatory lipid-related genes present in the DSS-gene signature.</caption><thead><tr><th>Function</th><th>Gene Symbol</th><th>P-value</th><th>Var.<sup>a</sup> (%)</th><th>Disease</th><th>Ref</th></tr></thead><tbody><tr><td colspan="6"><strong>Lipid-laden Mo/Mac-related genes</strong></td></tr><tr><td>scavenger receptors of modified LDL in Mo/Mac</td><td>OLR1, CD36, MSR1</td><td>&lt;0.00001 to 0.00013</td><td>0.21 to 0.32</td><td>metabolic diseases</td><td>[60–64]</td></tr><tr><td>lipid nuclear receptor/signalisation by lipids</td><td>PPARG, PPARA</td><td>0.00007 to 0.00732</td><td>0.21 to 0.34</td><td>metabolic diseases</td><td>[65]</td></tr><tr><td>efflux of modified cholesterol from Mo/Mac</td><td>NPC1</td><td>0.00005</td><td>0.32</td><td>Niemann-Pick disease, atherosclerosis</td><td>[66,67]</td></tr><tr><td></td><td>ABCA10</td><td>0.00016</td><td>0.14</td><td>none</td><td></td></tr><tr><td>migrating Mo/resident Mac chemokine receptors</td><td>CCR2, CX3CR1</td><td>0.00001 to 0.00099</td><td>0.22 to 0.40</td><td>atherosclerosis</td><td>[57]</td></tr><tr><td>other lipid-laden-related Mo/Mac genes</td><td>FABP4, SOCS6, RETN, IRS2</td><td>&lt;0.00001 to 0.00092</td><td>0.20 to 0.26</td><td>metabolic diseases</td><td>[68–72]</td></tr><tr><td></td><td>CHIT1</td><td>&lt;0.00001</td><td>0.48</td><td>Gaucher's disease, atherosclerosis</td><td>[73]</td></tr><tr><td></td><td>PCSK9</td><td>0.00001</td><td>0.42</td><td>familial hypercholesterolemia</td><td>[74]</td></tr><tr><td></td><td>SPP1</td><td>&lt;0.00001</td><td>0.49</td><td>metabolic and inflammatory diseases</td><td>[75,76]</td></tr><tr><td>anti-oxydant enzymes</td><td>LCAT, PAFAH2</td><td>0.00196 to 0.00461</td><td>21 to 26</td><td>metabolic diseases</td><td>[77]</td></tr><tr><td colspan="6"><strong>Arachidonic acid pathway-related genes</strong></td></tr><tr><td>phospholipase</td><td>PLA2G4A</td><td>0.00003</td><td>0.21</td><td>rheumatoid arthritis</td><td>[78]</td></tr><tr><td>eicosanoid synthesis enzymes</td><td>PTGES, LTA4H, PTGDS, TBXAS1, PTGDR</td><td>&lt;0.00001 to 0.00123</td><td>0.22 to 0.63</td><td>metabolic and inflammatory diseases, asthma, cancer</td><td>[79–82]</td></tr><tr><td>leukotrienes conversion enzyme</td><td>MGST2</td><td>0.00003</td><td>0.32</td><td>none</td><td></td></tr><tr><td>leukotriene transporter</td><td>SLCO2B1</td><td>0.00010</td><td>0.31</td><td>asthma</td><td>[83]</td></tr><tr><td>lipid oxidation</td><td>ALOX15B</td><td>0.00011</td><td>0.33</td><td>atherosclerosis</td><td>[84]</td></tr><tr><td>cytochrome P450 superfamily enzymes</td><td>CYP1B1, CYP2U1, CYP51A1</td><td>&lt;0.00001 to 0.00686</td><td>10 to 32</td><td>Vascular inflammation</td><td>[85]</td></tr></tbody></table>

HUGO gene names are indicated. When genes were represented by several clones on the microarray, p-value and variance medians were calculated. Genes in regular and bold are respectively under- and over-expressed in dengue shock syndrome patients.

demonstrating the classifying capability of this gene signature using unsupervised hierarchical clustering and SVM/leave-one-out methods [35,39,40]. Based on unsupervised hierarchical clustering, DHF grades I/II patients' expression profiles appear very close or indistinguishable from those of DF patients at the same time of disease evolution, while they group into two heterogeneous sub-groups (1 and 2b, Figure 1), of which significance should be investigated. Altogether, the present results highlight the inadequacy of the 1997 WHO classification of dengue clinical forms [19–21], that considers DF and DHF grades I/II as two separate disease phenotypes, and support the recently proposed classification [90].

Two important questions arise about the DSS-associated transcriptional profile: are the observed modifications of genes expression the cause or the consequence of the pathology, and could these modifications have a predictive value? We cannot definitively answer these questions from the present study, in part because blood samples were collected at the onset of shock (14 out of the 19 DSS patients) or after (5 patients). Functional study of each individual pathway will be required to fully understand the role of each gene in a complex network of molecular interactions.

The ability of some genes transcripts or genes products to accurately predict progression to DSS should be evaluated by multivariate regression models [91] using blood samples collected before the onset of shock, while this proves to be difficult in the context of dengue outbreaks [28]. In the present study, we chose to focus on those of the identified molecular mechanisms that made

the more sense to DSS pathophysiology and systemic vascular dysfunction, referring to recent findings on the role of innate immunity in systemic inflammatory processes leading to shock, multi-organ dysfunction syndromes or other pejorative clinical outcomes.

Second, while present results confirm some putative DSS-related biomarkers, it also reveals unreported alterations that make sense to hypovolemic shock pathophysiology. This reinforces the ability of a global and “open mind” approach to identify molecular processes relevant to the studied pathology. Blood cells transcriptional profiles clearly reveal alterations of different immune responses and the activation of a large pro-inflammatory response. A significant proportion of genes of which expression is modified are related to host innate immunity, lymphocyte functions and lipid metabolism in particular. This genome-wide expression analysis also confirms the over-expression of individual biomarkers previously associated with severe dengue, such as the acute phase pentraxin-related protein PTX3, the pro-inflammatory IL-18 cytokine or the anti-inflammatory IL-10 cytokine (Table S2) [11,12,15], providing a more comprehensive overview of their implication in the pathophysiology of DSS.

Our results differ however from those reported by Long *et al* in a genome-wide expression profiling study comparing DSS children with uncomplicated paediatric patients [28]. This study concluded on a global “benign” and “muted” immune transcriptional response but a decreased expression of genes involved in IL-10 and IFN type I-related pathways in DSS children blood cells [28].

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

<footnote>ªpercentage of variance associated to disease phenotype.</footnote>

<footnote>doi:10.1371/journal.pone.0011671.t005</footnote>
Differences in study design, size of cohorts and time of blood sampling from patients in the course of dengue disease may explain these differences. Indeed, in our study, two DSS children had gene expression profiles close to those of uncomplicated DF and DHF, and clustered within the DF/DHF cluster. Both proved to be the children from whom blood was sampled three days after the onset of shock, while the three DSS children sampled two days after shock onset still exhibited a typical DSS gene expression profile. This suggests that a shift from a “severe” to an “uncomplicated” transcriptional profile may occur within a very short time, and could explain the “uncomplicated” and benign gene immune transcriptional responses reported by Long *et al* [28]. Differences in strategies and methods used to filter genes differentially expressed between patients’ groups could also explain the finding that few IFN type I-related genes but a large diversity of other pathways were identified in the present study compared to other transcriptomic studies of DHF or DSS patients. Here, genes were selected considering only their statistical significance and their association with the disease phenotype. Differently from others [27,28,42], no “fold change” cut-off filter was applied since this non-statistically-motivated criteria selects preferentially genes prompt to high variations such as the IFN type I-induced genes [92,93], thus excluding from subsequent bio-informatic analysis a diversity of transcripts exhibiting more subtle variations but strong associations and biological relevance with the considered disease phenotype.

Third, unsuspected mechanisms identified in DSS patients could contribute importantly to the pathophysiology of this severe syndrome, as supported by similarities between those DSS-related alterations and other critical syndromes. Interestingly, a number of immune, repair-remodeling and metabolic-related related pathways are simultaneously altered in the blood cells of DSS children at the onset of shock. In particular, T and NK lymphocyte transcriptional responses are globally impaired while genes implicated in compensatory anti-inflammatory and repair/remodeling responses and in innate immune responses are over-expressed. This highlights the complexity of biological responses at the time of dengue shock syndrome, and points out similarities between DSS and other critical syndromes such as severe sepsis, or post-trauma SIRS that are similarly characterized by depressed T lymphocyte responses but over-expressed innate immunity [94,95].

Reduced abundance of a number of T-lymphocyte related transcripts at the time of DSS may reflect a feed-back mechanism aimed at limiting an initial early T lymphocyte activation, reported to occur in patients who further progress to severe dengue [96,97]. Such a negative feed-back may be sustained by the over-expression of a diversity of anti-inflammatory transcripts in DSS patients’ blood cells at the time of shock. In particular, the two potent immunomodulating factors prostaglandin E synthase and VSIG4, which dampen both T and NK lymphocyte responses [43–45] and have both a strong statistical association with the DSS phenotype, could have such a negative effect. Based on those observations and previous clinical reports, the benefit of corticotherapy in DSS patients might thus be questioned [98,99].

Over-expression in the blood of DSS children of several repair and remodeling genes encoding extracellular matrix proteins, vasoactive mediators and matrix metalloproteases such as the MMP9, likely reflects a compensatory response to inflammatory insults, and a number of those genes products are now considered putative biomarkers in systemic inflammatory syndromes such as severe sepsis [100]. Most proteins encoded by those genes are indeed secreted by activated immune cells such as monocytes/macrophages. They may have adverse effects

towards the vascular endothelium when produced in excess, since they may increase immobilization of inflammatory mediators at the surface of endothelial cells [101], permeability of capillaries [102], or induce direct damage to endothelial tissues [103]. Recently, one of them, MMP9 has been proposed as a putative candidate in the occurrence of plasma leakage during dengue infection [47].

While previous transcriptional studies failed to identify pro-inflammatory gene patterns in the blood cells of DSS patients [27,28,42], our study is the first one to report that a diversity of pro-inflammatory transcriptional responses at the interface of innate immunity, inflammation and host lipid metabolism are activated at the time of cardiovascular failure. Since those mechanisms are considered pathogenic in other systemic inflammatory diseases where systemic vascular dysfunction does occur, we suggest that they may altogether contribute to DSS pathophysiology.

Activation of a pro-inflammatory defence gene pattern in DSS patients’ blood cells (Table 4) has relevance to the pathophysiology of systemic vascular dysfunction, since most microbicidal peptides and enzymes have recognized pro-inflammatory and pathogenic effects towards vascular endothelial tissues [50]. Among them, the neutrophil microbicidal peptides alpha defensins and the highly pro-inflammatory calgranulins proteins S100A8/A9 and S100A12 are now considered putative pathogenic factors in sepsis, cardiovascular diseases, rheumatoid arthritis or atherosclerosis [51,56,104]. While neutrophils are considered the main source of those defence molecules, this cellular origin cannot be established from the present study due to the cellular complexity of unfractionated whole blood samples and to the possibility that other circulating cell types may express a neutrophil-like inflammatory repertoire under pathologic conditions [105]. A putative neutrophil origin of this gene expression pattern is however supported by the over-expression in DSS patients blood cells of transcripts encoding other neutrophil-related molecules such as the MMP8 matrix metalloprotease, CEACAM-6, CEACAM-8 and CD99L2 adhesion molecules (Tables 4 and S2) involved in the recruitment of neutrophils to vascular endothelia. Functional studies should confirm whether those first line-defence immune cells which produce an array of pro-inflammatory mediators highly damaging to host tissues and vascular endothelia [50,106], and poorly regarded in dengue [107], definitively play a role in DSS pathophysiology.

Alteration of a gene pattern related to homeostasis of cholesterol in monocytes/macrophages (Mo/Mac) in the blood cells of DSS children (Table 5) was an unexpected finding, while it should be considered regarding recent knowledge on the role of monocytes as a pivotal link between inflammation, innate immunity and host lipid metabolism [108,109]. Indeed, under physiological conditions, monocytes maintain cholesterol homeostasis by clearing modified LDL such as oxidized LDL (ox-LDL) from plasma. Under pathological conditions, balance between uptake and efflux of those modified cholesterol molecules may be altered [59,60,62,66,67,86], resulting in the intracellular accumulation of modified cholesterol. This turns classical monocytes towards a pro-inflammatory phenotype, lipid-laden monocytes/macrophages (Mo/Mac), a sub-type of pro-inflammatory immune cells initially identified in vascular lesions of chronic inflammatory metabolic diseases [110].

Recent knowledge has shown that those atypical monocytes produce a large array of pro-inflammatory mediators such as ROS, metalloproteases, eicosanoids, and pro-inflammatory adi-pokines, making these cells potent contributors to vascular damages, systemic inflammation and major metabolic changes

<figure><img src="image_15.png" /></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 4. Hypothesis of a second inflammatory amplification loop in dengue shock syndrome.** After induction of a first inflammatory and anti-viral response to dengue virus, disease resolution generally occurs around time of defervescence for most dengue-infected patients. Some patients however progress towards a life-threatening dengue shock syndrome. Results obtained in this study suggest that in those patients, a second inflammatory amplification loop, which involves a diversity of pro-inflammatory responses related to innate immunity, occurs and leads to a major inflammatory systemic syndrome and to vascular homeostasis breakdown. The putative role of different markers identified in vascular endothelial dysfunction is indicated. Thin black arrow, release of; Bold black arrow, interaction between; Punctuated black arrow, chemotactic effect; Thin red arrow, biological activity; Bold red arrow, direct activity on endothelium. DAMPs, danger-associated molecular pattern; GAG, glycosaminoglycane; ROI, reactive oxygen intermediates; TLR, Toll-like receptor.
doi:10.1371/journal.pone.0011671.g004</figcaption></figure>

<figure><img src="image_17.png" /></figure>
such as insulin-resistance [109,111], which altogether characterize systemic inflammatory syndromes such as DSS or severe sepsis. While the existence of functional lipid-laden Mo/Mac during DSS should be established by functional studies, such a molecular mechanism could explain decrease of circulating sub-fractions and total cholesterol previously reported in DSS [112–115] and in other critically ill patients where low cholesterol levels are associated with poor clinical outcome [116,117].

Altered homeostasis of cholesterol in blood cells from DSS patients could also favour replication of dengue viruses into host cells [118], thus contributing to increased viremia in patients with severe dengue infection [119,120], while this could not be evaluated in this study since part of patients had undetectable viremia at the time of blood sampling.

The factors contributing to altered homeostasis of cholesterol in the blood cells of DSS children at time of shock are numerous. Increased lipid peroxydation activity [113,114] insufficiently compensated by anti-oxydant mechanisms [121], as supported by related altered gene patterns identified in this study (Table 5), may result in high levels of circulating ox-LDL contributing to altered cholesterol metabolism. Differences in nutritional status [122–124] or host genetics may also contribute to altered homeostasis of cholesterol gene pattern in the blood cells of DSS patients. Interestingly, transcripts encoding molecules considered candidates to diseases characterized by impairment of cholesterol homeostasis such as *NPC1*, *PCSK9* and *PPARG* [66,67,74], have significant altered abundance in the blood cells of DSS children (Table 5). Further investigations should consider possible associations between DSS and allelic variants of such genes. Whatever the determinants of cholesterol metabolism alterations in DSS patients, our results reinforce interest in considering sub-fractions and total cholesterol as putative biomarkers of DSS [115]. They also suggest that drugs used to treat metabolic disorders such as atherosclerosis should deserve further attention for the control of such a pro-inflammatory process in dengue-infected patients, now proposed for other critical illnesses [117].

Transcriptional activation of the lipid-related arachidonic acid pathway in the whole blood cells of DSS children at the time of shock was another pro-inflammatory mechanism relevant to the pathophysiology of DSS [81]. Activation of this lipid metabolic pathway in innate cells such as neutrophils or lipid-laden monocytes during inflammatory process or infection [125], results in the production of eicosanoid lipid mediators, that are not only physiological regulators of vascular tone and permeability [81] but also potent pro-inflammatory mediators involved in a number of pathologies such as asthma [81]. Interestingly, formation of lipid bodies where eicosanoid synthesis takes place, can be induced by ox-LDL through activation of the PPARγ nuclear-lipid receptor [126], thus suggesting a direct link between the three pro-inflammatory pathways identified in DSS children and a contribution of arachidonic pathway-related inflammatory lipids and oxidative enzymes to the systemic vascular dysfunction leading to DSS.

Fourth, DAMPs and TLRs could be a link from primary to secondary inflammation, leading to DSS. Occurrence of DSS in only some patients at the late phase of infection is likely due to an inadequate control or an amplification of the primary inflammatory response aimed at fighting infection. The pro-inflammatory molecular responses activated in the blood cells of DSS children at time of shock involve a diversity of innate immune mediators that may amplify a first-line inflammatory response mediated by TNF, IL-6 or IL-1, thus contributing to a secondary inflammatory loop. Indeed, a number of repair/remodeling and of defence gene products over-expressed in DSS patients blood cells are considered endogenous danger signals or Danger-Associated Molecular

Patterns (DAMPs) (Tables 3 and 4; molecules with DAMP activity are indicated) capable to trigger secondary systemic inflammatory responses through direct interaction with surface or intracellular receptors such as TLRs or NODs expressed in endothelial or innate immune cells [127]. DAMPs include a diversity of molecules without structural similarity either actively produced by immune cells in the context of an infection or passively secreted by damaged tissues [128,129], now considered key inducers of secondary systemic inflammation in a number of acute inflammatory syndromes [130,131] or chronic diseases [132].

Amplification of inflammation during DSS through direct signalling by molecules harbouring DAMP activity via TLRs, is also supported by the increased abundance of DAMP-induced transcripts as those encoding the pro-inflammatory IL-18 cytokine or the NLRC4/CARD12 intracellular sensor [55]. Interestingly, association of allelic polymorphisms of TLR4 with DSS suggested by De Kruif and colleagues [46] suggests that differential signalling through TLRs may contribute to the severity of dengue disease, as suspected for other pathologies [133]. Accordingly, anti-inflammatory drugs targeting Toll-like receptors are now under development for a number of inflammatory pathologies where innate immunity and TLRs play a central pathogenic role [134].

## DSS pathophysiology: a secondary inflammatory loop hypothesis

To summarize, we report the identification of a specific gene expression profile in the blood cells of DSS children at time of shock, characterizing DSS as a unique entity at the transcriptional level whatever the immunological status of children regarding primary or secondary infection. Major immunological alterations identified at the time of shock are characterized by an altered balance between depressed T lymphocyte responses and exacerbated compensatory and pro-inflammatory innate immune responses that may, finally, be detrimental to the host [135–137], while functional studies should confirm the contribution of those molecular mechanisms to DSS pathophysiology.

Based on recent knowledge on molecular mechanisms altered in other systemic inflammatory diseases, DSS may result from a complex pro-inflammatory network involving a diversity of innate immune effectors sustaining a secondary systemic inflammatory loop, leading in turn to vascular homeostasis breakdown and systemic microcirculatory failure characterizing DSS (Figure 4).

We suggest that drugs available to treat metabolic and other systemic chronic inflammatory diseases could be considered for the treatment of dengue-infected patients before shock occurs, and that a number of bio-markers found altered in DSS patients blood cells should be evaluated, as putative predictive markers of progression to DSS.

## Supporting Information

**Figure S1** Validation of microarray results by RT-PCR. Pearson's correlation was calculated between microarray expression signals (horizontal axis) and Delta Ct values from real-time PCR (vertical axis) for nine genes highly associated to dengue shock syndrome. \*\* Correlation is significant at 0.01.
Found at: doi:10.1371/journal.pone.0011671.s001 (4.94 MB TIF)

**Table S1** Clinical and biological characteristics of each DF, DHF and DSS patient
Found at: doi:10.1371/journal.pone.0011671.s002 (0.04 MB XLS)

**Table S2** List of the 3515 clones corresponding to the 2959 genes differentially expressed between DF, DHF and DSS

<figure><img src="image_18.png" /></figure>
patients, identified using the multi-way ANOVA at a false discovery rate of 10. Clones corresponding to the 2959 genes are listed according to their association to DSS, the first one being the gene of which expression level variance is the most influenced by the clinical phenotype. HUGO gene names are indicated. The variation is the one related to the DSS group relatively to DF and DHF. ANOVA, analysis of variance; DF, dengue fever; DHF, dengue hemorrhagic fever; DSS, dengue shock syndrome; NA, not available. a percentage of variance associated to disease phenotype.

Found at: doi:10.1371/journal.pone.0011671.s003 (0.90 MB XLS)

## Acknowledgments

We greatly thank Pr Y. Buisson for supporting this program, I. Drouet and H. Puggelli for technical support and help in preparation of field work.

## References

1. Gubler DJ (2002) Epidemic dengue/dengue hemorrhagic fever as a public health, social and economic problem in the 21st century. *Trends Microbiol 10:* 100–103.

2. Khun S, Manderson L (2008) Poverty, user fees and ability to pay for health care for children with suspected dengue in rural Cambodia. *Int J Equity Health* 7: 10.

3. Peters KG (1998) Vascular endothelial growth factor and the angiopoietins: working together to build a better blood vessel. *Circ Res* 83: 342–343.

4. Basu A, Chaturvedi UC (2008) Vascular endothelium: the battlefield of dengue viruses. *FEMS Immunol Med Microbiol* 53: 287–299.

5. Pang T, Cardosa MJ, Guzman MG (2007) Of cascades and perfect storms: the immunopathogenesis of dengue haemorrhagic fever-dengue shock syndrome (DHF/DSS). *Immunol Cell Biol* 85: 43–45.

6. Green S, Rothman A (2006) Immunopathological mechanisms in dengue and dengue hemorrhagic fever. *Curr Opin Infect Dis* 19: 429–436.

7. Lin CF, Wan SW, Cheng HJ, Lei HY, Lin YS (2006) Autoimmune pathogenesis in dengue virus infection. *Viral Immunol* 19: 127–132.

8. Murgue B (2009) Severe dengue: questioning the paradigm. *Microbes Infect* 12: 113–118.

9. Fink J, Gu F, Vasudevan SG (2006) Role of T cells, cytokines and antibody in dengue fever and dengue haemorrhagic fever. *Rev Med Virol* 16: 263–275.

10. Libraty DH, Acosta LP, Tallo V, Segubre-Mercado E, Bautista A, et al. (2009) A prospective nested case-control study of Dengue in infants: rethinking and refining the antibody-dependent enhancement dengue hemorrhagic fever model. *PLoS Med* 6: e1000171.

11. Green S, Vaughn DW, Kalayanarooj S, Nimmannitya S, Suntayakorn S, et al. (1999) Elevated plasma interleukin-10 levels in acute dengue correlate with disease severity. *J Med Virol* 59: 329–334.

12. Mustafa AS, Elishbishi EA, Agarwal R, Chaturvedi UC (2001) Elevated levels of interleukin-13 and IL-18 in patients with dengue hemorrhagic fever. *FEMS Immunol Med Microbiol* 30: 229–233.

13. Juffrie M, Meer GM, Hack CE, Haasnoot K, Sutaryo, et al. (2001) Inflammatory mediators in dengue virus infection in children: interleukin-6 and its relation to C-reactive protein and secretory phospholipase A2. *Am J Trop Med Hyg* 65: 70–75.

14. Koraka P, Murgue B, Deparis X, Van Gorp EC, Setiati TE, et al. (2004) Elevation of soluble VCAM-1 plasma levels in children with acute dengue virus infection of varying severity. *J Med Virol* 72: 445–450.

15. Mairuhu AT, Peri G, Setiati TE, Hack CE, Koraka P, et al. (2005) Elevated plasma levels of the long pentraxin, pentraxin 3, in severe dengue virus infections. *J Med Virol* 76: 547–552.

16. Cardier JE, Marino E, Romano E, Taylor P, Liprandi F, et al. (2005) Proinflammatory factors present in sera from patients with acute dengue infection induce activation and apoptosis of human microvascular endothelial cells: possible role of TNF-alpha in endothelial cell damage in dengue. *Cytokine* 30: 359–365.

17. Tseng CS, Lo HW, Teng HC, Lo WC, Ker CG (2005) Elevated levels of plasma VEGF in patients with dengue hemorrhagic fever. *FEMS Immunol Med Microbiol* 43: 99–102.

18. Bozza FA, Cruz OG, Zagne SM, Azeredo EL, Nogueira RM, et al. (2008) Multiplex cytokine profile from dengue patients: MIP-1beta and IFN-gamma as predictive factors for severity. *BMC Infect Dis* 8: 86.

19. Deen JL, Harris E, Wills B, Balmaseda A, Hammond SN, et al. (2006) The WHO dengue classification and case definitions: time for a reassessment. *Lancet* 368: 170–173.

20. Bandyopadhyay S, Lum LC, Kroeger A (2006) Classifying dengue: a review of the difficulties in using the WHO case classification for dengue haemorrhagic fever. *Trop Med Int Health* 11: 1238–1255.

We thank Dr J. Desplans for help in using the Agilent microarrays platform.

We are also indebted to all doctors, nurses, patients and their families who participated to this study at the hospital of Kampong Cham, Cambodia, and to all technicians from the virology department of the Institut Pasteur in Cambodia who carried out dengue diagnosis assays.

## Author Contributions

Conceived and designed the experiments: SD CS VD PB PCP. Performed the experiments: SD CS PCP. Analyzed the data: SD CS AB PR PCP. Contributed reagents/materials/analysis tools: VD SO PTL NC SN PB PCP. Wrote the paper: SD PCP. Contributed to obtainment of funding: HT.

21. Ng CF, Lum LC, Ismail NA, Tan LH, Tan CP (2007) Clinicians' diagnostic practice of dengue infections. *J Clin Virol* 40: 202–206.

22. Cobb JP, Mindrinos MN, Miller-Graziano C, Calvano SE, Baker HV, et al. (2005) Application of genome-wide expression analysis to human health and disease. *Proc Natl Acad Sci U S A* 102: 4801–4806.

23. Feezor RJ, Cheng A, Paddock HN, Baker HV, Moldawer LL (2005) Functional genomics and gene expression profiling in sepsis: beyond class prediction. *Clin Infect Dis* 41 Suppl 7: S427–435.

24. Tang BM, McLean AS, Dawes IW, Huang SJ, Lin RC (2007) The use of gene-expression profiling to identify candidate genes in human sepsis. *Am J Respir Crit Care Med* 176: 676–684.

25. Tian L, Greenberg SA, Kong SW, Altschuler J, Kohane IS, et al. (2005) Discovering statistically significant pathways in expression profiling studies. *Proc Natl Acad Sci U S A* 102: 13544–13549.

26. Laudanski K, Miller-Graziano C, Xiao W, Mindrinos MN, Richards DR, et al. (2006) Cell-specific expression and pathway analyses reveal alterations in trauma-related human T cell and monocyte pathways. *Proc Natl Acad Sci U S A* 103: 15564–15569.

27. Simmons CP, Popper S, Dolocek C, Chau TN, Griffiths M, et al. (2007) Patterns of host genome-wide gene transcript abundance in the peripheral blood of patients with acute dengue hemorrhagic fever. *J Infect Dis* 195: 1097–1107.

28. Long HT, Hibberd ML, Hien TT, Dung NM, Van Ngoc T, et al. (2009) Patterns of gene transcript abundance in the blood of children with severe or uncomplicated dengue highlight differences in disease evolution and host response to dengue virus infection. *J Infect Dis* 199: 537–546.

29. World Health Organization (1997) Dengue haemorrhagic fever: diagnosis, treatment, prevention and control. *Geneva: World Health Organization*. 84 p.

30. Buchy P, Vo VL, Bui KT, Trinh TX, Glaziu P, et al. (2005) Secondary dengue virus type 4 infections in Vietnam. *Southeast Asian J Trop Med Public Health* 36: 178–185.

31. Reyes JM, Ong S, Mey C, Ngan C, Hoyer S, et al. (2003) Improved molecular detection of dengue virus serotype 1 variants. *J Clin Microbiol* 41: 3864–3867.

32. Edgar R, Domrachev M, Lash AE (2002) Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. *Nucleic Acids Res* 30: 207–210.

33. Smyth GK, Michaud J, Scott HS (2005) Use of within-array replicate spots for assessing differential expression in microarray experiments. *Bioinformatics* 21: 2067–2075.

34. Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. *Bioinformatics* 19: 185–193.

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

36. 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.

37. Eisen MB, Spellman PT, Brown PO, Botstein D (1998) Cluster analysis and display of genome-wide expression patterns. *Proc Natl Acad Sci U S A* 95: 14863–14868.

38. Saldanha AJ (2004) Java Treeview-extensible visualization of microarray data. *Bioinformatics* 20: 3246–3248.

39. Brown MP, Grundy WN, Lin D, Cristianini N, Sugnet CW, et al. (2000) Knowledge-based analysis of microarray gene expression data by using support vector machines. *Proc Natl Acad Sci U S A* 97: 262–267.

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

<figure><img src="image_19.png" /></figure>
41. Pavlidis P (2003) Using ANOVA for gene selection from microarray studies of the nervous system. Methods 31: 282-289.

42. Ubol S, Masrinoul P, Chaijaruwanich J, Kalayanooj S, Charoensirisuthikul T, et al. (2008) Differences in global gene expression in peripheral blood mononuclear cells indicate a significant role of the innate responses in progression of dengue fever but not dengue hemorrhagic fever. J Infect Dis 197: 1459-1467.

43. Stobo JD, Kennedy MS, Goldyne ME (1979) Prostaglandin E modulation of the mitogenic response of human T cells. Differential response of T-cell subpopulations. J Clin Invest 64: 1188-1203.

44. Joshi PC, Zhou X, Cuchens M, Jones Q (2001) Prostaglandin E2 suppressed IL-15-mediated human NK cell function through down-regulation of common gamma-chain. J Immunol 166: 885-891.

45. Vogt L, Schmitz N, Kurrer MO, Bauer M, Hinton HI, et al. (2006) VSIG4, a B7 family-related protein, is a negative regulator of T cell activation. J Clin Invest 116: 2817-2826.

46. de Kruif MD, Setiati TE, Mairuhu AT, Koraka P, Aberson HA, et al. (2008) Differential gene expression changes in children with severe dengue virus infections. PLoS Negl Trop Dis 2: e215.

47. Luplertlop N, Misse D, Bray D, Deleuze V, Gonzalez JP, et al. (2006) Dengue-virus-infected dendritic cells trigger vascular leakage through metalloproteinase overproduction. EMBO Rep 7: 1176-1181.

48. Predescu D, Predescu S, Shimizu J, Miyawaki-Shimizu K, Malik AB (2005) Constitutive cNOS-derived nitric oxide is a determinant of endothelial junctional integrity. Am J Physiol Lung Cell Mol Physiol 289: L371-381.

49. Yang D, Biragyn A, Hoover DM, Lubkowski J, Oppenheim JJ (2004) Multiple roles of antimicrobial defensins, cathelicidins, and eosinophil-derived neurotoxin in host defense. Annu Rev Immunol 22: 181-215.

50. DiStasi MR, Ley K (2009) Opening the flood-gates: how neutrophil-endothelial interactions regulate permeability. Trends Immunol 30: 547-556.

51. Foell D, Wittkowski J, Vogl T, Roth J (2007) S100 proteins expressed in phagocytes: a novel group of damage-associated molecular pattern molecules. J Leukoc Biol 81: 28-37.

52. Vogl T, Tenbrock K, Ludwig S, Leukert N, Ehrhardt C, et al. (2007) Mrp8 and Mrp14 are endogenous activators of Toll-like receptor 4, promoting lethal, endotoxin-induced shock. Nat Med 13: 1042-1049.

53. Yang D, Chen Q, Rosenberg HF, Rybak SM, Newton DL, et al. (2004) Human ribonuclease A superfamily members, eosinophil-derived neurotoxin and pancreatic ribonuclease, induce dendritic cell maturation and activation. J Immunol 173: 6134-6142.

54. Sur S, Glitz DG, Kita H, Kujawa SM, Peterson EA, et al. (1998) Localization of eosinophil-derived neurotoxin and eosinophil cationic protein in neutrophilic leukocytes. J Leukoc Biol 63: 715-722.

55. Pedra JH, Cassel SL, Sutterwala FS (2009) Sensing pathogens and danger signals by the inflammasome. Curr Opin Immunol 21: 10-16.

56. Foell D, Frosch M, Sorg C, Roth J (2004) Phagocyte-specific calcium-binding S100 proteins as clinical laboratory markers of inflammation. Clin Chim Acta 344: 37-51.

57. Gordon S (2007) Macrophage heterogeneity and tissue lipids. J Clin Invest 117: 89-93.

58. Barlic J, Murphy PM (2007) Chemokine regulation of atherosclerosis. J Leukoc Biol 82: 226-236.

59. Mosig S, Rennert K, Buttner P, Krause S, Lutjohann D, et al. (2008) Monocytes of patients with familial hypercholesterolemia show alterations in cholesterol metabolism. BMC Med Genomics 1: 60.

60. de Winther MP, van Dijk KW, Havekes LM, Hofker MH (2000) Macrophage scavenger receptor class A: A multifunctional receptor in atherosclerosis. Arterioscler Thromb Vasc Biol 20: 290-297.

61. Febbraio M, Hajjar DP, Silverstein RL (2001) CD36: a class B scavenger, receptor involved in angiogenesis, atherosclerosis, inflammation, and lipid metabolism. J Clin Invest 108: 785-791.

62. Pennings M, Meurs I, Ye D, Out R, Hoekstra M, et al. (2006) Regulation of cholesterol homeostasis in macrophages and consequences for atherosclerotic lesion development. FEBS Lett 580: 5588-5596.

63. Chen XP, Zhang TT, Du GH (2007) Lectin-like oxidized low-density lipoprotein receptor-1, a new promising target for the therapy of atherosclerosis? Cardiovasc Drug Rev 25: 146-161.

64. Renie G, Maingrette F, Li L (2007) Diabetic vasculopathy and the lectin-like oxidized low-density lipoprotein receptor-1 (LOX-1). Curr Diabetes Rev 3: 103-110.

65. Szanto A, Roszer T (2008) Nuclear receptors in macrophages: a link between metabolism and inflammation. FEBS Lett 582: 106-116.

66. Ory DS (2004) The niemann-pick disease genes; regulators of cellular cholesterol homeostasis. Trends Cardiovasc Med 14: 66-72.

67. Zhang JR, Coleman T, Langmade SJ, Scherrer DE, Lane L, et al. (2008) Niemann-Pick C1 protects against atherosclerosis in mice via regulation of macrophage intracellular cholesterol trafficking. J Clin Invest 118: 2281-2290.

68. Sesti G, Federici M, Hribal ML, Lauro D, Sbraccia P, et al. (2009) Defects of the insulin receptor substrate (IRS) system in human metabolic disorders. Faseb J 15: 2099-2111.

69. White MF (2002) IRS proteins and the common path to diabetes. Am J Physiol Endocrinol Metab 283: E413-422.

70. Makowski L, Hotamisligil GS (2005) The role of fatty acid binding proteins in metabolic syndrome and atherosclerosis. Curr Opin Lipidol 16: 543-548.

71. Howard JK, Flier JS (2006) Attenuation of leptin and insulin signaling by SOCS proteins. Trends Endocrinol Metab 17: 365-371.

72. Qatanani M, Szwergold NR, Greaves DR, Ahima RS, Lazar MA (2009) Macrophage-derived human resistin exacerbates adipose tissue inflammation and insulin resistance in mice. J Clin Invest.

73. Boot RG, van Achterberg TA, van Aken BE, Renkema GH, Jacobs MJ, et al. (1999) Strong induction of members of the chitinase family of proteins in atherosclerosis: chitotriosidase and human cartilage gp-39 expressed in lesion macrophages. Arterioscler Thromb Vasc Biol 19: 687-694.

74. Horton JD, Cohen JC, Hobbs HH (2009) PCSK9: a convertase that coordinates LDL catabolism. J Lipid Res 50 Suppl: S172-177.

75. Scatena M, Liaw L, Giachelli CM (2007) Osteopontin: a multifunctional molecule regulating chronic inflammation and vascular disease. Arterioscler Thromb Vasc Biol 27: 2302-2309.

76. Cho HJ, Cho HJ, Kim HS (2009) Osteopontin: a multifunctional protein at the crossroads of inflammation, atherosclerosis, and vascular calcification. Curr Atheroscler Rep 11: 206-213.

77. Forte TM, Subbanagounder G, Berliner JA, Blanche PJ, Clermont AO, et al. (2002) Altered activities of anti-atherogenic enzymes LCAT, paraoxonase, and platelet-activating factor acetylhydrolase in atherosclerosis-susceptible mice. J Lipid Res 43: 477-485.

78. Gilroy DW, Newson J, Sawmynaden P, Willoughby DA, Croxtall JD (2004) A novel role for phospholipase A2 isoforms in the checkpoint control of acute inflammation. Faseb J 18: 489-498.

79. Homaidan FR, Chakroun I, Haidar HA, El-Sabban ME (2002) Protein regulators of eicosanoid synthesis: role in inflammation. Curr Protein Pept Sci 3: 467-484.

80. Khanapure SP, Garvey DS, Janero DR, Letts LG (2007) Eicosanoids in inflammation: biosynthesis, pharmacology, and therapeutic frontiers. Curr Top Med Chem 7: 311-340.

81. Harizi H, Corcuff JB, Gualde N (2008) Arachidonic-acid-derived eicosanoids: roles in biology and immunopathology. Trends Mol Med 14: 461-469.

82. Rossi A, Cuzzocrea S, Sautebin L (2009) Involvement of leukotriene pathway in the pathogenesis of ischemia-reperfusion injury and septic and non-septic shock. Curr Vasc Pharmacol 7: 185-197.

83. Lima JJ, Blake KV, Tantisira KG, Weiss ST (2009) Pharmacogenetics of asthma. Curr Opin Pulm Med 15: 57-62.

84. Wittwer J, Hersberger M (2007) The two faces of the 15-lipoxygenase in atherosclerosis. Prostaglandins Leukot Essent Fatty Acids 77: 67-77.

85. Deng Y, Theken KN, Lee CR (2009) Cytochrome P450 epoxigenases, soluble epoxide hydrolase, and the regulation of cardiovascular inflammation. J Mol Cell Cardiol 48: 331-341.

86. Attie AD, Kastelein JP, Hayden MR (2001) Pivotal role of ABCA1 in reverse cholesterol transport influencing HDL levels and susceptibility to atherosclerosis. J Lipid Res 42: 1717-1726.

87. Hennuyer N, Tailleux A, Torpier G, Mezdour H, Fruchart JC, et al. (2005) PPARalpha, but not PPARgamma, activators decrease macrophage-laden atherosclerotic lesions in a nondiabetic mouse model of mixed dyslipidemia. Arterioscler Thromb Vasc Biol 25: 1897-1902.

88. Samuelsson B, Morgenstern R, Jakobsson PJ (2007) Membrane prostaglandin E synthase-1: a novel therapeutic target. Pharmacol Rev 59: 207-224.

89. Cipollone F, Fazia M, Iezzi A, Ciabattoni G, Pini B, et al. (2004) Balance between PGD synthase and PGE synthase is a major determinant of atherosclerotic plaque instability in humans. Arterioscler Thromb Vasc Biol 24: 1259-1265.

90. World Health Organization (2009) Dengue guidelines for diagnosis, treatment, prevention and control: new edition. Geneva: World Health Organization. 147 p.

91. Dyrskjot L, Zieger K, Real FX, Malats N, Carrato A, et al. (2007) Gene expression signatures predict outcome in non-muscle-invasive bladder carcinoma: a multicenter validation study. Clin Cancer Res 13: 3545-3551.

92. Whitney AR, Diehn M, Popper SJ, Alizadeh AA, Boldrick JC, et al. (2003) Individuality and variation in gene expression patterns in human blood. Proc Natl Acad Sci U S A 100: 1896-1901.

93. Eady JJ, Wortley GM, Wormstone YM, Hughes JC, Astley SB, et al. (2005) Variation in gene expression profiles of peripheral blood mononuclear cells from healthy volunteers. Physiol Genomics 22: 402-411.

94. Remick DG (2007) Pathophysiology of sepsis. Am J Pathol 170: 1435-1444.

95. Lenz A, Franklin GA, Cheadle WG (2007) Systemic inflammation after trauma. Injury 38: 1336-1345.

96. Chaturvedi UC, Shrivastava R, Tripathi RK, Nagar R (2007) Dengue virus-specific suppressor T cells: current perspectives. FEMS Immunol Med Microbiol 50: 285-299.

97. Green S, Vaughn DW, Kalayanooj S, Nimmannitya S, Suntayakorn S, et al. (1999) Early early activation in acute dengue illness is related to development of plasma leakage and disease severity. J Infect Dis 179: 755-762.

98. Panpanich R, Sornchai P, Kanjanaratanakorn K (2006) Corticosteroids for treating dengue shock syndrome. Cochrane Database Syst Rev 3: CD003488.

99. Rajapakse S (2009) Corticosteroids in the treatment of dengue illness. Trans R Soc Trop Med Hyg 103: 122-126.

100. Lorente L, Martin MM, Sole-Violan J, Blanquer J, Paramo JA (2010) Matrix metalloproteinases and their inhibitors as biomarkers of severity in sepsis. Crit Care 14: 402.

<figure><img src="image_20.png" /></figure>
101. Wight TN, Merrilees MJ (2004) Proteoglycans in atherosclerosis and restenosis: key roles for versican. *Circ Res* 94: 1158-1167.

102. Galley HF, Webster NR (2004) Physiology of the endothelium. Br J Anaesth 93: 105-113.

103. Raffetto JD, Khalil RA (2008) Matrix metalloproteinases and their inhibitors in vascular remodeling and vascular disease. *Biochem Pharmacol* 75: 346-359.

104. Ehrchen JM, Sunderkotter C, Foell D, Vogl T, Roth J (2009) The endogenous Toll-like receptor 4 agonist S100A8/S100A9 (calprotectin) as innate amplifier of infection, autoimmunity, and cancer. *J Leukoc Biol* 86: 557-566.

105. Borregaard N, Sorensen OE, Theilgaard-Monch K (2007) Neutrophil granules: a library of innate immunity proteins. *Trends Immunol* 28: 340-345.

106. Quinn K, Henriques M, Parker T, Slutsky AS, Zhang H (2008) Human neutrophil peptides: a novel potential mediator of inflammatory cardiovascular diseases. *Am J Physiol Heart Circ Physiol* 295: H1817-1824.

107. Butthep P, Bunyaratvej A, Bhamarapravati N (1993) Dengue virus and endothelial cell: a related phenomenon to thrombocytopenia and granulocytopenia in dengue hemorrhagic fever. *Southeast Asian J Trop Med Public Health* 24 Suppl 1: 246-249.

108. Michelsen KS, Doherty TM, Shah PK, Arditi M (2004) TLR signaling: an emerging bridge from innate immunity to atherogenesis. *J Immunol* 173: 5901-5907.

109. Wellen KE, Hotamisligil GS (2005) Inflammation, stress, and diabetes. *J Clin Invest* 115: 1111-1119.

110. Schmitz G, Grandl M (2008) Lipid homeostasis in macrophages - implications for atherosclerosis. *Rev Physiol Biochem Pharmacol* 160: 93-125.

111. Yagmur E, Trautwein C, Gressner AM, Tacke F (2006) Resistin serum levels are associated with insulin resistance, disease severity, clinical complications, and prognosis in patients with chronic liver diseases. *Am J Gastroenterol* 101: 1244-1252.

112. van Gorp EC, Suharti C, Mairuhu AT, Dolmans WM, van Der Ven J, et al. (2002) Changes in the plasma lipid profile as a potential predictor of clinical outcome in dengue hemorrhagic fever. *Clin Infect Dis* 34: 1150-1153.

113. Lee CY, Seet RC, Huang SH, Long LH, Halliwell B (2008) Different patterns of oxidized lipid products in plasma and urine of dengue fever, stroke and Parkinsons disease patients. *Cautions in the use of biomarkers of oxidative stress. Antioxid Redox Signal*.

114. Soundravally R, Sankar P, Bobby Z, Hoti SL (2008) Oxidative stress in severe dengue viral infection: association of thrombocytopenia with lipid peroxidation. *Platelets* 19: 447-454.

115. Suvarna JC, Rane PP (2009) Serum lipid profile: a predictor of clinical outcome in dengue infection. *Trop Med Int Health* 14: 576-585.

116. Marik PE (2006) Dyslipidemia in the critically ill. *Crit Care Clin* 22: 151-159, viii.

117. Kruger PS (2009) Forget glucose: what about lipids in critical illness? *Crit Care Resusc* 11: 305-309.

118. Rothwell C, Lebreton A, Young Ng C, Lim JY, Liu W, et al. (2009) Cholesterol biosynthesis modulation regulates dengue viral replication. *Virology* 389: 8-19.

119. Vaughn DW, Green S, Kalayanarooj S, Innis BL, Nimmannitya S, et al. (2000) Dengue viremia titer, antibody response pattern, and virus serotype correlate with disease severity. *J Infect Dis* 181: 2-9.

120. Murgue B, Roche C, Chungue E, Deparis X (2000) Prospective study of the duration and magnitude of viraemia in children hospitalised during the 1996-1997 dengue-2 outbreak in French Polynesia. *J Med Virol* 60: 432-438.

121. Ray G, Kumar V, Kapoor AK, Dutta AK, Batra S (1999) Status of antioxidants and other biochemical abnormalities in children with dengue fever. *J Trop Pediatr* 45: 4-7.

122. Kalayanarooj S, Nimmannitya S (2005) Is dengue severity related to nutritional status? *Southeast Asian J Trop Med Public Health* 36: 378-384.

123. Nguyen TH, Nguyen TL, Lei HY, Lin YS, Le BL, et al. (2005) Association between sex, nutritional status, severity of dengue hemorrhagic fever, and immune status in infants with dengue hemorrhagic fever. *Am J Trop Med Hyg* 72: 370-374.

124. Schaible UE, Kaufmann SH (2007) Malnutrition and infection: complex mechanisms and global impacts. *PLoS Med* 4: e115.

125. Pacheco P, Bozza FA, Gomes RN, Bozza M, Weller PF, et al. (2002) Lipopolysaccharide-induced leukocyte lipid body formation in vivo: innate immunity elicited intracellular Loci involved in eicosanoid metabolism. *J Immunol* 169: 6498-6506.

126. de Assis EF, Silva AR, Caiado LF, Marathe GK, Zimmerman GA, et al. (2003) Synergism between platelet-activating factor-like phospholipids and peroxisome proliferator-activated receptor gamma agonists generated during low density lipoprotein oxidation that induces lipid body formation in leukocytes. *J Immunol* 171: 2090-2098.

127. Castellheim A, Brekke OL, Espevik T, Harboe M, Mollnes TE (2009) Innate immune responses to danger signals in systemic inflammatory response syndrome and sepsis. *Scand J Immunol* 69: 479-491.

128. Oppenheim JJ, Yang D (2005) Alarmins: chemotactic activators of immune responses. *Curr Opin Immunol* 17: 359-365.

129. Bianchi ME (2007) DAMPs, PAMPs and alarmins: all we need to know about danger. *J Leukoc Biol* 81: 1-5.

130. Cinel I, Opal SM (2009) Molecular biology of inflammation and sepsis: a primer. *Crit Care Med* 37: 291-304.

131. Claus RA, Otto GP, Deigner HP, Bauer M (2010) Approaching clinical reality: markers for monitoring systemic inflammation and sepsis. *Curr Mol Med* 10: 227-235.

132. Gill R, Tsung A, Billiar T (2010) Linking oxidative stress to inflammation: Toll-like receptors. *Free Radic Biol Med* 48: 1121-1132.

133. Mockenhaupt FP, Cramer JP, Hamann L, Stegemann MS, Eckert J, et al. (2006) Toll-like receptor (TLR) polymorphisms in African children: Common TLR-4 variants predispose to severe malaria. *Proc Natl Acad Sci U S A* 103: 177-182.

134. O'Neill LA (2003) Therapeutic targeting of Toll-like receptors for inflammatory and infectious diseases. *Curr Opin Pharmacol* 3: 396-403.

135. Kim KD, Zhao J, Auh S, Yang X, Du P, et al. (2007) Adaptive immune cells temper initial innate responses. *Nat Med* 13: 1248-1252.

136. Barton GM (2008) A calculated response: control of inflammation by the innate immune system. *J Clin Invest* 118: 413-420.

137. Zhao J, Kim KD, Yang X, Auh S, Fu YX, et al. (2008) Hyper innate responses in neonates lead to increased morbidity and mortality after infection. *Proc Natl Acad Sci U S A* 105: 7528-7533.

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