Once bitten, twice shy: A modeling framework for incorporating heterogeneous mosquito biting into transmission models

The risk and intensity of mosquito-borne disease outbreaks are tightly linked to the frequency at which mosquitoes feed on blood, also known as the biting rate. However, standard models of mosquito-borne disease transmission inherently assume that mosquitoes bite only once per reproductive cycle-an assumption commonly violated in nature. Drivers of multiple biting, such as host defensive behaviors or climate factors, also affect the mosquito gonotrophic cycle duration (GCD), a quantity customarily used to estimate the biting rate. Here, we present a novel framework for incorporating more complex mosquito biting behaviors into transmission models. This framework can account for heterogeneity in and linkages between mosquito biting rates and multiple biting. We provide general formulas for the basic offspring number, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>N</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> , and basic reproduction number, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>R</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> , threshold measures for mosquito population and pathogen transmission persistence, respectively. To exhibit its flexibility, we expand on specific models derived from the framework that arise from empirical, phenomenological, or mechanistic modeling perspectives. Using the gonotrophic cycle duration as a standard quantity to make comparisons among the models, we show that assumptions about the biting process strongly affect the relationship between GCD and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>R</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> . While under the standard assumption of one bite per reproductive cycle, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>R</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> is an increasing linear function of the inverse of the GCD, alternative models of the biting process can exhibit saturating or concave relationships. Critically, from a mechanistic perspective, decreases in the GCD can lead to substantial decreases in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>R</mml:mi> <mml:mn>0</mml:mn></mml:msub> </mml:math> . Through sensitivity analysis of the mechanistic model, we determine that parameters related to probing and ingesting success are the most important targets for disease control. This work highlights the importance of incorporating the behavioral dynamics of mosquitoes into transmission models and provides a method for evaluating how individual-level interventions against mosquito biting scale up to determine population-level mosquito-borne disease risk.