Previous Article | Next Article ![]()
Infection and Immunity, October 2008, p. 4518-4529, Vol. 76, No. 10
0019-9567/08/$08.00+0 doi:10.1128/IAI.00525-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Laboratory of Malaria and Vector Research,1 Infectious Disease Pathogenesis Section, Comparative Medicine Branch, Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health,4 Division of Emerging and Transfusion Transmitted Diseases, Center for Biologics Evaluation and Research, Food and Drug Administration, Rockville,5 Emerging Infectious Diseases Program, Uniformed Services University of the Health Sciences,2 National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, Maryland3
Received 29 April 2008/ Returned for modification 5 June 2008/ Accepted 7 July 2008
|
|
|---|
1 proteins by Western blot analysis than mice unaffected by ECM, suggesting the possible utility of these molecules as prognostic biomarkers of CM in humans. We suggest that the higher expression of hemoglobin
1 in the brain may be associated with ECM and could be a source of excess heme, a molecule that is considered to trigger the pathogenesis of CM. Our studies greatly enhance the repertoire of host molecules for use as diagnostics and novel therapeutics in CM. |
|
|---|
Susceptibility to CM is believed to result from a combination of parasite and host factors. The ability of the parasite to sequester itself within the deep vascular beds of brain tissue is thought to be the dominant parasite characteristic required for CM (34). In previous studies, associations between a number of host molecules and susceptibility to CM have been established (65). The majority of these genes, including the tumor necrosis factor alpha (TNF-
) and gamma interferon (IFN-
) genes, encode proteins of the immune system. The development of a prognostic test that could be used to predict the outcome of P. falciparum infection in a susceptible host will require a complete knowledge of the expression pattern of the pathogenesis factors that are specifically triggered during a clinical episode of CM. To obtain a more accurate and complete picture of the pathogenesis of CM, we have identified the repertoire of genes whose levels are altered by measuring their expression in whole-brain tissue in real time and determined the biological associations at the protein level of a few selected molecules in the Plasmodium berghei ANKA murine model of experimental cerebral malaria (ECM).
Like P. falciparum, P. berghei ANKA, a strain causing murine malaria, possesses the ability to sequester itself within the microvasculature; it is the parasite of choice for in vivo studies of CM in mice (17). Infection with P. berghei ANKA parasites results in death for 100% of mice. However, the cause of death differs greatly based on the genetic background of the mouse, and different strains of mice are broadly categorized as susceptible or resistant to CM (17). In susceptible strains (e.g., C57BL/6 and CBA/J), approximately 60 to 100% of mice develop CM 6 to 10 days postinfection, and mice die with relatively low levels of parasitemia (below 10%). The features of ECM in these mice resemble those of P. falciparum clinical disease in humans: hemi- or paraplegia, deviation of the head, ataxia, distressed breathing, and convulsions. In contrast, resistant mice (e.g., BALB/c mice) do not develop any clinical or pathological features resembling those of CM (17) but instead succumb to anemia and hyperparasitemia approximately 15 to 21 days after infection with P. berghei ANKA. Susceptible mice that do not develop ECM follow a clinical course of disease similar to that of resistant mice and eventually die of hyperparasitemia and severe anemia. Furthermore, previous studies have shown that loss of function of CD8 or perforin molecules converts susceptible C57BL/6 mice to resistant mice whose resistance to ECM is comparable to that observed for the BALB/c strain (38, 47, 69). These studies suggest a role for the cell-cell interaction pathways dependent on these molecules in augmenting CM.
The goal of this study was to identify novel host molecules in the brain that are affected in the course of the pathological manifestation of CM. In particular, we wanted to identify novel host molecules that could serve as prognostic biomarkers of CM. Using high-density microarray chips, we compared the whole-brain tissue expression profiles of susceptible C57BL/6 mice with ECM (moribund mice) to those of (i) infected C57BL/6 mice without ECM (nonmoribund mice) and (ii) infected BALB/c mice (resistant mice) during the effector phase of ECM, the phase most pertinent to the management and therapeutic treatment of CM. Furthermore, to determine the mechanisms of CD8- and perforin-mediated pathogenesis of ECM, we also compared the expression profiles of infected C57BL/6 wild-type (WT) mice with ECM to those of (i) infected CD8 gene knockout (CD8-KO) mice and (ii) infected perforin gene knockout (PFP-KO) mice on the C57BL/6 background, both of which are resistant to ECM. Based on differential expression in brain tissue samples, we identified a wide range of host genes that were strongly associated with the manifestation of ECM. Bioinformatic analyses were performed using several computational tools for sequence analysis and contextual functional inference, and genes with significantly altered expression were assigned to several distinct biological functional categories. Several of these molecules belong to families of known immunological relevance that could either serve as novel cytoadherence molecules or contribute directly to the pathogenesis of CM.
|
|
|---|
Preparation of whole-brain RNA. Brain tissue samples were collected simultaneously from each group and stored at –80°C until use. For RNA preparation, tissue samples were pulse-homogenized in Tri-Reagent (Molecular Research Center, Cincinnati, OH), and RNA was isolated following two chloroform extractions, isopropanol precipitation, and resuspension in diethyl pyrocarbonate-treated water. The quantity of RNA was determined by optical densitometry, and its quality was evaluated by agarose gel electrophoresis.
Microarray and gene expression analysis. The global gene expression profiles of moribund C57BL/6 WT mice were compared with those of both nonmoribund (C57BL/6 WT) and resistant (BALB/c, C57BL/6 CD8-KO, and C57BL/6 PFP-KO) mice to determine the array of specific host genes induced and repressed during the expression of ECM. Microarray expression profiles were determined from RNA samples isolated from five to six individual mice per group. cDNA was synthesized from 30 µg of RNA extracted from host brain tissue and then labeled with Cy3 or Cy5 as described previously (40). Briefly, cDNA was labeled with 50 nmol of dUTP-Cy3 or dUTP-Cy5 during a reverse transcription reaction at 42°C for 90 min in a sample mixture containing 300 mM dithiothreitol; 15 mM each dATP, dCTP, dGTP, and dTTP; a mixture of poly(T) and random hexamer primers; and 300 U Superscript II. After incorporation of the label, RNA was degraded by the addition of NaOH, and labeled cDNA was purified and concentrated by ultrafiltration through a Vivaspin 500 column (Sartorius, Goettingen, Germany). Labeled cDNA was then hybridized to a murine oligonucleotide chip containing 16,600 oligonucleotide probes (Qiagen, Valencia, CA), and the slide was scanned by a GenePix microarray scanner. Microarray data were analyzed with GenePix Pro software (version 6.0; Axon Instruments, Inc., Union City, CA) and filtered using the NIAID microarray database tools (http://madb-niaid.cit.nih.gov), and extracted spots were normalized to the precalculated 50th percentile (median).
Western blot analysis.
Protein samples from brain tissue were prepared as a 10% brain homogenate, and tissue-specific p21, metallothionein, and hemoglobin
1 proteins were detected using antibodies specific for recombinant human p21WAF1 protein (Biosource, Camarillo, CA), metallothionein (Abcam Inc.), or hemoglobin
1 (Novus Biologicals Inc.) and a commercially obtained chemiluminescence-linked Western blotting kit (Western Light Tropix, Bedford, MA). Protein bands were visualized following incubation with enhanced chemiluminescence detection reagents, and the integrated optical densities (IOD) for each lane were measured using MetaMorph software, version 6.1.
Gene annotation, ontology, and identification of biological pathways. Protein sequences were obtained from the nonredundant database of protein sequences (National Center for Biotechnology Information, NIH, Bethesda, MD). Profile searches were conducted using the PSI-BLAST program (3) with either a single sequence or an alignment used as the query, and with a default profile inclusion expectation (E) value threshold of 0.01 (unless specified otherwise), and the searches were iterated until convergence. The multiple alignments used were constructed using the T_Coffee program (39), followed by manual correction based on the PSI-BLAST results. For all searches involving membrane-spanning domains or low-complexity sequences, we used a statistical correction for compositional bias to reduce false-positive results due to the general hydrophobicity of these proteins. A library of a large set of alignments of conserved protein domains, including those from the PFAM database (http://www.sanger.ac.uk/Software/Pfam/index.shtml) as well as an additional set of unpublished conserved domains, was used for domain searches with the hidden Markov models (HMMER package) or with PSI-BLAST (position-specific score matrices). Transmembrane regions were predicted in individual proteins using the TMHMM program, version 2.0, with default parameters. Signal peptides were predicted using the SIGNALP program (www.cbs.dtu.dk/services/SignalP-2.0/). More than 300 profiles of known domains were used to search the mouse proteins by the method described above. The domains collected from these searches were assembled using the scripts from the in-house TASS package. Functions were assigned using the contextual information from the domain architectures and published experimental results. Contextual analysis was performed using the protein-protein and genetic interaction data stored in the BioGrid and VisANT databases. All large-scale sequence analysis procedures were carried out using the TASS package (V. Anantharaman, S. Balaji, and L. Aravind, unpublished data).
|
|
|---|
Alterations in gene expression profiles related to susceptibility or resistance to ECM.
We compared the global expression profiles of host genes in the brains of individual mice with ECM (C57BL/6 WT moribund mice) to those for mice without ECM (C57BL/6 WT nonmoribund, C57BL/6 CD8-KO resistant, C57BL/6 PFP-KO resistant, and BALB/c resistant mice), in order to identify novel host molecules/biomarkers associated with the pathogenesis of ECM, by performing microarray hybridizations with individual RNAs from brain samples isolated from five to six mice per group. For the comparison of moribund C57BL/6 mice with nonmoribund C57BL/6 mice and resistant BALB/c mice, altered expression of an individual gene was defined as significant only if two criteria were met: (i) moribund mice exhibited an average
2-fold increase (upregulation) or decrease (downregulation) in gene expression relative to expression for both nonmoribund and resistant mice (P < 0.05 by a two-tailed Student t test) (criterion 1) and (ii) moribund mice exhibited a
2-fold increase (upregulation) or decrease (downregulation) in three of the five arrays per group (criterion 2).
When moribund and nonmoribund C57BL/6 mice were compared, the input data were five arrays and 34,928 arrayed sequences corresponding mainly to protein-coding genes. A total of 34,144 sequences were excluded by criterion 1, and 172 additional genes were excluded by criterion 2. When moribund C57BL/6 mice were compared to resistant BALB/c mice, the input data were five arrays and 34,826 arrayed sequences. In this comparison, 33,543 sequences were excluded by criterion 1 and 169 additional sequences were excluded by criterion 2. A total of 612 genes were considered to show significantly altered expression when moribund C57BL/6 mice were compared to nonmoribund C57BL/6 mice, and 1,114 genes were considered to show significantly altered expression when moribund C57BL/6 mice were compared to resistant BALB/c mice. Of these 1,726 genes, 210 genes (66% upregulated, 34% downregulated) showed significantly altered expression in both the nonmoribund and resistant groups relative to expression in the moribund group. All microarray gene expression data are available in Table S1 in the supplemental material.
For microarray analysis comparing moribund WT to resistant CD8-KO and PFP-KO mice (six mice per group), data were analyzed using the significant analysis of microarrays (SAM) program, and a false-discovery rate (FDR) was chosen that would yield a set of genes that were both statistically significant and also of a size manageable for in-depth functional analysis. For comparison with CD8-KO mice, genes with an FDR of <2.2% and an average twofold increase (upregulation) or decrease (downregulation) in expression were included, yielding a data set of 166 genes (92% upregulated, 8% downregulated). For PFP-KO mice, genes with an FDR of <4.2% and an average twofold increase (upregulation) or decrease (downregulation) in expression were included, yielding a data set of 9 genes (78% upregulated, 22% downregulated).
In order to understand the significance of the altered gene expression patterns, we carried out a number of comparisons of these patterns between different sets of mice. We first compared the expression patterns seen in the nonmoribund C57BL/6 and resistant BALB/c mice to determine if the similar phenotypes of avoidance of ECM in mice with disparate major histocompatibility complex alleles have similar molecular profiles. Second, we compared the expression profiles of these nonmoribund or resistant mice with those of the resistance phenotypes associated with CD8-KO and PFP-KO mice in order to identify potential shared and different pathways leading to resistance in these mice. Comparison of the nonmoribund C57BL/6 and resistant BALB/c mice showed that the two groups have similar patterns of expression changes with respect to the moribund mice (Pearson correlation,
0.83). At least 60 genes whose expression patterns in the resistant CD8-KO group were altered relative to those in WT mice with ECM were also present in the list of genes with altered expression in the C57BL/6 nonmoribund and resistant BALB/c mice. This commonality in gene expression patterns suggests that the distinct pathways utilized by these mice to avoid ECM share certain common mechanisms. The resistance phenotype of PFP-KO mice was associated with significant alterations in only nine genes. At face value, this would suggest that the resistance phenotype in PFP-KO mice occurs via a relatively restricted downstream mechanism compared to the other resistance phenotypes. These results should be considered hardly surprising, since the resistance to ECM resulted from the loss of a single genetic locus in the mouse genome.
To glean the functional significance of the genes with altered expression, we performed a comprehensive computational sequence analysis of their protein products, identified all conserved domains present in those products, and predicted subcellular targeting signals such as signal peptides and transmembrane helices. We then combined this with a survey of relevant literature on all genes for which there were published references and with contextual information from high-throughput interaction mapping studies using the Biogrid and VisANT databases (22). As a result, we were able to confidently classify the majority of genes into various functional guilds (Table 1). Using this functional classification and the inferred interactions of the gene products, we pieced together several biological processes that are affected by ECM. We discuss some key examples below.
|
View this table: [in a new window] |
TABLE 1. Biomarkers of ECM categorized by function
|
- and β-chain genes was considerably downregulated in the brain during ECM (Table 2) in comparison to that for nonmoribund and resistant BALB/c mice. Because reticulocytes are the sole source of hemoglobin mRNA, this finding suggested a possible reduction in the circulating-reticulocyte count. This possibility is consistent with previous reports of impaired erythropoiesis in mice infected with lethal murine malaria strains (7, 54) and in clinical cases of severe malaria (27). To further investigate the possible implications of the reduction in hemoglobin chain expression in ECM, we measured the expression of hemoglobin
1 chain protein in the brain tissue. In contrast to transcriptional levels, we found that hemoglobin
1 protein was overexpressed 6.3-fold in moribund mice relative to expression in nonmoribund mice (Fig. 1A). These counterintuitive results suggest the distinct possibility of an imbalance in the number of reticulocytes with respect to mature erythrocytes (which contain hemoglobin protein but not RNA) during ECM. This imbalance, in turn, might be related to the elevated expression of heme oxygenase 1 (Hmox1) that we observed for moribund mice with ECM. Recent studies have shown that deletion of the Hmox1 gene or inhibition of the enzyme results in an elevated incidence of ECM (44), suggesting a role for free heme liberated from lysed erythrocytes in pathogenesis. Taken together, our observations might imply that the increase in the number of mature infected erythrocytes in the brain during ECM might result in an excessive release of heme, which, in turn, could trigger elevated Hmox1 levels as a defensive response. Given the identification of a human microsatellite polymorphism in the Hmox1 promoter that is associated with susceptibility to CM (61), such a mechanism might also be relevant in human CM. In this context, it is of interest that lipocalin 2, an iron-sequestering protein that is present at the blood-brain barrier and prevents bacterial entry into the brain (33), is overexpressed in ECM. It remains to be seen if lipocalin 2 might play a role in sequestering the iron released from ruptured erythrocytes. |
View this table: [in a new window] |
TABLE 2. Biological functions of a subset of ECM-related genes in our data set
|
![]() View larger version (31K): [in a new window] |
FIG. 1. Biomarkers of ECM as detected by Western blotting. The expression of murine biomarkers of ECM was measured in brain tissue samples from individual moribund (n = 5), nonmoribund (n = 5), CD8-KO (n = 3), and healthy (uninfected) (n = 3) C57BL/6 mice by enhanced chemiluminescence-based Western blot analysis. Expression levels were measured based on the intensities of protein bands for each lane by using MetaMorph software (version 6.1) and are represented as average IOD units. The actual IOD units are the values shown multiplied by 1,000. Details of the specific secondary antibodies and Western blot reagents are provided in Materials and Methods. (A) Expression of murine hemoglobin 1 protein; (B) expression of metallothionein protein; (C) expression of the p21 molecule.
|
Additionally, in agreement with a recent report published by Wiese et al. (67), in our study, mice displaying symptoms of ECM had increased transcriptional and translational expression of the genes encoding metallothioneins I and II (MT-I and -II) (Table 2 and Fig. 1B). MT-I and MT-II belong to a superfamily of low-molecular-weight, cysteine-rich metalloproteins and are reported to participate in a wide array of activities, including immunoregulation, cell survival, and neuroprotection (8, 9, 26, 45).
Alterations in apoptotic and counterapoptotic pathways. We observed that the expression of the gene encoding the kinase inhibitor of the Cip/Kip family, p21, is significantly upregulated during ECM (Table 2). p21 has been shown to be important in the regulation of apoptosis. Earlier studies have shown that expression of p21 increases when cells are damaged, and this increased expression might be a key element for cellular survival during oxidative stress (43). We examined the expression of p21 at the level of translation in order to further establish the role of p21 during ECM and to determine whether a correlation exists between the synthesis of mRNA and the synthesis of protein. Our microarray results show that the expression of p21 in moribund mice is 5.5-fold greater than its expression in nonmoribund mice at the level of transcription. Quantification of p21 protein in brain tissue by Western blotting demonstrated a 14.5-fold increase in expression at the level of translation (Fig. 1C). Thus, the expression pattern of the p21 gene reflects its elevated production, which potentially represents a counterapoptotic response in face of strong oxidative stress.
At least 10 other genes with potential roles in the regulation of apoptosis also showed altered expression in mice with ECM. In particular, we found that the levels of transcripts of two members of the TNF family, namely, TRAIL/CD253, were lower in mice with ECM. Furthermore, in the experiments with CD8-KO mice, we found that the expression of Fas (a receptor of ligands of the TNF family) was upregulated in moribund mice relative to that in resistant mutant mice. This finding suggests that multiple signals involving the TNF pathway may be involved in regulating cell death during ECM. Interestingly, we also observed that levels of two lipoxygenases, Aloxe3 and Alox12b, involved in leukotriene production (18) were elevated in mice with ECM. Given the observed role of leukotrienes in neural stem cell survival (66), it is conceivable that the implied increase in leukotriene production might have an antiapoptotic role. However, the leukotrienes could also mediate inflammatory responses in the brain that account in part for the pathogenesis of ECM (18). Other molecules with antiapoptotic roles that were uncovered in the comparisons between WT mice with ECM and resistant CD8-KO mice were Bcl2a1a and Bcl2a1b. Both of these could potentially act in synergy with p21 in inhibiting apoptosis. In terms of neural surface molecules, of particular interest was the downregulation in ECM of Unc5a, a netrin receptor that has previously been shown to be downregulated in response to neural injury (1). This, again, might represent a neural regenerative mechanism responding to the tissue damage triggered by the pathogenesis.
Expression changes suggestive of repair responses to cell damage in ECM.
A total of 17 genes involved in cytoskeleton organization and remodeling showed altered expression in ECM; the majority of them showed at least a twofold upregulation. These include β-actin and associated microfilament binding proteins such as coronin, cofilin-1, kelch-like-15, and transgelin, microtubular components such as tubulin-
, intermediate filament proteins such as nestin and lamin-A, and the tight-junction component claudin-11. These hint at a repair response to the tissue damage triggered by ECM pathogenesis. Importantly, the repair process mediated by these cytoskeleton genes might also be central to the observed rapid clinical reversibility of, and recovery from, CM. At least two molecules, β-actin and claudin-11, are components of intercellular tight junctions. Claudin-11, also known as oligodendrocyte-specific protein (OSP), contributes 7% of total protein in purified myelin (6). Studies with OSP-null mice suggest that this gene plays a pivotal role in the formation of intramyelinic tight junctions in the central nervous system (CNS). These mice also exhibit hind limb weakness and slowed conduction velocities in the CNS (19), characteristics that appear to be consistent with the neurological effects observed in CM.
Inference of regulatory pathways and cell interactions associated with ECM.
Our systematic analysis of the genes with altered expression in the course of ECM in the brain suggested that several regulatory pathways were being activated or modulated. These included transcription factors (TFs) suggestive of secondary-transcription-level regulatory cascades, as well as signaling proteins from diverse signaling cascades. Additionally, we also found numerous secreted and cell surface proteins involved in cell-cell and cell-matrix interactions whose genes showed altered expression. Indeed, in terms of the number of genes with altered expression patterns, the latter class was overrepresented among all functional classes (P < 0.001 by the
2 test).
(i) TFs.
Overall, the expression of 30 TFs and 11 chromatin proteins was found to be altered in mice with ECM in our various experiments. In moribund mice with ECM, expression of 11 TFs was modulated, and 7 of these were significantly upregulated relative to their levels in nonmoribund and resistant BALB/c mice. Compared to resistant CD8-KO mice, WT mice with ECM showed upregulation of 12 TFs and downregulation of a single TF. The downregulation of Stat1, a proinflammatory and proapoptotic TF, in resistant CD8-KO mice suggests that this TF might be involved in the tissue damage observed in ECM. Interestingly, several of the other TFs with altered mRNA levels, such as MafK, MafF, FosB, JunD1, Myc, Runx1, Rbpsuh, GATA2, and I
B-
(an NF-
B inhibitor), are related to hematopoiesis or lymphocyte-specific transcriptional responses. Specifically, Runx1 and Rbpsuh are critical for T-cell development and specialization, whereas FosB and JunD1 dimerize to constitute the AP1 transcription factor, required for downstream responses to T-cell-receptor-mediated signaling and interleukin-2 activation (24). Hence, the overexpression of the latter two genes might be related to the escalation of inflammation in ECM. The upregulation of GATA2 and MafK, both positive regulators of erythropoiesis, and the downregulation of Hox5a, an inhibitor of erythropoiesis (13), during ECM might be responses to the reduction in reticulocyte levels inferred above. The overexpression of MafK, a TF subunit important for normal platelet production (57), in ECM may also correlate with the previously characterized thrombocytopenia that occurs during ECM (2).
(ii) Signal transduction. The expression of a number of signaling proteins, including polypeptide first messengers, diverse receptors, intracellular enzymes such as kinases, phosphatases, and ubiquitin E3 ligase subunits, and adaptor proteins with peptide-binding and protein-interacting domains such as SH2, SH3, and SAM domains, appears to be altered in our experiments. We systematically analyzed these to identify potential pathways activated during ECM. First, we observed that resistant CD8-KO mice exhibited downregulation of four chemokines, namely, Ccl3, Ccl9, Ccl5, and Ccl25, in contrast to moribund WT mice. This suggested that part of the resistance to ECM in CD8-KO mice may have emerged from the inability to activate these chemokines, which probably trigger a neural inflammatory response (23). Our hypothesis is supported by an earlier study that showed elevated expression of Ccl5 (RANTES) in postmortem brain tissue samples from CM patients (53). Similarly, the overexpression in ECM of the adipokine Angptl4 is consistent with recent reports that it is induced upon brain injury, leading to hypoxia/ischemia (68). This first messenger might represent an as yet unexplored pathway involved in the response to oxygen deprivation in the morbid brain.
At the level of the cell membrane, P2Y12, one of two ADP receptors expressed on platelets, was significantly downregulated (–2.6- to –3-fold) during ECM. We were particularly interested in this molecule due to its role in platelet activation. In vessel injury, platelets are activated and play an important role in the arrest of bleeding (20). Grau and colleagues have conducted a series of experiments suggesting an important role for platelets and microparticles in the pathogenesis of ECM. In these studies, expression of P-selectin on endothelial cells (but not platelets) (10) and microparticle formation (11) are both associated with the development of ECM. Downregulation of P2Y12 could disrupt the formation of microparticles, and future studies exploring the effect of P2Y12 expression on microparticle formation during ECM may throw further light on the role of microparticles in the pathogenesis of ECM. In contrast, five other G protein-coupled seven-transmembrane-domain receptors, namely, Gprc5a, Olfr705, Bdkrb1, Npy6r, and V1rf3, are upregulated in mice with ECM. Of these, the action of the bradykinin receptor Bdkrb1 and the neuropeptide Y receptor Npy6r might be of particular interest in terms of understanding the action of neuropeptides on leukocyte and thrombocyte function, especially in relation to inflammatory responses (12). Strikingly, we observed that Toll-like receptor 11 (Tlr11) has lowered expression in mice with ECM relative to that in nonmoribund C57BL/6 and resistant BALB/c mice. It has recently been proposed that Tlr11 is a key mediator of immunity against apicomplexan parasites (70); hence, its lowered expression might have significance for the pathogenesis of ECM.
Genes of nine kinases and one phosphatase, as well as two regulatory phosphatase subunits, show altered expression in our data sets. The majority of kinases are upregulated in ECM, suggesting the activation of diverse phosphorylation responses. Among these, of interest is casein kinase 2 (CK2), a ubiquitous and conserved nuclear and cytoplasmic serine/threonine protein kinase involved in cell growth and proliferation and in the suppression of apoptosis. Upregulation of CK2 is consistent with activation of other antiapoptotic factors (see above) and points to a concerted host response via multiple independent pathways to minimize additional damage to brain tissue during ECM. The GTPase GBP2 is upregulated in mice with ECM. This molecule belongs to the family of IFN-inducible GTPases and is induced in macrophages by IFN-
(32, 37), which is consistent with the upregulation of the IFN-
gene observed in moribund mice compared to resistant CD8-KO mice. The strong overexpression of the candidate tumor suppressor gene SASH1, which encodes an adaptor molecule with SAM and SH3 domains, during ECM suggests that this molecule might be an uncharacterized component of a lymphocyte signaling pathway, consistent with previous observations on Btk– B cells (28).
|
|
|---|
Recently, several laboratories have utilized microarray technology to identify the host molecules associated with the pathogenesis of ECM in mice (15, 31, 35, 55). The major biomarkers of ECM identified by transcriptional analysis included immune response genes and molecules associated with apoptosis pathways (15, 31, 35). In one study that measured transcriptional alterations based on changes during the early stages of P. berghei infection in CM-susceptible and CM-resistant mice, detection of potent inflammatory responses, expression of molecules associated with brain injury repair, and early disturbances of brain metabolic energy metabolism were reported (14). In spite of these efforts, our knowledge regarding the molecules and molecular pathways that are linked with the pathology of CM remains incomplete, and further studies are needed to fill the critical lacuna. In this communication, we report the identification of several previously unknown host molecules and associated biochemical pathways that are potential mediators of the pathogenesis of ECM in the murine P. berghei ANKA model. These molecules were identified by using a combination of approaches, including comparative global gene expression profiling, bioinformatic analysis of protein function, and pathway reconstruction; for a few of the molecules considered biologically relevant, we determined a possible association between protein expression (by Western blotting) in brain tissue and real-time expression of ECM (for moribund versus nonmoribund mice) in susceptible WT C57BL/6 and resistant CD8-KO C57BL/6 mice.
The unique strength of our approach relied on determining the differential expression of ECM-associated genes in brain samples from moribund and nonmoribund C57BL/6 mice and from BALB/c mice, which are inherently resistant to ECM. In addition, to define the molecular mechanism underlying the pathogenic role of CD8 T cells and the perforin molecule during ECM, we compared gene expression profiles for ECM-susceptible C57BL/6 mice with ECM to those for resistant CD8-KO C57BL/6 mice and PFP-KO C57BL/6 mice during the cerebral phase of P. berghei ANKA infection.
Our comparative microarray analyses revealed novel and highly specific host responses; we found that 210 genes (66% upregulated, 34% downregulated) showed significantly altered expression in moribund WT C57BL/6 mice relative to that in both nonmoribund WT C57BL/6 mice and resistant BALB/c mice. In comparison, the resistance phenotypes of CD8-KO and PFP-KO mice were associated with alterations in only 166 and 9 genes, respectively. Our comprehensive computational sequence analysis of the genes with altered expression, followed by functional classification and identification of cellular interactions based on reconstruction of biological pathways, revealed that several critical biological processes were affected during the period when mice were experiencing the pathology of ECM. Notable among the host genes with ECM-associated responses are several molecules of the immune system, a wide variety of antiapoptotic factors, TFs, receptors, and proteins involved in signal transduction and cell-cell interaction. Several of the ECM-associated alterations in molecular pathways are reminiscent of septic-shock-like immunological responses and of cellular responses to severe stress and trauma-related injury (25, 29, 36). To offset these sudden and dramatic events, the host cells appear to mount vigorous cellular responses, particularly in the form of antiapoptotic and neuroprotectant molecules. Some of these important cellular responses related to ECM are discussed below.
A major finding of our study is the observation that there is a significant downregulation in the transcription of the hemoglobin
1 gene (Table 2) but a significant overexpression of the hemoglobin
1 protein in moribund compared to nonmoribund mice (Table 2; Fig. 1). We posit that a marked reduction in the reticulocyte counts along with a simultaneous influx of mature erythrocytes, possibly infected with Plasmodium parasites, may account for the increased levels of hemoglobin
1 protein in the brain during ECM. Importantly, this increase in hemoglobin
1 protein levels observed in our study may be a consequence of an increase in brain hemorrhaging and/or increased blood deposition in the brain resulting from parasite sequestration and an inflammation-induced decrease in blood flow. Hemoglobin
1 is further cannibalized by malaria parasites to produce heme, a molecule recently implicated as a major "host toxin" that triggers the inflammatory responses leading to the clinical syndromes of ECM (44). We also find that a significant portion of the host response is directed toward countering the damaging cellular events triggered due to the pathology of ECM. This contention is supported by results demonstrating that mice undergoing the clinical symptoms of ECM had increased transcriptional and translational expression of the MT-I and -II molecules (Table 2 and Fig. 1). In the brain, astrocytes are the main source of MT-I and -II (67), which might play an active role in reducing damage to the CNS by at least two mechanisms. First, MT-I and -II reduce the activation and recruitment of macrophages and T cells and also inhibit the production and secretion of inflammatory cytokines, including TNF-
and lymphotoxin (58). Second, MT-I and -II are powerful scavengers of free radicals (63) and therefore may moderate oxidative stress during CM. Importantly, exogenous treatment with MT-I and -II is also neuroprotective, is well tolerated, and has been shown to inhibit experimental autoimmune encephalomyelitis (46), the murine model of multiple sclerosis. These results suggest that the CNS response to damage inflicted by CM may bear similarities to neuroprotection induced by external stimuli or chronic illness.
Further evidence of vigorous host responses mounted in response to the deleterious effects of ECM was provided by significant alterations in the expression of the host genes involved in the regulation of both proapoptotic and antiapoptotic pathways (Table 2; see also Results). In our data set, levels of both proapoptotic and antiapoptotic genes were altered during ECM. Two recent studies also reported alterations in the expression of proapoptotic pathway genes during the cerebral phase in mice (31, 35). We focused our attention on p21 (the kinase inhibitor of the Cip/Kip family), a molecule important in the regulation of apoptosis. Cytoplasmic p21 has been shown to prevent apoptosis during neuronal differentiation (62) by inhibiting the catalytic activity of stress-activated protein kinase (56), and mitochondrial p21 prevents Fas-mediated apoptosis by inhibiting the activation of caspase-3 through direct interactions with its precursor, procaspase-3 (59, 60).
Compared to that in nonmoribund mice, expression of p21 in moribund mice is 5.5-fold increased at the level of transcription and 14.5-fold increased at the level of translation (Fig. 1). We think that the greatly enhanced expression of the p21 molecule during ECM represents a counterapoptotic response in the face of strong oxidative stress. Furthermore, differential expression of p21 in mice with ECM indicates that this molecule might serve as an important host biomarker of ECM. While it remains to be determined whether p21 expression is elevated in patients exhibiting the symptoms of CM, this molecule has been used to prognosticate certain cancers (52).
In summary, we believe that directed exploration of the data set generated for the ECM-susceptible and ECM-resistant mice and in-depth experimental analysis of predicted candidate molecules has improved our understanding of the repertoire of host molecules that mediate the pathology of CM. Further validation studies using brain biopsy specimens from patients who have succumbed to CM will be needed in order to determine the relevance of the identified molecules as predictors of CM in humans. Nonetheless, the presentation of this data set along with the accompanying "proof of concept" by Western blot studies will provide the research community with novel candidates for the development of diagnostics and therapeutics against CM, especially for young children in the areas of high P. falciparum transmission in Africa.
The views and opinions expressed here are those of the authors and should not be construed as the official opinion of the Food and Drug Administration.
Published ahead of print on 21 July 2008. ![]()
Supplemental material for this article may be found at http://iai.asm.org/. ![]()
|
|
|---|
and apoptotic cell death during experimental autoimmune encephalomyelitis (EAE). Exp. Neurol. 170:1-14.[CrossRef][Medline]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»