Remarkable stability in patterns of blood-stage gene expression during episodes of non-lethal Plasmodium yoelii malaria
© Cernetich-Ott et al.; licensee BioMed Central Ltd. 2012
Received: 24 January 2012
Accepted: 22 July 2012
Published: 6 August 2012
Microarray studies using in vitro cultures of synchronized, blood-stage Plasmodium falciparum malaria parasites have revealed a ‘just-in-time’ cascade of gene expression with some indication that these transcriptional patterns remain stable even in the presence of external stressors. However, direct analysis of transcription in P. falciparum blood-stage parasites obtained from the blood of infected patients suggests that parasite gene expression may be modulated by factors present in the in vivo environment of the host. The aim of this study was to examine changes in gene expression of the rodent malaria parasite, Plasmodium yoelii 17X, while varying the in vivo setting of replication.
Using P. yoelii 17X parasites replicating in vivo, differential gene expression in parasites isolated from individual mice, from independent infections, during ascending, peak and descending parasitaemia and in the presence and absence of host antibody responses was examined using P. yoelii DNA microarrays. A genome-wide analysis to identify coordinated changes in groups of genes associated with specific biological pathways was a primary focus, although an analysis of the expression patterns of two multi-gene families in P. yoelii, the yir and pyst-a families, was also completed.
Across experimental conditions, transcription was surprisingly stable with little evidence for distinct transcriptional states or for consistent changes in specific pathways. Differential gene expression was greatest when comparing differences due to parasite load and/or host cell availability. However, the number of differentially expressed genes was generally low. Of genes that were differentially expressed, many involved biologically diverse pathways. There was little to no differential expression of members of the yir and pyst-a multigene families that encode polymorphic proteins associated with the membrane of infected erythrocytes. However, a relatively large number of these genes were expressed during blood-stage infection regardless of experimental condition.
Taken together, these results indicate that 1) P. yoelii gene expression remains stable in the presence of a changing host environment, and 2) concurrent expression of a large number of the polymorphic yir and pyst-a genes, rather than differential expression in response to specific host factors, may in itself limit the effectiveness of host immune responses.
KeywordsMalaria Plasmodium yoelii Blood-stage parasites DNA microarrays In vivo gene expression profiles
Of the four species of Plasmodium that infect humans, Plasmodium falciparum is most commonly associated with severe disease and mortality, especially in young children. Immunity to this species does develop, but requires repeated exposure over many years. Sequencing and annotation of the P. falciparum genome revealed ~5,300 genes, many of which are predicted to encode proteins involved in host cell invasion and immune evasion, but the function of a large proportion remains unknown[2, 3]. A detailed microarray study of gene expression using blood-stage P. falciparum parasites cultured in vitro showed that a large percentage of the genome is expressed during the asexual stage, and that there is a surprising ‘continuous cascade’ of gene expression. Most genes are expressed only once during the asexual cycle and groups of genes involved in similar processes are active at the same time[4, 5]. This highly coordinated expression profile implies tightly controlled regulation, although as relatively few transcription factors have been identified in Plasmodium, precisely how gene regulation occurs remains incompletely understood. In these transcriptional studies, in vitro cultured P. falciparum parasites were tightly synchronized to evaluate gene expression at distinct points during the life cycle, and growth conditions were controlled to avoid introducing unwanted stressors and other confounding factors.
Transcriptional evaluation using malaria parasites isolated from infected individuals is more technically challenging and is influenced by the diversity of both host and parasite, their interaction and the ensuing immune response. Additional environmental factors such as body temperature, nutritional status and hormone levels also vary considerably between individuals and may influence parasite growth, gene expression and ultimately disease outcome. In one study, a transcriptional analysis of ring-stage parasites isolated directly from human subjects revealed distinct transcriptional profiles, thought by the authors to occur in response to the in vivo environment and not necessarily detectible using parasites cultured in vitro. There has been some debate as to the reasons for the observed differences. More such studies are necessary to fully examine parasite gene expression patterns in vivo and to determine how they change in response to host factors and immune pressures.
The genomes of Plasmodium parasites contain a number of multi-gene families, including the relatively well-studied var genes of P. falciparum which encode ~60 antigenic variants of P. falciparum erythrocyte membrane protein-1 (Pf EMP-1). Pf EMP-1 is expressed on the surface of the infected red blood cell (iRBC) membrane, and mediates binding to a variety of host endothelial cell receptors in tissues such as brain, placenta, lung and kidney. Interestingly, the var genes appear to be unique to P . falciparum and possibly the chimpanzee parasite Plasmodium reichenowi[10, 11] and are not present in other plasmodial species[2, 12–15]. However, the genomes of these species of malarial parasites also contain multi-gene families thought to encode iRBC surface proteins, although the function of these proteins remains largely unknown. It is thought that altering expression of members of a multi-gene family functions as an immune evasion strategy while maintaining the essential function of the encoded protein. The Plasmodium Interspersed Repeats (PIR) multi-gene family is one such family that is well conserved in the human malaria parasite Plasmodium vivax (vir, n = 245), the monkey parasite Plasmodium knowlesi (kir, n = 68) and the rodent malaria parasites Plasmodium berghei (bir, n = 245), Plasmodium chabaudi (cir, n = 135) and Plasmodium yoelii (yir, n = 838)[12, 14–17]. Unlike the clonal expression of P. falciparum var genes[18, 19], many pir genes are concurrently transcribed during blood-stage infection[20–22], and in the case of the yir family, this expression appears relatively stable during a primary infection. Sequencing of the P. yoelii genome revealed another multi-gene family designated pyst-a (n = 140), which are homologous to P. chabaudi glutamate-rich proteins and to a single hypothetical protein in P. falciparum (PF14_0604), P. vivax (Pvx_117290,) and P. knowlesi (PKH_124210,). The glutamate-rich proteins of P. chabaudi, such as Pc90 (also known as Pc(em)93, Pc(em)96 and Pch105/RESA) are thought to localize to the cytoplasmic face of the RBC membrane and have been shown to be immunogenic[23–29]. To date, the function of the PYST-A family of proteins during blood-stage infection in P. yoelii has not been studied.
Rodent models of malaria present a unique opportunity to examine patterns of gene expression in vivo, allowing for use of cloned parasite strains replicating in inbred mice. These conditions may reduce some variability seen when using clinical isolates of P. falciparum parasites and are technically much more manageable. These in vivo model systems also allow for examination of gene expression at various time points post-infection and make it possible to manipulate the host immune system under controlled conditions. In this study, the P. yoelii 17X murine model was utilized to analyse how parasite gene expression patterns differ between individual animals infected with the same cloned parasite and how expression may change across independent infections initiated with different parasite populations. Gene expression at various time points post-infection and in the presence and absence of host immune pressure was also analyzed. A genome wide analysis was the primary focus but yir and pyst-a gene expression under these conditions was also examined. Results indicate that P. yoelii gene expression was remarkably stable supporting the notion that in vivo, external stressors and environmental factors have limited influence on gene transcription in Plasmodium blood-stage parasites.
All animal studies were reviewed, approved and conducted in compliance with the Institutional Animal Care and Use Committee (IACUC) of Drexel University College of Medicine (protocol approval ID #18874). This IACUC operates with Public Health Service approval (Animal Welfare Assurance Number: A-322-01).
Mice and parasites
Five to six week old male BALB/cByJ mice were purchased from The Jackson Laboratory (Bar Harbor, Maine, USA). Male BALB/cJ and B-cell-deficient JHD mice on a BALB/cJ background were purchased from Taconic Farms Inc. (Germantown, NY, USA). All animals were housed in the Animal Care Facility of Drexel University College of Medicine under specific pathogen-free conditions. Food and water were provided ad libitum and the room was maintained on a 12-hour light–dark cycle. The lethal and non-lethal 17XL and 17X strains of P. yoelii were originally obtained from Dr William P Weidanz (University of Wisconsin, Madison, WI, USA) and maintained as cryopreserved stabilates. Blood from P. yoelii -infected mice contains a mixture of ring, trophozoite and schizont stage parasites as replication is asynchronous. During infection with these P. yoelii lines, gametocytes are not observed and parasitized blood is not infective for mosquitoes.
To establish cloned lines of P. yoelii 17X, groups of three to four male BALB/cByJ mice were infected i.v. with 25, 10 or one P. yoelii 17X iRBC(s). Blood parasitaemia was monitored by the enumeration of parasitized erythrocytes in thin tail-blood smears stained with Giemsa (EMD Chemicals, Inc., Gibbstown, NJ, USA). Parasites from the lowest infective dose were selected as clones.
Py RMP preparation and generation of anti-Py RMP sera
Male BALB/cByJ mice were infected with 1x106P. yoelii 17X iRBCs and blood was obtained when parasitaemia was ~30%. Infected RBCs were separated from uninfected RBCs by density gradient centrifugation on a 70% Percoll gradient (GE Healthcare, Uppsala, Sweden). Recovered iRBCs were treated with PBS-0.05% saponin and erythrocyte membranes were separated from intact parasites by differential centrifugation. Membranes from uninfected RBCs were also isolated as a control. Protein concentrations were determined using the bicinchoninic acid protein assay (Pierce Chemical Company, Rockford, IL, USA). The P. yoelii reticulocyte membrane protein (Py RMP) fraction, uninfected membrane proteins and parasite-associated antigen were separated by SDS-PAGE (15 μg/lane) on a 10% polyacrylamide gel and stained with Coomassie Blue. For the generation of polyclonal anti-Py RMP sera, mice were immunized subcutaneously three times at three-week intervals with 50 μg/dose of the Py RMP preparation formulated with 25 μg Quil A adjuvant (Accurate Chemical and Scientific Corporation, Westbury, NY, USA). Two weeks following the third immunization, sera was collected.
The Py RMP fraction, uninfected membrane proteins and parasite-associated antigen were separated by SDS-PAGE as above and electroblotted onto nitrocellulose membranes. Membranes were then blocked with 5% (w/v) non-fat milk in TBS (25 mM Tris–HCl, pH 8.0, 150 mM NaCl) and probed with normal rabbit serum (1:20,000), polyclonal rabbit antiserum raised against recombinant P. yoelii MSP-8 (1:20,000), normal mouse serum (1:1000), or mouse anti-Py RMP (1:1000) diluted in TBS containing 0.1% (v/v) Tween 20 and 1% (w/v) BSA. Bound antibodies were detected by chemiluminescence using horseradish peroxidase conjugated protein A (Pierce Chemical Company) or rabbit anti-mouse IgG (Invitrogen, Carlsbad, CA, USA) and the SuperSignal West Pico Substrate (Pierce Chemical Company).
Blood was collected from P. yoelii 17X infected mice, washed and resuspended at ~25% haematocrit in PBS containing 1% (w/v) gelatin. Thin blood films were prepared, air dried and fixed in acetone:methanol (1:1) for 20 min at −20°C. Fixed smears were incubated for 30 min at 37°C in a humidified chamber with anti-Py RMP sera or normal mouse sera diluted 1:200 in PBS. Bound antibody was detected using tetramethyl rhodamine (TRITC)-conjugated goat anti-mouse immunoglobulin G (IgG) (Invitrogen) diluted 1:250 in PBS. Slides were stained with Hoechst 33258 diluted 1:47,000 in 1X PBS, washed, and mounted with AntiFade reagent (Invitrogen Molecular Probes, Eugene OR, USA). Images (1,000X magnification) were obtained using an Olympus BX60 fluorescent microscope (Olympus America Inc., Melville, NY, USA) and a SPOT RT Slider Digital Camera System (Diagnostic Instruments, Sterling Heights, MI, USA).
Immunizations and infections
Normal and immunodeficient naïve mice were infected by intraperitoneal injection of 1x105P. yoelii 17X infected RBCs. Infected RBCs for RNA isolation were obtained early in infection when parasitaemia was ~15% (day 10–11), at peak parasitaemia (~35-40%, day 14) or during descending parasitaemia (~15%, day 18). Alternatively, mice were immunized subcutaneously three times at three week intervals with 50 μg of the Py RMP fraction formulated with 25 μg of Quil A adjuvant. Control groups received Quil A alone. Two weeks following the final immunization, animals were challenged with 1x105P. yoelii 17X iRBCs. Infected RBCs for RNA isolation were obtained from Py RMP-immunized and adjuvant control mice on day 10–12 of infection.
Plasmodium yoelii DNA microarrays
Plasmodium yoelii DNA microarrays were produced in the Molecular Genomics Core Facility, Drexel University College of Medicine, under the direction of LWB. Each array contained 65-base oligonucleotides, spotted in duplicate, representing ~6,700 coding regions predicted from the P. yoelii genome sequence analysis and present in the current PlasmoDB and NCBI databases. Each quadrant of an array contained a pool of P. yoelii oligos (~2200), spotted in triplicate, which served as a positive control. Blood from P. yoelii -infected mice was obtained from individual animals or pooled from groups of animals (five to 10) as needed. Infected RBCs were isolated as above by Percoll density gradient centrifugation. Infected RBCs were saponin lysed and pelleted P. yoelii parasites were resuspended directly in the TRIzol Reagent (Invitrogen). Total RNA was extracted, precipitated and purified using an RNeasy RNA isolation kit (Qiagen, Inc., Valencia, CA, USA).
For all microarrays, gene expression in P. yoelii 17X blood-stage parasites was evaluated relative to a standard comparator of purified P. yoelii 17XL total RNA. P. yoelii RNA (5 μg/sample) was amplified in the presence of aminoallyl-dUTP using the Amino Allyl MessageAmp II aRNA amplification kit (Ambion, Inc., Austin, TX, USA). aRNA was then fluorescently labelled by reaction with monofunctional, NHS-activated Cy3 or Cy5 dyes (Amersham Biosciences Inc., Piscataway, NJ, USA). Cy dye labelled aRNA was purified with yield and specific activity of each probe determined by absorption spectroscopy. Pairs of Cy3 and Cy5 labelled aRNA probes (1 μg/probe) were pooled, fragmented (Ambion RNA Fragmentation Reagent) and hybridized to the P. yoelii microarrays for 14–16 hours at 65°C. Following hybridization and washing, slides were scanned using a GenePix 4000A microarray laser scanner (Axon Instruments, Inc., Union City, CA, USA) and the fluorescence intensity of each DNA feature was determined at 532 nm (Cy3) and 635 nm (Cy5). Data for each gene was obtained from replicate features on each array. In addition, replicate arrays and standard dye flips were also hybridized for each comparison. Fluorescence data were acquired and initially analysed using GenePixPro 5.1 Software (Axon Instruments, Inc.). Irregular and missing features flagged during image acquisition, features with a diameter of ≤60 μm and features with a signal to noise ratio (SNR) of <2 on the 532 nm channel and SNR <2 on the 635 nm channel were removed from the analysis.
In order to evaluate differential gene expression in P. yoelii 17X parasites replicating in various in vivo environments, Cy3 and Cy5 fluorescence intensity data were processed and analysed using limma version 2.18.2[31, 32] through limmaGUI version 1.20.0 in the R programming environment The normexp method with an offset value of 16 was utilized for background correction[35, 36]. Data were normalized using the print-tip loess method. After fit of a linear model to the expression data, genes differentially expressed by at least two-fold were identified (log2-fold change >1 or < −1). Such changes in gene expression were considered significant based on a p -value ≤ 0.01 obtained from the moderated t-statistic and adjusted for multiple testing using the method of Benjamini and Hochberg. The false discovery rate was 1%. An initial data analysis revealed relatively small numbers of differentially expressed genes across conditions. As such, the final statistical analysis was completed by multiple pairwise comparisons of groups, accepting the increased possibility of a type I error - erroneously detecting differences in gene expression.
The number of yir and pyst-a genes expressed in a given population of P. yoelii 17X parasites was estimated based on normalized signal intensity relative to the entire P. yoelii gene set. Fluorescence data were captured and filtered as described above and analysed using the Acuity 3.1 Microarray Informatics Software package (Axon instruments, Inc). On each array, Cy3 and Cy5 signal intensities were normalized such that the ratio of the median fluorescence intensity of 192 positive control features equalled 1. A normalized fluorescence intensity for all genes was calculated by setting the fluorescence intensity of the positive control to 3,000. Mean fluorescence intensity (MFI) was then calculated for each gene across an array set (i.e., replicate arrays, dye swaps) and percentiled. Based on mapping data available through PlasmoDB, 464/859 yir gene oligos and 74/140 pyst-a gene oligos on the arrays detect unique transcripts (Additional file1). Non-discriminating yir and pyst-a oligos were eliminated from the analysis. Expression of ‘unique’ genes was categorized as below detection (0-50th percentile), low (51-75th percentile), moderate (76-90th percentile) or high (>90th percentile).
Plasmodium yoelii gene/protein categories
Little variation in gene expression in Plasmodium yoelii 17X blood-stage parasites obtained from individual mice challenged with the same inoculum of iRBCs
Differential gene expression in Plasmodium yoelii 17X blood-stage parasites marginally increases over time
Gene expression profiles in Plasmodium yoelii 17X blood-stage parasites remain relatively stable in independently infected mice
Plasmodium yoelii gene expression is influenced by parasitaemia and/or reticulocyte levels
To further evaluate changes in parasite gene expression during the course of infection, P. yoelii 17X RNA was isolated on day 10 (19% parasitaemia, 18.6% reticulocytes) and on day 18 (14% parasitaemia, 70.2% reticulocytes). Although parasite load is similar at these time points, parasitaemia is ascending on day 10 and target host cells are relatively scarce. By day 18 post-infection, parasitaemia is declining due to immune-mediated mechanisms and reticulocytes are abundant. Comparing day 10 and day 18, there were 192 differentially expressed P. yoelii genes, including 67 genes of unknown function (Figure6 and Additional file6). Of those differentially expressed, 94 were up-regulated and 98 were down-regulated on day 18 versus day 10. This set of differentially expressed genes represents the most variable expression seen across all experimental conditions. Of interest, genes encoding proteins associated with the cytoplasmic translation machinery tended to be down-regulated on day 18 (n = 11) versus day 10, and included several ribosomal proteins. In addition, a small set of genes (n = 24) encoding proteins involved with schizont rupture and/or merozoite invasion of RBCs tended to be up-regulated on day 18 and these included glidosome- and rhoptry-associated genes, MSP1 and MSP7. Finally, there were 17 genes encoding apicoplast/mitochondrion-associated proteins that were differentially expressed, with up and down-regulated genes equally represented. As before, members the yir and pyst-a multigene families were not differentially expressed while the number of genes in each family expressed on day 18 remained stable (Figure3, 155 yir genes, 37 pyst-a genes).
Changes in immune pressure do not significantly alter Plasmodium yoelii 17X gene expression during a primary infection
To specifically evaluate whether the presence or absence of host immune pressure affects parasite gene expression in vivo, two experimental approaches were taken. In the first approach, all antibody-dependent immune pressure was removed by use of B cell deficient mice. Naïve, immunocompetent BALB/c mice and B cell deficient, JHD mice on a BALB/c background were concurrently infected with P. yoelii 17X iRBCs obtained from the same donor animal. P. yoelii gene expression in the two hosts was evaluated when ascending parasitaemia reached ~15%. Surprisingly, changes in gene expression were negligible with only seven differentially expressed genes, five of which were of unknown function (Figure6 and Additional file6). No yir or pyst-a genes were differentially expressed between the populations of P. yoelii 17X parasites.
Transcription in P. falciparum and P. vivax parasites appears to be tightly regulated, resulting in a ‘continuous cascade’ of gene expression during blood-stage development[4, 5, 45]. There is some evidence that this transcriptional programme continues unabated even in the presence of external stresses. In one study, anti-folate treatment of cultures of drug-sensitive P. falciparum blood-stage parasites did not lead to increased expression of the defined drug target, dihydrofolate reductase-thymidylate synthatase. Additionally, there were no genes differentially regulated greater than two-fold between control and drug treated cultures. Similarly, Gunasekera et al. found little variation in transcription patterns between untreated and chloroquine-treated P . falciparum parasites growing asynchronously in culture. These results were unexpected, as genome-wide studies in organisms such as Saccharomyces cerevisiae, Mycobacterium tuberculosis and Candida albicans revealed differential expression of hundreds of genes in response to drug treatment, and these changes often confirmed or revealed associated drug targets and affected pathways[48–51]. However, other studies in Plasmodium examining the effects of drug treatment or manipulation of culture conditions have found more substantial transcriptional changes (~300-400 differentially expressed genes), especially in genes predicted to encode proteins exported to the host cell cytoplasm and the RBC membrane[52, 53]. The work shown here extends these efforts significantly by focusing on malaria parasites replicating in vivo, in a more complex setting. With P. yoelii 17X blood-stage parasites, global gene expression profiles were remarkably consistent even in the presence of a changing host environment.
Initially, three genetically identical, age-matched BALB/cByJ mice were infected with the same inoculum of P. yoelii 17X iRBCs obtained from a donor mouse. P. yoelii RNA was isolated from each animal after a 10–11 day period of replication in vivo. As the in vivo environment is expected to be comparable in this setting, only a small set of differentially expressed genes was anticipated. Parasite gene expression patterns in each of the three animals were virtually identical with only five to six differentially expressed genes. Between the initiation of infection and harvest of iRBCs for RNA isolation, P. yoelii parasites in this study completed 10–12 replication cycles in vivo. Over this time, the number of differentially expressed genes was expected to increase somewhat relative to the population of P. yoelii iRBCs used to initiate the infection. This was the case, but the increase was relatively small involving only 40–60 genes out of the set represented by ~ 6,700 oligonucleotides on the arrays. A small cluster of genes expressed in late-stage parasites and potentially involved in merozoite invasion of host cells were differentially expressed between the donor mouse and recipients. This may not be surprising as iRBCs were harvested early during ascending parasitaemia when P. yoelii 17X parasites can be found in varying degrees in both normocytes and reticulocytes. This change in gene expression may reflect utilization of multiple invasion pathways. Finally, a number of potentially influencing variables were further increased by comparing gene expression profiles of P. yoelii 17X blood-stage parasites harvested from four independent infections initiated with iRBCs obtained from separate donor mice. Remarkably, the number of differentially expressed genes in pair-wise comparisons remained low (10–50 genes). These represent genes that are involved in diverse or unknown functions and at present, the associated changes do not appear to reflect a biologically significant response to the host environment. This, combined with a false discovery rate of 1% associated with this analysis suggest that detection of a low number of differentially expressed but largely unrelated genes could simply represent ‘noise’ inherent to large microarray analysis.
Non-lethal P. yoelii 17X parasites replicate preferentially in reticulocytes but do invade and develop within mature RBCs when reticulocytes are limiting. Parasites harvested on day 10 of P. yoelii 17X infection will be present in both reticulocytes and normocytes. As a result of an influx of new reticulocytes into circulation, P. yoelii 17X parasites will be found almost exclusively in reticulocytes by day 14 post-infection. As such, substantial differences in patterns of P. yoelii 17X gene expression in parasites harvested on day 10 versus day 14 were expected. In addition to the shift in host cell tropism, the ongoing infection increases parasite burden, host stress and immune pressure that could also influence P. yoelii 17X expression profiles. While the number of differentially expressed genes in this comparison increased to 95, these were not clustered based on related function or biological process and the majority were of unknown function. In an early study, utilizing lethal P. yoelii 17XL parasites, a similar number of genes altered in association with a shift in host cell preference from normocytes to reticulocytes were identified. The function of these genes will be revisited as the annotation of both P. falciparum and P. yoelii genomes progresses. To further examine the role of host cell availability and immune pressure, gene expression in parasites isolated on day 10 and day 18 of infection was compared. Parasitaemia at these two time points is comparable. However on day 18, 70-75% of RBCs in circulation are reticulocytes and P. yoelii iRBCs are rapidly being cleared from circulation by immune-mediated mechanisms. Nearly 200 genes were differentially expressed in P. yoelii iRBCs on day 18 relative to day 10 of infection; the largest set in this study. On day 18, several genes encoding components of the cytoplasmic translation machinery were down-regulated while a second set associated with merozoite invasion were up-regulated. In P. falciparum, transcription of genes encoding the cytoplasmic translation machinery generally peaks in ring-stage parasites at ~12 hours post-invasion, followed by a marked down-regulation in trophozoite- and schizont-stage parasites. In contrast, transcription of genes involved in schizont rupture and merozoite invasion is relatively low through most of the asexual cycle, peaking at ~42 hours post-invasion[4, 5]. Differences noted in gene expression in P. yoelii parasites on day 18 versus day 10 of infection may partially be explained by a modest shift noted in distribution of parasite stages in circulation toward increased schizonts and decreased trophozoites. It is also clear that infected blood isolated on day 18 contains damaged host cells and dying parasites, as well as free extracellular parasites. These factors likely contributed to a mixed alteration in expression profiles at this late time point during infection.
To focus on the role of the host immune response on parasite gene expression, two approaches to actively alter host immune pressure were taken. In the first approach, all antibody mediated immune pressure was removed by infecting B cell deficient JHD mice. It was expected that in the absence of antibodies, P. yoelii growth would be unchecked with an increase in the diversity of parasite populations and gene expression profiles. This did not occur as only seven genes were differentially expressed when P. yoelii parasites isolated from B cell deficient JHD mice were compared with those from immunocompetent control mice. In the second approach, immune pressure on P. yoelii 17X parasites was increased by immunizing mice with a preparation of parasite proteins associated with the membrane of infected reticulocytes. Here, immunization-induced antibodies were expected to select for the growth of parasites with altered gene expression patterns that included changes in members of the yir and pyst-a multigene families that encode proteins localized to the erythrocyte membrane. Unexpectedly, this did not occur as only 10 genes were differentially expressed. Even in these situations where the host immune response was experimentally manipulated, P. yoelii 17X gene expression proceeded unaltered. This is very reminiscent of the ‘hard-wired’ programme of gene expression observed with in vitro cultures of P. falciparum[4, 5, 46].
In addition to addressing global changes in gene expression, the expression of members of two multi-gene families in P. yoelii, the yir and pyst-a were of particular interest. Previous studies have shown that at a population level, many yir genes are transcribed during blood-stage infection, in a seemingly random order. Transcription at the level of individual parasites appears to be tightly controlled, with each parasite transcribing between one and three yir genes. In the present study, significant differences in the expression of yir family members across experimental conditions were not observed, including the comparison of P. yoelii 17X day 10 and day 18 parasites. This is in agreement with the study by Cunningham et al. examining the expression of a subset of the yir repertoire on days 12 and 18 post-infection. Of note, expression of a fairly large number of the yir genes (89–155) was consistently observed using signal intensity greater than the 50th percentile on each array as a benchmark. These expressed yir genes were distributed across five previously reported phylogenetic groups, with no obvious bias toward any one group.
The expression of pyst-a genes followed a similar pattern to the yir family, with little to no differential expression between conditions. It appears that a large proportion of pyst-a family members are expressed in a mixed blood-stage infection. Whether pyst-a gene expression is more restricted in individual parasites has not been determined. No yir or pyst-a members were differentially expressed in wild-type versus JHD animals, in agreement with a study examining expression of a subset of yir members in wild-type versus Rag2 knockout mice, which lack mature B and T cells. A more extensive analysis of yir and pyst-a expression during primary and secondary infections using a microarray approach would be of interest. It has been suggested that the polymorphic erythrocyte membrane antigens encoded by the yir and/or pyst-a multigene families may function as part of a parasite immune evasion strategy. For immune evasion, the current data suggest that changing the pattern of gene expression may be less important than showering the host immune system with a large repertoire of polymorphic antigens at any given time. In fact, prior immunization of mice with the Py RMP preparation may have further impeded the development of protective immune responses following challenge infection as these animals developed a persistent infection which was difficult to clear. In agreement with the current findings, members of the cir multi-gene family do not seem to be differentially expressed during the course of blood-stage P. chabaudi infection, but there is some indication that cir gene expression varies in parasites localized to different host tissues. However, it is still possible that much like the P. falciparum EMP1 and the var gene family, the yir and/or pyst-a encoded proteins may possess specific functions required for parasite growth and development in vivo.
Overall, in vivo, P. yoelii 17X gene expression did not appear to be appreciably influenced by the host environment. These data in the P. yoelii 17X model differ from that reported by Daily et al. who detected distinct patterns of gene expression in P. falciparum parasites isolated from malaria-infected patients. In light of these data, changes in the expression of genes encoding a set of mitochondrial proteins (n = 148), across all experimental conditions was examined. Significant changes that would be consistent with distinct physiological states or a response to environmental stress were not observed (Additional file3). Although parasite gene expression was assessed in vivo, the P. yoelii model is still less complex than with P. falciparum infected human subjects. Here, a single cloned line of P. yoelii 17X that does not produce gametocytes was utilized to simultaneously infect genetically identical mice housed in environmentally controlled, specific pathogen free conditions. In some regards, this lack of gametocyte-stage parasites is advantageous, as detection of differential gene expression due to varying gametocyte levels across samples can be ruled out. On the other hand, the repeated passage of this P. yoelii 17X line in the vertebrate host and/or the lack of exposure to the mosquito vector could have altered regulatory mechanisms that control expression of certain gene subsets (i.e. yir or pyst-a genes). It is also possible that the greater diversity in human hosts (genetic, environmental) may have a greater influence on parasite gene expression than we observed in the P. yoelii model. Finally, each isolate of P. falciparum may exhibit unique elements of a ‘hard-wired’ programme of gene expression that can be detected when comparing isolates obtained from individual malaria patients. Co-infection of human subjects with such distinct P. falciparum clones will increase diversity in the overall P. falciparum gene expression profile detected in a single host and may allow preferential growth of P. falciparum clones in specific in vivo environments.
The present study shows that gene expression patterns of blood-stage P. yoelii 17X parasites replicating in vivo are very stable and do not vary appreciably even at ascending and peak parasitaemia and in the presence or absence of host immune pressure. These findings support the notion of a ‘hard wired’ gene expression pattern in plasmodial species. Additionally, although a relatively large number of yir and pyst-a are expressed across the conditions studied here, there was no evidence of significant differential expression within these multi-gene families. Concurrent expression of a large number of polymorphic antigens in blood-stage malaria parasites may be a strategy to impede the development of protective immune responses.
Conserved, amino acid/carbohydrate metabolism
Conserved, apicoplast, mitochondrion
Conserved, cell division, DNA replication
Conserved, cytoplasmic translation machinery
Conserved, intracellular transport/vesicles
Conserved, protein folding, modification, degradation
Conserved, ribonucleotide synthesis
Conserved, schizont, merozoite, invasion
Conserved, structural protein
Conserved, transcription machinery
Conserved, unknown function
Rodent malaria parasites, multigene family
Rodent malaria parasites, unknown function
P. yoelii specific, pyst-a family
P. yoelii specific, yir family
P. yoelii specific, unknown function.
This work was supported by NIH-NIAID grant AI069147 (JMB). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
- Baird JK: Host age as a determinant of naturally acquired immunity to Plasmodium falciparum. Parasitol Today. 1995, 11: 105-111. 10.1016/0169-4758(95)80167-7.View ArticlePubMedGoogle Scholar
- Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S, Paulsen IT, James K, Eisen JA, Rutherford K, Salzberg SL, Craig A, Kyes S, Chan MS, Nene V, Shallom SJ, Suh B, Peterson J, Angiuoli S, Pertea M, Allen J, Selengut J, Haft D, Mather MW, Vaidya AB, Martin DM, Fairlamb AH, Fraunholz MJ, Roos DS, Ralph SA, McFadden GI, Cummings LM, Subramanian GM, Mungall C, Venter JC, Carucci DJ, Hoffman SL, Newbold C, Davis RW, Fraser CM, Barrell B: Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002, 419: 498-511. 10.1038/nature01097.View ArticlePubMedGoogle Scholar
- Aurrecoechea C, Brestelli J, Brunk BP, Dommer J, Fischer S, Gajria B, Gao X, Gingle A, Grant G, Harb OS, Heiges M, Innamorato F, Iodice J, Kissinger JC, Kraemer E, Li W, Miller JA, Nayak V, Pennington C, Pinney DF, Roos DS, Ross C, Stoeckert CJ, Treatman C, Wang H: PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res. 2009, 37: D539-543. 10.1093/nar/gkn814.PubMed CentralView ArticlePubMedGoogle Scholar
- Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 2003, 1: e5-PubMed CentralView ArticlePubMedGoogle Scholar
- Llinas M, Bozdech Z, Wong ED, Adai AT, DeRisi JL: Comparative whole genome transcriptome analysis of three Plasmodium falciparum strains. Nucleic Acids Res. 2006, 34: 1166-1173. 10.1093/nar/gkj517.PubMed CentralView ArticlePubMedGoogle Scholar
- LeRoux M, Lakshmanan V, Daily JP: Plasmodium falciparum biology: analysis of in vitro versus in vivo growth conditions. Trends Parasitol. 2009, 25: 474-481. 10.1016/j.pt.2009.07.005.View ArticlePubMedGoogle Scholar
- Daily JP, Scanfeld D, Pochet N, Le Roch K, Plouffe D, Kamal M, Sarr O, Mboup S, Ndir O, Wypij D, Levasseur K, Thomas E, Tamayo P, Dong C, Zhou Y, Lander ES, Ndiaye D, Wirth D, Winzeler EA, Mesirov JP, Regev A: Distinct physiological states of Plasmodium falciparum in malaria-infected patients. Nature. 2007, 450: 1091-1095. 10.1038/nature06311.View ArticlePubMedGoogle Scholar
- Lemieux JE, Gomez-Escobar N, Feller A, Carret C, Amambua-Ngwa A, Pinches R, Day F, Kyes SA, Conway DJ, Holmes CC, Newbold CI: Statistical estimation of cell-cycle progression and lineage commitment in Plasmodium falciparum reveals a homogeneous pattern of transcription in ex vivo culture. Proc Natl Acad Sci U S A. 2009, 106: 7559-7564. 10.1073/pnas.0811829106.PubMed CentralView ArticlePubMedGoogle Scholar
- Kraemer SM, Smith JD: A family affair: var genes, PfEMP1 binding, and malaria disease. Curr Opin Microbiol. 2006, 9: 374-380. 10.1016/j.mib.2006.06.006.View ArticlePubMedGoogle Scholar
- Trimnell AR, Kraemer SM, Mukherjee S, Phippard DJ, Janes JH, Flamoe E, Su XZ, Awadalla P, Smith JD: Global genetic diversity and evolution of var genes associated with placental and severe childhood malaria. Mol Biochem Parasitol. 2006, 148: 169-180. 10.1016/j.molbiopara.2006.03.012.View ArticlePubMedGoogle Scholar
- Bull PC, Buckee CO, Kyes S, Kortok MM, Thathy V, Guyah B, Stoute JA, Newbold CI, Marsh K: Plasmodium falciparum antigenic variation. Mapping mosaic var gene sequences onto a network of shared, highly polymorphic sequence blocks. Mol Microbiol. 2008, 68: 1519-1534. 10.1111/j.1365-2958.2008.06248.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Carlton JM, Angiuoli SV, Suh BB, Kooij TW, Pertea M, Silva JC, Ermolaeva MD, Allen JE, Selengut JD, Koo HL, Peterson JD, Pop M, Kosack DS, Shumway MF, Bidwell SL, Shallom SJ, van Aken SE, Riedmuller SB, Feldblyum TV, Cho JK, Quackenbush J, Sedegah M, Shoaibi A, Cummings LM, Florens L, Yates JR, Raine JD, Sinden RE, Harris MA, Cunningham DA, Preiser PR, Bergman LW, Vaidya AB, van Lin LH, Janse CJ, Waters AP, Smith HO, White OR, Salzberg SL, Venter JC, Fraser CM, Hoffman SL, Gardner MJ, Carucci DJ: Genome sequence and comparative analysis of the model rodent malaria parasitePlasmodium yoelii yoelii. Nature. 2002, 419: 512-519. 10.1038/nature01099.View ArticlePubMedGoogle Scholar
- Hall N, Karras M, Raine JD, Carlton JM, Kooij TW, Berriman M, Florens L, Janssen CS, Pain A, Christophides GK, James K, Rutherford K, Harris B, Harris D, Churcher C, Quail MA, Ormond D, Doggett J, Trueman HE, Mendoza J, Bidwell SL, Rajandream MA, Carucci DJ, Yates JR, Kafatos FC, Janse CJ, Barrell B, Turner CM, Waters AP, Sinden RE: A comprehensive survey of the Plasmodium life cycle by genomic, transcriptomic, and proteomic analyses. Science. 2005, 307: 82-86. 10.1126/science.1103717.View ArticlePubMedGoogle Scholar
- Pain A, Bohme U, Berry AE, Mungall K, Finn RD, Jackson AP, Mourier T, Mistry J, Pasini EM, Aslett MA, Balasubrammaniam S, Borgwardt K, Brooks K, Carret C, Carver TJ, Cherevach I, Chillingworth T, Clark TG, Galinski MR, Hall N, Harper D, Harris D, Hauser H, Ivens A, Janssen CS, Keane T, Larke N, Lapp S, Marti M, Moule S, Meyer IM, Ormond D, Peters N, Sanders M, Sanders S, Sargeant TJ, Simmonds M, Smith F, Squares R, Thurston S, Tivey AR, Walker D, White B, Zuiderwijk E, Churcher C, Quail MA, Cowman AF, Turner CM, Rajandream MA, Kocken CH, Thomas AW, Newbold CI, Barrell BG, Berriman M: The genome of the simian and human malaria parasite Plasmodium knowlesi. Nature. 2008, 455: 799-803. 10.1038/nature07306.PubMed CentralView ArticlePubMedGoogle Scholar
- Carlton JM, Adams JH, Silva JC, Bidwell SL, Lorenzi H, Caler E, Crabtree J, Angiuoli SV, Merino EF, Amedeo P, Cheng Q, Coulson RM, Crabb BS, Del Portillo HA, Essien K, Feldblyum TV, Fernandez-Becerra C, Gilson PR, Gueye AH, Guo X, Kang'a S, Kooij TW, Korsinczky M, Meyer EV, Nene V, Paulsen I, White O, SA R, Ren Q, Sargeant TJ, Salzberg SL, Stoeckert CJ, Sullivan SA, Yamamoto MM, Hoffman SL, Wortman JR, Gardner MJ, Galinski MR, Barnwell JW, Fraser-Liggett CM: Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature. 2008, 455: 757-763. 10.1038/nature07327.PubMed CentralView ArticlePubMedGoogle Scholar
- Janssen CS, Barrett MP, Turner CM, Phillips RS: A large gene family for putative variant antigens shared by human and rodent malaria parasites. Proc Biol Sci. 2002, 269: 431-436. 10.1098/rspb.2001.1903.PubMed CentralView ArticlePubMedGoogle Scholar
- Janssen CS, Phillips RS, Turner CM, Barrett MP: Plasmodium interspersed repeats: the major multigene superfamily of malaria parasites. Nucleic Acids Res. 2004, 32: 5712-5720. 10.1093/nar/gkh907.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Q, Fernandez V, Sundstrom A, Schlichtherle M, Datta S, Hagblom P, Wahlgren M: Developmental selection of var gene expression in Plasmodium falciparum. Nature. 1998, 394: 392-395. 10.1038/28660.View ArticlePubMedGoogle Scholar
- Scherf A, Hernandez-Rivas R, Buffet P, Bottius E, Benatar C, Pouvelle B, Gysin J, Lanzer M: Antigenic variation in malaria: in situ switching, relaxed and mutually exclusive transcription of var genes during intra-erythrocytic development in Plasmodium falciparum. EMBO J. 1998, 17: 5418-5426. 10.1093/emboj/17.18.5418.PubMed CentralView ArticlePubMedGoogle Scholar
- Cunningham DA, Jarra W, Koernig S, Fonager J, Fernandez-Reyes D, Blythe JE, Waller C, Preiser PR, Langhorne J: Host immunity modulates transcriptional changes in a multigene family (yir) of rodent malaria. Mol Microbiol. 2005, 58: 636-647. 10.1111/j.1365-2958.2005.04840.x.View ArticlePubMedGoogle Scholar
- Ebbinghaus P, Krucken J: Characterization and tissue-specific expression patterns of the Plasmodium chabaudi cir multigene family. Malar J. 2011, 10: 272-10.1186/1475-2875-10-272.PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandez-Becerra C, Pein O, de Oliveira TR, Yamamoto MM, Cassola AC, Rocha C, Soares IS, de Braganca Pereira CA, del Portillo HA: Variant proteins of Plasmodium vivax are not clonally expressed in natural infections. Mol Microbiol. 2005, 58: 648-658. 10.1111/j.1365-2958.2005.04850.x.View ArticlePubMedGoogle Scholar
- Gabriel JA, Holmquist G, Perlmann H, Berzins K, Wigzell H, Perlmann P: Identification of a Plasmodium chabaudi antigen present in the membrane of ring stage infected erythrocytes. Mol Biochem Parasitol. 1986, 20: 67-75. 10.1016/0166-6851(86)90143-X.View ArticlePubMedGoogle Scholar
- Wanidworanun C, Barnwell JW, Shear HL: Protective antigen in the membranes of mouse erythrocytes infected with Plasmodium chabaudi. Mol Biochem Parasitol. 1987, 25: 195-201. 10.1016/0166-6851(87)90008-9.View ArticlePubMedGoogle Scholar
- Wunderlich F, Brenner HH, Helwig M: Plasmodium chabaudi malaria: protective immunization with surface membranes of infected erythrocytes. Infect Immun. 1988, 56: 3326-3328.PubMed CentralPubMedGoogle Scholar
- Wunderlich F, Helwig M, Schillinger G, Speth V, Wiser MF: Expression of the parasite protein Pc90 in plasma membranes of erythrocytes infected with Plasmodium chabaudi. Eur J Cell Biol. 1988, 47: 157-164.PubMedGoogle Scholar
- Wiser MF, Leible MB, Plitt B: Acidic phosphoproteins associated with the host erythrocyte membrane of erythrocytes infected with Plasmodium berghei and P. chabaudi. Mol Biochem Parasitol. 1988, 27: 11-21. 10.1016/0166-6851(88)90020-5.View ArticlePubMedGoogle Scholar
- Wiser MF, Giraldo LE, Schmitt-Wrede HP, Wunderlich F: Plasmodium chabaudi: immunogenicity of a highly antigenic glutamate-rich protein. Exp Parasitol. 1997, 85: 43-54. 10.1006/expr.1996.4126.View ArticlePubMedGoogle Scholar
- Giraldo LE, Grab DJ, Wiser MF: Molecular characterization of a Plasmodium chabaudi erythrocyte membrane-associated protein with glutamate-rich tandem repeats. J Eukaryot Microbiol. 1998, 45: 528-534. 10.1111/j.1550-7408.1998.tb05112.x.View ArticlePubMedGoogle Scholar
- Chen J, Trounstine M, Alt FW, Young F, Kurahara C, Loring JF, Huszar D: Immunoglobulin gene rearrangement in B cell deficient mice generated by targeted deletion of the JH locus. Int Immunol. 1993, 5: 647-656. 10.1093/intimm/5.6.647.View ArticlePubMedGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: 3-Google Scholar
- Smyth GK, Michaud J, Scott HS: Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005, 21: 2067-2075. 10.1093/bioinformatics/bti270.View ArticlePubMedGoogle Scholar
- Wettenhall JM, Smyth GK: limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics. 2004, 20: 3705-3706. 10.1093/bioinformatics/bth449.View ArticlePubMedGoogle Scholar
- R: A Language and Environment for Statistical Computing. http://www.R-project.org.
- Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.View ArticlePubMedGoogle Scholar
- Silver JD, Ritchie ME, Smyth GK: Microarray background correction: maximum likelihood estimation for the normal-exponential convolution. Biostatistics. 2009, 10: 352-363.PubMed CentralView ArticlePubMedGoogle Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Statist Soc Ser B (Methodological). 1995, 57: 289-300.Google Scholar
- Hu G, Cabrera A, Kono M, Mok S, Chaal BK, Haase S, Engelberg K, Cheemadan S, Spielmann T, Preiser PR, Gilberger TW, Bozdech Z: Transcriptional profiling of growth perturbations of the human malaria parasite Plasmodium falciparum. Nat Biotechnol. 2010, 28: 91-98. 10.1038/nbt.1597.View ArticlePubMedGoogle Scholar
- Liew KJ, Hu G, Bozdech Z, Peter PR: Defining species specific genome differences in malaria parasites. BMC Genomics. 2010, 11: 128-10.1186/1471-2164-11-128.PubMed CentralView ArticlePubMedGoogle Scholar
- Mather MW, Henry KW, Vaidya AB: Mitochondrial drug targets in apicomplexan parasites. Curr Drug Targets. 2007, 8: 49-60. 10.2174/138945007779315632.View ArticlePubMedGoogle Scholar
- Gene Expression Omnibus Database. www.ncbi.nlm.nih.gov/geo.
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.View ArticlePubMedGoogle Scholar
- Minimal information about a microarray experiment (MIAMI). http://www.mged.org/Workgroups/MIAME/miame.html.
- Bozdech Z, Mok S, Hu G, Imwong M, Jaidee A, Russell B, Ginsburg H, Nosten F, Day NP, White NJ, Carlton JM, Preiser PR: The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites. Proc Natl Acad Sci U S A. 2008, 105: 16290-16295. 10.1073/pnas.0807404105.PubMed CentralView ArticlePubMedGoogle Scholar
- Ganesan K, Ponmee N, Jiang L, Fowble JW, White J, Kamchonwongpaisan S, Yuthavong Y, Wilairat P, Rathod PK: A genetically hard-wired metabolic transcriptome in Plasmodium falciparum fails to mount protective responses to lethal antifolates. PLoS Pathog. 2008, 4: e1000214-10.1371/journal.ppat.1000214.PubMed CentralView ArticlePubMedGoogle Scholar
- Gunasekera AM, Myrick A, Le Roch K, Winzeler E, Wirth DF: Plasmodium falciparum: genome wide perturbations in transcript profiles among mixed stage cultures after chloroquine treatment. Exp Parasitol. 2007, 117: 87-92. 10.1016/j.exppara.2007.03.001.View ArticlePubMedGoogle Scholar
- De Backer MD, Ilyina T, Ma XJ, Vandoninck S, Luyten WH, Vanden Bossche H: Genomic profiling of the response of Candida albicans to itraconazole treatment using a DNA microarray. Antimicrob Agents Chemother. 2001, 45: 1660-1670. 10.1128/AAC.45.6.1660-1670.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Bammert GF, Fostel JM: Genome-wide expression patterns in Saccharomyces cerevisiae: comparison of drug treatments and genetic alterations affecting biosynthesis of ergosterol. Antimicrob Agents Chemother. 2000, 44: 1255-1265. 10.1128/AAC.44.5.1255-1265.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Boshoff HI, Myers TG, Copp BR, McNeil MR, Wilson MA, Barry CE: The transcriptional responses of Mycobacterium tuberculosis to inhibitors of metabolism: novel insights into drug mechanisms of action. J Biol Chem. 2004, 279: 40174-40184. 10.1074/jbc.M406796200.View ArticlePubMedGoogle Scholar
- Wilson M, DeRisi J, Kristensen HH, Imboden P, Rane S, Brown PO, Schoolnik GK: Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization. Proc Natl Acad Sci U S A. 1999, 96: 12833-12838. 10.1073/pnas.96.22.12833.PubMed CentralView ArticlePubMedGoogle Scholar
- Natalang O, Bischoff E, Deplaine G, Proux C, Dillies MA, Sismeiro O, Guigon G, Bonnefoy S, Patarapotikul J, Mercereau-Puijalon O, Coppee JY, David PH: Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate. BMC Genomics. 2008, 9: 388-10.1186/1471-2164-9-388.PubMed CentralView ArticlePubMedGoogle Scholar
- Oakley MS, Kumar S, Anantharaman V, Zheng H, Mahajan B, Haynes JD, Moch JK, Fairhurst R, McCutchan TF, Aravind L: Molecular factors and biochemical pathways induced by febrile temperature in intraerythrocytic Plasmodium falciparum parasites. Infect Immun. 2007, 75: 2012-2025. 10.1128/IAI.01236-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi Q, Cernetich A, Daly TM, Galvan G, Vaidya AB, Bergman LW, Burns JM: Alteration in host cell tropism limits the efficacy of immunization with a surface protein of malaria merozoites. Infect Immun. 2005, 73: 6363-6371. 10.1128/IAI.73.10.6363-6371.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Cunningham D, Fonager J, Jarra W, Carret C, Preiser P, Langhorne J: Rapid changes in transcription profiles of the Plasmodium yoelii yir multigene family in clonal populations: lack of epigenetic memory?. PLoS One. 2009, 4: e4285-10.1371/journal.pone.0004285.PubMed CentralView ArticlePubMedGoogle Scholar
- Fonager J, Cunningham D, Jarra W, Koernig S, Henneman AA, Langhorne J, Preiser P: Transcription and alternative splicing in the yir multigene family of the malaria parasite Plasmodium y. yoelii: identification of motifs suggesting epigenetic and post-transcriptional control of RNA expression. Mol Biochem Parasitol. 2007, 156: 1-11. 10.1016/j.molbiopara.2007.06.006.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.