Previous Article | Next Article ![]()
Infection and Immunity, September 2005, p. 6091-6100, Vol. 73, No. 9
0019-9567/05/$08.00+0 doi:10.1128/IAI.73.9.6091-6100.2005
Divisions of Communicable Disease and Immunology,1 Retrovirology, Walter Reed Army Institute of Research, Silver Spring, Maryland 20910,3 Division of Emerging and Transfusion Transmitted Diseases, Center for Biologics Evaluation and Research, Food and Drug Administration, Rockville, Maryland 208522
Received 15 February 2005/ Returned for modification 25 March 2005/ Accepted 13 April 2005
|
|
|---|
|
|
|---|
Animal models of malaria are particularly useful in unraveling such protective mechanisms, since they provide flexibility in terms of experimental design and sample collection. The murine malaria parasite, Plasmodium yoelii, has numerous similarities to human P. falciparum and P. vivax and is a useful model for the comparison of differences in gene expression differences between lethal and nonlethal infections. The parasite has evolved two distinct strains, presumably due to the evolution of different pathophysiologic pathways. During the asexual blood stage infection, P. yoelii XNL parasites invade primarily reticulocytes and the infection is self-limiting. The sister strain, P. yoelii XL, invades all forms of murine erythrocytes, and the resultant high parasitemias cause a lethal infection (13, 17, 36). A consequence of infection is a profound anemia, the cause of which is not fully understood but involves both the destruction of parasitized erythrocytes and a dyserythropoiesis which impedes the production of sufficient numbers of newly formed red cells (18, 28).
During murine malaria infection, the spleen becomes the primary site of erythrocyte production and is involved in the removal of both dead parasites and malaria-infected red cells. As red cells are destroyed by cycles of parasite invasion and erythrocyte rupture, the host response to the evolving anemia differs in animals infected with either the lethal or nonlethal P. yoelii strain. Previous histologic studies have shown that the spleens of animals infected with 17XNL parasites erect a cellular barrier that functions to isolate newly produced reticulocytes away from ongoing infection. In contrast, a competent blood-spleen barrier fails to develop in animals infected with the lethal 17XL strain, resulting in phagocytosis of parasitized erythrocytes and late-stage erythroblasts (30-32). The spleen plays critical roles in establishing immune responses to infection, particularly the suppression of B-cell proliferative responses in animals infected with nonlethal forms of P. yoelii (29).
The primary objective of this study was to analyze the transcriptional changes in animals infected with either a nonlethal or lethal variant of P. yoelii malaria in order to discover molecular processes and pathways that correlate changes in gene expression with outcome of infection. Two independent but complementary approaches were used: (i) we studied transcriptional patterns during the course of a nonlethal 17XNL infection as the parasitemia increased and then decreased after reaching a peak parasite density of 50%, and (ii) we compared gene expression changes between the lethal and nonlethal P. yoelii infections at the identical stage of parasite density, which revealed specific gene ontology (GO) functional groups that differed with respect to clinical outcome of infection. We have identified three distinct patterns of global genetic signatures in early infection that are related to the host's response to regulatory and target gene expression of erythropoiesis in response to anemia, metabolic perturbations in the glycolytic enzyme pathway, and B-cell immune responses that distinguish lethal from nonlethal P. yoelii infections.
|
|
|---|
Experimental design. Frozen stocks of 17XNL and 17XL P. yoelii-infected erythrocytes were thawed and used to infect donor mice. On day 0, 24 mice were infected intravenously with 1 x 106 17XNL parasites and 16 mice with 1 x 106 17XL parasites. Six mice were injected with phosphate-buffered saline, and the mean gene expression value from each spleen served as the control baseline reference for all other samples. Parasite density in each animal was checked daily by direct microscopic examination of thin blood films and recorded as the percent infected erythrocytes. Mice infected with the 17XL strain generally succumb to infection 7 to 10 days postinfection, while mice infected with the 17XNL strain reach maximum parasitemia (50%) at approximately around day 10, followed by a rapid clearance of parasites and complete resolution of infection by day 17. For mice infected with 17XNL parasites, six mice were sacrificed at each successive time interval: at 24 h postinfection and on days 4 (5% parasitemia), 9 (25% parasitemia), 10 (50% parasitemia), 14 (25% parasitemia), and 17 (<1% parasitemia). For animals infected with 17XL parasites, four mice were sacrificed at 24 h postinfection and on days 3 (5% parasitemia), 5 (25% parasitemia), and 7 (50% parasitemia) (Fig. 1A). Six noninfected control animals were sacrificed 24 h after all the mice received either the 17XNL or 17XL infection. Spleens were removed, immediately snap frozen, and stored at 70°C until samples were processed.
![]() View larger version (26K): [in a new window] |
FIG. 1. Course of infection and differential gene expression in P. yoelii infection. (A) Kinetics of parasite replication in 17XNL and 17XL infection in mice. Error bars indicate standard deviations. (B) Number of differentially expressed genes (P < 0.005) during specific intervals during infection relative to mean gene expression from six uninfected mice. (C and D) Three-dimensional PCA of average expression from six 17XNL-infected mice at six time intervals (C) and from four 17XL-infected mice at 4 time intervals (D) during infection is shown in relation to gene expression from six uninfected mice (yellow circle). Significant differentially expressed genes were collapsed into three-dimensional vectors to reveal differences in transcript abundance as a function of parasite density.
|
Sample preparation and GeneChip analysis. Preparation of cDNA, in vitro transcription, staining, and scanning of Affymetrix U74Av2 GeneChips containing 12,489 probe sets and 8,305 genes (Affymetrix, Santa Clara, CA) were carried out essentially as described previously (27). Each spleen sample from uninfected and infected mice at each parasite density level was processed and hybridized individually to yield a data set of 58 total GeneChips for analysis. Scanned images were analyzed with Affymetrix MAS 5.0 to create .cel and Pivot data files. Two criteria were used to determine chip quality: the scaling factor determined by Affymetrix MAS 5.0 with the target signal intensity being set to its default (500) and the array outlier percentage determined by dChip version 1.3. GeneChip .cel files were normalized at the probe level by using the robust multichip average (RMA) method (12).
Differentially expressed genes were identified by a Student's t test comparing groups of mice at specified parasite density levels relative to spleen samples from six uninfected control mice, using a cutoff significance level set at a P value of <0.005. This significance value was chosen to minimize the false discovery rate (<1% by significance analysis of microarrays) (26) but to be generous enough to include a larger set of genes with functionally distinct gene expression programs. The data for each gene probe on the array were expressed as the log2 ratio of normalized fluorescence intensity of the sample and the average expression intensity from six uninfected control mice. Hierarchical cluster analysis of microarray data sets was performed using Cluster and Treeview software (http://rana.lbl.gov/EisenSoftware.htm) (10). Self-organizing maps and principal-component analysis (PCA) were performed with GeneLinker Gold software (Predictive Patterns Software, Kingston, Ontario, Canada). To assign gene ontology-annotated terms to differentially regulated sets of statistically significant genes or to groups of coexpressed genes detected in hierarchical clusters, the web-based tools Onto-Express (http://vortex.cs.wayne.edu:8080/index.jsp) (9) and dCHIP (http://biosun/harvard.edu/complab/dChip/) (15) were used.
|
|
|---|
Patterns of gene expression in 17XL and 17XNL infection varied with the intensity of parasitemia (Fig. 1B). PCA was used to reveal large differences in transcript abundance between groups of samples. PCA is an unsupervised clustering method that reduces dimensionality of data from which a large number of variables (genes) are interrelated. PCA identifies trends in the data that define distinct clusters by performing a covariance analysis between factors. Three-dimensional PCA plots often place objects that are similar next to each other and are reflective of the differences in gene expression in mice at various parasite density levels. The greatest change in gene expression relative to baseline in both the lethal and nonlethal infections occurred as the parasite density rose during two consecutive time intervals (from 0 to 5% and from 5 to 25% parasitemia), while significant changes in transcription continued to occur as the parasitemia increased to 50% in 17XL-infected mice (Fig. 1C and D).
Regardless of the ultimate outcome of an infection, transcriptional changes in innate immune genes occur within hours after a host senses an intact or degraded parasite. Within 24 h of infection and prior to the time parasites are detected on blood smears, a number of genes and associated gene ontology functional groups involved in innate immunity, inflammation, defense response, and cytokine signaling are induced in mice infected with either the 17XNL or 17XL strain (Fig. 2A). Such genes include that for tumor necrosis factor alpha (TNF-
) and genes involved in antigen presentation and T-cell activation. A larger set of differentially expressed genes, including cytokine genes (for interleukin-1ß, interleukin-6, gamma interferon, and TNF-
), chemokine genes, interferon-induced genes, and genes critical for adaptive immune responses were similarly induced by both lethal and nonlethal parasites as the parasite density rose to 5% (Fig. 2B). Many transcriptional changes precede the action of inflammatory cytokines such as TNF-
and gamma interferon and have been shown to vary with parasite density and return to baseline upon resolution of infection (16, 19, 25). The changes in expression in genes associated with the defense response in murine malaria were similarly observed in an early infection in nonhuman primates infected with the nonlethal parasite Plasmodium cynomolgi (35). As the infection progressed and the parasitemia rose in 17XL- and 17XNL-infected mice, albeit at different rates, identical patterns of expression for hundreds of differentially expressed genes were observed in both infections, which were independent of whether animals ultimately died or survived their infection (see Fig. S1 in the supplemental material).
![]() View larger version (45K): [in a new window] |
FIG. 2. Concordance in gene expression of immune response genes between lethal 17XL and nonlethal 17XNL infection. The fold change in intensity of differentially expressed genes is depicted as a heat map at 24 h postinfection (A) and at 5% parasitemia (B) for P. yoelii (17XL and 17XNL)-infected mice compared to six uninfected mice. Statistically significant gene ontology groups derived by Onto-Express from sets of up-regulated genes at 24 h postinfection and at 5% parasitemia are indicated.
|
![]() View larger version (29K): [in a new window] |
FIG. 3. Functional annotation of differentially expressed genes in 17XNL infection. (A and C) Pie charts depict the representative number of overrepresented GO terms from genes up-regulated (A) or down-regulated (C) during infection. The size of each slice represents the number of differentially expressed genes assigned to the particular GO functional classification. (B) Gene expression for one representative GO term (heat shock protein activity) is tightly regulated in nonlethal infection and varies directly with intensity of infection. (D) Expression of genes with chemotaxis/chemokine activity is repressed during 17XNL infection and remains suppressed even after resolution of parasitemia. The heat maps indicate the average fold change in gene expression greater than (red) or less than (blue) that in uninfected control spleens at each time point or stage of infection.
|
Differential expression of genes in erythropoiesis. Anemia is a severe consequence of experimental and human malaria infection, and the mechanism comprises both elements of red cell destruction and ineffective erythropoiesis. As the spleen is a major site for erythrocyte production during murine malaria infection, the kinetics of red cell production vis-a-vis parasite clearance is an important element in the host's response to anemia (6, 30). The expression of a number of transcription factors and their target genes associated with erythrocyte function, iron metabolism, heme biosynthesis, and erythrocyte transcriptional regulation was examined. Early in 17XNL and 17XL infection there was a down-regulation of genes that encode erythrocyte membrane proteins (glycophorin A, band 3), the erythropoietin receptor, iron metabolism proteins (transferrin receptors), and heme binding proteins (Fig. 4A). However, as the parasitemia increased from 5% to 25%, the transcript abundance of genes involved in erythropoiesis was consistently and dramatically induced in all mice infected with the nonlethal P. yoelii strain, while gene expression remained at or below basal levels in lethally infected mice relative to uninfected controls (Fig. 4A). Remarkably, all eight of the genes encoding enzymes involved in the biosynthesis of heme, including the rate-limiting enzyme aminolevulinic acid synthase, followed an identical expression profile that was observed for genes encoding red cell membrane receptors and iron metabolism proteins (Fig. 4B). The mean transcriptional activity for mice infected with the 17XNL strain for each of the genes involved in the heme biosynthetic pathway was significantly induced (mean 2.5-fold increase) compared to that for mice infected with the 17XL strain.
![]() View larger version (37K): [in a new window] |
FIG. 4. Differential gene expression in erythropoiesis. (A) Centroid plot of mean expression of 13 genes involved in iron transport and erythrocyte membrane proteins, with heat maps showing fold change in transcript abundance (columns represent individual mice). (B) Relative gene expression for genes representing the enzymes of the heme biosynthetic pathway at 25% parasitemia in nonlethal (NL) and lethal (L) malaria. (C) Cartoon representation of gene expression for regulatory transcription factors critical for erythropoiesis (adapted from reference 20). Red or blue colored boxes and ovals represent significantly induced or repressed gene expression, respectively. Gray objects represent no change in gene expression relative to uninfected control spleens.
|
Differential expression of genes regulating glycolysis. Metabolic acidosis resulting from the accumulation of lactate is an important prognostic indicator for the severity of malaria and has been implicated in the pathogenesis of human P. falciparum and murine Plasmodium berghei infection (8). The transcriptome of glycolytic pathway genes in 17XNL- and 17XL-infected mice at different stages of infection was examined. Each of the 10 genes encoding enzymatic reactions involved in the glycolytic pathway was significantly up-regulated early in 17XL and 17XNL infection (Fig. 5), and mice infected with the lethal 17XL strain showed increased gene expression for all the enzymes throughout infection. Conversely, in the nonlethal infection, a transient initial increase in expression was followed by a dramatic suppression of glycolytic enzyme transcription in the interval of time that the parasitemia rose from 5% to 25% (Fig. 5).
![]() View larger version (15K): [in a new window] |
FIG. 5. Expression of glycolytic enzyme genes. The log2 fold changes of the mean expression of glycolysis genes from 0 to 25% parasitemia for the 17XNL (n = 6) (Fig. 4A) and 17XL (n = 4) (Fig. 4B) infections are shown. The genes are listed, from top to bottom, in the order of the glycolytic enzyme pathway.
|
![]() View larger version (42K): [in a new window] |
FIG. 6. Expression of B-cell proliferation and plasma cell differentiation genes. (A) Hierarchical clustering of genes significantly differentially expressed between 17XNL and 17XL infection. The red-green matrix represents the normalized expression pattern for each gene across the samples. The intensity of the red or green shaded boxes represents the highest or lowest relative expression, respectively. Overrepresented GO terms identified by dCHIP in coexpressed clusters of genes are indicated. NL, nonlethal; L, lethal. (B) Comparisons in gene expression ratios for B-cell activation and proliferation genes. The color matrix plot shows relative expression across spleens from 17XNL and 17XL infection at 25% parasitemia. (C) Color matrix plot of a subset of immunoglobulin genes from the hierarchical cluster at 25% parasitemia, showing relative overexpression during nonlethal infection.
|
|
|
|---|
Biological and molecular processes linking sets of up-regulated genes to gene ontology-annotated terms were discovered from gene clusters identified through hierarchical clustering and self-organizing maps. The overrepresented GO terms assigned to sets of differentially expressed genes during specific phases of infection suggested that the functions of coexpressed genes may also be coregulated. The networks of cis-regulatory elements and transcription factors that regulate induction or silencing of gene expression during malaria infection play critical but as-yet-undefined roles in coordinating the activation of gene expression and may help explain why transcriptional activation continues after the host successfully clears the infection. The coordination of gene expression with gene regulation is tightly controlled and influences the expression of many genes involved in the pathological manifestations of disease. In particular, the differential expression of transcription factors and their target genes in regulating erythropoiesis illustrates how expression profiling differentiates between groups of animals that die (24; this paper) or recover from infection and associates that outcome with the host's response to severe malaria anemia. Such transcription factors and cis-regulatory elements, singly or in molecular complexes, regulate red cell proliferation and expression of erythrocyte genes for band 4.2, uroporphyrinogen III synthase, transferrin receptor 2, glycophorin A, and hemoglobin (3, 5, 20, 34).
As erythropoiesis is not limited to the spleen, future investigations should include genome-wide expression profiling of tissues such as the bone marrow and liver, in addition to correlating the expression patterns of regulatory proteins (transcription factors) with their target genes to improve our understanding of the molecular interactions that may predict prognostic signatures for severe malarial anemia in young children.
Physiologic adaptations to malaria infection have indicated that lactic acidosis resulting from the conversion of pyruvate to lactate is implicated in the pathogenesis of human P. falciparum and lethal murine P. berghei infection and is a prognostic indicator of severe disease (8). Sexton et al. have recently reported on the up-regulation of nine glycolytic pathway genes in the spleens of mice infected with P. berghei that favor the conversion of pyruvate to lactate and the down-regulation of counterregulatory genes which shunt or divert products from this pathway (24). Other reports have shown that acidosis correlates with increased parasitemia, progression of disease, and increased glycolytic activity (7, 8, 23). In the study presented here, mice infected with lethal 17XL P. yoelii parasites showed increased gene expression of all the enzymes present in the glycolytic pathway throughout infection. Conversely, in the nonlethal infection, transient increases in gene expression were followed by a dramatic reversal in gene transcription in the interval of time when the parasitemia rose from 5% to 25%. Whether the divergence in patterns of expression between lethal and nonlethal infections reflects the cause or effect of mechanisms that influence death or survival is unknown. However, the observation that the reversal in transcriptional activation of glycolytic pathway genes in nonlethal infection precedes the peak of maximal parasite density suggests that expression patterns may have prognostic significance as predictors of clinical outcome.
Signature patterns of gene expression in immune response genes that distinguished 17XNL from 17XL infection were associated primarily with B-cell proliferation and immunoglobulin production. Antibody responses and plasma cell production in lethal P. berghei and P. yoelii infection are suppressed (21, 30), while nonlethal P. chabaudi and P. yoelii infections induce a large expansion of plasma cells within the red pulp of the spleen (1, 2, 32). The suppression of gene expression of CXCL13, a chemokine important for B-cell homing to germinal centers (14), observed here coincided with the suppression observed for other B-cell proliferation genes, such as Lsp1 and Dock2, necessary for cytoskeletal rearrangements and lymphocyte motility and proliferation (4, 11), while genes involved in plasma cell differentiation were induced. Critically, the relative contribution of antibody-mediated immune mechanisms that appear to control the 17XNL infection studied here in an immunologically intact animal may mask or impede the transcription of genes that are crucial for resolving infections primarily by cell-mediated immune mechanisms (22). In addition, care must be taken not to overinterpret the differences in B-cell and plasma cell gene expression in animals infected with a lethal compared to a nonlethal infection, since the rate of 17XL parasite replication may outstrip the host's ability to elicit a timely and effectual immunologic response. The immune response in animals infected with the lethal strain of P. yoelii may in fact be normal (as shown by the early induction of gene ontology functional groups associated with immune response and cytokine and chemokine signaling pathways), but these animals die prior to the acquisition of antiparasite antibodies.
In conclusion, expression profiling is an extraordinarily useful tool to explore the molecular features of the host's response to infection and can provide insights into transcriptional regulatory mechanisms that influence both the pathogenesis of disease and the host's recovery from infection. While immune responses in human P. falciparum and P. vivax malaria may share many similar features of the global gene expression program observed in murine malaria, important differences in expression profiles in humans infected clinically or experimentally with malaria will depend heavily on the type of tissue (peripheral blood, bone marrow, spleen, or brain) and the stage of infection (early asymptomatic versus clinical malaria) that is studied.
Supplemental material for this article may be found at http://iai.asm.org/ ![]()
|
|
|---|
, with lethality of Plasmodium yoelii 17XL, a rodent model of cerebral malaria. Am. J. Trop. Med. Hyg. 59:852-858.[Abstract]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»