Next Article in Journal
SARS-CoV-2 Variants and Clinical Outcomes of Special Populations: A Scoping Review of the Literature
Next Article in Special Issue
Increased Virulence of Culicoides Midge Cell-Derived Bluetongue Virus in IFNAR Mice
Previous Article in Journal
Poly(rC)-Binding Protein 2 Does Not Directly Participate in HCV Translation or Replication, but Rather Modulates Genome Packaging
Previous Article in Special Issue
Reference Material Production and Milk Protein Concentration as Elements to Improve Bluetongue Serological Diagnosis in Bulk Tank Milk
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling the Influence of Climate and Vector Control Interventions on Arbovirus Transmission

by
Emma L. Fairbanks
1,*,
Janet M. Daly
2 and
Michael J. Tildesley
1
1
The Zeeman Institute for Systems Biology & Infectious Disease Epidemiology Research, Mathematics Institute and School of Life Sciences, University of Warwick, Coventry CV4 7AL, UK
2
One Virology—Wolfson Centre for Global Virus Research, School of Veterinary Medicine and Science, University of Nottingham, Loughborough LE12 5RD, UK
*
Author to whom correspondence should be addressed.
Viruses 2024, 16(8), 1221; https://doi.org/10.3390/v16081221
Submission received: 13 June 2024 / Revised: 16 July 2024 / Accepted: 18 July 2024 / Published: 30 July 2024

Abstract

:
Most mathematical models that assess the vectorial capacity of disease-transmitting insects typically focus on the influence of climatic factors to predict variations across different times and locations, or examine the impact of vector control interventions to forecast their potential effectiveness. We combine features of existing models to develop a novel model for vectorial capacity that considers both climate and vector control. This model considers how vector control tools affect vectors at each stage of their feeding cycle, and incorporates host availability and preference. Applying this model to arboviruses of veterinary importance in Europe, we show that African horse sickness virus (AHSV) has a higher peak predicted vectorial capacity than bluetongue virus (BTV), Schmallenberg virus (SBV), and epizootic haemorrhagic disease virus (EHDV). However, AHSV has a shorter average infectious period due to high mortality; therefore, the overall basic reproduction number of AHSV is similar to BTV. A comparable relationship exists between SBV and EHDV, with both viruses showing similar basic reproduction numbers. Focusing on AHSV transmission in the UK, insecticide-treated stable netting is shown to significantly reduce vectorial capacity of Culicoides, even at low coverage levels. However, untreated stable netting is likely to have limited impact. Overall, this model can be used to consider both climate and vector control interventions either currently utilised or for potential use in an outbreak, and could help guide policy makers seeking to mitigate the impact of climate change on disease control.

1. Introduction

Vectorial capacity, defined as the total number of potentially infectious bites that would arise from all the vectors biting a single infectious host on a single day [1], is a concept used in epidemiology to measure the potential ability of a vector species to transmit a pathogen to a susceptible host population. There are many models for vectorial capacity, each considering different bionomic parameters [2,3,4,5].
The rate at which vectors complete the gonotrophic cycle, defined as the time required for locating a host, blood feeding, egg maturation and oviposition, has been found to be temperature-dependent [6,7,8,9,10]. The extrinsic incubation period (EIP), defined as the time required for a vector to become infectious after consumption of an infected blood meal, is also temperature-dependent [11,12]. Multiple models consider how vectorial capacity changes with different climatic conditions [13,14,15,16,17,18,19,20,21]. These are usually used to assess variations in risk geographically and/or for climate change scenarios. Brand and Keeling [13] suggested a model for vectorial capacity of bluetongue virus (BTV), a Culicoides-borne orbivirus of veterinary significance, which considered the temperature dependence of the gonotrophic cycle length, vector mortality, and the EIP.
Chitnis et al. [22] suggested a vectorial capacity model for malaria transmission derived from a discrete-time entomological model of the Anopheles feeding cycle. Updated versions of this model have been utilised to predict the community-level effect of using a new vector tool on Plasmodium transmission [23,24,25]. The model considers the effects of vector-control interventions at different stages during the mosquito feeding cycle. Although most of the parameter values of the model are measurable, several, such as the ovary sac proportion (proportion of mosquitoes with uncontracted ovary sacs) require technical time-consuming experiments; therefore, there are limited data available on these parameters. One advantage of this model is that it incorporates host selection, which is influenced by the availability of competent and non-competent hosts and the vector’s preferences for these hosts. This consideration is quantified using the human blood index, defined as the proportion of vectors that have fed on humans relative to the total number of blood-fed vectors analysed.
We suggest an updated version of the Brand and Keeling [13] model that considers the influence of climate on vectors, as well as features from the Chitnis et al. [22] model, namely host selection and the impact of vector control tools at various stages of the gonotrophic cycle. This model considers both personal protection through reduction in feeding, as well as community protection through vector mortality and disarming (feeding inhibition due to intoxication). This model is then parameterised for Culicoides-borne arboviruses of veterinary importance in Europe: BTV, African horse sickness virus (AHSV), Schmallenberg virus (SBV), and epizootic haemorrhagic disease virus (EHDV). This model is then applied to temperature data from the United Kingdom, which has previously seen BTV and SBV outbreaks. The potential reduction in vectorial capacity for AHSV is then calculated for insecticide-treated stable netting, dependent on behaviour of owners.
Feeding multiple times per gonotrophic cycle is observed to occur naturally for some vectors (for example, Aedes mosquitoes), increasing opportunities for disease transmission [26]. This behaviour has not been investigated for Culicoides. In other cases, when vectors are disturbed during a blood meal, they may seek more blood from the same or another host [1]. This model considers the influence of the average number of blood meals per gonotrophic cycle on the predicted vectorial capacity and the impact of vector-control tools.

2. Methods

2.1. Model Development

Parameters for the model are given in Table 1.
This model, similarly to the Brand and Keeling [13] model, allows for several variables to be dependent on climate. Climatic-dependent variables will be considered as vectors, with each entry representing a predicted daily value, calculated from climatic data. Included in these variables are the vector gonotrophic cycle completion rate, vector mortality rate, and the pathogen incubation rate. Alternatively, following the approach of the Chitnis et al. [22] model, these variables can be treated as constants in regions with relatively homogeneous climates or when data are insufficient. Under these conditions, all elements of the vectors would be identical.
Given that a pathogen is transmitted to a vector on day t, the probabilities that the vector survives until and is infectious on day τ are
P μ ( t , τ ) = exp i = t τ μ ^ ( i ) and
P σ ( t , τ ) = Γ k j = t + 1 τ σ ^ ( j )
respectively, where Γ k is the cumulative distribution function of a gamma distribution with mean 1 and shape k. For the survival probability (Equation (1)), the accumulated rate includes the infection day, because the vector must survive that day to transmit the pathogen to other hosts. For the probability of the vector completing the EIP (Equation (3)), the accumulated rate starts on the day after infection. This is because most bites occur at sunset and twilight [27]; therefore, we assume the EIP only depends on climatic conditions during days after the bite has occurred. It is important to note that bites occurring in the early hours of the morning are typically counted as part of the previous day’s activity. The temperature during the day affects whether Culicoides will be active and bite that night. In this study, we assume that k = 1 ; therefore, we have
P σ ( t , τ ) = 1 exp j = t + 1 τ σ ^ ( j ) ,
which is a Markovian model.
In the model, when vectors encounter a host, they either bite, are killed preprandially (before feeding), are disarmed, or seek another host. If they bite, they may be killed postprandially (after feeding). Here, we refer to biting, preprandial mortality, and disarming as host-encountering events. The probability of a host-countering event on day s is considered to be Markovian, calculated as
P α ( s ) = 1 exp ϵ α ^ ( s ) χ ψ ( 1 ϕ ( π κ M κ D ) ) + ( 1 ψ ) + ( 1 χ ) ,
where parameters are defined in Table 1. Here, the impact of vector-control tools is scaled by their usage ( ϕ ), coverage ( ψ ) and the proportion of the vectors that would bite the host type in the absence of the vector-control tool. In this study, we assume that vector control tools are applied only to competent hosts. Vectors that only feed once during their gonotrophic cycle stop host seeking after host-encountering events.
The probability a bite with transmission of a pathogen from a host to a vector and vector to host on days t and τ are
P h v ( t ) = P α ( t ) ρ h v χ ψ ( 1 ϕ π ) + ( 1 ψ ) χ ψ ( 1 ϕ ( π κ M κ D ) ) + ( 1 ψ ) + ( 1 χ )
and
P v h ( t , τ ) = P α ( τ ) P σ ( t , τ ) ρ v h χ ψ ( 1 ϕ π ) + ( 1 ψ ) χ ψ ( 1 ϕ ( π κ M κ D ) ) + ( 1 ψ ) + ( 1 χ ) ,
respectively. Transmission does not occur if the vector is preprandially killed or disarmed.
The rates of the preprandial and postprandial killing effects of tools on any given day s are calculated as
γ ^ p r e ( s ) = P α ( s ) χ ψ ϕ κ M χ ψ ( 1 ϕ ( π κ M κ D ) ) + ( 1 ψ ) + ( 1 χ ) and
γ ^ p o s t ( s ) = P α ( s ) χ ψ ( 1 ϕ π ) ξ χ ψ ( 1 ϕ ( π κ M κ D ) ) + ( 1 ψ ) + ( 1 χ ) ,
respectively. The probability of these events are considered to be Markovian; therefore, the probabilities of a vector being killed preprandially or postprandially by a tool before a host-encountering event on day τ given it feeds on day t are
P κ M ( t , τ ) = 1 exp m = t + 1 τ γ ^ p r e ( m ) and
P ξ ( t , τ ) = 1 exp m = t τ 1 γ ^ p o s t ( m ) ,
respectively. Here, the vector can be killed preprandially while encountering the host on day τ .
The vectorial capacity on day t is then calculated as the probability the vector feeds on a host on day t and survives until the end of the day, multiplied by the probability it infects another host every day after, given it survives until that day:
V C ( t ) = P h v ( t ) P μ ( t , t ) 1 P ξ ( t , t ) × x = t + 1 P v h ( t , x ) P μ ( t , x ) 1 P κ M ( t , x ) 1 P ξ ( t , x 1 ) .

2.2. Model Application: Arbovirus Transmission and Control in the UK

2.2.1. Parameterising the Model for Culicoides-Borne Viruses

We parameterise the model for BTV, AHSV, EHDV, and SBV transmission in the UK facilitated by C. obsoletus, the same vector considered to transmit BTV in the Brand and Keeling [13] model. Here, the model will have a daily time step. Since we consider transmission by the same vector for all viruses, the parameters for the gonotrophic cycle completion rate and vector mortality rate are the same. The rate of gonotrophic cycle completion was parameterised by Mullens et al. [9] as
α ^ ( s ) = max 0 , 1.9 × 10 4 T ( s ) ( T ( s ) 3.7 ) ( 41.9 T ( s ) ) 1 / 2.7 ,
where T is the temperature in degrees Celsius on day s. For the rate of mortality, we use the reciprocal of the vector lifespan
μ ^ ( s ) = 1 111.84 exp ( 0.1547 T ( s ) ) ,
parameterised by Gerry and Mullens [28].
Carpenter et al. [12] parameterised a model for the EIP completion rate for BTV, AHSV and EHDV, as well as the probability of transmission from host to vector. This model assumes that
σ ^ ( s ) = max 0 , α ( T ( s ) T m i n ) ,
where T m i n is the minimum threshold temperature that virus replication occurs and α is the rate of viral replication per degree-day above this threshold. The model was parameterised for five BTV studies. The model was only parameterised for C. sonorensis for AHSV and EHDV, and since there is no parameterisation for C. obsoletus, we use the mean of the three C. sonorensis experiments. These parameter values, along with parameters for the other virus, are given in Table 2.
Gubbins et al. [29] estimated α , T m i n and ρ h v for SBV using statistical inference for a within-farm model, utilising seroprevalence data for cattle and sheep farms in Belgium and The Netherlands. In this study, the probability of transmission from host to vector was ρ h v = 0.14 . More recently, the susceptibility of field caught Culicoides in the UK was investigated by Barber et al. [30]. This study found that 12% of Culicoides contained significant quantities of SBV RNA within their heads 8 days post-feeding on viraemic blood through an artificial membrane. However, this value depends on the cycle threshold selected when analysing polymerase chain reaction data. For a more conservative cut-off level, this estimation decreased to 7%. We used ρ h v = 0.12 , reflecting the less conservative value since it is closer to the Gubbins et al. [29] value, which was jointly estimated with other model parameters, and the values for other viruses in this study.
The parameters for the probability of transmission from vector to host ( p v h ) and infectious period for each virus and host species are given in Table 3. There is limited data to inform this parameter for EHDV. Ruder et al. [31] experimentally infected two calves with EHDV by allowing infectious C. sonorensis to feed on them. As both animals became infectious, we assume that the probability of vector-to-host transmission is high, and assign the value estimated for BTV ( p h v 0.9), as this is the highest value estimated amongst the other viruses in this study.

2.2.2. Climate Data

The mean and maximum daily temperatures for the years 1973–2022 were extracted from the Hadley Centre Central England Temperature (HadCET) dataset [37], which are representative of a triangular region of the UK enclosed by Lancashire, London, and Bristol. The temperature-dependent rates were calculated for each day using the mean daily temperature.

2.2.3. Comparing Vectorial Capacity of Viruses

First, we will consider scenarios with no vector-control tools present. In this case, the parameters π , κ M , κ D and ξ are set to 0. In this scenario, the probability the vector feeds is
P α ( τ ) = 1 exp ( ϵ α ^ ( τ ) ) .
The equation for vectorial capacity also simplifies to
V C ( t ) = P h v ( t ) P μ ( t , t ) × x = t + 1 P v h ( x ) P μ ( t , x ) ,
which, if the blood index is set to 1, reproduces results from the Brand and Keeling [13] model.
Host selection is informed by a systematic literature search for the sources of blood meals in Culicoides [38]. Here, we use the average percentage of blood meals which were from the competent host, given the host was present at a location. For the C. obsoletus complex, the average for these values in Europe were 0.42, 0.56, and 0.41, for cattle, horse, and deer, respectively. This average was calculated as the mean weighted by the number of Culicoides caught per location.
The temperature data are divided into 10-year blocks: 1973–1982, 1983–1992, 1993–2002, 2003–2012 and 2013–2022. The vectorial capacity is calculated for each year. We then compare how the vectorial capacity has changed as the climate has changed according to the mean for each of the 10-year blocks.
For the 2013–2022 block, the mean number of days above the minimum threshold temperature that virus replication occurs ( T m i n ) is calculated using the minimum and maximum daily temperatures for each virus.

2.2.4. Uncertainty Analysis

The model simulated varying parameters individually to allow for the effects of perturbations in these parameters to be considered. Parameters considered varied in this analysis; their fixed values while varying other parameters and the ranges for which they are varied are given in Table 4. Here, the model is simulated for constant temperatures between 10 °C and 45 °C.

2.2.5. Comparing R 0 of Viruses

Brand and Keeling [13] suggested a formula for R 0 based on the vectorial capacity. Here, for each day t, assuming a host becomes infectious on day t, R 0 is calculated as the sum of the vectorial capacity over the expected infectious period multiplied by the expected number of vectors per host. However, Brand and Keeling [13] used estimations for the expected number of vectors per host derived from field locations with no host present [39]. These values may not reflect the number of vectors feeding on each host per day.
Due to a lack of knowledge on the number of vectors per host, in this study, we do not attempt to calculate R 0 . Instead, we calculate just the sum of the vectorial capacity over a host’s infectious period, without also multiplying by the expected number of vectors per host. The reciprocal of this value would predict the number of vectors per host required such that R 0 would be larger than 1. On each day t, we calculate the cumulative vectorial capacity from day t to day (t + host infectious period − 1). Here, it is important to note that we use vectorial capacity estimations that are scaled according to host selection; therefore, we are calculating the number of host-seeking vectors per competent host required, rather than the number of vectors required to bite each competent host.
Although the numbers of vectors in locations where each host type is found may vary; this method allows us to compare potential values of R 0 between viruses.

2.2.6. Feeding Multiple Times per Gonotrophic Cycle

We simulate the model for the years 2013–2022, assuming that the vector feeds 1, 2 or 3 times per gonotrophic cycle. Here, we consider host selection in the model.

2.2.7. Vector Control

Baker et al. [40] showed that untreated netting could reduce the number of Culicoides in the C. obsoletus group entering a horse stable by 75% ( π = 0.75 ). Mortality during World Health Organisation (WHO) cone bioassays across trialled net treatments varied substantially between days post-treatment and the treatments trialled. Effectiveness of interventions declined to less than 40% mortality after 7 days in most and 14 days for almost all interventions studied. We therefore simulate the model for 20% and 40% mortality induced by contact with the net. This is assumed to occur both on entrance and exit of the stable, and therefore is both a preprandial and postprandial killing effect.
We simulate the model for AHSV with 25%, 50% and 75% coverage ( ψ = 0.25 ,   0.50 ,   0.75 ) and 50% and 90% usage ( ϕ = 0.50 ,   0.90 ) for each proposed netting scenario in Table 5. As well as the varied mortality rates, we also simulate the model for an untreated net and a net which both kills and disarms. Previous studies have shown that as some insecticides age, their modes of action change, and they can disarm rather than kill [25].

3. Results

3.1. Model Development

Our study introduces a novel mathematical model for vectorial capacity that integrates both climatic factors and vector control interventions. This model builds upon existing frameworks by incorporating several key innovations. Firstly, it considers the impact of climate on various aspects of vector biology, including the gonotrophic cycle completion rate, vector mortality rate, and pathogen incubation rate. Secondly, it explicitly accounts for the effects of vector control measures at different stages of the vector feeding cycle, allowing for a more comprehensive assessment of intervention strategies. Thirdly, the model incorporates host availability and vector preferences, providing a more realistic representation of disease transmission dynamics. Finally, it allows for the consideration of vectors feeding multiple times per gonotrophic cycle, a behaviour observed in some vector species that can significantly impact disease transmission. This integrated approach enables a more comprehensive analysis of arbovirus transmission potential under varying climatic conditions and control scenarios, as demonstrated in the following results.

3.2. Comparing Vectorial Capacity of Viruses

Figure 1 shows the influence of the change in climate on the vectorial capacity. We observe that, for all diseases, the predicted vectorial capacity tends to be greater for more recent years. For EHDV, our calculations show that the vectorial capacity is only non-negligible for the most recent 2013–2022 block.
We observe that considering host selection reduces the vectorial capacity. Therefore, not considering this in models may lead to overestimating the ability of vectors to transmit a disease.
Overall, AHSV consistently has the largest vectorial capacity, compared to the other diseases, for both the non-scaled and scaled models. When comparing AHSV to BTV and SBV, this is due to the probability of transmission from host to vector being much more likely for AHSV. This probability is largest for EHDV; however, this virus has a much larger minimum threshold temperature for replication than the other viruses in this study. This is why there is a relatively short time-span where the vectorial capacity of EHDV is above 0 in the UK. However, here we have used the mean daily temperature for our calculations. There are likely to be days which are therefore assumed to have no virus replication because the mean temperature does not exceed the minimum threshold for virus replication, whereas in fact the maximum temperature does exceed this value. Table 6 shows that, for all viruses, there are more days with potential replication when the maximum daily temperature is considered, compared with the mean. This is particularly the case for EHDV, for which on average there were nearly seven times more days when the maximum temperature exceeded the replication threshold, compared to the mean temperature.

3.3. Uncertainty Analysis

Figure 2 shows the uncertainty analysis results. We observe that the model is not sensitive to small perturbations in the parameters describing the EIP completion rate ( T m i n and α ). Similarly, the model is not very sensitive to small changes in the probability of transmission from host to vector or vector to host ( ρ h v or ρ v h ). Therefore, the model’s outputs are likely to remain robust, even with minor inaccuracies in these parameter estimations. Although the same hold for moderate uncertainty around the blood index of susceptible hosts, the model is much more sensitive to this parameter. Within the range of the means for hosts in this study (0.41–0.56), small inaccuracies will not significantly effect the model output. However, this highlights the importance of considering the presents of other non-susceptible hosts when considering the capacity for transmission of these viruses.

3.4. Comparing R 0 of Viruses

Although the predicted vectorial capacity is found to be largest for AHSV when compared with BTV and EHDV, it has a relatively short host infectious period (due to high mortality rates). Figure 3 shows that when we look at the sum of the vectorial capacity over the infectious period of a host, given it became infectious on day t, the cumulative vectorial capacity over this period is comparable for AHSV and BTV. Results suggest that EHDV could also have the same cumulative vectorial capacity as SBV for a short period during the summer months in the UK.
If we consider the reciprocal of the sum of the vectorial capacity over the infectious period of the host, we note that these predictions suggest that around three host seeking vectors per host type of interest would result in the peak R 0 being above 1 for BTV and AHSV, whereas for SBV and EHDV, the prediction would be around 20 host-seeking vectors per host.

3.5. Feeding Multiple Times per Gonotrophic Cycle

Figure 4 demonstrates that, if vectors are feeding on hosts more than one time per gonotrophic cycle, the vectorial capacity of these viruses could be much larger. This increase is non-linear, since there is an increased probability each day of the vector both acquiring the virus and transmitting the virus once it becomes infectious.

3.6. Vector Control

Figure 5 shows the predicted reduction in vectorial capacity due to the four types of stable netting described in Table 5. For each netting type, the reduction in vectorial capacity is greater with higher coverage and usage/adherence. Predictions suggest that untreated netting could reduce vectorial capacity by approximately 70% with high coverage (75%) and usage (90%); however, all treated netting options are competitive with this, even at low coverage (25%) and usage (50%). With coverage or usage below 50%, untreated netting does not substantially reduce vectorial capacity.
Netting that induces 40% mortality on contact reduces vectorial capacity by over 80% for all coverage and usage scenarios. Netting that disarms 20% of vectors and kills 20% on contact does not significantly reduce the vectorial capacity compared to interventions that only kill 20% of vectors on contact. This suggests that the mode of action transitioning to disarming from mortality as the intervention ages, or as vectors become resistant, would not support continual reductions in vectorial capacity.
For vectors that feed multiple times per feeding cycle, netting with modes of action beyond repelling have increased reductions in vectorial capacity. This is because the vector potentially comes into contact with the netting more times per feeding cycle, increasing the probability of mortality before it becomes, and once it is, infectious.

4. Discussion

We developed a model for vectorial capacity that considers both climate and host selection behaviour of vectors and vector controls. This model can be applied to more accurately predict vectorial capacity with current, or potential future, vector control strategies. Although temperature was the only climatic effect considered in this study, this model allows for additional climatic variables to be considered. This study focused on the vectorial capacity of Culicoides-borne viruses; however, it can be utilised for other pathogens and vectors.
We calculated the vectorial capacity based on a singular estimate for mean temperature. However, temperature varies according to location. Even within relatively small areas, microclimates can affect the temperature and vectorial capacity [21,41,42]. Previous studies indicate that landscape can significantly influence disease seroprevalence [43]. However, this variation in disease seroprevalence across different landscapes could also be influenced by the varying abundance of vectors. This, in turn, may depend on the thermal preferences of the vectors [44].
When considering the mean daily temperature, the number of days above the minimum threshold temperature for virus replication within the vector can be much lower than for the daily maximum temperature, especially for EHDV. Models using the mean temperature may therefore underestimate the vectorial capacity and the potential for pathogen transmission. Future work could explore more nuanced methods to estimate the duration of potential virus transmission, such as accumulated degree days above the threshold temperature, to incorporate within-day temperature variation. However, as finer-scale temperature data are not commonly reported, future studies could focus on developing and validating proxy measures or statistical methods to estimate accumulated degree days from available daily minimum, maximum, and mean temperatures. Such approaches could significantly improve our ability to predict virus transmission potential across varying climatic conditions.
We showed that considering host selection, driven by host availability and preferences of vectors, in models can significantly reduce vectorial capacity predictions. Therefore, the availability and preferences for competent and non-competent hosts, the availability of which varies by location, should be considered when assessing risk and deciding where to implement vector control.
In this study, R 0 was not calculated due to uncertainty of the vector-landscape. Two viruses may have the same vectorial capacity but different values of R 0 due to a different number of vectors per host, which will be variable across location and host type [43,45]. Brand and Keeling [13] used the number of Culicoides caught in light-traps without hosts present to parameterise the model. However, light-traps have been shown to not be good surrogate hosts for Culicoides [45,46].
The peak vectorial capacity summed over the infectious period for AHSV and EHDV were similar to BTV and SBV, respectively. Given that the UK has previously seen outbreaks of BTV and SBV, this study suggests that AHSV and EHDV outbreaks could also be a threat. Recently, in 2024, recommencement of direct horse movements from South Africa to the European Union has been approved [47]. Our results suggest that an importation of an infected case could lead to further transmission. This study provided evidence that vectorial capacity of the viruses the model was applied to in this study has been increasing in the UK. This highlights the importance of preparedness for the potential emergence of these diseases in the UK and the rest of Europe.
In dividing the vectorial capacity estimates into ten-year blocks and calculating averages for each period, we aimed to observe long-term trends in vectorial capacity. However, this methodological choice introduces the potential for discontinuities between the blocks. Future analyses might benefit from applying smoothing techniques. Despite this, the model predicted that the vectorial capacity of EHDV was only significantly above 0 for the most recent block, ending in 2022. The first reported case of EHDV in Europe was 2022, and outbreaks have continued to occur [48,49].
Brand and Keeling [13] compared three models of vectorial capacity: deterministic, Markovian, and intermediate. The intermediate model requires knowledge of the distribution of the EIP by setting a shape parameter for the gamma distribution (k). Here, we simulated the Markovian model (by assuming k = 1 in Equation (2)), as k is not known for all viruses in the study. However, the model developed allows the EIP be modelled using the deterministic or intermediate model from Brand and Keeling [13]. If these other models are used, it is important to consider that the vectorial capacity predictions will generally be lowest for the deterministic model and largest for the Markovian model, with the intermediate model falling in between. Therefore, direct comparisons can not be made between results for different climates/viruses if a different model is used.
This study focused on adulticide vector-control tools. However, other tools, such as larvicide, have been shown to be effective against Culicoides [50]. Previously, the killing of larvae has been integrated into models of vectorial capacity by multiplying the vectorial capacity by the coverage of larvicide; however, larvicide does not affect vectorial capacity. By definition, vectorial capacity is the number of cases which arise from a single day of feeding on an infectious host, therefore vectors must have made it to the adult stage to influence vectorial capacity. Instead, larviciding impacts the number of vectors per host, hence reducing R 0 . If the Brand and Keeling [13] method is used to calculate R 0 from the vectorial capacity, proportionally reducing vectorial capacity or the number of vectors per host by larvicide coverage would yield the same result.
Usage/adherence is sometimes not considered in vector-control studies; however, this study shows that they can be just as important as coverage for reducing vectorial capacity. This is expected, since if a tool is unused, it will not have an effect. A drawback is that this parameter can be hard to quantify. Usage describes a reduction in efficacy due to behaviours such as not stabling the horse every night or whether the horse is outside the stable for some of the time vectors are host seeking. However, vector feeding behaviours are not consistent between nights or over the duration of a night [51].
The data used to parameterise the stable netting control strategy was gathered using WHO cone bioassays and light-traps within stables. Light traps have been observed to not effectively estimate protective efficacy for spatial repellents [52]. Considering this, and the relatively small number of vectors caught in the stable (n = 25), the estimation for the reduction in biting may be inaccurate. WHO cone bioassays have also been shown to induce more mortality compared to when the interventions are evaluated in semi-field or experimental-hut conditions [53,54]. In this study, we estimate the reduction in vectorial capacity for treated netting that disarms mosquitoes. However, it is important to note that this disarming behaviour was not directly observed during the trials that were used to parameterise the model. Despite these drawbacks, we were able to provide a framework for evaluating the potential implication of transmission for tools with a range of properties. We found that whilst untreated stable netting could need to have high coverage and usage to substantially reduce the vectorial capacity of AHSV in the UK, insecticide treated netting could potentially decrease vectorial capacity, even at low coverage and usage levels.
Other adulticides have been been shown to be effective against Culicoides for use on host species of relevance to the viruses in this study [55]. Robin et al. [56] found that there was no reduction in the percentage of Culicoides caught in traps which were blood-fed next to a horse treated with topical deltamethrin, compared to an untreated horse. However, the difference in the total number or number of blood-fed Culicoides was not reported. Killing the vectors before they feed (preprandial killing) would likely increase the percentage of blood-fed vectors caught, while disarming them would decrease it. In contrast, using repellents might not change this percentage, as the insects are still able to escape after encountering the repellent.
This study highlighted the potential added efficacy of vector control tools against mosquitoes that feed multiple times per gonotrophic cycle, compared to those who just feed once. There is evidence that successive blood meals after an infectious blood meal can increase vector competence and shorten the EIP of vectors [57,58,59,60,61]. Therefore, this behaviour could increase the infectiousness of vectors, which was not considered in this study. If feeding on multiple hosts were to occur due to disruption when feeding, the size, and therefore viraemic dose, of the blood meal may be less, which may also reduce the probability the vector becomes infectious [60].
We have established a novel model for vectorial capacity that considers the influence of both climate and vector control at each stage of the vector feeding cycle. We demonstrated how this model can be used to compare the vectorial capacity of viruses. This model was used to assess the potential for outbreaks by comparing the total vectorial capacity during a host’s infectious period between viruses that were previously unobserved in a region and those that had been observed there before. Extensions to this work could include global sensitivity analysis to explore the influence of tools across a range of vector bionomics or temperatures, using methods similar to [62] to identify target product profiles.
Our research provides valuable insights into how the strategic use of insecticide-treated netting, even at low coverage levels, can substantially reduce the transmission potential of AHSV. This model offers a powerful tool for policymakers and health professionals, aiding in the formulation of more effective vector management strategies that could mitigate the impact of these diseases, especially in the context of changing global climates.

Author Contributions

E.L.F.: Conceptualisation, Methodology, Software, Validation, Formal analysis, Investigation, Data Curation, Writing—Original Draft, Writing—Review and Editing, Visualisation, Funding acquisition; J.M.D.: Conceptualisation, Writing—Review and Editing, Funding acquisition; M.J.T.: Conceptualisation, Methodology, Writing—Review and Editing, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

The authors were supported by the Horserace Betting Levy Board (vet/prj/809). MJT is funded on a joint BBSRC/EEID grant (BB/T004312/1).

Institutional Review Board Statement

The authors confirm that the ethical policies of the journal, as noted on the journal’s author guidelines page, have been adhered to. No ethical approval was required.

Data Availability Statement

R code to run the model is available at https://github.com/emmafairbanks/EntoModels (DOI:10.5281/zenodo.12785165).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Garrett-Jones, C. The human blood index of malaria vectors in relation to epidemiological assessment. Bull. World Health Organ. 1964, 30, 241. [Google Scholar] [PubMed]
  2. Macdonald, G. The epidemiology and control of malaria. In The Epidemiology and Control of Malaria; Oxford University Press: Oxford, UK, 1957. [Google Scholar]
  3. Smith, D.; Ellis McKenzie, F. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar. J. 2004, 3, 13. [Google Scholar] [CrossRef] [PubMed]
  4. Brady, O.; Godfray, H.; Tatem, A.; Gething, P.; Cohen, J.; McKenzie, F.; Perkins, T.; Reiner, R.; Tusting, L.; Sinka, M.; et al. Vectorial capacity and vector control: Reconsidering sensitivity to parameters for malaria elimination. Trans. R. Soc. Trop. Med. Hyg. 2016, 110, 107–117. [Google Scholar] [CrossRef]
  5. Catano-Lopez, A.; Rojas-Diaz, D.; Laniado, H.; Arboleda-Sánchez, S.; Puerta-Yepes, M.; Lizarralde-Bejarano, D. An alternative model to explain the vectorial capacity using as example Aedes aegypti case in dengue transmission. Heliyon 2019, 5, e02577. [Google Scholar] [CrossRef] [PubMed]
  6. Alto, B.; Juliano, S. Precipitation and temperature effects on populations of Aedes albopictus (Diptera: Culicidae): Implications for range expansion. J. Med. Entomol. 2001, 38, 646–656. [Google Scholar] [CrossRef]
  7. Wittmann, E.; Mellor, P.; Baylis, M. Effect of temperature on the transmission of orbiviruses by the biting midge, Culicoides sonorensis. Med. Vet. Entomol. 2002, 16, 147–156. [Google Scholar] [CrossRef]
  8. Löwenberg Neto, P.; Navarro-Silva, M. Development, longevity, gonotrophic cycle and oviposition of Aedes albopictus Skuse (Diptera: Culicidae) under cyclic temperatures. Neotrop. Entomol. 2004, 33, 29–33. [Google Scholar] [CrossRef]
  9. Mullens, B.; Gerry, A.; Lysyk, T.; Schmidtmann, E. Environmental effects on vector competence and virogenesis of bluetongue virus in Culicoides: Interpreting laboratory data in a field context. Vet. Ital. 2004, 40, 160–166. [Google Scholar]
  10. Gunay, F.; Alten, B.; Ozsoy, E. Estimating reaction norms for predictive population parameters, age specific mortality, and mean longevity in temperature-dependent cohorts of Culex quinquefasciatus Say (Diptera: Culicidae). J. Vector Ecol. 2010, 35, 354–362. [Google Scholar] [CrossRef]
  11. Tjaden, N.; Thomas, S.; Fischer, D.; Beierkuhnlein, C. Extrinsic incubation period of dengue: Knowledge, backlog, and applications of temperature dependence. PLoS Negl. Trop. Dis. 2013, 7, e2207. [Google Scholar] [CrossRef]
  12. Carpenter, S.; Wilson, A.; Barber, J.; Veronesi, E.; Mellor, P.; Venter, G.; Gubbins, S. Temperature dependence of the extrinsic incubation period of orbiviruses in Culicoides biting midges. PLoS ONE 2011, 6, e27987. [Google Scholar] [CrossRef] [PubMed]
  13. Brand, S.; Keeling, M. The impact of temperature changes on vector-borne disease transmission: Culicoides midges and bluetongue virus. J. R. Soc. Interface 2017, 14, 20160481. [Google Scholar] [CrossRef] [PubMed]
  14. Jetten, T.; Martens, W.; Takken, W. Model simulations to estimate malaria risk under climate change. J. Med. Entomol. 1996, 33, 361–371. [Google Scholar] [CrossRef] [PubMed]
  15. Jetten, T.; Focks, D. Potential changes in the distribution of dengue transmission under climate warming. Am. J. Trop. Med. Hyg. 1997, 57, 285–297. [Google Scholar] [CrossRef] [PubMed]
  16. Bomblies, A.; Eltahir, E. Assessment of the impact of climate shifts on malaria transmission in the Sahel. EcoHealth 2009, 6, 426–437. [Google Scholar] [CrossRef] [PubMed]
  17. Ceccato, P.; Vancutsem, C.; Klaver, R.; Rowland, J.; Connor, S. A vectorial capacity product to monitor changing malaria transmission potential in epidemic regions of Africa. J. Trop. Med. 2012, 2012, 595948. [Google Scholar] [CrossRef] [PubMed]
  18. Yamana, T.; Eltahir, E. Projected impacts of climate change on environmental suitability for malaria transmission in West Africa. Environ. Health Perspect. 2013, 121, 1179–1186. [Google Scholar] [CrossRef] [PubMed]
  19. Shapiro, L.; Whitehead, S.; Thomas, M. Quantifying the effects of temperature on mosquito and parasite traits that determine the transmission potential of human malaria. PLoS Biol. 2017, 15, e2003489. [Google Scholar] [CrossRef]
  20. Viennet, E.; Mincham, G.; Frentiu, F.; Jansen, C.; Montgomery, B.; Harley, D.; Flower, R.; Williams, C.; Faddy, H. Epidemic potential for local transmission of Zika virus in 2015 and 2016 in Queensland, Australia. PLoS Curr. 2016, 8. [Google Scholar] [CrossRef]
  21. Wimberly, M.; Davis, J.; Evans, M.; Hess, A.; Newberry, P.; Solano-Asamoah, N.; Murdock, C. Land cover affects microclimate and temperature suitability for arbovirus transmission in an urban landscape. PLoS Negl. Trop. Dis. 2020, 14, e0008614. [Google Scholar] [CrossRef]
  22. Chitnis, N.; Smith, T.; Steketee, R. A mathematical model for the dynamics of malaria in mosquitoes feeding on a heterogeneous host population. J. Biol. Dyn. 2008, 2, 259–285. [Google Scholar] [CrossRef] [PubMed]
  23. Briët, O.; Impoinvil, D.; Chitnis, N.; Pothin, E.; Lemoine, J.; Frederic, J.; Smith, T. Models of effectiveness of interventions against malaria transmitted by Anopheles albimanus. Malar. J. 2019, 18, 263. [Google Scholar] [CrossRef] [PubMed]
  24. Denz, A.; Njoroge, M.; Tambwe, M.; Champagne, C.; Okumu, F.; van Loon, J.; Hiscox, A.; Saddler, A.; Fillinger, U.; Moore, S.; et al. Predicting the impact of outdoor vector control interventions on malaria transmission intensity from semi-field studies. Parasites Vectors 2021, 14, 64. [Google Scholar] [CrossRef]
  25. Fairbanks, E.; Saeung, M.; Pongsiri, A.; Vajda, E.; Wang, Y.; McIver, D.; Richardson, J.; Tatarsky, A.; Lobo, N.; Moore, S.; et al. Inference for entomological semi-field experiments: Fitting a mathematical model assessing personal and community protection of vector-control interventions. Comput. Biol. Med. 2024, 168, 107716. [Google Scholar] [CrossRef] [PubMed]
  26. Scott, T.; Clark, G.; Lorenz, L.; Amerasinghe, P.; Reiter, P.; Edman, J. Detection of multiple blood feeding in Aedes aegypti (Diptera: Culicidae) during a single gonotrophic cycle using a histologic technique. J. Med. Entomol. 1993, 30, 94–99. [Google Scholar] [CrossRef] [PubMed]
  27. Meiswinkel, R.; Elbers, A. The dying of the light: Crepuscular activity in Culicoides and impact on light trap efficacy at temperate latitudes. Med. Vet. Entomol. 2016, 30, 53–63. [Google Scholar] [CrossRef] [PubMed]
  28. Gerry, A.; Mullens, B. Seasonal abundance and survivorship of Culicoides sonorensis (Diptera: Ceratopogonidae) at a southern California dairy, with reference to potential bluetongue virus transmission and persistence. J. Med. Entomol. 2000, 37, 675–688. [Google Scholar] [CrossRef] [PubMed]
  29. Gubbins, S.; Turner, J.; Baylis, M.; Van der Stede, Y.; van Schaik, G.; Abrahantes, J.; Wilson, A. Inferences about the transmission of Schmallenberg virus within and between farms. Prev. Vet. Med. 2014, 116, 380–390. [Google Scholar] [CrossRef] [PubMed]
  30. Barber, J.; Harrup, L.; Silk, R.; Veronesi, E.; Gubbins, S.; Bachanek-Bankowska, K.; Carpenter, S. Blood-feeding, susceptibility to infection with Schmallenberg virus and phylogenetics of Culicoides (Diptera: Ceratopogonidae) from the United Kingdom. Parasites Vectors 2018, 11, 116. [Google Scholar] [CrossRef]
  31. Ruder, M.; Mead, D.; Stallknecht, D.; Kedmi, M.; Klement, E.; Brown, J.; Carter, D.; Howerth, E. Experimental infection of Holstein cows and calves with EHDV-7 and preliminary evaluation of different inoculation methods. Vet. Ital. 2015, 51, 289–299. [Google Scholar] [CrossRef]
  32. Jones, A.; Turner, J.; Caminade, C.; Heath, A.; Wardeh, M.; Kluiters, G.; Diggle, P.; Morse, A.; Baylis, M. Bluetongue risk under future climates. Nat. Clim. Change 2019, 9, 153–157. [Google Scholar] [CrossRef]
  33. Gubbins, S.; Carpenter, S.; Baylis, M.; Wood, J.; Mellor, P. Assessing the risk of bluetongue to UK livestock: Uncertainty and sensitivity analyses of a temperature-dependent model for the basic reproduction number. J. R. Soc. Interface 2007, 5, 363–371. [Google Scholar] [CrossRef] [PubMed]
  34. Backer, J.; Nodelijk, G. Transmission and control of African horse sickness in The Netherlands: A model analysis. PLoS ONE 2011, 6, e23066. [Google Scholar] [CrossRef]
  35. Ruder, M.; Allison, A.; Stallknecht, D.; Mead, D.; McGraw, S.; Carter, D.; Kubiski, S.; Batten, C.; Klement, E.; Howerth, E. Susceptibility of white-tailed deer (Odocoileus virginianus) to experimental infection with epizootic hemorrhagic disease virus serotype 7. J. Wildl. Dis. 2012, 48, 676–685. [Google Scholar] [CrossRef]
  36. Ruder, M.; Stallknecht, D.; Allison, A.; Mead, D.; Carter, D.; Howerth, E. Host and potential vector susceptibility to an emerging orbivirus in the United States: Epizootic hemorrhagic disease virus serotype 6. Vet. Pathol. 2016, 53, 574–584. [Google Scholar] [CrossRef] [PubMed]
  37. Parker, D.; Legg, T.; Folland, C. A new daily central England temperature series, 1772–1991. Int. J. Climatol. 1992, 12, 317–342. [Google Scholar] [CrossRef]
  38. Fairbanks, E.; Tildesley, M.; Daly, J. A systematic review quantifying host feeding patterns of Culicoides species responsible for pathogen transmission. bioRxiv 2024. [Google Scholar] [CrossRef]
  39. Sanders, C.; Shortall, C.; Gubbins, S.; Burgin, L.; Gloster, J.; Harrington, R.; Reynolds, D.; Mellor, P.; Carpenter, S. Influence of season and meteorological parameters on flight activity of Culicoides biting midges. J. Appl. Ecol. 2011, 48, 1355–1364. [Google Scholar] [CrossRef]
  40. Baker, T.; Carpenter, S.; Gubbins, S.; Newton, R.; Lo Iacono, G.; Wood, J.; Harrup, L. Can insecticide-treated netting provide protection for Equids from Culicoides biting midges in the United Kingdom? Parasites Vectors 2015, 8, 604. [Google Scholar] [CrossRef]
  41. Haider, N.; Kirkeby, C.; Kristensen, B.; Kjær, L.; Sørensen, J.; Bødker, R. Microclimatic temperatures increase the potential for vector-borne disease transmission in the Scandinavian climate. Sci. Rep. 2017, 7, 8175. [Google Scholar] [CrossRef]
  42. Haider, N.; Cuellar, A.; Kjær, L.; Sørensen, J.; Bødker, R. Microclimatic temperatures at Danish cattle farms, 2000–2016: Quantifying the temporal and spatial variation in the transmission potential of Schmallenberg virus. Parasites Vectors 2018, 11, 128. [Google Scholar] [CrossRef] [PubMed]
  43. Boyer, T.; Ward, M.; Singer, R. Climate, landscape, and the risk of orbivirus exposure in cattle in Illinois and western Indiana. Am. J. Trop. Med. Hyg. 2010, 83, 789–794. [Google Scholar] [CrossRef] [PubMed]
  44. Hug, D.; Stegmayer, R.; Blanckenhorn, W.; Verhulst, N. Thermal preference of adult mosquitoes (Culicidae) and biting midges (Ceratopogonidae) at different altitudes in Switzerland. Med. Vet. Entomol. 2023, 37, 562–573. [Google Scholar] [CrossRef] [PubMed]
  45. Elbers, A.; Gonzales, J. Culicoides (Diptera: Ceratopogonidae) Abundance Is Influenced by Livestock Host Species and Distance to Hosts at the Micro Landscape Scale. Insects 2023, 14, 637. [Google Scholar] [CrossRef] [PubMed]
  46. Gerry, A.; Monteys, V.; Vidal, J.; Francino, O.; Mullens, B. Biting rates of Culicoides midges (Diptera: Ceratopogonidae) on sheep in northeastern Spain in relation to midge capture using UV light and carbon dioxide-baited traps. J. Med. Entomol. 2014, 46, 615–624. [Google Scholar] [CrossRef] [PubMed]
  47. Equine Infectious Disease Surveilence. Equine Quarterly Surveillance Report (Q1 2024). 2024. Available online: https://equinesurveillance.org/landing/resources/reports/dsr20241.pdf (accessed on 15 July 2024).
  48. Department for Environment, Food and Rural Affairs. Epizootic Haemorrhagic Disease in Europe: Preliminary Outbreak Assessment. 2022. Available online: https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1124406/Epizootic_haemorrhagic_disease_in_europe.pdf (accessed on 15 July 2024).
  49. Department for Environment, Food and Rural Affairs. Epizootic Haemorrhagic Disease in Europe: Updated Outbreak Assessment #2. 2023. Available online: https://assets.publishing.service.gov.uk/media/6500653f1886eb000d9771eb/EHDv_Europe_30_Aug_23.pdf (accessed on 15 July 2024).
  50. Meloni, G.; Cossu, M.; Foxi, C.; Vento, L.; Circosta, S.; Burrai, E.; Masala, S.; Goffredo, M.; Satta, G. Combined larvicidal and adulticidal treatments to control Culicoides biting midges (Diptera: Ceratopogonidae): Results of a pilot study. Vet. Parasitol 2018, 257, 28–33. [Google Scholar] [CrossRef]
  51. Elbers, A.; Meiswinkel, R. Culicoides (Diptera: Ceratopogonidae) and livestock in the Netherlands: Comparing host preference and attack rates on a Shetland pony, a dairy cow, and a sheep. J. Vector Ecol. 2015, 40, 308–317. [Google Scholar] [CrossRef] [PubMed]
  52. Swai, J.; Kibondo, U.; Ntabaliba, W.; Ngoyani, H.; Makungwa, N.; Mseka, A.; Chura, M.; Mascari, T.; Moore, S. CDC light traps underestimate the protective efficacy of an indoor spatial repellent against bites from wild Anopheles arabiensis mosquitoes in Tanzania. Malar. J. 2023, 22, 141. [Google Scholar] [CrossRef]
  53. Gebremariam, B.; Birke, W.; Zeine, W.; Ambelu, A.; Yewhalaw, D. Evaluation of Long-Lasting Insecticidal Nets (DuraNet®) Under laboratory and Semi-Field Conditions Using Experimental Huts Against Anopheles Mosquitoes in Jimma Zone, Southwestern Ethiopia. Environ. Health Insights 2021, 15, 1178630220974730. [Google Scholar] [CrossRef] [PubMed]
  54. Bayili, K.; Ki, H.; Bayili, B.; Sow, B.; Ouattara, A.; Small, G.; Hien, A.; Dabire, R.; Diabate, A. Laboratory and experimental hut trial evaluation of VECTRON™ T500 for indoor residual spraying (IRS) against insecticide resistant malaria vectors in Burkina Faso. Gates Open Res. 2022, 6, 57. [Google Scholar] [CrossRef]
  55. Mullens, B.; Velten, R.; Gerry, A.; Braverman, Y.; Endris, R. Feeding and survival of Culicoides sonorensis on cattle treated with permethrin or pirimiphos-methyl. Med. Vet. Entomol. 2000, 14, 313–320. [Google Scholar] [CrossRef] [PubMed]
  56. Robin, M.; Archer, D.; McGowan, C.; Garros, C.; Gardes, L.; Baylis, M. Repellent effect of topical deltamethrin on blood feeding by Culicoides on horses. Vet. Rec. 2015, 176, 574. [Google Scholar] [CrossRef] [PubMed]
  57. Dong, S.; Balaraman, V.; Kantor, A.; Lin, J.; Grant, D.; Held, N.; Franz, A. Chikungunya virus dissemination from the midgut of Aedes aegypti is associated with temporal basal lamina degradation during bloodmeal digestion. PLoS Negl. Trop. Dis. 2017, 11, e0005976. [Google Scholar] [CrossRef] [PubMed]
  58. Serafim, T.; Coutinho-Abreu, I.; Oliveira, F.; Meneses, C.; Kamhawi, S.; Valenzuela, J. Sequential blood meals promote Leishmania replication and reverse metacyclogenesis augmenting vector infectivity. Nat. Microbiol. 2018, 3, 548–555. [Google Scholar] [CrossRef] [PubMed]
  59. Armstrong, P.; Ehrlich, H.; Magalhaes, T.; Miller, M.; Conway, P.; Bransfield, A.; Misencik, M.; Gloria-Soria, A.; Warren, J.; Andreadis, T.; et al. Successive blood meals enhance virus dissemination within mosquitoes and increase transmission potential. Nat. Microbiol. 2020, 5, 239–247. [Google Scholar] [CrossRef] [PubMed]
  60. Rozo-Lopez, P.; Londono-Renteria, B.; Drolet, B. Impacts of infectious dose, feeding behavior, and age of Culicoides sonorensis biting midges on infection dynamics of vesicular stomatitis virus. Pathogens 2021, 10, 816. [Google Scholar] [CrossRef] [PubMed]
  61. Johnson, R.; Stopard, I.; Byrne, H.; Armstrong, P.; Brackney, D.; Lambert, B. Investigating the dose-dependency of the midgut escape barrier using a mechanistic model of within-mosquito dengue virus population dynamics. PLoS Pathog. 2024, 20, e1011975. [Google Scholar] [CrossRef]
  62. Wang, Y.; Chitnis, N.; Fairbanks, E. Optimizing malaria vector control in the Greater Mekong Subregion: A systematic review and mathematical modelling study to identify desirable intervention characteristics. Parasites Vectors 2024, 17, 162. [Google Scholar] [CrossRef]
Figure 1. The predicted vectorial capacity for bluetongue virus (BTV) in cattle, Schmallenberg virus (SBV) in cattle, African horse sickness virus (AHSV) in horses, and epizootic haemorrhagic disease virus (EHDV) in deer; (Top) not scaled for host selection, (Bottom) scaled for host selection.
Figure 1. The predicted vectorial capacity for bluetongue virus (BTV) in cattle, Schmallenberg virus (SBV) in cattle, African horse sickness virus (AHSV) in horses, and epizootic haemorrhagic disease virus (EHDV) in deer; (Top) not scaled for host selection, (Bottom) scaled for host selection.
Viruses 16 01221 g001
Figure 2. The vectorial capacity estimates using fixed parameter values and under constant temperatures (top left) and the uncertainty analysis outputs for each parameter analysed, as described in Table 4.
Figure 2. The vectorial capacity estimates using fixed parameter values and under constant temperatures (top left) and the uncertainty analysis outputs for each parameter analysed, as described in Table 4.
Viruses 16 01221 g002
Figure 3. Mean vectorial capacity predictions summed over the duration of the infectious period for bluetongue virus (BTV) in cattle, Schmallenberg virus (SBV) in cattle, African horse sickness virus (AHSV) in horses, and epizootic haemorrhagic disease virus (EHDV) in deer.
Figure 3. Mean vectorial capacity predictions summed over the duration of the infectious period for bluetongue virus (BTV) in cattle, Schmallenberg virus (SBV) in cattle, African horse sickness virus (AHSV) in horses, and epizootic haemorrhagic disease virus (EHDV) in deer.
Viruses 16 01221 g003
Figure 4. Mean vectorial capacity predictions for 2013–2022, considering host selection, for bluetongue virus (BTV), Schmallenberg virus (SBV), African horse sickness virus (AHSV), and epizootic haemorrhagic disease virus (EHDV), assuming vectors feed on average 1, 2, or 3 times per gonotrophic cycle.
Figure 4. Mean vectorial capacity predictions for 2013–2022, considering host selection, for bluetongue virus (BTV), Schmallenberg virus (SBV), African horse sickness virus (AHSV), and epizootic haemorrhagic disease virus (EHDV), assuming vectors feed on average 1, 2, or 3 times per gonotrophic cycle.
Viruses 16 01221 g004
Figure 5. Mean vectorial capacity predictions for 2013–2022, considering host selection, for bluetongue virus (BTV), Schmallenberg virus (SBV), African horse sickness virus (AHSV), and epizootic haemorrhagic disease virus (EHDV), assuming various usage and coverage for four types of netting. Top: vectors are assumed to feed once per gonotrophic cycle. Bottom: vectors are assumed to feed an average of three times per gonotrophic cycle.
Figure 5. Mean vectorial capacity predictions for 2013–2022, considering host selection, for bluetongue virus (BTV), Schmallenberg virus (SBV), African horse sickness virus (AHSV), and epizootic haemorrhagic disease virus (EHDV), assuming various usage and coverage for four types of netting. Top: vectors are assumed to feed once per gonotrophic cycle. Bottom: vectors are assumed to feed an average of three times per gonotrophic cycle.
Viruses 16 01221 g005
Table 1. Model parameter definitions. Parameters denoted with α ^ are vectors with elements of daily calculated variables dependent on climatic conditions. Above and below the filled line gives fundamental and composite parameters, respectively.
Table 1. Model parameter definitions. Parameters denoted with α ^ are vectors with elements of daily calculated variables dependent on climatic conditions. Above and below the filled line gives fundamental and composite parameters, respectively.
ParameterDefinition
α ^ ( s ) Rate of gonotrophic cycle completion on day s.
μ ^ ( s ) Vector mortality rate on day s.
σ ^ ( s ) Rate of pathogen EIP completion on day s.
χ Proportion of blood meals from competent hosts, referred to as the blood index.
ϵ Average number of bites per vector per gonotrophic cycle.
π Reduction in the rate of vector biting due to the presence of a vector-control tool.
κ M Increase in the rate of vector mortality before biting due to the presence of a vector-control tool, relative to the rate of biting without the vector-control tool.
κ D Rate of vector disarming due to the presence of a vector-control tool, relative to the rate of biting without the vector-control tool.
ξ Increased probability of vector mortality after biting due to the presence of a vector control tool.
ψ Proportion of the target hosts with access to the vector-control tool, referred to as coverage.
ϕ Adherence to using the vector-control tool, referred to as usage.
ρ h v Probability of transmission from host to vector, given that the host is infectious.
ρ v h Probability of transmission from vector to host, given that the vector is infectious.
P μ ( t , τ ) Probability a vector which fed on day t survives until day τ .
P σ ( t , τ ) Probability that a vector which fed on day t is infectious day τ , given transmission occurred during the bite on day t.
P α ( s ) Probability of a host-encountering event on day s (defined as feeding or preprandial mortality or disarming due to the presence of a tool).
P h v ( t ) Probability of transmission of a pathogen from host to vector on day t.
P v h ( τ ) Probability of transmission of a pathogen from vector to host on day τ .
γ ^ p r e ( s ) Rate of preprandial killing on day s.
γ ^ p o s t ( s ) Rate of postprandial killing on day s.
P κ M ( t , τ ) Probability a vector is killed preprandially by day τ given it being fed on day t.
P ξ ( t , τ ) Probability a vector is killed postprandially by day τ given it being fed on day t.
V C ( t ) The vectorial capacity on day t.
Table 2. The parameters for the minimum threshold temperature that virus replication occurs ( T m i n ), the rate of viral replication per degree-day above this threshold ( α ), and the probability of transmission from host to vector ( ρ h v = 0.14 ) for bluetongue virus (BTV), African horse sickness virus (AHSV), Schmallenberg virus (SBV), and epizootic haemorrhagic disease virus (EHDV).
Table 2. The parameters for the minimum threshold temperature that virus replication occurs ( T m i n ), the rate of viral replication per degree-day above this threshold ( α ), and the probability of transmission from host to vector ( ρ h v = 0.14 ) for bluetongue virus (BTV), African horse sickness virus (AHSV), Schmallenberg virus (SBV), and epizootic haemorrhagic disease virus (EHDV).
Virus T min α ρ h v Ref.
BTV12.60.0190.13 [12]
AHSV12.60.0170.52 [12]
SBV12.350.030.12 [29,30]
EHDV19.50.0840.92 [12]
Table 3. The parameters for the transmission probability from vector to host ( ρ h v ) and host infectious period for bluetongue virus (BTV), African horse sickness virus (AHSV), Schmallenberg virus (SBV), epizootic haemorrhagic disease virus (EHDV), and their host species. †: Assumed to be the same as BTV.
Table 3. The parameters for the transmission probability from vector to host ( ρ h v ) and host infectious period for bluetongue virus (BTV), African horse sickness virus (AHSV), Schmallenberg virus (SBV), epizootic haemorrhagic disease virus (EHDV), and their host species. †: Assumed to be the same as BTV.
VirusHost Species ρ v h Infectious PeriodRef.
BTVCattle0.920.6 [32,33]
AHSVHorse0.774.4 [34]
SBVCattle0.763.04 [29]
EHDVDeer0.9 †11.7 [35,36]
Table 4. Parameters considered in the uncertainty analysis, their fixed value while varying other parameters, and the ranges for which they are varied.
Table 4. Parameters considered in the uncertainty analysis, their fixed value while varying other parameters, and the ranges for which they are varied.
ParameterFixed ValueRange
ρ h v 0.500–1
ρ v h 0.750–1
T m i n 1510–20
α 0.050.01–0.10
χ 0.500–1
Table 5. Modes of action of stable netting scenarios simulated. Full descriptions of parameters are given in Table 1.
Table 5. Modes of action of stable netting scenarios simulated. Full descriptions of parameters are given in Table 1.
ScenarioBite Rate Reduction ( π )Relative Preprandial Mortality Rate Increase ( κ M )Relative Disarming Rate Increase ( κ D )Postprandial Mortality Increase ( ξ )
Untreated0.750.000.000.00
20% mortality0.750.200.000.20
40% mortality0.750.400.000.40
20% mortality & 20% disarming0.750.200.200.20
Table 6. The mean number of days per year for 2013–2022 with temperature above the minimum temperature threshold for virus replication for bluetongue virus (BTV), Schmallenberg virus (SBV), African horse sickness virus (AHSV), and epizootic haemorrhagic disease virus (EHDV), and the mean and maximum daily temperatures across Central England.
Table 6. The mean number of days per year for 2013–2022 with temperature above the minimum temperature threshold for virus replication for bluetongue virus (BTV), Schmallenberg virus (SBV), African horse sickness virus (AHSV), and epizootic haemorrhagic disease virus (EHDV), and the mean and maximum daily temperatures across Central England.
VirusThreshold TemperatureDays above Threshold Temperature
Mean TemperatureMaximum Temperature
BTV12.6136.4213.7
SBV12.4140.4217.5
AHSV12.6136.4213.7
EHDV19.511.479.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fairbanks, E.L.; Daly, J.M.; Tildesley, M.J. Modelling the Influence of Climate and Vector Control Interventions on Arbovirus Transmission. Viruses 2024, 16, 1221. https://doi.org/10.3390/v16081221

AMA Style

Fairbanks EL, Daly JM, Tildesley MJ. Modelling the Influence of Climate and Vector Control Interventions on Arbovirus Transmission. Viruses. 2024; 16(8):1221. https://doi.org/10.3390/v16081221

Chicago/Turabian Style

Fairbanks, Emma L., Janet M. Daly, and Michael J. Tildesley. 2024. "Modelling the Influence of Climate and Vector Control Interventions on Arbovirus Transmission" Viruses 16, no. 8: 1221. https://doi.org/10.3390/v16081221

APA Style

Fairbanks, E. L., Daly, J. M., & Tildesley, M. J. (2024). Modelling the Influence of Climate and Vector Control Interventions on Arbovirus Transmission. Viruses, 16(8), 1221. https://doi.org/10.3390/v16081221

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop