In this article we present a spatially-explicit stochastic metapopulation model to simulate the spread and control of footand-mouth disease (FMD) in an endemic setting. We parameterise and validate the model using detailed outbreak data from the Republic of Turkey, 2001-2012. Subsequently, we assess the efficacy of ring vaccination, mass vaccination, and livestock movement restrictions with regards to incidence-reduction and likelihood of eradication. Key findings show that ring vaccination and mass vaccination are most effective for controlling endemic FMD; livestock movement controls do not lead to eradication on average. Combined biannual mass vaccination and 10 km ring vaccination around detected farms emerged as the optimal approach to maximise the probability of FMD elimination within a few years. The probability of disease detection, radius of ring vaccination, and coverage of mass vaccination are most critical to optimize for this policy. Our results suggest countries wishing to control the disease within their borders should focus on comprehensive surveillance and vaccination campaigns as their main policy goals. In summary, vaccination-based policies are more effective than movement restrictions in the endemic context.
African horse sickness is an equine orbivirus transmitted by Culicoides Latreille biting midges. In the last 80 years, it has caused several devastating outbreaks in the equine population in Europe, the Far and Middle East, North Africa, South-East Asia, and sub-Saharan Africa. The disease is endemic in South Africa; however, a unique control area has been set up in the Western Cape where increased surveillance and control measures have been put in place. A deterministic metapopulation model was developed to explore if an outbreak might occur, and how it might develop, if a latently infected horse was to be imported into the control area, by varying the geographical location and months of import. To do this, a previously published ordinary differential equation model was developed with a metapopulation approach and included a vaccinated horse population. Outbreak length, time to peak infection, number of infected horses at the peak, number of horses overall affected (recovered or dead), re-emergence, and Rv (the basic reproduction number in the presence of vaccination) were recorded and displayed using GIS mapping. The model predictions were compared to previous outbreak data to ensure validity. The warmer months (November to March) had longer outbreaks than the colder months (May to September), took more time to reach the peak, and had a greater total outbreak size with more horses infected at the peak. Rv appeared to be a poor predictor of outbreak dynamics for this simulation. A sensitivity analysis indicated that control measures such as vaccination and vector control are potentially effective to manage the spread of an outbreak, and shortening the vaccination window to July to September may reduce the risk of vaccine-associated outbreaks.
Background Culicoides biting midges exhibit a global spatial distribution and are the main vectors of several viruses of veterinary importance, including bluetongue (BT) and African horse sickness (AHS). Many environmental and anthropological factors contribute to their ability to live in a variety of habitats, which have the potential to change over the years as the climate changes. Therefore, as new habitats emerge, the risk for new introductions of these diseases of interest to occur increases. The aim of this study was to model distributions for two primary vectors for BT and AHS (Culicoides imicola and Culicoides bolitinos) using random forest (RF) machine learning and explore the relative importance of environmental and anthropological factors in a region of South Africa with frequent AHS and BT outbreaks. Methods Culicoides capture data were collected between 1996 and 2022 across 171 different capture locations in the Western Cape. Predictor variables included climate-related variables (temperature, precipitation, humidity), environment-related variables (normalised difference vegetation index-NDVI, soil moisture) and farm-related variables (livestock densities). Random forest (RF) models were developed to explore the spatial distributions of C. imicola, C. bolitinos and a merged species map, where both competent vectors were combined. The maps were then compared to interpolation maps using the same capture data as well as historical locations of BT and AHS outbreaks. Results Overall, the RF models performed well with 75.02%, 61.6% and 74.01% variance explained for C. imicola, C. bolitinos and merged species models respectively. Cattle density was the most important predictor for C. imicola and water vapour pressure the most important for C. bolitinos. Compared to interpolation maps, the RF models had higher predictive power throughout most of the year when species were modelled individually; however, when merged, the interpolation maps performed better in all seasons except winter. Finally, midge densities did not show any conclusive correlation with BT or AHS outbreaks. Conclusion This study yielded novel insight into the spatial abundance and drivers of abundance of competent vectors of BT and AHS. It also provided valuable data to inform mathematical models exploring disease outbreaks so that Culicoides-transmitted diseases in South Africa can be further analysed.
Bluetongue virus (BT) is a vector-borne virus that causes a disease, called bluetongue, which results in significant economic loss and morbidity in sheep, cattle, goats and wild ungulates across all continents of the world except Antarctica. Despite the geographical breadth of its impact, most BT epidemiological models are informed by parameters derived from the 2006-2009 BTV-8 European outbreak. The aim of this study was to develop a highly adaptable model for BT which could be used elsewhere in the world, as well as to identify the parameters which most influence outbreak dynamics, so that policy makers can be properly informed with the most current information to aid in disease planning.To provide a framework for future outbreak modelling and an updated parameterisation that reflects natural variation in infections, a newly developed and parameterised two-host, two-vector species ordinary differential equation model was formulated and analysed. The model was designed to be adaptable to be implemented in any region of the world and able to model both epidemic and endemic scenarios. It was parameterised using a systematic literature review of host-to-vector and vector-to-host transmission rates, host latent periods, host infectious periods, and vaccine protection factors. The model was demonstrated using the updated parameters, with South Africa as a setting based on the Western Cape’s known cattle and sheep populations, local environmental parameters, and Culicoides spp. presence data.The sensitivity analysis identified that the duration of the infectious period for sheep and cows had the greatest impact on the outbreak length and number of animals infected at the peak of the outbreak. Transmission rates from cows and sheep to C. imicola midges greatly influenced the day on which the peak of the outbreak occurred, along with the duration of incubation period, and infectious period for cows. Finally, the protection factor of the vaccine had the greatest influence on the total number of animals infected. This knowledge could aid in the development of control measures.Due to gradual climate and anthropological change resulting in alterations in vector habitat suitability, BT outbreaks are likely to continue to increase in range and frequency. Therefore, this research provides an updated BT modelling framework for future outbreaks around the world to explore transmission, outbreak dynamics and control measures.
Infectious disease dynamics operate across biological scales: pathogens replicate within hosts but transmit among populations. Functional changes in the pathogen-host interaction thus generate cascading effects across organizational scales. We investigated within-host dynamics and among-host transmission of three strains (SAT-1, -2, -3) of foot- and-mouth disease viruses (FMDVs) in their wildlife host, African buffalo. We combined data on viral dynamics and host immune responses with mathematical models to ask the following questions: How do viral and immune dynamics vary among strains? Which viral and immune parameters determine viral fitness within hosts? And how do within-host dynamics relate to virus transmission? Our data reveal contrasting withinhost dynamics among viral strains, with SAT-2 eliciting more rapid and effective immune responses than SAT-1 and SAT-3. Within-host viral fitness was overwhelmingly determined by variation among hosts in immune response activation rates but not by variation among individual hosts in viral growth rate. Our analyses investigating across-scale linkages indicate that viral replication rate in the host correlates with transmission rates among buffalo and that adaptive immune activation rate determines the infectious period. These parameters define the virus's relative basic reproductive number (R0), suggesting that viral invasion potential may be predictable from within-host dynamics.
We collected >40,000 mosquitoes from 5 provinces in South Africa during 2011–2018 and screened for zoonotic flaviviruses. We detected West Nile virus in mosquitoes from conservation and periurban sites and potential new mosquito vectors; Banzi virus was rare. Our results suggest flavivirus transmission risks are increasing in South Africa.
Transboundary livestock diseases are a high priority for policy makers because of the serious economic burdens associated with infection. In order to make well informed preparedness and response plans, policy makers often utilize mathematical models to understand possible outcomes of different control strategies and outbreak scenarios. Many of these models focus on the transmission between herds and the overall trajectory of the outbreak. While the course of infection within herds has not been the focus of the majority of models, a thorough understanding of within-herd dynamics can provide valuable insight into a disease system by providing information on herd-level biological properties of the infection, which can be used to inform decision making in both endemic and outbreak settings and to inform larger between-herd models. In this study, we develop three stochastic simulation models to study within-herd foot and mouth disease dynamics and the implications of different empirical data-based assumptions about the timing of the onset of infectiousness and clinical signs. We also study the influence of herd size and the proportion of the herd that is initially infected on the outcome of the infection. We find that increasing herd size increases the duration of infectiousness and that the size of the herd plays a more significant role in determining this duration than the number of initially infected cattle in that herd. We also find that the assumptions made regarding the onset of infectiousness and clinical signs, which are based on contradictory empirical findings, can result in the predictions about when infection would be detectable differing by several days. Therefore, the disease progression used to characterize the course of infection in a single bovine host could have significant implications for determining when herds can be detected and subsequently controlled; the timing of which could influence the overall predicted trajectory of outbreaks.
The prevalence and distribution of African alphaviruses such as chikungunya have increased in recent years. Therefore, a better understanding of the local distribution of alphaviruses in vectors across the African continent is important. Here, entomological surveillance was performed from 2014 to 2018 at selected sites in north-eastern parts of South Africa where alphaviruses have been identified during outbreaks in humans and animals in the past. Mosquitoes were collected using a net, CDC-light, and BG-traps. An alphavirus genus-specific nested RT-PCR was used for screening, and positive pools were confirmed by sequencing and phylogenetic analysis. We collected 64,603 mosquitoes from 11 genera, of which 39,035 females were tested. Overall, 1462 mosquito pools were tested, of which 21 were positive for alphaviruses. Sindbis (61.9%, N = 13) and Middelburg (28.6%, N = 6) viruses were the most prevalent. Ndumu virus was detected in two pools (9.5%, N = 2). No chikungunya positive pools were identified. Arboviral activity was concentrated in peri-urban, rural, and conservation areas. A range of Culicidae species, including Culex univittatus, Cx. pipiens s.l., Aedes durbanensis, and the Ae. dentatus group, were identified as potential vectors. These findings confirm the active circulation and distribution of alphaviruses in regions where human or animal infections were identified in South Africa.
Many infectious pathogens are shared through social interactions, and examining host connectivity has offered valuable insights for understanding patterns of pathogen transmission across wildlife species. African buffalo are social ungulates and important reservoirs of directly-transmitted pathogens that impact numerous wildlife and livestock species. Here, we analyzed African buffalo social networks to quantify variation in close contacts, examined drivers of contact heterogeneity, and investigated how the observed contact patterns affect pathogen invasion likelihoods for a wild social ungulate. We collected continuous association data using proximity collars and sampled host traits approximately every 2 months during a 15-month study period in Kruger National Park, South Africa. Although the observed herd was well connected, with most individuals contacting each other during each bimonthly interval, our analyses revealed striking heterogeneity in close-contact associations among herd members. Network analysis showed that individual connectivity was stable over time and that individual age, sex, reproductive status, and pairwise genetic relatedness were important predictors of buffalo connectivity. Calves were the most connected members of the herd, and adult males were the least connected. These findings highlight the role susceptible calves may play in the transmission of pathogens within the herd. We also demonstrate that, at time scales relevant to infectious pathogens found in nature, the observed level of connectivity affects pathogen invasion likelihoods for a wide range of infectious periods and transmissibilities. Ultimately, our study identifies key predictors of social connectivity in a social ungulate and illustrates how contact heterogeneity, even within a highly connected herd, can shape pathogen invasion likelihoods.
Download This Paper Open PDF in Browser Add Paper to My Library Share: Permalink Using these links will ensure access to this page indefinitely Copy URL Copy DOI
Experimental data related to manuscript "Consilience in disease ecology: FMDV transmission dynamics reflect viral growth and immune response rates."SAA and Haptoglobin used as measures of innate response [micrograms/mL], FMDV viral data [log10 genome copies/mL], Virus neutralization titre used as measure of adaptive response. Interferon alpha and gamma.
Mathematical models of vector-borne infections, including malaria, often assume age-independent mortality rates of vectors, despite evidence that many insects senesce. In this study we present survival data on insecticide-resistant Anopheles gambiae s.l . from experiments in Côte d’Ivoire. We fit a constant mortality function and two age-dependent functions (logistic and Gompertz) to the data from mosquitoes exposed (treated) and not exposed (control) to insecticide-treated nets (ITNs), to establish biologically realistic survival functions. This enables us to explore the effects of insecticide exposure on mosquito mortality rates, and the extent to which insecticide resistance might impact the effectiveness of ITNs. We investigate this by calculating the expected number of infectious bites a mosquito will take in its lifetime, and by extension the vectorial capacity. Our results show that the predicted vectorial capacity is substantially lower in mosquitoes exposed to ITNs, despite the mosquitoes in the experiment being highly insecticide-resistant. The more realistic age-dependent functions provide a better fit to the experimental data compared to a constant mortality function and, hence, influence the predicted impact of ITNs on malaria transmission potential. In models with age-independent mortality, there is a great reduction for the vectorial capacity under exposure compared to no exposure. However, the two age-dependent functions predicted an even larger reduction due to exposure, highlighting the impact of incorporating age in the mortality rates. These results further show that multiple exposures to ITNs had a considerable effect on the vectorial capacity. Overall, the study highlights the importance of including age dependency in mathematical models of vector-borne disease transmission and in fully understanding the impact of interventions.
By mid-May 2020, cases of COVID-19 in the UK had been declining for over a month; a multi-phase emergence from lockdown was planned, including a scheduled partial reopening of schools on 1 June 2020. Although evidence suggests that children generally display mild symptoms, the size of the school-age population means the total impact of reopening schools is unclear. Here, we present work from mid-May 2020 that focused on the imminent opening of schools and consider what these results imply for future policy. We compared eight strategies for reopening primary and secondary schools in England. Modifying a transmission model fitted to UK SARS-CoV-2 data, we assessed how reopening schools affects contact patterns, anticipated secondary infections and the relative change in the reproduction number, R. We determined the associated public health impact and its sensitivity to changes in social distancing within the wider community. We predicted that reopening schools with half-sized classes or focused on younger children was unlikely to push R above one. Older children generally have more social contacts, so reopening secondary schools results in more cases than reopening primary schools, while reopening both could have pushed R above one in some regions. Reductions in community social distancing were found to outweigh and exacerbate any impacts of reopening. In particular, opening schools when the reproduction number R is already above one generates the largest increase in cases. Our work indicates that while any school reopening will result in increased mixing and infection amongst children and the wider population, reopening schools alone in June 2020 was unlikely to push R above one. Ultimately, reopening decisions are a difficult trade-off between epidemiological consequences and the emotional, educational and developmental needs of children. Into the future, there are difficult questions about what controls can be instigated such that schools can remain open if cases increase. This article is part of the theme issue 'Modelling that shaped the early COVID-19 pandemic response in the UK'.
Foot-and-mouth disease (FMD) is a virulent and economically important disease of livestock, still endemic in many areas of Asia and sub-Saharan Africa. Transmission from persistently infected livestock, also known as carriers, has been proposed as a mechanism to support the persistence of FMD in endemic regions. However, whether carrier livestock can infect susceptible animals is controversial; recovered virus is infectious and there are claims of field transmission, but it remains undemonstrated experimentally. Alternate hypotheses for persistence include the movement of livestock within and between regions, and fomite contamination of the environment. Using a stochastic compartmental ordinary differential equation (ODE) model, we investigate the minimum rates of carrier transmission necessary to contribute to the maintenance of FMD in a region, and compare this to the alternate mechanism of persistence through cattle shipments. We find that carrier transmission can theoretically support persistence even at transmission rates much lower than the highest realistic rates previously proposed, and that the parameters with the most effect on the feasibility of carrier-mediated persistence are the average duration of both the carrier phase and natural immunity. However, shipment-mediated persistence remains a viable alternate mechanism for persistence without carrier transmission.
Shuni virus is associated with neurologic and febrile illness in animals and humans. To determine potential vectors, we collected mosquitoes in South Africa and detected the virus in species of the genera Mansonia, Culex, Aedes, and Anopheles. These mosquitoes may be associated with Shuni virus outbreaks in Africa and emergence in other regions.
Background The rising incidence of visceral leishmaniasis due to Leishmania infantum requires novel methods to control transmission by the sand fly vector. Indoor residual spraying of insecticide (IRS) against these largely exophilic / exophagic vectors may not be the most effective method. A synthetic copy of the male sex-aggregation pheromone of the key vector species Lutzomyia longipalpis in the Americas, was co-located with residual pyrethroid insecticide, and tested for its effects on vector abundance, hence potential transmission, in a Brazilian community study. Methods Houses within eight defined semi-urban blocks in an endemic municipality in Brazil were randomised to synthetic pheromone + insecticide or to placebo treatments. A similar number of houses located >100m from each block were placebo treated and considered as “True Controls” (thus, analysed as three trial arms). Insecticide was sprayed on a 2.6m 2 surface area of the property boundary or outbuilding wall, co-located within one metre of 50mg synthetic pheromone in controlled-release dispensers. Vector numbers captured in nearby CDC light traps were recorded at monthly intervals over 3 months post intervention. Recruited sentinel houses under True Control and pheromone + insecticide treatments were similarly monitored at 7–9 day intervals. The intervention effects were estimated by mixed effects negative binomial models compared to the True Control group. Results Dose-response field assays using 50mg of the synthetic pheromone captured a mean 4.8 (95% C.L.: 3.91, 5.80) to 6.3 (95% C.L.: 3.24, 12.11) times more vectors (female Lu . longipalpis ) than using 10mg of synthetic pheromone. The intervention reduced household female vector abundance by 59% (C.L.: 48.7, 66.7%) (IRR = 0.41) estimated by the cross-sectional community study, and by 70% (C.L.: 56.7%, 78.8%) estimated by the longitudinal sentinel study. Similar reductions in male Lu . longipalpis were observed. Beneficial spill-over intervention effects were also observed at nearby untreated households with a mean reduction of 24% (95% C.L.: 0.050%, 39.8%) in female vectors. The spill-over effect in untreated houses was 44% (95% C.L.: 29.7%, 56.1%) as effective as the intervention in pheromone-treated houses. Ownership of chickens increased the intervention effects in both treated and untreated houses, attributed to the suspected synergistic attraction of the synthetic pheromone and chicken kairomones. The variation in IRR between study blocks was not associated with inter-household distances, household densities, or coverage (proportion of total households treated). Conclusions The study confirms the entomological efficacy of the lure-and-kill method to reduce the abundance of this important sand fly vector in treated and untreated homesteads. The outcomes were achieved by low coverage and using only 1–2% of the quantity of insecticide as normally required for IRS, indicating the potential cost-effectiveness of this method. Implications for programmatic deployment of this vector control method are discussed.
Mosquitoes in the Aedes and Culex genera are considered the main vectors of pathogenic flaviviruses worldwide. Entomological surveillance using universal flavivirus sets of primers in mosquitoes can detect not only pathogenic viruses but also insect-specific ones. It is hypothesized that insect-specific flaviviruses, which naturally infect these mosquitoes, may influence their vector competence for zoonotic arboviruses. Here, entomological surveillance was performed between January 2014 and May 2018 in five different provinces in the northeastern parts of South Africa, with the aim of identifying circulating flaviviruses. Mosquitoes were sampled using different carbon dioxide trap types. Overall, 64,603 adult mosquitoes were collected, which were screened by RT-PCR and sequencing. In total, 17 pools were found positive for insect-specific Flaviviruses in the mosquito genera Aedes (12/17, 70.59%) and Anopheles (5/17, 29.41%). No insect-specific viruses were detected in Culex species. Cell-fusing agent viruses were detected in Aedes aegypti and Aedes caballus. A range of anopheline mosquitoes, including Anopheles coustani, An. squamosus and An. maculipalpis, were positive for Culex flavivirus-like and Anopheles flaviviruses. These results confirm the presence of insect-specific flaviviruses in mosquito populations in South Africa, expands their geographical range and indicates potential mosquito species as vector species.
Extremely contagious pathogens are a global biosecurity threat because of their high burden of morbidity and mortality, as well as their capacity for fast-moving epidemics that are difficult to quell. Understanding the mechanisms enabling persistence of highly transmissible pathogens in host populations is thus a central problem in disease ecology. Through a combination of experimental and theoretical approaches, we investigated how highly contagious foot-and-mouth disease viruses persist in the African buffalo, which serves as their wildlife reservoir. We found that viral persistence through transmission among acutely infected hosts alone is unlikely. However, the inclusion of occasional transmission from persistently infected carriers reliably rescues the most infectious viral strain from fade-out. Additional mechanisms such as antigenic shift, loss of immunity, or spillover among host populations may be required for persistence of less transmissible strains.