Geospatial analysis of malaria burden in Kagera region, North-western Tanzania using health facility and community survey data

Daniel A. Petro, Nyimvua Shaban, Sijenunu Aaron, Frank Chacky, Samwel Lazaro, Maciej F. Boni, Deus S. Ishengoma
doi: https://doi.org/10.1101/2024.06.11.24308780

Abstract

Generating evidence of malaria burden in areas with confirmed artemisinin partial resistance (ART-R) is the first step before developing a response strategy to the resistance. In this study, we assessed the burden of malaria in Kagera region (with recently confirmed ART-R) by geospatial analysis using data from health facilities (HFs) and community surveys from 2017 to 2023. In 2983717/8124363 patients from HFs, rapid diagnostic test positivity rate was 36.7% (range: 0-80%) and was similar in patients aged under-fives (33.1%, range: 0-79%) and ≥5 years (33.7%, range: 3-80%). The prevalence of malaria was 10.0% (range: 0-40.5%, n=84999/ 853761) in pregnant women but it was 26.1% (3409/13065) in school children (range: 0-78.4%). We identified hotspots and coldspots, and persistently high burden in 69/192 wards, school children, and patients aged ≥5 years, providing evidence for planning and developing an ART-R response for the region.

Tanzania is one of the priority countries defined by the World Health Organization (WHO) as a high-burden high-impact country; it accounted for an estimated 3.2% of all malaria cases and 4.4% of deaths reported worldwide in 2022 (1). In Tanzania, about 96% of all cases are caused by Plasmodium falciparum, and the remaining are due to other species (P. malariae and P. ovale) individually or in mixed infections (2–4). The National Malaria Control Programme’s (NMCP) strategy is to reduce the malaria burden and progress to substantially control and eliminate using different effective interventions. The main strategies include integrated malaria vector control (mainly with insecticide-treated nets), case management, preventive therapies, and other supportive strategies (5). However, progress to malaria elimination is under key biological threats such as the emergence and spread of insecticide (6) and antimalarial drug resistance (7), emergence of parasites with histidine-rich protein 2/3 (hrp2/3) gene deletions (8), and invasive vector species (Anopheles stephensi) in some malaria-endemic countries in Africa (1).

Malaria transmission in Tanzania is heterogeneous with 93% of the population living in areas where transmission occurs (5). Regions with very low and low transmission risk are located in the central, northeastern, and southwestern parts surrounded on both sides by moderate transmission regions and high transmission in the northwestern, southern, and western parts (9). Various factors including the types of mosquito vectors, environmental conditions, and other human-related factors are responsible for the current pattern of malaria transmission in Tanzania (10–12).

Over the past four decades, malaria-endemic countries in sub-Saharan Africa have reported resistance of P. falciparum to all major antimalarial drugs, including chloroquine, sulfadoxine-pyrimethamine (SP) (7), and artemisinin derivatives and their currently used partner drugs (13,14). Nevertheless, artemisinin-based combination therapies (ACTs) are still highly effective against parasites of African origin (7). However, artemisinin partial resistance (ART-R) has been confirmed in Rwanda (15,16), Uganda (17), Tanzania (18) and Eritrea (19), and one particular genotype has only been detected in the Horn of Africa, in Ethiopia (20), Somalia (21) and Sudan (22,23). In addition, the mutation that originated in Rwanda has been also reported in Kenya (24). In Tanzania, little was known about ART-R until in 2021 when a country-wide survey through the project on molecular surveillance of malaria in Tanzania (MSMT) reported a focus in Kagera region only, with high frequency (0.077) of parasites with mutations in Kelch13 gene (K13) and as high as 0.22 in Karagwe District Council (DC) which are linked to ART-R in Africa (25). Based on WHO’s criteria (7), a follow-up therapeutic efficacy study conducted in 2022 confirmed the presence of ART-R in the region although the tested drugs still had high efficacy exceeding 98.0% (18).

As per WHO’s strategy to respond to ART-R in Africa, when resistance is confirmed, new strategies to contain, prevent, and mitigate the spread and consequences of ART-R must be considered urgently (26). As the first step towards the development of a response strategy, WHO recommends mapping and providing a detailed profile of the disease burden and its spread in different areas (27). Therefore, this study aimed to investigate the temporal trends and spatial patterns of malaria burden in Kagera region and identify areas at the highest or lowest risk of malaria transmission as an initial stage towards the development of a response to ART-R in Tanzania. The findings from this study provide the government and other stakeholders with evidence to support the designing, deployment, and implementation of interventions to monitor, prevent, and contain the spread of ART-R within the region and to other areas of Tanzania.

Material and Methods

Study design and site

This study utilized aggregated secondary data collected retrospectively from health facilities (HFs) and community surveys in Kagera region of north-western Tanzania. The region consists of 8 councils, 192 wards, and 734 villages covering an area of approximately 35868 square kilometers. A description of the region’s geographical boundaries and the source of its economy is provided in Appendix A.

Data collation, management, and analysis

The data for this were downloaded from the District Health Information System – version 2 (DHIS2), a web-based software platform for reporting, analyzing, storing, and disseminating health data. This (DHIS2) is an electronic format for the paper-based Health Management Information System of the Tanzania Ministry of Health which is used by HFs to collect and report their routine data (5). The description of how the information is entered and maintained as well as quality control of the data is given elsewhere (9,28). The data used in this study included monthly malaria testing of pregnant women attending first antenatal care (ANC) clinics, symptomatic malaria patients seeking care at HFs (aged under-fives and ≥5 years) both from 2017 to 2023 and school children through school malaria parasitological surveys (SMPS) of 2015, 2017, 2019, and 2021 from the 37 selected wards in all councils of Kagera region.

Using Excel (Office 2013) and R software version 4.3.2: (www.r-project.org), data cleaning and preprocessing were done to ensure accuracy and consistency (Appendix B). Monthly data of total malaria tests performed and positive tests (by RDTs) reported by all HFs in the region were aggregated to provide annualized estimates. Malaria test-positivity rate (TPR) and prevalence were defined as the proportion of positive malaria tests among all malaria tests in patients aged under-fives and ≥5 years, and pregnant women, respectively.

For school children, the average prevalence for each council and ward was used. Using R software, data analysis was done for overall, per year and study group, council-wise, and at ward level, and the results presented in tables, figures or text. Python plotting tools (pandas v2.2.0, matplotlib v3.8.2, numpy v1.26.3) (Python.org) were used to show temporal trends of TPR and prevalence in the councils.

Kruskal-Wallis test was used to assess for any significant differences in malaria burden among councils, wards, and across years. Data on administrative boundaries of Tanzania at lower levels and its neighboring countries were downloaded from various sources including the Tanzanian National Bureau of Statistics (NBS), gadm.org, and naturalearthdata.com accessed via the R package rnaturalearth v1.0.1. These were used to create the study site and used for choropleth mapping, spatial autocorrelation (SA), and hotspot analysis. The spatial distribution of malaria burden in the region was mapped using R packages ggplot2 v3.5.1 and sf v1.0.15.

To investigate the spatial and temporal patterns of malaria burden in the region, the Global Moran’s Index was calculated for each group and year. The value of this index ranges from −1 (negative SA) to +1 (positive SA) (29,30). For hotspot analysis, local spatial statistical methods were applied through a spdep package v1.3.3 found in R software. The description of wards and HFs included in the analysis is found in Appendix B and Figure 1, while the cut-offs applied to characterize malaria prevalence were; <1% (very low), 1 - <5% (low), 5 -<30% (moderate), and ≥30% (high) as previously described (31,32). For malaria TPR, the cut-offs were <5% (very low), 5-<15% (low), 15-<30% (moderate), and ≥30% (high) (Appendix C and Table 1). The formulae and illustration of how SA and hotspot analysis were conducted can be found elsewhere (29,30) and in Appendix D.

Figure 1.
Figure 1.

A map showing the 26 regions of Mainland Tanzania (white) and Kagera region (red), Lake Victoria and Indian Ocean (lightblue), and neighboring countries (gray).

View this table:
Table 1. Results of global spatial autocorrelation of malaria test positivity rates in patients aged under-fives and ≥ 5 years in Kagera region

Results

Reported malaria cases from Health Facilities (HFs)

In the HFs in Kagera region, 36.7% of patients (range: 0-80%; n = 2983717/8124363) with malaria-like symptoms seeking care from 2017 to 2023, had RDT-positive results. The TPRs were similar in patients aged under-fives (33.1%, range: 0-79%) and ≥5 years (33.7%, range: 3-80%), with clear geographic and temporal variations (p<0.001). Among under-fives, most of the wards in the high stratum (TPR ≥30%) were from rural areas while wards in the very low stratum (TPR <5%) came from urban areas (Appendix Table 2). Furthermore, five councils had >10% of the wards with consistently high TPR (≥30%) throughout the years. These included 6/17 (35.3%) wards in Biharamulo DC, 12/43 (27.9%) in Muleba DC, 5/22 (22.7%) in Ngara DC, 4/23 (17.4%) in Karagwe DC and 4/29 (13.8%) in Bukoba DC. In Kyerwa DC, only 2/24 (8.3%) wards had consistently higher TPRs while only one ward (Kashai in Bukoba MC) recorded the lowest TPR (<5%) (Appendix Table 4). Among patients aged ≥5 years, the pattern was similar to that of under-fives, with a large number of wards observed in the high stratum throughout the years (Appendix Table 3).

View this table:
Table 2. Year-wise results of global spatial autocorrelation of malaria prevalence in Kagera region
View this table:
Table 3. Number of wards per council classified as malaria hotspots in patients aged under-fives and ≥5 years, and pregnant women in Kagera region at 95 and 99% confidence levels
View this table:
Table 4. Number of wards per council classified as malaria coldspots in patients aged under-fives and ≥5 years, and pregnant women

Similar to under-fives, the same five councils had >10% of the wards with consistently high TPRs (≥30%) over the years and these included, 7/17 (41.2%) wards in Biharamulo DC, 7/ 3 (30.4%) in Karagwe DC, 11/43 (25.6%) in Muleba DC, 5/22 (22.7%) in Ngara DC, and 4/29 (13.8%) wards in Bukoba DC. Only 1/24 (4.2%) ward in Kyerwa DC had consistently higher TPR (≥30%) (Appendix Table 4) (Figure 2B).

Figure 2.
Figure 2.

Spatial distribution of malaria test positivity rates (TPRs) in (A) under-fives (B) patients aged ≥5 years at ward level in Kagera region from 2017 to 2023. White areas are wards excluded from the analysis

Monthly variations of test positivity rates (TPRs) at health facility level

Monthly malaria TPR was higher (≥30%) in patients aged ≥5 years compared to under-fives in all councils (in Muleba DC, Ngara DC, Biharamulo DC, Bukoba DC, and Karagwe DC) and the pattern was similar in all years, except in Bukoba MC which recorded the lowest TPR (<5%) in both age groups. In all councils except Bukoba MC and Missenyi DC, TPRs less than 20% were rarely observed (Figure 3).

Figure 3.
Figure 3.

Monthly test positivity rates (TPRs) at health facility level in eight councils of Kagera region in patients aged under-fives (green) and ≥5 years (black) from 2017 to 2023.

Malaria prevalence among pregnant women

From 2017 to 2023, the prevalence of malaria in pregnant women was 10.0% (range: 0-40.5%, n=84999/ 853761), with significant differences among wards and across years (p<0.001). The prevalence was consistently higher in Biharamulo DC, Karagwe DC, Muleba DC, and Kyerwa DC, and low in Bukoba MC (<1%) throughout the years (Figure 4).

Figure 4.
Figure 4.

Trends of malaria prevalence in pregnant women by council, in Kagera region from 2017 to 2023. The pink-shaded regions show the Clopper-Pearson 95% confidence bands.

The largest number of wards (76.8%, n=116/151) were located in the moderate stratum while the lowest (0.6%, n=1/159) was in the high-risk stratum of prevalence (Appendix Table 5). In addition, only 3/8 (37.5%) councils had ≥5% of the wards with higher prevalence (≥30%) (including 3/17 (17.6%) wards in Biharamulo DC, 4/43 (9.3%) in Muleba DC and 2/23 (8.7%) wards in Karagwe DC) while Kyerwa DC (n=1/24, 4.2%) and Bukoba DC (n=1/29, 3.4%) had <5% of the wards with high prevalence. At least 5% of the wards from Bukoba DC (n=2/29, 6.9%) and Missenyi DC (n = 1/20, 5%) consistently recorded the lowest prevalence (Appendix Table 6) (Figure 5).

Figure 5.
Figure 5.

Spatial distribution of malaria prevalence in pregnant women at ward level in Kagera region from 2017 to 2023. White areas are wards excluded from the analysis.

Malaria prevalence in school children

The prevalence of malaria in school children was 26.1% (range: 0-78.4%; n=3409/1306) (Appendix E). A large number of wards (19/34, 55.9%) were located in the high stratum (≥30%) and few wards (1/34, 2.9%) came from areas with very low stratum (<1%) (Appendix Table 7). Furthermore, two wards from Ngara DC and two from Bukoba DC had persistently higher prevalence (≥30%) throughout the study period, while no ward recorded the lowest prevalence (<1%) for all years (Appendix Table 8) (Figure 6).

Figure 6.
Figure 6.

Spatial distribution of malaria prevalence in school children at ward level in Kagera region for 2015, 2017, 2019, and 2021 surveys. White areas are wards not selected for the survey or excluded from the analysis.

Overall malaria test positivity rates at ward level

Among patients aged between under-fives and ≥5 years, a declining trend was observed in under-fives and started to shift to older patients at the same time (Appendix Table 9). In both age groups, an increasing trend was observed in 2023 (Appendix Figure 2).

Overall malaria prevalence at ward level

Malaria prevalence was higher in school children than in pregnant women, a trend that was similar across years (Appendix Figure 3). Furthermore, prevalence in pregnant women showed a decreasing trend with minor variability before 2022 but increased by approximately 15% in 2023 (Appendix Table 9).

Global spatial autocorrelation analysis of malaria burden

Malaria burden in the region was clustered in under-fives, patients aged ≥5 years, and pregnant women. Councils/wards with high/low burden were surrounded by other areas with high/low values as well for all years and the observed patterns were statistically significant (p<0.05) in all groups except school children (Tables 1 and 2).

Hotspot analysis of malaria burden

Malaria hotspots were concentrated in the southern, central, and northwestern parts of Kagera region. Four councils (Ngara DC, Biharamulo DC, Muleba DC, and Bukoba DC) had >17% of the wards classified as malaria hotspots in patients aged under-fives and ≥5 years, and pregnant women while Karagwe DC and Kyerwa DC had <17% of wards as malaria hotspots in the same groups. The number of wards classified as hotspots of malaria in each group at 95% and 99% confidence levels (CLs) is presented in Table 3 (Appendix Table 10).

On the other hand, malaria coldspots were situated in the northeast and few in the remaining parts of the region. Three councils (Bukoba MC, Missenyi DC, and Bukoba DC) had >20% of the wards classified as malaria coldspots in patients aged under-fives, ≥5 years, and pregnant women while Kyerwa DC and Muleba DC had <20% of the wards. The number of wards classified as malaria coldspots in each group at 95% and 99% CLs is presented in Table 4 (Appendix Table 11).

The locations of malaria hotspots and coldspots at 95% and 99% CLs in under-fives, ≥5 years, and pregnant women respectively are indicated in Appendix Figures S4 to S6.

Discussion

This study was conducted as an initial stage to respond to ART-R as recommended by WHO (26) and aimed to undertake an assessment to generate, and provide evidence of the current burden of malaria in Kagera region where ART-R was recently confirmed (18). We used spatial statistical and disease mapping techniques to assess the local burden of malaria and identify potential malaria hotspots in Kagera region to support NMCP’s efforts of developing a response strategy against ART-R. Using aggregated secondary data from multiple sources, we showed that over 35% (69/192) of all wards in Kagera region, mostly from rural areas, had a high malaria burden, with these wards being clustered across years. The TPR was higher in patients aged ≥5 years than under-fives and this is comparable to previous studies conducted in Tanzania (11,12). The findings reported here provide important evidence to support the designing and implementation of mitigation and response strategies to recently confirmed ART-R in Kagera region and Tanzania in general.

Consistent with previous studies (12), this study reported high temporal and spatial heterogeneity of malaria burden at council and ward levels within the region. Most of the wards, especially from rural areas, recorded a high malaria burden for most of the years; while some wards mainly in urban and mixed/rural areas recorded the lowest burden (see Appendix Tables 4,6, and 8). This finding was consistent across years and study groups possibly due to conditions supporting malaria transmission in rural compared to urban areas as previously reported (33–35). The higher burden of malaria in rural areas is normally attributed to several factors including the presence of favorable conditions for malaria transmission in rural areas such as stagnant water bodies that tend to create desirable breeding sites for mosquito vectors (35). In addition, a large proportion of poor people tend to reside in rural areas and normally live in poorly constructed houses which support higher mosquito biting rates and transmission of malaria (12,36,37). To reduce the burden in these areas, various interventions such as distribution of insecticide-treated nets (ITNs), house improvement strategies, bio-larviciding in areas that are surrounded by water bodies, and behavioral change and communication strategies need to be implemented, as these have been reported to significantly reduce malaria transmission and the burden (38,39). On the other hand, the low burden of malaria in urban areas observed in this study is consistent with previous studies (33,35). This is possibly due to the high socio-economic status of the urban population and, the availability and access to malaria treatment and control measures when compared to rural areas (35).

A significant clustered spatial pattern of malaria burden reported in patients aged under-fives and ≥5 years, and pregnant women within the region aligns with other studies conducted elsewhere (29,30). The reports of areas with persistently high malaria burden in the southern, central, and north-western parts of the region as reported in this study may be attributed to the presence of similar risk factors such as environmental conditions, socio-economic status, healthcare infrastructures, and the presence of suitable mosquito breeding sites that tend to create foci of malaria transmission (10,12,40). Hence, future interventions in the region to reduce the local disease burden need to be designed based on these results to target the areas that have been identified to have a persistently high burden of malaria.

Malaria prevalence in Kagera region was high in school children compared to pregnant women, reaching as high as 78.4% in some wards. In previous studies, it was reported that school children spend most of their time in the evening playing outside and less often use bed nets, which tends to increase their chances of mosquito bites and hence increase the risk of malaria infections (41–43). To reduce malaria transmission in school children, various initiatives have been undertaken by NMCP and other stakeholders to protect this risky group including the distribution of ITNs through different channels like the school net program (5,44). However, the prevalence of malaria infections in this asymptomatic group is still high (12,45). Future studies are needed to explore more reasons and factors for this recurring observation in the region and the whole country and whether this will support the spread of the reported ART-R. On the other hand, pregnant women in Tanzania are normally given two or more doses of SP for intermittent preventive treatment as well as ITNs during their first ANC visits in HFs, which might have significantly contributed to reducing the risk of malaria infections in this group (28). However, the increasing trend of malaria prevalence in pregnant women observed in 2023 needs further follow-up to uncover all possible reasons for this change.

The overall malaria TPR at the ward level was high (≥30%) in both under-fives and patients aged ≥5 years with minor differences. These findings are comparable to what has been reported by previous studies in Tanzania and elsewhere (11,12). The year-wise comparison of TPR revealed a shift of TPR from under-fives to ≥5 years, similar to previous studies (46). This may be attributed to various efforts that have been undertaken by NMCP to implement effective malaria interventions in vulnerable groups like under-fives while these interventions may not have been widely considered in older age groups (11,28). This might have led to the shift of malaria infections to individuals aged ≥5 years, with the peak of malaria prevalence reported in children aged between 5 and 15 years (11). In addition, malaria TPR in under-fives and ≥5 years showed an increasing trend in 2023. An in-depth analysis to uncover possible reasons for these epidemiological changes in the region is required so that appropriate measures can be planned, now that ART-R has been confirmed.

Malaria transmission in Kagera region was observed throughout the year with some of the months recording high values of TPR compared to others. This indicates that malaria transmission in the region was not seasonal; an observation that can signify the presence of other factors like population movements, changes in malaria control, or unexpected increases/decreases in breeding sites or mosquito populations (12,47,48). Causal effects for these factors are in general difficult to establish, future studies are needed to explore and uncover the key drivers of persistent malaria transmission in this malaria-endemic region in the country.

The findings presented here had some limitations that may be addressed in future studies. Since this study included HFs with a reporting rate of 50% or more, future studies may consider including all HFs during analysis and use different geostatistical methods to account for monthly missing data, differences in utilization of HFs, treatment-seeking behaviors, and distance to HFs. In addition, the HF data did include personal information such as age and sex, and this limited additional analysis. For instance, previous studies showed a high burden of malaria in school children (aged 5 - 15 years) but this could not be teased out because all patients aged ≥5 years were put in one group. Future studies may also include populations at risk and use incidence to track changes in the disease burden within the region.

Overall, malaria infections in the region were higher in most of the wards from rural areas compared to urban areas, among school children compared to pregnant women, and among patients aged ≥5 years than in under-fives. The observed heterogeneity of malaria burden underscores prioritization of Kagera region for additional malaria control mainly due to the currently circulating P. falciparum parasites with K13 561H genotypes in the region (18). As early and preemptive action is preferred to late response (49), identifying specific resistance containment strategies (50) for 2024 and 2025 – in addition to general malaria control strategies– will be critical for slowing down the spread of artemisinin-resistant parasites from Kagera region to other parts of Tanzania.

Biographical Sketch

Mr. Petro is a doctoral student at the University of Dar es Salaam, Dar es Salaam. His research interests focus on malaria epidemiology.

Data Availability

Data used in this study are not publicly available and were obtained with a request from MoH through NMCP. Restrictions apply to the availability of this data and permission can be obtained with a reasonable request from the Permanent Secretary - Ministry of Health of Mainland Tanzania.

Appendices

Appendix A. Description of the study site

Appendix B. Administrative units and health facilities in Kagera region that were included in the analysis

Appendix C: Stratification of malaria burden in Kagera region

Appendix D: Description of spatial autocorrelation analysis

Appendix E: Description of wards and schools selected for school malaria parasitological surveys

Appendix Figure 1. Descriptive summary of the wards and health facilities included in the analysis

Appendix Figure 2 Comparison of malaria test positivity rates in Kagera region from 2017 to 2023 based on the data from health facilities

Appendix Figure 3. Comparison of malaria prevalence in Kagera region from 2017 to 2023 based on the data from pregnant women and school children

Appendix Figure 4 Locations of hotspots and coldspots of malaria in under-fives at ward level in Kagera region from 2017 to 2023

Appendix Figure 5 Locations of hotspots and coldspots of malaria in patients aged ≥5 years at ward level in Kagera region from 2017 to 2023

Appendix Figure 6. Locations of hotspots and coldspots of malaria in pregnant women at ward level in Kagera region from 2017 to 2023

Appendix Table 1. Cut-offs adopted to categorize councils/wards in Kagera region into different risk strata

Appendix Table 2. Number of wards by year in each stratum of malaria test positivity rates in under-fives in Kagera region from 2017 to 2023

Appendix Table 3. Number of wards by year in each stratum of malaria test positivity rates in patients aged ≥ 5 years in Kagera region from 2017 to 2023

Appendix Table 4. Ward names by council that had a high, moderate, low, and very low malaria test positivity rates in under-fives and ≥5 years from 2017 to 2023, Kagera region

Appendix Table 5. Number of wards by year in each stratum of malaria prevalence in pregnant women in Kagera region from 2017 to 2023

Appendix Table 6. Ward names by council that had a high, moderate, low, and very low malaria prevalence in pregnant women from 2017 to 2023, Kagera region

Appendix Table 7. Number of wards by year in each stratum of malaria prevalence in school children in Kagera region from 2017 to 2023

Appendix Table 8. Ward names by council that had a high, moderate, low, and very low malaria prevalence in school children from 2017 to 2023, Kagera region

Appendix Table 9. Overall malaria burden at ward level in Kagera region by years and study groups from 2015 to 2023

Appendix Table 10. Names of hotspots of malaria burden at ward level in Kagera region at 95% and 99% confidence levels from 2017 to 2023

Appendix Table 11. Names of coldspots of malaria burden at ward level in Kagera region at 95% and 99% confidence levels from 2017 to 2023

Acknowledgment

The authors would like to thank the MoH through NMCP, and all members of the Department of Mathematics, University of Dar es Salaam, and MSMT project team for their participation and invaluable support to make this work complete.

This study was partially funded by the National Institute for Medical Research (NIMR) under the MSMT project. The MSMT project is funded in whole, by the Bill & Melinda Gates Foundation [grant number INV 02202]. Under the grant conditions of the Foundation, a Creative Commons Attribution 4.0 Generic License has already been assigned to the author of the accepted Manuscript version that might arise from this submission.

PAP and DSI designed the study, and SA, FC and SL supported the acquisition and cleaning of the data. DAP conducted the initial analysis under the supervision of DSI and NS. DSI, NS, and MFB were involved in the interpretation of the results. DAP and MFB performed a final analysis of the data under the guidance of DSI, and DAP drafted the initial manuscript. All co-authors contributed to the revision of the initial draft and subsequent versions of the manuscript. All authors read and approved the final manuscript, and declared no competing interest.

Footnotes

  • Figure 6 updated to include data from school children survey conducted in 2021 and the manuscript has been formatted to Emerging Infectious Diseases Journal

References

  1. 1.↵
    World Health Organization. World Malaria Report 2023 [Internet]. World Health Organization; 2023. Available from: https://www.who.int/publications/i/item/9789240086173
  2. 2.↵
    Popkin Hall ZR, Seth MD, Madebe RA, Budodo R, Bakari C, Francis F, et al. Malaria species positivity rates among symptomatic individuals across regions of differing transmission intensities in Mainland Tanzania. J Infect Dis [Internet]. 2023 Nov 22; Available from: doi:10.1093/infdis/jiad522
    OpenUrlCrossRef
  3. 3.
    Popkin-Hall ZR, Seth MD, Madebe RA, Budodo R, Bakari C, Francis F, et al. Prevalence of non-falciparum malaria infections among asymptomatic individuals in four regions of Mainland Tanzania. Parasit Vectors [Internet]. 2024;17(1):153. Available from: doi:10.1186/s13071-024-06242-4
    OpenUrlCrossRef
  4. 4.↵
    Sendor R, Mitchell CL, Chacky F, Mohamed A, Mhamilawa LE, Molteni F, et al. Similar Prevalence of Plasmodium falciparum and Non-P. falciparum Malaria Infections among Schoolchildren, Tanzania1. Emerg Infect Dis [Internet]. 2023 Jun;29(6):1143–53. Available from: doi:10.3201/eid2906.221016
    OpenUrlCrossRef
  5. 5.↵
    National Malaria Control Programme (NMCP). “National Malaria Strategic Plan 2021-2025: Transitioning to Malaria Elimination Phase”, Tanzania Ministry of Health, Community Development, Gender, Elderly and Children. Dodoma; 2020 p. 1–179.
  6. 6.↵
    World Health Organization. Global-report-on-insecticide-resistance-in-malaria-vectors-2010-2016 [Internet]. 2018. Available from: https://fundacionio.com/wp-content/uploads/2020/11/Global-report-on-insecticide-resistance-in-malaria-vectors-2010-2016.pdf
  7. 7.↵
    World Health Organization. Report on antimalarial drug efficacy, resistance and response: 10 years of surveillance (2010-2019). World Health Organization; 2020.
  8. 8.↵
    Watson OJ, Tran TNA, Zupko RJ, Symons T, Thomson R, Visser T, et al. Global risk of selection and spread of Plasmodium falciparum histidine-rich protein 2 and 3 gene deletions. medRxiv. 2023;
  9. 9.↵
    Thawer SG, Golumbeanu M, Munisi K, Aaron S, Chacky F, Lazaro S, et al. The use of routine health facility data for micro-stratification of malaria risk in mainland Tanzania. Malar J [Internet]. 2022 Nov;21(1):345. Available from: doi:10.1186/s12936-022-04364-7
    OpenUrlCrossRef
  10. 10.↵
    Rumisha SF, Shayo EH, Mboera LEG. Spatio-temporal prevalence of malaria and anaemia in relation to agro-ecosystems in Mvomero district, Tanzania. Malar J. 2019;18(1):1–14.
    OpenUrlCrossRef
  11. 11.↵
    Ishengoma DS, Mmbando BP, Mandara CI, Chiduo MG, Francis F, Timiza W, et al. Trends of Plasmodium falciparum prevalence in two communities of Muheza district North-eastern Tanzania: correlation between parasite prevalence, malaria interventions and rainfall in the context of re-emergence of malaria after two decades of progressive. Malar J [Internet]. 2018 Jul;17(1):252. Available from: doi:10.1186/s12936-018-2395-1
    OpenUrlCrossRefPubMed
  12. 12.↵
    Mandai SS, Francis F, Challe DP, Seth MD, Madebe RA, Petro DA, et al. High prevalence and risk of malaria among asymptomatic individuals from villages with high rates of artemisinin partial resistance in Kyerwa district, North-western Tanzania. medRxiv. 2023;2010–23.
  13. 13.↵
    Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361(5):455–67.
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359(24):2619–20.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    Uwimana A, Legrand E, Stokes BH, Ndikumana JLM, Warsame M, Umulisa N, et al. Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda. Nat Med. 2020;26(10):1602–8.
    OpenUrlCrossRefPubMed
  16. 16.↵
    Uwimana A, Umulisa N, Venkatesan M, Svigel SS, Zhou Z, Munyaneza T, et al. Association of Plasmodium falciparum kelch13 R561H genotypes with delayed parasite clearance in Rwanda: an open-label, single-arm, multicentre, therapeutic efficacy study. Lancet Infect Dis. 2021;21(8):1120–8.
    OpenUrlCrossRefPubMed
  17. 17.↵
    Balikagala B, Fukuda N, Ikeda M, Katuro OT, Tachibana SI, Yamauchi M, et al. Evidence of artemisinin-resistant malaria in Africa. N Engl J Med. 2021;385(13):1163–71.
    OpenUrlCrossRefPubMed
  18. 18.↵
    Ishengoma DS, Mandara CI, Bakari C, Fola AA, Madebe RA, Seth MD, et al. Evidence of artemisinin partial resistance in North-western Tanzania: clinical and drug resistance markers study. medRxiv [Internet]. 2024 Feb 1; Available from: doi:10.1101/2024.01.31.24301954
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    Mihreteab S, Platon L, Berhane A, Stokes BH, Warsame M, Campagne P, et al. Increasing Prevalence of Artemisinin-Resistant HRP2-Negative Malaria in Eritrea. N Engl J Med [Internet]. 2023 Sep 28;389(13):1191–202. Available from: doi:10.1056/NEJMoa2210956
    OpenUrlCrossRef
  20. 20.↵
    Fola AA, Feleke SM, Mohammed H, Brhane BG, Hennelly CM, Assefa A, et al. Plasmodium falciparum resistant to artemisinin and diagnostics have emerged in Ethiopia. Nat Microbiol [Internet]. 2023 Oct;8(10):1911–9. Available from: doi:10.1038/s41564-023-01461-4
    OpenUrlCrossRef
  21. 21.↵
    Jalei AA, Na-Bangchang K, Muhamad P, Chaijaroenkul W. Monitoring antimalarial drug-resistance markers in Somalia. Parasites Hosts Dis [Internet]. 2023 Feb;61(1):78–83. Available from: doi:10.3347/PHD.22140
    OpenUrlCrossRef
  22. 22.↵
    Mohamed AO, Hussien M, Mohamed A, Suliman A, Elkando NS, Abdelbagi H, et al. Assessment of Plasmodium falciparum drug resistance molecular markers from the Blue Nile State, Southeast Sudan. Malar J [Internet]. 2020 Feb 18;19(1):78. Available from: doi:10.1186/s12936-020-03165-0
    OpenUrlCrossRef
  23. 23.↵
    Molina-de la Fuente I, Sagrado Benito MJ, Ousley J, Gisbert F de B, García L, González V, et al. Screening for K13-Propeller Mutations Associated with Artemisinin Resistance in Plasmodium falciparum in Yambio County (Western Equatoria State, South Sudan). Am J Trop Med Hyg [Internet]. 2023 Nov 1;109(5):1072–6. Available from: doi:10.4269/ajtmh.23-0382
    OpenUrlCrossRef
  24. 24.↵
    Owoloye A, Olufemi M, Idowu ET, Oyebola KM. Prevalence of potential mediators of artemisinin resistance in African isolates of Plasmodium falciparum. Malar J [Internet]. 2021 Dec 2;20(1):451. Available from: doi:10.1186/s12936-021-03987-6
    OpenUrlCrossRef
  25. 25.↵
    Juliano JJ, Giesbrecht DJ, Simkin A, Fola AA, Lyimo BM, Pereus D, et al. Country wide surveillance reveals prevalent artemisinin partial resistance mutations with evidence for multiple origins and expansion of high level sulfadoxine-pyrimethamine resistance mutations in northwest Tanzania. medRxiv [Internet]. 2023 Nov 30; Available from: doi:10.1101/2023.11.07.23298207
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    World Health Organization. Strategy to respond to antimalarial drug resistance in Africa 2022 [Internet]. Vol. 04. Geneva, Switzerland: World Health Organization; 2022. Available from: https://www.who.int/publications/i/item/9789240060265%0A https://cdn.who.int/media/docs/default-source/malaria/who-antimalarial-drug-resistance-strategy-for-consultation.pdf?sfvrsn=9d4eaa0_6
  27. 27.↵
    Zupko RJ, Nguyen TD, Ngabonziza JCS, Kabera M, Li H, Tran TNA, et al. Modeling policy interventions for slowing the spread of artemisinin-resistant pfkelch R561H mutations in Rwanda. Nat Med [Internet]. 2023 Nov;29(11):2775–84. Available from: doi:10.1038/s41591-023-02551-w
    OpenUrlCrossRef
  28. 28.↵
    National Malaria Control Programme (NMCP). Ministry Of Health, Community Development, Gender, Elderly and Children, “National Guidelines For Malaria Diagnosis, Treatment And Preventive Therapies.” Dodoma; 2020.
  29. 29.↵
    Tewara MA, Mbah-Fongkimeh PN, Dayimu A, Kang F, Xue F. Small-area spatial statistical analysis of malaria clusters and hotspots in Cameroon;2000–2015. BMC Infect Dis [Internet]. 2018 Dec 7;18(1):636. Available from: doi:10.1186/s12879-018-3534-6
    OpenUrlCrossRef
  30. 30.↵
    Saxena R, Nagpal BN, Das MK, Srivastava A, Gupta SK, Kumar A, et al. A spatial statistical approach to analyze malaria situation at micro level for priority control in Ranchi district, Jharkhand. Indian J Med Res [Internet]. 2012 Nov;136(5):776–82. Available from: https://www.ncbi.nlm.nih.gov/pubmed/23287124
    OpenUrl
  31. 31.↵
    Ministry of Health Community Development Gender Elderly & Children, National Malaria Control Programme. Supplementary malaria midterm strategic plan 2018–2020. 2019.
  32. 32.↵
    Ministry of Health, Community Development, Gender, Elderly and Children, National Malaria Control Programme. Consultative Malaria Expert Meeting Report 2018. 2018.
  33. 33.↵
    Nundu SS, Culleton R, Simpson SV, Arima H, Muyembe JJ, Mita T, et al. Malaria parasite species composition of Plasmodium infections among asymptomatic and symptomatic school-age children in rural and urban areas of Kinshasa, Democratic Republic of Congo. Malar J. 2021;20:1–13.
    OpenUrlCrossRefPubMed
  34. 34.
    Sultana M, Sheikh N, Mahumud RA, Jahir T, Islam Z, Sarker AR. Prevalence and associated determinants of malaria parasites among Kenyan children. Trop Med Health. 2017;45:1–9.
    OpenUrlCrossRef
  35. 35.↵
    Mbah CE, Ambe LA, Ngwewondo A, Kidzeru EB, Akwah L, Mountchissi C, et al. A Comparative Study of Asymptomatic Malaria in a Forest Rural and Depleted Forest Urban Setting during a Low Malaria Transmission and COVID-19 Pandemic Period. Biomed Res Int. 2022;2022.
  36. 36.↵
    Bofu RM, Santos EM, Msugupakulya BJ, Kahamba NF, Swilla JD, Njalambaha R, et al. The needs and opportunities for housing improvement for malaria control in southern Tanzania. Malar J. 2023;22(1):69.
    OpenUrl
  37. 37.↵
    Tusting LS, Rek J, Arinaitwe E, Staedke SG, Kamya MR, Cano J, et al. Why is malaria associated with poverty? Findings from a cohort study in rural Uganda. Infectious Diseases of Poverty [Internet]. 2016 Aug;5(1):78. Available from: doi:10.1186/s40249-016-0164-3
    OpenUrlCrossRef
  38. 38.↵
    McCann RS, van den Berg H, Diggle PJ, van Vugt M, Terlouw DJ, Phiri KS, et al. Assessment of the effect of larval source management and house improvement on malaria transmission when added to standard malaria control strategies in southern Malawi: study protocol for a cluster-randomized controlled trial. BMC Infect Dis [Internet]. 2017 Sep 22;17(1):639. Available from: doi:10.1186/s12879-017-2749-2
    OpenUrlCrossRef
  39. 39.↵
    Orkoh E, Efobi U. Effects of Behaviour Change Communication on Knowledge and Prevention of Malaria Among Women in Ghana. Eval Rev [Internet]. 2023 Aug 11;193841X231194565. Available from: doi:10.1177/0193841X231194565
    OpenUrlCrossRef
  40. 40.↵
    Dabaro D, Birhanu Z, Negash A, Hawaria D, Yewhalaw D. Effects of rainfall, temperature, and topography on malaria incidence in elimination targeted district of Ethiopia. Malar J. 2021;20(1):1–10.
    OpenUrlCrossRefPubMed
  41. 41.↵
    Finda MF, Moshi IR, Monroe A, Limwagu AJ, Nyoni AP, Swai JK, et al. Linking human behaviors and malaria vector biting risk in south-eastern Tanzania. PLoS One. 2019;14(6):e0217414.
    OpenUrl
  42. 42.
    Kihwele F, Gavana T, Makungu C, Msuya HM, Mlacha YP, Govella NJ, et al. Exploring activities and behaviors potentially increases school-age children’s vulnerability to malaria infections in south-eastern Tanzania. Malar J. 2023;22(1):293.
    OpenUrlCrossRefPubMed
  43. 43.↵
    Kihwele F, Gavana T, Makungu C, Msuya HM, Govella NJ, Mlacha YP, et al. Outdoor activities and behaviors increase the vulnerability of school-age children to malaria infections in south-eastern Tanzania [Internet]. Research Square. 2023. Available from: https://www.researchsquare.com/article/rs-2832075/v1
  44. 44.↵
    Yukich J, Stuck L, Scates S, Wisniewski J, Chacky F, Festo C, et al. Sustaining LLIN coverage with continuous distribution: the school net programme in Tanzania. Malar J [Internet]. 2020;19:158. Available from: doi:10.1186/s12936-020-03222-8
    OpenUrlCrossRefPubMed
  45. 45.↵
    Chacky F, Runge M, Rumisha SF, Machafuko P, Chaki P, Massaga JJ, et al. Nationwide school malaria parasitaemia survey in public primary schools, the United Republic of Tanzania [Internet]. Vol. 17, Malar J. 2023. Available from: doi:10.1186/s12936-018-2601-1
    OpenUrlCrossRef
  46. 46.↵
    Kigozi SP, Kigozi RN, Epstein A, Mpimbaza A, Sserwanga A, Yeka A, et al. Rapid shifts in the age-specific burden of malaria following successful control interventions in four regions of Uganda. Malar J [Internet]. 2020 Mar 30;19(1):128. Available from: doi:10.1186/s12936-020-03196-7
    OpenUrlCrossRef
  47. 47.↵
    Dida GO, Anyona DN, Abuom PO, Akoko D, Adoka SO, Matano AS, et al. Spatial distribution and habitat characterization of mosquito species during the dry season along the Mara River and its tributaries, in Kenya and Tanzania. Infectious diseases of poverty. 2018;7(1):1–16.
    OpenUrl
  48. 48.↵
    Fornace KM, Diaz AV, Lines J, Drakeley CJ. Achieving global malaria eradication in changing landscapes. Malar J. 2021;20(1):1–14.
    OpenUrlCrossRefPubMed
  49. 49.↵
    Boni MF, White NJ, Baird JK. The Community As the Patient in Malaria-Endemic Areas: Preempting Drug Resistance with Multiple First-Line Therapies. PLoS Med [Internet]. 2016 Mar;13(3):e1001984. Available from: doi:10.1371/journal.pmed.1001984
    OpenUrlCrossRefPubMed
  50. 50.↵
    Boni MF. Breaking the cycle of malaria treatment failure. Front Epidemiol [Internet]. 2022 Dec 14;2:1041896. Available from: doi:10.3389/fepid.2022.1041896
    OpenUrlCrossRef
Back to top
Posted June 23, 2024.
Download PDF

Supplementary Material

Data/Code
Email
Geospatial analysis of malaria burden in Kagera region, North-western Tanzania using health facility and community survey data
Daniel A. Petro, Nyimvua Shaban, Sijenunu Aaron, Frank Chacky, Samwel Lazaro, Maciej F. Boni, Deus S. Ishengoma
medRxiv 2024.06.11.24308780; doi: https://doi.org/10.1101/2024.06.11.24308780
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area