Non-falciparum malaria infections are as prevalent as P. falciparum among Tanzanian schoolchildren

Rachel Sendor, Cedar L. Mitchell, Frank Chacky, Ally Mohamed, Lwidiko E. Mhamilawa, Fabrizio Molteni, Ssanyu Nyinondi, Bilali Kabula, Humphrey Mkali, Erik J. Reaves, Naomi Serbantez, Chonge Kitojo, Twilumba Makene, Thwai Kyaw, Meredith Muller, Alexis Mwanza, Erin L. Eckert, Jonathan B. Parr, Jessica T. Lin, Jonathan J. Juliano, Billy Ngasala
doi: https://doi.org/10.1101/2022.06.07.22275625

Abstract

Efforts to achieve malaria elimination need to consider both falciparum and non-falciparum infections. The prevalence and geographic distribution of four Plasmodium species were determined by real-time PCR using dried blood spots collected during the 2017 School Malaria Parasitological Survey of eight regions of Tanzania. Among 3,456 schoolchildren, 22% had P. falciparum, 24% P. ovale spp., 4% P. malariae, and 0.3% P. vivax. Ninety-one percent of P. ovale infections had very low parasite densities, based on amplification at later cycle thresholds. Sixty-four percent of P. ovale infections were single-species, and 35% of these were detected in low malaria endemicity regions. P. malariae infections were predominantly co-infections with P. falciparum (73%). P. vivax was largely detected in northern and eastern regions. Overall, 43% of children with P. falciparum were co-infected with at least one non-falciparum species. A large, previously under-appreciated burden of P. ovale spp. infection exists among Tanzanian schoolchildren.

Article Summary A previously unrecognized burden of non-falciparum malaria infections was detected among Tanzanian schoolchildren in a 2017 cross-sectional study, with P. ovale spp. prevalence comparable to P. falciparum, and low-level prevalences of P. malariae and P. vivax detected.

Background

Sub-Saharan Africa harbors 95% of the global malaria burden1. While national surveys conducted by ministries of health throughout Africa regularly assess Plasmodium falciparum burden2, little is known about the prevalence and geographic distribution of non-falciparum malaria species, including Plasmodium malariae, Plasmodium vivax, and Plasmodium ovale spp.3–8. Although the clinical burden of non-falciparum malaria in sub-Saharan Africa is dwarfed by P. falciparum9, these species can still elicit symptomatic disease. P. malariae has been associated with an increased risk of anemia10, as well as clinical complications such as chronic nephrotic syndrome11,12. P. vivax causes severe anemia, pregnancy-related complications, and mortality following recurrent infections, although infections in sub-Saharan Africa are infrequent13–15. Clinical consequences of P. ovale spp. infection are least well-described, though well-described in travelers, and have been associated with severe infection in case reports16.

Recent findings suggest declining P. falciparum prevalence in East Africa may be associated with increases in non-falciparum infections17–20. However, comprehensive surveys of non-falciparum malaria in sub-Saharan Africa are infrequent because detection of these species remains challenging11,17. Existing field diagnostic methods such as microscopy and pan-Plasmodium spp. lactate dehydrogenase (LDH) or histidine-rich protein 2 (HRP2)-based rapid diagnostic tests (RDTs) lack sensitivity to detect non-falciparum species11,17. Non-falciparum parasite densities are often low, and the majority of those infected may not seek care. Mixed-infections with P. falciparum can also complicate detection of concurrent infection by non-falciparum species3,17. Molecular detection methods can sensitively detect non-falciparum species, but remain largely confined to research use.

In Tanzania, malaria burden is high, with the country accounting for 4.1% of malaria deaths globally in 20201. While approximately 93% of the mainland Tanzanian population remains at risk of infection, transmission throughout the country is highly heterogeneous21. Malaria transmission patterns are largely driven by geographic features of the country: across the arid highlands and in urban centers, malaria transmission is low, unstable, and seasonal; in southern, northern, and northwestern areas, transmission is moderate and varies seasonally; and along the coastal, lake, and southern regions, malaria transmission is high and perennial21,22.

Decades of concentrated malaria control interventions have helped lower the national prevalence from 18% in 2008 to 7% in 201723. The vast majority of malaria cases reported in Tanzania are attributed to P. falciparum9,21, although recent studies have also identified evidence of P. malariae, P. vivax, and P. ovale transmission4,18,24,25. Given widespread use of P. falciparum-specific HRP2-based RDTs for malaria diagnosis, the propensity for missed detection or misclassification of non-falciparum species in Tanzania is high, and large-scale, geographically representative studies to assess spatial distributions of these species are lacking.

This study deployed molecular methods in the context of a national survey to comprehensively characterize non-falciparum malaria epidemiology among schoolchildren in Tanzania. We demonstrate a high prevalence of P. ovale spp. infection, on par with P. falciparum prevalence. These infections vary geographically and do not always correlate with P. falciparum prevalence in a given region, with many P. ovale spp. infections detected in regions considered at low epidemiological risk for malaria. P. malariae was also broadly detected, whereas P. vivax was detected in limited geographical regions.

Methods

Study Design

The 2017 School Malaria Parasitological Survey (SMPS) was a cross-sectional study conducted in 2017 among children aged 5–16 years enrolled in public primary schools in mainland Tanzania. Methods for site selection and survey design mirrored the 2015 SMPS and have been previously summarized in detail22. Regions included in the study were selected through a multi-stage sampling scheme to maintain geographic representation and reflect the heterogeneity of malaria transmission across Tanzania22,26. The number of schools randomly selected per region was proportional to each region’s respective population22,26. Within each school, an average of 100 students were randomly selected for screening, and consented students were interviewed for demographic and clinical characteristics, had a malaria RDT performed, and provided a dried blood spot (DBS) sample22,26. Surveys largely coincided with each region’s rainy season. Among students who provided a DBS, a stratified random sample was selected to construct a sub-population for non-falciparum malaria testing, with students sampled in equal proportion to the regional proportion of students within the broader survey population to maintain representativeness.

Malaria detection was conducted using CareStart™ Malaria Pf/PAN (HRP2/pLDH) Ag Combo RDTs (Product code: RMRM 02571; Access Bio, Inc., Somerset, NJ, USA) detecting P. falciparum-specific HRP2 and pan-pLDH antigens; RDTs were considered positive if either band was positive. Schools and councils were grouped into epidemiological strata of malaria transmission risk based on P. falciparum prevalences in children estimated from the 2014-15 Tanzania SMPS22,26. P. falciparum prevalence was considered “very low” if <5%, “low” 5 to <10%, “moderate” 10 to <50%, and “high” ≥50%22,26. Dried blood spots (DBS) were collected on Whatman filter paper (Whatman plc, Cytiva, Buckinghamshire, UK) and shipped to the University of North Carolina, Chapel Hill for molecular testing.

Molecular Detection

DNA was extracted from three 6mm DBS punches using Chelex methods27, and samples underwent real-time PCR targeting the 18S ribosomal subunit for malaria testing, according to previously published protocols28 (Supplementary Table 1). Assays for each Plasmodium species (falciparum, ovale spp., malariae, and vivax) were run independently with appropriate controls. Positive controls for P. falciparum detection were made using whole blood and 3D7 cultured parasites (MRA-102; BEI Resources) to create mocked DBS, with DNA extracted using Chelex27 and serially diluted. Positive controls for non-falciparum species were based on standard dilutions of plasmid DNA from MR4 (MRA-180, MRA-179, MRA-178; BEI Resources), with semi-quantitative parasitemia estimated using six 18S rRNA gene copies per parasite28. P. malariae and P. falciparum assays were run for 40 cycles. P. ovale spp. and P. vivax assays were run for 45 cycles to enable detection of very-low density infections28. This approach was previously validated using 390 negative controls (comprising water [n=22] and human DNA [n=368]) and >170 positive controls of decreasing parasite densities per non-falciparum species; no false-positives were detected28. Assay specificity was also assessed previously by testing against other Plasmodium species29. In addition, the UNC laboratory participates in the World Health Organization malaria molecular quality assurance scheme, identifying and determining species of Plasmodium parasites in blinded samples every six months, with consistent high marks for assay performance across species. In the present study, among 20 negative controls per each species-specific assay, no false positive amplification was detected (Supplementary Table 2). A subset of positive P. ovale samples underwent further real-time PCR testing to distinguish wallikeri and curtisi species30,31. To evaluate potential bias due to differences in PCR cycles between species, a sensitivity analysis was conducted in a random sample of 750 students stratified by malaria transmission intensity, detecting P. falciparum and P. malariae infections by semi-quantitative, real-time PCR of the 18S RNA gene performed to 45 cycles.

Analysis

Malaria species-specific prevalences were calculated overall, and by single vs. mixed-species infections. Prevalences were not adjusted with sampling weights as samples for the non-falciparum analysis population were selected randomly and in equal proportion to the broader survey sample.

Descriptive statistics were computed for student-level characteristics by Plasmodium species. Differences between P. falciparum and non-falciparum species mono-infections were tested using Pearson’s chi-square analysis and Kruskal-Wallis rank sum testing in the presence of non-normality; Fisher’s exact test was applied where cell counts were <5. Similar analyses were performed comparing PCR–positive vs. PCR–negative students by species. Missing data were excluded from statistical analysis.

Spatial Mapping

Regional variation in the prevalence of each species was assessed through geospatial mapping of infections by council and region. Number of infections and number of students tested were aggregated by council, and prevalence estimates by council were estimated and mapped.

Scaled prevalence estimates were calculated by taking the proportion of each council’s prevalence out of the highest council prevalence for each species, as follows: Embedded Image where P is the prevalence for a given council, di.

Differences between scaled non-falciparum prevalence and scaled P. falciparum prevalence for each council were calculated and mapped to determine differential scaled prevalences. This method allowed comparison of prevalence estimates between each non-falciparum species with P. falciparum, while accounting for differences in the absolute burden of each species.

Analyses and mapping were performed using R (version 4.0.2), including eulerr32 and sf (version 0.9-7)33 packages. Shapefiles were sourced from GADM.org and elevation measurements were collected from the Shuttle Radar Topography Mission. Informed consent was obtained from students and their legal guardians prior to survey or blood sample collection and ethical clearance was obtained from the Tanzania National Institute for Medical Research. Analysis of de-identified samples was approved by the Institutional Review Board of the University of North Carolina, Chapel Hill (IRB #19-1495).

Results

Study Population

A total of 3,456 students from 180 schools across eight geographically-representative regions were selected for non-falciparum testing from among 17,131 students in the SMPS who had a DBS sample available. No significant differences were detected between the sub-population for non-falciparum testing and the SMPS DBS population (Supplementary Table 3). Median (IQR) student age in the non-falciparum analysis population was 11 (9–13) years, with an approximately equal distribution of males and females (51% vs. 49%). Malaria dual-antigen RDTs were positive in 20% of students. The majority of students attended schools in regions classified as ‘high’ (51%) or ‘moderate’ (13%) malaria transmission risk (Table 1).

View this table:
Table 1.

Student infection characteristics by malaria species*

Species Prevalence by PCR

P. falciparum infection was identified in 22% (n=755; 95% CI: 21%–23%) of students, P. ovale spp. in 24% (n=814; 95% CI: 22%–25%), P. malariae in 4% (n=136; 95% CI: 3%–5%), and P. vivax in 0.3% (n=11; 95% CI: 0.2%–0.6%), including single and mixed-species infections (Supplemental Table 4).

The majority of P. ovale spp. infections identified were single-species (64%, n=519), with 28% (n=224) mixed with P. falciparum only. Conversely, most P. malariae infections were co-infected with both P. ovale spp. and P. falciparum (40%, n=55), or with P. falciparum only (32%, n=44) [Figure 1]. Thirty-six percent (n=4) of P. vivax infections were single-species. Forty-three percent (n=326) of P. falciparum infections were co-infected with at least one non-falciparum species.

Figure 1.
Figure 1.

Distribution of malaria infection by species among Tanzanian schoolchildren

Note: Pf: P. falciparum; Pm: P. malariae; Po: P. ovale spp.; Pv: P. vivax. N=1 Pv + Pm co-infection is not depicted. Prevalence estimates by species are: P. falciparum: 21.8% (n=755; 95% CI: 20.5% – 23.3%); P. ovale spp.: 23.6% (n=814; 95% CI: 22.2% – 25.0%); P. malariae: 3.9% (n=136; 95% CI: 3.3% – 4.6%); P. vivax: 0.3% (n=11; 95% CI: 0.2% – 0.6%).

To evaluate the impact of different PCR cycles between assays, a sensitivity analysis testing P. falciparum and P. malariae within 45 PCR cycles was performed. A 25% (95% CI: 21%–29%) prevalence of P. falciparum and 3% (95% CI: 2%–5%) prevalence of P. malariae was observed, weighted to the distribution of students in the total non-falciparum population by transmission intensity strata. Within this subset, 0.5% (n=4) of P. falciparum and 0.3% (n=2) of P. malariae infections were detected at Cts between 40 to 45 (Supplemental Table 5). Thus, >99% of P. falciparum and P. malariae infections were detected by the primary 40-cycle assays utilized in the study.

Differences in subject-level characteristics between species are depicted in Table 1 (and Supplemental Table 6). Compared to P. falciparum, P. ovale spp. mono-infections were detected more frequently in slightly younger (median 11 vs. 12 years; p<0.001) and female (54% vs. 45%; p=0.009) students. In terms of sensitivity by RDT compared to PCR, 8% (n=40) of P. ovale spp. mono-infections were RDT-positive for any band, whereas 33% (n=8) of P. malariae and 69% (n=295) of P. falciparum mono-infections were RDT-positive. Non-falciparum-P. falciparum co-infections were positive by RDT in 78% (n=253/325) of infections detected by PCR.

Though only 3% (n=13) of P. falciparum mono-infections and no P. malariae or P. vivax mono-infections were detected in low-risk strata, 35% (n=181) of P. ovale spp. mono-infections occurred in regions considered to have ‘low’ or ‘very low’ epidemiological risk for malaria. ‘High” epidemiologic risk strata harbored the majority of mono-infections across all four malaria species, as well as mixed-infections with P. falciparum.

Parasite Density

Malaria parasitemia estimated by semi-quantitative PCR was low across non-falciparum species (Figure 2). The median (IQR; min–max) parasite densities for P. ovale spp. and P. malariae were comparable (1.8 [0.3 – 6.2; 0.02 – 42,149]) ovale parasites (p)/μL and 2.9 [0.7 – 13.7; 0.1 – 304] malariae p/μL, respectively). P. vivax parasite densities were an order of magnitude lower, around 0.1 (0.1 – 0.2; 0.02 – 2.0) P. vivax p/μL. Whereas a small proportion of P. malariae infections had a parasite density >50 p/μL (12%, n=16), this was rarely observed for P. ovale spp. (1%, n=11), and never the case for P. vivax, where the highest parasite density detected was 2.0 p/μL. Median (IQR; min-max) P. falciparum parasite density was also low at 13.1 (2.6 – 55.9; 0.1 – 8,248) p/μL, although 27% (n=203) of P. falciparum cases had a parasite density >50 p/μL, and 3% (n=24) >500 p/μL.

Figure 2.
Figure 2.

Estimated parasite density (p/μL) distributions by malaria species

Note: Parasite densities displayed on log scale. Mixed infections include P. falciparum and non-falciparum co-infections. P. ovale and P. vivax parasite densities were detected by including isolates amplifying at later PCR cycle thresholds. A) P. falciparum median (IQR) parasite density: (single-species [n=429]: 11.4 [2.5-54.7] p/μL; mixed [n=326]: 16.5 [3.5-56.9] p/μL; p=0.117). B) P. ovale spp. median (IQR) parasite density: (single-species [n=519]: 3.4 [0.3-7.5] p/μL; mixed [n=295]: 0.8 [0.3-2.8] p/μL; p<0.001). C) P. malariae median (IQR) parasite density: (single-species [n=24]: 4.0 [0.9-40.9] p/μL; mixed [n=112]: 2.8 [0.6-13.5] p/μL; p=0.169). D) P. vivax median (IQR) parasite density: (single-species [n=4]: 0.1 [0.1-0.2] p/μL; mixed [n=7] 0.2 [0.1-0.2] p/μL; p=0.571).

Median (IQR) parasite density among P. ovale spp. mixed-infections was lower than among P. ovale spp. mono-infections (0.8 [0.3-2.8] ovale p/μL vs. 3.4 [0.3-7.5] p/μL; p<0.001), whereas densities were similar between single and mixed-species infections among all other malaria species (Figure 2).

Plasmodium ovale Species Determination

Among 814 samples positive for P. ovale spp. infection, 60 (7%) of the highest parasite density samples were selected for real-time PCR detection of P. ovale wallikeri and P. ovale curtisi species. Species determination was successful in 35% (n=21) of samples, with P. ovale curtisi detected in 17 samples and P. ovale wallikeri in 9. Five students were positive for both P. ovale curtisi and P. ovale wallikeri infection. Further characterization by demographic and clinical characteristics was not performed due to limited sample size.

Geographic Distribution

P. ovale spp. were detected across all eight regions sampled in the study, indicating widespread distribution across Tanzania (Table 2, Figure 3). Prevalence of P. ovale spp. was highest within the northern Kagera region (34%, n=273) and the central region of Tabora (17%, n=139). P. ovale curtisi infections were detected in six of eight study regions (all but Arusha and Rukwa), and P. ovale wallikeri was detected in five (Kagera, Mara, Tabora, Tanga, and Iringa). High prevalence of P. malariae was also detected in Kagera (29%, n=39), similar to P. ovale spp., and in the southernmost Mtwara region (28%, n=38). P. vivax infections were predominantly distributed along the northern-western border of Tanzania in Kagera (55%, n=6), with select, isolated cases of single and mixed-species P. vivax infections also detected in southern and eastern regions. The Arusha region in the northeast and Iringa in the southwest did not have any cases of P. malariae or P. vivax infections, and had the lowest frequencies of P. ovale spp. (3%, n=23; 4%, n=30) and P. falciparum infections (0.4%, n=3; 0.1%, n=1).

View this table:
Table 2.

Student infections by Plasmodium species and school characteristics across Tanzania

Figure 3.
Figure 3. School council-level malaria prevalence distributions by species

Malaria infections across all species were predominantly detected at elevations <1500 meters (m), including 85% (n=693) of P. ovale spp., 94% (n=707) of P. falciparum, 95% (n=129) of P. malariae, and 100% (n=11) of P. vivax infections (Table 2). In comparison, 77% (n=2667) of students enrolled in the study were derived from schools at elevations <1500m.

Among those enrolled at higher elevations, P. ovale spp. infections were detected most frequently at 15% (n=121), compared to 6% (n=48) of students with P. falciparum, 1% (n=7) P. malariae, and 0 P. vivax.

Comparing scaled prevalence estimates for non-falciparum species with P. falciparum, we identified areas where prevalences were higher than expected for P. ovale spp. and P. malariae, based on P. falciparum burden (Figure 4); P. vivax infections were too infrequent for comparison. In the southern and south-western highlands and north-western lake regions (Iringa, Rukwa, Tabora, and Kagera) scaled P. ovale spp. prevalences were higher relative to P. falciparum. Scaled prevalence of P. malariae was notably higher than P. falciparum in the Karagwe council in Kagera and the Mtwara municipal council in Mtwara. In most other areas, scaled prevalence of P. malariae was similar or lower to P. falciparum prevalence.

Figure 4.
Figure 4.

Differential scaled prevalences of P. malariae and P. ovale spp. compared to P. falciparum at the school council-level

Note: scaled prevalence of P.f–P.v. is not depicted given the low number of P. vivax infections biasing the scaled measure. A) Blue shading indicates councils where P. falciparum scaled prevalence is greater than P. malariae scaled prevalence in that council. Purple indicates regions where P. malariae scaled prevalence is greater. B) Blue shading indicates councils where P. falciparum scaled prevalence is greater than P. ovale spp. scaled prevalence in that council. Green indicates regions where P. ovale spp. scaled prevalence is greater.

Discussion

This study represents the largest nationally representative molecular survey of non-falciparum malaria epidemiology across Tanzania to our knowledge. Real-time PCR was employed to estimate non-falciparum infection prevalences within school-aged children in eight regions of the country selected to maintain geographic diversity and heterogeneity of malaria transmission risk. One-quarter of schoolchildren (24%) harbored P. ovale spp. parasites, comparable to the detected P. falciparum prevalence (22%) in the population, with more than half (64%) constituting single-species infections. P. malariae was observed in 4% of students, with the majority co-infected with other malaria species. P. vivax infections were rare but present at a 0.3% prevalence.

The high prevalence of P. ovale spp. may be attributed to several factors. First, we extended the cycling time of our P. ovale PCR assay to 45 cycles to enable detection of low-density infections, which comprised 91% of all P. ovale spp. infections identified. There is precedent for this approach17,25,31, as low density parasitemia is characteristic of P. ovale spp. infections, making them challenging to detect. The use of 40-cycle P. ovale PCR assays yielded a lower 0.8% prevalence estimate in our prior work in the Democratic Republic of Congo (DRC)5. The prevalence of P. ovale spp. infections positive at <40 cycles in the present study was 9% (n=75), confirming the majority of infections occurred at very-low parasite densities. Second, many large-scale molecular surveys of non-falciparum malaria have focused on adults or all-age cohorts, whereas school-aged children are increasingly recognized as the main contributors to asymptomatic and infectious malaria reservoirs34–36. Finally, the high prevalence of P. ovale spp. in our study may reflect increasing P. ovale spp. transmission in the face of malaria control efforts targeting P. falciparum. Increasing or persistent transmission of P. ovale and P. malariae amid P. falciparum decline has been observed in molecular surveys from Tanzania and nearby Kenya and Uganda, including in symptomatic cases17,18,24,37. It’s unclear whether this is due to hypnozoite-induced relapses from P. ovale spp. infections not treated by artemisinin-based combination therapies (ACTs), day-biting or outdoor vectors that evade bed nets, or other factors.

In contrast to findings from other studies11,37–40, P. ovale spp. more commonly occurred as single-species infections compared to the other non-falciparum species, though this finding may be due in part to the increased sensitivity of the P. ovale PCR assay employed. These single-species infections were rarely detected by RDTs, rendering them more difficult to detect and treat. Single-species P. ovale infections were also largely the only species identified within regions categorized as low risk for malaria transmission, suggesting an unexpected transmission risk in areas where prevention measures may be less frequently used and P. falciparum risk is not a large concern. Our scaled differential prevalence map similarly highlights several councils where P. ovale spp. and P. malariae prevalence are proportionately higher than would be expected based on P. falciparum burden. Taken together, these characteristics indicate a hidden burden of P. ovale spp. in the country.

Detection of P. vivax in this study is notable given the control challenges this species poses. Infections were predominately detected in the northwest/Lake zone regions of Tanzania and in the east, where several studies have also observed low-level P. vivax prevalences4,24,41. The P. malariae prevalence of 4% detected in this study aligns with recent research in the region identifying a low prevalence of infection (2.5% in Malawi, 4.1% in DRC, 3.3% and 5.3% [symptomatic and asymptomatic] in western Kenya)12,28,40. Estimated parasite densities were low across all non-falciparum species, as expected. P. falciparum parasite densities were relatively low as well (median 13.1 p/μL), likely due to the survey of a predominantly asymptomatic population. Additionally, prevalence mapping confirmed low or non-existent prevalence of non-falciparum malaria within the northern highlands of Arusha and southern highlands/midlands of Iringa.

This study has several limitations. The use of different PCR cycling times by species introduces ascertainment bias. Because P. malariae and P. falciparum assays were run at 40 rather than 45 cycles, their relative prevalence compared to P. ovale spp. may be underestimated. However, a sensitivity analysis to quantify this bias indicated we would expect to detect only an additional 0.5% of P. falciparum and 0.3% of P. malariae infections if testing to 45 cycles, suggesting minimal underestimation of reported P. falciparum and P. malariae prevalences, and no meaningful impact on overall conclusions. Weighting sensitivity analysis results to the total study population yielded a P. falciparum prevalence of 25% if run to 45 cycles, compared to the observed prevalence of 22%. Additionally, this study did not sample all geographic regions in Tanzania, and findings cannot be extrapolated to other age groups with differing malaria risk profiles. School-based sampling likely underestimates prevalence of infection from symptomatic or severe malaria infections in school-aged children, as children may be absent due to illness. Lastly, the cross-sectional survey design reveals little about clinical implications of prevalent non-falciparum infections, especially given substantial, non-random missingness in fever data, or to what extent these infections represent chronic carriage versus transient parasitemia.

The overall high prevalence and broad geographic distribution of P. ovale spp. and to a lesser extent, P. malariae, with more focal distribution of P. vivax, detected in this study underscore an urgent need to understand clinical burden and transmission patterns of these species to inform malaria control programs in Tanzania. Current Tanzanian treatment protocols do not regularly treat for hypnozoite liver-stage P. ovale spp. infection, and relapses are expected after blood-stage clearance by ACTs42. With accumulating evidence of a previously unappreciated non-falciparum malaria burden present in sub-Saharan Africa39, new detection and treatment strategies may be required for continued progress toward malaria control and elimination.

Data Availability

Data are available upon reasonable request to the authors.

Funding

This study was in part funded by the National Institute of Health (K24AI134990 and R01TW010870 to JJJ; T32AI070114 to CLM; R21AI152260 to JTL; R21 AI148579 to JTL and JBP; R01 AI139520 to JBP, RS, and CLM). This work was also made possible through the Global Fund support to conducting the survey and the President’s Malaria Initiative via the United States Agency for International Development (USAID) Okoa Maisha Dhibiti Malaria (OMDM) (Cooperative Agreement Number: 72062118CA-00002) implemented by RTI, under the terms of an inter-agency agreement with Centers for Disease Control and Prevention (CDC) on data management and facilitation of the initial processing and exporting of the blood samples. Funding sources had no role in the study design, analysis, or writing of the manuscript.

Disclaimer

The findings and conclusions in this report are those of the author(s) and do not necessarily represent the official position of the U.S. Centers for Disease Control and Prevention, the President’s Malaria Initiative via the U.S. Agency for International Development, or other employing organizations or sources of funding.

Data Availability

Data are available upon reasonable request to authors.

Acknowledgments

We thank the 2017 SMPS study administrators and staff for their tireless work implementing the survey and we are grateful for the participation of all individuals in the study. The following reagents were obtained through BEI Resources, NIAID, NIH: Diagnostic Plasmid Containing the Small Subunit Ribosomal RNA Gene (18S) from Plasmodium malariae, MRA-179, Plasmodium ovale, MRA-180, and Plasmodium vivax, MRA-178, contributed by Peter A. Zimmerman. The following reagent was obtained through BEI Resources, NIAID, NIH: Plasmodium falciparum, Strain 3D7, MRA-102, contributed by Daniel J. Carucci.

Footnotes

  • ↵‡ Co-first authors

  • ↵* Co-senior authors

References

  1. 1.↵
    World Health Organization. World malaria report 2021. Accessed December 13, 2021. https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2021
  2. 2.↵
    Alegana VA, Macharia PM, Muchiri S, et al. Plasmodium falciparum parasite prevalence in East Africa: Updating data for malaria stratification. PLOS Global Public Health. 2021;1(12):e0000014. doi:10.1371/journal.pgph.0000014
    OpenUrlCrossRef
  3. 3.↵
    Oriero EC, Amenga-Etego L, Ishengoma DS, Amambua-Ngwa A. Plasmodium malariae, current knowledge and future research opportunities on a neglected malaria parasite species. Critical Reviews in Microbiology. 2021;47(1):44–56. doi:10.1080/1040841X.2020.1838440
    OpenUrlCrossRef
  4. 4.↵
    Kim MJ, Jung BK, Chai JY, et al. High malaria prevalence among schoolchildren on Kome Island, Tanzania. Korean J Parasitol. 2015;53(5):571–574. doi:10.3347/kjp.2015.53.5.571
    OpenUrlCrossRef
  5. 5.↵
    Mitchell CL, Brazeau NF, Keeler C, et al. Under the radar: epidemiology of Plasmodium ovale in the Democratic Republic of the Congo. The Journal of Infectious Diseases. 2021;223(6):1005–1014. doi:10.1093/infdis/jiaa478
    OpenUrlCrossRef
  6. 6.
    Björkman A, Shakely D, Ali AS, et al. From high to low malaria transmission in Zanzibar— challenges and opportunities to achieve elimination. BMC Med. 2019;17:14. doi:10.1186/s12916-018-1243-z
    OpenUrlCrossRef
  7. 7.
    Bousema JT, Drakeley CJ, Mens PF, et al. Increased Plasmodium falciparum gametocyte production in mixed infections with P. malariae. Am J Trop Med Hyg. 2008;78(3):442–448.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    Sutherland CJ. Persistent Parasitism: The adaptive biology of malariae and ovale malaria. Trends Parasitol. 2016;32(10):808–819. doi:10.1016/j.pt.2016.07.001
    OpenUrlCrossRef
  9. 9.↵
    World Health Organization. World malaria report, 2020. Accessed August 29, 2021. https://www.who.int/publications-detail-redirect/9789240015791
  10. 10.↵
    Langford S, Douglas NM, Lampah DA, et al. Plasmodium malariae infection associated with a high burden of anemia: a hospital-based surveillance study. PLoS Negl Trop Dis. 2015;9(12):e0004195. doi:10.1371/journal.pntd.0004195
    OpenUrlCrossRef
  11. 11.↵
    Mueller I, Zimmerman PA, Reeder JC. Plasmodium malariae and Plasmodium ovale--the “bashful” malaria parasites. Trends Parasitol. 2007;23(6):278–283. doi:10.1016/j.pt.2007.04.009
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    Lo E, Nguyen K, Nguyen J, et al. Plasmodium malariae prevalence and csp gene diversity, Kenya, 2014 and 2015. Emerg Infect Dis. 2017;23(4):601–610. doi:10.3201/eid2304.161245
    OpenUrlCrossRef
  13. 13.↵
    Twohig KA, Pfeffer DA, Baird JK, et al. Growing evidence of Plasmodium vivax across malaria-endemic Africa. PLoS Negl Trop Dis. 2019;13(1):e0007140. doi:10.1371/journal.pntd.0007140
    OpenUrlCrossRefPubMed
  14. 14.
    Howes RE, Battle KE, Mendis KN, et al. Global epidemiology of Plasmodium vivax. Am J Trop Med Hyg. 2016;95(6 Suppl):15–34. doi:10.4269/ajtmh.16-0141
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    Price RN, Commons RJ, Battle KE, Thriemer K, Mendis K. Plasmodium vivax in the era of the shrinking P. falciparum map. Trends Parasitol. 2020;36(6):560–570. doi:10.1016/j.pt.2020.03.009
    OpenUrlCrossRefPubMed
  16. 16.↵
    Groger M, Fischer HS, Veletzky L, Lalremruata A, Ramharter M. A systematic review of the clinical presentation, treatment and relapse characteristics of human Plasmodium ovale malaria. Malar J. 2017;16(1):112. doi:10.1186/s12936-017-1759-2
    OpenUrlCrossRef
  17. 17.↵
    Akala HM, Watson OJ, Mitei KK, et al. Plasmodium interspecies interactions during a period of increasing prevalence of Plasmodium ovale in symptomatic individuals seeking treatment: an observational study. The Lancet Microbe. 2021;2(4):e141–e150. doi:10.1016/S2666-5247(21)00009-4
    OpenUrlCrossRef
  18. 18.↵
    Yman V, Wandell G, Mutemi DD, et al. Persistent transmission of Plasmodium malariae and Plasmodium ovale species in an area of declining Plasmodium falciparum transmission in eastern Tanzania. PLoS Negl Trop Dis. 2019;13(5):e0007414. doi:10.1371/journal.pntd.0007414
    OpenUrlCrossRef
  19. 19.
    Doctor SM, Liu Y, Anderson OG, et al. Low prevalence of Plasmodium malariae and Plasmodium ovale mono-infections among children in the Democratic Republic of the Congo: a population-based, cross-sectional study. Malaria Journal. 2016;15:350. doi:10.1186/s12936-016-1409-0
    OpenUrlCrossRef
  20. 20.↵
    Taylor SM, Messina JP, Hand CC, et al. Molecular malaria epidemiology: mapping and burden estimates for the Democratic Republic of the Congo, 2007. PLOS ONE. 2011;6(1):e16420. doi:10.1371/journal.pone.0016420
    OpenUrlCrossRefPubMed
  21. 21.↵
    President’s Malaria Initiative. Tanzania Malaria Operational Plan FY 2015. https://d1u4sg1s9ptc4z.cloudfront.net/uploads/2021/03/fy-2015-tanzania-malaria-operational-plan.pdf
  22. 22.↵
    Chacky F, Runge M, Rumisha SF, et al. Nationwide school malaria parasitaemia survey in public primary schools, the United Republic of Tanzania. Malaria Journal. 2018;17(1):452. doi:10.1186/s12936-018-2601-1
    OpenUrlCrossRef
  23. 23.↵
    Ministry of Health, Community Development, Gender, Elderly and Children (MoHCDGEC) [Tanzania, Mainland], Ministry of Health (MoH) [Zanzibar], National Bureau of Statistics (NBS), Office of the Chief, Government Statistician (OCGS), and ICF. Tanzania Malaria Indicator Survey 2017. Dar es Salaam, Tanzania and Rockville, Maryland, USA. https://dhsprogram.com/pubs/pdf/MIS31/MIS31.pdf.
  24. 24.↵
    Cook J, Xu W, Msellem M, et al. Mass screening and treatment on the basis of results of a Plasmodium falciparum-specific rapid diagnostic test did not reduce malaria incidence in Zanzibar. The Journal of Infectious Diseases. 2015;211(9):1476–1483. doi:10.1093/infdis/jiu655
    OpenUrlCrossRefPubMed
  25. 25.↵
    Tarimo BB, Nyasembe VO, Ngasala B, et al. Seasonality and transmissibility of Plasmodium ovale in Bagamoyo District, Tanzania. Parasit Vectors. 2022;15:56. doi:10.1186/s13071-022-05181-2
    OpenUrlCrossRef
  26. 26.↵
    Mitchell CL, Ngasala B, Janko MM, et al. Evaluating malaria prevalence and land cover across varying transmission intensity in Tanzania using a cross-sectional survey of school-aged children. Malar J. 2022;21(1):80. doi:10.1186/s12936-022-04107-8
    OpenUrlCrossRef
  27. 27.↵
    Teyssier NB, Chen A, Duarte EM, Sit R, Greenhouse B, Tessema SK. Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples. Malar J. 2021;20:116. doi:10.1186/s12936-021-03630-4
    OpenUrlCrossRef
  28. 28.↵
    Gumbo A, Topazian HM, Mwanza A, et al. Occurrence and distribution of nonfalciparum malaria parasite species among adolescents and adults in Malawi. The Journal of Infectious Diseases. 2021;225(2):257–268. doi:10.1093/infdis/jiab353
    OpenUrlCrossRef
  29. 29.↵
    Brazeau NF, Mitchell CL, Morgan AP, et al. The epidemiology of Plasmodium vivax among adults in the Democratic Republic of the Congo. Nat Commun. 2021;12(1):4169. doi:10.1038/s41467-021-24216-3
    OpenUrlCrossRef
  30. 30.↵
    Perandin F, Manca N, Calderaro A, et al. Development of a real-time PCR assay for detection of Plasmodium falciparum, Plasmodium vivax, and Plasmodium ovale for routine clinical diagnosis. J Clin Microbiol. 2004;42(3):1214–1219. doi:10.1128/JCM.42.3.1214-1219.2004
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Mueller I
    Calderaro A, Piccolo G, Gorrini C, et al. A new real-time PCR for the detection of Plasmodium ovale wallikeri. Mueller I, ed. PLoS ONE. 2012;7(10):e48033. doi:10.1371/journal.pone.0048033
    OpenUrlCrossRefPubMed
  32. 32.↵
    Larsson J. Eulerr: area-proportional euler and venn diagrams with ellipses. R package, version 6.1.0; 2020. https://cran.r-project.org/package=eulerr.
  33. 33.↵
    Pebesma E. Simple features for R: standardized support for spatial vector data. The R Journal. 2018;10(1):439–446. doi:10.32614/RJ-2018-009
    OpenUrlCrossRef
  34. 34.↵
    Abdulraheem MA, Ernest M, Ugwuanyi I, et al. High prevalence of Plasmodium malariae and Plasmodium ovale in co-infections with Plasmodium falciparum in asymptomatic malaria parasite carriers in southwestern Nigeria. International Journal for Parasitology. 2021;52(1):23–33. doi:10.1016/j.ijpara.2021.06.003
    OpenUrlCrossRef
  35. 35.
    Andolina C, Rek JC, Briggs J, et al. Sources of persistent malaria transmission in a setting with effective malaria control in eastern Uganda: a longitudinal, observational cohort study. Lancet Infect Dis. 2021;21(11):1568–1578. doi:10.1016/S1473-3099(21)00072-4
    OpenUrlCrossRef
  36. 36.↵
    Walldorf JA, Cohee LM, Coalson JE, et al. School-age children are a reservoir of malaria infection in Malawi. PLoS One. 2015;10(7):e0134061. doi:10.1371/journal.pone.0134061
    OpenUrlCrossRefPubMed
  37. 37.↵
    Betson M, Clifford S, Stanton M, Kabatereine NB, Stothard JR. Emergence of nonfalciparum Plasmodium infection despite regular artemisinin combination therapy in an 18-month longitudinal study of Ugandan children and their mothers. The Journal of Infectious Diseases. 2018;217(7):1099–1109. doi:10.1093/infdis/jix686
    OpenUrlCrossRef
  38. 38.
    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(7):e2775. doi:10.1371/journal.pone.0002775
    OpenUrlCrossRefPubMed
  39. 39.↵
    Hawadak J, Nana RRD, Singh V. Global trend of Plasmodium malariae and Plasmodium ovale spp. malaria infections in the last two decades (2000–2020): a systematic review and meta-analysis. Parasites & Vectors. 2021;14:1–14. doi:10.1186/s13071-021-04797-0
    OpenUrlCrossRef
  40. 40.↵
    Parr JB, Kieto E, Phanzu F, et al. Analysis of false-negative rapid diagnostic tests for symptomatic malaria in the Democratic Republic of the Congo. Sci Rep. 2021;11(1):6495. doi:10.1038/s41598-021-85913-z
    OpenUrlCrossRef
  41. 41.↵
    Baltzell KA, Shakely D, Hsiang M, et al. Prevalence of PCR detectable malaria infection among febrile patients with a negative Plasmodium falciparum specific rapid diagnostic test in Zanzibar. Am J Trop Med Hyg. 2013;88(2):289–291. doi:10.4269/ajtmh.2012.12-0095
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    Groger M, Veletzky L, Lalremruata A, et al. Prospective clinical and molecular evaluation of potential Plasmodium ovale curtisi and wallikeri relapses in a high-transmission setting. Clinical Infectious Diseases. 2019;69(12):2119–2126. doi:10.1093/cid/ciz131
    OpenUrlCrossRef
Back to top
Posted June 08, 2022.
Download PDF

Supplementary Material

Data/Code
Email
Non-falciparum malaria infections are as prevalent as P. falciparum among Tanzanian schoolchildren
Rachel Sendor, Cedar L. Mitchell, Frank Chacky, Ally Mohamed, Lwidiko E. Mhamilawa, Fabrizio Molteni, Ssanyu Nyinondi, Bilali Kabula, Humphrey Mkali, Erik J. Reaves, Naomi Serbantez, Chonge Kitojo, Twilumba Makene, Thwai Kyaw, Meredith Muller, Alexis Mwanza, Erin L. Eckert, Jonathan B. Parr, Jessica T. Lin, Jonathan J. Juliano, Billy Ngasala
medRxiv 2022.06.07.22275625; doi: https://doi.org/10.1101/2022.06.07.22275625
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area