# Both consumptive and non-consumptive effects of predators impact mosquito populations and have implications for disease transmission

1. Marie C Russell marie.clare.russell@gmail.com 1
2. Catherine M Herzog 2
3. Zachary Gajewski 3
4. Chloe Ramsay 4
5. Fadoua El Moustaid 3
6. Michelle V Evans 5 6
7. Trishna Desai 7
8. Nicole L Gottdenker 5 8
9. Sara L Hermann 9
10. Alison G Power 10
11. Andrew C McCall 11

1. Department of Life Sciences, Imperial College London, Silwood Park Campus
2. Center for Infectious Disease Dynamics, Pennsylvania State University
3. Department of Biological Sciences, Virginia Polytechnic Institute and State University
4. Department of Biological Sciences, University of Notre Dame
5. Odum School of Ecology & Center for Ecology of Infectious Diseases, University of Georgia
6. MIVEGEC, IRD, CNRS, Université Montpellier
7. Nuffield Department of Population Health, University of Oxford
8. Department of Veterinary Pathology, University of Georgia College of Veterinary Medicine
9. Department of Entomology, Pennsylvania State University
10. Department of Ecology & Evolutionary Biology, Cornell University
11. Biology Department, Denison University

2022-01-19

- Research Article

- Ecology
- Epidemiology and Global Health

- vector ecology
- predation
- disease dynamics
- meta-analysis
- Mosquito

- publisher-id 71503
- doi 10.7554/eLife.71503
- elocation-id e71503

## Abstract

Predator-prey interactions influence prey traits through both consumptive and non-consumptive effects, and variation in these traits can shape vector-borne disease dynamics. Meta-analysis methods were employed to generate predation effect sizes by different categories of predators and mosquito prey. This analysis showed that multiple families of aquatic predators are effective in consumptively reducing mosquito survival, and that the survival of Aedes, Anopheles, and Culex mosquitoes is negatively impacted by consumptive effects of predators. Mosquito larval size was found to play a more important role in explaining the heterogeneity of consumptive effects from predators than mosquito genus. Mosquito survival and body size were reduced by non-consumptive effects of predators, but development time was not significantly impacted. In addition, Culex vectors demonstrated predator avoidance behavior during oviposition. The results of this meta-analysis suggest that predators limit disease transmission by reducing both vector survival and vector size, and that associations between drought and human West Nile virus cases could be driven by the vector behavior of predator avoidance during oviposition. These findings are likely to be useful to infectious disease modelers who rely on vector traits as predictors of transmission.

## Introduction

While it is well known that predation reduces vector populations through consumptive effects, non-consumptive effects of predators can also greatly impact prey demographics 106 Preisser et al.2004,12,31. Mosquitoes are vectors of a variety of debilitating and deadly diseases, including malaria, lymphatic filariasis, and arboviruses, such as chikungunya, Zika, and dengue 147 Weaver and Reisen 2009,12,31 148 WORLD HEALTH ORGANIZATION 2019,12,31. Consequently, there is motivation from a public health perspective to better understand the different drivers of variation in mosquito traits that can ultimately impact vector population growth and disease transmission. In addition, recent work has suggested that incorporation of vector trait variation into disease models can improve the reliability of their predictions 23 Cator et al.2019,12,31. In this study, systematic review and meta-analysis methods are used to synthesize a clearer understanding of the consumptive and non-consumptive effects of predators on mosquito traits, including survival, oviposition, development, and size.

Mosquito insecticide resistance is recognized as a growing problem 54 Hancock et al.2017,12,31 57 Hemingway and Ranson 1999,12,31 82 Liu 2014,12,31 leading some to suggest that control efforts should rely more heavily on ‘non-insecticide based strategies’ 14 Benelli et al.2015,12,31. The consumptive effects of predators on mosquitoes have previously been harnessed for biocontrol purposes. Past biocontrol efforts have used predators such as cyclopoid copepods 69 Kay et al.2001,12,31 87 Marten 1989,12,31 120 Russell et al.1995,12,31 142 Veronesi et al.2014,12,31 and mosquitofish 107 Pyke 2007,12,31 121 Seale 1916,12,31 to target the mosquito’s aquatic larval stage. The strength of the consumptive effects of these predators on mosquitoes can be influenced by multiple factors, including predator-prey size ratio and temperature. Predator-prey body size ratios tend to be higher in freshwater habitats than other types of habitats 18 Brose et al.2005,12,31, and attack rate tends to increase with temperature 68 Kalinoski and DeLong 2015,12,31 29 Dam and Peterson 1987,12,31, although other studies suggest a unimodal response to temperature 139 Uiterwaal and Delong 2019,12,31 41 Englund et al.2010,12,31.

Predators can also have non-consumptive effects on prey 103 Peacor and Werner 2000,12,31, and these effects are thought to be more pronounced in aquatic ecosystems than in terrestrial ecosystems 106 Preisser et al.2004,12,31. Non-consumptive effects of predators are the result of the prey initiating anti-predator behavioral and/or physiological trait changes that can aid in predator avoidance 59 Hermann and Landis 2016,12,31 80 Lima and Dill 1989,12,31. Such plasticity in certain prey traits may also result in energetic costs 81 Lima 1997,12,31. Predator detection is key for these trait changes to occur and can be mediated by chemical, tactile, and visual cues 58 Hermann and Thaler 2013,12,31. In mosquitoes, exposure to predators is known to affect a variety of traits including behavior, size, development, and survival 10 Arav and Blaustein 2005,12,31 17 Bond et al.2004,12,31 114 Roberts 2011,12,31 119 Roux et al.2014,12,31 151 Zuharah et al.2012,12,31. Experimental observations of predator effects on mosquito size and development are inconsistent and results sometimes vary by mosquito sex. For example, exposure to predation was found to increase the size of Culex pipiens mosquitoes 2 Alcalay et al.2017,12,31 but decrease the size of Culiseta longiareolata null 128 Stav et al.2004,12,31. In addition, female Aedes triseriatus exhibited shorter development times when exposed to predation at high nutrient availability 100 Ower and Juliano 2018,12,31, but male C. longiareolata had longer development times in the presence of predators 128 Stav et al.2004,12,31. In some cases, a shared evolutionary history between predator and prey organisms can strengthen the non-consumptive effects of predators on mosquitoes 21 Buchanan et al.2016,12,31 124 Sih 1985,12,31.

This investigation assesses the consumptive and non-consumptive effects of predators on mosquito traits and describes how these effects could impact disease transmission. The roles of vector genus, predator family, mosquito larval instar (an indicator of prey size), and temperature are also examined as potential moderators of predator effects. Non-consumptive effects of predators are expected to cause a smaller reduction in mosquito survival than consumptive effects because, in practice, measures of consumptive effects always include both consumptive and non-consumptive effects. Based on previous findings, larger predators are more likely to consumptively reduce mosquito survival 74 Kumar et al.2007,12,31. In addition, Aedes mosquito larvae may be more vulnerable to consumption than other genera because of the high degree of motility observed in this genus 35 Dieng et al.2002,12,31 89 Marten and Reid 2006,12,31 126 Soumare and Cilek 2010,12,31. The oviposition response to predation is expected to be weakest among Aedes species that oviposit above the water line, due in part to their delayed-hatching eggs 146 Vonesh and Blaustein 2009,12,31. Predation is predicted to reduce mosquito size and lengthen development time, consistent with the reduced growth response observed in other insect systems 59 Hermann and Landis 2016,12,31. Certain non-consumptive effects of predation, particularly oviposition site selection and decreased vector size, are likely to play important roles in the dynamics of mosquito-borne disease.

## Materials and methods

### Literature screening

A systematic search was conducted for studies on predation of mosquitoes that were published between 1970 and July 1, 2019 using both PubMed and Web of Science search engines, according to the PRISMA protocol 93 Moher et al.2008,12,31. Mosquito vectors of the Anopheles and Aedes genera were specifically highlighted in our search terms because these genera contain the vector species that transmit malaria, yellow fever, and dengue – the three most deadly mosquito-borne diseases worldwide 62 Hill et al.2004,12,31. Searches included 18 combinations of three vector predation terms (mosquito predat*, Anopheles predat*, Aedes predat*) and six trait terms (survival, mortality, development, fecundity, dispers*, host preference). Abstracts from the 1136 studies were each screened by two different co-authors, using the ‘metagear’ package in R 76 Lajeunesse 2015,12,31 108 R Development Core Team 2019,12,31. If either screener thought the study had information relevant to predation of mosquitoes, or both screeners thought the abstract was ambiguous, the study was read in full. This resulted in 306 studies that were fully reviewed to determine if any predation data could be extracted (Figure 1).

Figure 1.

#### Flowchart demonstrating the literature search, screening process, data exclusions, and the resulting seven different vector trait data subsets.

### Study exclusion criteria

Data were extracted from studies that collected data on non-consumptive and/or consumptive effects of predators on mosquitoes. Studies were required to have a mean, error measurement, and at least two replicates for both control and predator treatments. The control treatment was required to have all the same conditions as the predator treatment, such as prey density and type of water, without the predators. Studies that were not published in English and studies that did not differentiate between predators of multiple families were excluded. Studies were also excluded if oviposition by free-flying female mosquitoes could have interfered with observing the consumptive effects of predators on vector survival. The final database comprised data extracted from 60 studies (Supplementary file 1). The data included observations from laboratory experiments, as well as semi-field experiments, in which mesocosms of different treatments were observed in outdoor settings.

### Data extraction

Variables related to the publication, the vector, the predator, and the effect size (Table 1) were extracted from each study. Data from tables and text were recorded as they were published, and data from figures were extracted using WebPlotDigitizer 116 Rohatgi 2019,12,31. Error measurements that were not originally presented as standard deviations were converted to standard deviations prior to the effect size calculation.

Table 1.

#### Variables extracted from included studies.

Variable | Description
--- | ---
Publication data: |
Title | Full study title
Journal | Name of journal that published the study
Year | Year of publication
Study environment | Environment where the experiment took place: lab or semi-field
Vector data: |
Order, Family, Genus, Species | Taxonomic identification
Trait | Outcome that was measured (e.g. survival, development, etc.)
Stage | Life stage: egg, larva, pupa, or adult
Larval instar | Early (1st and 2nd instars), late (3rd and 4th instars), both, or NA (eggs, pupae, or adults)
Sex | Male or female
Predator data: |
Phylum, Class, Order, Family, Genus, Species | Taxonomic identification
Starved | Whether the predator was starved: yes or no
Time starved | Amount of time that the predator was starved (in minutes)
Predation effect | Consumptive or non-consumptive
Effect size data: |
Units | Units of extracted data
Control mean | Average of the outcome measured among the controls
Control standard deviation | Standard deviation of the outcome measured in the controls
Control number of replicates | Number of control replicates
Predation mean | Average of the outcome measured in the predator treatment
Predation standard deviation | Standard deviation of the outcome measured in the predator treatment
Predation number of replicates | Number of predation replicates
Experiment ID | Alphabetic assignment to mark observations sharing a control group or representing the same prey individuals as originating from the same experiment
Additional data: |
Experiment time (days) | Duration of the experiment in days
Data source | Graph or text
Number of predators | Number of predators with access to prey, or ‘cue’ if there are no predators with direct access to prey
Number of prey (vectors) | Number of mosquito prey that are exposed to predation
Arena volume (mL) | Volume of the arena where prey encounter predators
Time exposed to predator(s) | Amount of time (in days) when the predator has direct access to the mosquito prey
Temperature (°C) | Temperature during the predation interaction
Type of predator cue | Predator cues, or cues from both predator(s) and dying conspecifics; NA for observations with a consumptive predation effect

### Data exclusions

A PRISMA plot of literature inclusion and exclusion is provided in Figure 1. Observations where insecticide was used were excluded because insecticides are known to interfere with consumptive and non-consumptive effects of predators 32 Delnat et al.2018,12,31 63 Janssens and Stoks 2011,12,31. In addition, observations from experiments with mosquito prey of two or more species were excluded because it was not possible to account for effects from apparent competition or prey-switching. Observations of vector fecundity, vector competence, behavioral traits other than oviposition, as well as observations where the vector trait was marked as ‘other’ were not analyzed because each of these traits were only recorded from three or fewer studies.

Due to protandry, the earlier emergence of males to maximize their reproductive success, mosquitoes respond to sex-specific selective forces that influence their development time and body size 71 Kleckner et al.1994,12,31. Under low resource conditions, female mosquitoes are likely to maximize body mass by extending their development time, whereas males tend to minimize their development time at the expense of lower body mass 71 Kleckner et al.1994,12,31. Observations of mosquito development time and body size in our database that were not sex-specific were excluded so that these vector traits could be analyzed while controlling for sex. In addition, among the observations of development time and body size, some predator means did not necessarily represent an evenly weighted average of the replicates. For example, if a total of 20 mosquitoes from three different predator replicates survived to adulthood, the mean development time and size of those 20 individuals may have been reported. To represent an evenly weighted average of the replicates, it is necessary to first calculate summary statistics among multiple individuals that emerge from the same replicate, and then report the average of the replicate-specific means. Observations that might have been influenced by uneven representation of replicates were excluded to prevent pseudo-replication from altering later meta-analyses.

For consumptive observations where life stage-specific survival was reported after more than 10 days of predator exposure, only data on survival marked by adult emergence were included for analysis. Effects observed among immature vector stages after such a long period of predator exposure were not analyzed because they could have resulted from a combination of non-consumptive effects on development, and consumptive effects on survival. Development time observations that were reported as the inverse of development time (units of days–1) were excluded because although their means could be converted to units of days, their standard deviations could not be converted to match units of days. In cases where multiple body sections of the same mosquitoes were measured to produce multiple size observations, only the wing measurement was included in the analysis to prevent pseudo-replication. Observations in which both the control and the predator treatments had standard deviations of zero were excluded because the meta-analysis methods did not support non-positive sampling variances.

### Exclusions and data substitutions for predator treatment means of zero

One study that was included in our database reported egg survival data as the hatch rate of field collected Culex pervigilans rafts 151 Zuharah et al.2012,12,31. However, mosquitoes have been shown to lay eggs independent of mating 101 O’Meara 1978,12,31, and hatch rates of zero have previously been observed in rafts laid by Culex females that were held separately from males 131 Su and Mulla 1996,12,31. Thus, hatch rates of zero were excluded from further analysis because these values may represent unfertilized egg rafts, rather than a strong impact of predators on survival. Twenty of the 187 consumptive survival observations had a predation mean of zero, and each of these zeros resulted from experiments that began with a specified number of live larvae. Consumptive survival zeros were each replaced with 0.5% of the starting number of mosquito prey to avoid undefined effect sizes. In addition, there was one zero out of the 36 oviposition predation means; this value had units of ‘number of egg rafts laid’ and was replaced with 0.5 rafts. Similar methods for replacing zero values in the treatment mean with small non-zero values have previously been employed 134 Thapa et al.2017,12,31.

The final analysis dataset included seven subsets: consumptive effects on survival, non-consumptive effects on survival, oviposition, development (female and male), and size (female and male). The data included 187 observations from 34 studies of consumptive survival, 24 observations from seven studies of non-consumptive survival, 36 observations from 12 studies of oviposition, 14 observations from seven studies of female development, 14 observations from seven studies of male development, 27 observations from 10 studies of female size, and 18 observations from nine studies of male size (Figure 1). These observations covered seven different classes of predator families (Figure 2).

Figure 2.

#### Mosquito predator classes (bold font) and families (italicized font) included in the database and the vector traits that they may influence (in parentheses); predator images not to scale, and placed randomly with respect to the different mosquito life stages.

Image sources: phylopic.org (CC BY 3.0 or public domain): Actinopterygii (creator: Milton Tan), Arachnida (creators: Sidney Frederic Harmer & Arthur Everett Shipley, vectorized by Maxime Dahirel), Branchiopoda (creator: Africa Gomez), and Insecta (creator: Marie Russell). BioRender.com: Amphibia, Hexanauplia, and Malacostraca class silhouettes; mosquito larval instars, pupa, and blood-feeding adult. Trishna Desai: mosquito egg raft.

### Data analysis

    options(stringsAsFactors = F, scipen = 999)

    library("metafor")
    library("tidyverse")

    ## Full database:
    database = read.csv("Predator_effects_on_mosquitoes_database.csv")

    length(unique(database$Study_ID))

    ## Data exclusions & zero substitutions:
    clean_data = database %>%
      #Insecticide used with predators
      filter(Pesticide!="YES") %>%
      #Multiple prey species
      filter(Mixed_prey!="YES") %>%
      #Traits and behavior categories that didn't have enough data for meta-analysis
      filter(Vec_trait!="OTHER" & Vec_trait!="FECUNDITY" & Vec_trait!="COMPETENCE") %>%
      filter(!(grepl("BEHAVIOR", Vec_trait) & !grepl("oviposition", Vec_trait))) %>%
      #Size or development data that were not specific to sex
      filter(!(Vec_trait=="DEVELOPMENT" & (Sex=="BOTH" | is.na(Sex)))) %>%
      filter(!(Vec_trait=="SIZE" & (Sex=="BOTH" | is.na(Sex)))) %>%
      #Pseudo-replication due to uneven representation of reps
      filter(!(grepl("individual means", Units))) %>%
      #Immature life stage-specific survival reported after more than ten days
      filter(!(Pred_effect=="CONSUMPTIVE" & Study_ID=="U106" & grepl("larval", Units))) %>%
      filter(!(Pred_effect=="CONSUMPTIVE" & Study_ID==1175 & is.na(Sex))) %>%
      filter(!(Pred_effect=="CONSUMPTIVE" & Study_ID=="U105" & grepl("larvae", Units))) %>%
      filter(!(Pred_effect=="CONSUMPTIVE" & Study_ID=="U105" & grepl("pupae", Units))) %>%
      #Inverse dev time
      filter(!(Vec_trait=="DEVELOPMENT" & grepl("1/median", Units))) %>%
      #Pseudo-replication due to measuring different body parts of the same mosquitoes
      filter(!(Study_ID=="1224" & grepl("thorax", Units))) %>%
      filter(!(Study_ID=="1226" & grepl("thorax", Units))) %>%
      filter(!(Study_ID=="U166" & grepl("femur", Other_trait))) %>%
      filter(!(Study_ID=="U166" & grepl("thorax", Other_trait))) %>%
      #Both the control sd and the predation sd are zeros. This will result in non-positive sampling variances for the ROM measure of effect.
      filter(!(Con_sd==0 & Pred_sd==0)) %>%
      #For oviposition observations with a predation treatment mean of zero, replace zero egg rafts with 0.5 egg rafts.
      mutate(Pred_mean = ifelse(Pred_mean==0 & (grepl("oviposition", Vec_trait)), 0.5, Pred_mean)) %>%
      #Delete the egg stage survival data if the predation treatment has a mean of zero. These could represent unfertilized egg rafts.
      filter(!(Vec_trait=="SURVIVAL" & Vec_stage=="EGG" & Pred_mean==0)) %>%
      #Replace larval survival observations that have zero predator means with 0.5% of the starting larvae.
      mutate(Pred_mean = ifelse(Pred_mean==0 & (grepl("%", Units) | grepl("percent", Units)), 0.5,
                                ifelse(Pred_mean==0 & (grepl("number", Units) | grepl("#", Units) | grepl("abundance", Units)), (0.005*Num_prey),
                                       ifelse(Pred_mean==0 & grepl("proportion", Units), 0.005, Pred_mean)))) %>%
      #Add the Effect ID column.
      mutate(Effect_ID = 1: n())

    ###############################################################################
    #########################  CONSUMPTIVE SURVIVAL  ##############################
    ###############################################################################

    consumptive_survival = clean_data %>%
      filter(Pred_effect == "CONSUMPTIVE")

    length(unique(consumptive_survival$Study_ID))

    #Calculate the Ratio of Means (ROM) effect size estimate.
    ROM_consumptive_survival = escalc(measure = "ROM",
                                      n1i = Pred_n,
                                      n2i = Con_n,
                                      m1i = Pred_mean,
                                      m2i = Con_mean,
                                      sd1i = Pred_sd,
                                      sd2i = Con_sd,
                                      data = consumptive_survival,
                                      var.names = c("ROM_est", "ROM_sv"),
                                      append=T)

    #Random effects model:
    ROM.uni_consumptive_survival <- rma.uni(yi = ROM_est,
                                            vi = ROM_sv,
                                            data = ROM_consumptive_survival,
                                            method = "REML")
    ROM.uni_consumptive_survival[["version"]]=NULL

    summary(ROM.uni_consumptive_survival)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_consumptive_survival)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_consumptive_survival, model = "rma", predictor="sei")

    #Multilevel mixed effects models:
    #No moderators
    ROM.no_mods <- rma.mv(yi = ROM_est,
                          V = ROM_sv,
                          random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                          data = ROM_consumptive_survival,
                          method = "REML")
    ROM.no_mods[["version"]]=NULL
    summary(ROM.no_mods)

    #Vector genus as a moderator
    ROM.vec <- rma.mv(yi = ROM_est,
                      V = ROM_sv,
                      mods = ~Vec_genus-1,
                      random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                      data = ROM_consumptive_survival,
                      method = "REML")
    ROM.vec[["version"]]=NULL
    summary(ROM.vec)

    #Predator family as a moderator
    ROM.pred_fam <- rma.mv(yi = ROM_est,
                           V = ROM_sv,
                           mods = ~Pred_family-1,
                           random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                           data = ROM_consumptive_survival,
                           method = "REML")
    ROM.pred_fam[["version"]]=NULL
    summary(ROM.pred_fam)

    #Predator family and vector genus as moderators
    ROM.pred_and_vec <- rma.mv(yi = ROM_est,
                               V = ROM_sv,
                               mods = ~Pred_family+Vec_genus -1,
                               random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                               data = ROM_consumptive_survival,
                               method = "REML")
    ROM.pred_and_vec[["version"]]=NULL
    summary(ROM.pred_and_vec)

    #Predator family and vector genus interacting as moderators
    ROM.pred_vec_interact <- rma.mv(yi = ROM_est,
                                    V = ROM_sv,
                                    mods = ~Pred_family*Vec_genus -1,
                                    random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                    data = ROM_consumptive_survival,
                                    method = "REML")
    ROM.pred_vec_interact[["version"]]=NULL
    summary(ROM.pred_vec_interact)

    #182 out of 187 consumptive observations, 97%, are of the larval stage.
    #Testing moderators on observations of larval survival only:

    larval_stage = consumptive_survival %>%
      filter(!(is.na(Larval_instar)))

    length(unique(larval_stage$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_larval_stage = escalc(measure = "ROM",
                              n1i = Pred_n,
                              n2i = Con_n,
                              m1i = Pred_mean,
                              m2i = Con_mean,
                              sd1i = Pred_sd,
                              sd2i = Con_sd,
                              data = larval_stage,
                              var.names = c("ROM_est", "ROM_sv"),
                              append=T)

    #Random effects model:
    ROM.uni_larval_stage <- rma.uni(yi = ROM_est,
                                    vi = ROM_sv,
                                    data = ROM_larval_stage,
                                    method = "REML")
    ROM.uni_larval_stage[["version"]]=NULL
    summary(ROM.uni_larval_stage)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_larval_stage)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_larval_stage, model = "rma", predictor="sei")

    #Multilevel mixed effects models:
    #No moderators
    ROM.larval_no_mods <- rma.mv(yi = ROM_est,
                                 V = ROM_sv,
                                 random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                 data = ROM_larval_stage,
                                 method = "REML")
    summary(ROM.larval_no_mods)

    #Larval instar as a moderator
    ROM.larval_instar <- rma.mv(yi = ROM_est,
                                V = ROM_sv,
                                mods = ~Larval_instar-1,
                                random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                data = ROM_larval_stage,
                                method = "REML")
    ROM.larval_instar[["version"]]=NULL
    summary(ROM.larval_instar)

    #Vector genus as a moderator
    ROM.larval_vec <- rma.mv(yi = ROM_est,
                             V = ROM_sv,
                             mods = ~Vec_genus-1,
                             random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                             data = ROM_larval_stage,
                             method = "REML")
    ROM.larval_vec[["version"]]=NULL
    summary(ROM.larval_vec)

    #Predator family as a moderator
    ROM.larval_pred <- rma.mv(yi = ROM_est,
                              V = ROM_sv,
                              mods = ~Pred_family-1,
                              random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                              data = ROM_larval_stage,
                              method = "REML")
    ROM.larval_pred[["version"]]=NULL
    summary(ROM.larval_pred)

    #Larval instar and predator family as moderators
    ROM.larval_instar_and_pred <- rma.mv(yi = ROM_est,
                                         V = ROM_sv,
                                         mods = ~Larval_instar+Pred_family-1,
                                         random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                         data = ROM_larval_stage,
                                         method = "REML")
    ROM.larval_instar_and_pred[["version"]]=NULL
    summary(ROM.larval_instar_and_pred)

    #Vector genus and predator family as moderators
    ROM.larval_pred_and_vec <- rma.mv(yi = ROM_est,
                                      V = ROM_sv,
                                      mods = ~Vec_genus+Pred_family-1,
                                      random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                      data = ROM_larval_stage,
                                      method = "REML")
    ROM.larval_pred_and_vec[["version"]]=NULL
    summary(ROM.larval_pred_and_vec)

    #Larval instar and predator family interacting as moderators
    ROM.larval_instar_pred_interact <- rma.mv(yi = ROM_est,
                                              V = ROM_sv,
                                              mods = ~Larval_instar*Pred_family-1,
                                              random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                              data = ROM_larval_stage,
                                              method = "REML")
    ROM.larval_instar_pred_interact[["version"]]=NULL
    summary(ROM.larval_instar_pred_interact)

    #Vector genus and predator family interacting as moderators
    ROM.larval_pred_vec_interact <- rma.mv(yi = ROM_est,
                                           V = ROM_sv,
                                           mods = ~Vec_genus*Pred_family-1,
                                           random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                           data = ROM_larval_stage,
                                           method = "REML")
    ROM.larval_pred_vec_interact[["version"]]=NULL
    summary(ROM.larval_pred_vec_interact)

    #Test moderators only on observations that included temperature.
    #The temperature range is from 16 to 28.25 degrees C.

    temp = consumptive_survival %>%
      filter(!(is.na(Temp))) %>%
      mutate(Temp_squared = Temp*Temp)

    length(unique(temp$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_temp = escalc(measure = "ROM",
                      n1i = Pred_n,
                      n2i = Con_n,
                      m1i = Pred_mean,
                      m2i = Con_mean,
                      sd1i = Pred_sd,
                      sd2i = Con_sd,
                      data = temp,
                      var.names = c("ROM_est", "ROM_sv"),
                      append=T)

    #Random effects model:
    ROM.uni_temp <- rma.uni(yi = ROM_est,
                            vi = ROM_sv,
                            data = ROM_temp,
                            method = "REML")
    ROM.uni_temp[["version"]]=NULL
    summary(ROM.uni_temp)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_temp)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_temp, model = "rma", predictor="sei")

    #Multilevel mixed effects models:
    #No moderators
    ROM.temp_no_mods <- rma.mv(yi = ROM_est,
                               V = ROM_sv,
                               random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                               data = ROM_temp,
                               method = "REML")
    ROM.temp_no_mods[["version"]]=NULL
    summary(ROM.temp_no_mods)

    #Temp as a moderator
    ROM.temp <- rma.mv(yi = ROM_est,
                       V = ROM_sv,
                       mods = ~Temp,
                       random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                       data = ROM_temp,
                       method = "REML")
    ROM.temp[["version"]]=NULL
    summary(ROM.temp)

    #Temp_squared as a moderator
    ROM.temp_squared <- rma.mv(yi = ROM_est,
                               V = ROM_sv,
                               mods = ~Temp_squared,
                               random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                               data = ROM_temp,
                               method = "REML")
    ROM.temp_squared[["version"]]=NULL
    summary(ROM.temp_squared)

    ###############################################################################
    ######################  NON-CONSUMPTIVE SURVIVAL  #############################
    ###############################################################################

    non_con_survival = clean_data %>%
      filter(Pred_effect=="NONCONSUMPTIVE") %>%
      filter(Vec_trait=="SURVIVAL")

    length(unique(non_con_survival$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_non_con_survival = escalc(measure = "ROM",
                                  n1i = Pred_n,
                                  n2i = Con_n,
                                  m1i = Pred_mean,
                                  m2i = Con_mean,
                                  sd1i = Pred_sd,
                                  sd2i = Con_sd,
                                  data = non_con_survival,
                                  var.names = c("ROM_est", "ROM_sv"),
                                  append=T)

    #Random effects model:
    ROM.uni_non_con_survival <- rma.uni(yi = ROM_est,
                                        vi = ROM_sv,
                                        data = ROM_non_con_survival,
                                        method = "REML")
    ROM.uni_non_con_survival[["version"]]=NULL
    summary(ROM.uni_non_con_survival)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_non_con_survival)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_non_con_survival, model = "rma", predictor="sei")

    #There is significant publication bias, so the "trimfill" function is used to impute values from missing studies.
    trimfill(ROM.uni_non_con_survival, estimator = "L0", side = "left")
    trimfill(ROM.uni_non_con_survival, estimator = "L0", side = "right")

    #Visualization of the three imputed values on the left side of the funnel plot:
    funnel(trimfill(ROM.uni_non_con_survival, estimator = "L0", side = "left"))

    ###############################################################################
    ##############################  OVIPOSITION  ##################################
    ###############################################################################

    oviposition = clean_data %>%
      filter(Vec_trait=="BEHAVIOR: oviposition")

    length(unique(oviposition$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_oviposition = escalc(measure = "ROM",
                             n1i = Pred_n,
                             n2i = Con_n,
                             m1i = Pred_mean,
                             m2i = Con_mean,
                             sd1i = Pred_sd,
                             sd2i = Con_sd,
                             data = oviposition,
                             var.names = c("ROM_est", "ROM_sv"),
                             append=T)

    #Random effects model:
    ROM.uni_oviposition <- rma.uni(yi = ROM_est,
                                   vi = ROM_sv,
                                   data = ROM_oviposition,
                                   method = "REML")
    ROM.uni_oviposition[["version"]]=NULL
    summary(ROM.uni_oviposition)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_oviposition)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_oviposition, model = "rma", predictor="sei")

    #Multilevel mixed effects models:
    #No moderators
    ROM.oviposition_no_mods <- rma.mv(yi = ROM_est,
                                      V = ROM_sv,
                                      random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                      data = ROM_oviposition,
                                      method = "REML")
    ROM.oviposition_no_mods[["version"]]=NULL
    summary(ROM.oviposition_no_mods)

    #Vector genus as a moderator
    ROM.oviposition_vec <- rma.mv(yi = ROM_est,
                                  V = ROM_sv,
                                  mods = ~Vec_genus-1,
                                  random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                  data = ROM_oviposition,
                                  method = "REML")
    ROM.oviposition_vec[["version"]]=NULL
    summary(ROM.oviposition_vec)

    #Predator family as a moderator
    ROM.oviposition_pred <- rma.mv(yi = ROM_est,
                                   V = ROM_sv,
                                   mods = ~Pred_family-1,
                                   random = list(~1|Effect_ID, ~1|Study_ID/Experimental_ID),
                                   data = ROM_oviposition,
                                   method = "REML")
    ROM.oviposition_pred[["version"]]=NULL
    summary(ROM.oviposition_pred)

    ###############################################################################
    ##########################  FEMALE DEVELOPMENT  ###############################
    ###############################################################################

    development_females = clean_data %>%
      filter(Vec_trait=="DEVELOPMENT") %>%
      filter(Sex=="FEMALE")

    length(unique(development_females$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_development_females = escalc(measure = "ROM",
                                     n1i = Pred_n,
                                     n2i = Con_n,
                                     m1i = Pred_mean,
                                     m2i = Con_mean,
                                     sd1i = Pred_sd,
                                     sd2i = Con_sd,
                                     data = development_females,
                                     var.names = c("ROM_est", "ROM_sv"),
                                     append=T)

    #Random effects model:
    ROM.uni_development_females <- rma.uni(yi = ROM_est,
                                           vi = ROM_sv,
                                           data = ROM_development_females,
                                           method = "REML")

    ROM.uni_development_females[["version"]]=NULL
    summary(ROM.uni_development_females)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_development_females)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_development_females, model = "rma", predictor="sei")

    #There is significant publication bias, so the "trim-fill" function is used to impute values from missing studies.
    trimfill(ROM.uni_development_females, estimator = "L0", side = "left")
    trimfill(ROM.uni_development_females, estimator = "L0", side = "right")

    ###############################################################################
    ###########################  MALE DEVELOPMENT  ################################
    ###############################################################################

    development_males = clean_data %>%
      filter(Vec_trait=="DEVELOPMENT") %>%
      filter(Sex=="MALE")

    length(unique(development_males$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_development_males = escalc(measure = "ROM",
                                   n1i = Pred_n,
                                   n2i = Con_n,
                                   m1i = Pred_mean,
                                   m2i = Con_mean,
                                   sd1i = Pred_sd,
                                   sd2i = Con_sd,
                                   data = development_males,
                                   var.names = c("ROM_est", "ROM_sv"),
                                   append=T)

    #Random effects model:
    ROM.uni_development_males <- rma.uni(yi = ROM_est,
                                         vi = ROM_sv,
                                         data = ROM_development_males,
                                         method = "REML")
    ROM.uni_development_males[["version"]]=NULL
    summary(ROM.uni_development_males)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_development_males)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_development_males, model = "rma", predictor="sei")

    #There is significant publication bias, so the "trimfill" function is used to impute values from missing studies.
    trimfill(ROM.uni_development_males, estimator = "L0", side = "left")
    trimfill(ROM.uni_development_males, estimator = "L0", side = "right")

    ###############################################################################
    ##############################  FEMALE SIZE  ##################################
    ###############################################################################

    size_females = clean_data %>%
      filter(Vec_trait=="SIZE") %>%
      filter(Sex=="FEMALE")

    length(unique(size_females$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_size_females = escalc(measure = "ROM",
                              n1i = Pred_n,
                              n2i = Con_n,
                              m1i = Pred_mean,
                              m2i = Con_mean,
                              sd1i = Pred_sd,
                              sd2i = Con_sd,
                              data = size_females,
                              var.names = c("ROM_est", "ROM_sv"),
                              append=T)

    #Random effects model:
    ROM.uni_size_females <- rma.uni(yi = ROM_est,
                                    vi = ROM_sv,
                                    data = ROM_size_females,
                                    method = "REML")

    ROM.uni_development_males[["version"]]=NULL
    summary(ROM.uni_size_females)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_size_females)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_size_females, model = "rma", predictor="sei")

    #There is significant publication bias, so the "trimfill" function is used to impute values from missing studies.
    trimfill(ROM.uni_size_females, estimator = "L0", side = "left")
    trimfill(ROM.uni_size_females, estimator = "L0", side = "right")

    #Visualization of the two imputed values on the right side of the funnel plot:
    funnel(trimfill(ROM.uni_size_females, estimator = "L0", side = "right"))

    ###############################################################################
    ###############################  MALE SIZE  ###################################
    ###############################################################################

    size_males = clean_data %>%
      filter(Vec_trait=="SIZE") %>%
      filter(Sex=="MALE")

    length(unique(size_males$Study_ID))

    #Calculate the Ratio of Means (ROM).
    ROM_size_males = escalc(measure = "ROM",
                            n1i = Pred_n,
                            n2i = Con_n,
                            m1i = Pred_mean,
                            m2i = Con_mean,
                            sd1i = Pred_sd,
                            sd2i = Con_sd,
                            data = size_males,
                            var.names = c("ROM_est", "ROM_sv"),
                            append=T)

    #Random effects model:
    ROM.uni_size_males <- rma.uni(yi = ROM_est,
                                  vi = ROM_sv,
                                  data = ROM_size_males,
                                  method = "REML")
    ROM.uni_size_males[["version"]]=NULL

    summary(ROM.uni_size_males)

    #Assess publication bias visually with a funnel plot.
    funnel(ROM.uni_size_males)

    #Use Egger's regression test to quantitatively assess publication bias.
    regtest(ROM.uni_size_males, model = "rma", predictor="sei")

#### Measuring effect sizes and heterogeneity

All analyses were conducted in R version 4.0.2 108 R Development Core Team 2019,12,31. For each subset of trait data (Figure 1), the ratio of means (ROM) measure of effect size was calculated using the ‘escalc’ function from the ‘metafor’ package; this effect measure is equal to a log-transformed fraction, where predation mean is the numerator and control mean is the denominator 143 Viechtbauer 2009,12,31. Random effects models, using the ‘rma.uni’ function, were run with the ROM effect sizes as response variables; each model had a normal error distribution and a restricted maximum likelihood (REML) estimator for τ2, the variance of the distribution of true effect sizes 143 Viechtbauer 2009,12,31. Although these random effects models could not account for multiple random effects or moderators, they provided overall estimates of the ROM effect sizes and estimates of the I2 statistics. Each I2 statistic represented the percentage of total variation across studies due to heterogeneity 60 Higgins et al.2002,12,31. If the I2 statistic was equal to or greater than 75%, the heterogeneity was considered to be high 60 Higgins et al.2002,12,31, and high heterogeneity has previously motivated further testing of moderators 144 Vincze et al.2016,12,31.

#### Assessing publication bias

Publication bias was assessed by visually inspecting funnel plots and conducting Egger’s regression test (‘regtest’ function) with standard error as the predictor 129 Sterne and Egger 2000,12,31 143 Viechtbauer 2009,12,31. If the Egger’s regression test showed significant evidence of publication bias based on funnel plot asymmetry, the ‘trim and fill’ method (‘trimfill’ function) was used to estimate how the predation effect size might change after imputing values from missing studies 37 Duval and Tweedie 1999,12,31 38 Duval and Tweedie 1999,12,31 143 Viechtbauer 2009,12,31. The trim and fill method has previously been recommended for testing the robustness of conclusions related to topics in ecology and evolution 64 Jennions and Møller 2001,12,31. Of the two trim and fill estimators, Rnull and Lnull, that were originally recommended 37 Duval and Tweedie 1999,12,31 38 Duval and Tweedie 1999,12,31, the Lnull estimator was used in this study because it is more appropriate for smaller datasets 123 Shi and Lin 2018,12,31.

#### Testing moderators

Data subsets that had high heterogeneity, observations from at least 10 studies, and no evidence of publication bias according to Egger’s regression results were analyzed further using multilevel mixed effects models with the ‘rma.mv’ function 143 Viechtbauer 2009,12,31 61 Higgins et al.2019,12,31. All multilevel mixed effects models had normal error distributions, REML estimators for τ2, and accounted for two random factors: effect size ID, and experiment ID nested within study ID. Moderators, such as predator family, vector genus, larval instar (directly correlated to prey size), and temperature, were tested within each data subset to determine if they affected the observed heterogeneity in ROM effect sizes. For categorical moderators, the intercept of the multilevel mixed effects model was removed, allowing an analysis of variance (ANOVA) referred to as the ‘test of moderators’ to indicate if any of the categories had an effect size different than zero. For data subsets with observations from 10 to 29 studies, only one moderator was tested at a time to account for sample size constraints. For subsets with observations from a higher number of studies (30 or more), up to two moderators were tested at once, and interaction between moderators was also tested. The small sample corrected Akaike Information Criterion (AICc) was used to compare multilevel mixed effects models and to select the model of best fit within each data subset; differences in AICc greater than two were considered meaningful 22 Burnham and Anderson 2003,12,31.

## Results

### Random effects models

Each data subset (Figure 1) had an I2 statistic of greater than 75%, indicating high heterogeneity 60 Higgins et al.2002,12,31. Random effects model results showed that predators consumptively decreased mosquito survival with an effect size of –1.23 (95% CI −1.43,–1.03), p-value < 0.0001, and non-consumptively reduced survival with a smaller effect size of –0.11 (95% CI −0.17,–0.04), p-value = 0.0016. In addition, predators non-consumptively reduced oviposition behavior with an effect size of –0.87 (95% CI −1.31,–0.42), p-value = 0.0001, and mosquito body size was non-consumptively reduced by predators in both males and females; the female effect size was –0.13 (95% CI −0.19,–0.06), p-value = 0.0002, and the male effect size was –0.03 (95% CI −0.06,–0.01), p-value = 0.0184. There was not a significant non-consumptive effect of predators on either male or female development time; the female effect size was –0.01 (95% CI –0.09, 0.07), p-value = 0.7901, and the male effect size was –0.04 (95% CI –0.12, 0.04), p-value = 0.3273.

The Egger’s regression test results showed that the non-consumptive survival subset, both development time subsets (male and female), and the female size subset exhibited funnel plot asymmetry indicative of publication bias. The ‘trim and fill’ procedure identified missing studies in the non-consumptive survival subset and the female size subset, but the procedure did not identify any missing studies in either of the development time subsets. Three studies were estimated to be missing from the non-consumptive survival data, and accounting for imputed values from missing studies resulted in a shift in the predation effect size from –0.11 (95% CI −0.17,–0.04), p-value = 0.0016, to -0.13 (95% CI −0.20,–0.07), p-value < 0.0001. Two studies were estimated to be missing from the female size data, and accounting for imputed values from these missing studies shifted the predation effect size from –0.13 (95% CI −0.19,–0.06), p-value = 0.0002, to -0.10 (95% CI −0.17,–0.03), p-value = 0.0083. Shifts in effect size estimates due to the trim and fill procedure were minor and did not cause any of the observed effects of predators to change direction or become insignificant.

### Multilevel mixed effects models

The consumptive survival and oviposition data subsets met the criteria of high heterogeneity, observations from at least 10 studies, and no evidence of publication bias. Therefore, these data subsets were tested for moderators using multilevel mixed effects models. Predator families that decreased mosquito survival included Cyprinidae: –3.44 (95% CI −5.79,–1.09), p-value = 0.0042; Poeciliidae: –1.42 (95% CI −2.67,–0.16), p-value = 0.0270; Ambystomatidae: –5.18 (95% CI −7.94,–2.42), p-value = 0.0002; Aeshnidae: –2.93 (95% CI −4.80,–1.07), p-value = 0.0020; and Notonectidae: –2.14 (95% CI −3.07,–1.21), p-value < 0.0001 (Figure 3a). Vector genera that experienced significant decreases in survival due to consumptive effects of predators included Aedes: –1.23 (95% CI −1.81,–0.65), p-value < 0.0001; Anopheles: –1.34 (95% CI −2.01,–0.66), p-value = 0.0001; and Culex: –1.41 (95% CI −1.96,–0.86), p-value < 0.0001 (Figure 3b). Among all 187 consumptive survival observations from 34 studies, the best model fit, according to AICc value, was achieved when an interaction between predator family and vector genus was included in the model (Table 2). However, among the 163 larval stage consumptive survival observations from 30 studies, adding an interactive term between larval instar (an indicator of prey size) and predator family had a greater improvement on model fit than adding an interactive term between vector genus and predator family (Figure 3c, Table 3). Temperature did not affect the heterogeneity of consumptive survival data, either as a linear moderator: –0.01 (95% CI –0.10, 0.07), p-value = 0.7559, or a quadratic moderator: 0.00 (95% CI 0.00, 0.00), p-value = 0.8184. The best oviposition model fit, according to AICc value, was achieved when vector genus was added as a moderator (Table 4). The mean oviposition effect size was not significantly different than zero for Aedes: 0.32 (95% CI –2.14, 2.79), p-value = 0.7970, or Culiseta: –0.61 (95% CI –1.83, 0.62), p-value = 0.3329, but for Culex mosquitoes, oviposition was significantly decreased by predator presence: –1.69 (95% CI −2.82,–0.56), p-value = 0.0033 (Figure 4).

Table 2.

#### Candidate multilevel mixed effects models of consumptive effects from predators on mosquito survival, fitted to dataset of effect sizes (n = 187 from 34 studies), and ranked by corrected Akaike’s information criterion (AICc).

Moderator(s) | Test of moderators(degrees of freedom, p-value) | AICc | ΔAICc
--- | --- | --- | ---
Predator family x vector genus | 28, < 0.0001 | 500.5 | null
Predator family | 19, < 0.0001 | 507 | 6.5
Predator family + vector genus | 23, < 0.0001 | 508.1 | 7.6
Vector genus | 5, < 0.0001 | 573 | 72.5
None | ---- | 576.5 | 76

Table 3.

#### Candidate multilevel mixed effects models of consumptive effects from predators, fitted to dataset of effect sizes where larval instar is not missing (n = 163 from 30 studies), and ranked by corrected Akaike’s information criterion (AICc).

Moderator(s) | Test of moderators(degrees of freedom, p-value) | AICc | ΔAICc
--- | --- | --- | ---
Predator family x larval instar | 25, < 0.0001 | 429.2 | null
Predator family + larval instar | 19, < 0.0001 | 443.5 | 14.3
Predator family x vector genus | 25, < 0.0001 | 455 | 25.8
Predator family | 17, < 0.0001 | 456.8 | 27.6
Predator family + vector genus | 21, < 0.0001 | 458.4 | 29.2
Larval instar | 3, < 0.0001 | 503.1 | 73.9
Vector genus | 5, < 0.0001 | 504.7 | 75.5
None | ---- | 508.5 | 79.3

Table 4.

#### Candidate multilevel mixed effects models of non-consumptive effects of predators on mosquito oviposition behavior, fitted to dataset of effect sizes (n = 36 from 12 studies), and ranked by corrected Akaike’s information criterion (AICc).

Moderator(s) | Test of moderators(degrees of freedom, p-value) | AICc | ΔAICc
--- | --- | --- | ---
Vector genus | 3, 0.0149 | 122.1 | null
None | ---- | 125.2 | 3.1
Predator family | 12, 0.8855 | 167.9 | 45.8

Figure 3.

#### Effect sizes and 95 % confidence intervals for consumptive effects of predators, for different categories of moderators (with number of studies in parentheses).

(a) predator family with predator class in the right-hand column, (b) vector genus, and (c) larval instar.

Figure 4.

#### Oviposition effect sizes and 95 % confidence intervals for different categories of vector genus (with number of studies in parentheses).

## Discussion

In this study, laboratory and semi-field empirical data were obtained through a systematic literature review and used to conduct a meta-analysis that assessed consumptive and non-consumptive effects of predators on mosquito prey. Some results agree with previously observed trends, such as greater consumptive effects from larger predators 74 Kumar et al.2007,12,31 104 Peters 1982,12,31 and no oviposition response to predator cues among container-breeding Aedes mosquitoes 146 Vonesh and Blaustein 2009,12,31. However, this meta-analysis revealed additional trends. Mosquito larval instar had an important role in moderating consumptive effects of predators, likely because of its direct correlation to prey size. Furthermore, a small, but significant, decrease in mosquito survival due to non-consumptive effects of predators was observed, suggesting that mosquitoes can be ‘scared to death’ by predators 106 Preisser et al.2004,12,31. Both male and female body sizes were also reduced among mosquitoes that had been exposed to predators, and predator avoidance during oviposition was observed among female Culex mosquitoes. Effects of predators on different vector traits, particularly survival, body size, and oviposition behavior, have the potential to influence infectious disease dynamics.

### Consumptive effects of predators on survival

Several larger predators reduced mosquito survival, including freshwater fish (Cyprinidae and Poeciliidae), salamander larvae (Ambystomatidae), dragonfly larvae (Aeshnidae), and backswimmers (Notonectidae) (Figure 3a). This finding is consistent with a previous analysis which showed a positive linear relationship between predator body mass and ingestion rate across taxa 104 Peters 1982,12,31. In addition, more effect size heterogeneity in the consumptive survival data was explained by an interaction between predator family and larval instar than was explained by an interaction between predator family and vector genus (Table 3). This result suggests that the relative sizes of predator and prey groups could play a more important role in determining consumptive mosquito survival than variations in predator responses to different behaviors of prey genera, which are likely to be shaped by the degree of shared evolutionary history between trophic levels 21 Buchanan et al.2016,12,31. Larval instar is an indicator of mosquito size, and previous modeling work has provided evidence of prey size selection by predators to maximize energetic gain 92 Mittelbach 1980,12,31. While smaller cyclopoid copepods are more effective against early instar mosquito larvae 34 Dieng et al.2001,12,31, larger predators including tadpoles, giant water bugs, dragonfly larvae, fish, and backswimmers are more effective against late instar larvae 75 Kweka et al.2010,12,31.

### Non-consumptive effects of predators on survival

Exposure to predation cues significantly lowered mosquito survival, and this non-consumptive effect has also been observed in dragonfly larvae prey (Leucorrhinia intacta) that were exposed to caged predators 90 McCauley et al.2010,12,31. The reduction in mosquito survival from non-consumptive effects of predators was significantly smaller than the reduction that was observed from consumptive effects. This is partially due to the practical constraints of most experimental designs, which cause consumptive and non-consumptive effects of predators on survival to be grouped together and reported as consumptive effects. The greater impact of combined consumptive and non-consumptive effects, in comparison to only non-consumptive effects, has previously been observed in pea aphids (Acyrthosiphon pisum) 96 Nelson et al.2003,12,31.

### Non-consumptive effects of predators on body size

While predators did not significantly impact mosquito development time through non-consumptive effects in either sex, mosquito body size was decreased by the non-consumptive effects of predators in both sexes. Smaller body size is associated with lower reproductive success in mosquitoes because smaller females lay fewer eggs 16 Blackmore and Lord 1999,12,31 84 Lyimo and Takken 1992,12,31 98 Oliver and Howard 2010,12,31 130 Styer et al.2006,12,31 138 Tsunoda et al.2009,12,31, and smaller males produce less sperm 55 Hatala et al.2017,12,31 105 Ponlawat and Harrington 2006,12,31. These effects suggest that predation could non-consumptively reduce mosquito population growth. The smaller size of mosquitoes exposed to predators could also limit disease transmission. Vector lifespan contributes disproportionately to disease transmission because older vectors are more likely to have been exposed to pathogens, more likely to already be infectious after having survived the extrinsic incubation period, and more likely to survive long enough to bite subsequent hosts 23 Cator et al.2019,12,31. It is well-established that smaller mosquito body size is associated with shorter mosquito lifespan 9 Araújo et al.2011,12,31 56 Hawley 1984,12,31 111 Reisen et al.1983,12,31 112 Reiskind and Lounibos 2008,12,31 149 Xue et al.2009,12,31. Therefore, non-consumptive effects of predators may limit the transmission of mosquito-borne diseases.

### Non-consumptive effects of predators on oviposition behavior

Predator presence also non-consumptively reduced oviposition behavior in adult female mosquitoes. Meta-regression results showed that Culex females significantly avoid oviposition sites that contain predators or predator cues, but Aedes and Culiseta females do not avoid these sites, despite a slight non-significant trend toward predator avoidance in Culiseta (Figure 4). Both Culex and Culiseta mosquitoes have an ‘all-or-none’ oviposition strategy 66 Johnson and Fonseca 2013,12,31, in which they lay hundreds of rapidly hatching eggs in rafts on the water’s surface 30 Day 2015,12,31. Such an oviposition strategy is conducive to evolving predator avoidance behaviors, and a previous meta-analysis showed significant predator avoidance in both Culex and Culiseta during oviposition 146 Vonesh and Blaustein 2009,12,31. Conversely, it is likely that an oviposition response to predation is not particularly advantageous for Aedes because the delayed hatching of their eggs 30 Day 2015,12,31 can prevent the level of predation risk at the time of oviposition from matching the level of predation risk present in the eventual larval environment 146 Vonesh and Blaustein 2009,12,31. The predator avoidance response in Aedes species that lay their eggs above the water’s edge in containers has previously been described as ‘non-existent’ 146 Vonesh and Blaustein 2009,12,31. Both Aedes species included in this study’s oviposition data subset, Ae. albopictus and Ae. aegypti, meet the criterion of ovipositing above water in containers 67 Juliano 2008,12,31. Predator avoidance during oviposition has previously been found to increase the mosquito population size at equilibrium 127 Spencer et al.2001,12,31. However, this study’s results and those of a previous meta-analysis 146 Vonesh and Blaustein 2009,12,31 suggest that models of oviposition site selection, such as those using parameters from Notonectidae predators and Culiseta prey 70 Kershenbaum et al.2011,12,31, are not generalizable to Aedes vectors.

### Implications for West Nile Virus disease dynamics

Predator avoidance during oviposition by Culex mosquitoes (Figure 4) may be of particular importance to West Nile virus (WNV) disease dynamics. Previous work has shown that Cx. pipiens, Cx. restuans, and Cx. tarsalis all avoid predator habitats 146 Vonesh and Blaustein 2009,12,31, and that Cx. pipiens is the primary bridge vector of WNV responsible for spill-over transmission from avian reservoir hosts to humans 44 Fonseca et al.2003,12,31 51 Hamer et al.2007,12,31 72 Kramer et al.2007,12,31 8 Andreadis 2011,12,31. Cx. pipiens mosquitoes can live in permanent aquatic environments, such as ground pools 3 Amini et al.2019,12,31 11 Barr 1966,12,31 33 Dida et al.2017,12,31 133 Sulesco et al.2014,12,31, ponds 83 Lühken et al.2014,12,31, stream edges 3 Amini et al.2019,12,31, and lake edges 145 Vinogradova 1999,12,31 that are more common in rural areas, but Cx. pipiens are also found in urban and suburban residential areas, where they typically breed in artificial containers 133 Sulesco et al.2014,12,31, including tires 83 Lühken et al.2014,12,31 97 Nikookar et al.2016,12,31 140 Verna 2014,12,31, rainwater tanks 136 Townroe and Callaghan 2013,12,31, and catch basins 45 Gardner et al.2011,12,31. Small artificial containers, such as discarded tires, are generally unlikely to harbor larger predators, including freshwater fish (Cyprinidae and Poeciliidae), salamander larvae (Ambystomatidae), dragonfly larvae (Aeshnidae), and backswimmers (Notonectidae), because temporary aquatic environments cannot support the relatively long development times of these organisms. The mean dispersal distance of adult Culex mosquitoes is greater than one kilometer 27 Ciota et al.2011,12,31 53 Hamer et al.2013,12,31, and female Cx. pipiens have exhibited longer dispersal distances after developing in the presence of a fish predator 2 Alcalay et al.2017,12,31. Therefore, predator avoidance during oviposition may cause Cx. pipiens populations to disperse from permanent aquatic environments in more rural areas to artificial container environments in urbanized areas, where the risk of human WNV infection is higher 19 Brown et al.2007,12,31.

Predator cue levels may be altered by climate conditions, and these changes in cue levels can impact WNV transmission to humans. Drought has previously been associated with human WNV cases 65 Johnson and Sukhdeo 2012,12,31 86 Marcantonio et al.2014,12,31 115 Roehr 2011,12,31 122 Shaman et al.2004,12,31 42 Epstein and Defilippo 2000,12,31 102 Paull et al.2016,12,31, but the association has thus far lacked a clear underlying mechanism. Under drought conditions, the density of aquatic organisms increases and predation pressures can intensify due to compressed space and high encounter rates 4 Amundrud et al.2018,12,31. A previous study of a stream ecosystem found that impacts of fish predation are more severe during the dry season 36 Dudgeon 1992,12,31. In addition, reductions in water volume can facilitate consumption of mosquito larvae by crane fly larvae (Tipulidae), whereas mosquito consumption by tipulids was not observed at a higher water level 4 Amundrud et al.2018,12,31. Laboratory and semi-field studies have shown that mosquitoes respond to a gradient of predator cues 118 Roux et al.2013,12,31 125 Silberbush and Blaustein 2010,12,31. The frequency of larval anti-predator behavior is correlated with the concentration of predator cues 118 Roux et al.2013,12,31, and adult female mosquitoes prefer oviposition sites with lower predator densities 125 Silberbush and Blaustein 2010,12,31. Therefore, as predator cue levels increase due to drought, permanent aquatic habitats are likely to transition from suitable oviposition sites for one generation of female mosquitoes, to unsuitable oviposition sites for the next generation.

When suitable oviposition sites are absent, females retain their eggs until sites become available 15 Bentley and Day 1988,12,31. Cx. pipiens females can retain their eggs for up to five weeks, allowing them enough time to find container sites with low predation risk, often located in residential areas 66 Johnson and Fonseca 2013,12,31. The movement of gravid female Cx. pipiens to residential areas increases the risk of WNV spill-over to humans because these vectors are likely to have already blood-fed at least once 28 Clements 1991,12,31, suggesting that they have a higher risk of WNV infection, relative to non-gravid mosquitoes. This is consistent with studies that have reported associations between drought and WNV-infected mosquitoes in urban and residential areas 65 Johnson and Sukhdeo 2012,12,31 102 Paull et al.2016,12,31. In addition, vertical transmission of WNV from gravid females to their progeny may occur during oviposition 117 Rosen 1987,12,31, when the virus is transmitted by an accessory gland fluid that attaches eggs to one another 95 Nelms et al.2012,12,31. Because the rate of vertical transmission in Cx. pipiens increases with the number of days following WNV infection 7 Anderson et al.2007,12,31, extended searches for oviposition sites due to drought could increase the frequency of vertical transmission. However, the impact of vertical transmission on WNV epidemics is thought to be minimal because when transmission to an egg raft did occur, only 4.7% of the progeny were found to be infected as adults 7 Anderson et al.2007,12,31, and only about half of those infected adults are estimated to be female. In summary, the movement of Cx. pipiens females toward more residential areas, combined with potential limited WNV amplification from increased vertical transmission, suggests that the vector trait of predator avoidance during oviposition can serve as a plausible explanation for associations between drought and human WNV cases.

Another theory for the association between drought and human WNV cases is based on the hypothesis that increased contact between mosquito vectors and passerine reservoir hosts occurs during drought conditions 102 Paull et al.2016,12,31 122 Shaman et al.2004,12,31. The proposed aggregation of bird and mosquito populations during drought was originally thought to occur in humid, densely vegetated hammocks – a type of habitat that is specific to southern Florida 122 Shaman et al.2004,12,31, but WNV incidence is more consistently clustered in other regions of the US, particularly the Northern Great Plains 24 CENTERS FOR DISEASE CONTROL AND PREVENTION 2020,12,31 132 Sugumaran et al.2008,12,31. Northern cardinals (Cardinalis cardinalis), American robins (Turdus migratorius), and house sparrows (Passer domesticus) were among the bird species that most frequently tested seropositive for WNV antibodies in 2005 and 2006 in Chicago, where high numbers of human cases were reported 52 Hamer et al.2007,12,31, and these passerine species are more abundant in residential areas, regardless of precipitation patterns 6 Anderson 2005,12,31 13 Beddall 1962,12,31 79 Lepczyk et al.2007,12,31. Apart from drought, landowners’ participation in supplemental bird feeding, providing bird houses, gardening, and maintaining vegetation can strongly influence passerine abundance in residential areas 78 Lepczyk et al.2003,12,31. Furthermore, as terrestrial foragers that can obtain hydration from their diet of insects, fruits, and other plant material 5 Anderson 2005,12,31 20 Brzek et al.2008,12,31 85 Malmborg and Willson 1987,12,31 113 Renne et al.1999,12,31, passerine reservoir hosts of WNV are less likely to move in response to drought than the mosquito vectors of WNV, which have obligate aquatic life stages.

While hatch-year birds are more vulnerable to mosquito biting, and thus contribute to the amplification of WNV 52 Hamer et al.2007,12,31, it is illogical to expect an increased abundance of hatch-year birds during drought conditions. However, some have argued that in cases where drought decreases the abundance of juvenile birds, the ratio of mosquitoes to birds increases, and this could lead to higher WNV prevalence in the mosquito population 102 Paull et al.2016,12,31. Although reductions in both hatching success 46 George et al.1991,12,31 and survival of recently fledged birds 150 Yackel Adams et al.2005,12,31 have been observed during drought conditions, the impact of drought on avian abundance varies widely by species 141 Verner and Purcell 1998,12,31. In particular, synanthropic species, such as those likely to harbor WNV, are less negatively affected by drought 1 Albright et al.2008,12,31. Additionally, the droughts that impact avian abundance often occur over much longer periods of time than the seasonal droughts that predict WNV transmission to humans. For example, avian abundance has been modeled based on precipitation metrics spanning 32 weeks, and house wren (Troglodytes aedon) abundance has been predicted by precipitation averages spanning four years 141 Verner and Purcell 1998,12,31. Finally, birds with higher levels of stress hormones are more likely to be fed on by mosquitoes, and certain factors associated with residential areas, such as road noise, light pollution, and pesticide exposure, can cause avian stress 47 Gervasi et al.2015,12,31. Therefore, elevated avian stress hormones in these habitats may contribute to WNV prevalence in the mosquito population, independent of drought conditions.

### Implications for mosquito-borne disease modeling

Although the aquatic phase of the mosquito life cycle is often overlooked in mathematical models of mosquito-borne pathogen transmission 110 Reiner et al.2012,12,31, vector survival at immature stages plays an important role in determining mosquito population abundance, which is an essential factor for predicting disease transmission 12 Beck-Johnson et al.2012,12,31. The results of this study show that mosquito survival decreases among the Aedes, Anopheles, and Culex genera due to consumptive effects of predators (Figure 3b), and that there is also a reduction in mosquito survival due to non-consumptive effects. Other studies have demonstrated that aquatic predators dramatically impact mosquito survival and abundance. For example, a biocontrol intervention relying on the application of copepod predators eliminated Aedes albopictus from three communes in Nam Dinh, Vietnam, where dengue transmission was previously detected, and reduced vector abundance by 86–98% in three other communes 69 Kay et al.2001,12,31. Conversely, the annual abundance of Culex and Anopheles mosquitoes was observed to increase 15-fold in semi-permanent wetlands in the year following a drought, likely because the drought eliminated aquatic predators from wetlands that dried completely, and mosquitoes were able to re-colonize newly formed aquatic habitats more quickly than their most effective predators 26 Chase and Knight 2002,12,31.

While relationships between temperature and different vector traits, such as fecundity and lifespan, have been incorporated into models of temperature effects on mosquito population density 40 El Moustaid and Johnson 2018,12,31, models of predator effects on vector borne disease transmission have focused primarily on the impacts of predation on vector survival. Previous models have shown that predators of vector species can decrease or eliminate pathogen infection in host populations as vector fecundity increases 94 Moore et al.2009,12,31. The findings of this meta-analysis suggest that predators also decrease vector fecundity through non-consumptive effects on vector body size. In addition, the entomological inoculation rate (EIR) is likely to be reduced by effects of predators on mosquito fecundity and lifespan, as well as effects of predators on mosquito survival. The EIR has been defined as the product of three variables: (m) the number of mosquitoes per host, (a) the daily rate of mosquito biting, and (s) the proportion of mosquitoes that are infectious 12 Beck-Johnson et al.2012,12,31. Based on this study’s findings, predators are likely to decrease the number of mosquitoes per host by reducing mosquito survival through both consumptive and non-consumptive effects, and by reducing mosquito fecundity through non-consumptive effects on body size. In addition, predators are likely to decrease the proportion of mosquitoes that are infectious by shortening the vector lifespan through non-consumptive effects on body size. The relationship between mosquito body size and biting rate is unclear, with some studies showing higher biting rates among larger mosquitoes 9 Araújo et al.2011,12,31 50 Gunathilaka et al.2018,12,31, and others reporting higher biting rates among smaller mosquitoes 43 Farjana and Tuno 2012,12,31 77 Leisnham et al.2007,12,31. The links between factors that influence the EIR and observed effects of predators on mosquito prey demonstrate the necessity of including both consumptive and non-consumptive effects of predators in models of mosquito-borne disease.

### Conclusion

This meta-analysis on mosquito predation demonstrates that predators not only play an important role in directly reducing mosquito populations, but also have non-consumptive effects on surviving mosquitoes that may ultimately reduce further population growth and decrease disease transmission. While families of larger sized predators were effective in reducing mosquito survival, other factors, such as impacts on native species, as well as the economic cost of mass-rearing and field applications 73 Kumar and Hwang 2005,12,31 107 Pyke 2007,12,31, should be carefully considered before selecting a predator as a suitable biocontrol agent. Predictive disease models are likely to be more reliable when the non-consumptive effects of predation are incorporated. Although exposure of mosquito larvae to predators is commonplace in outdoor field settings, it remains rare in most laboratory-based assessments of vector traits. Therefore, mosquitoes observed in nature are likely to have smaller body sizes than those observed under optimal laboratory conditions. It is important for disease modelers to recognize these impacts of predation on vector traits as they can reduce mosquito population growth and limit disease transmission due to shorter vector lifespans. Within the WNV disease system, consideration of the oviposition behavioral response to predation cues by Culex vectors can improve current understanding of the association between drought and human cases. This study provides general estimates of the effects of predators on selected mosquito traits for use in predictive disease models.

### Future directions

Modeling efforts that aim to optimize the application of biocontrol predators should also consider incorporating predator effects on vector survival, fecundity, and lifespan. These additions to predictive models of various biocontrol interventions are likely to help public health officials choose the most cost-effective strategies for limiting disease transmission. In the 60-study database that was compiled, only one study was designed to directly measure the effect of larval-stage predation on vector competence 119 Roux et al.2014,12,31. Therefore, future efforts to assess the impact of predators on mosquito-borne disease transmission should prioritize experimental studies in which infected mosquito larvae are observed throughout an initial period of aquatic exposure to predators, followed by a period of blood-feeding in the adult stage.

Two studies from the compiled database examined the compatibility of predators with Bacillus thuringiensis var. israelensis (Bti), a commonly used bacterial biocontrol agent 25 Chansang et al.2003,12,31 99 Beeck et al.2015,12,31. Previous studies have supported the simultaneous application of cyclopoid copepod predators and Bti null 88 Marten et al.1992,12,31 135 Tietze et al.1993,12,31, but additional analyses are needed on the use of Bti with other families of mosquito predators. Populations of other insect pests, such as the southern green stink bug (Nezara viridula), are known to be regulated by both predators and parasites 39 Ehler 2001,12,31. The literature search conducted for this meta-analysis returned studies on water mite parasites 109 Rajendran and Prasad 1993,12,31 and nematode parasitoids 31 Valdez 2005,12,31 of mosquitoes, and ascogregarine parasites have previously been evaluated as biocontrol agents against Aedes mosquitoes 137 Tseng 2006,12,31. A more thorough review of the impacts of parasites and parasitoids on vector traits, such as survival, fecundity, and lifespan, is needed before incorporating these potential biocontrol agents into integrated vector control plans.

Three studies in the 60-study database included experiments where two mosquito prey species were made available to the predator species 48 Grill and Juliano 1995,12,31 49 Griswold and Lounibos 2004,12,31 91 Micieli et al.2001,12,31. In these cases, the effect size measurement for each mosquito species could be influenced by interspecific competition, or a preference of the predator species for a certain prey species. Hetero-specific prey observations were excluded from this meta-analysis, but future analyses centered on the concepts of interspecific competition or predator preferences might further evaluate these data. In addition, this meta-analysis investigated consumptive and non-consumptive effects of predators separately. More research is needed to determine how models should combine these different types of predator effects to accurately reflect predation interactions as they occur in natural environments.

## References

1. 1. TP Albright
   2. AM Pidgeon
   3. CD Rittenhouse
   4. MK Clayton
   5. CH Flather
   6. PD Culbert
   7. BD Wardlow
   8. VC Radeloff
   2008 Effects of drought on avian community structure 16 Global Change Biology 2158 2170
2. 1. Y Alcalay
   2. I Tsurim
   3. O Ovadia
   2017 Female mosquitoes disperse further when they develop under predation risk 29 Behavioral Ecology 1402 1408
3. 1. M Amini
   2. AA Hanafi-Bojd
   3. AA Aghapour
   4. AR Chavshin
   2019 Larval habitats and species diversity of mosquitoes (Diptera: Culicidae) in West Azerbaijan Province, Northwestern Iran 20 BMC ecology
4. 1. SL Amundrud
   2. SA Clay-Smith
   3. BL Flynn
   4. KE Higgins
   5. MS Reich
   6. DRH Wiens
   7. DS Srivastava
   2018 Drought alters the trophic role of an opportunistic generalist in an aquatic ecosystem 189 Oecologia 733 744
5. 1. TR Anderson
   2005 Chapter 6: Foraging Behavior and Food. Biology of the Ubiquitous House Sparrow: From Genes to Populations Oxford University Press
6. 1. TR Anderson
   2005 Chapter 10: Human Commensalism and Pest Management. Biology of the Ubiquitous House Sparrow: From Genes to Populations Oxford University Press
7. 1. JF Anderson
   2. AJ Main
   3. K Delroux
   4. E Fikrig
   2007 Extrinsic incubation periods for horizontal and vertical transmission of West Nile virus by Culex pipiens pipiens (Diptera: Culicidae)45 Journal of medical entomology 445 451
8. 1. TG Andreadis
   2011 The contribution of Culex pipiens complex mosquitoes to transmission and persistence of West Nile virus in North America 28 Journal of the American Mosquito Control Association 137 151
9. 1. M Araújo
   2. LHS Gil
   3. A Silva
   2011 Larval food quantity affects development time, survival and adult biological traits that influence the vectorial capacity of Anopheles darlingi under laboratory conditions 11 Malaria Journal
10. 1. D Arav
    2. L Blaustein
    2005 Effects of pool depth and risk of predation on oviposition habitat selection by temporary pool dipterans 43 Journal of Medical Entomology 493 497
11. 1. AR Barr
    1966 Occurrence and distribution of the Culex pipiens complex 37 Bulletin of the World Health Organization 293 296
12. 1. LM Beck-Johnson
    2. WA Nelson
    3. KP Paaijmans
    4. AF Read
    5. MB Thomas
    6. ON Bjørnstad
    2012 The effect of temperature on Anopheles mosquito population dynamics and the potential for malaria transmission 8 PLOS ONE
13. 1. BG Beddall
    1962 Range expansion of the cardinal and other birds in the northeastern states 75 The Wilson Bulletin 140 158
14. 1. G Benelli
    2. CL Jeffries
    3. T Walker
    2015 Biological Control of Mosquito Vectors: Past, Present, and Future 7 Insects
15. 1. MD Bentley
    2. JF Day
    1988 Chemical ecology and behavioral aspects of mosquito oviposition 34 Annual review of entomology 401 421
16. 1. MS Blackmore
    2. CC Lord
    1999 The relationship between size and fecundity in Aedes albopictus 25 Journal of Vector Ecology 212 217
17. 1. JG Bond
    2. JI Arredondo-Jimenez
    3. MH Rodriguez
    4. H Quiroz-Martinez
    5. T Williams
    2004 Oviposition habitat selection for a predator refuge and food source in a mosquito 30 Ecological Entomology 255 263
18. 1. U Brose
    2. T Jonsson
    3. EL Berlow
    4. P Warren
    5. C Banasek-Richter
    6. L-F Bersier
    7. JL Blanchard
    8. T Brey
    9. SR Carpenter
    10. M-FC Blandenier
    11. L Cushing
    12. HA Dawah
    13. T Dell
    14. F Edwards
    15. S Harper-Smith
    16. U Jacob
    17. ME Ledger
    18. ND Martinez
    19. J Memmott
    20. K Mintenbeck
    21. JK Pinnegar
    22. BC Rall
    23. TS Rayner
    24. DC Reuman
    25. L Ruess
    26. W Ulrich
    27. RJ Williams
    28. G Woodward
    29. JE Cohen
    2005 Consumer-resource body-size relationships in natural food webs 87 Ecology 2411 2417
19. 1. HE Brown
    2. JE Childs
    3. MA Diuk-Wasser
    4. D Fish
    2007 Ecological factors associated with West Nile virus transmission, northeastern United States 14 Emerging Infectious Diseases 1539 1545
20. 1. P Brzek
    2. K Kohl
    3. E Caviedes-Vidal
    4. WH Karasov
    2008 Developmental adjustments of house sparrow (Passer domesticus) nestlings to diet composition 212 The Journal of experimental biology 1284 1293
21. 1. AL Buchanan
    2. SL Hermann
    3. M Lund
    4. Z Szendrei
    2016 A meta-analysis of non-consumptive predator effects in arthropods: the influence of organismal and environmental characteristics 126 Oikos
22. 1. KP Burnham
    2. DR Anderson
    2003 Multimodel Inference: Understanding AIC and BIC in Model Selection 33 Sociological Methods & Research 261 304
23. 1. LJ Cator
    2. LR Johnson
    3. EA Mordecai
    4. FE Moustaid
    5. TRC Smallwood
    6. SL LaDeau
    7. MA Johansson
    8. PJ Hudson
    9. M Boots
    10. MB Thomas
    11. AG Power
    12. S Pawar
    2019 The Role of Vector Trait Variation in Vector-Borne Disease Dynamics 8 Frontiers in ecology and evolution
24. 1. CENTERS FOR DISEASE CONTROL AND PREVENTION
    2020 West Nile virus: Final Cumulative Maps & Data for 1999-2019 https://www.cdc.gov/westnile/statsmaps/cumMapsData.html
25. 1. UR Chansang
    2. A Bhumiratana
    3. P Kittayapong
    2003 Combination of Mesocyclops thermocyclopoides and Bacillus thuringiensis var. israelensis: a better approach for the control of Aedes aegypti larvae in water containers 29 Journal of Vector Ecology 218 226
26. 1. JM Chase
    2. TM Knight
    2002 Drought-induced mosquito outbreaks in wetlands 6 Ecology Letters 1017 1024
27. 1. AT Ciota
    2. CL Drummond
    3. MA Ruby
    4. J Drobnack
    5. GD Ebel
    6. LD Kramer
    2011 Dispersal of Culex mosquitoes (Diptera: Culicidae) from a wastewater treatment facility 49 Journal of medical entomology 35 42
28. 1. AN Clements
    1991 The Biology of Mosquitoes: Development, Nutrition and Reproduction Chapman & Hall
29. 1. HG Dam
    2. WT Peterson
    1987 The effect of temperature on the gut clearance rate constant of planktonic copepods 10 Journal of Experimental Marine Biology and Ecology
30. 1. JF Day
    2015 Mosquito Oviposition Behavior and Vector Control 7 Insects
31. 1. MRW Valdez
    2005 Parasitoid-induced behavioral alterations of Aedes aegypti mosquito larvae infected with mermithid nematodes (Nematoda: Mermithidae)31 Journal of vector ecology 344 354
32. 1. V Delnat
    2. TT Tran
    3. L Janssens
    4. R Stoks
    2018 Resistance to a chemical pesticide increases vulnerability to a biopesticide: Effects on direct mortality and mortality by predation 216 Aquatic toxicology
33. 1. GO Dida
    2. DN Anyona
    3. PO Abuom
    4. D Akoko
    5. SO Adoka
    6. A-S Matano
    7. PO Owuor
    8. C Ouma
    2017 Spatial distribution and habitat characterization of mosquito species during the dry season along the Mara River and its tributaries, in Kenya and Tanzania 7 Infectious diseases of poverty
34. 1. H Dieng
    2. M Boots
    3. N Tuno
    4. Y Tsuda
    5. M Takagi
    2001 A laboratory and field evaluation of Macrocyclops distinctus, Megacyclops viridis and Mesocyclops pehpeiensis as control agents of the dengue vector Aedes albopictus in a peridomestic area in Nagasaki, Japan 16 Medical and veterinary entomology 285 291
35. 1. H Dieng
    2. M Boots
    3. N Tuno
    4. Y Tsuda
    5. M Takagi
    2002 Life history effects of prey choice by copepods: implications for biocontrol of vector mosquitoes 19 Journal of the American Mosquito Control Association 67 73
36. 1. D Dudgeon
    1992 The effects of spate-induced disturbance, predation and environmental complexity on macroinvertebrates in a tropical stream 30 Freshwater Biology 189 197
37. 1. S Duval
    2. R Tweedie
    1999 Trim and fill: A simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis 56 Biometrics 455 463
38. 1. S Duval
    2. R Tweedie
    1999 A Nonparametric “Trim and Fill” Method of Accounting for Publication Bias in Meta-Analysis 95 Journal of the American Statistical Association 89 98
39. 1. LE Ehler
    2001 An evaluation of some natural enemies of >Nezara viridula in northern California 47 BioControl 309 325
40. 1. F El Moustaid
    2. LR Johnson
    2018 Modeling Temperature Effects on Population Density of the Dengue Mosquito Aedes aegypti 10 Insects
41. 1. G Englund
    2. G Ohlund
    3. CL Hein
    4. S Diehl
    2010 Temperature dependence of the functional response 14 Ecology letters 914 921
42. 1. PR Epstein
    2. C Defilippo
    2000 West Nile Virus and Drought 2 Global Change and Human Health 105 107
43. 1. T Farjana
    2. N Tuno
    2012 Multiple blood feeding and host-seeking behavior in Aedes aegypti and Aedes albopictus (Diptera: Culicidae)50 Journal of medical entomology 838 846
44. 1. DM Fonseca
    2. N Keyghobadi
    3. CA Malcolm
    4. C Mehmet
    5. F Schaffner
    6. M Mogi
    7. RC Fleischer
    8. RC Wilkerson
    2003 Emerging vectors in the Culex pipiens complex 303 Science 1535 1538
45. 1. AM Gardner
    2. GL Hamer
    3. AM Hines
    4. CM Newman
    5. ED Walker
    6. MO Ruiz
    2011 Weather variability affects abundance of larval Culex (Diptera: Culicidae) in storm water catch basins in suburban Chicago 49 Journal of medical entomology 270 276
46. 1. TL George
    2. AC Fowler
    3. RL Knight
    4. LC McEwen
    1991 Impacts of a Severe Drought on Grassland Birds in Western North Dakota 2 Ecological applications 275 284
47. 1. SS Gervasi
    2. N Burkett-Cadena
    3. SC Burgan
    4. AW Schrey
    5. HK Hassan
    6. TR Unnasch
    7. LB Martin
    2015 Host stress hormones alter vector feeding preferences, success, and productivity 283 Proceedings. Biological Sciences
48. 1. CP Grill
    2. SA Juliano
    1995 Predicting Species Interactions Based on Behaviour: Predation and Competition in Container-Dwelling Mosquitoes 65 The Journal of Animal Ecology
49. 1. MW Griswold
    2. LP Lounibos
    2004 Does differential predation permit invasive and native mosquito larvae to coexist in Florida?30 Ecological Entomology 122 127
50. 1. N Gunathilaka
    2. H Upulika
    3. L Udayanga
    4. D Amarasinghe
    2018 Effect of Larval Nutritional Regimes on Morphometry and Vectorial Capacity of Aedes aegypti for Dengue Transmission 2019 BioMed Research International 1 11
51. 1. GL Hamer
    2. UD Kitron
    3. JD Brawn
    4. SR Loss
    5. MO Ruiz
    6. TL Goldberg
    7. ED Walker
    2007 Culex pipiens (Diptera: Culicidae): a bridge vector of West Nile virus to humans 45 Journal of Medical Entomology 125 128
52. 1. GL Hamer
    2. ED Walker
    3. JD Brawn
    4. SR Loss
    5. MO Ruiz
    6. TL Goldberg
    7. AM Schotthoefer
    8. WM Brown
    9. E Wheeler
    10. UD Kitron
    2007 Rapid amplification of West Nile virus: the role of hatch-year birds 8 Vector Borne and Zoonotic Diseases 57 67
53. 1. GL Hamer
    2. TK Anderson
    3. DJ Donovan
    4. JD Brawn
    5. BL Krebs
    6. AM Gardner
    7. MO Ruiz
    8. WM Brown
    9. UD Kitron
    10. CM Newman
    11. TL Goldberg
    12. ED Walker
    2013 Dispersal of adult culex mosquitoes in an urban west nile virus hotspot: a mark-capture study incorporating stable isotope enrichment of natural larval habitats 8 PLOS Neglected Tropical Diseases
54. 1. PA Hancock
    2. A Wiebe
    3. KA Gleave
    4. S Bhatt
    5. E Cameron
    6. A Trett
    7. D Weetman
    8. DL Smith
    9. J Hemingway
    10. M Coleman
    11. PW Gething
    12. CL Moyes
    2017 Associated patterns of insecticide resistance in field populations of malaria vectors across Africa 115 PNAS 5938 5943
55. 1. AJ Hatala
    2. LC Harrington
    3. EC Degner
    2017 Age and Body Size Influence Sperm Quantity in Male Aedes albopictus (Diptera: Culicidae) Mosquitoes 55 Journal of medical entomology 1051 1054
56. 1. WA Hawley
    1984 The Effect of Larval Density on Adult Longevity of a Mosquito, Aedes sierrensis: Epidemiological Consequences 54 The Journal of Animal Ecology
57. 1. J Hemingway
    2. H Ranson
    1999 Insecticide resistance in insect vectors of human disease 45 Annual review of entomology 371 391
58. 1. SL Hermann
    2. JS Thaler
    2013 Prey perception of predation risk: volatile chemical cues mediate non-consumptive effects of a predator on a herbivorous insect 176 Oecologia 669 676
59. 1. SL Hermann
    2. DA Landis
    2016 Scaling up our understanding of non-consumptive effects in insect systems 20 Current opinion in insect science 54 60
60. 1. JPT Higgins
    2. SG Thompson
    3. JJ Deeks
    4. DG Altman
    2002 Measuring inconsistency in meta-analyses 327 BMJ 557 560
61. 1. JPT Higgins
    2. J Thomas
    3. J Chandler
    4. M Cumpston
    5. T Li
    6. MJ Page
    7. VA Welch
    2019 Cochrane Handbook for Systematic Reviews of Interventions John Wiley & Sons, Inc
62. 1. CA Hill
    2. FC Kafatos
    3. SK Stansfield
    4. FH Collins
    2004 Arthropod-borne diseases: vector control in the genomics era 3 Nature reviews. Microbiology 262 268
63. 1. L Janssens
    2. R Stoks
    2011 How does a pesticide pulse increase vulnerability to predation? Combined effects on behavioral antipredator traits and escape swimming 110 Aquatic toxicology 91 98
64. 1. MD Jennions
    2. AP Møller
    2001 Publication bias in ecology and evolution: an empirical assessment using the “trim and fill” method 77 Biological Reviews of the Cambridge Philosophical Society 211 222
65. 1. BJ Johnson
    2. MVK Sukhdeo
    2012 Drought-induced amplification of local and regional West Nile virus infection rates in New Jersey 50 Journal of Medical Entomology 195 204
66. 1. B Johnson
    2. DM Fonseca
    2013 The effects of forced-egg retention on the blood-feeding behavior and reproductive potential of Culex pipiens (Diptera: Culicidae)66 Journal of insect physiology 53 58
67. 1. SA Juliano
    2008 Species interactions among larval mosquitoes: context dependence across habitat gradients 54 Annual review of entomology 37 56
68. 1. RM Kalinoski
    2. JP DeLong
    2015 Beyond body mass: how prey traits improve predictions of functional response parameters 180 Oecologia 543 550
69. 1. BH Kay
    2. VS Nam
    3. TV Tien
    4. NT Yen
    5. TV Phong
    6. VTB Diep
    7. TU Ninh
    8. A Bektas
    9. JG Aaskov
    2001 Control of aedes vectors of dengue in three provinces of Vietnam by use of Mesocyclops (Copepoda) and community-based methods validated by entomologic, clinical, and serological surveillance 66 The American journal of tropical medicine and hygiene 40 48
70. 1. A Kershenbaum
    2. M Spencer
    3. L Blaustein
    4. JE Cohen
    2011 Modelling evolutionarily stable strategies in oviposition site selection, with varying risks of predation and intraspecific competition 26 Evolutionary Ecology 955 974
71. 1. CA Kleckner
    2. WA Hawley
    3. WE Bradshaw
    4. CM Holzapfel
    5. IJ Fisher
    1994 Protandry in aedes sierrensis - the significance of temporal variation in female fecundityROTANDRY IN AEDES SIERRENSIS - THE SIGNIFICANCE OF TEMPORAL VARIATION IN FEMALE FECUNDITY 76 Ecology 1242 1250
72. 1. LD Kramer
    2. LM Styer
    3. GD Ebel
    2007 A global perspective on the epidemiology of West Nile virus 53 Annual review of entomology 61 81
73. 1. R Kumar
    2. JS Hwang
    2005 Larvicidal efficiency of aquatic predators: A perspective for mosquito biocontrol 45 Zoological Studies 447 466
74. 1. R Kumar
    2. P Muhid
    3. HU Dahms
    4. LC Tseng
    5. JS Hwang
    2007 Potential of three aquatic predators to control mosquitoes in the presence of alternative prey: a comparative experimental assessment 59 Marine and Freshwater Research
75. 1. EJ Kweka
    2. G Zhou
    3. TM Gilbreath
    4. Y Afrane
    5. M Nyindo
    6. AK Githeko
    7. G Yan
    2010 Predation efficiency of Anopheles gambiae larvae by aquatic predators in western Kenya highlands 4 Parasites & Vectors
76. 1. MJ Lajeunesse
    2015 Facilitating systematic reviews, data extraction and metaanalysis with the metagear package for R 7 Methods in Ecology and Evolution 323 330
77. 1. PT Leisnham
    2. LM Sala
    3. SA Juliano
    2007 Geographic variation in adult survival and reproductive tactics of the mosquito Aedes albopictus 45 Journal of Medical Entomology 210 221
78. 1. CA Lepczyk
    2. AG Mertig
    3. J Liu
    2003 Assessing landowner activities related to birds across rural-to-urban landscapes 33 Environmental management 110 125
79. 1. CA Lepczyk
    2. CH Flather
    3. VC Radeloff
    4. AM Pidgeon
    5. RB Hammer
    6. J Liu
    2007 Human impacts on regional avian diversity and abundance 22 Conservation biology 405 416
80. 1. SL Lima
    2. LM Dill
    1989 Behavioral decisions made under the risk of predation: a review and prospectus 68 Canadian Journal of Zoology 619 640
81. 1. SL Lima
    1997 Nonlethal Effects in the Ecology of Predator-Prey Interactions: What are the ecological effects of anti-predator decision-making 48 Bioscience 25 34
82. 1. N Liu
    2014 Insecticide resistance in mosquitoes: impact, mechanisms, and research directions 60 Annual Review of Entomology 537 559
83. 1. R Lühken
    2. S Steinke
    3. M Leggewie
    4. E Tannich
    5. A Krüger
    6. S Becker
    7. E Kiel
    2014 Physico-Chemical Characteristics of Culex pipiens sensu lato and Culex torrentium (Diptera: Culicidae) Breeding Sites in Germany 52 Journal of medical entomology 932 936
84. 1. EO Lyimo
    2. W Takken
    1992 Effects of adult body size on fecundity and the pre-gravid rate of Anopheles gambiae females in Tanzania 7 Medical and veterinary entomology 328 332
85. 1. PK Malmborg
    2. MF Willson
    1987 Foraging ecology of avian frugivores and some consequences for seed dispersal in an illinois woodlot 90 The Condor 173 186
86. 1. M Marcantonio
    2. A Rizzoli
    3. M Metz
    4. R Rosà
    5. G Marini
    6. E Chadwick
    7. M Neteler
    2014 Identifying the environmental conditions favouring West Nile Virus outbreaks in Europe 10 PLOS ONE
87. 1. GG Marten
    1989 Elimination of Aedes albopictus from tire piles by introducing Macrocyclops albidus (Copepoda, Cyclopidae)6 Journal of the American Mosquito Control Association 689 693
88. 1. GG Marten
    2. W Che
    3. ES Bordes
    1992 Compatibility of cyclopoid copepods with mosquito insecticides 9 Journal of the American Mosquito Control Association 150 154
89. 1. GG Marten
    2. JW Reid
    2006 Cyclopoid copepods 23 Journal of the American Mosquito Control Association 65 92
90. 1. SJ McCauley
    2. L Rowe
    3. M-J Fortin
    2010 The deadly effects of “nonlethal” predators 92 Ecology 2043 2048
91. 1. MV Micieli
    2. G Marti
    3. JJ García
    2001 Laboratory evaluation of Mesocyclops annulatus (Wierzejski, 1892) (Copepoda: Cyclopidea) as a predator of container-breeding mosquitoes in Argentina 97 Memorias do Instituto Oswaldo Cruz 835 838
92. 1. GG Mittelbach
    1980 Foraging Efficiency and Body Size: A Study of Optimal Diet and Habitat Use by Bluegills 62 Ecology 1370 1386
93. 1. D Moher
    2. A Liberati
    3. J Tetzlaff
    4. DG Altman
    5. PRISMA Group
    2008 Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement 339 BMJ
94. 1. SM Moore
    2. ET Borer
    3. PR Hosseini
    2009 Predators indirectly control vector-borne disease: linking predator-prey and host-pathogen models 7 Journal of the Royal Society, Interface 161 176
95. 1. BM Nelms
    2. E Fechter-Leggett
    3. BD Carroll
    4. P Macedo
    5. S Kluh
    6. WK Reisen
    2012 Experimental and natural vertical transmission of West Nile virus by California Culex (Diptera: Culicidae) mosquitoes 50 Journal of medical entomology 371 378
96. 1. EH Nelson
    2. CE Matthews
    3. JA Rosenheim
    2003 Predators reduce prey population growth by inducing changes in prey behaviorREDATORS REDUCE PREY POPULATION GROWTH BY INDUCING CHANGES IN PREY BEHAVIOR 85 Ecology 1853 1858
97. 1. SH Nikookar
    2. M Fazeli-Dinan
    3. S Azari-Hamidian
    4. SN Mousavinasab
    5. M Arabi
    6. SP Ziapour
    7. J Shojaee
    8. A Enayati
    2016 Species composition and abundance of mosquito larvae in relation with their habitat characteristics in Mazandaran Province, northern Iran 107 Bulletin of entomological research 598 610
98. 1. J Oliver
    2. JJ Howard
    2010 Fecundity of wild-caught gravid Culiseta morsitans (Diptera: Culicidae)48 Journal of medical entomology 196 201
99. 1. L Beeck
    2. L Janssens
    3. R Stoks
    2015 Synthetic predator cues impair immune function and make the biological pesticide Bti more lethal for vector mosquitoes 26 Ecological Applications 355 366
100. 1. GD Ower
     2. SA Juliano
     2018 The demographic and life-history costs of fear: Trait-mediated effects of threat of predation on Aedes triseriatus 9 Ecology and evolution 3794 3806
101. 1. GF O’Meara
     1978 Variable expressions of autogeny in three mosquito species 1 International Journal of Invertebrate Reproduction 253 261
102. 1. SH Paull
     2. DE Horton
     3. M Ashfaq
     4. D Rastogi
     5. LD Kramer
     6. NS Diffenbaugh
     7. AM Kilpatrick
     2016 Drought and immunity determine the intensity of West Nile virus epidemics and climate change impacts 284 Proceedings. Biological Sciences
103. 1. SD Peacor
     2. EE Werner
     2000 The contribution of trait-mediated indirect effects to the net effects of a predator 98 PNAS 3904 3908
104. 1. RH Peters
     1982 The Ecological Implications of Body Size Cambridge University Press
105. 1. A Ponlawat
     2. LC Harrington
     2006 Age and body size influence male sperm capacity of the dengue vector Aedes aegypti (Diptera: Culicidae)44 Journal of medical entomology 422 426
106. 1. EL Preisser
     2. DI Bolnick
     3. MF Benard
     2004 Scared to death? the effects of intimidation and consumption in predator–prey interactions 86 Ecology 501 509
107. 1. GH Pyke
     2007 Plague Minnow or Mosquito Fish? A Review of the Biology and Impacts of Introduced Gambusia Species 39 Annual Review of Ecology, Evolution, and Systematics 171 191
108. 1. R Development Core Team
     2019 R: A language and environment for statistical computing R Foundation for Statistical Computing
109. 1. R Rajendran
     2. RS Prasad
     1993 A laboratory study on the life cycle and feeding behaviour of Arrenurus madaraszi (Acari: Arrenuridae) parasitizing Anopheles mosquitoes 88 Annals of tropical medicine and parasitology 169 174
110. 1. RC Reiner
     2. TA Perkins
     3. CM Barker
     4. T Niu
     5. LF Chaves
     6. AM Ellis
     7. DB George
     8. A Le Menach
     9. JRC Pulliam
     10. D Bisanzio
     11. C Buckee
     12. C Chiyaka
     13. DAT Cummings
     14. AJ Garcia
     15. ML Gatton
     16. PW Gething
     17. DM Hartley
     18. G Johnston
     19. EY Klein
     20. E Michael
     21. SW Lindsay
     22. AL Lloyd
     23. DM Pigott
     24. WK Reisen
     25. N Ruktanonchai
     26. BK Singh
     27. AJ Tatem
     28. U Kitron
     29. SI Hay
     30. TW Scott
     31. DL Smith
     2012 A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970-2010 10 Journal of the Royal Society, Interface
111. 1. WK Reisen
     2. MM Milby
     3. ME Bock
     1983 The effects of immature stress on selected events in the life-history of culex-tarsalisHE EFFECTS OF IMMATURE STRESS ON SELECTED EVENTS IN THE LIFE-HISTORY OF CULEX-TARSALIS 44 Mosquito News 385 395
112. 1. MH Reiskind
     2. LP Lounibos
     2008 Effects of intraspecific larval competition on adult longevity in the mosquitoes Aedes aegypti and Aedes albopictus 23 Medical and veterinary entomology 62 68
113. 1. IJ Renne
     2. SA Gauthreaux
     3. CA Gresham
     1999 Seed Dispersal of the Chinese Tallow Tree (Sapium sebiferum (L.) Roxb.) by Birds in Coastal South Carolina 144 The American Midland Naturalist 202 215
114. 1. D Roberts
     2011 Responses of three species of mosquito larvae to the presence of predatory dragonfly and damselfly larvae 145 Entomologia Experimentalis et Applicata 23 29
115. 1. B Roehr
     2011 US hit by massive West Nile virus outbreak centred around Texas 345 British Medical Journal
116. 1. A Rohatgi
     2019 WebPlotDigitizer https://automeris.io/WebPlotDigitizer
117. 1. L Rosen
     1987 Further observations on the mechanism of vertical transmission of flaviviruses by Aedes mosquitoes 39 The American journal of tropical medicine and hygiene 123 126
118. 1. O Roux
     2. A Diabaté
     3. F Simard
     2013 Divergence in threat sensitivity among aquatic larvae of cryptic mosquito species 83 The Journal of animal ecology 702 711
119. 1. O Roux
     2. A Vantaux
     3. B Roche
     4. KB Yameogo
     5. KR Dabiré
     6. A Diabaté
     7. F Simard
     8. T Lefèvre
     2014 Evidence for carry-over effects of predator exposure on pathogen transmission potential 282 Proceedings. Biological Sciences
120. 1. BM Russell
     2. LE Muir
     3. P Weinstein
     4. BH Kay
     1995 Surveillance of the mosquito Aedes aegypti and its biocontrol with the copepod Mesocyclops aspericornis in Australian wells and gold mines 10 Medical and veterinary entomology 155 160
121. 1. A Seale
     1916 The Mosquito Fish, Gambusia affinis (Baird and Girard), in the Philippine Islands 12 Philippine Journal of Science 177 187
122. 1. J Shaman
     2. JF Day
     3. M Stieglitz
     2004 Drought-induced amplification and epidemic transmission of West Nile virus in southern Florida 42 Journal of medical entomology 134 141
123. 1. L Shi
     2. L Lin
     2018 The trim-and-fill method for publication bias: practical guidelines and recommendations based on a large database of meta-analyses 98 Medicine
124. 1. A Sih
     1985 Antipredator Responses and the Perception of Danger by Mosquito Larvae 67 Ecology 434 441
125. 1. A Silberbush
     2. L Blaustein
     2010 Mosquito females quantify risk of predation to their progeny when selecting an oviposition site 25 Functional Ecology 1091 1095
126. 1. MKF Soumare
     2. JE Cilek
     2010 The effectiveness of Mesocyclops longisetus (Copepoda) for the control of container-inhabiting mosquitoes in residential environments 27 Journal of the American Mosquito Control Association 376 383
127. 1. M Spencer
     2. L Blaustein
     3. JE Cohen
     2001 Oviposition habitat selection by mosquitoes (culiseta longiareolata) and consequences for population size 83 Ecology 669 679
128. 1. G Stav
     2. L Blaustein
     3. Y Margalit
     2004 Individual and interactive effects of a predator and controphic species on mosquito populations 15 Ecological Applications 587 598
129. 1. JA Sterne
     2. M Egger
     2000 Funnel plots for detecting bias in meta-analysis: guidelines on choice of axis 54 Journal of clinical epidemiology 1046 1055
130. 1. LM Styer
     2. MA Meola
     3. LD Kramer
     2006 West Nile virus infection decreases fecundity of Culex tarsalis females 44 Journal of medical entomology 1074 1085
131. 1. T Su
     2. MS Mulla
     1996 Physiological aspects of autogeny in Culex tarsalis (Diptera: Culicidae): influences of sugar-feeding, mating, body weight, and wing length 22 Journal of Vector Ecology 115 121
132. 1. R Sugumaran
     2. SR Larson
     3. JP Degroote
     2008 Spatio-temporal cluster analysis of county-based human West Nile virus incidence in the continental United States 8 International Journal of Health Geographics
133. 1. TM Sulesco
     2. LG Toderas
     3. IG Uspenskaia
     4. IK Toderas
     2014 Larval Habitats Diversity and Distribution of the Mosquito (Diptera: Culicidae) Species in the Republic of Moldova 52 Journal of medical entomology 1299 1308
134. 1. R Thapa
     2. SB Mirsky
     3. KL Tully
     2017 Cover Crops Reduce Nitrate Leaching in Agroecosystems:A Global Meta-Analysis 47 Journal of environmental quality 1400 1411
135. 1. NS Tietze
     2. PG Hester
     3. KR Shaffer
     4. SJ Prescott
     5. ET Schreiber
     1993 Integrated management of waste tire mosquitoes utilizing Mesocyclops longisetus (Copepoda: Cyclopidae), Bacillus thuringiensis var. israelensis, Bacillus sphaericus, and methoprene 10 Journal of the American Mosquito Control Association 363 373
136. 1. S Townroe
     2. A Callaghan
     2013 British container breeding mosquitoes: the impact of urbanisation and climate change on community composition and phenology 9 PLOS ONE
137. 1. M Tseng
     2006 Ascogregarine parasites as possible biocontrol agents of mosquitoes 23 Journal of the American Mosquito Control Association 30 34
138. 1. T Tsunoda
     2. A Fukuchi
     3. S Nanbara
     4. M Takagi
     2009 Effect of body size and sugar meals on oviposition of the yellow fever mosquito, Aedes aegypti (Diptera: Culicidae)35 Journal of vector ecology 56 60
139. 1. SF Uiterwaal
     2. JP Delong
     2019 Functional responses are maximized at intermediate temperatures 101 Ecology
140. 1. TN Verna
     2014 Species Composition and Seasonal Distribution of Mosquito Larvae (Diptera: Culicidae) in Southern New Jersey, Burlington County 52 Journal of medical entomology 1165 1169
141. 1. J Verner
     2. KL Purcell
     1998 Fluctuating Populations of House Wrens and Bewick’s Wrens in Foothills of the Western Sierra Nevada of California 101 The Condor 219 229
142. 1. R Veronesi
     2. M Carrieri
     3. B Maccagnani
     4. S Maini
     5. R Bellini
     2014 Macrocyclops albidus (Copepoda: cyclopidae) for the Biocontrol of Aedes albopictus and Culex pipiens in Italy 31 Journal of the American Mosquito Control Association 32 43
143. 1. W Viechtbauer
     2009 Conducting Meta-Analyses in R with the metafor Package 36 Journal of Statistical Software 1 48
144. 1. E Vincze
     2. G Seress
     3. M Lagisz
     4. N Dinemanse
     5. P Sprau
     2016 Does Urbanization Affect Predation of Bird Nests? A Meta-Analysis 5 Frontiers in Ecology and Evolution
145. 1. EB Vinogradova
     1999 Culex Pipiens Pipiens Mosquitoes: Taxonomy, Distribution, Ecology, Physiology, Genetics, Applied Importance and Control Pensoft
146. 1. JR Vonesh
     2. L Blaustein
     2009 Predator-Induced Shifts in Mosquito Oviposition Site Selection: A Meta-Analysis and Implications for Vector Control 56 Israel Journal of Ecology and Evolution 263 279
147. 1. SC Weaver
     2. WK Reisen
     2009 Present and future arboviral threats 85 Antiviral research 328 345
148. 1. WORLD HEALTH ORGANIZATION
     2019 Vector-borne diseases https://www.who.int/news-room/fact-sheets/detail/vector-borne-diseases
149. 1. R-D Xue
     2. DR Barnard
     3. GC Muller
     2009 Effects of body size and nutritional regimen on survival in adult Aedes albopictus (Diptera: Culicidae)47 Journal of medical entomology 778 782
150. 1. AA Yackel Adams
     2. SK Skagen
     3. JA Savidge
     2005 Modeling post-fledging survival of Lark Buntings in response to ecological and biological factors 87 Ecology 178 188
151. 1. WF Zuharah
     2. N Fadzly
     3. PJ Lester
     2012 Lethal and sublethal impacts of predaceous backswimmer Anisops wakefieldi (Hemiptera: Notonectidae) on the life-history traits of the New Zealand mosquito Culex pervigilans (Diptera: Culicidae)50 Journal of Medical Entomology 1014 1024