Next Article in Journal
Early Olfactory Environment Influences Antennal Sensitivity and Choice of the Host-Plant Complex in a Parasitoid Wasp
Previous Article in Journal
Diversity and Ecosystem Services of Trichoptera
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Life History Effects Linked to an Advantage for wAu Wolbachia in Drosophila

1
Pest and Adaptation Research Group, School of BioSciences, Bio21 Institute, The University of Melbourne, Melbourne, VIC 3010, Australia
2
Institute of Plant and Environmental Protection, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
3
College of Life & Environmental Science, Shanghai Normal University, Shanghai 200234, China
*
Authors to whom correspondence should be addressed.
Insects 2019, 10(5), 126; https://doi.org/10.3390/insects10050126
Submission received: 8 April 2019 / Revised: 24 April 2019 / Accepted: 30 April 2019 / Published: 2 May 2019

Abstract

:
Wolbachia endosymbiont infections can persist and spread in insect populations without causing apparent effects on reproduction of their insect hosts, but the mechanisms involved are largely unknown. Here, we test for fitness effects of the wAu infection of Drosophila simulans by comparing multiple infected and uninfected polymorphic isofemale lines derived from nature. We show a fitness advantage (higher offspring number) for lines with the wAu Wolbachia infection when breeding on grapes, but only where there was Talaromyces and Penicillium fungal mycelial growth. When breeding on laboratory medium, the wAu infection extended the development time and resulted in larger females with higher fecundity, life history traits, which may increase fitness. A chemical associated with the fungi (ochratoxin A) did not specifically alter the fitness of wAu-infected larvae, which developed slower and emerged with a greater weight regardless of toxin levels. These findings suggest that the fitness benefits of Wolbachia in natural populations may reflect life history changes that are advantageous under particular circumstances, such as when breeding occurs in rotting fruit covered by abundant mycelial growth.

1. Introduction

Endosymbionts can have a variety of fitness effects on their invertebrate hosts, many of which are beneficial to the hosts, particularly through the provision of nutritional resources. For instance, in tsetse flies, Wigglesworthia glossinidia provides nutrients essential for female fertility [1], while in aphids, Buchnera aphidicola synthesizes essential amino acids [2]. Similarly, in ticks, Francisella bacteria synthesize B vitamins, which are lacking in their blood meals [3]. Not all invertebrate endosymbionts necessarily have beneficial effects on their hosts. For instance, aphids carry a pea aphid secondary symbiont (PASS), which suppresses growth and reproduction when Buchnera symbionts are present, although PASS also benefits aphids in the absence of Buchnera [4].
Of the different endosymbionts found in insects and other invertebrates, Wolbachia infections are the most widely studied because of their widespread distribution across invertebrates [5,6]; their diverse effects on host fitness, including various reproductive effects [7]; and their capacity to interfere with viral transmission in invertebrate hosts and thereby impact on human diseases associated with arboviruses [8,9,10,11]. Wolbachia strains can affect the reproduction of hosts in such a way as to enhance their spread and persistence in natural populations. This is particularly associated with cytoplasmic incompatibility (CI), the reduction in embryo and offspring viability when females that lack Wolbachia mate with males that carry Wolbachia, ensuring a fitness advantage of infected females over uninfected females [12] and leading to a rapid spread of Wolbachia in natural populations [13].
Although much of the focus of Wolbachia research has been on infections that spread via CI or that cause other reproductive effects, like feminization and male killing, many Wolbachia infections from natural populations also provide a range of benefits to their hosts [6]. This includes nutritional benefits in bedbugs [14] as well as the mutualistic associations of Wolbachia in filarial nematodes [15]. In Drosophila, Wolbachia may also provide nutritional benefits [16] although fitness effects on hosts often remain unclear and have only been identified indirectly [17,18].
One of the first Wolbachia infections to be described with uncertain host effects was the wAu infection of Drosophila simulans, identified from natural Australian populations [18], where it showed a potential to increase in frequency even though it had no impact on CI or other host reproductive effects [19]. Recent population cage experiments where flies were raised on fermenting fruits indicated a substantial fitness advantage of the wAu infection of around 20% per generation [20]. This infection is particularly interesting because after transfer to Aedes mosquitoes it has the potential to strongly inhibit Zika and dengue virus transmission [8]. If the wAu infection is associated with positive effects on fitness under some situations, this might assist in its spread from a low frequency in native hosts, which can normally be difficult to achieve when Wolbachia strains only influence host reproduction [21,22].
In this paper, we use a novel design to investigate the fitness advantages that might be associated with the wAu infection in Drosophila simulans by undertaking experiments that build on earlier work [19,20], which pointed to fitness advantages of this infection, but without identifying the relevant causes. We initially focused on the fitness of flies carrying wAu breeding in rotting fruit with different levels of Talaromyces and Penicillium fungal growth. We test whether a common compound associated with the predominant fungi (ochratoxin A) might mediate the fitness advantage of wAu in these breeding sites. However, we instead find that wAu-infected larvae had a delayed development time, but larger emergence size regardless of breeding conditions, and we suspect these changes are critical for the infection to persist and spread in sites with abundant fungal growth. Thus, an overall life history impact of this Wolbachia infection that does not affect host reproduction seems to form the basis for a fitness advantage.

2. Methods

2.1. Lines and Infection Status

Multiple lines of D. simulans (70 infected, 70 uninfected) were established for these experiments, each from 20 to 30 offspring of different field-collected females. Each isofemale line was confirmed to be either wAu-infected or uninfected based on a molecular assay (below). The lines were initiated with the offspring of field-caught females collected at the same time around Roleystone near Perth, Western Australia, where the infection frequency was around 60% based on data collected from the current study. Offspring of the field females forming a line were reared at a low density in 150 mL vials with a dead yeast-treacle- agar fly medium [19] for two generations before the three experiments (Table 1) started. Only lines reared on this medium were used in experiments. Isofemale lines were scored again for Wolbachia status after the experiments to confirm they were infected.
Flies were assayed for Wolbachia infection status and strain type using a real-time PCR/high resolution melt (RT-PCR/HRM) method designed to amplify a fragment of the wsp gene following Chelex extraction of DNA from individual flies using wsp primers (wsp_validation_Fwd: 5′-TTGGTTACAAAATGGACGACATCAG-3′; wsp_validation_Rev: 5′-CGAAATAACGAGCTCCAGCATAAAG-3′) [19,23]. These primers produce different melting temperatures of the wsp products. Separate D. simulans specific primers (Dsim RpS6 Fwd: 5′-CCAGATCGCTTCCAAGGAGGCTGCT-3′; DsimRpS6 Rev: 5′-GCCTCCTCGCGCTTGGCCTTAGAT-3′) were used as an internal control to test for the quality of the DNA extraction procedure, confirm species identification, and ensure that the correct PCR conditions were present. Individual flies from wAu-infected and wRi-infected D. simulans mass bred laboratory lines were used as positive and strain type controls for each PCR plate, while individual flies from an uninfected mass bred line and “no template” PCR reactions were used as negative controls. This assay can distinguish wAu from wRi infected individuals as well as uninfected individuals, which is critical because the wRi infection has recently spread through D. simulans populations in eastern Australia where only the wAu infection had previously been detected [19]. The wRi infection has not yet been found in Western Australia [20].

2.2. Comparison of Lines on Grapes with Fungal Growth

Kriesner and Hoffmann [20] showed that wAu had a fitness advantage when flies were held in population cages that were maintained on fruit, and that conditions provided by fruit may have accounted for the substantial fitness advantage of wAu, which resulted in an increase in infection frequency. We tested this by using organic table grapes (Red Globe) not contaminated by pesticides. Pairs of male and female adults (from each of the 23 wAu-infected and 23 uninfected lines, 3 pairs per isofemale line) were set up when D. simulans were five days post eclosion.
Adults were provided with bottles (250 mL) containing grapes, which were cut in half and placed on vermiculite. To accustom the flies to the grapes and conditions, flies were initially held on one grape per bottle (cut in half) for three days. Bottles were placed in a shade house, where they were exposed to fluctuating ambient conditions (temperatures in ranged from 13.9 to 31.6 °C, natural daylight). Adult pairs were then transferred (or “tipped”) to new bottles after three days of conditioning on the initial grapes. There were 6 grape halves (rather than 2) per bottle in these new bottles to ensure the fruit was not limiting for offspring production. Adults were left in these bottles for 3 days and then transferred to a new set of bottles with the same amount of grape resource for an additional three days. Emerging numbers of adults in bottles from this second and third set of bottles (offspring per pair of adults) were recorded after collecting flies 18 days after bottles had been set up (about 3 days after the first flies were observed). Where any male or female parent had died in a bottle (14.1% of cases), bottles were not used in life history assessments.
We suspected that growth of fungi, bacteria, or yeast might be an important factor influencing the fitness of the Wolbachia infection. We therefore compared the offspring number between wAu-infected and uninfected lines, focusing on the amount of fungal growth, which was easy to assess visually through the amount of mycelial mass. Fungus infestation levels on grapes were evaluated on the 18th day after flies were set up, with 3 levels recognized: Level 1, less than 1 piece of 6 pieces of grapes was covered by visible mycelial growth; level 2, 2 to 4 pieces of grapes were covered by visible growth; level 3, 5 to 6 pieces of grapes were covered by visible growth (Supplementary Figure S1). We compared the offspring number between wAu-infected and uninfected lines (infection as a fixed effect) with ANOVAs (IBM SPSS Statistics version 25) that also included fungus level as a fixed effect after testing for normality with a Kolmogorov-Smirnov test. Note that there were at least 9 replicate bottles (and a maximum of 38 bottles) in an infection-fungal growth level category. Because there were interactions between the fungal category and whether or not lines were infected (see Results below), we further considered the effect of the infection at each fungal level with t tests and showed patterns with kernel density (violin) plots, which show how trait values are distributed (R software package vioplot https://cran.r-project.org/web/packages/vioplot/index.html).
To examine fly weight, females and males from the same bottle were separated into glass vials and dried at 50 °C for 24 h in an oven. Weight of females and males from each bottle was obtained using an analytical balance (Sartorius Type 1712, Sartorius AG, Gottingen, Germany) with accuracy up to 0.01 mg. We excluded bottles with fewer than 4 individuals. Average weight of females and males from each bottle (expressed as average weight per fly in a bottle) were compared between wAu-infected and uninfected lines with an ANCOVA, which included adult parental age (successive tips) and fungal growth as fixed effects; emergence number was included as a covariate because adult weight may depend on larval density, although this effect was not significant. Sex ratio (female/(female + male)) was estimated by adding up all the offspring produced by a single pair of flies averaged across the two observation periods. A nested factor due to the isofemale line was included in initial models for this experiment and the other experiments did not appreciably affect the significance of infection effects, which is the focus of the current study rather than nuclear-based fecundity differences.
Morphology of fungi was examined under a microscope to assess species. Hyphae from the grapes were collected and DNA extracted using a Chelex based method. The ITS region of fungi was amplified and sequenced using primers ITS1/ITS4 (ITS1, 5′-TCC GTA GGT GAA CCT GCG G-3′; ITS4, 5′-TCC TCC GCT TAT TGA TAT GC-3′). The sequences were compared against GenBank database with BLAST (Basic Logical Alignment Search Tool). Fungal growths were identified as Talaromyces and Penicillium based on both ITS (Internal Transcribed Spacer) sequencing and morphology. Other microbes, including yeasts, also grow on grapes, but were not identified.

2.3. Comparisons of Lines on Laboratory Medium

We set up male and female pairs from isofemale lines 3 generations after they were established to score life history performance of wAu-infected flies on an agar-treacle-dead yeast laboratory medium (c.f. [24]). Pairs of males and females (2 per line, 92 in total) when flies were five days post eclosion were set up on spoons containing the laboratory medium brushed with live yeast paste to encourage egg-laying. Spoons were held at 25 °C and replaced every day over three days (i.e., tested on days 5, 6, and 7 post eclosion). Eggs on spoons were counted to estimate fecundity, and eggs were then transferred to vials with 20 mL of medium to estimate emergence rate, sex ratio, and development time (estimated by scoring emergence daily). Females and males from the same vial were separated into two glass vials and dried at 50 °C for 24 h in an oven. The weight of females and males from each vial was obtained as described above, with vials with less than 4 individuals not weighed.
We compared fecundity, development time, and offspring weight between wAu-infected and uninfected lines for each parental age class (5–7 days post eclosion). Data are again presented as violin plots (R software package vioplot). Fecundity, emergence rate, and offspring weight (which were normally distributed as tested by Kolmogorov-Smirnov tests) were compared among wAu-infected and uninfected lines by ANOVA, with infection and adult parental age as fixed effects. Rank-based tests with the raov function from the package, Rfit, was used to test the overall effect of infection, adult age, and their interaction on the emergence rate, given that the emergence rate was not normally distributed by a Komogorov-Smirnov test. Development time was compared among wAu-infected and uninfected lines with an ANCOVA, which included infection status and adult parental age as fixed effects and emergence number as a covariate because vials, with more larvae having the potential to have a delayed development time due to crowding. Finally, we compared sex ratio between offspring from wAu-infected and uninfected lines (pooled across parental age) with Wilcoxon-tests (given that the sex ratio was not normally distributed).

2.4. Toxin Effects

Penicillium and related fungi growing on grapes produce Ochratoxin A (OTA), which is one of the most-abundant food-contaminating mycotoxins [25]. Eggs from isofemale lines that had been established for 6 generations were tested for toxin responses by rearing larvae in two conditions where the toxin Ochratoxin A (Sigma-Aldrich, St. Louis, MO, USA) had been added to 12 mL of laboratory medium in vials (0.05 ng/mL OTA and 0.5 ng/mL OTA) in medium as well as untreated controls. The lower concentration is relevant to field conditions [26] while the higher concentration was included as an extreme condition. For each line, 30 eggs obtained from spoons were added to medium, and emergence rate, development time, as well as weight of adults were scored.
As in the previous experiment, given that emergence data were not normally distributed, rank-based tests with the raov function from the package, Rfit, were used to test the overall effect of infection, medium type, and their interaction on the emergence rate. We also compared the emergence rate between wAu-infected and uninfected lines of the adults with Wilcoxon-tests and provided kernel density (violin) plots (R software package, vioplot). We compared development time (which was normally distributed) between wAu-infected and uninfected lines with t-tests and used a two-way ANOVA to test for the effects of infection and medium type as fixed effects. Note that because the initial larval density in vials was constant (assuming that the hatchability of the eggs did not vary), we did not include density as a covariate as in the previous experiment. Dry weight was compared between wAu-infected and uninfected lines and medium type (fixed effects) with an ANOVA, while t-tests were used to compare infection status within medium type.

3. Results

3.1. Comparison of Lines on Grapes with Fungal Growth

There was an interaction between Wolbachia infection and fungus category for the number of offspring produced by 9-day-old parents in the ANOVA (Table 2). Offspring numbers were not significantly different between wAu-infected and uninfected flies feeding on low and medium levels of fungus infestations (Figure 1). However, the mean offspring number of wAu-infected flies (mean = 39.3, SD = 28.4) was substantially (132%) higher than the mean of the uninfected lines (mean = 16.6, SD = 16.0) when feeding on severely fungus-infested grapes, accounting for the interaction effect (Figure 1). It seems that the wAu lines produced larger numbers of offspring than uninfected lines due to higher fecundity and/or larval survival, but only when there was abundant mycelial growth.
For the 13-day-old parents, an interaction of Wolbachia infection and fungus was only marginally significant, and the number of offspring produced differed significantly by level of fungus infestation (Table 2). As for the 9-day-old parental comparison, the infected and uninfected parents differed in offspring production; wAu-infected flies averaged 43.8 (SD = 34.5) offspring compared to 24.4 (SD = 24.4) offspring for the uninfected flies, reflecting an advantage in offspring production over uninfected lines of 79.5% when grapes were heavily infested (Figure 1). A heavy infestation decreased offspring number.
For the average weight of females and males, the ANOVAs showed a significant difference among fungal levels due to a decrease of 35.8% (females) and 22.7% (males) in weight at higher fungal infestation levels compared to low levels. However, no significant difference was found for female and male weight between wAu-infected and uninfected lines (Table 2). The average sex ratio of wAu-infected offspring did not differ from that of uninfected offspring (Wilcoxon-test, p = 0.637).

3.2. Comparison of Lines on Laboratory Medium

Data were obtained from 40 fly pairs from the wAu-infected lines and 42 fly pairs from uninfected lines once some replicates without a complete set of parents were removed. Fecundity measured through the number of eggs laid on spoons was significantly different between wAu-infected and uninfected lines in 5-day-old parents and 7-day-old parents, with a mean fecundity in infected lines of 37.2 (SD = 7.2) eggs in 5-day-old parents and 27.3 (SD = 12.5) eggs in 7-day-old parents, compared to values of 31.5 (SD = 14.0) and 20.7 (SD = 8.6) for uninfected lines, respectively. For 6-day-old parents, mean fecundities of infected and uninfected pairs differed in the same direction, but this was not significant (Figure 2a). An ANOVA showed no interaction of infection status and parental age. There was a significant effect of infection (Table 3, increase of 20.5% overall in wAu-infected flies), and also for parental age (Table 3).
Sex-ratio was not significantly different in offspring produced by the wAu-infected and uninfected flies (Wilcoxon-test, p = 0.637). For the emergence rate, the rank-based test (Table 3) also indicated that wAu-infected flies produced offspring with similar emergence rates than uninfected flies.
For development time and offspring weight, the ANCOVAs did not show any interaction between the parental age and the infection status, but showed a significant effect of parental age, and a difference between wAu-infected and uninfected flies (Table 3). Development time was slowed in offspring from 7-day-old parents by an average of 2% when compared to the offspring of 5-day-old parents, while development time in the wAu-infected offspring was also slowed by 1.3% (Figure 2b) compared to the uninfecteds with a mean development time of 10.71 (SD = 0.41) days for the infected offspring compared to 10.60 (SD = 0.40) days for the uninfected offspring. Note that the covariate (emergence number) was significantly (p < 0.05) and positively (r = 0.124) associated with development time, reflecting a longer development when larval density was relatively higher. The average weight of female offspring from wAu infected offspring was 0.257 (SD = 0.025) mg, which was significantly heavier (Table 3) than the mean weight of uninfected line offspring (mean = 0.251, SD = 0.024), although the increase in weight in wAu-infected offspring was relatively small (2.7%). For males, the mean weight of infected offspring was 0.190 (SD = 0.015) mg compared to a mean of 0.186 (SD = 0.0186) mg for uninfected offspring, a significant difference (Table 3) with an increase of 2.5% to also produce extended development times (Figure 2b–d).

3.3. Toxin Effects

Emergence rate was lower when the toxin was present and decreased with concentration (Figure 3a), resulting in a significant effect of the OTA condition on this variable (Table 4). There was no interaction between infection and the OTA condition (Table 4), but a marginally significant overall effect of the infection due to the somewhat higher emergence of the infected lines (89.4%) compared to the control lines (86.9%).
For development time, the ANOVA (Table 4) showed a significant difference between wAu-infected and uninfected individuals (Table 4). The wAu-infected larvae had a longer development time, averaging 11.04 (SD = 0.098) days compared to the uninfecteds, which averaged 10.94 (SD = 0.13) days, representing an increase of 0.86% for the infected larvae. This increase was particularly evident for normal medium and for medium with 0.05ng/mL OTA added (Figure 3b).
The weight of female adults infected with wAu was higher than for uninfected females (Figure 3c,d) regardless of whether there was OTA in the medium; the effects of the infection and toxin were significant in the ANOVA and there was also an interaction effect (Table 4). On average, the female weight of infected flies (mean = 0.362, SD = 0.029 mg) was higher than that of the uninfected flies (mean = 0.336, SD = 0.040 mg) with an overall difference of 0.076%. The toxin presence decreased female weight, particularly at the highest concentration tested (Figure 3). There was no impact of the toxin on male weight, although there was a marginally significant effect of the infection on male weight (Table 4).

4. Discussion

The current work points to fitness advantages as well as costs associated with the wAu infection of D. simulans. The fitness advantages relate to a larger size (particularly in females) and a higher fecundity, likely to be a reflection of this larger size, while the cost is associated with an extended development time, again likely reflecting the fact that wAu-infected individuals would have pupated at a larger larval size (although this was not measured directly). These fitness effects probably contributed to the advantage that wAu-infected hosts showed under conditions where flies were reared on fruit contaminated with mycelial growth where a marked difference in productivity was detected. In particular, we anticipated that the higher fecundity of the infected lines observed on laboratory medium increased productivity although it is not obvious why this was only evident when there was mycelial growth. Perhaps there is an additional advantage stemming from a slower development time when mycelial growth is present (or there are other concomitant changes in microbes including yeasts), which allows the larvae to develop under the stressful fungal conditions (stress being evident from the small size of the emerging flies and lower emergence rate). Perhaps infected individuals somehow promote the growth of fungi. Although we were unable to link effects to one of the chemicals produced by the fungi (OTA), other compounds associated with microbial growth might still be important. To further investigate the potential effects of microbes on Wolbachia fitness, additional experiments are required, such as culturing specific strains of fungi, yeasts, and so on under sterile conditions, where nutrition can be carefully controlled [27].
Previous results from Kriesner et al. [20] highlighted the strong fitness effect of wAu in population cages, where flies were held with laboratory medium at a very large population size. These authors included two experimental treatments, where medium was “conditioned” to increase potential microbial activity by being exposed to fly populations to allow for oral or faecal transmission, as well as a situation where fruit was not conditioned before being used as a breeding site. However, these treatments did not influence the relative advantage of the wAu infection in D. simulans as reflected by the rate of increase of the infection in population cages. Instead, we suspect that based on the current results, there was an advantage of wAu associated with some of the fitness effects detected in the current study, although the advantage of wAu-infected hosts under conditions of abundant fungal mycelial growth would not have been a component of that study.
Our findings contrast with some previous work to identify fitness differences associated with the wAu-infected strain. In the initial paper on this strain, Hoffmann et al. [18] suggested that wAu acted as a neutral variant in natural D. simulans populations, being transmitted maternally at a high incidence near 100%, but not having any notable deleterious fitness effects (and not causing CI). The fitness effect measured in that study was fecundity, and it was measured by comparing an infected population with a derived population cured by exposure to tetracycline for a generation. However, the “cured” line was only given one generation without tetracycline to recover whereas it is now known that recovery from antibiotic treatment can take several generations [28], although the life history effect of tetracycline treatment unrelated to Wolbachia is unknown. The apparent absence of wAu fitness effects is also contradicted by the rapid increases in wAu frequency seen in natural populations of D. simulans in Australia prior to the arrival of the wRi infection, which has subsequently spread through all eastern Australian populations due to CI [24]. The relatively high frequency of wAu we observed here in a Western Australian population also confirms previous observations [24] that wAu has increased from a low frequency in Australia where it was first detected [18] to being relatively common in populations.
In the current study, a different experimental design was used to measure fitness effects associated with Wolbachia than is normally applied in fitness comparisons. We did not use tetracycline curing, which has been commonly applied, with the advantage that this approach controls for differences in nuclear background at least in the absence of genetic drift in lines after curing. In our design, the comparison of infected and uninfected lines involved multiple independently-derived isofemale lines obtained from the same field sample. This means that the nuclear genetic background of the infected and uninfected lines should both be representative of what is found in the field (given that each isofemale line captures a substantial fraction of genetic variation from the field population [29]), and randomized with respect to the infection given that nuclear-Wolbachia disequilibrium is not expected in populations [30]. While we averaged the effects of the infection across a sample of 20 or more nuclear backgrounds maintained identically as uninfected isofemale lines, fitness effects associated with Wolbachia infections are often undertaken through a comparison of one “cured” and uncured line or a comparison of one naturally infected or uninfected population. Comparisons based on a limited set of lines can include repeated backcrossing to homogenize lines to a common stock or reciprocal crosses to increase the likelihood that any detected effects are linked to Wolbachia, but the current design has the added advantage that inferences of infection effects are obtained on randomized nuclear backgrounds that are representative of those in the sampled field population. Obviously, such a design is much simpler to achieve in Drosophila, where a large number of lines can be readily reared as compared to in (say) mosquitoes, where line culture across multiple populations is much more challenging.
Within the limitations of generating meaningful fitness estimates, other studies have pointed to life history advantages of Wolbachia. In particular, the wRi infection of D. simulans was originally considered deleterious, but has evolved to increase fecundity [31]. The Wolbachia infection in D. suzukii is also expected to be beneficial given that it has increased in frequency, but does not cause CI [17] and there is evidence for a fitness advantage for Wolbachia associated with D. mauritiana [32]. Fitness benefits of Wolbachia and other endosymbionts have been well established in other organisms, particularly through the provision of nutrition. These include beneficial effects of Cardinium and Wolbachia on fecundity and survival of spider mites [33] and of Wolbachia on the growth and reproduction of bed bugs [14]. In some cases, the basis for fitness benefits have been established, such as the provisioning of riboflavin by Wolbachia in bedbugs [14] and the provisioning of amino acids by Buchnera in aphids [2,34]. However, there is also a level of inconsistency in reported fitness effects associated with Wolbachia infections, particularly in mosquitoes [35].
Apart from showing that life history traits are altered by the wAu infection, our experiments also tested for environment-specific fitness effects following the results linked to mycelial growth. In this respect, the toxin experiment failed to show specific impacts under certain conditions. In discussing fitness effects associated with Wolbachia, there is often a tendency in the literature to immediately attribute fitness associations to environment-specific components. In particular, the reduction in viral transmission associated with many Wolbachia strains in Drosophila [36,37] is seen as a way in which Wolbachia infections can be favoured in conditions where virus is present. However, at least in Drosophila, it is not yet apparent whether this applies; although several Wolbachia strains block viruses under laboratory conditions, natural viruses from field-collected flies and lines have not yet been linked to the presence of Wolbachia, at least with respect to the Wolbachia infection from D. melanogaster [38] that is persistently polymorphic in natural populations [19,39].
Our results also point to effects of parental age on the life history traits. This includes a reduction in offspring number when flies were reared on grapes, and a decrease in fecundity, increase in development time, and inconsistent effect of offspring weight in older parents. Female age effects on reproduction have been reported previously (e.g., [40]), and in our case, may reflect an exhaustion of eggs in older females, while the longer development time and weight changes may reflect changes in provisioning of eggs in older females. However, apparent age effects should be treated cautiously because the flies were transferred to new media/fruit, which may have differed from the initial batch used.

5. Conclusions

In conclusion, we found life history effects associated with the wAu infection that provide an advantage and potential cost to flies under at least one set of laboratory conditions, expressed under a field-relevant nuclear background. The fitness effects may be partly mediated through changes in fly size, a trait which has previously been related to the wMel infection in some (but not all) field collections [41]. Knowledge of such fitness effects is important in predicting the ability of infections to invade natural populations [22], particularly in the absence of cytoplasmic incompatibility. The advantage of wAu over uninfected flies may account for the increase in the frequency of this infection both in population cages and under field conditions seen in Australian populations [20,24]. The generality of these findings remains to be tested with other infections and hosts, but are likely to differ depending on Wolbachia density, tropism, and other factors [35]. The genomic basis of these effects also remains to be determined; the wAu infection lacks a wMel ortholog involved in cytoplasmic incompatibility, but also has numerous other differences [42] that may be linked to life history effects.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/10/5/126/s1, Figure S1: Fungus-infestation levels on grapes with three levels recognized: (a) Level 1, less than 1 piece of 6 pieces of grapes was covered by visible mycelial growth; (b) level 2, 2 to 4 pieces of grapes were covered by visible growth; (c) level 3, 5 to 6 pieces of grapes were covered by visible growth.

Author Contributions

Conceptualization, A.A.H., L.-J.C., and W.J.; methodology, L.-J.C., W.J. and A.A.H.; formal analysis, L.-J.C., W.J., and A.A.H.; investigation, L.-J.C. and W.J.; resources, A.A.H.; writing and editing, A.A.H., L.-J.C., and W.J.; visualization, L.-J.C. and W.J.; funding acquisition, A.A.H.

Funding

This research was funded by National Institutes of Health: 5R01GM104325; National Health and Medical Research Council: 1118640; Australian Research Council: LE150100083; Australian Research Council: DP120100916; Beijing Academy of Agriculture and Forestry Sciences: GJHZ2017-5.

Acknowledgments

This work was supported by a grant from the National Institutes of Health (5R01GM104325), a fellowship from the National Health and Medical Research Council (1118640), an infrastructure grant and discovery grant from the Australian Research Council (LE150100083, DP120100916) and grants from Beijing Academy of Agriculture and Forestry Sciences (GJHZ2017-5) and from Shanghai Normal University (17DZ2252700). Thanks to Kelly Richardson and Perran Ross for their assistance with laboratory work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Pais, R.; Lohs, C.; Wu, Y.N.; Wang, J.W.; Aksoy, S. The obligate mutualist Wigglesworthia glossinidia influences reproduction, digestion, and immunity processes of its host, the tsetse fly. Appl. Environ. Microbiol. 2008, 74, 5965–5974. [Google Scholar] [CrossRef] [PubMed]
  2. Douglas, A.E. Nutritional interactions in insect-microbial symbioses: Aphids and their symbiotic bacteria Buchnera. Annu. Rev. Entomol. 1998, 43, 17–37. [Google Scholar] [CrossRef]
  3. Duron, O.; Morel, O.; Noel, V.; Buysse, M.; Binetruy, F.; Lancelot, R.; Loire, E.; Menard, C.; Bouchez, O.; Vavre, F.; et al. Tick-bacteria mutualism depends on B vitamin synthesis pthways. Curr. Biol. 2018, 28, 1896–1902. [Google Scholar] [CrossRef] [PubMed]
  4. Koga, R.; Tsuchida, T.; Fukatsu, T. Changing partners in an obligate symbiosis: A facultative endosymbiont can compensate for loss of the essential endosymbiont Buchnera in an aphid. Proc. R. Soc. B Biol. Sci. 2003, 270, 2543–2550. [Google Scholar] [CrossRef] [PubMed]
  5. Sazama, E.J.; Bosch, M.J.; Shouldis, C.S.; Ouellette, S.P.; Wesner, J.S. Incidence of Wolbachia in aquatic insects. Ecol. Evol. 2017, 7, 1165–1169. [Google Scholar] [CrossRef] [PubMed]
  6. 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. [Google Scholar] [CrossRef]
  7. O’Neill, S.L.; Hoffmann, A.A.; Werren, J.H. Influential Passengers: Inherited Microorganisms and Arthropod Reproduction; Oxford University Press: Oxford, UK, 1997. [Google Scholar]
  8. Ant, T.H.; Herd, C.S.; Geoghegan, V.; Hoffmann, A.A.; Sinkins, S.P. The Wolbachia strain wAu provides highly efficient virus transmission blocking in Aedes aegypti. PLoS Pathog. 2018, 14, e1006815. [Google Scholar] [CrossRef] [PubMed]
  9. Walker, T.; Johnson, P.H.; Moreira, L.A.; Iturbe-Ormaetxe, I.; Frentiu, F.D.; McMeniman, C.J.; Leong, Y.S.; Dong, Y.; Axford, J.; Kriesner, P.; et al. The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations. Nature 2011, 476, 450–453. [Google Scholar] [CrossRef] [PubMed]
  10. Moreira, L.A.; Iturbe-Ormaetxe, I.; Jeffery, J.A.; Lu, G.J.; Pyke, A.T.; Hedges, L.M.; Rocha, B.C.; Hall-Mendelin, S.; Day, A.; Riegler, M.; et al. A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya, and Plasmodium. Cell 2009, 139, 1268–1278. [Google Scholar] [CrossRef] [PubMed]
  11. Bian, G.W.; Xu, Y.; Lu, P.; Xie, Y.; Xi, Z.Y. The endosymbiotic bacterium Wolbachia induces resistance to dengue virus in Aedes aegypti. Plos Pathogens 2010, 6, e1000833. [Google Scholar] [CrossRef] [PubMed]
  12. Caspari, E.; Watson, G.S. On the evolutionary importance of cytoplasmic sterility in mosquitoes. Evolution 1959, 13, 568–570. [Google Scholar] [CrossRef]
  13. Turelli, M.; Hoffmann, A.A. Rapid spread of an inherited incompatibility factor in California Drosophila. Nature 1991, 353, 440–442. [Google Scholar] [CrossRef] [PubMed]
  14. Moriyama, M.; Nikoh, N.; Hosokawa, T.; Fukatsu, T. Riboflavin provisioning underlies Wolbachia’s fitness contribution to its insect host. Mbio 2015, 6, e01732-15. [Google Scholar] [CrossRef]
  15. Taylor, M.J.; Bandi, C.; Hoerauf, A. Wolbachia bacterial endosymbionts of filarial nematodes. In Advances in Parasitology, Volume 60; Baker, J.R., Muller, R., Rollinson, D., Eds.; Elsevier Academic Press Inc.: San Diego, CA, USA, 2005; Volume 60, pp. 245–284. [Google Scholar]
  16. Brownlie, J.C.; Cass, B.N.; Riegler, M.; Witsenburg, J.J.; Iturbe-Ormaetxe, I.; McGraw, E.A.; O’Neill, S.L. Evidence for metabolic provisioning by a common invertebrate endosymbiont, Wolbachia pipientis, during periods of nutritional stress. Plos Pathog. 2009, 5, e1000368. [Google Scholar] [CrossRef]
  17. Hamm, C.A.; Begun, D.J.; Vo, A.; Smith, C.C.R.; Saelao, P.; Shaver, A.O.; Jaenike, J.; Turelli, M. Wolbachia do not live by reproductive manipulation alone: Infection polymorphism in Drosophila suzukii and D. subpulchrella. Mol. Ecol. 2014, 23, 4871–4885. [Google Scholar] [CrossRef] [PubMed]
  18. Hoffmann, A.A.; Clancy, D.; Duncan, J. Naturally-occurring Wolbachia infection in Drosophila simulans that does not cause cytoplasmic incompatibility. Heredity 1996, 76, 1–8. [Google Scholar] [CrossRef] [PubMed]
  19. Kriesner, P.; Conner, W.R.; Weeks, A.R.; Turelli, M.; Hoffmann, A.A. Persistence of a Wolbachia infection frequency cline in Drosophila melanogaster and the possible role of reproductive dormancy. Evolution 2016, 70, 979–997. [Google Scholar] [CrossRef] [PubMed]
  20. Kriesner, P.; Hoffmann, A.A. Rapid spread of a Wolbachia infection that does not affect host reproduction in Drosophila simulans cage populations. Evolution 2018, 72, 1475–1487. [Google Scholar] [CrossRef] [PubMed]
  21. Dobson, S.L.; Marsland, E.J.; Rattanadechakul, W. Mutualistic Wolbachia infection in Aedes albopictus: Accelerating cytoplasmic drive. Genetics 2002, 160, 1087–1094. [Google Scholar]
  22. Fenton, A.; Johnson, K.N.; Brownlie, J.C.; Hurst, G.D.D. Solving the Wolbachia paradox: Modeling the tripartite interaction between host, Wolbachia, and a natural enemy. Am. Nat. 2011, 178, 333–342. [Google Scholar] [CrossRef]
  23. Lee, S.F.; White, V.L.; Weeks, A.R.; Hoffmann, A.A.; Endersby, N.M. High-throughput PCR assays to monitor Wolbachia infection in the dengue mosquito (Aedes aegypti) and Drosophila simulans. Appl. Environ. Microbiol. 2012, 78, 4740–4743. [Google Scholar] [CrossRef]
  24. Kriesner, P.; Hoffmann, A.A.; Lee, S.F.; Turelli, M.; Weeks, A.R. Rapid sequential spread of two Wolbachia variants in Drosophila simulans. PLoS Pathog. 2013, 9, e1003607. [Google Scholar] [CrossRef] [PubMed]
  25. Alkhayyat, F.; Yu, J.-H. Upstream regulation of mycotoxin biosynthesis. In Advances in Applied Microbiology; Sariaslani, S., Gadd, G.M., Eds.; Academic Press: New York, NY, USA, 2014; Volume 86, pp. 251–278. [Google Scholar]
  26. Mutlu, A.G. Increase in mitochondrial DNA copy number in response to Ochratoxin A and methanol-induced mitochondrial DNA damage in Drosophila. Bull. Environ. Contam. Toxicol. 2012, 89, 1129–1132. [Google Scholar] [CrossRef]
  27. Begg, M.; Robertson, F.W. The nutritional requirements of Drosophila melanogaster. J. Exp. Biol. 1950, 26, 380–387. [Google Scholar] [CrossRef]
  28. Ballard, J.W.O.; Melvin, R.G. Tetracycline treatment influences mitochondrial metabolism and mtDNA density two generations after treatment in Drosophila. Insect Mol. Biol. 2007, 16, 799–802. [Google Scholar] [CrossRef] [PubMed]
  29. Hoffmann, A.A.; Parsons, P.A. The analysis of quantitative variation in natural populations with isofemale strains. Genet. Sel. Evol. 1988, 20, 87–98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Hoffmann, A.A.; Turelli, M. Cytoplasmic incompatibility in insects. In Influential Passengers: Microorganisms and Invertebrate Reproduction; O’Neill, S., Hoffmann, A., Werren, J., Eds.; Oxford University Press: Oxford, UK, 1997; pp. 42–80. [Google Scholar]
  31. Weeks, A.R.; Turelli, M.; Harcombe, W.R.; Reynolds, K.T.; Hoffmann, A.A. From parasite to mutualist: Rapid evolution of Wolbachia in natural populations of Drosophila. PloS Biol. 2007, 5, e114. [Google Scholar] [CrossRef] [PubMed]
  32. Fast, E.M.; Toomey, M.E.; Panaram, K.; Desjardins, D.; Kolaczyk, E.D.; Frydman, H.M. Wolbachia enhance Drosophila stem cell proliferation and target the germline stem cell niche. Science 2011, 334, 990–992. [Google Scholar] [CrossRef] [PubMed]
  33. Xie, R.R.; Sun, J.T.; Xue, X.F.; Hong, X.Y. Cytoplasmic incompatibility and fitness benefits in the two-spotted spider mite Tetranychus urticae (red form) doubly infected with Wolbachia and Cardinium. Syst. Appl. Acarol. 2016, 21, 1161–1173. [Google Scholar] [CrossRef]
  34. Vala, F.; Egas, M.; Breeuwer, J.A.J.; Sabelis, M.W. Wolbachia affects oviposition and mating behaviour of its spider mite host. J. Evol. Biol. 2004, 17, 692–700. [Google Scholar] [CrossRef]
  35. Hoffmann, A.A.; Ross, P.A.; Rasic, G. Wolbachia strains for disease control: Ecological and evolutionary considerations. Evol. Appl. 2015, 8, 751–768. [Google Scholar] [CrossRef] [PubMed]
  36. Martinez, J.; Longdon, B.; Bauer, S.; Chan, Y.S.; Miller, W.J.; Bourtzis, K.; Teixeira, L.; Jiggins, F.M. Symbionts commonly provide broad spectrum resistance to viruses in insects: A comparative analysis of Wolbachia strains. Plos Pathog. 2014, 10. [Google Scholar] [CrossRef] [PubMed]
  37. Osborne, S.E.; San Leong, Y.; O’Neill, S.L.; Johnson, K.N. Variation in antiviral protection mediated by different Wolbachia strains in Drosophila simulans. PLoS Pathog. 2009, 5, e1000656. [Google Scholar] [CrossRef] [PubMed]
  38. Shi, M.; White, V.L.; Schlub, T.; Eden, J.-S.; Hoffmann, A.A.; Holmes, E.C. No detectable effect of Wolbachia w Mel on the prevalence and abundance of the RNA virome of Drosophila melanogaster. Proc. R. Soc. B Biol. Sci. 2018, 285, 20181165. [Google Scholar] [CrossRef] [PubMed]
  39. Hoffmann, A.A. Partial cytoplasmic incompatibility between two Australian populations of Drosophila melanogaster. Entomol. Exp. Appl. 1988, 48, 61–67. [Google Scholar] [CrossRef]
  40. Koch, R.E.; Phillips, J.M.; Camus, M.F.; Dowling, D.K. Maternal age effects on fecundity and offspring egg-to-adult viability are not affected by mitochondrial haplotype. Ecol. Evol. 2018, 8, 10722–10732. [Google Scholar] [CrossRef] [PubMed]
  41. Hoffmann, A.A.; Hercus, M.; Dagher, H. Population dynamics of the Wolbachia infection causing cytoplasmic incompatibility in Drosophila melanogaster. Genetics 1998, 148, 221–231. [Google Scholar] [PubMed]
  42. Sutton, E.R.; Harris, S.R.; Parkhill, J.; Sinkins, S.P. Comparative genome analysis of Wolbachia strain wAu. BMC Genom. 2014, 15, 928–942. [Google Scholar] [CrossRef]
Figure 1. Violin plots for distribution of offspring number of wAu-infected and uninfected (wt) lines at different levels of fungal growth. (a) Offspring number of 9-day-old parents; (b) offspring number of 13-day-old parents. Levels 1 to 3 indicate the levels of fungus infestation on grapes as described in the text. Medians are indicated by central points, and the interquartile range by the solid bar, while the dotted line indicates the 95% confidence interval.
Figure 1. Violin plots for distribution of offspring number of wAu-infected and uninfected (wt) lines at different levels of fungal growth. (a) Offspring number of 9-day-old parents; (b) offspring number of 13-day-old parents. Levels 1 to 3 indicate the levels of fungus infestation on grapes as described in the text. Medians are indicated by central points, and the interquartile range by the solid bar, while the dotted line indicates the 95% confidence interval.
Insects 10 00126 g001
Figure 2. Violin plots for (a) fecundity (number of eggs), (b) development time, (c) female weight, and (d) male weight of wAu-infected and uninfected lines reared on laboratory medium with parents transferred across three days to new medium (5-, 6-, and 7-day-old parents). Medians are indicated by central points, and the interquartile range by the solid bar, while the dotted line indicates the 95% confidence interval.
Figure 2. Violin plots for (a) fecundity (number of eggs), (b) development time, (c) female weight, and (d) male weight of wAu-infected and uninfected lines reared on laboratory medium with parents transferred across three days to new medium (5-, 6-, and 7-day-old parents). Medians are indicated by central points, and the interquartile range by the solid bar, while the dotted line indicates the 95% confidence interval.
Insects 10 00126 g002
Figure 3. Violin plots for (a) emergence rate, (b) development time, (c) female weight, and (d) male weight of wAu-infected and uninfected (wt) larvae reared on three medium conditions (no OTA (Ochratoxin), 0.05 ng/mL OTA, and 0.5 ng/mL OTA in medium). Medians are indicated by central points, and the interquartile range by the solid bar, while the dotted line indicates the 95% confidence interval.
Figure 3. Violin plots for (a) emergence rate, (b) development time, (c) female weight, and (d) male weight of wAu-infected and uninfected (wt) larvae reared on three medium conditions (no OTA (Ochratoxin), 0.05 ng/mL OTA, and 0.5 ng/mL OTA in medium). Medians are indicated by central points, and the interquartile range by the solid bar, while the dotted line indicates the 95% confidence interval.
Insects 10 00126 g003
Table 1. Outline of experiments undertaken.
Table 1. Outline of experiments undertaken.
ExperimentNumber of Isofemale Lines UsedGeneration Since Isofemale Lines EstablishedTraits Scored
Comparison of lines on grapes with fungal growth23 wAu-infected and 23 uninfected lines2Number of offspring produced, dry weight, sex ratio
Comparisons of lines on laboratory mediumAs above3Fecundity, egg-to-adult emergence, sex ratio, development time, dry weight
Toxin effectsAs above6Egg-to-adult emergence, development time, dry weight
Table 2. ANOVAs (MS = mean square, F ratio, probability) for offspring number and offspring weight when flies were allowed to produce offspring on grapes with fungal growth; parental flies produced offspring in one set of bottles when they were 9 to 11 days post-eclosion, and a second set of bottles when they were 12 to 14 days post-eclosion. Significant results (p < 0.05) are indicated in bold.
Table 2. ANOVAs (MS = mean square, F ratio, probability) for offspring number and offspring weight when flies were allowed to produce offspring on grapes with fungal growth; parental flies produced offspring in one set of bottles when they were 9 to 11 days post-eclosion, and a second set of bottles when they were 12 to 14 days post-eclosion. Significant results (p < 0.05) are indicated in bold.
VariableInfection (df = 1)Fungal Level
(df = 2)
Infection by Fungal Level (df = 2)Error MS (df)
MSF (P)MSF (P)MSF (P)
Offspring number (parents 9–11 days old)24314.708 (0.032)23554.561 (0.012)18753.631 (0.030)516 (118)
Offspring number (parents 12–14 days old)25193.579 (0.061)1092115.52 (<0.001)18842.676 (0.074)704 (104)
Weight (mg, female)0.00000.001 (0.974)0.157737.99 (<0.001)0.00120.292 (0.747)0.0042 (191)
Weight (mg, male)0.00080.687 (0.408)0.025921.12 (<0.001)0.00050.386 (0.680)0.0012 (193)
Table 3. ANOVAs (MS = mean square or MRD = mean square for emergence based on ranked data, F ratio, probability) for fecundity and emergence rate, as well as an ANCOVA performed on preimaginal development when flies were reared on laboratory medium and transferred daily for three days to produce different parental ages. Significant results (p < 0.05) are indicated in bold. The ANCOVA on development time included the emergence number as a covariate.
Table 3. ANOVAs (MS = mean square or MRD = mean square for emergence based on ranked data, F ratio, probability) for fecundity and emergence rate, as well as an ANCOVA performed on preimaginal development when flies were reared on laboratory medium and transferred daily for three days to produce different parental ages. Significant results (p < 0.05) are indicated in bold. The ANCOVA on development time included the emergence number as a covariate.
VariableInfection (df = 1)Age (df = 2)Infection by Age (df = 2)Error MS (df)
MS/MRDF (P)MS/MRDF (P)MS/MRDF (P)
Fecundity1778.915.96 (<0.001)2206.119.80 (<0.001)39.00.350 (0.705)111.4 (240)
Emergence rate (ranked data)0.000150.0012(0.972)0.05290.441 (0.644)0.04950.413 (0.662)- (234)
Development time in days 11.1435 7.393 (0.007)0.80825.225 (0.007)0.1354 0.875 (0.418)0.1585 (229)
Weight (mg, female)0.00264.532 (0.034)0.00315.524 (0.005)0.00111.929 (0.148)0.5530 (208)
Weight (mg, male)0.00114.103 (0.044)0.00155.693 (0.004)0.00031.094 (0.337) 0.0003 (196)
1 ANCOVA with the emergence number included in the analysis. The number of flies emerged had a significant effect on development time (F(1,229) = 6.373, p = 0.012).
Table 4. ANOVAs (MS = mean square or MRD = mean square for emergence based on ranked data, F ratio, probability) for emergence rate, adult weight, and preimaginal development when flies were reared on medium with different OTA (Ochratoxin) levels. Significant results (p < 0.05) are indicated in bold.
Table 4. ANOVAs (MS = mean square or MRD = mean square for emergence based on ranked data, F ratio, probability) for emergence rate, adult weight, and preimaginal development when flies were reared on medium with different OTA (Ochratoxin) levels. Significant results (p < 0.05) are indicated in bold.
VariableInfection (df = 1)OTA (df = 2)Infection by OTA (df = 2)Error MS (df)
MS/MRDF (P)MS/MRDF (P)MS/MRDF (P)
Emergence rate (ranked data)0.12894.733 (0.032)0.688525.277 (<0.001)0.00760.280 (0.756)- (120)
Development time in days0.280125.419 (<0.001)0.10279.318 (<0.001)0.02892.619 (0.077)0.0110 (120)
Weight (mg, female)0.02119.4 (<0.001)0.0097.999 (<0.001)0.0044.042 (0.020)0.0011 (120)
Weight (mg, male)0.0000.359 (0.550)0.0012.846 (0.062)0.0001.075 (0.344) 0.0003 (120)

Share and Cite

MDPI and ACS Style

Cao, L.-J.; Jiang, W.; Hoffmann, A.A. Life History Effects Linked to an Advantage for wAu Wolbachia in Drosophila. Insects 2019, 10, 126. https://doi.org/10.3390/insects10050126

AMA Style

Cao L-J, Jiang W, Hoffmann AA. Life History Effects Linked to an Advantage for wAu Wolbachia in Drosophila. Insects. 2019; 10(5):126. https://doi.org/10.3390/insects10050126

Chicago/Turabian Style

Cao, Li-Jun, Weibin Jiang, and Ary A. Hoffmann. 2019. "Life History Effects Linked to an Advantage for wAu Wolbachia in Drosophila" Insects 10, no. 5: 126. https://doi.org/10.3390/insects10050126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop