Abstract
Background Burkitt lymphoma (BL) is one of the most common childhood cancers in sub-Saharan Africa and aetiologically linked to malaria. However, evidence for an effect of malaria interventions on BL is limited. The aim of this study was to investigate the potential population-level association between large-scale rollout of insecticide-treated bednets (ITNs) in sub-Saharan Africa in the 2000s and BL incidence.
Methods We conducted a systematic review in Embase, Global Health and MEDLINE to assemble all available data on BL incidence rates in children aged 0-15 years in malaria-endemic sub-Saharan African countries between Jan 1, 1990 and Feb 27, 2023. We calculated pooled estimates of BL incidence in sub-Saharan Africa for the time periods before and after ITN introduction. We used model estimates of sub-national ITN use to assess the association between average population ITN use and childhood BL incidence rates in a negative binomial regression model.
Results We included 23 studies comprising 66 datapoints on BL incidence based on 5,226 cancer cases from locations with large-scale ITN use in 17 countries. BL rates were 44% (95% CI 12-64%) lower in the time period after ITN introduction compared to before. The pooled incidence rate of BL was 1.36 (95% CI 0.88-2.10) and 0.76 (95% CI 0.50-1.16) per 100,000 person-years before and after introduction of ITNs, respectively. After adjusting for potential confounders, a 1% increase in average ITN use in the population in the 10 years before BL data collection was associated with a 2% (95% CI 1-4%) reduction in BL incidence.
Conclusion These findings suggest that large-scale rollout of ITNs in the 2000s was associated with a reduction in childhood BL burden in sub-Saharan Africa. Although published data may not be representative of all incidence rates across sub-Saharan Africa, our study highlights a potential additional benefit of malaria control programmes.
Introduction
Burkitt lymphoma (BL) is an aggressive type of non-Hodgkin lymphoma and among the most common childhood cancers in sub-Saharan Africa (1). Infection with the ubiquitous Epstein-Barr virus (EBV) is an established risk factor for BL, and there is substantial epidemiological and biological evidence for an aetiological relationship between Plasmodium falciparum malaria and BL in sub-Saharan Africa (1, 2). Malaria still poses an important public health problem, with over 200 million clinical cases annually (3). Malaria parasite prevalence and the cumulative number of infections correlate with BL incidence on the population level (4, 5), and individual-level associations between malaria infection or anti-malarial antibody titres and BL risk have been reported (6).
This combined evidence raises the possibility that interventions against malaria might have downstream beneficial effects on reducing the burden of BL. In an uncontrolled malaria suppression trial in 1977-1982 in Tanzania, a substantial decline in BL incidence rates was observed during a mass chloroquine distribution campaign, starting before the trial years and subsequently rising to pre-trial levels (7). Two hospital-based case-control studies suggested vector control interventions could potentially decrease the risk of developing BL in African children [odds ratio (OR) 0.2 (95% CI 0.03-0.9) and 0.3 (95% CI 0.2–0.7) for children using mosquito nets or living in a household with insecticide use, respectively] (8, 9). In contrast, a recent large population-based case-control study reported an increased risk of BL in children having used a mosquito net the previous night (OR 2.4, 95% CI 1.8-3.2) and in children living in households exposed to indoor residual insecticide spraying in the past year (OR 2.3, 95% CI 1.8-3.1) (10). These mixed results may indicate a difficulty in accounting for residual confounding from malaria burden in the study of such a rare disease, and for the temporal relationship between intervention exposure and BL development (2, 10).
On the population level, the question of whether BL incidence has declined over time following scale-up of malaria interventions has not been investigated. The 2000s marked a period of unprecedented donor funding and the first time malaria interventions were rolled out in sub-Saharan Africa at a large scale (3). Mass campaigns focused on distributing insecticide-treated nets (ITNs) and artemisinin-based combination therapies (3). These control efforts reduced P. falciparum infection prevalence in sub-Saharan Africa by an estimated 50% between 2000 and 2015, with most of this decline attributable to ITNs (11).
In this study, we aimed to investigate the potential association between large-scale rollout of ITNs in the early 2000s and BL incidence in sub-Saharan Africa using an ecological study design. First, we conducted a systematic review to assemble all publicly available data on BL incidence rates in sub-Saharan African children since 1990. We calculated pooled estimates of BL incidence across malaria-endemic sub-Saharan African countries before and after large-scale ITN introduction and compared time trends in incidence in individual locations. Second, we used model estimates of sub-national ITN use (12) to investigate the association between average ITN use in the population and BL incidence in children.
Methods
Systematic review
The systematic review was conducted on cancer registry publications and in literature databases. The International Agency for Research on Cancer (IARC) collates data from population-based cancer registries globally; we assembled BL incidence data from publications identified through a search of
IARC-associated websites. We also searched EMBASE, Global Health and Medline databases up to 27th February 2023 using the Ovid platform for all studies on BL incidence rates in sub-Saharan Africa published since 1st January 1990 without language restrictions. Search terms included keywords for sub-Saharan African countries, BL and non-Hodgkin lymphoma and the outcome of interest (Table S1).
Epidemiological studies reporting the incidence rate of BL in children aged 0-15 years in malaria-endemic sub-Saharan African countries were included. Both population-based and hospital-based studies were considered for identification of cancer cases, but incidence rates had to apply to the general population. The outcome of interest was incident cancer cases diagnosed as BL in a defined population during a specified time period with a midpoint after 1990. Studies were eligible for inclusion if they reported the number of BL cases and either the person-time at risk or the crude incidence rate.
Articles were deduplicated and screened in the Covidence software. Screening was performed independently by two reviewers blinded to each other’s decision, and disagreements were resolved through discussion. One article reported BL incidence in 0-19 year while otherwise meeting inclusion criteria, but unpublished data in the 0-15 year group could not be obtained on request (4). Data was extracted independently by the two reviewers into a pre-validated Microsoft Office Excel template, with any differences resolved in a final assessment. Harmonised extracted data included study location, time, study design, data collection and ascertainment methods, age range, number of cases, person-years at risk, and incidence rate. Studies identified from the published literature were cross-checked against the dataset extracted from cancer registry publications for duplicates and overlap in data from the same registry. Their quality was scored based on three predefined criteria (data collection, case ascertainment, and calculation of person-time at risk).
This systematic review was registered in the PROSPERO register (reference CRD42023411961) and is reported according to PRISMA guidelines (13). Full methodological details can be found in the Supplementary Materials.
Statistical analysis
Crude BL incidence rates in 0-15 year olds were converted to units of cases per 100,000 person-years. Due to the rarity of the cancer, BL incidence is usually reported for aggregate time periods. In the analysis, rates were assumed to apply to the calendar midpoint of the data collection period.
Time trends in BL incidence
To generate pooled incidence rates for sub-Saharan Africa in the period before and after large-scale ITN introduction, we estimated the year of ITN introduction for each location with BL data. The Malaria Atlas Project (MAP) generates comparable model estimates of yearly ITN use since 2000 at a fine spatial resolution (12). Based on information about the covered population in each study, we linked each BL datapoint to its corresponding first administrative level (admin1) unit. The urbanicity status of the population covered by the study was used to split ITN use estimates into urban and rural areas within each unit using a threshold density of 1500 people per square kilometre. We calculated the mean ITN use weighted by the population at risk in each location over time.
ITN use has generally increased since 2000, but temporal patterns vary between locations (Figure S1). To capture the beginning of large-scale population-wide ITN campaigns, we defined the year of ITN introduction as the first year with at least 5% usage and an increasing trend. Several southern African countries had nationwide BL data available but did not introduce ITNs at a large scale (mean ITN use below 1%) (14) (Figure S1). These countries employed other malaria interventions (e.g. artemisininbased combination therapies and indoor residual spraying) (14); therefore, we considered their data separately.
BL incidence rates in each group (time period before and after ITN introduction, and in countries without large-scale ITN use) were analysed using a negative binomial regression model with person-time at risk as an offset and including location-level random intercepts. The pooled BL incidence rate with 95% confidence intervals was predicted for each group and for individual countries. Additionally, for the subset of locations with at least one BL datapoint before and after ITN introduction, we calculated the rate ratio for BL incidence after ITN introduction compared to before to assess trends in BL incidence for each location individually.
Population-level association between ITN use and BL incidence
To investigate the association between average ITN use in the population and BL incidence in children, we considered different time lags for population-level exposure to ITNs. The latency period between malaria exposure and BL onset is not fully understood: most children in highly endemic areas have their first exposure to both EBV and malaria in the first few years of life (15), while BL is most common in children aged 5-10 years (1). This observation led to the hypothesis that repeated malaria infections during childhood increase the risk of BL (5). Given this and the variable patterns in ITN scale-up (Figure S1), we hypothesised that ITN use in a particular year would be a poor predictor of BL incidence and instead used the mean population ITN use in the 10 years before BL data collection as measure of exposure. In sensitivity analyses, we also tested models with mean population ITN use in the 5 or 15 years before BL data collection, as well as concurrent ITN use. We assumed no exposure to ITN use in the population before the year 2000 (3).
We identified several a priori potential confounders for the population-level relationship between ITN use and BL incidence: admin1-specific baseline malaria prevalence in 2000 (before implementation of interventions) and concurrent HIV prevalence in adults, the country-specific human development index (HDI) and the urbanicity status of the population covered by the cancer registry or study (Table 1).
We calculated crude and adjusted rate ratios for these predictors using a negative binomial model with person-time at risk as an offset and including geographical location as a random intercept to account for potential correlation in datapoints collected in the same area at different timepoints (Supplementary Materials). We identified influential outliers in the regression model based on a Cook’s distance greater than 4/(total number of studies). In sensitivity analyses, we removed each influential outlier in turn from the model and restricted the dataset to locations with BL data before and after ITN introduction to evaluate the robustness of the findings.
All analyses were conducted in R statistical software (v4.2.2) using the glmmTMB (v1.1.7) and ggeffects (v1.2.3) packages. MAP estimates of ITN use were extracted from the foresite package (v0.1.0, (16)).
Results
Study characteristics
The literature search returned a total of 2333 articles and 116 records identified on IARC-associated websites (Figure 1). After screening and deduplication of the two sets of data sources, a total of 13 publications from literature databases and 10 IARC reports were included (1, 5, 17-37).
76 location- and time-specific BL incidence rates based on a total of 5,700 cases were extracted from the 23 publications. Included data covered the 1990-2017 time period, with slightly more datapoints available after 2005 (Figure 2A). 80% of datapoints came from population-based cancer registries. Among the 13 studies identified in the literature search, quality scores ranged from 1 to 4 (highest score; lowest risk of bias). Common limitations included passive case finding and inconsistent or poorly defined confirmation of diagnosis (Table S2).
There was a moderate positive correlation between the mean reported BL incidence rate and the number of datapoints in a country, but no relationship between the number of datapoints and the HDI (Figure 2B). BL incidence rates came from 32 locations in 21 countries, with 0-8 per country (Figure 2C). Many countries in Central and West Africa had no data on BL incidence, while countries in East and Southern Africa had better coverage. The actual catchment area covered in these countries was usually very small (Figure 2D). Only The Gambia in West Africa, and South Africa, Eswatini, Botswana and Namibia in Southern Africa had cancer registries with national-level coverage. However, these Southern African countries generally had the lowest BL rates. They were also considered separately, since analyses relating to ITN timing and use were restricted to the 66 other datapoints from locations with large-scale ITN use, covering 17 countries.
Time trends in BL incidence
The median year of ITN introduction in the included locations was 2005 (range 2004-2011). BL incidence was estimated to be 44% (95% CI 12-64%) lower in the time period after large-scale ITN introduction compared to before. Table 2 shows pooled incidence rates for each country for the period before (years 1990-2008) and after (years 2007-2017) introduction of ITNs, and across the whole 1990-2017 time period for countries without large-scale ITN use. Across countries, these pooled incidence rates were 1.36 (95% CI 0.88-2.10) and 0.76 (95% CI 0.50-1.16) per 100,000 person years for the periods before and after ITN introduction, respectively. The pooled rate was much lower in countries without ITNs [0.18 (95% CI 0.07-0.46) per 100,000 person-years]. However, there was also large variation in the estimated incidence rates between countries, ranging from 0.09 per 100,000 (95% CI 0.02-0.37) in Mali to 5.98 per 100,000 (95% CI 3.82-9.37) in Malawi before ITN introduction, and from 0.07 per 100,000 (95% CI 0.01-0.43) in Ethiopia to 8.52 per 100,000 (95% CI 2.75-26.33) in Malawi after ITN introduction. 65% of all BL cases in the dataset were recorded in three countries (Kenya, Uganda and Tanzania).
In locations with at least one datapoint before and after ITN introduction, rate ratios were highly variable and time trends were not always consistent across locations (Figure 3, Figure S2). The incidence rate appeared to be lower in the period after ITN introduction in 7 out of 12 locations, but confidence intervals overlapped in 3 of these. The largest reductions were seen in Maputo in Mozambique [92% (95% CI 69-98)], Brazzaville in Congo [88% (95% CI 67-96)] and in The Gambia [71% (95% CI 33-88)]. These locations also had among the highest mean ITN use since 2000 (over 75th percentile at 18%, 15% and 19%, respectively). In Bamako in Mali, the incidence rate in the 2005-2017 period was significantly higher than in 1988-1997. In the two registries without large-scale ITN use (Namibia and South Africa), rates were slightly higher after 2005 than before, but there was no significant difference between the two time periods.
Population-level association between ITN use and BL incidence
Table 3 shows the crude and adjusted associations of BL incidence rates with average population ITN use over 10 years and potential confounders. Univariate and multivariate showed comparable results. There was evidence for a population-level association between ITN use and BL incidence rates after adjusting for baseline malaria prevalence, adult HIV prevalence, country HDI level and population urbanicity (p=0.005). A 1% increase in average ITN use was associated with a 2% (95% CI 1-4%) reduction in BL incidence, as illustrated in Figure 4. Baseline malaria prevalence was positively associated with BL incidence after adjusting for all other covariates, but there was no evidence for an association with the other variables (Table 3).
Comparing the main result to different exposure time periods for ITN use, the association was similar for average ITN use in the 5 and 15 years before the BL datapoint [rate ratio 0.98 (95% CI 0.97-0.99) and 0.97 (95% CI 0.95-0.99), respectively]. In line with our hypothesis, concurrent ITN use in the population was not strongly associated with BL incidence (rate ratio 0.99, 95% CI 0.98-1.00) (Table S3). In additional sensitivity analyses, removing influential outlier studies and restricting the analysis to locations with BL incidence data before and after ITN introduction did not affect conclusions about the association between ITN use and BL incidence, though the strength of the association was lower in some cases (Table S4).
Discussion
In this study, we investigated changes in BL incidence in sub-Saharan African children in relation to large-scale introduction of ITNs in the early 2000s. In a systematic review, we identified 66 datapoints on BL incidence from locations with large-scale ITN use, comprising 5,226 cancer cases and covering 17 countries between 1990 and 2017. We found BL rates were 44% (95% CI 12-64%) lower in the time period after ITN introduction compared to before. Pooled incidence rates were estimated at 1.36 (95% CI 0.88-2.10) and 0.76 (95% CI 0.50-1.16) per 100,000 person-years before and after ITN introduction, respectively. In an ecological analysis, we described an association between increasing ITN use in the population and decreases in childhood BL incidence after adjusting for potential confounders [rate ratio 0.98 (95% CI 0.96-0.99) for average ITN use in the 10 years before BL data collection].
Previous studies and cancer registries have also noted a decreasing temporal trend in incidence, both for BL (5, 17, 38, 39) and childhood non-Hodgkin lymphoma more generally (39). Given existing evidence for an aetiological role of malaria infection in BL development, the population-level association between ITN use and BL described in this study is plausible. ITNs offer both individual and community-wide protection from malaria infection due to the protective physical barrier and the killing effect of insecticides on mosquitoes, thereby reducing malaria prevalence (40). As previously suggested for malaria infection, our results also point towards past cumulative rather than concurrent exposure to ITNs as a predictor for lower BL incidence on the population level (5, 8).
However, several alternative factors could explain changes in BL incidence over time. Firstly, studies on the effect of antiretroviral treatment (ART) coverage on incidence of HIV-related cancers suggest that HIV trends in children could explain declines in incidence (39, 41). However, concurrent HIV prevalence in adults was not associated with BL incidence in our study, and BL rates can also remain constant in people living with HIV in the ART era (42).
Secondly, the variable time trends seen in individual registries suggest other influences on observed BL incidence, which limit comparability of data from different time periods, locations and sources. The coverage, completeness and diagnostic methods of a registry can vary significantly over time. In Mali, BL incidence rates appeared to be higher in the period after ITN introduction, but the rate before ITN introduction was based on only three cases detected over a 10-year period. Underdiagnosis of childhood cancers was suspected at the time (29), with subsequent improvement in case finding procedures reflected in total cancer numbers (31). Studies of temporal trends in population-based cancer registries in Africa have also reported increases in incidence for other cancer types in adults (43-45). Even though access to diagnosis might generally be expected to improve over time and lead to higher rates, other factors such as intermittent political instability or funding disruption could also contribute to variation in rates (31, 44, 46).
Despite improvements, cancer surveillance in sub-Saharan Africa is still of low quality overall (46). While the historical data used in this study must be interpreted in the context of these limitations, improved cancer surveillance is essential for assessing current and future cancer burden and mortality. Our estimate of BL incidence in the period after ITN introduction can inform planning and targeting of resources in cancer control programmes. It provides further support for the rate reported in a previous study [0.64 per 100,000 person-years in 2018], which extrapolated national estimates of non-Hodgkin lymphoma incidence in countries without primary BL data (47).
Burden estimates are of particular importance because access to effective treatment for childhood cancer faces numerous challenges and remains a low priority for healthcare systems in sub-Saharan Africa (48). Less than a third of children diagnosed with cancer are estimated to survive compared to 80% in high-income countries (49). The possibility that BL incidence may be decreasing in sub-Saharan Africa following successful malaria control programmes therefore has positive implications in terms of a reduced burden on patients and the health system. Given the comparatively higher political commitment and funding for infectious disease programmes like malaria and HIV, a multidisciplinary approach to addressing childhood cancer and infectious diseases has previously been suggested (49, 50). Collaboration with existing malaria programmes could raise awareness, attract funding and increase political commitment to BL prevention, treatment and research (50). Nevertheless, with plateauing funding and rising insecticide and treatment resistance, continued progress in malaria control also needs to remain a priority.
Our study used all publicly available data accrued over almost 30 years and, to the best of our knowledge, provides the first evidence for a population-level decline in childhood Burkitt lymphoma burden in malaria-endemic sub-Saharan Africa associated with large-scale rollout of ITNs for malaria control. These findings provide support for the existing observational evidence base for a potential role of malaria interventions in preventing BL in a context where conducting a randomised controlled trial would be unfeasible and unethical.
Beyond the considerations already outlined, there were several additional limitations. Firstly, BL incidence rates were assembled from various locations, and only a subset of these had data at more than one timepoint. It is not known how representative this published data would be of all incidence rates across malaria-endemic sub-Saharan Africa, and unknown systematic differences between the locations included in the period before ITN introduction compared to after could have biased the results. Secondly, our inclusion criteria were relatively broad to make use of all available data, but different methods in recording of cancer incidence between locations and over time could have affected the results, especially as confirmation of diagnosis was frequently suboptimal (Table S2). The small case numbers in some studies also suggest substantial underdiagnosis. Lastly, since an ecological study design was employed, the association between ITN use and BL should not be interpreted to apply on the individual level; both ITN use and potential confounders only represent proxies for the average exposure in the population, and were in many cases based on model estimates.
In conclusion, this study suggests that childhood BL burden in sub-Saharan Africa may be declining as ITN use in the population increases. Analysis of all available data on BL incidence in malaria-endemic African countries indicates that incidence is lower now than before the introduction of mass ITN distribution programmes around 2005. However, data came from various sources and was often of low quality, underlining the need for improved cancer surveillance in sub-Saharan Africa to allow robust investigations of trends in cancer incidence.
Contributors
NS conceived the study. NS and JK conducted the systematic review and were responsible for data curation. EKA provided domain knowledge for the study and contributed to study design. NS conducted the analysis and wrote the first draft of the manuscript. All authors contributed to interpreting the results, reviewing and revising the manuscript, and approved the final version for submission.
Competing interests
We declare no competing interests.
Data availability statement
All data reported in this study were extracted from published studies. All covariate estimates used in the analysis are publicly available from the sources cited in the article.
Funding
NS acknowledges funding from a Wellcome Trust grant awarded to Prof Azra Ghani and Prof Katharina Hauck [reference 220900/Z/20/Z] and from the MRC Centre for Global Infectious Disease Analysis [reference MR/R015600/1], jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth & Development Office (FCDO), under the MRC/FCDO Concordat agreement and is also part of the EDCTP2 programme supported by the European Union. EKA is supported by a Cancer Research UK (CRUK) grant awarded to Dr Tsilidis K [reference PPRCPJT\100005] and acknowledges funding from a return grant from the International Agency for Research on Cancer-World Health Organization (IARC-WHO). For the purpose of open access, the authors have applied a ‘Creative Commons Attribution’ (CC BY) license to any Author Accepted Manuscript version arising from this submission.