# Progeny fitness determines the performance of the parasitoid Therophilus javanus, a prospective biocontrol agent against the legume pod borer

Djibril Aboubakar Souna, Aimé Hippolyte Bokonon-Ganta, Marc Ravallec, Mesmin Alizannon, Ramasamy Srinivasan, Barry Robert Pittendrigh, Anne-Nathalie Volkoff, Manuele Tamò

► To cite this version:

Djibril Aboubakar Souna, Aimé Hippolyte Bokonon-Ganta, Marc Ravallec, Mesmin Alizannon, Ramasamy Srinivasan, et al. Progeny fitness determines the performance of the parasitoid Therophilus javanus, a prospective biocontrol agent against the legume pod borer. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-88644-3. http://hal-03259271

HAL Id: hal-03259271

https://hal.inrae.fr/hal-03259271

Submitted on 14 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Progeny fitness determines the performance of the parasitoid *Therophilus javanus*, a prospective biocontrol agent against the legume pod borer

Djibril Aboubakar Souna<fnref n="1" />, Aimé Hippolyte Bokonon-Ganta², Marc Ravallec³, Mesmin Alizannon<fnref n="1" />, Ramasamy Srinivasan⁴, Barry Robert Pittendrigh⁵, Anne-Nathalie Volkoff³ & Manuele Tamò<fnref n="1" />✉

*Therophilus javanus* (Bhat & Gupta) is an exotic larval endoparasitoid newly imported from Asia into Africa as a classical biological control agent against the pod borer *Maruca vitrata* (Fabricius). The parasitoid preference for the five larval instars of *M. vitrata* and their influence on progeny sex ratio were assessed together with the impact of larval host age at the time of oviposition on development time, mother longevity and offspring production. In a choice situation, female parasitoids preferred to oviposit in the first three larval instars. The development of immature stages of the parasitoid was observed inside three-day-old hosts, whereby the first two larval instars of *T. javanus* completed their development as endoparasites and the third larval instar as ectoparasite. The development time was faster when first larval instars (two- and three-day-old) of the host caterpillars were parasitized compared to second larval instar (four-day-old). The highest proportion of daughters (0.51) was observed when females were provided with four-day-old hosts. The lowest intrinsic rate of increase (r) (0.21 ± 0.01), the lowest rate of increase (λ) (1.23 ± 0.01), and the lowest net reproductive rate (Ro) (35.93 ± 6.51) were recorded on four-day-old hosts. These results are discussed in the light of optimizing mass rearing and release strategies.

*Therophilus javanus* (Bhat & Gupta) (Hymenoptera: Braconidae) is a solitary endoparasitoid attacking larval stages of the cowpea pod borer *Maruca vitrata* (Fabricius) (Lepidoptera: Crambidae)<fnref n="1" />. This parasitoid was recently imported from Taiwan into Benin as a candidate biological control agent against the cowpea pod borer in West Africa². Knowledge about life history traits—i.e., the main biological parameters that can affect the reproductive function and survivorship of an organism—is fundamental for understanding how organisms adapt to their environment. In the case of a parasitoid, this knowledge is a prerequisite prior to its introduction in a novel environment as a biological control agent³. Already, the biological potential of this newly introduced natural enemy has been investigated with regard to reproductive biology⁴ and the ability of this natural enemy to localize a plant-feeding *M. vitrata* caterpillar³. However, life traits estimation could provide additional and critical information on quantitative and ecological function of natural enemies⁶, and population growth can be estimated by assessing fertility life table parameters⁷. These parameters, when known for both a pest and its natural enemy, may also help to plan the outcome of their interaction⁸. Although *T. javanus* has been known in Asia for a long time², there has been no attempt to date to the authors' knowledge to study *T. javanus* life history traits in general.

Different factors can shape parasitoid life history trait: host quality (size, age, or stage)⁹, parasitoid feeding strategies during larval stage¹⁰,¹¹ and adult feeding behavior¹². Because *T. javanus* parasitizes a single host, its suitability and the amount of nutritious resources it provides, will be a determinant for parasitoid survival and

<footnote n="1">¹International Institute of Tropical Agriculture, Benin Research Station (IITA-Benin), 08 BP 0932 Tri Postal, Cotonou, Benin. ²Department of Crop Production, Faculty of Agronomic Sciences (FSA), University of Abomey-Calavi (UAC), 03 BP 2819 Cotonou, Benin. ³UMR DGIMI University of Montpellier, INRAE, Place Eugène Bataillon, 34 095 Montpellier, France. ⁴World Vegetable Center (WorldVeg), Shanhua, Tainan 74151, Taiwan. ⁵Department of Entomology, Michigan State University (MSU), East Lansing, USA. ✉email: m.tamo@cgiar.org</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** *Therophilus javanus* (Bhat & Gupta), female, Goergen Georg/IITA Benin. Scale bar = 5 mm.</figcaption></figure>

development. Thus, some parasitoids should maximize food consumption during larval stage and ultimately favor adult fitness¹³. Along this line, two different larval feeding strategies have been reported in koinobiont endoparasitoids¹⁴. Some complete part of their larval stages inside the host body, but after egression from the host at a late larval stage destructively consume all (or most) of the host tissues before formation of cocoon ("tissue feeders")¹⁰. Others complete their larval stages inside the host by feeding nondestructively on the host hemolymph or fat body and then egress from the host to form a cocoon on, or close, to the host, which may remain active and live temporally before dying ("hemolymph feeders")¹⁵. Overall, hemolymph-feeder parasitoids seem to be more advantaged because they can exploit a wide range of host stages and can use the dying host as a bodyguard against predators and hyperparasitoids after egression¹⁶. However, in tissue feeder parasitoids, feeding after egression can also enhance adult fitness¹⁰. Consequently, for both tissue feeders and hemolymph feeders, the quality of the host parasitized by the mother parasitoid can affect progeny fitness, including immature survival, development, adult longevity, and fecundity¹⁷,¹⁸.

Host quality at the time of oviposition and parasitoid larval feeding strategies are known to significantly influence adult fitness¹⁹,²⁰. The objectives of the present study were to determine the larval feeding behavior in *T. javanus*, investigate the impact of the host age on total development time, mother female longevity and offspring production, sex-ratio of the progeny, and estimate the life table parameters under laboratory conditions.

## Results

**Brief description of Therophilus javanus.** *Therophilus javanus* is a braconid parasitoid of Agathidinae subfamily with a long ovipositor sheath (Fig. 1). The mesonotum of the parasitoid is dark orange, the metasoma is black, while the first to third metasomal sternites are whitish (see²¹ for more details).

**Host stage preference and progeny sex ratio.** Host stage preference and the sex ratio—i.e. the ratio of males to females in the progeny — of the parasitoid were significantly affected by the larval instar of *M. vitrata* caterpillar (Table 1). Parasitization attacks and subsequent oviposition occurred more frequently on first and second larval stages, and less often on third larval instars (GLM: $\chi^2 = 132.3$, df = 4, *p* < 0.001). The female did not oviposit on the fourth and the fifth larval instars during the observation. The sex ratio was significantly influenced by the host larval instar chosen by the female at the time of oviposition ($\chi^2 = 7.0343$, df = 2, *p* = 0.03). More males emerged from *M. vitrata* caterpillars parasitized at the first larval instar ($\chi^2 = 5.7647$, df = 1, *p* = 0.02). However, there were no significant difference between the number of males and females when the host was parasitized at the second or the third larval instar ($\chi^2 = 0.80645$, df = 1, *p* = 0.4) and ($\chi^2 = 0.72727$, df = 1, *p* = 0.4), respectively.

**Larval feeding behavior and developmental time of T. javanus immature stages.** At 26 ± 1.1 °C, *T. javanus* egg development lasted an average 1.94 ± 0.16 days (n = 45). Larval development included three larval instars. The duration of the first larval instar was 3.04 ± 0.21 days (n = 137). The second and third larval instars lasted 2.04 ± 0.21 days (n = 92) and 1.84 ± 0.37 days (n = 83), respectively.
<table><thead><tr><td>Host larval instar</td><td>Oviposition</td><td>Sex ratio</td></tr></thead><tbody><tr><td>First</td><td>0.40 ± 0.10a</td><td>2.4a</td></tr><tr><td>Second</td><td>0.40 ± 0.14a</td><td>0.72b</td></tr><tr><td>Third</td><td>0.20 ± 0.08b</td><td>0.69b</td></tr><tr><td>Fourth</td><td>0</td><td>-</td></tr><tr><td>Fifth</td><td>0</td><td>-</td></tr></tbody></table>

**Table 1.** Percentage of successful oviposition (mean ± SE) and sex ratio of *T. javanus* in choice situation among the five instars of *M. vitrata* caterpillars. Means ± SE of oviposition followed by different letters in the same column are significantly different among treatments (GLM and Tukey's HSD test, *p* < 0.05). Sex ratio followed by different letters in the same column are significantly different among treatments (chi-square test followed by 'pairwise nominal independence').

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Larva of *T. javanus*, (arrow indicates the head): (a) first instar, (b) second instar, (c) third instar. Scale bar = 5 mm.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Third instar *T. javanus* larva egressed from *M. vitrata* caterpillar. (a) *T. javanus* larva, white colored, feeding on *M. vitrata* caterpillar within the *M. vitrata* cocoon indicated by an arrow; (b) arrow indicating the egression hole on *M. vitrata* caterpillar (picture taken after host cocoon has been removed). Scale bar = 5 mm.</figcaption></figure>

**Endoparasitic feeding.** First instar larvae were translucent, polypodeiform, with thirteen body segments and a distinct sclerotized head capsule harboring a pair of prominent mandibles (Fig. 2a). Each body segment except the last one presented two pairs of ventral processes of uniform size. The last segment was prolonged by a small, ventrally curved tail. The gut was visible, and content changed in color from translucent in early first instar to white or green in the late first instar of *T. javanus* reared on *M. vitrata*. The second instar was of hymenopteriform type (Fig. 2b). The head was hemispherical with no apparent mandibles. At this stage, there were no more visible paired ventral processes and the ventrally curved tail was lost. The larva was opaque white, and the gut was green. The early third instar was also of hymenopteriform type. The body was opaque white, but the gut was green and contained white granules (Fig. 2c).

**Ectoparasitic feeding.** About eight days after parasitism, the third instar larva egressed from the host to become ectoparasitic (Fig. 3a). The body color remained opaque white and the gut turned from green to yellow during this last part of the larval development. At this stage, the larval body was cylindrical, narrowed at both ends, and the anterior end was curved. Egressed larvae fed on the host larva by sucking all the host contents through the egression hole (Fig. 3b).

**Cocoon and pupa formation.** After the egressed larva emptied the host contents, the late third instar spun a cocoon (Fig. 4) before entering in pre-pupal stage. The pre-pupa instar displayed a dirty-white color and was identified when the head, thorax, and abdomen of the future adult started to be distinguishable (Fig. 5a); its duration was 1.11 ± 0.32 days (n = 50). The color of the newly formed pupa turned into ivory-white within a cocoon after expelling the brown meconium contained in the gut (Fig. 5b). During pupation, pigmentation was progressive and started from the head region, then reached the thorax, the abdomen, the legs, the antenna and, finally, the ovipositor (for the females) (Fig. 5c, d). Total duration of pupal stage was 5.65 ± 0.82 days (n = 41).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Morphology of *T. javanus* cocoon. (a) *T. javanus* spins the cocoon within the silken cocoon of *M. vitrata*; (b) view of the newly formed *T. javanus* cocoon (right) dissected from the cocoon of *M. vitrata* (left); (c) view of the fully formed cocoon of *T. javanus*. Scale bar = 5 mm.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Progressive body changes and pigmentation during pupation of female *T. javanus*. (a) pre-pupa, with arrow indicating the head; (b) hyaline pupa, with arrow indicating the ovipositor; c-d) progressive body pigmentation during pupation. Scale bar = 5 mm.</figcaption></figure>

**Impact of host age on adult parasitoid development time.** The total development time (from egg to adult) was found to be significantly influenced by the age of the host at the time of oviposition (GLM: $\chi^2=0.32489$, df=2, *p*<0.001). *Therophilus javanus* development time was faster when two- and three-day-old (first larval instars) *M. vitrata* caterpillars were parasitized [15.35±0.78 days (n=85) and 15.43±0.91 days (n=74), respectively] compared to four-day-old (second larval instar) *M. vitrata* caterpillars (16.73±0.82 days; n=64).

**Impact of host age on mother parasitoid longevity.** The host age at the time of oviposition did not significantly impact the mean longevity of adult females (GLM: $\chi^2=0.81308$, df=2, *p*=0.07). Adult females lived $5.46\pm2.81$ days (n=30), $5.63\pm1.54$ days (n=30) and $4.53\pm1.569$ days (n=30) when exposed to hosts that were two-, three-, and four-days-old, respectively. Survival of parasitoid females exposed to four-day-old hosts declined sharply after four days and reached zero after nine days, with that of three-day-old host displaying a similar pattern shifted by one day. While parasitoid females exposed to two-day-old hosts had a similar initial survival to those exposed to three-day-old host, their survival rate was the longest, reaching zero only after fourteen days (Fig. 6).

**Host age impact on offspring production.** The pre-oviposition period was similar whether females were exposed to two-, three-, or four-day-old hosts (*p*>0.05) with $0.46\pm0.32$, $0.56\pm0.12$, and $0.3\pm8.38$ days, respectively. The host age at the time of oviposition significantly affected the mean number of offspring produced per female per day (GLM: $\chi^2=350.27$, df=2, *p*<0.001), with a higher mean number of offspring produced on two- and three-day-old hosts (Fig. 7).
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Survival rates for adult females of *T. javanus* exposed to two-, three-, and four-day-old *M. vitrata* caterpillars.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Mean offspring produced per female *T. javanus* exposed to two-, three-, and four-day-old *M. vitrata* caterpillars. Bars with the same letters indicate not significantly different (GLM and Tukey's HSD test, *p* < 0.05).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Dynamics of offspring produced daily per adult *T. javanus* female exposed to two-, three- and four-day-old *M. vitrata* caterpillars.</figcaption></figure>

Females laid eggs until 12, 8, and 7 days post emergence, respectively, on two-, three-, or four-day-old *M. vitrata* caterpillars (Fig. 8). Mean number of offspring produced daily per female varied significantly (GLM: $\chi^2=879.8$, df=26, *p* < 0.001), with the curves corresponding to offspring production per female showing that within the first seven days after emergence the mean number of offspring produced daily in four-day-old hosts was lower than on two- and three-day-old hosts (Fig. 8). Maximum average numbers of offspring produced daily
<table><caption>**Table 2.** Life table parameters of *T. javanus* reared on two-, three- and four-day-old *M. vitrata* caterpillars (mean ± SE). Means ± SE followed by different letters in the same column are significantly different among treatments using the paired bootstrap test at 5% significance level.</caption><thead><tr><th rowspan="2">Host age</th><th colspan="4">Parameters</th></tr><tr><th>r (day<sup>-1</sup>)</th><th>λ (day<sup>-1</sup>)</th><th>Ro (offspring)</th><th>DT (day)</th></tr></thead><tbody><tr><td>Two-day-old</td><td>0.24 ± 0.01a</td><td>1.28 ± 0.01a</td><td>66.63 ± 6.62b</td><td>2.79 ± 0.07b</td></tr><tr><td>Three-day-old</td><td>0.26 ± 0.01a</td><td>1.29 ± 0.01a</td><td>92.03 ± 9.87a</td><td>2.66 ± 0.06b</td></tr><tr><td>Four-day-old</td><td>0.21 ± 0.01b</td><td>1.23 ± 0.01b</td><td>35.93 ± 6.51c</td><td>3.26 ± 0.16a</td></tr></tbody></table>

per adult *T. javanus* female were 16.57 ± 2.09, 18.19 ± 2.36, and 11.2 ± 1.50, observed at day 2, 5, and 2 on two-, three- and four-day-old hosts, respectively (Fig. 8).

**Impact of host age on offspring sex-ratio.** The sex-ratio was affected by the host age at the time of oviposition (χ²=261.09, df=2, *p*<0.001), and was male-biased on two- (χ²=197.92, df=1, *p*<0.001) and three- (χ²=739.6, df=1, *p*<0.001) but not four-day-old hosts (χ²=0.53432, df=1, *p*>0.05). Females exposed to three-day-old hosts produced the lowest proportion of females (0.24) compared to females exposed to two- (0.34) (χ²=58.216, df=1, *p*<0.001) and four-day-old hosts (0.51) (χ²=259.63, df=1, *p*<0.001). Similarly, females exposed to two-day-old hosts produced a lower proportion of females compared to females exposed to four-day-old hosts (χ²=81.996, df=1, *p*<0.001).

**Impact of host age on the life table parameters.** The life table parameters of the parasitoid were significantly affected by the host age at the time of oviposition (Table 2). Higher average values of the intrinsic rate of increase (*r*) were recorded on two- (*p*<0.01) and three-day-old hosts (*p*<0.001) compared to four-day-old hosts. Similarly, higher average values of the rate of increase (*λ*), the net reproductive rate (*Ro*), and shorter doubling time (*TD*) were recorded on two- and three-day-old hosts compared to four-day-old hosts (*p*<0.001).

Discussion

This study investigated the host stage preference of *T. javanus* adult females, it characterized the endoparasitic and ectoparasitic feeding behavior of the parasitoid larva, and evaluated the effect of host age on life history parameters of the parasitoid. *Therophilus javanus* preferred to oviposit in the first three larval instars of the host. Host larval instar and host age at the time of oviposition influenced offspring production and sex allocation. The higher intrinsic rate of increase (*r*) and the net reproductive rate (*Ro*) were recorded on three-day-old hosts, suggesting that the parasitoid population would increase more rapidly when provided three-day-old hosts. Inversely, the same three-day-old hosts produced the lowest proportion of daughter offspring.

In a choice situation, *T. javanus* preferred to oviposit in the first three larval instars of *M. vitrata*. In most cases, the two older instars resisted physically to the parasitoid and started to spin quickly a transparent silken web as a barrier to parasitoid attack (D. Aboubakar Souna, personal observation).

Life history observations of *T. javanus* have revealed three distinct larval stages, similar to the number of larval instars reported in other Braconidae wasps such as *Cardiochiles nigriceps* Viereck (Hymenoptera: Braconidae) and *Toxoneuron nigriceps* (Viereck) (Hymenoptera: Braconidae), larval endoparasitoids of *Heliothis virescens* (Fabricius) (Lepidoptera: Noctuidae)²²,²³. The first two larval stages of *T. javanus* feed inside the parasitized host larva, which continues to nourish itself and grow, while the third larval instar egresses from the host and continues host-feeding on it through the hole from which it egressed. Because parasitized *M. vitrata* caterpillars continue to feed and grow, we hypothesize that during the first two larval stages, *T. javanus* may have adapted to reduce the damage induced to vital organs of its host, such as feeding on the hemolymph and/or fat body. Unexpectedly, the third larval instar of *T. javanus* larva egressed destructively from the host body and continued to feed on it. This feeding behavior is similar to the one displayed by “tissue feeder” parasitoids¹⁴, and is utilized by some parasitoids to exploit a wide range of host stages¹⁴. However, after egression, some of them need to feed externally on the host tissue to acquire additional nutrients for successfully completing cocoon and adult stages¹¹. As an example, *T. nigriceps* feeds on *H. virescens* as a “tissue feeder” to regulate the host physiology, successfully completing its larval development after host egression¹¹. Subsequent studies¹⁰ have confirmed that *T. nigriceps* exploits the dual larval feeding behavior to increase food acquisition and enhance its adult fitness. *T. javanus* may also use a similar adaptation strategy to regulate the host physiology for successfully completing its development. While under current rearing conditions the total larval duration of *M. vitrata* takes on average ten days²⁴, the total development time of the larva until parasitoid egression was up to fourteen days, depending on the age of the larva at parasitization. This suggests that the developing parasitoids larva can actively manipulate the host metabolism to slow down its development and allow successful completion until egression.

To maximize the chance for oviposition in a suitable host, several parasitoids have developed host discrimination strategies for selecting suitable host stages that can increase their population fitness (c.f., ‘mother knows best’²⁵). On the other hand, some endoparasitoids can ensure the development of their progeny through maternal factors such as the venom fluid injected in the host at the time of parasitism²⁶. Previous work has highlighted the presence of a venom gland attached to the reproductive tract of *T. javanus* females⁴ and hypothesized that the venom fluid might be secreted in the host during parasitism²⁷. Although the venom function has not been
elucidated in more detail to date, we speculate it could regulate metabolic processes in parasitized *M. vitrata* caterpillars in order to affect its development and hence ensure complete development of the parasitoid larva²⁸.

The intrinsic rate of increase (r) obtained from the estimation of the life-table parameters indicate that the population of *T. javanus* would increase faster if reared on two and three-day-old *M. vitrata* caterpillars. However, three-day-old *M. vitrata* caterpillars induced the highest net reproductive rate (*Ro*), suggesting that oviposition by *T. javanus* adult females and subsequent feeding by developing immature stages on this host stage may engender a higher production of offspring and hence contribute positively to the overall population growth. Our results support previous observations reporting a higher suitability of two and three-day-old *M. vitrata* caterpillar to *Apanteles taragamae* (Viereck) (Hymenoptera: Braconidae)²⁴. In the same vein, the second instar of *Spodoptera exigua* (Hübner) (Lepidoptera: Noctuidae) was identified as the most suitable host stage for *Microplitis similis* (Lyle) (Hymenoptera: Braconidae), based on the high values of *r* and *Ro*²⁹. It was shown that parasitizing early host stages can reduce the development time in koinobionts parasitoids¹⁶; however, adult fitness of the offspring may be affected, in terms of short life expectancy, small size eggs, and reduced oviposition period³⁰. Moreover, it has already been assessed that the development in early *M. vitrata* larval instars reduced the offspring potential fecundity in *T. javanus*⁴. Our observations further show that the sex-ratio of the offspring population was male-biased when stemming from three-day-old hosts, while the parasitoid produced more females when the oviposition occurred in four-day-old *M. vitrata* hosts. Similar occurrences have been reported in other parasitoids, e.g., the sex-ratio of *Diadegma mollipla* (Holmgren) (Hymenoptera: Ichneumonidae) was male-biased when the parasitoid developed in the early instars (L1) compared to the older instars (L2, L3 and L4) for the Diamondback Moth *Plutella xylostella* (Linnaeus) (Lepidoptera: Plutellidae)³¹. Environmental and genetic factors, as well as maternal decision can regulate sex allocation in Braconid parasitoids³². The genetic basis has been studied for species such as *Cotesia glomerata*, *Asobara tabida* and *Alysia manducator* (Hymenoptera: Braconidae)³³,³⁴. A more direct involvement of female parasitoids in sex allocation was proposed for *Diachasmimorpha longicaudata* (Ashmead) (Hymenoptera: Braconidae), observed to exploit physical stimuli and, perhaps more importantly, volatiles emitted by larvae of *Anastrepha fraterculus* (Wiedemann) (Diptera: Tephritidae)³⁵. Similarly, during the oviposition process, female *T. javanus* may be able to use the same approach for discriminating host stages providing the best resources, maximizing chances for optimal progeny development and production of daughter offspring. However, the precise mechanisms for sex allocation in *T. javanus* will need to be elucidated by further investigations.

## Conclusions

Overall, the results provide key insights into the biology of *T. javanus* that can be useful to guide field releases of this parasitoid to control *M. vitrata* in West Africa. For optimizing the mass rearing of *T. javanus* in the laboratory, using three-day-old *M. vitrata* caterpillars (corresponding to first instar larvae) could accelerate the rate of colony production because of its higher net reproductive rate (*Ro*). However, knowing that female offspring was favored in older instars of *M. vitrata* caterpillar hosts, we suggest that a mix of first, second, and third instar larvae would be the best option for maximizing the total output of the mass rearing. Additional investigations are still required to assess the foraging behavior of the parasitoid in order to determine how females will detect and discriminate the most suitable stage or stages of *M. vitrata* caterpillars in the field.

## Material and methods

**Insect rearing.** Maruca vitrata caterpillars were collected from the existing colony at International Institute of Tropical Agriculture (IITA) Benin laboratories, while the *T. javanus* colony was issued from individuals received from The World Vegetable Center (WorldVeg), Taiwan, Republic of China. Insect colonies were reared under laboratory conditions at the IITA-Benin research station following the methods described in⁴. Briefly, *M. vitrata* caterpillars were reared in large cylindrical plastic boxes (11 cm height × 16.5 cm diameter) containing sprouting cowpea grains. Colonies of *T. javanus* were maintained on *M. vitrata* first instar (three-day-old) caterpillars under confined conditions. For parasitism, *M. vitrata* caterpillars were submitted to three-day-old *T. javanus* adult mated females. Parasitized caterpillars were then reared until pupal stage on sprouting cowpea grains. Following emergence, *T. javanus* adults were fed on honey solution. Insects were maintained under laboratory conditions at 26 ± 1.1 °C, 76 ± 7% with 12:12 L:D photoperiod.

**Influence of host stage on mother preference and progeny sex ratio.** Five *M. vitrata* caterpillars (one each for larval instar 1 to 5) were exposed in a petri dish (90 mm in diameter × 15 mm in depth) to two-days-old female *T. javanus*, mated and with no previous oviposition experience. Before the experiment, the female parasitoids were provided with a drop of honey. The behavior of the parasitoid was observed during one hour and behavioral data were collected using Pocket Observer software (version 3.2; Noldus, The Netherlands). Females which did not make any oviposition attempt during the observation period were considered as non-responding and removed from the analysis. After each observation, larvae, female parasitoid, and the petri dish were replaced by the new ones. At the end of the observation period, each stung caterpillar was separated and placed individually in the small cup (3 cm diameter × 3.5 cm height), and offered fresh sprouting cowpea grains as feeding substrate until the formation of parasitoid cocoon. The collected cocoons were kept individually in the small cup (3 cm diameter × 3.5 cm height) until adult emergence and the sex ratio was noted in progeny population per host stage. A total of 72 responding females were considered for this observation.

**Larval feeding behavior and developmental time of immature stages.** At first, three-day-old *M. vitrata* caterpillars were individually submitted to parasitism by three-day-old *T. javanus* mated females. Only a single stinging was allowed per caterpillar to ensure that a single parasitoid egg has been deposited. A total
of forty-five (45) parasitized *M. vitrata* caterpillars were dissected under a macro zoom microscope (Olympus MVX10) at 24 h, 36 h and 48 h after parasitism, then daily until pre-pupae formation. *Therophilus javanus* immature stages present in the host body were removed and their image taken using an Olympus XC50 camera. The shape of larval body and mandibles were used to distinguish the three larval instars of *T. javanus*. The development of forty-one (41) pupae was observed until adult emergence.

**Impact of host age on adult parasitoid development time.** To determine the impact of the host age on *T. javanus* development time from egg to adult, a total of one hundred sixty (160) two-, three- or four-day-old *M. vitrata* caterpillars were individually submitted to parasitism. One-time stinging was allowed per caterpillar. Parasitized caterpillars were reared until cocoon formation and adult emergence as previously described. Each cocoon was separated in individual small cup (3 cm diameter × 3.5 cm height) containing a honey drop. Cups were observed daily until adult emergence.

**Impact of host age on the life table parameters.** Thirty (30) newly emerged *T. javanus* couples were transferred to a rearing unit consisting of a plastic cup (diameter: 9 cm at the base and 12 cm at the top; height: 4.5 cm) containing sprouting cowpea grains and thirty (30) two-, three- or four-day-old *M. vitrata* caterpillars, respectively. Because the mating behavior had not been studied in this parasitoid, couples were allowed to mate all the time during the experiment. Each couple was allowed to parasitize a given host age during 24 h, and this, every day until female death. If the male died, it was replaced by another male. After exposure to parasitism, *M. vitrata* caterpillars remained in the rearing unit feeding on sprouting cowpea grains until parasitoid cocoon formation. For each couple and each day of exposure, the cocoons formed were regularly collected and placed in another plastic cup containing a drop of honey. Offspring emergence was recorded daily, and the sex of each progeny was noted. Life table parameters were estimated as described in³⁶: the adult female survival rate, the finite rate of increase (λ), the intrinsic rate of increase (r), the net reproductive rate (Ro), and doubling time (DT).

**Statistical analysis.** The Shapiro-Wilk normality test was realized to test the normal distribution of data³⁷. A general linear model (GLM) with a gamma distribution was used to investigate the effect of host age at the time of oviposition on the longevity of mother parasitoids and the development time of offspring. GLM with Poisson errors and log-link functions, corrected for over-dispersion, were used to investigate the influence of host age on the mean offspring produced per day, and the mean offspring produced daily per female. General linear models (GLM) with binomial errors distribution and logit-link function corrected for over-dispersion were used to analyze the proportion of oviposition per host larval instars. Means were compared using the 'glht' function in the "multcomp" package of the R software³⁸. The impact of host larval instar and host age at the time of oviposition on offspring sex allocation was tested by a chi-square test followed by a post-hoc test using the 'pairwise nominal independence' function in the "rcompanion" package of the R software³⁹. All these tests were conducted with R software package R 3.5.3⁴⁰. The life table parameters were calculated following standard methodology⁴¹,⁴² using the TWOSEX-MSChart software⁴³. The impact of host age at the time of oviposition on the duration of the oviposition periods and each life table parameter was accessed using the paired bootstrap test at P = 5%³⁶. To generate stable estimates of standard errors of the duration of the oviposition periods and each life table parameter we used the bootstrap value of 100,000⁴⁴.

**Data availability**

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

Received: 17 April 2020; Accepted: 7 April 2021

Published online: 26 April 2021

**References**

1. Dung, D. T., Thi, L., Phuong, H. & Long, K. D. Insect parasitoid composition on soybean, some eco-biological characteristics of the parasitoid, xanthopimpla punctata fabricius on soybean leaf folder omiodes indicata (Fabricius) in Hanoi Vietnam. *ISSAAS J.* **17**, 58–69 (2011).

2. Srinivasan, R., Yule, S., Lin, M. Y. & Khumsuwan, C. Recent developments in the biological control of legume pod borer (*Maruca vitrata*) on yard-long bean. *Acta Hort.* **1102**, 143–149 (2015).

3. Van Lenteren, J. C. et al. Environmental risk assessment of exotic natural enemies used in inundative biological control. *Biocontrol* **48**, 3–38 (2003).

4. Aboubakar Souna, D. et al. An insight in the reproductive biology of *Therophilus javanus* (hymenoptera, braconidae, and agathidinae), a potential biological control agent against the legume pod borer (Lepidoptera, Crambidae). *Psyche* (London) **2017**, 1–8 (2017).

5. Aboubakar Souna, D. et al. Volatiles from *Maruca vitrata* (Lepidoptera, Crambidae) host plants influence olfactory responses of the parasitoid *Therophilus javanus* (Hymenoptera, Braconidae, Agathidinae). *Biol. Control* **130**, 104–109 (2019).

6. Bellows, T. S. & Van Driesche, R. G. Life Table Construction and Analysis for Evaluating Biological Control Agents. in *Handbook of Biological Control* 199–223 (Elsevier, 1999). https://doi.org/10.1016/b978-012257305-7/50055-2.

7. Maia, A. D. H., Luiz, A. J. & Campanhola, C. Statistical inference on associated fertility life table parameters using jackknife technique: computational aspects. *J. Econ. Entomol.* **93**(2), 511–518 (2000).

8. Roy, M., Brodeur, J. & Cloutier, C. Effect of temperature on intrinsic rates of natural increase (rm) of a coccinellid and its spider mite prey. *Biocontrol* **48**, 57–72 (2003).

9. Harvey, J. A. Factors affecting the evolution of development strategies in parasitoid wasps: The importance of functional constraints and incorporating complexity. *Entomol. Exp. Appl.* **117**, 1–13 (2005).

10. Henderson, R. E., Kuriachan, I. & Vinson, S. B. Postregression feeding enhances growth, survival, and nutrient acquisition in the endoparasitoid *Toxoneuron nigriceps* (Hymenoptera: Braconidae). *J. Insect Sci.* **15**, 51 (2015).
11. Kuriachan, I., Henderson, R., Laca, R. & Vinson, S. B. Post-egression host tissue feeding is another strategy of host regulation by the koinobiont wasp toxoneuron nigriceps. *J. Insect Sci.* **11**, 1–11 (2011).

12. Benelli, G. et al. The impact of adult diet on parasitoid reproductive performance. *J. Pest. Sci.* **90**, 807–823 (2017).

13. Harvey, J. A. & Strand, M. R. The developmental strategies of endoparasitoid wasps vary with host feeding ecology. *Ecology* **83**, 2439–2451 (2002).

14. Harvey, J. A., Bezemer, T. M., Gols, R., Nakamatsu, Y. & Tanaka, T. Comparing the physiological effects and function of larval feeding in closely-related endoparasitoids (Braconidae: Microgastrinae). *Physiol. Entomol.* **33**, 217–225 (2008).

15. Harvey, J. A. et al. Development of a solitary koinobiont hyperparasitoid in different instars of its primary and secondary hosts. *J. Insect Physiol.* **90**, 36–42 (2016).

16. Harvey, J. A. & Malcicka, M. Nutritional integration between insect hosts and koinobiont parasitoids in an evolutionary framework. *Entomol. Exp. Appl.* **159**, 181–188 (2016).

17. Gols, R., Ros, V. I. D., Ode, P. J., Vyas, D. & Harvey, J. A. Varying degree of physiological integration among host instars and their endoparasitoid affects stress-induced mortality. *Entomol. Exp. Appl.* **167**, 424–432 (2019).

18. Vieira, L. J. P., Franco, G. M. & Sampaio, M. V. Host preference and fitness of *Lysiphlebus testaceipes* (Hymenoptera: Braconidae) in different instars of the aphid Schizaphis graminum. *Neotrop. Entomol.* **48**, 391–398 (2019).

19. Harvey, J. A. Dynamic effects of parasitism by an endoparasitoid wasp on the development of two host species: Implications for host quality and parasitoid fitness. *Ecol. Entomol.* **25**, 267–278 (2000).

20. Harvey, J. A., Kadash, K. & Strand, M. R. Differences in larval feeding behavior correlate with altered developmental strategies in two parasitic wasps: Implications for the size-fitness hypothesis. *Oikos* **88**, 621–629 (2000).

21. van Achterberg, C. & Long, K. D. Revision of the Agathidinae (Hymenoptera, Braconidae) of Vietnam, with the description of forty-two new species and three new genera. *Zookeys* https://doi.org/10.3897/zookeys.54.475 (2010).

22. Kuriachan, I., Consoli, F. L. & Vinson, S. B. In vitro rearing of *Toxoneuron nigriceps* (Hymenoptera: Braconidae), a larval endoparasitoid of *Heliothis virescens* (Lepidoptera: Noctuidae) from early second instar to third instar larvae. *J. Insect Physiol.* **52**, 881–887 (2006).

23. Pennacchio, F., Vinson, S. B. & Tremblay, E. Growth and development of Cardiochiles nigriceps viereck (hymenoptera, braconidae) larvae and their synchronization with some changes of the hemolymph composition of their host, Heliothis virescens (F.) (Lepidoptera, Noctuidae). *Arch. Insect Biochem. Physiol.* **24**, 65–77 (1993).

24. Dannon, E. A., Tamò, M., van Huis, A. & Dicke, M. Functional response and life history parameters of *Apanteles taragamae*, a larval parasitoid of *Maruca vitrata*. *Biocontrol* **55**, 363–378 (2010).

25. Henry, L. M., Gillespie, D. R. & Roitberg, B. D. Does mother really know best? Oviposition preference reduces reproductive performance in the generalist parasitoid *Aphidius ervi*. *Entomol. Exp. Appl.* **116**, 167–174 (2005).

26. Moreau, S. J. M. & Asgari, S. Venom proteins from parasitoid wasps and their biological functions. *Toxins (Basel)* **7**, 2385–2412 (2015).

27. Aboubakar Souna, D. Assessing the potential of *Therophilus javanus*, a biological control candidate against the cowpea pod borer *Maruca vitrata* in West Africa. (University of Montpellier (France); University of Abomey-Calavi (Benin), 2018).

28. Beckage, N. E. & Gelman, D. B. Wasp parasitoid disruption of host development: Implications for new biologically based strategies for insect control. *Annu. Rev. Entomol.* **49**, 299–330 (2004).

29. Qiu, B., Zhou, Z. & Xu, Z. Age Preference and Fitness of *Microplitis manilae* (Hymenoptera: Braconidae) Reared on *Spodoptera exigua* (Lepidoptera: Noctuidae). *Florida Entomol.* **96**, 602–609 (2013).

30. Mayhew, P. J. Comparing parasitoid life histories. *Entomol. Exp. Appl.* **159**, 147–162 (2016).

31. Sithole, R., Chinwada, P. & Lohr, B. L. Effects of host larval stage preferences and diet on life history traits of *Diadegma mollipla*, an African parasitoid of the Diamondback Moth. *Biocontrol Sci. Technol.* **28**, 172–184 (2018).

32. Boulton, R. A., Collins, L. A. & Shuker, D. M. Beyond sex allocation: The role of mating systems in sexual selection in parasitoid wasps. *Biol. Rev.* https://doi.org/10.1111/brv.12126 (2015).

33. Beukeboom, L. W., Ellers, J. & Van Alphen, J. J. M. Absence of single-locus complementary sex determination in the braconid wasps *Asobara tabida* and *Alysia manducator*. *Heredity* (Edinb). **84**, 29–36 (2000).

34. Zhou, Y., Gu, H. & Dorn, S. Single-locus sex determination in the parasitoid wasp *Cotesia glomerata* (Hymenoptera: Braconidae). *Heredity* (Edinb). https://doi.org/10.1038/sj.hdy.6800829 (2006).

35. Van Nieuwenhove, G. A. & Ovruski, S. M. Influence of *Anastrepha fraterculus* (Diptera: Tephritidae) larval instars on the production of *Diachasmimorpha longicaudata* (Hymenoptera: Braconidae) progeny and their sex ratio. *Florida Entomol.* **94**, 863–868 (2011).

36. Huang, Y. B. & Chi, H. Life tables of *Bactrocera cucurbitae* (Diptera: Tephritidae): With an invalidation of the jackknife technique. *J. Appl. Entomol.* **137**, 327–339 (2013).

37. Shapiro, A. S. S. & Wilk, M. B. Biometrika trust an analysis of variance test for normality (Complete Samples) Published by: Oxford University Press on behalf of Biometrika Trust Stable. *Biometrika* **52**, 591–611 (1965).

38. Bretz, F., Hothorn, T. & Westfall, P. H. Multiple comparisons using R. 187 (2011). https://doi.org/10.1128/AAC.03728-14.

39. Mangiafico, S. S. rcompanion: Functions to Support Extension Education Program Evaluation. R package version 2.0.0. *Cran R* (2016). https://doi.org/10.1016/j.AMJSURG.2007.06.026.

40. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria (2019).

41. Chi, H. & Su, H.-Y. Age-stage, two-sex life tables of *Aphidius gefiuisis* (Ashmead) (Hymenoptera: Braconidae) and its host *Myzus persicae* (Sulzer) (Homoptera: Aphididae) with mathematical proof of the relationship between female fecundity and the net reproductive rate. *Environ. Entomol.* **35**, 10–21 (2006).

42. Ning, S., Zhang, W., Sun, Y. & Feng, J. Development of insect life tables: Comparison of two demographic methods of *Delia antiqua* (Diptera: Anthomyiidae) on different hosts. *Sci. Rep.* **7**, 4821 (2017).

43. Chi, H. TWOSEX-MSChart: a computer program for the age stage, two-sex life table analysis. *National Chung Hsing University, Taichung, Taiwan. 2015b* Available: http://140.120.197.173/Ecology/ (access (2015)).

44. Akca, I., Ayvaz, T., Yazici, E., Smith, C. L. & Chi, H. Demography and population projection of *Aphis fabae* (Hemiptera: Aphididae): with additional comments on life table research criteria. *J. Econ. Entomol.* **108**, 1466–1478 (2015).

## Acknowledgements

We thank the Cooperation and Cultural Action Service (SCAC) of the French Embassy in Cotonou for partial support of DAS (N 898392E). Co-funding was also provided by the CGIAR Research Program on Grain Legumes and Dryland Cereals (GLDC) (DAS, MT and BRP), and the Bill and Melinda Gates Foundation (BMGF) for partial co-funding of the study (OPP1082463 to MT and BRP) and Open Access publication support. We also thank Georg Goergen for the habitus illustration of *T. javanus*, and Basile Dato for technical assistance with cowpea field installation at the International Institute of Tropical Agriculture (IITA), Benin Station.
**Author contributions**

Conceptualization and methodology: D.A.S., A.H. B.-G., M.R., R.S., A.-N.V., M.T. Data collection: D.A.S., M.A.

Data analysis: D.A.S. Writing-original draft preparation: D.A.S. Reviewing and editing: M.T., A.-N.V., A.H.B.-G., R.S., B.R.P. Supervision: A.H.B.-G., A.-N.V., M.T. Resource mobilization: D.A.S., B.R.P., A.-N.V., R.S., M.T.

**Competing interests**

The authors declare no competing interests.

**Additional information**

Correspondence and requests for materials should be addressed to M.T.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021