www.nature.com/scientificreports OPEN Progeny fitness determines the performance of the parasitoid Therophilus javanus, a prospective biocontrol agent against the legume pod borer Djibril Aboubakar Souna1, Aimé Hippolyte Bokonon‑Ganta2, Marc Ravallec3, Mesmin Alizannon1, Ramasamy Srinivasan4, Barry Robert Pittendrigh5, Anne‑Nathalie Volkoff3 & Manuele Tamò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)1. This parasitoid was recently imported from Taiwan into Benin as a candidate biological control agent against the cowpea pod borer in West A frica2. 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 a gent3. Already, the biological potential of this newly introduced natural enemy has been investigated with regard to reproductive b iology4 and the ability of this natural enemy to localize a plant-feeding M. vitrata caterpillar5. However, life traits estimation could provide additional and critical information on quantitative and ecological function of natural enemies6, and population growth can be estimated by assessing fertility life table parameters7. These parameters, when known for both a pest and its natural enemy, may also help to plan the outcome of their i nteraction8. Although T. javanus has been known in Asia for a long time2, 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)9, parasitoid feeding strategies during larval stage10,11 and adult feeding behavior12. 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 1International Institute of Tropical Agriculture, Benin Research Station (IITA-Benin), 08 BP 0932 Tri Postal, Cotonou, Benin. 2Department of Crop Production, Faculty of Agronomic Sciences (FSA), University of Abomey-Calavi (UAC), 03 BP 2819 Cotonou, Benin. 3UMR DGIMI University of Montpellier, INRAE, Place Eugène Bataillon, 34 095 Montpellier, France. 4World Vegetable Center (WorldVeg), Shanhua, Tainan 74151, Taiwan. 5Department of Entomology, Michigan State University (MSU), East Lansing, USA. *email: m.tamo@cgiar.org Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 1 Vol.:(0123456789) www.nature.com/scientificreports/ Figure 1. Therophilus javanus (Bhat & Gupta), female, Goergen Georg/IITA Benin. Scale bar = 5 mm. development. Thus, some parasitoids should maximize food consumption during larval stage and ultimately favor adult fi tness13. Along this line, two different larval feeding strategies have been reported in koinobiont e ndoparasitoids14. 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”)10. 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”)15. 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 e gression16. However, in tissue feeder parasitoids, feeding after egression can also enhance adult fitness10. 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 f ecundity17,18. Host quality at the time of oviposition and parasitoid larval feeding strategies are known to significantly influence adult fitness19,20. 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 off- spring 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 ( see21 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: χ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 influ- enced by the host larval instar chosen by the female at the time of oviposition (χ2 = 7.0343, df = 2, p = 0.03). More males emerged from M. vitrata caterpillars parasitized at the first larval instar (χ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 (χ2 = 0.80645, df = 1, p = 0.4) and (χ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. Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 2 Vol:.(1234567890) www.nature.com/scientificreports/ Host larval instar Oviposition Sex ratio First 0.40 ± 0.10a 2.4a Second 0.40 ± 0.14a 0.72b Third 0.20 ± 0.08b 0.69b Fourth 0 – Fifth 0 – 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 2. Larva of T. javanus, (arrow indicates the head): (a) first instar, (b) second instar, (c) third instar. Scale bar = 5 mm. 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. 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 hyme- nopteriform 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). Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 3 Vol.:(0123456789) www.nature.com/scientificreports/ 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. 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. 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: χ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: χ2 = 0.81308, df = 2, p = 0.07). Adult females lived 5.46 ± 2.81 days (n = 30), 5.63 ± 1.54 days (n = 30) and 4.53 ± 1.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 ± 0.32, 0.56 ± 0.12, and 0.3 ± 8.38 days, respectively. The host age at the time of oviposition significantly affected the mean number of offspring produced per female per day (GLM: χ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). Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 4 Vol:.(1234567890) www.nature.com/scientificreports/ Figure 6. Survival rates for adult females of T. javanus exposed to two-, three-, and four-day-old M. vitrata caterpillars. 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). Figure 8. Dynamics of offspring produced daily per adult T. javanus female exposed to two-, three- and four- day-old M. vitrata caterpillars. 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: χ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 Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 5 Vol.:(0123456789) www.nature.com/scientificreports/ Parameters Host age r (day-1) λ(day-1) Ro (offspring) DT (day) Two-day-old 0.24 ± 0.01a 1.28 ± 0.01a 66.63 ± 6.62b 2.79 ± 0.07b Three-day-old 0.26 ± 0.01a 1.29 ± 0.01a 92.03 ± 9.87a 2.66 ± 0.06b Four-day-old 0.21 ± 0.01b 1.23 ± 0.01b 35.93 ± 6.51c 3.26 ± 0.16a 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. 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 (χ2 = 261.09, df = 2, p < 0.001), and was male-biased on two- (χ2 = 197.92, df = 1, p < 0.001) and three- (χ2 = 739.6, df = 1, p < 0.001) but not four-day-old hosts (χ2 = 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) (χ2 = 58.216, df = 1, p < 0.001) and four-day-old hosts (0.51) (χ2 = 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 (χ2 = 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, sug- gesting 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)22,23. 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” p arasitoids14, and is utilized by some parasitoids to exploit a wide range of host stages14. However, after egression, some of them need to feed externally on the host tissue to acquire additional nutrients for successfully completing cocoon and adult s tages11. 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 egression11. Subsequent studies10 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 days24, 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 discrimi- nation strategies for selecting suitable host stages that can increase their population fitness (c.f., ’mother knows best’25). 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 parasitism26. Previous work has highlighted the presence of a venom gland attached to the reproductive tract of T. javanus females4 and hypothesized that the venom fluid might be secreted in the host during parasitism27. Although the venom function has not been Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 6 Vol:.(1234567890) www.nature.com/scientificreports/ 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 l arva28. 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)24. 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 Ro29. It was shown that parasitizing early host stages can reduce the development time in koinobionts p arasitoids16; however, adult fitness of the offspring may be affected, in terms of short life expectancy, small size eggs, and reduced oviposition p eriod30. Moreover, it has already been assessed that the development in early M. vitrata larval instars reduced the offspring potential fecundity in T. javanus4. 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 ovipo- sition 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)31. Environmental and genetic factors, as well as maternal decision can regulate sex allocation in Braconid p arasitoids32. The genetic basis has been studied for species such as Cotesia glomerata, Asobara tabida and Alysia manducator (Hymenoptera: Braconidae)33,34. 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)35. 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 labora- tory, 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 Insti- tute 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 i n4. Briefly, M. vitrata caterpillars were reared in large cylindrical plastic boxes (11 cm height × 16.5 cm diameter) contain- ing 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 labo- ratory 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 caterpil- lars (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 preview 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 popula- tion 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 Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 7 Vol.:(0123456789) www.nature.com/scientificreports/ 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 imma- ture 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 develop- ment 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 i n36: 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 d ata37. 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 software38. 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 s oftware39. All these tests were conducted with R software package R 3.5.340. The life table parameters were calculated following standard methodology41,42 using the TWOSEX-MSChart software43. 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%36. 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,00044. 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 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 leaffolder 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 agath- idinae), 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/1 0. 1016/ b978- 012257 305-7/ 50055-2. 7. Maia, A. D. H., Luiz, A. J. & Campanhola, C. Statistical inference on associated fertility life table parameters using jackknife tech- nique: 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. Postegression feeding enhances growth, survival, and nutrient acquisition in the endoparasitoid Toxoneuron nigriceps (Hymenoptera: Braconidae). J. Insect Sci. 15, 51 (2015). Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 8 Vol:.(1234567890) www.nature.com/scientificreports/ 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://d oi. org/ 10. 3897/ zooke ys.5 4. 475 (2010). 22. Kuriachan, I., Consoli, F. L. & Vinson, S. B. In vitro rearing of Toxoneuron nigriceps (Hymenoptera: Braconidae), a larval endopara- sitoid 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, braconi- dae) 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 per- formance 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://d oi. 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://d oi. org/ 10. 1038/ sj. hdy. 68008 29 (2006). 35. Van Nieuwenhove, G. A. & Ovruski, S. M. Influence of Anastrepha fraterculus (Diptera: Tephritidae) larval instars on the production of Diachasmimorpha longicaudata (Hymneoptera: 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-1 4. 39. Mangiafico, S. S. rcompanion: Functions to Support Extension Education Program Evaluation. R package version 2.0.0. Cran R (2016). https://d oi. org/ 10.1 016/J.A MJSUR G.2 007. 06.0 26. 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 gifuensis (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.1 97. 173/ Ecolo gy/ (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. Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 9 Vol.:(0123456789) www.nature.com/scientificreports/ 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. 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:// creat ivecom mons. org/ licens es/b y/4. 0/. © The Author(s) 2021 Scientific Reports | (2021) 11:8990 | https://doi.org/10.1038/s41598-021-88644-3 10 Vol:.(1234567890)