Immune tolerance caused by repeated P. falciparum infection against SE36 malaria vaccine candidate antigen and the resulting limited polymorphism

The call for second generation malaria vaccines needs not only the identification of novel candidate antigens or adjuvants but also a better understanding of immune responses and the underlying protective processes. Plasmodium parasites have evolved a range of strategies to manipulate the host immune system to guarantee survival and establish parasitism. These immune evasion strategies hamper efforts to develop effective malaria vaccines. In the case of a malaria vaccine targeting the N-terminal domain of P. falciparum serine repeat antigen 5 (SE36), now in clinical trials, we observed reduced responsiveness (lowered immunogenicity) which may be attributed to immune tolerance/immune suppression. Here, immunogenicity data and insights into the immune responses to SE36 antigen from epidemiological studies and clinical trials are summarized. Documenting these observations is important to help identify gaps for SE36 continued development and engender hope that highly effective blood-stage/multi-stage vaccines can be achieved.

Analysis of genome instability and implications for the consequent phenotype in Plasmodium falciparum containing mutated MSH2-1 (P513T)

Malarial parasites exhibit extensive genomic plasticity, which induces the antigen diversification and the development of antimalarial drug resistance. Only a few studies have examined the genome maintenance mechanisms of parasites. The study aimed at elucidating the impact of a mutation in a DNA mismatch repair gene on genome stability by maintaining the mutant and wild-type parasites through serial in vitro cultures for approximately 400 days and analysing the subsequent spontaneous mutations. A P513T mutant of the DNA mismatch repair protein PfMSH2-1 from Plasmodium falciparum 3D7 was created. The mutation did not influence the base substitution rate but significantly increased the insertion/deletion (indel) mutation rate in short tandem repeats (STRs) and minisatellite loci. STR mutability was affected by allele size, genomic category and certain repeat motifs. In the mutants, significant telomere healing and homologous recombination at chromosomal ends caused extensive gene loss and generation of chimeric genes, resulting in large-scale chromosomal alteration. Additionally, the mutant showed increased tolerance to N-methyl-N'-nitro-N-nitrosoguanidine, suggesting that PfMSH2-1 was involved in recognizing DNA methylation damage. This work provides valuable insights into the role of PfMSH2-1 in genome stability and demonstrates that the genomic destabilization caused by its dysfunction may lead to antigen diversification.

Plasmodium falciparum infection coinciding with the malaria vaccine candidate BK-SE36 administration interferes with the immune responses in Burkinabe children

Background

A vaccine targeting the erythrocyte stages of Plasmodium falciparum could play a role in preventing clinical disease. BK-SE36 is a promising malaria vaccine candidate that has shown a good safety profile and immunological responses during field evaluations. It was observed that repeated natural infections could result in immune tolerance against SE36 molecule.

Methods

The primary trial was conducted to assess the safety and immunogenicity of the BK-SE36 in two cohorts of children aged 25-60 months (Cohort 1) and 12-24 months (Cohort 2). Immunization was at full dose (1.0 mL) administered at 0, 1, and 6 months. Blood samples were collected before each vaccination for immunological assessments and detection of Plasmodium falciparum infection by microscopy. Blood samples were further collected one month post each vaccination to evaluate immunogenicity.

Results

Of seventy-two (72) subjects that have received BK-SE36 vaccination, 71 had available blood smears during vaccination days. One month post Dose 2, the geometric mean of SE36 antibodies was 263.2 (95% CI: 178.9-387.1) in uninfected individuals compared to 77.1 (95% CI: 47.3-125.7) in infected participants. The same trend was observed one-month post booster dose. Participants uninfected at the time of booster vaccination had significantly higher GMTs compared to those who were infected (424.1 (95% CI: 301.9-595.8) vs. 92.8 (95% CI: 34.9-246.6), p = 0.002. There was a 14.3 (95% CI: 9.7-21.1) and 2.4 (95% CI: 1.3-4.4) fold-change, respectively, in uninfected and infected participants between one-month post Dose 2 and booster. The difference was statistically significant (p < 0.001).

Conclusion

Concomitant infection by P. falciparum during BK-SE36 vaccine candidate administration is associated with reduced humoral responses. However, it is to be noted that the BK-SE36 primary trial was not designed to investigate the influence of concomitant infection on vaccine-induced immune response and should be interpreted cautiously.

Trial registration

WHO ICTRP, PACTR201411000934120.

African-specific polymorphisms in Plasmodium falciparum serine repeat antigen 5 in Uganda and Burkina Faso clinical samples do not interfere with antibody response to BK-SE36 vaccination

BK-SE36, based on Plasmodium falciparum serine repeat antigen 5 (SERA5), is a blood-stage malaria vaccine candidate currently being evaluated in clinical trials. Phase 1 trials in Uganda and Burkina Faso have demonstrated promising safety and immunogenicity profiles. However, the genetic diversity of sera5 in Africa and the role of allele/variant-specific immunity remain a major concern. Here, sequence analyses were done on 226 strains collected from the two clinical trial/follow-up studies and 88 strains from two cross-sectional studies in Africa. Compared to other highly polymorphic vaccine candidate antigens, polymorphisms in sera5 were largely confined to the repeat regions of the gene. Results also confirmed a SERA5 consensus sequence with African-specific polymorphisms. Mismatches with the vaccine-type SE36 (BK-SE36) in the octamer repeat, serine repeat, and flanking regions, and single-nucleotide polymorphisms in non-repeat regions could compromise vaccine response and efficacy. However, the haplotype diversity of SERA5 was similar between vaccinated and control participants. There was no marked bias or difference in the patterns of distribution of the SE36 haplotype and no statistically significant genetic differentiation among parasites infecting BK-SE36 vaccinees and controls. Results indicate that BK-SE36 does not elicit an allele-specific immune response.

Safety and immunogenicity of BK-SE36 in a blinded, randomized, controlled, age de-escalating phase Ib clinical trial in Burkinabe children

Background

A blood-stage vaccine targeting the erythrocytic-stages of the malaria parasite Plasmodium falciparum could play a role to protect against clinical disease. Antibodies against the P. falciparum serine repeat antigen 5 (SE47 and SE36 domains) correlate well with the absence of clinical symptoms in sero-epidemiological studies. A previous phase Ib trial of the recombinant SE36 antigen formulated with aluminum hydroxyl gel (BK-SE36) was promising. This is the first time the vaccine candidate was evaluated in young children below 5 years using two vaccination routes.

Methods

Safety and immunogenicity of BK-SE36 was assessed in a double-blind, randomized, controlled, age de-escalating phase Ib trial. Fifty-four Burkinabe children in each age cohort, 25–60 or 12–24 months, were randomized in a 1:1:1 ratio to receive three doses of BK-SE36 either by intramuscular (BK IM) or subcutaneous (BK SC) route on Day 0, Week 4, and 26; or the control vaccine, Synflorix^®via IM route on Day 0, Week 26 (and physiological saline on Week 4). Safety data and samples for immunogenicity analyses were collected at various time-points.

Results

Of 108 subjects, 104 subjects (96.3%) (Cohort 1: 94.4%; Cohort 2: 98.1%) received all three scheduled vaccine doses. Local reactions, mostly mild or of moderate severity, occurred in 99 subjects (91.7%). The proportion of subjects that received three doses without experiencing Grade 3 adverse events was similar across BK-SE36 vaccines and control arms (Cohort 1: 100%, 89%, and 89%; and Cohort 2: 83%, 82%, and 83% for BK IM, BK SC, and control, respectively). BK-SE36 vaccine was immunogenic, inducing more than 2-fold change in antibody titers from pre-vaccination, with no difference between the two vaccination routes. Titers waned before the third dose but in both cohorts titers were boosted 6 months after the first vaccination. The younger cohort had 2-fold and 4-fold higher geometric mean titers compared to the 25- to 60-month-old cohort after 2 and 3 doses of BK-SE36, respectively.

Conclusion

BK-SE36 was well tolerated and immunogenic using either intramuscular or subcutaneous routes, with higher immune response in the younger cohort.

Clinical Trial Registration

https://pactr.samrc.ac.za/TrialDisplay.aspx?TrialID=934, identifier PACTR201411000934120.

Characteristic features of the SERA multigene family in the malaria parasite

Serine repeat antigen (SERA) is conserved among species of the genus Plasmodium. Sera genes form a multigene family and are generally tandemly clustered on a single chromosome. Although all Plasmodium species encode multiple sera genes, the number varies between species. Among species, the members share similar sequences and gene organization. SERA possess a central papain-like cysteine protease domain, however, in some members, the active site cysteine residue is substituted with a serine. Recent studies implicate this gene family in a number of aspects in parasite biology and induction of protective immune response. This review summarizes the current understanding on this important gene family in several Plasmodium species. The Plasmodium falciparum (Pf)-sera family, for example, consists of nine gene members. Unlike other multigene families in Plasmodium species, Pf-sera genes do not exhibit antigenic variation. Pf-sera5 nucleotide diversity is also low. Moreover, although Pf-sera5 is highly transcribed during the blood stage of malaria infection, and a large amount is released into the host blood following schizont rupture, in malaria endemic countries the sero-positive rates for Pf-SERA5 are low, likely due to Pf-SERA5 binding of host proteins to avoid immune recognition. As an antigen, the N-terminal 47 kDa domain of Pf-SERA5 is a promising vaccine candidate currently undergoing clinical trials. Pf-SERA5 and Pf-SERA6, as well as P. berghei (Pb)-SERA3, and Pb-SERA5, have been investigated for their roles in parasite egress. Two P. yoelii SERA, which have a serine residue at the protease active center, are implicated in parasite virulence. Overall, these studies provide insight that during the evolution of the Plasmodium parasite, the sera gene family members have increased by gene duplication, and acquired various functions that enable the parasite to survive and successfully maintain infection in the host.

Phase 1b randomized trial and follow-up study in Uganda of the blood-stage malaria vaccine candidate BK-SE36

Background

Up to now a malaria vaccine remains elusive. The Plasmodium falciparum serine repeat antigen-5 formulated with aluminum hydroxyl gel (BK-SE36) is a blood-stage malaria vaccine candidate that has undergone phase 1a trial in malaria-naive Japanese adults. We have now assessed the safety and immunogenicity of BK-SE36 in a malaria endemic area in Northern Uganda.

Methods

We performed a two-stage, randomized, single-blinded, placebo-controlled phase 1b trial (Current Controlled trials ISRCTN71619711). A computer-generated sequence randomized healthy subjects for 2 subcutaneous injections at 21-day intervals in Stage1 (21–40 year-olds) to 1-mL BK-SE36 (BKSE1.0) (n = 36) or saline (n = 20) and in Stage2 (6–20 year-olds) to BKSE1.0 (n = 33), 0.5-mL BK-SE36 (BKSE0.5) (n = 33), or saline (n = 18). Subjects and laboratory personnel were blinded. Safety and antibody responses 21-days post-second vaccination (Day42) were assessed. Post-trial, to compare the risk of malaria episodes 130–365 days post-second vaccination, Stage2 subjects were age-matched to 50 control individuals.

Results

Nearly all subjects who received BK-SE36 had induration (Stage1, n = 33, 92%; Stage2, n = 63, 96%) as a local adverse event. No serious adverse event related to BK-SE36 was reported. Pre-existing anti-SE36 antibody titers negatively correlated with vaccination-induced antibody response. At Day42, change in antibody titers was significant for seronegative adults (1.95-fold higher than baseline [95% CI, 1.56–2.43], p = 0.004) and 6–10 year-olds (5.71-fold [95% CI, 2.38–13.72], p = 0.002) vaccinated with BKSE1.0. Immunogenicity response to BKSE0.5 was low and not significant (1.55-fold [95% CI, 1.24–1.94], p = 0.75). In the ancillary analysis, cumulative incidence of first malaria episodes with ≥5000 parasites/µL was 7 cases/33 subjects in BKSE1.0 and 10 cases/33 subjects in BKSE0.5 vs. 29 cases/66 subjects in the control group. Risk ratio for BKSE1.0 was 0.48 (95% CI, 0.24–0.98; p = 0.04).

Conclusion

BK-SE36 is safe and immunogenic. The promising potential of BK-SE36, observed in the follow-up study, warrants a double-blind phase 1/2b trial in children under 5 years.

Trial Registration

Controlled-Trials.com ISRCTN71619711 ISRCTN71619711

Plasmodium cynomolgi genome sequences provide insight into Plasmodium vivax and the monkey malaria clade

Plasmodium cynomolgi, a malaria parasite of Asian Old World monkeys, is the sister taxon of Plasmodium vivax, the most prevalent human malaria species outside Africa. Since P. cynomolgi shares many phenotypic, biologic and genetic characteristics of P. vivax, we generated draft genome sequences of three P. cynomolgi strains and performed comparative genomic analysis between them and P. vivax, as well as a third previously sequenced simian parasite, Plasmodium knowlesi. Here we show that genomes of the monkey malaria clade can be characterized by CNVs in multigene families involved in evasion of the human immune system and invasion of host erythrocytes. We identify genome-wide SNPs, microsatellites, and CNVs in the P. cynomolgi genome, providing a map of genetic variation for mapping parasite traits and studying parasite populations. The P. cynomolgi genome is a critical step in developing a model system for P. vivax research, and to counteract the neglect of P. vivax.

Geographic differentiation of polymorphism in the Plasmodium falciparum malaria vaccine candidate gene SERA5

SERA5 is regarded as a promising malaria vaccine candidate of the most virulent human malaria parasite Plasmodium falciparum. SERA5 is a 120 kDa abundantly expressed blood-stage protein containing a papain-like protease. Since substantial polymorphism in blood-stage vaccine candidates may potentially limit their efficacy, it is imperative to fully investigate polymorphism of the SERA5 gene (sera5). In this study, we performed evolutionary and population genetic analysis of sera5. The level of inter-species divergence (kS=0.076) between P. falciparum and Plasmodium reichenowi, a closely related chimpanzee malaria parasite is comparable to that of housekeeping protein genes. A signature of purifying selection was detected in the proenzyme and enzyme domains. Analysis of 445 near full-length P. falciparum sera5 sequences from nine countries in Africa, Southeast Asia, Oceania and South America revealed extensive variations in the number of octamer repeat (OR) and serine repeat (SR) regions as well as substantial level of single nucleotide polymorphism (SNP) in non-repeat regions (2562 bp). Remarkably, a 14 amino acid sequence of SERA5 (amino acids 59-72) that is known to be the in vitro target of parasite growth inhibitory antibodies was found to be perfectly conserved in all 445 worldwide isolates of P. falciparum evaluated. Unlike other major vaccine target antigen genes such as merozoite surface protein-1, apical membrane antigen-1 or circumsporozoite protein, no strong evidence for positive selection was detected for SNPs in the non-repeat regions of sera5. A biased geographical distribution was observed in SNPs as well as in the haplotypes of the sera5 OR and SR regions. In Africa, OR- and SR-haplotypes with low frequency (<5%) and SNPs with minor allele frequency (<5%) were abundant and were mostly continent-specific. Consistently, significant genetic differentiation, assessed by the Wright's fixation index (Fst) of inter-population variance in allele frequencies, was detected for SNPs and both OR- and SR-haplotypes among almost all parasite populations. The exception was parasite populations between Tanzania and Ghana, suggesting frequent gene flow in Africa. The present study points to the importance of investigating whether biased geographical distribution for SNPs and repeat variants in the OR and SR regions affect the reactivity of human serum antibodies to variants.

Serologic markers in relation to parasite exposure history help to estimate transmission dynamics of Plasmodium vivax

Plasmodium vivax infection has been gaining attention because of its re-emergence in several parts of the world. Southeastern Turkey is one of the places in which persistent focal malaria caused exclusively by P. vivax parasites occurs. Although control and elimination studies have been underway for many years, no detailed study has been conducted to understand the mechanisms underlying the ineffective control of malaria in this region. Here, for the first time, using serologic markers we try to extract as much information as possible in this region to get a glimpse of P. vivax transmission. We conducted a sero-immunological study, evaluating antibody responses of individuals living in Sanliurfa to four different P. vivax antigens; three blood-stage antigens (PvMSP1₁₉, PvAMA1-ecto, and PvSERA4) and one pre-erythrocytic stage antigen (PvCSP). The results suggest that a prior history of malaria infection and age can be determining factors for the levels and sustainability of naturally acquired antibodies. Significantly higher antibody responses to all the studied antigens were observed in blood smear-negative individuals with a prior history of malaria infection. Moreover, these individuals were significantly older than blood smear-negative individuals with no prior history of infection. These data from an area of sole P. vivax-endemic region may have important implications for the global malaria control/elimination programs and vaccine design.

Clues to evolution of the SERA multigene family in 18 Plasmodium species

SERA gene sequences were newly determined from 11 primate Plasmodium species including two human parasites, P. ovale and P. malariae, and the evolutionary history of SERA genes was analyzed together with 7 known species. All have one each of Group I to III cysteine-type SERA genes and varying number of Group IV serine-type SERA genes in tandem cluster. Notably, Group IV SERA genes were ascertained in all mammalian parasite lineages; and in two primate parasite lineages gene events such as duplication, truncation, fragmentation and gene loss occurred at high frequency in a manner that mimics the birth-and-death evolution model. Transcription profile of individual SERA genes varied greatly among rodent and monkey parasites. Results support the lineage-specific evolution of the Plasmodium SERA gene family. These findings provide further impetus for studies that could clarify/provide proof-of-concept that duplications of SERA genes were associated with the parasites' expansion of host range and the evolutionary conundrums of multigene families in Plasmodium.

Limited polymorphism of the Plasmodium vivax merozoite surface protein 1 gene in isolates from Turkey

The 200-kD merozoite surface protein of Plasmodium vivax (PvMSP-1) is one of the leading vaccine candidates against P. vivax malaria. However, the gene encoding PvMSP-1 (pvmsp1) is highly polymorphic and is a major obstacle to effective vaccine development. To further understand polymorphism in pvmsp1, we obtained 30 full-length pvmsp1 sequences from southeastern Turkey. Comparative analysis of sequences from Turkey and other areas showed substantially limited polymorphism. Substitutions were found at 280 and 162 amino acid sites in samples from other regions and those from Turkey, respectively. Eight substitutions were unique to Turkey. In one of them, D/E at position 1706 in the C-terminal 19-kD region, the K/E change at 1709 was the only polymorphism previously known. Limited diversity was also observed in microsatellites. Data suggest a recent population bottleneck in Turkey that may have obscured a signature for balancing selection in the C-terminal 42-kD region, which was otherwise detectable in other areas.

Concatenated mitochondrial DNA of the coccidian parasite Eimeria tenella

Apicomplexan parasites of the genus Plasmodium, pathogens causing malaria, and the genera Babesia and Theileria, aetiological agents of piroplasmosis, are closely related. However, their mitochondrial (mt) genome structures are highly divergent: Plasmodium has a concatemer of 6-kb unit and Babesia/Theileria a monomer of 6.6- to 8.2-kb with terminal inverted repeats. Fragmentation of ribosomal RNA (rRNA) genes and gene arrangements are remarkably distinctive. To elucidate the evolutionary origin of this structural divergence, we determined the mt genome of Eimeria tenella, pathogens of coccidiosis in domestic fowls. Analysis revealed that E. tenella mt genome was concatemeric with similar protein-coding genes and rRNA gene fragments to Plasmodium. Copy number was 50-fold of the nuclear genome. Evolution of structural divergence in the apicomplexan mt genomes is discussed.

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites

Background

The 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. Plasmodium vivax, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to P. vivax-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity.

Results

We performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (msp1) from P. vivax and nine P. vivax-related simian malaria parasites. The inferred phylogenetic tree of msp1 significantly differed from that of the mitochondrial genome, with a striking displacement of P. vivax from a position close to P. cynomolgi in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to P. inui and P. hylobati and the other leading to P. vivax, P. fieldi and P. cynomolgi. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of msp1 polymorphisms between P. vivax, P. inui and P. cynomolgi revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on msp1.

Conclusions

The present results indicate that the msp1 locus of P. vivax and related parasite species has lineage-specific unique evolutionary history with positive selection. P. vivax and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as msp1.

Divergence of the mitochondrial genome structure in the apicomplexan parasites, Babesia and Theileria

Mitochondrial (mt) genomes from diverse phylogenetic groups vary considerably in size, structure, and organization. The genus Plasmodium, causative agent of malaria, of the phylum Apicomplexa, has the smallest mt genome in the form of a circular and/or tandemly repeated linear element of 6 kb, encoding only three protein genes (cox1, cox3, and cob). The closely related genera Babesia and Theileria also have small mt genomes (6.6 kb) that are monomeric linear with an organization distinct from Plasmodium. To elucidate the structural divergence and evolution of mt genomes between Babesia/Theileria and Plasmodium, we determined five new sequences from Babesia bigemina, B. caballi, B. gibsoni, Theileria orientalis, and T. equi. Together with previously reported sequences of B. bovis, T. annulata, and T. parva, all eight Babesia and Theileria mt genomes are linear molecules with terminal inverted repeats (TIRs) on both ends containing three protein-coding genes (cox1, cox3, and cob) and six large subunit (LSU) ribosomal RNA (rRNA) gene fragments. The organization and transcriptional direction of protein-coding genes and the rRNA gene fragments were completely conserved in the four Babesia species. In contrast, notable variation occurred in the four Theileria species. Although the genome structures of T. annulata and T. parva were nearly identical to those of Babesia, an inversion in the 3-kb central region was found in T. orientalis. Moreover, the T. equi mt genome is the largest (8.2 kb) and most divergent with unusually long TIR sequences, in which cox3 and two LSU rRNA gene fragments are located. The T. equi mt genome showed little synteny to the other species. These results suggest that the Theileria mt genome is highly diverse with lineage-specific evolution in two Theileria species: genome inversion in T. orientalis and gene-embedded long TIR in T. equi.

Identification of Plasmodium malariae, a human malaria parasite, in imported chimpanzees

It is widely believed that human malaria parasites infect only man as a natural host. However, earlier morphological observations suggest that great apes are likely to be natural reservoirs as well. To identify malaria parasites in great apes, we screened 60 chimpanzees imported into Japan. Using the sequences of small subunit rRNA and the mitochondrial genome, we identified infection of Plasmodium malariae, a human malaria parasite, in two chimpanzees that were imported about thirty years ago. The chimpanzees have been asymptomatic to the present. In Japan, indigenous malaria disappeared more than fifty years ago; and thus, it is most likely inferred that the chimpanzees were infected in Africa, and P. malariae isolates were brought into Japan from Africa with their hosts, suggesting persistence of parasites at low level for thirty years. Such a long term latent infection is a unique feature of P. malariae infection in humans. To our knowledge, this is the first to report P. malariae infection in chimpanzees and a human malaria parasite from nonhuman primates imported to a nonendemic country.

Recent independent evolution of msp1 polymorphism in Plasmodium vivax and related simian malaria parasites

The Plasmodium MSP-1 is a promising malaria vaccine candidate. However, the highly polymorphic nature of the MSP-1 gene (msp1) presents a potential obstacle for effective vaccine development. To investigate the evolutionary history of msp1 polymorphism in P. vivax, we construct phylogenetic trees of msp1 from P. vivax and related monkey malaria parasite species. All P. vivax msp1 alleles cluster in the P. vivax lineage and are not distributed among other species. Similarly, all P. cynomolgi msp1 alleles cluster in the P. cynomolgi lineage. This suggests that, in contrast to presumed ancient origin of P. falciparum msp1 polymorphism, the origin of P. vivax msp1 polymorphism is relatively recent. We observed positive selection in the P. vivax lineage but not in P. cynomolgi. Also, positive selection acts on different regions of msp1 in P. vivax and P. falciparum. This study shows that the evolutionary history of msp1 differs greatly among parasite lineages.

Phylogeny and evolution of the SERA multigene family in the genus Plasmodium

The serine repeat antigen gene family of Plasmodium falciparum (Pf-SERA) consists of nine gene members. By sequence similarity search, 45 genes were identified to be homologous to the Pf-SERA genes in the ongoing seven Plasmodium genome sequencing project databases for the species: P. reichenowi, P. vivax, P. knowlesi, P. yoelii, P. berghei, P. chabaudi, and P. gallinaceum. In combination with additional PCR-based sequencing, we found that almost all SERA genes in each species were aligned in a tandem cluster and sandwiched between two conserved hypothetical protein genes, except for P. reichenowi, which could not be confirmed. The minimum and maximum numbers of clustered genes were 2 and 12 for P. gallinaceum and P. vivax, respectively. The best tree of the maximum likelihood analysis demonstrated that all Plasmodium SERA homologues, except for SERA1 of P. gallinaceum (Pg-SERA1), can be classified into four groups, represented by Pf-SERA5, Pf-SERA6, Pf-SERA7, and Pf-SERA8. Genes in the Pf-SERA8 group, although highly divergent and distantly related to the sequences of other groups, were not pseudogenes. P. berghei SERA5, the counterpart of Pf-SERA8, was expressed in the mosquito stage. P. gallinaceum lacks the orthologues to Pf-SERA5, Pf-SERA6, and Pf-SERA7, suggesting that P. gallinaceum diverged from a common ancestor of all eight Plasmodium species examined before gene duplication(s) occurred to generate these paralogous groups. Here, we reveal an evolutionary trail of SERA gene cluster in the genus Plasmodium and discuss a phylogeny of Plasmodium species from the viewpoint of the evolution of a multigene family.

Allelic dimorphism-associated restriction of recombination in Plasmodium falciparum msp1

Allelic dimorphism is a characteristic feature of the Plasmodium falciparum msp1 gene encoding the merozoite surface protein 1, a strong malaria vaccine candidate. Meiotic recombination is a major mechanism for the generation of msp1 allelic diversity. Potential recombination sites have previously been mapped to specific regions within msp1 (a 5' 1-kb region and a 3' 0.4-kb region) with no evidence for recombination events in a central 3.5-kb region. However, evidence for the lack of recombination events is circumstantial and inconclusive because the number of msp1 sequences analysed is limited, and the frequency of recombination events has not been addressed previously in a high transmission area, where the frequency of meiotic recombination is expected to be high. In the present study, we have mapped potential allelic recombination sites in 34 full-length msp1 sequences, including 24 new sequences, from various geographic origins. We also investigated recombination events in blocks 6 to 16 by population genetic analysis of P. falciparum populations in Tanzania, where malaria transmission is intense. The results clearly provide no evidence of recombination events occurring between the two major msp1 allelic types, K1-type and Mad20-type, in the central region, but do show recombination events occurring throughout the entire gene within sequences of the Mad20-type. Thus, the present study indicates that allelic dimorphism of msp1 greatly affects inter-allelic recombination events, highlighting a unique feature of allelic diversity of P. falciparum msp1.

Pathological role of Toll-like receptor signaling in cerebral malaria

Toll-like receptors (TLRs) recognize malaria parasites or their metabolites; however, their physiological roles in malaria infection in vivo are not fully understood. Here, we show that myeloid differentiation primary response gene 88 (MyD88)-dependent TLR signaling mediates brain pathogenesis of severe malaria infection, namely cerebral malaria (CM). A significant number of MyD88-, but not TIR domain containing adaptor-inducing IFN-beta (TRIF)-deficient or wild-type (WT) mice survived CM caused by Plasmodium berghei ANKA (PbA) infection. Although systemic parasitemia was comparable, sequestration of parasite and hemozoin load in the brain blood vessels was significantly lower in MyD88-deficient mice compared with those in TRIF-deficient or WT mice. Moreover, brain-specific pathological changes were associated with MyD88-dependent infiltration of CD8+, CCR5+ T cells and CD11c+ dendritic cells, including CD11c+, NK1.1+ and B220+ cells, and up-regulation of genes such as Granzyme B, Lipocalin 2, Ccl3 and Ccr5. Further studies using mice lacking various TLRs suggest that TLR2 and TLR9, but not TLR4, 5 and 7, were involved in CM. These results strongly suggest that TLR2- and/or TLR9-mediated, MyD88-dependent brain pathogenesis may play a critical role in CM, the lethal complication during PbA infection.

Comparative analysis of the ribosomal components of the hydrogenosome-containing protist, Trichomonas vaginalis

The ribosomes of the amitochondriate but hydrogenosome-containing protist lineage, the trichomonads, have previously been reported to be prokaryotic or primitive eukaryotic, based on evidence that they have a 70S sedimentation coefficient and a small number of proteins, similar to prokaryotic ribosomes. In order to determine whether the components of the trichomonad ribosome indeed differ from those of typical eukaryotic ribosomes, the ribosome of a representative trichomonad, Trichomonas vaginalis, was characterized. The sedimentation coefficient of the T. vaginalis ribosome was smaller than that of Saccharomyces cerevisiae and larger than that of Escherichia coli. Based on two-dimensional PAGE analysis, the number of different ribosomal proteins was estimated to be approximately 80. This number is the same as those obtained for typical eukaryotes (approximately 80) but larger than that of E. coli (approximately 55). N-Terminal amino acid sequencing of 18 protein spots and the complete sequences of 4 ribosomal proteins as deduced from their genes revealed these sequences to display typical eukaryotic features. Phylogenetic analyses of the five ribosomal proteins currently available also clearly confirmed that the T. vaginalis sequences are positioned within a eukaryotic clade. Comparison of deduced secondary structure models of the small and large subunit rRNAs of T. vaginalis with those of other eukaryotes revealed that all helices commonly found in typical eukaryotes are present and conserved in T. vaginalis, while variable regions are shortened or lost. These lines of evidence demonstrate that the T. vaginalis ribosome has no prokaryotic or primitive eukaryotic features but is clearly a typical eukaryotic type.

Root of the Eukaryota Tree as Inferred from Combined Maximum Likelihood Analyses of Multiple Molecular Sequence Data

Extensive studies aiming to establish the structure and root of the Eukaryota tree by phylogenetic analyses of molecular sequences have thus far not resulted in a generally accepted tree. To re-examine the eukaryotic phylogeny using alternative genes, and to obtain a more robust inference for the root of the tree as well as the relationship among major eukaryotic groups, we sequenced the genes encoding isoleucyl-tRNA and valyl-tRNA synthetases, cytosolic-type heat shock protein 90, and the largest subunit of RNA polymerase II from several protists. Combined maximum likelihood analyses of 22 protein-coding genes including the above four genes clearly demonstrated that Diplomonadida and Parabasala shared a common ancestor in the rooted tree of Eukaryota, but only when the fast-evolving sites were excluded from the original data sets. The combined analyses, together with recent findings on the distribution of a fused dihydrofolate reductase-thymidylate synthetase gene, narrowed the possible position of the root of the Eukaryota tree on the branch leading to Opisthokonta or to the common ancestor of Diplomonadida/Parabasala. However, the analyses did not agree with the position of the root located on the common ancestor of Opisthokonta and Amoebozoa, which was argued by Stechmann and Cavalier-Smith [Curr. Biol. 13:R665-666, 2003] based on the presence or absence of a three-gene fusion of the pyrimidine biosynthetic pathway: carbamoyl-phosphate synthetase II, dihydroorotase, and aspartate carbamoyltransferase. The presence of the three-gene fusion recently found in the Cyanidioschyzon merolae (Rhodophyta) genome sequence data supported our analyses against the Stechmann and Cavalier-Smith-rooting in 2003.

Sequence heterogeneity of the small subunit ribosomal RNA genes among blastocystis isolates

Genes encoding small subunit ribosomal RNA (SSUrRNA) of 16 Blastocystis isolates from humans and other animals were amplified by the polymerase chain reaction, and the corresponding fragments were cloned and sequenced. Alignment of these sequences with the previously reported ones indicated the presence of 7 different sequence patterns in the highly variable regions of the small subunit ribosomal RNA. Phylogenetic reconstruction analysis using Proteromonas lacertae as the outgroup clearly demonstrated that the 7 groups with the different sequence patterns are separated to form independent clades, 5 of which consisted of the Blastocystis isolates from both humans (B. hominis) and other animals. The presence of 3 higher order clades was also clearly supported in the phylogenetic tree. However, a relationship among the 4 groups including these 3 higher order clades was not settled with statistical confidence. The remarkable heterogeneity of small subunit ribosomal RNAs among different Blastocystis isolates found in this study confirmed, with sequence-based evidence, that these organisms are genetically highly divergent in spite of their morphological identity. The highly variable small subunit ribosomal RNA regions among the distinct groups will provide useful information for the development of group-specific diagnostic primers.

Mitochondrial-type hsp70 genes of the amitochondriate protists, Giardia intestinalis, Entamoeba histolytica and two microsporidians

Genes encoding putative mitochondrial-type heat shock protein 70 (mit-hsp70) were isolated and sequenced from amitochondriate protists, Giardia intestinalis, Entamoeba histolytica, and two microsporidians, Encephalitozoon hellem and Glugea plecoglossi. The deduced mit-hsp70 sequences were analyzed by sequence alignments and phylogenetic reconstructions. The mit-hsp70 sequence of these four amitochondriate protists were divergent from other mit-hsp70 sequences of mitochondriate eukaryotes. However, all of these sequences were clearly located within a eukaryotic mitochondrial clade in the tree including various type hsp70 sequences, supporting the emerging notion that none of these amitochondriate lineages are primitively amitochodrial, but lost their mitochondria secondarily in their evolutionary past.

Phylogenetic position of Blastocystis hominis and of stramenopiles inferred from multiple molecular sequence data

Blastocystis hominis, a parasite of the human intestine, has recently been positioned within stramenopiles by the small subunit rRNA phylogeny. To further confirm its phylogenetic position using multiple molecular sequence data, we determined the nucleotide sequences putatively encoding small subunit ribosomal RNA, cytosolic-type 70-kDa heat shock protein, translation elongation factor 2, and the non-catalytic 'B' subunit of vacuolar ATPase of B. hominis (HE87-1 strain). Moreover, we determined the translation elongation factor 2 sequence of an apicomplexan parasite, Plasmodium falciparum, that belongs to alveolates. The maximum likelihood analyses of small subunit rRNA and cytosolic-type 70-kDa heat shock protein clearly demonstrated that B. hominis (HE87-1 strain) is positioned within stramenopiles, being congruent with the previous small subunit rRNA analysis, including the sequences of B. hominis (Nand strain) and a Blastocystis isolate from guinea pig. Although no clear resolution among major eukaryotic groups was obtained by the individual phylogenies based on the four molecules analyzed here, a combined analysis of various molecules, including these, clearly indicated that Blastocystis/stramenopiles are the closest relatives of alveolates.

The phylogenetic position of the pelobiont Mastigamoeba balamuthi based on sequences of rDNA and translation elongation factors EF-1alpha and EF-2

The taxonomic position and phylogenetic relationships of the Pelobionta, an amitochondriate amoeboflagellate group, are not yet completely settled. To provide more information, we obtained sequences for the large subunit rDNA gene, the gene for translation elongation factor 1alpha, and for a large part of the gene encoding translation elongation factor 2 from a representative of this group, Mastigamoeba balamuthi (formerly Phreatamoeba balamuthi). The gene for the large subunit rDNA was unusually large compared to those of other protists, a phenomenon that had previously been observed for the gene encoding the small subunit rDNA. Phylogenetic reconstruction using a maximum likelihood method was performed with these sequences, as well as the gene encoding the small subunit rDNA. When evaluated individually, the M. balamuthi genes for the small and large subunit rDNAs and elongation factor 1alpha had a most recent common ancestor with either the Mycetozoa (slime molds) or with Entamoeba histolytica. A clade formed by M. balamuthi, E. histolytica, and Mycetozoa was not rejected statistically for any of the sequences. A combined maximum likelihood analysis using 3,935 positions from all molecules suggested that these three taxonomic units form a robust clade. We were unable to resolve the closest group to this clade using the combined analysis. These findings support the notion, which had previously been proposed primarily on cytological evidence, that both M. balamuthi and E. histolytica are closely related to the Mycetozoa and that these three together represent a major eukaryotic lineage.