Skip to main content

Diet and density dependent competition affect larval performance and oviposition site selection in the mosquito species Aedes albopictus (Diptera: Culicidae)



Oviposition-site choice is an essential component of the life history of all mosquito species. According to the oviposition-preference offspring-performance (P-P) hypothesis, if optimizing offspring performance and fitness ensures high overall reproductive fitness for a given species, the female should accurately assess details of the heterogeneous environment and lay her eggs preferentially in sites with conditions more suitable to offspring.


We empirically tested the P-P hypothesis using the mosquito species Aedes albopictus by artificially manipulating two habitat conditions: diet (measured as mg of food added to a container) and conspecific density (CD; number of pre-existing larvae of the same species). Immature development (larval mortality, development time to pupation and time to emergence) and fitness (measured as wing length) were monitored from first instar through adult emergence using a factorial experimental design over two ascending gradients of diet (2.0, 3.6, 7.2 and 20 mg food/300 ml water) and CD (0, 20, 40 and 80 larvae/300 ml water). Treatments that exerted the most contrasting values of larval performance were recreated in a second experiment consisting of single-female oviposition site selection assay.


Development time decreased as food concentration increased, except from 7.2 mg to 20.0 mg (Two-Way CR ANOVA Post-Hoc test, P > 0.1). Development time decreased also as conspecific density increased from zero to 80 larvae (Two-Way CR ANOVA Post-Hoc test, P < 0.5). Combined, these results support the role of density-dependent competition for resources as a limiting factor for mosquito larval performance. Oviposition assays indicated that female mosquitoes select for larval habitats with conspecifics and that larval density was more important than diet in driving selection for oviposition sites.


This study supports predictions of the P-P hypothesis and provides a mechanistic understanding of the underlying factors driving mosquito oviposition site selection.


Oviposition-site choice is an essential component of the life history of all mosquito species impacting offspring survival, juvenile development and growth, intra- and inter-specific interactions, predator avoidance and, ultimately, offspring phenotype and fitness [1]. Laboratory and field studies have shown that oviposition site selection by most mosquito species is not random [25]. Yet, egg-laying behavior may be complex and highly heterogeneous. This is perhaps because selecting particular habitats is based on several aspects including physical [6] and chemical [710] cues, as well as on interspecific and intraspecific density within the habitat [5, 11]. Although it is hypothesized that oviposition site selection maximizes offspring fitness by releasing immature stages from predation or competition with other mosquito larvae or organisms over food resources [5, 12, 13], few studies have empirically estimated the association between mosquito habitat quality, resource competition, and reproductive fitness in mosquitoes [11, 1417].

Though traditionally applied to food foraging behavior, optimal foraging models provide a conceptual framework to consider oviposition behavior as well. The ideal free distribution (IFD) theory postulates that habitats differ in their suitability to support a given species, and that organisms would preferably exploit those habitats that maximize their overall fitness [18]. Moreover, habitat suitability may change over time due to crowding and density-dependence, impacting future patch exploitation dynamics [18]. The IFD theory can be applied to mosquito oviposition site selection under the hypothesis that a female will allocate the majority of her eggs in the most offspring-suitable habitats available until the pressures from conspecific competition become too strong and she is obligated to lay her eggs in the next most suitable habitats (also known as the oviposition preference-offspring performance (P-P) hypothesis [11, 1517, 19]). Under the P-P hypothesis, there is a functional relationship between the suitability of an oviposition site, the density of potential competitors and the concentration of food resources.

Aedes spp. mosquitoes, as skip ovipositors [20], distribute individual eggs in several oviposition sites rather than as clustered egg rafts [21, 22] facilitating experimental evaluation of the P-P hypothesis. Skip oviposition entails high energy consumption for oviposition site detection relative to oviposition of entire clutches in one location, potentially increasing the risk of maternal mortality. Despite this cost, it may reduce sibling competition and ensure greater distribution of progeny, thereby optimizing larval development as well as offspring fitness [20, 23]. In this study we tested the P-P hypothesis through two interrelated laboratory experiments using the Asian tiger mosquito, Aedes albopictus, as our study system. This species not only feeds on a wide array of hosts (including humans) [24] and is an efficient vector of Chikungunya virus (CHIKV) [2527], but also has demonstrated the ability to colonize diverse artificial habitats such as bird baths, flowerpots catch basins or rainwater storage containers [27, 28]. In the first experiment, we evaluated Ae. albopictus larval performance (i.e., mortality rate, time to pupation and time to adult emergence) and fitness in response to two gradients of habitat conditions: diet level and density of conspecific larvae. By recreating those two-factor treatments with the most contrasting outcomes (from the first experiment), the second experiment assessed whether Ae. albopictus females select oviposition sites following the predictions of the P-P hypothesis. Given that immature mosquito survival is density-dependent [2931] and that potential oviposition sites vary in factors like pre-existing conspecific presence and nutrient levels, we predict that a female’s preference for oviposition sites will occur as a tradeoff between both trophic dimensions: habitats with high food content and low intraspecific competition pressure should be preferred as oviposition sites.


Mosquito colonies

All experimental Ae. albopictus larvae were F3 obtained from a laboratory colony originated from wild caught eggs collected in Raleigh, NC (Apperson C., personal communication). Mosquito colonies were maintained at 28 °C and 70% relative humidity on a 12:12 (L:D) light schedule. The diet used in this study was a mixture of solid tuna meal, low-heat desiccated, non-defatted solid Argentinian beef liver powder (Now Foods, Bloomingdale, IL), and solid vitamin mix (BioServ, Frenchtown, NJ) mixed at the ratio 2:2:1 respectively (Benedict M.Q., personal communication). The solid mixture was thoroughly mixed in pure water to make a 2% diet aqueous solution for efficient diet-transfer purposes.

Experiment 1: Assessing performance and fitness over a gradient of habitat conditions

We quantified Ae. albopictus immature performance and fitness in artificial containers across gradients of food concentrations (diet) and conspecific densities (CD). Performance was characterized by measuring mortality rate (dead individuals/10), time to pupation (days since egg-hatch), and time to adult emergence (days since egg-hatch). Four diet and four CD levels were considered to generate a four-by-four experimental design with 16 treatments encompassing all possible food-conspecific combinations (Table 1). A total of 10 experimental larvae were exposed to each treatment, keeping experimental larvae separate from the conspecifics by a mesh screen divider (Additional file 1: Figure S1). The experimental larvae were carefully monitored throughout their development until adult emergence, whereas the conspecific larvae were used solely to provide the CD conditions within the container treatments. Each treatment was replicated three times (total of 48 containers) under constant temperature (28%, relative humidity (80%), and light (diel length 12:12 L:D, and with two crepuscular hours mimicking sunrise and sunset illumination) regimes.

Table 1 Distribution of treatments per experiment

The CD conditions, in ascending order, were: zero larvae, 20 larvae, 40 larvae, and 80 larvae (Table 1). Conspecific larvae were added half as first-instars and half as third-instars to simulate a structured age distribution compatible with field observations as well as to maintain conspecific presence for the duration of the experiment. Diet treatments used in this study included a daily dose of 2.0 mg, 3.6 mg, 7.2 mg, and 20.0 mg of the diet mixture in 300 ml of deionized water.

Larval rearing containers were made out of 473 ml white, plastic, cylindrical food containers (Bauman Paper Co., Lexington, KY) divided by a white mesh screen positioned perpendicularly to the bottom of the container (Additional file 1: Figure S1). The mesh screen was permeable to the diet solution but impermeable to larvae movement. Twenty-four hours prior to the beginning of the experiment, each container was filled with 300 mL of purified water with the appropriate diet mixture to allow bacterial development (a food source for Ae. albopictus[32]) and settling of food particles. Following this settling period, 10 F3 generation, recently emerged (~24-hours post hatching) Ae. albopictus, I-instar larvae were added to one side of the container (experimental larvae) with conspecific larvae added to the other side of the mesh division at densities described in Table 1. The additions of experimental and conspecific larvae marked the beginning of the experiment.

Experimental larval development was monitored daily by counting the number of molts to each life stage starting from II-instar to pupal stage (after counting, molts were removed from the container). When treatment larvae reached the pupal stage, the pupae were removed from the container (they can be extracted without affecting per-capita food consumption of the remaining larvae). When experimental larvae pupated, the pupae were kept inside the container and the molt(s) removed. Each emerged adult from the experimental larvae was sexed and all emerged females were individually stored at −80° C for further wing-length measurements.

Wing length is considered a proxy of female body size and overall mosquito fitness [33, 34]. Female adult body size is important to assess fitness because larger females lay more eggs [35, 36], thus increasing the probability of successful progeny. Wing-length measurements were conducted for a subset consisting of two randomly selected replicates of each treatment. Wings were separated from the body of the adult females under a Leica MS5 stereozoom microscope (Leica Microsystems Incorporated, Bannockburn, IL). Following dissection, photographs of the wings were taken using a MagnaFire 2.1C CCD camera (Optronics, Goleta, CA) under an Olympus BX60F-3 microscope (Olympus Optical Co., Ltd., Japan) at 4X magnification. Wing length measurements were conducted using the image-processing program ImageJ [37].

Experiment 2: Empirical estimation of oviposition site preference

Twelve different containers with unique diet-CD combinations were used in single-female oviposition assays (Table 1). As before, half of the conspecific larvae were added as I-instars and half as III-instars. Coffee filter paper was lined along the walls of each container to collect mosquito eggs. The 12 containers were randomly positioned in an enclosed insect tent (75 × 75 × 115 cm; MegaView Science Co., Ltd., Taichung, Taiwan) and left for 48 hours to allow for bacterial growth.

Adult F3 mosquitoes were used for the oviposition assays. Upon adult emergence, a female was placed in a 20x20x20 cm acrylic chamber with male mosquitoes for 24 hours to mate, and provided access to a sedated rabbit on the following day. Use of laboratory animals for blood-feeding mosquitoes was approved by the Centers for Disease Control and Prevention Animal Use Protocols (IACUC protocol # 2245BARRABB). One day after feeding, gravid females were released into the experimental tent containing the array of oviposition sites and left for three days to oviposit the initial batch of eggs. We hypothesized that the initial oviposition would show the strongest correlation with oviposition site selection, given the absence of preexisting eggs, which can act as oviposition attractants [23, 38]. Environmental conditions were set as described in the first experiment. Upon completion of each 3-day assay, females were removed from the tents, filter papers were removed from each container and the number of eggs counted. The oviposition trials were replicated four times.

Data analysis

Experiment 1

Daily larval survival under different diet and CD combinations was analyzed using Kaplan-Meier Survival Analysis tests. Mantel-Cox Log-Rank tests were performed to determine whether increases in diet and/or CD significantly affected mosquito survival. Two-way Completely Randomized (CR) Analysis of Variance (ANOVA) tests compared the mean time to pupation and time to adult emergence (log10-transformed) across different treatments. The mean of the 10 experimental larvae for both time to pupation and time to emergence was computed for each container before the CR ANOVA. These computed means represented independent observations characterizing each container and thus validated the use of a CR ANOVA as opposed to a repeated/related measures ANOVA [39]. Percent pupation and emergence success rate were highly skewed and thus analyzed using the Kruskal-Wallis test.

A Generalized Additive Mixed Model (GAMM) was applied to determine the association between diet, CD levels and adult mosquito wing length [10, 40]. The full model had the form:

Y wing = α + f 1 Diet + f 2 C D + f 3 Diet C D + Z i + ε j

where, ε j N Ο , σ 2 . The parameter Zi represents a random effects term associated with each observation within a replicate. Given the fairly Normal wing length distribution we used a Gaussian link function to parameterize the model. GAMM models allow for non-linear relationships between the response variable and multiple explanatory variables by incorporating a smoothing function representing the additive component [40]. We fitted by applying a penalized cubic spline function to the data [40]. Four models were developed: one with each variable alone, one model with both variables non-linearly related with each other (expressed as Diet ~ CD) and the full model including all possible combinations of variables. We assessed each model’s ability to fit the data by comparing their Akaike Information Criterion (AIC) scores after a correction for small sample sizes (AICc). A model with ΔAICc = 2 or more units lower than any other model was considered the best, and used to predict the association between wing length, diet and CD. Once the best model was identified, we plotted each predicted f i as either a curve (single terms) or a surface (diet ~ CD) and assessed the model’s overall fit by reporting pseudo-R2 and deviance parameters.

Experiment 2

The IFD theory predicts that the number of eggs oviposited will be equally distributed among patches, or cups, with the same intrinsic quality. In habitats that differ in ‘quality’ (measured as the concentration of food and conspecifics), it is predicted that female mosquitoes will deposit eggs in high quality containers first. We developed a container quality index (QI), by first ranking food and CD from low to high and then assigning them values from 1 to 5 (Table 1). By adding the ranked values of food and CD to each treatment we developed the QI index (which ranges from 2 to 10). At the low end of the scale are treatments with no or low diet and high density. The high end of the scale represents treatments with high diet and no or low density (Table 1).

To quantify the association between diet, CD and number of eggs per container a general additive mixed model GAMM was applied. The full model had the form:

Y eggs = α + f 1 Diet + f 2 C D + F 3 Diet C D + Z i + ε j

where, ε j N Ο , σ . The parameter Zi represents a random effects term associated with each observation within a replicate. Given the distribution of eggs per container was skewed (i.e., many containers had no eggs, whereas a few containers had large number of eggs) we used a negative binomial link function to parameterize the model. The GAMM model was tested as before, and included the number of eggs per container on each of the four trials as observations. GAMM models were run in R statistical software’s package mgcv [41]. All other statistical analyses were performed in SPSS 17.0.0 [42].


Experiment 1: Assessing performance and fitness over a gradient of habitat conditions

Of the 480 experimental larvae, 35 (7.3%) died before emergence. Diet levels alone (i.e., in the absence of conspecific larvae) did not significantly affect immature lifespan (log rank Mantel-Cox, X2= 0.492, P > 0.1) (Figure 1A). Similarly, conspecific density (CD) alone did not significantly affect lifespan per larva (log rank, Mantel-Cox, X2= 2.858, P > 0.1) (Figure 1B). Overall pupation rate was 93.3% and mean time to pupation was 7.4 days (SD ± 1.93 days). The container treatment resulting in experimental larvae with the longest time to pupation was 2.0 mg diet and 80 CD (see Table 1 for full description of treatments), pupating at 11.4 days (SD ± 0.88 days), whereas the treatment with the shortest time to pupation was 20.0 mg and 0 CD, pupating at 5.5 days (SD ± 0.14 days) (Figure 2). Diet and CD (Two-Way CR ANOVA, diet: F (d.f. =3) = 57.6, P < 0.001; CD: F (3) = 3.6, P < 0.05) but not their interaction (F (9) = 1.5, P > 0.1) were significantly associated with time to pupation. Post-Hoc tests indicated that time to pupation was significantly shorter with each increasing diet level except from 7.2 mg to 20.0 mg (Two-Way CR ANOVA Post-Hoc test, P > 0.1) (Figure 2, Table 2). In addition, time to pupation significantly increased when CD increased from zero larvae to 80 larvae (Two-Way CR ANOVA Post-Hoc test, P < 0.05).

Figure 1

Impact of diet and CD on larval survival. (A). Survival Analysis Curve using diet (expressed in mg/300 ml) as a potential predictor of mean lifespan (per container type) across all replicates (N = 3). The Kaplan Meier test showed no difference in survival across diet treatments (P = 0.921). (B) Survival Analysis Curve using conspecific density (CD) as a potential predictor of mean lifespan (per container type) across all replicates (N = 3). The Kaplan Meier test showed no significant difference in treatments across CD (P = 0.414).

Figure 2

Development time to larval pupation across ascending diet levels by conspecific density. (A) zero conspecifics. (B) 20 conspecifics. (C) 40 conspecifics. (D) 80 conspecifics. Mean values are indicated in boxes. Diet is indicated in mg/300 ml of water. Each treatment was replicated three times.

Table 2 Association between time to pupation, diet and CD

A total of 445 larvae (92.7%) emerged as adults, comprising 98.0% of the total number of larvae that successfully pupated. Overall, the treatment with the longest time to adult emergence was 2.0 mg and 80 CD, at 13.5 days [SD ± 0.97 days], whereas the treatment with the shortest time to emergence was 20.0 mg and 0 CD, 7.7 days [SD ± 0.10 days]. These trends were qualitatively the same as the ones found for time to pupation (not shown). Diet and CD were both significant predictors of time to emergence (Two-Way CR ANOVA, diet: F (3) = 49.0, P < 0.001; CD: F (3) = 3.0, P < 0.05). The interaction between diet and CD was not significantly associated with the time to emergence (F (9) = 0.9, P > 0.1). Post-Hoc tests indicated that time to emergence was significantly affected with each increasing diet level except from 7.2 mg to 20.0 mg (Two-Way CR ANOVA Post-Hoc test, P > 0.1) (Table 3). In addition, time to emergence only increased significantly from CD of zero larvae to 80 larvae (Two-Way ANOVA Post-Hoc test, P < 0.05). Increasing CD levels from zero to 40 larvae did not significantly affect development time, suggesting that at least 80 conspecific larvae/300 ml are needed to induce intraspecific crowding and resource competition intense enough to significantly impact larval performance.

Table 3 Association between time to emergence, diet and CD Results from post-hoc comparisons performed after a Two-Way completely randomized analysis of variance (ANOVA) of time to emergence (days post-egg hatch), diet and conspecific density

A total of 131 females were dissected for wing length measurements. Wing length was positively and significantly associated with the container quality index (linear regression coefficient, b = 0.07;SE = 0.007; P < 0.001; Figure 3); at higher QI values mosquito fitness was predicted to be highest. The best GAM model (ΔAICc = 7 in comparison to the second best model) included the non-linear association between diet and CD (diet ~ CD) as its sole and significant term (Table 4). This model had an adjusted R2 of 0.60 and explained 62% of the deviance. Figure 4A shows the non-linear association between diet, CD, and wing length. Diet had a higher impact on wing length than CD; the number of conspecifics negatively affected wing length only when food levels were low (<5 mg). At intermediate diet levels (>10 mg), increases in CD did not exert any significant change in wing length (Figure 4A). Figure 4B shows the effect of varying diet and CD on adult wing length which, given the association between wing length and reproductive fitness, can be interpreted as an Ae. albopictus fitness surface. The diet and CD combination that produced the highest mosquito fitness were >11 mg and < 60 conspecifics/300 ml, respectively (Figure 4B).

Figure 3

Association between wing length and the container quality index. Boxplot indicating median (horizontal line), interquartile (box) and range (dotted lines). Dots represent outlier observations.

Table 4 Summary of the best GAMM model estimating the association between wing length, diet and conspecific density
Figure 4

Association between wing length, diet and CD. (A) Graphic representation of the association between diet, conspecific density and wing length (parameter Diet ~ CD in from the best GAM model). (B) Surface representing the association between diet and conspecific density (surface colors) and wing length (red lines). Numbers on top of lines indicate wing length values whereas colors interpolated wing length values.

Experiment 2: Empirical estimation of oviposition site preference

Because larval performance (i.e. time to pupation and adult emergence) differed significantly only between diets 2.0, 3.6, and 7.2 mg, these levels plus a level of no food were used as the experimental diets in the oviposition assays. We also concluded in the first experiment that CD only had an effect on larval performance in the range 0 to 80 conspecific larvae. Also, CD only had an effect on fitness (i.e. wing length) when the diet was <5 mg. Therefore, the CD levels used to test oviposition site selection were 0, 10, and 80 conspecific larvae.

The total numbers of eggs per replicate were: 40, 81, 101 and 122. The observed numbers of eggs laid in each container did not follow a pattern consistent with the expectation of the P-P hypothesis: that increases in the container quality index will be associated with increases in the number of eggs (Figure 5). Instead, most females laid eggs at intermediate container quality indices (Figure 5). Females whose total egg batch was low (~40-80 eggs distributed among all containers), oviposited heavily in containers harboring fewer conspecific larvae (i.e. CD of 10 larvae), whereas females with large egg batches (~100-125 eggs), oviposited preferentially in the high CD containers (i.e. 80 larvae). The best GAM model (ΔAICc = 25) included the non-linear diet ~ CD term as the sole and significant predictor of the number of eggs per container (Table 4, Figure 6). This model had an adjusted R2 of 0.1 and explained 23% of the deviance. CD had a larger impact on oviposition site selection than diet (Figure 6); the number of eggs was negatively associated with the number of conspecifics in a container at intermediate (20–50 per 300 ml container) larval densities. Figure 6 shows the importance of the simultaneous consideration of diet and CD when assessing the oviposition preference of a skip ovipositor such as Ae. albopictus.

Figure 5

Results from the oviposition site selection trials. Association between the number of eggs per container and the container quality index. Each symbol represents a replicate experiment in which a single female was given the choice to oviposit in containers with different combinations of diet and conspecific density (CD) of oviposition sites.

Figure 6

Association between oviposition, diet and CD. Surface representing the association between diet and conspecific density and number of laid eggs (surface colors) on a given container (term Diet ~ CD from the best GAM model). Numbers on top of lines indicate the natural logarithm of the number of eggs laid on a given container. Colors indicate interpolated number of eggs per container.


The P-P hypothesis proposes that females should oviposit in suitable habitats to minimize resource-mediated interactions while maximizing offspring fitness [11, 15, 16, 19]. This study confirmed that oviposition site selection by Ae. albopictus followed the predictions of the P-P hypothesis for oviposition site selection: 1) females allocated most of their eggs in an array of habitats that conferred high offspring performance and fitness, and 2) conspecific larvae presence and density were the most important factors influencing site selection within a medium to low diet gradient.

Previous studies have described that high diet conditions correlate with shorter Ae. albopictus development times [5, 31]. In the present study, larval development time decreased with increasing diet level, and such findings were not significantly influenced by the presence of conspecific individuals at medium diet concentrations. The lack of a scaling of development time with diet at high food concentrations indicates that diet only affected larval performance within the range in which competition for food resources operated (i.e., high diets provided more resources than the 10 experimental larvae could compete for) [43], a pattern similarly observed for Anopheles arabiensis mosquitoes [44].

Density-dependent competition for food during the early larval stages is considered to be one of the most important factors affecting mosquito population dynamics [31, 44, 45]. Overcrowding of mosquito larval habitats generally results in retarded growth, high mortality, small and non-uniform sized adults, and decreased fecundity. Although mortality rates differed between treatments, prolonged development time was exclusively observed in high-density treatments (80 conspecific larvae). Such observed density-dependent effect could have resulted from crowding effects (i.e., the direct impact of larger numbers of organisms on access to food and, ultimately, feeding success; [4, 5, 46]), increased food partitioning (lower food per capita; [46]), or toxins induced by such sources of stress [47]. By separating conspecific from experimental larvae and keeping the experimental larvae at low densities (10 larvae/150 ml of water), this study minimized any possible occurrence of crowding effects focusing on the effect of food partitioning on larval productivity. In our study, the estimated food available per capita (per larva) ranged from 0.00067 mg/mL (2 mg/300 mL per 10 individuals) to 0.0024 mg/mL (7.2 mg/300 mL per 10 individuals), which is comparable to previous studies testing diet gradients ranging from severely limiting to surplus levels [43, 44]. High conspecific larval density is associated with high levels of toxin production [47] and depletion of haemolymph glucose levels [38]. Thus, in our study stunting or prolonging development time to the adult stage could have been a product of both resource partition and sub-products of overcrowded larval habitats. By studying the chemical ecology of oviposition sites in association with mosquito population dynamics, future research may be able to further identify the biological processes underlying oviposition site selection by Ae. albopictus.

Adult female body size (i.e. quantified as wing length) was measured to estimate maternal fitness resulting from the tested diet and CD conditions. Female fitness is particularly important because larger females lay more eggs, contributing more to a population’s reproductive fitness [35, 36]. As predicted, wing length increased proportionally with increases in diet with the highest diet level (20.0 mg) resulting in the longest mean wing length. That the number of conspecific larvae negatively affected wing length only when food levels were low (<5 mg) suggests that shorter development and thus smaller body size outweigh the costs associated with high CD under nutrient limitations. Intense competition (most likely elicited under <11 mg of diet and >60 conspecific larvae [Figure 4B]) may decrease an individual larva’s probability of emergence and survival; therefore, larger body size may need to be sacrificed to ensure adult emergence and thereby maintain (or increase) overall adult population size [48].

In the oviposition assays, females tended to lay eggs in sites where conspecific larvae were present. This finding was counter to expectations following Experiment 1, which concluded that high competition resulting from low food availability negatively affects larval performance and fitness. It would seem advantageous for ovipositing females to avoid sites containing potential competitors to their own progeny. However, the benefit of using conspecifics as a cue of habitat suitability may counteract the potential drawback of conspecific competition [7, 49]. Wong et al. [50] hypothesized that the presence of Ae. aegypti conspecifics may serve as a signal that the site experiences infrequent water turnover and desiccation, and contains adequate food, two critical conditions for larval development. Whether females select oviposition sites due to their adequate larval conditions or other variables naturally found in oviposition sites (e.g. different types of diet, site color, presence of other mosquito and non-mosquito species) remains to be determined.

This study presented several limitations. The estimation of wing measurements from only two replicates affected the fit of the GAM models to the data (as observed in the low deviance explained by the model), which could have been improved with a larger sample size. Also, the CD treatments selected in this study were based on field observations performed in a residential neighborhood of Atlanta, GA, (MY unpublished data) and may not necessarily reflect the full range of naturally-occurring CDs in which Ae. albopictus larvae could develop. Also, the use of a water-based food formulation instead reduced the amount of debris and homogenized the color of the water between containers. By performing single female assays we were able to easily quantify oviposition site selection. However, the inclusion of more than one female could have provided additional insights about the role of female density on oviposition site selection. It is unknown whether the presence of other gravid females can act as an additional oviposition cue, and such behavioral mechanisms deserve further investigation.


Current control strategies have shown transient or no perceivable impacts in Ae. albopictus populations, and point to the need for a better understanding of the biological and ecological factors regulating natural populations [27]. Our findings expand current knowledge on the ecology of Ae. albopictus by showing that conspecific density may be an important driver in oviposition site selection. Further understanding of the mechanisms facilitating successful Ae. albopictus oviposition behavior will aid in the development of improved control strategies, particularly at the reproductive level.


  1. 1.

    Resetarits WJ: Oviposition site choice and life history evolution. Am Zool. 1996, 36: 205-215.

    Article  Google Scholar 

  2. 2.

    Chaves LF, Keogh CL, Vazquez-Prokopec GM, Kitron UD: Combined sewage overflow enhances oviposition of Culex quinquefasciatus (Diptera: Culicidae) in urban areas. J Med Entomol. 2009, 46 (2): 220-226. 10.1603/033.046.0206.

    Article  PubMed  Google Scholar 

  3. 3.

    Yap HH, Lee CY, Chong NL, Foo AE, Lim MP: Oviposition site preference of Aedes albopictus in the laboratory. J Am Mosq Control Assoc. 1995, 11 (1): 128-132.

    CAS  PubMed  Google Scholar 

  4. 4.

    Zahiri N, Rau ME: Oviposition attraction and repellency of Aedes aegypti (Diptera: Culicidae) to waters from conspecific larvae subjected to crowding, confinement, starvation, or infection. J Med Entomol. 1998, 35 (5): 782-787.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Zahiri N, Rau ME, Lewis DJ: Oviposition responses of Aedes aegypti and Ae. atropalpus (Diptera: Culicidae) females to waters from conspecific and heterospecific normal larvae and from larvae infected with Plagiorchis elegans (Trematoda: Plagiorchiidae). J Med Entomol. 1997, 34 (5)): 565-568.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Kitron UD, Webb DW, Novak RJ: Oviposition behavior of Aedes triseriatus (Diptera: Culicidae): prevalence, intensity, and aggregation of eggs in oviposition traps. J Med Entomol. 1989, 26 (5): 462-467.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Bentley MD, Day JF: Chemical ecology and behavioral aspects of mosquito oviposition. Ann Rev Entomol. 1989, 34: 401-421. 10.1146/annurev.en.34.010189.002153.

    CAS  Article  Google Scholar 

  8. 8.

    Li J, Deng T, Li H, Chen L, Mo J: Effect of water color and chemical compounds on the oviposition behavior of gravid Culex pipiens pallens females under laboratory conditions. J Agricult Urban Entomol. 2009, 26: 23-30. 10.3954/1523-5475-26.1.23.

    Article  Google Scholar 

  9. 9.

    Trexler JD, Apperson CS, Zurek L, Gemeno C, Schal C, Kaufman M, Walker E, Watson DW, Wallace L: Role of bacteria in mediating the oviposition responses of Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2003, 40 (6): 841-848. 10.1603/0022-2585-40.6.841.

    Article  PubMed  Google Scholar 

  10. 10.

    Peterkova-Koci K, Robles-Murguia M, Ramalho-Ortigao M, Zurek L: Significance of bacteria in oviposition and larval development of the sand fly Lutzomyia longipalpis. Parasit Vectors. 2012, 5: 145-10.1186/1756-3305-5-145.

    PubMed Central  Article  PubMed  Google Scholar 

  11. 11.

    Ellis AM: Incorporating density dependence into the oviposition preference-offspring performance hypothesis. J Anim Ecol. 2008, 77 (2): 247-256. 10.1111/j.1365-2656.2007.01345.x.

    Article  PubMed  Google Scholar 

  12. 12.

    Benzon GL, Apperson CS: Reexamination of chemically mediated oviposition behavior in Aedes aegypti (L.) (Diptera: Culicidae). J Med Entomol. 1988, 25 (3): 158-164.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Silver JB: Mosquito Ecology: Field Sampling Methods. 2008, New York: Springer, 3

    Google Scholar 

  14. 14.

    Heard SB: Imperfect Oviposition Decisions by the Pitcher Plant Mosquito (Wyeomyia-Smithii). Evol Ecol. 1994, 8 (5): 493-502. 10.1007/BF01238254.

    Article  Google Scholar 

  15. 15.

    Nufio CR, Papaj DR: Superparasitism of larval hosts by the walnut fly, Rhagoletis juglandis, and its implications for female and offspring performance. Oecologia. 2004, 141 (3): 460-467.

    Article  PubMed  Google Scholar 

  16. 16.

    Valladares G, Lawton JH: Host-Plant Selection in the Holly Leaf-Miner - Does Mother Know Best. J Anim Ecol. 1991, 60 (1): 227-240. 10.2307/5456.

    Article  Google Scholar 

  17. 17.

    Wong J, Morrison AC, Stoddard ST, Astete H, Chu YY, Baseer I, Scott TW: Linking Oviposition Site Choice to Offspring Fitness in Aedes aegypti: Consequences for Targeted Larval Control of Dengue Vectors. PLoS Negl Trop Dis. 2012, 6 (5): e1632-10.1371/journal.pntd.0001632.

    PubMed Central  Article  PubMed  Google Scholar 

  18. 18.

    Fretwell SD, Lucas HL: On territorial behavior and other factors influencing habitat distribution in birds I. Theoretical development. Acta Biotheoretica. 1970, 19: 16-36.

    Article  Google Scholar 

  19. 19.

    Rausher MD: Ecology of host-selection behaviour in phytophagous insects. 1983, New York: Academic

    Google Scholar 

  20. 20.

    Mogi M, Mokry J: Distribution of Wyeomyia smithii (Diptera: Culicidae) eggs in pitcher plants in Newfoundland, Canada. Trop Med. 1980, 22: 1-12.

    Google Scholar 

  21. 21.

    Apostol BL, Black WC, Reiter P, Miller BR: Use of randomly amplified polymorphic DNA amplified by polymerase chain reaction markers to estimate the number of Aedes aegypti families at oviposition sites in San Juan, Puerto Rico. Am J Trop Med Hyg. 1994, 51 (1): 89-97.

    CAS  PubMed  Google Scholar 

  22. 22.

    Chadee DD, Corbet PS: Seasonal incidence and diel patterns of oviposition in the field of the mosquito, Aedes aegypti (L.) (Diptera: Culicidae) in Trinidad, West Indies: a preliminary study. Ann Trop Med Parasitol. 1987, 81 (2): 151-161.

    CAS  PubMed  Google Scholar 

  23. 23.

    Chadee DD, Corbet PS, Greenwood JJD: Egg-laying yellow fever mosquitoes avoid sites containing eggs laid by themselves or by conspecifics. Ent Exp Appl. 1990, 57: 295-298. 10.1111/j.1570-7458.1990.tb01442.x.

    Article  Google Scholar 

  24. 24.

    Kamgang B, Nchoutpouen E, Simard F, Paupy C: Notes on the blood-feeding behavior of Aedes albopictus (Diptera: Culicidae) in Cameroon. Parasit Vectors. 2012, 5: 57-10.1186/1756-3305-5-57.

    PubMed Central  Article  PubMed  Google Scholar 

  25. 25.

    Black WC, Bennett KE, Gorrochotegui-Escalante N, Barillas-Mury CV, Fernandez-Salas I, de Lourdes Munoz M, Farfan-Ale JA, Olson KE, Beaty BJ: Flavivirus susceptibility in Aedes aegypti. Arch Med Res. 2002, 33 (4): 379-388. 10.1016/S0188-4409(02)00373-9.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Delatte H, Gimonneau G, Triboire A, Fontenille D: Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009, 46 (1): 33-41. 10.1603/033.046.0105.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Paupy C, Delatte H, Bagny L, Corbel V, Fontenille D: Aedes albopictus, an arbovirus vector: from the darkness to the light. Microbes Infect. 2009, 11 (14–15): 1177-1185.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Braks MAH, Honorio NA, Lounibos LP, Lourenco-De-Oliveira R, Juliano SA: Interspecific competition between two invasive species of container mosquitoes, Aedes aegypti and Aedes albopictus (Diptera: Culicidae), in Brazil. Ann Entomol Soc Am. 2004, 97 (1): 130-139. 10.1603/0013-8746(2004)097[0130:ICBTIS]2.0.CO;2.

    Article  Google Scholar 

  29. 29.

    Arrivillaga J, Barrera R: Food as a limiting factor for Aedes aegypti in water-storage containers. J Vector Ecol. 2004, 29 (1): 11-20.

    PubMed  Google Scholar 

  30. 30.

    Southwood TR, Murdie G, Yasuno M, Tonn RJ, Reader PM: Studies on the life budget of Aedes aegypti in Wat Samphaya, Bangkok, Thailand. Bull World Health Organ. 1972, 46 (2): 211-226.

    PubMed Central  CAS  PubMed  Google Scholar 

  31. 31.

    Walsh RK, Facchinelli L, Ramsey JM, Bond JG, Gould F: Assessing the impact of density dependence in field populations of Aedes aegypti. J Vector Ecol. 2011, 36 (2): 300-307. 10.1111/j.1948-7134.2011.00170.x.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Merritt RW, Dadd RH, Walker ED: Feeding-Behavior, Natural Food, and Nutritional Relationships of Larval Mosquitos. Ann Rev Entomol. 1992, 37: 349-376. 10.1146/annurev.en.37.010192.002025.

    CAS  Article  Google Scholar 

  33. 33.

    Clements A: The Biology of Mosquitoes. Vol. 1: Development, Nutrition, and Reproduction. 1992, London: Chapman & Hall

    Google Scholar 

  34. 34.

    Nasci RS: The Size of Emerging and Host-Seeking Aedes-Aegypti and the Relation of Size to Blood-Feeding Success in the Field. J Am Mosquito Contr. 1986, 2 (1): 61-62.

    CAS  Google Scholar 

  35. 35.

    Briegel H: Metabolic Relationship between Female Body Size, Reserves, and Fecundity of Aedes-Aegypti. J Insect Physiol. 1990, 36 (3): 165-172. 10.1016/0022-1910(90)90118-Y.

    Article  Google Scholar 

  36. 36.

    Steinwascher K: Relationship between Pupal Mass and Adult Survivorship and Fecundity for Aedes-Aegypti. Environ Entomol. 1982, 11 (1): 150-153.

    Article  Google Scholar 

  37. 37.

    Rasband MD: Bethesda, Maryland. ImageJ. 1983, USA: U.S: National Institutes of Health

    Google Scholar 

  38. 38.

    Zahiri N, Dunphy GB, Rau ME: Serum composition of Aedes aegypti (Diptera: Culicidae) larvae and the production of an oviposition repellent are influenced by infection with the entomopathogenic digenean Plagiorchis elegans (Trematoda: Plagiorchiidae), starvation, and crowding. J Med Entomol. 1998, 35 (2): 162-168.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Chaves LF: An Entomologist Guide to Demystify Pseudoreplication: Data Analysis of Field Studies With Design Constraints. J Med Entomol. 2010, 47 (3): 291-298. 10.1603/ME09250.

    Article  PubMed  Google Scholar 

  40. 40.

    Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM: Mixed Effects Models and Extensions in Ecology with R. 2009, New York: Springer

    Google Scholar 

  41. 41.

    Wood SN: Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J R Stat Soc B. 2011, 73: 3-36. 10.1111/j.1467-9868.2010.00749.x.

    Article  Google Scholar 

  42. 42.

    IBM SS: Chicago. SPSS 17.0. 2008, Illinois, USA: SPSS Inc

    Google Scholar 

  43. 43.

    Washburn JO, Hartmann EU: Could Aedes albopictus (Diptera: Culicidae) become established in California tree holes?. J Med Entomol. 1992, 29 (6): 995-1005.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Gilles JR, Lees RS, Soliban SM, Benedict MQ: Density-dependent effects in experimental larval populations of Anopheles arabiensis (Diptera: Culicidae) can be negative, neutral, or overcompensatory depending on density and diet levels. J Med Entomol. 2011, 48 (2): 296-304. 10.1603/ME09209.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Moore CG, Fisher BR: Competition in mosquitoes. Density and species ratio effects on growth, mortality, fecundity, and production of growth retardant. Ann Entomol Soc Am. 1969, 62 (6): 1325-1331.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Dye C: Competition Amongst Larval Aedes-Aegypti - the Role of Interference. Ecol Entomol. 1984, 9 (3): 355-357. 10.1111/j.1365-2311.1984.tb00859.x.

    Article  Google Scholar 

  47. 47.

    Kuno G, Moore CG: Production of Larval Growth-Retardant in Axenic Cultures of Aedes aegypti. Mosq News. 1975, 35 (2): 199-201.

    Google Scholar 

  48. 48.

    Dieng H, Mwandawiro C, Boots M, Morales R, Satho T, Tuno N, Tsuda Y, Takagi M: Leaf litter decay process and the growth performance of Aedes albopictus larvae (Diptera: Culicidae). J Vector Ecol. 2002, 27 (1): 31-38.

    PubMed  Google Scholar 

  49. 49.

    Stav G, Blaustein L, Margalith J: Experimental evidence for predation risk sensitive oviposition by a mosquito, Culiseta longiareolata. Ecol Entomol. 1999, 24: 202-207. 10.1046/j.1365-2311.1999.00183.x.

    Article  Google Scholar 

  50. 50.

    Wong J, Stoddard ST, Astete H, Morrison AC, Scott TW: Oviposition site selection by the Dengue vector Aedes aegypti and its implications for Dengue Control. PLoS Negl Trop Dis. 2011, 5: e1015-10.1371/journal.pntd.0001015.

    PubMed Central  Article  PubMed  Google Scholar 

Download references


We gratefully thank Paul Howell for helping with mosquito bloodfeeding and colony maintenance. This study was supported by funds from Emory University (Department of Environmental Studies James G. Lester Research Grant) and the Scholarly Inquiry and Research at Emory (SIRE) program through a student fellowship to M.Y. Development of the ideas presented here was assisted by support from Research and Policy for Infectious Disease Dynamics (RAPIDD) program of the Science and Technology Directorate, U.S. Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.

Author information



Corresponding author

Correspondence to Gonzalo M Vazquez-Prokopec.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors participated in the conceptualization of the study. MY, FK, JC and JM conducted laboratory experiments. MY and GVP analyzed the data. MY, and GVP drafted the manuscript. All authors participated in the revision of the manuscript and approved the submitted version. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Yoshioka, M., Couret, J., Kim, F. et al. Diet and density dependent competition affect larval performance and oviposition site selection in the mosquito species Aedes albopictus (Diptera: Culicidae). Parasites Vectors 5, 225 (2012).

Download citation


  • Habitat selection
  • Oviposition
  • Diet
  • Conspecific density
  • Mosquito