Previous Article | Next Article ![]()
Infection and Immunity, July 2006, p. 3904-3911, Vol. 74, No. 7
0019-9567/06/$08.00+0 doi:10.1128/IAI.02073-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Thomas Lavstsen,2,
Joseph Paschal Mugasa,3
Mirjam Kaestli,1
Anja T. R. Jensen,2
Dania Müller,1
Thor Theander,2 and
Hans-Peter Beck1*
Swiss Tropical Institute, Socinstrasse 57, Postfach 4002 Basel, Switzerland,1 Centre for Medical Parasitology at Institute for Medical Microbiology and Immunology, University of Copenhagen, Panum Institute 24-2, Blegdamsvej 3, 2200 Copenhagen N, Denmark,2 Ifakara Health Research and Development Centre, P.O. Box 53, Ifakara, Morogoro, Tanzania3
Received 23 December 2005/ Returned for modification 28 January 2006/ Accepted 4 April 2006
|
|
|---|
|
|
|---|
It has been shown that parasites from patients with severe malaria express a different subset of surface antigens that are more frequently recognized by sera from malaria-exposed individuals, including young children, than parasite antigens from older children with mild malaria (5, 23). It has also been shown that this subset of surface antigens is serologically conserved among different geographical regions (24), and it is therefore crucial to identify the molecular phenotype of such a subset to develop a disease-ameliorating vaccine.
Linking PfEMP1 expression to disease outcome is inherently difficult due to the extensive inter- and intragenomic variation in var genes. Previous studies predominantly relied on reverse transcription-PCR using degenerate primers, with subsequent cloning and sequencing (4, 17). While these studies have shown that the transcription of certain DBL1
domains is associated with either severe malaria or rosette formation, they have been unable to identify a clear correlation between var gene groups and disease outcomes.
Sequencing of the 3D7 genome revealed that P. falciparum parasites contain 50 to 60 var genes that can be grouped into three major groups, A, B, and C, and the single-copy intergenomic conserved var genes, var1 and var2csa, according to sequence similarities in both noncoding and coding sequences (12, 18, 22, 37). Evidence is emerging for the existence of subgroups of var group A, namely, type 3 var (18) and type 4 var (12, 15) genes, referred to here as var3 and var4, respectively. The functional relevance of this genetic structuring is indicated by the fact that CIDR domains of group B and C PfEMP1 variants bind to CD36, in contrast to CIDR domains of group A variants (26). Parasites selected for chondroitin sulfate A and human bone marrow endothelial cell binding in vitro dominantly express var2csa and var4, respectively (15, 29), also supporting the notion of functional genetic substructuring.
The genetic organization of var genes was exploited to design primers targeting the conserved regions defining var1-3 genes and group A, B, and C var genes. Using these primers in quantitative reverse transcription-PCR, the transcript abundances of var genes were measured in parasites collected from P. falciparum-infected Tanzanian children with asymptomatic malaria (AM), uncomplicated malaria (UM), and severe malaria (SM). Our data demonstrate an increase in transcript abundance for group A and B var genes in parasites causing severe malaria compared to that in parasites causing uncomplicated malaria.
|
|
|---|
Of all children seeking medical treatment, 40% were admitted to the hospital due to infection with P. falciparum. Malaria is reported to account for a case fatality rate of 2.4% in this hospital (30). Samples were collected from children (aged 4 to 59 months) presenting with malaria at the hospital. Severe malaria cases were defined according to the World Health Organization criteria for severe malaria (38). Uncomplicated malaria was defined as the presence of asexual P. falciparum, an axillary temperature of >37.5°C, or symptoms of headache or myalgia but no other signs of severe malaria. Exclusion criteria were confirmed coinfection, malnutrition (mid-upper-arm circumference [MUAC] of <12 cm), or antimalarial treatment during the last 14 days. Asymptomatic patients (presence of P. falciparum, axillary temperature of <37.5°C, and no other symptoms) were age-matched as closely as possible to the patients with severe cases by convenience sampling in the same area in January 2005. P. falciparum infection was determined by Giemsa-stained blood films, and parasitemia was counted as parasites per 200 white blood cells. Ethical clearance for this study was obtained from the Ifakara Health Research and Development Centre's scientific review board and the Medical Research Coordinating Committee of the National Institute for Medical Research in Tanzania.
Blood samples. After obtaining written informed consent from parents, 1 to 2 ml of venous blood from children was collected in EDTA tubes. Erythrocytes were separated from serum by centrifugation and washed with 40 ml phosphate-buffered saline, 5 volumes of TRIzol reagent (Invitrogen) was added to the erythrocyte pellet, and the sample was frozen at 70°C. Samples were collected from 52 patients with SM, 56 patients with UM, and 19 AM children. 3D7 parasite lines were generated as described elsewhere (33).
DNA, RNA, and cDNA. Genomic DNAs were isolated from infected red blood cells (100 µl blood was frozen after adding 2 volumes of 6 M guanidine HCl, 50 mM Tris, pH 8.0, 20 mM EDTA) with a QIAamp blood kit (QIAGEN), and total RNAs were extracted by using TRIzol reagent (Invitrogen) twice, as recommended by the manufacturer, and treated with DNase I (Invitrogen) for 30 min at 37°C. The absence of DNA in RNA samples was confirmed by stable base fluorescence after 40 cycles of real-time PCR with seryl-tRNA synthetase primers (29). Reverse transcription was performed using Superscript II (Invitrogen) and random hexamer primers in a total volume of 40 µl according to the manufacturer's recommendations. Hereafter, two different real-time PCR methods were used to quantify var transcript abundance, namely, a SYBR green-based assay and a minor groove binder (MGB) probe-based assay.
Validation of SYBR green method and quantification of var gene transcript abundance.
Quantitative real-time PCR using Quantitect SYBR green PCR master mix (QIAGEN) was performed on a Rotorgene thermal cycler system (Corbett Research) as previously described, using the seryl-tRNA synthetase (primer pair p90) and fructose-bisphosphate aldolase (primer pair p61) genes as endogenous controls (23). Primers targeting var gene groups A, B, and C were designed based on sequence similarities in the 3D7 var repertoire (22). The majority of var genes are flanked by conserved upstream region upsA, upsB, or upsC. Primer pairs B1 and B2 target the conserved upstream region of var B genes, whereas C1 and C2 target the upstream region of var C genes. Attempts to design primer pairs targeting upsA with sufficient amplification in the SYBR green assay were unsuccessful. In strain 3D7, grouping into groups A, B, and C is also maintained in the coding sequences for DBL
and ATS. This was exploited to design primer pairs A1, targeting DBL
of group A var genes; A2 and A3, targeting exon 2 of group A var genes; and BC1 and BC2, targeting exon 2 of both B and C var genes. Additional primers were designed to target conserved regions of var1 and var2csa 5' untranslated regions (UTRs) and var3 coding regions (pvar1utr, pvar2utr, and pvar3). Primers designed to target var4 genes in quantitative reverse transcriptase PCR amplified fragments with the expected melting temperature (Tm) for 3/20 genomic DNAs from field samples only, and this analysis was therefore left out subsequently. Seven annotated genes in 3D7 were predicted not to be targeted by any of these primer pairs. Primers are shown in Table 1, and the 3D7 genes expected to be amplified by the respective primers are listed in Table 2. The primers were validated as follows. Initially, all primers were tested on 10-fold dilutions of 3D7 genomic DNA (gDNA). All primers amplified fragments of the expected size and Tm, and sequencing of one PCR clone from each amplification reaction revealed an expected target sequence. All primers had amplification efficiencies (E) between 1.85 and 2 [E = 10(1/slope of 10-fold-dilution gDNA standard curve)] (data not shown). Cycle threshold (CT) values for primers targeting multiple genes were compared to those obtained with the control primer p90, targeting a single-copy gene. The observed
) was compared to the expected value estimated from the number of predicted genes targeted in 3D7. A
CT value reduction of 1 represents a duplication of targeted gene copies, and the estimated
= log(no. of targeted genes)/log(2) (Table 2). Most primer pairs amplified as expected; the primers were then tested on gDNAs from 20 field isolates, and the amplification results were compared to the p90 CT values (
) (Table 2). All primer pairs yielded fragments with the sizes and Tms expected from the amplifications of 3D7 gDNA, and all primers amplified with similar efficiencies (95% confidence interval for µ(n = 20)[
], <0.8) from all field isolates, with no significant difference (t test) between DNAs from severe and uncomplicated malaria cases. The amplification efficiency for field isolate gDNA was similar to that for 3D7 gDNA for most primer pairs, except C1, C2, and BC2, indicating that these primers targeted fewer genes in field isolates than in 3D7 and might reflect that a larger variation in var group C gene copy numbers exists among parasites from naturally infected individuals. However, none of the primer pairs showed any significant difference in amplification efficiencies between DNAs from severe and uncomplicated malaria cases.
|
View this table: [in a new window] |
TABLE 1. SYBR green primers
|
|
View this table: [in a new window] |
TABLE 2. Technical characteristics of var group primers for SYBR green real-time PCRa
|
![]() View larger version (10K): [in a new window] |
FIG. 1. Differences in var transcript abundance (x-fold changes) between 3D7UM and antibody-selected 3D7SM parasites (15). Transcript abundance was measured by using primers targeting var gene groups (white) or primers targeting single var genes (black) and are summarized corresponding to var gene groups. A twofold change in var transcript abundance (dashed lines) was arbitrarily defined as the cutoff for biologically significant changes in var transcript abundance.
|
(µ,
) = 21.2 and 2.6;
(µ,
) = 19.0 and 2.5]. As expected, there was a negative association between parasitemia and control gene p90 CT values (representing the overall amount of cDNA) (regression coefficient[log2(parasitemia/blood sample volume):
], 0.58; P = 0.005).
Validation of MGB method and quantification of var transcript abundance.
Quantitative real-time PCR using MGB probes was performed using an ABI PRISM 7200 sequence detection system (Applied Biosystems) as described by Kaestli et al. (16), with few modifications. Briefly, cDNAs were synthesized from total RNA, and a primary PCR with 16 cycles over the var 5' UTR-DBL1
target sequence was performed prior to real-time PCR for var groups A, B, and C. Primers and probes targeting upsA, -B, and -C for the MGB probe assay were described by Kaestli et al. (16) (upsA-probe, upsB-probe, and upsC-probe). The seryl-tRNA synthetase internal control gene was used for relative quantification without prior amplification. The primers and probe were designed using Primer Express software 2.0 (Applied Biosystems) and had the following sequences: primer p90Probe_for, 5'-ACCTCAGAACAACCATTATGTGCTT-3'; primer p90Probe_rev, 5'-TGTGCCCCTGCTTCTTTTCTAA-3'; and p90Probe, 5'-6-carboxyfluorescein-AGGTTACCACTCAAATACGCTGGATTCTCATCTTG- 6-carboxytetramethylrhodamine-6-carboxyfluorescein-3'.
Data analysis.
After all samples had been subjected to real-time PCR, the data set was cleaned for subsequent statistical analysis. Data points were not considered if the Tm diverged more than 1°C from the expected value or if the CT value was above 30. Transcript abundances were compared between clinical groups after normalization to internal controls (yielding
CT values) (Fig. 2). Based on these, x-fold changes were calculated by the 
CT method (see Table 4). Comparisons between groups were made with one-way analysis of variance and Intercooled Stata 8.0 analysis software.
![]() View larger version (22K): [in a new window] |
FIG. 2. Transcript abundances of var gene groups in parasites from children with AM and from children suffering from UM or SM. Transcript abundances are shown relative to the average abundances in uncomplicated cases ( CT values). Panel A shows var group A transcript abundances measured with primers A1-3 in quantitative PCR and with upsA-probe in the MGB assay. Similarly, panels B and C show the transcript abundances of group B and C var genes, respectively. Panel D shows the transcript abundances measured with primers BC1 and BC2, targeting group B and C genes. Panel E shows transcript abundances measured with primers targeting the conserved var genes var1-3. Boxes outline 25th to 75th percentiles, with medians indicated as a line inside each box and whiskers illustrating the 5th and 95th percentiles. Horizontal lines with asterisks below the plots indicate statistically significant differences in transcript abundance between groups (one-way analysis of variance; P < 0.05 after Bonferroni correction).
|
|
View this table: [in a new window] |
TABLE 4. Comparison of var group transcript abundances in parasites from patients with SM or UM
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Clinical and parasitological details of subjects
|
The BC1 primer pair was predicted to predominantly amplify fragments of B var genes, whereas BC2 primers were expected to amplify a smaller subset of genes consisting of group B, BC, and C var genes (Table 2). The BC1 primer results showed that the targeted genes were transcribed at higher levels by SM parasites than by UM or AM parasites. No significant changes were observed in transcription of genes targeted by the BC2 primer pair.
The primer pairs targeting the var1 and var2csa (Table 1) gene family showed no significant difference in transcript abundance between the cohorts, although a trend of increased var2csa transcript abundance with increased severity of disease was observed. In contrast, primers targeting the var3 family showed significantly higher transcript abundance in SM than in AM samples and a trend of higher transcript abundance with increased severity of disease (Fig. 2). var3 belongs to the group A var genes, and a correlation between the transcript abundances of var3 and group A var genes would be expected. The strongest correlation of var3 transcript abundance was found with transcripts measured by A2 (Rpvar3:A2 = 0.482; P = 0.0012 [Spearman rank]).
High transcript abundances of group A and B var genes are associated with severe disease. The association between transcript abundance and clinical presentation of malaria was tested in logistic regression models in which the dependent variable was the clinical presentation (SM or UM) and the independent variables were age, MUAC, parasitemia, sex, and transcript abundance measured by the respective primer pair. The logistic regression models were built for primer sets exhibiting statistically significant differences in transcript abundance (Table 4) and showed that young age and increased MUAC significantly increased the risk of severe disease. No significant association was found for parasitemia and sex (data not shown). According to the model, the risk of severe malaria is increased 20 to 61% with a twofold increased var group A or var group B transcript abundance (Table 5).
|
View this table: [in a new window] |
TABLE 5. Logistic regression model showing the risk of severe malaria for a twofold increase in transcript abundance of specific var gene groups after correcting for the effects of age, MUAC, parasitemia, and sex
|
3 (data not shown). Parasites from these children showed a trend towards a larger abundance of var group A transcripts than those for all other clinical cases (for upsA-probe, P = 0.1660; for A2 primer, P = 0.0077). Primer pairs A1 and A3 showed no difference in transcript abundance between the groups (data not shown).
Associations of var transcript abundance with other clinical features.
Linear regression models showed that increased var group B transcript abundances measured by both primer pairs B1 and B2 were positively associated with parasitemia. This was also the case when corrected for age [R(parasites/200 leukocytes):
CT(B1/B2) = 226/223; P = 0.045/0.019; an increase of parasitemia by 226 parasites/200 leukocytes resulted in a twofold increase in var group B transcript abundance (1 CT value decrease)]. In contrast, there was a nonsignificant trend for var group C transcript abundance, as measured by C1 or C2 primers, to be negatively correlated with parasitemia [R(parasites/200 leukocytes):
CT(C1/C2) = 224/207; P = 0.091/0.074). There was also no association between var group A transcript abundance and parasitemia, but var group A transcript abundance tended to decrease with age among the UM cases (Rage (months):
CT(A1/A2/A3) = 1.58/1.67/1.72; P = 0.050/0.058/0.081). No such trends were found with var group B or C.
|
|
|---|
PfEMP1 variants of group A were also associated with severe malaria in a study with Brazilian children. The dominant DBL1
transcripts were determined by their amplification frequencies upon reverse transcription-PCR using degenerate primers (17). However, it has been shown to be difficult to reliably reproduce results using this approach (10). Bull et al. (4) used a similar approach to study var gene expression in 12 clinical isolates from Tanzania. Although they were able to identify unique short DBL
sequence markers for var group A that correlated with the formation of rosettes, no association between var group expression and disease outcome was found. Recently, Kaestli et al. (16) analyzed parasites from Papua New Guinean children with asymptomatic, uncomplicated, or severe malaria infections by using quantitative real-time PCR to investigate changes in the proportion of var A, B, or C transcripts, using primers targeting the corresponding upstream regions. This study showed a significant increase in proportions of var group B transcripts in clinical cases, whereas var group C transcript levels were increased in asymptomatic cases. No particular involvement of var group A was reported.
The above-mentioned primers and probes were also applied in the present study, but to allow for relative comparisons to internal control genes, a new set of primers with specificity for the three main groups (A, B, and C) and var1-3 was designed. By using these primers, var transcript levels were measured in cDNAs from 106 Tanzanian children with asymptomatic, uncomplicated, or severe P. falciparum infections.
Similar to the results of the study by Kaestli et al. (16), a larger transcript abundance of var group B genes was observed with an increasing degree of disease severity. Importantly, a similar pattern of transcript abundance was found for the var genes of group A (including var3). Conversely, var genes of group C were found to be transcribed at the same level in all sample groups. The differential transcript abundance patterns determined with the BC1 and BC2 primers corresponded to the transcript abundance patterns measured by B and C group-specific primers, respectively. These conclusions were supported by data generated by two quantitative PCR methods with primers targeting both the 5' and 3' ends of the genes. Although we cannot exclude the possibility that the difference in PfEMP1 expression patterns in asymptomatic children and younger symptomatic children is due to the marked age difference or the severity of infection, the data suggest that var group C genes are not involved in severe childhood malaria. The fact that no var gene group was detected at higher transcription levels in AM and UM than in SM samples is puzzling. This could be explained in several ways. Firstly, the categorization of patients into AM, UM, or SM is operational, and other host and parasite factors, including other variant surface antigens, might play equally important roles in disease outcome. Secondly, since the current knowledge of the global var sequence repertoire is limited, unknown transcripts not targeted by our primers may be excluded from the analysis. Finally, there is a possibility that the present primers result in biased amplification of a subset of predicted target genes. In particular, the last case might be true for the var group C primer pairs C2 and BC2, as indicated by the relatively large differences between estimated and observed CT values for the 3D7 genomic DNA amplifications. Alternatively, group A or B var genes might be expressed in relatively larger abundance than var group C genes. This explanation would assume that not only the type of adhesion ligand, but also the amount of ligand, determines the adhesion phenotype. Future studies will test this hypothesis.
In logistic regression models, there was a statistically significant association between the risk of developing severe malaria and the transcript abundance of group A or B var genes. Thus, a twofold increase in var group A or B transcript abundance was associated with an increase of 20 to 61% in the risk of developing severe malaria. This observation does not indicate any difference between var group A and B transcript abundances in relation to disease severity. However, some indications of functional differences in group A and B var genes may be found in the linear regression models of var transcript abundance and clinical features. The association of increased parasitemia with var group B transcript abundance suggests that the parasites expressing these genes caused the severe infections in these children. In contrast, the lack of association between var group C transcription and parasitemia or clinical presentation supports a previous finding which suggested that group C var genes may be involved in establishing chronic infections (16). Similar to the case for var group B transcripts, a positive correlation between var group A transcript abundance and parasitemia would be expected in these data and in previous findings indicating the involvement of group A var genes in severe malaria (3, 15, 17). The lack of such an association might be explained if group A PfEMP1 variants confer the strongest cytoadhesion in naïve individuals only and if group A variants only dominate in first malaria infections. Since the average age of the children enrolled in this study was 29 months, most children would have undergone several, sometimes severe, infections and would have developed some immunity against PfEMP1 variants of group A. This is supported by a trend towards a lower var group A transcript abundance with increasing age in UM cases.
The observed trend of more abundant var group A transcripts in cerebral malaria cases than in all other cases might indicate that these var genes play a specific role in cerebral malaria. var2csa has been identified as the main chondroitin sulfate A binding ligand in pregnancy-associated malaria (28). The difference seen in var2csa transcript abundance between the disease groups was therefore unexpected. However, for all samples, the var2csa transcript abundance was >100-fold lower (data not shown) than that reported for placental parasites (35) or parasites selected in vitro on CSA (29). In addition, while var2csa transcription appears to be controlled by similar mechanisms to those controlling group A var genes (11), the translation of var2csa transcripts, unlike that of other var transcripts, seems to be controlled by translation of an upstream open reading frame (22; K. W. Deitsch, personal communication). Thus, var2csa is most likely not responsible for the disease outcomes of these children.
var1-like genes are unique among var genes because they are highly conserved between parasite genomes and appear to be controlled by a unique 5' region (36), which might indicate a specialized function for var1 products similar to the var2csa gene in pregnancy-associated malaria. However, the similar abundances of var1 transcripts in all three groups gave no indications of the function of var1 products. This, together with the observed constitutive var1 transcription throughout the intraerythrocytic stages in isogenic but phenotypically different parasite lines (20, 21, 29), leaves the function of var1 molecules enigmatized.
In conclusion, the data presented here show an association between disease outcomes and the transcription of var subtypes in African areas where malaria is endemic. Of specific importance, the association between severe malaria in young children and var group A and B transcription is demonstrated and supports the notion that a vaccine based on selected PfEMP1 molecules might be feasible.
This study was supported by the Swiss National Science Foundation (grant number 3100 AO 104043/1) and by grants from the Danish Medical Research Council (no. 22-02-0220) and the Commission of the European Communities (grant no. QLK2-CT-1999-01293) (EUROMALVAC).
M.R. and T.L. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»