On-going malaria transmission in The Gambia despite high coverage of control interventions: a nationwide cross-sectional survey

As indicators of burden of malaria have substantially decreased in The Gambia, reaching a pre-elimination status may be attainable. Achieving this goal requires in-depth understanding of the current burden of Plasmodium falciparum infection. A nationwide cross-sectional survey was conducted in 2012 to determine the prevalence of P.falciparum infection, and to describe its heterogeneity and associated risk factors. Finger-prick blood samples were collected for microscopy, species-specific PCR and haemoglobin measurement. A total of 9,094 participants were included and median age was 11.9 years (IQR 5, 28). Overall prevalence of P. falciparum was 16.01 % with marked heterogeneity between sites (4.32–36.75 %) and within villages in each site (1.63–49.13 %). Across all sites, 51.17 % (745/1,456) of infections were asymptomatic and 35.61 % (448/1,258) were sub-microscopic. The odds of P. falciparum infection were higher in older children; 5–15 years (OR = 1.90; 95 % CI 1.60–2.26), adults (OR = 1.48; 95 % CI 1.24–1.78) and participants with moderate anaemia (OR = 1.62; 95 % CI 1.32–1.99). The current malaria control interventions are not sufficient to interrupt transmission in The Gambia as malaria prevalence is still relatively high in the eastern part of the country. New interventions aiming at interrupting transmission are needed and should be urgently evaluated.


Background
The last decade has witnessed major progress in the fight against malaria as the burden has substantially reduced, even in sub-Saharan Africa, although malaria still contributes significantly to both morbidity and mortality [1,2]. The Gambia is one of the African countries where a substantial decline has been observed [3,4]. This decline resulted from the scaling-up of malaria control intervention [5], which includes: increased availability and access to long-lasting insecticide bed nets (LLINs), integrated vector control interventions such as indoor residual spraying (IRS), strengthened case management with rapid diagnostic tests (RDTs), and artemisinin combination therapy (ACT). Despite the scaling-up of these interventions, malaria transmission, which is highly seasonal, has not been interrupted. Earlier studies in schoolchildren [6,7] showed marked heterogeneity in malaria prevalence across the country and significant seasonal variation [8]. In addition, high resistance to DDT and pyrethroids observed in some areas raises concerns about the impact of current vector control interventions [9,10].

Open Access
*Correspondence: udalessandro@mrc.gm 1 Medical Research Council Unit, PO Box 273, Banjul, The Gambia Full list of author information is available at the end of the article As the country aims for pre-elimination status, gaining a detailed understanding of the current distribution of malaria infection, characterizing its heterogeneity across the country and further identifying asymptomatic carriers is critical as these are important elements for disease surveillance and targeting control efforts. Whereas health facility data have previously been used to describe malaria morbidity trends [3], this information is greatly limited by incompleteness, differing malaria case definitions and the lack of data on asymptomatic carrier status. This cross-sectional survey was, therefore, conducted to better characterize the dynamics of malaria transmission in The Gambia and to identify the determinants of its heterogeneity.

Methods
The Gambia is divided into five administrative regions, namely West Coast (WC), North bank region (NBR), Lower River (LRR), Central River (CRR) and Upper River (URR) regions ( Fig. 1). It is characterized by a long dry season from mid-October to mid-June followed by a single short rainy season from June to September when the average daily temperature is at 28 °C, creating a favourable environment for Anopheles gambiae mosquitoes. Malaria transmission occurs almost exclusively during the rainy season and immediately afterwards, until December-January. The entomological inoculation rate varies from 0 to 166 infective bites per person per year [11]. For this study, the country's five geographical regions were used; WC, NBR, LRR-south bank, CRRnorth bank and subdivided the URR into the north and south bank. In each of these regions, one primary school with the highest seroprevalence of anti-malarial antibodies (MSP1 19 ) among school children was identified using findings from a previous 2012 dry season nationwide seroprevalence survey [7]. In November 2012, six neighbouring villages around each of these schools, with populations between 100 and 500 were selected for inclusion in this survey. In villages with less than 300 inhabitants, the whole population was included while in larger villages the first 350 participants were consecutively enrolled. In all the villages participants were requested to gather at a central location for the survey, usually the village central meeting area. Infants less than 6 months of age or individuals who had stayed less than 4 weeks in the village were excluded. Village locations were mapped using a hand-held global positioning system (GPS) (Garmin eTrex ® 10).
Each participant was allocated a unique study ID code and individual interviews were conducted for older children (12-17 years) and adults, while parents/guardians completed the questionnaire for younger children (6 months-11 years). Individual participant information on demographic and general health indicators such as history of fever, sleeping habits (including outdoor or out in the field) and travel outside the village in the past month and IRS were collected using a structured questionnaire. Individual participants were asked if their house had received IRS in the last 1 year and in addition if they had slept under a LLIN the previous night. Axillary temperature was measured using a digital thermometer. A finger-prick blood sample was collected from each participant on a blood slide for microscopy and spotted on filter paper (Whatman 3 Corporation, Florham Park, NJ, USA) for molecular analysis. Haemoglobin (Hb) was measured with the HaemoCue 301 machine (Ängelholm, Sweden) following the manufacturer's instructions.

Laboratory procedures
To determine parasite density, microscopy was read on all samples that tested positive by PCR. Slides were stained with 2.5 % buffered Giemsa (pH 7.2) for 30 min and double read by trained microscopists who were blinded to the survey data of the participants. Slide smears were declared negative if no parasite was seen after examining 200 high power fields (HPF). If positive, parasites were counted against 500 white blood cells (WBCs and parasite densities were calculated assuming 8,000 WBC per μl of blood. A 20 % error check was used to identify discrepancies between slide readers. All discordant results were read by a senior microscopist who was also blinded to participant survey data and these results were taken as the final read. For diagnostic PCR, three 6-mm dried filter-paper blood spots (DBS) were punched into a 96-well plate. DNA extraction was carried out using the automated QIAxtractor robot (Qiagen). Negative and positive (3D7) controls were included to control for cross contamination and DNA extraction efficiency, respectively. The DBS were lysed by incubating them in tissue digest buffer at 60 °C for 1 h and digested eluates were applied onto capture plates, washed, and the DNA eluted into 80 µl. The extracted DNA (4 µl) was used in a nested PCR, amplifying the multicopy Plasmodium ribosomal RNA gene sequences using genus and species specific primers [12]. All PCR products were run using the QIAxcel capillary electrophoresis system (Qiagen), using the screening cartridge and 15-1,000 bp-alignment marker. Results were exported and double scored using both the QIAxcel binary scoring function and manually by visualization of the gel images and discrepancies were scored by a third independent reader. All readers were blinded to participant survey data.

Statistics
The sample size was based on the desired precision of the prevalence of infection by village; 300 individuals allowed an estimation of the prevalence with the following precision: 1 % (95 % CI 0.2-3 %); 5 % (3-8 %); 12 % (8-16 %), and 20 % (16-25 %). Data from questionnaires were verified for completeness and all data were double entered and verified using Microsoft Access database (Microsoft Corp., Redmond, USA). Errors in data entry were detected by running consistency checks before exportation to STATA 13.0 (StataCorp, College Station, Texas, USA) for analysis. Age groups were categorized into <5 years, ≥5 to ≤15 years and >15 years [13,14]; anaemia was defined according to the WHO Haemoglobin concentration for diagnosis of anaemia and assessment of severity; as severe (<8.0 g/dl), moderate (8.0-10.9 g/ dl), or mild (11.0-12.9 g/dl) and no anaemia (≥12.9 g/ dl) [15]. Carriage of Plasmodium falciparum infection was defined as parasitaemia detected by PCR at screening [16]. LLIN coverage was reported as the proportion of the population that slept under an LLIN the previous night.
Descriptive statistics are presented for continuous variables (median and interquartile range), and proportions for categorical variables and all point estimates presented with 95 % confidence intervals calculated using the Wilson score method with continuity correction. Uni-and multi-variable logistic regression analyses were performed to determine the independent predictors of P. falciparum infection. To account for missing data, complete case analyses were presented along with multiple imputation models (50 iterations), which allowed for the inclusion of partially observed cases. For multiple imputation analysis, a fully conditional model specification [17,18] was adopted in which a logistic model was used for binary variables (gender, slept under a LLIN the previous night, sleep outdoors at night, IRS, sleep in the fields, and used artemether-lumefantrine the previous week) and an ordinal logistic model was performed for age and anaemia categorical variables. There was an allowance for clustering between and within the six sites in the uni-multivariable models. For age and anaemia, categorical variables, tests for trend and significance of the associations were performed in the multivariable models and, a goodness of fit test performed for the final model. For the analysis of microscopy, the median asexual and sexual parasite densities were calculated and the proportions of parasitaemia determined. The geolocations and P. falciparum prevalence data of each village were imported from Microsoft access and transferred separately onto a layer containing the base map of The Gambia generated using Quantum GIS 2.0.1 (QGIS 2.0.1). The village GPS points with the corresponding variable attributes were then displayed on the new layer and the maps were exported as jpeg, 800dpi file format ( Fig. 1).
This study was approved by the Gambia Government/ MRC Joint Ethics Committee (SCC1318). Verbal consent for participation was first obtained during the village sensitization meetings. Written informed consent was obtained for all participants, with parents/guardians providing written consent for children less than 12 years. Assent was also obtained from children aged between 12 and 17 years.

Results
Over 3 weeks in November 2012, a total of 10,191 participants were enrolled from 36 villages. Of these, 9,094 (89.24 %) participants with PCR results and complete questionnaire data were included in this analysis. The median age was 11.9 years (IQR 5, 28) and the majority [59.57 % (5,417/9,094); 95 % CI 58.55-60.58 %] were females ( Table 1). The mean (SD) haemoglobin was 12.7 g/dl (2.26)  In the multivariable analyses, with clustering with-in and between sites the odds of P. falciparum infection were significantly higher among older children (OR = 1.90; 95 % CI 1.60-2.26; P ≤ 0.001) and adults (≥15 years) (OR = 1.48; 95 % CI 1.24-1.78; P ≤ 0.001) compared to younger children (<5 years). In addition, females had a significantly lower odds of parasitaemia (OR = 0.87; 95 % CI 0.76-1.00; P = 0.05) compared to males. Plasmodium falciparum infection was also strongly associated with anaemia, with the odds of parasitaemia increasing with severity of anaemia. The odds of infection were three-fold higher among participants with severe anaemia (moderate: OR = 1.62; 95 % CI 1.32-1.99; P ≤ 0.001; severe: OR = 3.11; 95 % CI 2.27-4.26; P ≤ 0.001) compared to participants with no anaemia (test for trend anaemia P = 0.002). Sleeping under a LLIN the previous night was protective against malaria infection (OR = 0.58; 95 % CI 0.41-0.82; P ≤ 0.001). In NBR, sleeping outdoors at night was associated with a higher odds of P. falciparum infection (OR = 1.98, 95 % CI 1.16-3.41; P = 0.01). However, the association was not statistically significant in the overall uni-multivariable analyses. "Sleeping outdoors at night", "indoor residual spraying" and "sleeping out in the fields" were not predictors of P. falciparum infection. There were no differences in the covariates predicting malaria infection when compared to the multilevel, multivariable, parsimonious model, suggesting missing data did not affect the observed trends in association (Table 4).

Discussion
Whereas the burden of malaria in The Gambia has decreased significantly over the last 10 years [3,4], probably because of the moderately high coverage of preventive interventions (LLINs and IRS) and free access to ACT, malaria transmission is ongoing. The prevalence of P. falciparum infection in villages in the NBR (site 2) was similar to that reported from community-based surveys in 2009 [8], implying that despite a substantial reduction, malaria transmission is still ongoing. The observed prevalence of infection close to 50 % in some villages in the eastern part of the country is comparable to the prevalence reported from high transmission areas in Tanzania [13], Gabon [19] or Burkina Faso [20]. In the rest of the country, prevalence was much lower and heterogeneous, with a five-to eight-fold difference between neighbouring villages. Although the selection of villages for this survey was done according to documented high prevalence of anti-malarial antibodies in school children, the results reported are likely to be a fair representative of the malaria transmission in The Gambia as geographical representation was taken into account. The difference in transmission between the eastern part and the rest of the country cannot be explained by different coverage of LLINs, as their reported use is relatively high throughout the country and LLIN use reduced odds infection of P. falciparum infection. However, the high resistance to DDT and pyrethroids observed in eastern Gambia [9,10] and in Senegal [21] could be responsible for the higher malaria prevalence. Other factors to consider are increased outdoor biting of the local vectors in response to environmental changes as has been reported previously [22,23], the differences in vector composition across the regions could also contribute to the higher prevalence of P. falciparum in the eastern region. A. gambiae s.s and A. arabiensis are the predominant and in general most efficient vectors in the east, while other areas of the country have less efficient vectors, such as A. melas and A. coluzzii. In addition, this region has more difficult access to health care and prompt treatment. For the latter, the eastern part of The Gambia, where the highest malaria prevalence has been observed, is the least developed and the farthest from the capital.

West Coast North bank region LRR-south bank CRR-north bank URR-south bank URR-north bank
The low malaria prevalence observed in the western part of the country is consistent with earlier reports in 2008 [8] indicating that 4 years later and despite high coverage of preventive interventions, the human reservoir of infection has not decreased substantially. A significant proportion of infections in this region were also sub-microscopic, implying that they would be missed by routinely used diagnostic tools and would therefore go untreated, hence contributing to maintain the human reservoir of infection and transmission. Conversely in the URR-south bank and NBR, whereas the majority of infections were detected by microscopy, a significant proportion of the participants were asymptomatic and thus unlikely to attend health facilities. Such individuals could therefore remain parasitaemic in the community and maintain residual transmission as well. Identifying individuals who constitute this reservoir of infection would be extremely important for the identification of 'hot spots' in these villages, and designing and targeting interventions aimed at further decreasing malaria transmission. Assuming such carriers are mostly the same individuals from one transmission season to the other, because of environmental factors, such as higher exposure to infectious bites, residing in households near breeding sites, or genetic factors [24], targeting them for treatment could interrupt transmission. An intense follow-up of some villages included in this survey and covering several transmission seasons is currently ongoing and should be able to prove or disprove this hypothesis.
As previously reported [25,26], about 36 % of PCRpositive samples were negative by microscopy and thus had sub-microscopic infections. The role of such infections in maintaining transmission remains unclear and urgently needs to be addressed [27]. Sub-microscopic infections may be important in maintaining transmission; an increasing body of evidence indicates that a large proportion of gametocyte carriers remain undetected by microscopy, and gametocytes may actually be present in the majority of infections, albeit at low densities [28]. When considering that the association between gametocyte density and mosquito infection rate is not very strong and is most variable at low gametocyte concentrations [28], the contribution of the carriers identified within this survey to the maintenance of malaria transmission is plausible. The challenge is to identify and treat them in the most efficient way. Active case detection strategies are increasingly popular across a number of countries although there is little evidence of their effectiveness [29]. Their main limitation is the sensitivity of the diagnostic screening tools used, which would miss an important proportion of sub-microscopic carriers unless field-adapted molecular diagnostic methods are used [30,31]. Other chemoprevention approaches, such as seasonal chemoprevention in children less than 5 years of age recently implemented in some regions of The Gambia and mass drug administration aimed at eliminating the reservoir of infection should be considered as additional interventions, especially in areas where transmission has not declined significantly.
The prevalence of asymptomatic parasitaemia was significantly higher in the older children and adults, possibly the result of stronger immunity, making them able to tolerate [32] but not eliminate a malaria infection. Therefore, infections acquired among older children and adults during the current transmission season may remain asymptomatic and untreated, resulting in the higher prevalence observed in these age groups. Nevertheless, such infections are not without consequences as infected individuals had a significantly higher risk of anaemia. This suggests that even if not associated with acute symptoms, malaria infections do have a deleterious effect on the health of infected individuals. A similar association was described in Ghana [33], Cameroon [34] and Tanzania [13] but not among Ugandan schoolchildren [35].

Conclusion
When considering the distribution of malaria infection across The Gambia, two main strata can be identified: the eastern part where transmission is still relatively intense, with villages in which a third to a half of the population infected, and the western part with low prevalence. Within each stratum, there is considerable heterogeneity between villages. It is remarkable to notice that even in the low endemicity stratum and despite high coverage of preventive and curative interventions, transmission has been maintained over the years. Obviously, current interventions are not sufficient to interrupt transmission in the current context of modified human and vector behaviour and new approaches need to be urgently evaluated. The low endemicity stratum in The Gambia offers an ideal setting to test new interventions aiming at interrupting malaria transmission.