Previous Article | Next Article ![]()
Infection and Immunity, May 2008, p. 1812-1824, Vol. 76, No. 5
0019-9567/08/$08.00+0 doi:10.1128/IAI.01650-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Nicholas H. Hunt,* and
Helen J. Ball
Molecular Immunopathology Unit, Bosch Institute, University of Sydney, Sydney, Australia
Received 11 December 2007/ Returned for modification 2 February 2008/ Accepted 18 February 2008
|
|
|---|
|
|
|---|
A number of significant changes in the murine brain that characterize CM have been described, including breakdown of the blood-brain barrier (62), redistribution and activation of glia (31, 32, 34), apoptosis of endothelial cells and astrocytes (46, 47), metabolic changes (44, 48), and distinct patterns of chemokine gene expression (38).
A recent approach to identifying novel changes within the brain during CM is to profile gene expression changes using microarray technology (9, 30). Microarray studies can generate a plethora of data that pinpoint involvement of unknown or poorly characterized genes, stimulating investigation of their biological functions. Some studies have used microarrays to examine global gene expression patterns during murine malaria, defining gene expression profiles that are associated with susceptibility or resistance to murine CM (26, 27), identifying genes in the brain that discriminate between different outcomes of Plasmodium infection (5), or describing trends in gene expression within the brain (55). Studies on the spleen have been performed to find genes that differentiate between cerebral and noncerebral (55) and between lethal and nonlethal (52) forms of murine malaria. Microarrays also have been used to examine Plasmodium infection in humans (14, 41) and monkeys (68). In these studies, the majority of genes identified have well defined roles in such functions as immunity, erythropoiesis, and metabolism.
In the current study we used microarrays to examine global gene expression patterns in the brains of Plasmodium berghei-infected mice. In contrast to previous studies, which examined the response of genetically divergent murine hosts to the parasite that causes FMCM, P. berghei ANKA (PbA), we investigated the response of genetically identical murine hosts to different, but related, parasite strains, PbA and P. berghei K173 (PbK), the latter being one that results in NCM. Furthermore, we used two different microarray technologies, oligonucleotide and Affymetrix arrays, and compared the outcomes. Using either technique, we found similar patterns, though different degrees, of gene expression in groups irrespective of progression to CM. These included genes belonging to interferon (IFN) signaling pathways and major histocompatibility complex (MHC) processing and presentation, as well as those associated with apoptotic, immunomodulatory, and antimicrobial processes. However, we found that very strong induction of these genes was associated with fatal outcome due to FMCM, as confirmed by real-time PCR. These results contribute to our understanding of the processes occurring in the brain during CM, and the roles of the genes identified deserve further attention in FMCM and other neurological diseases.
|
|
|---|
Tissue extraction and histopathology. Mice were anesthetized by isoflurane inhalation and euthanized at each stated time point. Mice were perfused with 10 ml cold sterile phosphate-buffered saline to clear the vessels of blood, and organs were removed. Formalin-fixed, paraffin-embedded sections (7 µm) were used to examine histopathological features, as described previously (47).
RNA extraction. Brain tissue from animals was collected in 1 ml TRIzol reagent (Invitrogen, Victoria, Australia) and immediately homogenized using Zirconia beads (Biospec Products, OK) and a Fast Prep homogenizer (Q-Biogene, CA). Total RNA was isolated according to the manufacturer's instructions (Invitrogen). The integrity of the RNA samples was tested using Bioanalyzer Nano chips (Agilent Technologies) according to the manufacturer's instructions. Samples were DNase treated using a DNA-free kit (Ambion, TX). The same RNA samples were hybridized to the oligonucleotide and Affymetrix arrays as were used in the real-time PCR analyses.
Experimental design of microarray studies. All experiments were designed to be compliant with Minimum Information About a Microarray Experiment (MIAME) standards. For oligonucleotide array experiments, control samples were pooled and used as a common reference for the test samples. Experimental groups of five mice were either inoculated with PbA and sacrificed for sample preparation at day 6 p.i. [PbA(6)] or inoculated with PbK and sacrificed for sample preparation at day 6 p.i. [PbK(6)] or day 14 p.i. [PbK(14)]. Five individual test samples were used per group and directly compared to the pooled (n = 5) control reference sample. For Affymetrix array experiments, biological replicates (n = 5) from all groups [control, PbA(6), PbK(6), and PbK(14)] were pooled and hybridized to an individual slide. These were compared with the pooled control sample in downstream computational analyses.
Oligonucleotide arrays: cDNA labeling, hybridization, scanning, and normalization.
For oligonucleotide arrays, indirect labeling was used to incorporate cyanine (Cy3 and Cy5) dyes into cDNA samples. Labeled probes were then applied to 22.5K Compugen mouse oligonucleotide arrays (
22,500 genes) (Ramaciotti Centre, University of New South Wales). Slides were washed and then scanned, and spot intensities were determined using an Axon scanner (GenePix 3.0; Molecular Devices, CA). Basic normalization procedures for oligonucleotide arrays were carried out in GenePix 3.0, including background subtraction to remove data generated from noise and removal of poor-quality spots due to hybridization artifacts or morphology. Data were then imported into Limma (R, Bioconductor project; http://www.bioconductor.org/) (12) for loess normalization. For full experimental details, including raw data sets, refer to the Gene Expression Omnibus website (http://www.ncbi.nlm.nih.gov/geo/index.cgi), accession number GSE9808.
Prior to normalization, the majority of the arrays revealed systemic biases (i.e., more spots with negative M values), as well as low-intensity spots (i.e., spots with low A values), which can affect ratios. Some arrays also showed fanning effects (i.e., spots with low intensity and high background), which can lead to very large but erroneous ratios. MA plots were used to visualize these spatial and intensity effects (see Fig. S1a in the supplemental material), while box plots were used to view the spread of M values (see Fig. S1c in the supplemental material). Many arrays showed variations across the M values, particularly across each print-tip group. These systemic biases suggested that local and not global methods of normalization are more appropriate in this situation. From the raw data, local normalization methods were employed to correct for technical effects that can contribute to differential gene expression. This was performed using loess normalization in the Limma program (57). Normalized data showed a better spread across A values, without extreme variation across M values, as viewed in MA plots (see Fig. S1b in the supplemental material). Similarly, in box plots, normalized M values were found to cluster more tightly around M = 0, with less spread than for raw data (see Fig. S1d in the supplemental material).
Affymetrix arrays: cDNA labeling, hybridization, scanning, and normalization.
For Affymetrix arrays, cRNA was labeled using the one-cycle target labeling protocol as described in the Affymetrix GeneChip Expression Analysis Manual (Affymetrix). Labeled cRNA samples were then hybridized to Affymetrix GeneChip Mouse Genome 430 2.0 arrays (
43,000 genes) (Ramaciotti Centre, University of New South Wales) before being scanned using the GeneArray Scanner (Affymetrix). Affymetrix arrays were initially assessed and single-array analyses carried out in the GeneChip Operating Software program according to the Affymetrix GeneChip Expression Analysis Data Analysis Fundamentals Manual (Affymetrix) on data generated from the image analyses. For full experimental details, including raw data sets, refer to the Gene Expression Omnibus website (http://www.ncbi.nlm.nih.gov/geo/index.cgi), accession number GSE9808.
Microarray data analysis. Following basic normalization procedures carried out using Limma and GeneChip Operating Software, data from both oligonucleotide arrays and Affymetrix arrays were imported into Genespring (Agilent Technologies) for further analyses. Within the Genespring program, genes are assigned individual P values, based on a one-sample Student t test and assuming a baseline of 0 on a log scale. For oligonucleotide arrays, genes were selected as differentially expressed if they had a normalized value greater than 2 (upregulated) or less than 0.5 (downregulated) on a log scale, with at least three slides in each group satisfying these criteria. Genes were then ranked according to their P values to generate lists of upregulated and downregulated genes. However, one-way analyses of variance were found to be too stringent when applied to these lists, as no genes were found to be differentially expressed, and thus were not used. For Affymetrix arrays, following normalization using robust multiarray analysis, arrays from P. berghei-infected mice were compared to the uninfected control array. Genes were selected as differentially expressed if they had a normalized value greater than 2 (upregulated) or less than 0.5 (downregulated) on a log scale. Genes were then ranked according to their fold change values to generate lists of upregulated and downregulated genes.
Functional annotation of genes. Functional annotation of genes was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (http://david.abcc.ncifcrf.gov/home.jsp) (6). Genes were annotated based on the three structured vocabularies (ontologies), describing gene products in terms of their "biological processes," "molecular functions," and "cellular compartments," which are considered independent of each other. DAVID has the advantage of clustering genes with similar gene ontology (GO) terms together, thus removing the high level of redundancy that exists in the ontologies and making interpretation easier.
Real-time PCR. Real-time PCR was performed as described previously (2), except that PCRs were performed with the Corbett Rotor-Gene (Corbett Research, New South Wales, Australia) and using Platinum Sybr green qPCR SuperMix UDG (Invitrogen). mRNA levels were normalized to a housekeeping gene, that for hypoxanthine guanine phosphoribosyltransferase (HPRT), and expressed relative to the mean for uninfected control samples using the delta threshold cycle method. The specificity of primers was checked using melting curve analyses, and the primer sets used (Table 1) had similar amplification efficiencies.
|
View this table: [in a new window] |
TABLE 1. Primer sequences used
|
Statistical analyses. Statistical analyses for gene expression data were performed using the Kruskal-Wallis test. Where statistical differences were observed, probabilities were then calculated using the Mann-Whitney test on selected comparisons. These analyses were carried out using GraphPad Prism version 4.00 for Windows (GraphPad Software, CA).
Microarray accession number. The Gene Expression Omnibus accession number for the microarray data reported in this paper, which also includes details of the MIAME-compliant experiments, is GSE9808.
|
|
|---|
![]() View larger version (43K): [in a new window] |
FIG. 1. Line graphs of array gene expression data as visualized in Genespring, for oligonucleotide arrays, showing upregulated genes (fold change of >2.0) in PbA(6) (A), PbK(6) (B), and PbK(14) (C) samples and downregulated genes (fold change of <0.5) in PbA(6) (D), PbK(6) (E), and PbK(14) (F) samples, and also for Affymetrix arrays, showing color-coded gene expression in PbA(6) (G), PbK(6) (H), and PbK(14) (I) samples (red, expression of >1.0 [upregulated]; green, expression of <1.0 [downregulated]).
|
-induced GTPase (IGTP), as well as poorly annotated ISGs such as IFN-induced protein with tetratricopeptide repeats 2 (IFIT2) and IFN-induced transmembrane protein 3 (IFITM3), IFN-induced gene 30 (Ifi30), Ifi35, and Ifi44 (see Tables S1 and S3 in the supplemental material). Genes associated with cellular compartments as well as unannotated genes were found on lists of downregulated genes (see Tables S2 and S4 in the supplemental material). Genes are differentially expressed between the PbA(6), PbK(6), and PbK(14) groups. On closer inspection, many genes were common to each group [i.e., PbA(6), PbK(6), or PbK(14)], indicating that gene expression patterns may overlap. In order to identify genes specific to each group, different lists were directly compared using Venn diagrams in the Genespring program (Fig. 2; see Tables S5 and S6 in the supplemental material).
![]() View larger version (13K): [in a new window] |
FIG. 2. Venn diagrams of oligonucleotide arrays, showing the intersection ( ) between different groups, where group A includes genes upregulated only in PbA(6) (n = 88), group B includes genes upregulated only in PbK(6) (n = 17), group C includes genes upregulated only in PbK(14) (n = 50), group D includes genes upregulated in PbA(6) PbK(6) (n = 4), group E includes genes upregulated in PbA(6) PbK(14) (n = 16), group F includes genes upregulated in PbK(6) PbK(14) (n = 6), and group G includes genes upregulated in all groups (n = 24). Group H includes genes downregulated only in PbA(6) (n = 65), group I includes genes downregulated only in PbK(6) (n = 168), group J includes genes downregulated only in PbK(14) (n = 81), group K includes genes downregulated in PbA(6) PbK(6) (n = 22), group L includes genes downregulated in PbA(6) PbK(14) (n = 11), group M includes genes downregulated in PbK(6) PbK(14) (n = 10), and group N includes genes downregulated in all groups (n = 9).
|
). Four genes were upregulated in both PbA(6) and PbK(6) (group D), 16 genes were upregulated in both PbA(6) and PbK(14) (group E), 6 genes were upregulated in both PbK(6) and PbK(14) (group F), and 24 genes were commonly upregulated in all three groups (group G). In these intersecting groups (i.e., D, E, F, and G), genes were commonly associated with immune functions (e.g., β2-m, CCL12, CXCL9, CXCL11, cytotoxic-T-lymphocyte-associated protein 2
[CTLA-2
], granzymes A and B, guanylate nucleotide binding 1, histocompatibility 2, Ii, IFIT1, IFIT3, IGTP, IRF1, and metallothionein 1) (Fig. 2; see Table S5 in the supplemental material). Similarly, many of these genes were found on lists of upregulated genes on the Affymetrix arrays following P. berghei infection (see Table S6 in the supplemental material).
Functional clustering reveals distinct groups of differentially expressed genes. Upregulated genes were clustered using the DAVID program, based on level 5 GO terms. For oligonucleotide arrays, some clusters were nonspecifically associated with P. berghei infection, including clusters of genes for biological processes related to immune functions (e.g., antigen processing and response to pathogen) as well as molecular functions associated with signaling (e.g., nucleotide binding). There were a few clusters of genes that were specific to PbA(6), including those for biological processes (e.g., catabolism, negative regulation of biological and cellular processes, and signal transduction) and cellular compartments (e.g., nuclear lumen and non-membrane-bound organelle). While no clusters were specific to PbK(6), a cluster for a cellular compartment (e.g., hemoglobin complex) was found to be specific to PbK(14) (see Table S7 in the supplemental material).
Clustering of Affymetrix arrays revealed similar groups of genes, particularly genes for biological processes related to immune functions (e.g., antigen processing and response to pathogen) and molecular functions associated with signaling (e.g., nucleotide binding and receptor binding). Again there were a few clusters of genes that were specific to PbA(6), including those for biological processes (e.g., cytokine biosynthesis, regulation of enzyme activity, and endocytosis) and molecular functions (e.g., ligase activity) (see Table S8 in the supplemental material).
Clustering groups were diverse in downregulated genes, and unannotated genes were prevalent (see Tables S7 and S8 in the supplemental material). There was no discernible pattern and little overlap in the annotations of downregulated genes during P. berghei infection, due in part to the green bias on oligonucleotide arrays. Moreover, since many of the genes found in these lists were not annotated, these were not explored further here.
Annotations and classes of differentially expressed genes. To gain a clearer overview of genes differentially expressed during P. berghei infection, gene lists from oligonucleotide and Affymetrix arrays were combined. Importantly, both oligonucleotide and Affymetrix array determinations were performed using the same RNA samples obtained from the same biological experiment. Genes were annotated using level 3 GO terms in DAVID and the percentage of functional classes determined, together with their statistical significance (see Table S9 in the supplemental material).
There was a large overlap in the annotations for genes that were upregulated during P. berghei infection, many of which were overrepresented. These were associated with GO terms from biological processes such as "immune response" [PbA(6), 110 genes, P = 1.98 x 10–49; PbK(6), 56 genes, P = 1.19 x 10–38; and PbK(14), 74 genes, P = 1.56 x 10–39], molecular functions such as "cytokine activity" [PbA(6), 16 genes, P = 2.63 x 10–4; PbK(6), 5 genes, P = 6.84 x 10–2; and PbK(14), 10 genes, P = 4.13 x 10–3], and cellular compartments such as "immunological synapse" [PbA(6), 13 genes, P = 3.45 x 10–8; PbK(6), 13 genes, P = 1.80 x 10–14; and PbK(14), 14 genes, P = 3.68 x 10–12] (see Table S9 in the supplemental material). Overall, there were more genes upregulated by PbA infection than by PbK infection. Genes belonging to selected annotation categories are shown in Table 2.
|
View this table: [in a new window] |
TABLE 2. Selected lists of functionally annotated genes that are upregulated during PbA infection
|
|
View this table: [in a new window] |
TABLE 3. Comparisons between selected gene expression data from oligonucleotide and Affymetrix arrays
|
Genes of interest that were upregulated during PbA infection from both oligonucleotide and Affymetrix arrays were then selected for real-time PCR validation. These included genes with known immune functions such as Ii, β2-m, and granzymes A and B, as well as some genes with interesting annotations, such as OAS-1b, ISG15, IGTP, cathepsin W (CTSW), and CTLA-2
. Many genes were upregulated during PbA infection but not PbK infection. Figure 3 shows representative graphs of genes upregulated at the mRNA level.
![]() View larger version (29K): [in a new window] |
FIG. 3. mRNA expression of OAS-1b (A), ISG15 (B), IGTP (C), Ii (D), CTSW (E), and CTLA-2 (F) in whole-brain homogenates of uninfected control mice and PbA(6)-, PbK(6)-, and PbK(14)-infected mice. Expression, measured by real-time PCR as described in Materials and Methods, is relative to the mean expression value in uninfected control samples. Values are means ± SEMs (n = 6). Samples were initially determined to be significantly different using the Kruskal-Wallis test (P < 0.001), and then the Mann-Whitney test was used: ++, P < 0.01 (values compared with uninfected controls); ##, P < 0.01 [values compared with PbA(6)]; **, P < 0.01 [values compared with PbK(6)]; *, P < 0.05 [values compared with PbK(6)].
|
There also were some genes that were predominantly upregulated at the end stage of the disease [PbA(6) and PbK(14)]. These included Ii mRNA (Fig. 3D), which was strongly induced in PbA(6) (54.1- ± 3.9-fold) and PbK(14) (58.4- ± 16.2-fold) compared with PbK(6) (19.8- ± 3.2-fold). CTSW mRNA (Fig. 3E) was upregulated in the PbA(6) (31.5- ± 9.1-fold) and PbK(14) (29.2- ± 11.4-fold) compared with PbK(6) (11.5- ± 0.7-fold). CTLA-2
mRNA (Fig. 3F) had a similar expression pattern, being upregulated in PbA(6) (28.0- ± 4.8-fold) and PbK(14) (44.2- ± 17.8-fold) compared with PbK(6) (4.0- ± 0.2-fold).
Overall, the real-time PCR data confirmed the upregulation of genes found on both the oligonucleotide and Affymetrix arrays and showed concordance between these two different platforms. However, the real-time PCR results also revealed that differential expression was not always easily identified by arrays, particularly by oligonucleotide arrays, which are less sensitive than Affymetrix arrays. By comparing the results from the arrays with the real-time PCR results, it is evident that, during hybridization of the samples, a threshold is reached on the oligonucleotide arrays that does not allow highly upregulated or downregulated genes to be distinguished from genes that are also differentially expressed at lower levels. Thus, a number of genes found to be upregulated by P. berghei infection on the arrays also were found to be differentially expressed between PbA(6), PbK(6), and PbK(14) upon validation with real-time PCR (Fig. 3).
There also appeared to be a significant problem associated with the detection of genes that are of low abundance. For example, tumor necrosis factor is known to be highly upregulated in the central nervous system during PbA infection (35) but was not identified as an upregulated gene on the arrays, despite being represented. With real-time PCR, tumor necrosis factor mRNA (data not shown) was found to be strongly upregulated in PbA(6) compared with PbK(6) and PbK(14). These findings suggest that a number of differentially expressed genes are not detected above background levels even though they may play critical roles during PbA infection.
Cellular expression of differentially expressed genes. The cellular sources of a number of these genes were then determined using laser capture microdissection. Vessels were separated from parenchymal cells of the brain (neurons and glia) in PbA-infected mice and real-time PCR performed to determine relative gene expression.
The majority of genes, including Ii, β2-m, OAS-1b, ISG15, and IGTP, were found to be expressed predominantly by vessel-associated cells compared with parenchymal cells, which expressed mRNAs for these genes at much lower levels. However, CTLA-2
and Fas were very highly expressed by vessel-associated cells but were not detectable in parenchymal cells (Table 4).
|
View this table: [in a new window] |
TABLE 4. Relative mRNA expression in microdissected cells during FMCMa
|
![]() View larger version (10K): [in a new window] |
FIG. 4. mRNA expression of OAS-1b (A), ISG15 (B), and IGTP (C) in whole-brain homogenates of uninfected control and PbA-infected mice over the course of the disease. Expression, measured by real-time PCR as described in Materials and Methods, is relative to the mean expression value in uninfected control samples. Values are means ± SEMs (n = 4 to 8). Samples were initially determined to be significantly different using the Kruskal-Wallis test (P < 0.001), and then the Mann-Whitney test was used: +++, P < 0.001 (values compared to uninfected controls, day 0); ++, P < 0.01 (values compared to uninfected controls, day 0); +, P < 0.05 (values compared to uninfected controls, day 0); **, P < 0.01 (values compared to day 4 p.i.); *, P < 0.05 (values compared to day 4 p.i.).
|
|
|
|---|
during the pathogenesis of FMCM (1, 13, 51). Some processes that were evident from gene expression patterns on the arrays included antigen processing and presentation, apoptosis, and antimicrobial pathways. While some of these genes previously have been established to play a role in the pathogenesis of FMCM, for example, intercellular adhesion molecule-1 (28) and granzyme B (46), there also were a number of upregulated genes, previously unreported, that represent interesting candidates for further investigation.
(i) IFN signaling pathways.
Microarray analysis demonstrated that genes for essential components of the IFN signaling pathway, including members of the signal transducers and activators of transcription family and the IRF family, were induced during P. berghei infection. IRFs are secondary transcription factors that are induced early following IFN-
stimulation and drive the transcription of IFN-
-responsive genes, including other ISGs (54). IRF1, a transcriptional activator, was upregulated on the arrays, whereas IRF2, a transcriptional repressor, was not. IRF1 is required for the development of a Th1-type immune response, and its absence leads to the induction of a Th2-type response (25, 59). Negative regulators of the IFN signaling pathway, for example, the family of suppressor of cytokine signaling proteins, also were upregulated (see Tables S1 and S3 in the supplemental material).
(ii) MHC processing and presenting pathways. Components of the MHC I and MHC II pathways were found by microarray analysis to be expressed during P. berghei infection (see Tables S1 and S3 in the supplemental material). This was not surprising given the upregulation of genes associated with IFN signaling pathways, particularly IRF1, and the importance of CD4+ and CD8+ T cells in malaria immunity and immunopathology (16, 46, 50, 58, 67). The expression of genes involved in antigen-processing and -presenting pathways known to be involved in the pathogenesis of murine malaria, for example, MHC peptides and costimulatory molecules (24), lends confidence to the validity of the upregulation of novel genes found by microarray analysis without further confirmation.
(iii) Apoptosis and cell proliferation. A number of genes associated with apoptotic pathways were expressed on the arrays following P. berghei infection. These included caspases, granzymes, cathepsins, calpains, protease inhibitors, OASs, and members of the p200 family (see Tables S1 and S3 in the supplemental material). Significantly, these were much more strongly induced in PbA infection than in PbK infection. Apoptosis occurs in the brain during human and murine CM in endothelial cells (45, 46, 65), astrocytes (47), and neurons (23, 53, 65). Both the death receptor-mediated pathway and the stress-induced pathway of apoptosis have been implicated in FMCM (21, 42, 45, 47). Both pathways converge on caspase 3, which also is involved in the pathogenesis of FMCM (23, 36, 45-47).
Other proteases, including granzymes, cathepsins, and calpains, also play a role in apoptosis, primarily by activating caspases. Some of these genes also have been implicated in neuronal death in neurological diseases (11, 17, 60), including CM (33, 56). A number of protease inhibitors were expressed on the arrays, including CTLA-2
(22), and may represent a host mechanism to limit the damage mediated by proteases, particularly since these proteolytic processes can lead to development of intracerebral hemorrhaging (63), which is an invariable finding in human and murine CM.
(iv) Immunomodulatory and antimicrobial genes. Genes involved in immune responses, such as Fc receptors and components of the complement system, were upregulated during P. berghei infection (see Tables S1 and S3 in the supplemental material), reflecting the important role of the immune system in this disease (19). Genes involved in adaptive immunity and other immunodulatory functions, including cytokines, chemokines, leukocyte surface antigens, and adhesion molecules, also were induced. There were a number of upregulated ISGs with interesting but poorly defined functions, including the family of IFIT genes and the family of IFITM genes. It remains to be seen whether these genes play other roles in host defense, particularly against parasites.
A number of interesting immunomodulatory genes expressed on the arrays have not yet been assigned a role in CM. For example, ISG15 is released by lymphocytes and monocytes (4, 20) and is inducible by IFN (20) but can itself induce the secretion of IFN-
from T cells and stimulate the activity of indoleamine 2,3-dioxygenase 1 (49), which has complex roles in malaria infection (18), as well as lymphocyte and NK proliferation, and can induce cytolytic activity (3). Other ISGs identified by arrays belong to the p47 and p65 GTPase (i.e., IGTP) families (29, 61). These GTPases have antimicrobial functions and are preferentially induced by IFN-
, a cytokine that is essential for the pathogenesis of FMCM through ill-understood mechanisms. While the roles of these genes are not well described, specific gene knockout mice have differences in their susceptibilities to infectious diseases, including those of the central nervous system (29, 61).
Real-time PCR Validation.
Selected genes were validated using real-time PCR (e.g., OAS-1b [Fig. 3A], ISG15 [Fig. 3B], IGTP [Fig. 3C], Ii [Fig. 3D], CTSW [Fig. 3E], and CTLA-2
[Fig. 3F]), confirming the pattern of gene expression seen in the arrays (i.e., upregulated in P. berghei infection). However, real-time PCR also revealed that the expression of some of these mRNAs was greater in PbA(6) and PbK(14) than in PbK(6).
Overall these results suggest that there are similar processes occurring in the brain at the end stages of P. berghei infection, independent of the occurrence of cerebral symptoms. These include those associated with MHC antigen-processing and -presenting pathways, although there was a striking difference in the mRNA upregulation of costimulatory molecules on the arrays. While PbA infection was associated with increased mRNAs for a number of costimulatory molecules (e.g., CD40 and CD86; see Table S3 in the supplemental material), this was not seen during PbK infection. In the absence of costimulation, T cells often undergo anergy and apoptosis instead of proliferation and differentiation. Since costimulation is an essential component of the antigen presentation pathway, this suggests that there is an absence of T-cell activation during PbK infection.
The location of gene expression in the brain clearly is important. The expression of Ii and β2-m mRNAs was greater in vessel-associated cells than in parenchymal cells (Table 4). This suggests that endothelial cells are most likely to express malaria antigens associated with either MHC I or MHC II. Indeed, brain endothelial cells are capable of expressing MHC I and MHC II, and their expression levels are associated with the development of FMCM (39). Microglia also can express MHC I, but not costimulatory molecules, during FMCM (43). Similarly, mRNAs of genes associated with apoptotic processes, such as Fas, OAS-1b, and CTLA-2
, were localized predominantly to vessel-associated cells compared with parenchymal cells (Table 4), suggesting that important interactions occur at the blood-brain barrier during FMCM. Furthermore, ISG15 and IGTP mRNAs also are predominantly localized to cells associated with cerebral vessels, and given the cytokine-like activities of ISG15 and its potential nonimmune roles in the brain, as well as the antimicrobial actions of IGTP, these are interesting candidates to examine in malaria infection.
These findings emphasize the key role played by cells of the blood-brain barrier, particularly endothelial cells, during FMCM. These cells are located at the key biological interface between the blood compartment, including the circulating intraerythrocytic malaria parasite, and the brain parenchyma. Changes to the blood-brain barrier during malaria during FMCM are well described (18, 19) and are thought to be both a cause and an effect of inflammatory processes. The prominent localization of gene induction to cells of the blood-brain barrier supports the idea that key interactions between the parasite and the host occur at the endothelium during FMCM, and as a consequence, pathological events are concentrated primarily within cerebral microvessels, with lesser and later involvement of cells from the brain parenchyma.
Concluding remarks. Overall, these data are supportive of an inflammatory reaction occurring in the brain during FMCM, which is typified by high circulating levels of inflammatory mediators, an activated endothelium, and adherent leukocytes (19). There was a high degree of concordance between results obtained from the oligonucleotide and Affymetrix arrays, which revealed similar patterns of gene expression.
Some gene expression changes typical of inflammatory responses also occurred at the end stage of PbK infection, in the absence of neurological signs. These changes likely result from high circulating levels of PRBCs and cytokines (37) but, overall, were small compared with changes due to PbA infection and presumably are insufficient for immunopathological processes. A number of genes, such as intercellular adhesion molecule-1 and CXCL9, were upregulated at the end stage of PbA and PbK infection and are largely indicative of an activated endothelium. Moreover, indoleamine 2,3-dioxygenase-1 is very strongly induced in endothelial cells during the late stages of PbK (15), perhaps as a tissue-protective response (18), suggesting that important responses also can occur in the brain during NCM.
These results support the view that leukocytes play a crucial role in the pathogenesis of FMCM. Leukocytes do not adhere to the microvasculature during noncerebral models of malaria, including late stages of PbK infection, and during PbA infection in genetically resistant mice, for example, those deficient in IFN-
or its receptor (1, 51), LT-
(10), and CXCR3 (38), despite continued expression of adhesion molecules by the endothelium. This suggests that leukocytes are unable to adhere to the endothelium during PbK infection, perhaps because of a lack of integrin activation or an absence of costimulation by antigen-presenting cells, thus leading to leukocyte anergy.
Overall, these data highlight the significant gene expression changes in the brain during PbA infection (FMCM). Furthermore, certain critical mediators that were upregulated in PbA were absent during PbK infection. In FMCM, several induced pathways were evident from the gene expression profiles on the arrays involving signaling and transcription, apoptosis, immunomodulation, and antimicrobial activities. Some of these have been identified previously in FMCM using arrays (5, 55), but without comparisons having been made with late stages of NCM. The results from this study contribute to our understanding of the complex pathogenic processes that occur during malaria.
This work was supported by grants from the National Health and Medical Research Council of Australia and the Sir Zelman Cowen Universities Fund to N.H.H. J.M. was supported by an Australian Postgraduate Award. H.J.B. was a Rolf Edgar Lake Research Fellow of the Faculty of Medicine, University of Sydney.
Published ahead of print on 25 February 2008. ![]()
Supplemental material for this article may be found at http://iai.asm.org/. ![]()
Present address: Department of Medicine, McGill University, Montreal, Quebec, Canada. ![]()
|
|
|---|
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»