- Open Access
Differences in PfEMP1s recognized by antibodies from patients with uncomplicated or severe malaria
Malaria Journalvolume 15, Article number: 258 (2016)
Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) variants are encoded by var genes and mediate pathogenic cytoadhesion and antigenic variation in malaria. PfEMP1s can be broadly divided into three principal groups (A, B and C) and they contain conserved arrangements of functional domains called domain cassettes. Despite their tremendous diversity there is compelling evidence that a restricted subset of PfEMP1s is expressed in severe disease. In this study antibodies from patients with severe and uncomplicated malaria were compared for differences in reactivity with a range of PfEMP1s to determine whether antibodies to particular PfEMP1 domains were associated with severe or uncomplicated malaria.
Parts of expressed var genes in a severe malaria patient were identified by RNAseq and several of these partial PfEMP1 domains were expressed together with others from laboratory isolates. Antibodies from Papuan patients to these parts of multiple PfEMP1 proteins were measured.
Patients with uncomplicated malaria were more likely to have antibodies that recognized PfEMP1 of Group C type and recognized a broader repertoire of group A and B PfEMP1s than patients with severe malaria.
These data suggest that exposure to a broad range of group A and B PfEMP1s is associated with protection from severe disease in Papua, Indonesia.
PfEMP1 is the immunodominant antigen of the malaria parasite Plasmodium falciparum expressed on the surface of the infected erythrocyte (IE). Adherence of this molecule to host receptors expressed on endothelial cells, uninfected erythrocytes and placental syncytiotrophoblasts facilitates sequestration of IE in vascular tissues, avoiding destruction in the spleen [1–3]. PfEMP1 molecules are encoded by the var multigene family [1–3]. Individual parasites have approximately 60 var gene variants and switching between single, transcribed var genes leads to changes in cytoadhesive phenotype as well as clonal antigenic variation and immune escape. var gene repertoires differ among isolates  and immunity to malaria is dependent on acquisition of antibodies to a range of PfEMP1 variants [5–8]. Immunity to both cerebral malaria  and non-cerebral, severe malaria  is acquired much more rapidly than immunity to uncomplicated malaria. Parasites that cause severe disease appear to express a conserved subset of variant antigens that are encountered earlier in life and that are thus more widely recognized by sera from semi-immune children than parasites causing uncomplicated disease [11, 12].
PfEMP1s contain combinations of Duffy binding-like domains (DBLα, β, γ, δ, ε, ζ and x) and cysteine rich inter-domain regions (CIDRα, β, γ and δ) . Some DBL and CIDR domain subtypes mediate adhesion to different host receptors (reviewed in [14, 15]), and some are organized in semi-conserved domain cassettes (DC) that are present in most parasites . var genes are also classified using their upstream sequence into groups A, B, C [16, 17] which comprise 20, 60 and 20 % respectively of the var gene repertoire ; the unique var gene called var2csa has a different upstream sequence (ups E) and is only involved in malaria during pregnancy .
The expression of particular subtypes of DBLα domains in severe malaria suggests severe disease may be preferentially caused by a restricted subset of var genes [19, 20]. Increased expression of group A and B var genes has been associated with clinical, but not specifically severe malaria in Papua New Guinea (PNG) [21, 22] and with severe malaria in Africa . Cerebral malaria in Africa was associated with increased expression of group A [20, 24, 25] or group B  var genes.
Consistent with its having a role in severe malaria, PfEMP1s encoded by group A and B var genes appear to be widely expressed by parasites that infect non- or semi-immune individuals. Antibodies from older children preferentially recognized PfEMP1s encoded by Group A var genes, indicating previous exposure . Group A and B var genes dominated infection of a naive individual , and more individuals develop antibodies to group A PfEMP1s than group B or C, and do so at a younger age .
Group A and B var genes also encode adhesion phenotypes associated with severe disease. In Africa the adhesion phenotype of rosetting is associated with severe malaria  and increased expression of group A var genes [19, 21, 25]. Some group A and B PfEMP1s can bind to intercellular adhesion molecule 1 (ICAM-1) [30, 31], and ICAM-1 expression was up-regulated in brain endothelium and co-localized with sequestered IEs in cerebral malaria patients . IE adhesion to ICAM-1 has variously been associated with cerebral malaria , clinical but not severe malaria  or inversely correlated with severe disease .
Another phenotype associated with severe disease is adhesion to endothelial protein C receptor, EPCR [36, 37]. Parasite isolates from African children with severe malaria bound EPCR and expressed DC8 or DC13 var genes [36, 38]. DC8 and DC13 PfEMP1s are primarily group B and A, respectively , and contain members of the subset of CIDRα1 domain types, which bind EPCR [36, 37]. Sera from African children with uncomplicated malaria recognize PfEMP1s containing DC8 and DC13 at higher levels than PfEMP1s without DC8 or 13, but it is unclear whether severe malaria specifically induces antibodies to DC8 and DC13 [39, 40].
DC5 PfEMP1s, which are nearly all group A, are recognized by sera from semi-immune children in a similar manner to other severe malaria associated isolates , and DC5 expression increased markedly during infection of a naive volunteer . African children acquired antibodies to DC5 more rapidly than to other PfEMP1 domains which is consistent with widespread expression of DC5 in non-immune individuals . Antibodies reactive with DC5 also correlated with protection from malaria episodes ; however, the evidence directly linking DC5 PfEMP1 expression and adhesion phenotype to severe malaria is less clear. High levels of DC5 sequence expression have been detected in severe malaria, but only together with expression of either DC8 or 13 sequences  and it is conceivable that these DC5 and DC8 or DC13 sequences were present on the same var genes. DC5 PfEMP1s bind platelet-endothelial cell adhesion molecule 1 (PECAM1)  and IE adhesion to PECAM1 has been implicated in cerebral malaria [44, 45]. However, IE adhesion to PECAM1 is also commonly found in samples from patients with uncomplicated malaria  and DC5 PfEMP1s were not expressed by parasites selected for adhesion to brain endothelium [39, 40]. A minority of var genes containing DC5 do encode EPCR-binding CIDRα1.5 domains although these domains are not part of the DC5 cassette .
Thus several promising candidates have emerged as members of the restricted population of PfEMP1s responsible for severe malaria, but the relative contributions of group A and B, and of DCs 8, 13 and 5 remain unclear. In particular very little is known about PfEMP1s in severe disease in the Asia Pacific region. Determining whether conserved PfEMP1sequences elicit protection from severe malaria disease globally is a priority for vaccine research. In this study plasmas from Papuan patients with severe or uncomplicated malaria were analysed for their reactivity with PfEMP1 polypeptides representative of the different groups and DCs including several expressed by parasites causing cerebral malaria.
Two ml of venous blood was collected from a 5-year old female patient with cerebral malaria in PNG. The IEs were separated from white blood cells by Plasmodipur column filtration (Europroxima) as per the manufacturer’s instructions. The IEs were lysed in TRIzol® (Life Technologies), incubated at 37 °C for 5 min and then stored at −80 °C. Venous samples were collected from patients with severe (n = 28) and uncomplicated (n = 35) malaria attending a healthcare facility in Timika, Papua Province, Indonesia. Plasma was separated from blood by centrifugation and stored at −20 °C. This area has unstable malaria transmission with estimated annual parasite incidence of 450 per 1000 population and symptomatic illness in all ages . Severe malaria was defined as peripheral parasitaemia with at least one modified World Health Organization (WHO) criterion of severity . Twenty-six of the 28 patients with severe malaria had parasitemias greater than 1000/μL, which is a previously-derived threshold that predicts clinical disease in northern Papua . Thus incidental parasitemia is unlikely in these 26 severe malaria patients but cannot be excluded in the two severe malaria patients with parasitemias less than 1000/μl.
Informed consent was provided by all participants. The study was approved in Indonesia by the Eijkman Institute Research Ethics Commission (project number 46), in Australia by the Melbourne Health Human Research Ethics Committee (project number 2010.284) and Human Research Ethics Committee of the NT Department of Health & Families and Menzies School of Health Research, Darwin, Australia (HREC 2010–1396), and in PNG by the Government of PNG Medical Research Advisory Committee (MRAC no. 11.12).
RNA extraction and RNAseq
Erythrocytes in TRIzol® were thawed at 37 °C, chloroform (1/5th of the TRIzol® volume) was added and vortexed 15 s, the solution was then subjected to centrifugation at 12,000×g for 30 min at 4 °C and the aqueous supernatant was aspirated and mixed with an equal volume of 70 % ethanol in RNase free water. This solution was then directly applied to RNeasy mini columns (QIAGEN), and RNA purification was performed with on-column DNasing as per the manufacturer’s instructions. Eluted RNA was oligo dT purified and used to generate a 65 bp paired end RNAseq library using the standard Illumina protocol. The library was sequenced on an Illumina GAII at Geneworks (Adelaide).
Alignments were performed using Subreads . Reads were first aligned to the 399 full length var genes that were previously described in seven P. falciparum isolates  and mapped reads were extracted. Reads were also aligned to the 3D7 P. falciparum genome and unmapped reads were extracted to include any var that may have been missed. The reads that mapped to var genes and the reads that did not map to the 3D7 genome were then merged, digitally normalized using khmer , and assembled using Oases . Assembled transcripts were subjected to a second round of de novo assembly using Cap3  to assemble contigs and extract consensus sequences from the contigs. Contigs were aligned to the domains from the 399 var genes with BLAST (version 2.2.25). The transcript abundance was determined by aligning the reads again to the assembled transcripts.
Proteins for Luminex assays were expressed in baculovirus-transfected insect cells as previously described . Proteins analysed by Luminex were IT4var02 protein 1 DBLγ12DBLδ5CIDRβ3DBLβ9, PF11_0008 DBLδ5CIDRβ4, IT4var02 protein 2 DBLδ5CIDRβ3, HB3var05 DC16 DBLα1.6CIDRδ1, PFI1820w var3 DBLα1.3 DBLε8, PFD0020c DBLβ12 DBLγ6, IT4var20 DBLα2CIDRα1.1DBLβ12DBLγ6DBLδ1CIDRβ1, PF11_0007 DBLα0.15CIDRα3.2DBLδ1CIDRβ1, PFL0020w DBLγ14DBLζ5DBLε4, IT4var13 DC 9 DBLα0.3CIDRα5DBLβ5DBLδ9CIDRγ9DBLγ11DBLζ4.
PfEMP1 domains identified by RNAseq as expressed by the parasites infecting the PNG CM patient were amplified from gDNA purified from the patient’s blood for use as ELISA antigens. The sequences were then cloned into the plasmid pEU-E01H-N1 (encoding an N terminal His-tagged protein; CellFree Sciences Matsuyama, Japan)  and expressed in the wheat germ cell-free expression system (CellFree Sciences) and purified on Ni-nitrilotriacetic acid agarose columns (Qiagen, Valencia, CA) as previously described [55, 56].
Antigens used for the ELISA included proteins one to seven that were encoded by var gene sequences transcribed by parasites infecting the cerebral malaria patient from PNG. Proteins 1 and 3 were encoded by 352 bp and 123 bp fragments of a contig with homology to the DC5 DBLβ7_D8 of Dd2var4 (ranked 33rd by transcript abundance) (Table 1). Protein 2 was encoded by a 861 bp orthologue of DC5 DBLδ5_D6_Dd2var4 that incorporated two contigs that were ranked 4th and 9th by abundance (Table 1). Protein 4 was encoded by a 697 bp fragment of a contig with homology to a DC8 DBLγ6_D3_Dd2var47 (ranked 5th by abundance). Protein 5 was encoded by a 186 bp fragment with homology to DBLγ10 of raj116 and was ranked 34th by abundance. Protein 6 was encoded by a 425 bp contig that incorporated two non-overlapping contigs with homology to NTSB3_DBLα0.11 of igh_var31 and DBLα0.15 of raj116_var34 (ranked respectively14th and 21st by abundance). Protein 7 was a 142 bp cloned orthologue of DBLδ1_D4_PFCLIN_var24 (ranked 96th by abundance).
Proteins 1, 2 and 3 were most closely related to the DC5 domain var gene Dd2var4. This gene has group A-like coding features, e.g. head structure and ATS, and is the only DC5 with a group B-like upstream sequence, the other 11 all being group A . Therefore, proteins 1, 2 and 3 were classified as group A for all subsequent analyses. These PNG derived proteins were supplemented by domains from HB3, 3D7 and ItG parasites that included two DBLδ domains from group C PfEMP1s, a DC8 CIDRα1.1, a DC13 CIDRα1.4 and a CIDRα3.1 (Fig. 2b).
The Luminex assay was performed as previously described . For each protein an eleven point standard curve was made using two-fold dilutions of pooled positive plasma starting with 1/40, which was assigned an arbitrary value of 1000 relative units (RU). Plasmas were diluted in 0.02 % Tween-20, 0.1 % BSA in PBS pH 7.4. Fluorescent intensities of patients’ plasmas were used to interpolate antibody concentrations in RU from the standard curves.
ELISA was performed as previously described . ELISA plates were blocked with 3 % (w/v) skim milk in PBS and all antibodies were diluted with 1 % (w/v) skim milk in PBS. Each plate included a pool of positive plasmas that was diluted two fold from 1/50 to 1/800 to generate a five point standard curve. The 1/50 dilution of pooled positive plasma was assigned an arbitrary value of 800 relative units (RU). All plasmas were tested at 1/50 and OD values interpolated from the standard curve for that plate. Any plasma that were below the curve were assigned the lowest value, any that were above the curve were re-tested at two-fold dilutions from 1/50 to 1/400. A pool of unexposed donor plasma at the same dilution as the test plasmas was included as a negative control in every Luminex and ELISA assay.
The association between disease severity with age and parasitaemia was assessed using a Mann–Whitney U-test and with gender using a Fisher’s exact test. RU values for individual proteins were compared by Mann–Whitney U-test. To compare between patients with severe and uncomplicated malaria for antibody responses to proteins belonging to a single domain cassette or PfEMP1 group, patients were categorized according to whether their plasma sample lay above or below the median concentration of RU for that antigen: those above or equal scoring 1, or 0 if below. To derive a single quantitative score for each plasma for all the proteins belonging to a single DC or PfEMP1 group, the plasma’s scores for each antigen in that DC or group that were determined by either Luminex or ELISA were summed.
Any individual plasma with a score of 1 for any protein within a DC or PfEMP1 group was classified as a responder to that DC or PfEMP1. Individual plasma samples with a combined score of 0 for all proteins within a DC or group were classified as non-responders.
Differences in the proportions of severe and uncomplicated malaria patients whose plasma responded to a DC or PfEMP1 group were compared by contingency tables using Fisher’s exact test. Differences between severe and uncomplicated malaria patients in the number of proteins within a DC or PfEMP1 group to which patients responded were compared only for patients who responded to at least one protein in the DC or PfEMP1 group using Mann–Whitney U- tests. This indicated differences in the breadth of the response to PfEMP1s within that group. The patients that did not respond to any protein within the group, i.e. had a score of zero were not included to remove any biases associated with large frequencies of zero values in non-parametric comparisons .
Var genes expressed in a Papua New Guinean cerebral malaria patient
To identify var sequences transcribed by parasites infecting this patient we used Illumina RNAseq to generate 65 bp paired-end short reads. Reads that mapped to the 399 full length var gene sequences available  were merged with reads that did not map to the P. falciparum 3D7 strain genome nor Homo sapiens, and the merged reads were subjected to two rounds of de novo assembly to generate 623 contigs that included 362 contigs with homology to var genes. The total length of assembled contigs was 200,158 bp, N50 423 bp, maximum contig length 4083 bp. The short contig length allowed assembly of only a few, full domains. However, the contigs could be aligned by BLASTN (E value <10−5) to individual domains from the 399 full-length var genes and orthologs of high identity to annotated domains and DC types were identified (Table 1).
We compared the percentage of reads from the patient that assembled in transcripts with homology to domain subtypes with the percentage of total var exon 1 sequence in the seven sequenced P. falciparum genomes that each domain subtype represented. DBLδ1, NTSB, CIDRβ1, NTSA and CIDRα3.1 were all abundant transcripts but also constituted a similar proportion of var transcripts as the proportion of total var exon1 sequences they constitute in the seven sequenced genomes (Fig. 1). Therefore, their abundance could represent random var gene transcription. Other abundant transcripts (present at more than 10 % of the level of total DBLδ1) but which were transcribed at more than three times their level of representation in the seven sequenced genomes included DBLγ6, 4, 9, and 10, DBLδ5, DBLβ6 and 7, CIDRα2.4, CIDRb2 and DBLε1 and 6 (Fig. 1).
The abundance of individual transcripts was compared using reads mapped per kb of assembled contiguous transcript (RPK) (see “Methods” section). The most abundant individual transcripts included a diverse range of CIDR and DBL domains. The first and eighth most abundant transcripts were CIDRα2.4 and CIDRα3.1 sequences that could potentially bind CD36  (Table 1). The second and third most abundant transcripts had no previous associations with severe disease or conserved domain cassettes. The fourth and the ninth most abundant were DC5 associated DBLδ5 sequences, 12 of the 13 published DBLδ5 sequences are present in DC5 genes. The 13th and 26th most abundant transcripts were also unique DC5 tandem-domain arrangements of DBLγ17-DBLδ5 and CIDRβ4-DBLβ3, respectively.
The most promising candidates for causing severe malaria are the EPCR-binding CIDRα1 domains of DC8 and DC13 PfEMP1s. The most abundant CIDRα1 capable of EPCR binding  in this patient was a CIDRα1.5-DBLβ7 contig ranked 50th overall by abundance. All described DC8s have a unique DBLα2-CIDRα1.1/1.6/1.8 tandem domain arrangement , but the most abundant DBLα2 was ranked 102nd by abundance. Four of the 50 most abundant transcripts had greatest homology to domains found in DC8 PfEMP1s (Table 1). Thus DC8 var genes may have been abundantly expressed by the parasites infecting this patient however these domains are also found at least as frequently in non-DC8 PfEMP1s. The DC13 is characterized by the tandem array of a DBLα1.7-CIDRα1.4 but neither of these domains were abundantly transcribed in this patient.
Var1csa sequences were the 7th, 12th and 45th most abundant transcripts; expression of this gene was not previously observed to be elevated in severe disease  and it is ubiquitously transcribed [59, 60], atypically late in the cell cycle after transcription of var genes encoding the adhesion phenotype [61, 62]. The 26th, 42nd and 35th most abundant transcripts were two DBLβ3 and a DBLβ5, respectively. DBLβ5 and some DBLβ3 including those in DC4 have been shown to bind ICAM-1 [63, 64]. DBLε10 from var2csa was the 43rd most abundant transcript. No other DCs were identified in the 50 most abundant transcripts.
Patterns of PfEMP1 antibody reactivity in severe and uncomplicated malaria in Papua
Antibody reactivity with PfEMP1 was assessed in plasma from 28 patients with severe malaria (median years 29, IQR 18.5-34; median P. falciparum parasites/µl 41,220, IQR 8260-334,273; 61 % male) and 35 patients with uncomplicated malaria (median years 22.5, IQR 18.0-25.5,; median P. falciparum parasites/µl 27,680, IQR 16,800-52,800; 54 % male). Patients with severe malaria tended to be older than those with uncomplicated malaria (p = 0.0599), but there was no significant difference in P. falciparum density (p = 0.7288) and gender (p = 0.7981). Twenty-three patients with severe malaria had a single diagnostic criterion (WHO) , including five with cerebral malaria, six with jaundice, eight with hyperparasitaemia, three with prostration, and one with acute renal failure. Five patients had two or more manifestations of severe malaria: one patient with jaundice and prostration, one with acute renal failure and acute respiratory distress syndrome, one with jaundice and hyperparasitaemia, and two with jaundice and acute renal failure.
To examine the reactivity of Papuan patient plasma with different PfEMP1 groups and DCs previously associated with severe or uncomplicated disease we tested 10 PfEMP1 DCs and the non-PfEMP1 proteins GLURP and MSP3 by a Luminex multiplexed bead assay; and twelve recombinant, partial PfEMP1 domains by ELISA. The proteins tested included groupings previously associated with severe disease (DC8, DC13, DC5, group A and group B) (Fig. 2). Individual datapoints for all plasmas and all antigens tested are presented in Additional file 1.
The DC domain constructs used in the Luminex assay were derived from the 3D7, ItG and HB3 isolates. Seven of the proteins used for ELISA were derived from sequences transcribed in the PNG cerebral malaria patient described above and included abundantly transcribed (by RPK) representatives of DC5 DBLβ7, DC5 DBLδ5, DC8 DBLγ6 and DBLγ10 (Table 1; Fig. 2b and see “Methods” section for details). These domain subtypes were all transcribed at more than threefold their level of representation in the seven sequenced genomes (Fig. 1). The most abundantly transcribed representative of the conserved NTS-DBLα arrangement was also included as well as a less abundant DBLd1 transcript (protein 7) (Table 1; Fig. 2b). By Luminex assay, levels of antibodies to MSP3 but not GLURP were higher in plasma from patients with uncomplicated malaria than in severe malaria, (Fig. 2a) (median MSP3 22 versus 11 RU, respectively, p = 0.055). The MSP3 data suggests that patients with uncomplicated malaria may have had more prior exposure to P. falciparum infection than the patients with severe malaria.
In the Luminex assay, patients with uncomplicated malaria generally had higher levels of antibody to individual PfEMP1s than patients with severe disease, the greatest difference being in a DC16 PfEMP1 (p = 0.0095) (Fig. 2a). Interestingly DC16 are group A PfEMP1s that have been shown not to be associated with severe disease . For four proteins, there was no significant difference (all p > 0.054) in antibody response between plasma from patients with severe and uncomplicated malaria. These proteins included three severe malaria-associated PfEMP1s, two DC5, and one DC8 (Fig. 2a).
In the ELISA, patients with uncomplicated malaria had significantly higher levels of antibody than patients with severe malaria to two group C PfEMP1s that were from lab isolates and a single group B DBLγ from the cerebral malaria patient (Fig. 2b; p < 0.05). A non-significant trend in the same direction (p < 0.10) was observed for a single DC5 and a single DC8.
The combined results for the ELISA and the Luminex assays revealed that a greater proportion of the patients with uncomplicated malaria than with severe malaria had antibodies to group C PfEMP1s (p = 0.004) and to PfEMP1s that were not DC5, nor 8, nor 13 (p = 0.008) (Fig. 3a). While the proportion of plasma reactive with severe malaria-associated proteins was higher in individuals with uncomplicated compared to severe malaria, the magnitude of the difference was smaller and not statistically significant (all p > 0.194). It is possible that analysing larger numbers of proteins in these groupings on a single platform may have detected differences.
Amongst patients with an antibody response to at least one protein of a cassette or group, plasma from patients with uncomplicated malaria recognized more PfEMP1s overall, and more group A and group B PfEMP1s than did plasma from patients with severe malaria (p = 0.031, p = 0.044 and p = 0.014, respectively, Fig. 3b). Thus although similar proportions of patients with severe and uncomplicated malaria had been exposed to at least one group A or group B PfEMP1 (Fig. 3a), the breadth of responses to the group A and group B PfEMP1s tested was lower in patients with severe malaria.
Sequencing the transcriptome of parasites causing CM in a single patient enabled assembly of a snapshot of the transcribed var repertoire in human malaria. Although entire genes could not be assembled, advances in the phylogenetics of var sequences  allowed the var contigs to be separated into useful classifications. Previous studies that implicated parasites expressing group A PfEMP1s [20, 23–25] and group B PfEMP1s  in severe disease provided no, or minimal, sequence data and were essentially restricted to classifying sequences to the groups defined by the var gene upstream sequences or by sequencing short DBLα tags.
The RNAseq of parasites causing CM in this PNG patient was consistent with previous studies of var genes in pathogenesis and abundantly expressed var genes identified included DC5 and possibly DC8, the DBLβ3 domain subtype and individual domains including CIDRα1.5 and DBLβ3. Recombinant proteins derived from the latter two domain subtypes have been shown to bind EPCR and ICAM-1 respectively [30, 36]. Other abundantly transcribed domain subtypes had not previously been identified in severe malaria. This limited study of a single patient indicates that RNAseq will be useful for identifying quantitative differences between transcribed var genes in severe disease in future studies.
The diversity of the transcribed var repertoire was consistent with a previous report of cerebral malaria in Africa . However, 45 % of the reads that assembled into var transcripts were in the 20 most abundant var contigs that between them represented 27 domains. Thus the quantitative nature of RNAseq revealed a hierarchy of var transcript abundance in this patient’s peripheral blood that would be difficult to detect using the non-quantitative, nested RT-PCR approaches available to this previous study . This suggests that the dominant var transcripts expressed by parasites causing cerebral malaria in a single patient are probably restricted in number.
Patients with uncomplicated malaria more commonly had antibodies to PfEMP1s that were from Group C or were not from DC5 nor DC8 nor DC13 than patients with severe malaria (Fig. 3a). In contrast, similar proportions of patients with severe and uncomplicated malaria had developed antibodies to the severe malaria associated PfEMP1s (group A and B, DC5, 8 and 13) (Fig. 3a), but the breadth of the response to group A and B PfEMP1s was greater in patients with uncomplicated than severe disease (Fig. 3b). Thus susceptibility to severe disease was associated with recognition of a narrower range of group A and B PfEMP1s and to an overall lack of antibodies to group C PfEMP1s and to PfEMP1s that were not DC5 nor DC8 nor DC13.
Overall, the serology findings in this Papuan adult population are consistent with existing models of infection in African children where parasites expressing severe malaria-associated group A and B PfEMP1s infect naive individuals and elicit antibodies [20, 66]. Susceptible, semi-immune individuals have antibody to some group A and B PfEMP1s , but protective immunity correlates with acquisition of antibodies recognizing a broader range of PfEMP1s [5–8]. Parasites expressing uncomplicated disease associated group C PfEMP1s, or PfEMP1s that were not DC5, 8 nor 13, would only dominate infections after parasites expressing severe disease associated PfEMP1s were controlled by acquisition of a broad antibody response. The alternative explanation is that group C PfEMP1s and PfEMP1s that were not DC5, 8 nor 13 were abundantly expressed by parasites causing acute, severe malaria but had not yet elicited antibodies. Although this cannot be excluded it is inconsistent with previous studies of var gene expression in both Africa and PNG [20–26].
In Papuan adults severe malaria is associated with a lack of antibodies to non-DC5 and 8 and group C PfEMP1s in general, and with antibodies to a narrower repertoire of group A and group B PfEMP1s than in patients with uncomplicated malaria. These findings from Papua are consistent with reports from Africa of elevated group A and B var gene expression in severe disease [20, 23–26] and of earlier development of antibodies to group A PfEMP1s in children [27, 29]. This study has also established the feasibility of performing RNAseq on patient isolates to identify expressed var gene sequences.
Baruch DI, Pasloske BL, Singh HB, Bi X, Ma XC, Feldman M, et al. Cloning the P. falciparum gene encoding PfEMP1, a malarial variant antigen and adherence receptor on the surface of parasitized human erythrocytes. Cell. 1995;82:77–87.
Smith JD, Chitnis CE, Craig AG, Roberts DJ, Hudson-Taylor DE, Peterson DS, et al. Switches in expression of Plasmodium falciparum var genes correlate with changes in antigenic and cytoadherent phenotypes of infected erythrocytes. Cell. 1995;82:101–10.
Su XZ, Heatwole VM, Wertheimer SP, Guinet F, Herrfeldt JA, Peterson DS, et al. The large diverse gene family var encodes proteins involved in cytoadherence and antigenic variation of Plasmodium falciparum-infected erythrocytes. Cell. 1995;82:89–100.
Rask TS, Hansen DA, Theander TG, Pedersen AG, Lavstsen T. Plasmodium falciparum erythrocyte membrane protein 1 diversity in seven genomes—divide and conquer. PLoS Comput Biol. 2010;6:e1000933.
Marsh K, Otoo L, Hayes RJ, Carson DC, Greenwood BM. Antibodies to blood stage antigens of Plasmodium falciparum in rural Gambians and their relation to protection against infection. Trans R Soc Trop Med Hyg. 1989;83:293–303.
Newbold CI, Pinches R, Roberts DJ, Marsh K. Plasmodium falciparum: the human agglutinating antibody response to the infected red cell surface is predominantly variant specific. Exp Parasitol. 1992;75:281–92.
Chan JA, Howell KB, Reiling L, Ataide R, Mackintosh CL, Fowkes FJ, et al. Targets of antibodies against Plasmodium falciparum-infected erythrocytes in malaria immunity. J Clin Invest. 2012;122:3227–38.
Marsh K, Howard RJ. Antigens induced on erythrocytes by P. falciparum: expression of diverse and conserved determinants. Science. 1986;231:150–3.
Gupta S, Snow RW, Donnelly C, Newbold C. Acquired immunity and postnatal clinical protection in childhood cerebral malaria. Proc Biol Sci. 1999;266:33–8.
Gupta S, Snow RW. C.A. D, Marsh K, Newbold C. Immunity to non-cerebral severe malaria is acquired after one or two infections. Nature Med. 1999;5:340–3.
Bull PC, Kortok M, Kai O, Ndungu F, Ross A, Lowe BS, et al. Plasmodium falciparum-infected erythrocytes: agglutination by diverse Kenyan plasma is associated with severe disease and young host age. J Infect Dis. 2000;182:252–9.
Nielsen MA, Staalsoe T, Kurtzhals JA, Goka BQ, Dodoo D, Alifrangis M, et al. Plasmodium falciparum variant surface antigen expression varies between isolates causing severe and nonsevere malaria and is modified by acquired immunity. J Immunol. 2002;168:3444–50.
Smith JD, Subramanian G, Gamain B, Baruch DI, Miller LH. Classification of adhesive domains in the Plasmodium falciparum erythrocyte membrane protein 1 family. Mol Biochem Parasitol. 2000;110:293–310.
Rowe JA, Claessens A, Corrigan RA, Arman M. Adhesion of Plasmodium falciparum-infected erythrocytes to human cells: molecular mechanisms and therapeutic implications. Expert Rev Mol Med. 2009;11:e16.
Smith JD, Rowe JA, Higgins MK, Lavstsen T. Malaria’s deadly grip: cytoadhesion of Plasmodium falciparum-infected erythrocytes. Cell Microbiol. 2013;15:1976–83.
Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, et al. Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002;419:498–511.
Lavstsen T, Salanti A, Jensen AT, Arnot DE, Theander TG. Sub-grouping of Plasmodium falciparum 3D7 var genes based on sequence analysis of coding and non-coding regions. Malar J. 2003;2:27.
Salanti A, Staalsoe T, Lavstsen T, Jensen AT, Sowa MP, Arnot DE, et al. Selective upregulation of a single distinctly structured var gene in chondroitin sulphate A-adhering Plasmodium falciparum involved in pregnancy-associated malaria. Mol Microbiol. 2003;49:179–91.
Normark J, Nilsson D, Ribacke U, Winter G, Moll K, Wheelock CE, et al. PfEMP1-DBL1alpha amino acid motifs in severe disease states of Plasmodium falciparum malaria. Proc Natl Acad. 2007;104:15835–40.
Warimwe GM, Keane TM, Fegan G, Musyoki JN, Newton CR, Pain A, et al. Plasmodium falciparum var gene expression is modified by host immunity. Proc Natl Acad Sci. 2009;106:21801–6.
Kaestli M, Cockburn IA, Cortes A, Baea K, Rowe JA, Beck HP. Virulence of malaria is associated with differential expression of Plasmodium falciparum var gene subgroups in a case-control study. J Infect Dis. 2006;193:1567–74.
Falk N, Kaestli M, Qi W, Ott M, Baea K, Cortes A, et al. Analysis of Plasmodium falciparum var genes expressed in children from Papua New Guinea. J Infect Dis. 2009;200:347–56.
Rottmann M, Lavstsen T, Mugasa JP, Kaestli M, Jensen AT, Muller D, et al. Differential expression of var gene groups is associated with morbidity caused by Plasmodium falciparum infection in Tanzanian children. Infect Immun. 2006;74:3904–11.
Kyriacou HM, Stone GN, Challis RJ, Raza A, Lyke KE, Thera MA, et al. Differential var gene transcription in Plasmodium falciparum isolates from patients with cerebral malaria compared to hyperparasitaemia. Mol Biochem Parasitol. 2006;150:211–8.
Warimwe GM, Fegan G, Musyoki JN, Newton CR, Opiyo M, Githinji G, et al. Prognostic indicators of life-threatening malaria are associated with distinct parasite variant antigen profiles. Sci Transl Med. 2012;4:129.
Kalmbach Y, Rottmann M, Kombila M, Kremsner PG, Beck HP, Kun JF. Differential var gene expression in children with malaria and antidromic effects on host gene expression. J Infect Dis. 2010;202:313–7.
Jensen AT, Magistrado P, Sharp S, Joergensen L, Lavstsen T, Chiucchiuini A, et al. Plasmodium falciparum associated with severe childhood malaria preferentially expresses PfEMP1 encoded by group A var genes. J Exp Med. 2004;199:1179–90.
Wang CW, Hermsen CC, Sauerwein RW, Arnot DE, Theander TG, Lavstsen T. The Plasmodium falciparum var gene transcription strategy at the onset of blood stage infection in a human volunteer. Parasitol Int. 2009;58:478–80.
Cham GK, Turner L, Lusingu J, Vestergaard L, Mmbando BP, Kurtis JD, et al. Sequential, ordered acquisition of antibodies to Plasmodium falciparum erythrocyte membrane protein 1 domains. J Immunol. 2009;183:3356–63.
Howell DP, Levin EA, Springer AL, Kraemer SM, Phippard DJ, Schief WR, et al. Mapping a common interaction site used by Plasmodium falciparum duffy binding-like domains to bind diverse host receptors. Mol Microbiol. 2008;67:78–87.
Oleinikov AV, Amos E, Frye IT, Rossnagle E, Mutabingwa TK, Fried M, et al. High throughput functional assays of the variant antigen PfEMP1 reveal a single domain in the 3D7 Plasmodium falciparum genome that binds ICAM1 with high affinity and is targeted by naturally acquired neutralizing antibodies. PLoS Pathog. 2009;5:e1000386.
Turner GDH, Morrison H, Jones M, Davis TME, Looareesuwan S, Buley ID, et al. An immunohistochemical study of the pathology of fatal malaria. Evidence for widespread endothelial activation and a potential role for intercellular adhesion molecule-1 in cerebral sequestration. Am J Pathol. 1994;145:1057–69.
Ochola LB, Siddondo BR, Ocholla H, Nkya S, Kimani EN, Williams TN, et al. Specific receptor usage in Plasmodium falciparum cytoadherence is associated with disease outcome. PLoS One. 2011;6:e14741.
Newbold C, Warn P, Black G, Berendt A, Craig A, Snow B, et al. Receptor-specific adhesion and clinical disease in Plasmodium falciparum. Am J Trop Med Hyg. 1997;57:389–98.
Rogerson SJ, Tembenu R, Dobano C, Plitt S, Taylor TE, Molyneux ME. Cytoadherence characteristics of Plasmodium falciparum-infected erythrocytes from Malawian children with severe and uncomplicated malaria. Am J Trop Med Hyg. 1999;61:467–72.
Turner L, Lavstsen T, Berger SS, Wang CW, Petersen JE, Avril M, et al. Severe malaria is associated with parasite binding to endothelial protein C receptor. Nature. 2013;498:502–5.
Lau CK, Turner L, Jespersen JS, Lowe ED, Petersen B, Wang CW, et al. Structural conservation despite huge sequence diversity allows EPCR binding by the PfEMP1 family implicated in severe childhood malaria. Cell Host Microbe. 2015;17:118–29.
Lavstsen T, Turner L, Saguti F, Magistrado P, Rask TS, Jespersen JS, et al. Plasmodium falciparum erythrocyte membrane protein 1 domain cassettes 8 and 13 are associated with severe malaria in children. Proc Natl Acad Sci. 2012;109:E1791–800.
Avril M, Tripathi AK, Brazier AJ, Andisi C, Janes JH, Soma VL, et al. A restricted subset of var genes mediates adherence of Plasmodium falciparum-infected erythrocytes to brain endothelial cells. Proc Natl Acad Sci. 2012;109:E1782–90.
Claessens A, Adams Y, Ghumra A, Lindergard G, Buchan CC, Andisi C, et al. A subset of group A-like var genes encodes the malaria parasite ligands for binding to human brain endothelial cells. Proc Natl Acad Sci. 2012;109:E1772–81.
Lavstsen T, Magistrado P, Hermsen CC, Salanti A, Jensen AT, Sauerwein R, et al. Expression of Plasmodium falciparum erythrocyte membrane protein 1 in experimentally infected humans. Malar J. 2005;4:21.
Magistrado PA, Lusingu J, Vestergaard LS, Lemnge M, Lavstsen T, Turner L, et al. Immunoglobulin G antibody reactivity to a group A Plasmodium falciparum erythrocyte membrane protein 1 and protection from P. falciparum malaria. Infect Immun. 2007;75:2415–20.
Berger SS, Turner L, Wang CW, Petersen JE, Kraft M, Lusingu JP, et al. Plasmodium falciparum expressing domain cassette 5 type PfEMP1 (DC5-PfEMP1) bind PECAM1. PLoS One. 2013;8:e69117.
Combes V, Taylor TE, Juhan-Vague I, Mege JL, Mwenechanya J, Tembo M, et al. Circulating endothelial microparticles in malawian children with severe falciparum malaria complicated with coma. JAMA. 2004;291:2542–4.
Faille D, Combes V, Mitchell AJ, Fontaine A, Juhan-Vague I, Alessi MC, et al. Platelet microparticles: a new player in malaria parasite cytoadherence to human brain endothelium. FASEB J. 2009;23:3449–58.
Heddini A, Chen Q, Obiero J, Kai O, Fernandez V, Marsh K, et al. Binding of Plasmodium falciparum-infected erythrocytes to soluble platelet endothelial cell adhesion molecule-1 (PECAM-1/CD31): frequent recognition by clinical isolates. Am J Trop Med Hyg. 2001;65:47–51.
Karyana M, Burdarm L, Yeung S, Kenangalem E, Wariker N, Maristela R, et al. Malaria morbidity in Papua Indonesia, an area with multidrug resistant Plasmodium vivax and Plasmodium falciparum. Malar J. 2008;7:148.
Yeo TW, Lampah DA, Kenangalem E, Tjitra E, Price RN, Anstey NM. Impaired skeletal muscle microvascular function and increased skeletal muscle oxygen consumption in severe falciparum malaria. J Infect Dis. 2013;207:528–36.
Tjitra E, Suprianto S, Currie BJ, Morris PS, Saunders JR, Anstey NM. Therapy of uncomplicated falciparum malaria: a randomized trial comparing artesunate plus sulfadoxine-pyrimethamine versus sulfadoxine-pyrimethamine alone in Irian Jaya, Indonesia. Am J Trop Med Hyg. 2001;65:309–17.
Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41:e108.
Crusoe M, Edvenson G, Fish J, Howe Adina, McDonald E, et al. The khmer software package: enabling efficient sequence analysis. Figshare. 2014. doi:10.6084/m9.figshare.979190.
Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012;28:1086–92.
Huang X, Madan A. CAP3: a DNA sequence assembly program. Genome Res. 1999;9:868–77.
Cham GK, Kurtis J, Lusingu J, Theander TG, Jensen AT, Turner L. A semi-automated multiplex high-throughput assay for measuring IgG antibodies against Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) domains in small volumes of plasma. Malar J. 2008;7:108.
Tsuboi T, Takeo S, Sawasaki T, Torii M, Endo Y. An efficient approach to the production of vaccines against the malaria parasite. Methods Mol Biol. 2010;607:73–83.
Tsuboi T, Takeo S, Iriko H, Jin L, Tsuchimochi M, Matsuda S, et al. Wheat germ cell-free system-based production of malaria proteins for discovery of novel vaccine candidates. Infect Immun. 2008;76:1702–8.
Delucchi KL, Bostrom A. Methods for analysis of skewed data distributions in psychiatric clinical studies: working with many zero values. Am J Psychiatry. 2004;161:1159–68.
Robinson BA, Welch TL, Smith JD. Widespread functional specialization of Plasmodium falciparum erythrocyte membrane protein 1 family members to bind CD36 analysed across a parasite genome. Mol Microbiol. 2003;47:1265–78.
Duffy MF, Caragounis A, Noviyanti R, Kyriacou HM, Choong EK, Boysen K, et al. Transcribed var genes associated with placental malaria in Malawian women. Infect Immun. 2006;74:4875–83.
Winter G, Chen Q, Flick K, Kremsner P, Fernandez V, Wahlgren M. The 3D7var5.2 (var COMMON) type var gene family is commonly expressed in non-placental Plasmodium falciparum malaria. Mol Biochem Parasitol. 2003;127:179–91.
Duffy MF, Brown GV, Basuki W, Krejany EO, Noviyanti R, Cowman AF, et al. Transcription of multiple var genes by individual, trophozoite-stage Plasmodium falciparum cells expressing a chondroitin sulphate A binding phenotype. Mol Microbiol. 2002;43:1285–93.
Kyes SA, Christodoulou Z, Raza A, Horrocks P, Pinches R, Rowe JA, et al. A well-conserved Plasmodium falciparum var gene shows an unusual stage-specific transcript pattern. Mol Microbiol. 2003;48:1339–48.
Janes JH, Wang CP, Levin-Edens E, Vigan-Womas I, Guillotte M, Melcher M, et al. Investigating the host binding signature on the Plasmodium falciparum PfEMP1 protein family. PLoS Pathog. 2011;7:e1002032.
Bengtsson A, Joergensen L, Rask TS, Olsen RW, Andersen MA, Turner L, et al. A novel domain cassette identifies Plasmodium falciparum PfEMP1 proteins binding ICAM-1 and is a target of cross-reactive, adhesion-inhibitory antibodies. J Immunol. 2013;190:240–9.
Montgomery J, Mphande FA, Berriman M, Pain A, Rogerson SJ, Taylor TE, et al. Differential var gene expression in the organs of patients dying of falciparum malaria. Mol Microbiol. 2007;65:959–67.
Cham GK, Turner L, Kurtis JD, Mutabingwa T, Fried M, Jensen AT, et al. Hierarchical, domain type-specific acquisition of antibodies to Plasmodium falciparum erythrocyte membrane protein 1 in Tanzanian children. Infect Immun. 2010;78:4653–9.
MFD conceived the study, performed molecular, serological and bioinformatic analyses and drafted the manuscript. RN helped design the study, designed and organized sample collection and performed molecular studies. TT and ET expressed recombinant proteins. BSF and LT organized sample collection and performed molecular studies. JM, DAL, RNP and NMA helped design the study and organized sample collection. FS performed molecular studies. LT, TL and TG helped design the study, expressed proteins and designed serological analyses. FJIF performed statistical analyses. ZPF performed bioinformatic analyses. SJR and GVB helped design the study and draft the manuscript. PS helped organize sample collection. ATP helped design the study and performed bioinformatic analyses. All authors read and approved the final manuscript.
We would like to thank all patients and their families from Papua Province, Indonesia, and Papua New Guinea for participating in the study. We acknowledge Laurens Manning for help in collecting patient samples. This work was supported by National Health and Medical Research Council of Australia [Grant 1007954, Grant 513782, Career Development Fellowship to ATP (1003856) and Practitioner Fellowship to NMA (1042072). The Timika Research Facility and Papuan Community Health Foundation were supported by AusAID (Australian Agency for International Development, Department of Foreign Affairs and Trade) and the National Health and Medical Research Council of Australia (Program Grant 1037304 awarded to RNP and NMA) and the Wellcome Trust (Senior Fellowship in Clinical Science 091625 awarded to RNP). The work benefitted from support by the Victorian State Government Operational Infrastructure Support and Australian Government NHMRC IRIISS.
The authors declare that they have no competing interests.