Risk assessment and spread of the potentially invasive Ceratitis rosa Karsch and Ceratitis quilicii De Meyer, Mwatawala & Virgilio sp. Nov. using life-cycle simulation models: Implications for phytosanitary measures and management

Distribution and host ranges of Ceratitis rosa and Ceratitis quilicii (Diptera: Tephritidae) in South Africa

Two fruit fly (Diptera: Tephritidae) species of economic importance: Ceratitis rosa Karsch and Ceratitis quilicii De Meyer, Mwatawala & Virgilio are present in South Africa. The two species were considered as one species prior to 2016, but were subsequently separated. In this study, the distribution and abundance of the two species were quantified in seven provinces in South Africa through trapping with Enriched Ginger Oil as an attractant. Trapping was conducted over three seasons across two years (2020 and 2021): late summer, autumn-winter, and spring-early summer. Host ranges of the two species were investigated by fruit sampling in and outside of trapping sites. Ceratitis quilicii was more widely distributed than C. rosa with the latter being recorded in only three north-eastern provinces. There were geographical limits for both species with no records of them in Northern Cape Province. Catches of C. quilicii were higher in summer with average temperatures varying from 15 to 27°C while for C. rosa, catches remained low and consistent between seasons. Ceratitis quilicii catches decreased at lower rates than those of C. rosa at temperatures below 15°C. The two species were reared from 13 plant species from nine families. Four of these hosts were infested by both C. quilicii and C. rosa in the same province where they occurred. Preferred hosts of the two species belonged to the Myrtaceae family. The characterisation of the distribution, abundance and host ranges of these pests will provide a baseline for pest status determination and implementation of management actions.

A generalized risk assessment index for forecasting insect population under the effect of temperature

Life history traits have been studied under various environmental factors, but the ability to combine them into a simple function to assess pest response to climate is still lacking complete understanding. This study proposed a risk index derived by combining development, mortality, and fertility rates from a stage-structured dynamic mathematical model. The first part presents the theoretical framework behind the risk index. The second part of the study is concerned with the application of the index in two case studies of major economic pest: the brown planthopper (Nilaparvata lugens) and the spotted wing drosophila (Drosophila suzukii), pests of rice crops and soft fruits, respectively. The mathematical calculations provided a single function composed of the main thermal biodemographic rates. This function has a threshold value that determines the possibility of population increase as a function of temperature. The tests carried out on the two pest species showed the capability of the index to describe the range of favourable conditions. With this approach, we were able to identify areas where pests are tolerant to climatic conditions and to project them on a geospatial risk map. The theoretical background developed here provided a tool for understanding the biogeography of Nilaparvata lugens and Drosophila suzukii. It is flexible enough to deal with mathematically simple (N. lugens) and complex (D. Suzukii) case studies of crop insect pests. It produces biologically sound indices that behave like thermal performance curves. These theoretical results also provide a reasonable basis for addressing the challenge of pest management in the context of seasonal weather variations and climate change. This may help to improve monitoring and design management strategies to limit the spread of pests in invaded areas, as some non-invaded areas may be suitable for the species to develop.

Fruit Flies: Challenges and Opportunities to Stem the Tide of Global Invasions

Global trade in fresh fruit and vegetables, intensification of human mobility, and climate change facilitate fruit fly (Diptera: Tephritidae) invasions. Life-history traits, environmental stress response, dispersal stress, and novel genetic admixtures contribute to their establishment and spread. Tephritids are among the most frequently intercepted taxa at ports of entry. In some countries, supported by the rules-based trade framework, a remarkable amount of biosecurity effort is being arrayed against the range expansion of tephritids. Despite this effort, fruit flies continue to arrive in new jurisdictions, sometimes triggering expensive eradication responses. Surprisingly, scant attention has been paid to biosecurity in the recent discourse about new multilateral trade agreements. Much of the available literature on managing tephritid invasions is focused on a limited number of charismatic (historically high-profile) species, and the generality of many patterns remains speculative.

Differential Cold Tolerance on Immature Stages of Geographically Divergent Ceratitis capitata Populations

Cold tolerance of adult medflies has been extensively studied but the effect of subfreezing temperatures on the immature stages remains poorly investigated, especially as far as different populations are regarded. In this study, we estimated the acute cold stress response of three geographically divergent Mediterranean fruit fly populations originating from Greece (Crete, Volos) and Croatia (Dubrovnik) by exposing immature stages (eggs, larvae, pupae) to subfreezing temperatures. We first determined the LT50 for each immature stage following one hour of exposure to different temperatures. Then eggs, larvae and pupae of the different populations were exposed to their respective LT50 for one hour (LT50 = -11 °C, LT50 = -4.4 °C, LT50 = -5 °C for eggs, larvae and pupae, respectively). Our results demonstrate that populations responded differently depending on their developmental stage. The population of Dubrovnik was the most cold-susceptible at the egg stage, whereas in that of Crete it was at the larval and pupal stage. The population of Volos was the most cold-tolerant at all developmental stages. The egg stage was the most cold-tolerant, followed by pupae and finally the 3rd instar wandering larvae. This study contributes towards understanding the cold stress response of this serious pest and provides data for important parameters that determine its successful establishment to unfavorable environments with an emphasis on range expansion to the northern, more temperate regions of Europe.

A comparative analysis of the chromosomes of three FARQ species complex members, Ceratitis rosa, C. quilicii, and C. fasciventris F2 (Diptera: Tephritidae)

The Ceratitis FARQ species complex consists of four highly destructive agricultural pests of Africa, namely C. fasciventris, C. anonae, C. rosa, and C. quilicii. The members of the complex are considered very closely related and the species limits among them are rather obscure. Their economic significance and the need for developing biological methods for their control makes species identification within the complex an important issue, which has become clear that can only be addressed by multidisciplinary approaches. Chromosomes, both mitotic and polytene, can provide a useful tool for species characterization and phylogenetic inference among closely related dipteran species. In the current study, we present the mitotic karyotype and the polytene chromosomes of C. rosa and C. quilicii together with in situ hybridization data. We performed a comparative cytogenetic analysis among the above two species and C. fasciventris, the only other cytogenetically studied member of the FARQ complex, by comparing the mitotic complement and the banding pattern of the polytene chromosomes of each species to the others, as well as by studying the polytene chromosomes of hybrids between them. Our analysis revealed no detectable chromosomal rearrangements discriminating the three FARQ members studied, confirming their close phylogenetic relationships.

Integrating temperature-dependent development and reproduction models for predicting population growth of the coffee berry borer, Hypothenemus hampei Ferrari

The coffee berry borer, Hypothenemus hampei Ferrari (Coleoptera: Curculionidae, Scolytinae), is the most devastating insect pest of coffee worldwide. It feeds on the beans inside the berries leading to significant crop losses and unmarketable products. This study aims to model the impact of temperature on H. hampei fecundity and population growth parameters, as a contribution to the prediction of infestation risk. The fecundity was assessed on fresh coffee beans at six constant temperatures in the range 15-30°C, with RH 80 ± 5% and photoperiod 12:12 L:D. Nonlinear models were fitted to the relationship between fecundity and temperature using the ILCYM software. The best fecundity model was combined to development models obtained for immature stages in a previous study in order to simulate life table parameters at different constant temperatures. Females of H. hampei successfully oviposited in the temperature range 15-30°C, with the highest fecundity observed at 23°C (106.1 offspring per female). Polynomial function 8 model was the best fitted to the relationship between fecundity and temperature. With this model, the highest fecundity was estimated at 23°C, with 110 eggs per female. The simulated net reproductive rate (R₀) was maximal at 24°C, with 50.08 daughters per female, while the intrinsic rate of increase (r_m) was the highest at 26°C, with a value of 0.069. Our results will help understand H. hampei population dynamics and develop an ecologically sound management strategy based on a better assessment of infestation risk.

Temperature-dependent phenology of the parasitoid Fopius arisanus on the host Bactrocera dorsalis

Fopius arisanus (Sonan), an egg-pupal parasitoid of numerous fruit fly species, was recently introduced into Africa for the control of the Oriental fruit fly, Bactrocera dorsalis (Hendel). In this study, life-table data of F. arisanus were generated under laboratory conditions at six constant temperatures (15, 20, 25, 30, 28 and 35 °C; 75% RH, L12:D12 photoperiod) and under variable conditions in a screenhouse on B. dorsalis. Several non-linear functions were fitted to model species development, mortality, longevity and oviposition using the Insect Life Cycle Modeling (ILCYM) software. The established phenology models were stochastically simulated at variable temperatures to estimate the life table parameters. Fopius arisanus completed development from egg to adult at all the temperatures tested except at 35 °C. Mean developmental time from egg to adult was inversely proportional to the temperature. The minimum temperature threshold (Tmin) from egg to adult was 8.15 °C, while the maximum temperature threshold (Tmax) was estimated at 34.2 °C. The optimal temperature for immature stages survival was predicted for 20-30 °C. The optimum fecundity estimated was 251 eggs/female at 22 °C. Following the stochastic simulations under natural conditions of the selected agro-ecological zones, it appears that the Humid Forest with Bimodal Rainfall provides a suitable thermal condition for potential population growth of F. arisanus. The present study shows the importance of temperature on the demographic parameters of F. arisanus. Implications of present findings on the biological control of B. dorsalis under climate change scenarios is discussed.

The Sterile Insect Technique: Success and Perspectives in the Neotropics

The sterile insect technique (SIT), an environmentally friendly means of control, is currently used against plant, animal, and human pests under the area-wide integrated pest management. It consists in the mass production, sterilization, and release of insects in an affected area where sterile males mate with wild females leading to no reproduction. Here, we review SIT in the Neotropics and focus on particular recent successful cases of eradication of the Mediterranean fruit fly, Ceratitis capitata (Wiedemann), as well as effective programs used against the Mexican fruit fly Anastrepha ludens (Loew), the New World screwworm fly Cochliomyia hominivorax (Coquerel)), and the Cactus moth Cactoblastis cactorum (Berg). We examine when SIT does not work and innovations that have made SIT more efficient and also highlight complimentary techniques that can be used in conjunction. We address potential candidate species that could be controlled through SIT, for example Philornis downsi Dodge & Aitken. Finally, we consider the impact of climate change in the context of the use of the SIT against these pests. Given the recent dramatic decline in insect biodiversity, investing in environmentally friendly means of pest control should be a priority. We conclude that SIT should be promoted in the region, and leadership and political will is needed for continued success of SIT in the Neotropics.

Phylogenomic resolution of the Ceratitis FARQ complex (Diptera: Tephritidae)

The Ceratitis FARQ complex (formerly FAR complex) includes four frugivorous tephritids, Ceratitis fasciventris, C. anonae, C. rosa and C. quilicii, the latter two causing important agricultural losses in Africa. Although FARQ species can be identified on the basis of subtle morphological differences, they cannot be resolved as monophyletic when trying phylogenetic tree reconstructions based on mitochondrial or nuclear gene fragments except for microsatellites. In this study, we used mitogenome and genome-wide SNPs to investigate the phylogenetic relationship within the complex as well as between all four Ceratitis subgenera. The analysis of 13 species supported the monophyly of the Ceratitis subgenera Ceratitis, Ceratalaspis, Pardalaspis, and recovered Pterandrus as paraphyletic but could not properly resolve species within the FARQ complex. Conversely, gene and species tree reconstructions based on 785,484 genome-wide SNPs could consistently resolve the FARQ taxa and provide insights into their phylogenetic relationships. Gene flow was detected by TreeMix analysis from C. quilicii to C. fasciventris, suggesting the existence of introgression events in the FARQ complex. Our results suggest that genome-wide SNPs represent a suitable tool for the molecular diagnosis of FARQ species and could possibly be used to develop rapid diagnostic methods or to trace the origins of intercepted samples.

Τhe complete mitochondrial genomes of Ceratitis rosa and Ceratitis quilicii, members of the Ceratitis FAR species complex (Diptera: Tephritidae)

Ceratitis FAR is an African species complex comprising insect pests of great economic interest and obscure species limits. Here, we report the mitochondrial genomes of two members of the FAR complex, namely Ceratitis rosa and the recently characterized Ceratitis quilicii. A phylogenetic analysis based on PCGs of available Tephritidae mitogenomes is presented. The current mitochondrial sequences from the FAR complex could contribute toward the resolution of phylogenetic relationships and species limits within this taxonomically challenging group, which is also an important issue for the development of environment-friendly and species-specific control methods, such as the sterile insect technique (SIT).

Temperature-based phenology model to predict the development, survival, and reproduction of the oriental fruit fly Bactrocera dorsalis

The oriental fruit fly Bactrocera dorsalis (Diptera: Tephritidae) is a major pest of fruit and vegetable production systems on several continents. The pest has invaded many countries, causing considerable impact on fruit production systems and commercialization. In this study we determined the relationship between temperature and development, survival and reproductive parameters of B. dorsalis on an artificial diet under laboratory conditions under 7 constant temperatures (10, 15, 20, 25, 30, 33 and 35 °C) with 70 ± 10% relative humidity and a photoperiod of L12:D12. We validated the laboratory results with a full life table analysis under semi-natural conditions in a screenhouse. We used the Insect Life Cycle Modeling (ILCYM) software for all mathematical models and simulations applied to all life history parameters. Bactrocera dorsalis completed its development at temperatures ranging between 15 and 33 °C with the mean developmental time of egg, larva, and pupa ranging between 1.46 and 4.31 days, 7.14-25.67 days, and 7.18-31.50 respectively. The models predicted temperatures ranging between 20 and 30 °C as favorable for development and survival, and 20 to 25 °C for optimal fecundity of B. dorsalis. Life table parameters showed the highest gross reproductive rate (GRR), net reproductive rate (R_o), intrinsic rate of increase (r_m), and finite rate of increase (λ) between 25 and 31 ᵒC while generation time (T) and doubling time (D_t) were low at this interval. The effects of future climate change on B. dorsalis life history parameters were further investigated and the outcome from this study will help in the management of B. dorsalis in different agroecologies in the context of ongoing climate change.

Exploring the role of host specialisation and oxidative stress in interspecific lifespan variation in subtropical tephritid flies

In herbivorous insects, the degree of host specialisation may be one ecological factor that shapes lifespan. Because host specialists can only exploit a limited number of plants, their lifecycle should be synchronised with host phenology to allow reproduction when suitable hosts are available. For species not undergoing diapause or dormancy, one strategy to achieve this could be evolving long lifespans. From a physiological perspective, oxidative stress could explain how lifespan is related to degree of host specialisation. Oxidative stress caused by Reactive Oxygen Species (ROS) might help underpin ageing (the Free Radical Theory of Aging (FRTA)) and mediate differences in lifespan. Here, we investigated how lifespan is shaped by the degree of host specialisation, phylogeny, oxidative damage accumulation and antioxidant protection in eight species of true fruit flies (Diptera: Tephritidae). We found that lifespan was not constrained by species relatedness or oxidative damage (arguing against the FRTA); nevertheless, average lifespan was positively associated with antioxidant protection. There was no lifespan difference between generalist and specialist species, but most of the tephritids studied had long lifespans in comparison with other dipterans. Long lifespan may be a trait under selection in fruit-feeding insects that do not use diapause.

Aged virgin adults respond to extreme heat events with phenotypic plasticity in an invasive species, Drosophila suzukii

Climate warming has increased the frequency of extreme heat events. Alien species usually invade new areas with a low-density population and often have limited mating opportunities due to the unsynchronized emergence of adults. Early-emerging virgin adults often have to wait to mate with later-emerging partners at the cost of aging, which reduces thermal tolerance. To understand the adaptive strategies of virgin males/females versus those of mated males/females in response to heat stress during aging, we conducted a fully factorial experiment to test the basal and plastic heat tolerance (CTmax, critical thermal maximum) of males and females with different mating statuses (virgin and mated) at different ages (5, 10, and 15 days after eclosion) after different acclimation regimes (null, rapid and developmental heat acclimation) in a well-known invasive species, Drosophila suzukii. We found that mating could change the heat tolerance of adults during aging. Mated females had higher basal heat tolerance than virgin females, while mated males had lower tolerance than virgin males. Mating could generally decrease the acclimation capacity (i.e., plasticity of heat tolerance) during aging. Aged virgin adults had a much higher acclimation capacity than aged mated adults. Our findings suggest that phenotypic plasticity of heat tolerance may be a main strategy used by virgin adults to cope with heat events. The phenotypic plasticity of thermal tolerance could increase the invasion success of alien species in new areas by allowing them to rapid respond to local temperature changes.

Integrating temperature-dependent life table data into Insect Life Cycle Model for predicting the potential distribution of Scapsipedus icipe Hugel & Tanga

Scapsipedus icipe Hugel and Tanga (Orthoptera: Gryllidae) is a newly described edible cricket species. Although, there is substantial interest in mass production of S. icipe for human food and animal feed, no information exists on the impact of temperature on their bionomics. Temperature-dependent development, survival, reproductive and life table parameters of S. icipe was generated and integrated into advanced Insect Life Cycle Modeling software to describe relative S. icipe population increase and spatial spread based on nine constant temperature conditions. We examined model predictions and implications for S. icipe potential distribution in Africa under current and future climate. These regions where entomophagy is widely practiced have distinctly different climates. Our results showed that S. icipe eggs were unable to hatch at 10 and 40°C, while emerged nymphs failed to complete development at 15°C. The developmental time of S. icipe was observed to decrease with increased in temperature. The lowest developmental threshold temperatures estimated using linear regressions was 14.3, 12.67 and 19.12°C and the thermal constants for development were 185.2, 1111.1- and 40.7-degree days (DD) for egg, nymph and pre-adult stages, respectively. The highest total fecundity (3416 individuals/female/generation), intrinsic rate of natural increase (0.075 days), net reproductive rate (1330.8 female/female/generation) and shortest doubling time (9.2 days) was recorded at 30°C. The regions predicted to be suitable by the model suggest that S. icipe is tolerant to a wider range of climatic conditions. Our findings provide for the first-time important information on the impact of temperature on the biology, establishment and spread of S. icipe across the Africa continent. The prospect of edible S. icipe production to become a new sector in food and feed industry is discussed.

Life table parameters of the red palm mite Raoiella indica (Acari: Tenuipalpidae) at various temperatures and for sexual and asexual reproduction

The red palm mite, Raoiella indica Hirst, is a threat to coconut, banana and native Arecaceae and Heliconiaceae in Brazil. This mite originated in the Eastern Hemisphere and was first reported in 2004 in the Americas, where the pest is spreading quickly and causing severe damage to its host plants. The objective of this work was to determine the life-history parameters of R. indica at constant temperatures, estimate its thermal requirements [threshold temperature (Tb) and thermal constant (K)] and also compare its life table parameters between sexual reproduction and parthenogenesis. The life tables were constructed on leaflets of Adonidia merrillii at 15, 20, 24, 27, 30 and 34 °C and 65% RH and a 12-h photoperiod. The longevity and the number of laid eggs of non-copulated adult females were evaluated at 27 °C. Raoiella indica had complete development, from egg to adult, only at 20, 24, 27 and 30 °C. At 15 °C, the eggs did not hatch, and at 34 °C, the mites survived only until the larval stage. For sexual reproduction, the optimal temperature was 27 °C, under which the reproductive parameters were higher. The reproductive parameters for sexual reproduction were higher than those for parthenogenesis. The Tb was 14.79 °C, and the thermal constant was 208.33 degree days. The life parameters estimated in this study can be used for modelling and predicting the population growth of R. indica in the field and consequently for improving their management strategies.

Effect of temperature on the biological parameters of the cabbage aphid Brevicoryne brassicae

The cabbage aphid, Brevicoryne brassicae, is a pest of many plants of the Brassicaceae family including cabbage, Brassica oleracea Linnaeus, 1753. We investigated the effect of temperature on the biological parameters of B. brassicae using different temperature-based models incorporated in the Insect Life Cycle Modelling (ILCYM) software. Nymphs of first stage were individually placed in the incubators successively set at 10°C, 15°C, 20°C, 25°C, 30°C, and 35°C; 75 ± 5% RH; and L12: D12-hr photoperiods. We found that first nymph reached the adult stage after 18.45 ± 0.04 days (10°C), 10.37 ± 0.26 days (15°C), 6.42 ± 0.07 days (20°C), 5.076 ± 0.09 days (25°C), and 5.05 ± 0.10 days (30°C), and failed at 35°C. The lower lethal temperatures for B. brassicae were 1.64°C, 1.57°C, 1.56°C, and 1.62°C with a thermal constant for development of 0.88, 0.87, and 0.08, 0.79 degree/day for nymphs I, II, III, and IV, respectively. The temperatures 10, 30, and 35°C were more lethal than 15, 20, and 25°C. Longevity was highest at 10°C (35.07 ± 1.38 days). Fertility was nil at 30°C and highest at 20°C (46.36 ± 1.73 nymphs/female). The stochastic simulation of the models obtained from the precedent biological parameters revealed that the life table parameters of B. brassicae were affected by the temperature. The net reproduction rate was highest at 20°C and lowest at 30°C. The average generation time decreased from 36.85 ± 1.5 days (15°C) to 6.86 ± 0.1 days (30°C); the intrinsic rate of increase and the finite rate of increase were highest at 25°C. In general, the life cycle data and mathematical functions obtained in this study clearly illustrate the effect of temperature on the biology of B. brassicae. This knowledge will contribute to predicting the changes that may occur in a population of B. Brassiace in response to temperature variation.

Threshold temperatures and thermal requirements of black soldier fly Hermetia illucens: Implications for mass production

Efforts to recycle organic wastes using black soldier fly (BSF) Hermetia illucens into high-nutrient biomass that constitutes a sustainable fat (biodiesel) and high-quality protein ingredient in animal feeds have recently gained momentum worldwide. However, there is little information on the most suitable rearing conditions for growth, development and survivorship of these flies, which is a prerequisite for mass production technologies. We evaluated the physiological requirements for growth and reproduction of H. illucens on two diets [spent grains supplemented with brewers’ yeast (D1) and un-supplemented (D2)]. Development rates at nine constant temperatures (10–42°C) were fitted to temperature-dependent linear and non-linear day-degree models. Thereafter, life history table parameters were determined within a range of favourable temperatures. The thermal maximum (TM) estimates for larval, pre-pupal and pupal development using non-linear model ranged between 37.2 ± 0.3 and 44.0 ± 2.3°C. The non-linear and linear day-degree model estimations of lower developmental temperature threshold for larvae were 11.7 ± 0.9 and 12.3 ± 1.4°C for D1, and 10.4 ± 1.7 and 11.7 ± 3.0°C for D2, respectively. The estimated thermal constant of immature life stages development of BSF was higher for the larval stage (250±25 DD for D1 and 333±51 for D2) than the other stages evaluated. Final larval wet weight was higher on D1 compared to D2. The population growth rate was most favourable at 30-degree celsius (°C) with higher intrinsic rate of natural increase (r_m = 0.127 for D1 and 0.122 for D2) and shorter doubling time (5.5 days for D1 and 5.7 days for D2) compared to the other temperatures. These results are valuable for the optimization of commercial mass rearing procedures of BSF under various environmental conditions and prediction of population dynamics patterns using computer simulation models.

Prediction of insect pest distribution as influenced by elevation: Combining field observations and temperature-dependent development models for the coffee stink bug, Antestiopsis thunbergii (Gmelin)

The antestia bug, Antestiopsis thunbergii (Gmelin 1790) is a major pest of Arabica coffee in Africa. The bug prefers coffee at the highest elevations, contrary to other major pests. The objectives of this study were to describe the relationship between A. thunbergii populations and elevation, to elucidate this relationship using our knowledge of the pest thermal biology and to predict the pest distribution under climate warming. Antestiopsis thunbergii population density was assessed in 24 coffee farms located along a transect delimited across an elevation gradient in the range 1000-1700 m asl, on Mt. Kilimanjaro, Tanzania. Density was assessed for three different climatic seasons, the cool dry season in June 2014 and 2015, the short rainy season in October 2014 and the warm dry season in January 2015. The pest distribution was predicted over the same transect using three risk indices: the establishment risk index (ERI), the generation index (GI) and the activity index (AI). These indices were computed using simulated life table parameters obtained from temperature-dependent development models and temperature data from 1) field records using data loggers deployed over the transect and 2) predictions for year 2055 extracted from AFRICLIM database. The observed population density was the highest during the cool dry season and increased significantly with increasing elevation. For current temperature, the ERI increased with an increase in elevation and was therefore distributed similarly to observed populations, contrary to the other indices. This result suggests that immature stage susceptibility to extreme temperatures was a key factor of population distribution as impacted by elevation. In the future, distribution of the risk indices globally indicated a decrease of the risk at low elevation and an increase of the risk at the highest elevations. Based on these results, we concluded with recommendations to mitigate the risk of A. thunbergii infestation.