| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Journal of Virology, August 2007, p. 7960-7973, Vol. 81, No. 15
0022-538X/07/$08.00+0 doi:10.1128/JVI.00536-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Institute of Human Virology, University of Maryland Biotechnology Center, Baltimore, Maryland 21201,1 Virginia Bioinformatics Institute at Virginia Tech, Blacksburg, Virginia 24061,2 Division of Pathology, Walter Reed Army Institute of Research, Silver Spring, Maryland 209103
Received 14 March 2007/ Accepted 6 May 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Our experiments on rhesus macaques were designed to model human disease following a lethal dose of virus. Needlestick injuries estimated to deliver approximately 102 to 104 infectious particles from patient blood that contains 105 to 107 infectious particles per milliliter have been described previously (62). In both human beings and macaques, intravenous exposure is followed approximately 1 week later by fever (5). In the monkey model, viremia is detectable in blood approximately 4 days after the initial exposure (22, 48). The blood later exhibits elevated liver enzymes and coagulation defects (18, 27). Additional disease signs include petechial rashes, pulmonary edema, cardialgia, renal failure, hemorrhage, and hypoxia, leading to death approximately 2 weeks after virus exposure (46, 62).
Although disease progression in the monkey model is often monitored by liver biopsy and bronchoalveolar lavage (47), venipuncture is less invasive and more practical for monitoring human beings. DNA microarray technology makes it possible to investigate the expression of thousands of genes in a blood sample simultaneously, and profiles of these genes may reveal prognostic biomarkers and insights into the molecular events causing viral hemorrhagic fever.
We characterized changes of gene expression in peripheral blood mononuclear cells (PBMC) isolated at different times from rhesus macaques that were uninfected, infected with the virulent strain LCMV-WE, or infected with the nonvirulent strain LCMV-Armstrong (LCMV-Arm). Although arenaviruses replicate well in periarterial spaces, liver, spleen, and peripheral lymphoid organs soon after inoculation (46, 65), they cannot infect most PBMC (primarily lymphocytes and immature monocytes), and hence viremia does not appear until day 4 in this model. Our goal was to identify molecular events in blood that could (i) serve as markers for disease progression and (ii) discriminate between virulent and nonvirulent infection. We characterized two major stages, previremic and viremic, and compared samples to uninfected samples or samples of blood from monkeys infected with LCMV that did not have disease. This analysis of the primate transcriptome after in vivo LCMV infection reveals critical molecular events that can be related to disease progression.
| MATERIALS AND METHODS |
|---|
|
|
|---|
LCMV infection of rhesus macaques is uniformly lethal at a dose of 103 PFU of LCMV-WE but shows no pathology after a dose of 103 to 108 PFU of LCMV-Arm (22, 46, 47, 68). Infection in macaques was determined by plaque assay of infectious particles in plasma, by infectious center assay of PBMC, and by reverse transcription-PCR (RT-PCR) to detect viral RNA in tissues as described previously (22, 43, 46, 48).
Flow cytometry analysis.
Fluorescence-activated cell sorter analysis was performed on monkey PBMC samples to identify percentages of cell subsets as described elsewhere (reference 68 and J. D. Rodas, C. Cairo, T. Ruckwardt, M. Djavani, J. C. Zapata, J. Bryant, C. D. Pauza, I. S. Lukashevich, and M. S. Salvato, submitted for publication). Overall, there were no dramatic changes in the circulating blood cell populations of rhesus macaques. Depletion of NK and 
T cells (less than 5% of blood cells) was detected in both LCMV-WE and LCMV-Arm samples.
RNA preparation and GeneChip hybridization. Total RNA was prepared from LCMV-infected or uninfected PBMC subjected to TRIzol (GIBCO-BRL) extraction and purification using an RNeasy system (QIAGEN, Valencia, CA) according to the manufacturer's instructions. A QIAGEN RNase-free DNase supplement kit was used to ensure that the RNA had no DNA contamination. RNA quantities were determined by measuring the absorbance at 260 nm. All RNA samples were checked for both quality and quantity on an Agilent 2100 BioAnalyzer system (Agilent Technologies, Palo Alto, CA). RNA that passed this initial quality control screen was then labeled according to the standard target labeling protocols provided by Affymetrix and hybridized to the GeneChip human genome U133 Plus 2.0 array (Affymetrix, Santa Clara, CA) as described by the manufacturer (http://www.affymetrix.com). This chip fully covers the human genome with 54,000 probe sets that represent approximately 22,000 genes.
Microarrays were washed by using a GeneChip fluidics station (Affymetrix). Staining with R-phycoerythrin-streptavidin mix (Molecular Probes) was followed by antibody amplification with a biotinylated streptavidin antibody (Vector Laboratories). Scanning was carried out by using a GeneArray scanner (Hewlett-Packard). The use of the Affymetrix human genome microarrays for monitoring transcriptome changes in rhesus macaque tissues has been validated by other studies (1, 69, 77).
Of the 24 hybridizations used for this analysis, 3 were repeated to determine that the chip-to-chip hybridization mean variation was less than 5%. This agrees with a study comparing different commercially available microarray platforms in which Affymetrix chips had a mean coefficient of variation of less than 5%, whereas the next closest commercial platform had a mean of 10% (71).
Microarray data analysis. Microarray data analysis was performed using the Array Data Analysis and Management System (ADAMS), currently being developed at VBI (http://pathport.vbi.vt.edu). The system uses publicly available tools for analysis of the data. Briefly, raw probe intensities were normalized and summarized using a robust multichip average of G+C content algorithm (gcRMA algorithm) (79). The detection calls (present, marginal, or absent) for each probe set were obtained using the mas5calls function in the Affy R package (29). For paired comparisons, only genes with at least one present call across the compared samples were included. A total of 24 samples were used to generate the microarray data (http://www.ncbi.nlm.nih.gov/geo/). The data from the 24 samples were grouped into four as follows to perform the analysis of variance (ANOVA): uninfected controls (11 samples), samples infected with the nonvirulent strain LCMV-Arm (5 samples), samples infected with the virulent strain LCMV-WE and showing the previremic response (4 samples), and samples infected with LCMV-WE and showing viremic response (4 samples). The P values were corrected using the Benjamini and Hochberg false-discovery-rate test (6).
Pairwise comparisons were also performed to calculate the changes (n-fold) using single pre- and postinfection samples from the same animals. A total of 6,416 genes with changes (n-fold) greater than or equal to two were used for further functional characterization of the genes. Among the 6,416 genes, we identified a total of 3,520 genes with P values from the ANOVA of
0.05 and 1,843 genes with P values of
0.01. The microarray data and the related experimental information from this study can be accessed at http://www.ncbi.nlm.nih.gov/geo/, platform number GSE5790.
Quantitative real-time reverse transcriptase PCR. Selected genes from the microarray analysis were subjected to validation by quantitative real-time RT-PCR to determine the extent to which gene expression was up- or downregulated as a result of infection. Monkey and human gene-specific primers were used to validate expression levels for selected host genes. The PCR primer pair sequences will be supplied upon request. Briefly, total RNA (1 µg) from each sample was used to generate a cDNA template using the iScript cDNA synthesis kit (Bio-Rad, Hercules, CA). The total reaction volume was 20 µl. The reaction was incubated as follows in a Tetrad thermocycler (MJResearch, Waltham, MA): 5 min at 25°C, 30 min at 42°C, 5 min at 85°C, and a hold at 4°C. cDNA products were diluted 1:10 in diethylpyrocarbonate-treated water. The incubations also were performed on controls with no RNA template and with the reverse transcriptase enzyme omitted. The PCR primer pairs were designed based on previously published sequences (GenBank) and using the Oligo 6 primer design software (Molecular Biology Insights, Cascade, CO).
The PCR used cDNA and Taq DNA polymerase (Invitrogen, Carlsbad, CA) and previously described conditions (21, 22). Each gene amplicon was purified by using a MiniElute PCR purification kit (QIAGEN). The purified amplicon for each gene was quantified on an agarose gel and with a GeneQuant Pro spectrophotometer (Amersham Biosciences, Piscataway, NJ). These purified amplicons were used further to optimize the real-time PCR conditions and to generate the standard curves in the real-time PCR assay. Primer concentrations and annealing temperatures were optimized for the iCycler iQ system (Bio-Rad) for each set of primers by use of the system's gradient protocol. PCR efficiencies were maintained at 100% for each primer set during optimization and also during the real-time PCR of sample cDNA.
RT-PCR products were analyzed by 1% agarose gel electrophoresis. The amplification plots were expressed as threshold cycle (CT) values. The CT is the reaction cycle at which PCR products or amplicons reach a threshold level of detection: the lower the CT value, the more abundant the substrate. CT values were normalized by using GAPDH (glyceraldehyde-3-phosphate dehydrogenase) or rRNA 18S amplicons as standards. Dissociation analysis of the PCR products was used to confirm specificity.
Protein determination using SearchLight proteome arrays. Plasma samples were collected from uninfected and LCMV-infected rhesus macaques and spun down to separate the debris, and then the supernatants were stored in 0.5-ml aliquots at 70°C and shipped frozen for analysis on sandwich enzyme-linked immunosorbent assays (ELISAs) on SearchLight proteome arrays (Pierce, Boston, MA).
Proteomic array data analysis. Proteomic array data were calculated as mean (± standard error of the mean) or median (range) values. Statistical significance between groups was determined using a one-way repeated-measures ANOVA. This accounts for individual differences among animals, thus removing animal variability from the analysis. After ANOVA, the Student t test was used to determine statistical significance between groups, with a maximum ANOVA P value of 0.05 considered as significant.
Ingenuity Pathway Analysis. The Ingenuity Pathway Analysis software tool was used to analyze 400 genes that discriminated between virulent and benign infections with changed P values of <0.05 in order to get a better idea of the most likely downstream events.
| RESULTS |
|---|
|
|
|---|
12%) were identified as differentially expressed genes (in Table SI in the supplemental material, these genes are arranged in the order of magnitude of differential expression). We use the word "genes" instead of "probe sets" herein for simplicity, even though each probe set actually identifies a subportion of a gene. We chose to use PBMC rather than whole blood because we wanted our results to be comparable to results from in vitro studies using PBMC.
|
Changes in circulating blood cell populations did not significantly affect blood gene expression profiles.
Using flow cytometry (Rodas et al., submitted), we analyzed fluctuations of CD8+ (Tc) cells, CD4+ (TH) cells, CD20+ (B) cells, CD3-CD16+ (NK) cells, CD3+ CD16+ (NK T) cells, and CD3 
TcR+ cells (
T lymphocytes). By day 2, we detected a 20 to 70% decrease in NK and 
T cells that lasted through day 7, without significant fluctuations in the other cell subsets (Rodas et al., submitted). We examined the current transcriptome data set for fluctuations in CD16 and 
TcR gene expression and saw that they remained below our twofold cutoff throughout the study (see Table SI in the supplemental material). This is reasonable, since the NK and 
T-cell populations constitute less than 5% of the total circulating PBMC, and their depletion during infection would not be expected to have a large impact on gene expression.
Gene expression profiling that discriminates between virulent and nonvirulent virus infection. After all results were normalized to those for uninfected controls, PBMC gene expression profiles showed significant differences between virulent and nonvirulent infections (Fig. 2). We determined that the levels of approximately 478 genes discriminated between nonvirulent and virulent LCMV infections (Fig. 2). Although data for individual animals are shown, statistical significance for the 478 genes was achieved by pooling results from eight animals (LCMV-WE infected) for comparison with results from three animals (LCMV-Arm infected). Genes that were downregulated during virulent LCMV-WE virus replication were relatively unchanged or upregulated with nonvirulent LCMV-Arm infection and vice versa (genes up with LCMV-WE were unchanged or down with LCMV-Arm). The 478 genes that differ most between virulent and nonvirulent infections are listed in Tables 1, 2, and 3 and in Tables SI to SV in the supplemental material. Tables 1, 2, and 3 are grouped temporally as previremic (Table 1), upregulated viremic (Table 2) and downregulated viremic (Table 3). The expression levels for some of these genes were verified by real-time PCR or by SearchLight ELISA (Table 4).
|
|
|
|
|
Functional grouping of genes regulated by LCMV infection. Genes regulated by LCMV were first grouped based on their cellular functions. Over 65% of the 6,416 differentially regulated genes were categorized as immune, proinflammatory, or stress response genes (Fig. 3; also see Tables SI to SV in the supplemental material for additional online data). A number of genes encoding transcription control factors (approximately 50) were downregulated, while 25 of them were upregulated in LCMV-WE-infected PBMC. The transcription factor human homeobox gene (HHEX) (see Table SIV in the supplemental material), also known as proline-rich homeobox (PRH), was not differentially expressed in PBMC, contrary to what we described previously for liver (22).
|
Temporal grouping of genes: early host response to the virulent strain of LCMV. Infection by the virulent LCMV-WE strain at the earliest (previremic [days 1 to 3]) time points revealed a number of downregulated genes greater than that of upregulated genes in macaque PBMC (Fig. 1 and 2; also see Table SI in the supplemental material). Prominent early downregulated genes included PTGS2, G protein-coupled receptor (GPR), EGR1 and -2, nuclear receptor subfamily 4 group A (NR4A2), chemokine (C-X-C motif) receptor 4 (CXCR4), insulin receptor substrate 2 (IRS2), and response gene to complement 32 (RGC32). Prominent genes that were upregulated during the previremic phase of the virulent infection included interferon (IFN)-inducible genes (e.g., CIG5, IFI44, IFIT4, GBP1, OAS2, OAS-L, and MX2), signal transducer and activator of transcription 1 (STAT1), dehydrogenase/reductase 9 (DHRS9), myristoylated alanine-rich protein kinase C substrate (MARCKS), and S100 calcium binding protein P (S100P). Profiling of the virulent infection demonstrated significant differences from the nonvirulent LCMV-Arm infection even before the onset of viremia (Table 1; also see Table SV in the supplemental material). Most of the differentially expressed genes remained unchanged after day 4 of infection (see Table SI in the supplemental material).
Viremic stage of response to the virulent strain of LCMV.
The host transcriptional response to LCMV-WE viremia was dominated by immune response, stress response, and proinflammatory genes. We also observed specific upregulation of genes involved in antiviral responses, including IL-15 (4-fold), as well as in a number of type I and II IFN-responsive genes, such as GBP1 (>37-fold), GBP2 (>12-fold), G1P2 (>15-fold), MX1 (>15-fold), and MX2 (>16-fold) (see Table SI in the supplemental material). Mx antiviral proteins belong to the family of GTPases and are induced preferentially by IFN-
/ß (25). Viperin (cig5), encoding an antiviral protein, is also upregulated by LCMV-WE infection but is downregulated with nonvirulent LCMV-Arm infection. Viperin is an IFN-inducible protein and is also upregulated by human cytomegalovirus and vesicular stomatitis virus infections (7, 15). Viperin and Mx proteins have been shown to impair influenza virus, vesicular stomatitis virus, rabies virus, and hepatitis C virus replication (15, 33, 40, 80). Increases in the expression of genes encoding neutrophil and/or monocyte/macrophage-targeted chemokines and chemokine receptors were observed by day 6 of the viremic stage. Interestingly, the gene which proved to be most strongly induced (119-fold) among all 6,416 differentially regulated transcripts was the IFN-stimulated chemokine CXCL10, also known as IP-10, which has a critical role in recruiting T-helper cells and other lymphocytes (see Table SI in the supplemental material). Expression of CCL3 (also known as MIP-1
), CCL4 (also known as MIP-1ß), CCL20, and the CCL3/CCL4 receptor CCR5 were all increased by 6 days postinfection. By contrast, the related CCL5 (RANTES) was downregulated in the previremic stage and was then unremarkable in the viremic stage. Cytokines, chemokines, and chemokine receptors taken together represented the strongest induction among all functional subclasses of immune/stress/inflammatory response genes (Tables 2 and 3).
Confirming microarray data by RT-PCR. In previous publications, we confirmed viral modulation of mRNA levels by use of quantitative real-time RT-PCR as well as measures of protein expression (21, 22). Real-time PCR was carried out using cDNA from LCMV-Arm and -WE-infected animals and from preinfection control animals to determine gene expression levels. Modulation of the specific mRNA for 16 genes relative to the expression of GAPDH was assessed at various times after infection. These genes were selected from different categories that contained the most highly up- or downregulated genes. Microarray gene expression levels were compared to those obtained by real-time PCR, using the same monkey blood samples that were used in the microarray hybridization. The microarray data were correlated with those obtained by quantitative real-time RT-PCR (Table 4), and sustained up- or downregulation was often also seen by PCR. The magnitude of up- or downregulation was usually more extreme in PCR than in microarray, since PCR data are accurate over a broader range.
Plasma chemokine and cytokine protein expression levels in LCMV-infected rhesus macaques.
To determine whether plasma chemokine and cytokine concentrations were significantly different between virulent LCMV and nonvirulent virus throughout the experiment, samples from infected rhesus macaques were analyzed by multiplexed sandwich ELISA. Plasma protein levels of LCMV-infected animals were compared to those of the control preinfection animals. The levels of plasma cytokines remained unchanged from the preinfected levels during the previremic stage of infection. After day 4 of infection (viremic stage), IFN-
, IP-10, and IL-8 levels were markedly upregulated in macaque plasma and remained elevated up to day 7 (Fig. 4A).
|
-inducible protein 10 (IP-10 or CXCL10) is a CXC chemokine that binds to the CXCR3 receptor on activated T lymphocytes and NK cells (73). We measured plasma levels of IP-10 in animals infected with virulent or nonvirulent LCMV strains. IP-10 was elevated markedly in plasma and remained elevated between day 4 and day 7, and levels were positively correlated with IFN-
concentrations and plasma viremia beginning 4 days after infection (Fig. 4B). Monocyte chemoattractant protein 1 (MCP-1), known for its ability to loosen endothelial cell junctions (41) and to participate in the breakdown of the blood-brain barrier (75), had a similarly high plasma profile in samples from virulently infected monkeys but not in samples from uninfected or LCMV-Arm-infected monkeys (Fig. 4B). Analysis of secreted proteins in plasma showed decreased levels of VEGF, IL-1ß, and tumor necrosis factor alpha (TNF-
) (Fig. 4C to E). Since VEGF is one of the principal factors responsible for inducing blood mononuclear cell COX-2 expression (34), it is reasonable that the downregulation of VEGF (Fig. 4E) should be concomitant with the downregulation of COX-2. A pathway analysis revealed the connections between a few important gene products, such as the transcription factors EGR1 and -2 and the COX-2 gene PTGS2 (Fig. 5). Additional nodes in this pathway included two gene products (IL-1R and epidermal growth factor receptor [EGFR]) that were prominent in a "kinomics" study by Bowick et al. (10). The severe downregulation of EGR1 and EGR2 during the virulent infection is directly linked to the downregulation of PTGS2, VEGF, IL-1R, and EGFR.
|
| DISCUSSION |
|---|
|
|
|---|
During days 2 to 4, the global impact of LCMV-WE infection was to reduce host gene expression. With the onset of viremia (after day 4), the trend was reversed as host genes were turned on in response to the rising viral burden (Fig. 1). This finding resembles an observation made by the Shenk laboratory about human cytomegalovirus in tissue culture: gene expression was predominantly downregulated in infected cells, but UV-inactivated virus caused upregulation of a large number of mRNAs, many of which contained IFN-responsive genes and encoded proinflammatory cytokines (11). Abundant viral replication, including the production of noninfectious particles, appears to drive the expression of IFN-responsive genes. A study of influenza-exposed cell cultures from the Katze laboratory revealed that a substantial proportion of the transcriptome response to viral infection can be driven by viral antigens in the absence of replication (30). Thus, the broad pattern of gene upregulation during the viremic phase of infection may reflect the host response to viral particles with an ensuing activation of IFN to control the spreading infection.
A recent profiling study with influenza virus-infected macaques demonstrated that gene expression arrays can be used to detect infection after exposure but prior to the onset of symptoms or virus shedding (3). In our study, we identified 91 genes expressed before the onset of viremia that differed significantly between virulent and nonvirulent infections (Table 1). Several of these early host responses have also been noted for murine LCMV infections. For example, CCR2, the MCP-1 receptor that recruits monocytes to inflamed sites, was upregulated in LCMV-WE-infected monkeys and in LCMV-infected mice but not in vesicular stomatitis virus-infected mice (57). Murine brains persistently infected with LCMV upregulated STAT1 and IFN-stimulated genes (39) that are prominent early in the monkey virulent infection; however, the murine tissues upregulated the major histocompatibility complex and ß2-microglobulin genes, but these are downregulated in monkeys. The monkey fibronectin gene was another gene upregulated early; similarly, in murine LCMV hepatitis, fibronectin deposits accompany viral antigen in hepatic sinusoids (42). Early growth response genes EGR1 and EGR2 are downregulated early in virulent monkey infections. In acute murine LCMV infections, EGR genes are critical for the development of a vigorous cytotoxic-T-lymphocyte response (72), and early downregulation of EGR1 and -2 in monkeys could be a mechanism for limiting the development of LCMV-specific cell-mediated immunity. Heat shock protein gp96 is involved in the stable antigen presentation of mature LCMV NP (4), and the gp96 gene, TRA1, was downregulated early in infected monkeys, indicating another mechanism for limiting antiviral immune responses. The resemblance between murine responses to LCMV infection and previremic primate responses supports the role of these transcripts as early indicators of infection.
The IFN response pathways constitute an important antiviral response during LCMV infection. Many IFN-inducible genes, including the MX2 gene, were highly induced by day 2 in LCMV-WE-infected PBMC (Fig. 3D). Murine studies showed that LCMV triggers early and high levels of IFN production in nonplasmacytoid dendritic cells (20). However, in our monkey studies plasma levels of type I IFN were unremarkable, though type II IFN levels did rise with viremia. The IFN-inducible genes are notable as a group because they are more readily upregulated by inactivated viral particles than by infectious particles (references 11 and 51 and M. S. Salvato, unpublished data). Reduced expression of IFN-inducible genes in the LCMV-Arm-infected monkeys is due either to lower levels of viral replication or to the successful suppression of the IFN pathway by LCMV-Arm.
TNF-
/NF-
B-responsive genes were suppressed during smallpox infection and strongly upregulated during the viremic phase of Ebola virus infection in a monkey transcriptome study by Rubins et al. (69). In our study, many of these same genes, categorized as IFN responsive or inflammatory, were also upregulated (with a few notable exceptions) during the viremic stage of virulent infection and were relatively diminished in expression during nonvirulent infection. The exceptions include IL-8, VEGF, IL-1ß, and CXCR4, which were downregulated, while most inflammatory genes (CXCL10, CCR5, CCR2, and TLR8) were upregulated in a virulent (lethal) infection. We do not know how these genes would behave in an animal that survives virulent disease. In a study of people infected with Lassa fever virus, survivors had higher levels of the proinflammatory cytokines (IL-8 and IP-10) than the fatalities (49), indicating that these cytochemokines could be essential components of host defense.
The notion that activated inflammatory pathways can avert lethal viral disease has been pursued in several animal studies. In a mouse model for Ebola, inoculation routes that elicited earlier and higher-level type I IFN promoted survival (50). Treatment of Marburg virus-infected monkeys with IFN-
prolonged survival but could not save the animals (38). In mice, the use of CpG oligonucleotides to stimulate type I IFN via TLR9 reduced arenavirus choriomeningitis (58). Unfortunately, type I IFNs not only induce antiviral responses but also have deleterious immune-suppressive activities. For example, type I and type III IFN responses synergize in suppressing the CD4 proliferative responses to respiratory syncytial virus infection (14). Initially, cell proliferation is blocked by IFN-mediated delays to cell cycling, and once cells have proliferated, IFN promotes apoptosis of activated lymphocytes (24). Our current hypothesis is that early expression of inflammatory genes stimulates acquired immunity but that high expression of inflammatory genes during the viremic phase contributes to the immune suppression of virus-specific host responses.
In support of this hypothesis, we note that innate effectors such as TLR1, TLR2, TLR4, and TLR8 are upregulated late during the viremic phase. TLR10, which is expressed in B cells, respiratory epithelia, and dendritic cells (32), is downregulated by late viremia along with the majority of HLA genes, which are downregulated in LCMV-WE-infected tissues compared to what is seen for LCMV-Arm-infected tissues. Reduced HLA expression likely causes a defect in antigen presentation that could contribute to pathogenesis.
Cell culture studies have been extremely useful in identifying host responses that are directly connected to the functions of viral gene products. For example, the NS proteins of respiratory syncytial virus and the NP protein of LCMV suppress IRF-3 and block NF-
B-promoted expression of IFN-ß, IFN-
, and proinflammatory chemokines like IL-8 and TNF-
(53, 74, 76). However, the IFN-suppressive activity of LCMV NP is found in both virulent and benign LCMV isolates, so its connection to pathogenesis is unclear. Cell culture studies with the LCMV Z gene product also failed to find disease-related mutations but did identify important virus-host interactions (12, 16, 22, 59). Our working hypothesis is that the toxic activities of viral gene products like NP or Z must reach a threshold that is dependent on polymerase activity to cause disease.
At first glance, cell culture studies that compare virulent and nonvirulent viruses do not support our findings about gene expression in vivo. A comparison of Lassa virus and Mopeia virus in cell culture showed that TNF-
and IL-8 are downregulated by the virulent virus in relation to the nonvirulent virus (45). In monkeys, although IL-8 is downregulated, TNF superfamily genes are upregulated, and several other NF-
B-driven responses rise during the viremic stage of the virulent disease. Comparison of virulent and nonvirulent arenavirus isolates, Pichinde virus isolates P18 and P2, respectively, in a monocytic cell line predicted the upregulation of gene products (Nurr1, Fez, Striatin, and PTEN) (9); however, the mRNAs for those products were downregulated in our study. A study using liver cell cultures infected with Marburg and Ebola viruses related virulence to increased suppression of the type I IFN and inflammatory responses, downregulation of coagulation genes, and upregulation of immune response genes (35). In our monkey studies, these three gene categories during the viremic phase of LCMV-WE disease look just the opposite of how they appear during primate infection with virulent Ebola virus (69). Our view is that the cell culture studies accurately depict virus infection in the microcosms of cells that are initially infected in vivo, but those results are drowned out by the responses from uninfected cells. Thus, the previremic stage in vivo reflects a small amount of virus in an overwhelmingly large population of uninfected cells, whereas the viremic stage, similar to what is seen in cell culture studies with high virus input, is driven by virus but dominated by host efforts to suppress viral replication.
A recent paper by Bowick et al. (10) compares guinea pigs infected with virulent and benign arenaviruses, Pichinde virus isolates P18 and P2, respectively. Using "kinomics" to assess kinase activity in peritoneal macrophages, this group showed that virulent infection significantly downregulated kinase activity as well as the activity of several signaling pathways. Their findings at the protein level are corroborated by our study, in which several mRNAs with products involved in signal transduction are downregulated in the virulent infection (Fig. 3F). Their most prominent findings, that virulent infection downregulates epidermal growth factor receptor (EGFR) and IL-1 receptor (IL-1R), are corroborated in a pathway analysis of our most prominent early gene downregulations (Fig. 5). Similar to our study, their study found major early responses that are not obviously related to results from virus-infected cell cultures (9). For example, we saw early downregulations of PTGS2 (encoding COX-2), NR4A2 (encoding nuclear receptor NURR1), and GPR109B (encoding G-protein coupled receptor 109B) that were not seen in PBMC cultures exposed to LCMV-WE or LCMV-Arm (J. C. Zapata, M. M. Djavani, and M. S. Salvato, unpublished observations).
Our observation of decreased PTGS2 (encoding COX-2) is remarkable in these studies. COX-2 is the key enzyme in prostaglandin biosynthesis; it is upregulated by many viruses (13, 17, 36, 37, 66, 67, 78), and cytomegalovirus is so dependent on COX-2 for replication that it encodes a COX-2 homologue (70). COX-2 expression can be upregulated by TNF-
, IL-1ß, or a variety of growth factors like VEGF and EGRF (2), but virulent arenaviruses in cell culture usually downregulate such products of NF-
B/RBP-mediated transcription (22, 26, 45). In fact, the genes for IL-1ß, VEGF, and IL-8 are all downregulated in virulent monkey infections compared to what is seen for nonvirulent monkey infections. Downregulation of COX-2 occurs during overdose of drugs such as rofecoxib (Vioxx), which are known to deplete prostacyclins, destroy clotting, and cause hemorrhage (63). Our findings with COX-2 corroborate studies by Fisher-Hoch and Cummins et al. that show a decrease in prostacyclin in monkeys with Lassa fever (18, 28). The decrease in PTGS2 expression during virulent disease is not directly attributable to viral gene functions but can provide a mechanism for the coagulation defects (18, 55), cardiac dysfunction (61, 64), and capillary leakage (55) that are hallmarks of arenavirus hemorrhagic fever.
In summary, microarray analysis of blood from rhesus monkeys revealed how infection with a virulent arenavirus regulates gene expression in an animal model. The early stage of the host response to virulent LCMV infection is dominated by innate immune responses based on complement, IFN, and proinflammatory processes. The alteration in the transcriptome during viremia strongly suggests that the higher level of replication of the virulent virus contributes to the pathogenicity of LCMV-WE-associated hemorrhagic fever. Of 54,000 probe sets in a full human genome array, 12% were differentially regulated in blood during primate infection with the virulent LCMV-WE. Roughly 400 of these genes differed in expression between LCMV-WE and LCMV-Arm infections, and of these genes, 91 differed during the early or previremic stage and could potentially be biomarkers of disease that are capable of classifying the disease prior to the onset of viremia. The characterization of gene expression changes provides a dynamic view of the host-pathogen interaction and corroborates a recent publication (10) on the guinea pig model for Lassa fever.
| ACKNOWLEDGMENTS |
|---|
We are grateful for Susan Martino-Catt and Clive Evans for helpful comments and technical support.
| FOOTNOTES |
|---|
Published ahead of print on 23 May 2007. ![]()
Supplemental material for this article may be found at http://jvi.asm.org/. ![]()
| REFERENCES |
|---|
|
|
|---|
B and RBP-J
by arenavirus infection of macrophages in vitro and in vivo. J. Virol. 76:1154-1162.
gene expression. J. Med. Virol. 59:552-560.[CrossRef][Medline]