How do biting disease vectors behaviourally respond to host availability?

BACKGROUND: Ecological theory predicts a diverse range of functional responses of species to resource availability; but in the context of human blood consumption by disease vectors, a simplistic, linear response is ubiquitously assumed. A simple and flexible model formulation is presented that extends the Holling's Types to account for a wider range of qualitatively distinct behaviours, and used to examine the impact of different vector responses to the relative availability of multiple blood-host species.
RESULTS: Epidemiological models of falciparum malaria, Chagas disease and Lyme disease demonstrate that the standard, often implicit, assumption of a linear functional response can lead to spurious under- or over-estimates in disease transmission potential, across a full range of pathogen life-cycles. It is shown how the functional response in vector biting can augment disease intervention outcomes. Interactions between vector biting behaviour and uneven pathogen transmission probabilities between alternative hosts, as is the case for Chagas disease, can render infection more resilient to control.
CONCLUSIONS: Both the novel response formula and the nested vector-borne disease structure offer a flexible framework that can be applied to other vector-borne diseases in assessing the role of this newly identified aspect of biting behavioural ecology.

Background

How species respond to availability in resources is highly variable and has fostered considerable interest among ecologists for decades. Seminal papers written by Holling describe the different forms of functional response that predators exhibit to prey density [1-3]. A term first coined by Solomon [4], ‘functional response’ refers to the influence of resource availability on the rate of its consumption. Three qualitatively distinct functional responses were originally described: Type I responses depict resource consumption as a linear function of availability; Type II responses depict resource consumption as a decelerating function of availability (convex-up); and Type III responses depict resource consumption as an initially accelerating but then decelerating function of availability (s-shaped curve). Although Holling initially drew strong delineation between the types of predator species and the functional responses that they exhibit (e.g. invertebrates are Type II whereas mammals are Type III), subsequent ecological studies have generalised this phenomenon beyond predator-prey interactions (to account for all manner of resource consumption) and recast various species across a continuous spectrum of Types [5].

Real [6] devised a general formula that enabled flexible characterisation of a Hollings Type I, II or III response. This constituted an important contribution because different response Types could be assessed simply through re-parameterisation of the same underlying model, and this type of nesting is an advantage in direct comparison between models in estimating best fit to data [7]. The formulation proposed by Real [6] is as follows:where, C is the resource consumption rate, N is the density (or availability) of the resource, w is the handling time of the resource (i.e. the amount of time taken between identifying the prey and consuming it), and α and β are shape parameters. In the case that w = 0 and β = 1, C = αN (a linear Type I response); in the case that w > 0 and β = 1, a formulation which is equivalent to the Michaelis-Menten equation of enzyme kinetics results and is characterised by a decelerating consumption that eventually saturates (Type II); in the case that w > 0 and β > 1, a sigmoidal relationship results (Type III).

However, the extent of analogous developments in the context of vector-borne diseases is very limited. Many disease vectors are obligatorily haematophagous, meaning blood is a critical resource for survival and/or reproduction [8]. An important distinction to make here is that host death would rarely be expected as a direct consequence of blood-feeding by a vector (an exception being neurotoxin-mediated tick paralysis), but instead may result from coincident pathogen transmission. Thus the bidirectional effects on both consumer and resource species population dynamics implicit to predator/prey (and parasitoid/host) systems do not necessarily hold here. A further necessary consideration for systems of haematophagous disease vectors is that, as a rule, blood can be (and is) sourced from multiple host species. Even the most discerning of vector species thrive on the blood of multiple hosts [9, 10]. For example, Anopheles gambiae (sensu stricto) famously shows extreme preference for human blood [11] but adapts to environments with low human availability by sourcing blood from alternative mammals [12] with reportedly little-to-no effect on its resulting fecundity [13]. This is in stark contrast to parasitoid systems which are typically highly specialised and to most predator-prey theoretical and empirical studies (although, see multi-species models of Abrams & Matsuda [14] and Rueffler et al. [15]). To avoid conflation, the biting response of disease vectors to host availability shall be referred to as ‘behavioural’ instead of ‘functional’.

Development in the understanding of behavioural responses in disease vector biting is timely as this field is anticipated to accelerate rapidly following recent advances in molecular approaches; novel biological fingerprinting methods allow for the inexpensive, rapid and sensitive identification of host species from vector blood meals [16]. In the context of human disease control, there has been a recently rekindled interest in targeting vectors that bite alternative (often domestic) host species [17-19]. However, the full potential of these advances in informing vector-borne disease epidemiology will only be realised with their appropriate interpretation through a more developed ecological theory. Of the few ecological and epidemiological studies that consider host preference among vectors, linearity between alternative host availability and vector response is typically assumed [17, 20, 21]. In other words, a doubling in the availability of a particular host species relative to all potential hosts doubles the proportion of bites taken on that species. Here, the behavioural response of haematophagous arthropod disease vectors to the availability of alternative hosts and the resulting consequences for vector-borne disease transmission are explored.

Methods

Modelling the behavioural response in biting disease vectors

The proportion of blood-meals taken from the host species of interest (here, humans) compared to alternative hosts is generally assumed to increase as a direct proportion of increasing relative human availability. A model was sought to relax this key assumption of a century’s worth of vector-borne disease (VBD) models in order to explore the epidemiological impact of non-linear vector biting behavioural responses to host availabilities. While the nested model of Real (1977) (Equation 1) offers a concise and flexible framework for exploring different qualitative responses, there were VBD-specific scenarios that could not be resolved using this existing framework. For example, in the case of a zoonotic VBD which has spilled over into a local human population, a vector with strong zoophilic speciation may only opportunistically bite humans when their preferred host becomes vanishingly rare - this is the entomological/ epidemiological situation reported in Louisiana, USA, where the local kissing bugs (Triatoma sanguisuga) were observed to only start biting humans and infecting them with Trypanosoma cruzi when the local armadillo population collapsed [22]. This behaviour is not described by a Type I, II or III response. A new, flexible formula was developed to account for a wider range of vector responses to host availability:where, p is the proportion of all blood meals that are derived from the species of interest, for humans, this metric has been termed the ‘human blood index’ (HBI) [20]; Q is the availability of the host species of interest relative to all potential hosts; α and β are parameters that shape the behavioural response. Figure 1 demonstrates the effects that α and β have on the functional response with this new formulation. In addition to a linear (Type I), convex-up (Type II) and sigmoidal (Type III) response, the new formula also allows for responses that are convex-down (hereon referred to as ‘Type IV’) and atypical sigmoidal (classic s-shape reflected in the y = x, ‘Type V’). The qualitatively distinct feeding behaviours which can be characterised as combinations of intrinsic (genetic) host preferences and vector phenotypic response to local conditions [23] are described in Table 1.

Fig. 1

The behavioural response in human blood index of a disease vector to varying levels of human host availability (relative to all potential blood sources). Distinct qualitative forms (denoted ‘I’ to ‘V’) are shaped by parameters α and β as described in Equation 2

Table 1

The qualitatively different behavioural responses (parameterisation and associated vector behaviours) described by the new formula

Response type	Ecological equivalent	Parametric conditions	Vector behaviour
Type I	Analogous to Holling’s Type I	α = 1β = 1	Indiscriminate; or vector biting that is consistent (proportionate) across relative availabilities of alternative hosts.
Type II	Analogous to Holling’s Type II	α < 1β ≥ 1	The HBI of an anthropophilic vector saturates whereby even when humans and non-humans have similar availability, almost all blood meals are secured from humans.
Type III	Analogous to Holling’s Type III	α ≥ 1β > 1	Similar to a Type II response, the HBI saturates, but at low levels of human availability vectors are uninclined to bite them. Corresponding with the analogous Holling’s Type, this could be associated with a learned behaviour with an increased rate of human encounters.
Type IV	Inversion of Holling’s Type II	α > 1β ≤ 1	A zoophilic vector is uninclined to bite humans until they constitute all but the only available blood source.
Type V	Inversion of Holling’s Type III	α ≤ 1β < 1	HBI saturates and becomes relatively invariant when humans and non-human hosts are at similar availability. This is analogous to ‘negative prey switching’ whereby the ‘predator’ consumes disproportionately less of the more available ‘prey’ [41]. Eventually, when non-humans become vanishingly rare, the HBI is forced to increase sharply to unity.

The impact of these qualitatively distinct behavioural responses are assessed for different classes of VBDs. Infectious disease agents can be categorised across a spectrum according to their transmission potential to humans relative to non-human species [24]; and a natural, human-centric stratification is to consider the transmission potential to humans of pathogens that are either strict-anthroponotic (where non-human species are incompetent reservoirs), generalist (where humans and non-human species are both competent reservoirs) or strict-zoonotic (where humans are incompetent reservoirs). To exemplify the epidemiological impact of vector biting responses, these three strata are respectively represented by models of falciparum malaria, Chagas disease and Lyme disease.

Nested ecological-epidemiological models

The transmission dynamics of falciparum malaria, Chagas disease and Lyme disease are all nested within the following general VBD framework:

The epidemiological categories of the vector population and the populations of different host species (subscript H refers to humans; N refers to non-human hosts) are tracked as proportions. Disease transmission is assumed frequency-dependent to maintain convention with almost all VBD models.

Susceptible hosts (S) become infected (I) following a bite from an infectious vector (Z). Infected hosts can either revert to susceptible at rate γ, or they can benefit from temporary (τ > 0 and/or θ > 0) or permanent (τ = θ = 0) immunity. Alternatively, hosts can become asymptomatically infected (A) directly progressing from symptomatic infection (π > 0) or following on from recovery and subsequent reinfection (θ > 0). Asymptomatic infection may have the same transmission potential to vectors as symptomatic infections (σ = 1) or different transmission potential to vectors (σ ≠ 1), and can either be lifelong (κ = 0) or the pathogen can be completely cleared and hosts recover at rate κ. Susceptible vectors (X) become infected (Y) following a bite from an infectious host, and after the extrinsic incubation period, become infectious (Z). Total vectors V = X+Y+Z. Vectors typically outnumber hosts and, following convention, the ratio of vector-to-hosts (all blood-source host species) is denoted m. Typically, they also live shorter lives than hosts and so the relatively long extrinsic incubation period is explicitly included, as are vector demographics. Here, a stable vector population is assumed whereby births are set to balance deaths (μ). The impact of biting behaviour is assessed for pathogens with markedly different aetiologies through the following specifications and parameterisations.

Model specification

Humans are considered the only intermediate host for Plasmodium falciparum (c.f. knowlesi malaria for which mixed-host species models now exist [25]). Therefore, transmission terms between vectors and non-human hosts are assumed to equal zero. Following convention of previously published malaria models, host recovery without imparting some level of immunity does not occur; nor does asymptomatic chronicity following initial infection and so γ and π equal zero respectively [26]. The resulting compartmental framework is equivalent to an SIRS model but with possibility of maintained asymptomatic infection status following temporally proximal sequential infections [27].

Trypanosoma cruzi is a more generalist pathogen, infecting marsupials, primates, bats, armadillos and rodents, among other species [28]. In highly endemic human communities of Latin America, domestic animals are the key infection source and dogs are the primary parasite reservoir [29]. Infection dynamics are therefore tracked between dogs, the kissing bug vector and humans. (However, Equations 3–13 can be extended to greater numbers of blood-hosts and this will constitute important future work). Similarly for dogs and humans, infection is not cleared and acute infection invariably leads to chronic, asymptomatic infection. Hence, ε, γ, τ, θ and κ equal zero and the transitions are described by an SIA model. A key difference between these two blood-hosts is that while asymptomatically infected dogs can continue to transmit the parasite to vectors (σN > 0), chronically infected humans do not constitute parasite reservoirs [30].

Humans are dead-end hosts of Lyme disease. Therefore, human host transmission to the vector (b) equals zero. No temporary or lasting immunity has been documented for humans and so dynamics are described by an SIS model; whereas infection is generally chronic and symptomless in amplification host species (e.g. white-footed field mice, deer) and described by an SIA model [31]. These substructures can be achieved as with the other disease examples by setting the redundant rates to zero. Rates of change between the remaining epidemiological categories for all infection models are described in full in Table 2.

Table 2

Parameterisation for vector-borne disease models

	Definition	Plasmodium falciparum	Trypanosoma cruzi	Borrelia burgdorferi
b i	Transmission coefficient (vectors→hosts) = bite rate x transmission probability	0.1 = 1/3 × 0.3 (humans) [42]; 0 (non-humans)	2 × 10-5 = ¼ × 8 × 10-5 (humans) [43, 44]; 2.5 × 10-4 = ¼ × 0.001 (non-humans) [45]	0.003=1365×1.0humansandnon‐humans [46]
b Vi	Transmission coefficient (hosts→vectors) = bite rate x transmission probability	0.007=13×0.02humans [47]; 0(non-humans)	0.015 = ½ × 0.03 (humans); 0.25 = ½ × 0.49 (non-humans) [48]	0(humans); 0.003=1365×1.0non‐humans [46]
γ	Recovery rate (no immunity)	0 (humans and non-humans)	0 (humans and non-humans)	1/28 (humans)a; 0 (non-humans) [31]
ε	Clearance rate of symptomatic infection	1/200 (humans) [49]; 0 (non-humans)	0 (humans and non-humans)	0 (humans and non-humans)
κ	Clearance rate of asymptomatic infection	1/200 (humans) [49]; 0 (non-humans)	0 (humans and non-humans)	0 (humans and non-humans)
π	Asymptomatic primary infection rate	0 (humans and non-humans)	1/40 (humans and non-humans) [50, 51]	0 (humans); 1/28 (non-humans) [31]
θ	Asymptomatic secondary infection rate	0.5 (assumed for humans); 0 (non-humans)	0 (humans and non-humans)	0 (humans and non-humans)
τ	Full susceptibility reversion rate	1/1000 (humans) [52]; 0 (non-humans)	0 (humans and non-humans)	0 (humans and non-humans)
μ	Birth (or maturation) and death rate of vectors (i.e. stable population)	1/10 [53]	1/365 [54]	1/365 [55]
σ	Adjustment factor for asymptomatic transmissibility to vector	0.25 (humans) [56]; 0 (non-humans)	≈0 humans [30];≈1 non-humans [57]b	0 (humans);≈1 (non-humans)
ζ	Rate of parasite development within vector	1/10 [58]	1/10 [59]	1/365 [60]c

aClassically, Lyme disease infection dynamics are of an SIS form whereby the pathogen is assumed to be cleared by the host’s immune system. However, Nadelman & Wormser [31] review several studies demonstrating that an SIA form is more appropriate for non-human hosts

bA longitudinal study of domestic dogs (a principal Chagas disease reservoir) demonstrated persistent infectiousness but it was unclear whether this was a result of repeat infections

cParasite development is assumed to correspond with the developmental delays between life stages of the tick (whereby the tick will take its blood-meal from a different host species)

Results

Control strategies focusing on reducing vector biting rates were simulated for models representing each of the three pathogen life-cycle strata. Specifically, the impact on control of different behavioural responses in vector biting behaviour was assessed separately for falciparum malaria, Chagas disease and Lyme disease (Fig. 2). Due to their long duration, rates of chronic infection with vector-borne diseases are slow to change following control and so results are presented in terms of the impact of vector behavioural response on rates of acute (symptomatic) infections which better characterise disease incidence.

Fig. 2

Disease control efficacy is contingent on the behavioural response of biting vectors to the availability of alternative blood-hosts. The parameters α and β determine the shape of the behavioural response as described in Equation 2. The human proportion of all blood-hosts is indicated in the top-right of each plot. For Plasmodium falciparum (top row), the region above the contours corresponds with controlled transmission, but for Trypanosoma cruzi (middle row) and Borrelia burgdorferi (bottom row), the regions below the contours correspond with controlled transmission. A special case is shown in the left plot for T. cruzi whereby the high α/ low β region (above broken line) delimits a second parameter space for controlled transmission (see text). The contour labels correspond with the percentage reduction in bite rate required to achieve control. (These models are all deterministic and so a 90 % reduction in the acute infections relative to the maximum level in the absence of control is used to infer controlled transmission)

Intuitively, acute malaria infection was easier to control when vectors had strong zoophilic preference across a wide range of host availabilities (high α, low β, i.e. Type IV) and most resilient to control when vectors were highly anthropophilic across a wide range of host availabilities (low α, high β i.e. Type II). Also expected for this parasite that requires passage through humans for its propagation, was the result that when humans constituted a larger proportion of the total blood-host availability, infection was more difficult to control (from left to right on top row of Fig. 2).

Acute human infection with both Lyme and Chagas disease had qualitatively similar relationships with the vectors’ behavioural response. Control was easiest when vectors were highly anthropophilic (low α, high β, i.e. Type II), and hardest when vectors were zoophilic (high α, low β i.e. Type IV). While this was an expected result for Lyme disease (whereby humans are dead-end hosts and thereby detract from the pathogen transmission cycle), the qualitatively similar result for Chagas was less intuitive and arose from the disproportionately high contribution of non-human hosts to the force of infection (dogs not only remain infectious for much longer than humans but are substantially more infectious to the vectors, Table 2). Host preference that was weighted against humans (high α, low β, i.e. Type IV) was most resilient to human infection control for all Lyme and Chagas disease scenarios, with one exception. For a generalist pathogen (T. cruzi), there are some circumstances when a more intermediate biting behaviour (less definitively anthropophilic/zoophilic) can give rise to a human infection that is less amenable to control (left subplot, middle row of Fig. 2). Here, a Type II response (low α, high β) reduces the proportion of bites on the more competent host (dogs), in turn reducing the prevalence of disease and subsequent force of infection on humans. However, when humans only constitute a minority of available blood-hosts, the exaggerated proportion of bites on dogs under a Type IV response (high α, low β) depletes the opportunities for pathogen spread to humans more than the reduced prevalence of disease that would result for a less zoophilic vector. Therefore, both extremes of vector biting behaviour are more conducive to control.

Discussion

Models of diseases spread by haematophagous arthropods are increasingly popular tools for understanding the spread of vector-borne diseases and strategizing their control. Persisting among even the most complex of contemporary models is the widespread assumption that vector bites are distributed in a directly proportionate manner on alternative hosts according to their relative availability, sometimes adjusted according to a constant intrinsic host preference as described by Bailey [32]. Critically, a linear (Type I) functional response is not only atypical for arthropods, it is without precedent [33].

Unfortunately, examples of studies that account for non-linear effects of human host availability on vector biting behaviour are scant. While parallel advances have been made in the context of predator-prey [34] and host-parasitoid systems [35], these have not translated to corresponding developments in VBD understanding. Antonovics et al. [36] used a general Type II model to simulate VBD transmission and described how a vector that is highly restricted in its movement could well have its bite rate limited by low host density (invalidating the general assumption of frequency-dependence); and, frequency-dependence and density-dependence are well known to generate markedly different transmission dynamics [37]. This concept has been recently built upon by Kershenbaum et al. [38] who showed in a one-vector two-host (one competent, one incompetent) model that parameter spaces exist whereby a reduction in competent host availability risks exacerbating disease prevalence through attenuated competition for limited vector feeding sites. The current study does not account for this potential transmission bottleneck; and a unified approach to account for not only the proportional distribution of bites across different alternative host availabilities but also a vector biting rate that can be influenced by host densities constitutes an important future endeavour.

A further aspect of the current work that justifies future development is the exploration of how different behavioural responses might impact pathogen transmission dilution and amplification [39]. Miller & Huppert [40] also used a Type II formulation to describe vector bites split between multiple host species and showed disease transmission is intensified (or ‘amplified’) when more host species are included in a system if the vector prefers the host with the highest transmission ability; otherwise, the addition of more host species dilutes transmission [40]. Although the current study is restricted to a vector shared between only two alternative hosts, the methods described are easily adaptable to more hosts for the analysis of species diversity and pathogen persistence.

Conclusions

This study makes the following contributions: it provides a new, two-parameter function that can be used to understand functional response ecology and that extends the qualitative Types achievable with earlier (three-parameter) models of Holling and Real; it introduces a general framework into which this vector behaviour can be incorporated to explore its consequences on infectious diseases with diverse epidemiology; it highlights an aspect of vector behaviour that is almost completely neglected; and demonstrates how this response of the vector to alternative host availabilities can drastically alter efforts to mitigate transmission. Data derived from laboratory or semi-field conditions (where the relative availabilities of alternative blood hosts can more easily be manipulated) are a high research priority in the field of vector-borne disease. Diverse aetiologies and idiosyncratic epidemiology severely limit opportunities for scientific discoveries that can potentially impact the whole gamut of vector-borne diseases. The ecology of this behaviour in biting disease vectors potentially offers one of the last largely unexplored avenues of generally applicable vector-borne disease research.