The effects of adaptation to urea on feeding rates and growth in Drosophila larvae

Abstract A collection of forty populations were used to study the phenotypic adaptation of Drosophila melanogaster larvae to urea‐laced food. A long‐term goal of this research is to map genes responsible for these phenotypes. This mapping requires large numbers of populations. Thus, we studied fifteen populations subjected to direct selection for urea tolerance and five controls. In addition, we studied another twenty populations which had not been exposed to urea but were subjected to stress or demographic selection. In this study, we describe the differentiation in these population for six phenotypes: (1) larval feeding rates, (2) larval viability in urea‐laced food, (3) larval development time in urea‐laced food, (4) adult starvation times, (5) adult desiccation times, and (6) larval growth rates. No significant differences were observed for desiccation resistance. The demographically/stress‐selected populations had longer times to starvation than urea‐selected populations. The urea‐adapted populations showed elevated survival and reduced development time in urea‐laced food relative to the control and nonadapted populations. The urea‐adapted populations also showed reduced larval feeding rates relative to controls. We show that there is a strong linear relationship between feeding rates and growth rates at the same larval ages feeding rates were measured. This suggests that feeding rates are correlated with food intake and growth. This relationship between larval feeding rates, food consumption, and efficiency has been postulated to involve important trade‐offs that govern larval evolution in stressful environments. Our results support the idea that energy allocation is a central organizing theme in adaptive evolution.


| 9517
BITNER ET al. Recent research has shown inferring gene-phenotype relationships requires a large number of independent populations to achieve adequate statistical power (Baldwin-Brown et al., 2014;Mueller et al., 2018). Most experimental evolution studies utilize just 6-10 independent populations. In this study, we have assembled a collection of forty independent populations to measure important phenotypes. We assessed the phenotypic differentiation that occurs in populations adapted to urea. Our results allowed us to determine the extent to which populations that have been selected for stress resistance or age-at-reproduction show correlated changes in the same phenotypes. Future work will collect genomic data from these forty populations to make gene-phenotype inferences.
The 40 study populations can be broken into three different categories. (a) Fifteen populations studied have been directly selected for tolerance to high levels of urea in their larval food. An additional five populations served as their direct controls. This group of twenty populations will be referred to as the urea selection populations.
(b) Another 10 populations were derived from populations selected for late reproduction, five were subjected to desiccation stress and five served as controls under mild starvation stress. This collection of 10 populations will be referred to as the stress-selected populations. (c) Finally, 10 populations were selected for late reproduction. These populations will be referred to as the demographically selected populations.
Selection for reproduction at later ages results in increased longevity as well as changes in a number of correlated traits (Partridge & Fowler, 1992;Service et al., 1985). Service et al. (1985) showed that populations selected for reproduction at later ages also showed increases in adult starvation resistance, desiccation resistance, and ethanol tolerance. Further evidence of these correlations came from populations directly selected for desiccation resistance and starvation resistance that also exhibited correlated increases in longevity (Rose et al., 1992). The stress-selected populations may show high levels of starvation resistance and desiccation resistance due to their derivation from populations selected for postponed selection or their direct history of desiccation selection. The demographically selected populations would be expected to show stress resistance due to selection on these correlated traits. It is unclear whether larval stress selection will affect adult stress traits. We test that possibility here with the urea-adapted populations.
Evolutionary studies using the model organism Drosophila melanogaster have often focused on the evolution of adult traits like behaviors (Mery & Kawecki, 2003), heat tolerance (Gilchrist & Huey, 1999), desiccation tolerance (Gibbs et al., 1997), starvation tolerance (Chippindale et al., 1996), and age-specific survival and fecundity (Rose, 1984). A number of studies have examined the evolution of larval phenotypes like competitive ability (Mueller, 1988a), urea and ammonia tolerance (Borash et al., 2000), parasite tolerance (Fellowes et al., 1999), and low nutrition (Kolss et al., 2009). A common adaptation among the studies of larval adaptation is the evolution of larval feeding rates. Thus, crowded larval conditions result in the increase in larval feeding rates (Joshi & Mueller, 1988; although see also Nagarajan et al., 2016), while stresses like food laced with ammonia or urea (Borash et al., 2000) and larval parasites (Fellowes et al., 1999) result in decreased feeding rates.
Feeding rates have been shown to be correlated with larval competitive ability (Burnet et al., 1977). However, Drosophila larvae can evolve increased competitive ability without changing their feeding rate in certain environments (Nagarajan et al., 2016). There is also evidence that the rate at which food passes through the larval alimentary tract is proportional to the feeding rate (Burnet et al., 1977).
Furthermore, larvae with very high feeding rates are less efficient than slower feeding larvae-that is they require more food to successfully pupate (Joshi & Mueller, 1996;Mueller, 1990). These findings suggested that energy trade-offs may be driving the evolution of feeding rates (Mueller & Barter, 2015). An important foundation for theories of life-history evolution is the concept of trade-offs (Cody, 1966;Dudley & Schmitt, 1996;Gadgil & Bossert, 1970;Sinervo, 1990;Stearns, 1992). Cody (1966) originally developed the idea that an organism limited energy budget determined the currency of these trade-offs. This idea was sufficiently powerful that entire fields, like optimal foraging, are based on the idea of maximizing energy intake (Charnov, 1976;MacArthur & Pianka, 1966).
Adaptation of Drosophila larvae to stress offers rich experimental system to illustrate this idea. Mueller and Barter (2015) suggested that in toxic food environments, a reduction in feeding rates would reduce the intake of toxins while increasing the efficiency of energy extraction from the food consumed. This would allow the larvae to route the extra energy they get from enhanced efficiency to detoxification. This theory hinges on the implied relationship between feeding rates and food consumption. This relationship has been called into question by Kaun et al. (2007) who found a weak positive but nonsignificant correlation between feeding rates and food consumption (see also Brown et al., 2019). In normal, nontoxic food, slowly feeding larvae should be below their optimum feeding rate, as measured by energy gained per unit time, and thus would be expected to grow more slowly than faster feeding larvae. We tested that prediction by measuring feeding rates and larval growth rates of all forty populations.

| Populations
Forty populations were used in this study from three different categories. The urea selection populations consisted of UX, RUX (reverse-selected UX), and UTB selected for urea resistance and the control AUC (urea control). The stress-selected populations consisted of the TDO (desiccation selected) populations which were originally selected for desiccation resistance in 1988 and the control TSO (desiccation controls) populations (Rose et al., 1992). The TSO populations experienced mild starvation selection. In 2005, the TSO and TDO populations were relaxed from selection and maintained on a 21-day culture regime. At the time of these experiments, the TSO and TDO populations had been removed from selection for about 243 generations (Figure 1). The demographically selected populations consisted of the five CO and five nCO (new CO) populations.
All populations are derived from an ancestral "IV" population collected from South Amherst, MA in 1975 by Philip Ives (Rose, 1984, Figure 1). After four and a half years of laboratory culture, the B 1-5 populations (baseline, 14-day generation cycle) and O 1-5 populations (70-day generation cycle) were derived from the single IV population in 1980 (Rose, 1984;Rose et al., 2004). The CO 1-5 populations (28-day generation cycle) were derived from the O 1-5 populations in 1989 (Rose et al., 1992). The nCO populations are a more recent creation of the CO selection regime.
All populations were maintained on a banana molasses food (Rose, 1984) at 25°C (24 hr light), at an uncontrolled humidity, and having a generation time of approximately 3-4 weeks. All selection regimes are fivefold replicated, uncrowded as larva (60-80 eggs per 8-dram vial), with emergent adults kept at a low density of approximately 50-60 flies per 8-dram vial, and transferred to a cage environment after 14 days of development from egg. Fresh food was provided in these cages about every other day for approximately 1 week. Effective population sizes in our experimental populations have an estimated range from 700-1,000 (Mueller et al., 2013) and are maintained at large breeding population sizes (≥1,000), with discrete generations.
The UX and the UTB populations were subjected to selection for increased larval tolerance to the presence of toxic levels of urea in the food for 382 and 98 generations respectively. Both selection regimes ultimately trace their ancestry back to the B populations ( Figure 1).
The level of urea was increased every few generations, when it was observed that a great proportion of larvae were surviving to adulthood. The final urea levels were 0.27 M. The urea-tolerant (UX) and unselected controls (AUC) were derived in the fall of 1996. Both sets of populations were derived from the five UU populations, which had a 3-week generation time, and were reared at low larval and adult densities ( Figure 1). The UU populations were derived in 1990, from the Rose B populations (Chippindale et al., 1994(Chippindale et al., , 1996Rose, 1984).
The RUX is a reverse-selected line of the UX. The RUX populations have been in a control environment for 100 generations (Figure 1).
The UTB line was created in October 2013 ( Figure 1). The AUC, RUX, UX, and UTB lines are all 21-day cycle flies. The TSO, TDO, nCO, and CO selection regimes did not include urea and were developed independently of the urea-selected lines and their control.
For every assay, eggs were collected from the eight population types and passed through two generations of common, standard conditions-low larval density, 1,000 adult density, discrete generation times, 24-hr light cycle, and regular banana molasses food (Bitner et al., 2020).

| Feeding rates
Feeding rates were collected on all 40 populations. Eggs were collected from adults who had undergone a two-generation standardization procedure. To measure the feeding rate, individual larvae around 48 hr old were gently moved onto a 3% agar coated with a 10% live yeast suspension. The larva was given 60 s to adjust to the new surroundings, and the sclerite retractions were recorded for 60 s. A total of twenty-five larvae per population were individually measured this way. This procedure for measuring sclerite rates is similar to Sewell et al. (1975) and described in Joshi and Mueller (1988) F I G U R E 1 The phylogeny of laboratory-selected populations used in this study. A selection regime is noted by a name like nCO or CO and each of these consists of five independent replicates. The numbers indicate the number of generations each selection regime has passed through at the time of the experiment. The feeding rates were measured in two blocks. The TSO populations were measured in each block and thus served as a reference.
Feeding rates were standardized in each block as the difference from that block's TSO mean feeding rate. With this standardization, we let y ijk be the feeding rate difference for selection regime i (i = 1,…,7, since TSO is now used as a reference), population j (j = 1,…,35), and individual k (k = 1,…,25). The linear mixed effects model for y ijk is, where δ i = 0, if i = 1, and 1 otherwise, B j and ε ijk are independent normally distributed random variables with zero mean and variances 2 B and 2 . Random variation due to drift, founder effects, and handling of individual populations is reflected in B j , random variation between individual feeding rate measurements is measured by ε ijk . The model was analyzed with the linear mixed effects R function lme (R Core Team, 2017).
Pairwise differences were evaluated after adjusting for multiple comparisons using Tukey's method implemented by the R function lsmeans.
The UX, RUX, and UTB populations all have a current or past history of being raised as larvae in food with urea. Previous research has demonstrated that feeding rates decline as populations adapt to urea (Borash et al., 2000). We were interested in testing the feeding rates of these three selection regimes as a group, and our initial statistical analysis showed no significant differences between UX, RUX, and UTB feeding rates. Thus, we pooled the feeding rate measurements for all 15 UX, RUX, and UTB populations into one urea population and repeated the analysis as described above comparing the pooled urea-selected populations to the AUC control.

| Starvation and desiccation
Individual female flies were collected from the eight selection regimes, a total of 30 flies per replicate per assay (150 flies per selection regime for one assay). Each fly was placed into a straw with pipette tips on both ends. The straws were wide enough for the flies to be able to move from one end to the other. The flies selected for the starvation assay were placed into straws with 3% agar while the flies selected for the desiccation assay were placed into straws with desiccant separated by cheese cloth. The cheese cloth functioned to prevent direct contact between the flies and desiccant. Flies subjected to desiccation were checked hourly, and flies undergoing starvation were checked every 4 hr.
Time to death by starvation and desiccation was collected on a total of eight selection regimes. We can let the desiccation (or starva-

| Viability versus food type
The viability experiment started with 50 first instar larvae. However, due to the size of this experiment and technical difficulty of counting out exactly 50 larvae, the actual number of larvae put in each vial is more properly thought of as a random variable which could be both higher or lower than 50. Thus, in this analysis, we have analyzed the total number of larvae that survived to become adults under the assumption that the mean number of input larvae was the same in all treatments. We also developed a statistical analysis of the relative level of adaptation to urea for each selection regime. Let y ureajkm be the number of survivors in urea food, selection regime j (TSO, AUC, CO, nCO, RUX, UX, and UTB), population k (k = 1,…,35), and replicate m (m = 1,…,10). Then we analyzed the differences, Δ jkm = y controljk − y ureajkm where y controljk is the mean survival in control food in selection regime j and population k.

| Developmental time versus food type
This assay looked at differences in the development time of the larvae. The measurements were made from the time first instar larvae (1) were placed in vials. For each population, 10 vials were set up with 50 freshly hatched larvae each in a banana molasses environment.
Another 10 vials were set up with urea food and 50 freshly hatched larvae. Eclosing adults were collected every 6 hr, separated by sex, and recorded.
The methods for the analysis of the development time experiment were the same as the viability analysis except there was one additional fixed effect, sex. The first analysis was on the difference between the urea development time and the control development time. Here, the mean control development time for each selection/ sex/population/rep sample was calculated and then subtracted from the corresponding development time in urea.
To test for differences due to selection among urea populations, we analyzed blocks (i) and (iii) using UX-urea as the standard. Thus, in block (i), we computed the mean UX-urea development time, and this mean was subtracted from the other control and urea observations. The same was done for block (iii) using the UX-urea mean from block (iii). To test the stress and demographically selected populations we analyzed blocks (ii) and (iii) using UTB as the standard. Thus, in block (ii) we computed the mean UTB-urea development time, and this mean was subtracted from the other control and urea observations. The same was done for block (iii) using the UTB-urea mean from block (iii).

| Larval growth rate
Pairs of populations, matched by the replicate number, were assayed together in a randomized block design and the assay followed the same procedures as mentioned in Santos et al. (1997). With a two-generation lead in for each population, 45 newly hatched first instar larvae were collected with a fine paint brush and placed onto non-nutritive agar petri dishes with 3 ml of yeast paste (188 g of yeast in 500 ml of DI water) and placed randomly into a 25°C incubator with 24-hr lighting.
For each population, a separate petri dish of 45 larvae was created for every time sample in this experiment. There were 13 different "hour numbers" at which larvae were sampled: 24,30,36,42,48,54,60,66,72,78,84,90, and 105 hr after the larvae were added to the petri dish.
At the designated hour, 30 larvae were washed with DI water and then allowed to air dry. The wet weight, to the nearest 0.01 mg, of the 30 larvae was recorded before they were then placed into an 80 degree C drying oven, after which their dry weights were recorded. Empirically, larval growth follows a logistic trajectory (Santos et al., 1997). We used a three-parameter logistic function to model basic growth dynamics. Under this model, the size of individual larvae after t-hours of growth is given by, In model (2) the asymptotic maximum size is equal to φ 1 , and the time to reach half the maximum size is φ 2 . As t → ∞, exp 2 − t ∕ 3 → 0, and the size approaches the asymptotic value φ 1 . The speed of this approach is governed by φ 3 . The smaller φ 3 the faster the approach to the asymptote. With this model, we let y ijkt be the average size of a larva from selection regime i (i = 1 (AUC), 2 (CO), 3 (nCO), 4 (RUX), 5 (UTB), 6 (UX), 7 (TSO), and 8(TDO)), population j (j = 1,…,40), and block k at time t. Random variation arises due to both population effects (random genetic drift), block effects, and individual variation. Consequently, the size of larvae from selection regime i and population j at time t is y ijkt = f(φ ijk ,t) + ε ijkt , and, where δ i = 0, if i = 1 and 1 otherwise. The within-population variation, ε, is assumed to be normally distributed with a zero mean. This variation increases as the larvae get larger so we assumed that Var(ε) = 2 |t| 2Δ where ∆ is estimated from the data. Population variation, b, and block variation, c, is assumed to only affect parameter φ 2 . We tested models with population variation in the other parameters, and the model with variation in φ 2 was chosen due to having the lowest Akaike and Bayesian information criterion (Pinheiro & Bates, 2000, chapter 8). The population variation is assumed independent of the within-population variation and has a normal distribution with zero mean and variance, 2 b . (3)  When displaying the predicted larval size from Equations (2 and 3), we also calculated 95% confidence intervals. With eight different selection regimes, we have 24 maximum likelihood parameters estimates and their covariance matrix estimates, ̂ = ̂ 1 ,̂ 2 ,̂ 3 ,̂ 12 ,̂ 22 ,̂ 32 , …,̂ 38 and Σ . These were assumed to have a t-distribution. From these distributions, we drew samples of the parameter vectors, ̃ k , (k = 1,…,m). For each sampled parameter vector, we made size predictions for each selection regime for ages, 42, 48, 54, and 60 hr. At a specific age, let the kth (out of m) prediction for selection regime i be ỹ ki . From these m predictions, we generated order statistics, Δ s (ỹ ki ), where Δ 1 (ỹ ki ) is the smallest predicted value at t and Δ m (ỹ ki ) is the largest. From the order statistics, we then used Δ l (ỹ ki ) as the lower confidence limit and Δ u (ỹ ki ) as the upper confidence limit. In our simulations, we set m = 5,000. Therefore, a 95% confidence interval corresponds to Δ l (ỹ ki ) = Δ 125 (ỹ ki ) and Δ u (ỹ ki ) = Δ 4876 (ỹ ki ).

Parameters of Equation
One hypothesis of interest was whether there was a relationship between the larval feeding rates and the growth of larvae. To test this, we fit a line to the larval size versus feeding rate observations at the four larval ages around 48 hr-the age our feeding rates were (3) measured. A significant positive slope for these lines was taken as evidence consistent with our hypothesis.
An alternative method of analysis would have been to use the actual observed size of larvae for each selection regime and larval age. However, the advantage of using the predictions is that a much larger sample is utilized to estimate the three-parameter model (Equation 2). The logistic model is an excellent empirical model for this biological process (see Figure 12). Alternatively, using individual estimates of larval size for each selection regime would rest on only five replicates for some selection regimes.

| Adult size
The adult size was collected from the same larvae collected for the larval growth rate assay (Section 2.6). Following Santos et al. (1997), at hour 105, 30 pupae were collected and placed into non-nutritive agar vials to allow for their development. When they had eclosed, FlyNap was used to anesthetize the flies to record their wet weights.
Afterward, they were placed into a drying oven at 80 degrees Celsius for 48 hr. The dry weight of the adult flies was then recorded after the 48 hr had passed. The statistics for analyzing the adult size of flies are the same as the development time differences.

| Larval feeding rates
The RUX, UX, and UTB populations fed at similar rates and are not significantly different from each other (Table 1, Figure 2). A significantly slower feeding rate was observed in the urea-adapted lines-UX, UTB, and RUX-compared to the stress line TDO (p < .0001) and the demographic line CO (p ≤ .0013). RUX and UX fed significantly slower than nCO (p ≤ .024) but UTB did not.
The pooled UX, RUX, and UTB showed slower feeding rates than the AUC controls (p = .014; Figure 3).

| Starvation and desiccation resistance
Two of the selection regimes, TSO and TDO, had undergone selec-

| Viability in urea
We compared the viability (number of larvae surviving to adult stage) of each population in the control environment to the urea environment by computing the difference in the two (Table 4, Figure 6).
If this difference is positive and significantly different from zero, it indicates sensitivity to the toxic effects of urea. The TSO, AUC, CO, and nCO populations show significant sensitivity to urea ( Figure 6).
The viability difference was not significantly different from zero in the urea-adapted populations, RUX, UX, and UTB populations.
The viability of the AUC populations in the control environment was not significantly greater than UX and RUX populations, but was F I G U R E 2 Larval feeding rates relative to the TSO populations measured as sclerite contractions per minute. The horizontal bars mark groups that are not significantly differentiated form each other. Between the separate groups are significant differences, with TDO feeding the fastest and RUX feeding the slowest significantly greater than UTB (p = .022, Figure 7). In the urea environment, AUC viability was significantly less than UX (p < .001), RUX (p < .001), and UTB (p < .001, Figure 7). In the control environment, the TSO selection regime had lower viability than the CO regime (p = .04) but not the nCO regime. In the urea environment, the TSO selection regime had lower viability than the CO regime (p < .0001) and the nCO regime (p < .0001). There were no significant differences between the CO and nCO regimes in either environment ( Figure 8).

| Development time in urea
The developmental time ( In the control environment, the UX development time difference is −24.4 (females) and −23.75 (males). We found that in the urea environment, AUC females (20.5 hr, p < .0001) and males (21.0 hr, p < .0001) had significantly longer development times than their respective UTB sexes, but there were no significant differences between the UX or RUX selection regimes ( Figure 10). In the control environment, the development time of the AUC females was not significantly greater than UTB, RUX, or UX, while AUC male development times were significantly less than RUX (9.6 hr, p = .045), but not different than UX or UTB.
Next, we analyzed the effects of urea on the stress and demographically selected populations. The actual mean development times for the relevant UTB populations, which were used as the standard, were as follows: Females: UTB-urea (block-ii) 301.8 hr, UTB-urea  were no significant differences between the nCO and CO flies.
The pupae that were collected at the 105-hr mark were weighed and dried. A separate collection of pupae was collected for the adult weights. The mean dry weights of each selection regime and the fitted growth equation (Equation 2) are shown in Figure 12. The asymptote seen in all populations between hours 90 and 105 represents the transition of larvae from an actively feeding stage to a wandering stage during which they search for a pupation site. The cessation of feeding thus leads to the plateau in larval size.
The predicted sizes of larvae from Equation (2) were used to estimate the linear relationship between larval size and feeding rate for different aged larvae ( Figure 13). The order of the eight selection regimes in Figure 13 from slowest feeding to fastest was RUX, UX, UTB, AUC, nCO, TSO, CO, and TDO. A significant correlation was seen at hours 48 (p = .004), 54 (p = .002), and 60 (p = .002) but not 42 hr. The 46-60-hr results support the conclusion that a slower feeding rate results in a slower growth rate.

| D ISCUSS I ON
The three selection regimes that have a history of exposure to urea in the larval environment (RUX, UX, and UTB) all show adaption to that environment in the form of increased viability in urea-laced food. They also show reduced feeding rates as has been previously documented (Borash et al., 2000). One difference between RUX, UX, and UTB populations is their development times in urea-laced food. Populations from the UX and UTB selection regimes show reduced development times relative to the AUC control, although only the UTB difference was significant, while populations from the RUX selection regime do not. The RUX populations were kept in control environments for 100 generations prior to these experiments but have continued to retain some of their adaptations to urea. One reasonable explanation for lack of convergence is the fixation of alleles at important loci which affect both survival in urea and larval feeding rates. Upon return to the control environment, the fixed alleles do not revert to the ancestral frequency.
Occasionally selection on one life-stage has resulted in phenotypic differentiation in another life stage. Chippindale et al. (1994) show that adult selection on age-at-reproduction may affect eggto-adult development time. Mery and Kawecki (2003)  CO populations in urea food and normal food is less than the AUC populations ( Figure 6). As noted in the introduction, demographic selection has resulted in increased stress resistance in previous experiments (Service et al., 1985). ins nor extreme competition. In crowded environments, larvae may increase their competitive ability by increasing their feeding rates (Guo et al., 1991;Mueller, 1988aMueller, , 1988bMueller et al., 2005). In addition, increased feeding rates appear to be correlated with increased foraging path lengths which in crowded environments may assist larvae in finding less crowded food resources (Mueller et al., 2005;Sokolowski et al., 1997). We note that there are complications to this simple view of the effects of crowding on feeding rates, but this does represent one possible outcome of evolution (Borash et al., 1998;Nagarajan et al., 2016). None of the populations used in this study have been subject to crowding so hence this particular selective pressure is not believed to have impacted any of the study populations.
Environments with toxins, like urea or ammonia, result in declining feeding rates (Borash et al., 2000;Mueller et al., 2005). Reduced feeding rates will reduce the intake of toxic compounds which the larva must then detoxify. Reduced feeding rates also increase the efficiency of energy extraction (Mueller, 1990), thus facilitating the energy-consuming process of detoxification (Mueller & Barter, 2015).
In both environments, crowding (Mueller & Joshi, 2000, figure 6.30) and toxins (Borash et al., 2000), we see a rapid change in feeding rates, at least in laboratory conditions, suggesting strong selection for these phenotypic changes in feeding rates. Of course, the RUX, UX, and UTB populations have all been subject to the selective pressures of urea detoxification and as seen in this study and elsewhere have all evolved slower feeding rates.
The third environment, with no or little competition and no toxins, is the one experienced by the AUC, nCO, CO, TSO, and TDO Our results have shown that larvae between 48 and 60 hr show a positive relationship between feeding rates and larval dry weight although the result at 42 hr is not significant. Thus, at 60 hr an increase in feeding rate of 10 retractions per minute will increase dry weight from 5.9% (at the highest feeding rate) to 8.1% (at the lowest feeding rate). The feeding rates are measured using 48-hr-old larvae. This conclusion is different from Kaun et al. (2007) who found no statistically significant correlation between feeding rates and food ingestion. There are three possible explanations for the difference between Kaun et al.'s results and our own listed as follows: (i) no physiological relationship exists between feeding rates and ingestion or growth-the observed statistical relationship is artifactual-(ii) a positive relationship exists but the Kaun et al. experiments lacked the power to detect that relationship, or (iii) no relationship exists between feeding rates and ingestion but there is a positive correlation between feeding rates and larval growth rates.
The third explanation is very difficult to explain especially since it stands to reason that faster feeding larvae would require more energy for this activity. So, we exclude explanation (iii) as an unreasonable alternative. We believe the second explanation is correct which we explain next.
Kaun et al. studied a total of 73 individuals whose feeding rate varied from 48 to 75 retractions per minute. Their data showed a positive relationship but with a p-value of only .11. In the present study, feeding rates of 1,375 individuals were measured with a range of 10 to 169. Selection regimes used in our study had mean feeding rates that varied from 79 to 141. Thus, not only was our sample size larger by a factor of almost 20, but also the range of feeding rates was much larger. Thus, we believe our analysis has much greater power to detect a correlation between feeding rates and growth due to both the increased sample size and greater range of feeding rates. To test this idea, we generated computer samples to recreate Figure 13. We restricted our database to feeding rates between 114 and 141. This included all selection regimes except RUX. We sampled only 73 individuals from all populations except TSO to estimate mean feeding rates on the x-axis of Figure 13. We used all of the growth rate data to create the same y-axis values as used in Figure 13. The TSO population was always added to these computer-generated data since TSO was the standard. For each of the four time intervals shown in Figure 13 (42, 48, 54, and 60 hr), we created 100 independent samples to test for significant slopes.
The fraction of 100 tests that were significant was 8%, 0%, 0%, and 0% at hours 42, 48, 54, and 60 respectively. These results support our contention that small sample sizes and a limited range of feeding rates would obscure the trend we observed.
This study revealed extensive variation among the forty populations for larval feeding rates, viability and development time in urea, and starvation resistance. However, there was little variation for desiccation resistance. Computer simulations have shown that gene-phenotype relationships are more likely to be uncovered for phenotypes with large levels of genetic variation relative to environmental variation (Mueller et al., 2018). Thus, except for desiccation resistance, we believe genomic studies of these populations should yield insights into gene-phenotype connections.

| CON CLUS IONS
Populations subject to demographic selection show elevated resistance to starvation relative to populations selected for urea resistance. These demographic populations may also show partial resistance to urea in their larval food. Selection for urea resistance does not enhance adult stress resistance against starvation or desiccation. The results of this study provide convincing evidence of the relationship between larval feeding rates and larval growth rates.
This relationship suggests that energy allocation to activities like growth and detoxification of urea are central to the evolution of feeding rates as suggested by Mueller and Barter (2015). In addition, they point to a way of understanding the central role feeding rates have played in the evolution of Drosophila larvae to a wide array of challenging environments.

ACK N OWLED G M ENTS
We thank Professor Michael Rose for use of his selection lines and valuable comments and two referees for valuable comments. This research was supported by funds from the School of Biological Sciences, University of California, Irvine.

F I G U R E 1 3
The predicted individual larval dry weight, with 95% confidence intervals, as a function of feeding rate. The p-values are for tests of the hypothesis that the slope of the line is zero. The order of the eight selection regimes along the x-axis is from slowest feeding to fastest: RUX, UX, UTB, AUC, nCO, TSO, CO, and TDO. The x-axis coordinates are offset to improve readability

CO N FLI C T O F I NTE R E S T
None declared. Writing-review & editing (lead).

DATA AVA I L A B I L I T Y S TAT E M E N T
The data can be accessed at: https://doi.org/10.7280/D1VT3G.