Mosquito and human characteristics influence natural Anopheline biting behavior and Plasmodium falciparum transmission ====================================================================================================================== * Christine F Markwalter * Zena Lapp * Lucy Abel * Emmah Kimachas * Evans Omollo * Elizabeth Freedman * Tabitha Chepkwony * Mark Amunga * Tyler McCormick * Sophie Bérubé * Judith N Mangeni * Amy Wesolowski * Andrew A Obala * Steve M Taylor * Wendy P O’Meara ## Abstract The human infectious reservoir of *Plasmodium falciparum* malaria parasites is governed by the efficiency of parasite transmission during vector human contact as well as mosquito biting preferences. Understanding mosquito biting bias in a natural setting can help inform precise targeting of interventions to efficiently interrupt transmission. In a 15-month longitudinal cohort study in a high transmission setting in western Kenya, we investigated human and mosquito factors associated with differential mosquito biting by matching human DNA in single- and multi-source *Anopheles* bloodmeals to the individuals they bit. We employed risk factor analyses and econometric models of probabilistic choice to assess mosquito biting behavior with respect to both human-to-mosquito transmission and mosquito-to-human transmission. We observed that *P. falciparum*-infected school-age boys accounted for 50% of bites potentially leading to onward transmission to mosquitoes and had an entomological inoculation rate 6.4x higher than any other group, that infectious mosquitoes were 2.8x more likely to bite cohort members harboring *P. falciparum* parasites compared to noninfectious mosquitoes, and that this preference to feed on infected people was enhanced by the presence of higher sporozoite loads in the mosquito head-thorax. Taken together, these results suggest that school-age boys disproportionately contribute to the *P. falciparum* transmission cycle and that *P. falciparum* sporozoites modify mosquito biting preferences to favor feeding on infected people. **Significance** The malaria parasite transmission cycle is doubly-dependent on mosquito-human contact rates. To efficiently deploy transmission-reducing interventions, it is important to understand how mosquito biting preferences shape the human infectious reservoir. Here, we match human DNA in mosquito bloodmeals to the people they bit to quantify mosquito biting preferences and understand how mosquito characteristics shape these preferences. We observed that school-age boys were bitten the most and contributed most to onward transmission to mosquitoes. We also observed that, compared to non-infectious mosquitoes, mosquitoes harboring infectious *Plasmodium falciparum* were more likely to bite *P. falciparum*-infected people. These observations increase our understanding of malaria parasite transmission and evolution and provide a foundation for developing effective transmission-reducing interventions. ## Introduction *Plasmodium falciparum* parasites are transmitted between human host and mosquito vector to maintain a complex and robust transmission cycle. The number of human *P. falciparum* infections arising from a single infected person (*R0*) is governed by the rate of vector-host contact and by the efficiency of parasite transmission during contact (1). Therefore, both nonrandom vector-host contact and transmission efficiency can drastically alter transmission patterns and influence intervention effectiveness (2). To accurately estimate transmission patterns and efficiently deploy transmission-reducing interventions, it is important to understand who is being bitten by *P. falciparum* vectors, who is transmitting to mosquitoes, and what mosquito characteristics shape biting behavior. Contacts between humans and mosquitoes are driven by mosquito vector abundance and biting preferences (3). Laboratory and semi-field studies have demonstrated high inter-individual variability in attractiveness to mosquitoes between people (4, 5) and that *Plasmodium* infection in the mosquito influences mosquito biting preferences and behaviors (6, 7). Studies of human and mosquito samples collected in natural settings have also demonstrated strong heterogeneity in mosquito bite exposure (8), as well as in which people carry gametocytes and are likely to transmit to mosquitoes (9–12). While these works have provided valuable insight into the human *P. falciparum* infectious reservoir, they are experimental (10–16), cross-sectional (8, 17–19), short-term (20), or semi-longitudinal studies (21). Few studies quantify variation in mosquito-human contact in a natural setting across a well-defined host population over an entire transmission season, and, to our knowledge, none have incorporated multisource bloodmeals, which can lead to increased biting rates and amplified transmission (22). Here, we explored mosquito biting behavior and human-to-mosquito transmission in a community-based cohort in Western Kenya using weekly collections of indoor-resting mosquitoes from 75 households over 15 months of observation. To precisely understand who was bitten by an individual mosquito, we matched human short tandem repeat (STR) genotypes between single- and multi-source mosquito bloodmeals and cohort members. We hypothesized that mosquitoes bite certain groups of people more than others, and that mosquito characteristics, including species and *P. falciparum* infection status, influence their biting behavior. ## Methods ### Ethical statement We obtained written informed consent from all heads of households and all participants or their parent for those <18 years; the latter also provided verbal assent if >8 years. The study was approved by the ethical review boards of Moi (2017/36) and Duke (Pro00082000) Universities. ### Study population and data collection We analyzed samples and data collected from July 2020 to September 2021 from a previously described ongoing cohort (23). Briefly, in Bungoma County, Kenya, all residents above 1 year of age were included from 75 households, which were selected by radial sampling in 5 villages of varying malaria intensity. We collected basic demographic and household data yearly. We collected dried blood spots (DBS) from fingerpricks from all participants every month and recorded bednet use (active collection). We captured interval episodes of malaria by providing testing by malaria rapid diagnostic test (RDT) upon request owing to self-reported symptoms; those with positive RDTs were referred for treatment with Artemether-Lumefantrine. A DBS was also collected at these visits (passive collection). One morning each week for three weeks each month, we collected indoor-resting mosquitoes from each household using vacuum aspiration with Prokopaks. Mosquitoes were collected twice per month for bloodmeal analysis and once per month for rearing. Mosquitoes collected for rearing were collected on the same day as monthly DBS collections. ### Entomological methods Collected mosquitoes either underwent bloodmeal analysis or were reared. For bloodmeal analysis, female *Anopheles spp.* were immediately processed after collection, identified morphologically, and their abdominal status graded. Each was transected to separate the head-thorax, wings, and abdomen. Abdomens graded as bloodfed were pressed onto filter paper, and other parts were stored in tubes. Mosquitoes collected for rearing were released into household-specific cages inside an insectary maintained at 27±2⁰C and 80±10% relative humidity. Mosquitoes fed on cotton saturated with 10% sucrose in distilled water for 7 days prior to sacrificing with chloroform. All female *Anopheles spp.* were morphologically identified, transected, and stored as described above. ### Molecular methods #### Detection of P. falciparum by qPCR gDNA was extracted using Chelex-100 from all i) DBS, ii) head-thoraces of mosquitoes processed for bloodmeal analysis, and iii) abdomens of reared mosquitoes. gDNA from each DBS and mosquito abdomen was tested in duplicate in a TaqMan real-time PCR assay targeting the *P. falciparum pfr364* motif as previously described (24). Samples were defined as *P. falciparum* positive if: i) both replicates amplified *P. falciparum* with Ct values <40 or ii) 1 replicate amplified *P. falciparum* with Ct value <38. Parasite densities were estimated on each reaction plate using standard curves, which consisted of extracts from *P. falciparum* 3D7 parasites across a range of parasites/uL (for DBS) or of *P. falciparum* 3D7 gDNA across a range of ng/uL (for sporozoites). The latter was converted to parasite density using a genome size of 22.8 Mbps. #### P. falciparum csp amplicon deep sequencing All *P. falciparum*-positive DBS and reared mosquito abdomens were genotyped across a variable segment of the parasite circumsporozoite gene (*csp*) as previously described (25) and detailed in the Supplementary Methods. #### Resolution of Anopheles gambiae sibling species by PCR gDNA was extracted by a hotshot technique (26) from a single wing of female *Anopheles* mosquitoes that were morphologically identified as *An. gambiae s.l.* Sibling species were resolved using a multiplex assay (27) based on the intergenic spacer region (IGS) which gives bands of differing sizes for each of 5 species in the *An. gambaie* complex. #### Human genotyping by short tandem repeats We genotyped human blood in both a DBS from each participant and in each mosquito abdomen classified as bloodfed for which a blood spot was available (1064/1563). We extracted gDNA from each sample type using Chelex-100, and genotyped using the Promega Geneprint 10 assay as described in Lapp et al (28). ### Matching mosquitoes and people #### Matching bloodmeals to human hosts Mosquito bloodmeals were matched to cohort participants using the bistro R package (v0.2.0; (28)), which can provide matches for incomplete STR profiles and bloodmeals derived from multiple human sources. Briefly, matches based on weight-of-evidence likelihood ratios were determined for each mosquito-human pair and individual mosquito-based thresholds. The estimated number of contributors (NOC) to a bloodmeal is computed as ⌈max(𝑎) /2⌉, where *a* is the number of alleles at a locus. We define single-source bloodmeals as those with NOC=1 and multisource bloodmeals as those with NOC>1. #### Matching infected mosquito abdomens to infected humans We matched *P. falciparum* infections in reared Anopheles abdomens to people that might have infected them by first subsetting to *P. falciparum*-infected individuals present in the family compound on the day the mosquito was collected. Among these mosquito-human pairs, we define a match as any pair that shared at least one *csp* haplotype. ### Estimating the maximum probability of a mosquito becoming infected after biting an infected individual If we define 𝐵*kI* as the mosquito biting *k* infected people, then our goal is to compute 𝑃(𝑀*I*|𝐵1*I*). Briefly, we use the equation: ![Formula][1] Where 𝑃(𝑀*I*) is the probability of a mosquito becoming infected. We take 𝑃(𝑀*I*|𝐵*I*) to be zero, assuming all infections in the mosquito abdomen had to arise from the most recent biting event, and we assume that no mosquito bites more than two infected people in a given night, giving us: ![Formula][2] Given that the probability of a mosquito becoming infected after biting to infected people is greater than or equal to that after biting one infected person (i.e. 𝑃(𝑀*I*|𝐵2*I*) ≥ 𝑃(𝑀*I*|𝐵1*I*)), we assume 𝑃(𝑀*I*|𝐵2*I*) = 𝑃(𝑀*I*|𝐵1*I*) and solve for a maximum transmission efficiency 𝑃(𝑀*I*|𝐵1*I*). 𝑃(𝑀*I*) was considered to be the proportion of abdomens from reared mosquitoes that developed *P. falciparum* infections divided by the proportion of immediately processed mosquitoes that recently fed (i.e. freshly fed, half-gravid, or gravid). The probabilities that a mosquito bit one or two infected people were computed from the STR matching data using only mosquitoes with complete STR profiles. 95% bootstrapped confidence intervals (CIs) were computed. ### Calculating biting bias We used the gini function in the R package DescTools v0.99.48 (29) to determine the Gini index and 95% CI for mosquito biting bias based on each person’s mean nightly bites, calculated as the cumulative number of bloodmeals matched to an individual divided by the total person-nights in the study. We define random biting as each person having the same probability of receiving a mosquito bite per night at risk. To evaluate the expected Gini index under a random distribution of bites, we computed the Gini index for 1000 simulated random bite distributions where the probability of a person being bitten was weighted by the amount of time they were present during the study period and identified the 2.5%, 50%, and 97.5% percentiles. ### Multilevel model analysis We performed a risk factor analysis to assess the characteristics of people who were more likely to be bitten by mosquitoes. For each household-night with at least one STR-matched mosquito, we included in the sample anyone who slept in the household within 30 days of that night; people were considered to have received *n* bites, where *n* is the number of STR-typed mosquitoes that they matched to on that night. We used the glmmTMB function from the R package glmmTMB v1.1.6 (30) to run a negative binomial model (family = nbinom2) with a random effect on the person. The outcome was the number of bites a person received on a given night. Risk factors included gender, categorized age (<5, 5-15, and >15 years), whether the person slept under a net, and whether the person was infected with *P. falciparum.* Detailed description of risk factors, adjustors, sensitivity analyses and subsetted models are in the Supplemental Methods. ### Choice model analysis We computed a weighted discrete choice model for each human risk factor for being bitten using the R package mlogit v1.1.1 (31). We subset to include only *An. gambiae SS* and *An. funestus* mosquitoes with STR and sporozoite detection results available. The following three mosquito characteristics were included in each model: mosquito species, single- vs multisource meal, and *P. falciparum* sporozoite presence in the mosquito head-thorax (i.e. mosquito infectiousness). Each mosquito was identified as an individual, and all people present in the household when that mosquito was collected were included in the model. Values of the human characteristics were the alternatives, weighted by the number of people in the household in each group. An example of the data structure and analysis can be found in the supplementary information. Mosquitoes that matched to two people, and the people present in the household at the time, were included twice in the model, with each bitten person coded as being bitten in one of the entries. We also computed a choice model with mosquito sporozoite load in the head-thorax as a continuous variable; uninfected mosquitoes were assigned a sporozoite density of 0. Finally, we assessed whether the human characteristics associated with taking multisource meals were driven by shared sleeping spaces (See Supplementary Methods). ### Data analysis and visualization All analyses and visualizations were performed in RStudio v2023.06.1.524 (32) with R v4.2.1 (33) using the following packages: tidyverse v2.0.0 (34), haven v2.5.3 (35), Hmisc v5.0.1 (36), ggpubr v0.6.0 (37), geomtextpath v0.1.1 (38), DiagrammeR v1.0.9 (39), broom v1.0.5 (40), broom.mixed v0.2.9.4 (41), ggeffects v1.3.2 (42), gtsummary v1.7.0 (43), ggmosaic v0.3.3 (44), biostat3 v0.1.8 (45), modelr v0.1.11 (46). Code is available at [https://github.com/duke-malaria-collaboratory/imbibe_manuscript](https://github.com/duke-malaria-collaboratory/imbibe_manuscript). Data supporting these analyses are available on request. ## Results From July 2020 to September 2021, we collected 3660 female Anopheles mosquitoes across 75 households in 5 villages(**Figure S1**). Most (3038) were immediately processed to determine abdominal status and preserve the blood from freshly fed mosquitoes (1563) to match to cohort members. The other 622 were reared to enable identification of successful abdominal infections. We also collected 6414 human DBS through active and passive case detection, of which 1758 (27%) were positive for *P. falciparum* (**Figure S2**). ### Female Anophelines rest where they fed Among the 1563 freshly-fed female Anophelines, 1064 were available for STR typing (**Figure S3**). Of these, 777 (73%) returned STR profiles (**Figures 1A**, **S3**), among which 654 (84%) were single-source bloodmeals, and 123 (16%) had more than one contributor (i.e. were multi-source). We identified at least one cohort participant in 662 (85%) bloodmeals (n=550 single-source, n=112 multisource), and an additional cohort participant in 58/112 (52%) matched multi-source bloodmeals (**Figure 1B**). Among the 720 identified mosquito-human pairs, the mosquito was collected from the family compound where the person lived in 94% (n=676) of cases (**Figure 1C**). Additionally, for the 58 multisource bloodmeals where two cohort members were matched, the two people were from the same family compound in 88% (n=51) of cases (**Figure 1D**). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/01/25/2024.01.24.24301433/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/01/25/2024.01.24.24301433/F1) Figure 1: Mosquitoes usually matched to human hosts from the household where they were collected. (A) Number of immediately processed mosquitoes (y) by collection date (x), colored by whether they were STR-typed and whether they matched. (B) Matching by number of contributors (NOC, x) to the mosquito bloodmeal. Fill indicates the number of different humans matched to a mosquito bloodmeal. (C) Circle plot indicating whether the mosquito was collected in the same village, family compound, or building as the person it bit. (D) Circle plot indicating whether people that were bitten by the same mosquito came from the same village, family compound, building, or sleeping space. ### Some *P. falciparum* haplotypes in reared mosquito abdomens spatiotemporally match to those from infected humans Among the 622 female *Anopheles* that were reared following collection, 412 survived to day 6 or later, of which 34 (8.3%) harbored *P. falciparum* in their abdomen. We estimated that the maximum transmission efficiency, defined as the probability of a mosquito successfully becoming infected after biting an infected individual, was 0.16 (95% CI: 0.14-0.27). Of the infected mosquitoes, 30 (88%) infections were successfully sequenced at the *Pf-csp* locus (**Figure S4A**). We observed 28 distinct haplotype sequences in mosquitoes, 10 (36%) of which were found in at least one *An. gambiae SS* and one *An. funestus* abdomen (**Figure S4B**). Additionally, we successfully sequenced 37/57 (65%) *P. falciparum* positive human DBS from household-days that yielded a reared mosquito with a sequenced infection (**Figure S4C**). For 13/30 (43%) infected mosquitoes, we identified at least one haplotype match to an infected household member on the household-day of mosquito collection (**Figure S4D**). Given this small sample size, we used the STR profile matches for all subsequent analyses. ### Mosquito biting is non-random We next investigated the distribution of mosquito bites across the cohort and which cohort members had the highest rate of biting. For 382/588 (65%) participants, we did not observe any bites. Among the 206 participants who were bitten, the maximum number of bites observed for an individual in a single night was 12 (median = 1). Across all collections, mean nightly biting rates ranged from 0.03-1.5 bites per person-night at risk. Individual biting rates were more unequal (Gini index: 0.82, 95% CI: 0.80-0.85) than expected by random chance (Gini index: 0.51, 95% CI: 0.48-0.54), with 20% of people receiving 86% of bites per night at risk (**Figure 2A-B**). The highest biting rates were observed for males 5-15 years old who were infected with *P. falciparum* and didn’t sleep under a bed net, with an estimated 12 bites per month per person (**Figure 2C**). These trends held in a multivariate model (**Figure 2D**, **Table 1, Table S1**): we observed higher biting rates on males (biting rate ratio (BRR): 1.68; CI: 1.28-2.19), children 5-15 (BRR: 1.49; CI: 1.13-1.98), and individuals infected with *P. falciparum* (BRR: 1.25; CI: 1.01-1.55), while lower biting rates were observed on net users (BRR: 0.51; CI: 0.40-0.65). All of these trends held in sensitivity analyses, with some wider confidence intervals due to smaller sample sizes (**Table S2**). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/01/25/2024.01.24.24301433/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/01/25/2024.01.24.24301433/F2) Figure 2: Biting is heterogeneous and nonrandom. (A) Individual observed bite counts (color) for each member (y) of each household (x) showing between- and within-household heterogeneity. Grey indicates individuals with no observed bites. Only people present for more than 3 months during the study period were included. (B) The observed distribution of bites is more unequal than random chance. Lorenz curves demonstrating cumulative share of nightly bites (y) across participants (x, sorted by ascending share of bites). Dashed line indicates the null of evenly distributed bites, dotted line the null of simulated (1000x) randomly-distributed bites, and solid line the observed distribution. Overall, 20% of people received 86% of observed bites. (C) Estimated monthly bites per person (from observed average nightly bites; color and number) by 4 risk factors (net use, *P. falciparum* infection status, gender, and age). (D) Biting rate ratios and 95% confidence intervals for covariates in risk factor analysis. Adjusters include high transmission season, number of STR-typed mosquitoes in the household, number of household members, whether any household member tested positive for *P. falciparum* by a rapid diagnostic test in the previous month, and number of people in the sleeping space. View this table: [Table 1:](http://medrxiv.org/content/early/2024/01/25/2024.01.24.24301433/T1) Table 1: Person-night characteristics of bitten and not bitten people To estimate biting rates on potentially infectious individuals and thereby better understand how the observed differential biting by mosquitoes could influence human contributions to onward transmission, we restricted our multivariate model to individuals infected with *P. falciparum*. This model yielded biting rates (**Figure S5A**) similar to those observed among all people (**Figure 2D**). From this model we estimated that, in our cohort, males 5-15 receive about 50% of bites that could lead to onward *P. falciparum* transmission to mosquitoes (**Figure 3A**). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/01/25/2024.01.24.24301433/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/01/25/2024.01.24.24301433/F3) Figure 3: Estimated group-level biting rates with respect to potential *P. falciparum* transmission. (A) Human-to-mosquito *P. falciparum* transmission using a model including only infected individuals. Groups separated by gender on the x-axis, categorical age on the y-axis (and color), and net use by shade. Size of box and number displayed in box show the proportion of bites received by individuals infected with *P. falciparum* in the cohort, representing potential contribution to onward transmission. (B) Mosquito-to-human transmission using a model including only infectious mosquitoes. Color and shade indicate the categories as in panel A. Next, to better understand biting heterogeneity on the other side of the transmission cycle, from mosquitoes to humans, we investigated biting rates on people by infectious mosquitoes only (**Figure S5B**). We observed similar estimates for the BRRs compared to those among all mosquitoes except for individuals infected with *P. falciparum*, who had a higher BRR among infectious mosquitoes (2.26; CI: 1.40-3.66) compared to all mosquitoes. Estimating relative entomological inoculation rates (EIRs) for each category of individuals revealed that, in our cohort, males 5-15 infected with *P. falciparum* have an EIR 6.4x higher than any other group of people (**Figure 3B**). ### Mosquito characteristics influence their choice of bloodmeal We investigated associations between biting behavior and three mosquito characteristics: species, harboring a multi-source bloodmeal, and the presence of sporozoites in the head-thorax. *P. falciparum*-infected *An. funestus* took more multisource bloodmeals compared to uninfected *An. funestus* (p = 0.002), whereas no difference was observed in the proportion of multisource meals taken by infected and uninfected *An. gambiae SS* (p = 0.439) (**Figure S6A**). However, among mosquitoes positive for *P. falciparum* sporozoites, no differences were observed in sporozoite densities in the head-thoraces between mosquitoes that took single- vs multi-source bloodmeals for both *An. funestus* (p = 0.15) and *An. gambiae SS* (p = 0.35) (**Figure S6B**). We next employed a discrete choice model framework to explore how mosquito characteristics impact the relative risk of a bite for humans that differ in each of the risk factors from the previous analysis, assuming a human population similar to the cohort analyzed here (**Figure 4A, Table S3**). We observed that host selection was in general similar between species; however, biting preferences varied in relation to number of bloodmeals and mosquito infectiousness. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/01/25/2024.01.24.24301433/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/01/25/2024.01.24.24301433/F4) Figure 4: Mosquito characteristics influence their choice of bloodmeal. (A) Relative risk ratios estimated in choice models on the preference for humans with various characteristics. A separate multivariate choice model was run for each human characteristic. Dots indicate point estimates, and bars the 95% confidence intervals. (B) Probability of biting an infected person as a function of sporozoite density. Ref: Reference, y: years. In contrast with the observations using all bloodmeals, mosquitoes that bit more than one person were less likely to choose to bite a male host (RRR 0.50, 95% CI 0.33-0.77) and more likely to have bitten a child under 5 years (RRR 4.10, 95% CI 1.80-9.35) compared to an individual > 15 years. This could be driven, in part, because pairs of people who share a sleeping space are more likely to be found in a multisource bloodmeal compared to people who do not share a sleeping space (unadjusted RR 2.80, 95% CI 1.50-5.23). Females (vs. males) and children <5 years (vs individuals > 15 years) are more likely to share a sleeping space with someone (unadjusted RR females 1.09, 95% CI 1.00 - 1.20; unadjusted RR < 5 years 1.23, 95% CI 1.12-1.36). Amongst the investigated mosquito characteristics, biting behavior was most strongly associated with infection by *P. falciparum* sporozoites. Compared to uninfectious mosquitoes, infectious mosquitoes prefer males over females (RRR 2.02, 95% CI 1.28-3.18), children < 5 years (RRR 2.69, 95% CI 1.19-6.09) and 5-15 years (RRR 2.21, 95% CI 1.37-3.59) compared to >15 years, and individuals not sleeping under a net over those who did (RRR 2.22, 95% CI 1.15-4.26). Most strikingly, individuals infected with *P. falciparum* were nearly 3x more likely than uninfected individuals to be bitten by infectious mosquitoes (RRR 2.76, 95% CI 1.65-4.61). This association could result from inherent biases for mosquitoes to bite specific people; therefore, we tested for a dose-response relationship by repeating the discrete choice model using continuous sporozoite density in place of binary infectiousness. In this model, higher sporozoite densities were associated with a greater likelihood of biting people infected with *P. falciparum* than uninfected people (RRR per 100 sporozoites: 1.92, CI 1.23-2.98; **Figure 4B**, **Table S4**), regardless of species or number of human bloodmeal sources (**Figure S7**). This observation further supports the finding that the preference of Anopheline mosquitoes for biting people infected with *P. falciparum* is enhanced by the presence of sporozoites in the mosquito. ## Discussion In a 15-month longitudinal cohort study in a high-transmission setting in western Kenya, we investigated the human and mosquito factors associated with differential mosquito biting. To do so, we matched human DNA in single- and multi-source *Anopheles* bloodmeals to the individuals they bit and employed measures of inequality, risk factor analyses, and econometric models of probabilistic choice to assess mosquito biting behavior with respect to both human-to-mosquito transmission and mosquito-to-human transmission. We observed that school-age boys were bitten the most, that infectious mosquitoes were more likely to bite cohort members harboring *P. falciparum* parasites, and that this preference to feed on infected people was enhanced by the presence of higher sporozoite loads in the mosquito head-thorax. Taken together, these results suggest that school-age boys disproportionately contribute to onward *P. falciparum* transmission and that *P. falciparum* sporozoites modify mosquito biting preferences to favor feeding on infected people. Infectious mosquitoes were nearly 3x more likely to bite people harboring *P. falciparum* parasites (vs. uninfected people) compared to uninfected mosquitoes. This observation was consistent when expressed as a density of sporozoites, supporting a causal link between the presence of sporozoites and a preference to feed upon infected people. Note that, after a mosquito bites an infected person it takes about 7-10 days for the mosquito to become infectious; therefore, we can assign directionality to the relationship between mosquito infectiousness and human *P. falciparum* infection status here. To our knowledge, this is the first evidence that *P. falciparum* sporozoites may manipulate mosquito feeding to favor *P. falciparum*-infected humans. Several prior studies have provided experimental evidence that *Plasmodium spp* sporozoite infection modifies specific aspects of Anopheline feeding behavior, including enhancing probing time (47), feeding persistence (48), feeding success (49, 50), and attraction to human odors (6). Prior laboratory and semi-field studies have also reported that being infected with *Plasmodium spp* increases a mouse’s (51) or person’s attractiveness to mosquitoes (52–54), and our observations suggest that this increased attractiveness is further enhanced when the mosquito is likewise infected with *P. falciparum*. Such enhancement could result if sporozoite infection increases olfactory responsiveness (7, 55) to volatiles released during human *P. falciparum* infection. One question that remains unanswered is whether this heightened attractiveness is present among all *P. falciparum*-infected individuals, or differs by asexual parasite or gametocyte density. Regardless, how this phenomenon might influence overall transmission intensity is not clear: although modeling suggests that altered mosquito behavior due to *P. falciparum* infection could increase the force of infection by several-fold (56), it is possible that, given the high biting heterogeneity we observed, it may also concentrate the force of infection into a small group of perpetually infected people and therefore not increase population-level prevalence. Either way, our observations support the notion that, in this natural setting, *P. falciparum* parasites manipulate Anopheline vectors to favor feeding on humans who are already likewise infected with *P. falciparum*. We estimated biting rates for three distinct groups of mosquito-human pairs: (1) all mosquitoes and all cohort members, to assess overall mosquito biting preferences; (2) all mosquitoes and people harboring *P. falciparum* parasites, to assess what portion of the human population likely contributes most to transmission to mosquitoes; and (3) infectious mosquitoes and all participants, to assess the relative EIR for each human risk factor. The highest biting rates for all of these groups were on school-age boys. This demographic accounted for 50% of all bites potentially leading to parasite transmission to mosquitoes despite comprising only 19.2% of the population at risk. Previous studies have shown that this age group is also most infective to mosquitoes in feeding assays (10), and comprises the largest segment of the infectious reservoir, defined as infectiousness to mosquitoes projected over a population (10, 12, 57). Given that the rate of secondary mosquito infections is a product of an individual’s mosquito exposure and their infectiousness to mosquitoes (58), the disproportionate contribution of school-age boys to onward *P. falciparum* transmission reported here is likely underestimated. These observations, coupled with the understanding that school age children more often harbor asymptomatic or chronic infections (59–62), point to them playing a significant role in maintaining transmission and highlight the importance of targeting interventions toward this group as a means to efficiently reduce overall community transmission. Prior studies using STR typing in wild-caught bloodfed mosquitoes have also investigated features that make certain people more likely to be bitten by mosquitoes. These studies showed either no difference in biting rates by gender (17, 63) or higher biting rates for males compared to females (8, 18–20), consistent with our findings. They also observed that individuals aged >15 years tend to receive the most bites from malaria vectors, a pattern often attributed to higher body surface area (8, 18–20). However, we observed that school age boys have higher biting rates than other age groups, and others have observed that biting preferences by age varied between sites and seasons (64). One potential reason for this is that elements of human attractiveness beyond surface area, such as body odor (5), play a role in host choice. Furthermore, our observation that people sleeping in the same sleeping space were more likely to receive bites from mosquitoes that took multi-source meals also supports the notion that accessibility plays a role in host selection. Taken together, our observations highlight that mosquito biting choices likely result from a combination of host attractiveness, human behavior, and accessibility. Our study provides a unique opportunity to consider the potential evolutionary advantages resulting from the influence of mosquito behavior by parasites on parasite diversification. Genetic diversification arises from de novo mutation and from recombination during the sexual stages in the mosquito; the latter likely plays a major role in adaptation (63). For diversifying recombination to occur, the mosquito must be infected with multiple distinct *P. falciparum* strains during the span of a few hours (65). This can occur if the mosquito feeds on multiple infected individuals in a short timespan (mosquito super-infection), which occurred only 3.5% of the time in our dataset. Alternatively, polygenomic mosquito infections can also result when a mosquito feeds upon a person who has been co- or super-infected with multiple parasite strains, which is a common phenomenon in high-transmission settings like ours, where we have observed up to 17 strains in human infections (24). Our observation that infectious mosquitoes are more likely to bite infected individuals provides a plausible mechanism by which these genetically diverse infections accumulate within hosts, and thereby increases the possibility of advantageous recombination in a subsequent mosquito co-infection. This would allow parasites to explore a more diverse fitness landscape through recombination between different strains from different initial infections. Moreover, our observation that *An. gambiae SS* and *An. funestus* have relatively similar biting preferences and share haplotypes suggests that there is likely one recombining reservoir of *P. falciparum* strains shared amongst both of these species. Given the evidence that diverse blood-stage infections may be partially prevented by prevalent blood-stage *Plasmodium* infections (66), this phenomenon suggests a mechanism by which the parasite has adapted to generate opportunities for out-crossing and genetic diversification. Our study has several limitations. First, it was not possible to determine the initial abdominal status of reared mosquitoes. Furthermore, the small sample size of reared infected abdomens limited our ability to make inferences about differential onward *P. falciparum* transmission to mosquitoes. However, despite these limitations we were still able to estimate the maximum human-to-mosquito transmission efficiency to be between 0.14-0.27, which aligns with what has been observed in direct- and membrane-feeding experiments (11, 67, 68). Second, we could not match all mosquito STR profiles to people; however, our ability to incorporate a portion of incomplete and multisource profiles into our analysis increased our sample size by 37% (28), allowing us to more accurately use the matched bloodmeals and human infection data to estimate potential contributions to the transmission reservoir of different groups of people. Finally, these estimates assume equal transmission efficiency and gametocyte carriage amongst groups, which several studies have shown are also heterogeneous and nonrandom. Nonetheless, our large sample size of matched bloodmeals representing longitudinal observations of the same participants over several transmission seasons provide robust insights into mosquito biting behaviors and the potential impact on transmission and parasite diversity. In conclusion, we observed that school-age boys are disproportionately bitten and disproportionately contribute to onward transmission to mosquitoes, and that infectious mosquitoes are more likely to bite infected individuals. Our approach to measuring and quantifying mosquito biting bias in a natural setting contributes to accurate estimation of parasite transmission dynamics, and these observations provide insight into parasite diversity and evolution as well as the development of transmission-reducing interventions. ## Supporting information Supplemental Information [[supplements/301433_file03.pdf]](pending:yes) ## Data Availability All code for analyses is available at https://github.com/duke-malaria-collaboratory/imbibe_manuscript. Study data are available upon request. [https://github.com/duke-malaria-collaboratory/imbibe\_manuscript](https://github.com/duke-malaria-collaboratory/imbibe_manuscript) ## Funding This work was supported by NIAID (R01AI146849 to W.P.O. and S.M.T. and F32AI149950 and K01AI175527 to C.F.M.). ### Competing interests The authors declare no competing interests. ### Author contributions WPO, SMT, AAO, JNM, AW, and CFM acquired funding. JNM, AAO, WPO, LA performed project administration. EO, TC, MA, EF, and CFM performed the investigations. EK curated the data. CFM, ZL, SMT, and WPO conceptualized the project and wrote the original manuscript draft. CFM and ZL developed the methodology and software, and performed the formal analysis, and visualization. TK, AW, WPO, and SMT supervised the project. All authors reviewed and edited the manuscript. ## Acknowledgements We thank field technicians Ibrahim Khaoya, Lucas Marango, Ezna Mukeli, Eric Nalianya, Jane Nyongesa, Lilian Nukewa, Edith Wamalwa, and Aggrey Wekesa for their engagement with the study participants; Joseph Kipkoech Kirui, Sarah Laing, Julius Maiyo and Emily Robie for operational assistance and coordination; Kelsey Sumner, Erica Zeno, and Sophie Berube for scientific discussion; Thynn Thane, Jillian Grassia, Jenna Decurzio, Laura-Leigh Rowlette and Scott Langdon for sample processing; and Jamie Mills, Robert Rono, Francis Kithuku and Nikita Poujai for administrative support. Ultimately, we are indebted to the study household members for their participation in this study. ## Footnotes * # Co-first authors * * Co-senior authors * Received January 24, 2024. * Revision received January 24, 2024. * Accepted January 25, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. 1. D. L. Smith, F. E. McKenzie, R. W. Snow, S. I. Hay, Revisiting the Basic Reproductive Number for Malaria and Its Implications for Malaria Control. PLoS Biol. 5, e42 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.0050042&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17311470&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 2. 2. M. E. J. Woolhouse, et al., Heterogeneities in the transmission of infectious agents: Implications for the design of control programs. Proc. Natl. Acad. Sci. 94, 338–342 (1997). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czo4OiI5NC8xLzMzOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAxLzI1LzIwMjQuMDEuMjQuMjQzMDE0MzMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 3. 3. J. Martinez, A. Showering, C. Oke, R. T. Jones, J. G. Logan, Differential attraction in mosquito–human interactions and implications for disease control. Philos. Trans. R. Soc. B Biol. Sci. 376, 20190811 (2021). 4. 4. Y. T. Qiu, R. C. Smallegange, J. J. A. Van Loon, C. J. F. Ter Braak, W. Takken, Interindividual variation in the attractiveness of human odours to the malaria mosquito Anopheles gambiae s. s. Med. Vet. Entomol. 20, 280–287 (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2915.2006.00627.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17044878&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241155500003&link_type=ISI) 5. 5. D. Giraldo, et al., Human scent guides mosquito thermotaxis and host selection under naturalistic conditions. Curr. Biol. CB 33, 2367–2382.e7 (2023). 6. 6. R. C. Smallegange, et al., Malaria Infected Mosquitoes Express Enhanced Attraction to Human Odor. PLOS ONE 8, e63602 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0063602&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23691073&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 7. 7. N. M. Stanczyk, et al., Species-specific alterations in Anopheles mosquito olfactory responses caused by Plasmodium infection. Sci. Rep. 9, 3396 (2019). 8. 8. B. P. Gonçalves, et al., Examining the human infectious reservoir for Plasmodium falciparum malaria in areas of differing transmission intensity. Nat. Commun. 8, 1133 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-01270-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29074880&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 9. 9. J. T. Lin, et al., Microscopic Plasmodium falciparum Gametocytemia and Infectivity to Mosquitoes in Cambodia. J. Infect. Dis. 213, 1491–1494 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv599&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26667316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 10. 10. C. Andolina, et al., Sources of persistent malaria transmission in a setting with effective malaria control in eastern Uganda: a longitudinal, observational cohort study. Lancet Infect. Dis. 21, 1568–1578 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1473-3099(21)00072-4&link_type=DOI) 11. 11. A. Ahmad, et al., Infectivity of patent Plasmodium falciparum gametocyte carriers to mosquitoes: establishing capacity to investigate the infectious reservoir of malaria in a low-transmission setting in The Gambia. Trans. R. Soc. Trop. Med. Hyg. 115, 1462–1467 (2021). 12. 12. J. Rek, et al., Asymptomatic School-Aged Children Are Important Drivers of Malaria Transmission in a High Endemicity Setting in Uganda. J. Infect. Dis. 226, 708–713 (2022). 13. 13. T. Bousema, et al., Mosquito Feeding Assays to Determine the Infectiousness of Naturally Infected Plasmodium falciparum Gametocyte Carriers. PLOS ONE 7, e42821 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0042821&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22936993&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 14. 14. A. C. Gamage-Mendis, J. Rajakaruna, R. Carter, K. N. Mendis, Infectious Reservoir of Plasmodium Vivax and Plasmodium Falciparum Malaria in an Endemic Region of Sri Lanka. Am. J. Trop. Med. Hyg. 45, 479–487 (1991). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI0NS80LzQ3OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAxLzI1LzIwMjQuMDEuMjQuMjQzMDE0MzMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 15. 15. C. J. Drakeley, N. I. J. Akim, R. W. Sauerwein, B. M. Greenwood, G. A. T. Targett, Estimates of the infectious reservoir of Plasmodium falciparum malaria in The Gambia and in Tanzania. Trans. R. Soc. Trop. Med. Hyg. 94, 472–476 (2000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0035-9203(00)90056-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11132369&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 16. 16. A. Lin Ouédraogo, et al., Dynamics of the Human Infectious Reservoir for Malaria Determined by Mosquito Feeding Assays and Ultrasensitive Malaria Diagnosis in Burkina Faso. J. Infect. Dis. 213, 90–99 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv370&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26142435&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 17. 17. T. W. Scott, A. K. Githeko, A. Fleisher, L. C. Harrington, G. Yan, DNA profiling of human blood in anophelines from lowland and highland sites in western Kenya. Am. J. Trop. Med. Hyg. 75, 231–237 (2006). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI3NS8yLzIzMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAxLzI1LzIwMjQuMDEuMjQuMjQzMDE0MzMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 18. 18. S. Soremekun, et al., Measuring the efficacy of insecticide treated bednets: the use of DNA fingerprinting to increase the accuracy of personal protection estimates in Tanzania. Trop. Med. Int. Health 9, 664–672 (2004). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-3156.2004.01250.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15189456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222300500003&link_type=ISI) 19. 19. J. B. Keven, et al., Nonrandom Selection and Multiple Blood Feeding of Human Hosts by Anopheles Vectors: Implications for Malaria Transmission in Papua New Guinea. Am. J. Trop. Med. Hyg. 105, 1747–1758 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.21-0210&link_type=DOI) 20. 20. R. B. Mbewe, et al., Genotyping of Anopheles mosquito blood meals reveals nonrandom human host selection: implications for human-to-mosquito Plasmodium falciparum transmission. Malar. J. 22, 115 (2023). 21. 21. W. M. Guelbéogo, et al., Variation in natural exposure to anopheles mosquitoes and its effects on malaria transmission. eLife 7, e32625 (2018). 22. 22. L. C. Norris, C. M. Fornadel, W.-C. Hung, F. J. Pineda, D. E. Norris, Frequency of Multiple Blood Meals Taken in a Single Gonotrophic Cycle by Anopheles arabiensis Mosquitoes in Macha, Zambia. Am. J. Trop. Med. Hyg. 83, 33–37 (2010). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo3OiI4My8xLzMzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMjUvMjAyNC4wMS4yNC4yNDMwMTQzMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 23. 23. E. E. Zeno, et al., Risk of Malaria Following Untreated Subpatent Plasmodium falciparum Infections: Results Over 4 Years From a Cohort in a High-Transmission Area in Western Kenya. J. Infect. Dis., jiad398 (2023). 24. 24. K. M. Sumner, et al., Genotyping cognate Plasmodium falciparum in humans and mosquitoes to estimate onward transmission of asymptomatic infections. Nat. Commun. 12, 909 (2021). 25. 25. Z. Lapp, et al., Analytic optimization of Plasmodium falciparum marker gene haplotype recovery from amplicon deep sequencing of complex mixtures. 2023.08.17.23294237 (2023). 26. 26. G. E. Truett, et al., Preparation of PCR-quality mouse genomic DNA with hot sodium hydroxide and tris (HotSHOT). BioTechniques 29, 52, 54 (2000). 27. 27. E. E. Wilkins, P. I. Howell, M. Q. Benedict, IMP PCR primers detect single nucleotide polymorphisms for Anopheles gambiae species identification, Mopti and Savanna rDNA types, and resistance to dieldrin in Anopheles arabiensis. Malar. J. 5, 125 (2006). 28. 28. Z. Lapp, et al., bistro: An R package for vector bloodmeal identification by short tandem repeat overlap. Methods Ecol. Evol. 00, 1–9 (2023). 29. 29. A. Signorell, DescTools: Tools for Descriptive Statistics (2023). 30. 30. M. E. Brooks, et al., glmmTMB Balances Speed and Flexibility Among Packages for Zero-inflated Generalized Linear Mixed Modeling. R J. 9, 378–400 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.32614/RJ-2017-066&link_type=DOI) 31. 31. Y. Croissant, Estimation of Random Utility Models in R: The mlogit Package. J. Stat. Softw. 95, 1–41 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v095.i01&link_type=DOI) 32. 32.Posit team, RStudio: Integrated Development Environment for R (2023). 33. 33. R Core Team, R: A Language and Environment for Statistical Computing (2023). 34. 34. H. Wickham, et al., Welcome to the Tidyverse. J. Open Source Softw. 4, 1686 (2019). 35. 35. H. Wickham, E. Miller, D. Smith, haven: Import and Export “SPSS”, “Stata” and “SAS” Files (2023). 36. 36. F. E. Harrell Jr., Hmisc: Harrell Miscellaneous (2023). 37. 37. A. Kassambara, ggpubr: “ggplot2” Based Publication Ready Plots (2023). 38. 38. A. Cameron, T. van den Brand, geomtextpath: Curved Text in “ggplot2” (2022). 39. 39. R. Iannone, DiagrammeR: Graph/Network Visualization (2023). 40. 40. D. Robinson, A. Hayes, S. Couch, broom: Convert Statistical Objects into Tidy Tibbles (2023). 41. 41. B. Bolker, D. Robinson, broom.mixed: Tidying Methods for Mixed Models (2022). 42. 42.D. Lüdecke, ggeffects: Tidy Data Frames of Marginal Effects from Regression Models. J. Open Source Softw. 3, 772 (2018). 43. 43. D. D. Sjoberg, K. Whiting, M. Curry, J. A. Lavery, J. Larmarange, Reproducible Summary Tables with the gtsummary Package. R J. 13, 570–580 (2021). 44. 44. H. Jeppson, H. Hofmann, D. Cook, ggmosaic: Mosaic Plots in the “ggplot2” Framework (2021). 45. 45. A. Karlsson, M. Clements, biostat3: Utility Functions, Datasets and Extended Examples for Survival Analysis (2023). 46. 46. H. Wickham, modelr: Modelling Functions that Work with the Pipe (2023). 47. 47. J. W. Wekesa, R. S. Copeland, R. W. Mwangi, Effect of Plasmodium Falciparum on Blood Feeding Behavior of Naturally Infected Anopheles Mosquitoes in Western Kenya. Am. J. Trop. Med. Hyg. 47, 484–488 (1992). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI0Ny80LzQ4NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAxLzI1LzIwMjQuMDEuMjQuMjQzMDE0MzMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 48. 48. R. A. Anderson, J. C. Koellaf, H. Hurd, The effect of Plasmodium yoelii nigeriensis infection on the feeding persistence of Anopheles stephensi Liston throughout the sporogonic cycle. Proc. R. Soc. Lond. B Biol. Sci. 266, 1729–1733 (1999). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1999.0839&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10518321&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000082866800003&link_type=ISI) 49. 49. J. C. Koella, M. J. Packer, Malaria parasites enhance blood-feeding of their naturally infected vector Anopheles punctulatus. Parasitology 113, 105–109 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182000066348&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8760311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 50. 50. . J. c. Koella, F. l. Sørensen, R. a. Anderson, The malaria parasite, Plasmodium falciparum, increases the frequency of multiple feeding of its mosquito vector, Anopheles gambiae. Proc. R. Soc. Lond. B Biol. Sci. 265, 763–768 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1998.0358&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9628035&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000073580200004&link_type=ISI) 51. 51. C. M. De Moraes, et al., Malaria-induced changes in host odors enhance mosquito attraction. Proc. Natl. Acad. Sci. 111, 11079–11084 (2014). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTExLzMwLzExMDc5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMjUvMjAyNC4wMS4yNC4yNDMwMTQzMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 52. 52. R. Lacroix, W. R. Mukabana, L. C. Gouagna, J. C. Koella, Malaria Infection Increases Attractiveness of Humans to Mosquitoes. PLOS Biol. 3, e298 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.0030298&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16076240&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 53. 53. A. Robinson, et al., Plasmodium-associated changes in human odor attract mosquitoes. Proc. Natl. Acad. Sci. 115, E4209–E4218 (2018). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE1LzE4L0U0MjA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMjUvMjAyNC4wMS4yNC4yNDMwMTQzMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 54. 54. S. N. Emami, et al., A key malaria metabolite modulates vector blood seeking, feeding, and susceptibility to infection. Science 355, 1076–1080 (2017). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNTUvNjMyOS8xMDc2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMjUvMjAyNC4wMS4yNC4yNDMwMTQzMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 55. 55. M. Hajkazemian, et al., Mosquito host-seeking diel rhythm and chemosensory gene expression is affected by age and Plasmodium stages. Sci. Rep. 12, 18814 (2022). 56. 56. L. J. Cator, P. A. Lynch, A. F. Read, M. B. Thomas, Do malaria parasites manipulate mosquitoes? Trends Parasitol. 28, 466–470 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2012.08.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23044288&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310939800003&link_type=ISI) 57. 57. J. E. Coalson, et al., Simulation models predict that school-age children are responsible for most human-to-mosquito Plasmodium falciparum transmission in southern Malawi. Malar. J. 17, 147 (2018). 58. 58. W. Stone, B. P. Gonçalves, T. Bousema, C. Drakeley, Assessing the infectious reservoir of falciparum malaria: past and future. Trends Parasitol. 31, 287–296 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2015.04.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25985898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 59. 59. M. K. Mwandagalirwa, et al., Individual and household characteristics of persons with Plasmodium falciparum malaria in sites with varying endemicities in Kinshasa Province, Democratic Republic of the Congo. Malar. J. 16, 456 (2017). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 60. 60. J. E. Coalson, et al., High prevalence of Plasmodium falciparum gametocyte infections in school-age children using molecular detection: patterns and predictors of risk from a cross-sectional study in southern Malawi. Malar. J. 15, 527 (2016). 61. 61. J. A. Walldorf, et al., School-Age Children Are a Reservoir of Malaria Infection in Malawi. PLOS ONE 10, e0134061 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0134061&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26207758&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 62. 62. Z. Zhou, et al., Assessment of submicroscopic infections and gametocyte carriage of Plasmodium falciparum during peak malaria transmission season in a community-based cross-sectional survey in western Kenya, 2012. Malar. J. 15, 421 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-016-1482-4&link_type=DOI) 63. 63. A. Miles, et al., Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum. Genome Res. 26, 1288–1299 (2016). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjI2LzkvMTI4OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAxLzI1LzIwMjQuMDEuMjQuMjQzMDE0MzMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 64. 64. T. W. Scott, A. K. Githeko, A. Fleisher, L. C. Harrington, G. Yan, DNA profiling of human blood in anophelines from lowland and highland sites in western Kenya. Am. J. Trop. Med. Hyg. 75, 231–237 (2006). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI3NS8yLzIzMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAxLzI1LzIwMjQuMDEuMjQuMjQzMDE0MzMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 65. 65. L. A. Baton, L. C. Ranford-Cartwright, Spreading the seeds of million-murdering death. Trends Parasitol. 21, 573–580 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2005.09.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16236552&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234320300011&link_type=ISI) 66. 66. S. Portugal, et al., Host-mediated regulation of superinfection in malaria. Nat. Med. 17, 732–737 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.2368&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21572427&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) 67. 67. M. Diallo, et al., Evaluation and optimization of membrane feeding compared to direct feeding as an assay for infectivity. Malar. J. 7, 248 (2008). 68. 68. S. Bonnet, C. Gouagna, I. Safeukui, J.-Y. Meunier, C. Boudin, Comparison of artificial membrane feeding with direct skin feeding to estimate infectiousness of Plasmodium falciparum gametocyte carriers to mosquitoes. Trans. R. Soc. Trop. Med. Hyg. 94, 103– 106 (2000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0035-9203(00)90456-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10748913&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F25%2F2024.01.24.24301433.atom) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif