Skip to main content

Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication

Abstract

Background

The highly polyphagous Queensland fruit fly (Bactrocera tryoni Froggatt) expanded its range substantially during the twentieth century and is now the most economically important insect pest of Australian horticulture, prompting intensive efforts to develop a Sterile Insect Technique (SIT) control program. Using a “common garden” approach, we have screened for natural genetic variation in key environmental fitness traits among populations from across the geographic range of this species and monitored changes in those traits induced during domestication.

Results

Significant variation was detected between the populations for heat, desiccation and starvation resistance and wing length (as a measure of body size). Desiccation resistance was correlated with both starvation resistance and wing length. Bioassay data for three resampled populations indicate that much of the variation in desiccation resistance reflects persistent, inherited differences among the populations. No latitudinal cline was detected for any of the traits and only weak correlations were found with climatic variables for heat resistance and wing length. All three stress resistance phenotypes and wing length changed significantly in certain populations with ongoing domestication but there was also a strong population by domestication interaction effect for each trait.

Conclusions

Ecotypic variation in heat, starvation and desiccation resistance was detected in Australian Qfly populations, and these stress resistances diminished rapidly during domestication. Our results indicate a need to select source populations for SIT strains which have relatively high climatic stress resistance and to minimise loss of that resistance during domestication.

Background

Climate plays a major role in determining the geographical distribution of species, and their climate adaptability is considered a key component of their expansion and invasion potential [1, 2]. The Queensland fruit fly (Bactrocera tryoni Froggatt) (‘Qfly’) is native to subtropical regions of eastern Australia but has also established populations in the Northern Territory and down the east coast to temperate regions of Victoria, with some invasive populations also found in New Caledonia, French Polynesia, Pitcairn Island and Cook Island [3, 4]. As its range has expanded so has its pest status; it is now the major insect pest of Australian horticulture [5, 6]. The species’ wide climatic and geographical range suggests a high level of climate adaptability.

Latitudinal clines in climatic stress resistance and life history traits are well documented in insects, particularly in Drosophila melanogaster [7]. In some cases, such as body size variation in Drosophila, strong clinal signals are detected in different continents [8,9,10,11,12]. These parallel geographical patterns reflect repeatable adaptive genetic responses to climatic selection. In contrast, latitudinal patterns of starvation resistance in D. melanogaster differ among continents: a strong cline is present in India [13], but absent in South America [14] and Australia [15]. Furthermore, while linear clines are abundant in insects, non-linear patterns also exist (e.g., overwinter egg production and longevity in D. melanogaster [16]). In addition to these classical examples in D. melanogaster, clinal variations have also been observed in other insects such as D. birchii and D. serrata in eastern Australia [7, 17,18,19] and D. subobscura, Eldana saccharina, Glossina pallidipes, and Anopheles gambiae in other continents [20,21,22,23]. In several cases the phenotypic variation in traits such as climate stress resistance has been associated with gradients in specific climatic variables [17, 19,20,21,22, 24,25,26,27,28] and, in physiological terms, with differences in body size and sometimes also developmental rates [10, 12, 17, 18, 29, 30].

In Qfly, survival and reproduction are heavily influenced by temperature, moisture and availability of suitable host fruits [31,32,33,34,35], and Fletcher [36, 37] has suggested that the populations now persisting in the temperate regions do so in part by short distance movements between orchards and nearby water sources. Desiccation stress is therefore speculated to be a major determinant of Qfly distributions [32, 35, 38], although there is as yet no empirical data in the literature about variation in desiccation resistance amongst wild populations of the species. Not much is known about variation in Qfly thermal resistance either. Bateman [39] reported that northern populations of Qfly have higher fecundity and survival rates at higher temperatures (30 °C) than at the standard laboratory rearing conditions of 25 °C. However, he acknowledged some possible bias in his work through the inclusion of what is now recognised as a sibling species, B. neohumeralis (which he termed ‘variety neohumeralis’), in some of his collections. Furthermore, Meats [32] found no difference in cold resistance between northern and southern Qfly collections and he later reported rapid acclimation to low temperatures in several populations in the laboratory. Overall the current literature on geographic variation in the climatic stress resistances of Qfly is fragmentary and provides insufficient detail for understanding the role of such variation in the ecology and invasive potential of the species, and some of the literature is not based on current taxonomy [6, 38].

Better understanding of genetic variation in Qfly climate stress resistance has direct relevance to the Sterile Insect Technique (SIT) programmes now being implemented to suppress established populations and eliminate new outbreaks in key horticultural regions. These bisexual SIT programmes involve the mass rearing, gamma ray sterilisation and mass release of a domesticated Qfly strain into outbreak areas, resulting in wild flies wasting their reproductive effort on non-productive mating. These programmes are now being carried out in several regions of south-eastern Australia [40]. Climate stress resistances (along with other life history traits outside the scope of this paper such as mating competitiveness [41,42,43], lifespan [44] and predator evasion [45]) have been identified as important for the success of these programmes [46,47,48]. Some climate stress resistance is lost during the domestication process that occurs when the large numbers of Qflies that need to be produced for irradiation and release are mass reared in factory settings for multiple generations [46, 49, 50]. Changes in various traits during domestication have also been observed in strains of Mexican fruit fly Anastrepha ludens used in SIT programmes [51]. Given the harsh conditions into which the flies are often released, it is important to select strains for mass rearing from relatively stress resistant populations and to maintain as much of that resistance through multiple generations of mass rearing as possible.

The present study adopts a “common garden” approach to investigate genetic variation in heat, cold, desiccation and starvation resistance in Qfly populations of diverse geographical origins. We also explore how the variation in those traits relates to wing length (as a measure of body size), geographic origin (latitude and coastal vs inland) and weather variables. Additionally, we investigate how the climate stress resistances change during domestication and whether those changes differ between populations. Populations from three sites showing desiccation resistance differences in the initial survey were resampled in subsequent years and rescored for that trait to determine the genetic stability of the population differences observed. We discuss our results in relation to the ecology of this invasive pest and implications for SIT-based pest control programmes.

Results

Heat resistance

Our primary survey of 12 populations from widely separated sites (Fig. 1, Table 1, Additional file 1 Table S1) found significant population differences at G2/G3 (\( {\chi}_{11}^2 \) = 5.17, dispersion (ϕ) = 0.098, P < 0.001) in heat resistances as measured by knockdown time at 42 °C (Fig. 2). These differences were not explained by latitude (t8 = 1.72, P > 0.05), or coastal vs inland origin of the populations (t8 = 0.84, P > 0.05). There was a significant association between the differences in heat resistance and the geographical distance between populations (Mantel’s r = 0.46, P = 0.02), although it only explained 21% of the variance in resistance. The resistance differences were largely due to the relatively long knockdown times (i.e., higher resistances) of the Darwin, Alice Springs and Sydney populations (on average, 32% longer than those of the other populations).

Fig. 1
figure 1

Geographical distribution of Qfly populations studied. Climatic zones are as described in the Köppen climate classification scheme for the years 1980 to 2016 as updated by Beck et al. [85]. Hatched areas encompass the known distribution of Qfly across Australia. Note that Qfly populations in the arid inland regions of Australia are largely confined to townships and irrigated areas. The map was created with ‘ggplot2’ [109] in R using climate zones data from Beck et al. [85] under the Creative Commons Attribution 4.0 International Licence

Table 1 Populations studied in the primary and resampled surveys
Fig. 2
figure 2

Variation in heat resistance of wild and domesticated Qfly populations. Data are first presented as normalised median knockdown times in minutes on exposure to 42 °C for each population in the generations designated as wild (G2/G3) and domesticated (G11–15). The data for both the wild and domesticated generations are then also shown as ratios against the S06 control used in the same batch of bioassays. Finally, the ratios of the median knockdown times for the domesticated vs wild generations for each population are given. Standard errors for the knockdown times were on average 5.34 and 2.55 for the wild and domesticated populations respectively. Note that not all populations were rescored at G11–15

Regressing the median knockdown time against the five key weather variables yielded one significant association, a positive association with the minimum temperature of the coldest month, and this was marginally significant (t10 = 2.26, P = 0.047). As this variable was correlated with several other weather variables not used in the model (annual maximum temperature, annual minimum temperature, annual temperature, minimum temperature of the coldest month, and precipitation of the wettest month, Additional file 1 Fig. S1), it was unclear which aspect(s) of climate was causally involved.

The heat resistances of the Darwin, Alice Springs and Sydney populations at G2/G3 were all slightly higher than that of the long term domesticated S06 control, while the resistances of all the other populations at that point were slightly lower than that of S06 (Fig. 2). However, by G12–15 Darwin, Alice Springs and Sydney had all lost resistance and all the rescored populations then showed lower (20%) resistance than S06 (Fig. 2), although some populations showed no significant change during domestication. No significant population differences remained by G12–15 (\( {\chi}_8^2 \) = 1.72, ϕ = 0.19, P > 0.05).

Cold resistance

Cold resistance (as measured by shorter chill coma recovery time) was not correlated with heat resistance (Table 2) and, unlike heat resistance, there were no significant differences between populations in cold resistance at G2/G3 (Fig. 3; \( {\chi}_{11}^2 \) = 1.55, ϕ = 0.10, P > 0.05). Nor were there significant population differences among the G12–15 flies (\( {\chi}_8^2 \) = 0.87, ϕ = 0.14, P > 0.05), or domestication or population by domestication interaction effects (Fig. 3; \( {\chi}_8^2 \) = 1.12, ϕ = 0.11, P > 0.05, and \( {\chi}_8^2 \) = 0.48, ϕ = 0.11, P > 0.05, respectively). Neither the G2/G3 nor G12–15 data for the test populations showed consistent differences from the corresponding S06 data, although the G12–15 data were generally closer to S06 than were the G2/G3 data.

Table 2 Spearman correlations between stresses and wing length across wild (G2/G3) and domesticated (G10–15) samples
Fig. 3
figure 3

Variation in cold resistance of wild and domesticated Qfly populations. Data are presented as recovery times in minutes following exposure to 0 ± 0.5 °C for 16 h. Standard errors for recovery times were on average 18.6 and 6.1 for the wild and domesticated populations respectively. The format of the figure otherwise follows Fig. 2

Desiccation and starvation resistance

Desiccation resistance (as measured by longer survival time under desiccating conditions) was not correlated with either heat or cold resistance (Table 2). However, similar to heat resistance, we found significant variation in desiccation resistance among the 12 G2/G3 populations (Fig. 4; \( {\chi}_{11}^2 \) = 6.77, ϕ = 0.08, P < 0.001). These differences were not explained by latitude (t8 = 0.97, P > 0.05), the coastal vs inland origins of the populations (t8 = 0.04, P > 0.05), the weather variables (t8 = 1.95, P > 0.05) or geographical distances between them (Mantel’s r = − 0.11, P > 0.05). The main outlier populations were Sydney, again, and to a lesser extent Griffith and Cape Tribulation, which all showed higher levels of resistance than the other populations (on average, 44.3, 21.8 and 19.4% longer survival times respectively; Fig. 4). All 12 G2/G3 test populations showed higher desiccation resistance than S06 (> 2-fold for Sydney and an average of 55% across all populations).

Fig. 4
figure 4

Variation in desiccation and starvation resistance of wild and domesticated Qfly populations. Data are presented as survival times in hours in the presence (desiccation) and absence (starvation) of desiccant (silica gel beads). Standard errors of survival times were on average 11.30 and 5.37 for starvation, and 5.76 and 6.08 for desiccation, in the wild and domesticated populations respectively. The format of the figure otherwise follows Fig. 2

There was still significant variation in desiccation resistance between populations in the G10/G11 samples (\( {\chi}_7^2 \) = 9.65, ϕ = 0.11, P < 0.001), and most of these samples still had substantially higher resistance than S06 (average 32% longer survival times; Fig. 4). However, the differences between populations at G10/G11 were not the same as those at G2/G3 above, which was reflected in a highly significant population by domestication interaction term (\( {\chi}_7^2 \) = 10.33, ϕ = 0.08, P < 0.001). In particular, Sydney, which had the highest resistance at G2, had relatively low resistance at G10, at which point it was not significantly different from the S06 control. The three individual populations for which the G2/G3 and G10/G11 data differed significantly were Sydney and Alice Springs, both showing decreases over time (on average 54 and 68% shorter survival times respectively), and Narrabri, which showed an increase (on average 23% longer survival time; Additional file 1 Table S2).

Starvation resistance (as measured by longer survival time without access to food or water) was not correlated with heat or cold resistance but was positively correlated with desiccation resistance (Spearman’s correlation r = 0.62, P < 0.01; Table 2). Some correlation with the latter was expected since ‘desiccation’ was also measured in the absence of food and water, although the correlation observed implies that only 38% of the variance in the two measures was shared between them. There were significant population differences at G2/G3 in starvation resistance (\( {\chi}_{11}^2 \) = 11.32, ϕ = 0.17, P < 0.05), largely due to the higher resistance of the Cape Tribulation and Sydney populations (on average 75 and 53% longer survival, respectively). However, there was no association of resistance with latitude (t8 = 0.18, P > 0.05), coastal vs inland origin (t8 = 0.49, P > 0.05), weather variables (P > 0.05 for all t8 tests) or geographical distance between populations (Mantel’s r = − 0.02, P > 0.05). The G2/G3 samples all showed higher starvation resistances than S06 (up to 2.7-fold and an average of 73% longer survival) (Fig. 4).

Starvation resistances had decreased in all except the Batemans Bay population by G10/G11 but were still higher than S06 at that point (up to 1.7-fold) (Fig. 4). However, Sydney was the only population in which the difference (decrease) was statistically significant (Additional file 1 Table S2). As with desiccation resistance, there was a significant main effect of population at G10/G11, as well as a significant population by domestication interaction (\( {\chi}_7^2 \) = 4.91, ϕ = 0.16, P < 0.05 and \( {\chi}_7^2 \) = 4.18, ϕ = 0.16, P < 0.05, respectively).

Intermediate generations

An analysis including data for the intermediate generations (i.e., G4, G6 and G8) for eight populations again showed a positive correlation between desiccation and starvation resistances, albeit only explaining about 10% of variance (Spearman correlation r = 0.31, P < 0.01). Much of the correlation reflected a sharp drop over time in both measures in Sydney; changes in the other populations were less pronounced and not all were decreases (Fig. 5). Both desiccation and starvation resistances decreased over time in Alice Springs and Griffith, Brisbane and Utchee Creek showed drops in starvation but not desiccation resistance, Narrabri rose in desiccation but not starvation resistance, and Batemans Bay and Mareeba showed no consistent directional change in either desiccation or starvation resistance.

Fig. 5
figure 5

Trends over generations in the desiccation and starvation resistance in eight Qfly populations. The data are presented as normalised median survival times in hours in the presence (desiccation) and absence (starvation) of desiccant (silica gel beads). The G2 and G10+ data for each population are taken from Fig. 4

Resampled populations

The Sydney, Cape Tribulation and Alice Springs sites were recollected in a subsequent growing season and rescored for desiccation resistance. The initial G2 assays of these populations scored them as the highest, medium-high and medium-low respectively for desiccation resistance. The G2 assays for the recollected samples yielded resistance values that were highly correlated with those of the original collections (Fig. 6).

Fig. 6
figure 6

Correlation between the G2 desiccation resistances of the primary and resampled collections from Cape Tribulation, Alice Springs and Sydney. The standard errors were on average 9.02 and 6.17 for the primary and resampled collections respectively

Relationships to wing length

There were significant differences in wing length among G2/G3 populations (F11, 329 = 12.61, P < 0.001), mostly due to the longer wings of the Sydney and Narrabri samples (Fig. 7). These differences were not associated with latitude (t10 = 1.68, P > 0.05), coastal vs inland origins (t10 = 0.42, P > 0.05), or geographic distances (Mantel’s r = − 0.15, P > 0.05), but there was a significant association (t10 = 3.31, P < 0.05) with a weather variable, maximum temperature of the warmest month. As with the association of heat resistance with minimum temperature of the coldest month above, maximum temperature of the warmest month was also correlated with other weather variables not tested in this analysis (Additional file 1 Fig. S1) and we cannot determine from our results what aspect of climate was causally involved in the association.

Fig. 7
figure 7

Variation in body size in wild and domesticated Qfly populations. Data are presented as median wing length in millimetres for each population in the generations designated as wild (G2) and domesticated (G10/G11). The format of the figure otherwise follows Fig. 2. BB: Batemans Bay; MA: Mareeba; SY: Sydney. Standard errors were on average 0.087 and 0.094, for the wild and domesticated population groups respectively

We also found population differences in wing length in the G10/G11 samples (Fig. 7; F8, 259 = 41.09, P < 0.001), but they were generally not the same differences as at G2/G3 (F8, 588 = 16.51, P < 0.001) and there was no consistent direction to the changes seen: wing length increased in Mareeba, Griffith and Batemans Bay but decreased in Brisbane and Sydney.

Overall, neither heat, cold or desiccation resistance measurements were significantly correlated with wing length in either the G2/G3 or G10/G11 samples (all regression models with P > 0.05; Fig. 8). However, the change in starvation resistance in the period between the two sampling times did show a positive correlation with the change in wing length in this period (t6 = 2.91, P < 0.05), mostly due to the relatively large increases in both measures in Mareeba and Batemans Bay and relatively large decreases in Sydney (Fig. 8). While the corresponding correlation between the changes in desiccation resistance and wing length was not significant across all populations, it is notable that the Sydney population again changed from the highest score for both measures at G2/G3 to much lower scores for both at G10/G11 (Fig. 8).

Fig. 8
figure 8

Associations between changes in stress resistances and wing length during domestication. Ratios of heat, desiccation and starvation resistances in G10–15 vs G2/G3 for each population are plotted against the corresponding ratios for wing length. The resistance data are taken from Figs. 2 and 4. BB: Batemans Bay; MA: Mareeba; SY: Sydney

Discussion

We have found genetic differences in heat, desiccation and starvation but not cold resistance among the Qfly populations which were collected in 2016 and 2017 and scored at G2/G3 in a ‘common garden’ laboratory environment. The only population survey of any of these traits reported previously for Qfly was that of Meats [32] who, consistent with our findings, found no differences in cold resistance between populations sampled from Cairns (Queensland) to Nowra (NSW). Overall, the only significant correlations among the traits we measured involved desiccation resistance, which was positively correlated with both starvation resistance and body size (measured as wing length). However, one population, Sydney, stood out on several measures, being one of three populations with the highest levels of heat resistance (together with Alice Springs and Darwin) and desiccation and starvation resistance (along with Cape Tribulation and Griffith) and also being the largest in size. Importantly, we also found that the population differences in desiccation resistance were stable, inherited features of the populations, including Sydney, which we were able to resample in the subsequent year.

Unlike the situation with heat and cold resistance and wing size in D. melanogaster populations in eastern Australia [10, 19], we found no strong association between any of our stress resistance traits and latitude or weather variables. We acknowledge that our sample size (12 populations) was too small to detect small changes with latitude. Nevertheless, it appears that Qfly has more complex ecotypic variation in these variables in eastern Australia than does D. melanogaster.

Our findings concur with the climate modelling work of Yonow and Sutherst [35], which suggested that cold is less important than heat and desiccation stress in limiting the distribution of Qfly in eastern Australia. Desiccation stress may be particularly important; it has been suggested [35] that the uptake of irrigation systems has been key to the southward expansion of the species during the twentieth century, and that without irrigation this species would still be restricted to the north-eastern coast in Queensland. However, it is unclear why Sydney should be the major outlier in terms of heat, desiccation and starvation resistance and body size. Bateman [39] did include Sydney samples in his thermal stress experiments but, as noted in the Introduction, some of his samples included B. neohumeralis. Yu et al. [52] and Cameron et al. [53] did not find genetic differences between Sydney and Brisbane populations but their work was based on a small number of microsatellite markers which would be unlikely to influence climate stress resistance or body size. Follow-up work is needed on these traits, and their genetics, in populations between Sydney and the nearest populations we sampled to the north (Brisbane, Narrabri) and south (Batemans Bay).

While other relevant literature on Sydney samples may be lacking, in broader terms our evidence for geographic differentiation among Qfly populations is corroborated by several other lines of evidence. The microsatellite studies cited above did find differences between Northern Territory, particularly Alice Springs, populations and east coast samples [52, 53] and there is also some evidence for geographically restricted populations which are relatively less melanised than most populations (e.g. form melas in the Brisbane-Rockhampton area in Queensland, which is now considered a synonym species of Qfly) [4, 54, 55]. The existence of the very closely related Bactrocera aquilonis in northern Western Australia [54, 56] (but see also [53]) also indicates a propensity of the lineage to differentiate geographically. Similar to the ecotypic differences observed in other polyphagous fruit fly species [57, 58], our abiotic stress resistance profiling of Qfly populations also indicates the presence of genetically stable ecotypes in Australia.

We found the resistance of most populations to heat and starvation stress declined to values closer to those of the long domesticated S06 control population (~ 12 years in the laboratory at the time of testing) during the average of ten generations of domestication. Domestication effects were less consistent for desiccation resistance; decreases were observed for three populations, Sydney, Cape Tribulation and Griffith, which initially showed relatively high resistance, but others showed little change or, in a couple of cases, increases. The changes in body size were also variable and, in that case, did not show any convergence of extreme values towards the mean during domestication. The only previous study of the effects of domestication on any of these traits in Qfly was Weldon et al. [59] who found that the desiccation resistance of a recent Sydney collection was higher than that of one that had been kept in the laboratory for 20 generations. This concurs with our findings for the Sydney population. Reductions in starvation and desiccation resistance during domestication are also well documented in Drosophila subobscura [60, 61] and D. melanogaster [62].

Work on a variety of insects, including some tephritid fruit flies [63, 64], has shown that various other traits also change during domestication (Hoffmann and Ross [65] for a review). Many of these, for example reductions in pre-adult development time and sexual maturation time in adults and increased fecundity, involve an increase in fitness under laboratory conditions. Some of these are also known to occur in Qfly. For instance, Meats et al. [66] found age of mating decreased after just four generations of domestication, and higher fecundity also developed subsequently. Other traits demonstrated to change as the species adapts to the laboratory include thermal preferences [67, 68], protein consumption [66], increased abundance of volatiles released during courtship [69], locomotor activity [49], and reduced dispersal rates if released into the field [50]. Gilchrist et al. [70] also found a reduction in microsatellite variability during Qfly domestication. Given the number of characters affected, it seems likely that genetic changes during domestication would be widespread across the Qfly genome.

Few of these other studies compared the effects of domestication on different populations of the species in question as we have done. One exception is Simões et al. [61] who found significant differences between D. subobscura populations in the rate of the decline in starvation resistance. Various other Drosophila studies, using genetic markers ranging from chromosome inversions through to SNPs from genome-wide resequencing, have reported some convergent allele/haplotype frequency changes among different source populations brought into a common laboratory environment but also many cases of divergent trajectories, even among replicate populations from the same source [71,72,73]. We concur with Santos et al. [74] that genetic variation involved in field fitness but essentially irrelevant to the laboratory environment is likely to diverge rapidly in laboratory lines of necessarily small founding populations.

The implications of our findings for Qfly SIT programmes are twofold. Firstly, our results show that genetic changes with significant deleterious effects on the performance of flies released into the field could accumulate within the timeframe in which a strain is scaled up for mass rearing and release. It will clearly be important to minimise this timeframe, or the rate of change. Secondly our results show that source populations for candidate SIT strains can vary substantially in key traits, and in the changes in those traits during the course of domestication. Thus, it would be prudent to evaluate populations from multiple sources and to explore rearing conditions for candidate strains which minimise any deterioration in the key traits during domestication. Notable here is the finding of Gilchrist and Meats [50] that intercrossing among domesticated Qfly strains can ameliorate some aspects of field fitness losses, although the practicality of implementing such a crossing scheme in a factory setting may be challenging. Also notable is the demonstration in the Mexican fruit fly Anastrepha ludens that directional selection regimes can be developed that increase the level of desiccation resistance of candidate SIT strains in the laboratory whilst retaining their reproductive fitness [51, 75].

Molecular genetic analyses are now needed in order to elucidate both the genetic bases for the ecotypic variation in the resistance traits we have found and the changes in these traits that occur during domestication. Such analyses could also provide the molecular markers needed for breeding and selection programmes aimed at minimising loss of resistance during domestication. Importantly a draft genome sequence is now available for Qfly which would enable these analyses [76]. Precedent work with D. melanogaster suggests polygenic control will be found, albeit with a few genes of major effect, and also allowing the possibility of some epigenetic effects [77].

Also now needed is a screen for any differences in the ‘plastic’ responses of Qfly populations to the various stresses. Abundant literature on various Drosophila species indicates that prior exposure to milder climate stress can augment climate stress resistance of individuals and that the extent of the augmentation can vary in an inherited fashion between populations [30, 78,79,80,81]. Thus, future molecular genetic work should also consider the ‘plastic’ components of the resistance phenotypes.

Conclusions

We have found significant natural differences in heat, desiccation and starvation resistance between Qfly populations of diverse geographical origins. The fact that the differences were detected in a ‘common garden’ laboratory environment after two - three generations in the laboratory strongly suggests a genetic basis for the differences. No association with latitude or coastal vs inland origins was detected in any of the three stress-related traits, and only a relatively weak association with a climate variable was found for heat resistance and wing length. This contrasts with the strong latitudinal and climatic associations found for climate stress resistance among D. melanogaster and D. serrata populations also collected from eastern Australia and suggests complex ecotypic variation exists across the species’ range. Our data corroborate earlier evidence for ecotypic variation among Qfly populations based on microsatellite markers [52, 53, 82]. In general terms the data support the idea that the success of this pest in invading a relatively wide range of climatic environments is associated with genetic variation affecting its adaptation to different climatic conditions. We have also found significant losses of heat and desiccation resistances with ongoing domestication of some but not all of the tested populations. The differences between the populations in this respect bear out the complexity of the ecotypic variation. Both the geographic variation and the changes with domestication also have practical implications for the SIT programmes now being mounted to control Qfly [83, 84]. In particular they point to the importance of developing a strain from a relatively stress resistance base population and that procedures to preserve resistance need to be implemented from the outset of mass-rearing [42, 59, 70].

Methods

Fly stocks and husbandry for the primary survey

Twelve populations of Qfly were established in the laboratory from infested fruits collected from across the Australian distribution of the species during the fruit growing seasons (largely summer) of 2016/2017. The collections sites cover a broad range of the climatic conditions (following the Köppen climate classification [85]) within the species’ range (Fig. 1, Table 1 and Additional file 1 Table S1). Infested fruits were kept in the laboratory at 24–26 °C and third instar larvae were allowed to pupate in fine-grade vermiculite. The vermiculite was kept in drainage trays to avoid drenching the pupae with juice from the decomposing fruit, but in some cases it was also necessary to lightly moisten the vermiculite to reduce the risk of pupae dehydrating [86]. Pupae were sieved out of the vermiculite approximately every 10 days for three to 4 weeks after fruit collection. Adult flies were kept in Bugdorm rearing cages (31.5 × 31.5 × 31.5 cm BugDorm 43,030 Insect Rearing Cages, MegaView Science, Taiwan). Adult flies (less than 3 days old) identified as B. tryoni by taxonomic features detailed in Drew [54] were used to establish 1–4 replicate colonies for each population, using 80–172 flies of each sex per replicate cage (Table 1).

The colonies were then maintained for up to 15 generations under the following regimen. Adult flies were provided hydrolysed yeast (Yeast Hydrolysate Enzymatic, Cat No 103304, MP Biomedicals LLC, Australia), granulated white sugar, and water (separately ad libitum). Cages were maintained at 24-26 °C and 65 ± 5% relative humidity with 12 h light per day plus 1 h of simulated dusk and dawn. The simulated dusk (low light intensity, < 100 lux) conditions were used to induce mating [87]. Eggs were collected over 3 days from 19 to 24 days old females (6–10 h oviposition per day) using as an egg laying device, a translucent 250 mL plastic bottle (Décor, Sauce dispenser, Item No 128440) with ~ 80 small punctures on one side and a slice of apple (as an attractant) hanging inside (Additional file 1 Fig. S2). The bottle was rinsed with water to collect the eggs at the bottom (if few eggs were produced the water was strained through a nappy liner (Woolworths Homebrand) placed over the bottle mouth to collect as many eggs as possible). Up to ~ 300 eggs were transferred to 35 mL of gel diet [88] in a 90 mL round plastic container (Ø76 x 25H) (Cat. No. 01C2, Chanrol Pty Ltd.) for larval development, and the lid partially closed to prevent dehydration of the diet. The plastic container with the eggs was then placed on vermiculite in another 750–1000 mL rectangular plastic container. The lid of the container holding the diet and eggs was removed 5 days later to allow larvae to exit the diet to pupate in the vermiculite. During this time the external plastic container was covered with a nappy liner (Woolworths Homebrand) secured by a lid with an 8 cm × 14 cm hole cut out of the middle.

All populations were scored for heat, cold, desiccation and starvation resistance and wing length at generations 2 or 3 (G2 or G3) (considering the flies that emerged from the infested fruit as G0), with nine of those populations also scored in subsequent generations (Table 1). Eight of those nine populations were rescored for desiccation and starvation resistances in trials conducted every 2–3 generations until G10 or G11. The same eight domesticated populations were also rescored for heat and cold resistances and wing length once at G12–15. The ninth population, from Darwin, was rescored for heat and cold resistances only, and only at G15.

We used a domesticated line, S06, originally collected in Sydney in 2006 [70] and maintained under the same conditions as the colonies above as a control for batch effects in every resistance bioassay (see below). The phenotypic characteristics of this long-term domesticated line had stabilised well before our experiments and showed no systematic change during the course of our work.

The resistance data presented herein were based on bioassays of adult males because previous studies and our own preliminary work suggested the stress resistance phenotypes could be confounded by the reproductive history of the females, probably due to changes in their resource utilisation following mating [89, 90]. The males used were 19–24 days old when tested, at which age they would still be robust and at the peak of their fertility under laboratory conditions [91,92,93,94]. We used this age range to ensure all flies would be fully mature when tested, because wild-caught Qflies generally mature later than domesticated strains (noting that the domesticated SIT Qflies are mostly mature at 8 days) [66]. The early generation colonies were not highly productive, necessitating pooling of flies from different replicate colonies for the bioassays. However as fly fecundity increased over time [66, 70], from G6 onwards replicate colonies could be assayed independently for desiccation and starvation resistances.

Following well established precedents in the literature for various higher dipterans, we assessed thermal and desiccation resistances with heat knockdown and chill coma recovery time assays and desiccation survival time assays respectively [17, 19, 26,27,28,29, 95,96,97]. Our desiccation survival time assays closely followed a published method for Qfly, including the use of parallel ‘control’ assays of starvation survival times as well [46]. Our heat knockdown time and chill coma recovery time assays were based on the methods used to investigate latitudinal clines for those traits in eastern Australian D. melanogaster and Drosophila serrata [95,96,97].

Heat knockdown (heat resistance) assays

Heat resistance assays were based on the protocol of Sgrò et al. [96]. Thirty adult males from each population were placed individually in 5 mL polypropylene tubes (Ref. 60.992.523, Sarstedt, Australia), submerged in a transparent Perspex water bath preheated to 42 °C, and knockdown time recorded for each fly. Pilot experiments showed 42 °C knocked down the great majority of flies in less than 1.5 h whilst still giving a broad distribution of knockdown times within that interval (Additional file 1 Fig. S3). Knockdown time was defined as the time it took a fly to become immobilised and unresponsive to illumination by an LED flash light (Mini LED Torch, Cat No. TOR9LEDBK, Officeworks, Australia). Trials were stopped after 1.5 h, and any flies not knocked down at that point (< 1% of the total) were scored with the stoppage time. The same experimental set up was used to document changes in heat resistance in domesticated colonies after 12–15 generations of laboratory adaptation (Table 1).

Chill coma recovery (cold resistance) assays

To measure cold resistance thirty adult males from each population were placed individually in 2 mL polypropylene tubes. The tubes were then sandwiched between two identical 70-well Perspex plates (297 × 210 × 13 mm), held together by screws at the corners (Additional file 1 Fig. S4). Each unit was then sealed inside two waterproof bags and submerged in ice water (0 ± 0.5 °C) for 16 h. In a given trial, equal numbers of males of each population were randomly distributed into each unit. At the end of the 16 h the units were removed from the ice water and placed on a white topped bench at room temperature and lighting. The tubes were quickly removed, and the flies gently moved to the centre of their respective wells. A transparent, 3 mm thick Perspex sheet (297 × 210 × 3 mm) was then placed onto the top of each unit. Recovery time was defined as the time it took the flies to stand on their legs. Flies (< 1% of the total) that had not recovered at our final scoring time (300 mins after the 16 h cold treatment) were assigned the recovery time of 300 mins, following Gilbert et al. [97] and Hoffman et al. [98]. The Batemans Bay, Bega Valley, Canberra, Griffith and Mareeba populations were measured at G3, with the other populations tested at G2. The same experimental setup was used to monitor changes in cold resistance in domesticated colonies after 12–15 generations of laboratory adaptation (Table 1).

Desiccation and starvation resistance assays

Desiccation and starvation resistances were measured following the protocol of Weldon and Taylor [46], with minor modifications. Two cohorts of 30 males from each population or colony being tested in a given trial were transferred individually to 5 mL polypropylene tubes. For one of the cohorts, the desiccation treatment, the lid of each tube contained eight silica gel beads (Stock keeping unit (SKU): SBR-P, Silica Gel Australia) which were held in place and physically separated from the fly by cotton wool. The lids for the other cohort, the starvation treatment, had cotton wool but not silica gel beads. The tubes were then returned to standard rearing conditions and mortality scored every 8 h for the following 16 h, then every 2 h for a further 76 h, and finally every 4 h until the last fly died. Flies were scored as dead if no movement was observed after flicking the tubes. A red-filtered torch light was used to score flies during the night cycles.

Wing length measurement

Studies on a variety of insect species have found variation in climate stress resistance within the species is often correlated with body size, using body weight, wing length or wing centroid size as a proxy for body size [17, 18, 24, 46, 99, 100]. The structure and magnitude of our experiments prohibited us measuring any of these traits on treatment flies, so we measured wing length on frozen material as a surrogate, since wing length has been shown to correlate strongly with body size and weight in Qfly [101, 102]. Wing length was measured on 10 to 12 randomly selected males (frozen at the end of each generation scored, i.e., G2/G3 and G12-G15) per replicate cage for each population. The right wing of each fly was mounted onto a microscope slide. A Leica (M205 A) stereomicroscope fitted with a Leica DFC500 digital camera (Leica Microsystems) was used to photograph individual fly wings. Wing length was measured in mm from the proximal edge of the basal radial cell to the intersection of the costal and R4 + 5 veins on the margin of the wing. All measurements were done using the Leica Application Suite (V 4.12.0, Build 86) with the live interactive measurement module.

Resampling populations for desiccation resistance

Populations from three sites were resampled across two subsequent summers and rescored for desiccation resistance at G2. Colony establishment and husbandry were as described above, and desiccation resistance was also scored essentially as above, the only exception being some minor differences in desiccant presentation for the 2017/2018 collections, as detailed in Additional file 1 Table S3. Once again S06 was also scored in each trial to enable comparisons across trials.

Statistical analysis

As noted, resistance values for S06 flies were used to standardise the corresponding data from the experimental flies in each batch of bioassays. The structure of the experiment necessarily meant that not all populations and generations could be assayed in a single batch so, following an approach recommended by Kleynhans et al. [21] and Mitchell et al. [103], we incorporated in every batch of assays a long-established standard reference strain (S06 in our case) whose performance showed no systematic change over the course of the experiment. In our case standardisation entailed dividing the value for each tested individual by the ratio of the median of S06 values for that batch of assays over the overall median of S06 values across all batches.

All statistical analyses were carried out and figures generated using the R software and packages [104]. Relevant helper functions and packages used in the present study are listed in the Additional file 1 and additional scripts and raw data can be found at http://github.com/Angel-Popa/Qfly_abiotic_stress_resistance. Goodness of fit to the distributions assumed in the models was checked and confirmed with standard statistical diagnostics: residuals versus fitted values, distribution of standardised residuals, homogeneity of residual variance and Cook’s distances (Additional file 1 Fig. S5-S16).

Gamma generalized linear models (Gamma-GLM) with log-link transformations of the S06-standardised G2/G3 data for the four resistance variables were used to investigate population differences in those variables. The significance of population differences was then assessed by analysis of deviance. Given that the deviance statistic follows a χ2 distribution, the deviance values obtained could be tested for significance against a χ2 distribution for the appropriate degrees of freedom [105]. Significant effects were followed up by testing for specific effects of coastal vs inland collection sites, their latitudes, five climatic variables for the nearest weather stations (http://www.bom.gov.au) and geographic distances between the sites.

The five climatic variables used in these tests were the 5-year averages between 2013 and 2018 for annual rainfall, annual solar exposure, maximum temperature of the warmest month, minimum temperature of the coldest month, and rainfall for the driest month. Following Hoffmann et al. [19] these were chosen as the least correlated of a larger set of eleven variables (Additional file 1 Fig. S1). Bidirectional elimination was used to select the best regression model for each of the geographically varying resistance measures against the five key climate variables.

Mantel tests were used to test for any effect of geographic distance between collection sites on resistance differences between the populations. Euclidean distances for the distances between populations and between the median stress responses of these populations (Additional file 1 Table S4) were calculated using the base R [104] dist function, and the Mantel tests were carried out using the mantel function in the Vegan R packages, version 2.5–5 [106].

Gamma-GLM analyses with log-link transformations of S06-standardised G2/G3 and G12–15 data were used to investigate the effect of domestication on the four resistance variables. For these analyses a dummy variable was set up in which all the G2/G3 data were coded as ‘Wild’ and all the G12–15 data as ‘Domesticated’. Significant domestication by population interaction effects exposed in this analysis were followed up with post-hoc tests via the EMMEANS procedure from the EMMEANS R package version 1.3.2 [107].

For desiccation and starvation resistance, where data for three intervening generations were available, we also compared the regression slopes of the medians for each colony across generations using the EMTREND function from the EMMEANS package, version 1.3.2 [107].

Wing length data were subjected to analysis of variance (ANOVA) to test for effects of population and domestication, with Tukey post-hoc testing using EMMEANS as above. Associations between changes in resistance variables and wing length were tested by linear regression analysis.

Spearman correlation coefficients among various data sets were calculated using the Hmisc package, version 4.2–0 [108].

Availability of data and materials

The data presented here are available in the Variation in stress resistance in Queensland fruit fly repository https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.9751634.v2

Abbreviations

G10/G11:

Generation 10 or 11

G12-G15:

Generation 12 to 15

G2/G3:

Generation two or three

Gamma-GLM:

Gamma generalized linear models

Qfly:

Queensland fruit fly

S06:

Qfly collected in Sydney in 2006

SIT:

Sterile Insect Technique

SNPs:

Single nucleotide polymorphisms

References

  1. Lee CE. Evolutionary genetics of invasive species. Trends Ecol Evol. 2002;17(8):386–91.

    Article  Google Scholar 

  2. Bowler K, Terblanche JS. Insect thermal tolerance: what is the role of ontogeny, ageing and senescence? Biol Rev. 2008;83:339–55.

    Article  PubMed  Google Scholar 

  3. Gilchrist AS, Dominiak B, Gillespie PS, Sved JA. Variation in population structure across the ecological range of the Queensland fruit fly, Bactrocera tryoni. Aust J Zool. 2006;54:87–95.

    Article  Google Scholar 

  4. Clarke AR, Powell KS, Weldon CW, Taylor PW. The ecology of Bactrocera tryoni (Diptera: Tephritidae): what do we know to assist pest management? Ann Appl Biol. 2011;158:26–54.

    Article  Google Scholar 

  5. NFFS Implementation committee. National Fruit Fly Strategy: Implementation Action Plan. Plant Health Australia. 2010. Available from: http://www.planthealthaustralia.com.au/national-programs/fruit-fly/national-fruit-fly-strategy-implementation-action-plan/.

  6. Dominiak BC, Daniels D. Review of the past and present distribution of Mediterranean fruit fly (Ceratitis capitata Wiedemann) and Queensland fruit fly (Bactrocera tryoni Froggatt) in Australia. Aust J Entomol. 2012;51:104–15.

    Article  Google Scholar 

  7. Hoffmann AA, Weeks AR. Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia. Genetica. 2007;129(2):133–47.

    Article  PubMed  Google Scholar 

  8. Coyne JA, Beecham E. Heritability of two morphological characters within and among natural populations of Drosophila melanogaster. Genetics. 1987;117:727–37.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Imasheva AG, Bubli OA, Lazebny OE. Variation in wing length in Eurasian natural populations of Drosophila melanogaster. Heredity. 1994;72(5):508–14.

    Article  PubMed  Google Scholar 

  10. James AC, Azevedo RBR, Partridge L. Cellular basis and developmental timing in a size cline of Drosophila melanogaster. Genetics. 1995;140(2):659–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. van’T Land J, Van Putten P, Zwaan B, Kamping A, Van Delden W. Latitudinal variation in wild populations of Drosophila melanogaster: Heritabilities and reaction norms. J Evol Biol. 1999;12:222–32.

    Article  Google Scholar 

  12. Loeschcke V, Bundgaard J, Barker JSF. Variation in body size and life history traits in Drosophila aldrichi and D. buzzatii from a latitudinal cline in eastern Australia. Heredity. 2000;85(5):423–33.

    Article  PubMed  Google Scholar 

  13. Karan D, Dahiya N, Munjal AK, Gibert P, Moreteau B, Parkash R, et al. Desiccation and starvation tolerance of adult Drosophila: opposite latitudinal clines in natural populations of three different species. Evolution. 1998;52(3):825–31.

  14. Robinson SJW, Zwaan B, Partridge L. Starvation resistance and adult body composition in a latitudinal cline of Drosophila melanogaster. Evolution. 2000;54(4):1819–24.

    Article  CAS  PubMed  Google Scholar 

  15. Hoffmann AA, Hallas R, Sinclair C, Partridge L, Mitrovski P. Levels of variation in stress resistance in Drosophila among strains, local populations, and geographic regions: patterns for desiccation, starvation, cold resistance, and associated traits. Evolution. 2001;55(8):1621–30.

    Article  CAS  PubMed  Google Scholar 

  16. Mitrovski P, Hoffmann AA. Postponed reproduction as an adaptation to winter conditions in Drosophila melanogaster: evidence for clinal variation under semi-natural conditions. Proc R Soc B Biol Sci. 2001;268(1787):2163–8.

    Article  CAS  Google Scholar 

  17. Hallas R, Schiffer M, Hoffmann AA. Clinal variation in Drosophila serrata for stress resistance and body size. Genet Res. 2002;79(2):141–8.

    Article  PubMed  Google Scholar 

  18. Griffiths JA, Schiffer M, Hoffmann AA. Clinal variation and laboratory adaptation in the rainforest species Drosophila birchii for stress resistance, wing size, wing shape and development time. J Evol Biol. 2005;18(1):213–22.

    Article  CAS  PubMed  Google Scholar 

  19. Hoffmann AA, Anderson A, Hallas R. Opposing clines for high and low temperature resistance in Drosophila melanogaster. Ecol Lett. 2002;5:614–8.

    Article  Google Scholar 

  20. Gilchrist GW, Jeffers LM, West B, Folk DG, Suess J, Huey RB. Clinal patterns of desiccation and starvation resistance in ancestral and invading populations of Drosophila subobscura. Evol Appl. 2008;1(3):513–23.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Kleynhans E, Mitchell KA, Conlong DE, Terblanche JS. Evolved variation in cold tolerance among populations of Eldana saccharina (Lepidoptera: Pyralidae) in South Africa. J Evol Biol. 2014;27:1149–59.

    Article  CAS  PubMed  Google Scholar 

  22. Terblanche JS, Clusella-Trullas S, Deere JA, Van Vuuren BJ, Chown SL. Directional evolution of the slope of the metabolic rate-temperature relationship is correlated with climate. Physiol Biochem Zool. 2009;82(5):495–503.

    Article  PubMed  Google Scholar 

  23. Cheng C, White BJ, Kamdem C, Mockaitis K, Costantini C, Hahn MW, et al. Ecological genomics of Anopheles gambiae along a latitudinal cline: a population-resequencing approach. Genetics. 2012;190(4):1417–32.

  24. Weldon CW, Nyamukondiwa C, Karsten M, Chown SL, Terblanche JS. Geographic variation and plasticity in climate stress resistance among southern African populations of Ceratitis capitata (Wiedemann) (Diptera: Tephritidae). Sci Rep. 2018;8:9849.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Nyamukondiwa C, Kleynhans E, Terblanche JS. Phenotypic plasticity of thermal tolerance contributes to the invasion potential of Mediterranean fruit flies (Ceratitis capitata). Ecol Entomol. 2010;35(5):565–75.

    Article  Google Scholar 

  26. Castañeda LE, Lardies MA, Bozinovic F. Interpopulational variation in recovery time from chill coma along a geographic gradient: a study in the common woodlouse, Porcellio laevis. J Insect Physiol. 2005;51(12):1346–51.

    Article  PubMed  CAS  Google Scholar 

  27. David JR, Gibert P, Pla E, Petavy G, Karan D, Moreteau B. Cold stress tolerance in Drosophila: analysis of chill coma recovery in D. melanogaster. J Therm Biol. 1998;23(5):291–9.

    Article  Google Scholar 

  28. Karl I, Janowitz SA, Fischer K. Altitudinal life-history variation and thermal adaptation in the copper butterfly Lycaena tityrus. Oikos. 2008;117(5):778–88.

    Article  Google Scholar 

  29. Zeilstra I, Fischer K. Cold tolerance in relation to developmental and adult temperature in a butterfly. Physiol Entomol. 2005;30(1):92–5.

    Article  Google Scholar 

  30. Kellermann V, van Heerwaarden B. Terrestrial insects and climate change: adaptive responses in key traits. Physiol Entomol. 2019;44(2):99–115.

    Article  Google Scholar 

  31. Bower CC. The Ecology and Behaviour of the Larvae of the Queensland Fruit Fly, Dacus tryoni (Frogg.) PhD Thesis. University of Sydney; 1975.

    Google Scholar 

  32. Meats A. The bioclimatic potential of the Queensland fruit fly, Dacus tryoni, in Australia. Proc Ecol Soc Aust. 1981;11:151–61.

    Google Scholar 

  33. Meats A. Thermal constraints to successful development of the Queensland fruit fly in regimes of constant and fluctuating temperature. Entomol Exp Appl. 1984;36:55–9.

    Article  Google Scholar 

  34. Meats A, Fitt GP. Survival of repeated frosts by the Queensland fruit fly, Dacus tryoni: experiments in laboratory simulated climates with either step or ramp fluctuations of temperature. Entomol Exp Appl. 1987;45:9–16.

    Article  Google Scholar 

  35. Yonow T, Sutherst RW. The geographical distribution of the Queensland fruit fly, Bactrocera (Dacus) tryoni, in relation to climate. Aust J Agric Res. 1998;49:935–53.

    Article  Google Scholar 

  36. Fletcher BS. The ecology of a natural population of the Queensland fruit fly, Dacus tryoni. IV the immigration and emigration of adults. Aust J Zool. 1973;21:541–65.

    Article  Google Scholar 

  37. Fletcher BS. The ecology of a natural population of the Queensland fruit fly, Dacus tryoni V. the dispersal of adults. Aust J Zool. 1974;22:189–202.

    Article  Google Scholar 

  38. Dominiak BC, Mapson R. Revised distribution of Bactrocera tryoni in eastern Australia and effect on possible incursions of Mediterranean fruit fly: development of Australia’s eastern trading block. J Econ Entomol. 2017;110(6):2459–65.

    Article  PubMed  Google Scholar 

  39. Bateman MA. Adaptations to temperature in geographic races of the Queensland fruit fly, Dacus (strumeta) tryoni. Aust J Zool. 1967;15:1141–61.

    Article  Google Scholar 

  40. Horticulture Innovation Australia. Fact sheet: Sterile Insect Technique Plus (SITPlus) Program. [accessed 2020 May 15]; Available from: https://hortfrontiers.com.au/fruit-fly-fund/.

  41. Rull J, Brunel O, Mendez ME. Mass rearing history negatively affects mating success of male Anastrepha ludens (Diptera: Tephritidae) reared for sterile insect technique programs. J Econ Entomol. 2009;98(5):1510–6.

    Article  Google Scholar 

  42. Schutze MK, Dammalage T, Jessup A, Vreysen MJB, Wornoayporn V, Clarke AR. Effects of laboratory colonization on Bactrocera dorsalis (Diptera, tephritidae) mating behaviour: ‘what a difference a year makes.’. Zookeys. 2015;2015(540):369–83.

    Article  Google Scholar 

  43. Pereira R, Silva N, Quintal C, Abreu R, Andrade J, Dantas L. Sexual performance of mass reared and wild mediterranean fruit flies (Diptera: Tephritidae) from various origins of the Madeira Islands. Florida Entomol. 2007;90(1):10–4.

    Article  Google Scholar 

  44. Hernández E, Toledo J, Artiaga-López T, Flores S. Demographic changes in Anastrepha obliqua (Diptera: Tephritidae) throughout the laboratory colonization process. J Econ Entomol. 2009;102(2):542–51.

    Article  PubMed  Google Scholar 

  45. Hendrichs MA, Wornoayporn V, Katsoyannos B, Hendrichs J. Quality control method to measure predator evasion in wild and mass-reared mediterranean fruit flies (Diptera: Tephritidae). Florida Entomol. 2007;90(1):64–70.

    Article  Google Scholar 

  46. Weldon CW, Taylor PW. Desiccation resistance of adult Queensland fruit flies Bactrocera tryoni decreases with age. Physiol Entomol. 2010;35:385–90.

    Article  Google Scholar 

  47. Raphael KA, Shearman DCA, Gilchrist AS, Sved JA, Morrow JL, Sherwin WB, et al. Australian endemic pest tephritids: genetic, molecular and microbial tools for improved sterile insect technique. BMC Genet. 2014;15:1–13.

  48. Sherwin WB, Frommer M, Sved JA, Raphael KA, Oakeshott JG, Shearman DCA, et al. Tracking invasion and invasiveness in Queensland fruit flies: From classical genetics to ‘omics.’. Curr Zool. 2015;61(3):477–87.

  49. Weldon CW, Prenter J, Taylor PW. Activity patterns of Queensland fruit flies (Bactrocera tryoni) are affected by both mass-rearing and sterilization. Physiol Entomol. 2010;35(2):148–53.

    Article  Google Scholar 

  50. Gilchrist AS, Meats AW. Factors affecting the dispersal of large-scale releases of the Queensland fruit fly, Bactrocera tryoni. J Appl Entomol. 2012;136(4):252–62.

    Article  Google Scholar 

  51. Tejeda MT, Arredondo J, Liedo P, Pérez-Staples D, Ramos-Morales P, Díaz-Fleischer F. Reasons for success: rapid evolution for desiccation resistance and life-history changes in the polyphagous fly Anastrepha ludens. Evolution. 2016;70(11):2583–94.

    Article  PubMed  Google Scholar 

  52. Yu H, Frommer M, Robson MK, Meats AW, Shearman DC, Sved J. a. Microsatellite analysis of the Queensland fruit fly Bactrocera tryoni (Diptera: Tephritidae) indicates spatial structuring: implications for population control. Bull Entomol Res. 2001;91(2):139–47.

    CAS  PubMed  Google Scholar 

  53. Cameron EC, Sved JA, Gilchrist AS. Pest fruit fly (Diptera: Tephritidae) in northwestern Australia: one species or two? Bull Entomol Res. 2010;100(02):197–206.

    Article  CAS  PubMed  Google Scholar 

  54. Drew RAI. The tropical fruit flies (Diptera: Tephritidae: Dacinae) of the Australasian and Oceanian regions. Mem Queensl Museum. 1989;26:1–521.

    Google Scholar 

  55. Vargas RI, Piñero JC, Leblanc L. An overview of pest species of Bactrocera fruit flies (Diptera: Tephritidae) and the integration of biopesticides with other biological approaches for their management with a focus on the Pacific region. Insects. 2015;6(2):297–318.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Drew RAI, Lambert DM. On the specific status of Dacus (Bactrocera) aquilonis and D. (Bactrocera) tryoni (Diptera: Tephritidae). Ann Entomol Soc Am. 1986;79(6):870–8.

    Article  Google Scholar 

  57. Duyck PF, David P, Quilici S. Climatic niche partitioning following successive invasions by fruit flies in La Réunion. J Anim Ecol. 2006;75:518–26.

    Article  PubMed  Google Scholar 

  58. Terblanche JS, Kleynhans E. Phenotypic plasticity of desiccation resistance in Glossina puparia: are there ecotype constraints on acclimation responses? J Evol Biol. 2009;22(8):1636–48.

    Article  CAS  PubMed  Google Scholar 

  59. Weldon CW, Yap S, Taylor PW. Desiccation resistance of wild and mass-reared Bactrocera tryoni (Diptera: Tephritidae). Bull Entomol Res. 2013;103(6):690–9.

    Article  CAS  PubMed  Google Scholar 

  60. Matos M, Rose MR, Rocha Pité MT, Rego C, Avelar T. Adaptation to the laboratory environment in Drosophila subobscura. J Evol Biol. 2000;13(1):9–19.

    Article  Google Scholar 

  61. Simões P, Santos J, Fragata I, Mueller LD, Rose MR, Matos M. How repeatable is adaptive evolution? The role of geographical origin and founder effects in laboratory adaptation. Evolution. 2008;62(8):1817–29.

    Article  PubMed  Google Scholar 

  62. Hoffmann AA, Hallas R, Sinclair C, Partridge L. Rapid loss of stress resistance in Drosophila melanogaster under adaptation to laboratory culture. Evolution. 2001;55(2):436–8.

    Article  CAS  PubMed  Google Scholar 

  63. Kuba H, Koyama J. Mating behavior of the melon fly, Dacus cucurbitae Coquilett (Diptera: Tephritidae): comparative studies of one wild and two laboratory strains. Appl Entomol Zool. 1982;17(4):559–69.

    Article  Google Scholar 

  64. Ekesi S, Nderitu PW, Chang CL. Adaptation to and small-scale rearing of invasive fruit fly Bactrocera invadens (Diptera: Tephritidae) on artificial diet. Ann Entomol Soc Am. 2007;100(4):562–7.

    Article  Google Scholar 

  65. Hoffmann AA, Ross PA. Rates and patterns of laboratory adaptation in (mostly) insects. J Econ Entomol. 2018;111(2):501–9.

    Article  PubMed  Google Scholar 

  66. Meats A, Holmes HM, Kelly GL. Laboratory adaptation of Bactrocera tryoni (Diptera: Tephritidae) decreases mating age and increases protein consumption and number of eggs produced per milligram of protein. Bull Entomol Res. 2004;94:517–24.

    Article  CAS  PubMed  Google Scholar 

  67. Lynch KE, White TE, Kemp DJ. The effect of captive breeding upon adult thermal preference in the Queensland fruit fly (Bactrocera tryoni). J Therm Biol. 2018;78:290–7.

    Article  PubMed  Google Scholar 

  68. Inskeep JR, Weldon CW, Taylor PW, Rempoulakis P. Differences in the microhabitat distribution and behaviour of sterile mass-reared and wild wild Queensland fruit flies, Bactrocera tryoni (Froggatt). Poster Present 10th Int Symp Fruit Flies Econ Importance. Tapachula; 2018.:23–7.

  69. Pérez J, Park SJ, Taylor PW. Domestication modifies the volatile emissions produced by male Queensland fruit flies during sexual advertisement. Sci Rep. 2018;8:16503.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Gilchrist AS, Cameron EC, Sved JA, Meats AW. Genetic consequences of domestication and mass rearing of pest fruit fly Bactrocera tryoni (Diptera: Tephritidae). J Econ Entomol. 2012;105(3):1051–6.

    Article  CAS  PubMed  Google Scholar 

  71. Simões P, Fragata I, Seabra SG, Faria GS, Santos MA, Rose MR, et al. Predictable phenotypic, but not karyotypic, evolution of populations with contrasting initial history. Sci Rep. 2017;7:913.

  72. Simões P, Pascual M, Coelho MM, Matos M. Divergent evolution of molecular markers during laboratory adaptation in Drosophila subobscura. Genetica. 2010;138(9):999–1009.

    Article  PubMed  CAS  Google Scholar 

  73. Fragata I, Simões P, Lopes-Cunha M, Lima M, Kellen B, Bárbaro M, et al. Laboratory selection quickly erases historical differentiation. PLoS One. 2014;9(5):e96227.

  74. Santos M, Fragata I, Santos J, Simões P, Marques A, Lima M, et al. Playing Darwin. Part B. 20 years of domestication in Drosophila subobscura. Theory Biosci. 2010;129(2–3):97–102.

  75. Tejeda MT, Arredondo-Gordillo J, Orozco-Dávila D, Quintero-Fong L, Díaz-Fleischer F. Directional selection to improve the sterile insect technique: survival and sexual performance of desiccation resistant Anastrepha ludens strains. Evol Appl. 2017;10:1020–30.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Gilchrist AS, Shearman DCA, Frommer M, Raphael KA, Deshpande NP, Wilkins MR, et al. The draft genome of the pest tephritid fruit fly Bactrocera tryoni: Resources for the genomic analysis of hybridising species. BMC Genomics. 2014;15:1153.

  77. Telonis-Scott M, Sgró CM, Hoffmann AA, Griffin PC. Cross-study comparison reveals common genomic, network, and functional signatures of desiccation resistance in Drosophila melanogaster. Mol Biol Evol. 2016;33(4):1053–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Sgrò CM, Terblanche JS, Hoffmann AA. What can plasticity contribute to insect responses to climate change? Annu Rev Entomol. 2016;61(1):433–51.

    Article  PubMed  CAS  Google Scholar 

  79. Valladares F, Matesanz S, Guilhaumon F, Araujo MB, Balaguer L, Benito-Garz NM, et al. The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change. Ecol Lett. 2014;17(11):1351–64.

  80. Chown SL, Slabber S, McGeoch MA, Janion C, Leinaas HP. Phenotypic plasticity mediates climate change responses among invasive and indigenous arthropods. Proc R Soc B Biol Sci. 2007;274(1625):2531–7.

    Article  Google Scholar 

  81. Manenti T, Loeschcke V, Moghadam NN, SØrensen JG. Phenotypic plasticity is not affected by experimental evolution in constant, predictable or unpredictable fluctuating thermal environments. J Evol Biol. 2015;28(11):2078–87.

    Article  CAS  PubMed  Google Scholar 

  82. Wang Y, Yu H, Raphael K, Gilchrist AS. Genetic delineation of sibling species of the pest fruit fly Bactrocera (Diptera: Tephritidae) using microsatellites. Bull Entomol Res. 2003;93(4):351–60.

    Article  CAS  PubMed  Google Scholar 

  83. Horticulture Innovation Australia. Establishment of the Qfly SITplus facility in southern Australia (HG14035). 2018. Available from: https://www.horticulture.com.au/growers/help-your-business-grow/research-reports-publications-fact-sheets-and-more/hg14035/.

  84. Horticulture Innovation Australia. SITplus: Port Augusta Qfly SIT factory pilot operation (FF15000). 2019. Available from: https://www.horticulture.com.au/growers/help-your-business-grow/research-reports-publications-fact-sheets-and-more/ff15000/.

  85. Beck HE, Zimmermann NE, McVicar TR, Vergopolan N, Berg A, Wood EF. Data Descriptor: Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci Data. 2018;5:180214.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Hulthen AD, Clarke AR. The influence of soil type and moisture on pupal survival of Bactrocera tryoni (Froggatt) (Diptera: Tephritidae). Aust J Entomol. 2006;45:16–9.

    Article  Google Scholar 

  87. Tychsen PH. Mating behaviour of the Queensland fruit fly, Dacus tryoni (Diptera: Tephritidae), in field cages. Aust J Entomol. 1977;16:459–65.

    Article  Google Scholar 

  88. Moadeli T, Taylor PW, Ponton F. High productivity gel diets for rearing of Queensland fruit fly, Bactrocera tryoni. J Pest Sci. 2017;90:507–20.

    Article  Google Scholar 

  89. Rush B, Sandver S, Bruer J, Roche R, Wells M, Giebultowicz J. Mating increases starvation resistance and decreases oxidative stress resistance in Drosophila melanogaster females. Aging Cell. 2007;6(5):723–6.

    Article  CAS  PubMed  Google Scholar 

  90. Yap S, Fanson BG, Taylor PW. Mating reverses actuarial aging in female Queensland fruit flies. PLoS One. 2015;10(7):e0132486.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  91. Gilchrist AS, Meats AW. An evaluation of outcrossing to improve mass-reared strains of the Queensland fruit fly Bactrocera tryoni. Int J Trop Insect Sci. 2014;34(1):35–44.

    Article  Google Scholar 

  92. Weldon CW, Pérez-Staples D, Taylor PW. Feeding on yeast hydrolysate enhances attraction to cue-lure in Queensland fruit flies, Bactrocera tryoni. Entomol Exp Appl. 2008;129(2):200–9.

    Article  Google Scholar 

  93. Pérez-Staples D, Weldon CW, Taylor PW. Sex differences in developmental response to yeast hydrolysate supplements in adult Queensland fruit fly. Entomol Exp Appl. 2011;141(2):103–13.

    Article  Google Scholar 

  94. Akter H, Adnan S, Morelli R, Rempoulakis P, Taylor PW. Suppression of cuelure attraction in male Queensland fruit flies provided raspberry ketone supplements as immature adults. PLoS One. 2017;12(8):e0184086.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  95. Colinet H, Lee SF, Hoffmann AA. Functional characterization of the frost gene in Drosophila melanogaster: importance for recovery from chill coma. PLoS One. 2010;5(6):e10925.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  96. Sgrò CM, Overgaard J, Kristensen TN, Mitchell KA, Cockerell FE, Hoffmann AA. A comprehensive assessment of geographic variation in heat tolerance and hardening capacity in populations of Drosophila melanogaster from eastern Australia. J Evol Biol. 2010;23:2484–93.

    Article  PubMed  Google Scholar 

  97. Gilbert P, Moreteau B, Petavy G, Karan D, David JR. Chill-coma tolerance, a major climatic adaptation among Drosophila species. Evolution. 2001;55(5):1063–8.

    Article  Google Scholar 

  98. Hoffmann AA, Sørensen JG, Loeschcke V. Adaptation of Drosophila to temperature extremes: bringing together quantitative and molecular approaches. J Therm Biol. 2003;28(3):175–216.

    Article  Google Scholar 

  99. Weldon CW, Boardman L, Marlin D, Terblanche JS. Physiological mechanisms of dehydration tolerance contribute to the invasion potential of Ceratitis capitata (Wiedemann) (Diptera: Tephritidae) relative to its less widely distributed congeners. Front Zool. 2016;13:15.

    Article  PubMed  PubMed Central  Google Scholar 

  100. Tejeda MT, Arredondo J, Pérez-Staples D, Ramos-Morales P, Liedo P, Díaz-Fleischer F. Effects of size, sex and teneral resources on the resistance to hydric stress in the tephritid fruit fly Anastrepha ludens. J Insect Physiol. 2014;70:73–80.

    Article  CAS  PubMed  Google Scholar 

  101. Reynolds OL, Orchard BA. Effect of adult chill treatments on recovery, longevity and flight ability of Queensland fruit fly, Bactrocera tryoni (Froggatt) (Diptera: Tephritidae). Bull Entomol Res. 2011;101(1):63–71.

    Article  CAS  PubMed  Google Scholar 

  102. Prabhu V, Pérez-Staples D, Taylor PW. Protein: carbohydrate ratios promoting sexual activity and longevity of male Queensland fruit flies. J Appl Entomol. 2008;132(7):575–82.

    Article  CAS  Google Scholar 

  103. Mitchell KA, Boardman L, Clusella-Trullas S, Terblanche JS. Effects of nutrient and water restriction on thermal tolerance: a test of mechanisms and hypotheses. Comp Biochem Physiol -Part A. 2017;212:15–23.

    Article  CAS  Google Scholar 

  104. R Core Team. R: A language and environment for statistical computing. https://www.r-project.org/. Vienna, Austria; 2018. Available from: https://www.r-project.org/.

  105. Lane PW. Generalized linear models in soil science. Eur J Soil Sci. 2002;53:241–51.

    Article  Google Scholar 

  106. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community ecology package. R package version 2.5–3. Available from: https://cran.r-project.org/package=vegan.

  107. Lenth R. emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.4. Available from: https://cran.r-project.org/package=emmeans.

  108. Harrell Jr. FE. Hmisc: Harrell Miscellaneous. R package version 4.2–0. Available from: https://cran.r-project.org/package=Hmisc.

  109. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag. New York; 2016. VIII, 213. Available from: https://ggplot2.tidyverse.org.

Download references

Acknowledgements

We thank Dr. Stuart Gilchrist for providing the S06 control strain, and Freya Robinson, Valentina Colombo, Jenna Ryu and Tim Hogarty for help with the maintenance of stocks and bioassays. We also thank the following colleagues for the assistance in Qfly sampling: Chris Chilcott (CSIRO), Magen Pettit (CSIRO), Brian Thistleton (NT DPIF) and Mike Neal (NT DPIF) in Darwin, Geoff Miers, Richard Tucker (NT DPIF) and Chris Pavey (CSIRO) in Alice Springs, the Cape Tribulation Farm and Daintree Ice Cream Co in Cape Tribulation, Stefano De Faveri (QDAFF) and Jodie Cheesman (QDAFF) in Mareeba, David Warmington (Cairns Botanical Garden), Sybilla Oczkowicz (QDAFF), Peter Leach (QDAFF), Yan Diczbalis (QDAFF) and Peter Salleras (Fruit Forest Farm) in Cairns and Utchee Creek, Brendan Missenden (QDAFF) and Jane Royer (QDAFF) in Brisbane, Simone Heimoana (CSIRO) and Sharon Downes (CSIRO) in Narrabri, Andrew Creek (NSW-DPI) and Ottorino Cappello (NSW-DPI) in Griffith, Michael Calkovics (Canberra Organic Growers Society), Carol Hartley (CSIRO), Nigel England (CSIRO) and Tom Walsh (CSIRO) in Canberra, Peter Gow (SAGE community garden Moruya) in Batemans Bay; and the Bega Valley Seed Savers in the Bega Valley region.

About this supplement

This article has been published as part of BMC Genetics Volume 21 Supplement 2, 2020: Comparing rearing efficiency and competitiveness of sterile male strains produced by genetic, transgenic or symbiont-based technologies. The full contents of the supplement are available online at https://0-bmcgenet-biomedcentral-com.brum.beds.ac.uk/articles/supplements/volume-21-supplement-2.

Funding

This research was conducted as part of the SITplus collaborative fruit fly programmes. It is a component of the project Raising Q-fly Sterile Insect Technique to World Standard (HG14033) funded by the Hort Frontiers Fruit Fly Fund, part of the Hort Frontiers strategic partnership initiative developed by Hort Innovation, with co-investment from Macquarie University and contributions from the Australian Government. We also acknowledge CSIRO for its generous support of this project. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. This study benefitted from discussions at International Atomic Energy Agency funded meetings for the Coordinated Research Project “Comparing Rearing Efficiency and Competitiveness of Sterile Male Strains Produced by Genetic, Transgenic or Symbiotic-based Technologies”. Publication costs are funded by the Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture, IAEA (CRP No.: D4.20.16) Vienna, Austria.

Author information

Authors and Affiliations

Authors

Contributions

ADPB: study design, data collection, analysis and interpretation, co-wrote the manuscript; SFL: study design, data collection, analysis and interpretation, co-wrote the manuscript; HLY: data collection and analysis, co-wrote the manuscript; SSP: data collection; MS: data collection; RGM: data collection; CCV: study design and interpretation; ORE: study design and interpretation; PWT: study design and interpretation; JGO: study design, data analysis and interpretation, co-wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ángel-David Popa-Báez.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Climatic variables from the Qfly collection sites. Table S2. Individual populations for which the wild (G2/G3) and domesticated (G10–15) bioassays results differed significantly. Contrast is calculated for the estimated mean response variable for each population by looking at the differences of the domesticated over the wild populations. The estimated mean of the contrast is calculated on the log-transformed data for the response variables. Table S3. Methodological differences between the standard desiccation resistance and that used for the resampled 2017/2018 collection. Table S4. Euclidean distance between site’s geographical coordinates. Fig. S1. Correlation among 11 climatic variables. Correlation values are presented together with asterisks indicating significance values for each correlation. ‘*’ P < 0.05; ‘**’ P < 0.01; ‘**’ P < 0.001. mean.max = Annual maximum temperature; mean.min = Annual minimum temperature; mean.rain = Annual rainfall; mean.solar = Annual solar exposure; annual.temp = Annual temperature; max.high.temp = Maximum temperature of the warmest month; min.high.temp = Minimum temperature of the warmest month; min.low.temp = Minimum temperature of the coldest month; max.low.temp = Maximum temperature of the coldest month; ppt.dry.month = Precipitation of the driest month; ppt.wet.month = Precipitation of the wettest month. Fig. S2. Egging device used in present study. Fig. S3. Results of the pilot experiment on heat knock down recovery time. Fig. S4. Cold resistance apparatus used in present study. Fig. S5. Diagnostic plots Gamma-GLM heat resistance in wild populations of the Queensland fruit fly. Fig. S6. Diagnostic plots Gamma-GLM heat resistance in domesticated populations of the Queensland fruit fly. Fig. S7. Diagnostic plots Gamma-GLM heat resistance change during domestication. Fig. S8. Diagnostic plots Gamma-GLM cold resistance in wild populations of the Queensland fruit fly. Fig. S9. Diagnostic plots Gamma-GLM cold resistance in domesticated populations of the Queensland fruit fly. Fig. S10. Diagnostic plots Gamma-GLM cold resistance change during domestication. Fig. S11. Diagnostic plots Gamma-GLM desiccation resistance in wild Qfly populations. Fig. S12. Diagnostic plots Gamma-GLM desiccation resistance in domesticated populations of the Queensland fruit fly. Fig. S13. Diagnostic plots Gamma-GLM desiccation resistance change during domestication. Fig. S14. Diagnostic plots Gamma-GLM starvation resistance in wild Qfly populations. Fig. S15. Diagnostic plots Gamma-GLM starvation resistance in domesticated populations of the Queensland fruit fly. Fig. S16. Diagnostic plots Gamma-GLM starvation resistance change during domestication.

Rights and permissions

Open Access This is an open access article distributed under the terms of the Creative Commons Attribution IGO License (https://creativecommons.org/licenses/by/3.0/igo/) which permits unrestricted use, distribution, and reproduction in any medium, provided appropriate credit to the original author(s) and the source is given.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Popa-Báez, ÁD., Lee, S.F., Yeap, H.L. et al. Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication. BMC Genet 21 (Suppl 2), 135 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12863-020-00935-2

Download citation

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12863-020-00935-2

Keywords