The influence of larval competition on Brazilian Wolbachia-infected Aedes aegypti mosquitoes
© Dutra et al. 2016
Received: 28 December 2015
Accepted: 2 May 2016
Published: 16 May 2016
With field releases starting in Brazil, particular interest must be given to understanding how the endosymbiotic bacterium Wolbachia pipientis affects Aedes aegypti mosquitoes with a Brazilian genetic background. Currently, there is limited information on how the bacterium affects phenotypic traits such as larval development rate, metabolic reserves and morphometric parameters in Ae. aegypti. Here, we analyze for the first time, the effect of Wolbachia on these key phenotypes and consider how this might impact the potential of the bacterium as a disease control agent in Brazil.
We examined the influence of the wMel strain of Wolbachia in laboratory Ae. aegypti with a Brazilian genetic background, reared under different larval densities. Pupae formation was counted daily to assess differences in development rates. Levels of metabolic reserves and morphometric parameters were assessed in adults resulting from each larval condition.
wMel infection led to more rapid larval development at higher densities for both males and females, with no effect under less crowded conditions in females. Infection also led to reduced body size at both high and low density, but not at intermediate density, although the scale of this difference was maintained regardless of larval density, in comparison to uninfected individuals. Wing shape also varied significantly between infected and uninfected mosquitoes due to larval density. Glycogen levels in uninfected mosquitoes decreased under higher larval density, but were consistently high with Wolbachia infection, regardless of larval density.
We demonstrate that the wMel Wolbachia strain can positively influence some important host fitness traits, and that this interaction is directly linked to the conditions in which the host is reared. Combined with previously published data, these results suggest that this Wolbachia strain could be successfully used as part of the Eliminate Dengue Program in Brazil.
Mosquitoes of the family Culicidae are considered the most important group of insects involved in disease transmission in humans. The mosquito Aedes aegypti is primarily responsible for the transmission of dengue, the most prevalent mosquito-borne viral disease , which has seen a 30-fold increase in incidence over the last half century, dramatically increasing the burden on human health . The transmission cycle of dengue and other vector-borne diseases depends on the existence of a sophisticated tripartite relationship between the pathogen and invertebrate and vertebrate hosts. Vector competence in mosquitoes is dependent on a wide variety of extrinsic and intrinsic factors [3, 4].
Dengue transmission is therefore strongly linked to host fitness and physiology. As such, environmental conditions during larval development can affect many important life history traits during adulthood, including fecundity, fertility, immune response and host seeking [5–7]. These environmental conditions include the availability of nutritious particles dissolved in the water of the breeding sites that larvae use as food, the influence of temperature, and crowding condition itself . Competition for these resources between larvae can detrimentally affect development time  and the levels of key adult energy reserves such as lipids or glycogen , the primary carbohydrate reserves in insects . Glycogen is the principal energetic component utilized by mosquitoes to fuel flight , particularly over short distances, as undertaken during host-seeking by female mosquitoes . The conditions experienced during larval development directly influence adult mosquito size and fitness, with smaller mosquitoes experiencing impairments in blood-feeding , flight  and mating performance . Smaller mosquitoes also display a high susceptibility to infection, increased dissemination rate of dengue virus (DENV) infection  and have altered host gene expression , all of which affect mosquito vector competence . As such, larval competition, commonly present in nature [19–21], plays a major role in shaping the vectorial capacity of a mosquito species.
Wolbachia pipientis (hereafter Wolbachia) is amongst the most widespread endosymbiotic bacteria known . This species is currently being evaluated as a novel biocontrol agent to reduce dengue transmission, by exploiting some of the many phenotypic changes the bacteria make to their host’s biology . Wolbachia’s effects can range from fitness costs [24–27] to fitness benefits [28, 29]. Wolbachia is not naturally present in Ae. aegypti , but several strains have been artificially introduced into this mosquito via the process of transinfection [31, 32]. Critically and encouragingly from the perspective of disease control, Wolbachia limits infection by key pathogens in the mosquito host, including dengue virus [32–34]. Consequently, it represents a promising biocontrol agent that is currently being used in population replacement strategies around the world [32, 35–37].
There are many physiological changes associated with Wolbachia infection in adult Ae. aegypti, some of which have an associated fitness cost [38, 39]. However, there is relatively little information on the role played by Wolbachia in the development of immature mosquitoes and how this might contribute to adult fitness. Any potential fitness effects from altered development in Wolbachia-infected mosquitoes must be considered when using the bacterium as a biocontrol agent, as it could affect their ability to compete against wild type mosquitoes in the field, and thus hamper bacterial spread . In Australian Ae. aegypti mosquitoes, the wMelPop Wolbachia strain causes infected larvae to develop faster than uninfected individuals under low nutritional conditions and high larval density, whereas the opposite effect occurred when food availability was high and larval density low . Both the wMel and wMelPop strains caused a delay in larval development time compared to uninfected larvae in mixed cohorts of infected and uninfected individuals .
Wolbachia infection can affect levels of key nutritional reserves, as seen with the native wFlu strain in the mosquito Aedes fluviatilis, where levels of glycogen were higher in infected mosquitoes . Wolbachia infection can also lead to altered adult body size, with wMelPop causing a decrease and wMel an increase in comparison to uninfected mosquitoes ; however, this was specifically linked to the conditions in which the experiments were performed, given that a further set of experiments saw no effect of Wolbachia in the size of mosquito larvae reared under extremely nutritionally deprived conditions .
Releases of wMel-infected Ae. aegypti have recently commenced in Brazil as part of the Eliminate Dengue Program. These mosquitoes possess a Brazilian genetic background, different to that of mosquitoes used in previous releases in Australia . It is possible that different host genetic backgrounds could influence the fitness effects caused by Wolbachia, hindering the spread of this bacterium in the field, or induce fitness benefits that could actually enhance the chance of successful population invasion [29, 44]. So far, the effect of wMel on fecundity, maternal transmission and cytoplasmic incompatibility (CI) in Brazilian Ae. aegypti have been characterized, where the wMel strain, similarly to what was observed in the Australian genetic background, causes strong CI, a high rate of maternal transmission and has no evident detrimental effect on host fecundity or fertility . However, the effects of larval competition on the fitness of Wolbachia-infected mosquitoes are unknown. In this study, we explored how changes in larval density conditions affected both larval and adult stages of Brazilian Wolbachia-infected Ae. aegypti. We measured the effect of wMel infection on larval development time at different densities, quantified levels of the key energetic reserve metabolite glycogen in individual adult females, and then performed geometric morphometric analyses to assess the impact of Wolbachia infection and larval density on adult size and wing shape. Our results offer new insight into the role that host genetic background plays in Wolbachia-infected mosquitoes.
Mosquito colony maintenance
All Ae. aegypti mosquitoes used in our experiments were maintained in a climate controlled insectary under previously described conditions . All experiments involved two previously described mosquito lines; the first (wMel_Br) was generated by introducing the wMel Wolbachia strain into a Brazilian genetic background. The second line (wMel_BrTET) was cured of its wMel infection by tetracycline treatment .
Larval rearing and dietary conditions
In order to analyze the effect of crowding on phenotypic traits of wMel-infected Ae. aegypti, we reared larvae under three different larval densities: either 10 (low density condition), 50 (intermediate density condition) or 250 (high density condition) larvae per tray. Eggs from both wMel_Br and wMel_BrTET lineages were hatched synchronously in separate trays containing filtered, dechlorinated water for 1 h without any source of food. After eclosion, first instar larvae were separated into small, black flowerpots (12 × 8.8 × 10 cm) containing 150 ml of filtered, dechlorinated water, according to their respective larval density conditions. In order to obtain around 60–100 adult mosquitoes per treatment per lineage for phenotypic traits, we raised larvae in 40, 10 or 5 trays, for the low, intermediate and the most crowded conditions, respectively. For larval development time, we recorded pupae formation in all trays in order to analyze similar number among treatments, i.e. around 350 individuals per biological replicate. For adult phenotypic traits, we randomly selected a specific number of individuals within a pool made out of each treatment (range of individuals tested is described in each specific experimental section below). Relevant to state, is the fact that we did our analysis on the larval level rather than the tray level, thus ignoring the issue of variation between trays, i.e. we compared a similar sample size.
Larvae were fed Tetramin Tropical Fish Food flakes (Tetra) as a food source. The highest level of food provided was 0.25 mg of food per larvae, each day, as previously described . All trays received 2.5 mg of food per day, with the low density condition functioning as the control treatment, with these larvae receiving the highest amount of food. The other two conditions represented increasing stress due to larval competition for food, with the intermediate and high density treatments receiving 5 and 25 times less food than the low density condition, respectively. We measured larval development for each density condition by recording the time from egg hatch to pupation. Trays were examined every 24 h, where pupae were removed and sexed by visual analysis of their terminalia using a stereomicroscope.
To determine if Wolbachia infection and larval density affected mosquito morphometric traits, we randomly collected a total of 455 adult female mosquitoes within 24 h of emergence. These samples were collected across all treatments and from two independent biological replicates (60–100 mosquitoes measured per treatment) although not all individuals were analyzed. These mosquitoes were stored at room temperature in 70 % ethanol for later analysis, where we randomly selected approximately the same number of female mosquitoes as stated above, separately, from each biological replicate. The right wings from all mosquitoes were detached and the scales were removed manually with a paintbrush, as scales can interfere with landmark data acquisition during wing measurement . Wings were mounted on microscope slides and photographed at 30.2× magnification using a Carl Zeiss AxioCam MRc camera coupled to Carl Zeiss Stemi SV 6 stereomicroscope (Zeiss, Oberkochen, Germany) using AxioVision version 22.214.171.124 image capturing software (Zeiss).
In order to account for the potential effect of allometry (the influence of size in shape) in our analyses, we performed a multivariate regression of Procrustes coordinates against centroid size in all conditions analyzed. From the shape variables generated, wing shape differences between individuals were analyzed based on the Mahalanobis distances (MD), a method for measuring how similar a specific set of conditions is to an ideal set of conditions (defined as the square root of the distances squared between the superimposed individual to the mean shape that are standardized by the covariance matrix of the distance variables), coupled to a permutation test with 10,000 randomizations , where smaller MD numbers between two groups means that they are more closely related to each other.
User measurement error in data acquisition was taken into account through repeatability test procedures, described by the Pearson correlation coefficient between two measurements. This was performed by double landmarking 30 randomly chosen samples from each group using statistical tests previously described . All repeatability values were above 0.99 indicating that data acquisition was accurate.
To determine if Wolbachia infection and larval density affected adult reserves, the total glycogen content of mosquitoes was quantified as previously described . Briefly, a total of 300 female Ae. aegypti adult mosquitoes were collected from the different larval rearing conditions (c.40–60 mosquitoes per treatment; these mosquitoes belonged to a different cohort to that used for previous analyses of morphometric traits and development time), within 24 h of emergence and without the opportunity to feed, thus maintaining the nutritional reserves that were carried over to adulthood after pupation. These mosquitoes were homogenized individually in 200 mM sodium acetate (pH 4.8). The supernatant was removed and incubated with 1 U/ml of α-amyloglucosidase (Sigma, St Louis, MO, USA) at 40 °C for 4 h, in triplicate for each sample. This reaction produced glucose as a final product, and this was used in a colorimetric glucose quantification assay using Glucox 500 (Doles, Goiânia, Brazil). Samples were incubated for 30 min at 37 °C and then quantified using a spectrophotometer at 510 nm (Shimadzu UV-1240). As a control for our reactions, we used free glucose samples obtained from reactions without α-amyloglucosidase where the resulting values were subtracted from the experimental values in order to obtain the final amount of glycogen (i.e. reaction with enzyme - reaction without enzyme = actual value of glycogen). The glycogen content was determined using a standard curve run in parallel. Total protein levels in each mosquito were quantified using the Bradford assay. Bovine serum albumin was used to make standard curves with known concentrations of protein, which we ran in parallel with the experimental samples . The glycogen content in each mosquito was normalized to total protein content, with the final glycogen data presented as a ratio of total glycogen to total protein.
Statistics and data analysis
All assays were repeated twice on biologically independent groups of mosquitoes. The development time data were first combined for all trays in each treatment and log-probit transformed to obtain the rate of pupae formation. These data were analyzed by regression with a 95 % confidence interval in order to compare the development dynamics between groups, where dynamics refers to difference in the rate of pupation over time (peak pupae formation time). Independently, we analyzed the median development time from both groups as well as through pairwise comparisons (described below). The effects of crowding conditions and Wolbachia infection on development time, centroid size and glycogen levels were considered using generalized linear models of regression (SPSS V17, IBM). Pairwise effects were compared using Mann-Whitney U tests followed by multiple test correction using a false discovery rate of 0.05. For non-parametric data, Kruskal-Wallis one-way analysis of variance, followed by pairwise comparison using Dunn’s tests. Plots were made using GraphPad Prism version 6.0 g for Macintosh (GraphPad Software, San Diego, CA, USA). The morphometric statistical analyses of wing centroid size and shape (Generalized Procrustes and Canonical Variate Analysis) were performed using TpsUtil 1.60 and TpsRelw 1.53 (James Rohlf) and MorphoJ 1.06d (Flywings).
Analysis of the effects of Wolbachia infection and larval density on pupal development time status using generalized linear models
Female development data
Male development data
Wolbachia × Crowding
Wolbachia × Crowding
In males, the development dynamics were different for infected and uninfected individuals in the low density condition (F (0.05) = 5.899, P = 0.0008) (Fig. 2d), however, in general, both took a median of 7 days to develop (Additional file 1: Figure S1b). Developmental dynamics also differed for the intermediate density (F (0.05) = 2.778, P = 0.0457) (Fig. 2e), where there was a median pupation time of 9 and 9.5 days for infected and uninfected individuals, respectively (Additional file 1: Figure S1b). For the high-density condition, there was no difference in development dynamics (F (0.05) = 2.072, P = 0.1054), but the median pupation time was 23 days for infected and 26 days for uninfected individuals (Additional file 1: Figure S1b).
Analysis of the effect of Wolbachia infection and larval density on centroid size using generalized linear models
Centroid size data
Wolbachia × Crowding
Mahalanobis distances (MD) for different crowding conditions without allometry
The MD between different larval densities for wMel_Br and wMel_BrTET is depicted in Table 3. As expected, the MD showed that there was a pattern of proximity between larval densities as the crowding condition increased for both infected and uninfected groups. Infected individuals from the low larval density were closer to those reared at intermediate larval density (1.53) than to those from the high larval density condition (2.71). However, for uninfected mosquitoes, this pattern of distribution was reversed, with the lower larval density being more distant from the intermediate condition (2.80), which then was closely related to the high density (1.96). Overall, there was a greater difference between the highest and lowest densities for uninfected individuals (3.34) than for infected mosquitoes (2.71).
In an exploratory Principal Component Analysis (PCA), we found that a total of ten variables (dimensions) were required in order to explain 90.8 % of the variability associated with wing shape in our data. A large proportion of the total observed shape variation (62.3 %) was associated with the first three principal components (PC) based on the covariance matrix, which explained the majority of the observed differences between groups. PC1 (34.1 %) accounted for more variation than PC2 and PC3 combined (28.1 %). Deformation grids produced with the thin-plate spline showing wing shape temporal variations are depicted in Fig. 4c. For PC1, most of the variation can be explained by landmark 2 changing in the opposite direction to landmark 17 and 18. In PC2 the variation in wing shape can be attributed to landmark 2 moving into the opposite direction to the landmarks 12–16, whereas for PC3, it can be attributed to landmarks 12–16 moving towards the edge of the wing while landmarks 17–18 were positioned further towards the center.
Analysis of the effects of Wolbachia infection and larval density on adult glycogen levels using generalized linear models
Wolbachia × Crowding
We have examined how wMel strain of Wolbachia affects Ae. aegypti larval development under different crowding conditions, and how this in turn alters adult mosquito morphology and energetic reserves (measured as glycogen content). Our results show that wMel-infected larvae developed faster at higher densities for both males and females, with no effect under less crowded conditions in females. Infection also led to reduced body size at both high and low, but not intermediate density. Variation in wing shape was also observed between infected and uninfected mosquitoes due to larval density. Wolbachia infection also served to maintain stable glycogen levels at higher larval densities, whereas glycogen levels in uninfected mosquitoes were decreased.
Our observations of larval development dynamics corroborated previous reports [54, 55] showing a gradual increase in the time needed for both infected and uninfected larvae to reach pupal stage as the crowding conditions increased from low to high density.
In general, as larval density increased for both males and females, there were differences in the median development time due to infection at the intermediate and high densities, with infected individuals developing sooner than their uninfected counterparts. This result highlights the beneficial effect of Wolbachia, which acts as a mutualistic agent as conditions worsen and decreases development time for infected larvae. This could make the spread of the bacterium in the field much easier in places where field releases are taking place in Brazil , where a faster development rate makes juvenile individuals less susceptible to predators , a fact that can directly influence the population size in a given area through time.
These results differed from what was previously observed  on natural Wolbachia infections in Ae. albopictus where detrimental effects of infection on both males and females were observed as crowding increased. It is important to note that such differences may have occurred due to differences in strain biology or the fact that wMel in Ae. aegypti is a transinfection . In an Australian genetic background, larval density had no effect on development time in cohorts of wMel-infected larvae, but there was a delayed development associated with infection when uninfected and infected larvae were reared together . The conditions presented by these authors differ from what was tested here. Although they analyzed a broader spectrum of larval densities, ranging from 50 to 800 larvae per tray, this only represented a 16-fold increase in population density whereas here, we considered a 25-fold increase, which seems to be reasonable, considering the larval density of Ae. aegypti in the field [60, 61]. This highlights the importance of host genotype in host-symbiont interactions  and also the importance of taking into account different methodologies.
Previous studies have shown that Wolbachia infection can produce divergent phenotypic effects depending on the sex of the host, where Wolbachia-infected males of Ae. albopictus under high larval densities displayed a delayed development time compared to females [55, 58]. This is not unexpected given that Wolbachia biology is geared towards maternal transmission, with female hosts experiencing heavy selection pressure and male hosts left as little more than genetic dead-ends [63, 64]. As such, particular importance has to be given to sex-based analyses in Wolbachia-infected insects. Although females and males compete in the same environment, the latter given their reduced body mass, develop more quickly and utilize a fair amount of available resources, which indirectly influences the rate of population growth, potentially at the expense of the more epidemiologically important females .
When considering morphometric traits, we first analyzed wing size as a proxy of total body size . There is evidence that greater female body size correlates to an increased ability to transmit pathogens , as well as increased fecundity, flight capacity, host-seeking in order to blood feed  and also resistance to desiccation in the field . Previous work has shown that there is an inversely proportional correlation between wing size and larval density in Ae. aegypti, even when infected with Wolbachia  and this was similar to our results. We also observed that wMel-infection was an important factor affecting wing size, with infected individuals displaying reduced body size for both low and high larval density conditions. This contrasts with what was observed in Australian Ae. aegypti mosquitoes infected with the same Wolbachia strain, where there was a beneficial effect of infection that resulted in larger wings . Here, a large proportion of the differences between groups of infected and uninfected individuals could be attributed to a group of female mosquitoes displaying extreme values for body size (females that were either very small or very big), rather than the values around the actual median of the group.
Interestingly, we observed that wing size decreased to a similar degree for infected and uninfected mosquitoes as larval density increased (by 23.7 and 26 %, respectively): These data suggest that any density × Wolbachia effect affecting wing size and thus adult size, is likely to be only a minor effect. However, we did not compare this effect using a mixed cohort of Wolbachia-infected and uninfected individuals competing in the same tray. From the perspective of releasing wMel-infected Brazilian females into the field for mosquito control, these results do not suggest that wMel infection will serve as a significant disadvantage during competition against uninfected mosquitoes, as also shown elsewhere .
Variation in insect wing morphological traits such as the overall shape can affect dispersal  and the ability to find resources , as flight performance depends on thoracic mass associated with flight muscles, wing asymmetry and the ability to spatially and temporarily make use of variations in the prevailing winds, amongst other factors . As such, wing morphology together with wing size are important factors to be taken into consideration when planning to rear competitive laboratory mosquitoes for future releases in the field . Our results indicate that Wolbachia infection directly influences the wing shape of its host. Larval density also accounted for a degree of shape variation between treatments, with a pattern of proximity between groups as larval density increased (individuals from the low density were more closely related to those from the intermediate density, which were closer to females from the higher density treatment). Other studies support our observation that external factors can influence wing shape, as this variation was observed when comparing field mosquitoes against laboratory-reared mosquitoes. It was also observed when food availability was manipulated and also as an indicator of stress [70–72]. In our data, much of the variance in wing shape was associated with landmark 2 (costal notch) when analyzing the two main principal components that accounted for wing shape variation (PC1 and PC2). However, it is difficult to precisely say if this particular landmark can have direct implications in the fitness of field mosquitoes (given that this is a parameter influenced not only by just one variable) even though it was influenced by changes in food availability.
One important factor relating to wing size and its influence on host biology is the conditions in which the experiments were conducted (extrinsic factors), given the multifactorial nature of this trait. Some studies involved variation in temperature , relative humidity  and food availability [25, 71], while we considered different densities and food availability. Other authors analyzed the influence of body size on fecundity in Ae. aegypti  and other species . Fundamentally, all these studies will provide a different outcome on how the conditions tested influenced wing size. However, all small mosquitoes are not identical. Lack of nutrients and high temperatures can both produce smaller mosquitoes, however, these mosquitoes are likely to have very different metabolic and nutritional profiles. As such, we reiterate the importance of a multifactorial approach, as presented here, to compare the overall fitness of an individual, instead of being based only in a single factor, i.e. size [74, 75].
Insects in general are constantly expending energy; when they are not ingesting food, they are using their energy stocks in order to sustain many important bodily functions . Triglycerides and glycogen are the two main forms of energy reserves that are stored in animal adipocyte cells. These are responsible for controlling the synthesis and utilization of energy reserves, but are also involved in the synthesis of most of the hemolymph proteins and circulating metabolites . Previous work has shown that infection can increase levels of glycogen, as seen for Plasmodium-infected mosquitoes . We saw that the wMel-infected mosquitoes maintained stable levels of glycogen as larval density conditions became increasingly crowded, whereas uninfected mosquitoes experienced decreased levels of this energetic reserve as crowding stress increased. Higher glycogen levels under stressed conditions likely represents a fitness benefit due to wMel infection, as seen elsewhere, with Wolbachia becoming a metabolic provisioning agent under conditions of high dietary iron stress in infected Drosophila melanogaster  whereas other Wolbachia strains display a range of nutritional mutualisms [79, 80]. Together, these traits could assist competitiveness in the field, especially given that it has already been demonstrated that in female Ae. aegypti mosquitoes, feeding on carbohydrates providing higher glycogen reserves, allows extended flight activities and is essential for survival in the absence of blood meals . Interestingly, under high larval density, infected individuals presented reduced body size but higher glycogen levels. Previous studies with Anopheles gambiae have demonstrated that glycogen determines male mating success in a swarm, with the ability to initiate and sustain swarming being positively associated with carbohydrates reserves [10, 81] whereas for other species, there is a correlation between body size and mating competitiveness ; however this correlation is not always positive and there are cases where size is not a determinant of mating capacity .
Implications for mosquito control
With field releases of the wMel-infected Ae. aegypti mosquitoes just starting in Brazil, in areas that display great variability in their physical, human/mosquito population and demographic characteristics, it is critical to assess the effects of this bacterial strain in mosquitoes with a Brazilian genotype. Different Wolbachia strains can differently affect host physiology and this variability can be influenced by environmental conditions as well as host genetic background [27, 37, 83]. In the Brazilian genetic background, wMel infection causes almost complete CI, has a high rate of maternal transmission and has no detrimental effect on fecundity, when females where reared under optimal conditions in a laboratory setting . However, field conditions do not always provide the best-case scenario for insect development, as crowding and other factors can detrimentally affect fitness [19, 84–87]. As such, analysis of how Wolbachia affects mosquitoes reared under different density-dependent regimes will provide information that is vital to the planning stages of field releases of Wolbachia-infected mosquitoes in the field, along with an additional measure of infected mosquito competitiveness.
Our results indicate that the effects of wMel on its host are variable depending on larval density conditions. Perhaps the most important benefit is the comparatively decreased development time, as insects that develop faster tend to have an advantage in the field in terms of exposure to larval predation , mating competition  and food availability . Infected mosquitoes with faster development times also had elevated levels of glycogen, which would likely prove to be of great benefit in the field and could provide Wolbachia-infected mosquitoes with the fuel to fly further or longer and this could potentially aid in host or oviposition site location , as well as invasion. Both of these traits would likely improve adult competitiveness in the field.
In combination with the previously described physiological effects  it would be tempting to suggest that wMel has an entirely beneficial effect in the Brazilian background. However, our data indicate that there are infection-dependent decreases in wing size and alterations in shape. These changes could affect flight and behaviours such as host-seeking which depends on flight . The nature and extent of these effects are unclear and should be investigated further.
Nutritional benefits due to Wolbachia infection, as observed in our glycogen results are fairly common, but have not typically been associated with Wolbachia infection in Ae. aegypti. In terms of metabolic pathways, the wMel strain is genetically limited, suggesting a reliance on host resources . In Ae. aegypti, wMel has been previously shown to deplete host cholesterol levels . Our results indicate that wMel has a complex metabolic relationship with its host, varying between parasitic and mutualistic, with this being the first example of a nutritional mutualism for this strain in mosquitoes [90, 91]. The provision of useful resources such as glycogen may also contribute to the faster development rates than we observed, especially given that similar effects on glycogen have been observed with the native wFlu strain in Ae. fluviatilis eggs .
Heavy fitness costs due to detrimental effects on the biology of Wolbachia-infected mosquitoes could hinder the spread of the bacterium in the field . Our results indicate that such effects are not likely to occur for wMel infections in Brazilian Ae. aegypti, particularly under conditions of high larval density and competition. Faster development rates and higher glycogen levels in infected in comparison to uninfected mosquitoes, may lead to increased competitiveness associated with infection and such beneficial fitness effects are desirable for the successful deployment of Wolbachia in the field [10, 92]. However these benefits could potentially be balanced by observed decreases in mosquito wing size, suggesting smaller body size and altered wing shape, which could affect flight, although not necessarily detrimentally.
The phenotypic effects we have observed are distinct from those previously observed in Australian mosquitoes, where this bacterium was first deployed as a biocontrol agent [25, 43]. This highlights the importance of investigating the effects of different host genetic backgrounds on Wolbachia infection, particularly in advance of releases of mosquitoes for vector control purposes. Our data provide additional information that can be used to inform release planning and design in Wolbachia population replacement mosquito control strategies in Brazil, where mosquitoes are being released into areas with large local mosquito populations  and frequent dengue transmission . Specifically, they will be useful to developing an accurate picture of the fitness effects of Wolbachia infection, which is a requirement in the mathematical modeling used to predict whether Wolbachia can invade specific localities [56, 94]. These data could enable researchers to more accurately estimate the factors involved in Wolbachia release dynamics, such as the overall release period, number of mosquitoes releases and the distribution of releases within a field site, given our improved understanding of the competitiveness of wMel-infected Brazilian Ae. aegypti.
Our results suggest that the wMel strain of Wolbachia can actually have beneficial effects on host physiology under certain conditions. We observed beneficial effects through decreased time to pupation and increased glycogen content under different larval crowding scenarios; however these were offset by a mild decrease in body size of female mosquitoes and distinct wing shapes associated with infection. These results, in combination with previous work, indicate that the wMel strain of Wolbachia can likely be used successfully in open field releases in Brazil, as part of the Eliminate Dengue Program.
We are grateful to all members of the Mosquitos Vetores: Endossimbiontes e Interação Patógeno Vetor (MV - CPqRR/FIOCRUZ), especially Dr. Eric Caragata for critically reviewing the paper. Members of Biologia Molecular e Imunologia da Malária (BMIM - CPqRR/FIOCRUZ), Laboratório de Transmissores de Hematozoários (LATHEMA - FIOCRUZ) and Laboratório de Sanidade Animal, Laboratório de Química e Função de Proteínas e Peptídeos e Unidade de Experimentação Animal (UENF). We also extend our thanks to the Brazilian and Australian coordinators of the Eliminate Dengue program, in particular to Prof Scott L. O’Neill for donating the original wMel line.
This work was partially supported by FAPEMIG, CNPq, the Brazilian Ministry of Health (DECIT/SVS) as well as a grant to Monash University from the Foundation for the National Institutes of Health through the Vector-Based Transmission of Control: Discovery Research (VCTR) program of the Grand Challenges in Global Health Initiatives of the Bill and Melinda Gates Foundation.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496:504–7.View ArticlePubMedPubMed CentralGoogle Scholar
- World Health Organization. WHO | Global Strategy for dengue prevention and control, 2012–2020. Geneva: World Health Organization; 2012. http://www.who.int/denguecontrol/9789241504034/en/index.html (accessed Nov. 9, 2015).Google Scholar
- Hardy JL, Houk EJ, Kramer LD, Reeves WC. Intrinsic factors affecting vector competence of mosquitoes for arboviruses. Annu Rev Entomol. 1983;28:229–62.View ArticlePubMedGoogle Scholar
- Tjaden NB, Thomas SM, Fischer D, Beierkuhnlein C. Extrinsic incubation period of dengue: knowledge, backlog, and applications of temperature dependence. PLoS Negl Trop Dis. 2013;7:e2207.View ArticlePubMedPubMed CentralGoogle Scholar
- Beerntsen BT, James AA, Christensen BM. Genetics of mosquito vector competence. Microbiol Mol Biol Rev. 2000;64:115–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Bara J, Rapti Z, Cáceres CE, Muturi EJ. Effect of larval competition on extrinsic incubation period and vectorial capacity of Aedes albopictus for dengue virus. PLoS One. 2015;10:e0126703.View ArticlePubMedPubMed CentralGoogle Scholar
- Dye C. Vectorial capacity: must we measure all its components? Parasitol Today. 1986;2:203–9.View ArticlePubMedGoogle Scholar
- Couret J, Dotson E, Benedict MQ. Temperature, larval diet, and density effects on development rate and survival of Aedes aegypti (Diptera: Culicidae). PLoS One. 2014;9:e87468.View ArticlePubMedPubMed CentralGoogle Scholar
- Tun-Lin W, Burkot TR, Kay BH. Effects of temperature and larval diet on development rates and survival of the dengue vector Aedes aegypti in north Queensland, Australia. Med Vet Entomol. 2000;14:31–7.View ArticlePubMedGoogle Scholar
- Briegel H, Knüsel I, Timmermann SE. Aedes aegypti: size, reserves, survival, and flight potential. J Vector Ecol. 2001;26:21–31.PubMedGoogle Scholar
- Arrese EL, Soulages JL. Insect fat body: energy, metabolism, and regulation. Annu Rev Entomol. 2010;55:207–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Kaufmann C, Brown MR. Regulation of carbohydrate metabolism and flight performance by a hypertrehalosaemic hormone in the mosquito Anopheles gambiae. J Insect Physiol. 2008;54:367–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Takken W, Verhulst NO. Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2011;58:120928130709004.Google Scholar
- Benjamin SN, Bradshaw WE. Body size and flight activity effects on male reproductive success in the pitcherplant mosquito (Diptera: Culicidae). Ann Entomol Soc Am. 1994;87:331–6.View ArticleGoogle Scholar
- Ponlawat A, Harrington LC. Factors associated with male mating success of the dengue vector mosquito, Aedes aegypti. Am J Trop Med Hyg. 2009;80:395–400.PubMedGoogle Scholar
- Alto BW, Reiskind MH, Lounibos LP. Size alters susceptibility of vectors to dengue virus infection and dissemination. Am J Trop Med Hyg. 2008;79:688–95.PubMedPubMed CentralGoogle Scholar
- Price DP, Schilkey FD, Ulanov A, Hansen IA. Small mosquitoes, large implications: crowding and starvation affects gene expression and nutrient accumulation in Aedes aegypti. Parasit Vectors. 2015;8:1–14.View ArticleGoogle Scholar
- Alto BW, Lounibos LP, Mores CN, Reiskind MH. Larval competition alters susceptibility of adult Aedes mosquitoes to dengue infection. Proc Biol Sci. 2008;275:463–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Mueller LD. Theoretical and empirical examination of density-dependent selection. Annu Rev Ecol Syst. 1997;28:269–88.Google Scholar
- Moore CG, Whitacre DM. Competition in mosquitoes. Production of Aedes aegypti larval growth retardant at various densities and nutrition levels. Ann Entomol Soc Am. 1972;65:915–8.View ArticleGoogle Scholar
- 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:211–26.PubMedPubMed CentralGoogle Scholar
- Zug R, Hammerstein P. Still a host of hosts for Wolbachia: analysis of recent data suggests that 40 % of terrestrial arthropod species are infected. PLoS One. 2012;7:e38544.View ArticlePubMedPubMed CentralGoogle Scholar
- McGraw EA, O’Neill SL. Beyond insecticides: new thinking on an ancient problem. Nat Rev Microbiol. 2013;11:181–93.View ArticlePubMedGoogle Scholar
- Fleury F, Vavre F, Ris N, Fouillet P, Boulétreau M. Physiological cost induced by the maternally-transmitted endosymbiont Wolbachia in the Drosophila parasitoid Leptopilina heterotoma. Parasitology. 2000;121(Pt 5):493–500.View ArticlePubMedGoogle Scholar
- Ross PA, Endersby NM, Yeap HL, Hoffmann AA. Larval competition extends developmental time and decreases adult size of wMelPop Wolbachia-infected Aedes aegypti. Am J Trop Med Hyg. 2014;91:198–205.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhukova MV, Kiseleva E. The virulent Wolbachia strain wMelPop increases the frequency of apoptosis in the female germline cells of Drosophila melanogaster. BMC Microbiol. 2012;12 Suppl 1:S15.View ArticlePubMedPubMed CentralGoogle Scholar
- Baton LA, Pacidônio EC, Gonçalves DDS, Moreira LA. wFlu: characterization and evaluation of a native Wolbachia from the mosquito Aedes fluviatilis as a potential vector control agent. PLoS One. 2013;8:e59619.View ArticlePubMedPubMed CentralGoogle Scholar
- Dobson SL, Rattanadechakul W, Marsland EJ. Fitness advantage and cytoplasmic incompatibility in Wolbachia single- and superinfected Aedes albopictus. Heredity. 2004;93:135–42.View ArticlePubMedGoogle Scholar
- Dean MD. A Wolbachia-associated fitness benefit depends on genetic background in Drosophila simulans. Proc Biol Sci. 2006;273:1415–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Kittayapong P, Baisley KJ, Baimai V, O’Neill SL. Distribution and diversity of Wolbachia infections in southeast asian mosquitoes (Diptera: Culicidae). J Med Entomol. 2000;37:340–5.View ArticlePubMedGoogle Scholar
- McMeniman CJ, Lane RV, Cass BN, Fong AWC, Sidhu M, Wang Y-F, et al. Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science. 2009;323:141–4.View ArticlePubMedGoogle Scholar
- Walker T, Johnson PH, Moreira LA, Iturbe-Ormaetxe I, Frentiu FD, McMeniman CJ, et al. The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations. Nature. 2011;476:450–3.View ArticlePubMedGoogle Scholar
- Moreira LA, Iturbe-Ormaetxe I, Jeffery JA, Lu G, Pyke AT, Hedges LM, et al. A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya, and Plasmodium. Cell. 2009;139:1268–78.View ArticlePubMedGoogle Scholar
- Caragata EP, Dutra HLC, Moreira LA. Exploiting intimate relationships: controlling mosquito-transmitted disease with Wolbachia. Trends Parasitol. 2016;32:207–18.View ArticlePubMedGoogle Scholar
- Nguyen TH, Nguyen HL, Nguyen TY, Vu SN, Tran ND, Le TN, et al. Field evaluation of the establishment potential of wMelpop Wolbachia in Australia and Vietnam for dengue control. Parasit Vectors. 2015;8:563.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoffmann AA, Goundar AA, Long SA, Johnson PH, Ritchie SA. Invasion of Wolbachia at the residential block level is associated with local abundance of Stegomyia aegypti, yellow fever mosquito, populations and property attributes. Med Vet Entomol. 2014;28 Suppl 1:90–7.View ArticlePubMedGoogle Scholar
- Dutra HLC, dos Santos LMB, Caragata EP, Silva JBL, Villela DAM, Maciel-de-Freitas R, et al. From lab to field: the influence of urban landscapes on the invasive potential of Wolbachia in Brazilian Aedes aegypti mosquitoes. PLoS Negl Trop Dis. 2015;9:e0003689.View ArticlePubMedPubMed CentralGoogle Scholar
- Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008;6:741–51.View ArticlePubMedGoogle Scholar
- Yeap HL, Mee P, Walker T, Weeks AR, O’Neill SL, Johnson P, et al. Dynamics of the “popcorn” Wolbachia infection in outbred Aedes aegypti informs prospects for mosquito vector control. Genetics. 2011;187:583–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Hancock P, Sinkins S, Godfray H. Population dynamic models of the spread of Wolbachia. Am Nat. 2011;177:323–33.View ArticlePubMedGoogle Scholar
- Da Rocha Fernandes M, Martins R, Pessoa Costa E, Casagrande Pacidônio E, Araujo De Abreu L, Da Silva Vaz I, et al. The modulation of the symbiont/host interaction between Wolbachia pipientis and Aedes fluviatilis embryos by glycogen metabolism. PLoS One. 2014;9(6):e98966.Google Scholar
- Ross PA, Endersby NM, Hoffmann AA. Costs of three Wolbachia infections on the survival of Aedes aegypti larvae under starvation conditions. PLoS Negl Trop Dis. 2016;10:e0004320.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoffmann AA, Montgomery BL, Popovici J, Iturbe-Ormaetxe I, Johnson PH, Muzzi F, et al. Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission. Nature. 2011;476:454–7.View ArticlePubMedGoogle Scholar
- Reynolds KT, Thomson LJ, Hoffmann AA. The effects of host age, host nuclear background and temperature on phenotypic effects of the virulent Wolbachia strain popcorn in Drosophila melanogaster. Genetics. 2003;164:1027–34.PubMedPubMed CentralGoogle Scholar
- Lorenz C, Suesdek L. Short report: evaluation of chemical preparation on insect wing shape for geometric morphometrics. Am J Trop Med Hyg. 2013;89:928–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Lorenz C, Marques TC, Sallum MAM, Suesdek L. Morphometrical diagnosis of the malaria vectors Anopheles cruzii, An. homunculus and An. bellator. Parasit Vectors. 2012;5:257.View ArticlePubMedPubMed CentralGoogle Scholar
- Siegel JP, Novak RJ, Lampman RL, Steinly BA. Statistical appraisal of the weight-wing length relationship of mosquitoes. J Med Entomol. 1992;29:711–4.View ArticlePubMedGoogle Scholar
- Stoyan D, Bookstein FL. Morphometric tools for landmark data. Geometry and biology. Volume 35. New York: Cambridge University Press; 1993. p. 512.Google Scholar
- Rohlf FJ. Shape statistics: procrustes superimpositions and tangent spaces. J Classif. 2014;16:197–223.View ArticleGoogle Scholar
- Klingenberg CP, Mcintyre GS. Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with procrustes methods. Evolution. 1998;52:1363–75.View ArticleGoogle Scholar
- Arnqvist G, Martensson T. Measurement error in geometric morphometrics: empirical strategies to assess and reduce its impact on measures of shape. Acta Zool Acad Sci. 1998;44(1):73-96.Google Scholar
- Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.View ArticlePubMedGoogle Scholar
- Gayon J. History of the concept of allometry. Am Zool. 2000;40:748–58.Google Scholar
- Serpa L, Kakitani I, Voltolini J. Competition between Aedes aegypti and Aedes albopictus larvae in the laboratory. Rev da Soc Bras. 2008;41:479–84.View ArticleGoogle Scholar
- Gavotte L, Mercer D. Wolbachia infection and resource competition effects on immature Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2009;46:451–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hancock PA, Godfray HCJ. Modelling the spread of Wolbachia in spatially heterogeneous environments. J R Soc Interface. 2012;9(76):3045–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Alto BW, Malicoate J, Elliott SM, Taylor J. Demographic consequences of predators on prey: trait and density mediated effects on mosquito larvae in containers. PLoS One. 2012;7:e45785.View ArticlePubMedPubMed CentralGoogle Scholar
- Gavotte L, Mercer SR, Stoeckle JJ, Dobson SL. Costs and benefits of Wolbachia infection in immature Aedes albopictus depend upon sex and competition level. J Invertebr Pathol. 2010;105:341–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Islam M, Dobson S. Wolbachia effects on Aedes albopictus (Diptera: Culicidae) immature survivorship and development. J Med Entomol. 2006;43:689–95.View ArticlePubMedGoogle Scholar
- Peres RC, Rego R, Maciel-de-Freitas R. The use of the Premise Condition Index (PCI) to provide guidelines for Aedes aegypti surveys. J Vector Ecol. 2013;38:190–2.View ArticlePubMedGoogle Scholar
- Codeço CT, Lima AWS, Araújo SC, Lima JBP, Maciel-de-Freitas R, Honório NA, et al. Surveillance of Aedes aegypti: comparison of house index with four alternative traps. PLoS Negl Trop Dis. 2015;9:e0003475.View ArticlePubMedPubMed CentralGoogle Scholar
- Kondo N, Shimada M, Fukatsu T. Infection density of Wolbachia endosymbiont affected by co-infection and host genotype. Biol Lett. 2005;1:488–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y-K, Ding X-L, Rong X, Hong X-Y. How do hosts react to endosymbionts? A new insight into the molecular mechanisms underlying the Wolbachia-host association. Insect Mol Biol. 2015;24:1–12.View ArticlePubMedGoogle Scholar
- Correa CC, Ballard JWO. What can symbiont titres tell us about co-evolution of Wolbachia and their host? J Invertebr Pathol. 2014;118:20–7.View ArticlePubMedGoogle Scholar
- Armbruster P, Hutchinson R. Pupal mass and wing length as indicators of fecundity in Aedes albopictus and Aedes geniculatus (Diptera: Culicidae). J Med Entomol. 2002;39:699–704.View ArticlePubMedGoogle Scholar
- Mogi M, Miyagi I, Abadi K, Syafruddin. Inter- and intraspecific variation in resistance to desiccation by adult Aedes (Stegomyia) spp. (Diptera: Culicidae) from Indonesia. J Med Entomol. 1996;33:53–7.View ArticlePubMedGoogle Scholar
- Hassall C, Thompson DJ, Harvey IF. Latitudinal variation in morphology in two sympatric damselfly species with contrasting range dynamics (Odonata: Coenagrionidae). Eur J Entomol. 2008;105:939–44.View ArticleGoogle Scholar
- Hoffmann AA, Ratna E, Sgrò CM, Barton M, Blacket M, Hallas R, et al. Antagonistic selection between adult thorax and wing size in field released Drosophila melanogaster independent of thermal conditions. J Evol Biol. 2007;20:2219–27.View ArticlePubMedGoogle Scholar
- Breuker CJ, Brakefield PM, Gibbs M. The association between wing morphology and dispersal is sex-specific in the glanville fritillary butterfly Melitae cinxia (Lepidoptera: Nymphalidae). Eur J Entomol. 2007;104:445–52.View ArticleGoogle Scholar
- Yeap HL, Endersby NM, Johnson PH, Ritchie SA, Hoffmann AA. Body size and wing shape measurements as quality indicators of Aedes aegypti mosquitoes destined for field release. Am J Trop Med Hyg. 2013;89:78–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Jirakanjanakit N, Leemingsawat S, Thongrungkiat S, Apiwathnasorn C, Singhaniyom S, Bellec C, et al. Influence of larval density or food variation on the geometry of the wing of Aedes (Stegomyia) aegypti. Trop Med Int Health. 2007;12:1354–60.View ArticlePubMedGoogle Scholar
- Hoffmann AA, Woods RE, Collins E, Wallin K, White A, McKenzie JA. Wing shape versus asymmetry as an indicator of changing environmental conditions in insects. Aust J Entomol. 2005;44:233–43.View ArticleGoogle Scholar
- Briegel H. Metabolic relationship between female body size, reserves, and fecundity of Aedes aegypti. J Insect Physiol. 1990;36:165–72.Google Scholar
- Reiskind MH, Zarrabi AA. Is bigger really bigger? Differential responses to temperature in measures of body size of the mosquito, Aedes albopictus. J Insect Physiol. 2012;58:911–7.View ArticlePubMedGoogle Scholar
- Lanciani CA, Le TM. Effect of temperature on the wing length-body weight relationship in Anopheles quadrimaculatus. J Am Mosq Control Assoc. 1995;11:241–3.PubMedGoogle Scholar
- Caragata EP, Rancès E, O’Neill SL, McGraw EA. Competition for amino acids between Wolbachia and the mosquito host, Aedes aegypti. Microb Ecol. 2014;67:205–18.View ArticlePubMedGoogle Scholar
- Zhao YO, Kurscheid S, Zhang Y, Liu L, Zhang L, Loeliger K, et al. Enhanced survival of plasmodium-infected mosquitoes during starvation. PLoS One. 2012;7:1–8.Google Scholar
- Brownlie JC, Cass BN, Riegler M, Witsenburg JJ, Iturbe-Ormaetxe I, McGraw EA, et al. Evidence for metabolic provisioning by a common invertebrate endosymbiont, Wolbachia pipientis, during periods of nutritional stress. PLoS Pathog. 2009;5:e1000368.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor MJ, Bandi C, Hoerauf A. Wolbachia bacterial endosymbionts of filarial nematodes. Adv Parasitol. 2005;60:245–84.View ArticlePubMedGoogle Scholar
- Kaiser W, Huguet E, Casas J, Commin C, Giron D. Plant green-island phenotype induced by leaf-miners is mediated by bacterial symbionts. Proc Biol Sci. 2010;277:2311–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Huho BJ, Ng’habi KR, Killeen GF, Nkwengulila G, Knols BGJ, Ferguson HM. Nature beats nurture: a case study of the physiological fitness of free-living and laboratory-reared male Anopheles gambiae s.l. J Exp Biol. 2007;210:2939–47.View ArticlePubMedGoogle Scholar
- Charlwood JD, Pinto J, Sousa CA, Ferreira C, Do Rosário VE. Male size does not affect mating success (of Anopheles gambiae in São Tomé). Med Vet Entomol. 2002;16:109–11.View ArticlePubMedGoogle Scholar
- Brelsfoard CL, Dobson SL. Wolbachia effects on host fitness and the influence of male aging on cytoplasmic incompatibility in Aedes polynesiensis (Diptera: Culicidae). J Med Entomol. 2011;48:1008–15.View ArticlePubMedGoogle Scholar
- Walsh RK, Bradley C, Apperson CS, Gould F. An experimental field study of delayed density dependence in natural populations of Aedes albopictus. PLoS One. 2012;7:e35959.View ArticlePubMedPubMed CentralGoogle Scholar
- Yeap HL, Axford JK, Popovici J, Endersby NM, Iturbe-Ormaetxe I, Ritchie SA, et al. Assessing quality of life-shortening Wolbachia-infected Aedes aegypti mosquitoes in the field based on capture rates and morphometric assessments. Parasit Vectors. 2014;7:58.View ArticlePubMedPubMed CentralGoogle Scholar
- Mori A. Effects of larval density and nutrition on some attributes of immature and adult Aedes albopictus. Trop Med. 1979;21:85–103.Google Scholar
- 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:1325–31.View ArticlePubMedGoogle Scholar
- Bedhomme S, Agnew P, Sidobre C, Michalakis Y. Sex-specific reaction norms to intraspecific larval competition in the mosquito Aedes aegypti. J Evol Biol. 2003;16:721–30.View ArticlePubMedGoogle Scholar
- Wu M, Sun LV, Vamathevan J, Riegler M, Deboy R, Brownlie JC, et al. Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: a streamlined genome overrun by mobile genetic elements. PLoS Biol. 2004;2:E69.View ArticlePubMedPubMed CentralGoogle Scholar
- Zientz E, Dandekar T, Gross R. Metabolic interdependence of obligate intracellular bacteria and their insect hosts. Microbiol Mol Biol Rev. 2004;68:745–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Foster J, Ganatra M, Kamal I, Ware J, Makarova K, Ivanova N, et al. The Wolbachia genome of Brugia malayi: endosymbiont evolution within a human pathogenic nematode. PLoS Biol. 2005;3:e121.View ArticlePubMedPubMed CentralGoogle Scholar
- Levi T, Ben-Dov E, Shahi P, Borovsky D, Zaritsky A. Growth and development of Aedes aegypti larvae at limiting food concentrations. Acta Trop Elsevier BV. 2014;133:42–4.View ArticleGoogle Scholar
- Van Panhuis WG, Hyun S, Blaney K, Marques ETA, Coelho GE, Siqueira JB, et al. Risk of dengue for tourists and teams during the World Cup 2014 in Brazil. PLoS Negl Trop Dis. 2014;8(7):e3063.Google Scholar
- Ndii MZ, Hickson RI, Allingham D, Mercer GN. Modelling the transmission dynamics of dengue in the presence of Wolbachia. Math Biosci. 2015;262:157–66.View ArticlePubMedGoogle Scholar