- Open Access
High sensitivity detection of Plasmodium species reveals positive correlations between infections of different species, shifts in age distribution and reduced local variation in Papua New Guinea
Malaria Journal volume 8, Article number: 41 (2009)
When diagnosed by standard light microscopy (LM), malaria prevalence can vary significantly between sites, even at local scale, and mixed species infections are consistently less common than expect in areas co-endemic for Plasmodium falciparum, Plasmodium vivax and Plasmodium malariae. The development of a high-throughput molecular species diagnostic assay now enables routine PCR-based surveillance of malaria infections in large field and intervention studies, and improves resolution of species distribution within and between communities.
This study reports differences in the prevalence of infections with all four human malarial species and of mixed infections as diagnosed by LM and post-PCR ligase detection reaction – fluorescent microsphere (LDR-FMA) assay in 15 villages in the central Sepik area of Papua New Guinea.
Significantly higher rates of infection by P. falciparum, P. vivax, P. malariae and Plasmodium ovale were observed in LDR-FMA compared to LM diagnosis (p < 0.001). Increases were particularly pronounced for P. malariae (3.9% vs 13.4%) and P. ovale (0.0% vs 4.8%). In contrast to LM diagnosis, which suggested a significant deficit of mixed species infections, a significant excess of mixed infections over expectation was detected by LDR-FMA (p < 0.001). Age of peak prevalence shifted to older age groups in LDR-FMA diagnosed infections for P. falciparum (LM: 7–9 yrs 47.5%, LDR-FMA: 10–19 yrs 74.2%) and P. vivax (LM: 4–6 yrs 24.2%, LDR-FMA: 7–9 yrs 50.9%) but not P. malariae infections (10–19 yrs, LM: 7.7% LDR-FMA: 21.6%). Significant geographical variation in prevalence was found for all species (except for LM-diagnosed P. falciparum), with the extent of this variation greater in LDR-FMA than LM diagnosed infections (overall, 84.4% vs. 37.6%). Insecticide-treated bednet (ITN) coverage was also the dominant factor linked to geographical differences in Plasmodium species infection prevalence explaining between 60.6% – 74.5% of this variation for LDR-FMA and 81.8% – 90.0% for LM (except P. falciparum), respectively.
The present study demonstrates that application of molecular diagnosis reveals patterns of malaria risk that are significantly different from those obtained by standard LM. Results provide insight relevant to design of malaria control and eradication strategies.
Malaria is an infectious disease that shows considerable spatial heterogeneity on global [1, 2], regional [3–6] and local scales [7–10]. While it has been possible at some levels of resolution to relate variation in prevalence or incidence of infections with differences in environment and measures of transmission potential (such as entomological inoculation rates) on global and regional levels [4, 6, 11], it has been difficult to evaluate Plasmodium species prevalence and determine the causes of the often pronounced local heterogeneity in malarial prevalence .
As all human Plasmodium species share both host and vector, mixed species infections are commonly observed. An extensive review of the literature on mixed species malaria infections by Ritchie  concluded that there are geographic differences in the way that human malaria species interact and that these interactions may even change from year to year for a given location. In selected field surveys from areas with a Plasmodium falciparum, Plasmodium malariae and Plasmodium vivax mix, a deficit of mixed infections has generally been recorded [14, 15], suggesting heterologous suppression between those parasite species. In contrast, a surplus of mixed infections has been observed in areas with P. falciparum, P. malariae and Plasmodium ovale [14, 15]. However, there is generally great variability in results between studies and it cannot be ruled out that observed fluctuations in the frequency of mixed infections may be due to limitations in the sensitivity of light microscopic (LM) species identification . More recent studies using PCR techniques found much higher levels of mixed infections  and random distribution of different infections [18, 19]. Earlier studies of mixed infections in mostly asymptomatic people during cross-sectional surveys has not provided conclusive evidence for interactions between the different Plasmodium species.
Previous studies on spatial patterns of malarial infections were based exclusively on light microscopy diagnosis (LM) of infection, which has a significantly lower sensitivity for detecting malarial infections compared to PCR-based diagnostic assays [17, 20]. The more limited sensitivity is a particular problem in areas such as Papua New Guinea (PNG) where four human malaria species are endemic (P. falciparum, P. vivax, P. malariae and P. ovale), and mixed species infections are very common . In these circumstances LM has been observed to be particularly poor in detecting the less common P. malariae and P. ovale infections  as well as for accurate diagnosis of mixed species infections .
The developments of a high-throughput post-PCR, LDR-FMA species typing assay  now makes feasible the routine molecular diagnosis of malaria infections in large field and intervention studies. This assay evaluates P. falciparum, P. vivax, P. malariae and P. ovale simultaneously and is performed in 96-well plate format to ensure efficient sample processing. Overall high specificity and sensitivity of the assay provides more accurate assessment of minority species in mixed infections. In addition, good correlations between parasite density and median fluorescence intensity confirmed that the assay is semi-quantitative [23, 24]. In a recent comparison of malarial age group-based infection patterns detected by LM vs LDR-FMA methods in samples of over 1000 people (all ages), an increase in the prevalence of non-falciparum compared to falciparum infections, and a significant shift in distribution of infections to older age groups by LDR-FMA compared to LM diagnosis was observed . Moreover, we found that adolescents and adults very commonly harbour sub-microscopic malarial infections. Similarly, in a cohort of children 5–14 years we found a significantly lower incidence of new P. vivax compared to P. falciparum blood stage infections when diagnosed by LM but identical rates when PCR diagnosis was used .
It is, therefore, becoming clear that increased application of PCR-diagnosis significantly alters the understanding of the epidemiology of different Plasmodium species. In order to further evaluate the effect of the increased sensitivity of detection on patterns of Plasmodium species infection, here LM and LDR-FMA diagnostic methods were compared in a series of 15 cross-sectional population surveys at five distinct locations in the central Sepik region of East Sepik province, PNG.
Study site and design
The selected study sites in the Maprik and Wosera-Gawi districts cover a combined area of over 160 km2 that is characterized by low hills, plains and riverine plains with a wet tropical climate . The natural vegetation is lowland hill forest that has mostly been replaced by re-growth following cultivation. Extensive grasslands are common on the riverine plans near the Sepik River. The Wosera study site described previously  is situated near the center of this expanded study region. In the early 1990's malaria was found to be holoendemic with a peak prevalence of 77% in children 5–9 years of age in the central Wosera area , with entomological inoculation rates (EIR) estimated at 35.4 infective bites/person-year for P. falciparum, 12.1 for P. vivax, and 9.6 for P. malariae . By comparison, recent human studies in the Wosera found prevalence of P. falciparum to be significantly reduced .
For the present study five distinct areas surrounding local health facilities (i.e. Brukham, Burui, Ilaita, Ulupu and Wombisa HC) were selected (Figure 1). Brukham, Ilaita and Ulupu are situated in the foothills of the Prince Alexander range, Burui and Wombisa are in the Sepik River plain. In April and May 2005, 3 villages were surveyed in each of the five areas. In order to achieve a near-random sample of the village population a household-based sampling strategy was pursued. A number of households with a total population of 150–200 persons were included in each survey. The surveys included every member of the selected households who could be reached on the day of survey. Individual informed consent was obtained from all study participants (or consent granted by parents or guardians for inclusion of children) using study protocols evaluated and approved by the PNG Medical Research Advisory Committee (MRAC 00.26 & MRAC 05.20).
From each household a semi-structured questionnaire was administered to collect data on type of house, household assets, education of the parents, personal use of bednets, recent health facility attendance and use of antimalarial drugs. From each individual, a thick and thin film was prepared on a single slide and a 250 μl blood sample collected into a K+EDTA microtainer from a finger or heel prick. Haemoglobin measurements were made from this blood sample using a HemoCue 201+ Hb meter (Angholm, Sweden); remaining blood was preserved for extraction of DNA.
A socio-economic index was created using data on household ownership of consumer durables (i.e. bed, mattress, bednet, chairs, umbrella, clothing cupboard, kerosene pressure lamp, kerosene cooker, electric torch, radio, television, car). Households that owned < five consumer durables were classified as "low", > five consumer durables were classified as "medium", those with > 10 as "high" socio-economic status.
Blood smear preparation and examination
All smears were stained with 2.5% buffered Giemsa (pH 7.2) for 35 minutes and examined by LM. Slides were declared negative if no parasites were seen in 100 thick film fields. The parasite species in positive films were identified and densities were recorded as the number of parasites/200 WBC. Densities were calculated assuming 8,000 WBC/μl . Sexual and asexual stage parasite densities were reported separately for P. falciparum only. All slides with densities less than 200/μl, along with a randomly selected 20% of all blood films were routinely re-examined. If less than 80% concordance was achieved between LM evaluations, the entire batch of slides was re-read.
DNA extraction and molecular diagnostic assay
Parasite DNA was extracted from the cell pellets using QIAmp 96 DNA Blood kits (Qiagen, CA).
A semi-quantitative post-PCR, ligase detection reaction/fluorescent microsphere assay (LDR-FMA)  was used to determine the presence of infection by P. falciparum, P. vivax, P. malariae and P. ovale. The design and sensitivity of this assay has been described previously [23, 24, 28]. In short, this assay combines PCR amplification of the 18S ribosomal RNA gene (491–500 bp fragments) using genus specific primers, followed by a multiplex species-specific ligation detection reaction (LDR). The LDR products are hybridized to Luminex® FlexMAP™ "classification" bead sets (5'), and receive "reporter" labelling following incubation with streptavidin-R-phycoerythrin that binds to biotin (3'). Doubly labelled species-specific LDR complexes are detected using a Bio-Plex array reader (Bio-Rad Laboratories, Hercules, CA). Species-specific fluorescence data were collected by Bio-Rad software, Bio-Plex Manager 3.0 (Bio-Rad Laboratories, Hercules, CA). In order to guarantee maximum sensitivity for the detection of Plasmodium infections the PCR-cycle number was set at 35. Differentiation of negative from positive fluorescent signals was performed by comparing median fluorescent intensity (MFI) from study participants against values obtained from two negative controls on each 96-well plate. Cut-off values for positivity were set at the two standard deviations above the mean MFI of negative controls.
Differences in prevalence of malarial infections among villages and areas and univariate association between categorical explanatory variables were tested using χ2 test, multivariate predictors of infection risk were determined using logistic regression. Non-parametric Spearman's correlations were used to compare parasite density by light microscopy and mean fluorescence intensities, while parametric correlations were calculated for the association of prevalence of infection and mean haemoglobin at village level.
Within and between village variance in risk of malarial infections were estimated using a hierarchical Bayesian logistic model, i.e.
ri ~ Bernoulli(pi)
where logit (pi) = βa(i) + σj(i) + θj(i)k(i).
and "i" indexes the individuals, "a" the age-groups, "j" the area and "k" the village within area. "ri" indicates whether individual "i" is infected with a given Plasmodium species, "pi" denotes expected value, "βa(i)" is a fixed effect for the age group, "σj(i)" is a random effect for the area and the "θj(i)k(i)" is a random effect term for the village. In the models adjusting for bednet coverage, an additional linear effect is introduced into the model. All models were fitted using WinBugs 1.4  assuming normal priors, with the distributions of the random effects centred on zero.
Variance components were estimated for all species with and without adjusting for bednet coverage. The relative contribution of between-area variance was determined by dividing the area variance component with the total spatially structure variance, i.e σj(i)/(σj(i)+ θj(i)k(i)). The effect of difference in bednet coverage on spatially structured variance was determined by comparing the individual variance component from models with and without inclusion of the bednet covariate. Variance components, their relative sizes, and changes following adjustment for bednet coverage are all reported on a logit scale.
A total of 2,744 volunteers from 15 villages in five distinct geographical areas (Figure 1) participated in the cross-sectional surveys. Of these 121 (4.4%) were excluded because of lacking demographic or LM data, while insufficient finger-prick blood sample for LDR-FMA analysis led to exclusion of an additional 96 (3.5%) individuals. Overall 2,527 participants from 659 households with completed demographic information (median participation per household: 4, range (1–17)) were available for study. Of these 1,331 (52.7%) were female, 982 (38.9%) children < 10 yrs and 1,189 (47.1%) adults > = 20 yrs of age.
In this population, 1,121 (44.4%) and 1,844 (73.0%) participants had LM and LDR-FMA detectable Plasmodium species infections, respectively (Table 1). Diagnosis of infections by LM and LDR-FMA were highly significantly associated (X2 > = 493.6, p < 0.001) for all species except P. ovale which is rarely detected by LM in PNG . While the higher sensitivity of the LDR-FMA assay is reflected in a large number of infections that were only detected by this assay (P. falciparum: 667, P. vivax 561, P. malariae 251, P. ovale 121), very few LM-positive infections were not detected by the LDR-FMA assay (P. falciparum: 46, P. vivax 28, P. malariae: 12, P. ovale: 0). Fluorescence intensities of samples that were positive by both LM and LDR-FMA were significantly higher than those that were only LDR-FMA positive (P. falciparum:11,963 vs 4,970; P. vivax: 10,718 vs 5,109; P. malariae: 10,630 vs 6,422; all t-tests p < 0.001). For infections positive by both methods, fluorescence intensities were significantly correlated with densities by LM for P. falciparum (rho = 0.51, p < 0.001) and P. vivax (rho = 0.41, p < 0.001) but not for P. malariae (rho = -0.05, p = 0.67).
By LM, P. falciparum was the most common parasite with a prevalence of 29.9% (756/2,527) followed by P. vivax and P. malariae with a prevalence of 14.6% (368/2,527) and 3.9% (99/2,527), respectively (Table 1). No infections with P. ovale were detected by LM. Plasmodium infections were most frequent in the 7–9, and 4–6 year-old children (67.2%, CL95: 61.2–72.8 and; 60.4%, CL95: 54.9–65.8) respectively, followed by adolescents (10–19 years: 55.9%, CL95: 50.6–61.1), 2–3 year olds (54.7%, CL95: 47.7–61.5), adults (20–39 years: 33.9%, CL95: 30.5–37.4, > = 40 years: 29.7, CL95: 25.5–34.3)) and infants (0–1 years: 26.4%, CL95: 20.1–33.4). P. vivax reached peak prevalence in younger children (Figure 2: 4–6 years: 24.2%, C95: 22.1–26.3) than P. falciparum (7–9 years: 47.5%, CL95: 41.1–53.7) or P. malariae (10–19 years: 7.7%, CL95: 3.7–8.8).
A significantly higher proportion of individuals were positive for infection when diagnosed by LDR-FMA (Table 1). This was the case for all species but was most pronounced for the less common species; P. falciparum prevalence increased 1.8 fold (29.9% to 54.9%, p < 0.001), P. vivax 2.4-fold (14.6% to 35.7%, p < 0.001) and P. malariae 3.4-fold (3.9% to 13.4%, p < 0.001), while P. ovale was detected in samples from 4.8% (p < 0.001) of all participants. In parallel with these increases from LM to LDR-FMA, more mixed infections were detected (Table 2, 4.0% to 28.6%, p < 0.001) and mixed infections were detected in a greater proportion of the infected samples (9.0% to 39.2%, p < 0.001).
As in earlier studies [18, 19, 24], we used the multiple-kind lottery model to determine whether there was a random distribution of Plasmodium species in the survey population (Table 2). Among both LM and LDR-FMA diagnosed infections a significant departure from random distribution was observed, albeit in different ways. By LM diagnosis single species infections were observed at significantly higher frequencies than expected, and different mixed Plasmodium species infections were observed at frequencies that were equal to, or lower than the expected frequencies. Overall, mixed infections were found in only 4.0% of all samples compared with the expected 5.8% (p < 0.001). Fewer mixed species infections were observed in Burui and Wombisa compared to the other areas (Additional files 1 &2). In a given individual, infection with one species was negatively associated with the presence of any other species (Table 3). By LDR-FMA diagnosis, single and double species infections were observed at frequencies equal to or less than expected, whereas triple and quadruple infection assemblages (Table 1: 6.5% vs 3.8%, p < 0.001) and negative samples (27.0% vs 23.9%, p = 0.002) were significantly more common than expected. Among LDR-FMA diagnosed infections, however, positive associations were observed for all species pairs except P. falciparum/P. vivax (Table 3) in any given individual. While the general trend of negative associations between species in LM, and positive associations in LDR-FMA diagnoses are consistent across different age groups, the effects tended to be more pronounced in adults and adolescents (> 10 yrs) compared to children, except for LDR-FMA diagnosed P. falciparum, P. malariae and P. ovale which show a stronger association in children (Table 3).
For all Plasmodium species, the increase in apparent prevalence of infections detected by the LDR-FMA assay was proportionally larger in adults than in children (Figure 2). As a consequence the prevalence of infection appeared to peak at a later age for LDR-FMA-detectable infections than when diagnosis was performed by LM (Figure 2). Prevalence of infections peaked in adolescents (10–19 years: P. falciparum: 74.2%, CL95: 69.3–78.6; P. malariae: 21.6%, CL95: 17.5–26.3; P. ovale: 7.0%, CL95: 4.6–10.2; mixed infections: 43.0%, CL95: 37.9–48.3) in all species except P. vivax, where infections were most commonly found in children 7–9 years of age (50.9%, CL95: 44.8–57.1).
The risk of infections with malaria parasites was independently associated with a number of risk factors. The most important predictors were coverage of insecticide treated ITNs, educational level of the female household head and socio-economic status of family (Table 4). ITN coverage was significantly associated with a reduction in risk of infection with all different malaria species both by LM and LDR-FMA, although for P. vivax a decrease was only observed at high ITN coverage (> 90%). People living in households where the mother spoke English or had attended at least 3 years of schooling also had a significantly lower risk of being infected with P. falciparum, P. vivax (LDR-FMA diagnosed only) and P. malariae infections but not with P. ovale. Only risk of infection with P. falciparum was significantly reduced in households of higher socio-economic status, while difference in house construction was only associated with difference in risk of LDR-FMA diagnosed P. falciparum and P. malariae infection.
For all species and both diagnostic methods, there were significant differences in prevalence of Plasmodium species infections among villages (Table 1, p < 0.001 for all species, except P. ovale (LDR-FMA) p = 0.03), and areas (p < 0.001 for all species, except P. falciparum (LM) p = 0.06) and P. ovale (LDR-FMA) p = 0.01). By LM, infection prevalence for individual villages varied between 15.5% to 45.2% for P. falciparum (coefficient of variance (CV): 26.4%), 6.4% to 25.1% for P. vivax (CV: 39.8%), 0.6% to 15.4% for P. malariae (CV: 94.5%) and 32.0% to 67.6% for any Plasmodium infections (CV: 21.9%). P. vivax and P. malariae infections were found more commonly in the villages in the foot hills (Ulupu (P. vivax only), Brukham & Ilaita) compared to the villages situated on the Sepik River flood plain. The prevalence of LDR-FMA detectable infections were significantly associated with those by LM (P. falciparum: r2 = 0.65, p = 0.008, P. vivax: r2 = 0.67, p = 0.007, P. malariae: r2 = 0.85, p < 0.001). Although variation in infection prevalence between villages was less pronounced (P. falciparum: 35.0% to 78.7%, CV: 17.8%); P. vivax: 22.9% to 45.2%, CV: 16.8%; P. malariae: 2.5% to 30.9%, CV: 55.3%); P. ovale: 1.2% to 9.0%, CV: 45.3% and any infection: 59.1% to 85.6%, CV: 17.8%), broadly similar geographical trends were observed. Interestingly, prevalence of any malaria infection by both LM and LDR-FMA was significantly associated with a lower population mean Hb level (LM: r2 = -0.69, p = 0.004; LDR-FMA: r2 = -0.75, p = 0.001) and a higher incidence of moderate-to-severe anaemia (Hb < 8 g/dl, LM: r2 = 0.61, p = 0.016; LDR-FMA: r2 = 0.65, p = 0.009). When diagnosed by LM, the reductions in age- and gender-adjusted Hb levels were most strongly associated with the presence of a mixed infection (-1.2 g/dl, CL95: 0.9–1.4, p < 0.001), followed by P. malariae (-0.6 g/dl, CL95: 0.3–0.9, p < 0.001), P. falciparum (-0.5 g/dl, CL95: 0.4–0.7, p < 0.001) and P. vivax single infections (-0.3 g/dl, CL95: 0.1–0.5, p = 0.005). Among LDR-FMA diagnosed infections only mixed species (-0.4 g/dl, CL95: 0.2–0.5, p < 0.001) and P. falciparum single infections (-0.5 g/dl, CL95: 0.3–0.7, p < 0.001) were associated with a significant reduction in Hb levels. Hb reductions were larger in children than adults and adolescents.
The Bayesian estimation of among-area and village-within-area variance components revealed differences in scale of this geographical variation in risk of infection among the different Plasmodium species and between methods of diagnosis (Table 5). Among LM diagnosed infection, area-specific differences accounted for only 11.5% of variation in rates of infection with P. falciparum, but for 54.9% and 53.7% of variation in rates of P. vivax and P. malariae infections respectively. With LDR-FMA the proportion of location-structured variation increased to 53.8% for P. falciparum, 88.1% for P. vivax and 83.7% for P. malariae, 86.4% of variation in rates of P. ovale infections was observed between sites.
The prevalence of Plasmodium species infections was significantly reduced in villages with higher ITN coverage (Figure 3). Such a negative association was observed for all species when diagnosed by LDR-FMA (r2: -0.59 to -0.82, p < 0.05) as well as P. vivax (r2: -0.69, p < 0.01) and P. malariae (r2: -0.77, p < 0.001) but not P. falciparum (r2: -0.15) when diagnosed by LM. Accounting for differences in bednet coverage among villages consequently reduced total spatial variance by 56–85% for all species and both diagnostic methods, except for LM-diagnosed P. falciparum infections (-2%). In contrast, substantial differences were observed between the two diagnostic methods and between Plasmodium species in regard to the spatial level at which the effect of bednet coverage was most pronounced (Table 5). Among LDR-FMA diagnosed infections, bednet coverage was associated with a comparable reduction in among-area and village-within-area variances of P. malariae and P. ovale prevalence rates, but a stronger reduction in among-area variance was observed for P. falciparum and P. vivax. Among LM diagnosed P. malariae infections, reduction in among-area variance were comparable to those observed at village-level, while the increased reduction in among-area variance of P. vivax infection was even more pronounced than among LDR-FMA diagnosed infection. Adjusting for bednet coverage further increased the differences in the proportion of spatial variance observed at among-area level variance observed in LM (13–41%) compared to LDR-FMA diagnosed infection (34–90%).
With over 2,500 samples processed, the present study is the largest epidemiological study to-date that has used a post-PCR LDR-FMA assay for the diagnosis of malarial infection on all samples collected. The direct comparison of LM and LDR-FMA results confirmed that the LDR-FMA assay is both significantly more sensitive than LM and semi-quantitative [23, 24], with the increased sensitivity most notable for P. malariae and P. ovale infections.
The higher sensitivity of the LDR-FMA in detecting low density infection was most notable in the assessment of mixed species assemblages and in age-specific prevalence. In concordance with the observations by Kasehagen et al. , a significant deficit of mixed species infections was observed in LM diagnosis, but an excess of high-level mixed infections, and interestingly also of uninfected study participants, were detected by LDR-FMA diagnosis. Prevalence of infections detected by LDR-FMA peaked at a later age than did infections diagnosed by LM (Figure 2), confirming that adolescents and adults contribute significantly more to the total burden of infections in a population than appreciated when only LM diagnosis is used . Having observed the same patterns of infections in two independent surveys suggests that similar changes may be generally observed when PCR-based diagnosis is used in epidemiological studies.
Associations between different Plasmodium species have long been a subject of investigation, but results have been inconsistent and inconclusive (see the review by Richie ). Since the vectors of different Plasmodia are the same, it may be expected that the different parasites would be associated, unless there is immune-mediated competition, species-specific difference in red blood cell susceptibility. Meta-analysis of 19 studies by McKenzie and Bossert  found a tendency for statistically significant surplus of mixed infections by P. falciparum-P. malariae-P. ovale and for deficits in P. falciparum-P. malariae-P. vivax. The present analyses suggest that these deficits may be artifacts of the lower sensitivity of LM, since they were seen in the PNG data only when the assessment was made by LM. This is further highlighted by the weaker negative association in children under 10 that have a lower immune status and thus a more limited control of parasitaemia . A recent series of cross-sectional studies in Malawi  that used PCR for diagnosis also found an significant increase of mixed infection that we more pronounced in adults than children. These observations are consistent with a suggestion  that parasite densities may be suppressed by co-infections, but that co-infections do not prevent blood stage infections, or stimulate cross-species immunity capable of limiting infections.
With both diagnostic methods considerable variation in prevalence of infections was observed among villages for all Plasmodium species. While prevalence by LM and LDR-FMA was strongly correlated and overall geographical patterns comparable, LM-based prevalence showed a much higher degree of variability within individual areas for all Plasmodium species than LDR-FMA-based prevalence. In particular, for the non-falciparum species, almost all spatial variation in LDR-FMA-based rates of infection was observed as difference between areas. The proportionally higher within area variance is likely to be related to the lower sensitivity of LM (both for low-level infection and correct identification of infecting species). Differences in the quality of different slide batches, or between-reader variation may further contribute to within-area variation. In intervention trials where more than one slide is made and all slides have at least two independent reads these types of variability may be reduced. In large scale epidemiology studies or in monitoring of malaria interventions where thousands of people are screened, this is not necessarily common practice. If the trend for a decrease in survey-to-survey variation in PCR-based diagnosis holds true elsewhere, the heterogeneity of malaria risk at the local level may have been overestimated in earlier studies. Results of this nature would emphasize the importance of implementing molecular diagnostic methods.
Interestingly, the least amount of between area variation (by both LM and LDR-FMA) was observed in the most common parasites, i.e. P. falciparum. Whether this reflects a true difference in ecology of different Plasmodium species in the study area is difficult to determine from a single cross-sectional survey. In addition, multiple P. falciparum strain infections are commonly observed in PNG  and the true amount of geographical difference in transmission may be obscured if only presence or absence, but not the number of infecting strains were detected.
The coverage of ITNs in a given village was the best predictor of risk of malaria infections at individual and village levels. Both high coverage of ITNs in the village of residence as well as personal ITN use significantly reduced a participant's risk of being infected with any of the human malaria species. Interestingly, personal use was less strongly associated with protection (data not shown) than overall coverage in a given village. This suggests that if high coverage is achieved, ITN may have an effect on overall transmission and thereby also protect people who do not always sleep under a ITN . The protective effect of ITNs against P. falciparium is well described  and there increasing evidence for their effectiveness against P. vivax, although the effect may be smaller than that against P. falciparum [35–37]. To date there has been very little information on the effect of ITNs on P. malariae and P. ovale. While ITNs were highly effect in reducing P. malariae in a study in Burkina Faso, this may be the first description of an effect of ITNs on risk of P. ovale infections.
Differences in ITN coverage among villages also reduced total spatial variance of Plasmodium infections (measured on the logistic scale) by up to 85%. Similar spatial associations of coverage in untreated bednets with reduced (peak) prevalence and increased age at peak prevalence were observed P. falciparum, P. vivax and P. malariae in earlier studies in the Wosera . In addition, high coverage of untreated bednets led to significant reduction in sporozoite rates in mosquito vectors . In African studies high ITN coverage reduced P. falciparum transmission by up to 90% [11, 40]. Therefore, it is very likely that if sufficiently high coverage rates can be achieved, the ongoing PNG national bednet distribution program will not only reduce the burden of illness, but also have an impact of prevalence of infections. The very high reductions in between-area variation in P. malariae and P. ovale prevalence suggest that ITNs may be quite effective for the control of these species, while control of P. vivax may be more difficult as indicated by the smallest reduction in spatial variation (by LDR-FMA). Interestingly, the variation in prevalence of infection was closely linked to a difference in mean haemoglobin levels in the study villages, indicating that control of malaria through ITNs may also result in improved haemoglobin levels. Given the well known problems with spurious correlations in ecological comparisons, these associations with ITN coverage need to be interpreted with care. Direct evaluation of the ongoing PNG national bednet distribution program on the transmission of non-falciparum species is needed.
The only other factor consistently associated with a reduced prevalence of Plasmodium infections is the level of maternal education. It is well know that educated mothers have a better understanding of health, seek treatment more promptly, and use ITN more regularly [41–43] than uneducated mothers [44, 45]. Similarly, higher socioeconomic status is generally associated with better health outcomes. Interestingly, there was no consistent association of house-type with prevalence of malaria infections of any species. This suggests that different housing types are similarly attractive to the local mosquito vector species. The tendency of many PNG anophelines to bite outdoors readily in the early or late hours of the night  may further reduce any effect of different house types on risk of Plasmodium species infection.
The present study demonstrates that application of the high-throughput, PCR-based LDR-FMA diagnostic assay is feasible in large malaria field studies. Additionally, results show that epidemiological patterns of malaria risk detected by molecular diagnosis differ significantly from risk assessed by LM. Overall, results from this study enhance our understanding of the impact that low level asymptomatic infection, infection in older age groups, and infection by less common malaria parasites contributes to the burden and transmission of Plasmodium species infections. In combination with high-throughput Plasmodium strain typing assays [47, 48] and novel approaches to define the levels of gametocytemia, the accurate classification of Plasmodium infections in the field and the understanding of malaria epidemiology is significantly improved., Future large scale malaria field studies and, in particular, intervention trials would therefore benefit significantly by improving capacity for molecular diagnosis of Plasmodium infections. Increased investment into establishing molecular diagnosis capacity in malaria endemic countries should be considered a high priority.
Haworth J: The global distribution of malaria and the present control effort. Malaria – Principle and Practice of Malariology. Edited by: Wernsdorfer WH, McGregor IA. 1988, Edinburgh: Churchill Livingstone, 2: 1379-1420.
Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.
Briet OJ, Gunawardena DM, Hoek van der W, Amerasinghe FP: Sri Lanka malaria maps. Malar J. 2003, 2: 22-10.1186/1475-2875-2-22.
Gemperli A, Sogoba N, Fondjo E, Mabaso M, Bagayoko M, Briet OJ, Anderegg D, Liebe J, Smith T, Vounatsou P: Mapping malaria transmission in West and Central Africa. Trop Med Int Health. 2006, 11: 1032-1046. 10.1111/j.1365-3156.2006.01640.x.
Kazembe LN, Kleinschmidt I, Holtz TH, Sharp BL: Spatial analysis and mapping of malaria risk in Malawi using point-referenced prevalence of infection data. Int J Health Geogr. 2006, 5: 41-10.1186/1476-072X-5-41.
Omumbo JA, Hay SI, Snow RW, Tatem AJ, Rogers DJ: Modelling malaria risk in East Africa at high-spatial resolution. Trop Med Int Health. 2005, 10: 557-566. 10.1111/j.1365-3156.2005.01424.x.
Cattani JA, Moir JS, Gibson FD, Ginny M, Paino J, Davidson W, Alpers MP: Small-area variations in the epidemiology of malaria in Madang Province. P N G Med J. 1986, 29 (1): 11-17.
Incardona S, Vong S, Chiv L, Lim P, Nhem S, Sem R, Khim N, Doung S, Mercereau-Puijalon O, Fandeur T: Large-scale malaria survey in Cambodia: novel insights on species distribution and risk factors. Malar J. 2007, 6: 37-10.1186/1475-2875-6-37.
Kleinschmidt I, Sharp B, Mueller I, Vounatsou P: Rise in malaria incidence rates in South Africa: a small-area spatial analysis of variation in time trends. Am J Epidemiol. 2002, 155: 257-264. 10.1093/aje/155.3.257.
Sogoba N, Vounatsou P, Doumbia S, Bagayoko M, Toure MB, Sissoko IM, Traore SF, Toure YT, Smith T: Spatial analysis of malaria transmission parameters in the rice cultivation area of Office du Niger, Mali. Am J Trop Med Hyg. 2007, 76: 1009-1015.
Maxwell CA, Chambo W, Mwaimu M, Magogo F, Carneiro IA, Curtis CF: Variation of malaria transmission and morbidity with altitude in Tanzania and with introduction of alphacypermethrin treated nets. Malar J. 2003, 2: 28-10.1186/1475-2875-2-28.
Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, Swalm C, Keating J, Regens JL, Shililu JI: Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. Am J Trop Med Hyg. 2003, 68: 734-742.
Richie TL: Interactions between malaria parasites infecting the same vertebrate host. Parasitology. 1988, 96: 607-639.
McKenzie FE, Bossert WH: Mixed-species Plasmodium infections of humans. J Parasitol. 1997, 83: 593-600. 10.2307/3284229.
McKenzie FE, Bossert WH: Multispecies Plasmodium infections of humans. J Parasitol. 1999, 85: 12-18. 10.2307/3285692.
Zimmerman PA, Mehlotra RK, Kasehagen LJ, Kazura JW: Why do we need to know more about mixed Plasmodium species infections in humans?. Trends Parasitol. 2004, 20: 440-447. 10.1016/j.pt.2004.07.004.
Snounou G, Pinheiro L, Goncalves A, Fonseca L, Dias F, Brown KN, do R, V E: Identification of the four human malaria parasite species in field samples by polymerase chain reaction and detection of a high prevalence of mixed infections. Mol Biochem Parasitol. 1993, 58: 283-292. 10.1016/0166-6851(93)90050-8.
Mehlotra RK, Kasehagen LJ, Baisor M, Lorry K, Kazura JW, Bockarie MJ, Zimmerman PA: Malaria infections are randomly distributed in diverse holoendemic areas of Papua New Guinea. Am J Trop Med Hyg. 2002, 67: 555-562.
Mehlotra RK, Lorry K, Kastens W, Miller SM, Alpers MP, Bockarie M, Kazura JW, Zimmerman PA: Random distribution of mixed species malaria infections in Papua New Guinea. Am J Trop Med Hyg. 2000, 62: 225-231.
Coleman RE, Sattabongkot J, Promstaporm S, Maneechai N, Tippayachai B, Kengluecha A, Rachapaew N, Zollner G, Miller RS, Vaughan JA: Comparison of PCR and microscopy for the detection of asymptomatic malaria in a Plasmodium falciparum/vivax endemic area in Thailand. Malar J. 2006, 5: 121-10.1186/1475-2875-5-121.
Muller I, Bockarie M, Alpers M, Smith T: The epidemiology of malaria in Papua New Guinea. Trends Parasitol. 2003, 19: 253-259. 10.1016/S1471-4922(03)00091-6.
Mueller I, Zimmerman PA, Reeder JC: Plasmodium malariae and Plasmodium ovale – the "bashful" malaria parasites. Trends Parasitol. 2007, 23: 278-283. 10.1016/j.pt.2007.04.009.
McNamara DT, Kasehagen LJ, Grimberg BT, Cole-Tobian J, Collins WE, Zimmerman PA: Diagnosing infection levels of four human malaria parasite species by a polymerase chain reaction/ligase detection reaction fluorescent microsphere-based assay. Am J Trop Med Hyg. 2006, 74: 413-421.
Kasehagen LJ, Mueller I, McNamara DT, Bockarie MJ, Kiniboro B, Rare L, Lorry K, Kastens W, Reeder JC, Kazura JW, Zimmerman PA: Changing patterns of Plasmodium blood-stage infections in the Wosera region of Papua New Guinea monitored by light microscopy and high throughput PCR diagnosis. Am J Trop Med Hyg. 2006, 75: 588-596.
Michon P, Cole-Tobian JL, Dabod E, Schoepflin S, Igu J, Susapu M, Tarongka N, Reeder JC, Beeson JG, Schofield L, King CL, Mueller I: The risk of malarial infections and disease in Papua New Guinean children. Am J Trop Med Hyg. 2007, 76: 997-1008.
Genton B, Al Yaman F, Beck HP, Hii J, Mellor S, Narara A, Gibson N, Smith T, Alpers MP: The epidemiology of malaria in the Wosera area, East Sepik Province, Papua New Guinea, in preparation for vaccine trials. I. Malariometric indices and immunity. Ann Trop Med Parasitol. 1995, 89: 359-376.
Smith T, Hii JL, Genton B, Muller I, Booth M, Gibson N, Narara A, Alpers MP: Associations of peak shifts in age – prevalence for human malarias with bednet coverage. Trans R Soc Trop Med Hyg. 2001, 95: 1-6. 10.1016/S0035-9203(01)90314-1.
McNamara DT, Thomson JM, Kasehagen LJ, Zimmerman PA: Development of a multiplex PCR-ligase detection reaction assay for diagnosis of infection by four human malaria parasite species. J Clin Microbol. 2004, 42: 2403-2410. 10.1128/JCM.42.6.2403-2410.2004.
Lunn DJ, Thomas A, Best N, Spiegelhalter D: WinBUGS – a Bayesian modelling framework: concepts, structure, and extensibility. Statistics and Computing. 2000, 10: 325-337. 10.1023/A:1008929526011.
Michon P, Cole-Tobian JL, Dabod E, Schoepflin S, Igu J, Susapu M, Tarongka N, Zimmerman PA, Reeder JC, Beeson JG, Schofield L, King CL, Mueller I: The risk of malarial infections and disease in Papua New Guinean children. Am J Trop Med Hyg. 2007, 76: 997-1008.
Bruce MC, Macheso A, Kelly-Hope LA, Nkhoma S, McConnachie A, Molyneux ME: Effect of Transmission Setting and Mixed Species Infections on Clinical Measures of Malaria in Malawi. PLoS One. 2008, 3: e2775-10.1371/journal.pone.0002775.
Bruce MC, Galinski MR, Barnwell JW, Donnelly CA, Walmsley M, Alpers MP, Walliker D, Day KP: Genetic diversity and dynamics of Plasmodium falciparum and P. vivax populations in multiply infected children with asymptomatic malaria infections in Papua New Guinea. Parasitology. 2000, 121: 257-272. 10.1017/S0031182099006356.
Hawley WA, Phillips-Howard PA, ter Kuile FO, Terlouw DJ, Vulule JM, Ombok M, Nahlen BL, Gimnig JE, Kariuki SK, Kolczak MS, Hightower AW: Community-wide effects of permethrin-treated bed nets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003, 68 (4 Suppl): 121-127.
Lengeler C: Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004, 2: CD000363-
Luo D, Lu D, Yao R, Li P, Huo X, Li A, Wen L, Ge C, Zhang S, Huo H: Alphamethrin-impregnated bed nets for malaria and mosquito control in China. Trans R Soc Trop Med Hyg. 1994, 88: 625-628. 10.1016/0035-9203(94)90199-6.
Richards FO, Klein RE, Flores RZ, Weller S, Gatica M, Zeissig R, Sexton J: Permethrin-impregnated bed nets for malaria control in northern Guatemala: epidemiologic impact and community acceptance. Am J Trop Med Hyg. 1993, 49: 410-418.
Rowland M, Bouma M, Ducornez D, Durrani N, Rozendaal J, Schapira A, Sondorp E: Pyrethroid-impregnated bed nets for personal protection against malaria for Afghan refugees. Trans R Soc Trop Med Hyg. 1996, 90: 357-361. 10.1016/S0035-9203(96)90505-2.
Habluetzel A, Cuzin N, Diallo DA, Nebie I, Belem S, Cousens SN, Esposito F: Insecticide-treated curtains reduce the prevalence and intensity of malaria infection in Burkina Faso. Trop Med Int Health. 1999, 4: 557-564. 10.1046/j.1365-3156.1999.00455.x.
Hii JL, Smith T, Vounatsou P, Alexander N, Mai A, Ibam E, Alpers MP: Area effects of bednet use in a malaria-endemic area in Papua New Guinea. Trans R Soc Trop Med Hyg. 2001, 95: 7-13. 10.1016/S0035-9203(01)90315-3.
Gimnig JE, Vulule JM, Lo TQ, Kamau L, Kolczak MS, Phillips-Howard PA, Mathenge EM, ter Kuile FO, Nahlen BL, Hightower AW: Impact of permethrin-treated bed nets on entomologic indices in an area of intense year-round malaria transmission. Am J Trop Med Hyg. 2003, 68 (4 Suppl): 16-22.
Cochrane SH, Leslie J, O'Hara DJ: Parental education and child health: intracountry evidence. Health Policy Educ. 1982, 2: 213-250. 10.1016/0165-2281(82)90011-X.
Olaogun AA, Adebayo AA, Ayandiran OE, Olasode OA: Effects of mothers' socio-economic status on the management of febrile conditions in their under five children in a resource limited setting. BMC Int Health Hum Rights. 2006, 6: 1-10.1186/1472-698X-6-1.
Noor AM, Omumbo JA, Amin AA, Zurovac D, Snow RW: Wealth, mother's education and physical access as determinants of retail sector net use in rural Kenya. Malar J. 2006, 5: 5-10.1186/1475-2875-5-5.
Bicego GT, Boerma JT: Maternal education and child survival: a comparative study of survey data from 17 countries. Soc Sci Med. 1993, 36: 1207-1227. 10.1016/0277-9536(93)90241-U.
Wamani H, Tylleskar T, Astrom AN, Tumwine JK, Peterson S: Mothers' education but not fathers' education, household assets or land ownership is the best predictor of child health inequalities in rural Uganda. Int J Equity Health. 2004, 3: 9-10.1186/1475-9276-3-9.
Benet A, Mai A, Bockarie F, Lagog M, Zimmerman PA, Alpers MP, Reeder JC, Bockarie MJ: Polymerase chain reaction diagnosis and the changing pattern of vector ecology and malaria transmission dynamics in Papua New Guinea. Am J Trop Med Hyg. 2004, 71: 277-284.
Falk N, Maire N, Sama W, Owusu-Agyei S, Smith T, Beck HP, Felger I: Comparison of PCR-RFLP and Genescan-based genotyping for analyzing infection dynamics of Plasmodium falciparum. Am J Trop Med Hyg. 2006, 74: 944-950.
Cole-Tobian J, Zimmerman PA, King CL: High throughput identification of the predominant malaria parasite clone in complex blood stage infections using a multi-SNP molecular haplotyping assay. Am J Trop Med Hyg. 2007, 76: 12-19.
Karl S, David M, Moore L, Grimberg BT, Michon P, Mueller I, Zborowski M, Zimmerman PA: Enhanced detection of gametocytes by magnetic deposition microscopy predicts higher potential for Plasmodium falciparum transmission. Malar J. 2008, 7: 66-10.1186/1475-2875-7-66.
First and foremost, we could like to thank all members of study communities that participated in the surveys. We thank Lawrence Rare and Jim Wildig for assistance with the field work and Rex Ivivi and Jonah Iga for assistance with microscopy. Finally, we thank Dr. Charles King for critical comments leading to completion of this manuscript. This work as supported with funding from National Institutes of Health (AI063135, AI-46919, and TW007872) and AusAID.
The authors declare that they have no competing interests.
IM designed the study, analysed the data and wrote the first draft. SW, DM, SM, BK AS conducted the field work and read blood slides. SW and DM also assisted in data analyses. DTM and PAZ conducted LDR-FMA typing of all samples. TAS assisted with the spatial data analyses. All authors contributed to the interpretation of the data and the writing of manuscript.
Electronic supplementary material
Additional file 2: Table S2. Comparison of observed species mixtures in different regions by light microscopy. (PDF 61 KB)
About this article
Cite this article
Mueller, I., Widmer, S., Michel, D. et al. High sensitivity detection of Plasmodium species reveals positive correlations between infections of different species, shifts in age distribution and reduced local variation in Papua New Guinea. Malar J 8, 41 (2009). https://doi.org/10.1186/1475-2875-8-41
- Mixed Infection
- Parasite Density
- Plasmodium Species