key: cord-0001253-0epeljaf
authors: Zapata, Juan Carlos; Carrion, Ricardo; Patterson, Jean L.; Crasta, Oswald; Zhang, Yan; Mani, Sachin; Jett, Marti; Poonia, Bhawna; Djavani, Mahmoud; White, David M.; Lukashevich, Igor S.; Salvato, Maria S.
title: Transcriptome Analysis of Human Peripheral Blood Mononuclear Cells Exposed to Lassa Virus and to the Attenuated Mopeia/Lassa Reassortant 29 (ML29), a Vaccine Candidate
date: 2013-09-12
journal: PLoS Negl Trop Dis
DOI: 10.1371/journal.pntd.0002406
sha: 0efede5fe4185af68b00ac20b99399af133e80bb
doc_id: 1253
cord_uid: 0epeljaf

Lassa virus (LASV) is the causative agent of Lassa Fever and is responsible for several hundred thousand infections and thousands of deaths annually in West Africa. LASV and the non-pathogenic Mopeia virus (MOPV) are both rodent-borne African arenaviruses. A live attenuated reassortant of MOPV and LASV, designated ML29, protects rodents and primates from LASV challenge and appears to be more attenuated than MOPV. To gain better insight into LASV-induced pathology and mechanism of attenuation we performed gene expression profiling in human peripheral blood mononuclear cells (PBMC) exposed to LASV and the vaccine candidate ML29. PBMC from healthy human subjects were exposed to either LASV or ML29. Although most PBMC are non-permissive for virus replication, they remain susceptible to signal transduction by virus particles. Total RNA was extracted and global gene expression was evaluated during the first 24 hours using high-density microarrays. Results were validated using RT-PCR, flow cytometry and ELISA. LASV and ML29 elicited differential expression of interferon-stimulated genes (ISG), as well as genes involved in apoptosis, NF-kB signaling and the coagulation pathways. These genes could eventually serve as biomarkers to predict disease outcomes. The remarkable differential expression of thrombomodulin, a key regulator of inflammation and coagulation, suggests its involvement with vascular abnormalities and mortality in Lassa fever disease.

Lassa Fever (LF) is one of the most neglected tropical diseases and has the second highest global incidence of any viral hemorrhagic fever (VHF), after Dengue HF. The estimated ''at risk'' population in Sierra Leone, Guinea, and Nigeria, may be as high as 59 million. The number of symptomatic cases per year is around 300,000 with between 5,000 and 10,000 deaths [1] [2] [3] . LF is also one of the most prominent imported VHF affecting public health of non-endemic countries worldwide and has been classified as a category A bio-threat agent [4] , [5] , [2] , [6] .

Lassa virus (LASV) belongs to the Arenaviridae, a fast-growing family of rodent-borne viruses, currently including 22 enveloped viruses with bi-segmented, ambisense single-stranded RNA genomes divided into two groups [7] , [8] . The Old World (OW) group contains the prototypic lymphocytic choriomeningitis virus (LCMV) which can cause neuropathology in adults, fetal abnormalities in newborns, and a fatal LF-like disease in immunocompromised patients. OW viruses also include LASV; and non-pathogenic viruses, Mopeia (MOPV), Mobala (MOBV), and Ippy (IPPYV), as well as newly-found LuJo (pathogenic) [9] and Luna (unknown pathogenicity) [10] viruses. The New World group includes Junín, Machupo, Guanarito, Sabia, and Chapare viruses associated with severe VHF in South America [11] .

The disease mechanism for LF is still not clearly understood. As with other VHF, events that determine LF fatal outcome occur very early in the disease and depend on the host immune system [12] , [13] . Macrophages (MP) and dendritic cells (DC) are the initial targets of infection, followed by viremic seeding of a wide variety of cells and organs. High viremia and elevated levels of AST in plasma are two major risk factors predicting fatal outcome [14] . Although the most consistent pathologic finding is hepatocellular necrosis, the extent of liver damage is usually insufficient to implicate hepatic failure as the cause of death [14] , [2] . Functional damage of platelets and endothelial cells seems to play an important role [15] , [16] , [2] . However, contribution of proand anti-inflammatory stimuli to LF pathogenesis is not clearly understood.

Our laboratory [17] , [18] and others [19] , [20] , [21] , [22] showed that LASV infection is associated with suppression of proinflammatory responses contributing to unchecked viremia and fatal outcome. Study of LF pathogenesis is complicated by the absence of reliable and validated animal models. Using a surrogate model of LF in rhesus macaques inoculated with the WE strain of LCMV [23] we showed that early transcriptional changes in blood, mostly in genes involved in complement, interferon, and pro-inflammatory pathways, can discriminate virulent from nonvirulent infection [24] .

There is no approved LF vaccine and therapeutic options are limited to intravenous or oral ribavirin, which is often impractical in endemic areas. Vaccine candidate ML29 is a reassortant containing the L genomic segment, or replication genes, of MOPV and the S segment encoding the major structural proteins of LASV [25] , [26] . Eighteen mutations distinguish the ML29 genome from the parental arenaviruses and likely contribute to ML29 attenuation [27] . In animal models, ML29 is broadly cross-protective against diverse strains of LASV [28] , and was able to elicit LASV-specific immunity in a monkey model for AIDS [29] , [30] .

In the current study we have directly exposed human PBMC from healthy donors to LASV or to its derivative, ML29. We used PBMC profiling previously to show that ML29 elicited less inflammatory gene expression than MOPV [29] , recapitulating similar profiles of ML29 in naïve primates [26] , and in SIVinfected macaques [30] . Human PBMC exposed to virus is a model for the viremic stage of infection in vivo when much of the transcriptome response is due to non-infectious signal transduction by viral particles [31] [24] . PBMC were used here and in other studies [32] [33] [34] rather than cell lines [35] [36] to more closely simulate the situation in vivo in which circulating blood cells are exposed to virus-infected sites. Here we have shown that in vitro PBMC exposure to arenaviruses with different pathogenic potential, LASV and ML29, resulted in differential expression of ISG, apoptotic, NF-kB, and coagulation pathway genes.

Viruses and cells LASV (strain Josiah) was obtained from CDC (Atlanta, GA) and ML29 was previously described [25] , [26] , [27] .

Vero E6 cells (ATCC, CRL-1586) were cultured in Dulbecco's modified Eagle's medium (DMEM, GIBCO-BRL) supplemented with 10% fetal bovine serum (FBS, GIBCO-BRL), 1% penicillinstreptomycin, and L-glutamine (2 mM); incubated at 37uC in 5% CO 2 for 24 hours. Vero cell monolayers were infected with LASV and ML29 at a multiplicity of infection (MOI) of 0.01 and incubated 1 hour at 37uC in 5% CO 2 , washed with PBS and then covered in DMEM 2% FBS. Supernatants were collected at 48 and 72 hours post infection, titrated in Vero cells and stored at 270uC at 10 7 PFU/ml. Human PBMC were obtained by apheresis from healthy volunteers followed by Ficoll-hypaque isolation [37] and resuspended in RPMI 1640 medium with 10% heat-inactivated human AB serum (Sigma).

Monocytes were isolated from purified PBMC by plastic adherence and immature dendritic cells (DC) were generated from them after 5 days culture in supplemented RPMI. Briefly, the method for generating DC was to plate human PBMC at 1.67610 7 cells/2 ml/well in a 6-well plate in RPMI, 10% fetal calf serum. After 2 hrs, non-adherent cells were removed and the adherent monocytes were fed with 2.5 ml RPMI containing 20 ng/ml GM-CSF (Peprotech, Rocky Hill NJ), and 20 ng/ml IL-4 (Peprotech). After 5 days of incubation, cells were harvested, analyzed by flow cytometry, and processed for incubation with virus. On average, 22% of the cells were differentiated, monocytederived DC (DR+, CD11c+). [38] .

Human blood samples were used under an IRB exemption to Dr. Salvato for the use of de-identified human cells obtained through the Red Cross.

Human PBMC and DC exposure 1.2610 7 PBMC from healthy donors were exposed to 1 MOI of LASV or ML29 (BSL-4 facilities at SFBR, San Antonio, TX) for 45 minutes at 37uC in 5% CO 2 and plated in 6 well plates in duplicate. Supernatants and cells were harvested at 4, 8, and 24 hours post-exposure. Cells were re-suspended in Trizol and kept at 270uC for RNA isolation. 1610 7 DC from 3 healthy donors were exposed to 1 MOI of LASV, ML29 or LPS (as an inhibitor control: LPS is usually a positive control, here it was used to inhibit thrombomodulin) for 45 minutes at 37uC in 5% CO 2 (BSL-4 facilities at CDC, Atlanta, GA). Supernatants and cells were harvested at 24 hours post-exposure. Supernatants were kept at 270uC for ELISA analysis.

Gene expression profiling using cDNA microarrays Total RNA was isolated from PMBC samples using the Trizol method (Invitrogen, Carlsbad, CA) followed by a cleaning step with Rneasy mini kit (Qiagen, Valencia, CA), according to the manufacturer's instructions. Quality and quantity of all RNA samples were evaluated on the Agilent 2100 BioAnalyzer 116 system (Agilent Technologies, Palo Alto, CA) by looking at 18 & 28 s rRNA peaks and by the RIN (RNA integrity number). High quality RNA was labeled and hybridized according to Affymetrix protocols using the GeneChip human genome U133 Plus 2.0 array (Affymetrix, Santa Clara, CA) and as described previously [24] . This chip covers the whole human genome using 54,000 probe sets representing approximately 22,000 genes.

The virulent Lassa fever virus (LASV) and the nonpathogenic Mopeia virus (MOPV) infect rodents and, incidentally, people in West Africa. The mechanism of LASV damage in human beings is unclear. There is no licensed Lassa fever vaccine and therapeutic intervention is usually too late. The ML29 vaccine candidate derived from Lassa and Mopeia viruses protects rodents and primates from Lassa fever disease. Peripheral blood mononuclear cells from healthy human subjects were exposed to either LASV or ML29 in order to identify early cellular responses that could be attributed to the difference in virulence between the two viruses. Differential expression of interferon-stimulated genes as well as coagulation-related genes could lead to an explanation for Lassa fever pathogenesis and indicate protective treatments for Lassa fever disease.

Images from each hybridization, were inspected manually and percentage of present calls of each array was checked. Cluster [39] and Principal Component Analyses (PCA) [40] against all conditions using genes with normalized maximum value/minimum value .5 were also performed. Raw data from the arrays were normalized at probe level using a robust multichip average of G+C content algorithm (gcRMA) [41] , [42] and then log2 transformed. The detection calls (Present, Marginal, Absent) for each probeset were obtained using Affymetrix GeneChip Operating Software (GCOS) (http://www.affymetrix.com/browse/products. jsp?productId = 131429&navMode = 34000&navAction = jump& aId = productsNav#1_1). Only genes with at least one present call across all compared hybridizations were selected for further statistical analysis. All data have been submitted under series record GSE41300.

For the identification of differentially-expressed genes the Linear Models for Microarray Data (LIMMA) package was used [43] .

Due to high variation among donors the pairwise comparison Significance Analysis of Microarrays (SAM) was used to eliminate disturbance. PBMC from each donor were compared before and after exposure at different time points.

To identify the over-represented Gene Ontology (GO) terms in each list of differentially expressed genes from 3 donors and 2 repetitions, the GO::TermFinder was used [44] . GO terms with adjusted p value (false discovery rate or FDR) less than 0.05 were selected.

In order to identify over-represented Pathways from each list of differentially expressed genes from 3 donors and 2 repetitions, the Kegg's pathway [with adjusted p value (FDR) less than 0.05 for each list of genes] was used.

Due to the high number of interferon-related genes involved in the LASV and ML29 gene profiles, they were selected for validation by RT-PCR. 24 hour post-exposure (hpe) RNA samples, from 3 donors, and 3 time points used in the microarray test, were evaluated by the Interferon and Receptor PCR Array Kit (Human Interferon and Receptor array, Cat# No. PAHS-064, SABiosciences, Frederick, MD). This platform tested 84 interferon-related genes and 5 housekeeping genes for data normalization. Briefly, 1 mg of RNA from PBMC was reversetranscribed in a 20-mL reaction volume into first-strand cDNA using SuperArray's ReactionReady First Strand cDNA Synthesis Kit (Cat. No C-01, SABiosciences) containing random and oligo dT primers. After mixing the cDNA with RT2 SYBR Green/ ROX qPCR Master Mix (Cat. No PA-012, SABiosciences), 10 mL of the cDNA mixture was dispensed into a 96 well format plate RT 2 profiler PCR array (Cat. No PAHS-064, SAbioscience).

Real-time PCR was performed on an ABI Fast 7900 Real-Time PCR System (Applied Biosystems, Foster City, CA) using the following cycling parameters: 10 min at 95uC (heat activation step); 40 cycles of 15 sec at 95uC, 1 min at 60uC. Relative changes in gene expression were calculated using the DDC t (threshold cycle) method. This method first subtracts the ct (threshold cycle number) of the gene-average ct of the 5 house keeping genes on the array to normalize to the RNA amounts. Finally the DDct is calculated which involves subtracting the normalized average ct of LASV samples -normalized average ct of the Control Samples. Then this DDct is raised to the negative power of 2 in order to calculate the fold changes. The p value was calculated using a 2tailed Student's t test. These results were treated to SAM data analysis to avoid donor variation.

Whole blood from healthy donors was collected in heparin tubes, mixed with 20 ng/ml lipopolysaccharide (LPS Cat No-tlrlekpls, InvivoGen, San Diego, CA) as down-regulation control for thrombomodulin [45] , MEM as non-stimulated control and ML29 virus. Then this mix was brought to 3 ml with RPMI 10% FBS and incubated for 4, 8, or 24 hours. For each time-point cells were stained with phycoerytrhin-conjugated mouse anti-human thrombomodulin monoclonal antibody (PE-CD141, BD Biosciences, San Jose, CA) and fluorescein isothiocyanate mouse anti-human CD14+ antibody (FITC-CD14+ BD Biosciences, San Jose, CA). Cells were collected (10,000 events per sample) using a FACSCalibur cytometer (BD Biosciences) and the data were analyzed using FlowJo 8.8.4 software (Tree Star, San Carlos, CA). The monocyte population was gated and percentages of double positive cells (CD14+ and CD141+) were determined.

Dendritic cells were isolated from 3 different individuals for virus exposure and expression of soluble thrombomodulin. Immature DC were derived from human PBMC as described above in ''Viruses and cells''.

Soluble thrombomodulin was quantified by ELISA (Human Thrombomodulin ELISA kit. Cat. No ab46509, Abcam, Cambridge, MA), in DC supernatants after 24 hour of exposure to LPS, ML29, LASV or LASV plus LPS. Plates were read at 450 nm using SpectraMax M2 ELISA reader (Molecular Devices Corporation, Sunnyvale, CA). The estimated amount of thrombomodulin was calculated using the Softmax pro 4.8 program (Molecular Devices Corporation, Sunnyvale, CA) and presented in picograms per ml.

In a previous transcriptome profiling study using a surrogate LF model in monkeys we showed that significant changes in PBMC gene expression occurred during pre-viremic days 1, 2, and 3 [24] , clearly differentiating pathogenic versus non-pathogenic infections even before detectable virus replication. Clearly, a large part of the response to viruses is due to signal transduction of uninfected bystander cells.

In this study, PBMC from three healthy donors were exposed to LASV or ML29 viruses at 1 MOI in duplicate and were collected at 4, 8, and 24 hours post-exposure (hpe) for isolation of total RNA. cDNA was hybridized to Affymetrix chips covering the whole human genome. Cluster and Principal Component Analysis (PCA) showed that all samples clustered together by time (4, 8, and 24 hpe), by donor (3 healthy donors) and by treatment (unexposed, LASV-exposed, and ML29-exposed cells) ( Figure  S1 ). Since donor variation exceeded treatment variation, pairwise Significance Analysis of Microarrays (SAM) was used to compare RNA samples before and after exposure for each donor at different time points. Differentially-expressed genes were defined as those that had at least a .2-fold changes with respect to media controls.

A total of 122 genes were identified as being differentially expressed at tested time-points ( Figure 1 and Table S1 ). Gene expression changes in ML29 and LASV-exposed PBMC were detected at 4 and 8 hpe respectively, with all genes up-regulated in LASV in contrast to ML29 with 11 down-regulated genes (4 at 4 hpe and 7 at 8 hpe). At 24 hpe LASV-treated PBMC showed the maximum changes (122 genes) in comparison with ML29exposed samples.

Gene ontology analysis of global changes in LASV versus ML29-exposed PBMC In order to have a general view of the pathways affected by both virus exposures we performed gene ontology (GO) analysis of differentially expressed genes. The analyses revealed that the most affected genes were involved in type I IFN-stimulated pathways (21%-17 genes from 122), apoptosis (12.29%-15/122), and NFkappa B pathways (10.65%-13/122) ( Figure 2A) . Although, Several genes could be put into more than one category, the interferon-related genes had the higher percentage of presence. From those categories the immune response, defense response and response to viruses were selected for validation analysis ( Table 1 ). The most differentially expressed genes between LASV and ML29 are shown in Table 2 .

LASV exposure up-regulates genes involved in innate immunity: IFN-stimulated genes, transcription factors, and genes involved in apoptosis As shown in Figure 2A , the most striking differences involved gene expression in three different functional categories: IFNstimulated, apoptosis, and NF-kB pathways, during the late timepoint, 24 h after exposure.

To confirm the observed gene expression changes in PBMC exposed to LASV and ML29 viruses, the most affected probesets, the ISG, were selected for further validation. RNA samples from the 24-hr LASV-and ML29-exposed cells were assayed in duplicate using a human IFN and Receptors qRT-PCR Superarray kit (see Methods). Levels of expression of 17 gene-targets were expressed in fold changes (Table 1 and figure 2B ).

Exposure of human PBMC to LASV resulted in over-expression of all 17 IFN-stimulated genes by 2-to 100-fold over background. Exposure of PBMC to ML29 had relatively little effect on the 17 tested genes. Comparative analyses showed a good correlation (76%) between the two techniques, DNA array and IFN qRT-PCR array (Table 1) . For a few genes measured in the DNA array, ML29 exposure gave a higher background than quantitation by real-time PCR. This is because the housekeeping genes we used for normalizing the PCR data, including GAPDH, were not constantly expressed in the compared data sets, whereas probesets used for normalizing the microarray data varied less than 0.04% in all data sets (see Methods).

The most affected gene, IFI6, is the first IFN-stimulated gene known to inhibit apoptosis [46] . Among other most affected genes were IFI35, IFITM1, IRF-7, and SP110. While the function of the IFI35 gene is not well-defined, the IFITM proteins have roles in immune cell signaling and adhesion, cancer, germ cell physiology, and bone mineralization. IFITMs are also viral restriction factors that block entry of Influenza, Ebola, Dengue and West Nile viruses but not arenavirus-envelope pseudotyped particles [47] [48] [49] . IRF-7 is over-expressed in LASVexposed PBMC and others found that it could be up-regulated in 3-day LASV-infected monocyte-derived macrophages [50] . IRF-7 is also known to limit the burst size of LCMV in culture but is Table S1 in the supplemental material. Green color represents down-regulated genes (two-fold or more when compared with non-exposed control). Red color represents up-regulated genes (two-fold or more when compared with non-exposed control). doi:10.1371/journal.pntd.0002406.g001 not ultimately essential for LCMV clearance in vivo [51] . IRF-7 and IRF-3 are transcription factors which together with NF-kB and ATF-2/cJUN induce production of IFN-I. IRF-3 is required for the activation of IFN-b, which in turn, primes the expression of most IFN-a genes by IFN-induced IRF-7 through the STAT-1 pathway. Human SP110 plays an important role in resistance to intracellular pathogens such as human mycobacteria and intracellular cancers. Recently cellular proteins have been described that mediate the interaction of SP110 and viral proteins during viral infections [52] . Among apoptosis-related genes, the most striking difference was observed in the expression of matrix metalloproteinase, MMP-1 (interstitial collagenase), up-regulated .65-fold in LASV-exposed cells versus ML-29-exposed cells ( Table 2 ). Proteins of the MMP family normally breakdown extracellular matrix during embryonic development, reproduction, and tissue remodeling, as well as in disease processes, such as arthritis and metastasis. MMP-1 is known to interact with Tat, the HIV transactivating protein, and this interaction results in the degradation of Tat and depression of Tat-mediated neurotoxicity and transactivation [53] .

Two other apoptosis-related genes, SPP1 and INHBA, were upregulated 14-and 20-fold, respectively (Table 2) . SPP1, a secreted arginine-glycine-aspartate (RGD)-containing phosphoprotein 1 (previous names: osteopontin, bone sialoprotein I), enhances IFNc and IL-12 production and contributes to host defense, bone formation, and wound healing by stimulating macrophage migration as well as protecting against viral and bacterial infections through its pro-Th1 effect [54] . Inhibin-beta A (INHBA) joins Inhibin-alpha to form a pituitary follicle-stimulating hormone (FSH) secretion inhibitor. Inhibin has been shown to negatively-regulate gonadal stromal cell proliferation and to have tumor-suppressor activity. INHBA is also up-regulated by vaccinia and inhibits a diverse array of cytokines by enabling vaccinia virus E3 protein to antagonize several distinct but interlinked signaling cascades [55] .

Exposure of human PBMC to LASV was associated with strong up-regulation (22-25-fold) of CXCL6 and IL1A genes in comparison to ML29 exposure ( Table 2) . Among CXC chemokines, CXCL6, a granulocyte chemotactic protein-2, shares the highest functional homology with IL-8. Both CXCL6 and IL1A exhibit disparate regulatory effects on recruited neutrophil responses [56] [57] .

Due to the nature of LF disease, six LASV-upregulated genes related to the coagulation pathway were analyzed (Table 3) . Thrombomodulin (THBD), heparin-binding growth factor (HBGF), and integrin alpha M (ITGAM), were up-regulated after LASV exposure at all time-points in contrast to ML29 that downregulated them. Plasminogen activator urokinase (PLAU) and integrin beta3 (ITGB3) were mostly up-regulated with LASV exposure and remained unchanged or down-regulated with ML29. PD-ECGF (TYMP) was slightly up-regulated after exposure to both viruses. The analysis of those coagulation-related genes using the Gene Ontology (GO) software showed interaction between them and genes related to immune response ( Figure 3 ). In order to validate cDNA array results, THBD expression was measured on the surface of ML29 exposed-monocytes by flow cytometry, or quantified by ELISA in supernatants of DC exposed to LASV or ML29. There were no significant changes in thrombomodulin expression in CD14+ cells from PBMC treated with LPS or ML29 when compared with non-treated cells at 4 hpe (data not shown). However, as shown in figure 4A , there was a significant reduction in THBD expression on the surface of CD14+ cells after LPS (from ,74% control to 52% at 8 hpe and 47% at 24 hpe) or ML29 stimulus (from ,74% control to 57% at 8 hpe and 51% at 24 hpe), all in agreement with microarray results showing lower THBD RNA in ML29-exposed PBMC ( Figure 1 and Table S1 ). DC cultures exposed to LASV secreted significantly higher levels of THBD than DC exposed to ML29 (Figure 4B ), in agreement with transcriptome data. Control DC secreted 0.502 pg/ml THBD, whereas Lassa-exposed DC expressed 0.724 pg/ml THBD, significantly more than ML29-exposed DC (0.161 pg/ml). Treatment with LASV and LPS together diminished the THBD expression from 0.724 pg/ml after LASV exposure to 0.481 pg/ml. The other coagulation-related genes could also contribute to LF disease: Plasminogen activator urokinase (PLAU) expression could be contributing to the anti-clotting effect mediated by LASV exposure. Heparin-binding epidermal growth factor-like (HB-EGF) expression is involved in several blood vessel physiologies [58] , in kidney pathology [59] , [60] , inhibition of NF-kappa B activation [61] , in cardiac hypertrophy, and together with PD-ECGF promotes vascular maturation [62] , [63] . The ITGB3 (CD61, platelet marker) was 12-fold up-regulated in cells exposed to LASV at 24 h. Integrins are integral cell-surface proteins composed of a and b chains. A given chain may combine with multiple partners resulting in different integrins. The ITGB3 is found along with the aIIb chain in platelets working as a receptor for fibrinogen, von Willebrand factor, plasminogen, prothrombin, thrombospondin, fibronetin, osteospontin (SPP1), and vitronectin [64] , [65] , [66] . Integrins are known to participate in cell adhesion as well as cell-surface mediated signaling. High affinity LASV binding to cellular alpha-dystroglycan (a-DG) receptor perturbs the signaling cross- Cell surface THBD was detected by flow cytometry on gated CD14+ cells. This experiment was repeated with PBMC from 4 different healthy donors with similar results. After exposure to either LPS or ML29, cellsurface THBD was down-regulated by approximately 30%, at both 8 and 24 hours post exposure (hpe). According to a paired Student's t test, this down-regulation was significant (p,0.05) at both time-points. B). Soluble THBD measured by ELISA in dendritic cell supernatants exposed to LPS inhibitor, ML29, LASV and LASV plus LPS 24 hpe. Standard deviation of THBD concentrations are shown from three different experiments. THMD from ML29 cells is significantly depressed with respect to THBD from non-exposed control cells. Lassa-exposed cells express significantly more THBD than ML29-exposed cells. doi:10.1371/journal.pntd.0002406.g004 talk between DG and a6b1 integrins, shifting the normal signaling equilibrium towards inhibition of the MEK/ERK pathway [67] . The DG is a ubiquitous receptor for extracellular matrix (ECM) proteins, which cooperates with b1 integrins to control cell-matrix interactions.

When this paper was in post-review, Malhotra et al. published results of transcriptome profiling PBMC from cynomolgous macaques infected with LASV [68] . Their results are in good correlation with the current study and with our previous results in the LCMV-infected monkey model for LF [24] (Table S2) .

As would be expected from PBMC exposure to other RNA viruses, transcriptome analysis showed that the set of IFNstimulated genes (ISG) is the most affected after exposure to LASV or ML29. The antiviral-ISG, including GBP1, IFI44, IFIH1, IFIT1, MX1, and LY6E, were the most common group of up-regulated genes after influenza infection [69] , yellow fever vaccination [70] , and during febrile episodes of dengue hemorrhagic fever [71] , [72] .

In comparison to the strong up-regulation by LASV, ML29 has muted expression of INHBA (a negative regulator of IFN-c), IFI44, TNFSF10 (TRAIL), SPP1, and LY6E (RIG-E) genes in the interferon pathway, as well as other genes related to immune response such as integrin alpha-M/beta-2 (ITGAM). These genes were down-regulated at 4, 8 and, some, 24 hpe, suggesting a mechanism to avoid recognition by the innate immune system, and to interfere with essential signals in the pathway that leads to the development of type I immunity [73] , [74] , [75] . For example, interferon-alpha/beta inducible IFI44 is involved in the antiviral action of type I interferon [76] , [77] . ML29 down-regulation of IFI44 probably allows the virus to replicate early during the infection so that a strong immune response can control viral replication and spreading. In support of this notion, other antiviral genes such as RSAD2 (Viperin), and TRIM5 [78] , [79] , [80] were more expressed in LASV than in ML29-exposed PBMC.

LASV up-regulates CXCL6, an IL-8 homolog, as well as three other genes that control cell trafficking: monocyte-attracting chemokines CCL7 [81] and CCL8 (or monocyte chemotactic protein-2, MCP-2 [82] [83] , and the mannose receptor C type 1 (MRC1) [84] . These in vitro responses may be related to clinical observations of higher serum IL-8 levels in patients with acute non-fatal LF than in control subjects. Levels of IL-8 and other proinflammatory cyto/chemokines (IL-1b, IL-6, IL-10, IP-10, TNFa) were low or undetectable in patients with fatal LF [21] supporting our previous observations in primary human cells [18] . These observations were also confirmed in animal models of fatal LF [20] , [17] , [85] , [30] providing no evidence for the type of ''cytokine storm'' observed in filovirus HF [86] , [87] . The idea that an uncontrolled replication of LASV in targets cells, DC and macrophages, is associated with a marked suppression of innate and adaptive immune responses resulting in clinical progression and death of experimental animals and LF patients is a wellaccepted concept.

Bleeding is not a salient manifestation in most cases of LF. In contrast to South American arenaviral HF caused by Junín and Machupo viruses, disseminated intravascular coagulation (DIC) is rarely seen in LF patients. Nevertheless, virus-induced coagulation abnormalities and impaired vascular function is a hallmark of LF pathogenesis [3] , [86] . Platelets remain either normal or moderately low, but with significantly reduced function. Plasma from LF patients has inhibitory activity on platelets from healthy individuals, suggesting the presence of soluble factors causing platelets malfunction [88] . Our previous work supports the notion that LASV can replicate vigorously in endothelial cells without cytopathic effects but with reduced cell function, such as reduced production of interleukin IL-8 [18] .

Lymphopenia and impaired lymphocyte proliferative responses are also characteristic of severe Lassa fever [89] [88] . Consequently, in addition to hematopoietic effects of soluble factors, lymphocyte apoptosis due to LASV up-regulation of MMP-1 and INHBA may contribute to Lassa pathogenesis.

Coagulation defects associated with LF disease could potentially be affected by LASV-upregulation of six coagulation-related genes, the most notable being thrombomodulin (THBD). THBD is predominantly synthesized by vascular endothelial cells as a 60 kDa type I trans-membrane protein that inhibits thrombotic, inflammatory and redox related responses [90] . THBD forms a 1:1 molar complex with thrombin and significantly enhances the rate of thrombin inactivation by ATIII while accelerating activation of protein C (aPC). We observed levels of secreted THBD in dendritic cell (DC) culture medium that correlated with our transcriptome data. DC exposure to LASV up-regulated secretion of THBD, whereas exposure to ML29 had the opposite effect. Notably, LPS treatment of LASV-exposed cells resulted in lowering THBD as well. In the presence of thrombin, THBD generates aPC, which inhibits pro-coagulant and pro-inflammatory responses including: fibrinogen cleavage, factor V activation, platelet activation or local cytokine-induced chemotaxis for monocytes, neutrophils and up-regulation of leukocyte adhesion molecules [91] , [92] , [93] . It has been shown that THBD mRNA half-life is shortened by treatment of cells with IFN-c [94] . This observation in addition to the LASV-increased THBD expression seen here, and the absence of IFN responses seen in vitro, are in agreement with LF findings where there is no clear inflammatory response and no cellular infiltrates, in conjunction with the high viral loads seen in the presence of hemorrhagic symptoms. The absence of THBD up-regulation in the LCMV-WE infected model for LF may be a key failing of the model (see details in Table S2 ).

Our current hypothesis is that LASV infects endothelial and antigen-presenting cells and increases the expression of both membrane and soluble THBD forms, activating aPC or capturing HMGB1 with anti-inflammatory and anti-coagulation effects ( Figure 5 ). The up-regulation of THBD and its consequent anticoagulation effects could also explain why fibrin deposits are rarely seen in histological examination of tissues from LF patients and disseminated intravascular coagulation (DIC) is rare in LF, but frequently observed in other viral hemorrhagic fevers. Given these results, inhibition of THBD could be a potential therapeutic intervention for LF.

LASV-induced perturbation of the cross-talk between integrins and the extracellular matrix [67] may contribute to the functional alterations of epithelial and vascular endothelial cells that precede shock and death in LF patients [95] . It is also possible that, during high states of viremia, LASV binds directly to platelets inhibiting their function. After an arterial injury, leukocytes are recruited to the site and there is an increase in ITGAM expression allowing adhesion of platelets and fibrinogen to the vessel walls [96] . We propose that LASV induction of ITGAM could increase platelet adhesion to tissues and this may explain the decreased number of circulating platelets during hemorrhagic episodes. Platelet-derived endothelial cell growth factor (PD-ECGF) is involved in angiogenesis and wound repair [97] , [98] . The differential expression of all those genes may be contributing to the non-pathogenic phenotype of ML29 and to the pathogenic phenotype of LASV.

In LASV-infected primates, DC and macrophages were identified as prominent targets at early stages of infection, while Kupffer cells, hepatocytes, adrenal cortical cells, and endothelial cells were more frequently infected at late stages [99] . We showed previously that mononuclear cells in PBMC are not permissive for LASV and the virus can replicate only in differentiated monocytes, DC and macrophages [18] , representing less than 1% of circulating PBMC. Since LASV and ML29 replicate at the same rate in monocyte-derived cultures (ISL unpublished), the differential gene expression reported here is most likely due to differential signaling of uninfected bystander cells.

Both viruses, LASV and ML29, have identical surface-exposed GP1 responsible for interaction with a-DG; however ML29 and LASV GP2 differ by the lysine to glutamic acid (K272E) substitution located between the fusion domain and the RRLL motif, the cellular subtilase SKI1/S1P cleavage site [26] [7] . The amount of negative charge at that site in ML29 is greater than in other arenaviral GP2, and could potentially affect fusion and signaling by ML29.

In summary, exposure of human PBMC to viruses modeled the viremic stage of disease, and revealed clear differences between responses to LASV and the attenuated ML29. It remains to be shown whether the charge difference in LASV and ML29 GP2 proteins accounts for their different transcriptome profiles. ISG and genes involved in apoptotic and NF-kB pathways were among the most LASV-affected genes at 24 h after exposure. LASV and not ML29 also affected genes involved in coagulation pathways, e.g. LASV strongly up-regulates THBD suggesting a connection to the vascular abnormalities in fatal cases of LF. Based on our current knowledge about THBD involvement in pro-coagulant and pro-inflammatory pathways, THBD inhibition could be a promising drug target. Additionally, it is important to further investigate the immune response-related genes to elucidate the protective mechanisms that are undermined during LF pathogenesis. Figure S1 Cluster analysis showed that all samples clustered together by time (4, 8, and 24 hpe), by donor (3 healthy donors) and by treatment (unexposed, LASV-exposed, and ML29-exposed cells). Sample IDs are organized by Time-Donor-Stimulus-Sample number. Red color represents up-regulated genes and green color represents down-regulated genes. (TIF)

Table S1 A total of 122 genes were identified as being differentially expressed at tested time-points. 136 appear in the table due to the use of more than one probe to detect the same gene. Column A shows the Affymetrix gene name. Column B corresponds to the probe identification reference. Columns C to E are the comparison between LASV exposure and un-exposed cells at different time points. Columns F to H are the comparison between ML29 and un-exposed cells at different time points. Column I has the Unigene ID numbers, and column J has the GenBank accession numbers. The numbers indicate fold-changes in gene expression when compared with the control. Red color represents the more significant upregulated genes. Green color represents the more significant down-regulated genes. (TIF )   Table S2 SAM ANALYSIS of the RNA profiles of PBMC exposed to LASV and ML29 in vitro compared to RNA profiles of PBMC from WE/ARM infected monkeys. ,140 genes with common patterns and ,7 with different patterns, all at 24 hpe. The navy blue highlighting indicates the discrepancies in transcriptome between the LASV/ML29 pair and the LCMV-WE/ARM pair. Yellow highlighting indicates interferon pathway genes and coagulation-related genes. (TIF) Author Contributions

Lassa fever: epidemiology, therapy and vaccine development

Lassa fever

Lassa fever, a new virus disease of man from West Africa. I. Clinical description and pathological findings

A prospective study of the epidemiology and ecology of Lassa fever

Historical Lassa fever reports and 30-year clinical update

Lassa Virus Genome

Virus Taxonomy, Ninth Report of the International Committee on Taxonomy of Viruses

Genetic detection and characterization of Lujo virus, a new hemorrhagic feverassociated arenavirus from southern Africa

Molecular surveillance and phylogenetic analysis of Old World arenaviruses in Zambia

Phylogeny of the genus Arenavirus

Exotic emerging viral diseases: progress and challenges

Visualizing antigenspecific and infected cells in situ predicts outcomes in early viral infection

Lassa fever. Effective therapy with ribavirin

Physiological and immunologic disturbances associated with shock in a primate model of Lassa fever

Hematologic dysfunction in Lassa fever

Arenavirus-mediated liver pathology: acute lymphocytic choriomeningitis virus infection of rhesus macaques is characterized by high-level interleukin-6 expression and hepatocyte proliferation

Lassa and Mopeia virus replication in human monocytes/macrophages and in endothelial cells: different effects on IL-8 and TNF-alpha gene expression

Lassa virus infection of human dendritic cells and macrophages is productive but fails to activate cells

Early and strong immune responses are associated with control of viral replication and recovery in lassa virus-infected cynomolgus monkeys

Low levels of interleukin-8 and interferon-inducible protein-10 in serum are associated with fatal infections in acute Lassa fever

Cutting edge: impairment of dendritic cells and adaptive immunity by Ebola and Lassa viruses

Hemorrhagic fever occurs after intravenous, but not after intragastric, inoculation of rhesus macaques with lymphocytic choriomeningitis virus

Early blood profiles of virus infection in a monkey model for Lassa fever

Lukashevich IS (1992) Generation of reassortants between African arenaviruses

A live attenuated vaccine for Lassa fever made by reassortment of Lassa and Mopeia viruses

Molecular characterization of a reassortant virus derived from Lassa and Mopeia viruses

A ML29 reassortant virus protects guinea pigs against a distantly related Nigerian strain of Lassa virus and can provide sterilizing immunity

Safety, immunogenicity, and efficacy of the ML29 reassortant vaccine for Lassa fever in small non-human primates

An attenuated Lassa vaccine in SIV-infected rhesus macaques does not persist or cause arenavirus disease but does elicit Lassa virus-specific immunity

Use of differential display analysis to assess the effect of human cytomegalovirus infection on the accumulation of cellular RNAs: induction of interferon-responsive RNAs

Modeling of SEB-induced host gene expression to correlate in vitro to in vivo responses

Blood genomic profiles of exposures to Venezuelan equine encephalitis in Cynomolgus macaques (Macaca fascicularis)

Early indicators of exposure to biological threat agents using host gene profiles in peripheral blood mononuclear cells

Differential signaling networks induced by mild and lethal hemorrhagic fever virus infections

Analysis of gene expression in Lassa virusinfected HuH-7 cells

Separation of peripheral leukocytes by Ficoll density gradient centrifugation

Binding and uptake of agalactosyl IgG by mannose receptor on macrophages and dendritic cells

Cluster analysis and display of genome-wide expression patterns

Principal component analysis for clustering gene expression data

Summaries of Affymetrix GeneChip probe level data

A model based background adjustment for oligonucleotide expression arrays

Linear models and empirical bayes methods for assessing differential expression in microarray experiments

GO::TermFinderopen source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes

Lipopolysaccharide down-regulates the thrombomodulin expression of peripheral blood monocytes: effect of serum on thrombomodulin expression in the THP-1 monocytic cell line

G1P3, an interferon inducible gene 6-16, is expressed in gastric cancers and inhibits mitochondrial-mediated apoptosis in gastric cancer cell line TMK-1 cell

The IFITM proteins mediate cellular resistance to influenza A H1N1 virus, West Nile virus, and dengue virus

IFITM3 inhibits influenza A virus infection by preventing cytosolic entry

Distinct patterns of IFITM-mediated restriction of filoviruses, SARS coronavirus, and influenza A virus

Role of interferons in the control of Lassa virus replication in human dendritic cells and macrophages. Microbes and infection

Role of interferon regulatory factor 7 in T cell responses during acute lymphocytic choriomeningitis virus infection

Identification of proteins interacting with human SP110 during the process of viral infections

Interaction of HIV Tat and matrix metalloproteinase in HIV neuropathogenesis: a new host defense mechanism

Eta-1 (osteopontin): an early component of type-1 (cell-mediated) immunity

Vaccinia virus E3 suppresses expression of diverse cytokines through inhibition of the PKR, NF-kappaB, and IRF3 pathways

In vivo neutrophil recruitment by granulocyte chemotactic protein-2 is assisted by gelatinase B/MMP-9 in the mouse

IL-1alpha modulates neutrophil recruitment in chronic inflammation induced by hydrocarbon oil

Heparin binding EGF is necessary for vasospastic response to endothelin

Heparin-binding EGF-like growth factor contributes to reduced glomerular filtration rate during glomerulonephritis in rats

Heparin-binding EGF-like growth factor mRNA is upregulated in the peri-infarct region of the remnant kidney model: in vitro evidence suggests a regulatory role in myofibroblast transformation

Inhibition of NF-kappa B activation and its target genes by heparin-binding epidermal growth factor-like growth factor

Dual intracellular signaling by proteolytic cleavage of membrane-anchored heparin-binding EGF-like growth factor

Endothelial-derived PDGF-BB and HB-EGF coordinately regulate pericyte recruitment during vasculogenic tube assembly and stabilization

Glycoprotein Ic-IIa functions as an activation-independent fibronectin receptor on human platelets

The membrane glycoprotein Ia-IIa (VLA-2) complex mediates the Mg++-dependent adhesion of platelets to collagen

Structure and function of the platelet integrin alphaIIbbeta3

Binding of Lassa virus perturbs extracellular matrix-induced signal transduction via dystroglycan

Transcriptional profiling of the circulating immune response to lassa virus in an aerosol model of exposure

Global host immune response: pathogenesis and transcriptional profiling of type A influenza viruses expressing the hemagglutinin and neuraminidase genes from the 1918 pandemic virus

Systems biology approach predicts immunogenicity of the yellow fever vaccine in humans

Gene expression profiling during early acute febrile stage of dengue infection can predict the disease outcome

Gene expression patterns of dengue virus-infected children from nicaragua reveal a distinct signature of increased metabolism

Integrin alphavbeta3 is a coreceptor for human cytomegalovirus

Osteopontin as a means to cope with environmental insults: regulation of inflammation, tissue remodeling, and cell survival

Role of osteopontin in adhesion, migration, cell survival and bone remodeling

Induction of the human gene for p44, a hepatitis-C-associated microtubular aggregate protein, by interferon-alpha/beta

HIV-1 Tat reprograms immature dendritic cells to express chemoattractants for activated T cells and macrophages

TRIM8/ GERP RING finger protein interacts with SOCS-1

The cytoplasmic body component TRIM5alpha restricts HIV-1 infection in Old World monkeys

Promyelocytic leukemia protein mediates interferon-based anti-herpes simplex virus 1 effects

Human monocyte chemotactic protein-3 (MCP-3): molecular cloning of the cDNA and comparison with other chemokines

Human monocyte chemotactic proteins-2 and -3: structural and functional comparison with MCP-1

Monocyte chemotactic protein-2 activates CCR5 and blocks CD4/CCR5-mediated HIV-1 entry/replication

C-C chemokines, but not the C-X-C chemokines interleukin-8 and interferon-gamma inducible protein-10, stimulate transendothelial chemotaxis of T lymphocytes

Cytokine patterns in a comparative model of arenavirus haemorrhagic fever in guinea pigs

The role of the vascular endothelium in arenavirus haemorrhagic fevers

Pathogenesis of arenavirus hemorrhagic fevers

A plasma inhibitor of platelet aggregation in patients with Lassa fever

Effective vaccine for lassa fever

The roles of protein C and thrombomodulin in the regulation of blood coagulation

Local overexpression of thrombomodulin for in vivo prevention of arterial thrombosis in a rabbit model

Thrombomodulin-protein C-EPCR system: integrated to regulate coagulation and inflammation. Arteriosclerosis

Thrombomodulin improves early outcomes after intraportal islet transplantation

Thrombomodulin RNA is destabilized through its 39-untranslated element in cells exposed to IFN-gamma

Pathophysiology and treatment of Lassa fever

Decreased neointimal formation in Mac-1(2/2) mice reveals a role for inflammation in vascular repair after angioplasty

Organization and chromosomal localization of the human platelet-derived endothelial cell growth factor gene

The growth promoting effects of bFGF, PD-ECGF and VEGF on cultured postimplantation rat embryos deprived of serum fractions

Pathogenesis of Lassa fever in cynomolgus macaques