Assessing the benefits of five years of different approaches to treatment of urogenital schistosomiasis: A SCORE project in Northern Mozambique

Background In Mozambique, schistosomiasis is highly endemic across the whole country. The Schistosomiasis Consortium for Operational Research and Evaluation (SCORE) coordinates a five-year study that has been implemented in various African countries, including Mozambique. The overall goal of SCORE was to better understand how to best apply preventive chemotherapy with praziquantel (PZQ) for schistosomiasis control by evaluating the impact of alternative treatment approaches. Methods This was a cluster-randomised trial that compared the impact of different treatment strategies in study areas with prevalence among school children of ≥21% S. haematobium infection by urine dipstick. Each village was randomly allocated to one of six possible combinations of community-wide treatment (CWT), school-based treatment (SBT), and/or drug holidays over a period of four years, followed by final data collection in the fifth year. The most intense intervention arm involved four years of CWT, while the least intensive arm involved two years of SBT followed by two consecutive years of PZQ holiday. Each study arm included 25 villages randomly assigned to one of the six treatment arms. The primary outcome of interest was change in prevalence and intensity of S. haematobium among 100 children aged 9-to-12-years that were sampled each year in every village. In addition to children aged 9-to-12 years, 100 children aged 5–8 years in their first-year of school and 50 adults (aged 20–55 years) were tested in the first and final fifth year of the study. Prevalence and intensity of S. haematobium infection was evaluated by two filtrations, each of 10mL, from a single urine specimen. Principal findings In total, data was collected from 81,167 individuals across 149 villages in ten districts of Cabo Delgado province, Northern Mozambique. Overall PZQ treatment resulted in a significant reduction in the prevalence of S. haematobium infection from Year 1 to Year 5, where the average prevalence went from 60.5% to 38.8%, across all age groups and treatment arms. The proportion of those heavily infected also reduced from 17.6% to 11.9% over five years. There was a significantly higher likelihood of males being infected than females at baseline, but no significant difference between the sexes in their response to treatment. The only significant response based on a study arm was seen in both the 9-to-12-year-old and first-year cross sections, where two consecutive treatment holidays resulted in a significantly higher final prevalence of S. haematobium than no treatment holidays. When the arms were grouped together, four rounds of treatment (regardless of whether it was CWT or SBT), however, did result in a significantly greater reduction in S. haematobium prevalence than two rounds of treatment (i.e. with two intermittent or consecutive holiday years) over a five-year period. Conclusions Although PC was successful in reducing the burden of active infection, even among those heavily infected, annual CWT did not have a significantly greater impact on disease prevalence or intensity than less intense treatment arms. This may be due to extremely high starting prevalence and intensity in the study area, with frequent exposure to reinfection, or related to challenges in achieving high treatment coverage More frequent treatment had a greater impact on prevalence and intensity of infection when arms were grouped by number of treatments, however, cost efficiency was greater in arms only receiving two treatments. Finally, a significant reduction in prevalence of S. haematobium was seen in adults even in the SBT arms implying the rate of transmission in the community had been decreased, even where only school children have been treated, which has significant logistical and cost-saving implications for a national control programme in justifying CWT.


Background
In Mozambique, schistosomiasis is highly endemic across the whole country. The Schistosomiasis Consortium for Operational Research and Evaluation (SCORE) coordinates a five-year study that has been implemented in various African countries, including Mozambique. The overall goal of SCORE was to better understand how to best apply preventive chemotherapy with praziquantel (PZQ) for schistosomiasis control by evaluating the impact of alternative treatment approaches.

Methods
This was a cluster-randomised trial that compared the impact of different treatment strategies in study areas with prevalence among school children of !21% S. haematobium infection by urine dipstick. Each village was randomly allocated to one of six possible combinations of community-wide treatment (CWT), school-based treatment (SBT), and/or drug holidays over a period of four years, followed by final data collection in the fifth year. The most intense intervention arm involved four years of CWT, while the least intensive arm involved two years of SBT followed by two consecutive years of PZQ holiday. Each study arm included 25 villages randomly assigned to one of the six treatment arms. The primary outcome of interest was change in prevalence and intensity of S. haematobium among 100 children aged 9-to-12-years that were sampled each year in every village. In addition to children aged 9-to-12 years, 100 children aged 5-8 years in their first-year of school and 50 adults (aged 20-55 years) were tested in the first and final fifth year of the study. Prevalence and intensity of S. haematobium infection was evaluated by two filtrations, each of 10mL, from a single urine specimen.

Principal findings
In total, data was collected from 81,167 individuals across 149 villages in ten districts of Cabo Delgado province, Northern Mozambique. Overall PZQ treatment resulted in a significant reduction in the prevalence of S. haematobium infection from Year 1 to Year 5, where the average prevalence went from 60.5% to 38.8%, across all age groups and treatment arms. The proportion of those heavily infected also reduced from 17.6% to 11.9% over five years. There was a significantly higher likelihood of males being infected than females at baseline, but no significant difference between the sexes in their response to treatment. The only significant response based on a study arm was seen in both the 9-to-12-year-old and first-year cross sections, where two consecutive treatment holidays resulted in a significantly higher final prevalence of S. haematobium than no treatment holidays. When the arms were grouped together, four rounds of treatment (regardless of whether it was CWT or SBT), however, did result in a significantly greater reduction in S. haematobium prevalence than two rounds of treatment (i.e. with two intermittent or consecutive holiday years) over a five-year period.

Conclusions
Although PC was successful in reducing the burden of active infection, even among those heavily infected, annual CWT did not have a significantly greater impact on disease prevalence or intensity than less intense treatment arms. This may be due to extremely high starting prevalence and intensity in the study area, with frequent exposure to reinfection, or related to challenges in achieving high treatment coverage More frequent treatment had a greater impact on prevalence and intensity of infection when arms were grouped by number of treatments, however, cost efficiency was greater in arms only receiving two treatments. Finally, a significant reduction in prevalence of S. haematobium was seen in adults even in the SBT arms implying the rate of transmission in the community had been decreased, even where only school children have been treated, which has significant logistical and cost-saving implications for a national control programme in justifying CWT.

Author summary
Urogenital schistosomiasis is highly endemic in Mozambique. This study was part of a multi-country trial, including Mozambique, designed to understand the impact of different schistosomiasis treatment strategies involving community-wide treatment (CWT), school-based (SBT), and treatment holidays over a five-year period. Results from Mozambique showed that although preventive chemotherapy was successful in reducing the prevalence of Schistosoma haematobium over five-years, the most intense treatment approach, annual CWT, did not have a significantly greater impact than less intense treatment strategies, such as bi-annual SBT. Infection rates were higher among males, but there was no difference in response to treatment by gender. Four rounds of treatment (regardless of whether it was given in the community or school) did result in a significantly greater

Introduction
Schistosomiasis is a major yet neglected public health problem, second only to malaria in terms of parasite-induced human morbidity and mortality worldwide [1]. Estimates show that globally at least 218 million people required preventive treatment in 2015 [2][3][4], and at least 20 million suffer from severe and debilitating forms of the disease [5]. Schistosomiasis is a major public health problem in Mozambique, as shown by an epidemiological survey of schistosomiasis and soil-transmitted helminthiasis among school children carried out between 2005 and 2007 [6]. The mean estimated prevalence of urogenital schistosomiasis, Schistosoma haematobium, was 47% while that of intestinal schistosomiasis, S. mansoni, was much lower (around 1%) across all of Mozambique. In Cabo Delgado province, the area where this study took place, the prevalence of S. haematobium was 57.9%, ranging from 8.8% on the coast to 93% inland [6].
In 2001, the World Health Organisation (WHO) endorsed preventive chemotherapy (PC) as the global strategy to control morbidity due to schistosomiasis through regular treatment with praziquantel (PZQ) [7]. Treating all school-aged children (SAC) is more cost-effective than a test-and-treat approach [8,9], which is achieved by allocating a geographic area (typically the district) to a recommended WHO treatment category for schistosomiasis based on infection prevalence [7]. As a result, most programmes control morbidity associated with schistosomiasis through school-based deworming, as it is highly cost-effective [10]. However, the limitation of this approach is poor coverage among out-of-school children, particularly older children and those living far from schools, and neglecting children younger than four years and adults [11]. Studies have highlighted the advantages of a community-wide treatment (CWT) approach in addition to a school-based treatment (SBT) in reducing prevalence and intensity of schistosomiasis infection, particularly in communities where not all children attend school [11][12][13][14][15]. Furthermore, questions remain about optimal frequency of PZQ treatment for infection and morbidity control. It has been observed in some studies that more frequent dosing can improve reductions in worm burden and increase parasitological cure, particularly where schistosomiasis is highly endemic [16]. Whereas others have shown that a single dose of PZQ results in sustained low transmission of S. haematobium for two years [17].
The Schistosomiasis Consortium for Operational Research and Evaluation (SCORE) was established in 2008 to answer strategic questions about schistosomiasis control [18]. It includes multi-country field studies that aim to understand the benefits and costs of alternative approaches to PC involving CWT, SBT and "drug holidays" (i.e. years without PC). The gaining and sustaining SCORE study protocol and baseline characteristics have been described elsewhere [18][19][20][21]. The primary research question presented here is which strategy for PC provided the greatest reduction in prevalence and intensity of S. haematobium infection among 9-to-12-year olds after four years of intervention in Cabo Delgado, Northern Mozambique. In addition, the impact of treatment on first-year students and adults sampled in each village was also assessed. These findings will provide an evidence base for the Mozambique national control programme for schistosomiasis to address strategic questions about schistosomiasis treatment and potentially shift from morbidity control to interruption of transmission, and subsequently elimination.

Study arms
This study was a parallel cluster-randomised, intervention trial with six study arms (Fig 1). Communities received various combinations of CWT, SBT or drug holiday over a four-year period, with the final round of data collection carried out 12 months after Year 4 treatment. All communities received CWT directly after Year 5 data collection, as they were no longer participating in the trial. The most intensive intervention arm involves four years of CWT (Arm 1), whilst the least intensive treatment strategy was two years of SBT PC followed by two consecutive years of drug holidays (Arm 5).

Study area and population
The area for this study was 10 out of 17 districts of Mozambique's northern most province, Cabo Delgado. The area was chosen as previous mapping had shown a high prevalence of S. haematobium and a low prevalence of S. mansoni, as one of the country selection criteria was no mixed infections to ensure that species type did not confound response to a particular treatment strategy [18].
Communities to be included in the study were determined by convenience sampling based on three criteria: (i) selecting only villages that had a primary school (so that it could be randomised to SBT), (ii) village and school had no history of PC using PZQ against schistosomiasis, (iii) and the village has a school attended by a minimum of 100 children aged 9-to-12-years. Given that children who test positive must be treated, the eligibility survey was carried out among fifty 13-to-14-year olds in each village, therefore treating those infected did not affect subsequent study results, especially where prevalence was high [18]. The starting prevalence of S. haematobium was evaluated by reagent strip testing for microhaematuria on a single midday urine (communities were eligible only if !21% by urine dipstick was detected) [22]. A total of 150 communities found eligible were randomly assigned using a computer-based randomisation procedure, without stratification, to one of the six study arms. Since the rainy season lasts from November to April, data were collected between July and October each year to omit seasonal variation. A door-to-door census was conducted prior to treatment in Year 1 in all selected communities to provide a baseline for estimating coverage.

Sample size
The aim was to enrol 25 villages per study arm, and to monitor each village's prevalence and intensity of S. haematobium infection among 100 school children aged 9-to-12-years. In each school, 50 boys and 50 girls were selected using systematic random sampling from those in the second, third and fourth classes. A total of 150 villages were enrolled, with 15,000 children tested each year. Parasitological data were collected prior to treatment from Year 1 through to Year 5. To evaluate the effect of treatment holidays, currently recommended by WHO PC guidelines whereby moderate-risk communities are treated once every two years [7], no testing of children was conducted during holiday years as children who were found to be infected would have needed to be treated. Sample size calculations assumed treatment interventions would reduce schistosomiasis prevalence in high endemicity areas from 50% in Year 1 to 15% by the end of the study in the most intense treatment (Arm 1). Power computations used generalized estimating equations to fit a logistic regression model that included treatment arm and time effects and treatment-by-time interaction, and assumed an over-dispersion parameter of φ = 5.0. The minimum effect size was computed so that a difference over time could be detected with 90% power for a 2-sided α = 0.05 level test. Negligible correlation was assumed between variables in Year 1 and at study conclusion. Based on these assumptions, the calculations estimated that studying 25 villages per arm and 100 children per village were sufficient to detect an overall difference between arms of 11.4% prevalence at the conclusion of the study. In addition to children aged 9-to-12-years, systematic random sampling was used in each village to select 100 first-year students (aged 5-to-8-years old) and a convenience sample of 50 adults (aged 20-to-55-years) in the first and fifth years only.

Urine examination
In each school, selected children were given a plastic pot labelled with an adhesive barcode with unique identification numbers and asked to provide a single urine specimen. The sex, age, and school class for each child were recorded. All samples were collected between 10am and 2pm and examined immediately in the school. Two filtrations of 10mL each were carried out on a single urine sample using a nylon filter with a pore size of 20µm (Sefar AG; Heiden, Switzerland) [23,24]. Both filters were examined for S. haematobium eggs under a light microscope and number of S. haematobium eggs were counted. All parasitological examinations were performed by trained laboratory technicians. For quality control, around 10% of all microscope slides were re-examined by a senior technician. The intensity of infection was expressed as the number of eggs per 10mL of urine filtered. For specimens of less than 10 mL, the volume of urine filtered was measured and the number of eggs per 10mL calculated. If estimated counts were above 1,000 eggs per 10mL, they were truncated at 1,000. The arithmetic mean of two filtrations taken from a single urine specimen was calculated to be the egg count of the child [25]. A child was deemed egg positive if one or more eggs were found in any of the slides examined.

Treatment
Praziquantel tablets were delivered at approximately 40mg/kg using a WHO dose pole [26]. All treatment was carried out in the village on the same day as sampling whereby everyone was treated regardless of their parasitological results, which were anonymous. During SBT, directly-observed treatment with PZQ was administered by trained teachers to all children attending school. Efforts were made to treat non-school attendees through community sensitisation and mobilisation efforts. CWT involved providing treatment to the entire eligible population using community drug distributors either door-to-door or at fixed points in the study community, which only excluded children under 4 years of age or under 94 cm in height. Everyone who had taken the medication in SBT was monitored for adverse events for four hours after treatment, since the analysis was carried out in the school. In the CWT, since the drug distributor lives in the village they were available in the community.
At baseline (Year 1) a census was carried out in all selected communities of the total population and proportion of all SAC (aged 5-14-years) recorded in the population to provide a baseline for estimating coverage.

Ethical statement
Informed consent was obtained from all individuals !18 years of age and from parents or legal guardians of children less than 18 years of age. The purpose of the study was explained to all schoolchildren and verbal assent was obtained from the children. Permission was also obtained from school headmasters. Ethical clearance was obtained from the National Bio-ethical Committee for Health of Mozambique (NBCHM), and the survey was conducted according to NBCHM guidelines (reference no. IRB00002657). The trial is registered with the International Standard Randomised Controlled Trial registry under ISRTC number 14117624 for Mozambique. The study protocol was also approved by Imperial College London (ICREC_10_2_2).

Data handling and analysis
Demographic data were collected on smartphones and uploaded to a dedicated database maintained on a central server (EpiCollect at Imperial College London). Laboratory data were collected on paper forms and entered onto the smartphone retrospectively using the barcode with unique identification number to link with the individual's demographic information. The electronic data capture system involved entering data in the field, synching the data to a central server and downloading the data retrospectively for cleaning and analysis. Data cleaning and management was carried out by biostatisticians at the Schistosomiasis Control Initiative (SCI). Data analysis was performed using R version 3.2 (R Core Team (2013): R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria).
Since 9-to-12-year olds are the only age group to be consistently sampled from baseline to Year 5, the primary research question presented here is which strategy for PC provided the greatest reduction in prevalence and intensity of S. haematobium infection among 9-to-12-year olds after four years of intervention. In addition, the impact of treatment on first-year students and adults have also been reviewed and analysed.
Individuals without data on age, sex, and presence or absence of eggs on at least one slide were not included in the study. For each year of data, three infection categories for S. haematobium infection were created: no infection, light infection (<50 eggs/10mL of urine), and heavy infection (!50 eggs/10mL of urine) [7]. Indicator variables were created for infection status (0 = not infected, 1 = infected) and heavy infection status (0 = egg count <50, 1 = egg count> 50). These were averaged across each village to produce a village level prevalence and village level prevalence of heavy infected, respectively. 95% confidence intervals were calculated for village and study arm level prevalence values. Survey package [27] was used to consider survey design effects. Prevalence of infection and heavy infection was calculated as arithmetic means of the infection categories, aggregated by the relevant factors (e.g. age, sex). The arithmetic mean of infection intensity was calculated using both egg positive and negative, aggregated by grouping factors (village, study arm and gender). Although, mean intensity was underestimated as egg counts were capped at 1000 eggs/10mL, this was done consistently over the years therefore the trend in intensity change is believed to be valid. The egg reduction rate was calculated as the reduction in the intensity of infection assessed indirectly, using egg count via the following formula: % egg = 10mL reduction (1−arithmetric mean of eggs/10mL urine after treatment in Year 5) x 100 arithmetic mean of eggs = 10mL urine before treatment at Baseline.
To assess overall population-level programmatic impact on the prevalence of S. haematobium between study arms, as per the standard analysis plan, agreed prior to analysis, Generalised Estimated Equations (GEE) were used to estimate the differences between arms in year 5 only; unadjusted estimates using only village and study arm as covariates and adjusted estimates including sex, age and weighting for number of children who provided data were modelled using SAS software 9.1.2. (SAS Institute Inc., Cary, NC, USA). The models were run on the three cross-sections of data: 9-to-12-year olds, first-year students (5-to-8-year olds), and adults (20-to-55-year olds). Prevalence was modelled using a binomial GEE with logit link function with Village IDs treated as the repeated measure and 'Lsmestimate' used to test prespecified differences between arms. In some villages, the enrollment was less than the target study population size. To account for this, a village-weight term was added to the GEE model to weight results according to numbers of children tested per village. Mean number of eggs per 10mL urine (mean intensity of infection) was modelled with negative-binomial GEE and log link function. To avoid issues of multiple testing we focused on specific arm comparisons. These comparisons were: Arm 1 (cccc) vs Arm 2 (ccss), Arm 1 (cccc) vs Arm 3 (cchh), Arm 1 (cccc) vs Arm 4 (ssss), Arm 4 (ssss) vs Arm 5 (sshh), and Arm 4 (ssss) vs Arm 6 (shsh) (Fig 1).
To assess the significance of the observed reduction in prevalence of S. haematobium in adults in the SBT arms (4, 5, and 6) we performed a binomial Generalised Linear Mixed Model (GLMM) with logit link function with infection status (0/1) as the response, age (standardised), sex, study arm and study year as fixed effects and village as a random effect on a subset of the parasitological data including only adults in the school based study arms for Year 1 and Year 5.
Nominal treatment coverage stipulated in the study protocol was 75%. In practice this was not achieved in many places. Despite additional data collection and considerable analytical effort, we were unable to adequately validate coverage; for this reason further analysis incorporating treatment coverage would be inappropriate. As the issues appear to be consistent throughout the study villages, the randomised allocation of villages to study arms means between arm comparisons could still have validity.
Kulldorff's space-time scan statistics were used to detect high-risk spatial clusters for schistosomiasis infection using the SatScan software version 9.1.1. (Information Management Services, Inc., Boston, MA, USA; www.satscan.org) [28,29]. To detect the high-risk spatial clusters of cases, the statistics use a moving elliptical window scanning the study area, and the maximum size of which is no more than 50% of the total population. We reported statistically significant clusters with an indicated p value <0.05. The SaTScan program was used to obtain reported observed cases, expected cases, relative risk of infection, and locations of specific cluster.

Results
Data were collected in the dry season between July and October each year from 2011 until 2015 in ten districts of Cabo Delgado province, Northern Mozambique. It was intended to sample 105,000 people from a total of 150 villages over the five years, but due to poor school enrolment and absence from many villages by children and adults working in the fields and mines, our final analysis included 81,167 individuals (77.3% participation). One village dropped out of the study prior to data collection due to political issues, therefore data is available from 149 villages.
Change in village-level prevalence from baseline to Year 5 is shown in Fig 2 and the SatScan in S1 Fig. Overall, forty-three percent of the population was female. Prevalence of infection decreased from baseline (60.5%) to Year 5 (38.8%), with prevalence of heavy infection also falling from 17.6% at baseline to 11.9% in Year 5 ( Table 1) Table 1 for gender, prevalence and intensity of infection across all age groups from baseline to Year 5.

9-to-12-year old cross-section
Over five years, 42,731 independent observations were taken of children aged 9-to-12 years from 149 communities randomly allocated across six study arms. See S1 Table for a summary of S. haematobium infection rates for each study year and study arm (9-to-12-year olds only). The results in Table 2 and Fig 2 demonstrate the impact of treatment strategy on the prevalence and proportion of heavy infection of S. haematobium for the 9-to-12-year olds at baseline and Year 5. The initial prevalence of S. haematobium infection was 66.7% across all study arms for this age group, with 22.8% prevalence of heavy intensity of infection. Our findings show a reduction in Year 5 where S. haematobium infection decreased to 42.5% over all treatment arms and the proportion of those heavily infected reduced to 13.4%. The arithmetic mean intensity of infection at the village level, reduced over the course of the study from 69.1 eggs per 10ml urine to 58.1 eggs per 10ml, with changes in mean intensity between the study arms. The egg reduction rate from baseline to Year 5 was 13.5% (Table 2).   Prevalence by gender demonstrated a higher prevalence of infection, as well as a greater proportion of those heavily infected, among males (p<0.001) ( Table 3). Response to treatment was similar by gender, with both males and females showing reduced infection over time across all arms, except for Arms 3 and 5 (Fig 4). Thus, prevalence and intensity of infection appears to have increased after the treatment hiatus. The effects of treatment were analysed using Arm 1, the most frequent and wide treatment approach, as the reference case. Arm 1 (annual CWT), had a comparable impact on prevalence as the less intense Arm 6 (SBT every two years). Arm 1 had a greater reduction on prevalence than Arm 2, 3 and 4 across all years, although the difference was not statistically significant (Table 3). Arm 1 did have a significantly larger impact on S. haematobium infection than Arm 5 in Year 5 (p<0.024).
Significance testing to assess whether the prevalence of S. haematobium varied over time between study arms showed that gender (p<0.001) and age (p<0.001) had a significant effect on infection status with higher prevalence and proportion of those heavily infected seen in males and among the 9-to-12-year age group in Year 5 (Table 3 and S2 Table). Across study year, prevalence of infection significantly decreases over the five years (p<0.001) with no statistically significant effect across the study arms. For the SBT arms, Arm 4 (annual SBT) was used as the reference arm to assess the impact of treatment holidays on prevalence. There was no significant difference seen following holidays. The impact of four versus two treatments was also assessed with the GEE model, where four treatments over four years (for both CWT and SBT) had a significantly greater impact in reducing prevalence of infection (p<0.024)

Variable
Arm 1 (cccc)  than two treatments over four years. There was no significant impact between the rounds of treatment and intensity of infection (S3 Table).

Other age cross sections
In Year 1 and 5 of the study an additional 31,542 observations were recorded for first-year students (aged 5-to-8-years) and 50 adults (aged 20-to-55-years) ( Table 4). At baseline, examination revealed that of the 7,456 first-year students, 62.8% were infected with S. haematobium and out of 4,253 adults, 44.5% were found to be infected. The proportion of heavily infected first-year students was higher compared to adults, with 19.3% and 7.1%, respectively. In Year 5 prevalence had fallen to 39.3% and 26.4% among first-year students and adults, respectively, and proportion of heavily infected had decreased to 12.8% and 5.6%, respectively. Significance testing to assess whether the prevalence of S. haematobium differed between study arms across the first-year students and adults in Year 5 are shown in Table 3 and S4  Table. From the GEE prevalence models adjusted for age and sex, gender was a significant Impact of different treatment approaches on urogenital schistosomiasis effect in first-students only (p<0.001) with males displaying a higher probability of being infected than females. Increasing age was also significantly associated with infection among first-year students, with the oldest respondents (8 years) being more likely to be infected than seven and six-year olds (p<0.001). There was no statistical association between age or sex with risk of infection in adults. For first-year students and adults, both the adjusted and unadjusted prevalence models show no significant difference in prevalence between the study arms at the end of the study.
Of particular interest is the statistically significant reduction in prevalence of infection for adults from Year 1 to 5 in the SBT arms (Fig 5). Adults did not receive treatment in these arms, however, the GEE prevalence model shows no significant difference in prevalence for adult's focal comparison arms. This would suggest that SBT alone has an impact on disease transmission.

Treatment coverage by strategy over time
Drug coverage was defined as the proportion of individuals who have ingested the PZQ [7]. The denominator is the estimated target population. The SBT coverage was determined by the Table 3 Impact of different treatment approaches on urogenital schistosomiasis percentage of school-aged children (5-12-years) that were treated, including those not enrolled in school. CWT coverage was determined by the total population treated using a denominator of the whole population based on the census carried out by the SCORE study in 2011. Coverage by study arm for each year is outlined in Table 5.

Estimate (95% CI) Pr > |z| Estimate (95% CI) Pr > |z| Estimate (95% CI) Pr > |z| Estimate (95% CI) Pr
The nominal treatment coverage stipulated by the study protocol was coverage of more than 75% for both SBT and CWT [7], which did not occur in many communities due to suspected inaccuracies in denominator values. The results shown in Table 5 shows that SBT coverage remained consistently low over time across all arms. CWT coverage, however, did increase over time across all arms, with no significant difference seen between arms.
When coverage was put into the multivariate model looking at impact on changes of prevalence over time, there was no clear trend or significant relationship with impact on infection.

Financial costs of different treatment strategies
Delivery costs (defined as the cost per person treated minus the drug costs) was $0.31 for SBT and $0.36 for CWT. The relative costs by each activity and input are summarised in Table 6. Our cost-analysis does not consider unpaid days of labour for individuals such as teachers, or other costs related to the SCORE field work. When the financial costs were broken down by activity, there was substantial cost saving for SBT for supervision and monitoring as this was carried out by the SCORE field teams during the prevalence survey. For CWT on the other hand, additional costs were incurred as supervisors were sent back to the villages post-survey to supervise and collect treatment registers.
According to input costs, there was a greater cost in per diems due to additional costs of community drug distributors salary and incentives (where one distributor was recruited for every 500 people in the village), whereas teachers were not paid extra in the schools but given incentives (a t-shirt) only. Fuel costs were also higher for CWT since the teams needed to revisit the communities a second time to collect any left-over drugs and the treatment registers. In SBT the treatment was done on the same day as the prevalence survey and the registers taken. The additional material costs for CWT were incurred due to additional treatment registers printed for larger numbers being treated.

Discussion
In Cabo Delgado, Mozambique, S. haematobium infection was found prevalent throughout all 149 schools in the study, where many schools had 90% baseline prevalence (Fig 2). The primary outcome of the study was to examine the change in infection among 9-to-12-year olds between specified study arms. Our findings showed a significant decrease in S. haematobium infection in this age group from a mean prevalence of 66.4% in Year 1 (of which 22.5% were heavily infected) to 42.5% in Year 5 (when 13.4% were heavily infected). When comparing the impact of treatment strategy each year, the only significant response to study arm was seen between annual SBT (arm 4) and two years of SBT followed by two consecutive treatment holidays (arm 5), which resulted in a higher final prevalence than when there were no treatment Table 4  holidays. There were no significant differences in Year 5 prevalence for the other treatment comparisons. The general uptrend in prevalence and intensity of S. haematobium infection shown after two subsequent holiday years, compared with one-year post-treatment, demonstrates rebound of S. haematobium if PC is interrupted. These findings of significant risk for re-emergence of infection if treatment is suspended supports previous schistosome infection modelling [30,31].

Mean infection intensity
Other studies have shown that the highest prevalence and intensities of infection occurred in males and among the 9-to-12-year age group, with infection decreasing in adulthood [32][33][34][35][36]. Despite a higher starting prevalence in adolescent males, there was no statistically significant difference in response to treatment by gender or age across the study arms. There is a need to explore further whether boys in these regions are more susceptible to subtle morbidity and impaired growth, mental or physical performance than girls.
The limitations of the current diagnostic methods for urogenital schistosomiasis have been demonstrated elsewhere, where the current 'gold standard' of urine filtration has been shown across all study arms. Mean intensity of infection decreased in first-year students and adults in arm 1 and 2 only. In the remaining study arms (3)(4)(5)(6), mean intensity levels increase from baseline to Year 5 in both first-year students and adults. It appears therefore that although the number of individuals infected reduces over time, the burden of infection in infected individuals seem to be heavier. It is not clear whether individuals harbouring a high parasite burden in such highly endemic populations are being re-infected each year, are more likely to be non-compliant in terms of taking the treatment, or if the parasites have developed resistance to treatment.
https://doi.org/10.1371/journal.pntd.0006061.g005 to be insensitive, especially in areas with low endemicity, and egg secretion varies from day-today [37,38]. However, in the context of Cabo Delgado where the prevalence and intensity of infection even in Year 5 remained high, it is likely that infections were not underestimated. Furthermore, aim of this study is to look at relative reduction in infection over time and given that microscopy was done consistently over the years, the trend in change is believed to be valid. As shown in the SatScan spatial scan, there is a decrease in the number of high prevalence spatial clusters from two western clusters at baseline to a single and less scattered central cluster in Year 5. A possible explanation for this is that as the prevalence of infection decreases over time, heterogeneities in the force of transmission by geography or social group become more visible. Such geographical clustering has also been seen in other SCORE projects in Kenya [39]. Further research is therefore necessary to understand how transmission dynamics change by space and social group as prevalence reduces, particularly among more resilient places/groups, and therefore how to adapt interventions accordingly.
There is a growing body of evidence demonstrating the burden of infection and morbidity in adults, as well as their potential role in sustaining transmission [14,30,31,[40][41][42]. Indeed, the results here show that in arms 3-6 mean intensity levels increase from baseline to Year 5 in adults. It appears that although the number of individuals infected reduces over time, the burden of infection in those infected seem to be heavier. It is not clear whether individuals harbouring a high parasite burden in such highly endemic populations are being re-infected each year, are non-compliant in taking the treatment, or if the parasites have developed resistance to treatment. This implies a greater need to include adults in schistosomiasis PC treatment programmes in some endemic settings. The SCORE study is therefore crucial as it investigates the impact of CWT in reaching more adults. The findings from Mozambique, however, have Impact of different treatment approaches on urogenital schistosomiasis demonstrated a significant reduction in prevalence of schistosomiasis in adults even in the SBT arms, implying transmission in the community has been decreased, even where only school children have been treated. These results support previous findings that demontrates the potential benefit of switching to community-wide mass treatment is highly variable in different contexts [43,44]. The delivery cost for SBT and CWT in this SCORE study is summarised in the Table 6 where $0.31 was calculated for SBT and $0.36 for CWT per person treated. Although this difference is small on an individual basis, in 2016 the Mozambique Ministry of Health treated around 8.5 million people and therefore these results therefore have significant logistical and cost-saving implications for a national control programme as they challenge the justification for CWT. Furthermore, a crude analysis of the relative cost per person treated by percentage reduction in infection for each study arm implies that even though two treatments over five years reduces prevalence less than four treatments, two treatments has a greater percentage reduction by relative cost. Village-level data, both in communities and at the school, were collected about water sources, sanitation, and hygiene (WaSH) issues through key-informant interviews. The baseline results from the WaSH survey have been published elsewhere [21]. In brief, among our study population, use of open surface water was widely reported for washing and bathing, and although pit latrines were reported to exist in most communities, open urination was widely practised. Furthermore, nearly all the communities in this study rely on subsistence farming whereby a large proportion of villagers spend months at a time sleeping in tented accommodation whilst working in the 'machambas' (fields), with no access to sanitation. Despite, such reported high-risk behaviour when the WaSH data was incorporated into the multivariate model, there were no significant findings associated with change in prevalence and intensity of infection across the study arms.
The most plausible explanation for why more intense treatment approaches did not result in a better impact on infection reduction was most likely due to poor treatment coverage. The nominal treatment coverage stipulated by the protocol was more than 75% for both SBT and CWT, which did not occur in many communities. A qualitative Knowledge Attitudes and Practices (KAP) survey was carried out in Year 3 to understand the challenges to achieving sufficient coverage, this showed contextual issues extraneous to the study that included extremely low school attendance (around 30% in majority of communities), high proportion of individuals (including school-aged children) working in fields or mines during the day, acute flooding which proved challenging to the CWT, severe political unrest and subsequent mistrust in any health-related programmes following cholera outbreaks, and poor study sensitisation in some areas. Although several days were spent by the survey team in each community to mobilise inhabitants, many communities were sparsely populated for several weeks whilst people were living in their fields or away in the mines for extended periods of time, which also highlighted the issue of unreliable denominator data. Nevertheless, there was an increase in CWT coverage over time that was due raised awareness of the programme and less trepidation in taking the treatment as the programme progressed, distribution of t-shirts among key community individuals in Year 3 and in Year 4 food was provided free-of-charge with treatment to help increase treatment coverage. Other studies have shown that provision of a pre-treatment snack can result in reduced side-effects as well as decreased prevalence of schistosomiasis [45]. These strategies proved effective with regards to community treatment coverage, but less impact on school treatment coverage. Such challenges in attaining good coverage is a reality for national control programmes, particularly in the light of poor school attendance [46]. There it is even harder to achieve than in a study setting such as ours, where there are generally more time and resources allocated to sensitising and motivating communities to participate.

Conclusions
The results obtained show that Cabo Delgado province in Northern Mozambique is highly endemic for urinary schistosomiasis. Although PC was successful in reducing the burden of active infection in all study arms, a more intense treatment approach, such as annual CWT, did not have as much impact as expected. This may be due to extremely high starting prevalence and intensity in the study area, with frequent exposure to reinfection [47]. It also may be related to poor treatment coverage and the challenge of performing a long-term study in a resource-poor setting, which may not be valid in other locations. These findings highlighted the need to complement mass drug administration with extensive public health education and behavioural modifications to reduce the risk of reinfection, and to consider the use of vector control through mollusciciding in areas where compliance to treatment proves to be challenging [38,48]. There is a need to explore further those individuals harbouring a high parasite burden in such highly endemic populations to understand whether they are being re-infected each year, are more likely to be non-compliant in terms of taking the treatment, or if the parasites might have developed resistance to treatment. The results of this study have also highlighted the challenge of achieving high treatment coverage especially with such highly migrant communities, specifically acquiring accurate population figures, and the importance of independent coverage surveys that do not rely on denominator values to accurately evaluate the impact of a treatment programme. Furthermore, more frequent treatment had a greater impact on prevalence of infection when arms were grouped by number of treatments, however, the costs efficiency was greater in arms only receiving two treatments. Finally, the significant reduction in prevalence of schistosomiasis among adults even in the SBT Arms implying transmission in the community has been decreased, has significant logistical and cost-saving implications for a national control programme as they challenge the justification for CWT.  Table. Summary of S. haematobium infection rates for each study year and study arm (9-to-12-year-old cross section only). Village-level S. haematobium prevalence, prevalence of heavy infection, and mean intensity for each study year from baseline to Year 5, by study arm among 9-12-year old cross section only. (TIFF) S2 Table. GEE prevalence-and intensity-adjusted parameter estimates (9-12-year olds only). (TIFF) S3 Table. GEE logistic model to assess differences in prevalence and intensity of S. haematobium infection between four versus two rounds of treatment in 9-12-year olds. (TIFF) S4 Table. GEE logistic model to assess differences in prevalence and intensity of S. haematobium infection in first-year students only. (TIFF)