Abstract
Schistosomiasis affects more than 230 million people worldwide. Control and elimination of this parasitic infection is based on mass drug administration of praziquantel (PZQ), which has been in use for several decades. Because of the limitations of the efficacy of PZQ especially against juvenile worms, and the threat of the emergence of resistance, there is a need to consider alternative formulations or delivery methods, or new drugs that could be more efficacious. We use an individual-based stochastic model of parasite transmission to investigate the effects of possible improvements to drug efficacy. We consider an increase in efficacy compared to PZQ, as well as additional efficacy against the juvenile life stage of schistosome parasites in the human host, and a slow-release formulation that would provide long-lasting efficacy for a period of time following treatment.
Analyses suggest a drug with a high efficacy of 99%, or with efficacy lasting 24 weeks after treatment, are the two most effective individual improvements to the drug profile of PZQ. A drug with long lasting efficacy is most beneficial when MDA coverage is low. However, when prevalence of infection has already been reduced to a low level, a high efficacy is the most important factor to accelerate interruption of transmission. Our results indicate that increased efficacy against juvenile worms can only result in modest benefits, but the development of a new drug formulation with higher efficacy against adult worms or long-lasting efficacy would create an improvement to the community impact over the currently used formulation.
Author summary The World Health Organization has set the target of elimination of schistosomiasis as a public health problem by 2030. Currently, praziquantel is the sole drug used in mass drug administration (MDA) for schistosomiasis, raising concerns about the potential evolution of resistance and whether an improved drug profile would be necessary to achieve elimination and the more ambitious goal of interruption of transmission. We present the results from an individual-based stochastic mathematical model that simulates schistosome transmission and the impact of MDA. Three potential improvements to the properties of the drug are considered: improved efficacy, long lasting efficacy post a single treatment and improved efficacy against the juvenile life stage of the schistosome. Our findings reveal that with good coverage, an improved drug efficacy is best at reducing prevalence and achieving interruption of transmission. However, when MDA quality is compromised due to low coverage, infrequent treatment or high rates of non-adherence, then a long-lasting efficacy performs best. There is only a modest benefit of improved drug efficacy against juvenile schistosomes. These results highlight the importance of good MDA program and inform future drug development aims.
Introduction
The human schistosome parasites are parasitic flatworms that cause schistosomiasis, a disease that affects more than 230 million people worldwide (1). There are three main species responsible for human infection: Schistosoma mansoni, Schistosoma haematobium and Schistosoma japonicum. Each species cover different geographic areas, S. mansoni is endemic in several countries in Africa and the Caribbean, S. haematobium is endemic in Africa and Middle East, and S. japonicum is endemic in Asia (2). Like several other neglected tropical diseases (NTDs), mass drug administration (MDA) is a key component in the control and elimination programmes of schistosomiasis. Currently, praziquantel (PZQ) is the only drug available for use against schistosomes and has been in use for over 40 years (3).
Because of the high disease burden, schistosomiasis has been targeted for elimination as a public health problem by the World Health Organization in several countries by 2030 (4). However, attaining this goal will likely require very high levels of coverage of MDA (5), as well as improvements in access to safe water, sanitation and hygiene (WASH) and possibly some form of snail control (World Health Organization, 2022). As there is the potential for Schistosoma species to develop resistance against PZQ, which also has other limitations such as a low efficacy against the juvenile stage in the schistosome life cycle, alternative drugs have been considered (8,9). An ideal drug would be efficacious against all stages of the life cycle and include other benefits such as a slow-release mechanism that could act to provide a longer duration of protection post treatment.
The schistosome life cycle is complex, going through multiple development stages (10,11). When a free-living cercaria infects a human host by penetrating the skin, it goes on a journey through the body, via the vascular system. Juvenile worms travel to the lungs and the liver, before reaching the intestine (S. mansoni and S. japonicum) or bladder (S. haematobium) where they typically remain for 5-10 years as adult worms (12).
Detailed examination of the life cycle within a human host is difficult but estimates of the length of time between infection and egg production in mice hosts has been estimated to be 36-38 days (13). In that time, the proportion of juvenile worms surviving to maturity has been found to be approximately 40% (10,14). Studies of the human immune response have found a similar time span of 4 to 6 weeks between infection and detection of antibodies targeting adult worms (15).
PZQ has been found to be ineffective against juvenile worms, possibly due to their location within the body at the time of treatment (8,16,17). For example, one study found that the concentration of PZQ required to reduce the worm burden by 50% was 30 times higher when treating after 4 weeks, compared to 7 weeks (14). Alternatives to PZQ have been suggested for their efficacy against juveniles (8,18), however the impact that this would have at a population transmission level is unknown.
An alternative approach to targeting juveniles at the time of treatment, is to instead administer a slow- release drug, which will continue to kill adult worms for several weeks post pill swallowing. Research has shown that the slow-release delivery of drugs is effective against schistosomes (19,20). If such a treatment were applied during MDA, it could theoretically target juvenile worms as they mature. If the period of slow release is sufficiently long (longer than the typical maturation period of juveniles), or if the new drug was efficacious against juveniles as well, there is also the added benefit of giving lasting protection against new infections for a period after treatment.
Intuition suggests that coverage and efficacy are the two most important factors in an MDA program and the product of the two will give a good reflection on its effectiveness. However, intuition is less reliable when additional parameters are included, such as the length of time that efficacy will last after treatment. In cases such as this, mathematical models are valuable tools for helping to formulating policy and in drug design, by facilitating quantitative analyses of the impact of different improvements at the population level in the drug profile of action on the parasite.
In this paper, we investigate the impact of drug profile improvements to the current practice of using PZQ in MDA, employing an individual based stochastic mathematical model of the transmission and mass drug treatment of S. mansoni and S. haematobium. We assess via simulations additional efficacy against juvenile worms as well as adult worms, long lasting efficacy that gives protection against adult worms for several weeks post treatment, and an improvement in drug efficacy against adult worms only. To measure the effectiveness of each of these improvements, we show the average time taken to reach interruption of transmission based on 200 repeated simulations, i.e. zero prevalence in the entire population.
Methods
Model formulation
We use a stochastic individual based model developed at Imperial College London (the Imperial model) that is based on the deterministic model that has been described in previous publications (21,22). To include the juvenile stage in a human we adjust the model of mean adult worm burden as a function of age and time, M(a, t), to include the mean juvenile worm burden J(a, t), so that in the deterministic version,
where L represents the concentration of the infectious material in the environment, β(a) is the age- related contact rate, μ is the maturation rate of juvenile worms, and σj and σm are the natural death rates of juvenile and adult worms respectively.
In the Imperial model, infection and maturation of juvenile worms and death of juvenile and adult worms are discrete events, modelled as a Poisson process. All other details are as described in previous publications (23–25); we assume the worms are monogamous and have a negative binomial distribution which arises from a gamma distributed predisposition to infection amongst the individual human hosts. Individuals also have an age-dependent contact rate with infectious material, and the same age-dependence is also included in the output of eggs, which is density-dependent (26).
We also replicate some of the simulations with the SCHISTOX model, developed at the University of Oxford (27,28).
MDA treatment
Starting in a closed, untreated population, we simulate biannual MDA, with a coverage of 75% in the community (everyone over the age of 2). We assume that treatment is delivered at random at each round. We first consider a ‘baseline’ approach, whereby MDA of PZQ occurs with an efficacy of 86.3% against S. mansoni, or 94% against S. haematobium (29). The effect of the drug is assumed to be immediate, so that any adult worm in a host that is treated will have a probability of dying at the time of treatment.
When considering improvements to treatment that include efficacy against juvenile worms, they are also killed at the time of treatment, with the same efficacy as for adult worms. In the slow-release treatment, we assume that MDA still occurs as normal, but there is a fixed period of time following treatment whereby a juvenile worm that reaches maturation will die with the same probability as the original drug efficacy. We consider efficacious periods of 8, 16 and 24 weeks. When there is a slow-release treatment that also has efficacy against juveniles, the worm will be killed at the time of infection, rather than the time of maturation.
Initially we investigate the impact of improving drug efficacy to 99%, with different levels of coverage and at different baseline prevalences (R0 values) and record the length of time taken to reach interruption of transmission (IOT) up to a maximum of 20 years after the start of the MDA programme. We define this as the time taken for 90% of the 200 simulation runs to reach a 0% prevalence of infection, using the Imperial and SCHISTOX models (22,27). Later, we consider the impacts of additional improvements to the drug profile by adapting the Imperial model. Here we restrict our focus to a scenario with high levels of transmission (R0 = 1.4-4, k=0.24), with an initial baseline prevalence of >50%. For each improvement to drug efficacy, we show the mean prevalence and the mean time to reach different prevalence thresholds, including IOT (0% prevalence). We also consider alterations to the MDA whereby treatment occurs annually instead of biannually, or a small proportion of the population (the ‘never treated’, set at a range of values between 1 to 5%) are permanently excluded from MDA and never receive treatment. The parameter values used in model simulations are documented in Table 1 for S mansoni and in Table 2 for S haematobium.
Results
The effect of drug efficacy on the ability to achieve IOT within 20 years from different baseline prevalences is shown in Table 3, along with the effect of MDA coverage levels and the proportion of the population that is never treated. For both age profiles used to model S. mansoni, where IOT can be reached within 20 years, an increased 99% efficacy can reduce the time taken by up to four years. The benefits are smaller with S. haematobium but note that we assume the efficacy of PZQ is already 94% for this parasite so the increase in efficacy is smaller. Most importantly, in both cases, if the coverage with PZQ is 60%, a bigger decrease in time to IOT can be attained with an increase in coverage compared to an increase in efficacy. The largest impact on reaching IOT is the proportion of the population that are systematically never treated. In scenarios with a high baseline prevalence (high R0 value), it is impossible to reach IOT if any of the population are never treated. In moderate prevalences it is possible in some scenarios when the never treated is 1% and in low prevalence settings (small R0 values but >1) it is possible but can require up to 10 years of additional MDA. If the proportion never treated is as high as 5% then it is not possible to achieve IOT for any scenario we considered.
Number of years of MDA required for interruption of transmission when varying the baseline prevalence, MDA coverage, the proportion of the population systematically never treated and the drug efficacy. We model S. mansoni with a low and high adult burden and S. haematobium. MDA is annual when the baseline prevalence is low or moderate and biannual if the baseline prevalence is high.
To examine the effects of additional improvements to a potential new drug formulation, we now restrict our attention to an area of high baseline prevalence (high values of the basic reproductive number R0 – R0 = 3.3-4 for S mansoni, R0 = 2 for S haematobium). As shown in Table 4, in this scenario the entire community must be available for treatment (0% never treated) to achieve IOT, no matter the improvements to the drug profileTable 4. If this is the case, a new drug with an efficacy against juvenile and adult worms, lasting 24 weeks can reach a 0% prevalence of S. mansoni faster than the current drug efficacy profile by up to four years, depending on the age intensity profile, equivalent to six rounds of MDA. The same drug would only be two years faster when used against S. haematobium, probably due to the already high efficacy of PZQ against this species. This drug with the combination of efficacy against juvenile worms and a slow-release period is the fastest approach to achieve IOT in this scenario but is only marginally better than a slow-release drug without efficacy against juveniles, or one that has a near-perfect 99% efficacy against adult worms. Little benefit is seen from a drug that has efficacy against juvenile worms without any additional improvements in the drug properties profile.
Number of years of MDA required for interruption of transmission, starting at a high baseline prevalence (>50%). The baseline efficacy is 86.3% for S. mansoni and 94% for S. haematobium.
These results are broken down further for the S. mansoni, low adult burden case in Figure 1, showing the reduction of prevalence, the probability of achieving elimination over time and the time to reach different prevalence thresholds. This shows how the best performing drug, with 99% efficacy against adult worms, quickly reduces the prevalence in just over 6 years on average, compared with the baseline case (86.3% efficacy against adult worms), which takes just under 8 years in total. Note that results from the bar chart differ from those in Table 4, as the former was calculated using the mean, whereas the latter shows the time for 90% of simulations to reach elimination.
Results of biannual MDA with 75% community coverage against S. mansoni with a low adult burden, starting with a high baseline prevalence. (A) Prevalence, based on worm counts during MDA with a range of drug improvements over time. (B) The probability of reaching interruption of transmission with each drug improvement over time. (C) Years of MDA required to reach interruption of transmission (bar height) and the time taken to reach different prevalence thresholds (coloured blocks).
Interestingly, most improvements to the speed of the reduction of prevalence are realised early during MDA, when the prevalence is still high. For example, a drug with a 24-week efficacy against adult and juvenile worms reduces the prevalence from 70% to 50% in 0.8 years and 50% to 25% in 1 year, compared to the baseline case which takes 1.1 and 1.7 years to reach those intervals respectively. However, the time taken to reduce the prevalence from 10% to 1% and from 1% to 0%, takes 1.6 and
1.2 years with a 24-week efficacy period, compared to 1.8 and 1.1 years for the baseline case. The only improvement that outperforms the baseline case consistently throughout MDA is a drug with the high 99% efficacy. A 99% efficacy drug outperforms all other drug characteristics at high (>50%) and low (<10%) prevalences, but it is not as effective as a drug with efficacy lasting at least 16 weeks at the moderate prevalences (10%-50%).
If MDA is administered annually rather than biannually, as shown in Figure 2, drugs with long lasting efficacy (slow release) perform better than any other. If efficacy were to last for 24 weeks, the time taken to reach elimination would be 7 years quicker than the baseline case. In this scenario, a drug with 99% efficacy performs similarly to one which has 86% efficacy against juvenile worms, which achieve elimination approximately two years quicker than the baseline case. However, once a 10% prevalence threshold is reached, a drug with 99% efficacy is the better option, taking 4 years to achieve elimination compared to approximately 5 years in all other cases.
As in Figure 1, except MDA is administered annually instead of biannually. (A) Prevalence, based on worm counts during MDA with a range of drug improvements over time. (B) The probability of reaching interruption of transmission with each drug improvement over time. (C) Years of MDA required to reach interruption of transmission (bar height) and the time taken to reach different prevalence thresholds (coloured blocks).
A similar effect is seen when the frequency is restored to biannual MDA, but community coverage is reduced to 40%, as shown in Figure 3. The time to achieve IOT is approximately 25 years but can be reduced to approximately 16 years with a drug that has efficacy lasting 24 weeks. In both cases where fewer treatments are administered, the reduction in time to reach each prevalence threshold and achieve IOT from an improved drug characteristic is greater than when coverage and frequency is high. This is also the case when the efficacy of PZQ is lower than the 86.3% we have otherwise assumed, as shown in Figure 4 where the baseline efficacy is 60%. If this were the case then the benefit of a new drug with 99% efficacy is obvious, but a slow-release formulation lasting 24 weeks would decrease the time to achieve IOT, even with the same 60% efficacy.
As in Figure 1, except MDA coverage is 40% of the community. (A) Prevalence, based on worm counts during MDA with a range of drug improvements over time. (B) The probability of reaching interruption of transmission with each drug improvement over time. (C) Years of MDA required to reach interruption of transmission (bar height) and the time taken to reach different prevalence thresholds (coloured blocks).
As in Figure 1, except the efficacy of PZQ is assumed to be 60%. (A) Prevalence, based on worm counts during MDA with a range of drug improvements over time. (B) The probability of reaching interruption of transmission with each drug improvement over time. (C) Years of MDA required to reach interruption of transmission (bar height) and the time taken to reach different prevalence thresholds (coloured blocks).
Figure 5 illustrates the importance of good adherence amongst the population eligible for treatment. Here, MDA is administered biannually again but 1% of the population is systemically excluded from treatment. As in Table 4, this shows that IOT is difficult to achieve when this is the case. While some drug improvements decrease the time needed to reduce prevalence to 10%, the time taken to reach elimination from a 10% or 1% prevalence is much larger than when everyone adheres to treatment during the course of MDA. Once the 10% prevalence threshold is passed, there are no benefits from improvements to treatment and the lines in Figure 5A and B converge, while the variation in bar heights in Figure 5C is very limited.
As in Figure 1, except 1% of the population are systematically excluded from MDA and are therefore never treated. (A) Prevalence, based on worm counts during MDA with a range of drug improvements over time. (B) The probability of reaching interruption of transmission with each drug improvement over time. (C) Years of MDA required to reach interruption of transmission (bar height) and the time taken to reach different prevalence thresholds (coloured blocks).
Although the above results show the potential benefits of improved drug characteristics, it is clear from comparing the baseline in Figure 1 with Figures 2-5 that delivering frequent MDA with good coverage with no proportion never treated is more important than the characteristics of the drug. To illustrate this, we compare the effect of PZQ with good coverage and assuming an efficacy of 86.3%, against an ideal new drug that has a 99% efficacy against adults and juveniles and a 24-week slow release period, with suboptimal MDA. The results in Figure 6 show that while the ideal drug is very effective at bringing down prevalence and achieving IOT within 5 years (panels A and E), when the frequency of MDA is reduced to annual administration the results are very similar to the effects of PZQ with biannual MDA (panels B and F). When coverage is reduced to 40% (panels C and G), the ideal drug will initially reduce prevalence faster than PZQ with ideal MDA but is much slower to achieve IOT. This is also the case when 5% of the population are systematically never treated (panels D and H), where the prevalence falls much quicker with the ideal drug, but IOT is impossible to achieve within 15 years and prevalence is lower when using PZQ with ideal treatment after approximately 5 years.
A comparison between PZQ with ideal MDA, defined as 75% community coverage, biannual, 0% NT (red lines), and an ideal proved drug with 99% efficacy against juvenile and adult worms, lasting for 24 weeks. Individual plots show prevalence (A-D) and robability of interruption of transmission (E-H) for 15 years after starting MDA. (A, E) Comparison between PZQ (red) and an ideal rug (blue), both with ideal coverage. (B, F) The ideal drug is given annually instead of biannually (cyan). (C, G) Coverage of the ideal rug is reduced to 40% (yellow). (D, H) 5% of the population are permanently excluded from MDA with the ideal drug (magenta).
Discussion
Use of PZQ has long been the principal approach to reducing the prevalence of schistosome parasite infection and controlling morbidity due to high burdens of worms. However, the ambition to achieve the elimination of transmission in defined settings has led to research into new drugs, alongside other interventions such as improvements in water, sanitation and hygiene, plus snail control (the intermediate host in the life cycle) (8). Many of these ideas have yet to be tested at a population level. Past discussion about what might be a sensible set of properties for any new drug to support schistosomiasis control have been based largely on what technically might be possible (e.g. activity against juvenile worms and slow-release technologies). The research described in this paper sets out to provide a more quantitative background for such discussions by investigating via simulation what the population level impact of different drug properties would be based on the predictions of a mathematical model of the transmission of S. mansoni and S. haematobium and MDA impact. The magnitude of the improvements we consider are not based on products currently in development. However, the results provide guidance on the effect that each improvement will have at the community level.
Although a lack of efficacy against juvenile worms is often highlighted as a problem with PZQ, we find that at a population level this has little impact on reducing prevalence. This is due to the time it takes to mature in the human host being small relative to adult worm life expectancy. We have assumed that the juvenile maturation period is 5 weeks on average compared to the adult life span of 5.7 years for S. mansoni and 4 years for S. haematobium. This means that the number of juveniles as a proportion of the total number of worms in the population is small, approximately 1.5-2.5%. Therefore, even the impact of a drug with 100% efficacy against juvenile worms would only have the effect of improving efficacy against all worms by 2.5% at best.
It is therefore unsurprising that a drug with high efficacy against adult worms achieves better outcomes than one with an efficacy against juveniles but no improvement in efficacy against adults. A drug with 99% efficacy against adult worms is effective at decreasing the time to achieve IOT, particularly when coverage is also good. Similarly, it is unsurprising that most of the improvements don’t show substantially better results for S. haematobium compared to the baseline (see Table 4: Number of years of MDA required for interruption of transmission, starting at a high baseline prevalence (>50%). The baseline efficacy is 86.3% for S. mansoni and 94% for S. haematobium., because the efficacy is already quite high against this schistosome species. The estimate for the efficacy of PZQ used in this model is one that is commonly referred to (23,24), however most estimates of egg reduction rate range from 60-100% (40). If resistance is being developed, the efficacy of PZQ in the field could be much lower. If this is the case, or if efficacy is less than assumed here, the benefits of a new drug with higher efficacy are obviously highly desirable, but a slow-release period would also be beneficial, even if efficacy remained low. The caveat to this is that a slow-release formulation may affect the way resistance develops.
What is surprising, perhaps, is the result that when a low prevalence is already reached, there appears to be no benefit to drug improvements, with the exception of a drug with a very high efficacy. This is counter-intuitive, as one might expect that there would be value in killing juveniles when approaching elimination. However, at lower prevalences, the rate of new infections will have already slowed, and thus the number of juveniles will be small, whereas there will still be many adult worms in the population particularly if not everyone eligible has received treatment, or the efficacy against adults is not sufficiently high. A drug with high efficacy against adult worms only is also the most effective at reducing prevalence from a high baseline level as almost everyone who is treated will be free from infection. However, at moderate prevalences, drugs with a long-lasting efficacy are most effective as they provide protection from reinfection and reduce the immediate bounce back in prevalence after MDA.
In this paper we mostly focus on the goal of interruption or elimination of transmission and for the purpose of the model, define this as a complete absence of worms in the population. Although this can be achieved in a closed community, there is always the potential for infection to be reintroduced and prevalence to bounce back if infected people move into or visit areas that have cleared infection by high coverage of MDA over many years (41). In this scenario where prevalence is low, high efficacy against adult worms and good coverage of MDA, with no proportion never treated, are more important than the other potential improvements.
The improvements in the properties of an anti-schistosome infection drug were chosen to illustrate how they might impact infection at a community level. The main conclusion arising from the simulation studies is that high coverage and high efficacy are key to success, but long-lasting efficacy can also be very effective. This is particularly the case at moderate and high prevalences, and when MDA coverage is not ideal. Activity against juvenile worms is of limited additional benefit. When coverage is low, or some of the population are systematically excluded from MDA, even an ideal drug with very high efficacy that lasts for 24 weeks is not as effective as PZQ with good coverage. For many helminth infections WHO has recommended MDA coverage targets of around 75% for school aged children. We believe this is too low. Coverage targets should be around 80-90% of the total population in moderate or high transmission settings if transmission interruption is the goal. As has already been shown elsewhere, good coverage over the course of MDA year by year, with a very low proportion of never treated, is vital (5,28). WASH interventions are also strongly recommended as a key supplement to MDA and for long lasting general improvements in health (7).
Data Availability
All relevant data are within the manuscript. References
References
