Author Correction: Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Evolutionary analysis of the most polymorphic gene family in falciparum malaria

The var gene family of the human malaria parasite Plasmodium falciparum encode proteins that are crucial determinants of both pathogenesis and immune evasion and are highly polymorphic. Here we have assembled nearly complete var gene repertoires from 2398 field isolates and analysed a normalised set of 714 from across 12 countries. This therefore represents the first large scale attempt to catalogue the worldwide distribution of var gene sequences We confirm the extreme polymorphism of this gene family but also demonstrate an unexpected level of sequence sharing both within and between continents. We show that this is likely due to both the remnants of selective sweeps as well as a worrying degree of recent gene flow across continents with implications for the spread of drug resistance. We also address the evolution of the var repertoire with respect to the ancestral genes within the Laverania and show that diversity generated by recombination is concentrated in a number of hotspots. An analysis of the subdomain structure indicates that some existing definitions may need to be revised From the analysis of this data, we can now understand the way in which the family has evolved and how the diversity is continuously being generated. Finally, we demonstrate that because the genes are distributed across the genome, sequence sharing between genotypes acts as a useful population genetic marker.

Repeated clinical malaria episodes are associated with modification of the immune system in children

BACKGROUND

There are over 200 million reported cases of malaria each year, and most children living in endemic areas will experience multiple episodes of clinical disease before puberty. We set out to understand how frequent clinical malaria, which elicits a strong inflammatory response, affects the immune system and whether these modifications are observable in the absence of detectable parasitaemia.

METHODS

We used a multi-dimensional approach comprising whole blood transcriptomic, cellular and plasma cytokine analyses on a cohort of children living with endemic malaria, but uninfected at sampling, who had been under active surveillance for malaria for 8 years. Children were categorised into two groups depending on the cumulative number of episodes experienced: high (≥ 8) or low (< 5).

RESULTS

We observe that multiple episodes of malaria are associated with modification of the immune system. Children who had experienced a large number of episodes demonstrated upregulation of interferon-inducible genes, a clear increase in circulating levels of the immunoregulatory cytokine IL-10 and enhanced activation of neutrophils, B cells and CD8⁺ T cells.

CONCLUSION

Transcriptomic analysis together with cytokine and immune cell profiling of peripheral blood can robustly detect immune differences between children with different numbers of prior malaria episodes. Multiple episodes of malaria are associated with modification of the immune system in children. Such immune modifications may have implications for the initiation of subsequent immune responses and the induction of vaccine-mediated protection.

Optimized whole-genome amplification strategy for extremely AT-biased template

Pathogen genome sequencing directly from clinical samples is quickly gaining importance in genetic and medical research studies. However, low DNA yield from blood-borne pathogens is often a limiting factor. The problem worsens in extremely base-biased genomes such as the AT-rich Plasmodium falciparum. We present a strategy for whole-genome amplification (WGA) of low-yield samples from P. falciparum prior to short-read sequencing. We have developed WGA conditions that incorporate tetramethylammonium chloride for improved amplification and coverage of AT-rich regions of the genome. We show that this method reduces amplification bias and chimera formation. Our data show that this method is suitable for as low as 10 pg input DNA, and offers the possibility of sequencing the parasite genome from small blood samples.

A comprehensive evaluation of assembly scaffolding tools

Background

Genome assembly is typically a two-stage process: contig assembly followed by the use of paired sequencing reads to join contigs into scaffolds. Scaffolds are usually the focus of reported assembly statistics; longer scaffolds greatly facilitate the use of genome sequences in downstream analyses, and it is appealing to present larger numbers as metrics of assembly performance. However, scaffolds are highly prone to errors, especially when generated using short reads, which can directly result in inflated assembly statistics.

Results

Here we provide the first independent evaluation of scaffolding tools for second-generation sequencing data. We find large variations in the quality of results depending on the tool and dataset used. Even extremely simple test cases of perfect input, constructed to elucidate the behaviour of each algorithm, produced some surprising results. We further dissect the performance of the scaffolders using real and simulated sequencing data derived from the genomes of Staphylococcus aureus, Rhodobacter sphaeroides, Plasmodium falciparum and Homo sapiens. The results from simulated data are of high quality, with several of the tools producing perfect output. However, at least 10% of joins remains unidentified when using real data.

Conclusions

The scaffolders vary in their usability, speed and number of correct and missed joins made between contigs. Results from real data highlight opportunities for further improvements of the tools. Overall, SGA, SOPRA and SSPACE generally outperform the other tools on our datasets. However, the quality of the results is highly dependent on the read mapper and genome complexity.

REAPR: a universal tool for genome assembly evaluation

Methods to reliably assess the accuracy of genome sequence data are lacking. Currently completeness is only described qualitatively and mis-assemblies are overlooked. Here we present REAPR, a tool that precisely identifies errors in genome assemblies without the need for a reference sequence. We have validated REAPR on complete genomes or de novo assemblies from bacteria, malaria and Caenorhabditis elegans, and demonstrate that 86% and 82% of the human and mouse reference genomes are error-free, respectively. When applied to an ongoing genome project, REAPR provides corrected assembly statistics allowing the quantitative comparison of multiple assemblies. REAPR is available at http://www.sanger.ac.uk/resources/software/reapr/.

A post-assembly genome-improvement toolkit (PAGIT) to obtain annotated genomes from contigs

Genome projects now produce draft assemblies within weeks owing to advanced high-throughput sequencing technologies. For milestone projects such as Escherichia coli or Homo sapiens, teams of scientists were employed to manually curate and finish these genomes to a high standard. Nowadays, this is not feasible for most projects, and the quality of genomes is generally of a much lower standard. This protocol describes software (PAGIT) that is used to improve the quality of draft genomes. It offers flexible functionality to close gaps in scaffolds, correct base errors in the consensus sequence and exploit reference genomes (if available) in order to improve scaffolding and generating annotations. The protocol is most accessible for bacterial and small eukaryotic genomes (up to 300 Mb), such as pathogenic bacteria, malaria and parasitic worms. Applying PAGIT to an E. coli assembly takes ∼24 h: it doubles the average contig size and annotates over 4,300 gene models.

Genome wide adaptations of Plasmodium falciparum in response to lumefantrine selective drug pressure

The combination therapy of the Artemisinin-derivative Artemether (ART) with Lumefantrine (LM) (Coartem®) is an important malaria treatment regimen in many endemic countries. Resistance to Artemisinin has already been reported, and it is feared that LM resistance (LMR) could also evolve quickly. Therefore molecular markers which can be used to track Coartem® efficacy are urgently needed. Often, stable resistance arises from initial, unstable phenotypes that can be identified in vitro. Here we have used the Plasmodium falciparum multidrug resistant reference strain V1S to induce LMR in vitro by culturing the parasite under continuous drug pressure for 16 months. The initial IC₅₀ (inhibitory concentration that kills 50% of the parasite population) was 24 nM. The resulting resistant strain V1S_LM, obtained after culture for an estimated 166 cycles under LM pressure, grew steadily in 378 nM of LM, corresponding to 15 times the IC₅₀ of the parental strain. However, after two weeks of culturing V1S_LM in drug-free medium, the IC₅₀ returned to that of the initial, parental strain V1S. This transient drug tolerance was associated with major changes in gene expression profiles: using the PFSANGER Affymetrix custom array, we identified 184 differentially expressed genes in V1S_LM. Among those are 18 known and putative transporters including the multidrug resistance gene 1 (pfmdr1), the multidrug resistance associated protein and the V-type H+ pumping pyrophosphatase 2 (pfvp2) as well as genes associated with fatty acid metabolism. In addition we detected a clear selective advantage provided by two genomic loci in parasites grown under LM drug pressure, suggesting that all, or some of those genes contribute to development of LM tolerance – they may prove useful as molecular markers to monitor P. falciparum LM susceptibility.

Iterative Correction of Reference Nucleotides (iCORN) using second generation sequencing technology

Motivation: The accuracy of reference genomes is important for downstream analysis but a low error rate requires expensive manual interrogation of the sequence. Here, we describe a novel algorithm (Iterative Correction of Reference Nucleotides) that iteratively aligns deep coverage of short sequencing reads to correct errors in reference genome sequences and evaluate their accuracy.

Results: Using Plasmodium falciparum (81% A + T content) as an extreme example, we show that the algorithm is highly accurate and corrects over 2000 errors in the reference sequence. We give examples of its application to numerous other eukaryotic and prokaryotic genomes and suggest additional applications.

Availability: The software is available at http://icorn.sourceforge.net

Contact:tdo@sanger.ac.uk; cnewbold@hammer.imm.ox.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Functional identification of the Plasmodium centromere and generation of a Plasmodium artificial chromosome

The artificial chromosome represents a useful tool for gene transfer, both as cloning vectors and in chromosome biology research. To generate a Plasmodium artificial chromosome (PAC), we had to first functionally identify and characterize the parasite's centromere. A putative centromere (pbcen5) was cloned from chromosome 5 of the rodent parasite P. berghei based on a Plasmodium gene-synteny map. Plasmids containing pbcen5 were stably maintained in parasites during a blood-stage infection with high segregation efficiency, without drug pressure. pbcen5-containing plasmids were also stably maintained during parasite meiosis and mitosis in the mosquito. A linear PAC (L-PAC) was generated by integrating pbcen5 and telomere into a plasmid. The L-PAC segregated with a high efficiency and was stably maintained throughout the parasite's life cycle, as either one or two copies. These results suggest that L-PAC behaves like a Plasmodium chromosome, which can be exploited as an experimental research tool.

New insights into the blood-stage transcriptome of Plasmodium falciparum using RNA-Seq

Recent advances in high-throughput sequencing present a new opportunity to deeply probe an organism's transcriptome. In this study, we used Illumina-based massively parallel sequencing to gain new insight into the transcriptome (RNA-Seq) of the human malaria parasite, Plasmodium falciparum. Using data collected at seven time points during the intraerythrocytic developmental cycle, we (i) detect novel gene transcripts; (ii) correct hundreds of gene models; (iii) propose alternative splicing events; and (iv) predict 5' and 3' untranslated regions. Approximately 70% of the unique sequencing reads map to previously annotated protein-coding genes. The RNA-Seq results greatly improve existing annotation of the P. falciparum genome with over 10% of gene models modified. Our data confirm 75% of predicted splice sites and identify 202 new splice sites, including 84 previously uncharacterized alternative splicing events. We also discovered 107 novel transcripts and expression of 38 pseudogenes, with many demonstrating differential expression across the developmental time series. Our RNA-Seq results correlate well with DNA microarray analysis performed in parallel on the same samples, and provide improved resolution over the microarray-based method. These data reveal new features of the P. falciparum transcriptional landscape and significantly advance our understanding of the parasite's red blood cell-stage transcriptome.

ABACAS: algorithm-based automatic contiguation of assembled sequences

Summary: Due to the availability of new sequencing technologies, we are now increasingly interested in sequencing closely related strains of existing finished genomes. Recently a number of de novo and mapping-based assemblers have been developed to produce high quality draft genomes from new sequencing technology reads. New tools are necessary to take contigs from a draft assembly through to a fully contiguated genome sequence. ABACAS is intended as a tool to rapidly contiguate (align, order, orientate), visualize and design primers to close gaps on shotgun assembled contigs based on a reference sequence. The input to ABACAS is a set of contigs which will be aligned to the reference genome, ordered and orientated, visualized in the ACT comparative browser, and optimal primer sequences are automatically generated.

Availability and Implementation: ABACAS is implemented in Perl and is freely available for download from http://abacas.sourceforge.net

Contact:sa4@sanger.ac.uk

Genome-wide discovery and verification of novel structured RNAs in Plasmodium falciparum

We undertook a genome-wide search for novel noncoding RNAs (ncRNA) in the malaria parasite Plasmodium falciparum. We used the RNAz program to predict structures in the noncoding regions of the P. falciparum 3D7 genome that were conserved with at least one of seven other Plasmodium spp. genome sequences. By using Northern blot analysis for 76 high-scoring predictions and microarray analysis for the majority of candidates, we have verified the expression of 33 novel ncRNA transcripts including four members of a ncRNA family in the asexual blood stage. These transcripts represent novel structured ncRNAs in P. falciparum and are not represented in any RNA databases. We provide supporting evidence for purifying selection acting on the experimentally verified ncRNAs by comparing the nucleotide substitutions in the predicted ncRNA candidate structures in P. falciparum with the closely related chimp malaria parasite P. reichenowi. The high confirmation rate within a single parasite life cycle stage suggests that many more of the predictions may be expressed in other stages of the organism's life cycle.

Plasmodium falciparum var gene expression is developmentally controlled at the level of RNA polymerase II‐mediated transcription initiation

The Plasmodium falciparum var gene family codes for a major virulence factor in this most lethal of human malaria parasites. A single var protein variant type is expressed on each infected red blood cell, with antigenic variation allowing progeny parasites to escape host immune detection. The control of mutually exclusive var gene expression in the parasite relies on in situ epigenetic changes. Whether control of expression occurs at transcription initiation or post transcription, however, remains to be established. Recent evidence supports existence of a unique var transcription site at the nuclear periphery containing the dominantly expressed var gene, although silent var genes can colocalize to the same region. We demonstrate here that exclusive var gene expression is controlled at the level of transcription initiation during ring stages and that var genes are transcribed by RNA polymerase II. This represents another example where P. falciparum differs from the paradigm for antigenic variation, Trypanosoma brucei.

Genome variation and evolution of the malaria parasite Plasmodium falciparum

Infections with the malaria parasite Plasmodium falciparum result in more than 1 million deaths each year worldwide. Deciphering the evolutionary history and genetic variation of P. falciparum is critical for understanding the evolution of drug resistance, identifying potential vaccine candidates and appreciating the effect of parasite variation on prevalence and severity of malaria in humans. Most studies of natural variation in P. falciparum have been either in depth over small genomic regions (up to the size of a small chromosome) or genome wide but only at low resolution. In an effort to complement these studies with genome-wide data, we undertook shotgun sequencing of a Ghanaian clinical isolate (with fivefold coverage), the IT laboratory isolate (with onefold coverage) and the chimpanzee parasite P. reichenowi (with twofold coverage). We compared these sequences with the fully sequenced P. falciparum 3D7 isolate genome. We describe the most salient features of P. falciparum polymorphism and adaptive evolution with relation to gene function, transcript and protein expression and cellular localization. This analysis uncovers the primary evolutionary changes that have occurred since the P. falciparum-P. reichenowi speciation and changes that are occurring within P. falciparum.

Microarray-based comparative genomic analyses of the human malaria parasite Plasmodium falciparum using Affymetrix arrays

Microarray-based comparative genomic hybridization (CGH) provides a powerful tool for whole genome analyses and the rapid detection of genomic variation that underlies virulence and disease. In the field of Plasmodium research, many of the parasite genomes that one might wish to study in a high throughput manner are not laboratory clones, but clinical isolates. One of the key limitations to the use of clinical samples in CGH, however, is the miniscule amounts of genomic DNA available. Here we describe the successful application of multiple displacement amplification (MDA), a non-PCR-based amplification method that exhibits clear advantages over all other currently available methods. Using MDA, CGH was performed on a panel of NF54 and IT/FCR3 clones, identifying previously published deletions on chromosomes 2 and 9 as well as polymorphism in genes associated with disease pathology.

Transient cross-reactive immune responses can orchestrate antigenic variation in malaria

The malaria parasite Plasmodium falciparum has evolved to prolong its duration of infection by antigenic variation of a major immune target on the surface of the infected red blood cell. This immune evasion strategy depends on the sequential, rather than simultaneous, appearance of immunologically distinct variants. Although the molecular mechanisms by which a single organism switches between variants are known in part, it remains unclear how an entire population of parasites within the host can synchronize expression to avoid rapidly exhausting the variant repertoire. Here we show that short-lived, partially cross-reactive immune responses to parasite-infected erythrocyte surface antigens can produce a cascade of sequentially dominant antigenic variants, each of which is the most immunologically distinct from its preceding types. This model reconciles several previously unexplained and apparently conflicting epidemiological observations by demonstrating that individuals with stronger cross-reactive immune responses can, paradoxically, be more likely to sustain chronic infections. Antigenic variation has always been seen as an adaptation of the parasite to evade host defence: we show that the coordination necessary for the success of this strategy might be provided by the host.

Transcription of a subtelomerically located var gene variant in Plasmodium falciparum appears to require the truncation of an adjacent var gene

The Plasmodium falciparum R29 clone preferentially transcribes the R29var gene variant on rosette selection, unlike other isogenic clones from the same parasite lineage. Characterisation of the R29var gene locus revealed that this gene lies internal to, and is in a tail-to-tail orientation with, a second var gene variant (A4var) at one end of chromosome 13. In the R29 clone, a spontaneous deletion event between these two var variants deletes all of the A4var gene and the subtelomeric repetitive sequence arrays. We have previously shown that a simple disruption of the A4var gene is not sufficient to preferentially activate the R29var gene in rosette-selected parasites. We therefore hypothesised that the truncation of the chromosome end may be a key factor in predisposing the R29var variant to transcription under rosette selection conditions. Here, we have generated a panel of isogenic parasite clones with both intact and truncated A4var-R29var loci, and show that R29var transcription is only detected in rosette-selected clones with a truncated locus. Furthermore, we present provisional data describing the relative frequency with which this spontaneous deletion event occurs. These data have implications in our understanding of how spontaneous deletion events within subtelomeric var loci may affect transcription of these var gene variants.

Sequence of Plasmodium falciparum chromosomes 1, 3-9 and 13

Since the sequencing of the first two chromosomes of the malaria parasite, Plasmodium falciparum, there has been a concerted effort to sequence and assemble the entire genome of this organism. Here we report the sequence of chromosomes 1, 3-9 and 13 of P. falciparum clone 3D7--these chromosomes account for approximately 55% of the total genome. We describe the methods used to map, sequence and annotate these chromosomes. By comparing our assemblies with the optical map, we indicate the completeness of the resulting sequence. During annotation, we assign Gene Ontology terms to the predicted gene products, and observe clustering of some malaria-specific terms to specific chromosomes. We identify a highly conserved sequence element found in the intergenic region of internal var genes that is not associated with their telomeric counterparts.

Genome sequence of the human malaria parasite Plasmodium falciparum

The parasite Plasmodium falciparum is responsible for hundreds of millions of cases of malaria, and kills more than one million African children annually. Here we report an analysis of the genome sequence of P. falciparum clone 3D7. The 23-megabase nuclear genome consists of 14 chromosomes, encodes about 5,300 genes, and is the most (A + T)-rich genome sequenced to date. Genes involved in antigenic variation are concentrated in the subtelomeric regions of the chromosomes. Compared to the genomes of free-living eukaryotic microbes, the genome of this intracellular parasite encodes fewer enzymes and transporters, but a large proportion of genes are devoted to immune evasion and host-parasite interactions. Many nuclear-encoded proteins are targeted to the apicoplast, an organelle involved in fatty-acid and isoprenoid metabolism. The genome sequence provides the foundation for future studies of this organism, and is being exploited in the search for new drugs and vaccines to fight malaria.

Stage-specific promoter activity from stably maintained episomes in Plasmodium falciparum

Genomic DNA is organised at its simplest level within phased arrays of nucleosomes, a structure key to the correct transcriptional regulation of the encoded genes. Here we studied chromatin formation on DNA transfected into Plasmodium falciparum either as an episomal plasmid or following integration by homologous recombination. We show that stably maintained and replicated plasmid assembles phased arrays of nucleosomes and that a reporter gene is transcribed in an appropriate temporal manner. These data provide a key observation for the future investigation of promoter structure and function with transfected DNA in Plasmodium spp.

Cellular responses to Plasmodium falciparum erythrocyte membrane protein-1: use of relatively conserved synthetic peptide pools to determine CD4 T cell responses in malaria-exposed individuals in Benin, West Africa

BACKGROUND

Plasmodium falciparum erythrocyte membrane protein-1, a variant antigen of the malaria parasite, is potentially a target for the immune response. It would be important to determine whether there are CD4 T cells that recognise conserved regions. However, within the relatively conserved region, there is variation. It is not possible to test T cell responses from small field samples with all possible peptides.

METHODS

We have aligned sequences that are relatively conserved between several PfEMP1 molecules, and chosen a representative sequence similar to most of the PfEMP1 variants. Using these peptides as pools representing CIDRalpha, CIDRbeta and DBLbeta-delta domains, DBLalpha domain, and EXON 2 domain of PfEMP1, we measured the CD4 T cell responses of malaria-exposed donors from Benin, West Africa by a FACS based assay.

RESULTS

All the three peptide pools elicited a CD4 T cell response in a proportion of malaria-exposed and non-exposed donors. CD4 T cell proliferation occurs at a relatively higher magnitude to peptide pools from the DBLalpha and EXON 2 in the malaria-exposed donors living in Benin than in the UK malaria-unexposed donors.

CONCLUSIONS

These findings suggest that an immunological recall response to conserved peptides of a variant antigen can be measured. Further testing of individual peptides in a positive pool will allow us to determine those conserved sequences recognised by many individuals. These types of assays may provide information on conserved peptides of PfEMP1 which could be useful for stimulating T cells to provide help to P. falciparum specific B cells.

CD4 T cell responses to a variant antigen of the malaria parasite Plasmodium falciparum, erythrocyte membrane protein-1, in individuals living in malaria-endemic areas

Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP-1) is a variant antigen on the surface of malaria-infected red blood cells. Antibody responses to PfEMP-1 correlate with immunity, and, therefore, PfEMP-1 may be a good candidate for a malaria vaccine. However, the specificity of CD4 T cells required for a protective variant-specific antibody response is not known. We have measured the CD4 T cell response to 3 different regions that are relatively homologous among different PfEMP-1 variants. The response to the cysteine-rich interdomain region was unusual in that the majority of donors, whether malaria exposed or not, had positive CD4 T cell, interleukin-10, and interferon-gamma responses. The CD4 T cell response to the exon 2 and duffy binding-like domain proteins was significantly greater in malaria-exposed donors than in unexposed Europeans, which suggests that these regions contain peptides recognized by T cells, which thus may be useful as components of a vaccine.

Antigenic variation at the infected red cell surface in malaria

Many pathogens that either rely on an insect vector to complete their life cycle (e.g., Trypanosoma spp. and Borrelia spp.) or exist in a unique ecological niche where transmission from host to host is sporadic (e.g., Neisseria spp.) have evolved strategies to maintain infection of their mammalian hosts for long periods of time in order to ensure their survival. Because they have to survive in the face of a fully functional immune system, a common feature of many of these organisms is their development of sophisticated strategies for immune evasion. For the above organisms and for malaria parasites of the genus Plasmodium, a common theme is the ability to undergo clonal antigenic variation. In all cases, surface molecules that are important targets of the humoral immune response are encoded in the genome as multicopy, nonallelic gene families. Antigenic variation is accomplished by the successive expression of members of these gene families that show little or no immunological cross-reactivity. In the case of malaria parasites, however, some of the molecules that undergo antigenic variation are also major virulence factors, adding an additional level of complication to the host-parasite interaction. In this review, we cover the history of antigenic variation in malaria and then summarize the more recent data with particular emphasis on Plasmodium falciparum, the etiological agent of the most severe form of human malaria.

A functional analysis of a natural variant of intercellular adhesion molecule-1 (ICAM-1Kilifi)

Intercellular adhesion molecule-1 (ICAM-1) is involved in a range of interactions both within the host and between the host and a number of pathogens. Recently we described a mutation within the coding region of the first N-terminal immunoglobulin-like domain of ICAM-1, present at high frequency within African populations, which increased the risk of cerebral malaria. To understand the mechanism by which such a polymorphism might be maintained despite counter-selection by malaria, we have carried out functional assays using both forms of ICAM-1 as soluble Fc chimeric fusion proteins. ICAM-1Kilifi has reduced avidity for LFA-1 compared with ICAM-1ref and binding to soluble fibrinogen was completely abolished with the Kilifi variant. In Plasmodium falciparum adhesion assays, ITO4-A4u binding to ICAM-1Kilifi was reduced compared with binding to the reference form. These results allow for the possibility of balanced selection between the reference and Kilifi forms of ICAM-1 through modulation of inflammatory responses and indicate the existence of differences within ICAM-1-binding P. falciparum isolates which may be relevant to pathogenesis.

The complete nucleotide sequence of chromosome 3 of Plasmodium falciparum

Analysis of Plasmodium falciparum chromosome 3, and comparison with chromosome 2, highlights novel features of chromosome organization and gene structure. The sub-telomeric regions of chromosome 3 show a conserved order of features, including repetitive DNA sequences, members of multigene families involved in pathogenesis and antigenic variation, a number of conserved pseudogenes, and several genes of unknown function. A putative centromere has been identified that has a core region of about 2 kilobases with an extremely high (adenine + thymidine) composition and arrays of tandem repeats. We have predicted 215 protein-coding genes and two transfer RNA genes in the 1,060,106-base-pair chromosome sequence. The predicted protein-coding genes can be divided into three main classes: 52.6% are not spliced, 45.1% have a large exon with short additional 5' or 3' exons, and 2.3% have a multiple exon structure more typical of higher eukaryotes.

Cytoadherence, pathogenesis and the infected red cell surface in Plasmodium falciparum

The particular virulence of Plasmodium falciparum compared with the other malaria species which naturally infect humans is thought to be due to the way in which the parasite modifies the surface of the infected red cell. Approximately 16 hours into the asexual cycle, parasite encoded proteins appear on the red cell surface which mediate adherence to a variety of host tissues. Binding of infected red cells to vascular endothelium, a process which occurs in all infections, is thought to be an important factor in the pathogenesis of severe disease where concentration of organisms in particular organs such as the brain occurs. Binding to uninfected red cells to form erythrocyte rosettes, a property of some isolates, is linked to disease severity. Here we summarise the data on the molecular basis of these interactions on both the host and parasite surfaces and review the evidence for the involvement of particular receptors in specific disease syndromes. Finally we discuss the relevance of these data to the development of new treatments for malaria.

Immunity to non-cerebral severe malaria is acquired after one or two infections

In areas of stable transmission, clinical immunity to mild malaria is acquired slowly, so it is not usually effective until early adolescence. Life-threatening disease is, however, restricted to a much younger age group, indicating that resistance to the severe clinical consequences of infection is acquired more quickly. Understanding how rapidly immunity develops to severe malaria is essential, as severe malaria should be the primary target of intervention strategies, and predicting the result of interventions that reduce host exposure will require consideration of these dynamics. Severe disease in childhood is less frequent in areas where transmission is the greatest. One explanation for this is that infants experience increased exposure to infection while they are protected from disease, possibly by maternal antibody. They therefore emerge from this period of clinical protection with considerably more immunity than those who experience lower transmission intensities. Here we use this data, assuming a period of clinical protection, to estimate the number of prior infections needed to reduce the risk of severe disease to negligible levels. Contrary to expectations, one or two successful infective bites seem to be all that is necessary across a broad range of transmission intensities.

Acquired immunity and postnatal clinical protection in childhood cerebral malaria

By analysing data on the age distribution of cerebral malaria among sites of different transmission intensities, we conclude that the most plausible explanation for the epidemiological patterns seen is that (i) cerebral malaria is caused by a distinct set of Plasmodium falciparum antigenic types; (ii) these antigenic types or 'CM strains' are very common and induce strong strain-specific immunity; and (iii) the postnatal period of protection against cerebral malaria is much longer than the period of protection against other forms of severe disease. The alternative hypothesis that cerebral malaria may be caused by any 'strain' of P. falciparum is compatible with the data only if a single exposure is sufficient to protect against further episodes. This is not consistent with observations on the history of exposure of patients with cerebral malaria. Finally, it is clear that although the delayed peak in incidence of cerebral malaria (with age) can be generated by assuming that subsequent exposures carry a higher risk of disease, such an explanation is not compatible with the observation that severe disease rates are low among infants and young children in areas of high transmissibility.

Receptor-specific adhesion and clinical disease in Plasmodium falciparum

One important factor in the virulence of infections with Plasmodium falciparum is the adherence of infected erythrocytes to small vessel endothelium. In infections that lead to serious, life-threatening disease accumulation of large numbers of infected cells in particular organs is thought to lead to organ dysfunction or failure. This is of particular relevance when the affected organ is the brain, leading to the development of cerebral malaria. Many different endothelial receptors for infected red blood cells have been identified. Some receptors such as CD36 and thrombospondin are used by all parasite isolates, whereas others such as intercellular adhesion molecule-1 (ICAM-1) or vascular cell adhesion molecule (VCAM) are used by a subset of field and laboratory isolates. While it has been speculated that the ability to bind or affinity of binding to a particular endothelial receptor may be related to the pattern of disease, only studies with limited numbers of patients have been carried out to date and these have been in general inconclusive. Here we have taken parasite isolates from 150 patients with defined clinical syndromes as well as isolates from 50 healthy but parasitized community controls and quantitatively assessed their binding to purified endothelial receptors in vitro. Our results show that disregarding the level of adhesion, all parasites bind to CD36, most bind to ICAM-1, few bind to VCAM, and almost none bind to E-selectin. In assessing the degree of binding we show that 1) binding to all receptors was reduced in parasites taken from severely anemic patients; 2) binding to CD36 is identical in parasites from cerebral malaria patients and community controls but slightly elevated in parasites from nonsevere cases; and 3) binding to ICAM-1 is highest in cerebral malaria patients. Because rosette formation by uninfected cells has also been a phenotype associated with disease severity and one that may interfere in vitro with receptor binding, we also assessed rosette formation in all isolates. In this study the highest level of rosette-forming parasites was found in the anemic group and not the cerebral malaria group. Stratifying the data for the frequency of rosette formation showed that the above results were not significantly altered by this phenomenon. Our data are not consistent with a role for binding to CD36 in the development of severe disease but show an association between the degree of binding to ICAM-1 and clinical illness in nonanemic patients.

Limited spatial clustering of individual Plasmodium falciparum alleles in field isolates from coastal Kenya

We describe Plasmodium falciparum genetic diversity in coastal Kenya, typing S-antigen and the merozoite surface proteins 1 and 2 (MSP-1 and MSP-2) in field isolates by the polymerase chain reaction (PCR). Malaria in coastal Kenya is characterized by low seasonal transmission, and a relatively high incidence of severe disease, which tends to occur in time-space clusters. We chose the highly polymorphic S-antigen as a marker for localized parasite diversity because it has been shown to vary in serotype prevalence in time and space. A total of 261 children (up to nine years of age) in two neighboring locations with different transmission rates were sampled for blood-stage parasites in cross-sectional surveys before and after the main transmission period in 1991, and also in a concomitant one-year longitudinal survey tracing clinical infections. Six major sequence types of S-antigen were identified, which were subdivided into 70 alleles; however, only 50% of isolates were typeable. The S-antigen sequence types varied qualitatively between locations, over time, and between asymptomatic and clinical disease infections, but not between different age groups. The MSP-1 and MSP-2 sequence type prevalences, in contrast, did not differ in any of these comparisons. We describe the use of the Mantel test for assessing clustering of individual parasite alleles at the household level, and demonstrate low-level clustering of MSP-1 and MSP-2 alleles and S-antigen sequence types, at the end of a long period of low transmission.

Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa

BACKGROUND

Malaria remains a major cause of mortality and morbidity in Africa. Many approaches to malaria control involve reducing the chances of infection but little is known of the relations between parasite exposure and the development of effective clinical immunity so the long-term effect of such approaches to control on the pattern and frequency of malaria cannot be predicted.

METHODS

We have prospectively recorded paediatric admissions with severe malaria over three to five years from five discrete communities in The Gambia and Kenya. Demographic analysis of the communities exposed to disease risk allowed the estimation of age-specific rates for severe malaria. Within each community the exposure to Plasmodium falciparum infection was determined through repeated parasitological and serological surveys among children and infants. We used acute respiratory-tract infections (ARI) as a comparison.

FINDINGS

3556 malaria admissions were recorded for the five sites. Marked differences were observed in age, clinical spectrum and rates of severe malaria between the five sites. Paradoxically, the risks of severe disease in childhood were lowest among populations with the highest transmission intensities, and the highest disease risks were observed among populations exposed to low-to-moderate intensities of transmission. For severe malaria, for example, admission rates (per 1000 per year) for children up to their 10th birthday were estimated as 3.9, 25.8, 25.9, 16.7, and 18.0 in the five communities; the forces of infection estimated for those communities (new infections per infant per month) were 0.001, 0.034, 0.050, 0.093, and 0.176, respectively. Similar trends were noted for cerebral malaria and for severe malaria anaemia but not for ARI. Mean age of disease decreased with increasing transmission intensity.

INTERPRETATION

We propose that a critical determinant of life-time disease risk is the ability to develop clinical immunity early in life during a period when other protective mechanisms may operate. In highly endemic areas measures which reduce parasite transmission, and thus immunity, may lead to a change in both the clinical spectrum of severe disease and the overall burden of severe malaria morbidity.

Periodicity and space-time clustering of severe childhood malaria on the coast of Kenya

Traditionally malaria epidemiology has focused on factors such as parasite rates and vector dynamics without specific reference to disease. There are limited comprehensive data on malaria as a life-threatening event in African children. We have identified, through hospital surveillance, 581 episodes of severe malaria in residents of a defined area on the Kenya coast over a period of 3 years. This represents an absolute minimum risk of developing severe malaria by the fifth birthday of 1 in 15. The presentation of severe malaria showed marked seasonality, but the timing and magnitude of these fluctuations varied considerably between years. A satellite navigational system was used to define the exact location of the home of each severe malaria case. Space-time clustering of severe malaria was evident in this community. Seasonal peaks in incidence of severe malaria may comprise discrete mini-epidemics. In contrast, parasite rates in the community varied little during the course of the surveillance. The monitoring of disease, as opposed to parasitization, in children may result in more effective targeting of intervention resources.

Rheological analysis of the adhesive interactions of red blood cells parasitized by Plasmodium falciparum

Adhesion of parasitized red blood cells (RBCs) to vascular endothelium is thought to be a key factor in the pathology of falciparum malaria. However, quantitative analyses of the intercellular forces and of the effects of flow on adhesion have been lacking. We have characterized cytoadhesion of RBCs parasitized by the strains ITO4 (which can bind to receptors ICAM-1 or CD36) and FCR3A2 (which can bind to CD36 only) using micropipette manipulation and flow chamber techniques. Target cells were unfixed or glutaraldehyde-fixed human umbilical vein endothelial cells (HUVEC, bearing ICAM-1 only) or human amelanotic melanoma cells (C32, bearing CD36 and ICAM-1). In the static, micropipette assay, 60% to 70% of parasitized cells would adhere when tested at up to three successive sites. The percentage of cells adhering and the force required for their detachment (approximately 10(-10) N) were similar for each combination of parasite strain and adhesion target (ITO4/HUVEC, ITO4/C32, FCR3A2/C32). In the flow chamber, efficiency of initial adhesion of parasitized cells was essentially constant (at about 1%) up to a stress of 0.1 Pa, and then decreased rapidly with increasing stress. Either receptor (ICAM-1 or CD36) could immobilize flowing cells at a physiologic flow stress (0.1 Pa), but the numbers of cells adhering varied for the different combinations (ITO4/C32 greater than ITO4/HUVEC greater than FCR3A2/C32). When flow was increased in steps, adhered cells were gradually washed off but many could withstand stresses at which they would not initially adhere. The force for detachment estimated in this way was similar to the pipette value, and again, was similar for the different combinations of strains and targets. Adhesion from flow depends on the affinity between surfaces being above a critical level, and once adhesion is established, the fracture energy determines resistance to disruption of adhesion. The results show that the fracture energy is greater than the affinity (ie, that adhesion becomes stabilized after it is initially established) and that the ratio of affinity to fracture energy is different for different receptor/ligand pairs, with ICAM-1 appearing to be the more efficient immobilizing receptor. Also, static and flow-based assays of adhesion clearly differ; the affinity is less critical in the static situation, so that most parasitized cells were capable of adhering in a static assay, but fewer did so under flow. Adhesiveness varied markedly from cell to cell, both for targets and parasitized cells.(ABSTRACT TRUNCATED AT 400 WORDS)