Previous Article | Next Article ![]()
Infection and Immunity, December 2007, p. 5565-5574, Vol. 75, No. 12
0019-9567/07/$08.00+0 doi:10.1128/IAI.00405-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Susanne K. Miller,1,2,
Tobias Sargeant,1,2
Robert T. Good,1
Anthony N. Hodder,1
Terence P. Speed,1
Tania F. de Koning-Ward,1 and
Brendan S. Crabb1*
The Walter and Eliza Hall Institute of Medical Research, Melbourne, Victoria 3050, Australia,1 Department of Medical Biology, University of Melbourne, Melbourne, Australia2
Received 18 March 2007/ Returned for modification 14 May 2007/ Accepted 4 September 2007
|
|
|---|
|
|
|---|
Accordingly, the P. falciparum serine repeat antigens (SERAs) also have potential as drug targets. These proteins are expressed in the parasitophorous vacuole of mature, erythrocytic-stage parasites (8, 13, 20), and all of them possess a papain protease-like central domain flanked by two novel and relatively cysteine-rich domains. Several SERA proteins, termed the "serine-type" SERAs, possess an atypical serine residue in the active-site position of the enzyme domain usually occupied by a canonical cysteine. The enzyme domain from at least one serine-type SERA, SERA5, retains a degree of catalytic function; however, it is uncertain if this proteolytic activity has biological relevance (12). The P. falciparum genome contains a total of nine SERA genes, six of which encode serine-type SERA proteins and three of which encode more classical cysteine-type SERAs (2, 11, 20). These two types of SERA proteins form distinct phylogenetic groups not just within P. falciparum but also across the Plasmodium genus, suggesting different biological roles for members of each type (3, 12). The precise nature and relevance of these roles remain to be determined.
While it has been established that all members of the P. falciparum SERA gene family, eight of which are tandemly arrayed genes on chromosome 2 (SERA1 to SERA8), and a ninth gene on chromosome 9 (SERA9), are expressed in a coregulated manner late in the blood-stage cycle, it is clear that each member is expressed at very different levels (2, 4, 18-20). This expression pattern is similar across all parasite lines examined to date, and multiple SERAs are expressed in the one parasite, suggesting that family members are not differentially expressed. Proteomic data, although not strictly quantitative, suggest a strong concordance between mRNA and protein levels, with SERA5 peptides detected with extraordinary frequency (16). SERA5 ranked fourth (159 peptides) in a list of 712 trophozoite/schizont proteins detected. SERA4 peptides were also well represented in this analysis, ranking 79th (22 peptides), together with SERA6 (208th; nine peptides), SERA3 (209th; nine peptides), SERA7 (380th; three peptides), and SERA9 (457th; two peptides). A separate proteomic analysis of P. falciparum schizonts yielded almost identical results (17). It is important that some SERA family members that are expressed weakly in the erythrocytic cycle may be predominantly expressed at other stages.
Several lines of evidence suggest that a possible role for some members of this family is to promote merozoite egress from infected erythrocytes. Apart from their appropriate location in the parasitophorous vacuole, it has been shown that antibodies specific for SERA5 appear to impede schizont rupture (23). Related to this, the targeted deletion of a P. falciparum cysteine-type SERA homologue (SERA8) in P. berghei, which is predominantly expressed in sporozoites, results in the arrest of sporozoite egress from oocysts in the mosquito (1). While it is attractive to contemplate a generic role in parasite egress for this family of proteins at different life stages, it remains to be demonstrated if this is the case. Moreover, the reason for so many SERA genes, and for a variable number between species (3), remains a mystery.
In order to probe the biological role of the SERAs during the blood stage, we previously attempted to disrupt each of the SERA genes clustered on chromosome 2 using a single-crossover recombination approach (7). Using this technology, we disrupted the peripherally located genes (SERA2, SERA3, SERA7, and SERA8) but were unsuccessful in targeting the peripherally located gene SERA1 or any of the more strongly expressed genes, SERA4, SERA5, and SERA6 (20). Since then, another more powerful gene-targeting technology has been developed for P. falciparum. This approach involves selecting double-crossover integrants using a combination of both positive and negative selectable markers (9). In this study, we used this method as well as repeating the single-crossover approach in an attempt to disrupt the two most strongly expressed members of the SERA family, SERA4 and SERA5, and the remaining two serine-type SERAs, SERA1 and SERA9. While the SERA5 gene remained refractory to deletion, we were able to disrupt SERA1, SERA4, and SERA9. SERA4-null parasites displayed a twofold increase in the expression of SERA5 mRNA. We also performed a comprehensive analysis of genome sequences of Plasmodium species parasites and demonstrate that, in contrast to the cysteine-type members, the serine-type SERAs form a monophyletic group that have expanded postspeciation. Together, these data provide evidence for a key, probably single, role for serine-type SERAs in the erythrocytic cycle.
|
|
|---|
Parasite culture and transfection.
P. falciparum asexual erythrocytic-stage parasites (D10 or 3D7 lines) were cultivated at 37°C, 5% CO2, 1% N2, and 1% O2 (27) in 4% hematocrit using human O-positive erythrocytes and fed every second day with complete culture medium (RPMI-HEPES supplemented with 0.2% NaHCO3, 5% heat inactivated human serum, and 5% Albumax). Enriched trophozoite and schizont preparations used in Western blot analyses were obtained by Percoll purification (29). Ring-stage parasites (
5% parasitemia) were transfected with 100 µg of purified plasmid DNA (plasmid maxi kit; Qiagen) (6) using modified electroporation conditions (10), and drug cycling commenced according to methods described previously by Crabb et al. and Duraisingh et al. (7, 9). Briefly, parasites were cultured in a 10-cm culture dish for 48 h prior to selection with 2.5 nM WR99210. Transfected parasites were visible after 3 weeks of continuous culture. Parasites containing integrated forms of the plasmid were obtained by drug cycling, where drug was removed for 2 to 3 weeks before being reapplied.
Southern blotting. Total parasite DNA was isolated from saponin (0.15%)-lysed trophozoite- and schizont-infected erythrocytes at 5% parasitemia. Buffer A (50 mM NaAc [pH 5.2], 100 mM NaCl, 1 mM EDTA) and sodium dodecyl sulfate (SDS) (3% [wt/vol]) were added to the saponin-lysed pellet and mixed by inversion. The resultant mixture was extracted using phenol-chloroform (1:1) (Progen). The aqueous phase was ethanol precipitated as described previously (26), and the DNA pellets were air dried and resuspended in 500 µl of Tris-EDTA. Manipulation of recombinant DNA and analysis of nucleic acids by Southern blot hybridization were carried out using standard procedures (26).
SDS-polyacrylamide gel electrophoresis and Western blotting. Parasite proteins extracted from trophozoite and schizont preparations that had been purified by Percoll density gradient (64%) were separated by SDS-polyacrylamide gel electrophoresis (15) on 4 to 20% polyacrylamide gradient Tris-HEPES-sodium dodecyl sulfate-polyacrylamide gel electrophoresis gels under reducing conditions and transferred onto polyvinylidene difluoride membranes (Amersham) using the Bio-Rad Mini-Protean II electrophoresis cell and a Mini Tran-Blot electrophoretic transfer cell apparatus according to the manufacturer's instructions. After transfer, polyvinylidene difluoride membranes were blocked in 5% (wt/vol) skim milk in phosphate-buffered saline overnight at 4°C. Membranes were probed with rabbit anti-SERA-specific polyclonal antibodies. The antibodies and methodology used were described in a previous study (20). Horseradish peroxidase-conjugated sheep anti-rabbit or sheep anti-human immunoglobulin G (Silenus Laboratories) diluted 1/3,000 was used for detection.
Growth rate assay. Ring-stage wild-type and knockout parasites were synchronized by sorbitol lysis twice at 4-h intervals and plated in duplicate at 0.5% parasitemia into complete culture medium containing a 4% hematocrit in the absence of drug selection. Methanol-fixed and Giemsa-stained thin blood smears were made every 48 h for a total of 28 days. Fresh medium was added every 48 h, and cultures were maintained below a parasitemia of 5%.
Microarray. Microarrays comprised of 6,912 P. falciparum 70-mer oligonucleotides were printed in duplicate on polylysine-coated microscope slides (Operon Technologies) as described previously (4). Highly synchronous P. falciparum D10 parasite cultures were obtained by double-sorbitol treatment (5%; Sigma) and used to measure SERA gene expression levels between mature (trophozoite/schizont) and young (ring) D10 parasites. For this, parasitized erythrocytes were lysed by saponin lysis (0.15%), and the RNA was extracted using RNeasy Minicolumns (Qiagen). Total RNA (20 µg) was reverse transcribed overnight using SuperScript II (Gibco BRL). The reactive amino group of aa-dUTP was used to conjugate 2 µg of the purified cDNA preparations with 100 µg of the fluorescent dye esters N-hydroxysuccinimide-Cy3 (Amersham) and N-hydroxysuccinimide-Cy5 (Amersham). Microarrays were hybridized overnight at 42°C with 500 ng of each fluorescently labeled probe in a solution containing 25% (vol/vol) formamide, 5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), and 0.1% (wt/vol) SDS. Arrays were washed twice at room temperature with 1x SSC-0.2% SDS, followed by two stringent washes (0.1x SSC) prior to scanning on a ScanArray 4000 apparatus (GSI Lumonics). Data were extracted from the raw images and mined using GeneSpring (Silicon Genetics).
Quantitative reverse transcription (RT)-PCR.
Parasitized erythrocytes at 36 to 40 h postsynchronization were lysed by saponin lysis, and RNA was extracted using RNeasy Minicolumns (Qiagen). Total RNA (5 µg) was reverse transcribed using SuperScript II (Gibco BRL), and 500 ng of the resultant cDNA was added to the QuantiTect SYBR green PCR mixture (Qiagen) containing the gene-specific oligonucleotides specified in Table S1 in the supplemental material. PCR was performed using a LightCycler apparatus (Roche), and reaction conditions were as follows: 95°C (15 min), followed by 40 cycles of 95°C (15 s), 45°C (30 s), and 67°C (20 s). Serial dilutions of P. falciparum D10 genomic DNA were used as a standard reference control. The same preparation procedure was followed when performing quantitative PCR on parasite preparations isolated following a time course for wild-type D10 and D10-
SERA4 parasite lines. Parasites in each line were collected at 30 h postsynchronization and subsequent 4-h intervals for 20 h and prepared as described above.
Database searches and phylogenetic analysis. HMMER (http://hmmer.janelia.org) profiles were generated for the conserved protease and C-terminal domains of aligned P. falciparum 3D7 SERA proteins. These profiles were used to search a database of six-frame translations of genomic sequences as well as protein sequences of apicomplexan species. From the results, matches were extracted for the following species: P. falciparum (three strains), Plasmodium reichenowi, Plasmodium vivax, Plasmodium knowlesi, Plasmodium yoelii, Plasmodium berghei, Plasmodium chabaudi, Plasmodium gallinaceum, and Theileria annulata. In the case of T. annulata, only the protease domain was matched by the hidden Markov model (HMM) profiles generated. Alignment of the corresponding T. annulata protein sequence with P. falciparum cysteine-type SERA sequences indicated the presence of a divergent (and therefore unrecognized by the C-terminal HMM profile) but homologous C-terminal domain in which all seven cysteine residues were conserved. The C-terminal portion of the T. annulata sequence was extracted from the multiple alignment and used for further analyses.
A multiple alignment was constructed from concatenated protease and C-terminal domains using Hmmalign. This alignment consisted of 354 characters over 71 sequences. P. reichenowi SERA8 was removed because the obtained sequence was partial. For cysteine-type SERAs lacking a C-terminal domain (e.g., P. falciparum SERA8, P. yoelii SERA5, and P. vivax SERA12), gaps were inserted to pad the alignment. Phylogenetic analyses were carried out using MrBayes (24) with fixed Blosum62 transition probabilities and substitution rates of
(4) + I. In parallel, topologies were inferred using neighbor-joining and maximum likelihood methods as implemented by Phylip (Phylogeny Inference Package version 3.6; J. Felsenstein, Department of Genome Sciences, University of Washington, Seattle). Internal branch splits congruent between Markov chain Monte Carlo (MCMC), maximum likelihood, and neighbor-joining methods were recorded. In order to examine support for the Cys-Ser split in the inferred phylogeny, we removed all trees from the MCMC sampling that contained this branch partition and reanalyzed the remainder of the data.
|
|
|---|
![]() View larger version (34K): [in a new window] |
FIG. 1. Targeted disruption of four serine-type SERA genes. (A) Schematic representation detailing the outcome of expected integration events in SERA knockout parasite lines. Solid lines under the exons represent the location of the gene targeting regions used to clone fragments into pHH1 (SERA1), pTK SERA4, pTK SERA5, and pTK SERA9. The targeting regions of the coding sequence are indicated by the gray shaded regions. (B) Southern blots summarizing the result of the attempted targeted disruption of each SERA gene using restriction-digested genomic DNA from each transfected line following drug cycling and parental parasites. Genomic DNA was digested BstBI (B)/SphI (S) for SERA1, AflIII (A) for SERA4, XmnI (X) for SERA5, and BclI (Bc)/NotI/SwaI (Sw) for SERA9. The membranes were probed with each respective targeting sequence shown above (A). The positions of endogenous (E), plasmid (P), and integration-specific fragment (In) events that result from digestion with the respective enzymes and the probing of blots with the specific targeting region (gray) are indicated, and their molecular weights (in thousands) are given below each relevant blot. RE, restriction enzyme.
|
In the case of SERA1, Southern blot analysis identified a band of 4.6 kb corresponding to the endogenous SERA1 locus in wild-type 3D7. The disappearance of this fragment together with the appearance of fragments representing the targeted locus confirms successful single-crossover homologous recombination (Fig. 1). It should be noted that plasmids that integrate by single-crossover recombination often insert more than one copy into the locus, and if this has occurred, a band corresponding to that expected for the episomal plasmid would be observed. This was the case for SERA1, where, based on band intensity, one additional copy was inserted (Fig. 1).
With respect to the targeting of the SERA4 gene, hybridization with the FL1 and FL2 targeting sequences of SERA4 identified single bands of 5.4 kb and 1.8 kb in D10 corresponding to the endogenous gene. In the transfected D10-SERA4/C0 line, which had not undergone ganciclovir selection, the bands at 2.6 kb and 4.0 kb indicated the presence of plasmid pTK
SERA4. However, in D10-
SERA4 parasites, which had undergone negative selection pressure, single bands at 6.0 kb and 2.5 kb indicated that the plasmid had integrated into the SERA4 gene via the expected double-crossover event (Fig. 1). Therefore, although we were previously unable to disrupt SERA4 using single-crossover technology, a SERA4 "knockout" was obtained by using the negative selection procedure.
In an attempt to delete SERA5, plasmid pHTK
SERA5 was transfected into D10 parasites on four separate occasions, and parasites were cycled as described above. Parasites were analyzed by pulsed-field gel electrophoresis, and in all rounds of drug cycling hybridizations of both the hDHFR probe and the pGEM probe to D10-pHTK
SERA5, either parasites produced a smeared multiple banding pattern consistent with episomal maintenance of the plasmid or the hDHFR probe appeared to hybridize to a chromosome other than chromosome 2, suggesting that the plasmid integrated elsewhere into the genome, possibly by nonhomologous recombination (data not shown). Southern blot analysis revealed that all attempts to integrate the plasmid into the genome were unsuccessful (Fig. 1). Hence, our attempts to disrupt SERA5 in the D10 line made in this study, together with the multiple attempts in both D10 and 3D7 lines described in our previous study (21), suggest an important role for this protein in the erythrocytic cycle.
SERA9, which is located on chromosome 9 and is the only SERA gene not contained in the chromosome 2 locus, has not been previously targeted for deletion. Plasmid pHTK
SERA9 was transfected into D10 parasites and cycled in the same way as described above. Hybridization with the targeting sequence of SERA9 identified a single band of 2.0 kb in D10 corresponding to the endogenous gene. In the transfected D10-SERA9/C0 line that had not undergone negative selection, the bands at 2.0 kb and 9.6 kb indicate the presence of plasmid pTK
SERA9. However, in the D10-
SERA9 parasites, which had undergone negative selection pressure with ganciclovir, a single band at 5.4 kb indicated that the plasmid had integrated into the SERA9 gene via the expected event (Fig. 1).
In summary, we have successfully disrupted SERA1, SERA4, and SERA9, while SERA5 remained refractory to deletion. As we previously deleted SERA2 and SERA3, it is apparent that five of the six members of the serine-type SERA family are not essential to in vitro blood-stage growth, while SERA5 appears to be important for the maintenance of this life cycle stage.
Deletion of SERA4 leads to a transcriptional up-regulation of SERA5.
The phenotype of SERA4 knockout parasites was analyzed in detail, considering that it is the second-highest-expressing member of the SERA family (after SERA5) in the erythrocytic cycle (2, 20). Late-stage parental D10 and D10-
SERA4 parasite proteins were separated on a 10% polyacrylamide gel and subjected to Western blot analysis to confirm that the disruption of the SERA4 gene had specifically affected production of the SERA4 protein (Fig. 2A). Duplicate membranes were probed with either rabbit anti-SERA4 polyclonal antibodies or rabbit anti-Hsp70 polyclonal antibodies as a protein loading control. Full-length SERA4 was detected at approximately 100 kDa in D10 and D10-
SERA4/C0 parasites, a line containing an episomal form of the plasmid, but not in the two cloned D10-
SERA4 lines. Integration of plasmid pHTK
SERA4 into the SERA4 gene was expected to result in a truncated SERA4 gene predicted to encode the first 157 amino acids of the 963-amino-acid wild-type SERA4 protein. A truncated form of SERA4 was not detected by Western blot analysis of D10-
SERA4 parasites (data not shown).
![]() View larger version (57K): [in a new window] |
FIG. 2. Western blot analysis of SERA4-null parasites. (A) Western blots of total material from synchronized late-blood-stage parasites comparing the parental D10 wild type with two clones of D10- SERA4 were probed with either rabbit anti-Hsp70 polyclonal or rabbit anti-SERA4 polyclonal antibody. The antisera used to probe each membrane are indicated at the bottom of the gel. Molecular mass markers (kDa) are shown to the left of the membrane. (B and C) Western blots of total material from synchronized late-blood-stage parasites from wild-type D10 (B) or D10- SERA4 (clone 1) (C) probed with rabbit anti-SERA1 to anti-SERA9 polyclonal antibodies specific to the N-terminal region of each SERA protein (lanes 1 to 9, respectively). The antiserum used to probe each membrane strip is indicated at the top of the lane. The strip probed with the prebleed pool is indicated by P. Molecular mass markers (kDa) are shown to the left of each membrane.
|
SERA4 (Fig. 2C) parasite proteins extracted from late-stage parasites. Anti-SERA4, anti-SERA5, and anti-SERA6 antibodies reacted strongly with the expected species in the parental D10 extracts. However, as expected, only SERA5 and SERA6 were detected in the D10-
SERA4 knockout line. SERA proteins 1, 2, 3, 7, 8, and 9 were not obvious in either D10 or D10-
SERA4 parasites, suggesting that they remain at low levels of expression in both circumstances.
In vitro growth of the D10 parental line was compared with the D10-
SERA4 parasite line to determine whether the deletion of SERA4 had any deleterious effect on the ability of these parasites to propagate in blood-stage cultures (Fig. 3). Parasites were synchronized at ring stages, and growth and stage distribution were followed thereafter over a 28-day time period by examination of thin blood smears every 2 days. No detectable difference in growth rates was observed between the parental line and the mutant parasites at any stage over this period.
![]() View larger version (13K): [in a new window] |
FIG. 3. Disruption of SERA4 does not affect in vitro parasite growth. Shown are data from an extended growth rate assay comparing D10 wild-type parasites with the D10- SERA4 (D10- S4) parasite line over a 28-day period. Parasites were maintained below 5%. Error bars represent the standard deviations of mean parasitemias from three individual counts.
|
SERA4 line using a combination of microarray and quantitative PCR analyses. Competitive microarray analysis was performed on two occasions using 6,912 70-mer oligonucleotides printed in duplicate on polylysine-coated microscope slides (as described in reference 4). Interestingly, in both arrays, SERA5 was up-regulated by approximately threefold in the D10-
SERA4 knockout line in comparison to the parental D10 line (see Fig. S1 in the supplemental material). This was the only gene that showed consistent and obvious transcriptional changes using this method. In addition to this, the coregulated schizont-stage control genes MSP-2, MSP-4, and MSP-5 remained unchanged in the D10-
SERA4 parasite line.
To independently confirm an up-regulation of SERA5 mRNA in SERA4 knockout parasites, we used quantitative RT-PCR. To account for any differences in the synchrony of the parasite lines, we tested RNA extracted from parasites isolated from time courses at 4-h intervals. Expression levels were normalized to the expression of MSP-1, a gene that has a very similar transcriptional profile (Fig. 4). As expected, the level of SERA4 mRNA transcripts was dramatically reduced in the D10-
SERA4 parasite line. We consistently observed a twofold increase in SERA5 transcript in SERA4 knockout parasites, while mRNA levels of other SERA genes did not change. These data suggest that a compensatory up-regulation of SERA5 is required in the absence of SERA4. As mRNA and protein levels of SERA appear to be roughly concordant, we expect that such a transcriptional change leads to increased SERA5 protein in SERA4 knockout parasites. However, any change in protein levels was not sufficiently large to be detected by semiquantitative Western blotting (data not shown).
![]() View larger version (30K): [in a new window] |
FIG. 4. Higher levels of SERA5 mRNA are present in SERA4-null parasites. (A and B) SERA gene transcript expression levels in late-blood-stage D10 wild-type (A) andD10- SERA4 (B) parasites measured by quantitative RT-PCR. Parasites in each parasite line were collected at 30 h postsynchronization and subsequent 4-h intervals for 20 h. All expression levels (mRNA) were normalized to the coregulated gene MSP-1. The peak mRNA expression level observed for SERA5 is indicated on the panels to illustrate the difference between D10 and D10 S4. (C) Change (n-fold) in gene expression levels between D10 and D10- SERA4 parasite lines. This was calculated by dividing the mean expression levels of the individual values shown in B by those shown in A.
|
Beginning with HMM profiles generated from P. falciparum SERA sequences, we sought to locate and extract conserved domains of SERA genes from the genomes of sequenced apicomplexan species. By using this approach and by restricting ourselves to species for which whole-genome sequencing has been undertaken, we aimed to generate a map of SERA loci in the surveyed genomes that was as complete as possible. This map of the genomic arrangement (Fig. 5) together with the phylogeny (Fig. 6) constructed from the retrieved sequences shed light on the evolution of the SERA family and lend support to the hypothesis of a single function for serine-type SERAs.
![]() View larger version (34K): [in a new window] |
FIG. 5. SERA loci across the Plasmodium and Theileria genera. HMM hits for the protease and C-terminal SERA domains are indicated schematically on the contigs on which they were found. Breaks in species tracks indicate discontinuities at contig ends. Where they exist, gene annotations have been indicated. When HMM matches existed on multiple short contigs in a species, the contigs were ordered according to orthology relationships with related species inferred from the phylogeny. SERA genes that have been successfully knocked out in P. falciparum are marked in red, while those refractory to deletion are in black. Where orthology is clearly discernible from the phylogeny, it is marked by vertical blocks that join species. Orthologies of P. gallinaceum and T. annulata genes are ambiguous due to a postspecies duplication. This ambiguity is represented by a bifurcation in the orthology block. The phylogeny of Cys SERA genes indicates a duplication of SERA8 present only in P. gallinaceum as well as a duplication resulting in SERA6 and SERA7 in Plasmodium spp. (excluding P. gallinaceum). Ambiguities arising from these duplications are represented by bifurcations in the orthology blocks. Where it is clearly present, the orthology of flanking sequence is also noted. Protease domains marked in green do not contain active-site mutations, whereas protease domains marked in blue contain the cysteine-to-serine active-site mutation. The presence of secondary-site mutations has been marked below the corresponding protease domain. RT-PCR expression levels (this study) and protein abundance as a percentages of SERA5 are shown above the P. falciparum (3D7) track (16, 17). Similarly, relative copy number expression measurements for P. vivax as a percentage of P. vivax SERA4 are shown above the P. vivax track (22). KO, knockout.
|
![]() View larger version (68K): [in a new window] |
FIG. 6. SERA sequences belong to one of two major phyletic groups that separate according to their catalytic residues. Colored segments that make up the inner disc delineate active-site and secondary-site mutations in SERA sequences. Colored segments that are part of the outer ring mark the clade to which sequences belong. Thick lines in the phylogeny indicate branches shared by MCMC, maximum-likelihood, and neighbor-joining topologies. Thin lines indicate branches supported by MCMC and one other method. Dashed lines indicate branches existing only in the MCMC topology. MCMC branch probabilities are given on the appropriate branches where they are less than 1. Members of the genus Theileria appear to have a single-copy SERA, the only observation to date of a SERA ortholog outside the genus Plasmodium. Radiation in SERA sequences containing the serine active-site mutation tends to be species specific, whereas radiation of cysteine SERA sequences has occurred prior to speciation. Branch lengths for SERA6 and SERA7 tend to be shorter than those of serine-type SERAs and those of SERA8.
|
By examining the consensus tree constructed from MCMC samples not containing the branch dividing the cysteine- and serine-type SERA subgroups, we sought to discover whether there was a plausible evolutionary scenario in which the cysteine SERA subgroup was not monophyletic. The consensus tree constructed from these data maintained the broad grouping of sequences but moved the SERA6/SERA7 subtree to a point in the tree that would require a second Ser
Cys mutation to correct the original Cys
Ser mutation and would also imply that SERA6 and SERA7 have arisen relatively recently. Because the Cys-Ser division is supported by multiple methods of tree reconstruction and is more parsimonious than the alternative hypothesis, it is reasonable to assume monophyly of the cysteine-type SERAs.
Various numbers of serine-type SERA genes are observed in different species, ranging from none in P. gallinaceum to nine in P. vivax. The serine-type subgroup of the phylogeny divides neatly into three clades: P. falciparum/P. reichenowi, P. vivax/P. knowlesi, and P.yoelii/P.berghei/P.chabaudi (Fig. 6), and it therefore appears that the duplications that gave rise to multiple copies of the serine-type SERAs have occurred after the speciation events that gave rise to these clades.
The pattern of clade- and species-specific duplication among serine-type SERAs suggests that the expansion of the serine-type SERAs has been relatively recent and followed an ancient cysteine-type-to-serine-type SERA split in an ancestral Plasmodium species. The species-specific expansion of the subfamily, especially in the case of P. vivax and P. falciparum, which infect the same host, supports the hypothesis that a single function for the subfamily exists and that duplication is being driven by host selective pressure.
|
|
|---|
The presence of SERA5 in the parasitophorous vacuole and demonstrated (relatively weak) in vitro catalytic activity are suggestive of a role in schizont rupture or cleavage of invasion ligands. However, there remains no proof of this, and indeed, other possible roles at this stage can be envisaged, including possible noncatalytic regulatory functions.
While an important blood-stage role for the serine-type SERAs appears to be likely, the number of different functions that these proteins perform appears to be far fewer than the number of genes present. The comprehensive gene identification and phylogenetic analysis performed in this study confirm that the serine-type SERAs form a single monophyletic group arising from a single ancient duplication and active-site mutation of a cysteine-type SERA. This supports the notion that redundancy is possible between the serine-type SERAs, and indeed, it is possible that there is only one blood-stage role for all members of this group. In addition, the ability to disrupt all except the most strongly expressed of the serine-type SERAs indicates that these more weakly expressed members do not perform independent and distinct functions within the erythrocytic cycle. Moreover, the apparent compensatory increase in SERA5 expression in response to the disruption of SERA4, already discussed above in relation to a blood-stage function, is also consistent with functional redundancy among this group.
If all serine-type SERAs perform the same function, why is SERA5 refractory to genetic deletion? The most likely explanation for this is that the required compensatory change in the expression of another of the serine-type SERAs is too great to be readily achieved in culture. Formal proof of this and proof of the hypothesis of redundant roles for the serine-type SERAs require functional complementation experiments where the SERA5 coding sequence is replaced with that of another serine-type SERA. As SERA genes are large, this remains technically challenging but is theoretically achievable.
Given that we have now seen that the serine-type SERAs form a single monophyletic group, we propose that only one function is likely for members of the serine-type SERAs. This role need not be restricted to the blood stages, and indeed, mRNA profiling suggests that some members of the serine-type SERAs are expressed at other stages (19). With respect to the cysteine-type SERAs, phylogenetic analysis suggests that as many as three orthologous groups represented by P. falciparum members SERA6, SERA7, and SERA8 (Fig. 5 and 6). Some functional data to support this separation exist. The rodent malaria version of the "SERA8-like" SERAs functions in oocyst rupture (1), while "SERA6-like" members are most likely to function in the erythrocytic stage, given their strong expression and apparently essential role at this stage (20).
SERAs that function in the blood stage remain attractive drug targets, and indeed, SERA5 has even progressed considerably as a vaccine candidate. There is now a strong likelihood that interference with serine-type SERA function will interrupt normal erythrocytic development and promote parasite clearance from infected individuals. SERA5 may indeed be obligatory for blood-stage growth. However, if other family members such as SERA4 can compensate if expressed highly enough, it is important that interventions targeting SERA5 should be designed to also block the other serine-type SERAs.
This work was supported by the National Health and Medical Research Council of Australia and the Wellcome Trust, United Kingdom. J.E.M. and T.S. are the recipients of NHMRC Dora Lush postgraduate scholarships, S.K.M. was the recipient of a Melbourne Research scholarship (University of Melbourne), and B.S.C. is an International Research Scholar of the Howard Hughes Medical Institute.
Published ahead of print on 24 September 2007. ![]()
Supplemental material for this article may be found at http://iai.asm.org/. ![]()
J.E.M. and S.K.M. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»