Next Article in Journal
Soil Fungal Diversity and Ecology Assessed Using DNA Metabarcoding along a Deglaciated Chronosequence at Clearwater Mesa, James Ross Island, Antarctic Peninsula
Previous Article in Journal
New Psychoactive Substances Intoxications and Fatalities during the COVID-19 Epidemic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Friend or Foe: Symbiotic Bacteria in Bactrocera dorsalis–Parasitoid Associations

by
Rehemah Gwokyalya
1,2,*,
Christopher W. Weldon
2,
Jeremy Keith Herren
1,
Joseph Gichuhi
1,
Edward Edmond Makhulu
1,
Shepard Ndlela
1 and
Samira Abuelgasim Mohamed
1,*
1
International Centre of Insect Physiology and Ecology, Nairobi P.O. Box 30772-00100, Kenya
2
Department of Zoology and Entomology, University of Pretoria, Private Bag X20, Pretoria 0028, South Africa
*
Authors to whom correspondence should be addressed.
Biology 2023, 12(2), 274; https://doi.org/10.3390/biology12020274
Submission received: 9 November 2022 / Revised: 14 January 2023 / Accepted: 16 January 2023 / Published: 9 February 2023

Abstract

:

Simple Summary

Host-associated gut bacteria influence the eco-physiological functions of their hosts, for example, how they interact with natural enemies in their ecosystems, and hence the outcome of biological pest management regimes. Here, we assessed the effect of three common gut bacteria of the oriental fruit fly, Bactrocera dorsalis, on the biological control of this pest using parasitoid wasps. We found that some gut bacteria increased parasitoid emergence and positively impacted the size and fecundity of the emerging parasitoid offspring. We, therefore, conclude that some bacteria can be used as probiotics in the mass rearing of parasitoids to boost the biological control of the oriental fruit fly.

Abstract

Parasitoids are promising biocontrol agents of the devastating fruit fly, Bactrocera dorsalis. However, parasitoid performance is a function of several factors, including host-associated symbiotic bacteria. Providencia alcalifaciens, Citrobacter freundii, and Lactococcus lactis are among the symbiotic bacteria commonly associated with B. dorsalis, and they influence the eco-physiological functioning of this pest. However, whether these bacteria influence the interaction between this pest and its parasitoids is unknown. This study sought to elucidate the nature of the interaction of the parasitoids, Fopius arisanus, Diachasmimorpha longicaudata, and Psyttlia cosyrae with B. dorsalis as mediated by symbiotic bacteria. Three types of fly lines were used: axenic, symbiotic, and bacteria-mono-associated (Lactococcus lactis, Providencia alcalifaciens, and Citrobacter freundii). The suitable stages of each fly line were exposed to the respective parasitoid species and reared until the emergence of adult flies/parasitoids. Thereafter, data on the emergence and parasitoid fitness traits were recorded. No wasps emerged from the fly lines exposed to P. cosyrae. The highest emergence of F. arisanus and D. longicaudata was recorded in the L. lactis fly lines. The parasitoid progeny from the L. lactis and P. alcalifaciens fly lines had the longest developmental time and the largest body size. Conversely, parasitoid fecundity was significantly lower in the L. lactis lines, whereas the P. alcalifaciens lines significantly improved fecundity. These results elucidate some effects of bacterial symbionts on host–parasitoid interactions and their potential in enhancing parasitoid-oriented management strategies against B. dorsalis.

1. Introduction

The oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), is a major horticultural pest that infests numerous fruit and vegetable crops across Asia, Africa, Oceania, and the Americas, resulting in hefty losses to the horticultural sectors of affected countries [1,2,3,4]. The management of this invasive pest is largely dependent on chemical insecticides, with far-reaching consequences on human and environmental health [5]. Furthermore, the overreliance on these chemicals very often results in the development of pest resistance. To ameliorate these effects, integrated pest management (IPM) is slowly being embraced in various parts of the world [6,7,8]. One of the IPM tactics being deployed is the use of biological control agents such as parasitoids [9,10,11,12], which, when released in sufficient numbers under ideal conditions, can significantly suppress pest populations [11,12,13].
The egg-prepupal parasitoid Fopius arisanus (Sonan) (Hymenoptera: Braconidae) and the larval-prepupal parasitoid Diachasmimorpha longicaudata (Ashmead) (Hymenoptera: Braconidae) are co-evolved, natural enemies of B. dorsalis [14,15,16]. These parasitoids have been credited with the most outstanding successful biological control programs against tephritids in Hawaii [11,12] and French Polynesia [13,17]. Fopius arisanus and D. longicaudata were imported into Africa for the classical biological control of B. dorsalis [9,18,19], especially because the pest lacked efficient resident parasitoid species on the continent [20]. One of the native parasitoid species of Africa, Psyttalia cosyrae (Wilkinson) (Hymenoptera: Braconidae), has also been investigated for possible use in the management of this invasive pest alongside D. longicaudata and F. arisanus [8].
Previous research has shown great variations in parasitism success by these three parasitoids on B. dorsalis. Successful development has been recorded for D. longicaudata and F. arisanus [14,16,18,21] but none for P. cosyrae [15,22]. While these variations in parasitism success could be due to the ability of the parasitoids to surmount their hosts’ defense mechanisms [22], there has been an increasing appreciation of other factors such as host-associated symbionts and their roles in shaping the dynamics of host–parasitoid interactions [23,24,25,26,27,28].
Bacterial symbionts have been shown to alter the outcomes of host–parasitoid interactions, with the majority of the bacteria showing host-protective functions [23,26,29,30,31,32,33]. These symbiont-conferred protective phenotypes have been linked to the symbionts’ ability to synthesize endotoxins [34,35,36] or alter host immune and metabolic processes [37,38]. As a result, the development of the parasitoid offspring is hampered. In addition to their protective roles, symbiotic bacteria have been shown to affect host selection and oviposition cues by parasitoids, as well as the fitness parameters of parasitoid progenies, by influencing their developmental times, sex ratios, size, and survival [39,40,41,42,43,44]. The ability of bacterial symbionts to alter the dynamics of host–parasitoid interactions is significant for the ecological dynamics of the host insects and parasitoids. This relationship could be exploited to develop symbiont-based biocontrol management strategies for pests of agricultural importance, such as the invasive B. dorsalis.
Diverse communities of symbiotic bacteria are associated with B. dorsalis [4,45,46]. Some of these bacteria have been shown to alter the development, foraging behavior, mating, and immunity of their host [46,47,48,49,50]. Specifically, Lactococcus lactis, Providencia alcalifaciens, and Citrobacter freundii are among the most common gut bacteria associated with B. dorsalis [46,51,52]. These bacteria have been implicated in the detoxification of insecticides, performance of irradiated flies, host fly fitness, and susceptibility to entomopathogenic fungi [5,46,51]. Therefore, microbial symbionts are key mediators of the physiological and ecological roles of B. dorsalis. However, whether they affect the outcome of the interactions between B. dorsalis and its associated parasitoids has not been fully elucidated.
Since such bacterial symbionts present considerable physiological and ecological benefits to their host, B. dorsalis, it is of interest to determine whether they promote host-protective phenotypes against parasitic wasps and if they influence the fitness parameters of the parasitoid progeny. Hence, the current study aimed at determining whether L. lactis, P. alcalifaciens, and C. freundii, the symbionts commonly associated with B. dorsalis, (1) affect parasitoid oviposition behaviors, (2) confer general protection against the braconid parasitoids associated with B. dorsalis, and (3) alter the fitness parameters of D. longicaudata, F. arisanus, and P. cosyrae progeny.

2. Materials and Methods

2.1. Bactrocera dorsalis and Parasitoid Rearing

The Bactrocera dorsalis flies used in this study were reared on ripe, insecticide-free mangoes, as previously described [22], in the insectary of the International Centre of Insect Physiology and Ecology (icipe). The flies were derived from the eggs collected from a 15th-generation, laboratory-reared colony, which was initiated with field collections and was regularly infused with field collections and maintained, as previously described [53]. The eclosed adult flies were fed on yeast and water ad libitum [15]. The parasitoids used in this study, F. arisanus, D. longicaudata, and P. cosyrae, were also reared in the insectary at icipe at a temperature range of 25–27 °C, 60–70% relative humidity, and a 12:12 day:light photoperiod. Psyttalia cosyrae was reared on a laboratory colony of Ceratitis cosyra Walker (Diptera: Tephritidae) [54], while D. longicaudata and F. arisanus were reared on B. dorsalis as described by Mohamed et al. [9] and Mohamed et al. [18], respectively. The parasitoids D. longicaudata, F. arisanus, and P. cosyrae used in this study were at the 183rd, 164th, and 177th generations in the laboratory, respectively.

2.2. Bacteria Culturing

The bacteria P. alcalifaciens, L. lactis, and C. freundii used in this study were retrieved from glycerol stocks (stored at −80 °C) previously isolated from B. dorsalis [46]. Using a sterile tip, individual species, each 5 μL, were retrieved from the glycerol stocks, spread on a nutrient-rich medium, i.e., brain heart infusion (BHI) agar (Thermo Scientific™ Oxoid, Hampshire, UK), and incubated at 37 °C for 14 h. They were then re-cultured in BHI broth (Thermo Scientific™ Oxoid, Hampshire, UK) at 37 °C for 14 h in a shaking incubator.
The bacteria were harvested by centrifuging at 12879 RCF for five minutes, and the precipitate containing the bacteria was diluted in distilled water to 1 × 108 concentration.

2.3. Generation of Axenic, Symbiotic, and Mono-Associated Bacterial Lines

The generation of bacteria-free fly lines was performed as described by Gichuhi et al. [46]. Briefly, freshly laid B. dorsalis eggs were collected from perforated mango domes (used as oviposition sites), surface-sterilized in 70% ethanol, dechorionated in a 7% v/v sodium hypochlorite solution for three minutes and rinsed in distilled water thrice. To confirm the elimination of bacteria from the flies, the dechorionated embryos were inoculated in autoclaved liquid larval diets [55] (modified by excluding streptomycin and nipagen) and left to develop until the 2nd instar. Consequently, 10 2nd instar larvae were randomly selected and individually homogenized in 100 µL of phosphate-buffered saline (PBS), and the homogenate was plated on BHI agar plates. The plates were incubated at 37 °C for 14 h and observed for bacterial growth. There was no bacterial growth observed on the plates, and this confirmed the axenic state of the larvae.
From these dechorionated embryos, negative (axenic) and positive (symbiotic) controls, as well as the bacterial mono-associated lines (BMALs) (L. lactis or P. alcalifaciens or C. freundii), were generated. Axenic lines were generated by inoculating the dechorionated embryos on autoclaved artificial liquid larval diets [55] (modified by excluding streptomycin and nipagen) supplemented with distilled water only. The BMALs were generated by directly inoculating the dechorionated eggs on autoclaved artificial liquid larval diets [55] (modified by excluding streptomycin and nipagen) supplemented with 1000 μL of 1 × 108 CFU/mL inoculum of the respective bacterial isolates, individually. Symbiotic fly lines were generated by directly inoculating the dechorionated embryos on non-autoclaved larval diets previously fed upon by larvae with an intact microbiome, i.e., from a non-dechorionated embryo fly line reared under normal laboratory conditions. This unsterilized/recycled diet was selected because it would contain microbiota from the frass of the previously reared larvae. For each type of fly line, inoculations per replicate were carried out as previously described [46], with a slight modification in the size of tubes used: flat-based glass tube cylinders (75 × 25 mm). The tubes were capped with autoclaved cotton wool and maintained in a microbiological safety cabinet) at 25–27 °C, 60–70% relative humidity, and 12:12 day:light photoperiod. To minimize confounding effects, all fly lines were generated from the same batch of collected embryos that were dechorionated in the same batch before the respective inoculations and then reared under the exact same conditions.

2.4. Parasitic Wasp Infection of BMALs and Axenic and Symbiotic Fly Lines

For D. longicaudata and P. cosyrae parasitization assays, B. dorsalis larvae were used. Briefly, B. dorsalis embryos were dechorionated and raised as described above for each experimental group until the larvae reached the 2nd instar stage. Using soft forceps, a set of a hundred 2nd instar B. dorsalis larvae were randomly selected from each fly line and individually transferred to larval oviposition units containing autoclaved carrot diets [56] (modified by excluding nipagen). The carrot diets for the BMALs and axenic and symbiotic lines were supplemented with 500 μL (1 × 108) of the respective bacterial species, distilled water, or a liquid diet previously fed on by larvae with an intact microbiome, respectively. The oviposition units were transferred to a Perspex cage (12 × 12 × 12 cm) holding 10 D. longicaudata or P. cosyrae mated female parasitoids for parasitization. In a recent study, Gwokyalya et al. [22] demonstrated that D. longicaudata achieved oviposition above 90%, whereas P. cosyrae achieved similar oviposition rates when host exposure was carried out for six hours. Consequently, in this study, exposure to D. longicaudata was implemented for two hours, while exposure to P. cosyrae was carried out for six hours. Post-parasitization, the host larvae were retrieved from the oviposition units and transferred to the autoclaved carrot diets. The carrot diets were supplemented with 500 µL of the respective bacteria inoculum for the BMALs, a regular liquid larval diet for the symbiotic lines, or distilled water for the axenic lines. The host larvae were left to feed and subsequently pupariate. The puparia were transferred to four-liter transparent lunch boxes (18 × 11 × 15 cm) covered with a cotton mesh and monitored until fly and/or parasitoid emergence.
Fopius arisanus parasitization assays were conducted using B. dorsalis eggs. Briefly, 150 dechorionated embryos were transferred to a filter paper laid on an agar plate and moistened with a 1000 µL (1 × 108) inoculum of the respective bacteria species for each BMAL. For the axenic lines, distilled water was used to moisten the filter paper, whereas, for the symbiotic lines, 1000 µL of the liquid diet used to raise normal laboratory flies was used. The agar plates were incubated for one hour (to allow for efficient bacterial infection in mono-associated and symbiotic lines), after which the oviposition units were transferred to individual Perspex cages (12 × 12 × 12 cm) containing 10 mated F. arisanus females for parasitization, for six hours. Post-parasitization, B. dorsalis eggs were retrieved from the oviposition units and transferred to autoclaved liquid larval diets [55] (excluding streptomycin and nipagen) supplemented with 500 µL of the respective bacteria inoculum for the BMALs, a normal larval liquid diet for the symbiotic lines, or distilled water for the axenic lines and left to feed and subsequently pupariate. Consequently, the puparia were transferred to transparent 4-liter lunch boxes (18 × 11 × 15 cm) covered with a cotton mesh and monitored for fly and/or parasitoid emergence.
The parasitization experiments for each fly line and parasitoid were conducted individually in sterile conditions. All the wasps used in these assays underwent parasitization by exposing them to bacteria-free B. dorsalis larvae or eggs prior to parasitization assays. The wasps that did not exhibit oviposition behavior within 15 min were removed and replaced with new ones. All fly lines and parasitoids were kept at 25–27 °C, 60–70% humidity, and 12:12 day:light photoperiod.

2.5. Effect of Bacterial Symbionts on Host Acceptability for Oviposition by Parasitoid Females

To check for possible variations in acceptability by parasitoids across the different fly lines, host larvae and/or eggs from all fly lines were exposed to parasitoids, as described in Section 2.4, in no-choice assays and dissected to check for oviposition success. Diachasmimorpha longicaudata- and P. cosyrae-parasitized fly lines were dissected six hours post-parasitization, while those exposed to F. arisanus were dissected after four days. The host larvae containing parasitoid egg(s)/larva(e) were recorded and expressed as a percentage of the total number of exposed hosts to reflect host acceptability. These experiments were replicated five times.

2.6. Parasitoid Emergence and Fitness Assays

To investigate the effect of bacterial symbionts on parasitism success and fitness of the emerging parasitoids, the host fly lines were generated and parasitized, as described in Section 2.4. The parasitized fly lines were reared until pupation, and the pupae were collected, transferred to transparent four-liter lunch boxes (18 × 11 × 15 cm) covered with a cotton net cloth, and monitored for host fly and parasitoid emergence. The emerging host flies were collected, counted, and killed by freezing, while the parasitoids were transferred to separate transparent four-liter lunch boxes (18 × 11 × 15 cm) and fed on 50% honey solution ad libitum. The emerged wasps from each line for each parasitoid species were counted and sexed, and the number of each sex was recorded. On day 28 of post-exposure to parasitoid wasps, the host cadavers (unhatched host puparia) were dissected to check for the presence of dead parasitoid offspring. If found, these were recorded as parasitoid deaths.
The total number of emerging parasitoids was expressed as a percentage of the total recovered puparia to reflect parasitoid emergence. The total number of emerging host flies was expressed as a percentage of the total recovered puparia to reflect fly emergence. The total number of host puparia containing a dead parasitoid was expressed as a percentage of the total recovered puparia to reflect parasitoid death. Sex ratios were calculated as the proportion of female wasps out of all the emerged parasitoids.

2.7. Effect of Bacterial Symbionts Parasitoid Offspring’s Fitness Traits

2.7.1. Parasitoid Progeny Developmental Time

In a separate set of experiments, B. dorsalis larvae and eggs were exposed to D. longicaudata and F. arisanus for two and six hours, respectively, and treated in a similar manner as described in Section 2.4. The host puparia were transferred to four- liter lunch boxes and monitored until the onset of parasitoid emergence upon which the number of emerging male and female parasitoids was recorded daily until the last parasitoid emerged. These data were used to calculate the developmental time of the parasitoids, which was computed using the formula below [14]. These experiments were replicated 10 times.
Developmental   time = i = 1 n n D i n F i i = 1 n   n F i
where i signifies an individual of a total of n insects; nF i is the daily individual insect (parasitoid) emergence; and nD i is the duration in days for the development of the i th insect (parasitoid) from oviposition to adult emergence.

2.7.2. Parasitoid Progeny Body Size

The hind tibia of the F1 parasitoid progeny of D. longicaudata and F. arisanus were used as a proxy to measure the parasitoid progeny size. Briefly, one hind leg of each of the 10-day-old F1 parasitoids was removed, mounted on a microscope slide using entellanTM rapid mounting medium for microscopy (Merck, Darmstadt, Germany), and layered with coverslips. The hind tibiae (one tibia from each parasitoid) were measured under a stereomicroscope (Zeiss Stemi 508, Zeiss, Oberkochen, Germany). From each fly line, 10 individuals (five males and five females) of each parasitoid species were randomly selected for the assays, and this was replicated five times.

2.7.3. Longevity of Host-Deprived F1 Parasitoid Progeny

From the same groups of emerging parasitoids described in Section 2.7.1 above, five replicates were randomly selected for the parasitoid progeny survival analysis. From each selected replicate of each F1 parasitoid line, 10 male and 10 female parasitoids were randomly selected, transferred to sterile transparent four-liter lunch boxes (18 × 11 × 15 cm), and fed on 50% honey solution ad libitum, but they were not provided with hosts for oviposition. The parasitoids were monitored daily for mortality until the last parasitoid died. For each F1 progeny line within a parasitoid species, the survival data were obtained for 50 males and 50 females.

2.7.4. Fecundity of F1 Female Parasitoid Progeny

In a separate set of assays, 100 host larvae or 150 eggs from the axenic line, symbiotic line, and BMALs were exposed to 10 mated female D. longicaudata or F. arisanus for 2 and 6 hours, respectively, and incubated until parasitoid emergence, as described in Section 2.5. The emerging F1 parasitoids of each species from each fly line were used to investigate the symbiont-mediated effects on F1 parasitoid fecundity. Briefly, 10 randomly selected 5-day-old, mated females of D. longicaudata or F. arisanus from each fly line were transferred to four-liter lunch boxes (18 × 11 × 15 cm) and exposed to axenic larvae or eggs for two or six hours, respectively. The exposure was carried out once every day for 15 consecutive days. A 1:10 parasitoid: host ratio was maintained throughout the experiment. If a female parasitoid died, the number of host larvae/eggs offered daily was adjusted to maintain this ratio. The Diachasmimorpha longicaudata-parasitized larval groups were dissected 48 h after parasitization, while the F. arisanus-parasitized groups were dissected after four days, under a stereomicroscope (Zeiss Stemi 508, Zeiss, Oberkochen, Germany). The total number of parasitoid eggs and/or larvae in each host larva was recorded daily, and this was presented as the number of eggs laid per female parasitoid per day. These assays were performed in four independent replicates for each parasitoid species from each fly line.

2.8. Data Analysis

All analyses were performed using R software (version 4.0.1) [57].
The data obtained from host acceptability bioassays were analyzed using a generalized linear model (GLM) with quasi-Poisson distribution considering parasitoid species and host infection status as fixed factors. Data from the host suitability (parasitoid and fly emergence as well as parasitoid mortality) bioassay and the parasitoid progeny fitness trait (sex ratio, developmental time, body size, fecundity, and longevity) bioassay were analyzed separately for each parasitoid species. Since they were collected from the same batch of host fly lines, parasitoid and fly emergence, as well as parasitoid mortality data, were analyzed using multivariate analysis of variance (MANOVA), with host infection status as a fixed factor. In cases in which a significant effect of host infection status was detected, a univariate one-way analysis of variance (ANOVA) was conducted to check for significant differences within the response variables (parasitoid emergence, fly emergence, and parasitoid death), and this was followed by Student–Neuman–Keuls (SNK) post hoc multiple comparisons tests.
The sex ratios of D. longicaudata and F. arisanus F1 progeny were analyzed using a GLM with a binomial error distribution, considering host infection status as a fixed factor. Parasitoid progeny developmental time was analyzed using a GLM assuming a Gamma error distribution with host infection status as a fixed factor. Parasitoid progeny body size (hind tibia length) was analyzed using a linear mixed-effect model from the ‘lme4’ package (version 1.1.27.1) [58], considering host infection status and the sex of the parasitoid progeny as fixed and random factors, respectively. A generalized, linear mixed-effect model (using the ‘glmer’ function from the ‘lme4’ package [58]), with a Gamma error distribution, was fitted to analyze the fecundity of the female parasitoid progeny of F. arisanus and D. longicaudata, considering host infection status and replicate population as fixed and random factors, respectively. A Cox proportional hazards model using the “coxme” function from the “survival” package (version 3.2.13) [59] was fitted to determine the effect of bacterial symbionts on the longevity of the F1 progeny of both parasitoid species. In the model, both host infection status and replicate population were included as fixed and random factors, respectively. The significance of all these models was determined using analysis of deviance with chi-square tests, and mean separation was determined using Tukey’s multiple comparison tests at α ≤ 0.05.

3. Results

3.1. Effect of Bacterial Symbionts on Host Acceptability for Oviposition by Female Parasitoids

The no-choice assays revealed that host acceptance by all the parasitoid species was not affected by symbiotic bacteria (χ2 = 0.518, df = 3, p = 0.775), parasitoid species (χ2 = 0.246, df = 2, p = 0.769), or the interaction between these two factors (χ2 = 0.343, df = 6, p = 0.994) (Figure 1).

3.2. Effect of Bacterial Symbionts on Host Suitability for Development of the Immature Stages of the Parasitoids

We recorded a significant effect of symbiotic bacteria on parasitoid emergence, fly emergence, and parasitoid mortality in the F. arisanus-parasitized fly lines (MANOVA: F12,135 = 5.913, p < 0.001). Further, analysis using ANOVA showed significant effects of bacterial symbionts on parasitoid emergence (F4,45 = 24.122, p < 0.001; Figure 2a) as well as host fly emergence (F4,45 = 33.708, p < 0.001; Figure 2b) but not parasitoid mortality (F4,45 = 0.872, p = 0.488; Figure 2c). The flies inoculated with L. lactis yielded the highest percentage of parasitoid wasps and the lowest percentage of eclosed host flies, while the reverse was true for those inoculated with P. alcalifaciens (Figure 2a,b).
In general, the effect of bacterial symbionts on the three response variables in the D. longicaudata-parasitized host fly lines (MANOVA: F12,135 = 2.666, p = 0.003) showed a similar trend to that recorded for F. arisanus. Inoculating B. dorsalis with symbiotic bacteria significantly influenced parasitoid emergence (F4,45 = 8.246, p < 0.001; Figure 3a) and fly emergence (F4,45 = 11.527, p < 0.001; Figure 3b) but not parasitoid mortality (F4,45 = 0.304, p = 0.874; Figure 3c). Overall, the L. lactis-infected fly lines yielded the highest percentage of parasitoids when compared with other host fly lines (Figure 3a).
We did not record any parasitoid emergence from P. cosyrae-parasitized larval groups irrespective of the infection status. However, we found a significant multivariate effect of bacterial symbionts on host fly emergence and parasitoid mortality (MANOVA: F8,90 = 4.667, p < 0.001). Further analysis revealed that host fly emergence significantly varied across all the host fly lines (F4,45 = 13.419, p < 0.001; Figure 4a), with the highest percent fly emergence recorded in the symbiotic and P. alcalifaciens host fly lines (Figure 4a). Parasitoid mortality was not affected by symbiotic bacteria (F4,45 = 1.805, p = 0.145; Figure 4b).

3.3. Effect of Bacterial Symbionts on Parasitoid Offspring’s Fitness Traits

3.3.1. Parasitoid Progeny Sex Ratio

The bacterial symbionts of B. dorsalis did not alter the sex ratios of the progeny of F. arisanus (ꭓ2 = 0.048, df = 4, p = 1, Figure 5a), nor those of D. longicaudata (ꭓ2 = 0.078, df = 4, p = 0.999, Figure 5b). However, the sex ratio was female-biased for both parasitoid species (Figure 5a, b).

3.3.2. Parasitoid Progeny Developmental Time

Developmental times significantly varied among host fly lines for both parasitoid species (χ2 = 122.360, df = 4, p < 0.001; χ2 = 18.20, df = 4, p = 0.001, for F. arisanus and D. longicaudata, respectively), being the longest for the wasps emerging from P. alcalifaciens and L. lactis fly lines (Table 1). Overall, the average developmental time for F. arisanus was 23.56 ± 0.37 days, while it was 17.50 ± 0.21 days for D. longicaudata.

3.3.3. Parasitoid Progeny Body Size and Longevity

Like other fitness parameters, the body size (as measured by the hind tibia length) of the F1 progeny of both parasitoid species was significantly affected by bacterial symbionts (F. arisanus: χ2 = 56.220, df = 4, p < 0.001; D. longicaudata: χ2 = 345.750, df = 4, p < 0.001). Overall, the parasitoids emerging from L. lactis and P. alcalifaciens fly lines had longer tibia than those emerging from other fly lines for both parasitoids (Figure 6). On the other hand, there were no symbiont-mediated effects on the longevity of the parasitoid progeny of both parasitoid species (χ2 = 0.006, df = 4, p = 1.000 and χ2 = 0.096, df = 4, p = 0.998 for F. arisanus and D. longicaudata, respectively) (Figure 7a,b).

3.3.4. Fecundity of Female Parasitoid Progeny

The fecundity of the F1 progeny of both parasitoids was significantly affected by symbiotic bacteria (χ2 = 11.698, df = 4, p = 0.020 and χ2 = 60.834, df = 4, p < 0.001 for F. arisanus and D. longicaudata, respectively), with the parasitoid progeny emerging from the P. alcalifaciens-inoculated fly lines ovipositing more eggs (F. arisanus = 6.18 ± 0.29 and D. longicaudata = 6.74 ± 0.26 eggs/♀/day), whereas the least fecund parasitoid wasps were those reared from the L. lactis-inoculated fly lines (F. arisanus = 5.44 ± 0.16 and D. longicaudata = 5.30 ± 0.13 eggs/♀/day, Figure 8a,b).

4. Discussion

Augmentative biological control is a widely used approach for the management of fruit flies in many countries [11,12,13,16,20]. This approach does not only require the production of parasitoids on a large scale, but the produced parasitoid wasps also need to be of high quality in terms of fitness, to ensure their optimal performance in the field. Symbiotic bacteria could be exploited to enhance the quality of mass-reared parasitoids [41]. Indeed, symbiotic bacteria have been reported to influence the physiological quality of their hosts, inevitably influencing how they are perceived by parasitoids and their suitability for the development of parasitoid immatures and, eventually, the fitness of the emerging parasitoid progeny [39,41,42,60].
The comparable acceptability of the three parasitoid species (F. arisanus, P. cosyrae, and D. longicaudata) across the B. dorsalis fly lines reported in this study might be an indication that these symbiotic bacteria had no influence of the host cues used by the parasitoid females to locate their hosts. Indeed, Sochard et al. [27] reported that Hamiltonella defensa and Regiella insecticola did not alter the acceptability of Acyrthosiphon pisum (Harris) (Hemiptera: Aphididae) by its parasitoid, Aphidius ervi Haliday (Hymenoptera: Braconidae). Likewise, Xie et al. [32] showed that Spiroplasma sp. had no effect on the parasitism rates of Leptopilina heterotoma Förster (Hymenoptera: Figitidae) in Drosophila melanogaster Meigen (Diptera: Drosophilidae). However, the findings by Attia et al. [39] and Frago et al. [29] revealed that symbionts influence the host acceptability by parasitoid wasps, a phenotype that was attributed to the symbiont-mediated cues that influence host finding and ovipositing behavior by the parasitoids.
Nonetheless, despite the similar outcome in host acceptability, parasitoid emergence varied among the B. dorsalis fly lines, suggesting that bacterial symbionts had a profound effect on parasitoid development. Most facultative bacteria described to date possess protective roles where, in some instances, complete host protection against parasitoids is realized [30]. These protective phenotypes are largely attributed to the symbiont-encoded viruses that impair parasitoid development [35,36], the ability of the symbionts to alter the host immune defenses [37], or their ability to compete for nutrients, which are needed for the successful development of the parasitoids [38]. However, the high parasitism rates recorded in the L. lactis-inoculated fly lines in the current study suggest that L. lactis increases the susceptibility of B. dorsalis to attack by parasitoids. This susceptibility could be associated with the suppression of immune defense mechanisms of B. dorsalis and/or the L. lactis-mediated pathogenicity mechanism, previously highlighted in the B. dorsalisL. lactisMetarhizium anisopliae interaction [46]. However, this mechanism remains poorly understood.
Surprisingly, one evaluated symbiont, P. alcalifaciens, was found to confer protection to B. dorsalis against F. arisanus but not D. longicaudata, suggesting that bacteria-mediated protection against parasitoids is not universal across host–parasitoid systems. These two parasitoids attack different stages of their host, with F. arisanus using the eggs and D. longicaudata using the larvae as hosts, and this might have contributed to the discrepancy in P. alcalifaciens-mediated protection. Alternatively, it is possible that P. alcalifaciens produces toxins or engages pathogenic mechanisms to which only F. arisanus is susceptible. Similar symbiont-specific host-protective phenotypes have been demonstrated in other host–symbiont–parasitoid tripartite interactions [31,40,61,62].
Similar to earlier studies [40,42], we found that symbiotic bacteria did not affect the sex ratios of the progeny of F. arisanus and D. longicaudata. It has been demonstrated that parasitoids discriminate between good- and poor-quality hosts and, consequently, allocate more female offspring to good-quality hosts [63,64,65,66]. The results obtained in the current study showed that the F1 progeny of both parasitoid species from all fly lines were female-biased. The bias may be caused by the parasitoids perceiving hosts across all fly lines as good quality, thus choosing to oviposit more females in all of them. Furthermore, the female parasitoids used in this study were newly mated and in their early stages of oviposition, which may have contributed to this sex bias. Studies have shown that mated female parasitoids are likely to lay more females than males since most eggs at this stage are diploid (fertilized) eggs that develop into female progeny [67,68].
Earlier studies suggested that the parasitoid progeny’s developmental time and size are a function of host size [42,69], wherein the parasitoids emerging from smaller hosts take longer to develop, especially for koinobiont parasitoid species [70,71]. Indeed, the B. dorsalis pupae from the L. lactis and P. alcalifaciens fly lines were reported to have higher weights than those from the other fly lines [46]. It was, therefore, expected that the parasitoid progeny from these lines would emerge earlier. However, in this study, these parasitoids emerged later than those reared on the, C. freundii, axenic, and symbiotic host fly lines with associated lower puparia weights. This phenotype could be a symbiont-mediated manipulation of host physiology such that the parasitoids developing in these lines delay development and subsequent emergence until the hosts reach the desired physiological state. Alternatively, it is known that the parasitoid rate of the consumption of its host is inversely proportional to the size of the host, which may increase the development time of the parasitoids, especially in tissue-feeding endoparasitoids [69,70]. It is, therefore, possible that the parasitoids from the P. alcalifaciens and L. lactis fly lines had longer developmental times, as they had to completely consume their relatively larger hosts prior to emergence. These findings, however, contrast with a previous report that other Providencia species shorten the developmental time and alter the sex ratio of D. longicaudata progeny [41]. Differences in bacteria strains used, inoculation strategy, and host species may account for the contrast between our findings and the findings of that study, more so since in the latter, bacterial isolates were supplemented in non-axenic diets, thus accommodating additive effects from other native bacterial symbionts in the host larvae.
Moreover, the finding that the parasitoids emerging from the P. alcalifaciens and L. lactis host fly lines were larger than those reared on the axenic, symbiotic, and C. freundii counterparts conforms with the host–parasitoid size theory, according to which larger hosts are predicted to yield larger parasitoid progeny [69,70,71,72,73]. The larger body size of the parasitoid progeny from these lines could, therefore, be an additive effect of these two symbionts on the size of their hosts, a trait that could be attributed to symbiont-mediated effects on the nutritional resources of the hosts reared on these bacterial lines.
Previous studies suggested that the size of a female parasitoid is directly proportional to its fecundity and longevity, e.g., [72,74,75]. While this held true for the F1 parasitoid progeny emerging from the P. alcalifaciens fly lines, the opposite was recorded for those reared on the L. lactis fly lines, which, despite being larger than the parasitoids from the axenic, symbiotic, and C. freundii fly lines, were less fecund. This finding suggests the L. lactis-associated fitness costs in such a way that the female parasitoid progeny from this fly line allocate more resources towards growth and adult maintenance at the detriment of their reproductive potential. This symbiont may, therefore, have cryptic fitness costs on the reproductive potential of parasitoids that supersede the progeny size benefits. Similar symbiont-mediated tradeoffs in parasitoid fitness have been demonstrated in other host–symbiont–parasitoid systems [39,62,76].

5. Conclusions

This study reveals a mixed nature of interactions: a general synergism with the parasitoids with L. lactis and a more beneficial role to the host flies in the case of P. alcalifaciens, suggesting that symbiotic-bacteria could have significant consequences on host–parasitoid dynamics. In this regard, these findings underscore the previously unknown role of bacterial symbionts in mediating the host–parasitoid interactions in B. dorsalis by providing information vital for parasitoid mass rearing. For example, L. lactis can be harnessed as a probiotic in mass-rearing programs of D. longicaudata and F. arisanus, hence enhancing the performance of parasitoids in biological control programs targeting this pest. Nevertheless, future studies are needed to elucidate the mechanisms underlying these symbiont-mediated host–parasitoid interactions.

Author Contributions

Conceptualization, S.A.M., J.K.H., S.N. and R.G.; methodology and validation, R.G., J.G., S.N., C.W.W., J.K.H. and S.A.M.; investigation and data curation, R.G., E.E.M. and J.G.; formal analysis, R.G. and C.W.W.; writing—original draft preparation, R.G.; writing—review and editing, R.G., E.E.M., J.G., S.N., C.W.W., J.K.H. and S.A.M.; supervision, S.N., C.W.W., J.K.H. and S.A.M.; project administration and funding acquisition, S.A.M. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work received financial support from the International Development Research Centre (IDRC) and the Australian Centre for International Agricultural Research (ACIAR) grant number 109040; the Norwegian Agency for Development Cooperation (NORAD), the Section for Research, Innovation, and Higher Education grant number RAF-3058 KEN-18/0005; the Swedish International Development Cooperation Agency (Sida); the Swiss Agency for Development and Cooperation (SDC); the Australian Centre for International Agricultural Research (ACIAR); the Federal Democratic Republic of Ethiopia; and the Government of the Republic of Kenya. The views expressed herein do not necessarily reflect the official opinion of the donors.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author(s) on reasonable request.

Acknowledgments

The authors are grateful to the arthropod pathology unit, the SymbioVector team, as well as the African fruit fly project team at icipe for their technical support.

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. Ekesi, S.; Nderitu, P.W.; Rwomushana, I. Field Infestation, Life History and Demographic Parameters of the Fruit Fly Bactrocera invadens (Diptera: Tephritidae) in Africa. Bull. Entomol. Res. 2006, 96, 379–386. [Google Scholar] [CrossRef] [PubMed]
  2. Lux, S.A.; Copeland, R.S.; White, I.M.; Manrakhan, A.; Billah, M.K. A New Invasive Fruit Fly Species from the Bactrocera dorsalis (Hendel) Group Detected in East Africa. Int. J. Trop. Insect Sci. 2003, 23, 355–361. [Google Scholar] [CrossRef]
  3. Shi, Z.; Wang, L.; Zhang, H. Low Diversity Bacterial Community and the Trapping Activity of Metabolites from Cultivable Bacteria Species in the Female Reproductive System of the Oriental Fruit Fly, Bactrocera dorsalis Hendel (Diptera: Tephritidae). Int. J. Mol. Sci. 2012, 13, 6266. [Google Scholar] [CrossRef]
  4. Wang, H.; Jin, L.; Peng, T.; Zhang, H.; Chen, Q.; Hua, Y. Identification of Cultivable Bacteria in the Intestinal Tract of Bactrocera dorsalis from Three Different Populations and Determination of Their Attractive Potential. Pest Manag. Sci. 2014, 70, 80–87. [Google Scholar] [CrossRef]
  5. Cheng, D.; Guo, Z.; Riegler, M.; Xi, Z.; Liang, G.; Xu, Y. Gut Symbiont Enhances Insecticide Resistance in a Significant Pest, the Oriental Fruit Fly Bactrocera dorsalis (Hendel). Microbiome 2017, 5, 13. [Google Scholar] [CrossRef] [PubMed]
  6. Rwomushana, I.; Tanga, C.M. Fruit Fly Species Composition, Distribution and Host Plants with Emphasis on Mango-Infesting Species. In Fruit Fly Research and Development in Africa-Towards a Sustainable Management Strategy to Improve Horticulture; Ekesi, S., Mohamed, S.A., De Meyer, M., Eds.; Springer: Cham, Switzerland, 2016; pp. 71–106. ISBN 9783319432267. [Google Scholar]
  7. Heve, W.K.; Adjadeh, T.A.; Billah, M.K. Overview and Future Research Needs for Development of Effective Biocontrol Strategies for Management of Bactrocera dorsalis Hendel (Diptera: Tephritidae) in Sub-Saharan Africa. Pest Manag. Sci. 2021, 77, 4224–4237. [Google Scholar]
  8. Ndlela, S.; Niassy, S.; Mohamed, S.A. Important Alien and Potential Native Invasive Insect Pests of Key Fruit Trees in Sub-Saharan Africa: Advances in Sustainable Pre- and Post-Harvest Management Approaches. CABI Agric. Biosci. 2022 31 2022, 3, 1–46. [Google Scholar] [CrossRef]
  9. Mohamed, S.A.; Ekesi, S.; Hanna, R. Evaluation of the Impact of Diachasmimorpha Longicaudata on Bactrocera invadens and Five African Fruit Fly Species. J. Appl. Entomol. 2008, 132, 789–797. [Google Scholar] [CrossRef]
  10. Mohamed, S.A.; Ekesi, S.; Wharton, R.A.; Lux, S.A.; Overholt, W.A. Old and New Host-Parasitoid Associations: Parasitism of the Native African and Invasive Fruit Flies Species. In Proceedings of the Proceedings of the 5th International Symposium on Biological Control of Arthropods; Mason, P.G., Gillespie, D.R., Vincent, C., Eds.; CABI: Wallingford, UK, 2017; pp. 224–227. [Google Scholar]
  11. Vargas, R.; Stark, J.D.; Uchida, G.K.; Purcell, M. Opiine Parasitoids (Hymenoptera: Braconidae) of Oriental Fruit Fly (Diptera: Tephritidae) on Kauai Island, Hawaii: Islandwide Relative Abundance and Parasitism Rates in Wild and Orchard Guava Habitats. Environ. Entomol. 1993, 22, 246–253. [Google Scholar] [CrossRef]
  12. Vargas, R.; Leblanc, L.; Harris, E.J.; Manoukis, N.C. Regional Suppression of Bactrocera Fruit Flies (Diptera: Tephritidae) in the Pacific through Biological Control and Prospects for Future Introductions into Other Areas of the World. Insects 2012, 3, 727–742. [Google Scholar] [CrossRef] [PubMed]
  13. Vargas, R.I.; Leblanc, L.; Putoa, R.; Eitam, A. Impact of Introduction of Bactrocera dorsalis (Diptera: Tephritidae) and Classical Biological Control Releases of Fopius Arisanus (Hymenoptera: Braconidae) on Economically Important Fruit Flies in French Polynesia. J. Econ. Entomol 2007, 100, 670–679. [Google Scholar] [CrossRef] [PubMed]
  14. Nanga Nanga, S.; Hanna, R.; Gnanvossou, D.; Fotso Kuate, A.; Fiaboe, K.K.M.; Djieto-Lordon, C.; Schmidt-Jeffris, R. Fruit Preference, Parasitism, and Offspring Fitness of Fopius arisanus (Hymenoptera: Braconidae) Exposed to Bactrocera dorsalis (Diptera: Tephritidae) Infested Fruit Species. Environ. Entomol. 2019, 48, 1286–1296. [Google Scholar] [CrossRef]
  15. Ndlela, S.; Mohamed, S.A.; Azrag, A.G.A.; Ndegwa, P.N.; Ong’amo, G.O.; Ekesi, S. Interactions between Two Parasitoids of Tephritidae: Diachasmimorpha longicaudata (Ashmead) and Psyttalia cosyrae (Wilkinson) (Hymenoptera: Braconidae), under Laboratory Conditions. Insects 2020, 11, 671. [Google Scholar] [CrossRef]
  16. Gnanvossou, D.; Hanna, R.; Bokonon-Ganta, A.H.; Ekesi, S.; Mohamed, S.A. Release, Establishment and Spread of the Natural Enemy Fopius arisanus (Hymenoptera: Braconidae) for Control of the Invasive Oriental Fruit Fly Bactrocera dorsalis (Diptera: Tephritidae) in Benin, West Africa. In Fruit Fly Research and Development in Africa—Towards a Sustainable Management Strategy to Improve Horticulture; Springer: Cham, Switzerland, 2016; pp. 575–600. [Google Scholar] [CrossRef]
  17. Leblanc, L.; Vargas, R.; Putoa, R. From Eradication to Containment: Invasion of French Polynesia by Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) and Releases of Two Natural Enemies: A 17-Year Case Study. In Proceedings of the Proceedings of the Hawaiian Entomological Society; Hawaiian Entomological Society: Honolulu, HI, USA, 2013; pp. 31–43. [Google Scholar]
  18. Mohamed, S.A.; Ekesi, S.; Hanna, R. Old and New Host-Parasitoid Associations: Parasitism of the Invasive Fruit Fly Bactrocera invadens (Diptera: Tephritidae) and Five African Fruit Fly Species by Fopius arisanus, an Asian Opiine Parasitoid. Biocontrol Sci. Technol. 2010, 20, 183–196. [Google Scholar] [CrossRef]
  19. Vargas, R.; Leblanc, L.; McKenney, M.; Mackey, B.; Harris, E.; Badji, K. Rearing Fopius arisanus (Sonan) (Hymenoptera: Braconidae) on Mediterranean Fruit Fly and Its Introduction into Senegal against Oriental Fruit Fly (Diptera: Tephritidae). Proc. Hawaiian Entomol. Soc. 2016, 48, 85–94. [Google Scholar]
  20. Mohamed, S.A.; Ramadan, M.M.; Ekesi, S. In and out of Africa: Parasitoids Used for Biological Control of Fruit Flies. In Fruit Fly Research and Development in Africa-Towards a Sustainable Management Strategy to Improve Horticulture; Ekesi, S., Ed.; Springer International Publishing: New York, NY, USA, 2016; pp. 325–368. ISBN 9783319432267. [Google Scholar]
  21. da Silva Gonçalves, R.; Manoukis, N.C.; Nava, D.E. Effect of Fopius arisanus Oviposition Experience on Parasitization of Bactrocera dorsalis. BioControl 2017, 62, 595–602. [Google Scholar] [CrossRef]
  22. Gwokyalya, R.; Herren, J.K.; Weldon, C.W.; Khamis, F.M.; Ndlela, S.; Mohamed, S.A. Differential Immune Responses in New and Old Fruit Fly-Parasitoid Associations: Implications for Their Management. Front. Physiol. 2022, 13, 1663. [Google Scholar] [CrossRef]
  23. Oliver, K.M.; Russell, J.A.; Moran, N.A.; Hunter, M.S. Facultative Bacterial Symbionts in Aphids Confer Resistance to Parasitic Wasps. Proc. Natl. Acad. Sci. USA 2003, 100, 1803–1807. [Google Scholar] [CrossRef]
  24. Oliver, K.M.; Moran, N.A.; Hunter, M.S. Variation in Resistance to Parasitism in Aphids is Due to Symbionts Not Host Genotype. Proc. Natl. Acad. Sci. USA 2005, 102, 12795–12800. [Google Scholar] [CrossRef]
  25. Oliver, K.M.; Noge, K.; Huang, E.M.; Campos, J.M.; Becerra, J.X.; Hunter, M.S. Parasitic Wasp Responses to Symbiont-Based Defense in Aphids. BMC Biol. 2012, 10, 11. [Google Scholar] [CrossRef]
  26. Asplen, M.K.; Bano, N.; Brady, C.M.; Desneux, N.; Hopper, K.R.; Malouines, C.; Oliver, K.M.; White, J.A.; Heimpel, G.E. Specialisation of Bacterial Endosymbionts That Protect Aphids from Parasitoids. Ecol. Entomol. 2014, 39, 736–739. [Google Scholar] [CrossRef]
  27. Sochard, C.; Bellec, L.; Simon, J.-C.; Outreman, Y. Influence of “Protective” Symbionts throughout the Different Steps of an Aphid–Parasitoid Interaction. Curr. Zool. 2020, 67, 1–13. [Google Scholar] [CrossRef]
  28. Leclair, M.; Buchard, C.; Mahéo, F.; Simon, J.C.; Outreman, Y. A Link Between Communities of Protective Endosymbionts and Parasitoids of the Pea Aphid Revealed in Unmanipulated Agricultural Systems. Front. Ecol. Evol. 2021, 9, 187. [Google Scholar]
  29. Frago, E.; Mala, M.; Weldegergis, B.T.; Yang, C.; McLean, A.; Godfray, H.C.J.; Gols, R.; Dicke, M. Symbionts Protect Aphids from Parasitic Wasps by Attenuating Herbivore-Induced Plant Volatiles. Nat. Commun. 2017, 8, 1860. [Google Scholar] [CrossRef]
  30. Guay, J.F.; Boudreault, S.; Michaud, D.; Cloutier, C. Impact of Environmental Stress on Aphid Clonal Resistance to Parasitoids: Role of Hamiltonella defensa Bacterial Symbiosis in Association with a New Facultative Symbiont of the Pea Aphid. J. Insect Physiol. 2009, 55, 919–926. [Google Scholar] [CrossRef]
  31. McLean, A.H.C.; Godfray, H.C.J. Evidence for Specificity in Symbiont-Conferred Protection against Parasitoids. Proc. R. Soc. B Biol. Sci. 2015, 282, 1–8. [Google Scholar] [CrossRef]
  32. Xie, J.; Butler, S.; Sanchez, G.; Mateos, M. Male Killing Spiroplasma Protects Drosophila melanogaster against Two Parasitoid Wasps. Heredity (Edinb). 2014, 112, 399–408. [Google Scholar] [CrossRef]
  33. Xie, J.; Vilchez, I.; Mateos, M. Spiroplasma Bacteria Enhance Survival of Drosophila hydei Attacked by the Parasitic Wasp Leptopilina heterotoma. PLoS ONE 2010, 5, e12149. [Google Scholar] [CrossRef] [PubMed]
  34. Degnan, P.H.; Moran, N.A. Evolutionary Genetics of a Defensive Facultative Symbiont of Insects: Exchange of Toxin-Encoding Bacteriophage. Mol. Ecol. 2008, 17, 916–929. [Google Scholar] [CrossRef]
  35. Moran, N.A.; Degnan, P.H.; Santos, S.R.; Dunbar, H.E.; Ochman, H. The Players in a Mutualistic Symbiosis: Insects, Bacteria, Viruses, and Virulence Genes. Proc. Natl. Acad. Sci. USA 2005, 102, 16919–16926. [Google Scholar] [CrossRef]
  36. Oliver, K.M.; Degnan, P.H.; Hunter, M.S.; Moran, N.A. Bacteriophages Encode Factors Required for Protection in a Symbiotic Mutualism. Science 2009, 325, 992–994. [Google Scholar] [CrossRef]
  37. Fytrou, A.; Schofield, P.G.; Kraaijeveld, A.R.; Hubbard, S.F. Wolbachia Infection Suppresses Both Host Defence and Parasitoid Counter-Defence. Proc. R. Soc. B Biol. Sci. 2005, 273, 791–796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Paredes, J.C.; Herren, J.K.; Schüpfer, F.; Lemaitre, B. The Role of Lipid Competition for Endosymbiont-Mediated Protection against Parasitoid Wasps in Drosophila. MBio 2016, 7. [Google Scholar] [CrossRef]
  39. Attia, S.; Renoz, F.; Pons, I.; Louâpre, P.; Foray, V.; Piedra, J.-M.; Sanané, I.; Le Goff, G.; Lognay, G.; Hance, T. The Aphid Facultative Symbiont Serratia symbiotica Influences the Foraging Behaviors and the Life-History Traits of the Parasitoid Aphidius Ervi. Entomol. Gen. 2021. [Google Scholar] [CrossRef]
  40. Hopper, K.R.; Kuhn, K.L.; Lanier, K.; Rhoades, J.H.; Oliver, K.M.; White, J.A.; Asplen, M.K.; Heimpel, G.E. The Defensive Aphid Symbiont Hamiltonella defensa Affects Host Quality Differently for Aphelinus glycinis versus Aphelinus atriplicis. Biol. Control 2018, 116, 3–9. [Google Scholar]
  41. Koskinioti, P.; Ras, E.; Augustinos, A.A.; Beukeboom, L.W.; Mathiopoulos, K.D.; Caceres, C.; Bourtzis, K. The Impact of Fruit Fly Gut Bacteria on the Rearing of the Parasitic Wasp Diachasmimorpha longicaudata. Entomol. Exp. Appl. 2020, 168, 541–559. [Google Scholar] [CrossRef]
  42. Nyabuga, F.N.; Outreman, Y.; Simon, J.C.; Heckel, D.G.; Weisser, W.W. Effects of Pea Aphid Secondary Endosymbionts on Aphid Resistance and Development of the Aphid Parasitoid Aphidius ervi: A Correlative Study. Entomol. Exp. Appl. 2010, 136, 243–253. [Google Scholar] [CrossRef]
  43. Pons, I.; Renoz, F.; Noël, C.; Hance, T. New Insights into the Nature of Symbiotic Associations in Aphids: Infection Process, Biological Effects, and Transmission Mode of Cultivable Serratia symbiotica Bacteria. Appl. Environ. Microbiol. 2019, 85. [Google Scholar] [CrossRef]
  44. Schmid, M.; Sieber, R.; Zimmermann, Y.S.; Vorburger, C. Development, Specificity and Sublethal Effects of Symbiont-Conferred Resistance to Parasitoids in Aphids. Funct. Ecol. 2012, 26, 207–215. [Google Scholar] [CrossRef]
  45. Andongma, A.A.; Wan, L.; Dong, Y.C.; Li, P.; Desneux, N.; White, J.A.; Niu, C.Y. Pyrosequencing Reveals a Shift in Symbiotic Bacteria Populations across Life Stages of Bactrocera dorsalis. Sci. Rep. 2015, 5, 9470. [Google Scholar] [CrossRef]
  46. Gichuhi, J.; Khamis, F.; Van den Berg, J.; Mohamed, S.; Ekesi, S.; Herren, J.K. Influence of Inoculated Gut Bacteria on the Development of Bactrocera dorsalis and on Its Susceptibility to the Entomopathogenic Fungus, Metarhizium anisopliae. BMC Microbiol. 2020, 20, 321. [Google Scholar] [CrossRef] [PubMed]
  47. Damodaram, K.J.P.; Ayyasamy, A.; Kempraj, V. Commensal Bacteria Aid Mate-Selection in the Fruit Fly, Bactrocera dorsalis. Microb. Ecol. 2016, 72, 725–729. [Google Scholar] [CrossRef]
  48. Akami, M.; Ren, X.-M.; Qi, X.; Mansour, A.; Gao, B.; Cao, S.; Niu, C.-Y. Symbiotic Bacteria Motivate the Foraging Decision and Promote Fecundity and Survival of Bactrocera dorsalis (Diptera: Tephritidae). BMC Microbiol. 2019, 19, 1–13. [Google Scholar] [CrossRef] [Green Version]
  49. Akami, M.; Andongma, A.A.; Zhengzhong, C.; Nan, J.; Khaeso, K.; Jurkevitch, E.; Niu, C.-Y.Y.; Yuval, B. Intestinal Bacteria Modulate the Foraging Behavior of the Oriental Fruit Fly Bactrocera dorsalis (Diptera: Tephritidae). PLoS ONE 2019, 14, e0210109. [Google Scholar] [CrossRef]
  50. Khan, M.; Seheli, K.; Bari, M.A.; Sultana, N.; Khan, S.A.; Sultana, K.F.; Hossain, M.A. Potential of a Fly Gut Microbiota Incorporated Gel-Based Larval Diet for Rearing Bactrocera dorsalis (Hendel). BMC Biotechnol. 2019, 19, 94. [Google Scholar] [CrossRef]
  51. Khaeso, K.; Andongma, A.A.; Akami, M.; Souliyanonh, B.; Zhu, J.; Krutmuang, P.; Niu, C.-Y. Assessing the Effects of Gut Bacteria Manipulation on the Development of the Oriental Fruit Fly, Bactrocera dorsalis (Diptera; Tephritidae). Symbiosis 2017, 74, 97–105. [Google Scholar] [CrossRef]
  52. Zhang, Q.; Cai, P.; Wang, B.; Liu, X.; Lin, J.; Hua, R.; Zhang, H.; Yi, C.; Song, X.; Ji, Q.; et al. Manipulation of Gut Symbionts for Improving the Sterile Insect Technique: Quality Parameters of Bactrocera dorsalis (Diptera: Tephritidae) Genetic Sexing Strain Males After Feeding on Bacteria-Enriched Diets. J. Econ. Entomol. 2021, 114, 560–570. [Google Scholar] [CrossRef]
  53. Ekesi, S.; Mohamed, S.A. Mass Rearing and Quality Control Parameters for Tephritid Fruit Flies of Economic Importance in Africa. In Wide Spectra of Quality Control; Akyar, I., Ed.; IntechOpen: London, UK, 2011; pp. 387–410. ISBN 978-953-307-683-6. [Google Scholar]
  54. Mohamed, S.A.; Overholt, W.A.; Wharton, R.A.; Lux, S.A.; Eltoum, E.M. Host Specificity of Psyttalia cosyrae (Hymenoptera: Braconidae) and the Effect of Different Host Species on Parasitoid Fitness. Biol. Control 2003, 28, 155–163. [Google Scholar] [CrossRef]
  55. Chang, C.L. Fruit Fly Liquid Larval Diet Technology Transfer and Update. J. Appl. Entomol. 2009, 133, 164–173. [Google Scholar] [CrossRef]
  56. Ekesi, S.; Mohamed, S.A.; Chang, C.L. A Liquid Larval Diet for Rearing Bactrocera invadens and Ceratitis fasciventris (Diptera: Tephritidae). Int. J. Trop. Insect Sci. 2014, 34, S90–S98. [Google Scholar] [CrossRef]
  57. R Core Team R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. 2021. Available online: https://www.gbif.org/tool/81287/r-a-language-and-environment-for-statistical-computing (accessed on 8 November 2022).
  58. Bates, D.; Mächler, M.; Bolker, B.M.; Walker, S.C. Fitting Linear Mixed-Effects Models Using Lme4. J. Stat. Softw. 2015, 67. [Google Scholar] [CrossRef]
  59. Therneau, T.M. A Package for Survival Analysis in R. Version 2.38. 2015. Available online: https://CRAN.R-project.org/package=survival (accessed on 8 November 2022).
  60. Desneux, N.; Asplen, M.K.; Brady, C.M.; Heimpel, G.E.; Hopper, K.R.; Luo, C.; Monticelli, L.; Oliver, K.M.; White, J.A. Intraspecific Variation in Facultative Symbiont Infection among Native and Exotic Pest Populations: Potential Implications for Biological Control. Biol. Control 2018, 116, 27–35. [Google Scholar] [CrossRef]
  61. Cayetano, L.; Vorburger, C. Symbiont-Conferred Protection against Hymenopteran Parasitoids in Aphids: How General Is It? Ecol. Entomol. 2015, 40, 85–93. [Google Scholar] [CrossRef] [Green Version]
  62. Luo, C.; Gatti, J.L.; Monticelli, L.S.; Poirié, M.; Desneux, N.; Zhao, H.; Hu, Z. An Increased Risk of Parasitism Mediated by the Facultative Symbiont Regiella insecticola. J. Pest Sci. 2020, 93, 737–745. [Google Scholar] [CrossRef]
  63. Charnov, E.L.; Los-Den Hartogh, R.L.; Jones, W.T.; Van Den Assem, J. Sex Ratio Evolution in a Variable Environment. Nature 1981, 289, 27–33. [Google Scholar] [CrossRef]
  64. Charnov, E.L. The Theory of Sex Allocation; Princeton University Press: Princeton, NJ, USA, 1982; Volume 18, ISBN 9780691083124. [Google Scholar]
  65. King, B.H.; Napoleon, M.E. Using Effects of Parasitoid Size on Fitness to Test a Host Quality Model Assumption with the Parasitoid Wasp Spalangia endius. Can. J. Zool. 2006, 84, 1678–1682. [Google Scholar] [CrossRef]
  66. King, B.H. Host-Size-Dependent Sex Ratios among Parasitoid Wasps: Does Host Growth Matter ? Oecologia 1989, 78, 420–426. [Google Scholar]
  67. Ode, P.J.; Hardy, I.C.W. Parasitoid Sex Ratios and Biological Control. In Behavioral Ecology of Insect Parasitoids: From Theoretical Approaches to Field Applications; Wiley-Blackwell: Hoboken, NJ, USA, 2008; pp. 253–291. ISBN 9781405163477. [Google Scholar]
  68. Harvey, J.A.; de Haan, L.; Verdeny-Vilalta, O.; Visser, B.; Gols, R. Reproduction and Offspring Sex Ratios Differ Markedly among Closely Related Hyperparasitoids Living in the Same Microhabitats. J. Insect Behav. 2019, 32, 243–251. [Google Scholar] [CrossRef]
  69. Arthur, A.P.; Wylie, H.G. Effects of Host Size on Sex Ration, Development Time and Size of Pimpla turionellae (L.) (Hymenoptera: Ichneumonidae). Entomophaga 1959, 4, 297–301. [Google Scholar] [CrossRef]
  70. Harvey, J.A.; Kadash, K.; Strand, R. Micheal Differences in Larval Feeding Behavior Correlate with Altered Developmental Strategies in Two Parasitic Wasps: Implications for the Size-Fitness Hypothesis. Oikos 2000, 88, 621–629. [Google Scholar]
  71. Mackauer, M.; Sequeira, R. Patterns of Development in Insect Parasites. In Parasites and Pathogens of Insects; Beckage, N.E., Thompson, S.N., Federici, B.A., Eds.; Academic Press: San Diego, CA, USA, 1993; Volume 1, p. 12. ISBN 9780120844418. [Google Scholar]
  72. Godfray, H.C.J. Parasitoids: Behavioral and Evolutionary Ecology, 1st ed.Princeton University Press: Princeton, NJ, USA, 1994; Volume 67, ISBN 9780691000473. [Google Scholar]
  73. Wang, X.G.; Messing, R.H. Fitness Consequences of Body-Size-Dependent Host Species Selection in a Generalist Ectoparasitoid. Behav. Ecol. Sociobiol. 2004, 56, 513–522. [Google Scholar] [CrossRef]
  74. Gao, S.; Tang, Y.; Wei, K.; Wang, X.; Yang, Z.; Zhang, Y. Relationships between Body Size and Parasitic Fitness and Offspring Performance of Sclerodermus pupariae Yang et Yao (Hymenoptera: Bethylidae). PLoS ONE 2016, 11, e0156831. [Google Scholar] [CrossRef] [PubMed]
  75. Sagarra, L.A.A.; Vincent, C.; Stewart, R.K.K. Body Size as an Indicator of Parasitoid Quality in Male and Female Anagyrus kamali (Hymenoptera: Encyrtidae). Bull. Entomol. Res. 2001, 91, 363–367. [Google Scholar] [CrossRef] [PubMed]
  76. Luo, C.; Monticelli, L.; Meng, L.; Li, D.; Fan, J.; Zhao, H.; Hu, Z. Effect of the Endosymbiont Regiella insecticola on an Aphid Parasitoid. Entomol. Gen. 2017, 36, 301–308. [Google Scholar] [CrossRef]
Figure 1. Acceptability by three parasitoid species (Diachasmimorpha longicaudata, Fopius arisanus, and Psyttalia cosyrae) for Bactrocera dorsalis fly lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) (mean ± standard error of the mean (SEM)). n.s denotes non-significant differences (Tukey’s tests, α ≤ 0.05, n = 500).
Figure 1. Acceptability by three parasitoid species (Diachasmimorpha longicaudata, Fopius arisanus, and Psyttalia cosyrae) for Bactrocera dorsalis fly lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) (mean ± standard error of the mean (SEM)). n.s denotes non-significant differences (Tukey’s tests, α ≤ 0.05, n = 500).
Biology 12 00274 g001
Figure 2. Effects of symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the outcome of Bactrocera dorsalis and Fopius arisanus interaction. Within each category ((a) parasitoids, (b) host flies, and (c) parasitoid mortality), boxes capped with different letters indicate significant statistical differences; n.s denotes nonsignificant differences (Student–Neuman–Keuls (SNK) test, α ≤ 0.05, n = 1000). Within each box, horizontal bars denote the median values, and the ends of each boxplot whisker represent the minimum and maximum values of the data. Each data point represents the percent individuals (emerging parasitoids and host flies as well as dead parasitoids) recorded for a single replicate.
Figure 2. Effects of symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the outcome of Bactrocera dorsalis and Fopius arisanus interaction. Within each category ((a) parasitoids, (b) host flies, and (c) parasitoid mortality), boxes capped with different letters indicate significant statistical differences; n.s denotes nonsignificant differences (Student–Neuman–Keuls (SNK) test, α ≤ 0.05, n = 1000). Within each box, horizontal bars denote the median values, and the ends of each boxplot whisker represent the minimum and maximum values of the data. Each data point represents the percent individuals (emerging parasitoids and host flies as well as dead parasitoids) recorded for a single replicate.
Biology 12 00274 g002
Figure 3. Effects of symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the outcome of Bactrocera dorsalis and Diachasmimorpha longicaudata interaction. Within each category ((a) parasitoids, (b) host flies, and (c) parasitoid mortality), boxes capped with different letters indicate significant statistical differences; n.s denotes nonsignificant differences (Student–Neuman–Keuls (SNK) test, α ≤ 0.05, n = 1000). Within each box, horizontal bars denote the median values and the ends of each boxplot whisker represent the minimum and maximum values of the data. Each data point represents the percent individuals (emerging parasitoids and host flies as well as dead parasitoids) recorded for a single replicate.
Figure 3. Effects of symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the outcome of Bactrocera dorsalis and Diachasmimorpha longicaudata interaction. Within each category ((a) parasitoids, (b) host flies, and (c) parasitoid mortality), boxes capped with different letters indicate significant statistical differences; n.s denotes nonsignificant differences (Student–Neuman–Keuls (SNK) test, α ≤ 0.05, n = 1000). Within each box, horizontal bars denote the median values and the ends of each boxplot whisker represent the minimum and maximum values of the data. Each data point represents the percent individuals (emerging parasitoids and host flies as well as dead parasitoids) recorded for a single replicate.
Biology 12 00274 g003
Figure 4. Effects of symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the outcome of Bactrocera dorsalis and Psyttalia cosyrae interaction. Within each category ((a) host flies and (b) parasitoid mortality), boxes capped with different letters indicate significant statistical differences; n.s denotes nonsignificant differences (Student–Neuman–Keuls (SNK) test, α ≤ 0.05, n = 1000). Within each box, horizontal bars denote the median values and the ends of each boxplot whisker represent the minimum and maximum values of the data. Each data point represents the percent individuals (emerging host flies as well as dead parasitoids) recorded for a single replicate.
Figure 4. Effects of symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the outcome of Bactrocera dorsalis and Psyttalia cosyrae interaction. Within each category ((a) host flies and (b) parasitoid mortality), boxes capped with different letters indicate significant statistical differences; n.s denotes nonsignificant differences (Student–Neuman–Keuls (SNK) test, α ≤ 0.05, n = 1000). Within each box, horizontal bars denote the median values and the ends of each boxplot whisker represent the minimum and maximum values of the data. Each data point represents the percent individuals (emerging host flies as well as dead parasitoids) recorded for a single replicate.
Biology 12 00274 g004
Figure 5. Sex ratios (mean ± SEM) of (a) Fopius arisanus and (b) Diachasmimorpha longicaudata progeny emerging from different Bactrocera dorsalis lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic). n.s denotes non-significant differences among the sex ratios of the emerging parasitoids (Tukey’s tests, α = 0.05).
Figure 5. Sex ratios (mean ± SEM) of (a) Fopius arisanus and (b) Diachasmimorpha longicaudata progeny emerging from different Bactrocera dorsalis lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic). n.s denotes non-significant differences among the sex ratios of the emerging parasitoids (Tukey’s tests, α = 0.05).
Biology 12 00274 g005
Figure 6. Effect of Bactrocera dorsalis symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on F1 parasitoid body size (hind tibia lengths (mm; mean ± standard error of the mean (SEM)). Within each parasitoid species, scatter plots capped with different letters are significantly different (Tukey’s tests, α = 0.05, n = 50). Each data point represents the tibia length of a single adult parasitoid.
Figure 6. Effect of Bactrocera dorsalis symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on F1 parasitoid body size (hind tibia lengths (mm; mean ± standard error of the mean (SEM)). Within each parasitoid species, scatter plots capped with different letters are significantly different (Tukey’s tests, α = 0.05, n = 50). Each data point represents the tibia length of a single adult parasitoid.
Biology 12 00274 g006
Figure 7. Effect of Bactrocera dorsalis symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the survival of F1 progeny of (a) Fopius arisanus and (b) Diachasmimorpha longicaudata (Tukey’s tests, α = 0.05, n = 100). The data represent 100 (50 males and 50 females) F. arisanus and D. longicaudata F1 offspring kept in four-liter lunch boxes and fed with 50% honey solution ad libitum. Survival of both parasitoid species was recorded until the last parasitoid died.
Figure 7. Effect of Bactrocera dorsalis symbiotic bacteria (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) on the survival of F1 progeny of (a) Fopius arisanus and (b) Diachasmimorpha longicaudata (Tukey’s tests, α = 0.05, n = 100). The data represent 100 (50 males and 50 females) F. arisanus and D. longicaudata F1 offspring kept in four-liter lunch boxes and fed with 50% honey solution ad libitum. Survival of both parasitoid species was recorded until the last parasitoid died.
Biology 12 00274 g007
Figure 8. Cumulative fecundity of F1 female parasitoid progeny of (a) Fopius arisanus and (b) Diachasmimorpha longicaudata reared on different Bactrocera dorsalis fly lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) (Tukey’s tests, α = 0.05, n = 400). Error bars reflect the standard error of the mean (SEM).
Figure 8. Cumulative fecundity of F1 female parasitoid progeny of (a) Fopius arisanus and (b) Diachasmimorpha longicaudata reared on different Bactrocera dorsalis fly lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic) (Tukey’s tests, α = 0.05, n = 400). Error bars reflect the standard error of the mean (SEM).
Biology 12 00274 g008
Table 1. Effect of bacterial symbionts on the developmental times of Fopius arisanus and Diachasmimorpha longicaudata progeny emerging from the different Bactrocera dorsalis fly lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic). Within each column, means followed by the same letter are not significantly different (Tukey’s tests, α = 0.05, n = 1000). SEM = standard error of the mean.
Table 1. Effect of bacterial symbionts on the developmental times of Fopius arisanus and Diachasmimorpha longicaudata progeny emerging from the different Bactrocera dorsalis fly lines (axenic, Lactococcus lactis, Providencia alcalifaciens, Citrobacter freundii, and symbiotic). Within each column, means followed by the same letter are not significantly different (Tukey’s tests, α = 0.05, n = 1000). SEM = standard error of the mean.
Developmental Time Days (Mean ± SEM)
Parasitoid SpeciesFly LineMaleFemalePooled
Fopius arisanusAxenic19.96 ± 0.109 a21.39 ± 0.118 a20.76 ± 0.133 ab
L. lactis20.30 ± 0.147 ab21.86 ± 0.083 b21.05 ± 0.174 bc
P. alcalifaciens20.69 ± 0.077 b21.80 ± 0.074 b21.16 ± 0.210 c
C. freundii19.99 ± 0.115 a21.29 ± 0.122 a20.78 ± 0.086 ab
Symbiotic19.93 ± 0.078 a21.12 ± 0.084 a20.55 ± 0.118 a
Diachasmimorpha longicaudataAxenic14.97 ± 0.977 a17.65 ± 0.096 a16.97 ± 0.107 a
L. lactis16.91 ± 0.220 b18.57 ± 0.243 c17.89 ± 0.224 b
P. alcalifaciens16.75 ± 0.101 ab18.18 ± 0.048 bc17.61 ± 0.052 b
C. freundii15.81 ± 0.113 ab17.56 ± 0.096 a16.81 ± 0.106 a
Symbiotic16.08 ± 0.117 ab17.72 ± 0.097 ab16.99 ± 0.096 a
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gwokyalya, R.; Weldon, C.W.; Herren, J.K.; Gichuhi, J.; Makhulu, E.E.; Ndlela, S.; Mohamed, S.A. Friend or Foe: Symbiotic Bacteria in Bactrocera dorsalis–Parasitoid Associations. Biology 2023, 12, 274. https://doi.org/10.3390/biology12020274

AMA Style

Gwokyalya R, Weldon CW, Herren JK, Gichuhi J, Makhulu EE, Ndlela S, Mohamed SA. Friend or Foe: Symbiotic Bacteria in Bactrocera dorsalis–Parasitoid Associations. Biology. 2023; 12(2):274. https://doi.org/10.3390/biology12020274

Chicago/Turabian Style

Gwokyalya, Rehemah, Christopher W. Weldon, Jeremy Keith Herren, Joseph Gichuhi, Edward Edmond Makhulu, Shepard Ndlela, and Samira Abuelgasim Mohamed. 2023. "Friend or Foe: Symbiotic Bacteria in Bactrocera dorsalis–Parasitoid Associations" Biology 12, no. 2: 274. https://doi.org/10.3390/biology12020274

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