Introduction
Animal African Trypanosomiasis continues to be a serious veterinary problem, costing 37 African countries USD 4.5 billion annually (FAO, 2017). The trypanosome species of veterinary interest in eastern and southern Africa are Trypanosoma vivax and T. congolense. Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) used a catalytic function to model age-related changes in the prevalence of T. vivax infections in female G. pallidipes. For T. congolense, they suggested a slightly more complex model that allowed for a reduction in prevalence among older flies. Age was estimated using ovarian dissection, which can provide precise estimates of chronological age for flies in their first ovarian cycle – i.e., flies that have ovulated fewer than four times – but not, in general, for older flies (Challier, Reference Challier1965). Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) assumed for their sample, however, that negligible numbers of flies survived to a third ovarian cycle so that uncorrected ovarian category is identical to ovarian age – i.e., the number of times the fly has ovulated – and is thus a valid index of chronological age. This assumption was carried over into other modelling of the same data (Woolhouse and Hargrove, Reference Woolhouse and Hargrove1998; Lord et al., Reference Lord, Barnard, Day, Hargrove, McNamara, Paul, Trenholme and Woolhouse1999). The larger numbers dissected in the current study makes clear that the above assumption cannot be sustained. Accordingly, we develop improved models that account for the longer survival of adult flies.
In so doing, we also investigated the suggestion that the probability of a tsetse acquiring a trypanosome infection is higher when the fly takes its first blood meal – i.e., as a teneral fly – than at all subsequent feeds. This so-called ‘teneral effect’ is well established for T. brucei, being first described by Van Hoof et al. (Reference Van Hoof, Henrard and Peel1937) for the acquisition of T. b. gambiense by G. palpalis. Haines (Reference Haines2013) provides an excellent review of the subject and notes that refractoriness to T. brucei infections in non-teneral tsetse can be reduced among starved flies (Kubi et al., Reference Kubi, Van den Abbeele, De Deken, Marcotty, Dorny and Van den Bossche2006). In his general model for the African trypanosomiases, Rogers (Reference Rogers1988) assumed that there was no such teneral effect in the case of T. congolense, which causes disease in livestock, but not humans. Distelmans et al. (Reference Distelmans, D'Haeseleer, Kaufman and Rousseeuw1982) had, however, found that teneral G. p. palpalis were more readily infected with T. congolense than older flies. Later, Welburn and Maudlin (Reference Welburn and Maudlin1992) using laboratory G. m. morsitans, found the risk of acquiring a T. congolense infection was seven-times higher in tenerals than in all older flies. This led them to suggest that non-teneral flies do not play a significant part in the epidemiology of T. congolense. The motivation for our paper is thus two-fold: (i) to provide improved models for the dynamics of trypanosome infection in the field; (ii) to assess the importance of the teneral effect in the dynamics of infection of G. pallidipes females with T. congolense.
Methods
Study area
The study was carried out between September 1990 and April 1999 at Rekomitjie Research Station (16°10′S, 29° 25′E, altitude 520 m), Zambezi Valley, Zimbabwe. Daily maximum and minimum temperatures were recorded using mercury thermometers housed in a Stevenson screen at the Station, and rainfall was recorded using a gauge placed 4 m from the Stevenson screen.
Fly sampling
Female G. pallidipes were collected at various sites within 2 km of the Research Station; >99% of flies were captured in epsilon traps (Hargrove & Langley, Reference Hargrove and Langley1990) baited with artificial host-odour consisting of acetone, 1-octen-3-ol, 3-n-propyl phenol and 4-methyl phenol released at ~200, 0.4, 0.01 and 0.8 mg h−1, respectively. The cages in which flies were captured were wrapped in a heavy-duty black cotton cloth and kept in a polystyrene box; in about 90% of cases, traps were cleared at 30-minute intervals. Flies were dissected within 24 hours of capture.
Classification of flies by ovarian and wing fray category and trypanosome infection status
We followed the general dissection procedure described by Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993). We dissected the ovaries and uterus of female tsetse and assigned each to an ovarian category. The general technique was originally developed for estimating the age of female mosquitoes (Detinova, Reference Detinova1949, Reference Detinova1959, Reference Detinova1962). It was then applied to tsetse by Saunders (Reference Saunders1960, Reference Saunders1962) and improved by Challier (Reference Challier1965). Details of the age estimation procedure, as applied to female tsetse, are provided by the above authors, and the technique has been used and cited by many others since then (Rogers et al., Reference Rogers, Randolph and Kuzoe1984; Dransfield et al., Reference Dransfield, Brightwell, Kiilu, Chaudhury and Adabie1989; Randolph et al., Reference Randolph, Rogers and Kiilu1990; Hargrove, Reference Hargrove, Howard-Kitto, Kelleher and Ramesh1993, Reference Hargrove1995, Reference Hargrove1999a, Reference Hargrove1999b; Reference Hargrove2012; Hargrove et al., Reference Hargrove, English, Torr, Lord, Haines, van Schalkwyk, Patterson and Vale2019). Nonetheless, for convenience – and particularly given the complexities of the technique – we provide the following details of the age estimation procedure, based largely on the description provided by Hargrove (Reference Hargrove2012).
Figure 1 illustrates the disposition of the four oocytes within the paired ovaries at different stages of the life of a female tsetse fly. In the normal fly, oocytes develop and are ovulated in the strict order: right inner (labelled A), left inner (C), right outer (B) and left outer (D). Ovulation of all four oocytes comprises a single ovarian cycle. The relative size of the oocytes allows distinction between flies that have ovulated zero, one, two or three times. The oocytes of a fly that has ovulated four times have, however, the same conformation as a fly that has not yet ovulated (fig. 1). Flies that have ovulated 0 and 4 times can, however, be distinguished because ovulation results in the ovariole bearing a follicular relic of the previous ovulation. A fly that has ovulated four times has a relic on each ovariole, whereas a fly that has not yet ovulated for the first time has no relic on any ovariole. Similarly, since flies in ovarian categories 1, 2 and 3 have 1, 2 or 3 relics, respectively, they can be distinguished from flies that have ovulated 5, 6 or 7 times – all of which have four relics.
Diagrammatic representation of changes in the appearance of tsetse fly ovaries during successive ovulations. Symbols next to each diagram indicate the ovarian category: 0, 1, 2, 3, 4 + 4n, 5 + 4n, 6 + 4n, 7 + 4n. 0a and 0b indicate females that have not yet ovulated for the first time and where the largest oocyte is shorter, or longer, respectively, than about 0.6 mm.
Problems arise for flies that have ovulated more than seven times because, in tsetse – though not in mosquitoes (Detinova, Reference Detinova1962) – only one relic is ever seen on an ovariole. Further ovulations do not result in changes in the appearance of the relics. Since it is not therefore possible to distinguish a fly that has ovulated 4 times from one that has ovulated 8, 12, 16 … , etc., times, we define ovarian category 4 + 4n, n = 0, 1, 2, 3… as a category that includes flies of all these ovarian age groups. Ovarian categories 5 + 4n, 6 + 4n and 7 + 4n are similarly defined. In what follows we distinguish between a fly's ovarian category – as defined above – and its ovarian age group, being the number of times the fly has ovulated.
For flies in ovarian categories 0 to 3, only, the ovarian category and ovarian age group are identical. In principle n has no upper bound but the relationship between wing fray and age in G. pallidipes females suggests that <1% of G. pallidipes females survived to ovulate >11 times in the field at Rekomitjie (Hargrove, Reference Hargrove2020). Negligible numbers of females would thus survive long enough to ovulate more than 15 times, and the same conclusion can be reached from field estimates of the age-specific mortality of G. m. morsitans (Hargrove et al., Reference Hargrove, Ouifki and Ameh2011). These results suggest n takes a maximum value of 3 for our study situation and we have assumed this is the case. Flies are thus assumed to survive a maximum of four ovarian cycles; ovarian category 4 + 4n includes flies in ovarian age groups 4, 8 and 12 – and analogous statements apply to flies in ovarian categories 5 + 4n, 6 + 4n and 7 + 4n.
We used the classification of Jackson (Reference Jackson1946) to assign each fly to one of six wing fray categories. While ovarian dissection provides a more accurate estimate of a fly's chronological age for tsetse, this is only true for flies in ovarian age group <4. Wing fray, by contrast, provides a relative measure of age for all flies and allows us to gauge how trypanosome prevalences change over flies’ entire lives.
The mouthparts of each fly were dissected, and the labrum and hypopharynx were examined under a light microscope at 240 × magnification. If no trypanosomes were found in the mouthparts the fly was diagnosed as negative for all trypanosome infections and no further dissection was carried out. If trypanosomes were detected, the midgut was dissected out and screened for infection. Where trypanosomes were detected in both the midgut and mouthparts, the fly was diagnosed as being infected with T. congolense; if trypanosomes were found only in the mouthparts the fly was diagnosed as being infected with T. vivax. It is also possible for a fly to be infected with both trypanosome species.
The delay between the time when a fly becomes infected and when it is possible to detect trypanosomes in the fly by dissection is such that there is effectively zero probability of finding trypanosomes in a fly before it has ovulated for the first time. Thus, among 1228 flies dissected that were adjudged to have recently ovulated for the first time and were either 8 or 9 days old – 0 and 2 were found to harbour a T. vivax or a T. congolense infection, respectively. The estimated prevalences, with 95% confidence interval, were thus extremely close to zero: 0 (0–0.003) and 0.0016 (0.0002–0.006), respectively. Flies in ovarian category 0, which are by assumption younger than 8 days, will have expected prevalences even lower than the above values. Accordingly, flies found to be in ovarian category 0 were not dissected further to look for trypanosomes, and the trypanosome prevalence of such flies was assumed to be identically zero.
Modelling patterns of acquisition of trypanosome infections
Catalytic model: assuming no flies ovulate more than seven times
Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) used a catalytic model to explain age-related changes in the prevalence of T. vivax infections in female G. pallidipes, fitting their data using the function:
In equation 1, p(k) is the predicted trypanosome prevalence among flies of age k days. Equivalently p(k) can be interpreted as the probability that any fly will become infected by the time it reaches the age of k days. The parameter λ is the rate of infection and τ is the number of days between the acquisition of an infection and the time when trypanosomes can first be detected in the fly. Equation 1 predicts prevalence directly for flies of any known daily age, without requiring knowledge of fly mortality, and without using explicit predicted values for the numbers of infected and uninfected flies. We follow Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) in applying our models to data pooled over time, so that model parameter estimates are viewed as long-term averages.
To help interpret λ, one can replace the exponential function in equation 1 with its power series expansion. For fractional values of λ(k- τ) close to zero, as in our estimates below, the squared and higher-power terms of the power series are negligible in magnitude. Thus, equation 1 can be closely approximated by p(k) ≈ λ(k - τ). Since p(k + 1) –p(k) ≈ λ, the λ parameter of equation 1 can be interpreted as the increase in a fly's infection probability from any one day to the next. That is, λ is the approximate probability that any uninfected fly will become infected within the next 24 hours.
We applied the catalytic model to our data. Each dissected fly in our data was identified only to its ovarian category; to apply equation 1, we specified the range of daily ages spanned by each category. Following Hargrove (Reference Hargrove1995, Reference Hargrove2012) we allocate flies of age 0–7 days to category zero: thereafter we assume a 9-day pregnancy duration, so that successive ovarian age groups include successive 9-day intervals. For the catalytic model only, we follow Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) in assuming that no flies survived to ovulate more than seven times so that, for this model only, ovarian age group i and ovarian category i are identical for i = 1,7. Thus, ovarian category 1 includes days 8 to 16, category 2 includes days 17 to 25, and so forth. In general, category i includes the daily ages (i*9-1) to (i*9 + 7), with all lifespans terminating at (7*9 + 7) = 70 days.
We then assume that the predicted prevalence P(i), for flies of ovarian category i, i = 1,7 is equal to the mean of the daily prevalences predicted for all 9 days included in that ovarian category. That is,
As explained above, we assume P(0) = 0. In words, equation (2) says that we find the average prevalence for any given ovarian category by summing the prevalences for each day of pregnancy and dividing by 9.
Susceptible-infected SI model: assuming flies can survive to ovulate more than seven times
When fitting the catalytic model, no account was taken of fly survival – nor, explicitly, of the numbers of susceptible and infected flies found in each ovarian category. Instead, trypanosome observed and predicted prevalences were compared directly. As discussed earlier, however, it is clear that significant proportions of female tsetse do survive to ovulate more than seven times. Under these circumstances we need to take fly mortality into consideration and the catalytic model of Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) is not appropriate. Instead, we used an SI model to predict, separately, the numbers of infected (I), noninfected (susceptible, S), and total (N = S + I) flies, that would occur at each daily age, k. We assume that, as flies age, they progressively transition from the susceptible to the infected group, with equation 1 again predicting the probability that a surviving fly will become infected by the age of k. We pool the predicted daily values of S and I into their respective ovarian age groups and for, age groups >3, further pool these for ovarian categories 4 + 4n to 7 + 4n.
To model fly mortality, we assume that flies of all ages, and infection status, die at the instantaneous daily rate μ, so that an initial, total number, N(0) of newly-emerged flies will have decreased to N(k) = N(0) exp(- μk) by age k. The daily survival probability of any fly is φ = exp(- μ). With this mortality model, the numbers of surviving flies that are infected is given by: I(k) = N(0) exp(-μk)(1 - exp (-λ(k - τ)), where λ and τ are as defined for equation 1. The total numbers, and numbers infected, at each daily age are again pooled into 9-day ovarian age groups, now indexed by i = 1,2, … ,15. As before, ovarian age group i includes the daily ages (i*9-1) to (i*9 + 7), but all flies are now assumed to have died by the age of (15*9 + 7) = 142 days. For ovarian age groups i = 1 to 3, and daily age indexed by k, we sum the numbers predicted for each daily age within each age group, to get group totals:
For the ovarian age groups 1 to 3, prevalence is then predicted as
For older flies, we need to further pool the ovarian age groups to calculate the prevalence for the corresponding ovarian category, as is shown for the (4 + 4n) category:
and analogous equations are used for flies in ovarian categories 5 + 4n, 6 + 4n and 7 + 4n. Note that N(0) can be an arbitrary value because it divides out when calculating P(i).
Susceptible, exposed and infective (SEI) model
The fixed time delay (τ) in the SI model is biologically unrealistic in assuming a sharp boundary between the ages when tsetse are and are not infective. As an alternative, we consider a model where we partition flies into three states – Susceptible, Exposed and Infective (SEI). Instead of a fixed delay (τ) we model the proportions of flies transitioning from the S to E to I compartments during each age increment, again indexed by k. Susceptible flies are those that have never acquired a trypanosome infection. Exposed flies have been so infected, but the trypanosome infection has not developed to the point where the fly can pass on that infection – i.e., it is not infective. Infective flies can pass on the trypanosome infection; for both the SI and SEI models we assume that flies, once infected, stay infected for life. They thus only leave the infective state by dying. They are also the only flies in which trypanosome can be detected using dissection techniques; accordingly, only the infective flies are included in compartment I in the modelling. Our SEI model assumes the same extended set of 16 ovarian age groups (including category/group 0) as were used in the SI model, and also assumes the same mortality model.
We assume an arbitrary initial number S0 of susceptible flies. Since none has yet been exposed or infected, E0 = 0 and I0 = 0. Then the following transition equations apply for the SEI model as the flies age:
where, for each unit of age φ is again the probability that a fly will survive to the next age, regardless of its infection status. The parameter λ is the probability that susceptibles (Sk) move to the exposed class (Ek +1). Since all flies surviving the Exposed class ultimately become infected, λ is again the probability of becoming infected in a single age increment. The γ parameter is the probability of moving from the exposed class (Ek) to the infective class (Ik +1), between ages k and k + 1. The total count is given by Nk = Sk + Ek + Ik. When implementing equations 5 we used the same, separately estimated, value of φ as used in the SI model.
When applying this model, we changed the age increment from 1 to 3 days, assuming that the latter increment approximated the interval between successive feeds, and that the flies take an average of three meals in a 9-day pregnancy (Randolph et al., Reference Randolph, Williams, Rogers and Connor1992). The change enabled us to investigate the suggestion that the force of infection is higher at the first blood meal than at subsequent meals (Welburn and Maudlin, Reference Welburn and Maudlin1992). This was achieved by allowing that λ in equation 5 could take a different (presumed higher) value when the fly took its first blood meal, compared with all later meals. Notice that the first meal is taken during ovarian category 0 when no flies can yet be in category I and when, therefore, the numbers play no role in the model fitting. If, however, λ were indeed very much higher at the first meal, we should expect a surge of infectives during ovarian category 1.
With a 3-day age increment, there are 3 increments per 9-day pregnancy. We now begin ovarian age group 1 at day 7, rather than day 8, so that age group 0 includes the 3-day age increments k = 1 and 2, while group 1 includes increments k = 3, 4, and 5. In general, ovarian age group i includes the age increments from k = (3*i) to k = (3*i + 2), for i = 1, 2, …15. Thus, predicted total values of S, E, I and N for each of the 15 ovarian age groups are given by:
Finally, the values of I(i) and N(i) from equation 6 are inserted into equations 3 and 4, to calculate the SEI- predicted prevalences for the younger ovarian categories i = 1, 2, 3 and for the older, pooled ovarian categories 4 + 4n to 7 + 4n.
Statistical methods
Prevalence models: parameter estimation and model assessment
We estimated the parameters of all models using maximum likelihood (ML), as implemented in the R package maxlik (Henningsen and Toomet, Reference Henningsen and Toomet2011). For the catalytic, SI and SEI models, the predicted prevalence can also be interpreted as the predicted probability that a randomly sampled fly is infected. Thus, the log likelihood for each of the three models has the binomial form:
where nI(i) and nU(i) = n(i) – nI(i) are the observed counts of infected and uninfected flies, respectively, in the 7 ovarian categories i = 1 to 7 + 4n. The observed total count in each category is n(i).
The three prevalence models differ only in how they predict the P(i) values of equation 7. For the catalytic model, P(i) is given by equation 2 and, for the SI model, the predicted P(i) values are given by equations 3 and 4. For the SEI model, predicted values of P(i) are again given by equations 3 and 4, but with the components of P(i) (that is, I(i) and N(i)) now calculated from equations 6.
We report 95% confidence intervals for all estimated parameters. In addition, we used the Akaike information criterion (AIC), as well as graphical displays, to assess and compare how well the prevalence models fit the data (Burnham and Anderson, Reference Burnham and Anderson2002). AIC measures a model's quality of fit, with a penalty for the number of estimated parameters needed to achieve that quality. In calculating AIC for the SI and SEI models, we also included φ in the count of estimated parameters. As a rule of thumb, models having AIC values within 2 units of each other do not differ in their quality of fit (Burnham and Anderson, Reference Burnham and Anderson2002).
Estimation of survival probability
We also used ML to estimate the survival probability φ separately from the prevalence models and based only on the total counts in the ovarian categories. This estimate was then used as a fixed parameter value when estimating the other parameters of the SI and SEI models. To estimate φ, we assumed that the total counts, N(i), represent a sample from a stationary age distribution of live flies. Then, the number surviving to ovarian age group i is given by N(i) = N(0)φi, where N(0) is an arbitrary initial count of age 0 flies, and φ is now the 9-day survivorship probability. Assuming that no flies survive beyond 15 ovulations, the grand total count of all flies is NT = N(0)∑i φi, where the sum ranges over i = 1,2 … 15. Then, the ratio N(i)/NT expresses the age distribution as a multinomial distribution of the probabilities p(i) that a randomly-selected fly is of ovarian age group i.
For ovarian categories i = 1, 2, 3, the probabilities are:
The number of survivors in ovarian category 4 + 4n includes all flies that survived through ovarian age groups 4, 8 or 12, and similarly for the other pooled categories. Thus the survival probabilities for the grouped ovarian categories 4 + 4n and 7 + 4n are:
with analagous expressions for 5+6n and 6+6n. Finally, the log likelihood is given by
where n(i) are the observed total counts in ovarian categories 1, 2, … 7 + 4n, and the probabilities p(i) are given by equations 8 and 9. We maximized equation 10 directly. Equivalently, Hargrove (Reference Hargrove, Howard-Kitto, Kelleher and Ramesh1993) set the derivative of Equation 10 to zero, thus yielding a closed form for the maximum, which was then solved numerically.
As defined above, φ is the probability of survival for 9 days. We converted the estimate of φ into 3-day and 1-day survival probabilities, for use in the SEI and SI models, by calculating its cube root and ninth root, respectively.
Results
Temporal variations in temperature and rainfall, and in trypanosome prevalence
Rekomitjie Research Station has a hot tropical climate with temperatures increasing between July and November, leading up to the start of the single rainy season, largely restricted to November-February (fig. 2A). A total of 29,360 female G. pallidipes were dissected for determination of ovarian category and trypanosome infection status. Of these, the wing fray category was gauged in 27,951 flies. T. vivax and T. congolense infection prevalences – for all flies pooled on age – were 3.45 and 2.66%, respectively and declined during each year as temperatures increased between about July and December. Thereafter they increased sharply for 1–2 months, then more slowly to peak again in July (fig. 2B).
(A) Minimum, mean and maximum daily temperatures and mean monthly rainfall 1991–1999, Rekomitjie Research Station, Zambezi Valley, Zimbabwe. (B) Monthly prevalence of infections of T. vivax and T. congolense in female G. pallidipes captured and dissected between 1990 and 1999. Error bars indicate the 95% confidence intervals. Samples sizes are for each month, pooled on year.
Age distribution of all flies dissected, and of those found infected with trypanosomes
Among all flies dissected, there were similar numbers in ovarian categories 1 and 2, smaller numbers in ovarian category 3, and then a marked jump in numbers when moving to ovarian category 4 + 4n (fig. 3A) – in accord with the fact that ovarian category 4 + 4n includes flies that have ovulated 4, 8 or 12 times (see above). For infected flies, numbers increase between ovarian categories 1 and 3 and then take an even larger jump for those in ovarian category 4 + 4n – reflecting the fact that flies can accrue infections throughout life (fig. 3B). The numbers of infected flies only begin to diminish for flies of ovarian category 4 + 4n, reflecting the changing balance between reduced numbers of susceptible flies, and thus new infections, and deaths of those that are already infected. When wing fray was used as a cruder measure of age, the numbers dissected increased until wing fray category 4, and the numbers of infected flies until category 4 or 5 (fig. 3C, D). The difference reflects the fact that wing fray cannot decrease with age – and will naturally increase at different rates in different flies.
Distribution by ovarian or wing fray category of female G. pallidipes captured at Rekomitjie Research Station 1990-1999. Figures show numbers (n) dissected and numbers (i) infected. A., C. T. vivax; B., D. T. congolense by ovarian and wing fray categories, respectively. Notice the difference in scales on the left and right vertical axes on all graphs.
Trypanosome prevalence as a function of fly age and ambient mean temperature
Stepwise logistic regressions were carried out, using the prevalences of either T. vivax or T. congolense infections as the dependent variable. Ovarian, and wing fray, categories and mean ambient temperatures all accounted for significant fractions of the variance for both dependent variables. For both species of trypanosome, ovarian category removed the largest proportion of the variance, and the odds ratios increased steadily with increasing ovarian category (table 1). Once this variation was removed, the odds ratios for increases with increasing wing fray category were more modest.
Results of logistic regression of the prevalence of T. vivax (A) and T. congolense (B) in female G. pallidipes, as a function of the fly's ovarian, and wing fray, category and the mean daily temperature (tbar91) over the nine days prior to its capture. Sample size 27,951
Trypanosome prevalences, adjusted for ovarian and wing fray categories, declined significantly with increasing temperature (table 1) measured as the mean temperature averaged over the nine days prior to the day on which the fly was captured – this period approximating the duration of pregnancy. The result is in accord with the finding that trypanosome prevalences, pooled on age, generally decrease in the hottest months of the year (fig. 2B).
Fitting data using the catalytic function
The model defined in equation 1, based on the assumption that no flies ovulate more than seven times, provides a reasonable fit to the data for both species of trypanosome (fig. 4) – although there are clear trends in the residuals. For T. vivax, the best fit was obtained with λ = 0.00141 per day and τ = 9.7 days (table 2), very close to the values of λ = 0.00149 per day and τ = 14.2 days of Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993), using only data for flies dissected in 1990 and 1991. For T. congolense, using the same model, the values were λ = 0.00107 per day and τ = 9.8 days (table 2). There is considerable overlap of the 95% confidence intervals between species, for λ and for τ. Results obtained using our larger data set also suggest that there is no need to postulate a decline in T. congolense prevalence among older flies (fig. 4B, cf. Woolhouse et al., Reference Woolhouse, Hargrove and McNamara1993).
Fitting the function p(k) = 1 – exp [-λ(k - τ)] to ovarian age-specific prevalences for (A) T. vivax and (B) T. congolense in female G. pallidipes – where k is the daily age of the fly, λ the rate of infection and τ the delay between the time that a fly is infected and when trypanosomes can first be detected in the fly. Predicted daily prevalences are averaged within each ovarian age group to estimate the prevalence of each ovarian age. It is assumed that no flies survive >7 ovulations. Error bars on observed prevalences are exact 95% confidence intervals for a binomial proportion.
Maximum likelihood estimates for catalytic, SI, and SEI infection models, for T. vivax and T. congolense infections of female G. pallidipes
For the SEI(2) model the rate of infection (λ1) at the first meal can be different from that at all subsequent meals. For all other models λ is identical for all ages. Units are daily for all models. Figures in parentheses adjacent to the estimate provide the 95% confidence interval (CI). We assume a fixed survival probability of 0.976/day, previously estimated via ML, from total count data only. If a CI spans 0, then that parameter is not significantly different from 0, when that parameter is added to a model already containing the other parameters. Lower AIC values within each species denote higher quality of model fit.
* Akaike information criterion.
Examination of the raw data on which the prevalences were calculated shows, however, that the Woolhouse et al. (Reference Woolhouse, Hargrove and McNamara1993) assumption that a negligible fraction of flies survives to a third ovarian cycle cannot be sustained (fig. 3). The increase, between ovarian categories 3 and 4 + 4n, in the total numbers of flies caught and dissected – and the large jump in the numbers of infected flies – make it obvious that significant proportions of flies do survive to a third ovarian cycle: and we need to account for this fact. Note that a jump in observed trypanosome prevalence is not predicted by the fitted catalytic model (fig. 4).
Fitting data using a susceptible-infected (SI) model
The SI model defined by equations (3) and (4) requires knowledge of the survival probability. From the ovarian category distribution of the total number of adult female G. pallidipes, dissected to detect infections of T. vivax (fig. 3A), we estimated that φ = 0.976 and the daily mortality rate μ = -ln(φ) = 0.0242 (95% CI: 0.0238–0.0246). Using this estimate, and the procedures leading to equations 3 and 4, resulted in improved fits to the data for T. vivax, relative to the catalytic model (cf figs 4A and 5A; AIC values (table 2) are lower by about 20).
(A) Fitting age-specific percentage prevalences of A. T. vivax and B. T. congolense in female G. pallidipes, using a Susceptible-Infected (SI) model. Prevalences were calculated from the sums of the predicted numbers of infected and uninfected flies over successive 9-day pregnancy periods. It is assumed that negligible numbers of flies survive >15 ovulations. Error bars on observed prevalences are exact 95% confidence intervals for a binomial proportion.
For T. congolense, similarly, the fit was markedly improved (cf figs 4B and 5B), with the AIC values again lower by about 20, relative to the catalytic model. As is obvious from the inspection of fig. 5, however, the fit to the T. vivax data is markedly better than for T. congolense when the SI model is applied to the data for both species.
Fitting data using a susceptible-exposed-infected (SEI) model
The relatively poor fit of the SI model to T. congolense data results from the use of the value (τ) to allow for the time taken to mature an infection. Using the SEI model, described by equation 5, results in an improved fit to the data for T. congolense (cf figs 6A and 5B). The AIC for the SEI(1) model was 6.0 lower than that for the SI model (table 2, rows 6 and 7). For T. vivax the AIC criteria achieved using the SI and SEI(1) models differed by only 1.0 < 2.0 so that there was no ‘quality of fit’ difference between these two model forms (table 2, rows 2 and 3). Visually, the fit obtained using the SEI model (not shown) is almost indistinguishable from that for the SI model (fig. 5A).
(A) Fitting age-specific percentage prevalences of T. congolense in female G. pallidipes, using a Susceptible-Exposed-Infected (SEI(1)) model. (B) Predicted prevalence for ovarian age groups 1–15, calculated from the sums of the predicted numbers of infected and uninfected flies over successive 9-day pregnancy periods. It is assumed that negligible numbers of flies survive >15 ovulations. Error bars on observed prevalences are exact 95% confidence intervals for a binomial proportion.
In the SEI(2) model, our estimates of λ1 were not significantly different from zero, for both T. vivax and T. congolense (table 2). Thus, our data provided no support for a teneral effect, in which infection rates are higher for the first bloodmeal than for all subsequent bloodmeals. Our estimates of γ for both species, and of λ for T. congolense, were also not significant for the SEI(2) model (table 2). This suggests that our counts of infected and uninfected flies from only 8 ovarian categories cannot support reliable estimates of the 3 parameters of the SEI(2) model.
Note also that, in the single-lambda model, SEI(1), there was considerable overlap in CIs, between T. vivax and T. congolense, for the estimates of λ and γ (table 2, rows 3 and 7). Thus, while the point estimates of the two parameters indicate a higher probability of infection for T. vivax – and a shorter maturation time, resulting in higher predicted prevalences for this species (fig. 5B), the difference is not statistically significant.
Discussion
Recognition of the fact that significant proportions of female tsetse survive long enough in the field that they ovulate more than seven times means that it is not appropriate to use a catalytic model, which ignores fly mortality, to fit age-prevalence data for trypanosome infections. The more complex SI and SEI(1) models provide much improved fits to the data. Differences in the estimated infection rates between the models may appear small in absolute terms. However, the SI and SEI(1) values are 25 and 26% lower, respectively, than the values for the catalytic model in the case of T. vivax – and 29 and 16% lower in the case of T. congolense. This is as expected – because the catalytic model predicts the observed lifetime increase in prevalence over a lifespan of 70 days, in contrast to that same increase predicted over a lifespan of 142 days, for the SI and SEI(1) models.
The difficulty with the SI and SEI models is that they require knowledge of the mortality rate among adult female tsetse. This can be a problem because estimating this mortality from age distributions has proved to be extremely difficult (Hargrove and Ackley, Reference Hargrove and Ackley2015), particularly using samples collected over a short period of time. More believable estimates of mortality can be obtained, however, if samples are pooled over time periods of a year or more – such that problems associated with the seasonal variability of the population age structure are avoided (Van Sickle and Phelps, Reference Van Sickle and Phelps1988; Ackley and Hargrove, Reference Ackley and Hargrove2017). Accordingly, we have only applied our analysis to data pooled over the whole period of the study and have not attempted to estimate the parameters λ, λ1 and γ over shorter periods of time.
Coinfection rates
We assume in our analyses that the processes whereby female G. pallidipes acquire infections of the two trypanosome species, T. vivax and T. congolense, are completely independent of each other. In particular, being infected with one of the species should not alter the probability that a fly becomes infected with the other species. Coinfections, with more than one trypanosome species in a single fly, can be detected using techniques such as the polymerase chain reactions (PCR) that detect and identify, to species level, the DNA of different trypanosomes. The dissection process we used does not allow of such exact identification. Notice, however, that any fly that only has a mouthpart infection, can have only a T. vivax-type infection – because, if it also had a T. congolense infection, it would have trypanosomes in the midgut. The possibility of a vivax-congolense coinfection arises, therefore, only in flies where there are trypanosomes detected both in the mouthparts and the midgut. The following calculation suggests that the chances of this coinfection occurring for the female G. pallidipes in our study are negligible. Thus only 2.66 and 3.45% of all flies dissected were diagnosed to have a T. congolense-type or T. vivax-type infection, respectively. If a fly's chances of acquiring an infection of one of the trypanosome species is independent of whether it is already infected with the other species, then the probability of finding a fly that is positive for both species is thus 0.0266 × 0.0345 = 0.00092 or about 1 in 1000 flies. Coinfections would be even rarer if infection with one species of trypanosome reduced the probability of acquiring an infection with another trypanosome species. Coinfections could be more prevalent than indicated above only if an infection with one species of trypanosome increased the probability of acquiring an infection with the other species. We are not aware of any suggestion in the literature that this might happen.
T. congolense infections: the teneral effect
Welburn and Maudlin (Reference Welburn and Maudlin1992) found that laboratory G. m. morsitans were about 7-times more likely to acquire a T. congolense infection at their first meal than at any subsequent meal. By contrast, in G. pallidipes sampled in the field, we could not reject the null hypothesis of equality between the rate of infection for T. congolense at the first vs all later meals. Given the very large sample size available to us, and the excellent fit of the SEI(1) model to the data, we suggest that it might not even be feasible to separate λ1 and λ using field data.
The following calculations suggest that, in any case, a teneral effect of the order estimated by Welburn and Maudlin (Reference Welburn and Maudlin1992) would be of little importance even if it did occur – basically because of the longevity of tsetse, particularly females, and their need to feed every 2.5–3 days (Randolph et al., Reference Randolph, Williams, Rogers and Connor1992; Hargrove, Reference Hargrove1999a, Reference Hargrove1999b). Thus, suppose first that there were no teneral effect and that λ = 0.00271 at every meal, including the first (table 2, row 7). If, as we assume, a fly could live for 142 days (Hargrove et al., Reference Hargrove, Ouifki and Ameh2011), we expect it to take about 47 blood meals. Its chances of remaining uninfected with T. congolense, throughout its life, is then estimated as (1–0.00271)47 = 0.88, and the probability that it is infected is 1–0.88 = 12%, in reasonable accord with the observed data and the predictions of our model (fig. 6B). Suppose, alternatively, that there is a teneral effect of the order suggested by Welburn and Maudlin (Reference Welburn and Maudlin1992), such that λ1 = 7λ = 0.0190. Then the estimated probability of remaining uninfected for 142 days is (1–0.0190) × (1–0.00271)46 = 0.87 and the probability that the fly has become infected is 1–0.87 = 13%, only slightly higher than when there is no teneral effect. We thus find no support, at least in G. pallidipes, for Welburn and Maudlin's (Reference Welburn and Maudlin1992) suggestion that bloodmeals other than the first are of no importance epidemiologically and conclude, conversely, that the so-called teneral effect is of minor importance in the epidemiology of T. congolense as transmitted by female G. pallidipes.
Trypanosome prevalence in vertebrate hosts
The probability (λ) that tsetse will mature a trypanosome infection after taking a bloodmeal is the product λ = p1 × p2, where p1 is the probability that the bloodmeal they take is indeed infected with trypanosomes, and p2 is the conditional probability that a fly matures an infection given that it takes an infected meal. Results for the SEI(1) model (table 2, row 7) suggest λ ≈ 0.003 per 3-day period (i.e., per bloodmeal) for T. congolense infections of female G. pallidipes. From Welburn and Maudlin (Reference Welburn and Maudlin1992, their fig. 1) p1 for the first bloodmeal is about 0.7 and p2 for all later meals is about 1/7 of this value, i.e., p2 ≈ 0.1. If p2 takes a similarly low value for non-teneral female G. pallidipes at Rekomitjie, this suggests that the proportion p1 that take a meal from a host infected with T. congolense is p1 = λ/p2 = 0.003/0.1 = 0.03. That is to say, of the hosts fed on by G. pallidipes at Rekomitjie, only about 3% are harbouring sufficient T. congolense to ensure that tsetse feeding off them take an infected meal. If that is the case then the proportion of flies developing a mature infection after their first meal is λ = 0.03 × 0.7 = 0.021, or about 2.1% of the population. In fact, however, the 2% infection rate is only reached at the age of about 35 days (fig. 6A) – suggesting that p1 or p2, or both, are even lower than estimated.
For T. vivax our estimate for λ is also about 0.003 (table 2, line 3) but we do not have any estimate for p2. If it is markedly higher than for T. congolense then the proportion of animals with T. vivax infections is even lower than for T. congolense. For example, if p2 = 0.5 then p1 = λ2/p2 = 0.003/0.5 = 0.006 and we predict that only 0.6% of vertebrate hosts would be harbouring sufficient T. vivax to provide an infective meal. Such low levels of infection in wild hosts support our suggestion that the teneral effect is unlikely to be important in the epidemiology of either T. vivax of T. congolense; the chances of becoming infected by a single meal is simply too small to explain the observed age-related changes in trypanosome prevalence.
Sex and species differences
In the current study we have considered only females of G. pallidipes, because there is no method available for gauging the chronological age of males. Nonetheless, the degree of wing fray in males does provide estimates of relative age, albeit with less accuracy than is associated with the use of ovarian dissection in females. We also have available less extensive data for age-specific trypanosome prevalences in female G. m. morsitans sampled at Rekomitjie. Future work will investigate whether similar models, to those developed here for female G. pallidipes, are also appropriate for the pattern of acquisition of trypanosome infection for male G. pallidipes, and for G. m. morsitans. Given that the teneral effect was suggested based on trypanosome infections acquired by G. m. morsitans in the laboratory, we need to decide, particularly, whether the effect is important in field populations of this tsetse species.
Data
All data used, and materials generated, in this publication are available from the authors on reasonable request.
Acknowledgements
The authors acknowledge the support and continued interest of Dr William Shereni, Director of the Trypanosomiasis Control Department, Government of Zimbabwe. All experimental work was carried out at Rekomitjie Research Station in the Zambezi Valley. Our grateful thanks are due to the tsetse capture and dissection teams, ably led by Mr. Pio Chimanga. We acknowledge the contribution of Prof Mark Woolhouse, who trained the original field team to dissect flies to detect trypanosome infections.
Author contributions
JWH designed the study, collated the data and carried out preliminary analyses. JVS carried out detailed statistical analysis. Both authors contributed equally to the preparation of the manuscript and the preparation of all figures and tables.
Financial support
SACEMA receives core funding from the Department of Science and Innovation, Government of South Africa. The contents of this publication are the sole responsibility of the authors and do not reflect the views of any funding agency.
Conflict of interest
The authors declare that they have no competing interests.
Ethical standards
Not applicable; the study does not involve mammalian subjects.
Consent for publication
Both authors gave consent for this publication. No other consent required.
Open access
