The epidemiology of foot-and-mouth disease at the wildlife-livestock interface in northern Tanzania

Foot-and-mouth disease (FMD), a disease of cloven hooved animals caused by FMD virus (FMDV), is one of the most economically devastating diseases of livestock worldwide. The global burden of disease is borne largely by livestock-keepers in areas of Africa and Asia where the disease is endemic and where many people rely on livestock for their livelihoods and food-security. Yet, there are many gaps in our knowledge of the drivers of FMDV circulation in these settings. In East Africa, FMD epidemiology is complicated by the circulation of multiple FMDV serotypes (distinct antigenic variants) and by the presence of large populations of susceptible wildlife and domestic livestock. The African buffalo (Syncerus caffer) is the only wildlife species with consistent evidence of high levels of FMDV infection, and East Africa contains the largest population of this species globally. To inform FMD control in this region, key questions relate to heterogeneities in FMD prevalence and impacts in different livestock management systems and to the role of wildlife as a potential source of FMDV for livestock. To develop FMD control strategies and make best use of vaccine control options, serotype-specific patterns of circulation need to be characterised. In this study, the impacts and epidemiology of FMD were investigated across a range of traditional livestock-keeping systems in northern Tanzania, including pastoralist, agro-pastoralist and rural smallholder systems. Data were generated through field studies and laboratory analyses between 2010 and 2015. The study involved analysis of existing household survey data and generated serological data from cross-sectional livestock and buffalo samples and longitudinal cattle samples. Serological analyses included non-structural protein ELISAs, serotype-specific solid-phase competitive ELISAs, with optimisation to detect East African FMDV variants, and virus neutralisation testing. Risk factors for FMDV infection and outbreaks were investigated through analysis of cross-sectional serological data in conjunction with a case-control outbreak analysis. A novel Bayesian modeling approach was developed to infer serotype-specific infection history from serological data, and combined with virus isolation data from FMD outbreaks to characterise temporal and spatial patterns of serotype-specific infection. A high seroprevalence of FMD was detected in both northern Tanzanian livestock (69%, [66.5 - 71.4%] in cattle and 48.5%, [45.7-51.3%] in small ruminants) and in buffalo (80.9%, [74.7-86.1%]). Four different serotypes of FMDV (A, O, SAT1 and SAT2) were isolated from livestock. Up to three outbreaks per year were reported by households and active surveillance highlighted up to four serial outbreaks in the same herds within three years. Agro-pastoral and pastoral livestock keepers reported more frequent FMD outbreaks compared to smallholders. Households in all three management systems reported that FMD outbreaks caused significant impacts on milk production and sales, and on animals’ draught power, hence on crop production, with implications for food security and livelihoods. Risk factor analyses showed that older livestock were more likely to be seropositive for FMD (Odds Ratio [OR] 1.4 [1.4-1.5] per extra year) and that cattle (OR 3.3 [2.7-4.0]) were more likely than sheep and goats to be seropositive. Livestock managed by agro-pastoralists (OR 8.1 [2.8-23.6]) or pastoralists (OR 7.1 [2.9-17.6]) were more likely to be seropositive compared to those managed by smallholders. Larger herds (OR: 1.02 [1.01-1.03] per extra bovine) and those that recently acquired new livestock (OR: 5.57 [1.01 – 30.91]) had increased odds of suffering an FMD outbreak. Measures of potential contact with buffalo or with other FMD susceptible wildlife did not increase the likelihood of FMD in livestock in either the cross-sectional serological analysis or case-control outbreak analysis. The Bayesian model was validated to correctly infer from ELISA data the most recent serotype to infect cattle. Consistent with the lack of risk factors related to wildlife contact, temporal and spatial patterns of exposure to specific FMDV serotypes were not tightly linked in cattle and buffalo. In cattle, four serial waves of different FMDV serotypes that swept through southern Kenyan and northern Tanzanian livestock populations over a four-year period dominated infection patterns. In contrast, only two serotypes (SAT1 and SAT2) dominated in buffalo populations. Key conclusions are that FMD has a substantial impact in traditional livestock systems in East Africa. Wildlife does not currently appear to act as an important source of FMDV for East African livestock, and control efforts in the region should initially focus on livestock management and vaccination strategies. A novel modeling approach greatly facilitated the interpretation of serological data and may be a potent epidemiological tool in the African setting. There was a clear temporal pattern of FMDV antigenic dominance across northern Tanzania and southern Kenya. Longer-term research to investigate whether serotype-specific FMDV sweeps are truly predictable, and to shed light on FMD post-infection immunity in animals exposed to serial FMD infections is warranted.

Applied evolution: An integrated approach to studying life history traits in response to drug selection

The use of chemical control measures to reduce the impact of parasite and pest species has frequently resulted in the development of resistance. Thus, resistance management has become a key concern in human and veterinary medicine, and in agricultural production. Although it is known that factors such as gene flow between susceptible and resistant populations, drug type, application methods, and costs of resistance can affect the rate of resistance evolution, less is known about the impacts of density-dependent eco-evolutionary processes that could be altered by drug-induced mortality. The overall aim of this thesis was to take an experimental evolution approach to assess how life history traits respond to drug selection, using a free-living dioecious worm (Caenorhabditis remanei) as a model. In Chapter 2, I defined the relationship between C. remanei survival and Ivermectin dose over a range of concentrations, in order to control the intensity of selection used in the selection experiment described in Chapter 4. The dose-response data were also used to appraise curve-fitting methods, using Akaike Information Criterion (AIC) model selection to compare a series of nonlinear models. The type of model fitted to the dose response data had a significant effect on the estimates of LD50 and LD99, suggesting that failure to fit an appropriate model could give misleading estimates of resistance status. In addition, simulated data were used to establish that a potential cost of resistance could be predicted by comparing survival at the upper asymptote of dose-response curves for resistant and susceptible populations, even when differences were as low as 4%. This approach to dose-response modeling ensures that the maximum amount of useful information relating to resistance is gathered in one study. In Chapter 3, I asked how simulations could be used to inform important design choices used in selection experiments. Specifically, I focused on the effects of both within- and between-line variation on estimated power, when detecting small, medium and large effect sizes. Using mixed-effect models on simulated data, I demonstrated that commonly used designs with realistic levels of variation could be underpowered for substantial effect sizes. Thus, use of simulation-based power analysis provides an effective way to avoid under or overpowering a study designs incorporating variation due to random effects. In Chapter 4, I 3 investigated how Ivermectin dosage and changes in population density affect the rate of resistance evolution. I exposed replicate lines of C. remanei to two doses of Ivermectin (high and low) to assess relative survival of lines selected in drug-treated environments compared to untreated controls over 10 generations. Additionally, I maintained lines where mortality was imposed randomly to control for differences in density between drug treatments and to distinguish between the evolutionary consequences of drug treatment versus ecological processes affected by changes in density-dependent feedback. Intriguingly, both drug-selected and random-mortality lines showed an increase in survivorship when challenged with Ivermectin; the magnitude of this increase varied with the intensity of selection and life-history stage. The results suggest that interactions between density-dependent processes and life history may mediate evolved changes in susceptibility to control measures, which could result in misleading conclusions about the evolution of heritable resistance following drug treatment. In Chapter 5, I investigated whether the apparent changes in drug susceptibility found in Chapter 4 were related to evolved changes in life-history of C. remanei populations after selection in drug-treated and random-mortality environments. Rapid passage of lines in the drug-free environment had no effect on the measured life-history traits. In the drug-free environment, adult size and fecundity of drug-selected lines increased compared to the controls but drug selection did not affect lifespan. In the treated environment, drug-selected lines showed increased lifespan and fecundity relative to controls. Adult size of randomly culled lines responded in a similar way to drug-selected lines in the drug-free environment, but no change in fecundity or lifespan was observed in either environment. The results suggest that life histories of nematodes can respond to selection as a result of the application of control measures. Failure to take these responses into account when applying control measures could result in adverse outcomes, such as larger and more fecund parasites, as well as over-estimation of the development of genetically controlled resistance. In conclusion, my thesis shows that there may be a complex relationship between drug selection, density-dependent regulatory processes and life history of populations challenged with control measures. This relationship could have implications for how resistance is monitored and managed if life histories of parasitic species show such eco-evolutionary responses to drug application.