Characterization and tissue-specific expression patterns of the Plasmodium chabaudi cir multigene family
© Ebbinghaus and Krücken; licensee BioMed Central Ltd. 2011
Received: 14 June 2011
Accepted: 19 September 2011
Published: 19 September 2011
Variant antigens expressed on the surface of parasitized red blood cells (pRBCs) are important virulence factors of malaria parasites. Whereas Plasmodium falciparum erythrocyte membrane proteins 1 (PfEMP1) are responsible for sequestration of mature parasites, little is known about putative ligands mediating cytoadherence to host receptors in other Plasmodium species. Candidates include members of the pir superfamily found in the human parasite Plasmodium vivax (vir), in the simian pathogen Plasmodium knowlesi (kir) and in the rodent malarias Plasmodium yoelii (yir), Plasmodium berghei (bir) and Plasmodium chabaudi (cir). The aim of this study was to reveal a potential involvement of cir genes in P. chabaudi sequestration.
Subfamilies of cir genes were identified by bioinformatic analyses of annotated sequence data in the Plasmodium Genome Database. In order to examine tissue-specific differences in the expression of cir mRNAs, RT-PCR with subfamily-specific primers was used. In total, 432 cDNA clones derived from six different tissues were sequenced to characterize the transcribed cir gene repertoire. To confirm differences in transcription profiles of cir genes, restriction fragment length polymorphism (RFLP) analyses were performed to compare different host tissues and to identify changes during the course of P. chabaudi infections in immunocompetent mice.
The phylogenetic analysis of annotated P. chabaudi putative CIR proteins identified two major subfamilies. Comparison of transcribed cir genes from six different tissues revealed significant differences in the frequency clones belonging to individual cir gene subgroups were obtained from different tissues. Further hints of difference in the transcription of cir genes in individual tissues were obtained by RFLP. Whereas only minimal changes in the transcription pattern of cir genes could be detected during the developmental cycle of the parasites, switching to expression of other cir genes during the course of an infection was observed around or after peak parasitemia.
The tissue-specific expression of cir mRNAs found in this study indicates correlation between expression of CIR antigens and distribution of parasites in inner organs. Together with comparable results for other members of the pir superfamily this suggests a role of cir and other pir genes in antigenic variation and sequestration of malaria parasites.
Antigenic variation is a major characteristic of malaria parasites of the genus Plasmodium leading to severe and chronic infections in a variety of vertebrates. Malaria parasites have developed strategies to evade host immune responses by expressing a large and diverse repertoire of variant proteins on the surface of parasitized red blood cells (pRBCs) [1–4]. By rapidly switching between these antigens, the parasites avoid antibody-mediated immunity of the host thus enabling the parasites to proliferate in the host without being completely eliminated by the adaptive immune response. Furthermore, these surface antigens were proposed to be involved in adherence to endothelial cells causing sequestration of late trophozoites and schizonts in post-capillary venules in specific inner organs, which considerably contributes to severe pathology in the host [5–8].
In the most virulent human pathogen Plasmodium falciparum, the PfEMP1 protein family encoded by widely studied var genes was shown to be expressed at the surface of pRBC mediating the parasite binding to endothelial cells lining small blood vessels [9, 10]. Several in vitro assays with P. falciparum revealed a tight regulation of the expression of the individual variant antigens by silencing, activation and mutually exclusive expression, thus only one variant protein is expressed in parasites at any given time [10–13]. By switching the expression to another variant, antigenic properties of the surface of pRBCs changes and the parasites prevent their complete clearance [10, 14]. Typical waves of parasitaemia in persistent P. falciparum infections reflect the repeated switching between different members of the diverse repertoire of variant antigens .
Further variant surface antigens encoded by multigene families were found in several other Plasmodium species infecting humans, monkeys and rodents [4, 15–17]. However, homologues of the var genes are not present in Plasmodium vivax or other malaria species , with exception of several var-like sequences from the chimpanzee parasite Plasmodium reichenowi, a close relative of P. falciparum.
The largest multigene family in Plasmodium genomes was described to be formed by the Plasmodium interspersed repeat (pir) genes [20, 21]. The superfamily of pir genes constitutes the major variant surface antigen family in most Plasmodium species. They were found in the human pathogen P. vivax (vir) , in the simian parasite Plasmodium knowlesi (kir)  and in the rodent malaria species Plasmodium yoelii (yir), Plasmodium berghei (bir) and Plasmodium chabaudi (cir) [15, 17, 23]. Transcriptional changes of yir genes modulated by host immunity were reported in immunocompetent mice infected with P. yoelii. In contrast to the exclusive expression of individual var genes and consecutive activation of different genes in P. falciparum, transcriptional profiling analyses in P. yoelii showed a simultaneous expression of a broad range of different yir genes within different intra-erythrocytic developmental parasite stages . In a single parasite, however, only one to three yir transcripts were detectable. Cunningham and colleagues  concluded that antigenic variation in P. yoelii probably involves exposure of the immune system to many different YIR antigens and transcriptional switching takes place without any epigenetic memory. For the variant VIR proteins of the most widely distributed human malaria pathogen P. vivax, a similar differential expression in natural infections could also be detected  and for the cir genes of the rodent malaria species P. chabaudi antigenic switching in laboratory mice has already been described .
Although the PIR proteins were localised close to the surface of pRBCs infected with P. vivax, P. yoelii, P. berghei, and P. chabaudi[20, 22, 24, 27], less is known about the role of PIR proteins in host/parasite interactions. For instance, a supposed correlation of antigenic variation of the PIRs with sequestration in inner organs has not yet been analysed. In P. falciparum, adherence of pRBCs to different host receptor such as CD36, ICAM-1 or chondroitin-sulfate A (CSA) is mediated by the major variant protein PfEMP1 [6, 8, 28], but no homologues to these proteins were detected in most other malaria pathogens. Therefore, it is conceivable that the PIR proteins, which are also expressed on the surface membrane of pRBCs, are as well involved in adherence to endothelial cells and sequestration of pRBCs in the microvasculature of inner organs.
The complex phenomenon of sequestration is hitherto not completely understood as indicated by the common assumption that sequestration of mature pRBCs in the microvasculature of the host tissues, as found in P. falciparum infections, does not occur for P. vivax. Thus, it has been assumed that the human parasite P. vivax must have developed different strategies e.g. adherence to barrier cells in the spleen to avoid spleen clearance . Recently, however, in vitro assays provided new evidence for a cytoadherence of P. vivax pRBCs to endothelial cells and placental cryosections suggesting a cytoadherence comparable to that of P. falciparum pRBCs in vivo. Direct or indirect evidences for an involvement of VIR antigens in binding to host receptors are currently still missing. In the most studied rodent model P. berghei, real-time in vivo imaging of transgenic parasites revealed CD36-mediated sequestration of schizonts in adipose tissues and lung as well as an accumulation of schizonts in the spleen . Accumulation of other blood stages was also observed in different tissues including the brain and placenta of pregnant mice [32–34]. Nevertheless, identification of parasite ligands involved in binding to CD36 or other host receptors are still missing. In vitro adherence of P. chabaudi infected erythrocytes to purified human CD36 has been observed . Moreover, sequestration to microvascular endothelial cells was reported to appear in an organ specific manner in P. chabaudi infections, predominantly in liver, but also in brain and spleen thus resembling at least partially the sequestration pattern in P. falciparum infections.
In this study, the relationship between antigenic variation of the cir multigene family and accumulation of parasites localization in different host tissues was analysed at the molecular level. Bioinformatic characterization of the cir genes was used to identify different subfamilies within the multigene family. First indications for a tissue-specific expression of cir genes were obtained by amplification of a representative sample of transcribed cir genes from different host tissues. Transcriptional profiling analyses of cir genes in different host tissues using RFLP analysis of RT-PCR products provided further evidence of a tissue-specific expression of cir genes in P. chabaudi infections. These results provide first indications for a possible correlation of antigenic variation of the cir multigene family and sequestration in host tissues in P. chabaudi infections.
All annotated sequence data of putative cir genes in the Plasmodium Genome Database Resource version 8.0  (PlasmoDB; plasmodb.org) were compiled (Additional file 1). For confirmation of correct annotation, a CD-BLAST [37–39] was performed to identify and locate the conserved domain of the CIR-BIR-YIR superfamily [22, 23] (Pfam protein families database accession number [PF06022]). Sequences with a partial CIR-BIR-YIR conserved domain were excluded from all subsequent analyses.
For phylogenetic reconstruction, putative CIR proteins were aligned using ClustalW2  and a phylogenetic tree was calculated with PhyML 3.0  using maximum likelihood estimation and the JTT model  for amino acid substitution. The gamma shape parameter and the proportion of invariable sites were estimated and the number of substitution rate categories was set to four. The implemented BIONJ algorithm was used to build the starting tree. Resulting trees (Newick format) were visualised and processed with MEGA4 [43, 44].
For validation of structural protein motifs such as transmembrane domains, signal peptides and PEXEL motifs, analyses of the protein sequences were examined with Protscale , TMHMM 2.0 [46, 47] and SignalP . Results were compared with annotations extracted from PlasmoDB.
All experiments were performed with 5-8 weeks-old outbred female NMRI mice (Crl:NMRI(Han)) provided by Charles River (Sulzfeld, Germany). The mice were kept in cages with a maximum of five animals per cage and received food and water ad libitum. The experiments were planed according to all relevant guidelines for animal protection and approved by German authorities responsible for animal protection.
Infection of mice
A non-clonal line of P. chabaudi very similar but not identical to P. chabaudi AS was used [49–51]. Blood stages of P. chabaudi were weekly passaged in female NMRI mice by intraperitoneal injection (i.p.) of a droplet of tail vein blood diluted in PBS. Parasitaemia was evaluated in Giemsa-stained blood smears and total erythrocytes number was counted in a Neubauer chamber.
For the experiments, blood of an infected NMRI mouse was collected by cardiac puncture under anesthesia. For each transcriptional profiling experiment, six mice were infected i.p. with approximately 100 parasitized red blood cells (pRBCs). Organs and blood were collected at about 30% parasitaemia, i.e. just before peak parasitaemia.
For transcriptional profiling during the course of an infection, tail vein blood of mice infected i.p. with 100 pRBCs was passaged i.p. into naïve female NMRI mice at days 7 (early infection), 14 (around peak parasitaemia), 21 (7 days after peak parasitaemia) and 35 (21 days after peak parasitaemia). Expression of cir genes was analysed when parasitaemia reached about 30%.
RNA and DNA extraction
Blood samples of P. chabaudi infected mice were collected by cardiac puncture, rapidly frozen and stored at -80°C. Small pieces of liver, spleen, kidneys, lung and brain were transferred into RNA Later (Sigma Aldrich) and kept at -80°C for long-term storage. Total RNA was extracted using NucleoSpin® RNA II kit (Macherey-Nagel) according to the manufactures instructions. Genomic DNA extraction was performed with the NucleoSpin® Blood kit (Macherey-Nagel).
Verification of complete cir gene structures
To verify the complete gene structure of selected cir genes, RT-PCR and genomic PCR were performed in parallel. The primer pairs used for amplification and detailed information about PCR conditions can be found in Additional file 2.
Cloning of RT-PCR products and sequencing
Residual contaminating genomic DNA in total RNA preparations was digested with DNase I (Fermentas). First strand cDNA was synthesised using 1 μg RNA and the RevertAid™ Premium Reverse Transcriptase (Fermentas) with Oligo dT18 primers including reactions without reverse transcriptase as negative controls for amplification.
For transcription analysis, specific primers spanning the second exon encoding an essential part of the CIR-BIR-YIR domain (> 98% of its length) were designed for both major cir gene subfamilies. The primer pairs f1up (5'-AATACGCTATTTTATGGTTTAGTTATAAA-3')/f1lo (5'-TGAAATTCCTAAAATAATGGGTATTATTAAAA-3') and f2up (5'-TATGCTATTTTATGGTTAAGTTATATGCTA-3')/f2lo (5'-ATGAATACTTATAAGCAATTCCCAAGAAAA-3') were targeted to highly conserved sequence regions for amplification of a broad range of cir genes.
Following cDNA synthesis, amplification with the AccuPrime™ DNA polymerase (Invitrogen) was performed using subfamily-specific primers for cir subfamily 1 and 2. After an initial denaturation for 30 s at 94°C, 40 cycles of 10 s at 94°C, 30 s at 55°C, 30 s at 72°C followed by a final extension of 10 min at 72°C. RT-PCR products were then isolated from a 0.8% agarose gel and precipitated in the presence of glycogen.
In order to sample a first repertoire of transcribed P. chabaudi cir genes, RT-PCR products of each tissue (blood, liver, spleen, kidney, lung and brain) were gel-purified and cloned into the pCR™4-TOPO® TA vector (Invitrogen). Thirty-six clones for each tissue and subfamily were sequenced (GATC, Constance, Germany) resulting in sequences for 216 independent cDNA clones for cir subfamily 1 and 2.
Transcriptional profiling of RT-PCR products by RFLP
For transcriptional profiling, DNase digestion of RNA and cDNA synthesis was performed as described before and RT-PCR products were amplified with the Phusion® Hot Start II High-Fidelity DNA Polymerase (Fermentas) using subfamily-specific primers for cir subfamily 1 and 2. After an initial denaturation for 30 s at 98°C, 50 cycles of 10 s at 98°C, 30 s at 55°C, 30 s at 72°C were performed followed by a final extension for 10 min at 72°C.
Gel-purified RT-PCR products (150 ng) were digested with Alu I or Xap I (Fermentas) in 5 μl for 3 h at 37°C followed by enzyme inactivation for 20 min at 65°C. The restricted fragments (30 ng) were analysed with the Agilent 2100 Bioanalyzer using the DNA 1000 LabChip® kit (Agilent) following the manufactures instructions. In order to ensure reproducibility of the RFLP profiles, both RT-PCR and RFPL were usually repeated at least twice resulting in 4 replicates.
In order to evaluate whether clones representing certain cir subfamilies were significantly more frequently recovered from one tissue than from others, frequencies were compared using a Z-test.
Phylogenetic analyses of putative CIR proteins in PlasmoDB
For identification of the individual proteins in each branch, the phylogenetic tree of CIR and YIR proteins is given in Newick format in Additional file 3.
Localization motifs found in CIR proteins
For protein structural analyses, the 186 putative CIR proteins were examined in regard to subcellular localization characteristics, such as putative signal peptides, transmembrane domains, and the PEXEL motif (Plasmodium export element) [53, 54]. Results were then compared to annotations in PlasmoDB. Detailed information on protein motifs identified in all putative CIR proteins can be obtained from Additional file 1. In only one of the 186 sequences (PCHAS_114740) a signal peptide had been annotated, however, the re-analysis with SignalP revealed only a 1.6% probability for the presence of a signal peptide in this particular predicted CIR protein with a maximum cleavage site probability of 1.2% between amino acid position 26 and 27. Moreover, a cDNA sequence encoding a corresponding CIR protein with signal peptide could not be amplified by RT-PCR using primers derived from PCHAS_114740 in the non-clonal P. chabaudi line used throughout this study. Altogether, this indicates that the CIR proteins lack a signal peptide and presumably use alternative pathways for transport out of the parasite cell across the parasitophorous vacuole to the host cell membrane.
Finally, the putative CIR proteins were scanned for the conserved five-residue PEXEL motif Rx[L/I]x[E/Q/D] with the Protein Motif Pattern search tool in PlasmoDB. This short hydrophobic peptide has been proposed to be localised 16-24 amino acids downstream of the NH2-terminus and to mediate the export of proteins across the parasitophorous vacuole to the erythrocyte cytoplasm [53, 54]. Seventeen sequences of the putative CIR proteins were identified to have the specific PEXEL motif but none of these sequences exhibited the PEXEL motif at its canonical position. In the protein sequence PCHAS_114700 the motif is located at amino acid position 4-8 and in four sequences between position 30-70. In another 10 small CIR proteins and in two large CIR proteins, the PEXEL motif is found between the amino acids 90-300 or 640-660, much closer to the COOH- than to the NH2-terminus. Hence, all annotated putative CIR proteins only possess an Rx[L/I]x[E/Q/D] stretch in their sequence but lack a functioning PEXEL motif.
Verification of complete gene structures
RT-PCR of cir genes with subfamily-specific primers
For the subsequent transcriptional profiling analyses (see below), subfamily-specific primers for cir subfamily 1 and subfamily 2 were designed spanning the second exon, the essential part of the CIR-BIR-YIR conserved domain. The primer pairs were located at highly conserved sequence regions for amplification of a broad range of cir cDNAs (Figure 4A). The expected RT-PCR product sizes for subfamily 1 and subfamily 2 were about 600 bp (Figure 4B). In most PCRs, however, an additional smaller but very faint RT-PCR product (approx. 300 bp) appeared. Sequencing analyses indicated that this product is a notably smaller cir transcript with a deletion of about 250 bp within the second exon. This cDNA most likely also represents an alternatively spliced cir variant, in particular since this cDNA [GenBank:JF969288] has an uninterrupted open reading frame.
Cloning of cir cDNAs from individual tissues
To examine whether the expression of certain cir antigen mRNAs correlates with the localization of the parasites in particular internal organs, a first repertoire of transcribed cir genes was amplified from parasites using RT-PCR and subfamily-specific primers for subfamilies 1 and 2. For this purpose, a single female NMRI mouse was i.p. infected with 100 pRBCs and, at a parasitaemia of 30%, blood, liver, spleen, kidney, lung and brain were collected for RNA isolation. RT-PCR products of about 600 bp from each tissue were excised from agarose gels and directly cloned into the pCR4-TOPO vector. For each tissue, clones were obtained from at least three independent RT-PCR reactions. Thirty-six clones for each tissue and subfamily were sequenced resulting in 216 cir sequences for each cir subfamily. Information about individual sequences with GenBank® accession numbers can be found in Additional file 5. A broad range of different sequences from both subfamilies was amplified as shown by phylogenetic maximum likelihood analysis of the deduced protein sequences in comparison with the 186 putative conserved domains obtained from PlasmoDB (Additional file 6) indicating that the method is able to amplify a representative repertoire of both subfamilies. The facts that (i) the length of the amplified fragments differs only minimally (mean length and SD is 561 bp ± 2.5% for subfamily 1 and 523.7 bp ± 4.4% for subfamily 2) and (ii) the percent identity is very high (mean percent identity and SD is 80.8% ± 8.5% for subfamily 1 and 76.7% ± 10.9% for subfamily 2) a strong bias introduced by unequal amplification efficacies is at least unlikely.
Statistical analyses using a Z-test were done to identify significant differences in the frequency of different tissues represented within those individual subfamilies that are sufficiently represented within the total population. The method does not allow to detect expression differences for minor cir transcripts. Results are also indicated in Figure 5. Most remarkably, the cir cDNA belonging to subgroup 1.1 could be frequently recovered from liver (24.2% of all subgroup 1.1 clones), lung, (24.2%) blood (21.4%), spleen (17.1%) and kidney (12.9%) but not from brain (0%). The absence of this cir cDNA from our brain samples is statistically significant when compared to all other tissues (p < 0.01). In contrast, the cir cDNA sequences of subgroup 1.4 were not found in blood cDNA but significantly more often present in brain (34.6%; p < 0.01) and spleen (23.1%; p < 0.05) samples. Furthermore, a significant difference could be detected between brain and lung (7.7%) in subgroup 1.4 (p < 0.05). The cir cDNAs in subgroup 1.5 were absent from liver samples (0%) but a higher number (60%) of these clones came from brain samples than from all other tissues investigated (10% for each tissue). While the differences between brain and liver were statistically significant (p < 0.05), differences between brain and the other tissues (10%) did just not reach significance (p < 0.07). The cDNA sequences of subgroup 1.6 were absent from spleen, kidney and brain samples but were found in blood (45.5%), liver (27.3%) and lung (27.3%). The frequency of clones obtained from blood within the cir subgroup 1.6 were significantly higher when compared to spleen, kidney and brain (p < 0.01). Finally, for subgroup 1.8 no cir cDNAs could be detected in liver and kidney. In contrast, this subgroup was obtained with a high frequency of 43.8% from the spleen and with intermediate frequencies of 12.5% from blood, 18.8% from lung and 25.0% from brain. However, only the differences between the frequencies for spleen, liver and kidney were statistically significant (p < 0.05).
All these differences in tissue distribution of the transcribed cir cDNAs found in this cloning and sequencing study are first indications for a tissues specific expression of cir genes and were analysed in more detail in the subsequent experiments.
Expression profiling of cir genes in different host tissues
Due to these first differences in the expression of cir genes, further analyses of different transcription profiles in various host tissues during the infection were carried out using RFLP of RT-PCR products in order to confirm differences in expression of cir genes between different host tissues for a larger number of mice.
Comparing the RFLP profiles, considerable differences in the banding patterns of the individual organs for subfamily 1 and 2 were observed thus confirming differences in expression of cir mRNAs between different host tissues. Nevertheless, the prominent fragments observed in RFLP agreed in many, though not all cases well with the restriction pattern obtained in silico from the sequences of the cir clones described above (see Additional file 9). This similarity of experimental RFLP patterns and the in silico predicted restriction patterns of some cir cDNA sequences of the cloning and sequencing experiment suggested at least a similar repertoire of amplified cir genes in both infection experiments. On the other hand, none of the profile pattern was identical between the four mice infected with 100 P. chabaudi pRBCs taken from the same parental population. This indicates that random selection of 100 pRBCs as starter population from a large parental population severely influences the initial cir gene expression pattern. Despite the fact that patterns for organs differed within individual mice, there was no clear tissue-specific pattern detectable that was conserved between mice.
For subfamily 1, no transcriptional differences between the tissues were found in mouse A (Figure 7A). Moreover, the restriction pattern obtained for this mouse was highly different in comparison to the restriction profiles of mice B-D and the sizes of the individual fragments did not clearly refer to any in silico fragments of cir sequences obtained in the cloning study. This presence or absence of individual fragments in the RT-PCR RFLP profiles for subfamily 1 are further hints for a differential expression of cir genes in different host tissues. In particular, in mouse B, additional fragments in restriction patterns were observed e.g. in brain (~ 370 bp) or in kidney and lung (~ 340 bp) which were absent from samples of all other tissues (Figure 7A). Further differences within the transcriptional profiles of cir subfamily 1 were detected e.g. in the spleen of mouse C, showing a distinct banding pattern compared to the other five tissues of the same mouse.
Differences - though generally less pronounced - were also found in cir subfamily 2 (Figure 7B). An additional fragment of ~ 400 bp was present e.g. in the liver of mouse B which was absent in blood, spleen kidney, lung, and brain samples from the same mouse. Furthermore, two additional bands of cir fragments of approximately 100-200 bp appeared in the spleen of mouse B.
In conclusion, the differences in the transcriptional profiles of the different host tissues found in this experiment confirmed a correlation between parasite tissue localization on expression of specific cir genes.
Expression profiling of cir genes during the course of infection
Surprisingly, only two out of four immunocompetent mice (mouse 2 and 4) showed changes in expression of subfamily 1 cir genes in the blood between day 7 (early infection), day 14 (around peak parasitaemia) and day 21 (after peak parasitaemia) (Figure 8A) p.i. However, the virtually identical patterns produced here provide additional evidence for the reproducibility of the RFLP method. In mouse 2, a completely different restriction fragment pattern was observed at day 21 p.i., seven days after peak parasitaemia, compared to the profiles of early parasitaemia (day 7 p.i.) and peak parasitaemia (day 14 p.i.). In mouse 4 the first differences in the expression pattern of the cir genes were observed already at peak parasitaemia (day 14 p.i.) and a completely changed RFLP pattern was found at day 21 p.i. (after parasitaemia) compared to day 7 p.i. in the early phase of infection. Surprisingly, however, no changes in the expression pattern could be observed for subfamily 1 in mouse 1 and 3 (Figure 8A) and for all four mice in subfamily 2. The same RFLP pattern of amplified cir cDNAs could be identified throughout the course of an infection (Figure 8B). Since Alu I did apparently not cut the cir subfamily 2 members expressed in this experiment, an additional restriction digestion with Xap I was then performed (Figure 8B) confirming further that no changes in expression of cir subfamily 2 occurred in the course of an infection in these mice.
At day 35 p.i., all mice had completely resolved the infection irrespectively of changes in cir gene expression pattern before. Apparently, a chronic infection of the mice failed to occur after i.p. infection with 100 pRBC whether or not changes in cir gene expression pattern were detectable in the earlier infection.
Transcriptional profiling of cir genes throughout intraerythrocytic development
For RFLP analyses of transcriptional changes of the cir genes at different life cycle stages, 30 μl tail vein blood of female NMRI mice infected with 100 pRBCs was collected at 3 h (early trophozoites), 10 h (late trophozoites) and 17 h (mature trophozoites and early schizonts) after beginning of the light cycle on day 13 p.i. (parasitaemia about 30%).
In order to improve understanding of host/parasite interactions correlating antigenic variation and accumulation of parasite localization in different host tissues, expression of the P. chabaudi cir multigene family was analysed in vivo at the transcriptional level. Despite diverse in vitro or in vivo investigations regarding antigenic variation in Plasmodium species, the complete molecular mechanism of antigenic switching in the parasites is, hitherto, far away from being fully understood. The present study, therefore, examined P. chabaudi infections in immunocompetent mice to get further insight in the complex phenomenon of antigenic variation and sequestration.
In an initial phylogenetic analysis of the annotated repertoire of putative CIR proteins, two subfamilies and an unassigned group of very long CIRs could be identified. While the subfamily 1 and 2 exhibit the predicted primary structure of common PIR proteins, the unassigned CIR proteins show more divergence often due to large insertions within the conserved CIR-BIR-YIR domain.
It has already previously been proposed that different subsets within the pir superfamily exhibit different functions. In P. yoelii and P. berghei, for example, a stage specific role was assumed since not all pir genes are transcribed equally when different life cycle stages were compared [24, 55]. Indeed, different roles for P. chabaudi cir subfamily 1 (with changes in expression pattern) and subfamily 2 (with minor or no changes in expression pattern) genes are likely. The large insertions and additional transmembrane domains found within the group of unassigned CIR proteins might also be interpreted as hints for distinct functions of these proteins. Since none of them has been examined in any functional or even localization study, there is still no clue for the role these proteins might play in the Plasmodium life cycle or in the pathogens immune evasion strategy. For the VIR proteins in P. vivax diverse protein domain and structure predictions have been suggested to indicate different functions for different subsets of these proteins as for example those with PEXEL motif or those with additional transmembrane domains . However, the presence or absence of the PEXEL motifs might simply result in usage of alternative transport signals or pathways to the host cell membrane  and currently no functional role has been suggested for any of the PIR proteins with more than one transmembrane domain. Since bioinformatic analysis of the CIR proteins revealed not a single strong PEXEL motif, it can be assumed that PEXEL motifs are not an important feature in the function of CIR proteins and that at least the majority of PIR proteins is able to be transported to the pRBC surface without an obvious PEXEL translocation signal.
A common ancestry of the gene families forming the pir superfamily has been postulated due to e.g. conserved sequence motifs as well as structural predictions within the PIR sequences . By phylogenetic comparison of CIR and YIR proteins as performed here, this hypothesis can be further strengthened. The phylogenetic relationship of the CIR proteins with the widely studied YIRs in P. yoelii shows that the YIR protein sequences only share high similarity to CIR proteins of subfamily 2 and neither to subfamily 1 nor the group of unassigned CIR proteins. Comparison to the YIR antigens of the close relative P. yoelii suggests that most of the CIR sequences in the P. chabaudi genome have evolved by diversification after separation of P. yoelii and P. chabaudi from a common ancestor. Such a dispersing evolution of distinctly evolved subfamilies has also been demonstrated by phylogenetic analysis of PIR proteins encoded in P. yoelii and P. berghei genomes  and it is likely that the high variability of antigens in the individual malaria species distinctly evolved probably in response of the host immune pressure.
In this study a non-clonal P. chabaudi line very similar to the clonal P. chabaudi AS strain was used. The deduced protein sequences of the amplified cir transcripts of these parasites, however, were quite similar but not identical to the annotated putative CIR predicted from the genome sequence of the clonal P. chabaudi AS strain supposing a high variability of antigens even between closely related strains within the same species.
Cloning and sequencing of a first subset of cir genes demonstrated that a broad range of subfamily 1 and subfamily 2 cir genes is transcribed during a P. chabaudi infection in immunocompetent mice infected with a starter population of a minimal size. These findings are consistent with those described previously for other pir multigene families, e.g. the vir genes in P. vivax, the yir genes in P. yoelii and the bir genes in P. berghei[25, 26, 55]. In contrast to the mutually exclusive expression of only one var gene found in P. falciparum parasites cultured in vitro, it has been shown that in P. vivax and in several rodents Plasmodium species many different pir genes were transcribed in vivo within a parasite population in an individual host. Interestingly, examinations of natural P. falciparum infected human samples has shown that - in contrast to the limited transcription pattern in cultured P. falciparum parasites - many transcripts of var and stevor genes are also transcribed simultaneously in vivo[59, 60]. Such contradictory in vivo and in vitro findings make clear that in vivo models such as P. chabaudi or other rodent malaria parasites are essential for complex investigations of antigenic variation.
In accordance with previous studies of the yir genes of P. yoelii, switching of cir gene expression could be detected around peak parasitaemia in the course of infection suggesting that antigenic variation might be modulated by selective forces exercised by the host immune system. However, these transcriptional changes were not observed in all infected mice and, moreover, were not detectable to the same extent for all cir subfamilies. Most marked differences in mRNA expression patterns could be observed for cir subfamily 1 whereas no or obviously less transcriptional switching was detected for cir subfamily 2 during the infection.
Completely different cir gene expression patterns of the progenies derived from 100 pRBCs starter populations originating from the same parental population strongly suggest that a large parasite population can be extremely heterogeneous. In micromanipulated P. yoelii pRBCs, it has been shown that only one to three different yir transcripts were transcribed within an individual cell but many different transcripts were detected within a whole parasite population . In addition, rapid switching in the transcribed repertoires of yir genes between different clonal host parasites populations and parasite developmental stages has been described. Although no infections with single pRBCs were performed in the present study, the initial diversity of parasites was apparently kept at a minimum using the minimal infectious dose resulting in patent infections by intraperitoneal infection as revealed by the substantial differences in RT-PCR RFLP patterns between mice infected with 100 pRBCs. Therefore, the number of transcribed cir genes per pRBC is presumably also much lower than that found to be transcribed in a large population.
Previous studies have analysed switching between different yir genes in P. yoelii infected mice both in immunocompetent  and highly immunodeficient  mice. Efficient and frequent switching of expressed yir genes could be observed even in the absence of any selecting force of an adaptive immune system . The fact that cir gene subfamily 1 and subfamily 2 expression patterns apparently do not necessarily change in the course of a P. chabaudi infection in immunocompetent mice even within three weeks suggests lower on-off switching frequencies at least for these groups of cir genes. This is particularly surprising as adaptive immune responses are well known to effectively select for parasites expressing new variant antigens in other protozoan parasites [61–63].
Changes in the expression patterns of variant antigens appear to occur more frequently in cir subfamily 1 than in subfamily 2. No data are currently available for switching frequencies within the highly heterogeneous group of unassigned cir genes. Fonager et al already speculated that within the yir gene family those genes belonging to highly heterogeneous groups would be more important for antigenic variation than members of more conventional subfamilies. The phylogenetic analysis including members of all previously defined yir groups  shows that even the most divergent YIR proteins cluster together with the CIR subfamily 2. This suggests that subfamily 2 is ancient among the CIR proteins while subfamily 1 and the unassigned CIR proteins have probably evolved from ancient genes by diversification. The higher apparent switching frequency observed for cir subfamily 1 than subfamily 2 genes is a first experimental hint corroborating the hypothesis that those pir family members that evolved relatively recently might play a more prominent role for antigenic variation than those showing ancient properties . Future work in both P. yoelii and P. chabaudi should, therefore, no longer neglect these unusual pir members from the analyses but should pay special attention to them.
Analysis of changes in cir gene expression patterns in the course of one asexual round of multiplication in 24 h revealed only minor changes in gene expression for subfamily 1 and virtually no changes for subfamily 2. In particular, no reduction in the complexity of the cir genes in late parasite stages was observed in comparison to early rings as has been described for var genes . Since Cunningham et al found evidence that in P. yoelii even schizonts transcribe at least up to three different yir genes, a mutually exclusive expression of pir genes in late individual parasites or parasite populations appears to be highly unlikely. The same was also shown for vir genes in P. vivax with more than one antigen expressed in a single parasite and different expression patterns between parasites . Transcriptomic analysis of P. vivax intraerythrocytic developmental cycle stages also revealed that many vir genes are expressed early during the ring stage and turned off later while others are expressed predominantly in schizonts . These authors could not find any linkage between the position of the vir gene within the phylogenetic tree and its predominant expression time. In contrast to these results, only minor changes during the intraerythrocytic cycle were found in the present study. Possible explanations for this discrepancy include the fact that of course the RFLP analysis is less sensitive with regard to the detection of changes in expression of individual genes when compared to the microarray method used by Bozdech et al. In addition, only for about 60% of the vir genes consistent temporal expression patterns could be observed for three different P. vivax isolates. If the highly expressed genes show no temporal expression pattern, such a pattern would clearly not be detectable with the RFLP method used here since it is not able to detect minor cir transcripts at all. Finally, the different experimental designs with only 100 pRBCs as founders in the cir gene study and non-selected parasites from naturally infected patients for the P. vivax transcriptome study might well contribute to the observed differences.
Remarkable differences in the mRNA expression patterns between pRBCs in blood, liver, spleen, kidney, lung and brain could be observed suggesting a host-tissue specific expression of cir genes. It is believed that the PIR proteins, like the Pfemp1 protein family in P. falciparum, are possibly involved in adhesion to host receptors thus mediating accumulation and sequestration in different host tissues. For the BIR or YIR proteins, for example, an expression of these molecules on the surface of pRBCs has already been demonstrated [24, 27] but neither direct nor indirect evidence for a correlation of an adhesion of PIR proteins to host endothelial cells has yet been found. Accumulation of pRBCs in different tissues expressing different cir genes, however, can be considered to be a first experimental hint that CIR and maybe also other PIR proteins might indeed be involved in adhesion and sequestration. Whether rapid and dramatic changes in tissue distribution of P. chabaudi parasites at peak parasitaemia, i.e. exclusion from the splenic red pulp , has effects on cir gene expression patterns would also be interesting to analyse in the future. Since P. chabaudi is the only frequently used experimental malaria model with synchronous development and a strong sequestration phenotype, these results suggest that further studies of CIR proteins will provide important new insights into the interaction of non-P. falciparum malaria pRBC with host epithelia.
The present study has demonstrated for the first time that there are differences in the tissue-specific expression of some cir genes. These results suggest a possible correlation between the expression of CIR antigens and accumulation of parasites in inner organs of the host. The high agreement in results obtained for transcriptional switching and antigenic variation of the cir genes with that for other members of the pir superfamily indicate very well that the cir genes are promising candidates for further functional studies. For example parasites expressing certain CIRs constitutively as transgenes might improve understanding of the functional role of the pir superfamily in Plasmodium infections.
This study was supported by a research grant from the Deutsche Forschungsgemeinschaft (DFG) (KR 2245/5-1). The authors thank Georg von Samson-Himmelstjerna and Janina Demeler for support and critical revision of the manuscript.
- Brown KN, Brown IN: Immunity to malaria: antigenic variation in chronic infections of Plasmodium knowlesi. Nature. 1965, 208: 1286-1288. 10.1038/2081286a0.View ArticlePubMedGoogle Scholar
- McLean SA, Pearson CD, Phillips RS: Plasmodium chabaudi: antigenic variation during recrudescent parasitaemias in mice. Exp Parasitol. 1982, 54: 296-302. 10.1016/0014-4894(82)90038-8.View ArticlePubMedGoogle Scholar
- Barnwell JW, Howard RJ, Miller LH: Altered expression of Plasmodium knowlesi variant antigen on the erythrocyte membrane in splenectomized rhesus monkeys. J Immunol. 1982, 128: 224-226.PubMedGoogle Scholar
- Mendis KN, Ihalamulla RI, David PH: Diversity of Plasmodium vivax-induced antigens on the surface of infected human erythrocytes. Am J Trop Med Hyg. 1988, 38: 42-46.PubMedGoogle Scholar
- Franke-Fayard B, Fonager J, Braks A, Khan SM, Janse CJ: Sequestration and tissue accumulation of human malaria parasites: can we learn anything from rodent models of malaria?. PLoS Pathog. 2010, 6:Google Scholar
- Miller LH, Baruch DI, Marsh K, Doumbo OK: The pathogenic basis of malaria. Nature. 2002, 415: 673-679. 10.1038/415673a.View ArticlePubMedGoogle Scholar
- Mackintosh CL, Beeson JG, Marsh K: Clinical features and pathogenesis of severe malaria. Trends Parasitol. 2004, 20: 597-603. 10.1016/j.pt.2004.09.006.View ArticlePubMedGoogle Scholar
- Rasti N, Wahlgren M, Chen Q: Molecular aspects of malaria pathogenesis. FEMS Immunol Med Microbiol. 2004, 41: 9-26. 10.1016/j.femsim.2004.01.010.View ArticlePubMedGoogle Scholar
- Su XZ, Heatwole VM, Wertheimer SP, Guinet F, Herrfeldt JA, Peterson DS, Ravetch JA, Wellems TE: The large diverse gene family var encodes proteins involved in cytoadherence and antigenic variation of Plasmodium falciparum-infected erythrocytes. Cell. 1995, 82: 89-100. 10.1016/0092-8674(95)90055-1.View ArticlePubMedGoogle Scholar
- Dzikowski R, Templeton TJ, Deitsch K: Variant antigen gene expression in malaria. Cell Microbiol. 2006, 8: 1371-1381. 10.1111/j.1462-5822.2006.00760.x.View ArticlePubMedGoogle Scholar
- Deitsch KW, Calderwood MS, Wellems TE: Malaria. Cooperative silencing elements in var genes. Nature. 2001, 412: 875-876.View ArticlePubMedGoogle Scholar
- Voss TS, Healer J, Marty AJ, Duffy MF, Thompson JK, Beeson JG, Reeder JC, Crabb BS, Cowman AF: A var gene promoter controls allelic exclusion of virulence genes in Plasmodium falciparum malaria. Nature. 2006, 439: 1004-1008.PubMedGoogle Scholar
- Recker M, Buckee CO, Serazin A, Kyes S, Pinches R, Christodoulou Z, Springer AL, Gupta S, Newbold CI: Antigenic variation in Plasmodium falciparum malaria involves a highly structured switching pattern. PLoS Pathog. 2011, 7: e1001306-10.1371/journal.ppat.1001306.PubMed CentralView ArticlePubMedGoogle Scholar
- Deitsch KW, Moxon ER, Wellems TE: Shared themes of antigenic variation and virulence in bacterial, protozoal, and fungal infections. Microbiol Mol Biol Rev. 1997, 61: 281-293.PubMed CentralPubMedGoogle Scholar
- Fischer K, Chavchich M, Huestis R, Wilson DW, Kemp DJ, Saul A: Ten families of variant genes encoded in subtelomeric regions of multiple chromosomes of Plasmodium chabaudi, a malaria species that undergoes antigenic variation in the laboratory mouse. Mol Microbiol. 2003, 48: 1209-1223. 10.1046/j.1365-2958.2003.03491.x.View ArticlePubMedGoogle Scholar
- al-Khedery B, Barnwell JW, Galinski MR: Antigenic variation in malaria: a 3' genomic alteration associated with the expression of a P. knowlesi variant antigen. Mol Cell. 1999, 3: 131-141. 10.1016/S1097-2765(00)80304-4.View 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 parasite Plasmodium yoelii yoelii. Nature. 2002, 419: 512-519. 10.1038/nature01099.View ArticlePubMedGoogle Scholar
- Korir CC, Galinski MR: Proteomic studies of Plasmodium knowlesi SICA variant antigens demonstrate their relationship with P. falciparum EMP1. Infect Genet Evol. 2006, 6: 75-79. 10.1016/j.meegid.2005.01.003.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
- 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
- Cunningham D, Lawton J, Jarra W, Preiser P, Langhorne J: The pir multigene family of Plasmodium: antigenic variation and beyond. Mol Biochem Parasitol. 2010, 170: 65-73. 10.1016/j.molbiopara.2009.12.010.View ArticlePubMedGoogle Scholar
- del Portillo HA, Fernandez-Becerra C, Bowman S, Oliver K, Preuss M, Sanchez CP, Schneider NK, Villalobos JM, Rajandream MA, Harris D, Pereira da Silva LH, Barrell B, Lanzer M: A superfamily of variant genes encoded in the subtelomeric region of Plasmodium vivax. Nature. 2001, 410: 839-842. 10.1038/35071118.View 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
- 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
- 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
- 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
- Di Girolamo F, Raggi C, Birago C, Pizzi E, Lalle M, Picci L, Pace T, Bachi A, de Jong J, Janse CJ, Waters AP, Sargiacomo M, Ponzi M: Plasmodium lipid rafts contain proteins implicated in vesicular trafficking and signalling as well as members of the PIR superfamily, potentially implicated in host immune system interactions. Proteomics. 2008, 8: 2500-2513. 10.1002/pmic.200700763.View ArticlePubMedGoogle Scholar
- Sherman IW, Eda S, Winograd E: Cytoadherence and sequestration in Plasmodium falciparum: defining the ties that bind. Microbes Infect. 2003, 5: 897-909. 10.1016/S1286-4579(03)00162-X.View ArticlePubMedGoogle Scholar
- del Portillo HA, Lanzer M, Rodriguez-Malaga S, Zavala F, Fernandez-Becerra C: Variant genes and the spleen in Plasmodium vivax malaria. Int J Parasitol. 2004, 34: 1547-1554. 10.1016/j.ijpara.2004.10.012.View ArticlePubMedGoogle Scholar
- Carvalho BO, Lopes SC, Nogueira PA, Orlandi PP, Bargieri DY, Blanco YC, Mamoni R, Leite JA, Rodrigues MM, Soares IS, Oliveira TR, Wunderlich G, Lacerda MV, del Portillo HA, Araujo MO, Russell B, Suwanarusk R, Snounou G, Renia L, Costa FT: On the cytoadhesion of Plasmodium vivax-infected erythrocytes. J Infect Dis. 2010, 202: 638-647. 10.1086/654815.View ArticlePubMedGoogle Scholar
- Franke-Fayard B, Janse CJ, Cunha-Rodrigues M, Ramesar J, Buscher P, Que I, Lowik C, Voshol PJ, den Boer MA, van Duinen SG, Febbraio M, Mota MM, Waters AP: Murine malaria parasite sequestration: CD36 is the major receptor, but cerebral pathology is unlinked to sequestration. Proc Natl Acad Sci USA. 2005, 102: 11468-11473. 10.1073/pnas.0503386102.PubMed CentralView ArticlePubMedGoogle Scholar
- Spaccapelo R, Janse CJ, Caterbi S, Franke-Fayard B, Bonilla JA, Syphard LM, Di Cristina M, Dottorini T, Savarino A, Cassone A, Bistoni F, Waters AP, Dame JB, Crisanti A: Plasmepsin 4-deficient Plasmodium berghei are virulence attenuated and induce protective immunity against experimental malaria. Am J Pathol. 2010, 176: 205-217. 10.2353/ajpath.2010.090504.PubMed CentralView ArticlePubMedGoogle Scholar
- Neres R, Marinho CR, Goncalves LA, Catarino MB, Penha-Goncalves C: Pregnancy outcome and placenta pathology in Plasmodium berghei ANKA infected mice reproduce the pathogenesis of severe malaria in pregnant women. PLoS One. 2008, 3: e1608-10.1371/journal.pone.0001608.PubMed CentralView ArticlePubMedGoogle Scholar
- Marinho CR, Neres R, Epiphanio S, Goncalves LA, Catarino MB, Penha-Goncalves C: Recrudescent Plasmodium berghei from pregnant mice displays enhanced binding to the placenta and induces protection in multigravida. PLoS One. 2009, 4: e5630-10.1371/journal.pone.0005630.PubMed CentralView ArticlePubMedGoogle Scholar
- Mota MM, Jarra W, Hirst E, Patnaik PK, Holder AA: Plasmodium chabaudi-infected erythrocytes adhere to CD36 and bind to microvascular endothelial cells in an organ-specific way. Infect Immun. 2000, 68: 4135-4144. 10.1128/IAI.68.7.4135-4144.2000.PubMed CentralView 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
- Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Liebert CA, Liu C, Lu F, Lu S, Marchler GH, Mullokandov M, Song JS, Tasneem A, Thanki N, Yamashita RA, Zhang D, Zhang N, Bryant SH: CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res. 2009, 37: D205-210. 10.1093/nar/gkn845.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Bryant SH: CD-Search: protein domain annotations on the fly. Nucleic Acids Res. 2004, 32: W327-331. 10.1093/nar/gkh454.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: D225-229. 10.1093/nar/gkq1189.PubMed CentralView ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal × version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Kumar S, Nei M, Dudley J, Tamura K: MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9: 299-306. 10.1093/bib/bbn017.PubMed CentralView ArticlePubMedGoogle Scholar
- Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157: 105-132. 10.1016/0022-2836(82)90515-0.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: A model recognition approach to the prediction of all-helical membrane protein structure and topology. Biochemistry. 1994, 33: 3038-3049. 10.1021/bi00176a037.View ArticlePubMedGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.View ArticlePubMedGoogle Scholar
- Emanuelsson O, Brunak S, von Heijne G, Nielsen H: Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007, 2: 953-971. 10.1038/nprot.2007.131.View ArticlePubMedGoogle Scholar
- Wunderlich F, Stubig H, Konigk E: Development of Plasmodium chabaudi in mouse red blood cells: structural properties of the host and parasite membranes. J Protozool. 1982, 29: 60-66.View ArticlePubMedGoogle Scholar
- Wunderlich F, Dkhil MA, Mehnert LI, Braun JV, El-Khadragy M, Borsch E, Hermsen D, Benten WP, Pfeffer K, Mossmann H, Krücken J: Testosterone responsiveness of spleen and liver in female lymphotoxin beta receptor-deficient mice resistant to blood-stage malaria. Microbes Infect. 2005, 7: 399-409. 10.1016/j.micinf.2004.11.016.View ArticlePubMedGoogle Scholar
- Krücken J, Delic D, Pauen H, Wojtalla A, El-Khadragy M, Dkhil MA, Mossmann H, Wunderlich F: Augmented particle trapping and attenuated inflammation in the liver by protective vaccination against Plasmodium chabaudi malaria. Malar J. 2009, 8: 54-10.1186/1475-2875-8-54.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
- Marti M, Good RT, Rug M, Knuepfer E, Cowman AF: Targeting malaria virulence and remodeling proteins to the host erythrocyte. Science. 2004, 306: 1930-1933. 10.1126/science.1102452.View ArticlePubMedGoogle Scholar
- Hiss JA, Przyborski JM, Schwarte F, Lingelbach K, Schneider G: The Plasmodium export element revisited. PLoS One. 2008, 3: e1560-10.1371/journal.pone.0001560.PubMed CentralView 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
- Merino EF, Fernandez-Becerra C, Durham AM, Ferreira JE, Tumilasci VF, d'Arc-Neves J, da Silva-Nunes M, Ferreira MU, Wickramarachchi T, Udagama-Randeniya P, Handunnetti SM, Del Portillo HA: Multi-character population study of the vir subtelomeric multigene superfamily of Plasmodium vivax, a major human malaria parasite. Mol Biochem Parasitol. 2006, 149: 10-16. 10.1016/j.molbiopara.2006.04.002.View ArticlePubMedGoogle Scholar
- Spielmann T, Gilberger TW: Protein export in malaria parasites: do multiple export motifs add up to multiple export pathways?. Trends Parasitol. 2010, 26: 6-10. 10.1016/j.pt.2009.10.001.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
- Daily JP, Le Roch KG, Sarr O, Fang X, Zhou Y, Ndir O, Mboup S, Sultan A, Winzeler EA, Wirth DF: In vivo transcriptional profiling of Plasmodium falciparum. Malar J. 2004, 3: 30-10.1186/1475-2875-3-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Daily JP, Le Roch KG, Sarr O, Ndiaye D, Lukens A, Zhou Y, Ndir O, Mboup S, Sultan A, Winzeler EA, Wirth DF: In vivo transcriptome of Plasmodium falciparum reveals overexpression of transcripts that encode surface proteins. J Infect Dis. 2005, 191: 1196-1203. 10.1086/428289.PubMed CentralView ArticlePubMedGoogle Scholar
- Lopez-Rubio JJ, Riviere L, Scherf A: Shared epigenetic mechanisms control virulence factors in protozoan parasites. Curr Opin Microbiol. 2007, 10: 560-568. 10.1016/j.mib.2007.10.003.View ArticlePubMedGoogle Scholar
- Zambrano-Villa S, Rosales-Borjas D, Carrero JC, Ortiz-Ortiz L: How protozoan parasites evade the immune response. Trends Parasitol. 2002, 18: 272-278. 10.1016/S1471-4922(02)02289-4.View ArticlePubMedGoogle Scholar
- Vanhamme L, Pays E, McCulloch R, Barry JD: An update on antigenic variation in African trypanosomes. Trends Parasitol. 2001, 17: 338-343. 10.1016/S1471-4922(01)01922-5.View ArticlePubMedGoogle Scholar
- 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 USA. 2008, 105: 16290-16295. 10.1073/pnas.0807404105.PubMed CentralView ArticlePubMedGoogle Scholar
- Krücken J, Mehnert LI, Dkhil MA, El-Khadragy M, Benten WP, Mossmann H, Wunderlich F: Massive destruction of malaria-parasitized red blood cells despite spleen closure. Infect Immun. 2005, 73: 6390-6398. 10.1128/IAI.73.10.6390-6398.2005.PubMed CentralView 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.