Open Access

Characterization of the repertoire diversity of the Plasmodium falciparum stevor multigene family in laboratory and field isolates

  • Jane E Blythe1, 2,
  • Makhtar Niang1,
  • Kevin Marsh3,
  • Anthony A Holder2,
  • Jean Langhorne2 and
  • Peter R Preiser1Email author
Contributed equally
Malaria Journal20098:140

https://doi.org/10.1186/1475-2875-8-140

Received: 24 April 2009

Accepted: 26 June 2009

Published: 26 June 2009

Abstract

Background

The evasion of host immune response by the human malaria parasite Plasmodium falciparum has been linked to expression of a range of variable antigens on the infected erythrocyte surface. Several genes are potentially involved in this process with the var, rif and stevor multigene families being the most likely candidates and coding for rapidly evolving proteins. The high sequence diversity of proteins encoded by these gene families may have evolved as an immune evasion strategy that enables the parasite to establish long lasting chronic infections. Previous findings have shown that the hypervariable region (HVR) of STEVOR has significant sequence diversity both within as well as across different P. falciparum lines. However, these studies did not address whether or not there are ancestral stevor that can be found in different parasites.

Methods

DNA and RNA sequences analysis as well as phylogenetic approaches were used to analyse the stevor sequence repertoire and diversity in laboratory lines and Kilifi (Kenya) fresh isolates.

Results

Conserved stevor genes were identified in different P. falciparum isolates from different global locations. Consistent with previous studies, the HVR of the stevor gene family was found to be highly divergent both within and between isolates. Importantly phylogenetic analysis shows some clustering of stevor sequences both within a single parasite clone as well as across different parasite isolates.

Conclusion

This indicates that the ancestral P. falciparum parasite genome already contained multiple stevor genes that have subsequently diversified further within the different P. falciparum populations. It also confirms that STEVOR is under strong selection pressure.

Background

Plasmodium malaria is an infectious disease of global significance with more than 500 million people suffering from the disease and at least one million dying each year [1]. The P. falciparum genome project [2] has revealed the presence of several multicopy gene families which code for hypervariable antigens that are exported to the surface of the infected host erythrocyte and represent targets of naturally acquired immunity to malaria. To date var, rif as well as stevor have been described [24] as the major variant surface antigens (VSAs). These multigene families are predominantly situated at the sub-telomeric ends of chromosomes [2, 5, 6], where gene rearrangements are frequent [7, 8] leading to high rates of recombination, thus facilitating their rapid evolution and diversity.

While some rif and var genes are found in centrally located clusters, most stevor genes are found at chromosome ends [2]. The high sequence diversity of proteins encoded by these gene families may have evolved as an immune evasion strategy that enables the parasite to establish long lasting chronic infections [9].

The P. falciparum genome contains between 50 and 60 var genes [2] that code for the Erythrocyte Membrane Protein 1 (PfEMP1) [1012]. EMP1 is linked not only to immune evasion but also to parasite sequestration and pathology and is by far the best studied multigene family in P. falciparum [11, 13, 14].

Stevor and rif share structural similarities and this relationship is emphasized by the existence of a RIFIN-STEVOR family (PF02009) in the PFAM database [15]. Characteristically, they are defined as small polypeptides with a putative signal sequence followed by a semi-conserved domain, an HVR and a conserved C-terminal domain [16]. In both protein families, two predicted transmembrane-spanning regions flank the HVR which has been predicted to form a loop on the infected red blood cell (iRBC) surface, exposing it to the host immune system and probably resulting in distinct evolution into novel functional units [3, 16]. However, STEVOR and RIFIN differ in that the HVR of RIFIN is up to 300 bp longer than the equivalent region in STEVOR. Moreover, STEVOR are a distinct, more conserved, and lower copy number family [17] than RIFIN. The rif and stevor (150–200 copies and 30–33 copies respectively in the 3D7 genome) repertoire diversity has been studied in the 3D7 parasite clone, and in several other laboratory parasite lines from different geographic origins [18, 19]. Based on these studies, it is clear that the HVR of STEVOR has significant sequence diversity both within as well as across different P. falciparum lines. These studies did not address whether or not there are ancestral stevor that can be found in different parasites.

The entire genome of the P. falciparum 3D7 parasite clone was sequenced completely first as part of the malaria genome project [2]. A number of additional lines have now also been at least partially sequenced so that there is now considerable information on the VSA repertoire on these parasite species. The ancestors of parasites isolated worldwide are most likely to be from Africa [20]. A study of var gene sequences in 12 Kilifi parasite isolates has shown a comparable number of broadly classified types as the 3D7 genome [21]. Similarly, another study of var sequences from Amazonian parasite isolates indicated that the expected gene repertoire was again similar to that of the 3D7 genome and that the Amazonian var genes were most similar to var identified in other American isolates [18].

The 3D7 var genes have been grouped into three majors types based on sequence analysis of the intron and 5' and 3' un-translated regions (UTR) [2, 22, 23], and the Kilifi var had a similar type distribution [21]. Recently, phylogenetic approaches have enabled subdivision of the RIFIN protein family into two major A- and B-subgroups [24] associated with distinct subcellular location and developmental expression patterns [25].

Such important data are currently missing for the stevor multigene family and it is not known whether field isolates contain a similar number of stevor to that of the 3D7 clone. In addition, a more extensive study of stevor from isolates worldwide is required in order to estimate both the repertoire size and the extent of diversity. Comparative analysis of long-term laboratory-adapted parasites and fresh clinical isolates has shown that laboratory-adapted lines are not significantly more or less divergent than recent isolates [26]. However, the 3D7 clone was shown to have accumulated DNA in culture through insertions of +1nucleotide/kb [27], and this might suggest divergent sequences compared to fresh isolates.

In this study, DNA and RNA sequence analysis were combined as well as phylogenetic approaches to analyse the stevor sequence repertoire and diversity in laboratory lines and fresh isolates from Kilifi in Kenya. Conserved stevor genes (orthologues) were identified in different P. falciparum isolates from different global locations. Consistent with previous studies, the HVR of the stevor gene family was found to be highly divergent both within and between isolates. Importantly, phylogenetic analysis shows some clustering of stevor sequences both within a single parasite clone as well as across different parasite isolates. These data indicate that the ancestral P. falciparum parasite genome already contained multiple stevor genes that subsequently diversified further within the different P. falciparum populations.

Methods

Parasites

Parasites used in this study were either well-characterized laboratory lines or were isolated from blood samples obtained at the KEMRI-Wellcome Trust, Kilifi, Kenya. Parasite lines and their origin are listed (Table 1). The African parasites were collected at Kilifi district hospital, situated 60 km north of Mombasa on the Kenyan coast. More than 10% of the children under five years of age, who reside within the Kilifi district, are admitted into hospital each year. In this area, P. falciparum transmitted by mosquitoes of the Anopheles gambiae complex [28], is seasonal and usually occurs after the long and short rainy seasons.
Table 1

Origin of P. falciparum clones used for in vitro culture

Clone

Parent isolate

Geographical origin

Reference

3D7

NF54

Amsterdam airport -Netherlands (unknown origin)

(Walliker et al., 1987)

A4

IT04

Brazil

(Roberts et al., 1992)

C10

1776

Malaysia

(Ang et al., 1996)

D10

FCQ-27

Papua New Guinea

(Chen et al., 1980)

Dd2

W2-MEF

Indochina

(Oduola et al., 1988)

FcB1

Unknown

Colombia

(Schrevel et al., 1994)

T9/96

T9

Thailand

(Thaithong et al., 1984)

This study was approved by the Institutional Ethical Review Board of Nanyang Technological University, Singapore, and by the Kenyan Medical Research Institute National Ethics Committee

Culture of P. falciparum laboratory lines and Kilifi isolates

Laboratory parasite lines were cultured in fresh human RBC (supplied by National UK Blood Donation Service) in RPMI 1640 complete medium supplemented with Albumax and 2 mM glutamine (Gibco) as previously described [29]. Parasitaemia was monitored daily by microscopy examination of Giemsa stained thin blood smears.

Plasmodium falciparum parasite isolates from Kilifi were recovered from liquid nitrogen frozen stock and cultured in RPMI 1640 medium supplemented with 10% AB human serum (Kilifi General Hospital, Kenya), 2 mM L-glutamine (Gibco), 37.5 mM HEPES (Gibco), 20 mM Glucose, and 2 μg/ml Gentamycin (Gibco). No fresh red blood cells (RBCs) were added and the parasites were maintained in culture for approximately 30 hours, or until they reached late trophozoite/schizont-stages. Parasite samples were separated into two fractions for later processing of parasite genomic DNA (gDNA) and RNA isolation.

Preparation of parasite genomic DNA (gDNA) and RNA

Plasmodium falciparum iRBCs were pelleted at 560× g for 5 minutes. Resulting pellets were stored at -80°C in screw top cryovials. A phenol/chloroform extraction method was used for larger (greater than 500 μl) iRBC pellets as previously described [30]. The DNeasy® kit (Qiagen) protocols 1 and 2 were used for small (less than 500 μl) blood samples according to manufacturer's recommendations. Genomic DNA was resuspended in nuclease-free water.

For RNA isolation, iRBC pellets (approximately 100 μl) were lysed in 1 ml prewarmed (37°C) Trizol LS reagent (Gibco Invitrogen) and then stored at -80°C in screw top cryovials. RNA was extracted as described previously [31], using DNase I digestion (Invitrogen) to remove contaminating DNA. RNA was resuspended in nuclease-free water.

Primers designed for stevor specific PCR/RT-PCR

The primers were assessed using multiple sequence alignments to ensure that the chosen sequences were in regions conserved within the stevor gene family, but not in other genes (for example the rif gene family). A combination of two sets of primers in a nested PCR was used RepF1, RepF2 and RepR internal-PCR primers were designed as previously described [16]. The RepF1/2 primers were redesigned to remove an opal stop codon (TGA) replacing it with (TGC). The smf1/smr1 external primers were previously described [32]. All primers were synthesized by Eurogentec/Oswel (Southampton, UK) or Research Biolabs (Singapore). The use of a nested-PCR protocol for RT-PCR was problematic due to amplification of contaminating gDNA, therefore a stevor-specific single step amplification was designed, facilitated by the complete genome sequence [2], using the primers JBSTEVORF1 and JBSTEVORR1. Primers sequences are summarized in Table 2.
Table 2

Primers sequences

Primer name

Direction

Primer Sequence (5'-3')

Tm (°C)

Smf1

sense

GAYSCAGAACTCAARGAAATWATT

46.42

Smr1

antisense

GCAGMACCAAAGYWGYAATACC

47.9

JBSTEVORF1

sense

AGATGTACCCGTGGTATATGTTCTTGCA

54.96

JBSTEVORR1

antisense

GCTAATATAAGTAGAACCAAAGCTGCAATAC

54.44

RepF1

sense

AGATGAATTCGTGGTATATRTTYTTGYTC

55.84

RepF2

sense

AGATGAATTCGGGGTTTAGTTGTGTGTGC

63.86

RepR

antisense

TTTAGGATCCAGAACCAAARYTGCAATACC

62.21

DNA and RNA products were separated by electrophoresis in 1% multi purpose agarose (Roche) in 1× Tris buffer contained 0.25 μg/ml ethidium bromide to enable UV visualization. RNA gels contained 5 mM guanidine thiocyanate.

Small-scale preparation of plasmid DNA (minipreps)

Plasmid DNA was isolated using the QIAprep spin SNAP Miniprep® Kit (Invitrogen). The QIAgen Maxiprep® kit was used according to manufacturer's protocol to isolate larger quantities of pET-24a (+) vector.

Automated DNA Sequencing

Following DNA expansion, purification, and confirmation that a stevor insert was present; the DNA was then quantified by spectrophotometry, pelleted and dried before sequencing.

Three systems have been used for sequencing reactions:

a) The BigDye® Terminator Cycle Sequencing Ready Reaction kit (Perkin-Elmer, Applied Biosystems, Warrington, UK) was used according to manufacturer's instructions. The reaction was carried out using approximately 600 ng of pCR®2.1 or pCR® II DNA (vector plus insert) per sample with either M13 forward or reverse primers in the Eppendorf Mastercycler. Unincorporated dye terminators were removed by ethanol precipitation. Samples were resuspended in 3 μl loading dye 16.5% v/v Blue Dextran/EDTA (Applied Biosystems) in formamide (Amresco, Solon USA), followed by electrophoresis and data collection on the ABI PRISM™ 377 DNA sequencer (Perkin -Elmer).

b) The BigDye® reactions were also used with T7 primers and pET-24a (+) DNA, or TrxFus forward or T7 reverse primers with pET102/D-TOPO DNA. Sequencing was performed at the National Neuroscience Institute, Singapore.

c) The MegaBACE system (Amersham Biosciences) was used. The ABI PRISM® BigDye® Terminator v1.1 Cycle Sequencing kit (Applied Biosystems) cycle sequencing mix was used according to the manufacturer's protocols. Samples were ethanol precipitated before addition of MegaBACE™ loading solution (Amersham Biosciences) and loading onto a MegaBACE 96-capillary machine.

Sequence alignment analysis (HVR region)

A fragment of approximately 280 base pairs (bp) was successfully isolated through PCR and RT-PCR amplification. However, due to the variant nature of the HVR, sequences varied in length, coding for a maximum of 82 amino acids. Clones have been obtained from multiple PCR reactions to reduce PCR bias, and were sequenced in forward and reverse directions. P. falciparum 3D7, IT and Ghanaian isolates sequences were obtained from BLAST and Tools-oligo searches within PlasmoDB [33].

Editing of nucleotide and amino acid (translated) sequences was carried out using the Autoassembler®, Chromas®, DNAstar (2001) freeware programmes. Sequences from different parasite sources and different amplification reactions were aligned by using ClustalX® [34], and Bioedit7.0.1® softwares. Multiple alignments and their corresponding chromatograms were edited by eye for discrepancies and mistaken nucleotide assignment by the automated sequencing method. BLAST analysis, version 4.4 from PlasmoDB [33] was used in order to confirm that DNA sequences were the result of the amplification of stevor exon-2 hyper-variable loop locus.

Phylogenetic analysis

Phylogenetic trees showing genetic distance were constructed using the following programmes: Molecular Evolutionary Genetics Analysis (MEGA) 2.1 or 3.0®, PAUP version 4.0 [35, 36], and Tree view® [37]. These programmes enable the construction of trees using three-building models: Neighbor-Joining, Maximum Parsimony and Maximum likelihood, taking into account the observed nucleotide frequency, substitution rate and transition/transversion ratio, based on Nei genetic distance estimates [35]. The bootstrap analysis [38] was used in order to statistically assign support to hypotheses of phylogenetic relationship among sequences estimated by the above-mentioned methods.

Results

Identification of common stevor genes in laboratory lines and Kilifi field isolates

The P. falciparum genome sequencing project [2] as well as subsequent sequencing efforts has given some insight into the diversity of stevor in different P. falciparum lines. Still these data are limited and, therefore, additional sequences from both established laboratory lines and freshly isolated parasite samples would provide valuable additional information. For this purpose, genomic DNA was extracted from five different laboratory lines of various geographic origins: A4; FcB1; C10; T9/96 and 3D7 (Table 1) and also from 10 representative blood samples taken from patients in the Kilifi region of Kenya.

Genomic DNA was amplified using the previously described stevor-specific nested primers flanking the predicted hyper-variable loop region [32] (Table 2). PCR from all laboratory lines (Figure 1A) as well as Kilifi isolates (Figure 1B) gave amplification products of the expected size (600 bp and 300 bp) indicating that these primers were able to amplify stevor sequences from both the laboratory lines and Kilifi field isolates. Similar to what was found with the P. falciparum 3D7, more than one PCR product was detected in the expected size range, amplified by both the external (550–700 bp), and internal primers (250–350 bp) (Figure 1A). This is indicative of the amplification of multiple stevor sequences from both laboratory lines and field isolates. The 300 bp PCR products from the different samples were cloned and a number of clones were sequenced in both forward and reverse directions. The sequence data obtained here from genomic DNA was further expanded by the addition of the transcribed stevor sequences which had previously been obtained by RT-PCR from RNA extracted from the same original blood samples [39].
Figure 1

Agarose gel showing PCR products obtained from laboratory line and Kilifi isolate parasite genomic DNA. Gel electrophoresis of nested stevor PCR products from genomic DNA extracted from 5 representative P. falciparum laboratory lines: A4; FcB1; C10; T9/96; and 3D7 (A) and 10 representative Kilifi parasite isolates (B). 5 μl of each product was loaded onto a 2.5% agarose gel. Results show two amplicons in the expected size ranges, products carried over in the nested PCR-step from the external PCR (approximately 600 bp) with primers smf1 and smr1, and internal PCR products (approximately 300 bp) from primers RepF1/2 and RepR. 100 bp marker (Fermentas) is shown. dH2O = negative control containing sterile water instead of template DNA.

Despite the low throughput approach, a total of 213 genomic DNA and RNA products were successfully sequenced, 133 were unique DNA sequences [see Additional file 1] coding for 108 unique amino acid sequences. Comparison with known STEVOR sequences confirmed that the PCR and RT-PCR products indeed represented STEVOR. The number of times each sequence was obtained per isolate is shown (Figure 2A). Overall, six different stevor sequences were identified from genomic DNA in the reference 3D7 clone isolated in The Netherlands but of unknown geographic origin [40] that had been included as a positive control. Interestingly, two (3D7_2_41202 and 3D7_1_181202) of the sequences obtained from the 3D7 clone maintained at the NIMR in London were absent from the published 3D7 genome stevor repertoire (Figure 2A). All other stevor sequences present in the 3D7 genome are shown with their gene identification number for the purpose of comparison. Of the 133 unique stevor DNA sequences, twenty-one were found to be present in more than one parasite line or isolate (Figure 2A, asterisk), five (5) of these paralogues were identified in A4, C10, FcB1 and D10 (blue asterisk), six (6) were found in A4, C10 and FcB1 (green asterisk), three (3) were found in T9/96, C10 and A4 (orange asterisk), while two (2) were found in 3D7 and T9/96 (red asterisk), 3D7 and A4 (brown asterisk) as well as K11 and K14 (black asterisk). Finally a single sequence was found as part of the 3D7 genome and K1640 (Figure 2A, pink asterisk). The highly represented stevor sequences most likely represent "common" stevor that can be found in a variety of parasite isolates from around the world. For a number of these more common stevor transcription was also detected by RT-PCR indicating that these genes are of potential functional relevance.
Figure 2

Identification of stevor sequences per parasite isolate. A) Number of different stevor obtained from parasite DNA and RNA. Parasite isolates are shown in different colors below the graph. The sequence name is denoted along the x-axis, each sequence has a unique identifier consisting of the parasite isolate (+ if amplified from RNA), clone number and date. The number of clones per sequence is shown on the y-axis. P. falciparum 3D7 genomic stevor sequences are included for comparison (solid black). B) Representation of the different transcribed (+) stevor obtained in this study.

From this initial analysis it is clear that there are a number of stevor sequences that are common to a number of different laboratory lines; in contrast only a single sequence (K1640_1_BD141105) from a Kilifi isolate was found to be identical to a gene (PFL2610w) from the 3D7 genome. In contrast despite the large number of stevor sequences identified in the Kilifi isolates, only a single stevor sequence was found to be present in more than one isolate. This particular sequence was transcribed in both K11 (K11+_111104) and K14 (K14+_1_111104) parasites (Figure 2B).

The RT-PCR data indicate that stevor is transcribed in all the P. falciparum lines examined with the exception of the parasite clone A4, supporting previous observations that this parasite clone does not transcribe nor express STEVOR [4, 39]. This lack of STEVOR expression in the A4 parasite is not due to the absence of the genetic information as this clone contains at least 20 different stevor genes (Figure 2A). In line with published observations [21] that field isolates transcribe multiple var at the same time, numerous unique stevor transcripts were also detected in both laboratory and patient samples (Figure 2B).

Conservation of stevor across different parasite lines

The observation that identical HVR sequences could be found in different parasite lines raised the possibility that stevor could be grouped into distinct subgroups with different biological functions, similar to the rif and pir multigene family [25, 41, 42]. For this purpose the sequences identified here, together with the corresponding regions from the 3D7 stevor repertoire were clustered into similarity groups by constructing a Neighbor Joining as well as a Minimum Evolution Tree (Figure 3).
Figure 3

Phylogenetic trees of stevor DNA sequences. Neighbor-joining (A) and Minimum Evolution (B) methods were used for generation of phylogenetic trees of 152 HVR stevor alleles. Trees are based upon a nucleotide alignment of the HVRs and condensed to show only those branches with a bootstrap value of 95% significance and above. Sequences are color-coded according to originating isolate: Kilifi-yellow; T9/96-dark blue; 3D7-pink; A4 and C10-green; FcB1-light blue; D10-light green. Orange open-circle around a branch point indicates the major stevor sub-grouping, blue open-circles highlight T9/96 sub-groups, green open-circles highlight sub-groupings including FcB1, C10 and A4 sequences, yellow open-circles highlight Kilifi sequence sub-groupings.

The Neighbor Joining tree (Figure 3A) clustered the 152 HVR DNA sequences into two groups with one large sub-group containing approximately one-third of the stevor sequences (Figure 3A, orange branch). This subgroup contained sequences from all parasites except the D10 laboratory line and the K14, K11 and K6 Kilifi field isolates. At this stage it is not clear whether the absence of this subgroup of stevor in these parasites is a reflection of the limited sequence data available or truly represents a lack of this subgroup. It is also interesting to note that there is further sub-clustering of sequences within this branch, where clusters either contain sequences from the same line or alternatively from a number of different parasites. The second large group containing the remaining two-thirds of stevor sequences appears to contain numerous branch points of related sequences rather than containing a single defined cluster. These data would indicate that overall the sequence divergence is much higher in this group. Sequences obtained from a single parasite line also appear to form distinct clusters, with each cluster containing a number of different related sequences. This can be clearly seen for T9/96 stevor (Figure 3A, blue circles) that contains at least three distinct clusters. Similar clustering is also seen for sequences from A4, C10, FcB1 and D10 with the majority of sequences falling within just four sub-groups (Figure 3A, green circles) while the majority of Kilifi isolates stevor sequences were contained within seven clusters (Figure 3A, yellow circles). Importantly, these specific sequence clusters often also contain sequences from other parasite lines or isolates indicating that they do not result from a specific phylogenetic branch. Rather it appears that while there is a trend for some sequences from the same parasite to cluster together, most are distributed across the tree and cluster with sequences from other strains. To further validate these findings, the stevor sequences were clustered using the Minimum Evolution tree construction method at 95% significance bootstrap values (Figure 3B). This approach produced the same clusters as the Neighbor Joining method (compare Figure 3A with 3B). The 3D7 stevor repertoire is distributed throughout the trees (Figure 3A, B; pink points) with clusters containing a maximum of two 3D7 sequences.

To assess whether or not this clustering of stevor DNA sequences is maintained at the amino acid level, 108 unique STEVOR amino acid sequences were used to construct a Neighbor Joining as well as a Minimum Evolution Tree with a bootstrap value of 85% or higher (Figure 4A, B). The two approaches resulted in nearly identical trees with a number of broad clusters containing sequences from different lines and isolates. Using this analysis a single major sub-group was observed (highlighted by an orange circle) reflecting an apparent trend for sequences of the same geographic region to be more closely related (Figure 4A, B). More precisely, one obvious cluster was observed containing majority of Kilifi protein sequences (Figure 4A, orange circle, yellow points). However, as already identified in the DNA sequence trees, clusters are not completely exclusive and a number of sequences from laboratory lines are found within this cluster (Figure 4A, orange circle, green, pink, blue and light blue points). Conversely a number of Kilifi isolate sequences are found distributed across the other branches of the tree where they cluster both with other Kilifi sequences as well as sequences from a range of laboratory lines (Figure 4, yellow points).
Figure 4

Phylogenetic trees of stevor HVR protein sequences. Unrooted phylogenetic trees using neighbor-joining and minimum evolution methods showing all branches (A) and 85% condensed trees (B) using 108 amino acid sequences. Trees are based upon an alignment of the HVRs then condensed to show only those branches with a bootstrap value of 85% significance and above.

Similar to the outcome of the DNA sequence analysis, the T9/96 amino acid sequences clustered within three small sub-groups (Figure 4, blue circles) while sequences from laboratory lines A4 (green), C10 (green), FcB1 (light blue) could be found dispersed in different clusters on the tree (Figure 4). Interestingly, the significant clustering seen within the trees produced from the nucleotide sequences (Figure 3) was also maintained at the amino acid level. This clustering is consistent using either of two methods for generating the phylogenetic trees.

This analysis was further extended to include all STEVOR sequences including those from published databases (IT and Ghanaian isolates) and 75 South American region sequences obtained from Albrecht et al. (2006)[18]. A total of 226 HVR amino acid sequences were used for the analysis (Figure 5). However, no obvious cluster was observed, though a number of non-significant groups were present when the tree was condensed at 50% bootstrap values. The high divergence within the HVR encoding the predicted erythrocyte surface-exposed region of the STEVOR family is likely in response to host immunity.
Figure 5

50% condensed tree using 226 amino acid sequences. Sequences are color-coded according to originating isolate: Kilifi field isolates-yellow; T9/96-dark blue; 3D7-pink; A4 and C10-green; FcB1-light blue; D10-light green, IT-black, Ghanaian field isolate-red. Orange open-circle indicates a major stevor sub-group, blue open-circles highlight T9/96 sub-groups, green open-circles highlight sub-groupings including FcB1, C10 and A4 sequences, yellow open-circles highlight Kilifi sequence sub-groupings.

Discussion

The human malaria parasite P. falciparum is able to evade the host's immune response by expressing a range of variable antigens on the surface of the infected erythrocyte. Several genes are potentially involved in this process with the var, rif and stevor [2] multigene families being the most likely candidates and coding for rapidly evolving proteins. The rapid changes are believed to result from a strong immune-mediated diversifying selection [19, 43, 44].

In this study, the stevor sequence repertoire and diversity were analysed in laboratory lines and Kilifi field isolates. As expected the stevor family was found in all P. falciparum genomes and was transcribed in all field and laboratory parasites (exception for the A4 clone) with multiple transcripts being detected. This result is in line with previous studies in which numerous var as well as stevor transcripts were detected in single patient samples using a microarray [45]. From the RT-PCR and sequencing data it is clear that both in cultured laboratory lines and in patient samples a significant proportion of the stevor repertoire is transcribed at the same time. In isolate K5 a total of 16 different sequences were detected; and despite the close phylogenetic relationship between these sequences (Figure 3) at this stage it can not be ruled out that these transcripts arose from a multiclonal infection.

Recent phylogenetic studies have shown that the rif family form two clear subgroups [24] that potentially can be further divided into functionally distinct subgroups [46]. This is clearly not the case for stevor where the presence of distinct subgroups is not that strongly supported. Still there is sufficient evidence to support some clustering into different groups. Importantly, the identification of identical sequences in laboratory lines from different geographic regions as well as field samples would imply some functional constraints acting upon this sequence. In addition sequences within the Kilifi isolates appear to be more closely related to each other than to any of the laboratory parasites, and no common stevor were found between Kilifi except for a single stevor sequence shared between one Kilifi isolate and 3D7. The stevor sequence found in two Kilifi isolates (K11 and K14) may represent a common sequence detected in the two isolates and also seen in a number of laboratory lines.

Conclusion

Conserved genes were identified in different P. falciparum isolates from different global locations for example: A4, C10 and FcB1 in particular, contain several closely related stevor genes, suggesting that the ancestral P. falciparum parasite genome already had multiple stevor genes that may have subsequently diversified further within the different P. falciparum populations. The high sequence variability observed is consistent with the HVR of STEVOR being under strong selection pressure. On the other hand, there is only weak evidence that STEVOR like the RIFIN can be subdivided into distinct phylogenetic subgroups implying that STEVOR function has less diversified in P. falciparum.

Notes

Declarations

Acknowledgements

We like to thank Pete Bull and Vandana Thathy for critical reading of the manuscript. JEB was in receipt of a studentship from the Medical Research Council (UK). This work was supported by a Biomedical Research Council (BMRC) Singapore Grant (04/1/22/19/364 and 07/1/22/19/514) and the BioMalPar European Network of Excellence supported by a European grant (LSHP-CT-2004-503578) from Priority 1 "Life Sciences, Genomics and Biotechnology for Health" in the 6th Framework Programme.

Authors’ Affiliations

(1)
Nanyang Technological University, School of Biological Sciences
(2)
Division of Parasitology, National Institute for Medical Research
(3)
Kenya Medical Research Institute (KEMRI)

References

  1. Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434 (7030): 214-217. 10.1038/nature03342.PubMed CentralView ArticlePubMedGoogle Scholar
  2. 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
  3. Kyes SA, Rowe JA, Kriek N, Newbold CI: Rifins: a second family of clonally variant proteins expressed on the surface of red cells infected with Plasmodium falciparum. Proc Natl Acad Sci U S A. 1999, 96 (16): 9333-9338. 10.1073/pnas.96.16.9333.PubMed CentralView ArticlePubMedGoogle Scholar
  4. Niang M, Yan Yam X, Preiser PR: The Plasmodium falciparum STEVOR multigene family mediates antigenic variation of the infected erythrocyte. PLoS Pathog. 2009, 5 (2): e1000307-10.1371/journal.ppat.1000307.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Sam-Yellowe TY, Florens L, Johnson JR, Wang T, Drazba JA, Le Roch KG, Zhou Y, Batalov S, Carucci DJ, Winzeler EA, Yates JR: A Plasmodium gene family encoding Maurer's cleft membrane proteins: structural properties and expression profiling. Genome Res. 2004, 14 (6): 1052-1059. 10.1101/gr.2126104.PubMed CentralView ArticlePubMedGoogle Scholar
  6. Thompson JK, Rubio JP, Caruana S, Brockman A, Wickham ME, Cowman AF: The chromosomal organization of the Plasmodium falciparum var gene family is conserved. Mol Biochem Parasitol. 1997, 87 (1): 49-60. 10.1016/S0166-6851(97)00041-8.View ArticlePubMedGoogle Scholar
  7. Freitas-Junior LH, Bottius E, Pirrit LA, Deitsch KW, Scheidig C, Guinet F, Nehrbass U, Wellems TE, Scherf A: Frequent ectopic recombination of virulence factor genes in telomeric chromosome clusters of P. falciparum. Nature. 2000, 407: 1018-1022. 10.1038/35039531.View ArticlePubMedGoogle Scholar
  8. Hernandez-Rivas R, Hinterberg K, Scherf A: Compartmentalization of genes coding for immunodominant antigens to fragile chromosome ends leads to dispersed subtelomeric gene families and rapid gene evolution in Plasmodium falciparum. Mol Biochem Parasitol. 1996, 78: 137-148. 10.1016/S0166-6851(96)02618-7.View ArticlePubMedGoogle Scholar
  9. 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
  10. Baruch DI, Pasloske BL, Singh HB, Bi X, Ma XC, Feldman M, Taraschi TF, Howard RJ: Cloning the P. falciparum gene encoding PfEMP1, a malarial variant antigen and adherence receptor on the surface of parasitized human erythrocytes. Cell. 1995, 82 (1): 77-87. 10.1016/0092-8674(95)90054-3.View ArticlePubMedGoogle Scholar
  11. Smith JD, Chitnis CE, Craig AG, Roberts DJ, Hudson-Taylor DE, Peterson DS, Pinches R, Newbold CI, Miller LH: Switches in expression of Plasmodium falciparum var genes correlate with changes in antigenic and cytoadherent phenotypes of infected erythrocytes. Cell. 1995, 82: 101-110. 10.1016/0092-8674(95)90056-X.PubMed CentralView ArticlePubMedGoogle Scholar
  12. 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 (1): 89-100. 10.1016/0092-8674(95)90055-1.View ArticlePubMedGoogle Scholar
  13. Buffet PA, Gamain B, Scheidig C, Baruch D, Smith JD, Hernandez-Rivas R, Pouvelle B, Oishi S, Fujii N, Fusai T, Parzy D, Miller LH, Gysin J, Scherf A: Plasmodium falciparum domain mediating adhesion to chondroitin sulfate A: a receptor for human placental infection. Proc Natl Acad Sci USA. 1999, 96: 12743-12748. 10.1073/pnas.96.22.12743.PubMed CentralView ArticlePubMedGoogle Scholar
  14. Rowe JA, Moulds JM, Newbold CI, Miller LH: Plasmodium falciparum rosetting mediated by a parasite-variant erythrocyte membrane protein and complement-receptor 1. Nature. 1997, 388: 292-295. 10.1038/40269.View ArticlePubMedGoogle Scholar
  15. Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, Eddy SR, Sonnhammer EL, Bateman A: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, 34: D247-251. 10.1093/nar/gkj149.PubMed CentralView ArticlePubMedGoogle Scholar
  16. Cheng Q, Cloonan N, Fischer K, Thompson J, Waine G, Lanzer M, Saul A: stevor and rif are Plasmodium falciparum multicopy gene families which potentially encode variant antigens. Mol Biochem Parasitol. 1998, 97: 161-176. 10.1016/S0166-6851(98)00144-3.View ArticlePubMedGoogle Scholar
  17. Gardner MJ, Tettelin H, Carucci DJ, Cummings LM, Aravind L, Koonin EV, Shallom S, Mason T, Yu K, Fujii C, Pederson J, Shen K, Jing J, Aston C, Lai Z, Schwartz DC, Pertea M, Salzberg S, Zhou L, Sutton GG, Clayton R, White O, Smith HO, Fraser CM, Adams MD, Venter JC, Hoffman SL: Chromosome 2 sequence of the human malaria parasite Plasmodium falciparum. Science. 1998, 282: 1126-1132. 10.1126/science.282.5391.1126.View ArticlePubMedGoogle Scholar
  18. Albrecht L, Merino EF, Hoffmann EH, Ferreira MU, de Mattos Ferreira RG, Osakabe AL, Dalla Martha RC, Ramharter M, Durham AM, Ferreira JE, Del Portillo HA, Wunderlich G: Extense variant gene family repertoire overlap in Western Amazon Plasmodium falciparum isolates. Mol Biochem Parasitol. 2006, 150 (2): 157-165. 10.1016/j.molbiopara.2006.07.007.View ArticlePubMedGoogle Scholar
  19. Lavazec C, Sanyal S, Templeton TJ: Hypervariability within the Rifin, Stevor and Pfmc-2TM superfamilies in Plasmodium falciparum. Nucleic Acids Res. 2006, 34: 6696-6707. 10.1093/nar/gkl942.PubMed CentralView ArticlePubMedGoogle Scholar
  20. Wootton JC, Feng X, Ferdig MT, Cooper RA, Mu J, Baruch DI, Magill AJ, Su XZ: Genetic diversity and chloroquine selective sweeps in Plasmodium falciparum. Nature. 2002, 418: 320-323. 10.1038/nature00813.View ArticlePubMedGoogle Scholar
  21. Bull PC, Berriman M, Kyes S, Quail MA, Hall N, Kortok MM, Marsh K, Newbold CI: Plasmodium falciparum variant surface antigen expression patterns during malaria. PLoS Pathog. 2005, 1: e26-10.1371/journal.ppat.0010026.PubMed CentralView ArticlePubMedGoogle Scholar
  22. Mercereau-Puijalon O, Barale JC, Bischoff E: Three multigene families in Plasmodium parasites: facts and questions. Int J Parasitol. 2002, 32 (11): 1323-1344. 10.1016/S0020-7519(02)00111-X.View ArticlePubMedGoogle Scholar
  23. Voss TS, Thompson JK, Waterkeyn J, Felger I, Weiss N, Cowman AF, Beck HP: Genomic distribution and functional characterisation of two distinct and conserved Plasmodium falciparum var gene 5' flanking sequences. Mol Biochem Parasitol. 2000, 107: 103-115. 10.1016/S0166-6851(00)00176-6.View ArticlePubMedGoogle Scholar
  24. Joannin N, Abhiman S, Sonnhammer EL, Wahlgren M: Sub-grouping and sub-functionalization of the RIFIN multi-copy protein family. BMC Genomics. 2008, 9: 19-10.1186/1471-2164-9-19.PubMed CentralView ArticlePubMedGoogle Scholar
  25. Petter M, Haeggstrom M, Khattab A, Fernandez V, Klinkert MQ, Wahlgren M: Variant proteins of the Plasmodium falciparum RIFIN family show distinct subcellular localization and developmental expression patterns. Mol Biochem Parasitol. 2007, 156: 51-61. 10.1016/j.molbiopara.2007.07.011.View ArticlePubMedGoogle Scholar
  26. Volkman SK, Sabeti PC, DeCaprio D, Neafsey DE, Schaffner SF, Milner DA, Daily JP, Sarr O, Ndiaye D, Ndir O, Mboup S, Duraisingh MT, Lukens A, Derr A, Stange-Thomann N, Waggoner S, Onofrio R, Ziaugra L, Mauceli E, Gnerre S, Jaffe DB, Zainoun J, Wiegand RC, Birren BW, Hartl DL, Galagan JE, Lander ES, Wirth D: A genome-wide map of diversity in Plasmodium falciparum. Nat Genet. 2007, 39: 113-119. 10.1038/ng1930.View ArticlePubMedGoogle Scholar
  27. Jeffares DC, Pain A, Berry A, Cox AV, Stalker J, Ingle CE, Thomas A, Quail MA, Siebenthall K, Uhlemann AC, Kyes S, Krishna S, Newbold C, Dermitzakis ET, Berriman M: Genome variation and evolution of the malaria parasite Plasmodium falciparum. Nat Genet. 2007, 39: 120-125. 10.1038/ng1931.PubMed CentralView ArticlePubMedGoogle Scholar
  28. Mbogo CN, Snow RW, Kabiru EW, Ouma JH, Githure JI, Marsh K, Beier JC: Low-level Plasmodium falciparum transmission and the incidence of severe malaria infections on the Kenyan coast. Am J Trop Med Hyg. 1993, 49: 245-253.PubMedGoogle Scholar
  29. Trager W, Jensen JB: Human malaria parasites in continuous culture. Science. 1976, 193: 673-675. 10.1126/science.781840.View ArticlePubMedGoogle Scholar
  30. Beck HP: Extraction and purification of Plasmodium parasite DNA. Methods Mol Med. 2002, 72: 159-163.PubMedGoogle Scholar
  31. Kyes S, Pinches R, Newbold C: A simple RNA analysis method shows var and rif multigene family expression patterns in Plasmodium falciparum. Mol Biochem Parasitol. 2000, 105: 311-315. 10.1016/S0166-6851(99)00193-0.View ArticlePubMedGoogle Scholar
  32. Kaviratne M, Khan SM, Jarra W, Preiser PR: Small variant STEVOR antigen is uniquely located within Maurer's clefts in Plasmodium falciparum-infected red blood cells. Eukaryot Cell. 2002, 1: 926-935. 10.1128/EC.1.6.926-935.2002.PubMed CentralView ArticlePubMedGoogle Scholar
  33. Plamsmo DB: [http://www.plasmodb.org]
  34. Higgins DG, Sharp PM: CLUSTAL: a package for performing multiple sequence alignment on a microcomputer. Gene. 1988, 73: 237-244. 10.1016/0378-1119(88)90330-7.View ArticlePubMedGoogle Scholar
  35. Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
  36. Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (and other methods) 4.0 Beta. 2003, Sunderland, MA: Sinauer AssociatesGoogle Scholar
  37. Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
  38. Felsenstein J: Phylogenies from molecular sequences: inference and reliability. Annu Rev Genet. 1988, 22: 521-565. 10.1146/annurev.ge.22.120188.002513.View ArticlePubMedGoogle Scholar
  39. Blythe JE, Yam XY, Kuss C, Bozdech Z, Holder AA, Marsh K, Langhorne J, Preiser PR: Plasmodium falciparum STEVOR proteins are highly expressed in patient isolates and located in the surface membranes of infected red blood cells and the apical tips of merozoites. Infect Immun. 2008, 76: 3329-3336. 10.1128/IAI.01460-07.PubMed CentralView ArticlePubMedGoogle Scholar
  40. Walliker D, Quakyi IA, Wellems TE, McCutchan TF, Szarfman A, London WT, Corcoran LM, Burkot TR, Carter R: Genetic analysis of the human malaria parasite Plasmodium falciparum. Science. 1987, 236: 1661-1666. 10.1126/science.3299700.View ArticlePubMedGoogle Scholar
  41. 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, Ralph SA, 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
  42. 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
  43. Mu J, Joy DA, Duan J, Huang Y, Carlton J, Walker J, Barnwell J, Beerli P, Charleston MA, Pybus OG, Su XZ: Host switch leads to emergence of Plasmodium vivax malaria in humans. Mol Biol Evol. 2005, 22: 1686-1693. 10.1093/molbev/msi160.View ArticlePubMedGoogle Scholar
  44. 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
  45. 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
  46. Petter M, Bonow I, Klinkert MQ: Diverse expression patterns of subgroups of the rif multigene family during Plasmodium falciparum gametocytogenesis. PLoS ONE. 2008, 3: e3779-10.1371/journal.pone.0003779.PubMed CentralView ArticlePubMedGoogle Scholar

Copyright

© Blythe et al; licensee BioMed Central Ltd. 2009

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.

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Advertisement