Abstract
Background
The durations of untreated stage 1 (early stage, haemolymphatic) and stage 2 (late stage, meningoencephalitic) human African trypanosomiasis (sleeping sickness) due to Trypanosoma brucei gambiense are poorly quantified, but key to predicting the impact of screening on transmission. Here, we outline a method to estimate these parameters.
Methods
We first model the duration of stage 1 through survival analysis of untreated serological suspects detected during Médecins Sans Frontières interventions in Uganda and Sudan. We then deduce the duration of stage 2 based on the stage 1 to stage 2 ratio observed during active case detection in villages within the same sites.
Results
Survival in stage 1 appears to decay exponentially (daily rate = 0.0019; mean stage 1 duration = 526 days [95%CI 357 to 833]), possibly explaining past reports of abnormally long duration. Assuming epidemiological equilibrium, we estimate a similar duration of stage 2 (500 days [95%CI 345 to 769]), for a total of nearly three years in the absence of treatment.
Conclusion
Robust estimates of these basic epidemiological parameters are essential to formulating a quantitative understanding of sleeping sickness dynamics, and will facilitate the evaluation of different possible control strategies.
Background
Human African trypanosomiasis (sleeping sickness, HAT) has historically been a predominant parasitic infection in Africa, causing many millions of deaths in the late 1800s through early 1900s, and has reemerged in historical foci after breakdown of control programmes [1]. Two forms of HAT are recognised, due respectively to Trypanosoma brucei gambiense (a mainly human disease found mainly in western and central Africa) and Trypanosoma brucei rhodesiense (a zoonosis observed in eastern and southern Africa) [2]. Gambiense HAT has an insidious onset and progresses over years, while rhodesiense HAT is fulminant, causing death within a few months of infection [3]. In recent years, 10 000–15 000 HAT cases (mostly gambiense) are reported annually [4], but this is likely an underestimate of the true burden.
Fundamental aspects of the pathogenesis, clinical profile and epidemiology of HAT remain poorly understood. Infecting trypanosomes first colonise the haemolymphatic system, evading specific immunity through antigenic variation, but causing only mild and intermittent symptoms – this is known as stage 1 (early or haemolymphatic stage). Eventually, parasites cross the bloodbrain barrier and cause severe neurological consequences, systemic deterioration and, ultimately, death [5]. The period after bloodbrain barrier penetration is called stage 2 (late or meningoencephalitic stage).
The durations of both disease stages in gambiense HAT, key parameters in the determination of the reproductive ratio of HAT, are not well quantified, and understood only on the basis of informal analysis or anecdote. Both stage 1 and 2 are believed to last from months to years in gambiense disease. However, atypically fast or slow progressions are reported and gambiense HAT, just like rhodesiense, seems to exhibit a wide range of virulence, often following geographical patterns [3,6,7].
In a separate paper (Checchi et al., submitted), we review available evidence on the natural evolution of untreated gambiense HAT, and consider the possible phenomenon of trypanotolerance, including instances where patients might spontaneously clear their infections or become chronic, asymptomatic carriers. We find that most infections are likely to be pathogenic and fatal if untreated, but that chronic carriage cannot be ruled out based on available evidence, and would, even if occurring at low frequencies, play a crucial role in maintaining transmission even in HAT foci where intensive case detection campaigns are carried out. Trypanotolerance is a welldescribed phenomenon in cattle and other animal species [8], and has been postulated by some as an explanation for the mysterious persistence of various historical HAT foci [6,9].
In the present paper, we focus on pathogenic T.b.gambiense cases only, which probably constitute the majority of infections, and estimate indirectly the duration of both stage 1 and stage 2 from time of infection to death, in the absence of treatment, based on observations of untreated serological suspects, and the prevalence of the two stages in defined populations. Rigorous and robust estimation of these parameters is crucial to enable quantitative evaluation of the potential of interventions to reduce transmission by minimising the period of human infectiousness to flies [10]. Estimates for T.b.rhodesiense have been derived elsewhere [11,12] by analysing patientreported time from symptom onset to stage 2 presentation or death.
Methods
Overall study design
The progression of HAT infection from stage 1 to stage 2 and from stage 2 to death or any other outcome, in the absence of any treatment, can be modelled simply (Figure 1). Provided that model parameters specifying the progression rates from stage 1 (r_{1}) and stage 2 (r_{2}) can be estimated, the mean and median durations of both stages, and their confidence intervals, can easily be calculated. Furthermore, if the numbers of stage 1 and stage 2 cases are at steady state (equilibrium), the ratio of stage 1 to stage 2 prevalence (S_{1}/S_{2}) will also be constant over time. Assuming equilibrium, S_{1}/S_{2 }must be directly proportional to the relative durations of stage 1 and stage 2; thus, if one of the two durations is known, the other can be deduced.
Figure 1. Simple model for the progression of untreated HAT.
The challenge is to identify data from actual HAT foci that can inform model parameters; and in particular, to estimate r_{1 }and/or r_{2 }whilst maintaining the ethical obligation to treat parasitologically confirmed infections immediately. Here, we use data from Médecins Sans Frontières (MSF) HAT programmes to estimate (i) r_{1 }from a survival analysis of serological suspects kept under observation without treatment, and (ii) S_{1 }and S_{2 }based on the output of active screening sessions, in which communities are tested exhaustively over the space of a few days, yielding the point prevalence of each stage. We then deduce r_{2 }from the estimates of r_{1}, S_{1 }and S_{2}.
Analysis was performed with Stata 9.0 software (Stata Corporation, College Station, Texas, USA). Secondary analysis of these operational data was approved by the Ethics Committee of the London School of Hygiene and Tropical Medicine.
Estimation of stage 1 to stage 2 progression rate (r_{1})
We pooled serological screening and individual patient datasets from six MSF HAT control programmes (Moyo, Uganda, 1986–1993; Adjumani, Uganda, 1991–1996; Arua, Uganda, 1995–2002; Yumbe, Uganda, 2000–2002; Kiri (Kajo Keji), southern Sudan, 2000–2005; Maridi, southern Sudan, 1999–2005). Programmes and datasets have been described previously [13,14]. Briefly, cases are detected through a combination of passive (fixedcentre) and active (mass screening) case detection. Screening relies on the serological Card Agglutination Test for Trypanosomiasis (CATT, Institute of Tropical Medicine, Antwerp, Belgium). Persons positive for the CATT in whole blood go through a complex diagnostic algorithm consisting of a repeat CATT test at progressively higher dilutions, direct microscopy of aspirate from cervical lymph nodes (if present), microscopical examination of blood using one or more concentration techniques (the Quantitative Buffy Coat, Woo, and/or miniAnion Exchange Column tests), and a lumbar puncture for staging purposes (microscopic evidence of trypanosomes or a white blood cell [WBC] count >5/μL in cerebrospinal fluid [CSF] leads to stage 2 classification).
Frequently, individuals display a strong CATT reaction but no parasitological evidence of HAT infection and are thus classified as serological suspects. Because the CATT alone is poorly predictive (specificity is around 95% but infection prevalence is usually <5%, giving positive predictive values <50%) [15], treatment is withheld in such cases. Serological suspects are asked to return for repeat testing after three, six, nine, and twelve months, and are treated immediately if found to have parasitologically confirmed HAT, or released from further followup if the strong CATT reaction disappears.
Our analysis focussed on stage 1 serological suspects, defined as patients positive at a CATT dilution ≥1:4, but with no parasitological evidence of HAT in any fluids, and a CSF WBC count ≤5/μL. Among patients who were seen at least once after initial screening, possible endpoints at the last followup visit were defined as (i) noncase if the CATT became negative or positive at a dilution <1:4, in the absence of other parasitological evidence of HAT, (ii) persistent suspect if positivity at CATT dilution ≥1:4 remained but other tests were still negative, (iii) confirmed stage 1 if blood or gland aspirate became positive and CSF remained negative, or (iv) confirmed stage 2 if parasites were observed in the CSF, if parasites were observed only in blood or lymph fluid but CSF WBC count exceeded 5/μL, or if the CSF WBC count exceeded 20/μL in the presence of a CATT dilution ≥1:4, irrespective of other tests. Only patients who developed either stage 1 or 2 HAT during followup (endpoints (iii) and (iv) above) were retained in this study for further survival analysis. Among these, patients reaching endpoint (iv) were considered to have progressed to stage 2, and those with endpoint (iii) to have remained in stage 1 for the duration of followup.
Maximumlikelihood estimation was used to fit parametric models to the observed survival data. While estimating confidence intervals for model parameters, possible clustering of survival probabilities within project sites was adjusted for by considering each site as a cluster and taking into account withincluster correlation of observations. So as to select an optimal model, different hypotheses about the most closely fitting underlying distribution (exponential, Weibull, loglogistic, gamma, and Gompertz) were evaluated using Wald tests if models were nested (for example, the exponential distribution is a particular case of the Weibull), or the Akaike Information Criterion if they were not. Ordinarily, survival analysis assumes that events occur when they are observed (rightcensoring). In our case, because followup visits took place months apart, we decided against such an assumption, and instead treated observations of stage 2 progression as intervalcensored, i.e. considered simply that progression had occurred within a time interval defined by the visit when stage 2 was confirmed and the previous visit (when it was known that the patient had not yet progressed). This method implies likelihood expressions that are somewhat different from the right censoring case [16] and a userdefined Stata module for intervalcensored survival analysis was applied [17]. So as to provide some graphical assessment of the goodness of fit of the parametric models, we compared fits of these models to a nonparametric KaplanMeier curve of the survival data, with Greenwood 95% confidence intervals, censoring each progression event at the exact midpoint of the interval between the two visits between which it was known to have occurred (a close approximation to intervalcensoring).
The resulting model enables estimation of survival time in stage 1 after detection. However, upon detection patients will have already been infected for some time; the distribution of these times since infection is unknown, and is itself dependent on incidence trends in the predetection period. If incidence is stable, this distribution will be uniform; if incidence is increasing, proportionately more will have been infected recently. As illustrated in Figure 2, if survival is dependent on time since infection, the observed shape of the distribution of survival times after detection will reflect both the distribution of times since infection (itself a result of past incidence trends) and the true distribution of survival in stage 1 (i.e. r_{1}), as in scenarios B1, B2, B3, C1, C2 and C3. If survival is in reality timeindependent (constant rate), both the true and observed distributions will be exponential (scenarios A1, A2 and A3), and will share the same parameter value, irrespective of whether incidence is stable or not.
Figure 2. Expected shapes of the observed distribution of survival times after detection. Shaded cells A1C3 show the expected distribution of survival times after detection, given different scenarios for the true distribution of stage 1 survival (A, B and C) and incidence in the predetection period (1, 2 and 3). P = proportion of cases remaining in stage 1; t = time; * = exponential function , with mean stage 1 duration = 1/r_{1; }** = linear function P(t) = r_{1}t with fixed stage 1 duration = 1/r_{1}.
Estimation of S_{1 }and S_{2}
Data on active screening sessions were available from the Adjumani, Arua, and Kiri programmes, and were used to provide realistic values for S_{1 }and S_{2 }in situations unmodified by control. Screening sessions were retained for analysis only if they were the first to be conducted in the respective villages, reached a coverage (persons screened/total population) ≥70% (a commonly used target in HAT programmes), and yielded at least one case of HAT. The mean and median S_{1 }and S_{2}, and their ratio, were calculated after weighting for village population size. We also did sensitivity analysis to see whether passive case detection in the prescreening period affected the S_{1 }to S_{2 }ratio.
Estimation of stage 2 removal rate (r_{2})
Assuming that no treatment is available, the prevalence of cases in stage 2 at time t is given by the prevalence at the previous time point (t  δt), plus new progressions from stage 1, minus stage 2 patients removed, such that:
Assuming equilibrium conditions, S_{2}(t) = S_{2}(t  δt) = S_{2 }and S_{1}(t) = S_{1}(t  δt) = S_{1}. Thus, solving for r_{2 }gives
In an epidemiological scenario meeting the above assumptions, the S_{1}/S_{2 }ratio is thus linearly proportional to the relative duration of r_{2 }and r_{1}.
Results
Description of included and excluded patient cohorts
Out of 7742 patients listed as suspects in the six databases, 1079 were not eligible for the analysis: 385 (5.0%) had laboratory findings compatible with confirmed stage 1 or 2 HAT and were in fact treated; 577 (7.5%) were stage 2 serological suspects; and 117 (1.5%) had an unclear stage due to missing variables. This left 6663 (86.1%) patients meeting stage 1 suspect criteria.
Baseline characteristics (Table 1) differed significantly across sites (p < 0.001 for age; p < 0.001 for proportion with zero WBC in CSF, with high levels in Maridi). Followup rates were low everywhere (Table 1) with the exception of Arua. Of the 6663 patients eligible, 2250 (33.8%) attended one followup visit, 580 (8.7%) attended two, 232 (3.5%) three, 114 (1.7%) four, and two (0.03%) five visits. Overall, 47.7% of suspects returned for followup at least once. Attendance declined with time (however, the low followup rates at months 9 and 12 partly reflect the fact that by then many patients had been released from followup after becoming CATT negative). Among patients seen once, persontime under observation was roughly comparable across sites.
Table 1. Baseline characteristics and rates of followup of stage 1 serological suspects, by site
Marked differences in patient endpoints were evident across sites (Table 2). Most persons lost their strong CATT reactivity during followup (noncases), about a third remained suspects, while parasitologically confirmed HAT was noted in only a minority. Overall, 199 suspects were confirmed as stage 1 during followup, and 99 progressed to confirmed stage 2 (Table 2), leaving 298 "true" HAT cases among whom survival analysis was conducted.
Table 2. Endpoints during followup among stage 1 serological suspects who attended at least one control visit, by site
The 298 patients retained differed systematically from the remaining 6365 stage 1 suspects who were not confirmed as stage 1 or 2 during followup, and thus excluded: they were slightly younger (mean age 23.4 years versus 25.0 years, p = 0.029), a lower proportion had a WBC count of zero in their CSF (109 [36.6%] versus 2801 [44.0%], p = 0.011), and a lower proportion were screened actively (132 [44.3%] versus 3427 [53.8%], p = 0.001). Patients followed up at least once also differed systematically from those never seen again: more were female (1528/3178 or 48.1% versus 1565/3485 or 44.9%; p = 0.01); more had a zeroWBC count in CSF (1438/3178 or 45.2% versus 1472/3485 or 42.2%; p = 0.01); and considerably more were positive at a CATT dilution >1:4 (2165/2713 or 79.8% versus 1468/2305 or 63.7%; p < 0.001).
Stage 1 duration
After maximum likelihood estimation, a simple exponential model fit the data no worse (p > 0.20 for all comparisons among nested models) than any more complex alternatives (Weibull, loglogistic, Gamma, and Gompertz), as evidenced visually (Figure 3), and was thus retained. Because distribution of survival times post detection is exponential, the overall function for survival in stage 1 must also be exponential and have the same parameter (Figure 2). The model thus yields the following function for the proportion P of cases remaining in stage 1 up to time t:
Figure 3. Fitted intervalcensored stage 1 survival models. For comparison purposes, the figure also superimposes a KaplanMeier survival curve (with Greenwood 95% confidence intervals) based on censoring progression events at the midpoint between the two visits between which the event is known to have occurred.
where r_{1}, the daily rate of progression to stage 2, is 0.0019 (95%CI 0.0012 to 0.0028). Accordingly, the mean duration of stage 1 is 1/r_{1}, or 526 days (95%CI 357 to 833), corresponding to about one and a half years. The median duration is ln(2)/r_{1}, or 364 days (95%CI 248 to 578).
The estimate of r_{1 }is higher for patients detected though passive screening at the HAT centre (0.0021, 95%CI 0.0012 to 0.0038) than for patients detected actively in the community (0.0016, 95%CI 0.0013 to 0.0020). However, the difference is nonsignificant.
Stage 2 duration
Overall, 88 active screening sessions were eligible for analysis (Table 3). Infection prevalences were comparable in Adjumani and Kiri, but considerably lower in Arua. Nevertheless, the ratio of stage 1 and 2 prevalence was similar across sites, and the overall average was close to unity (1.04) considering the ratio of means, while a ratio of medians yielded 1.25.
Table 3. Details of eligible active screening sessions and stage 1 and 2 prevalences
Using Equation 2, the stage 2 removal rate r_{2 }was thus very close to r_{1}: 0.0020 per day (95%CI 0.0013 to 0.0029). Accordingly, the mean duration of stage 2 is 1/0.0020 or 500 days (95%CI 345 to 769), and the mean total duration of stage 1 and 2 combined is 1026 days (95%CI 702 to 1602), or approximately 2 years and 10 months. Similarly, the median stage 2 duration is 347 days (95%CI 239 to 533), giving a total median gambiense HAT duration of 711 days (95%CI 487 to 1111), or almost two years.
Sensitivity analysis to compare S_{1}/S_{2 }ratios according to the prescreening passive detection rate did not show significant differences (p = 0.41 for detection rate <1 vs ≥1 per 1000 personmonths, KruskalWallis test).
Discussion
Our estimate of the "natural" duration of gambiense HAT (almost three years, evenly split between stage 1 and 2) helps to refine commonly held views that both stages last 'months to years', and is broadly consistent with ranges provided by experts: "several months to two years" for stage 1 [7], and, for stage 2, "four months to one year" [7] or "from one to three years" [3]. Fèvre et al. did survival analysis on 96 gambiense patients reported in the historical literature, and estimated a median duration of three years [18].
To our knowledge this is the first attempt to quantify the duration of stage 1 and 2 specifically. Furthermore, the distribution of stage 1 survival post detection has a clear exponential shape. As shown in Figure 2, this observation is most consistent with a constant, i.e. timeindependent rate of progression, irrespective of incidence trend before detection (scenarios A1, A2 and A3), rendering our estimates entirely robust in this respect. Distributions approximating an exponential could also arise if incidence is stable, and the progression rate increases with time (scenario B1), or if incidence is decreasing (B3 or C3). However, the latter is implausible: in all sites, MSF intervened because incidence was credibly reported to be increasing. Our findings are least consistent with a timedependent progression rate in a scenario of increasing incidence, which would yield strongly convex distributions (B2 and C2).
The mathematical property of constant r_{1 }is crucial to interpreting observed variation in the virulence of HAT: exponential distributions allow for a tail of unusually long durations, reports of which abound in the literature (Checchi et al., submitted). For example, using our estimated survival function, 3% of cases would be expected to remain in stage 1 for >5 years. Such cases should probably no longer be considered extraordinary, but simply a result of the intrinsic pattern of disease variability. The implication for control is that a small but important proportion of cases can carry the infection for many years, and, if undetected by active screening rounds, could maintain a tiny but persistent infectious reservoir even when local elimination appears all but achieved. Mopup active screening rounds, even years after transmission in a HAT focus is considered to have been brought under control, could thus play a vital role in reducing the risk of epidemic resurgence, something passive case detection might by itself never achieve.
If stage 2 invariably progresses to death, as is widely assumed, then r_{2 }is the specific HAT death rate after stage 2 initiation, and 1/r_{2 }is the life expectancy of stage 2. To our knowledge only one other study has attempted to estimate this parameter: Jusot et al. reanalysed observations by Greggio on the survival of absconded HAT patients diagnosed at a Belgian Congo hospital between 1907 and 1915 [19]. They showed an exponential decay of survival, and estimated a daily death rate of 0.00235 (mean duration = 425 days), very close to our estimate.
Potential Biases
Our findings are subject to six main potential sources of bias, which mainly affect the estimates of stage duration:
1. A systematic delay in ascertainment of stage 2 progression, leading to an underestimation of stage 1 duration, could have occurred if patients had presented to the treatment centre only when CNS involvement had become symptomatic. Interval censoring removes part of this potential bias by considering only the interval during which the event occurred. Furthermore, only 29.3% (1458/4971) of followup visits in our entire cohort of 6663 patients occurred >30 days earlier or later than scheduled, suggesting that event ascertainment was mainly driven by the followup schedule, not symptom onset. There was a mean delay of 34 days between the scheduled date of followup and the date when the visit actually took place. This delay was greater when patients were found to have progressed to stage 2 at the visit (56 days versus 34 days for all other visits; p = 001 for difference of means; p = 0.070 for difference of medians), but this is a small difference considering HAT evolves over many months. The distribution of visit timing with respect to the scheduled date was similar whether or not progression to stage 2 was detected at the visit (Figure 4).
Figure 4. Timing of actual followup visits with respect to the scheduled date. Data are provided for visits at which progression to stage 2 was detected (n = 99), as well as all other visits (n = 4774).
2. About half of analysiseligible patients were followed up, each for an average of five months. At baseline, this group had a higher prevalence of strong CATT reactivity than defaulters, suggesting it harboured more true HAT cases. This difference should have been eliminated by virtue of including only confirmed stage 1 or 2 patients in the survival analysis. Stage 1 duration may nonetheless have been underestimated if patients progressing to stage 2 and thus feeling ill were more likely to present for followup; however, this pattern may be unlikely given that progression events occurred around scheduled visits. It is also possible that patients progressing may have stayed away from MSF HAT centres due to a perception that the programme had failed to treat them: in this case, duration would be overestimated.
3. Our analysis assumes that suspects have the same rate of progression to stage 2 as 'typical' confirmed stage 1 patients. In fact, CATT serological suspects are a complex population [20], probably composed of a mixture of false parasitological negatives and false serological positives. False negatives would be due to (i) testing during the prepatent phase (first 1–2 weeks after the tsetse bite); (ii) lowvirulence cases featuring a parasite density below the detection threshold of microscopy; (iii) testing when parasite density is at its minimum, as part of its welldescribed undulating pattern, and thus below the detection threshold; or (iv) laboratory error. False positives would result from (i) CATT crossreaction with other trypanosome species (some of which are known to cause transient and possibly even pathogenic infections), other parasites such as leishmania and Plasmodium spp., or rheumatoid factors [21]; (ii) past HAT infection (antibodies can survive for up to 3 years) [22]; or (iii) error in CATT performance. So as to remove as many of the false positives as possible and isolate a population most likely to comprise typical pathogenic infections, we restricted the denominator at risk to suspect cases that were indeed true HAT cases, as evidenced by stage 1 or 2 parasitological confirmation during followup. Arguably, we could have retained in the analysis individuals who remained suspects until the end of followup, and the survival function would have been markedly different, yielding a far longer stage 1 duration (r_{1 }= 0.00014, or about 20 years), which does not seem plausible. However, the vast majority of these cases have been shown in other studies to be false positives by both parasitological and PCR techniques [20,2325]. Were suspects progressing more slowly than other cases, the true stage 1 duration would be shorter than our estimate: however, any difference is unlikely to be great, since a much faster r_{1 }than the one we measured would yield durations of only a few months, consistent with empirical observations of rhodesiense, but not gambiense HAT, for which one to three years are the norm (Checchi et al, submitted).
4. Our estimate of r_{2 }should be regarded as less robust and potentially more biased than that of r_{1}, since it relies on an assumption of stable incidence, which is atypical in HAT (indeed, past models of HAT transmission have suggested unstable incidence and lack of equilibrium might be an inherent feature of HAT dynamics) [19,26,27]. A rising incidence, which is more plausible, would skew the S_{1}/S_{2 }ratio upwards, and result in an overestimation of r_{2}, i.e. underestimated stage 2 duration. This bias is probably not large, since all foci in this analysis featured relatively low prevalence (<2–3%), suggesting incidence may have been rising slowly. However, our deduced stage 2 duration should be probably regarded as a lowerend estimate.
5. So as to assemble a sufficiently large cohort of patients, we combined datasets from various HAT foci. We assumed that this cohort would be representative of the population of HAT patients in these foci. However, HAT foci might differ in several ways. Particularly, incidence trends might have been divergent (e.g. increasing in one site and decreasing in another), affecting the stable incidence assumption differently; and the profile of serological suspects, including the prevalence of true HAT stage 1 cases (see point 3 above) might vary due to heterogeneity in incidence, coverage and timeliness of case detection or diagnostic sensitivity across sites: unequal representation of the different foci in the overall cohort, due merely to the duration and coverage of the respective treatment programmes and to incidence during the programme durations, would bias results towards the patterns in foci that are overrepresented in the cohort.
6. Estimates here apply to strains of T. gambiense present in Sudan and Uganda. More or less virulent strains encountered elsewhere might feature different progression and removal rates.
Conclusion
Our quantitative understanding of the dynamics of HAT is poor relative to that of many other infectious diseases. This is partly a consequence of imprecise knowledge of the fundamental parameters that govern HAT disease dynamics. Here we have estimated two of these key parameters based on substantial datasets, providing estimates of the true duration of untreated infections. These parameter estimates are essential prerequisites to the development of formal approaches for the quantitative evaluation of different strategies to control this neglected disease. The duration of infection and thus infectiousness is a fundamental determinant of HAT's reproductive ratio, and present control strategies essentially work by reducing r_{1 }and r_{2}: mathematical models will only offer reliable predictions of the impact of different screening strategies if they incorporate realistic values for r_{1 }and r_{2 }[19]. Furthermore, the exponential property of stage 1 to stage 2 progression determines the distribution of asymptomatic infection durations, and therefore the probability that some longduration cases will fall through the net of case detection before they become symptomatic, and seed renewed HAT outbreaks once control is relaxed. These quantitative dynamic processes should inform any future policies to control and possibly eliminate gambiense HAT.
Abbreviations
CATT, Card Agglutination Test for Trypanosomiasis; CNS, Central nervous system; CSF, Cerebrospinal fluid; HAT, Human African trypanosomiasis; IQR, Interquartile range; MSF, Médecins Sans Frontières; PCR, Polymerase Chain Reaction; S_{1}, point prevalence of HAT cases in stage 1 of the disease; S_{2}, point prevalence of HAT cases in stage 2 of the disease; λ, incidence rate; r_{1}, rate of progression from stage 1 to stage 2; r_{2}, rate of removal from stage 2 (through death or spontaneous cure); P(t), proportion of HAT cases still in stage 1 at time t after infection; t, time; WBC, white blood cells.
Competing interests
The author(s) declare that they have no competing interests.
Authors' contributions
F Checchi designed the study, carried out analyses and wrote this paper. J Filipe, D Haydon and D Chandramohan assisted in data analysis and interpretation. F Chappuis was involved in study design and interpretation. All authors made substantial contributions to the manuscript.
Acknowledgements
We are grateful to MSF operational sections for sharing data used for this analysis, and to the Médecins Sans Frontières HAT Working Group for encouragement and support. Thanks in particular to Gerardo Priotto (Epicentre) and the Epicentre team in Uganda for help with data collection. Veerle Lejon provided advice on stage 1 and 2 definitions. This work was conducted without any specific source of funding.
References

Stich A, Barrett MP, Krishna S: Waking up to sleeping sickness.
Trends Parasitol 2003, 19(5):195197. PubMed Abstract  Publisher Full Text

Pepin J, Meda HA: The epidemiology and control of human African trypanosomiasis.
Adv Parasitol 2001, 49:71132. PubMed Abstract

Apted FIC: Clinical manifestations and diagnosis of sleeping sickness. In The African trypanosomiases. first edition. Edited by Mulligan HW. London: George Allen and Unwin Ltd; 1970:661682.

World Health Organization: Human African trypanosomiasis (sleeping sickness): epidemiological update.

Kennedy PG: Diagnostic and neuropathogenesis issues in human African trypanosomiasis.
Int J Parasitol 2006, 36(5):505512. PubMed Abstract  Publisher Full Text

Frezil JL, Carnevale P: Le probleme du reservoir du virus et du mantien des foyers de trypanosomiase humaine en Afrique Centrale.

Burri C, Brun R: Human African trypanosomiasis. In Manson's tropical diseases. 21st edition. Edited by Cook GC, Zumla A. London: Elsevier Science limited; 2003:13031323.

Murray M, Morrison WI, Whitelaw DD: Host susceptibility to African trypanosomiasis: trypanotolerance.
Adv Parasitol 1982, 21:168. PubMed Abstract

Sternberg JM: Human African trypanosomiasis: clinical presentation and immune response.
Parasite Immunol 2004, 26(11–12):469476. PubMed Abstract  Publisher Full Text

Gouteux JP, Artzrouni M: [Persistence and resurgence of sleeping sickness caused by Trypanosoma brucei gambiense in historic foci. Biomathematical approach of an epidemiologic enigma].
C R Acad Sci III 2000, 323(4):351364. PubMed Abstract  Publisher Full Text

Odiit M, Coleman PG, Liu WC, McDermott JJ, Fevre EM, Welburn SC, Woolhouse ME: Quantifying the level of underdetection of Trypanosoma brucei rhodesiense sleeping sickness cases.
Trop Med Int Health 2005, 10(9):840849. PubMed Abstract  Publisher Full Text

Odiit M, Kansiime F, Enyaru JC: Duration of symptoms and case fatality of sleeping sickness caused by Trypanosoma brucei rhodesiense in Tororo, Uganda.
East Afr Med J 1997, 74(12):792795. PubMed Abstract

Paquet C, Castilla J, Mbulamberi D, Beaulieu MF, Gastellu Etchegorry MG, Moren A: [Trypanosomiasis from Trypanosoma brucei gambiense in the center of northwest Uganda. Evaluation of 5 years of control (1987–1991)].
Bull Soc Pathol Exot 1995, 88(1):3841. PubMed Abstract

Chappuis F, Udayraj N, Stietenroth K, Meussen A, Bovier PA: Eflornithine is safer than melarsoprol for the treatment of secondstage Trypanosoma brucei gambiense human African trypanosomiasis.
Clin Infect Dis 2005, 41(5):748751. PubMed Abstract  Publisher Full Text

Chappuis F, Loutan L, Simarro P, Lejon V, Buscher P: Options for field diagnosis of human african trypanosomiasis.
Clin Microbiol Rev 2005, 18(1):133146. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Williamson JM, Satten GA, Hanson JA, Weinstock H, Datta S: Analysis of dynamic cohort data.
Am J Epidemiol 2001, 154(4):366372. PubMed Abstract  Publisher Full Text

INTCENS: Stata module to perform intervalcensored survival analysis [http://ideas.repec.org/c/boc/bocode/s453501.html] webcite

Fevre EM, Coleman PG, Welburn SC, Maudlin I: Reanalyzing the 1900–1920 sleeping sickness epidemic in Uganda.
Emerg Infect Dis 2004, 10(4):567573. PubMed Abstract  Publisher Full Text

Jusot JF, de Vlas SJ, van Oortmarssen GJ, De Muynck A: [Contribution of a mathematical model in the control of a parasitosis: the case of human African trypanosomiasis due to Trypanosoma brucei gambiense].
Ann Soc Belg Med Trop 1995, 75(4):257272. PubMed Abstract

Garcia A, Jamonneau V, Magnus E, Laveissiere C, Lejon V, N'Guessan P, N'Dri L, Van Meirvenne N, Buscher P: Followup of Card Agglutination Trypanosomiasis Test (CATT) positive but apparently aparasitaemic individuals in Cote d'Ivoire: evidence for a complex and heterogeneous population.
Trop Med Int Health 2000, 5(11):786793. PubMed Abstract  Publisher Full Text

Paquet C, Ancelle T, GastelluEtchegorry M, Castilla J, Harndt I: Persistence of antibodies to Trypanosoma brucei gambiense after treatment of human trypanosomiasis in Uganda.
Lancet 1992, 340(8813):250. PubMed Abstract  Publisher Full Text

Magnus E, Vervoort T, Van Meirvenne N: A cardagglutination test with stained trypanosomes (C.A.T.T.) for the serological diagnosis of T. B. gambiense trypanosomiasis.
Ann Soc Belg Med Trop 1978, 58(3):169176. PubMed Abstract

Chappuis F, Stivanello E, Adams K, Kidane S, Pittet A, Bovier PA: Card agglutination test for trypanosomiasis (CATT) enddilution titer and cerebrospinal fluid cell count as predictors of human African Trypanosomiasis (Trypanosoma brucei gambiense) among serologically suspected individuals in southern Sudan.
Am J Trop Med Hyg 2004, 71(3):313317. PubMed Abstract  Publisher Full Text

Koffi M, Solano P, Denizot M, Courtin D, Garcia A, Lejon V, Buscher P, Cuny G, Jamonneau V: Aparasitemic serological suspects in Trypanosoma brucei gambiense human African trypanosomiasis: a potential human reservoir of parasites?
Acta Trop 2006, 98(2):183188. PubMed Abstract  Publisher Full Text

Simarro PP, Ruiz JA, Franco JR, Josenando T: Attitude towards CATTpositive individuals without parasitological confirmation in the African Trypanosomiasis (T.b. gambiense) focus of Quicama (Angola).
Trop Med Int Health 1999, 4(12):858861. PubMed Abstract  Publisher Full Text

Rogers DJ: A general model for the African trypanosomiases.
Parasitology 1988, 97(Pt 1):193212. PubMed Abstract

Artzrouni M, Gouteux JP: Control strategies for sleeping sickness in Central Africa: a modelbased approach.
Trop Med Int Health 1996, 1(6):753764. PubMed Abstract
Prepublication history
The prepublication history for this paper can be accessed here: