Mechanistic movement models reveal ecological drivers of tick-borne pathogen spread

Identifying ecological drivers of tick-borne pathogen spread has great value for tick-borne disease management. However, theoretical investigations into the consequences of host movement behaviour on pathogen spread dynamics in heterogeneous landscapes remain limited because spatially explicit epidemiological models that incorporate more realistic mechanisms governing host movement are rare. We built a mechanistic movement model to investigate how the interplay between multiple ecological drivers affects the risk of tick-borne pathogen spread across heterogeneous landscapes. We used the model to generate simulations of tick dispersal by migratory birds and terrestrial hosts across theoretical landscapes varying in resource aggregation, and we performed a sensitivity analysis to explore the impacts of different parameters on the infected tick spread rate, tick infection prevalence and infected tick density. Our findings highlight the importance of host movement and tick population dynamics in explaining the infected tick spread rate into new regions. Tick infection prevalence and infected tick density were driven by predictors related to the infection process and tick population dynamics, respectively. Our results suggest that control strategies aiming to reduce tick burden on tick reproduction hosts and encounter rate between immature ticks and pathogen amplification hosts will be most effective at reducing tick-borne disease risk.     

In-roads to the spread of antibiotic resistance: regional patterns of microbial transmission in northern coastal Ecuador

The evolution of antibiotic resistance (AR) increases treatment cost and probability of failure, threatening human health worldwide. The relative importance of individual antibiotic use, environmental transmission and rates of introduction of resistant bacteria in explaining community AR patterns is poorly understood. Evaluating their relative importance requires studying a region where they vary. The construction of a new road in a previously roadless area of northern coastal Ecuador provides a valuable natural experiment to study how changes in the social and natural environment affect the epidemiology of resistant Escherichia coli. We conducted seven bi-annual 15 day surveys of AR between 2003 and 2008 in 21 villages. Resistance to both ampicillin and sulphamethoxazole was the most frequently observed profile, based on antibiogram tests of seven antibiotics from 2210 samples. The prevalence of enteric bacteria with this resistance pair in the less remote communities was 80 per cent higher than in more remote communities (OR = 1.8 [1.3, 2.3]). This pattern could not be explained with data on individual antibiotic use. We used a transmission model to help explain this observed discrepancy. The model analysis suggests that both transmission and the rate of introduction of resistant bacteria into communities may contribute to the observed regional scale AR patterns, and that village-level antibiotic use rate determines which of these two factors predominate. While usually conceived as a main effect on individual risk, antibiotic use rate is revealed in this analysis as an effect modifier with regard to community-level risk of resistance.     

Stochastic extinction and the selection of the transmission mode in microparasites.

Stochastic fluctuations in the transmission process of microparasites generate a risk of parasite extinction that cannot be assessed by deterministic models, especially in host populations of small size. While this risk of extinction represents a strong selection pressure for microparasites, it is usually not clearly separated from the deterministic ones. We suggest here that this stochastic selection pressure can affect the selection of the transmission mode of microparasites. To avoid extinction, parasites should maximize their inter-population transmission to ensure frequent reintroductions. Since the types of contacts may differ if congeners belong to the same or distinct populations, strains that are mainly transmitted through inter-population contacts might be selected. To examine this assumption, we analyse the issue of the competition between two strains differing in their transmission mode using a stochastic metapopulation model in which hosts may display different behaviours inside and outside their populations. We show that stochastic selection pressures may drive parasite evolution towards a transmission mode that maximizes the persistence of the parasite. We study the conditions under which stochastic selection pressures may surpass the deterministic ones. Our results are illustrated by the cases of feline immunodeficiency virus in cats and of sexually transmitted diseases in mammals.     

Incorporating demographic stochasticity into multi-strain epidemic models: application to influenza A

We develop mathematical models of the transmission and evolution of multi-strain pathogens that incorporate strain extinction and the stochastic generation of new strains via mutation. The dynamics resulting from these models is then examined with the applied aim of understanding the mechanisms underpinning the evolution and dynamics of rapidly mutating pathogens, such as human influenza viruses. Our approach, while analytically relatively simple, gives results that are qualitatively similar to those obtained from much more complex individually based simulation models. We examine strain dynamics as a function of cross-immunity and key transmission parameters, and show that introducing strain extinction and modelling mutation as a stochastic process significantly changes the model dynamics, leading to lower strain diversity, reduced infection prevalence and shorter strain lifetimes. Finally, we incorporate transient strain-transcending immunity in the model and demonstrate that it reduces strain diversity further, giving patterns of sequential strain replacement similar to that seen in human influenza A viruses.     

Campylobacter jejuni colonization and transmission in broiler chickens: a modelling perspective.

Campylobacter jejuni is one of the most common causes of acute enteritis in the developed world. The consumption of contaminated poultry, where C. jejuni is believed to be a commensal organism, is a major risk factor. However, the dynamics of this colonization process in commercially reared chickens is still poorly understood. Quantification of these dynamics of infection at an individual level is vital to understand transmission within populations and formulate new control strategies. There are multiple potential routes of introduction of C. jejuni into a commercial flock. Introduction is followed by a rapid increase in environmental levels of C. jejuni and the level of colonization of individual broilers. Recent experimental and epidemiological evidence suggest that the celerity of this process could be masking a complex pattern of colonization and extinction of bacterial strains within individual hosts. Despite the rapidity of colonization, experimental transmission studies exhibit a highly variable and unexplained delay time in the initial stages of the process. We review past models of transmission of C. jejuni in broilers and consider simple modifications, motivated by the plausible biological mechanisms of clearance and latency, which could account for this delay. We show how simple mathematical models can be used to guide the focus of experimental studies by providing testable predictions based on our hypotheses. We conclude by suggesting that competition experiments could be used to further understand the dynamics and mechanisms underlying the colonization process. The population models for such competition processes have been extensively studied in other ecological and evolutionary contexts. However, C. jejuni can potentially adapt phenotypically through phase variation in gene expression, leading to unification of ecological and evolutionary time-scales. For a theoretician, the colonization dynamics of C. jejuni offer an experimental system to explore these 'phylodynamics', the synthesis of population dynamics and evolutionary biology.     

Time and dose-dependent risk of pneumococcal pneumonia following influenza: a model for within-host interaction between influenza and Streptococcus pneumoniae

A significant fraction of seasonal and in particular pandemic influenza deaths are attributed to secondary bacterial infections. In animal models, influenza virus predisposes hosts to severe infection with both Streptococcus pneumoniae and Staphylococcus aureus. Despite its importance, the mechanistic nature of the interaction between influenza and pneumococci, its dependence on the timing and sequence of infections as well as the clinical and epidemiological consequences remain unclear. We explore an immune-mediated model of the viral–bacterial interaction that quantifies the timing and the intensity of the interaction. Taking advantage of the wealth of knowledge gained from animal models, and the quantitative understanding of the kinetics of pathogen-specific immunological dynamics, we formulate a mathematical model for immune-mediated interaction between influenza virus and S. pneumoniae in the lungs. We use the model to examine the pathogenic effect of inoculum size and timing of pneumococcal invasion relative to influenza infection, as well as the efficacy of antivirals in preventing severe pneumococcal disease. We find that our model is able to capture the key features of the interaction observed in animal experiments. The model predicts that introduction of pneumococcal bacteria during a 4–6 day window following influenza infection results in invasive pneumonia at significantly lower inoculum size than in hosts not infected with influenza. Furthermore, we find that antiviral treatment administered later than 4 days after influenza infection was not able to prevent invasive pneumococcal disease. This work provides a quantitative framework to study interactions between influenza and pneumococci and has the potential to accurately quantify the interactions. Such quantitative understanding can form a basis for effective clinical care, public health policies and pandemic preparedness.     

The control of sea lice in Atlantic salmon by selective breeding

Sea lice threaten the welfare of farmed Atlantic salmon and the sustainability of fish farming across the world. Chemical treatments are the major method of control but drug resistance means that alternatives are urgently needed. Selective breeding can be a cheap and effective alternative. Here, we combine experimental trials and diagnostics to provide a practical protocol for quantifying resistance to sea lice. We then combined quantitative genetics with epidemiological modelling to make the first prediction of the response to selection, quantified in terms of reduced need for chemical treatments. We infected over 1400 young fish with Lepeophtheirus salmonis, the most important species in the Northern Hemisphere. Mechanisms of resistance were expressed early in infection. Consequently, the number of lice per fish and the ranking of families were very similar at 7 and 17 days post infection, providing a stable window for assessing susceptibility to infection. The heritability of lice numbers within this time window was moderately high at 0.3, confirming that selective breeding is viable. We combined an epidemiological model of sea lice infection and control on a salmon farm with genetic variation in susceptibility among individuals. We simulated 10 generations of selective breeding and examined the frequency of treatments needed to control infection. Our model predicted that substantially fewer chemical treatments are needed to control lice outbreaks in selected populations and chemical treatment could be unnecessary after 10 generations of selection. Selective breeding for sea lice resistance should reduce the impact of sea lice on fish health and thus substantially improve the sustainability of Atlantic salmon production.     

Predicting the effect of climate change on African trypanosomiasis: integrating epidemiology with parasite and vector biology

Climate warming over the next century is expected to have a large impact on the interactions between pathogens and their animal and human hosts. Vector-borne diseases are particularly sensitive to warming because temperature changes can alter vector development rates, shift their geographical distribution and alter transmission dynamics. For this reason, African trypanosomiasis (sleeping sickness), a vector-borne disease of humans and animals, was recently identified as one of the 12 infectious diseases likely to spread owing to climate change. We combine a variety of direct effects of temperature on vector ecology, vector biology and vector–parasite interactions via a disease transmission model and extrapolate the potential compounding effects of projected warming on the epidemiology of African trypanosomiasis. The model predicts that epidemics can occur when mean temperatures are between 20.7°C and 26.1°C. Our model does not predict a large-range expansion, but rather a large shift of up to 60 per cent in the geographical extent of the range. The model also predicts that 46–77 million additional people may be at risk of exposure by 2090. Future research could expand our analysis to include other environmental factors that influence tsetse populations and disease transmission such as humidity, as well as changes to human, livestock and wildlife distributions. The modelling approach presented here provides a framework for using the climate-sensitive aspects of vector and pathogen biology to predict changes in disease prevalence and risk owing to climate change.     

Assessing the impact of adherence to anti-retroviral therapy on treatment failure and resistance evolution in HIV

The adherence of patients to therapy is a crucial factor for successful HIV anti-retroviral therapy. Imperfect adherence may lead to treatment failure, which can cause the emergence of resistance within viral populations. We have developed a stochastic model that incorporates compartments of latently infected cells and virus genotypes with different susceptibilities to three simultaneously used drugs. With this model, we study the impact of several key parameters on the probability of treatment failure, i.e. insufficient viral suppression, and the emergence of resistance. Specifically, we consider the impact of drug dosage, drug half-lives, fitness costs for resistance, different basic reproductive numbers of the virus and the influence of pre-existing mutations under various levels of adherence. Furthermore, we also investigate the influence of different temporal distributions of non-adherent days (drug holidays) during a treatment. Factors that promote resistance evolution include a high reproductive number, extended drug holidays and poor adherence. Pre-existing mutations only have a substantial effect if they confer resistance against more than one drug. Overall, our study highlights the importance of the interactions between imperfect adherence, pharmacodynamics, pharmacokinetics and latently infected cells for our understanding of drug resistance and therapy failure in HIV anti-retroviral therapy.     

Complexity and anisotropy in host morphology make populations less susceptible to epidemic outbreaks

One of the challenges in epidemiology is to account for the complex morphological structure of hosts such as plant roots, crop fields, farms, cells, animal habitats and social networks, when the transmission of infection occurs between contiguous hosts. Morphological complexity brings an inherent heterogeneity in populations and affects the dynamics of pathogen spread in such systems. We have analysed the influence of realistically complex host morphology on the threshold for invasion and epidemic outbreak in an SIR (susceptible–infected–recovered) epidemiological model. We show that disorder expressed in the host morphology and anisotropy reduces the probability of epidemic outbreak and thus makes the system more resistant to epidemic outbreaks. We obtain general analytical estimates for minimally safe bounds for an invasion threshold and then illustrate their validity by considering an example of host data for branching hosts (salamander retinal ganglion cells). Several spatial arrangements of hosts with different degrees of heterogeneity have been considered in order to separately analyse the role of shape complexity and anisotropy in the host population. The estimates for invasion threshold are linked to morphological characteristics of the hosts that can be used for determining the threshold for invasion in practical applications.     

Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity

Dengue, a mosquito-borne virus of humans, infects over 50 million people annually. Infection with any of the four dengue serotypes induces protective immunity to that serotype, but does not confer long-term protection against infection by other serotypes. The immunological interactions between serotypes are of central importance in understanding epidemiological dynamics and anticipating the impact of dengue vaccines. We analysed a 38-year time series with 12 197 serotyped dengue infections from a hospital in Bangkok, Thailand. Using novel mechanistic models to represent different hypothesized immune interactions between serotypes, we found strong evidence that infection with dengue provides substantial short-term cross-protection against other serotypes (approx. 1–3 years). This is the first quantitative evidence that short-term cross-protection exists since human experimental infection studies performed in the 1950s. These findings will impact strategies for designing dengue vaccine studies, future multi-strain modelling efforts, and our understanding of evolutionary pressures in multi-strain disease systems.     

A metapopulation modelling framework for gonorrhoea and other sexually transmitted infections in heterosexual populations

Gonorrhoea continues to be a public health problem in the UK, and is the second most common bacterial sexually transmitted infection (STI) after chlamydia. In the UK, gonorrhoea is disproportionately concentrated in epidemiologically distinct subpopulations, with much higher incidence rates in young people, some ethnic minorities and inner city subpopulations. The original model of STI transmission proposed by Hethcote and Yorke explained some of these features through the concept of the ‘core group’. Since then, several authors have modified the original model approach to include multiple sexual activity classes, but found this modelling approach to be inadequate when applied to low-prevalence settings such as the UK. We present a metapopulation framework for modelling gonorrhoea and other STIs. The model proposes that the epidemiology of gonorrhoea is largely driven by subpopulations with higher than average concentrations of individuals with high sexual risk activity. We show how this conceptualization of gonococcal epidemiology overcomes key limitations associated with some of the prior efforts to model gonorrhoea. We also use the model to explain several epidemiological features of gonorrhoea, such as its asymmetric distribution across subpopulations, and the contextual risk experienced by members of at-risk subpopulations. Finally, we extend the model to explain the distribution of other STIs, using chlamydia as an example of a more ubiquitous bacterial STI.     

Role of high-dose exposure in transmission hot zones as a driver of SARS-CoV-2 dynamics

Epidemiological data about SARS-CoV-2 spread indicate that the virus is not transmitted uniformly in the population. The transmission tends to be more effective in select settings that involve exposure to relatively high viral dose, such as in crowded indoor settings, assisted living facilities, prisons or food processing plants. To explore the effect on infection dynamics, we describe a new mathematical model where transmission can occur (i) in the community at large, characterized by low-dose exposure and mostly mild disease, and (ii) in so-called transmission hot zones, characterized by high-dose exposure that can be associated with more severe disease. The model yields different types of epidemiological dynamics, depending on the relative importance of hot zone and community transmission. Interesting dynamics occur if the rate of virus release/deposition from severely infected people is larger than that of mildly infected individuals. Under this assumption, we find that successful infection spread can hinge upon high-dose hot zone transmission, yet the majority of infections are predicted to occur in the community at large with mild disease. In this regime, residual hot zone transmission can account for continued virus spread during community lockdowns, and the suppression of hot zones after community interventions are relaxed can cause a prolonged lack of infection resurgence following the reopening of society. This gives rise to the notion that targeted interventions specifically reducing virus transmission in the hot zones have the potential to suppress overall infection spread, including in the community at large. Epidemiological trends in the USA and Europe are interpreted in light of this model.     

Escherichia coli O157 infection on Scottish cattle farms: dynamics and control

In this study, we parametrize a stochastic individual-based model of the transmission dynamics of Escherichia coli O157 infection among Scottish cattle farms and use the model to predict the impacts of both targeted and non-targeted interventions. We first generate distributions of model parameter estimates using Markov chain Monte Carlo methods. Despite considerable uncertainty in parameter values, each set of parameter values within the 95th percentile range implies a fairly similar impact of interventions. Interventions that reduce the transmission coefficient and/or increase the recovery rate of infected farms (e.g. via vaccination and biosecurity) are much more effective in reducing the level of infection than reducing cattle movement rates, which improves effectiveness only when the overall control effort is small. Targeted interventions based on farm-level risk factors are more efficient than non-targeted interventions. Herd size is a major determinant of risk of infection, and our simulations confirmed that targeting interventions at farms with the largest herds is almost as effective as targeting based on overall risk. However, because of the striking characteristic that the infection force depends weakly on the number of infected farms, no interventions that are less than 100 per cent effective can eradicate E. coli O157 infection from Scottish cattle farms, implying that eliminating the disease is impractical.     

Modelling the effects of past and future climate on the risk of bluetongue emergence in Europe

Vector-borne diseases are among those most sensitive to climate because the ecology of vectors and the development rate of pathogens within them are highly dependent on environmental conditions. Bluetongue (BT), a recently emerged arboviral disease of ruminants in Europe, is often cited as an illustration of climate''s impact on disease emergence, although no study has yet tested this association. Here, we develop a framework to quantitatively evaluate the effects of climate on BT''s emergence in Europe by integrating high-resolution climate observations and model simulations within a mechanistic model of BT transmission risk. We demonstrate that a climate-driven model explains, in both space and time, many aspects of BT''s recent emergence and spread, including the 2006 BT outbreak in northwest Europe which occurred in the year of highest projected risk since at least 1960. Furthermore, the model provides mechanistic insight into BT''s emergence, suggesting that the drivers of emergence across Europe differ between the South and the North. Driven by simulated future climate from an ensemble of 11 regional climate models, the model projects increase in the future risk of BT emergence across most of Europe with uncertainty in rate but not in trend. The framework described here is adaptable and applicable to other diseases, where the link between climate and disease transmission risk can be quantified, permitting the evaluation of scale and uncertainty in climate change''s impact on the future of such diseases.     

Non-random biodiversity loss underlies predictable increases in viral disease prevalence

Disease dilution (reduced disease prevalence with increasing biodiversity) has been described for many different pathogens. Although the mechanisms causing this phenomenon remain unclear, the disassembly of communities to predictable subsets of species, which can be caused by changing climate, land use or invasive species, underlies one important hypothesis. In this case, infection prevalence could reflect the competence of the remaining hosts. To test this hypothesis, we measured local host species abundance and prevalence of four generalist aphid-vectored pathogens (barley and cereal yellow dwarf viruses) in a ubiquitous annual grass host at 10 sites spanning 2000 km along the North American West Coast. In laboratory and field trials, we measured viral infection as well as aphid fecundity and feeding preference on several host species. Virus prevalence increased as local host richness declined. Community disassembly was non-random: ubiquitous hosts dominating species-poor assemblages were among the most competent for vector production and virus transmission. This suggests that non-random biodiversity loss led to increased virus prevalence. Because diversity loss is occurring globally in response to anthropogenic changes, such work can inform medical, agricultural and veterinary disease research by providing insights into the dynamics of pathogens nested within a complex web of environmental forces.     

Analytical methods for quantifying environmental connectivity for the control and surveillance of infectious disease spread

The sustained transmission and spread of environmentally mediated infectious diseases is governed in part by the dispersal of parasites, disease vectors and intermediate hosts between sites of transmission. Functional geospatial models can be used to quantify and predict the degree to which environmental features facilitate or limit connectivity between target populations, yet typical models are limited in their geographical and analytical approach, providing simplistic, global measures of connectivity and lacking methods to assess the epidemiological implications of fine-scale heterogeneous landscapes. Here, functional spatial models are applied to problems of surveillance and control of the parasitic blood fluke Schistosoma japonicum and its intermediate snail host Oncomelania haupensis in western China. We advance functional connectivity methods by providing an analytical framework to (i) identify nodes of transmission where the degree of connectedness to other villages, and thus the potential for disease spread, is higher than is estimated using Euclidean distance alone and (ii) (re)organize transmission sites into disease surveillance units based on second-order relationships among nodes using non-Euclidean distance measures, termed effective geographical distance (EGD). Functional environmental models are parametrized using ecological information on the target organisms, and pair-wise distributions of inter-node EGD are estimated. A Monte Carlo rank product analysis is presented to identify nearby nodes under alternative distance models. Nodes are then iteratively embedded into EGD space and clustered using a k-means algorithm to group villages into ecologically meaningful surveillance groups. A consensus clustering approach is taken to derive the most stable cluster structure. The results indicate that novel relationships between nodes are revealed when non-Euclidean, ecologically determined distance measures are used to quantify connectivity in heterogeneous landscapes. These connections are not evident when analysing nodes in Euclidean space, and thus surveillance and control activities planned using Euclidean distance measures may be suboptimal. The methods developed here provide a quantitative framework for assessing the effectiveness of ecologically grounded surveillance systems and of control and prevention strategies for environmentally mediated diseases.     

The dynamic influence of genetic variation on the susceptibility of sheep to gastrointestinal nematode infection.

The interaction between sheep and the nematode Teladorsagia circumcincta is one of the best understood of all host-parasite interactions. Following infection, there is considerable variation among lambs in the number of nematode eggs produced, the number of early fourth-stage larvae and the number of adult worms in the mucosa. These traits have a high variance to mean ratio (i.e. they are overdispersed or aggregated among hosts), they are skewed and approximately negative binomially distributed. The sources of overdispersion are differences among lambs in the ingestion of infective larvae and the immune response. Both forces can produce aggregation but their relative importance is unknown. The key components of variation can be identified by variance analysis. The sum of the average effects of polymorphic genes is known as additive genetic variation and this increases essentially from zero at one month of age to quite high values at six months of age. The major mechanism underlying genetic variation appears to be the differences among individuals in immune responses. Two of the major sources of variation in immune responses are differences in antigen recognition and differences in the type of cytokines produced. Genes that influence both these sources of variation are associated with differences in resistance to nematode infection. Therefore, much of the heterogeneity among animals in parasite transmission appears to be due to genetic variation in immune responsiveness.     

The effect of population structure on the emergence of drug resistance during influenza pandemics.

The spread of H5N1 avian influenza and the recent high numbers of confirmed human cases have raised international concern about the possibility of a new pandemic. Therefore, antiviral drugs are now being stockpiled to be used as a first line of defence. The large-scale use of antivirals will however exert a strong selection pressure on the virus, and may lead to the emergence of drug-resistant strains. A few mathematical models have been developed to assess the emergence of drug resistance during influenza pandemics. These models, however, neglected the spatial structure of large populations and the stochasticity of epidemic and demographic processes. To assess the impact of population structure and stochasticity, we modify and extend a previous model of influenza epidemics into a metapopulation model which takes into account the division of large populations into smaller units, and develop deterministic and stochastic versions of the model. We find that the dynamics in a fragmented population is less explosive, and, as a result, prophylaxis will prevent more infections and lead to fewer resistant cases in both the deterministic and stochastic model. While in the deterministic model the final level of resistance during treatment is not affected by fragmentation, in the stochastic model it is. Our results enable us to qualitatively extrapolate the prediction of deterministic, homogeneous-mixing models to more realistic scenarios.